Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ab initio phasing”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

First-Principles Study of Tritium Trapping by Point Defects in Fe-Al Aluminide Coating Phases

Density functional theory simulations have been carried out to investigate the potential for tritium trapping by metal vacancies in five different Fe-Al aluminide coating phases. It was found that tritiation of Fe and Ni vacancies is generally less favorable than the tritiation of Al vacancies. However, for the first tritiation, a trend in the defect formation energy can be obtained such that metal defects in the Fe 2 Al x family of materials (i.e., Fe 2 Al 4 , Fe 2 Al 5 , and Fe 2 Al 6 ) trap tritium species more favorably than metal vacancies in Fe 4 Al 13 and FeNiAl 5 . Further investigations using ab initio thermodynamics calculations confirmed that trend for a range of tritium partial pressure at a temperature of 700 K. Especially, it was shown that the energy difference between tritiated and non-tritiated metal vacancies is smaller and more favorable for the Fe 2 Al x family, followed by Fe 4 Al 13 , and FeNiAl 5 . While this study shows that tritium interacts differently in the various Fe-Al aluminide phases, it also suggests that tritium trapping and retention could be more efficient if metal defects are present in some Fe-Al phases.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Development of an interatomic potential for the Ta–Li system

A new interatomic potential for the Ta–Li system is introduced to facilitate the study of phase stability, mechanical properties and non-equilibrium dynamics after Li implantation in Ta. Here, this potential is based on a generalization of the embedded atom method (GEAM) and includes contributions from embedding energy, explicit two- and three-body interactions, and nonlocal many-body interaction terms. The parameters of the potential are optimized using energies and atomic forces for a wide range of configurations obtained from ab initio density functional theory (DFT) calculations. The potential is rigorously validated across a range of physical properties, including elastic constants, equations of state, phonon dispersion curves, point defect properties, and melting temperatures for different compositions. Although our potential is trained on a small dataset, its accuracy is comparable to that of available machine learning potentials for Li and Ta. Our simulations show that at temperatures below 500 K, Li atoms in Ta–Li alloys form clusters separated by Ta-rich domains, and we find no evidence of ordered phase formation. For Li concentrations below a few percent, Li atoms preferentially segregate to surfaces and grain boundaries. However, in alloys containing more than ~10% Li, the accumulation of Li in symmetric-tilt grain boundaries can lead to one of the following effects: formation of amorphous-like regions, changes in grain boundary structural units, or lateral movement of the grain boundary.

GEAM potential↗
DOE OSTI2021

Discrete quantum geometry and intrinsic spin Hall effect

We show that the quantum geometry of the Fermi surface can be numerically described by a three-dimensional discrete quantum manifold. This approach not only avoids singularities in the Fermi sea, but it also enables the precise computation of the intrinsic Hall conductivity resolved in spin, as well as any other local properties of the Fermi surface. The method assures numerical accuracy when the Fermi level is arbitrarily close to singularities, and it remains robust when Kramers degeneracy is protected by symmetry. Furthermore, the approach is demonstrated by calculating the anomalous Hall and spin Hall conductivities of a two-band lattice model of a Weyl semimetal and a full-band ab initio model of zinc-blende GaAs.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

A two-dimensional type I superionic conductor

Superionic conductors (SICs) possess liquid-like ionic diffusivity in the solid state, finding wide applicability from electrolytes in energy storage to materials for thermoelectric energy conversion. Type I SICs (e.g., AgI, Ag 2 Se, etc.) are defined by a first-order transition to the superionic state and have so far been found exclusively in three-dimensional crystal structures. Here, we reveal a two-dimensional type I SIC, α-KAg 3 Se 2 by scattering techniques and complementary simulations. Quasi-elastic neutron scattering and ab initio molecular dynamics simulations confirm that the superionic Ag + ions are confined to sub-nanometre sheets, with the simulated local structure validated by experimental X-ray powder pair-distribution-function analysis. Finally, we demonstrate that the phase transition temperature can be controlled by chemical substitution of the alkali metal ions that comprise the immobile charge-balancing layers. Our work thus extends the known classes of SICs and will facilitate the design of new materials with tailored ionic conductivities and phase transitions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Signatures of a liquid–liquid transition in an ab initio deep neural network model for water

Significance Water is central across much of the physical and biological sciences and exhibits physical properties that are qualitatively distinct from those of most other liquids. Understanding the microscopic basis of water’s peculiar properties remains an active area of research. One intriguing hypothesis is that liquid water can separate into metastable high- and low-density liquid phases at low temperatures and high pressures, and the existence of this liquid–liquid transition could explain many of water’s anomalous properties. We used state-of-the-art approaches in computational quantum chemistry, statistical mechanics, and machine learning and obtained evidence consistent with a liquid–liquid transition, supporting the argument for the existence of this phenomenon in real water.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Natural Charge-Transfer Analysis: Eliminating Spurious Charge-Transfer States in Time-Dependent Density Functional Theory via Diabatization, with Application to Projection-Based Embedding

