Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density functional calculations”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Developing Machine Learning Interatomic Potential for Fe-Cr-Ni Alloys

Accurate prediction of creep and fatigue behavior of stainless steel at elevated temperatures in hydrogen environment requires fundamental understanding of alloy-hydrogen interaction at cross-scale including bulk lattice and key defects such as vacancies, grain boundaries, surfaces, stacking faults, dislocations, and precipitates. This project aims to predict creep behavior of 347H stainless steel with H using machine learning interatomic potentials based on first-principles density functional theory simulations. The Moment Tensor Potentials platform is adopted for this work since it demonstrates a fine balance between model accuracy and computational efficiency. The potential is well trained based on large amount of high-fidelity density functional theory calculations. The validation is carried out by comparing various important properties including short range order, coefficient of thermal expansion, elastic properties, stacking fault energy, grain boundary energy, and surface energy. This work lays the foundation for reliable atomistic simulation of high temperature hydrogen attack of stainless steel.

density functional theory (DFT)↗
DOE OSTI2024

Cryogenic Photoelectron Spectroscopic and Theoretical Study of the Electronic and Geometric Structures of Undercoordinated Osmium Chloride Anions OsCl n – ( n = 3–5)

A series of anionic transition metal halides OsCl n - (n=3–5) have been investigated using a newly developed, home-constructed, cryogenic anion cluster photoelectron spectroscopy. Here, the target anionic species are generated through collision-induced dissociation in a two-stage ion funnel. The measured vertical detachment energies (VDEs) are 3.48 eV, 4.54 eV, and 4.81 eV for n = 3, 4, and 5, respectively. Density functional theory calculations at the B3LYP-D3(BJ)//aug-cc-pVTZ (-pp) level predict the lowest energy structures of OsCl n - (n = 3–5) to be a quintet triangle, quartet square, and quintet square-based pyramid, respectively. The CCSD(T)-calculated VDEs and corresponding adiabatic detachment energies (ADEs) agree well with our experimental measurements. Analysis of the corresponding frontier molecular orbitals (FMOs) and charge density differences suggests that the d-orbitals of the transition metal Os play a primary role in the single-photon detachment processes, and the detached electrons originating from different molecular orbitals are distinguishable.

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

Violation of the Cauchy–Born rule in multi-principal element alloys

Multi-principal element alloys are a novel class of materials that are formed by combining multiple elements in high concentrations and show exceptional properties compared to conventional alloys. These alloys have high configurational entropy due to inherent atomic disorder. The Cauchy–Born rule is a popular homogenization method for linking atomistic to continuum length scales. In this Letter, we use ab initio density functional theory calculations to report that the Cauchy–Born rule, which holds in conventional alloys in the absence of defects, is not valid in multi-principal element alloys. Further, the violation of the Cauchy–Born rule in these alloys is attributed to the presence of atomic disorder due to which the deformations are inhomogeneous. Our results also show the deviations in deformation are related to the magnitude of the stretch and shear by power laws.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Tunable magnetic order in low-symmetry SeO 3 ligand linked TM 3 (SeO 3 ) 3 H 2 O ( TM =Mn , Co, and Ni) compounds

The ability to design crystal lattices with specific sublattice geometries has long been sought after especially in the context of harnessing magnetic frustration to elicit emergent physics. One approach which has seen some success recently is to design magnetic sublattices chemically through the use of nonmagnetic linker ligands as scaffolding for the magnetic ions. In this work we report the magnetic properties of one such family of materials, the transition metal (TM) selenite hydrates with chemical formula TM 3 (SeO 3 ) 3 H 2 O. These materials link highly distorted TMO 6 octahedra via nonmagnetic [SeO 3 ] 2+ linkers. Studying members with TM = Mn, Co, and Ni we use magnetic susceptibility and neutron powder diffraction to identify antiferromagnetic order in all three compounds with moderate frustration indexes. A comparison of the magnetic structures suggests that while the overall structure of the SeO 3 scaffolding remains unchanged the TM effects changes in the magnetic properties through tuning the internal bonding parameters and the single-ion physics. Changing from Mn to Co is found to be particularly consequential manifesting changes in both the ordered-moment direction and in the direction of the ordering vector. Field-dependent measurements of the susceptibility and heat capacity reveal metamagnetic transitions in both Mn 3 (SeO 3 ) 3 H 2 O and Co 3 (SeO 3 ) 3 H 2 O indicating nearby magnetic ground states accessible under relatively small applied fields. Density functional theory calculations broadly confirm these results, showing both a sensitivity of the magnetic structure to the TM and its local environment. Although no spin liquid behavior is achieved, these results suggest the fruitfulness of such synthesis philosophies and encourage future work to engender higher frustration in these materials via doping, field, pressure, or larger linker ligands.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Predicting Hugoniot equation of state in erythritol with ab initio and reactive molecular dynamics

