Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “First principles”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

First-principles derivation of magnetic interactions in the triangular quantum spin liquid candidates KYb C h 2 ( C h = S , Se , Te ) and A YbSe 2 ( A = Na , Rb )

The AYb Ch 2 (A=alkali metal, Ch=chalcogen) delafossites are a family of crystals which have a triangular lattice of effective S=1/2 moments on their Yb atoms. They hold great promise for the realization of the triangular quantum spin liquid state because of their defect-minimized growth and the ability to interchange chemical constituencies among the family. Here we use ab initio computations to evaluate the exchange couplings of four realized (and one theoretical) rhombohedral delafossite structures and examine the influence of chemical substitution on promoting the development of a quantum spin liquid state. We find good agreement with experiment regarding the antiferromagnetic nearest-neighbor exchange J 1 , but our calculations underestimate J 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Learning Grain-Boundary Segregation: From First Principles to Polycrystals

The segregation of solute atoms at grain boundaries (GBs) can strongly impact the structural and functional properties of polycrystals. Yet, due to the limited availability of simulation tools to study polycrystals at the atomistic scale (i.e., interatomic potentials), there is a minimal understanding of the variation of solute segregation tendencies across the very complex space of GB microenvironments and the large range of alloys in which it can occur. Here, we develop an algorithmic framework that can directly learn the full spectrum of segregation energies for a metal solute atom in a metal polycrystal from ab initio methods, bypassing the need for alloy interatomic potentials. This framework offers a pathway to a comprehensive catalog of GB solute segregation with quantum accuracy, for the entire alloy space. As an initial demonstration in this pursuit, we build an extensive GB segregation database for aluminum-based alloys across the periodic table, including dozens of alloys for which there are substantially no prior data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Photoinduced electronic excitations drive polymerization of carbon monoxide: A first-principles study

Under pressure, carbon monoxide (CO) transforms into a polymer that can be recovered to ambient conditions. While this transformation can occur without additional stimuli, experimental observations have shown that laser irradiation can induce a similar transformation at reduced pressure. The resulting polymeric phase, which is metastable under ambient conditions, releases energy through decomposition into more stable configurations. Using time-dependent density functional theory and Born–Oppenheimer molecular dynamics simulations, we investigate the mechanism by which electronic excitation facilitates CO polymerization. Our calculations reveal that electronic excitation enhances carbon–carbon bonding, enabling polymerization at pressures significantly lower than those required by conventional compression methods. In conclusion, these findings suggest that a photo-assisted approach could be employed to synthesize novel, potentially energetic materials under less demanding pressure conditions.

Born-Oppenheimer molecular dynamics↗
DOE OSTI2022

Predicting Phonon-Induced Spin Decoherence from First Principles: Colossal Spin Renormalization in Condensed Matter

Developing a microscopic understanding of spin decoherence is essential to advancing quantum technologies. Electron spin decoherence due to atomic vibrations (phonons) plays a special role as it sets an intrinsic limit to the performance of spin-based quantum devices. Two main sources of phonon-induced spin decoherence—the Elliott-Yafet and Dyakonov-Perel mechanisms—have distinct physical origins and theoretical treatments. Here, we show calculations that unify their modeling and enable accurate predictions of spin relaxation and precession in semiconductors. We compute the phonon-dressed vertex of the spin-spin correlation function with a treatment analogous to the calculation of the anomalous electron magnetic moment in QED. We find that the vertex correction provides a giant renormalization of the electron spin dynamics in solids, greater by many orders of magnitude than the corresponding correction from photons in vacuum. Finally, our Letter demonstrates a general approach for quantitative analysis of spin decoherence in materials, advancing the quest for spin-based quantum technologies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

First-Principles Investigation of Near-Surface Divacancies in Silicon Carbide

The realization of quantum sensors using spin defects in semiconductors requires a thorough understanding of the physical properties of the defects in the proximity of surfaces. We report a study of the divacancy (V Si V C ) in 3C-SiC, a promising material for quantum applications, as a function of surface reconstruction and termination with -H, -OH, -F and oxygen groups. Here we show that a V Si V C close to hydrogen-terminated (2 x 1) surfaces is a robust spin-defect with a triplet ground state and no surface states in the band gap and with small variations of many of its physical properties relative to the bulk, including the zero-phonon line and zero-field splitting. However, the Debye-Waller factor decreases in the vicinity of the surface and our calculations indicate it may be improved by strain-engineering. Overall our results show that the V Si V C close to SiC surfaces is a promising spin defect for quantum applications, similar to its bulk counterpart.

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

