Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “first principles simulations”

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.

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

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Simulations of Attosecond Metallization in Quartz and Diamond Probed with Inner-Shell Transient Absorption Spectroscopy

When dielectrics are hit with intense infrared (IR) laser pulses, transient metalization can occur. The initial attosecond dynamics behind this metallization are not entirely understood. Therefore, simulations are needed to understand this process and to help interpret experimental observations of it, such as with attosecond transient absorption (ATA). In this paper, we present first-principles simulations of ATA based on bulk-mimicking clusters and real-time time-dependent density functional theory (RT-TDDFT), with Koopmans-tuned range-separated hybrid functionals and Gaussian basis sets. Our method gives good agreement with the experiment for the breakdown threshold in silica and diamond. This breakdown voltage corresponds to a Keldysh parameter of approximately one and thus involves a transition to a regime where the dynamics are driven by tunneling. Pumping at an amplitude just below this value causes a mixture of multiphoton and tunneling excitations across the band gap to occur. The computed extreme ultraviolet and X-ray attosecond transient spectra also agree well with the experiment and show a decrease in optical density due to the transient population of the conduction band from the IR field. First-principles approaches such as this are valuable for interpreting the complicated modulations in a spectrum and for guiding future attosecond experiments on solids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Resolving the dynamic correlated disorder in KTa 1- x Nb x O 3

Understanding the complex temporal and spatial correlations of ions in disordered perovskite oxides is critical to rationalize their functional properties. Here, in this study, we provide new insights into the longstanding controversy regarding the off-centering of transition metal ions in the archetypal ferroelectric alloy KTa 1-x Nb x O 3 (KTN). By mapping the full energy (E) and wavevector (Q) dependence of the dynamical structure factor S(Q, E) using neutron scattering, and rationalizing our observations with atomistic simulations leveraging machine learning, we fully resolve the static vs dynamic nature of diffuse scattering sheets, as well as their composition (x) and temperature dependence. Our first-principles simulations, extended with machine-learning molecular dynamics, reproduce both inelastic neutron spectra and diffuse features, and establish how dynamically-correlated transition metal off-centerings couple to phonons, unifying local and collective viewpoints. This study sheds new light into an exemplary ferroelectric systems and shows the importance of mapping the full S(Q, E) to reveal critical spatio-temporal correlations of atomic disorder from which functional properties emerge.

42 ENGINEERING↗
DOE OSTI2025

SciDAC ISEP: Integrated Simulation of Energetic Particles in Burning Plasmas

The objective of the SciDAC Center for Integrated Simulation of Energetic Particles in Burning Plasmas (ISEP) is to improve physics understanding of energetic particle (EP) confinement and EP interactions with burning thermal plasmas through large-scale simulations. The ISEP center will develop a multiscale and multiphysics ISEP framework for a predictive capability of EP physics and deliver an EP module incorporating both first-principles simulations and high fidelity reduced transport models to the fusion whole device modeling (WDM) project. The ISEP framework will enable us to perform long time, global kinetic simulations of EP physics in burning plasmas, by utilizing the full power of the next generation supercomputers. Our research and development activities will build on fruitful collaborations with computer scientists and applied mathematicians to offer enabling technologies for performance scalability, portability, solvers, coupling for integration with the fusion WDM project, and long-term preservation of data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Machine learning inversion of interatomic force constants from single-crystal inelastic neutron scattering

Atomic vibrations govern many macroscopic properties of materials, but experiments to comprehensively probe them remain challenging. Inelastic neutron scattering (INS) is a powerful technique to map phonon dispersions in crystals, especially when leveraging modern time-of-flight (ToF) spectrometers with large detectors. However, efficiently and robustly extracting interatomic force constants (FCs) parameterizing phonon dynamics from experimental spectra remains a bottleneck due to the complexity and high dimensionality of ToF INS datasets. Here, we present a machine learning approach for the direct inversion of FCs from single-crystal INS measurements. The framework leverages synthetic training data generated using universal machine-learned force fields and an efficient physics-based forward model. We benchmark two neural architectures–one emphasizing structured latent representation learning and the other direct, supervised spectral regression–across simulated datasets for two materials under idealized and noisy conditions. The latent-representation model is subsequently applied to experimental single-crystal INS data on germanium. The model is shown to reproduce FCs derived from both first-principles simulations and from iterative optimization, and furthermore achieves reliable inference even from sparse, single-orientation measurements representing short data acquisitions. Analysis of the learned latent space reveals semantically continuous and physically interpretable encodings that support strong cross-domain generalization. By bridging theoretical and experimental domains, we establish a path toward rapid inversion of experimental spectra and data-driven interpretation of temperature-dependent lattice dynamics.