Erythritol has been proposed as an inert surrogate for developing theoretical and computational models to study aging in energetic materials. In this work, we present a comparison of mechanical and shock properties of erythritol computed using the ReaxFF reactive force field and from ab initio calculations employing density functional theory (DFT). We screened eight different ReaxFF parameterizations, of which the CHO parameters developed for hydrocarbon oxidation provide the most accurate predictions of mechanical properties and the crystal structure of erythritol. Further validation of the applicability of this ReaxFF parameterization for modeling erythritol is demonstrated by comparing predictions of the elastic constants, crystal structure, vibrational density of states, and Hugoniot curves against DFT calculations. The ReaxFF predictions are in close agreement with the DFT simulations for the elastic constants and shock Hugoniot when the crystal is loaded along its c axis but show as much as 30% disagreement in the elastic constants in the ab plane and 12% difference in shock pressures when shocked along the a or b crystal axes. Last, we compare thermomechanical properties predicted from classical molecular dynamics with those calculated using the quasi-harmonic approximation and show that quantum mechanical effects produce large discrepancies in the computed values of heat capacity and thermal expansion coefficients compared with classical assumptions. Combining classical molecular dynamics predictions of mechanical behavior with phonon-based calculations of thermal behaviors, we show that predicted shock-induced temperatures for pressures up to 6.5 GPa do not exceed the pressure-dependent melting point of erythritol.

Hu, Jing↗
DOE OSTI2020

Tuning Fermi Levels in Intrinsic Antiferromagnetic Topological Insulators MnBi 2 Te 4 and MnBi 4 Te 7 by Defect Engineering and Chemical Doping

MnBi 2 Te 4 and MnBi 4 Te 7 are intrinsic antiferromagnetic topological insulators, offering a promising materials platform for realizing exotic topological quantum states. However, high densities of intrinsic defects in these materials not only cause bulk metallic conductivity, preventing the measurement of quantum transport in surface states, but may also affect magnetism and topological properties. In this paper, systematic density functional theory calculations reveal specific material chemistry and growth conditions that determine the defect formation and dopant incorporation in MnBi 2 Te 4 and MnBi 4 Te 7 . The large strain induced by the internal heterostructure promotes the formation of large-size-mismatched antisite defects and substitutional dopants. The results here show that the abundance of antisite defects is responsible for the observed n-type metallic conductivity. Additionally, a Te-rich growth condition is predicted to reduce the bulk free electron density, which is confirmed by experimental synthesis and transport measurements in MnBi 2 Te 4 . Furthermore, Na doping is proposed to be an effective acceptor dopant to pin the Fermi level within the bulk band gap to enable the observation of surface quantum transport. The defect engineering and doping strategies proposed here should stimulate further studies for improving synthesis and for manipulating magnetic and topological properties in MnBi 2 Te 4 , MnBi 4 Te 7 , and related magnetic topological insulators.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Covalent Functionalization of Nickel Phosphide Nanocrystals with Aryl-Diazonium Salts

Covalent functionalization of Ni 2 P nanocrystals was demonstrated using aryl-diazonium salts. Spontaneous adsorption of aryl functional groups was observed, with surface coverages ranging from 20 to 96% depending on the native reactivity of the salt as determined by the aryl substitution pattern. Increased coverage was possible for low reactivity species using a sacrificial reductant. Functionalization was confirmed using thermogravimetric analysis, Fourier transform infrared spectroscopy, and X-ray photoelectron spectroscopy. The structure and energetics of this nanocrystal electrocatalyst system, as a function of ligand coverage, were explored with density functional theory calculations. The Hammett parameter of the surface functional group was found to linearly correlate with the change in Ni and P core-electron binding energies and the nanocrystal's experimentally and computationally determined work function. The electrocatalytic activity and stability of the functionalized nanocrystals for hydrogen evolution were also improved when compared to the unfunctionalized material, but a simple trend based on electrostatics was not evident. Density functional theory was used to understand this discrepancy, revealing that H adsorption energies on the covalently functionalized Ni 2 P also do not follow the electrostatic trend and are predictive descriptors of the experimental results.