MoS 2 -supported Au 31 for CO hydrogenation: A first-principle study

While the basal plane of single-layer molybdenum disulphide (MoS 2 ) is inert, attempts have been made to functionalize it chemically through the creation of defects and/or addition of dopants. With nanoparticles as dopants, the authors present density functional theory-based calculations of the hydrogenation of CO on a 31-atom, bilayer Au cluster supported on single-layer MoS 2 (Au31/MoS 2 ). Not surprisingly, the adsorption and reaction of all species involved in the hydrogenation occur at the edges of the cluster—the regions at which the interaction between MoS 2 and the Au cluster is tracked to be the strongest. The authors find two possible mechanisms that lead to the formation of methanol: (1) CO* → CHO* → CH 2 O* → CH 3 O* → CH 3 OH* and (2) CO* → CHO* → HCOH* → CH2OH* → CH3OH*, where * indicates the adsorbed species. Detailed analysis of the reaction pathways, however, does not result in favoring one mechanism over the other. Rather, both mechanisms are facile. In addition, the rate-limiting step in each mechanism is found to be the formation of the formyl radical group (CO* + H* → CHO*).

2D materials↗
DOE OSTI2023

First-principles study of grain-boundary wetting in Fe-Σ5(013)[100] tilt boundary

Here, we report an ab-initio study of grain-boundary wetting (GBW) in three systems (Fe/Zn, Fe/Pb & Fe/Bi) on a Σ5(013)[100] tilt grain boundary. An energy-based approach is pursued and adapted to the study of liquid-metal/solid-metal interfaces. We provide a sound calculation of the wettability parameter S (i.e., defined as the minimum energy difference between the unbroken initial state and final state with liquid-metal/solid-metal interfaces). Solid/liquid interfaces are modelled by ab-initio molecular dynamics, also allowing a study of structural effects in the liquid phase near the solid interface. We show how surface energies for solid-metal/liquid-metal interfaces need a special treatment for the liquid-metal chemical potential due to the change in crystallinity near the interface. The computation of the wettability criterion for the three systems correctly predicts the GBW trend known experimentally for these systems. The computations are compared with small-scale TEM experiments on GBW of two systems (Fe/Pb and Fe/Bi) and it is shown to provide reasonable results. This approach offers new perspectives in predicting GBW trends across various systems.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Ferroelectric switching pathways and domain structure of SrBi 2 ⁢(Ta,Nb) 2 ⁢O 9 from first principles

Several families of layered perovskite oxide ferroelectrics exhibit a coupling between polarization and structural order parameters, such as octahedral rotation distortions. This coupling provides opportunities for novel electric field-based manipulation of material properties, and it also stabilizes complex domain patterns and domain-wall vortices. Among layered perovskites with such coupled orders, the Aurivillius-phase oxides SrBi 2 ⁢B 2 ⁢O 9 (B = Ta,Nb) are well-known for their excellent room-temperature ferroelectric performance. This work combines group theoretic analysis with density functional theory calculations to examine the ferroelectric switching processes of SrBi2⁢B2⁢O9. Low-energy two-step ferroelectric switching paths are identified, with polarization reversal facilitated by structural order-parameter rotations. Analysis of the domain structure reveals how the relative energetics of the coupled order parameters translates into a network of several distinct domain-wall types linked by domain-wall vortex structures. Comparisons are made between the ferroelectric switching and domain structure of SrBi 2 ⁢B 2 ⁢O 9 and those of the layered n = 2 Ruddlesden-Popper hybrid improper ferroelectrics. The results provide insight into how ferroelectric properties may be optimized by engineering the complex crystal structures of Aurivillius-phase oxides.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Impacts of point defects on shallow doping in cubic boron arsenide: A first principles study

Cubic boron arsenide (BAs) stands out as a promising material for advanced electronics, thanks to its exceptional thermal conductivity and ambipolar mobility. However, effective control of p- and n-type doping in BAs poses a significant challenge, mostly as a result of the influence of defects. In the present study, we employed density functional theory (DFT) to explore the impacts of the common point defects and impurities on p-type doping of Be B and Si As , and on n-type doping of Si B and Se As . We found that the most favorable point defects formed by C, O, and Si are C As , O B O As , Si As , C As Si B , and O B Si As , which have formation energies of less than 1.5 eV. While the O impurity detrimentally affects both p- and n-type dopings, C and Si impurities are harmful for n-type dopings, making n-type doping a potential challenge. Interestingly, the antisite defect pair A s B B As benefits both p- and n-type doping. Finally, the doping limitation analysis presented in this study can potentially pave the way for strategic development in the area of BAs-based electronics.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Understanding of three different polyvinylpyrrolidone (PVP) based battery binders blends on graphene surfaces from first principles via DFT simulations

