Wide-ranged multiphase equation of state for iron and model variations addressing uncertainties in high-pressure melting
Not Available
36 MATERIALS SCIENCE↗
SEARCH · Engineering Papers
Results for “first principles density functional theory”
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.
Half-magnetization plateau and the origin of threefold symmetry breaking in an electrically switchable triangular antiferromagnet
Not Available
2-dimensional systems↗
Exploring interlayer coupling in the twisted bilayer PtTe 2
We have investigated interlayer interactions in the bilayer PtTe 2 system, which influence the electronic energy bands near the Fermi level. Our diffusion Monte Carlo (DMC) calculations for the high-symmetry bilayer stackings (AA, AB, AC) manifest distinct interlayer binding characteristics among the stacking modes by revealing significantly different interlayer separations depending on the stacking, which is critical to understanding the interlayer coupling of the twisted bilayers consisting of various local stacking arrangements. Furthermore, a comparison between the interlayer separations obtained from DMC and density functional theory (DFT) shows that meta-generalized gradient approximation (GGA)-based van der Waals–DFT results agree with DMC for different layer stackings, including twisted bilayers, but only the ground-state AA stacking matches well with GGA-based DFT predictions. This underscores the importance of accurate exchange-correlation potentials even for capturing the stacking-dependent interlayer binding properties. We further show that the variability in DFT-predicted interlayer separations is responsible for the large discrepancy of band structures in the 21 . 79 ∘ twisted bilayer PtTe 2 , affecting its classification as metallic or insulating. These results demonstrate the importance of obtaining a correct description of stacking-dependent interlayer coupling in modeling delicate bilayer systems at finite twists. Published by the American Physical Society 2024
36 MATERIALS SCIENCE↗
First-principles equation of state of CHON resin for inertial confinement fusion applications
A wide-range (0 to 1044.0 g/ cm 3 and 0 to 10 9 K) equation-of-state (EOS) table for a CH 1.72 O 0.37 N 0.086 quaternary compound has been constructed based on density-functional theory (DFT) molecular-dynamics (MD) calculations using a combination of Kohn-Sham DFT MD, orbital-free DFT MD, and numerical extrapolation. The first-principles EOS data are compared with predictions of simple models, including the fully ionized ideal gas and the Fermi-degenerate electron gas models, to chart their temperature-density conditions of applicability. The shock Hugoniot, thermodynamic properties, and bulk sound velocities are predicted based on the EOS table and compared to those of C-H compounds. The Hugoniot results show the maximum compression ratio of the C-H-O-N resin is larger than that of CH polystyrene due to the existence of oxygen and nitrogen; while the other properties are similar between CHON and CH. Radiation hydrodynamic simulations have been performed using the table for inertial confinement fusion targets with a CHON ablator and compared with a similar design with CH. Here, the simulations show CHON outperforms CH as the ablator for laser-direct-drive target designs.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
High-precision predictions of properties of chemically disordered crystals
Multiple scattering theory (MST) combined with density functional theory (DFT) allows to predict properties of chemically disordered materials from the first principles. However, such predictions often suffer from the systematic errors, which depend on crystal geometry. Each computed property of a particular crystal structure typically has a relatively small random error and a larger systematic error. Cancellation of systematic errors allows more accurate predictions. We propose a computational methodology based on the subtraction of the systematic errors in MST. To exemplify it, we apply it to the precipitated alloys. Considering precipitation strengthening in Ni superalloys, we compute the relative enthalpies of the competing Ni_3(Al_{1-x}Ti_x)_1 crystal structures with a chemical disorder on the Al+Ti sublattice. Such predicted composition-structure-property dependencies are useful for the guided design of the next-generation alloys with improved strength. Our predictions are validated by comparison with the results of other DFT methods (having a higher computational cost) and with experiment.
density functional theory↗
Extreme sensitivity of higher-order interatomic force constants and thermal conductivity to the energy surface roughness of exchange-correlation functionals
In this Letter, we report that the fourth-order interatomic force constants (4th-IFCs) are significantly sensitive to the energy surface roughness of exchange-correlation (XC) functionals in density functional theory calculations. This sensitivity, which is insignificant for the second- (2nd-) and third-order (3rd-) IFCs, varies for different functionals in different materials and can cause misprediction of thermal conductivity by several times of magnitude. As a result, when calculating the 4th-IFCs using the finite difference method, the atomic displacement needs to be taken large enough to overcome the energy surface roughness, in order to accurately predict phonon lifetime and thermal conductivity. We demonstrate this phenomenon on a benchmark material (Si), a high-thermal conductivity material (BAs), and a low thermal conductivity material (NaCl). For Si, we find that the LDA, PBE, and PBEsol XC functionals are all smooth to the 2nd- and 3rd-IFCs but all rough to the 4th-IFCs. This roughness can lead to a prediction of nearly one order of magnitude lower thermal conductivity. For BAs, all three functionals are smooth to the 2nd- and 3rd-IFCs, and only the PBEsol XC functional is rough for the 4th-IFCs, which leads to a 40% underestimation of thermal conductivity. For NaCl, all functionals are smooth to the 2nd- and 3rd-IFCs but rough to the 4th-IFCs, leading to a 70% underprediction of thermal conductivity at room temperature. In conclusion, with these observations, we provide general guidance on the calculation of 4th-IFCs for an accurate thermal conductivity prediction.
42 ENGINEERING↗
Giant linear magnetoelectric effect at the morphotropic phase boundary of epitaxial Sr 0.5 Ba 0.5 MnO 3 films
First-principles calculations are conducted to compute linear magnetoelectric coupling coefficients in epitaxial (001) Sr 0.5 Ba 0.5 MnO 3 films. A large enhancement of different linear magnetoelectric elements is found in a strained-induced morphotropic phase boundary region. Such enhancement is demonstrated to originate from the behavior of the dielectric susceptibility, thanks to a simple phenomenological model that is presently shown to be relevant and accurate. Finally, this work can thus provide a promising approach towards designing highly desired single-phase multiferroic with a colossal magnetoelectric conversion.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Distinguishing erbium dopants in Y 2 O 3 by site symmetry: Ab initio theory of two spin-photon interfaces
Here, we present a first-principles study of the defect formation and electronic structure of erbium (Er)–doped yttria (Y 2 O 3 ). This is an emerging material for spin-photon interfaces in quantum information science due to the narrow-linewidth optical emission from Er dopants at standard telecommunication wavelengths and their potential for quantum memories and transducers. We calculate formation energies of neutral and negatively and positively charged Er dopants and find the charge-neutral configuration to be the most stable, consistent with experiment. Of the two substitutional sites of Er for Y, the C 2 (more relevant for quantum memories) and C 3i (more relevant for quantum transduction), we identify the former as possessing the lowest formation energy. The electronic properties are calculated using the Perdew-Burke-Ernzerhof functional along with the Hubbard U parameter and spin-orbit coupling, which yields a ~ 6 μ B orbital and a ~ 3 μ B spin magnetic moment, and 11 electrons in the Er 4 f shell, confirming the formation of charge-neutral Er 3+ . This standard density functional theory approach underestimates the band gap of the host and lacks a first-principles justification for U . To overcome these issues we performed screened hybrid functional calculations, including a negative U for the 4 f orbitals, with mixing (α) and screening (w) parameters. These produced robust electronic features with slight modifications in the band gap and the 4 f splittings depending on the choice of tuning parameters. We also computed the many-particle electronic excitation energies and compared them with experimental values from photoluminescence.
36 MATERIALS SCIENCE↗
First-Principles Study of Molecular Hydrogen Activation by Defects in Boron Nitride
Here, we used density functional theory simulations in combination with ab initio thermodynamics to determine the H 2 partial pressure (p H2 )-dependent energetics associated with H 2 activation and recovery at various defect sites in hexagonal boron nitride (h-BN). We found that some defects are very reactive with hydrogen, thereby definitely trapping hydrogen in defective h-BN. However, depending on hydrogen partial pressure, less reactive defect sites can be populated. Because of the lower binding capability of these sites, they would allow hydrogen to be recycled and recovered. For small defect sizes, we found that hydrogen preferentially binds to nitrogen sites by forming N–H bonds, and if no N sites are available then boron sites would be the next to bind hydrogen. Hydrogen dissociation via frustrated Lewis pair is found to be more favorable than forming only N–H bonds but only if the defect size is large enough to accommodate steric effects. For specific conditions such as T = 400 K, p H2 = 1 bar, and only considering one molecular H 2 per defect, three defects, namely, the N monovacancy, 3V(1B2N), and hexagonal 6V(3B3N) could play a role in both the activation and recycling of H 2 as they would be reacting enough to allow a favorable splitting of H 2 while not binding too strongly to allow its recovery. More broadly, a range of p H2 and hydrogen loading conditions were investigated for different types of defects and the finding suggests that p H2 could be used to fine-tune the Gibbs free energy of hydrogenation, thereby allowing several types of defects at different hydrogen loading contents to play a role in the activation/recovery process of H 2 in defective h-BN.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Mechanisms of Ethylene Epoxidation over Silver from Machine Learning-Accelerated First-Principles Modeling and Microkinetic Simulations
In this work we employed machine learning-augmented density functional theory (DFT) thermodynamic calculations to assess the stability of different AgO x structures under catalytic ethylene epoxidation reaction conditions. We found that there are multiple AgO x surface motifs that could co-exist under the relevant conditions. These included Ag surface oxides (e.g., AgO_p(4 × 4) and Ag 1.83 O) and atomic oxygen-covered Ag(111) surfaces. Furthermore, we employed DFT calculations to evaluate the energetics of different reaction mechanisms by which ethylene and oxygen can react on these surfaces. These studies revealed several energetically viable reaction pathways for ethylene epoxidation. Microkinetic modeling analyses, based on the DFT-calculated reaction pathways, showed that ethylene epoxidation can proceed on all surfaces and that multiple pathways, including those involving Langmuir–Hinshelwood and Eley–Rideal mechanisms, could be involved in selective and unselective reactions. The diversity of mechanisms that we discovered in the context of the relatively simple ethylene epoxidation reaction on Ag suggests that the richness and complexity of surface chemistry are most likely a rule rather than an exception in heterogeneous catalytic chemical transformations on metal surfaces and that the concept of a single or even a dominant mechanism and reaction intermediates might need to be revisited for many reactions.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
First-principles study of intrinsic point defects and Xe impurities in uranium monocarbide
Based on density functional theory (DFT) calculations, we perform an extensive investigation of intrinsic point defects and Xe impurities in uranium monocarbide (UC). The DFT calculations involve both the conventional generalized gradient approximation (GGA) and the GGA+U approach with the Hubbard parametric term (U), using up to 5×5×5 supercells. GGA calculations for the formation energy of intrinsic defects demonstrate the significant effect of using larger supercells than in previous studies. Results confirm that the ⟨111⟩ and ⟨100⟩ dumbbell interstitials are the most stable interstitial configurations for U and C, respectively. The interstitial mechanisms are favored for self-diffusion of both uranium and carbon and diffusion of Xe under equilibrium conditions. Calculations also reveal that the Xe substitutional defect at the C lattice site tends to adopt an off-site configuration, which can be interpreted as a Xe interstitial–C vacancy complex. We also utilize GGA+U to assess the impact of effective U parameter (Ueff) on the results. Moreover, we introduce a method to estimate the carbon chemical potential by fitting the phase diagram composition data and propose a selection Ueff=1.25 eV based on the experimental Xe diffusion activation energy. With this approach, GGA+U calculations reproduce the available experimental data for the formation energy of the carbon Frenkel pair and can explain the stepwise recovery of intrinsic properties and burst Xe release behavior in UC observed in annealing experiments.
Huang, Gui-Yang (ORCID:0000000301447376)↗
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↗
Electronic Structure of InAs and InSb Surfaces: Density Functional Theory and Angle-Resolved Photoemission Spectroscopy
The electronic structure of surfaces plays a key role in the properties of quantum devices. However, surfaces are also the most challenging to simulate and engineer. Here, in this work, the electronic structure of InAs(001), InAs(111), and InSb(110) surfaces is studied using a combination of density functional theory (DFT) and angle-resolved photoemission spectroscopy (ARPES). Large-scale first principles simulations are enabled by using DFT calculations with a machine-learned Hubbard U correction [npj Comput. Mater. 6, 180 (2020)]. To facilitate direct comparison with ARPES results, a “bulk unfolding” scheme is implemented by projecting the calculated band structure of a supercell surface slab model onto the bulk primitive cell. For all three surfaces, a good agreement is found between DFT calculations and ARPES. For InAs(001), the simulations clarify the effect of the surface reconstruction. Different reconstructions are found to produce distinctive surface states, which may be detected by ARPES with low photon energies. For InAs(111) and InSb(110), the simulations help elucidate the effect of oxidation. Owing to larger charge transfer from As to O than from Sb to O, oxidation of InAs(111) leads to significant band bending and produces an electron pocket, whereas oxidation of InSb(110) does not. The combined theoretical and experimental results may inform the design of quantum devices based on InAs and InSb semiconductors, for example, topological qubits utilizing the Majorana zero modes.
42 ENGINEERING↗
Revisiting the assignment of atomic charges in metal oxides based on core-level x-ray photoelectron spectra: The case of Ti in SrTiO3(001)
We demonstrate that assigning formal charges to transition metal (TM) cations based on core-level (CL) x-ray photoemission binding energies in oxides leads to physically inconsistent pictures of electronic structure. O 2p–TM 3d hybridization is well known to result in significant covalency in TM–O bonds, thereby reducing TM cation charges from their fully ionic values. However, the ionic bonding model remains the working paradigm for assigning TM CL features, and the resulting cation charges are often taken to be representative of the material under study. Here, we show that a more physically meaningful way to assign charges is to extract information about charge distributions utilizing Dirac–Hartree–Fock theory to calculate CL spectra from first principles and then use the resulting wave functions to determine charges based on orbital occupancies. TM cation charges can also be determined using density functional theory and Bader population analysis. We illustrate these two methods using the Ti 2p spectrum for SrTiO3(001) and show that the agreement between them is excellent. Significantly, the resulting Ti charge is considerably lower than the formal charge. The high degree of similarity between the Ti 2p spectrum for SrTiO3 and those for the rutile and anatase polymorphs of TiO2 suggests that the charge densities surrounding Ti in the latter materials are similar to that in SrTiO3. Taking a broader perspective, oxides containing other first-row transition metals also exhibit covalent character, leading to TM cation charges lower than the analogous fully ionic values in these materials as well.
Chambers, Scott A. (ORCID:000000025415043X)↗
Spin–spin interactions in defects in solids from mixed all-electron and pseudopotential first-principles calculations
Abstract Understanding the quantum dynamics of spin defects and their coherence properties requires an accurate modeling of spin-spin interaction in solids and molecules, for example by using spin Hamiltonians with parameters obtained from first principles calculations. We present a real-space approach based on density functional theory for the calculation of spin-Hamiltonian parameters, where only selected atoms are treated at the all-electron level, while the rest of the system is described with the pseudopotential approximation. Our approach permits calculations for systems containing more than 1000 atoms, as demonstrated for defects in diamond and silicon carbide. We show that only a small number of atoms surrounding the defect needs to be treated at the all-electron level, in order to obtain an overall all-electron accuracy for hyperfine and zero-field splitting tensors. We also present results for coherence times, computed with the cluster correlation expansion method, highlighting the importance of accurate spin-Hamiltonian parameters for quantitative predictions of spin dynamics.
36 MATERIALS SCIENCE↗
Thermodynamics of Plutonium Monocarbide from Anharmonic and Relativistic Theory
Thermodynamics of plutonium monocarbide is studied from first-principles theory that includes relativistic electronic structure and anharmonic lattice vibrations. Density-functional theory (DFT) is expanded to include orbital-orbital coupling in addition to the relativistic spin-orbit interaction for the electronic structure and it is combined with anharmonic, temperature dependent, lattice dynamics derived from the self-consistent ab initio lattice dynamics (SCAILD) method. The obtained thermodynamics are compared to results from simpler quasi-harmonic theory and experimental data. Formation enthalpy, specific heat, and Gibbs energy calculated from the anharmonic model are validated by direct comparison with a calculation of phase diagram (CALPHAD) assessment of PuC and sub-stochiometric PuC0.896. Overall, the theory reproduces CALPHAD results and measured data for PuC rather well, but the comparison is hampered by the sub-stoichiometric nature of plutonium monocarbide. It was found that a bare theoretical approach that ignores spin-orbit and orbital-orbital coupling (orbital polarization) of the plutonium 5f electrons promotes too soft phonons and Gibbs energies that are incompatible with that of the CALPHAD assessment of the experimental data. The investigation of PuC suggests that the electronic structure is accurately described by plutonium 5f electrons as “band like” and delocalized, but correlate through spin polarization, orbital polarization, and spin-orbit coupling, in analogy to previous findings for plutonium metal.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Thermodynamics Modeling for Actinide Monocarbides and Mononitrides from First Principles
The high-temperature thermodynamical properties for the actinide monocarbides and mononitrides ThC, ThN, UC, UN, PuC, and PuN are calculated from first-principles electronic-structure theory. The electronic structure is modeled with density-functional theory (DFT) and is fully relativistic, including the spin-orbit interaction. Furthermore, the DFT is extended to account for orbital–orbital interactions, by means of a parameter-free orbital-polarization (OP) technique, that has proven to be essential for the 5f electrons in plutonium. Strong anharmonicity and the temperature dependence of the lattice vibrations are captured with the self-consistent ab initio lattice dynamics (SCAILD) method. The calculated free energies and heat capacities are compared to published results from quasi-harmonic (QH) theory, and experiments, where available. For the uranium and plutonium compounds, we make use of CALPHAD assessments to help evaluate the theory. Generally, our anharmonic relativistic approach compares well with both CALPHAD and experiments. For the thorium compounds, our theory is in good accord with QH modeling of the free energy at lower temperatures but for the heat capacity the comparison is less favorable.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Electronic structure and the van Hove singularity scenario in high-T(sub c)H(g)Ba2CuO(4+delta) superconductors
The electronic structure and the hole concentrations in the high Tc superconductor HgBa2CuO(4+delta) (delta = O, 1) has been investigated by employing a first principles full potential self-consistent LMTO method with the local density functional theory. The scalar relativistic effects have been considered. The hole concentrations of the Cu-d and O-p(x,y) orbitals are seen to be larger for the HgBaCuO5 system than those of the HgBaCuO4 solid. However, the van Hove singularity (vHs) induced Cu-d and O-p peak which is seen to lie comparatively away and above the Fermi level in the delta = 1 system shifts towards the Fermi level in the delta = 0 system. Thus, the superconducting behavior appears to originate from the occurrence of the vHs peak at the Fermi level. The Fermi surface nesting area in the delta = 0 compound is seen to be larger than in the delta = 1 compound. The calculation reveals that the increase in pressure on the crystal enhances the hole concentrations but without showing any optimum value, On the other hand, the vHs peak approaches to-wards the Fermi level with pressure and crosses the Fermi surface near V/Vo approximately equals 0.625 (V and Vo are the crystal volumes at high and normal pressures, respectively). Our calculated value of the bulk modulus equal to 0.626 Mbar predicts the occurrence of this crossover at about 24 GPa which is in complete agreement with the experimental value. At this pressure the compound has maximum nesting area and self-doped behavior.
Agrawal, Bal K.↗