14 SOLAR ENERGY↗
DOE OSTI2021

Double sulfur vacancies by lithium tuning enhance CO2 electroreduction to n-propanol

Electrochemical CO 2 reduction can produce valuable products with high energy densities but the process is plagued by poor selectivities and low yields. Propanol represents a challenging product to obtain due to the complicated C 3 forming mechanism that requires both stabilization of *C 2 intermediates and subsequent C 1 –C 2 coupling. Herein, density function theory calculations revealed that double sulfur vacancies formed on hexagonal copper sulfide can feature as efficient electrocatalytic centers for stabilizing both CO* and OCCO* dimer, and further CO–OCCO coupling to form C 3 species, which cannot be realized on CuS with single or no sulfur vacancies. The double sulfur vacancies were then experimentally synthesized by an electrochemical lithium tuning strategy, during which the density of sulfur vacancies was well-tuned by the charge/discharge cycle number. The double sulfur vacancy-rich CuS catalyst exhibited a Faradaic efficiency toward n-propanol of 15.4 ± 1% at -1.05 V versus reversible hydrogen electrode in H-cells, and a high partial current density of 9.9 mA cm -2 at -0.85 V in flow-cells, comparable to the best reported electrochemical CO 2 reduction toward n-propanol. Our work suggests an attractive approach to create anion vacancy pairs as catalytic centers for multi-carbon-products.

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

Optical conductivity of overdoped cuprates from ab initio out-of-plane impurity potentials

Dopant impurity potentials determined by ab initio supercell density functional theory calculations are used to calculate the optical conductivity of overdoped La 2-x ⁢Sr x CuO 4 (LSCO) and Tl 2 ⁢Ba 2⁢ CuO 6+δ (Tl-2201) in the superconducting and normal states. Vertex corrections are included to account for the effect of forward scattering on two-particle properties. This approach was previously shown to provide good, semiquantitative agreement with measurements of superfluid density in LSCO. Here, in this study, we compare calculations of conductivity with measurements of terahertz conductivity on LSCO using identical impurity, band, and correlation parameters and find similarly good correspondence with experiment. In the process, we delineate the impact of the different disorder mechanisms on single-particle and transport relaxation processes. In particular, we reveal the critical role of apical oxygen vacancies in transport scattering and show that transport relaxation rates in LSCO are significantly reduced when apical oxygen vacancies are annealed out. These considerations are shown to be crucial for understanding the variability of experimental results for overdoped LSCO in samples of nominally identical doping but different types. Finally, we give predictions for Tl-2201 terahertz conductivity experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Precisely computing phonons via irreducible derivatives

Computing phonons from first principles is typically considered a solved problem, yet inadequacies in existing techniques continue to yield deficient results in systems with sensitive phonons. Here, in this study, we circumvent this issue using the lone irreducible derivative (LID) and bundled irreducible derivative (BID) approaches to computing phonons via finite displacements, where the former optimizes precision via energy derivatives and the latter provides the most efficient algorithm using force derivatives. A condition number optimized basis for BID is derived which guarantees the minimum amplification of error. Additionally, a hybrid LID-BID approach is formulated, in which select irreducible derivatives computed using LID replace BID results. We illustrate our approach on two prototypical systems with sensitive phonons: the shape memory alloy AuZn and metallic lithium. Comparing our resulting phonons in the aforementioned crystals to calculations in the literature reveals nontrivial inaccuracies. Our approaches can be fully automated, making them well suited for both niche systems of interest and high-throughput approaches.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Structure, stability, and superconductivity of N-doped lutetium hydrides at kbar pressures

Here, the structure of the material responsible for the room temperature and near ambient pressure superconductivity reported in an N-doped lutetium hydride has not been conclusively determined. Herein, density functional theory calculations are performed in an attempt to uncover what it might be. Guided by a range of strategies including crystal structure prediction and modifications of existing structure types, we present an array of Lu-N-H phases that are dynamically stable at experimentally relevant pressures. Although none of the structures found are thermodynamically stable, and none are expected to remain superconducting above ~17 K at 10 kbar, a number of metallic compounds with fcc Lu lattices – as suggested by the experimental X-ray diffraction measurements of the majority phase – are identified. The system whose calculated equation of states matches best with that measured for the majority phase is fluorite-type LuH 2 , whose 10 kbar superconducting critical temperature was estimated to be 0.1 K using the Allen-Dynes modified McMillan equation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Origin of insulating and nonferromagnetic SrRuO 3 monolayers

