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

Precipitation of dopants on acceptor-doped LaMnO3±δ revealed by defect chemistry from first principles

Perovskite oxides degrade at elevated temperatures while precipitating dopant-rich particles on the surface. A knowledge-based improvement of surface stability requires a fundamental and quantitative understanding of the dopant precipitation mechanism on these materials. We propose that dopant precipitation is a consequence of the variation of dopant solubility between calcination and operating conditions in solid oxide fuel cells (SOFCs) and electrolyzer cells (SOECs). To study dopant precipitation, we use 20% (D = Ca, Sr, Ba)-doped LaMnO3+δ (LDM20) as a model system. We employ a defect model taking input from density functional theory calculations. The defect model considers the equilibration of LDM20 with a reservoir consisting of dopant oxide (DO), peroxide (DO2), and O2 in the gas phase. The equilibrated non-stoichiometry of the A-site and B-site as a function of temperature, T, and oxygen partial pressure, p(O2), reveals three regimes for LDM20: A-site deficient (oxidizing conditions), A-site rich (atmospheric conditions), and near-stoichiometric (reducing conditions). Assuming an initial A/B non-stoichiometry, we compute the dopant precipitation boundaries in a p–T phase diagram. Our model predicts precipitation both under reducing (DO) and under highly oxidizing conditions (DO2). We found precipitation under anodic, SOEC conditions to be promoted by large dopant size, while under cathodic, SOFC conditions precipitation is promoted by initial A-site excess. The main driving forces for precipitation are oxygen uptake by the condensed phase under oxidizing conditions and oxygen release assisted by B-site vacancies under reducing conditions. Possible strategies for mitigating dopant precipitation under in electrolytic and fuel cell conditions are discussed.

Hess, Franziska (ORCID:000000024460643X)↗
DOE OSTI2021

Light–matter interaction of a molecule in a dissipative cavity from first principles

Cavity-mediated light–matter coupling can dramatically alter opto-electronic and physico-chemical properties of a molecule. Ab initio theoretical predictions of these systems need to combine non-perturbative, many-body electronic structure theory-based methods with cavity quantum electrodynamics and theories of open-quantum systems. Here, we generalize quantum-electrodynamical density functional theory to account for dissipative dynamics of the cavity and describe coupled cavity–single molecule interactions in the weak-to-strong-coupling regimes. Specifically, to establish this generalized technique, we study excited-state dynamics and spectral responses of benzene and toluene under weak-to-strong light–matter coupling. By tuning the coupling, we achieve cavity-mediated energy transfer between electronically excited states. This generalized ab initio quantum-electrodynamical density functional theory treatment can be naturally extended to describe cavity-mediated interactions in arbitrary electromagnetic environments, accessing correlated light–matter observables and thereby closing the gap between electronic structure theory, quantum optics, and nanophotonics.

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

First principles reactive simulation for equation of state prediction

The high cost of density functional theory (DFT) has hitherto limited the ab initio prediction of the equation of state (EOS). In this article, we employ a combination of large scale computing, advanced simulation techniques, and smart data science strategies to provide an unprecedented ab initio performance analysis of the high explosive pentaerythritol tetranitrate (PETN). Comparison to both experiment and thermochemical predictions reveals important quantitative limitations of DFT for EOS prediction and thus the assessment of high explosives. In particular, we find that DFT predicts the energy of PETN detonation products to be systematically too high relative to the unreacted neat crystalline material, resulting in an underprediction of the detonation velocity, pressure, and temperature at the Chapman–Jouguet state. The energetic bias can be partially accounted for by high-level electronic structure calculations of the product molecules. Furthermore we demonstrate a modeling strategy for mapping chemical composition across a wide parameter space with limited numerical data, the results of which suggest additional molecular species to consider in thermochemical modeling.

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

Understanding the interaction between carboxylates and coinage metals from first principles

Carboxylate groups have recently been explored as a new type of ligand to protect superatomic copper and silver nanoclusters, but little is known of the interfacial structure and bonding. Here, we employ density functional theory to investigate the interfaces of a model carboxylate group, CH3COO, on the coinage metal surfaces and clusters. We found that μ2-CH3COO is the most preferred binding mode on the three M(111) surfaces (M = Cu, Ag, and Au), while μ3-CH3COO is also stable on Cu(111) and Ag(111). The saturation coverage was found to be about seven CH3COO groups per nm2 for all surfaces. CH3COO has the strongest binding on Cu and weakest on Au. Moving from the flat surfaces to the icosahedral M13 clusters, we found that the eight-electron superatomic [M13(CH3COO)6]− nanoclusters also prefer the μ2-CH3COO mode on the surface. The icosahedral kernel in [Cu13(CH3COO)6]− and [Ag13(CH3COO)6]− was well maintained after geometry optimization, but a larger deformation was found in [Au13(CH3COO)6]−. Given the broad availability and variety of carboxylic acids including amino acids, our work suggests that carboxylate groups could be the next-generation ligands to further expand the universe of atomically precise metal clusters, especially for Cu and Ag.