42 ENGINEERING↗
DOE OSTI2025

OmniXAS: A universal deep-learning framework for materials x-ray absorption spectra

X-ray absorption spectroscopy (XAS) is a powerful characterization technique for probing the local chemical environment of absorbing atoms. However, analyzing XAS data presents significant challenges, often requiring extensive, computationally intensive simulations, as well as significant domain expertise. These limitations hinder the development of fast, robust XAS analysis pipelines that are essential in high-throughput studies and for autonomous experimentation. Here, we address these challenges with OmniXAS, a framework that contains a suite of transfer learning approaches for XAS prediction, each uniquely contributing to improved accuracy and efficiency, as demonstrated on the K-edge spectra database covering eight 3⁢d transition metals (Ti–Cu). The OmniXAS framework is built upon three distinct strategies. First, we use M3GNet [Nat. Comput. Sci. 2, 718 (2022)] to derive latent representations of the local chemical environment of absorption sites as input for XAS prediction, achieving significant improvements over conventional featurization techniques. Second, we employ a hierarchical transfer learning strategy, training a universal multitask model across elements before fine-tuning for element-specific predictions. Models based on this cascaded approach after elementwise fine-tuning outperform element-specific models by up to 69%. Third, we implement cross-fidelity transfer learning, adapting a universal model to predict spectra generated by simulation of a different fidelity with a much higher computational cost. This approach improves prediction accuracy by up to 11% over models trained on the target fidelity alone. Our approach significantly boosts the throughput of XAS modeling by orders of magnitude as compared to first-principles simulations and is extendable to XAS prediction for a broader range of elements. The proposed transfer learning framework is generalizable to enhance deep-learning models that target other properties in materials research.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Computationally guided experimental validation of divacancy defect formation in 4H-SiC

Recent research into solid-state qubits for quantum information science has focused on optically addressable spin defects such as the negatively charged nitrogen-vacancy center in diamond and the neutrally charged divacancy (VV) in 4H-SiC as scalable quantum sensors and networking qubits. Within this context, direct investigations of the structural origin and defect formation dynamics of a sub-set of the VV center in 4H-SiC remain lacking. Here, we take a systematic experimental approach guided by predictions from first-principles simulations to gain a thorough mechanistic understanding of the VV defect formation and control in 4H-SiC. We study the effect of annealing time and temperature on VV formation in high-purity semi-insulating 4H-SiC samples following electron irradiation. Three different temperatures (1123, 1273, and 1473 K) and annealing duration (from 0.5 to 72 h) are chosen to explore VV formation in different regions. We find that samples annealed at 1273 K give the highest VV-related photoluminescence (PL) intensities, in agreement with the prediction from first-principles calculations. Furthermore, the logarithmic dependence of VV-related PL intensities on the annealing duration at 1273 K indicates that 1273 K provides sufficient thermal energy for silicon vacancy migration but not for VV migration. Together, these results suggest that efficient VV formation occurs above the V Si migration temperature and below the VV migration threshold.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Molecular dynamics simulation of hydrodynamic transport coefficients in plasmas

Molecular dynamics simulations are used to calculate transport coefficients in a two-component plasma interacting through a repulsive Coulomb potential. The thermal conductivity, electrical conductivity, electrothermal coefficient, thermoelectric coefficient, and shear viscosity are computed using the Green–Kubo formalism over a broad range of Coulomb coupling strength, 0.01 ≤ Γ ≤ 140. Emphasis is placed on testing standard results of the Chapman–Enskog solution in the weakly coupled regime (Γ ≪ 1) using these first-principles simulations. As expected, the results show good agreement for Γ ≲ 0.1. However, this agreement is only possible if careful attention is paid to the definitions of linear constitutive relations in each of the theoretical models, a point that is often overlooked. For example, the standard Green–Kubo expression for thermal conductivity is a linear combination of thermal conductivity, electrothermal, and thermoelectric coefficients computed in the Chapman–Enskog formalism. Meaningful results for electrical conductivity are obtained over the full range of coupling strengths explored, but it is shown that potential and virial components of the other transport coefficients diverge in the strongly coupled regime (Γ ≫ 1). In this regime, only the kinetic components of the transport coefficients are meaningful for a classical plasma.

Electrical conductivity↗
DOE OSTI2026

When does global attention help: a unified empirical study on atomistic graph learning