For many types of vertical excitation energies, linear-response time-dependent density functional theory (LR-TDDFT) offers a useful degree of accuracy combined with unrivaled computational efficiency, although charge-transfer excitation energies are often systematically and dramatically underestimated, especially for large systems and those that contain explicit solvent. As a result, low energy electronic spectra of solution-phase chromophores often contain tens to hundreds of spurious charge-transfer states, making LR-TDDFT needlessly expensive in bulk solution. More nefariously, intensity borrowing by the low-energy charge-transfer states can affect intensities of the valence excitations even if those excitation energies are accurate. At higher excitation energies, it is difficult to distinguish spurious CT states from genuine charge-transfer-to-solvent (CTTS) excitations. In this work, we introduce an automated diabatization scheme that enables fast and effective screening of the CTTS acceptor space in bulk solution. Our procedure introduces the concept of “natural charge-transfer orbitals”, which provide a means to isolate characteristic pairs of orbitals that are most likely to participate in a CTTS excitation. The projection of these orbitals onto solvent-centered virtual orbitals provides a criterion for defining the most important solvent molecules in a given excitation. We apply this method to analyze an ab initio molecular dynamics (MD) trajectory of I-(aq) and report the lowest-energy CTTS band in the absorption spectrum. Our results are in excellent agreement with experimental measurements for bulk I-(aq), and only one-third of the water molecules in the I-(H2O)96 simulation cell need to be described with LR-TDDFT in order to obtain excitation energies that are converged to < 0.1 eV. The tools introduced herein will improve the accuracy, efficiency, and usability of LR-TDDFT in solution-phase environments

Carter-Fenk, Kevin D.↗
DOE OSTI2024

Emulating ab initio computations of infinite nucleonic matter

We construct efficient emulators for the computation of the infinite nuclear matter equation of state. These emulators are based on the subspace-projected coupled-cluster method for which we here develop a new algorithm called small-batch voting to eliminate spurious states that might appear when emulating quantum many-body methods based on a non-Hermitian Hamiltonian. The efficiency and accuracy of these emulators facilitate a rigorous statistical analysis within which we explore nuclear matter predictions for > 10 6 different parametrizations of a chiral interaction model with explicit Δ -isobars at next-to-next-to leading order. Constrained by nucleon-nucleon scattering phase shifts and bound-state observables of light nuclei up to He 4 , we use history matching to identify nonimplausible domains for the low-energy coupling constants of the chiral interaction. Within these domains we perform a Bayesian analysis using sampling and importance resampling with different likelihood calibrations and study correlations between interaction parameters, calibration observables in light nuclei, and nuclear matter saturation properties. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2022

Ab initio Screening of Refractory Nitrides and Carbides for High Temperature Hydrogen Permeation Barriers

Density functional theory was used to screen eleven refractory materials – two pure metals, six nitrides, and three carbides–as high-temperature hydrogen permeation barriers to prevent hydrogen embrittlement. Activation energies were calculated for atomic hydrogen (H) diffusion into the first subsurface layer from the lowest energy surface of the high-temperature phase of each candidate material. The candidate barrier materials with the highest activation energies are h-BN, c-BN, HfN, and ZrN with predicted barriers of 3.25 eV, 3.23 eV, 3.14 eV, and 2.76 eV, respectively. Strain energies, Bader charges, and density of states were calculated for the diffusing H at the relaxed initial state and the transition state to provide insight into contributing factors to high energy barriers. The diffusing H atom in materials with the highest predicted barriers are protic. In addition, interstitial H atoms induce mid-gap states in the density of states of both BN polymorphs. The nitrogen retention of each nitride material at high temperatures was predicted using nitrogen vacancy formation energies with respect to gaseous nitrogen. Experimental evaluation of nitrogen retention in h-BN, ZrN, and TiN confirmed their resistance to nitrogen loss at 1773 K. However, of these nitrides, TiN is predicted to be the least stable. This work identifies multiple promising materials that are predicted to be effective hydrogen barriers at high temperatures and that are stable at temperatures above 2700 K, with BN predicted to perform best.

Density Functional Theory (DFT)↗
DOE OSTI2020

Quantifying Uncertainties In Neutron-Alpha Scattering With Chiral Nucleon-Nucleon And Three-Nucleon Forces