Liu, Tongyu↗
DOE OSTI2021

Understanding copper diffusion in CuInSe 2 with first-principles based atomistic and continuum models

We investigate the diffusion of copper in CuInSe 2 using thermodynamic and kinetic models based on density functional theory calculations, attempting to reconcile large differences in reported experimental diffusivities. We find that observations of rapid chemical diffusion can be explained by large thermodynamic factors, which we calculate using a compositionally constrained model of intrinsic point defect formation. We further characterize how copper diffusion coefficients depend on material synthesis conditions and exhibit their variation across the CuInSe 2 secondary phase diagram. In doing so, we identify stable off-stoichiometries that are dominated by either vacancy- or interstitial-mediated diffusion mechanisms. These results are employed in the development of a continuum reaction–diffusion model, which we use to simulate experimental depth profiles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

A first-principles investigation of point defect structure and energetics in ThO 2

The structure and energetics of charged point defects in thorium dioxide (ThO 2 ) have been investigated using the density functional theory (DFT) and phonon simulations. DFT simulations were performed under both zero-pressure and constant volume conditions. Termed as the free volume change of the point defects, the change in volume of the supercell has been computed in the zero-pressure case. Supercell expansion was observed with the increase of the (nominal) charge state of anion (O) interstitials and cation (Th) vacancies from neutral to its maximum. On the contrary, contraction of the supercell has been observed with anion vacancies and cation interstitials as the defect charge increases. The supercell volume change with respect to the charge state has been correlated with the resulting defect energetics. It has been observed that, as the defect charge increased, the internal energy and entropy of defect formation of the cation vacancies and anion interstitials were found to increase, while that of the cation interstitials and anion vacancies decreased. Here, the temperature dependence of internal energy and entropy has also been examined. It was found that, as the temperature increases, the internal energies of the formation of cation vacancies and anion interstitials decrease, while those of the cation interstitials and anion vacancies increase. An opposite observation is seen for the entropies of formation defects when above room temperatures.

42 ENGINEERING↗
DOE OSTI2024

Native point defects in HgCdTe infrared detector material: Identifying deep centers from first principles

We investigate the native point defects in the long-wavelength infrared (LWIR) detector material Hg0.75Cd0.25Te using a dielectric-dependent hybrid density functional combined with spin–orbit coupling. Characterizing these point defects is essential as they are responsible for intrinsic doping and nonradiative recombination centers in the detector material. The dielectric-dependent hybrid functional allows for an accurate description of the bandgap (Eg) for Hg1−xCdxTe (MCT) over the entire compositional range, a level of accuracy challenging with standard hybrid functionals. Our comprehensive examination of the native point defects confirms that cation vacancies VHg(Cd) are the primary sources of p-type conductivity in the LWIR material given their low defect formation energies and the presence of a shallow acceptor level (−/0) near the valence-band maximum. In addition to the shallow acceptor level, the cation vacancies exhibit a deep charge transition level (2−/−) situated near the midgap, characteristic of nonradiative recombination centers. Our results indicate that Hg interstitial could also be a deep center in the LWIR MCT through a metastable configuration under the Hg-rich growth conditions. While an isolated Te antisite does not show deep levels, the formation of VHg–TeHg defect complex introduces a deep acceptor level within the bandgap.

Physics↗
DOE OSTI2024

Prediction challenge: First principles simulation of the ultrafast electron diffraction spectrum of cyclobutanone

Computer simulation has long been an essential partner of ultrafast experiments, allowing the assignment of microscopic mechanistic detail to low-dimensional spectroscopic data. However, the ability of theory to make a priori predictions of ultrafast experimental results is relatively untested. Herein, as a part of a community challenge, we attempt to predict the signal of an upcoming ultrafast photochemical experiment using state-of-the-art theory in the context of preexisting experimental data. Specifically, we employ ab initio Ehrenfest with collapse to a block mixed quantum–classical simulations to describe the real-time evolution of the electrons and nuclei of cyclobutanone following excitation to the 3s Rydberg state. The gas-phase ultrafast electron diffraction (GUED) signal is simulated for direct comparison to an upcoming experiment at the Stanford Linear Accelerator Laboratory. Following initial ring-opening, dissociation via two distinct channels is observed: the C3 dissociation channel, producing cyclopropane and CO, and the C2 channel, producing CH2CO and C2H4. Direct calculations of the GUED signal indicate how the ring-opened intermediate, the C2 products, and the C3 products can be discriminated in the GUED signal. We also report an a priori analysis of anticipated errors in our predictions: without knowledge of the experimental result, which features of the spectrum do we feel confident we have predicted correctly, and which might we have wrong?