Binders play a crucial role in binding electrodes, current collectors, and conductive agents during battery manufacturing. The selection of the binders and their blends considerably affect various properties such as mechanical durability, adhesion performance, ionic/electronic conductivities, and solid electrolyte interface stability for lithium-ion batteries. In this work, polyvinylidene di-fluoride (PVDF), polyacrylic acid (PAA), polyvinylpyrrolidone (PVP), and lithiated polyacrylic acid (Li-PAA) were selected as the model binders, and their blends were studied on graphene surfaces. Further, among them, PVDF is hardly soluble in water unlike the other binders, so in practice, the blend of PVP:PVDF is not thought to be a proper combination. However, the PVP:PVDF binder has also been included because the purpose of this work is to calculate several electronic properties such as binding energies, intermolecular interaction energies, bond critical points, electron density, and the Laplacian of electron density at critical bonding regions of all these blends at the molecular level via density functional theory (DFT) simulations in order to evaluate and compare how is the interaction strength and bonding type of three different lithium-ion batteries (LIBs) binder blends on graphene surfaces. Most stable binder pairs and their binding mechanisms on graphene surfaces were studied. The delocalization of lithium in Li-PAA was studied with the presence of water.

25 ENERGY STORAGE↗
DOE OSTI2022

Machine learning prediction of the mechanical properties of refractory multicomponent alloys based on a dataset of phase and first principles simulation

In this work, a dataset including structural and mechanical properties of refractory multicomponent alloys was developed by fusing computations of phase diagram (CALPHAD) and density functional theory (DFT). The refractory multicomponent alloys, also named refractory complex concentrated alloys (CCAs) which contain 2–5 types of refractory elements were constructed based on Special Quasi-random Structure (SQS). The phase of alloys was predicted using CALPHAD and the mechanical property of alloys with stable and single body-centered cubic (BCC) at high temperature (over 1,500°C) was investigated using DFT-based simulation. As a result, a dataset with 393 refractory alloys and 12 features, including volume, melting temperature, density, energy, elastic constants, mechanical moduli, and hardness, were produced. To test the capability of the dataset on supporting machine learning (ML) study to investigate the property of CCAs, CALPHAD, and DFT calculations were compared with principal components analysis (PCA) technique and rule of mixture (ROM), respectively. It is demonstrated that the CALPHAD and DFT results are more in line with experimental observations for the alloy phase, structural and mechanical properties. Furthermore, the data were utilized to train a verity of ML models to predict the performance of certain CCAs with advanced mechanical properties, highlighting the usefulness of the dataset for ML technique on CCA property prediction.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Performance of Heterostructural TaC/AlGaN Schottky Diodes Based on First Principles Electronic Structure Properties

Advances in ultra-wide bandgap materials, such as high Al-content AlxGa1-xN (AlGaN), are essential for next generation power electronics, but the requirement for lattice matched substrates is currently a significant obstacle. Recently, conductive TaC has emerged as a promising virtual substrate for AlGaN heteroepitaxy, with wurtzite (0001) AlxGa1-xN lattice-matched to rocksalt (111) TaC at x ~ 0.5. Thus, understanding and controlling the electronic properties of the TaC/AlGaN interface is key for developing technological applications based on TaC/AlGaN devices. Using density functional theory and electronic structure calculations, we here investigate TaC/Al0.5Ga0.5N interfaces, where we include explicit alloy models in the slab calculations. We predict the Schottky barrier height and the electric field discontinuity resulting from interface charges. Considering all possible combinations of (Ta, C) substrate termination, (Al/Ga, N) nucleation, and (Al/Ga, N) polarity, we construct a chemical potential phase diagram to identify the stable interfaces that can be accessed through variation of the synthesis conditions. The predicted interface electronic properties are implemented in device performance simulations to demonstrate a practical design for a strain-free, high-efficiency TaC/AlGaN Schottky diode with a low barrier height and without interface charges, underscoring the potential of TaC as a substrate for ultra-wide bandgap devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Solutes that reduce yield strength anisotropies in magnesium from first principles

