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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Observing and Modeling the Sequential Pairwise Reactions that Drive Solid‐State Ceramic Synthesis

Abstract Solid‐state synthesis from powder precursors is the primary processing route to advanced multicomponent ceramic materials. Designing reaction conditions and precursors for ceramic synthesis can be a laborious, trial‐and‐error process, as heterogeneous mixtures of precursors often evolve through a complicated series of reaction intermediates. Here, ab initio thermodynamics is used to model which pair of precursors has the most reactive interface, enabling the understanding and anticipation of which non‐equilibrium intermediates form in the early stages of a solid‐state reaction. In situ X‐ray diffraction and in situ electron microscopy are then used to observe how these initial intermediates influence phase evolution in the synthesis of the classic high‐temperature superconductor YBa 2 Cu 3 O 6+ x (YBCO). The model developed herein rationalizes how the replacement of the traditional BaCO 3 precursor with BaO 2 redirects phase evolution through a low‐temperature eutectic melt, facilitating the formation of YBCO in 30 min instead of 12+ h. Precursor selection plays an important role in tuning the thermodynamics of interfacial reactions and emerges as an important design parameter in planning kinetically favorable synthesis pathways to complex ceramic materials.

36 MATERIALS SCIENCE↗
DOE OSTI2026

The CP-PAW Code Package for First-Principles Calculations from a User’s Perspective

CP-PAW is a combined electronic structure and ab initio molecular dynamics code to perform mixed quantum and classical simulations of atomistic condensed phase systems, such as solids, liquids, and molecular systems. As the name suggests, the CP-PAW code unifies the all-electron projector augmented-wave (PAW) method with the Car–Parrinello (CP) approach to determine not only the electronic and nuclear ground states of condensed matter but also to study their properties and dynamics. In addition to briefly outlining the underlying theory, the focus will be on the unique aspects of CP-PAW and how to correctly employ them as a user. How to install CP-PAW using the new build system will also be briefly mentioned.