Chemistry↗
DOE OSTI2024

Relativistic Douglas–Kroll–Hess calculations of hyperfine interactions within first-principles multireference methods

A relativistic magnetic hyperfine interaction Hamiltonian based on the Douglas–Kroll–Hess (DKH) theory up to the second order is implemented within the ab initio multireference methods, including spin–orbit coupling in the Molcas/OpenMolcas package. This implementation is applied to calculate relativistic hyperfine coupling (HFC) parameters for atomic systems and diatomic radicals with valence s or d orbitals by systematically varying active space size in the restricted active space self-consistent field formalism with restricted active space state interaction for spin–orbit coupling. The DKH relativistic treatment of the hyperfine interaction reduces the Fermi contact contribution to the HFC due to the presence of kinetic factors that regularize the singularity of the Dirac delta function in the nonrelativistic Fermi contact operator. This effect is more prominent for heavier nuclei. As the active space size increases, the relativistic correction of the Fermi contact contribution converges well to the experimental data for light and moderately heavy nuclei. The relativistic correction, however, does not significantly affect the spin-dipole contribution to the hyperfine interaction. In addition to the atomic and molecular systems, the implementation is applied to calculate the relativistic HFC parameters for large trivalent and divalent Tb-based single-molecule magnets (SMMs), such as Tb(III)Pc2 and Tb(II)(CpiPr5)2 without ligand truncation using well-converged basis sets. In particular, for the divalent SMM, which has an unpaired valence 6s/5d hybrid orbital, the relativistic treatment of HFC is crucial for a proper description of the Fermi contact contribution. Even with the relativistic hyperfine Hamiltonian, the divalent SMM is shown to exhibit strong tunability of HFC via an external electric field (i.e., strong hyperfine Stark effect).

Chemistry↗
DOE OSTI2024

First-principles study of defects and doping limits in CaO

Calcium oxide (CaO) is a promising host for quantum defects because of its ultrawide bandgap and potential for long spin coherence times. Using hybrid functional calculations, we investigate the intrinsic point defects and how they limit Fermi-level positions and doping in CaO. We find calcium and oxygen vacancies to be the most common intrinsic defects, acting as compensating acceptors and donors, respectively. Oxygen interstitials are also prevailing under O-rich conditions and act as compensating donors. Due to compensation by these defects, O-poor conditions are required to dope CaO n-type, while O-rich conditions are required for p-type doping. We find that, at room temperature, intrinsic CaO can only achieve Fermi-level positions between 1.76 eV above the valence-band maximum (VBM) and 1.73 eV below the conduction-band minimum (CBM). If suitable shallow dopants are found, the allowed range of Fermi levels would increase to between VBM + 0.53 eV and CBM − 0.27 eV and is set by the compensating intrinsic defects. Additionally, we study hydrogen impurities, and show that hydrogen will not only limit p-type doping but can also act as shallow donor when substituting oxygen (HO defects).

Physics↗
DOE OSTI2025

Validation of a comprehensive first-principles-based framework for predicting the performance of future stellarators

This paper presents the validation of the GENE-KNOSOS-Tango framework for recovering both the steady-state plasma profiles in the considered radial domain and selected turbulence trends in a stellarator. This framework couples the gyrokinetic turbulence code GENE, the neoclassical transport code KNOSOS, and the transport solver Tango in a multi-timescale simulation feedback loop. Ion-scale kinetic-electron and electron-scale adiabatic-ion flux-tube simulations were performed to evolve the density and temperature profiles for four OP1.2b W7-X scenarios. The simulated density and temperature profiles showed good agreement with the experimental data using a reasonable set of boundary conditions. Equally important was the reproduction of observed trends for several turbulence properties, such as density fluctuations and turbulent heat diffusivities. Key effects were also touched upon, such as electron-scale turbulence and the neoclassical radial electric field shear. The validation of the GENE-KNOSOS-Tango framework enables credible predictions of physical phenomena in stellarators and reactor performance based on a given set of edge parameters.

Gyrokinetic simulations↗
DOE OSTI2023

A first principles analysis of potential-dependent structural evolution of active sites in Fe-N-C catalysts