Graph neural networks (GNNs) are widely used as surrogates for costly experiments and first-principles simulations to study the behavior of compounds at atomistic scale, and their architectural complexity is constantly increasing to enable the modeling of complex physics. While most recent GNNs combine more traditional message passing neural networks (MPNNs) layers to model short-range interactions with more advanced graph transformers (GTs) with global attention mechanisms to model long-range interactions, it is still unclear when global attention mechanisms provide real benefits over well-tuned MPNN layers due to inconsistent implementations, features, or hyperparameter tuning. We introduce the first unified, reproducible benchmarking framework–built on HydraGNN–that enables seamless switching among four controlled model classes: MPNN, MPNN with chemistry/topology encoders, GPS-style hybrids of MPNN with global attention, and fully fused localglobal models with encoders. Using seven diverse open-source datasets for benchmarking across regression and classification tasks, we systematically isolate the contributions of message passing, global attention, and encoder-based feature augmentation. Our study shows that encoder-augmented MPNNs form a robust baseline, while fused localglobal models yield the clearest benefits for properties governed by long-range interaction effects. We further quantify the accuracycompute trade-offs of attention, reporting its overhead in memory. Together, these results establish the first controlled evaluation of global attention in atomistic graph learning and provide a reproducible testbed for future model development.

Equivariant graph neural networks↗
DOE OSTI2025

Bridging Scales in Black Hole Accretion and Feedback: Relativistic Jet Linking the Horizon to the Host Galaxy

Simulating black hole (BH) accretion and feedback from the BH horizon to galactic scales is extremely challenging, as it involves a vast range of scales. Recently, our multizone method has successfully achieved global dynamical steady states of hot accretion flows in 3D general relativistic magnetohydrodynamic simulations by tracking the bidirectional interaction between a nonspinning BH and its host galaxy. In this paper, we present technical improvements to the method and apply it to spin a * = 0.9 BHs, which power relativistic jets. We first test the new multizone setup with a smaller Bondi radius, R B ≈ 400 r g , where r g is the gravitational radius. The strongly magnetized accretion launches a relativistic jet with an intermediate feedback efficiency η ∼ 30%, in between that of a prograde (η ∼ 100%) and retrograde (η ∼ 10%) torus. Interestingly, both prograde and retrograde simulations also eventually converge to the same intermediate efficiency when evolved long enough, as accumulated magnetic fields remove gas rotation. We then extend strongly magnetized simulations to larger Bondi radii, R B ≈ 2 × 10 3 , 2 × 10 4 , 2 × 10 5 r g . We find that the BH accretion rate $\dot{M}$ is suppressed with respect to the Bondi rate as $\dot{M}_{\textrm{B}}$ as $\dot{M}/\dot{M}_{\textrm{B}} ∝ R_{\textrm{B}}^{-1/2}$. However, despite some variability, the time-averaged feedback efficiency remains at η ∼ 30%, independent of R B . This suggests that BH feedback efficiency in hot accretion flows is mainly governed by the BH spin (a * ) rather than by the galactic properties (R B ). From these first-principles simulations, we provide a feedback subgrid prescription for cosmological simulations: $\dot{E}_{\textrm{fb}} = 2$ x $10^{-3}[R_{\textrm{B}}/(2$ x $10^5 r_g)]^{-1/2}$ $\dot{M}_{\textrm{B}}c^2$ for BH spin a * = 0.9.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Exploring dielectric properties in atomistic models of amorphous boron nitride

Abstract We report a theoretical study of dielectric properties of models of amorphous Boron Nitride, using interatomic potentials generated by machine learning. We first perform first-principles simulations on small (about 100 atoms in the periodic cell) sample sizes to explore the emergence of mid-gap states and its correlation with structural features. Next, by using a simplified tight-binding electronic model, we analyse the dielectric functions for complex three dimensional models (containing about 10.000 atoms) embedding varying concentrations of sp 1 , sp 2 and sp 3 bonds between B and N atoms. Within the limits of these methodologies, the resulting value of the zero-frequency dielectric constant is shown to be influenced by the population density of such mid-gap states and their localization characteristics. We observe nontrivial correlations between the structure-induced electronic fluctuations and the resulting dielectric constant values. Our findings are however just a first step in the quest of accessing fully accurate dielectric properties of as-grown amorphous BN of relevance for interconnect technologies and beyond.

Materials Science↗
DOE OSTI2024

Identification of Defects and the Origins of Surface Noise on Hydrogen–Terminated (100) Diamond