Blöchl, Peter E [Institute for Theoretical Physic↗
DOE OSTI2024

Unveiling the Electrocatalytic Hydrogen Evolution Reaction Pathway on RuP 2 through Ab Initio Grand Canonical Monte Carlo

In this study, the high catalytic reactivity of ruthenium phosphide (RuP 2 ) has been identified by first-principles density functional theory (DFT) calculations for the electrocatalytic hydrogen evolution reaction (HER). Complex surface reconstructions are considered by applying the ab initio grand canonical Monte Carlo (ai-GCMC) algorithm, efficiently providing a sufficient phase-space exploration of possible surfaces. Combined with surface-phase Pourbaix diagrams, we are able to identify the actual surfaces that obtained under specific experimental environments, thus leading to a more accurate understanding of the nature of the active sites and the binding strength of adsorbates. Specifically, through hundreds of surface reconstructions and hydrogenation states generated with ai-GCMC, we identify the most favorable surface phases of RuP 2 under aqueous acidic conditions. We discover that the HER activity is determined by multiple surfaces with different stoichiometries within a narrow electrode potential window. Low HER overpotential (η) has been found for each of the identified surfaces, as low as 0.04 V. High H-coverage reconstructed surfaces have been discovered under acidic conditions, and the surface Ru sites introduced by additional Ru adatoms or exposed by P-vacancies serve as the active sites for HER based on their nearly reversible H binding. Furthermore, this work provides atomistic insights into the origin of high HER activity on RuP 2 by exploring the dynamic surface phases of electrocatalysts and features a generalizable method to explore the reconstructed/hydrogenated surface space as a function of experimental conditions.

25 ENERGY STORAGE↗
DOE OSTI2025

Conformationally Adaptable Extractant Flexes Strong Lanthanide Reverse-Size Selectivity

Chemical selectivity is traditionally understood in the context of rigid molecular scaffolds with precisely defined local coordination and chemical environments that ultimately facilitate a given transformation of interest. By contrast, nature leverages dynamic structures and strong coupling to enable specific interactions with target species in otherwise complex media. Taking inspiration from nature, we demonstrate unconventional selectivity in the solvent extraction of light over heavy lanthanides using a conformationally flexible ligand called octadecyl acyclopa (ODA). This novel ligand forms pseudocyclic molecular complexes with lanthanide ions at organic/aqueous interfaces, revealed by vibrational sum frequency generation spectroscopy. These complexes are extracted into the organic phase, where femtosecond structural dynamics are probed by two-dimensional infrared spectroscopy and ab initio molecular dynamics simulations to mechanistically frame the macroscopic selectivity trends. We find larger-than-expected structural fluctuations and bond lengths for heavy Ln–ODA complexes that arise from an inability of ODA to contort around the smaller ions to satisfy all would-be bonding interactions, despite forming some individually strong bonds. This finding contrasts with the binding of ODA with lighter lanthanides where, despite individually weaker bonds, collective interactions manifest that minimize structural fluctuations and give rise to enhanced thermodynamic stability. Furthermore, these results point to a new paradigm where conformational dynamics and cumulative bonding interactions can be used to facilitate unconventional chemical transformations.

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

Thermal transport in phase-stabilized lithium zirconate phosphates

The thermal properties of yttrium-stabilized lithium zirconate phosphate [LZP: Li 1+x+y Y x Zr 2-x (PO 4 ) 3 with x = 0.15, -0.2 ≤ y ≤ 0.4 and with x = 0.0, y = 0.0] are presented over a wide temperature range from 30 to 973 K, elucidating the interplay between structural phase transformations and thermal properties in a solid state superionic conducting material. At room temperature, the thermal conductivity decreases by more than 75% as the stoichiometry is changed from lithium deficient to excess and increases with increasing temperature, indicative of defect-mediated transport in the spark plasma sintered materials. The phase transformations and their stabilities are examined by x-ray diffraction and differential scanning calorimetry and indicate that the Y3 + substitution of Zr 4+ is effective in stabilizing the ionically conductive rhombohedral phase over the entire temperature range measured, the mechanism of which is found through ab initio theoretical calculations. Finally, these insights into thermal transport of LZP superionic conductors are valuable as they may be generally applicable for predicting material stability and thermal management in the ceramic electrolyte of future all-solid-state-battery devices.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Smooth periodic gauge satisfying crystal symmetry and periodicity to study high-harmonic generation in solids

Intense lasers can easily drive nonadiabatic transitions of excited electron wave packets across the Brillouin zones, thus transition dipole moments (TDM) between energy bands of solids should be continuous, satisfying crystal symmetry, and periodic at zone boundaries. While current ab initio algorithms are powerful in calculating band structures of solids, they all introduced random phases into the eigenfunctions at each crystal momentum k. In this work, we show how to choose a “smooth-periodic” gauge where TDMs can be smooth versus k, preserving crystal symmetry, as well as maintaining periodic at boundaries. The symmetry properties of TDMs with respect to k ensure the absence of even-order harmonics from MgO with inversion symmetry, while the TDM in the “smooth-periodic” gauge for broken-symmetry ZnO is responsible for even harmonics that were underestimated in previous simulations. These results reveal the importance of correctly treating the complex TDMs that satisfy crystal symmetry and continuous across zone boundaries in nonlinear laser-solid interactions, which has been elusive in most theories so far.

36 MATERIALS SCIENCE↗
NASA NTRS2010

Theoretical Calculation of Jet Fuel Thermochemistry: Tetrahydrodicylopentadiene (JP10) Thermochemistry Using the CBS-QB3 and G3(MP2)//B3LYP Methods - 1

High-level ab initio calculations have been performed on the exo and endo isomers of gas-phase tetrahydrodicyclopentadiene (THDCPD), a principal component of the jet fuel JP10, using the Gaussian Gx and Gx(MPx) composite methods, as well as the CBS-QB3 method, and using a variety of isodesmic and homodesmotic reaction schemes. The impetus for this work is to help resolve large discrepancies existing between literature measurements of the formation enthalpy Delta (sub f)H deg (298) for exo-THDCPD. We find that use of the isodesmic bond separation reaction C10H16 + 14CH4 yields 12C2H6 yields results for the exo isomer (JP10) in between the two experimentally accepted values, for the composite methods G3(MP2), G3(MP2)//B3LYP, and CBS-QB3. Application of this same isodesmic bond separation scheme to gas-phase adamantane yields a value for Delta (sub f)H deg (298) within 5 kJ/mol of experiment. Isodesmic bond separation calculations for the endo isomer give a heat of formation in excellent agreement with the experimental measurement. Combining our calculated values for the gas-phase heat of formation with recent measurements of the heat of vaporization yields recommended values for Delta (sub f)H deg (298)liq of -126.4 and -114.7 kJ/mol for the exo and endo isomers, respectively.

Zehe, Michael J.↗
DOE OSTI2026

Predictive Chemical Kinetic Modeling: Where We Succeed, Where We Struggle, and What Comes Next

Chemical kinetic modeling plays a foundational role in fields ranging from energy to environmental science, pharmaceuticals, and advanced materials. The past two decades have seen remarkable progress, particularly in modeling gas-phase reactions for thermochemical processes, leading to impactful industrial applications such as steam cracking and air quality management. However, new challenges are emerging. The successful development of systematic methodologies for the description of gas-phase kinetics opens the possibility to apply the same approach to the study of more challenging systems. Here, we review recent advances, including ab initio transition state theory-based master equation estimation of elementary rates, automated mechanism generation, machine-learning-assisted kinetics, and uncertainty quantification, and discuss the advances needed to apply the same methodological approach in areas such as heterogeneous catalysis, electrochemistry, liquid-phase and solid-state reactivity, and multiscale model integration. We advocate for the development of targeted tools, especially methods that go beyond empirical tuning toward first-principles-based predictions. We highlight the need for accessible software and AIaugmented workflows to democratize modeling for industry and academia alike. In this perspective, we call attention to not only what has worked but also what remains unsolved, advocating to avoid overemphasizing successes in scientific works at the expense of realism. The next decade should focus on predictive capability, physical accuracy, and community infrastructure (e.g., databases and services) to enable innovation across diverse fields. We argue that kinetic modeling, properly equipped, can accelerate discovery far beyond its traditional domains.

ab initio calculations↗
DOE OSTI2025

Thermodynamic and Kinetic Mechanisms Governing the Synthesis of Nickel-Poor Cathodes

A deeper understanding of the thermodynamics and kinetics governing the lithiation and layering mechanisms of NMC cathode materials (LiNixMnyCozO 2 , where x + y + z = 1) offers valuable insights for enhancing synthesis methods and improving cathode performance. By employing atomistic and mesoscale approaches informed by in situ powder X-ray diffraction (PXRD) experiments, critical parameters for comprehending lithiation and layering processes and reaction rates were identified. The mesoscale approach captured the evolution of the phases and crystallite size observed in the in situ PXRD, revealing the differences in reaction rates with the use of different lithium salts and starting precursors. Ab initio molecular dynamics (AIMD) underscored the importance of vacancies and structural defects in promoting ion mobility and facilitating the nucleation of a layered domain. This nucleation disrupts the symmetry of disordered phases, ultimately creating a strained phase that serves as a buffer between layered and disordered regions. The lithiation and layering processes reflect a dynamic balance between the thermodynamic drive for a low-energy layered structure and the kinetic of diffusion, which is influenced by temperature and lithium vacancy concentration. Overall, reaction mechanisms are driven by the inherent defects of the intermediate phase that differ for NMC cathode materials. The lithium salts impact the rates of lithiation and layering, with a much slower process for Li 2 CO 3 .

25 ENERGY STORAGE↗
DOE OSTI2022

A high-dimensional neural network potential for molecular dynamics simulations of condensed phase nickel and phase transitions

A high-dimensional neural network interatomic potential was developed and used in molecular dynamics simulations of condensed phase Ni and Ni systems with liquid–solid phase coexistence. The reference data set was generated by sampling the potential energy surface over a broad temperature-pressure domain using ab initio MD simulations to train a unified potential. Excellent agreement was achieved between bulk face-centered cubic nickel thermal expansion simulations and relevant experimental data. The same potential also yields accurate structures and diffusivities in the liquid state. The phase transition between liquid and solid phases was simulated using the two-phase interface method. The predicted melting point temperature is within a few kelvins of the literature value. Here, the general methodology could be applied to describe crystals with much more complex phase behaviors.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Thermally frustrated phase transition at high pressure in B 2-ordered FeV

X-ray diffraction measurements of equiatomic B2-ordered FeV were performed in a diamond-anvil cell at room temperature at several pressure points up to 80 GPa that showed the cubic phase to be stable with no indication of structural phase transitions. Density functional theory at 0 K predicts Fermi surface nesting, an electronic topological transition, and a phonon dynamical instability within the experimentally investigated pressure range. Nevertheless, the instability is absent in phonon dispersion curves extracted from ab initio molecular dynamics simulations below the critical volume at temperatures as low as 10 K, indicating that thermal atomic displacements can frustrate the phase transition by renormalizing the phonon dispersion curves. Ferrimagnetism is critical for the stability of the cubic phase at low temperature, but thermal atomic displacements are enough to support the structure at and above the Néel temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Design of high polarization low switching barrier hybrid improper ferroelectric perovskite oxide superlattices

Hybrid improper ferroelectricity is a useful tool to design ABO 3 /A'BO 3 polar superlattices from non polar building blocks. In this study, we have designed high polarization-low switching barrier hybrid improper ferroelectric superlattices with efficient polarization, and polarization-magnetization switching properties above room temperature, using density functional theory and ab initio molecular dynamics simulations. Superlattices with a chemical formula of (AAlO 3 ) m /(A'AlO 3 ) n , where m/n = 1/1, 1/3, 3/1, 1/5 and 5/1, A, A' = Lanthanide and Y cations are considered to outline the design principles behind polarization switching and (LaFeO 3 ) 3 /(CeFeO 3 ) 1 is investigated for polarization-magnetization switching. We find that the unconventional switching paths via out-of-phase rotation Q R- (a 0 a 0 c - ) and tilt precession QTP always yield lower switching barrier compared to those via in-phase rotation Q R+ (a 0 a 0 c + ) and tilt QT (a - a - c 0 ) of BO 6 octahedra. Results from ab initio molecular dynamics simulations estimate the temperature at which the lowest energy barrier can be overcome. Furthermore, it is possible to tune the polarization switching barrier by tuning the tolerance factor, A,A' cation radius mismatch and super lattice periodicity. For switching via Q R- , the switching barrier varies exponentially with rotation angle, indicating how high switching barrier is expected for systems, away from cubic symmetry. We provide a recipe to overcome such a bottleneck by tuning superlattice periodicity. Finally, we have proposed the multiferroic device application concept through a proposed polarization-temperature hysteresis loop and magnetization switching.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Multi-slice electron ptychographic tomography for three-dimensional phase-contrast microscopy beyond the depth of focus limits