We report that modern ab initio theory combined with high-quality nucleon-nucleon (NN) and three-nucleon (3N) interactions from chiral effective field theory (EFT) can provide a predictive description of low-energy light-nuclei reactions relevant for astrophysics and fusion-energy applications. However, the high cost of computations has so far impeded a complete analysis of the uncertainty budget of such calculations. Starting from NN potentials up to fifth order (N 4 LO) combined with leading-order 3N forces, we study how the order-by-order convergence of the chiral expansion and confidence intervals for the 3N contact and contact-plus-one-pion-exchange low-energy constants (c E and c D ) contribute to the overall uncertainty budget of many-body calculations of neutron- 4 He (n-α) elastic scattering. We compute structure and reaction observables for three-, four-, and five-nucleon systems within the ab initio frameworks of the no-core shell model and no-core shell model with continuum. Using a small set of design runs, we construct a Gaussian process model (GPM) that acts as a statistical emulator for the theory. With this, we gain insight into how uncertainties in the 3N low-energy constants propagate throughout the calculation and determine the Bayesian posterior distribution of these parameters with Markov-Chain Monte Carlo. We find rapidly converging n-α phase shifts with respect to the chiral order. With the adopted leading-order 3N force, calculations based on the NN interaction at N 4 LO of Entem, Machleidt, and Nosyk are unable to reproduce the experimental phase shifts in the 3/2 - channel within the estimated chiral truncation errors. Closer agreement with empirical data is found when using an older parametrization of the NN interaction at order N 3 LO, and the position and width of the P-wave resonances can be used to reduce the uncertainty of the 3N low-energy constants. The present results point to a lack of spin-orbit strength when the newer parametrization of the chiral NN force up to fifth order is combined with the leading-order 3N force. The inclusion of higher-order 3N-force terms may be required to recover the missing strength. GPMs can act as fast and accurate emulators of ab initio many-body calculations of low-energy scattering and reactions of light nuclei, opening the way to a robust quantification of theoretical uncertainties grounded in the description of the underlying chiral Hamiltonian.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Generalizable machine learning potentials for quantum-accurate predictions of non-equilibrium behavior in 2D materials

Machine learning interatomic potentials (ML-IAPs) are emerging as transformative tools in materials modeling, promising quantum-level accuracy at a fraction of the computational cost. However, their ability to generalize beyond equilibrium configurations and to reliably capture defect- and temperature-driven behavior remains underexplored. Here, we develop and benchmark two state-of-the-art ML-IAPs, Spectral Neighbor Analysis Potential (SNAP) and Allegro, on a comprehensive dataset for monolayer MoSe₂. Using density functional theory (DFT) as the reference, we evaluate their performance in capturing stress–strain behavior, phase transition energetics, defect evolution, edge stability, and fracture toughness. Allegro, a deep equivariant neural network potential, surpasses both SNAP and the classical Tersoff potential in accuracy, efficiency, and transferability. Importantly, both ML potentials accurately reproduce experimental fracture measurements and ab initio predictions of inversion domain formation—phenomena well beyond their training sets. Our findings establish ML-IAPs as viable replacements for traditional force fields in the study of non-equilibrium mechanical phenomena, enabling large-scale, high-fidelity simulations in 2D materials and beyond. In conclusion, this work provides a broadly applicable framework for data-driven modeling of structural and functional transformations under extreme conditions.

2D materials↗
DOE OSTI2024

Real-Time KMC Simulation of Vacancy-Mediated Intermixing in Au@Ag Octahedral Core–Cubic Shell Nanocrystals with Ab Initio-Guided Kinetics

Utilization of core–shell rather than monometallic nanocrystals (NCs) facilitates fine-tuning of NC properties for applications. However, compositional evolution via intermixing can degrade these properties prompting recent experimental studies. We develop an atomistic-level stochastic model for vacancy-mediated intermixing exploiting a formalism which allows incorporation at an ab initio density functional theory level of not just the thermodynamics of vacancy formation, but also relevant diffusion barriers for a vast number of possible local environments (in the core and in the shell, at the interface, and in the intermixed phase). This facilitates a predictive treatment and comprehensive understanding of intermixing on the relevant time scale (e.g., 10 1 –10 3 s). In contrast, previous modeling at the atomistic level utilized only unrealistic generic prescriptions of barriers or employed simplified continuum treatments. For Au@Ag octahedral core–cubic shell NCs, our modeling not only captures the experimentally observed rate or time scale for intermixing of ~100 s at 450 °C for 60 nm NCs, but also elucidates the underlying rate controlling processes and the effective intermixing barrier.

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