Near-surface nitrogen vacancy centres are critical to many diamond-based quantum technologies such as information processors and nanosensors. Surface defects play an important role in the design and performance of these devices. The targeted creation of defects is central to proposed bottom-up approaches to nanofabrication of quantum diamond processors, and uncontrolled surface defects may generate noise and charge trapping which degrade shallow NV device performance. Surface preparation protocols may be able to control the production of desired defects and eliminate unwanted defects, but only if their atomic structure can first be conclusively identified. This work uses a combination of scanning tunnelling microscopy (STM) imaging and first-principles simulations to identify several surface defects on H:C(100)—2 × 1 surfaces prepared using chemical vapour deposition (CVD). The atomic structure of these defects is elucidated, from which the microscopic origins of magnetic noise and charge trapping are determined based on the modeling of their paramagnetic properties and acceptor states. Rudimentary control of these deleterious properties is demonstrated through STM tip-induced manipulation of the defect structure. Furthermore, the results validate accepted models for CVD diamond growth by identifying key adsorbates responsible for the nucleation of new layers.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Rigidity‐Driven Structural Isomers in the NaCl–Ga 2 S 3 System: Implications for Energy Storage

Alternative energy sources require the search for innovative materials with promising functionalities. Systems with unusual chemical properties represent an insufficiently explored domain, concealing unexpected features. Using diffraction and Raman spectroscopy over a wide temperature range, supported by first‐principles simulations, a rare phenomenon is unveiled: phase‐dependent chemical interactions between binary components in the NaCl–Ga 2 S 3 system. In this unique occurrence, previously intact binary crystalline species transform upon melting into mixed liquid structural isomers, forming bonds with new partners. The chemical combinatorics appears to be fully reversible for stable crystals and liquids. Despite this, rapidly frozen glasses out of thermodynamic equilibrium remain in a metastable isomeric state, offering remarkable properties, particularly a high room‐temperature Na + conductivity, comparable to the best sodium halide superionic conductors and therefore encouraging for sodium solid‐state batteries and energy applications. A rigidity paradigm is responsible for the observed phenomenon, as the extremely constrained Ga 2 S 3 crystal lattice does not survive viscous flow, breaking up at a short‐range level. The removal of rigidity constraints and dense packing leads to a significant increase in empty space, which is the origin of high sodium diffusivity. Broadly, the rigidity‐driven structural isomerism opens up an inspiring path to the discovery of atypical materials.

25 ENERGY STORAGE↗
DOE OSTI2024

Large Electrically and Chemically Tunable Rashba–Dresselhaus Effects in Ferroelectric CsGeX 3 (X = Cl, Br, I) Perovskites

Rashba–Dresselhaus effects, which originate from spin–orbit coupling and allow for spin manipulations, are actively explored in materials, following the pursuit of spintronics and quantum computing. However, materials that possess practically significant Rashba–Dresselhaus effects often contain toxic elements and offer little opportunity for the tunability of the effects. We used first-principles simulations to reveal that the recently discovered halide ferroelectrics in the CsGeX 3 (X = Cl, Br, I) family possess large and tunable Rashba-Dresselhaus effects. In particular, they give origin to the spin splitting of up to 171 meV in the valence band of CsGeI 3 . The value is chemically tunable and can decrease by 25% and 70% for CsGeBr 3 and CsGeCl 3 , respectively. Such chemical tunability could result in the engineering of desired values through a solid solution technique. Application of an electric field was found to result in structural changes that could decrease and increase spin splitting, leading to electrical tunability of the effect. In the vicinity of conduction and valence band extrema, the spin textures are mostly of the Rashba type, which is promising for spin-to-charge conversion applications. The spin directions are coupled with the polarization direction, leading to Rashba-ferroelectricity cofunctionality. Furthermore, our work identifies lead-free perovskite halides as excellent candidates for spin-based applications and is likely to stimulate further research in this direction.

Electric fields↗
DOE OSTI2025

Morphology, Deformations, and Photocatalytic Activity of Thermally Treated Brookite Titanium Dioxide Thin Films

Metastable states, in which the coupling between long-range lattice deformations and electronic properties can be controlled, provide a pathway to tailoring the behavior of photocatalytic materials by directing the flow of photoinduced charge carriers. Brookite is a metastable polymorph of earth-abundant TiO 2 that exhibits photocatalytic function and, due to its high energy relative to the anatase and rutile polymorphs, may serve as a precursor for the formation of transitional metastable structures. In this work, facile thermal annealing is employed to promote the formation of predominantly brookite-phase films, regulate the brookite lattice distortions, and determine the effect of these distortions on charge separation, ultimately directed at enhancing photocatalytic activity. Profile fitting of X-ray diffraction patterns and peak shifts in Raman spectra revealed structural distortions of the brookite lattice. Structural defects, including lattice gliding, dislocations, stacking faults, and twin boundaries, were observed using scanning transmission electron microscopy. First-principles simulations reveal how the lattice distortions associated with stacking faults induce band bending, thus increasing the photocatalytic activity of brookite. In conclusion, this study provides insight into the microstructural tuning of metastable phases to enhance their unique functionalities.