Electron ptychography is a powerful computational method for atomic-resolution imaging with high contrast for weakly and strongly scattering elements. Modern algorithms coupled with fast and efficient detectors allow imaging specimens with tens of nanometers thicknesses with sub-0.5 Ångstrom lateral resolution. However, the axial resolution in these approaches is currently limited to a few nanometers, limiting their ability to solve novel atomic structures ab initio. Here, we experimentally demonstrate multi-slice ptychographic electron tomography, which allows atomic resolution three-dimensional phase-contrast imaging in a volume surpassing the depth of field limits. We reconstruct tilt-series 4D-STEM measurements of a $\mathrm{Co_3O_4}$ nanocube, yielding 2 Å axial and 0.7 Å transverse resolution in a reconstructed volume of $\mathrm{(18.2\,nm)^3}$. Our results demonstrate a 13.5-fold improvement in axial resolution compared to multi-slice ptychography while retaining the atomic lateral resolution and the capability to image volumes beyond the depth of field limit. Multi-slice ptychographic electron tomography significantly expands the volume of materials accessible using high-resolution electron microscopy. We discuss further experimental and algorithmic improvements necessary to also resolve single weakly scattering atoms in 3D.

36 MATERIALS SCIENCE↗
DOE OSTI2023

High-throughput ab initio design of atomic interfaces using InterMatch

