SEARCH · Search NASA

Results for “Diffusion Couple”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Modeling thermal radiation waves in silica plasmas for the Mooncat NIF experiment

The Mooncat experiment on the National Ignition Facility uses a laser-driven hohlraum to create a thermal radiation wave in a titanium-doped silica plasma. The titanium dopant enables absorption spectroscopy measurements to infer the temperature of the wave as it propagates. This measurement can be used to constrain multi-physics simulation codes to better understand when simulations do not match an experiment. In this paper, we present radiation-hydrodynamics simulations of the thermal radiation wave in the first full-platform shots of the Mooncat experiment. We examine the important parameters of the simulation, focusing on the radiation temperature source, the material model of the silica plasma as it pertains to radiation transport, and lateral leakage through a beryllium tube enclosing the silica. We compare different simulation modeling strategies to an analytic model of diffusive radiation transport and find that the simulation agrees with the analytic model when it is sufficiently simplified. These simulations show how radiation energy couples to matter to develop a shock wave in a radiative heat wave, an important topic in astrophysics and nuclear fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Preliminary Plan to Inform Testing of a Heat Exchanger Test Article

This report presents a preliminary plan to guide the qualification testing of advanced heat exchanger (HX) components for nuclear-to-industrial heat transfer applications. The objective is to establish a defensible, physics-based methodology that integrates computational modeling, targeted experimentation, and in-service inspection considerations to demonstrate component performance and reliability under representative reactor conditions. The analysis identifies Sodium-cooled Fast Reactor (SFR) and High-Temperature Gas-cooled Reactor (HTGR) systems as reference configurations in terms of temperature, pressure, and chemical environment. Within these operating envelopes, dominant degradation mechanisms— including creep–fatigue interaction, flow-induced vibration, corrosion, and diffusion-bond deterioration—were evaluated to define test requirements. A comprehensive computationalexperimental framework is proposed to support life prediction and qualification activities. The framework couples high-fidelity structural-mechanics, thermal-hydraulic, and fluid-structure interaction models with accelerated degradation testing to produce a traceable linkage between microstructural evolution, mechanical performance, and remaining useful life (RUL). The approach adheres to established Verification, Validation, and Uncertainty Quantification (VVUQ) standards (ASME V&V 10/20; NUREG-2152) and incorporates a digital-twin architecture for continuous model refinement through data assimilation. The plan further outlines testing methodologies, including pre-test analyses, test-loop design parameters, and sensor placement strategies that maximize information yield while maintaining mechanistic fidelity. Complementary sections describe in-service inspection (ISI), on-line monitoring (OLM), and structural-health-monitoring (SHM) techniques applicable to compact HX geometries typical of advanced reactors. Collectively, these activities establish the technical foundation for demonstrating 40-60-year equivalent service life of advanced heat exchangers in support of the U.S. Department of Energy’s Advanced Reactor and Integrated Energy Systems programs. The forthcoming phase will execute the defined pre-test analyses, initiate hardware fabrication, and implement the integrated testing campaign to validate the proposed qualification methodology.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Corrosion Response of Ni-19Cr-5Fe in Molten NaCl-MgCl 2 Salt

Ni-based alloys are strong candidates for use in high-temperature molten salt reactors due to their superior corrosion resistance and mechanical stability compared to stainless steels. In this study, the static corrosion behavior of a Ni-19Cr-5Fe model alloy was systematically evaluated in a purified molten NaCl-MgCl 2 salt at 700°C for 30, 240, and 500 h. Post-exposure analyses were conducted to assess microstructural evolution, corrosion depth, and elemental depletion profiles. Corrosion rates, quantified by chromium depletion depth, followed an inverse power-law trend with increasing exposure time, indicating diffusion-limited kinetics. This trend is attributed to limitations in mass transport in the static salt and to the progressive local depletion of reactive chromium species. Coupled electron backscatter diffraction and energy dispersive X-ray spectroscopy analysis further revealed that the grain boundary character significantly influences corrosion susceptibility: high-angle grain boundaries exhibited pronounced Cr depletion and pitting, while low-angle and Σ3 boundaries remained comparatively resistant. These results offer valuable insight into the role of microstructure in corrosion processes and reinforce the importance of time-dependent material evaluation in molten salt environments relevant to advanced reactor designs.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Modeling Microwave-Enhanced Chemical Vapor Infiltration Process for Preventing Premature Pore Closure