Using Labusch-type solid solution strengthening models parameterized with DFT-computed solute-dislocation interaction energies, we perform a computational search for 63 solutes across the periodic table to find those that lower anisotropy ratios (non-basal to basal CRSS) of magnesium potentially increasing its ductility per the von Mises criterion. For this purpose, we compute changes in strength for solutes as a function of composition and temperature, and compute anisotropy ratios for solutes that include both rare earth and non-rare earth elements. Here we specifically focus on solute-dislocation interaction energies in the following DFT-optimized dislocations as representative of three non-basal plastic deformation modes: $\langle c + a \rangle$ edge, (10$\bar1$2) tension twinning edge, and the (10$\bar{1}$1) compression twinning edge. We find that solute-induced changes in non-basal deformation modes can be approximated using a second-order polynomial in the size misfit of the solutes, which permits rapid screening of solutes. Our approach to identify solutes known to improve strengthening incorporates solute solubility, and suggests other solutes that not have been previously explored for strengthening. The 8 rare-earth solutes that our method suggests as the best, ordered by increasing anisotropy ratios at their optimal concentrations, are: Gd, Tb, Dy, Nd, Ho, Er, Tm, and Yb. The 12 non-rare-earth solutes that our method suggests as the best, ordered by increasing anisotropy ratios, are: Y, Mn, Sc, Pb, Ca, Ag, Bi, Tl, Zn, Li, Ga, and Al. Of these, Gd, Nd, Er, Yb, Y, Mn, Ca, Zn, Li, and Al are used in commercial Mg alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Confinement effects on the solvation structure of solvated alkaline metal cations in a single-digit 1T-MoS 2 nanochannel: A first-principles study

Confinement plays an important role in determining ion transport in porous materials, which, in turn, may influence the performance of many energy storage and desalination devices. In this work, we combined density functional theory (DFT) with an implicit solvation model and ab initio molecular dynamics (AIMD) to investigate the effects of nanoconfinement on several solvated alkaline metal cations in a single-digit 1T-MoS 2 nanochannel. Our DFT calculations with a solvation model indicated that cations with stronger hydration energy introduce a higher number of co-intercalated water molecules into the channel, consistent with early experimental observation obtained for MXene (2D transition metal carbide) channels. The predicted optimal water numbers for the cations were then used for AIMD simulations that explicitly include the effects of the solvent. When compared with the cations in bulk solution, our simulations showed that the hydration structure and coordination number (CN) of the solvated cations confined in the MoS 2 channel can be significantly altered. Further, we found that larger cations with weaker hydration energy (K + , Rb + , and Cs + ) exhibited a distinctive CN decrease under confinement, while smaller cations (Li + and Na + ) retained a similar hydration shell as in the bulk solution. More specifically, the hydration shell of large cations (K + , Rb + , and Cs + ) in MoS 2 showed similar features of the coordination angle to the bulk, which suggests the partially broken hydration shell with no geometry change under confinement. Our simulations provided insights into the change of the hydration structure of alkaline metal cations under confinement, which may have important implications on their transport in the 1T-MoS 2 channel.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Effect of Ni deposition on CrI 3 monolayer: First-principles study

In our recent experimental work, we demonstrated that Ni thin film deposition can be used to tune the surface magnetic states in bulk CrI 3 . Here, we present the results of a computational study of Ni deposition on a CrI 3 monolayer. In particular, we compare structural, electronic, and magnetic properties of three different systems: a pristine CrI 3 monolayer, a CrI 3 monolayer with partial Ni coverage, and a CrI 3 monolayer with full Ni coverage. The corresponding cells used in our calculations are Cr 2 I 6 , Cr 2 I 6 Ni 1 , and Cr 2 I 6 Ni 3 . In Cr 2 I 6 Ni 1 , the Ni atom is positioned over the hollow site, which corresponds to the lowest energy configuration. In Cr 2 I 6 Ni 3 , the three Ni atoms are positioned over the single hollow site, and over the two Cr sites. Our calculations indicate that the partial coverage of CrI 3 monolayer with Ni retains the semiconducting nature of the parent compound, but results in a strong reduction of the energy band gap. At the same time, the full coverage of CrI 3 monolayer with Ni results in a metallic transition, with a highly spin-polarized (possibly half-metallic) electronic structure. The magnetic alignment of the CrI 3 monolayer with the deposited Ni is ferrimagnetic, due to the anti-aligned magnetic moments of Cr and Ni atoms. While the partial coverage of the CrI 3 monolayer with Ni has a moderate impact on the lattice parameter of the parent compound, full coverage by Ni results in a significant increase of the lattice constant. For both partial and full coverages, the lowest energy configuration corresponds to the Ni atoms being approximately aligned with the out-of-plane coordinate of the surface iodine layer. In addition, in the full coverage mode, the Cr layer exhibits a significant shift away from the Ni atoms, effectively “detaching” the Ni / I atomic layer from the rest of the cell.