Forming a hetero-interface is a materials-design strategy that can access an astronomically large phase space. However, the immense phase space necessitates a high-throughput approach for an optimal interface design. Here we introduce a high-throughput computational framework, InterMatch, for efficiently predicting charge transfer, strain, and superlattice structure of an interface by leveraging the databases of individual bulk materials. Specifically, the algorithm reads in the lattice vectors, density of states, and the stiffness tensors for each material in their isolated form from the Materials Project. From these bulk properties, InterMatch estimates the interfacial properties. We benchmark InterMatch predictions for the charge transfer against experimental measurements and supercell density-functional theory calculations. We then use InterMatch to predict promising interface candidates for doping transition metal dichalcogenide MoSe 2 . Finally, we explain experimental observation of factor of 10 variation in the supercell periodicity within a few microns in graphene/α-RuCl 3 by exploring low energy superlattice structures as a function of twist angle using InterMatch. We anticipate our open-source InterMatch algorithm accelerating and guiding ever-growing interfacial design efforts. Moreover, the interface database resulting from the InterMatch searches presented in this paper can be readily accessed online.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

First-Principles-Based Machine Learning Models for Phase Behavior and Transport Properties of CO 2

In this work, we construct distinct first-principles-based machine-learning models of CO 2 , reproducing the potential energy surface of the PBE-D3, BLYP-D3, SCAN and SCAN-rvv10 approximations of density functional theory. Here we employ the Deep Potential methodology to develop the models and consequently achieve a significant computational efficiency over ab initio molecular dynamics (AIMD) that allows for larger system sizes and time scales to be explored. Although our models are trained only with liquid phase configurations, they are able to simulate a stable interfacial system and predict vapor-liquid equilibrium properties, in good agreement with results from the literature. Because of the computational efficiency of the models, we are also able to obtain transport properties, such as viscosity and diffusion coefficients. We find that the SCAN-based model presents a temperature shift in the position of the critical point, while the SCAN-rvv10-based model shows improvement but still exhibits a temperature shift that remains approximately constant for all properties investigated in this work. We find that the BLYP-D3-based model generally performs better for liquid phase and vapor-liquid equilibrium properties, but the PBE-D3-based model is better suited for predicting transport properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2001