The chemical vapor infiltration (CVI) process involves infiltrating a porous preform with reacting gases that undergo chemical transformation at high temperatures to deposit the ceramic phase within the pores, ultimately leading to a dense composite. The conventional CVI process in composite manufacturing needs to follow an isothermal approach to minimize temperature differences between the external and internal surfaces of the preform, ensuring that reactive gases infiltrate internal pores before external surfaces seal. Here, this study addresses the challenge of premature pore closure in CVI processes through microwave heating. A frequency-domain microwave solver is developed in OpenFOAM to investigate volumetric heating mechanisms within the preform. Through numerical studies, we demonstrate the capability of microwave heating of creating an inside-out temperature inversion. This inversion accelerates reactions proximal to the preform center, effectively mitigating the risk of premature external pore closure and ensuring uniform densification. The results reveal a significant enhancement in temperature inversion when high-permittivity reflectors are incorporated to generate resonant waves. This microwave heating strategy is then coupled with high-fidelity direct numerical simulation (DNS) of reacting flow, enabling the analysis of resulting densification processes. The DNS includes detailed chemistry and realistic diffusion coefficients. The numerical results can be used to estimate the impact of microwave-induced temperature inversion on densification in productions.

42 ENGINEERING↗

DOE OSTI2025

Concurrent two-way coupling of global and local models across internal boundaries with non-matching discretizations

Coupling local and global models enables efficient simulation of multiscale systems, where global models capture large-scale behavior and local models, with enhanced physics, resolve finer details over a smaller region. Here, this paper presents a mathematically consistent method for coupling physics-based models of varying fidelity across adjacent, non-overlapping subdomains, even when discretizations do not match at the immersed interdomain interfaces. Incompressible Navier-Stokes equations (NSE) constitute the global model while residual-based turbulence model serves as the local high-fidelity model. In addition, a scalar advection-diffusion equation that models the convection of an active scalar field is appended to the turbulence model in the local domain. This scalar field does not have its complement in the global model, giving rise to unequal number of equations at the immersed boundary between local and global models. Interdomain coupling terms are derived via the Variational Multiscale Discontinuous Galerkin (VMDG) method with new developments in scale representation and efficient fine-scale estimation. While transient laminar flows modeled with NSE in the global domain can be resolved with relatively coarse mesh, turbulent flow calculations in the local model require much finer spatial discretizations as well as smaller time-step for appropriately resolving the turbulent flow physics. The proposed framework also accommodates non-matching meshes at the immersed boundaries. Test problems in 2D and 3D numerically showcase the concurrent two-way coupling of unknown fields across the immersed boundaries. The 3D test presents a case with an unequal number of equations, where the scalar field represents the convection of contaminant concentration. This provides more detailed physics in the local region and highlights its application in climate modeling and atmospheric sciences.