The electro-magnetic properties of ultrathin epitaxial ruthenate films have long been the subject of debate. Here we combine experimental with theoretical investigations of (SrTiO 3 ) 5 -(SrRuO 3 ) n -(SrTiO 3 ) 5 (STO 5 -SRO n -STO 5 ) heterostructures with n = 1 and 2 unit cells, including extensive atomic-resolution scanning-transmission-electron-microscopy imaging, electron-energy-loss-spectroscopy chemical mapping, as well as transport and magneto-transport measurements. The experimental data demonstrate that the STO 5 -SRO 2 -STO 5 heterostructure is nearly stoichiometric, metallic, and ferromagnetic with T C ~ 128 K, even though it lacks the characteristic bulk-SRO octahedral tilts and matches the cubic STO structure. In contrast, the STO 5 -SRO 1 -STO 5 heterostructure features Ru-Ti intermixing in the RuO 2 layer, also without octahedral tilts, but is accompanied by a loss of metallicity and ferromagnetism. Density-functional-theory calculations show that stoichiometric n = 1 and n = 2 heterostructures are metallic and ferromagnetic with no octahedral tilts, while non-stoichiometry in the Ru sublattice in the n = 1 case opens an energy gap and induces antiferromagnetic ordering. Furthermore, the results indicate that the observed non-stoichiometry is the cause of the observed loss of metallicity and ferromagnetism in the n = 1 case.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Interwoven atypical quantum states in CeLiBi 2

Here we report the discovery of CeLiBi 2 , the first example of a material in the tetragonal CeTX 2 ( T = transition metal; X = pnictogen) family wherein an alkali cation replaces the typical transition metal. Magnetic susceptibility and neutron powder diffraction measurements are consistent with a crystal-field Γ 6 ground-state Kramers doublet that orders antiferromagnetically below T N = 3.4 K with an incommensurate propagation wave vector k = ( 0, 0.0724(4), 0.5) that generates a nanometric modulation of the magnetic structure. The best model of the ordered state is an elliptical cycloid with Ce moments primarily residing in the ab plane. This is highly unusual, as all other Γ 6 CeTX 2 members order ferromagnetically. Further, we observe an atypical hard-axis metamagnetic transition at 2 T in magnetostriction, magnetization, and resistivity measurements. CeLiBi 2 is a rare example of a highly conductive material with dominant skew scattering leading to a large anomalous Hall effect. Quantum oscillations with five frequencies arise in magnetostriction and magnetic susceptibility data to T = 30 K and μ 0 H = 55 T, which indicate small Fermi pockets of light carriers with effective masses as low as 0.07 m e . Density functional theory calculations indicate that square-net Dirac-like Bi-p bands are responsible for these ultralight carriers. Together, our results show that CeLiBi 2 enables multiple atypical magnetic and electronic properties in a single clean material.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Dominant effect of anharmonicity on the equation of state and thermal conductivity of MgO under extreme conditions

Understanding heat transfer in the earth mantle region is of great scientific interest but highly challenging due to dissimilar conditions from our daily lives. Thermal conductivity of highly pressured and thermally activated materials can be completely different from that at ambient condition. In this study we calculate the equation of state and thermal conductivity of cubic crystalline MgO, a major component of the mantle, in such an environment. It is shown that the material properties are not accurately captured unless rigorous treatment of anharmonicity of phonon vibrational motion is explicitly applied. To predict renormalized phonon dispersion relations, interatomic force constants of high-order expansion are required, especially at elevated temperature. The anharmonicity plays an increasingly substantial role for regulating volume expansion as temperature rises, because the Gibbs free energy of MgO considerably deviates from what was predicted by the quasiharmonic approximation. We find that the lattice thermal conductivity of MgO and the heat transfer of the earth's mantle are reliably estimated only when both the temperature-dependent renormalization of phonon frequencies and four-phonon scattering effects are taken into consideration. We map the characteristics at Geotherm into contour diagrams in which equilibrium volumes and the thermal conductivities are described as functions of temperature and pressure. Furthermore, the thermal conductivity values of MgO in this study are smaller than those in previous theoretical studies due to the additional incorporation of four-phonon scatterings in the calculations. Specifically, deviation between the values in previous studies and in this study becomes larger as thermodynamic variables reach those near the core-mantle boundary. For instance, thermal conductivities of 52.15 and 23.67 Wm -1 K -1 are estimated at ambient and extreme conditions (T = 4000K and P = 135 GPa), respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Influence of External Conditions on the Black-to-Yellow Phase Transition of CsPbI 3 Based on First-Principles Calculations: Pressure and Moisture