2D magnetism↗
DOE OSTI2022

A combined first principles study of the structural, magnetic, and phonon properties of monolayer CrI3

The first magnetic 2D material discovered, monolayer (ML) CrI3, is particularly fascinating due to its ground state ferromagnetism. However, because ML materials are difficult to probe experimentally, much remains unresolved about ML CrI3’s structural, electronic, and magnetic properties. Here, we leverage Density Functional Theory (DFT) and high-accuracy Diffusion Monte Carlo (DMC) simulations to predict lattice parameters, magnetic moments, and spin–phonon and spin–lattice coupling of ML CrI3. We exploit a recently developed surrogate Hessian DMC line search technique to determine CrI3’s ML geometry with DMC accuracy, yielding lattice parameters in good agreement with recently published STM measurements—an accomplishment given the ~10% variability in previous DFT-derived estimates depending upon the functional. Strikingly, we find that previous DFT predictions of ML CrI3’s magnetic spin moments are correct on average across a unit cell but miss critical local spatial fluctuations in the spin density revealed by more accurate DMC. DMC predicts that magnetic moments in ML CrI3 are 3.62 μB per chromium and -0.145 μB per iodine, both larger than previous DFT predictions. The large disparate moments together with the large spin–orbit coupling of CrI3’s I-p orbital suggest a ligand superexchange-dominated magnetic anisotropy in ML CrI3, corroborating recent observations of magnons in its 2D limit. We also find that ML CrI3 exhibits a substantial spin–phonon coupling of ~3.32 cm-1. Our work, thus, establishes many of ML CrI3’s key properties, while also continuing to demonstrate the pivotal role that DMC can assume in the study of magnetic and other 2D materials.

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

Jet Noise Physics and Modeling Using First-principles Simulations

An extensive analysis of our jet DNS database has provided for the first time the complex correlations that are the core of many statistical jet noise models, including MGBK. We have also for the first time explicitly computed the noise from different components of a commonly used noise source as proposed in many modeling approaches. Key findings are: (1) While two-point (space and time) velocity statistics are well-fitted by decaying exponentials, even for our low-Reynolds-number jet, spatially integrated fourth-order space/retarded-time correlations, which constitute the noise "source" in MGBK, are instead well-fitted by Gaussians. The width of these Gaussians depends (by a factor of 2) on which components are considered. This is counter to current modeling practice, (2) A standard decomposition of the Lighthill source is shown by direct evaluation to be somewhat artificial since the noise from these nominally separate components is in fact highly correlated. We anticipate that the same will be the case for the Lilley source, and (3) The far-field sound is computed in a way that explicitly includes all quadrupole cancellations, yet evaluating the Lighthill integral for only a small part of the jet yields a far-field noise far louder than that from the whole jet due to missing nonquadrupole cancellations. Details of this study are discussed in a draft of a paper included as appendix A.

Freund, Jonathan B.↗
DOE OSTI2026

Vibrational and electronic properties of Np 2 O 5 from experimental spectroscopy and first principles calculations

High-valence actinide oxides are critical to understanding the behavior of 5f-electrons, yet their structural and electronic properties remain poorly understood due to challenges in synthesis and handling. We report the first Raman spectroscopic study of single-crystalline Np 2 O 5 and the first scanning tunneling spectroscopy (STS) measurement on any neptunium-containing material. Hydrothermally synthesized crystals were structurally verified by X-ray diffraction. Raman spectra revealed sharply resolved vibrational features, including previously unreported low-frequency modes. STS measurements revealed a band gap of 1.5 eV. Density functional theory (DFT) enables vibrational mode assignments, reveals neptunium-dominated low-frequency phonons, oxygen-dominated high-frequency modes, and predicts an indirect band gap of 1.68 eV. This predicted value is in excellent agreement with the experimentally measured STS gap. This combined Raman, DFT, and STS approach provides a robust framework for correlating lattice dynamics and electronic structure in actinide materials, providing benchmark data for Np 2 O 5 , and opening new avenues for probing structure–property relationships in complex f-electron materials.

36 - MATERIALS SCIENCE↗