Variational Multiscale Discontinuous Galerkin (VMD↗

DOE OSTI2025

Sensitivities of time-dependent temperature profile predictions for NSTX with the multi-mode model

The Multi-Mode Model (MMM) for turbulent transport was applied to a large set of well-analyzed discharges from the National Spherical Torus Experiment (NSTX) in order to evaluate its sensitivities to a wide range of plasma conditions. MMM calculations were performed for hundreds of milliseconds in each discharge by performing time-dependent predictive simulations with the 1.5D tokamak integrated modeling code TRANSP. A closely related study (Lestz et al 2025 Plasma Phys. Control. Fusion 67 105029) concluded that MMM predicted electron and ion temperature profiles that were in reasonable agreement with NSTX observations, generally outperforming a different reduced transport model, TGLF. This finding motivates the more thorough investigation of the characteristics of the MMM predictions conducted in this work. The simulations with MMM have electron energy transport dominated by electron temperature gradient modes in the examined discharges with relatively low plasma β (ratio of kinetic plasma pressure to magnetic field pressure) and high collisionality, transitioning to a mixture of different modes for higher β and lower collisionality. The thermal ion diffusivity predicted by MMM is much smaller than the neoclassical contribution, in line with previous experimental analysis of NSTX. Nonetheless, the electron and ion temperature profiles are coupled via collisional energy exchange and thus sensitive to which transport channels are predicted. The time-dependent simulations with MMM are robust to the simulation start time, converging to remarkably similar temperature profiles later during the discharge. MMM typically overpredicts confinement relative to NSTX observations, leading to the prediction of overly steep temperature profiles. Plasmas with spatially broader temperature profiles, higher plasma β, and longer energy confinement times tend to be predicted by MMM with better agreement with the experiment. As a result, these findings provide useful context for understanding the regime-dependent tendencies of MMM in anticipation of self-consistent, time-dependent predictive simulations of NSTX-U discharges with these same modeling tools.

MMM↗

DOE OSTI2024

Modeling Microwave-Enhanced Chemical Vapor Infiltration Process for Preventing Premature Pore Closure

The chemical vapor infiltration (CVI) process involves infiltrating a porous preform with reacting gases that undergo chemical transformation at high temperatures to deposit the ceramic phase within the pores, ultimately leading to a dense composite. The conventional CVI process in composite manufacturing needs to follow an isothermal approach to minimize temperature differences between the external and internal surfaces of the preform, ensuring that reactive gases infiltrate internal pores before external surfaces seal. This study addresses the challenge of premature pore closure in CVI processes through microwave heating. A frequency-domain microwave solver is developed in Open-FOAM to investigate volumetric heating mechanisms within the preform. Through numerical studies, we demonstrate the capability of microwave heating of creating an inside-out temperature inversion. This inversion accelerates reactions proximal to the preform center, effectively mitigating the risk of premature external pore closure and ensuring uniform densification. The results reveal a significant enhancement in temperature inversion when high-permittivity reflectors are incorporated to generate resonant waves. This microwave heating strategy is then coupled with high-fidelity direct numerical simulation (DNS) of reacting flow, enabling the analysis of resulting densification processes. The DNS simulation includes detailed chemistry and realistic diffusion coefficients. The numerical results can be used to estimate the impact of microwave-induced temperature inversion on densification in productions.

Ge, Wenjun↗

DOE OSTI2026

Processing-Dependent Structure and Poroelasticity of Nafion in Liquid Water

Ionomers act as the solid polymer electrolyte membrane in many modern electrochemical devices, yet the role of their nanostructure in modulating the poroelastic response remains poorly understood, especially in liquid water, where few techniques can measure simultaneous transport-mechanical properties. Poroelastic Relaxation Indentation (PRI) is uniquely suited for measuring time-dependent transport-mechanical properties of porous solids, specifically hydraulic diffusivity, elastic modulus, Poisson’s ratio, and intrinsic permeability, for porous solids. While ionomers such as Nafion are not porous in the typical sense, Nafion has a nanophase-segregated structure that, when fully swollen in liquid water, behaves as a poroelastic solid with a coupled mechanical-transport response. Using a poroelastic framework, we investigate how casting and pretreatment of Nafion membranes alter their poroelastic response in liquid environments. We characterize both extruded and dispersion-cast Nafion membranes pretreated in water at 25 or 100 °C and relate the mechanical-transport properties to the ionomer structure via hydrophilic and intercrystalline domain spacing measured using Small-Angle X-ray Scattering (SAXS). At 25 °C, dispersion-cast membranes exhibit distinctly lower hydraulic diffusivity and intrinsic permeability than extruded membranes but with comparable mechanical properties. Pretreatment at 100 °C increases hydrophilic domain spacing, improving transport but at the expense of mechanical stiffness. Dispersion-cast membranes respond more strongly to pretreatment than extruded membranes. Using the Carman-Kozeny pore network model and the hydrophilic domain-spacing, we estimate the pore radius and tortuosity to show how pretreatment reduces structure-related tortuosity differences between dispersion-cast and extruded membranes. Here, in this work, we show that nanophase-segregated materials such as Nafion can be rigorously characterized using poroelasticity, resulting in physically meaningful transport-mechanical properties. Coupling PRI with SAXS provides insights into the nanostructural features that govern the coupled mechanical-transport response. By establishing PRI for a nanophase-segregated material, this approach opens avenues for this technique’s application in other hydrated polymeric materials not typically considered to be poroelastic.

Shen, Margaret [University of California, Berkeley↗

DOE OSTI2025

Insights from designing an artificial cascade catalysis system using principles from substrate channeling in enzymes

Generalizing the key requirements of highly-selective, multi-step chemical conversions involving spatially separated reaction centers remains one of the grand challenges of chemistry. Much work towards this effort has focused on decomposing multi-step conversions into their constituent reactions, whose intermediates are successively upgraded in a chemical cascade via diffusion from center to center. This approach for synthesizing more complex molecules takes its cues from biochemical networks, where near-unit conversion of even complex carbohydrates is achieved by upgrading chemical precursors via enzymatic cascades. In this computational study we examine a simple cascade involving coupled Ag and Cu catalysts that sequentially converts CO2 to CO and then CO2 and CO to reduced products, generically named CO2Product and COProduct. The system architecture is inspired by the phenomenon of biological substrate channeling, and components are examined to evaluate their effects on conversion efficiency in the cascade. Aside from a substrate channel linking two reaction centers, we find efficient cascades must also incorporate directional substrate diffusion, compartmentalization of the reaction centers, and proper timing of substrate arrival at the active center. We make explicit linkages between these requirements and chemical conversion in known biological systems, revealing additional control elements that could be incorporated.

CO 2 reduction↗

DOE OSTI2024

Controlled patterning of crystalline domains by frontal polymerization

Materials with hierarchical architectures that combine soft and hard material domains with coalesced interfaces possess superior properties compared with their homogeneous counterparts. These architectures in synthetic materials have been achieved through deterministic manufacturing strategies such as 3D printing, which require an a priori design and active intervention throughout the process to achieve architectures spanning multiple length scales. Here we harness frontal polymerization spin mode dynamics to autonomously fabricate patterned crystalline domains in poly(cyclooctadiene) with multiscale organization. This rapid, dissipative processing method leads to the formation of amorphous and semi-crystalline domains emerging from the internal interfaces generated between the solid polymer and the propagating cure front. The size, spacing and arrangement of the domains are controlled by the interplay between the reaction kinetics, thermochemistry and boundary conditions. Small perturbations in the fabrication conditions reproducibly lead to remarkable changes in the patterned microstructure and the resulting strength, elastic modulus and toughness of the polymer. Furthermore, this ability to control mechanical properties and performance solely through the initial conditions and the mode of front propagation represents a marked advancement in the design and manufacturing of advanced multiscale materials. Drawing inspiration from biological systems in which structural complexity develops through dissipative reaction–diffusion processes, this study explores a transformative synthetic manufacturing strategy aimed at harnessing the principles underpinning morphogenic growth, unlocking new avenues for advanced materials design and fabrication. Synthetic coupled reaction-transport processes offer a versatile yet relatively underexplored method to manipulate the spatial attributes of synthetic materials10. Here we introduce an innovative manufacturing approach based on frontal ring-opening metathesis polymerization (FROMP) that draws parallels with morphogenic growth and development, enabling the formation of patterned microstructures within polymeric materials.

36 MATERIALS SCIENCE↗

DOE OSTI2025

From disorganized data to emergent dynamic models: Questionnaires to partial differential equations

Starting with sets of disorganized observations of spatially varying and temporally evolving systems, obtained at different (also disorganized) sets of parameters, we demonstrate the data-driven derivation of parameter dependent, evolutionary partial differential equation (PDE) models capable of generating the data. This tensor type of data is reminiscent of shuffled (multidimensional) puzzle tiles. The independent variables for the evolution equations (their “space” and “time”) as well as their effective parameters are all emergent , i.e. determined in a data-driven way from our disorganized observations of behavior in them. We use a diffusion map based questionnaire approach to build a smooth parametrization of our emergent space/time/parameter space for the data. This approach iteratively processes the data by successively observing them on the “space,” the “time” and the “parameter” axes of a tensor. Once the data become organized, we use machine learning (here, neural networks) to approximate the operators governing the evolution equations in this emergent space. Our illustrative examples are based (i) on a simple advection–diffusion model; (ii) on a previously developed vertex-plus-signaling model of Drosophila embryonic development; and (iii) on two complex dynamic network models (one neuronal and one coupled oscillator model) for which no obvious smooth embedding geometry is known a priori. This allows us to discuss features of the process like symmetry breaking, translational invariance, and autonomousness of the emergent PDE model, as well as its interpretability.

generative models↗

DOE OSTI2026

Conditional diffusion machine-learning framework for mapping valence electron distribution from convergent beam electron diffraction

Quantitative convergent beam electron diffraction (CBED) enables determination of aspherical valence electron distributions through refinement of low-order structure factors, which are highly sensitive to chemical bonding and charge density variations. However, conventional quantitative CBED (QCBED) requires solving a highly nonlinear inverse problem with many coupled parameters, and computationally intensive dynamical diffraction calculations, making it time-consuming and difficult to apply to complex systems. More broadly, reconstructing charge density and orbital electron distribution from diffraction data has long been a central challenge in both x-ray and electron crystallography. Here, in this study, we introduce an artificial-intelligence (AI)-based framework that replaces traditional refinement with a data-driven inverse solver. Using a large synthetic CBED dataset generated by Bloch-wave simulations, we train a conditional diffusion model to directly infer crystal structural parameters and multipole density formalism parameters, and hence valence electron distributions, from CBED patterns alone. By learning from forward simulations across realistic parameter space, the model effectively solves the inverse problem. Compared with direct regression approaches, the diffusion-based framework provides posterior parameter distributions for rigorous uncertainty quantification while preserving quantitative fidelity and reducing analysis time by orders of magnitude. By eliminating the need for external single-crystal x-ray diffraction data and complex nonlinear refinement, this approach enables practical, high-throughput, and in situ quantitative CBED, enabling real-time mapping of valence electron distributions and their correlation with functional responses in quantum and energy materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Modeling kinetic effects of charged vacancies on electromechanical responses of ferroelectrics: Rayleighian approach

Understanding the time-dependent effects of charged vacancies on the electromechanical responses of materials is at the forefront of research for designing materials exhibiting metal-insulator transitions and memristive behavior. A Rayleighian approach is used to develop a model for studying the nonlinear kinetics of the reaction leading to generation of vacancies and electrons via the dissociation of vacancy-electron pairs. Also, diffusion and elastic effects of charged vacancies are considered to model polarization-electric potential and strain-electric potential hysteresis loops. The model captures multiphysics phenomena by introducing couplings among polarization, the electric potential, stress, strain, and concentrations of charged (multivalent) vacancies and electrons (treated as classical negatively charged particles), where the concentrations can vary due to association-dissociation reactions. A derivation of coupled time-dependent equations based on the Rayleighian approach is presented. Three limiting cases of the governing equations are considered, highlighting the effects of (1) nonlinear reaction kinetics on the generation of charged vacancies and electrons, (2) Vegard's law (i.e., the concentration-dependent local strain) on asymmetric strain-electric potential relations, and (3) coupling between a fast component and the slow component of the net polarization on the polarization-electric-field relations. The Rayleighian approach discussed in this work should pave the way for developing a multiscale modeling framework in a thermodynamically consistent manner while capturing multiphysics phenomena in ferroelectric materials. Published by the American Physical Society 2025

Kumar, Rajeev (ORCID:0000000194943488)↗

DOE OSTI2026

A multiscale packed-bed reactor model for sustainable ethylene production via chemical looping oxidative coupling of methane

The rising global warming concerns and shale gas discovery have prompted research in the direction of greenhouse gas (GHG), such as methane, reduction and conversion. Oxidative coupling of methane (OCM) offers a pathway to low carbon-intense valorization of methane while producing ethylene, a chemical regarded as central to the petrochemical industry. Even after decades of OCM discovery, researchers keep understanding the process and underlying chemical reactions in a pursuit to achieve industrial viability for OCM. Here, in general, OCM suffers from low C 2 selectivity, yield and reactor temperature runaways due to highly exothermic nature of its reactions. Computational Fluid Dynamics (CFD) tools help analyze spatial gradients within the reactor to deeply understand the diffusion of species, mass and heat transfer phenomena. Furthermore, challenges associated with scaling up such as hot spot formation and parametric sensitivity can be addressed without having to expend on costly experiments. The current paper presents a multiscale packed-bed reactor CFD model coupled with a chemical kinetic model for the chemical looping OCM. The CFD model includes two scales i.e., macroscale for catalyst bed and microscale for individual pellets. Moreover, a chemical kinetic model based on 10 gas-phase reactions is integrated with the CFD model. An additional surface reaction for the formation of gas-phase oxygen from catalyst surface is added to account for the absence of feed oxygen. The model is calibrated against experimental results. The calibrated model captures trends in CH 4 conversion, C 2 selectivity and C 2 yield within a ± 4.35 % range across a temperature range of 700-900 °C. Moreover, model fidelity is evaluated by varying key computational parameters such as mesh resolution and time step size. The model is also verified by varying the inlet methane concentration and the gas hourly space velocity (GHSV) and comparing the results with literature. A sensitivity analysis and scale-up of the current model is undergoing.

Chemical looping↗

DOE OSTI2024

Fast Oxygen Redox Kinetics Induced by CoO 6 Octahedron With π –Interaction in P2–Type Sodium Oxides

Enhancing the kinetics of lattice oxygen redox (LOR) in P2-type layered sodium oxide cathodes is crucial for the advancement of sodium-ion batteries (SIBs) with superior energy and power densities. Electronic structure regulation stands out as a highly effective approach to address the inherent limitations of P2-type layered oxides with LOR, including sluggish kinetics, phase transitions, voltage hysteresis, and local structural distortion. In this work, a strategy involving the introduction of CoO 6 octahedra with π-interaction into Na 0.6 Li 0.1 Fe 0.3–x Co x Mn 0.6 O 2 (x = 0, 0.15, 0.3) cathodes to facilitate Na-ion transport is proposed. Furthermore, the impact of FeO 6 octahedra with σ-interaction in P2-type cathodes on electrochemical performance is comprehensively investigated. Through multimodal in-situ and ex-situ characterization techniques, it is revealed that Co–O with π-interaction effectively mitigates P2-OP4 phase transitions by strengthening Na–O, reduces voltage hysteresis, and stabilizes the local structure. Consequently, Na 0.6 Li 0.1 Co 0.3 Mn 0.6 O 2 demonstrates enhanced Na-ion diffusion kinetics, leading to improved rate performance and a reversible capacity of 55 mAh g –1 at 10 C, significantly outperforming cathodes with Fe–O σ-interaction. Moreover, when coupled with hard carbon, the full cell achieves a remarkable energy density of 395 Wh kg –1 (on cathode) at 0.1 C, with a capacity retention of 75% over 100 cycles at 1 C.

25 ENERGY STORAGE↗

DOE OSTI2025

Investigation of Corrosion Behavior of Carbide Phase Strengthened Nickel‐Based Alloy in Molten Fluoride Salt

The corrosion behavior of a high creep strength carbide-phase strengthened Ni-based alloy in molten FLiNaK (LiF-NaF-KF: 46.5-11.5-42 mol%) salt in the temperature range of 700°C–750°C has been investigated as a part of the development of structural alloys for fluoride salt-based molten salt reactors (MSRs). The alloy composition was designed based on Hastelloy-N, but with the goal of improving creep strength. Cr depletion depth, a measure of corrosion, was observed to be single micrometers after several hundred hours of corrosion testing. Sequential corrosion testing involving testing of pre-tested samples in fresh salt coupled with SEM-EDS, scanning transmission electron microscopy (STEM) examinations and thermodynamic and kinetic modeling, suggest that the corrosion rate at the alloy-salt interface is governed by diffusion of elements from the alloy bulk to the surface. The carbide phases in the corrosion-tested sample microstructure were identified to be largely M 6 C-type Mo-rich carbides and MC-type mixed carbides. Atom probe tomography (APT) showed some partitioning of Cr and Ti to the carbide phase and showed the carbide phases to be stable at the salt-facing surface.

Gordon, Ryan [Univ. of Wisconsin, Madison, WI (Uni↗

DOE OSTI2026

Contrasting structural reversibility and magnetic correlations in isostructural honeycomb magnets CrCl3 and 𝛼−RuCl3

We report a comparative neutron single crystal diffraction study of the structural and magnetic properties of layered halides CrCl3 and 𝛼−RuCl3. They host a honeycomb arrangement of transition metal ions with distinct electronic configurations and undergo a first-order structural transition between high-temperature 𝐶⁢2/𝑚 and low-temperature 𝑅⁢‾‾‾3. Both compounds show a step-like change in the 𝑐-lattice parameter across the structure transition. In contrast, the in-plane lattice response is quite different: 𝛼−RuCl3 exhibits an abrupt hysteretic change across the transition accompanied by progressive crystalline degradation upon thermal cycling, whereas CrCl3 shows a smooth in-plane lattice evolution and remains structurally robust. Magnetically, CrCl3 orders into an A-type antiferromagnetic structure at 𝑇𝑁=14 K and exhibits pronounced diffuse magnetic scattering extending up to about 40 K. 𝛼−RuCl3 shows no observable magnetic diffuse scattering above its zigzag antiferromagnetic ordering temperature 𝑇𝑁=7.6 K. These results suggest that the contrasting structural responses arise from an interplay between interlayer sliding energetics, stacking-strain coupling, and elastic accommodation of the stacking transition. The distinct chemical bonding and electronic configurations of the two compounds provide a microscopic basis for their different lattice responses to the structure transition and magnetic correlations.

Morgan, Zachary [ORNL] (ORCID:0000000243625911)↗

DOE OSTI2024

First-principles study of the Stark shift effect on the zero-phonon line of the NV center in diamond

Point defects in semiconductors are attractive candidates for quantum information science applications owing to their ability to act as spin-photon interface or single-photon emitters. However, the coupling between the change of dipole moment upon electronic excitation and stray electric fields in the vicinity of the defect, an effect known as Stark shift, can cause significant spectral diffusion in the emitted photons. In this work, using first principles computations, we revisit the methodology to compute the Stark shift of point defects up to the second order. The approach consists of applying an electric field on a defect in a slab and monitoring the changes in the computed zero-phonon line (i.e., difference in energy between the ground and excited state) obtained from constraining the orbital occupations (constrained-DFT). Here, we study the Stark shift of the negatively charged nitrogen-vacancy (NV) center in diamond using this slab approach. We discuss and compare two approaches to ensure a negatively charged defect in a slab and we show that converged values of the Stark shift measured by the change in dipole moment between the ground and excited states (Δ⁢μ) can be obtained. We obtain a Stark shift of Δ⁢μ = 2.68⁢D using the semilocal GGA-PBE functional and of Δ⁢μ = 2.23⁢D using the HSE hybrid functional. These values are in good agreement with experimental results. We also show that modern theory of polarization can be used on constrained-DFT to obtain Stark shifts in very good agreement with the slab computations.

36 MATERIALS SCIENCE↗