Fe-N-C (iron–nitrogen–carbon) electrocatalysts have emerged as potential alternatives to precious metal-based materials for the oxygen reduction reaction (ORR). However, the structure of these materials under electrochemical conditions is not well understood, and their poor stability in acidic environments poses a formidable challenge for successful adoption in commercial fuel cells. To provide molecular-level insights into these complex phenomena, we combine periodic density functional theory (DFT) calculations, exhaustive treatment of coadsorption effects for ORR reaction intermediates, including O and OH, and comprehensive analysis of solvation stabilization effects to construct voltage-dependent ab initio thermodynamic phase diagrams that describe the in situ structure of the active sites. These structures are further linked to activity and stability descriptors that can be compared with experimental parameters such as the half-wave potential for ORR and the onset potential for carbon corrosion and CO 2 evolution. The results indicate that pyridinic Fe sites at zigzag carbon edges, as well as other edge sites, exhibit high activity for ORR compared to sites in the bulk. However, edges neighboring the active sites are prone to instability via overoxidation and consequent site loss. The results suggest that it could be beneficial to synthesize Fe-N-C catalysts with small sizes and large perimeter edge lengths to enhance ORR activity, while voltage fluctuations should be limited during fuel cell operation to prevent carbon corrosion of overoxidized edges.

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

First-principles calculation of excited states of diatomic molecules: a benchmark for the Gutzwiller conjugate gradient minimisation method

We recently proposed the Gutzwiller conjugate gradient minimisation (GCGM) method for efficient and accurate calculation of the ground state total energy of molecular and bulk systems. The GCGM method is developed under the framework of Gutzwiller wave function but goes beyond the commonly adopted Gutzwiller approximation to improve the accuracy and flexibility in treating the correlation effects. In this conference proceeding, we benchmark the GCGM method with the calculation of excited state potential energy curves of three diatomic molecules, namely H 2 , N 2 , and O 2 . Our calculations demonstrate the flexibility and reasonable accuracy of the method.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Mg and native defects in cubic silicon carbide from first principles

The diffusion of Mg defects in 3C-SiC is studied using the density functional theory. Mg has the highest burn-in rate as a transmutant in 3C-SiC when it is placed in high-energy neutron irradiation environment of a fusion reactor. The presence and evolution of transmutant defects impact thermal and mechanical properties of this important structural material. This study is focused on understanding the structure, stability, and evolution of Mg defects and the interaction of Mg with native defects in 3C-SiC. Our calculations of diffusion coefficients for different Mg defects suggest that Mg is likely to diffuse faster in pristine 3C-SiC than in the damaged one, in agreement with earlier experimental observations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

First-principles calculation of gate-tunable ferromagnetism in magic-angle twisted bilayer graphene under pressure

Magic-angle twisted bilayer graphene (MATBG) is notable as a highly tunable platform for investigating strongly correlated phenomena such as unconventional superconductivity and quantum spin liquids, due to easy control of doping level through gating and sensitive dependence of the magic angle on hydrostatic pressure. Experimental observations of correlated insulating states, unconventional superconductivity and ferromagnetism in MATBG indicate that this system exhibits rich exotic phases. In this work, using density functional theory calculations in conjunction with the effective screening medium method, we find the MATBG under pressure at a twisting angle of 2.88 °and simulate how its electronic states evolve when doping level and electric field perpendicular to plane are tuned by gating. Our calculations show that, at doping levels between two electrons and four holes per moiré unit cell, a ferromagnetic (FM) solution with spin density localized at AA stacking sites is lower in energy than the nonmagnetic solution. The magnetic moment of this FM state decreases with both electron and hole doping and vanishes at four electrons/holes doped per moiré unit cell. Hybridization between the flat bands at the Fermi level and the surrounding dispersive bands can take place at finite doping. On increasing the out-of-plane electric field at zero doping, a transition from the FM state to the nonmagnetic one is seen. An investigation of impurity effects shows that both absorption of O 2 molecules and occurrence of Stone–Wales impurities suppress the FM state, and the mechanisms are understood from our calculations. We also analyze the interlayer bonding character due to flat bands via Wannier functions. Finally, we report trivial band topology of the flat bands in the FM state at a certain doping level.

Physics↗
DOE OSTI2026

First principles band structure of interacting phosphorus and boron/aluminum δ -doped layers in silicon

Silicon can be heavily doped with phosphorus in a single atomic layer (a δ layer), significantly altering the electronic structure of the conduction bands within the material. Recent progress has also made it possible to further dope silicon with acceptor-based δ layers using either boron or aluminum, making it feasible to create devices with interacting δ layers with opposite polarity. Using density functional theory, we calculate the electronic structure of a phosphorus-based δ layer interacting with a boron or aluminum δ layer, varying the distances between the δ layers. At separations 1 nm and smaller, the dopant potentials overlap and largely cancel each other out, leading to an electronic structure closely mimicking intrinsic silicon. At separations greater than 1 nm, the two δ layers behave independently of one another, with an equivalent electronic structure to a p–n diode with an intrinsic layer taking the place of the depletion region. One mechanism for charge transfer between δ layers at larger distances could be tunneling, where we see a tunneling probability exceeding what would be seen for a standard silicon 1.1 eV triangular barrier, indicating that the interaction between delta layers may enhance tunneling compared to a traditional junction.

delta layers↗