Active learning of reactive Bayesian force fields applied to heterogeneous catalysis dynamics of H/Pt

Abstract Atomistic modeling of chemically reactive systems has so far relied on either expensive ab initio methods or bond-order force fields requiring arduous parametrization. Here, we describe a Bayesian active learning framework for autonomous “on-the-fly” training of fast and accurate reactive many-body force fields during molecular dynamics simulations. At each time-step, predictive uncertainties of a sparse Gaussian process are evaluated to automatically determine whether additional ab initio training data are needed. We introduce a general method for mapping trained kernel models onto equivalent polynomial models whose prediction cost is much lower and independent of the training set size. As a demonstration, we perform direct two-phase simulations of heterogeneous H 2 turnover on the Pt(111) catalyst surface at chemical accuracy. The model trains itself in three days and performs at twice the speed of a ReaxFF model, while maintaining much higher fidelity to DFT and excellent agreement with experiment.

42 ENGINEERING↗
DOE OSTI2021

A Chemical Dynamics Study of the Reaction of the Methylidyne Radical (CH, X 2 Π) with Dimethylacetylene (CH 3 CCCH 3 , X 1 A 1g )

The gas-phase reaction of the methylidyne (CH; X 2 Π) radical with dimethylacetylene (CH 3 CCCH 3 ; X 1 A 1g ) was studied at a collision energy of 20.6 kJ mol -1 under single collision conditions with experimental results merged with ab initio calculations of the potential energy surface (PES) and ab initio molecule dynamics (AIMD) simulations. Here, the crossed molecular beams experiment reveals that the reaction proceeds barrierless via indirect scattering dynamics through long-lived C 5 H 7 reaction intermediate(s) ultimately dissociating to C 5 H 6 isomers along with atomic hydrogen with atomic hydrogen predominantly released from the methyl groups as verified by replacing the methylidyne with the D1-methylidyne reactant. AIMD simulations reveal that the reaction dynamics are statistical leading predominantly to p28 (1-methyl-3-methylene-cyclopropene, 13 %) and p8 (1-penten-3-yne, 81 %) plus atomic hydrogen with a significant amount of available energy being channeled into the internal excitation of the polyatomic reaction products.

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

Algorithm advances and applications of time‐dependent first‐principles simulations for ultrafast dynamics

Abstract Far from equilibrium phenomenon is a central theme of contemporary material research. Such phenomenon can exhibit itself in atomic structure and dynamics, but very often it also happens as non‐equilibrium phenomenon in the electronic structure. In ab initio material simulation, density functional theory (DFT) has played an essential role in studying electronic ground state problems. For excited states, besides many‐body perturbation theory, another powerful tool is the time dependent DFT (TDDFT) method. In particular, the real‐time TDDFT (rt‐TDDFT) method can be used to simulate many non‐equilibrium phenomena directly. Here we introduce our works on some algorithm advances based on our recently rt‐TDDFT method. This method uses the plane‐wave basis set, and significantly accelerates its efficiency by increasing the time step from 0.1–1 as in traditional methods to 0.2–0.5 fs. The noncollinear magnetic moments and spin–orbit coupling have also been included in our rt‐TDDFT method. Furthermore, a Boltzmann‐TDDFT algorithm has been developed to solve the hot carrier overheating problem in Ehrenfest dynamics, and a natural orbital branching algorithm has been developed to overcome the mean‐field approximation in Ehrenfest dynamics nuclear trajectory, thus allows stochastic multiple paths in chemical reactions. Utilizing these methods, we have studied the photoinduced ultrafast demagnetization, ultrafast phase transition, energy transfer between plasmon and hot carriers, as well as the high‐energy ion implantation and low‐energy atomic diffusion in semiconductors. We believe the tools as the ones introduced here can enable us to study a wide range of phenomena which are of great interest in modern day material research. This article is categorized under: Structure and Mechanism > Computational Materials Science Electronic Structure Theory > Ab Initio Electronic Structure Methods Electronic Structure Theory > Density Functional Theory

Liu, Wen‐Hao↗
DOE OSTI2023

Ab-initio General-relativistic Neutrino-radiation Hydrodynamics Simulations of Long-lived Neutron Star Merger Remnants to Neutrino Cooling Timescales