band bending↗
DOE OSTI2024

Nonsolvating Fluoroaromatic Cosolvent Enabled Long-Term Cycling of High-Voltage Lithium-Ion Batteries with Organosulfur Electrolytes

Here, the structure–activity relationships of nonsolvating cosolvents for organosulfur-based electrolyte systems were revealed. The performance of nonsolvating dilutant fluorobenzene (FB) was compared to various fluorinated ether dilutants in high-voltage electrolytes containing a concentration of 1.2 M LiPF 6 dissolved in fluoroethylene carbonate (FEC), ethyl methyl sulfone (EMS), and the dilutant. In a high-voltage and high-loading LiNi 0.8 Mn 0.1 Co 0.1 O 2 (NMC811) full cell configuration, the organosulfur-based electrolyte containing FB dilutant enabled superior electrochemical performance compared to the electrolytes using other nonsolvating fluorinated ether formulations. Moreover, the FB-containing electrolyte exhibited the highest ionic conductivity and lowest viscosity among all organosulfur-based electrolytes containing nonsolvating dilutant. These improvements are attributed to the enhanced physical properties of electrolyte and lithium-ion mobility. Furthermore, by employing first-principles simulations, the observed suppression of side reactions at high voltage is linked to FB’s lower reactivity toward singlet dioxygen, which is likely produced at the NMC interface. Overall, FB is considered an excellent diluent that does not impede cell operation by mass decomposition at the cathode.

25 ENERGY STORAGE↗
DOE OSTI2024

Quasi-one-dimensional hydrogen bonding in nanoconfined ice

Abstract The Bernal-Fowler ice rules stipulate that each water molecule in an ice crystal should form four hydrogen bonds. However, in extreme or constrained conditions, the arrangement of water molecules deviates from conventional ice rules, resulting in properties significantly different from bulk water. In this study, we employ machine learning-driven first-principles simulations to identify a new stabilization mechanism in nanoconfined ice phases. Instead of forming four hydrogen bonds, nanoconfined crystalline ice can form a quasi-one-dimensional hydrogen-bonded structure that exhibits only two hydrogen bonds per water molecule. These structures consist of strongly hydrogen-bonded linear chains of water molecules that zig-zag along one dimension, stabilized by van der Waals interactions that stack these chains along the other dimension. The unusual interplay of hydrogen bonding and van der Waals interactions in nanoconfined ice results in atypical proton behavior such as potential ferroelectric behavior, low dielectric response, and long-range proton dynamics.

Science & Technology - Other Topics↗
DOE OSTI2025

Nonlinear excitation of energetic particle driven geodesic acoustic mode by resonance overlap with Alfvén instability in ASDEX Upgrade

The Alfvén instability nonlinearly excited the energetic-particle-driven geodesic acoustic mode on the ASDEX-Upgrade tokamak, as demonstrated experimentally. The mechanism of the energetic-particle-driven geodesic acoustic mode excitation and the mode nonlinear evolution is not yet fully understood. In the present work, a first-principles simulation using the MEGA code investigated the mode properties in both the linear growth and nonlinear saturated phases. Here we show that the simulation successfully reproduced the excitation and coexistence of these two modes, and agreed with the experimental results well. Conclusive evidence showed that the resonance overlap is the excitation mechanism of the energetic-particle-driven geodesic acoustic mode. In the linear growth phase, energetic particles that satisfied different resonance conditions excited the Alfvén instability, which then caused energetic particle redistribution in phase space. These redistributed energetic particles caused resonance overlap, exciting the energetic-particle-driven geodesic acoustic mode in the nonlinear phase.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

First-principles study of the magneto-Raman effect in van der Waals layered magnets

Magneto-Raman spectroscopy has been used to study spin-phonon coupling in two-dimensional (2D) magnets. Raman spectra of CrI 3 show a strong dependence on the magnetic order within a layer and between the layers. Here we carry out the first systematic theoretical investigation of the magneto-Raman effect in 2D magnets by performing density functional theory calculations and developing a generalized polarizability model. Our first-principles simulations well reproduce experimental Raman spectra of CrI 3 with different magnetic states. The model reveals how the change of spin orientation in each layer is coupled to the layer’s vibration to induce or eliminate the spin-dependent anti-symmetric off-diagonal terms in the Raman tensor for altering the selection rules. We also uncover that the correlation between phonon modes and magnetic orders is a universal phenomenon, which should exist in other phonon modes and 2D magnets. Our predictive simulations and modeling are expected to guide the research in 2D magnets.

36 MATERIALS SCIENCE↗