All-inorganic CsPbI 3 perovskite has emerged as a promising candidate for next-generation solar cells. However, owing to the poor structural stability of black phase CsPbI 3 perovskite at room temperature, it spontaneously transforms to the yellow, photoinactive, nonperovskite phase at room temperature, limiting further development of CsPbI 3 perovskite solar cells. Here, to better understand the mechanism driving such undesirable phase transformation, we examine the thermodynamics and kinetics associated with γ-to-δ phase transformation using density functional theory calculation. A solid-state nudged elastic band method was employed to find minimum energy pathways assuming a concerted solid-state phase transformation between these two phases. The gas molecules representing atmospheric (H 2 O, O 2 , and N 2 ) and inert (Ar) ambient environments, as well as applied external pressure, were considered to examine the influence of external conditions on the black-to-yellow phase transition. Our calculation reveals that, with increasing pressure, the thermodynamic driving force for converting the γ-phase to the δ-phase increases, whereas a larger activation energy must be overcome for the phase transition to occur. We also investigate the moisture-induced phase transformation with an H 2 O molecule and its dissociated species (H + /OH – ) and demonstrate that the reaction energy barriers can be significantly lowered in the presence of H 2 O or OH – . On the other hand, other nonpolar species (O 2 , N 2 , and Ar) have a negligible effect on the phase transformation kinetics, while they might reduce the thermodynamic driving force for the phase transformation and suppress the undesirable phase transformation. Our theoretical prediction could support recent observations that the γ-to-δ phase transformation occurs rapidly and catalytically in the presence of moisture, whereas in dry argon and dry oxygen atmospheres, the γ-CsPbI 3 remains stable.

14 SOLAR ENERGY↗
DOE OSTI2024

Thermally induced structural evolution and nanoscale interfacial dynamics in Bi-Sb-Te layered nanostructures

Layered chalcogenides, including Bi-Sb-Te ternary alloys and heterostructures, are renowned as thermoelectric and topological insulators and have recently been highlighted as plasmonic building blocks beyond noble metals. Here, we conduct joint in situ transmission electron microscopy and density functional theory calculations to investigate the temperature-dependent nanoscale dynamics and interfacial properties, identifying the role of native defects and edge configurations in the anisotropic sublimation of Bi 2 Te 3 -Sb 2 Te 3 heterostructures and Sb 2-x Bi x Te 3 alloys. We report structural dynamics, including edge evolution, layer-by-layer sublimation, and the formation and coalescence of thermally induced polygonal nanopores. These nanopores are initiated by preferential dissociation of tellurium, reducing thermal stability in heterostructures. Triangular and quasi-hexagonal configurations dominate nanopore structures in heterostructures. Our calculations reveal antisite defects (Te Sb and Te Bi ) as key players in defect-assisted sublimation. These findings enhance our understanding of nanoscale dynamics and assist in designing tunable low-dimensional chalcogenides.

Bi2Te3-Sb2Te3 heterostructure↗
DOE OSTI2023

Machine learning the relationship between Debye temperature and superconducting transition temperature

Recently a relationship between the Debye temperature $Θ_D$ and the superconducting transition temperature $T_c$ of conventional superconductors has been proposed [Esterlis et al., npj Quantum Mater. 3, 59 (2018)]. The relationship indicates that $T_c$ ≤ $AΘ_D$ for phonon-mediated BCS superconductors, with $A$ being a prefactor of order ~ $0.1$. In order to verify this bound, we train machine learning (ML) models with 10 330 samples in the Materials Project database to predict $Θ_D$. Here, by applying our ML models to 9860 known superconductors in the NIMS SuperCon database, we find that the conventional superconductors in the database indeed follow the proposed bound. We also perform first-principles phonon calculations for $\mathrm{H_3S}$ and $\mathrm{LaH_{10}}$ at 200 GPa. The calculation results indicate that these high-pressure hydrides essentially saturate the bound of $T_c$ versus $Θ_D$.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