2026-10-20 | | Total: 55
Strong thermal gradients in tungsten plasma-facing components can bias the migration of irradiation-induced defects and thereby alter their redistribution, interaction, and retention. Here, we extend an atomistic jump-statistics framework previously developed for hydrogen to quantify the thermomigration of a single vacancy and interstitial He$_1$--He$_4$ clusters in body-centered-cubic tungsten. Non-equilibrium molecular-dynamics simulations performed over several thermal gradients are used to determine the directional jump ratio, while geometrically derived projected jump lengths convert the fitted directional-bias coefficient into the heat-of-transport coefficient $Q^*=αkT^2$, where $k$ is the Boltzmann constant, $T$ the temperature, and $α$ a constant. All investigated defects exhibit negative values of $α$, corresponding to net migration toward the hot region. He$_4$ has the largest response, $α=-2.2457\times10^{-2}$ K$^{-1}$, followed closely by the vacancy, whereas He$_3$ exhibits the weakest response despite its intermediate cluster size. Independent constant-temperature simulations over 1000--1600 K show that He$_3$ is also the least mobile cluster and has the largest effective migration energy, $E_m^{\mathrm{eff}}=0.271$ eV. In contrast, He$_4$ is more mobile than He$_3$ and has the smallest effective migration energy, 0.087 eV. Agreement between the He$_3$ concentration profile predicted from the jump statistics and that obtained from a long molecular-dynamics trajectory provides an independent validation of the method. These results show that thermomigration cannot be inferred from cluster size, diffusivity, or migration energy alone and provide defect-specific heat-of-transport parameters for higher-length-scale models of tungsten under fusion-relevant thermal gradients.
Piezoelectric nanogenerators (PENGs) are an attractive platform for self-powered mechanical-energy harvesting, but their output is often limited by depolarization fields, charge leakage, and inefficient stress transfer to the active dipoles, particularly in polymer-dominated composites. Here we report HoFeO3-BaTiO3 (HFO-BTO) nanofiber heterostructures prepared by alternating sol-gel electrospinning, in which orthorhombic HoFeO3 (Pbnm) and tetragonal BaTiO3 (P4mm) coexist in intimate contact. Rietveld-refined X-ray diffraction, high-resolution transmission electron microscopy, and Raman spectroscopy indicate interfacial lattice distortion. Polarization-electric field loops show that the heterostructure retains a polarization comparable to that of BaTiO3 (16.39 vs 15.80 microcoulombs per square centimeter at zero field) while its switching-current maximum moves to a higher field (about 7 vs about 3.5 kilovolts per centimeter), consistent with strain-stabilized ferroelectric domains. When the nanofibers are encapsulated in polydimethylsiloxane (PDMS) and excited under cyclic impact, the device produces a hybrid output: the piezoelectric response of the strained BaTiO3 is supplemented by a triboelectric contribution from contact electrification. The HFO-BTO device delivers open-circuit voltage peaks of about +6 volts and -8 volts and a power density of 24.31 microwatts per cubic centimeter at a 4.7 megaohm load. This is 2.4 times the value for an identically fabricated BaTiO3 device (10.00 microwatts per cubic centimeter) and 8.7 times that of the HoFeO3 device (2.78 microwatts per cubic centimeter). These results link interface-driven lattice distortion to improved electromechanical output and identify all-oxide nanofiber heterostructures as a promising route to lead-free energy harvesting.
Plasma-facing components (PFCs) in future fusion reactors must withstand extreme environments, making tungsten (W) a prime candidate material due to its high thermal conductivity, high melting temperature, low sputtering yield, low tritium absorption. However, neutron irradiation introduces complex nanostructural defects that severely degrade these critical properties. In this work, we use Non-Equilibrium Molecular Dynamics (NEMD) simulations coupled with the Creation Relaxation Algorithm (CRA) to investigate the lattice thermal conductivity (kappa) of W across a wide range of irradiation doses. Systems treated with different doses up to 1.8 canonical displacements per atom (cdpa) were simulated to generate and study diverse defect structures, including vacancies and self-interstitial atom (SIA) clusters. By correlating spatially-resolved temperature gradients with local defect concentrations, we show that even small changes in defect type and distribution significantly affect thermal transport. We employ Machine Learning (ML) models to develop a quantitative relationship between the defect features and the resulting thermal conductivity degradation. The ML analysis identifies vacancy density as the dominant structural factor controlling thermal transport. We found a significant, nonlinear suppression of kappa with increasing defect density. Predictive multi-defect mapping built upon our ML model resolved a step-like degradation transition at a critical vacancy threshold of 0.05 at.%, beyond which phonon scattering completely saturated and capped the maximum relative lattice thermal conductivity degradation at approximately 30 at.%. This integrated simulation and ML approach provides a powerful, predictive framework for quantifying the complex relationship between structure and thermal conductivity in irradiated fusion materials.
Leveraging our earlier formalism for charge-state decay [see R. Kuate Defo et al., 108 235208 (2023)], we investigate relative population partitioning in the presence of singly negatively charged substitutional boron impurities (B$^-_{\rm C}$) in diamond samples, where these diamond samples also contain substitutional nitrogen (N$_{\rm C }$) impurities and have ether-terminated surfaces. Effectively, given any two defects, one donor and one acceptor, such that their charge-transition levels, their relative populations, and the temperature are known, the timescale for charge transfer between the two defects is completely determined. We use this insight and experimental measurements for the timescale for charge transfer between NV$^-$ centers and N$_{\rm C}$ hybridized with the ether-terminated surface of diamond, as well as the temperature and the relative populations for the NV$^-$ centers, [see Z. Yuan et al., Phys. Rev. Res. 2, 033263 (2020)] to argue that an additional defect must be modulating the relative populations of N$_{\rm C}$ hybridized with the ether-terminated surface of diamond. The shifts in the local Fermi level for N$_{\rm C}$ hybridized with the ether-terminated surface of diamond are consistent with the additional defect being a substitutional boron impurity (B$_{\rm C}$). Ultimately, our results allow for the determination of local relative concentrations of ionized defects.
Metal nitridocarbonates containing the (NCN)2- ion are established for a multitude of metal atoms throughout the periodic table. The preparation of nitridocarbonate compounds containing a (CN3)5-, (C2N4)6-, or (C3N6)6- ion remains challenging, although corresponding species have been established. The new compound La14.33(CN3)8N was synthesized in a conventional solid-state metathesis reaction from a mixture of LaF3, Na2(CN2) and Li3N in an open system at 550 °C. The reaction conditions, including the formation and de-composition of La14.33(CN3)8N are analyzed on basis of thermoanalytic studies. The crystal structure of La14.33(CN3)8N was refined by single-crystal X-ray diffraction, and the mid-IR frequencies resulting from two crystallographically distinct (CN3)5- ions in the structure are reported. 23Na, 7Li and 13C NMR studies provide additional support for the identity of the compound. The crystal structure of La14.33(CN3)8N is isotypic to that of Ce15(BN3)8N.
Diamond is an important host material for functional defects such as solid-state spin qubits. During the defect creation by irradiation and subsequent thermal annealing, carbon-related defects, including vacancies and interstitials, are generated and can migrate to form new defect complexes. To investigate this kinetically driven process and reveal the formation mechanisms of the di-interstitials and di-vacancies in the neutral charge state, we fine-tuned a universal machine learning interatomic potential using a small dataset generated from density functional theory calculations. We obtained a barrier of 1.84 eV for mono-interstitial migration and found that two initially separated mono-interstitials only need to overcome a lower barrier to form di-interstitial configurations. Beyond the training data, the MLIP also identified the formation pathway for an additional di-interstitial configuration and a stepwise divacancy dissociation mechanism. These results reveal migration-assisted defect complex formation mechanisms and demonstrate how targeted MLIP training can identify reaction pathways.
Extensive studies of hidden spin polarization have established a broad class of hidden phenomena in quantum materials. Here, we introduce the hidden piezo-Hall effect, in which strain induces finite anomalous Hall conductivities (AHCs) in individual sectors while the global AHC remains zero. At zero strain, neither the full system nor its sectors exhibit an anomalous Hall effect (AHE). Under uniaxial strain, however, opposite sector AHCs emerge while their sum remains canceled. We establish this effect for a $d$-wave hidden altermagnet using a tight-binding model and first-principles calculations for the experimentally synthesized material $\mathrm{RbV_2Te_2O}$. Rotating the strain direction from the $x$ to the $y$ axis reverses the signs of the sector AHCs without generating a net Hall response. The same symmetry mechanism may also give rise to hidden piezomagneto-optical, piezo-thermal Hall, piezo-Nernst, and piezo-Ettingshausen effects. Our work extends hidden physics to strain-controlled transport and identifies $\mathrm{RbV_2Te_2O}$ as a promising candidate for experimental verification.
GaN-on-diamond heterostructures are promising for high-power and high-frequency electronics, however, direct epi- taxial growth on diamond (100) remains challenging due to severe lattice and thermal mismatches, which induce high threading dislocation densities. To overcome this challenge, we have used Plasma Assisted Molecular Beam Epitaxy (PAMBE) to grow wurtzite GaN layers on diamond(100) substrates using nanowire-assisted epitaxial lateral overgrowth (ELOG). Two nanowire templates (NW100, NW600) were benchmarked against a nanowire-free reference (NW000). High-resolution XRD, φ -scans, (0002) pole figures, and reciprocal space mapping show that NW100 attains the lowest dislocation density (5.51 x 10-4 nm-2) and microstrain (2.736 x 10-3), the narrowest (0002) reflection, and a c-axis lattice parameter (5.1821 Å) within 0.06% of the bulk value; NW000, in contrast, exhibits a contracted c-axis (5.1633 Å) and elevated dislocation density (6.26 x 10-4 nm-2), consistent with residual tensile stress. Cross-sectional TEM con- firms the intended diamond/AlN/nanowire/ELOG/planar-GaN layer sequence and shows threading dislocations bending into near-lateral trajectories, directly evidencing defect filtering by the nanowire template. Together, these results estab- lish nanowire-assisted ELOG as an effective strategy to suppress dislocation propagation and recover near-bulk lattice parameters in GaN-on-diamond heteroepitaxy, providing a structurally well-defined platform for device integration.
It is commonly believed that a bulk photovoltaic current requires broken spatial inversion symmetry. However, in this work, we demonstrate that dc photocurrents can arise in centrosymmetric solids from the bilinear coupling between the electric and magnetic components of light. This magnetoelectric bulk photovoltaic effect (MBPE) has a quantum-geometric origin in the dipole moments of the Zeeman quantum metric and Zeeman Berry curvature. We identify four distinct MBPE mechanisms, including time-reversal-even and time-reversal-odd injection-like and shift-like photocurrents driven by linearly or circularly polarized light. Our symmetry analysis shows that the MBPE is allowed in broad classes of two-dimensional centrosymmetric nonmagnetic and magnetic materials under normally incident light. As representative examples, we predict sizable MBPEs in centrosymmetric monolayer Pt$_2$HgSe$_3$ and altermagnetic RuF$_4$, with current densities on the order of $1~\mathrm{mA/m}$. Notably, the injection-like photocurrent in RuF$_4$ reverses sign upon switching either the light helicity or the Néel vector. Our results extend the bulk photovoltaic effect to centrosymmetric solids and open new opportunities for electrical readout of altermagnetic order.
Ferrimagnetic oxides such as rare-earth iron garnets (REIG) have emerged as essential components of spin-based devices owing to their high Curie temperatures, low Gilbert damping, and excellent insulating properties. Recently, REIGs have attracted immense interest for hosting robust chiral spin textures that can be potential information carriers in next-generation spintronic devices. Here we investigate the structural, chemical, and magnetic behaviour of single-crystal Lu3Fe5O12 (LuIG) thin films for different film thicknesses on a Gd3(Mg,Zr)Ga5O12(111) or SGGG(111) substrate. Our results demonstrate the presence of a 3 nm-thick Gd-rich LuIG layer at the LuIG/SGGG interface, which gives rise to a magnetic compensation temperature around 50 K. We further probe the anomalous Hall effect (AHE) in LuIG layers with an adjacent platinum ultrathin film. In addition to the magnetic compensation point, the AHE loops show pronounced nonmonotonic features, often attributed to the spin Hall topological-like Hall effect (SH-THE) and non-collinear spin textures arising from the Dzyaloshinskii-Moriya interaction (DMI). LuIG thickness-dependent AHE measurements reveal that the composition variation due to Gd diffusion breaks the inversion symmetry and could contribute to DMI. These findings shed light on fundamental understanding and engineering of the different interfaces of magnetic garnets for advanced spintronic device applications.
The structural behaviour of cubic (Fm-3m) K2OsCl6 under compression was investigated to 35.4 GPa using angle-dispersive synchrotron powder X-ray diffraction. Upon increasing pressure, K2OsCl6 undergoes a sequence of structural transformations from the ambient-pressure cubic phase to a tetragonal (P4/mnc) phase above 9.5(3) GPa, followed by a transition to a monoclinic (P21/n) phase beginning at approximately 14.2 GPa. The tetragonal and monoclinic phases coexist over the pressure range 14.2-17.2 GPa, while the monoclinic phase is observed as the sole phase at 17.9 GPa and remains stable up to the maximum pressure investigated. Analysis of the diffraction profiles and unit-cell parameters indicates that the intermediate tetragonal structure is likely associated with stress-induced symmetry lowering rather than a distinct thermodynamic phase. This distortion, which is adequately described by the tetragonal P4/mnc space group, is attributed to non-hydrostatic deviatoric stresses arising from the pressure-transmitting medium. Equation-of-state analysis reveals a pronounced reduction in compressibility across the structural transitions, with the bulk modulus increasing from 26(2) GPa for the cubic phase to 71(3) GPa for the monoclinic phase. The pressure-induced structural sequence is analogous to that reported for related compounds such as K2IrCl6, although the transitions in K2OsCl6 occur at lower pressures. All structural transformations are reversible upon decompression.
AI/ML models are useful for the prediction of materials' properties. Recently, there has been great progress in developing Neural Network models to predict the properties of solids, molecules, and crystals, avoiding expensive and time-consuming laboratory testing. This study presents a novel AI/ML-based multi-modal approach to predict material magnetic ordering using material crystal structure information as input to a Graph Neural Network model (GNN) and semantic text descriptions of the material's properties and structure as input to a Large Language Model (LLM), before performing information fusion of the GNN and LLM outputs to predict material magnetic ordering in four classes (Non-magnetic, Ferromagnetic, Antiferromagnetic, Ferrimagnetic). The model is trained, validated, and tested on a dataset which is a combination of Materials Project dataset, Magndata dataset, and a set of uranium compounds with magnetic ordering determined in laboratory testing, with an 80/10/10 (training/validation/testing) data split. The proposed GNN and LLM multi-modal approach achieves better overall accuracy than both the standalone GNN and LLM approaches with three fusion strategies.
Self-driving laboratories accelerate materials exploration through autonomous closed-loop integration of synthesis, characterization, and next-experiment selection. However, characterization does not always directly provide the target property; measured spectra or images often require automated analysis to extract quantitative descriptors based on physically motivated assumptions and predefined criteria. Such procedures can fail for measurements outside these criteria, requiring a decision on whether a reliable value can be recovered or no numerical value should be reported. Here, we implemented an autonomous adaptive spectral-analysis agent using a large language model (LLM) to supervise a deterministic numerical-analysis routine. For each spectrum, the agent executes the predefined analysis script and assesses the result using the measured spectrum and fitting diagnostics. Valid results are accepted unchanged, recoverable failures are reanalyzed through bounded adjustments of predefined fitting parameters, and NaN (not a number) is returned when no scientifically defensible value can be obtained. We evaluated the workflow on optical transmittance spectra from 203 epitaxial beta-Ga2O3 films grown during a previously reported self-driving thin-film growth campaign. The predefined optical analysis script returned Urbach energy, an optical-quality metric, for 183 spectra and failed for 20. The agent retained those 183 values and evaluated 20 exception cases, recovering 16 through bounded parameter adjustments and returning NaN for 4 when the spectral evidence did not support a reliable Urbach tail. The audit records captured spectral features, parameter changes, and rationale for each decision. These results demonstrate autonomous scientific exception handling by an LLM agent while leaving quantitative property extraction to deterministic numerical analysis.
A weak or vanishing circular optical signal need not imply weak helicity-sensitive transitions: it can be the residual of much larger opposing microscopic contributions. We resolve the finite-frequency interband Kubo response before summing over transitions and momentum, separating cancellation among transitions at fixed momentum from cancellation across the Brillouin zone. An exactly compensated $d$-wave model establishes the symmetry benchmark. In MnTe, first-principles Wannier calculations reveal a symmetry-dark Néel orientation with an extended sign-changing circular kernel whose bulk integral vanishes. Reorienting the Néel vector removes this extinction, yet the allowed response remains a small imbalance between large positive- and negative-momentum sectors; the imbalance reverses with photon energy and flips the macroscopic sign. Thus magnetic symmetry decides whether circular response is allowed, whereas microscopic compensation determines how much survives and with what sign. This distinction provides a general framework for interpreting weak, vanishing, and sign-changing magneto-optical signals in compensated magnets.
We present NEP-kappa 2.0, a unified framework that integrates density-functional theory, material-specific machine-learning interatomic potentials, and pretrained force models for advanced phonon calculations. A common workflow enables harmonic lattice dynamics, quasi-harmonic and stochastic self-consistent phonon treatments, three- and four-phonon transport, Wigner population--coherence transport, and time-dependent phonon dynamics. Applications to 3C-SiC, BAs, and Cs$_2$AgInCl$_6$ demonstrate consistent force-model comparisons, finite-temperature renormalization, higher-order scattering, and intermode coherence. Beyond steady-state transport, mode-resolved simulations of selectively excited 3C-SiC reveal that half of the injected optical-phonon energy reaches acoustic modes within approximately 3.75 ps, primarily through one-optical--two-acoustic scattering. Acoustic branches exhibit response times differing by a factor of five, with competing scattering channels producing a transient energy overshoot. Channel-resolved entropy production and persistent nonthermal occupations distinguish interbranch energy redistribution from intrabranch thermalization. The framework connects equilibrium phonon transport with microscopic nonequilibrium energy dynamics.
Five metals share a common sublattice in the (Hf,Ta,Ti,Nb,Zr)B$_2$ (HEB$_2$). Most high-entropy materials are produced as polycrystalline materials in the form of sintered ceramics. The question of which metal species are located at grain boundaries is still unknown. Here, we calculated the grain-boundary energy $γ_{GB}$ of four types of boundaries in five binary diborides and HEB$_2$ using density functional theory (DFT). We analyzed the bonding at these boundaries using COHP and sampled cation segregation at six boundaries of HEB$_2$, ordered and general, using Monte Carlo / molecular dynamics with a MACE potential that was fine-tuned for this system. According to our results $γ_{GB}$ splits by group: the group V diborides NbB$_2$ and TaB$_2$ form boundaries with about half the energy of the group IV at three of the four types of GBs, ranging from 0.8 to 2.7 compared to 2.5--4.8~J/m$^2$. COHP shows a significant difference at the Fermi level, namely the Me--B bonds of NbB$_2$ and TaB$_2$ are antibonding at Fermi energy ($E_F$) in the bulk but bonding at the GB, while their integrated strength hardly changes. In HEB$_2$ the segregation follows a similar pattern: Nb and Ta are enriched at five of the six grain boundaries, including general boundaries, up to 23--35~at.\% at 2000~K, although the size of the atoms would favor Zr and Hf in the expanded sites. A segregated $Σ193$ grain boundary can withstand 26 against 20~GPa for a random boundary in tension at 300~K. The grain boundaries of HEB$_2$ are composed of group V metals, as suggested by the boundary energies of binaries compounds and their bonding at the Fermi level. Grain boundaries of the (Hf,Ta,Ti,Nb,Zr)B$_2$ preferentially include Nb and Ta. This lowers the grain boundary energy and strengthens the boundaries. This consists of the lower boundary energy and changes Fermi level bonding in NbB$_2$ and TaB$_2$.
The ability to create and control nonequilibrium spin states is central to understanding ultrafast magnetism and to developing optical routes for manipulating magnetic materials beyond conventional thermal pathways.In metallic systems, photoexcited electronic distributions rapidly evolve toward quasi-thermal states through efficient electron-electron and electron-phonon scattering, making long-lived modification of spin-dependent interactions difficult to achieve. Here, we investigate coherent magnetization precession in an amorphous GdFeCo thin film driven by simultaneous excitation with two femtosecond laser pulses.Using a symmetric dual-pump geometry that continuously evolves from single-pulse to spatially periodic excitation, we separate effects associated with the excitation pattern from those arising from simultaneous dual-pump excitation. We identify a narrow excitation window in which the lifetime of the ferromagnetic-resonance-like mode increases strongly while its frequency exhibits a small but reproducible increase.Time-window analysis shows that the enhanced lifetime persists for several hundred picoseconds, substantially longer than transient magnetostatic effects associated with the periodic excitation pattern.The combined frequency and decay time changes cannot be explained by spatially structured heating alone and instead point to a long-lived nonequilibrium spin-dependent electronic configuration. Within a phenomenological Elliott-Yafet framework, the simultaneous increase in FMR decay time and frequency is consistent with a transient reduction of the effective spin-mixing parameter governing spin-flip scattering. Our results demonstrate that the excitation pathway, rather than excitation energy alone, can serve as a control parameter for spin relaxation and provide access to nonequilibrium spin states that are not readily reached by conventional single-pulse excitation.
Recently, $\text{EuCd}_2\text{P}_2$ has attracted increased attention for its colossal magnetoresistive (CMR) property, the mechanism of which seems to deviate from that of known canonical CMR compounds such as manganites. Many studies underscore the role of low density of hole carriers caused by Eu vacancies, particularly in stabilizing ferromagnetic (FM) order and modifying electrical transport and magnetoresistive responses in this compound. A magnetic polaron based mechanism rather than double exchange based mechanism is being proposed to explain the CMR effect for this compound. In this work, we explore the effect of elctron doping on the magnetic transition and CMR behaviour via Gd substitution at Eu sites with varying concentration (x) in Eu$_{1-x}$Gd$_{x}$Cd$_2$P$_2$. Our Density Functional Theory (DFT) and Monte Carlo-based calculations reveal an oscillating ground state magnetic order as a function of electron carrier density. The ground state is found to be A-AFM metal for x=0.056, 0.125 and 0.5 while it becomes ferromagnetic metal at x= 0.33, 0.67 and 1. Most importantly, we predict colossal anisotropic magnetoresistance in this compound for x=0.125 and 0.5. The underlying mechanism of anisotropic CMR in this compound looks similar to double exchange like due to the presence of a strong intra-site 4\textit{f}-5\textit{d} ferromagnetic exchange interaction and an itinerant electron in the hybridized conduction states. The oscillating ground state magnetic order and the enhancement of transition temperature (as high as 150K for x=0.67) as a function of electron carrier density offers a new avenue to tune the magnetotransport property in this compound and offers an excellent platform for designing future spintronics devices.
Electric hidden unconventional collinear magnetism, in which a perpendicular electric field (\(\mathcal{E}_z\)) lifts the \(\mathcal{PT}\)-enforced degeneracy and selectively releases the hidden crystallographic constraints in conventional 2D antiferromagnets to drive emergent magnetic phases, has remained largely unexplored. Here we establish a complete symmetry classification of all 88 collinear spin layer groups (cSLGs) compatible with conventional \(\mathcal{PT}\)-symmetric antiferromagnets, and show that the \(\mathcal{E}_z\) drives the system into three distinct classes of unconventional collinear magnetism: type-IV magnetism (Class-I, 23 cSLGs), altermagnetism (Class-II, 34 cSLGs), and fully compensated ferrimagnetism (Class-III, 31 cSLGs). We illustrate these categories using the most common and feasibly exfoliable 2D candidate materials. These include Sr$_2$CaCr$_2$P$_2$O$_9$ and Sr$_2$CaCr$_2$As$_2$O$_9$ (Class-I), BaMnGe and MnPS$_3$ (Class-II), and Sr$_3$Fe$_2$O$_5$Cl$_2$ and bilayer MnBi$_2$Te$_4$ (Class-III). We further demonstrate that the \(\mathcal{E}_z\) generates distinctly different spin-splitting fingerprints across the three classes, yielding distinct spin-dependent responses: the electric Hall effect in Class-I, spin-polarized transport in Class-II, and the electric reversible magneto-optical Kerr effect in Class-III. Our work not only closes a complete internal cycle from conventional to unconventional collinear magnetism in 2D antiferromagnets, but also establishes a viable all-electrical compensated magnetism platform and control scheme for future spintronic applications.
The extrinsic bulk photovoltaic effect from asymmetric electron-defect scattering can make an important contribution to the total photocurrent and is highly diagnostic of defect type. We present an ab initio Boltzmann framework combining optical excitation with electron-phonon and symmetric and skew electron-defect scattering, and we apply it to monolayer GeS. Because the defects that generate the current also relax the anisotropy that sustains it, the response peaks at a finite defect fraction. Within the scattering-limited model, cooling to $0$~K shifts the peak position toward lower defect fractions and yields a maximum enhancement of more than six orders of magnitude at fixed defect fraction. At selected low-temperature defect fractions, the total extrinsic current can exceed the intrinsic shift current by more than an order of magnitude, so that the defect fraction itself becomes a tunable parameter.
This article describes the application of HEXRD, an open-source Python library for X-ray diffraction analysis, and HEXRDGUI, its companion graphical application, to dynamic compression experiments. Built on a generic multi-detector geometric model, the toolkit supports arbitrary instrument configurations with independently positioned planar or cylindrical detector panels, arbitrary beam directions, and multiple X-ray sources. It calibrates instruments from powder rings, Laue spots or the rings of an unknown phase, fits whole patterns by Le Bail or Rietveld refinement, and imports data from the TARDIS, FIDDLE and PXRDIP platforms at the National Ignition Facility (NIF) and OMEGA EP. Although originally developed for high-energy diffraction microscopy at synchrotron sources, HEXRD has been adopted as an analysis tool for dynamic compression X-ray diffraction experiments at facilities including the NIF, the Matter in Extreme Conditions end-station at the Linac Coherent Light Source, the Dynamic Compression Sector (DCS) at the Advanced Photon Source, OMEGA EP at the Laboratory for Laser Energetics, and the High Energy Density (HED) instrument at the European XFEL. Calibrated from the same ring positions on a single planar detector, HEXRD and the independent code pyFAI predict the same scattering angle 2$θ$ to numerical precision. On a dual-source dataset from the NIF, calibrating the detectors against both X-ray sources at once fits the rings of both to an RMS of 0.18° in 2$θ$, about one image-plate pixel, whereas a fit to either source alone, with the plate tilts held at their nominal values, fits its own rings nearly as well but mispredicts the other's by about 0.6°. The software is freely available under a BSD license.
The thermal conductivity of natural diamond can be further improved by ~50% by eliminating isotope impurity to obtain pure 12C diamond. However, unlike isolated point defects, how color centers scatter phonons and correspondingly limit the thermal conductivity remain poorly understood. Here, we report the thermal conductivities of two isotope-purified diamond with different Nitrogen-vacancy (NV) color center concentrations, confirmed by photoluminescence (PL) and secondary-ion mass spectrometry (SIMS). Results show ~ 3200 W m-1 K-1 and ~ 2100 W m-1 K-1 for two samples of low and high NV concentrations, respectively. Modified Callaway model analyses show that first-principles-derived scattering parameters can reproduce room-temperature thermal conductivity of the high NV concentration sample under a vacancy concentration assumption (~30% of the total Nitrogen concentration), but fail to reproduce its temperature dependence. Compared to PL and SIMS results, which show ~7.1-fold NV concentration and ~92-fold difference in total Nitrogen concentration between two samples, suggesting the assumed 30% vacancy fraction may be greatly overestimated. Taken together, these results suggest that NV-color-center defects may scatter phonons more strongly than models treating Nitrogen and vacancy as independent point defects, providing new insights into color-center-induced phonon-scattering mechanisms.
Thermal conductivity ($κ$) measurement of a thin film is critical for developing thermal management materials in heterogeneous systems. A nanoscale point contact method based on scanning thermal microscopy (SThM) offers a unique approach to probe $κ$ with nanoscale spatial resolution in sub-100 nm thick films. Extracting $κ$ of a thin film from this method requires separating film, substrate, and interface contributions that are comparable in magnitude. Multiple unknowns make this extraction a high dimensional problem. In this work, we used dimensional analysis to reduce the number of variables to provide physical insights into the realistic bounds of $κ$ measurement. We show that the total tip-sample thermal resistance $R_{ss}$, originally a function of six variables, reduces to four independent dimensionless groups. These groups are used to construct a mathematical relation that captures the dependence of $R_{ss}$ on the governing parameters. The resulting analysis shows that the point contact tip-film-substrate problem is governed by a single dimensionless invariant, $Ω= R_{ss}κ_{\mathrm{sub}}t$, and two control parameters: the confinement ratio $t/b$ and the thermal contrast $κ_{\mathrm{sub}}/κ_{\mathrm{film}}$. From the invariant, we derive an amplification factor, set only by the contact geometry and the interfaces, that converts the uncertainty of a measurement into an uncertainty in the extracted conductivity. This factor scales with the ratio of the parasitic interface resistance to the film's sheet thermal resistance and bounds the film conductivity that a given contact can resolve. We apply the framework to a 14$\pm$1 nm Al film grown on SiO$_2$ to locate the measurement within the established bound. The framework then lets us map its thermal resistance with sub-100 nm spatial resolution and recover the local thermal conductivity.
Pretrained diffusion models provide powerful priors for generating new materials, including inorganic crystals, but generating useful materials requires satisfying specific chemical and physical properties, from mechanical stiffness to bandgaps. Existing approaches either retrain models for new target properties or steer only the atom coordinates and lattice during generation, leaving composition unchanged despite its fundamental role in determining physical properties. We present FK-ReMask, an inference-time, gradient-free method that steers both the continuous and discrete channels of a pretrained crystal diffusion model. During generation, FK-ReMask maintains multiple denoising trajectories and replicates promising candidates using Feynman-Kac particle selection. However, selection alone is insufficient: in the atom-type channel, committed elements cannot be revised, so cloning propagates existing compositions and limits the population's ability to explore new ones. We address this with an atom re-masking mechanism, which reopens committed sites during the generation process and allows alternative compositions to emerge. Across four electronic, photovoltaic, and mechanical objectives, our method improves the yield of SUN (stable, unique, novel) and property-on-target materials. We further demonstrate joint steering of three photovoltaic properties, showing that multiple physical objectives can be specified entirely at inference time without retraining.
In this paper, we introduce an Abaqus user-subroutine implementation of the phase-field model of fracture initiated by Kumar, Francfort, and Lopez-Pamies (\emph{J. Mech. Phys. Solids} 112, 523--551, 2018). The implementation is presented for the basic case of isotropic linear elastic brittle materials with a Drucker-Prager strength surface. The material inputs are therefore the Young's modulus $E$, the Poisson's ratio $ν$, the uniaxial tensile strength $\sts$, the uniaxial compressive strength $\scs$, and the fracture toughness $G_c$. The mathematical correspondence between the heat equation and the evolution equation for the phase field $v(\bfX,t)$ is exploited through a UMAT subroutine, wherein the driving force $c_{\texttt{e}}(\bfX,t)$ --- accounting for the material strength surface --- in the phase-field model plays the role of a heat source. To ensure that the equations for the displacement field $\bfu(\bfX,t)$ and the phase field $v(\bfX,t)$ are solved in a staggered manner --- a solution strategy that, as opposed to monolithic approaches, reliably yields the correct physical solution --- we make use of a combination of UEXTERNALDB, UAMP, MPC, and URDFIL subroutines. Within this approach, the MPC and URDFIL subroutines are additionally utilized to enforce the required bounding and irreversibility inequalities on $v(\bfX,t)$ via a primal-dual active set approach. The capabilities and robustness of the proposed implementation are demonstrated by simulating the ``Nine Circles of Elastic Brittle Fracture'' on a silicate glass --- a set of tests that span the entire range of well-settled experimental knowledge on fracture nucleation and propagation. These simulations also serve to illustrate the deployment of the code using several native Abaqus finite elements, including the simplicial elements CPS3, CAX3, and C3D4.