Abstract We perform the first 3D ab-initio general-relativistic neutrino-radiation hydrodynamics of a long-lived neutron star merger remnant spanning a fraction of its cooling timescale. We find that neutrino cooling becomes the dominant energy loss mechanism after the gravitational-wave dominated phase (∼20 ms postmerger). Electron flavor antineutrino luminosity dominates over electron flavor neutrino luminosity at early times, resulting in a secular increase of the electron fraction in the outer layers of the remnant. However, the two luminosities become comparable ∼20–40 ms postmerger. A dense gas of electron antineutrinos is formed in the outer core of the remnant at densities ∼10 14.5 g cm −3 , corresponding to temperature hot spots. The neutrinos account for ∼10% of the lepton number in this region. Despite the negative radial temperature gradient, the radial entropy gradient remains positive, and the remnant is stably stratified according to the Ledoux criterion for convection. A massive accretion disk is formed from the material squeezed out of the collisional interface between the stars. The disk carries a large fraction of the angular momentum of the system, allowing the remnant massive neutron star to settle to a quasi-steady equilibrium within the region of possible, stable, rigidly rotating configurations. The remnant is differentially rotating, but it is stable against the magnetorotational instability. Other MHD mechanisms operating on longer timescales are likely responsible for the removal of the differential rotation. Our results indicate the remnant massive neutron star is thus qualitatively different from a protoneutron stars formed in core-collapse supernovae.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Spin-orbit-controlled metal–insulator transition in Sr 2 IrO 4

In the context of correlated insulators, where electron–electron interactions (U) drive the localization of charge carriers, the metal–insulator transition is described as either bandwidth- or filling-controlled. Motivated by the challenge of the insulating phase in Sr 2 IrO 4 , a new class of correlated insulators has been proposed, in which spin–orbit coupling (SOC) is believed to renormalize the bandwidth of the half-filled j eff = 1/2 doublet, allowing a modest U to induce a charge-localized phase. Although this framework has been tacitly assumed, a thorough characterization of the ground state has been elusive. Furthermore, direct evidence for the role of SOC in stabilizing the insulating state has not been established, because previous attempts at revealing the role of SOC have been hindered by concurrently occurring changes to the filling. Here, we overcome this challenge by employing multiple substituents that introduce well-defined changes to the signatures of SOC and carrier concentration in the electronic structure, as well as a new methodology that allows us to monitor SOC directly. Specifically, we study Sr 2 Ir 1-x T x O 4 (T = Ru, Rh) by angle-resolved photoemission spectroscopy, combined with ab initio and supercell tight-binding calculations. This allows us to distinguish relativistic and filling effects, thereby establishing conclusively the central role of SOC in stabilizing the insulating state of Sr 2 IrO 4 . Most importantly, we estimate the critical value for SOC in this system to be λ c = 0.42 ± 0.01 eV, and provide the first demonstration of a spin–orbit-controlled metal–insulator transition.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Unconventional Anomalous Hall Effect in Hexagonal Polar Magnet Y 3 Co 8 Sn 4

A rare realization of unconventional anomalous Hall effect (UAHE) both below and above the magnetic transition temperature (T C ) in a hexagonal noncentrosymmetric magnet Y 3 Co 8 Sn 4 , using a combined experimental and ab initio calculations. Occurrence of such UAHE is mainly attributed to the reciprocal ($\mathscr{KS}$) topology (i.e., the presence of topological Weyl points at/near the Fermi level), along with some contribution from the topological magnetic texture (at low temperatures), as inferred from the measured field-dependent ac susceptibility. The effect of UAHE on the measured transport behavior however evolves differently with temperature above and below T C = 53 K, suggesting different physical mechanism responsible in the two phases. A unique planar ferrimagnetic ordering is found to be the most stable state with ab-plane as the easy plane below T C . The simulated net magnetization and the moment per Co atom agrees fairly well with the experimentally measured values. A reasonably large anomalous Hall conductivity (AHC) is also observed in both the phases (above and below T C ) of the present compound, which is again not so ubiquitous. In conclusion, the results underscore the family of R 3 Co 8 Sn 4 (R = rare earth) polar magnets as a compelling backdrop for exploring the synergy of magnetism and topological non-trivial electronic states, pivotal for spintronic applications.

Density functional theory calculations↗
DOE OSTI2023

CoTe 2 : A Quantum Critical Dirac Metal with Strong Spin Fluctuations

Abstract Quantum critical points separating weak ferromagnetic and paramagnetic phases trigger many novel phenomena. Dynamical spin fluctuations not only suppress the long‐range order, but can also lead to unusual transport and even superconductivity. Combining quantum criticality with topological electronic properties presents a rare and unique opportunity. Here, by means of ab initio calculations and magnetic, thermal, and transport measurements, it is shown that the orthorhombic CoTe 2 is close to ferromagnetism, which appears suppressed by spin fluctuations. Calculations and transport measurements reveal nodal Dirac lines, making it a rare combination of proximity to quantum criticality and Dirac topology.

36 MATERIALS SCIENCE↗