Ab-Initio Description and Prediction of Properties of Carbon-Based and Other Non-Metallic Materials

We have resolved the long-standing problem consisting of 30%-50% theoretical underestimates of the band gaps of non-metallic materials. We describe the Bagayoko, Zhao, and Williams (BZW) method that rigorously circumvents the basis-set and variational effect presumed to be a cause of these underestimates. We present ab-initio, computational results that are in agreement with experiment for diamond (C), silicon (Si), silicon carbides (3C-SiC and 4H-SiC), and other semiconductors (GaN, BaTiO3, AlN, ZnSe, ZnO). We illustrate the predictive capability of the BZW method in the case of the newly discovered cubic phase of silicon nitride (c-Si3N4) and of selected carbon nanotabes [(10,0), and (8,4)]. Our conclusion underscores the inescapable need for the BZW method in ab-initio calculations that employ a basis set in a variational approach. Current nanoscale trends amplify this need. We estimate that the potential impact of applications of the BZW method in advancing our understanding of nonmetallic materials, in informing experiment, and particularly in guiding device design and fabrication is simply priceless.

Bagayoko, D.↗
DOE OSTI2025

Charge order driven multiferroic behavior in Sr 4 Fe 6 O 12 : An ab-initio study

Here, in this paper, we report the structural, electronic, and ferroelectric properties of the layered mixed-valent transition-metal compound Sr 4 ⁢ Fe 6 ⁢O 12 (SFO). We demonstrate how SFO undergoes a phase transition from a high-temperature (𝑇) centrosymmetric tetragonal phase (𝑃⁢4 2 /𝑚⁢𝑛⁢𝑚) to a low-𝑇 polar orthorhombic phase (𝑃⁢𝑚⁢𝑛⁢2 1 ). The transition is primarily driven by charge ordering at tetrahedral Fe-layer creating Fe 3+ and Fe 2+ cations between two edge-sharing tetrahedra. This charge ordering induces electronic polarization, which is remarkably larger (3.5 times) in magnitude than ionic polarization and oppositely directed, giving a net polarization of 0.05 C/m 2 which is comparable to the state-of-the-art rare-earth nickelets and manganite perovskites. The direction of structural distortion, governed by the polar mode irrep $Γ^−_5$, depends sensitively on the type of magnetic ordering in the Fe-octahedral layer. Consequently, both the ionic and electronic polarization directions are influenced by magnetic ordering, suggesting the potential for multiferroic behavior with strong magnetoelectric coupling in this material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Mott vs Kondo: Influence of various density functional based methods on the Ce isostructural phase transition mechanism

The cerium iso-structural phase transition (γ–α) is dominated by f-electron localization changes that results in a magnetic ordering change and a volume collapse. Generally, these physics are difficult to capture with ab initio and first principles methods. However, previous works have shown various methods to be successful in predicting at least some of the physics of the γ to α phase transition. Therefore, here, we perform a broad survey of density functional based methods across three levels of theory and types of functions [GGA (generalized gradient approximation), MetaGGA, and hybrid functionals] and compare the results, focusing on hydrostatic compression across the phase boundary at zero kelvin. For the methods that best reproduce experimental results, we directly probe the predicted mechanisms and frame the results in the Mott/Kondo debate, assessing how the underlying methods and assumptions of different functionals can assess the physical drivers in the phase transition, providing insight into the governing processes of this unique phase transition.

36 MATERIALS SCIENCE↗