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 163 records · Page 9

High pressure stability analysis and chemical bonding of Ti{sub 1-x}Zr{sub x}N alloy: A first principle study

First-principles pseudo-potential calculations have been performed to analyze the stability of Ti{sub 1-x}Zr{sub x}N alloy under high pressures. The first order phase transition from B1 to B2 phase has been observed in this alloy at high pressure. The variation of lattice parameter with the change in concentration of Zr atom in Ti{sub 1-x}Zr{sub x}N is also reported in both the phases. The calculations for density of states have been performed to understand the alloying effects on chemical bonding of Ti-Zr-N alloy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Benchmarking the accuracy of coverage-dependent models: adsorption and desorption of benzene on Pt (111) and Pt 3 Sn (111) from first principles

Bimetallic catalysts have demonstrated properties favorable for upgrading biofuel through catalytic hydrodeoxygenation. However, the design and optimization of such bimetallic catalysts requires the ability to construct accurate, predictive models of these systems. To generate a model that predicts the kinetic behavior of benzene adsorbed on Pt (111) and a Pt 3 Sn (111) surface alloy (Pt 3 Sn (111)), the adsorption of benzene was studied for a wide range of benzene coverages on both surfaces using density functional theory (DFT) calculations. The adsorption energy of benzene was found to correlate linearly with benzene coverage on Pt (111) and Pt 3 Sn (111); both surfaces exhibited net repulsive lateral interactions. Through an analysis of the d-band properties of the metal surface, it was determined that the coverage dependence is a consequence of the electronic interactions between benzene and the surface. The linear coverage dependence of the adsorption energy allowed us to quantify the influence of the lateral interactions on the heat of adsorption and temperature programmed desorption (TPD) spectra using a mean-field model. A comparison of our simulated TPD to experiment showed that this mean-field model adequately reproduces the desorption behavior of benzene on Pt (111) and Pt 3 Sn (111). In particular, the TPD correctly exhibits a broadening desorption peak as the initial coverage of benzene increases on Pt (111) as well as and a low temperature desorption peak on Pt 3 Sn (111). Furthermore, due to the sensitivity of the TPD peak temperature to the desorption energy, precise alignment of experimental and theoretical TPD spectra demands an accurate calculation of the adsorption energy. Therefore, an analysis of the effect of the exchange-correlation functional on TPD modeling is presented. Through this work, we show the necessity of incorporating lateral interactions into theoretical models in order to correctly predict experimental behavior.

09 BIOMASS FUELS↗

Origin of structural analogies and differences between the atomic structures of GeSe{sub 4} and GeS{sub 4} glasses: A first principles study

First-principles molecular dynamics simulations based on density functional theory are employed for a comparative study of structural and bonding properties of two stoichiometrically identical chalcogenide glasses, GeSe{sub 4} and GeS{sub 4}. Two periodic cells of 120 and 480 atoms are adopted. Both glasses feature a coexistence of Ge-centered tetrahedra and Se(S) homopolar connections. Results obtained for N = 480 indicate substantial differences at the level of the Se(S) environment, since Ge–Se–Se connections are more frequent than the corresponding Ge–S–S ones. The presence of a more prominent first sharp diffraction peak in the total neutron structure factor of glassy GeS{sub 4} is rationalized in terms of a higher number of large size rings, accounting for extended Ge–Se correlations. Both the electronic density of states and appropriate electronic localization tools provide evidence of a higher ionic character of Ge–S bonds when compared to Ge–Se bonds. An interesting byproduct of these investigations is the occurrence of discernible size effects that affect structural motifs involving next nearest neighbor distances, when 120 or 480 atoms are used.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effect of Fe and Ti Substitution Doping on Magnetic Property of Monolayer CrSi{sub 2}: a First-Principle Investigation

First-principle calculations based on spin-polarized density functional theory were performed to investigate the effect of Fe and Ti substitution doping on magnetic property of monolayer CrSi{sub 2}. The electronic structures, binding energy, magnetic property, total and partial density of states, and spin density of monolayer CrSi{sub 2} are scientifically studied. Calculated binding energy reveals that Fe-doped monolayer CrSi{sub 2} is more stable than Ti-doped monolayer CrSi{sub 2}. The local magnetic moment of Fe and Ti atom all decrease compared with atomic moment in free gas phase due to variation of bond interaction and charge transfer. The density of states and spin-density results indicated that local magnetic moment of Fe atom is larger than Ti atom, leading to total magnetic moment of Fe-doped monolayer CrSi{sub 2} is bigger than Ti-doped monolayer CrSi{sub 2}.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Communication: Surface-to-bulk diffusion of isolated versus interacting C atoms in Ni(111) and Cu(111) substrates: A first principle investigation

Using first principle calculations, we study the surface-to-bulk diffusion of C atoms in Ni(111) and Cu(111) substrates, and compare the barrier energies associated with the diffusion of an isolated C atom versus multiple interacting C atoms. We find that the preferential Ni-C bonding over C–C bonding induces a repulsive interaction between C atoms located at diagonal octahedral voids in Ni substrates. This C–C interaction accelerates C atom diffusion in Ni with a reduced barrier energy of ∼1 eV, compared to ∼1.4-1.6 eV for the diffusion of isolated C atoms. The diffusion barrier energy of isolated C atoms in Cu is lower than in Ni. However, bulk diffusion of interacting C atoms in Cu is not possible due to the preferential C–C bonding over C–Cu bonding, which results in C–C dimer pair formation near the surface. The dramatically different C–C interaction effects within the different substrates explain the contrasting growth mechanisms of graphene on Ni(111) and Cu(111) during chemical vapor deposition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Near-infrared radiation absorption properties of covellite (CuS) using first-principles calculations

First-principles density functional theory was used to investigate the electronic structure, optical properties and the origin of the near-infrared (NIR) absorption of covellite (CuS). The calculated lattice constant and optical properties are found to be in reasonable agreement with experimental and theoretical findings. The electronic structure reveals that the valence and conduction bands of covellite are determined by the Cu 3d and S 3p states. By analyzing its optical properties, we can fully understand the potential of covellite (CuS) as a NIR absorbing material. Our results show that covellite (CuS) exhibits NIR absorption due to its metal-like plasma oscillation in the NIR range.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Phase Stability in U-6Nb Alloy Doped with Ti from the First Principles Theory

First-principles calculations within the density-functional-theory (DFT) approach are conducted in order to explore and explain the effect of small amounts of titanium on phase stability in the U-6Nb alloy. During rapid quenching from high to room temperature, metastable phases α' (orthorhombic), α" (monoclinic), and γ 0 (tetragonal) can form, depending on Nb concentration. Important mechanical properties depend on the crystal structure and, therefore, an understanding of the effect of impurities on phase stability is essential. Insights on this issue are obtained from quantum-mechanical DFT calculations. The DFT framework does not rely on any material-specific assumptions and is therefore ideal for an unbiased investigation of the U-Nb system.

36 MATERIALS SCIENCE↗

First-principles high-pressure unreacted equation of state and heat of formation of crystal 2,6-diamino-3, 5-dinitropyrazine-1-oxide (LLM-105)

We report dispersion-corrected density functional theoretical calculations of the unreacted equation of state (EOS) of crystal 2,6-diamino-3, 5-dinitropyrazine-1-oxide (LLM-105) under hydrostatic compression of up to 45 GPa. Convergence tests for k-points sampling in the Brillouin zone show that a 3 × 1 × 2 mesh is required to reproduce the X-ray crystal structure at ambient conditions, and we confirm our finding with a separate supercell calculation. Our high-pressure EOS yields a bulk modulus of 19.2 GPa, and indicates a tendency towards anisotropic compression along the b lattice vector due to molecular orientations within the lattice. We find that the electronic energy band gap decreases from a semiconductor type of 1.3 eV at 0 GPa to quasi-metallic type of 0.6 eV at 45 GPa. The extensive intermolecular hydrogen bonds involving the oxide (–NO) and dioxide (–NO{sub 2}) interactions with the amine (–NH{sub 2}) group showed enhanced interactions with increasing pressure that should be discernible in the mid IR spectral region. We do not find evidence for structural phase transitions or chemically induced transformations within the pressure range of our study. The gas phase heat of formation is calculated at the G4 level of theory to be 22.48 kcal/mol, while we obtain 25.92 kcal/mol using the ccCA-PS3 method. Density functional theory calculations of the crystal and the gas phases provided an estimate for the heat of sublimation of 32.4 kcal/mol. We thus determine the room-temperature solid heat of formation of LLM-105 to be −9.9 or −6.5 kcal/mol based on the G4 or ccCA-PS3 methods, respectively.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

First-principles studies of the local structure and relaxor behavior of Pb ( Mg 1 / 3 Nb 2 / 3 ) O 3 – PbTiO 3 -derived ferroelectric perovskite solid solutions

We have investigated the effect of transition-metal dopants on the local structure of the prototypical 0.75Pb(Mg 1/3 Nb 2/3 )O 3 –0.25PbTiO 3 relaxor ferroelectric. We find that these dopants give rise to very different local structure and other physical properties. For example, when Mg is partially substituted by Cu or Zn, the displacement of Cu or Zn is much larger than that of Mg and is even comparable to that of Nb. The polarization of these systems is also increased, especially for the Cu-doped solution, due to the large polarizability of Cu and Zn. As a result, the predicted maximum dielectric constant temperatures T m are increased. On the other hand, the replacement of a Ti atom with a Mo or Tc atom dramatically decreases the displacements of the cations and the polarization, and thus, the T m values are also substantially decreased. Here, the higher T m cannot be explained by the conventional argument based on the ionic radii of the cations. Furthermore, we find that Cu, Mo, or Tc doping increases the cation displacement disorder. Here, the effect of the dopants on the temperature dispersion ΔT m , which is the change in T m for different frequencies, is also discussed. Our findings lay the foundation for further investigations of unexplored dopants.

25 ENERGY STORAGE↗

First principles calculations of point defect diffusion in CdS buffer layers: Implications for Cu(In,Ga)(Se,S){sub 2} and Cu{sub 2}ZnSn(Se,S){sub 4}-based thin-film photovoltaics

We investigate point defects in CdS buffer layers that may arise from intermixing with Cu(In,Ga)Se{sub 2} (CIGSe) or Cu{sub 2}ZnSn(S,Se){sub 4} (CZTSSe) absorber layers in thin-film photovoltaics (PV). Using hybrid functional calculations, we characterize the migration barriers of Cu, In, Ga, Se, Sn, Zn, Na, and K impurities and assess the activation energies necessary for their diffusion into the bulk of the buffer. We find that Cu, In, and Ga are the most mobile defects in CIGS-derived impurities, with diffusion expected to proceed into the buffer via interstitial-hopping and cadmium vacancy-assisted mechanisms at temperatures ∼400 °C. Cu is predicted to strongly favor migration paths within the basal plane of the wurtzite CdS lattice, which may facilitate defect clustering and ultimately the formation of Cu-rich interfacial phases as observed by energy dispersive x-ray spectroscopic elemental maps in real PV devices. Se, Zn, and Sn defects are found to exhibit much larger activation energies and are not expected to diffuse within the CdS bulk at temperatures compatible with typical PV processing temperatures. Lastly, we find that Na interstitials are expected to exhibit slightly lower activation energies than K interstitials despite having a larger migration barrier. Still, we find both alkali species are expected to diffuse via an interstitially mediated mechanism at slightly higher temperatures than enable In, Ga, and Cu diffusion in the bulk. Our results indicate that processing temperatures in excess of ∼400 °C will lead to more interfacial intermixing with CdS buffer layers in CIGSe devices, and less so for CZTSSe absorbers where only Cu is expected to significantly diffuse into the buffer.

36 MATERIALS SCIENCE↗

The high pressure structure and equation of state of 2,6-diamino-3,5-dinitropyrazine-1-oxide (LLM-105) up to 20 GPa: X-ray diffraction measurements and first principles molecular dynamics simulations

Recent theoretical studies of 2,6-diamino-3,5-dinitropyrazine-1-oxide (C{sub 4}H{sub 4}N{sub 6}O{sub 5} Lawrence Livermore Molecule No. 105, LLM-105) report unreacted high pressure equations of state that include several structural phase transitions, between 8 and 50 GPa, while one published experimental study reports equation of state (EOS) data up to a pressure of 6 GPa with no observed transition. Here we report the results of a synchrotron-based X-ray diffraction study and also ambient temperature isobaric-isothermal atomistic molecular dynamics simulations of LLM-105 up to 20 GPa. We find that the ambient pressure phase remains stable up to 20 GPa; there is no indication of a pressure induced phase transition. We do find a prominent decrease in b-axis compressibility starting at approximately 13 GPa and attribute the stiffening to a critical length where inter-sheet distance becomes similar to the intermolecular distance within individual sheets. The ambient temperature isothermal equation of state was determined through refinements of measured X-ray diffraction patterns. The pressure-volume data were fit using various EOS models to yield bulk moduli with corresponding pressure derivatives. We find very good agreement between the experimental and theoretically derived EOS.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Bayesian prior construction for uncertainty quantification in first-principles statistical mechanics

First-principles statistical mechanics enables the prediction of thermodynamic and kinetic properties of materials, but is computationally expensive. Many approaches require surrogate models to calculate energies within Monte Carlo or molecular dynamics simulations. Inexpensive surrogates such as cluster expansions enable otherwise intractable calculations by interpolating data from higher accuracy methods, such as Density Functional Theory (DFT). Surrogate models introduce uncertainty into downstream calculations, in addition to any uncertainty inherent to DFT calculations. Bayesian frameworks address this by quantifying uncertainty and incorporating expert knowledge through priors. However, constructing effective priors remains challenging. This work introduces and describes practical strategies for building Bayesian cluster expansions, focusing on basis truncation, hyperparameter selection, and ground state replication. We analyze multiple basis truncation schemes, compare cross-validation to the evidence-approximation for hyperparameter optimization, and provide methods to find and enforce ground-state-preserving models through priors. Additionally, we compare the uncertainties between different approximations to DFT (LDA, PBE, SCAN) against the uncertainty introduced with the use of cluster expansion surrogate models. These approaches are demonstrated on the BCC Li x Mg 1-x and Li x Al 1-x alloys, which are both of interest for solid-state Li batteries. Our results provide guidelines for constructing and utilizing Bayesian cluster expansions, thereby improving the transparency of materials modeling. Furthermore, the approaches and insights developed in this work can be transferred to a wide range of cluster expansion surrogate models, including the atomic cluster expansion and related machine-learned interatomic potential architectures.

Alloy theory↗

Optical properties of vacancies in aluminum oxide (𝛼−Al 2 ⁢O 3 ) from first principles

We employ first-principles calculations based on hybrid density functional theory to investigate the structural and optical properties of the oxygen vacancy (𝑉 O ) and aluminum vacancy (𝑉 Al ) in 𝛼−Al 2 ⁢O 3 , the most stable (corundum) phase of alumina. Our calculations facilitate the identification of experimental excitation and luminescence spectra with specific electronic transitions at the vacancy sites. The absorption line shape for excitation of an electron at 𝑉$^0_O$ to the conduction-band minimum (CBM) is in excellent agreement with the 6.1 eV band detected by optical absorption spectroscopy, and we find that the 5.9 eV absorption band is generated by an internal electron transition at 𝑉$^0_O$. We confirm that the slowly decaying 3.0 eV emission band of the 𝐹 center is due to a triplet-singlet transition at 𝑉$^0_O$. Our calculations also reveal that the 4.8 eV/5.4 eV absorption and 3.8 eV emission bands assigned to the 𝐹 + center are generated by internal transitions at 𝑉$^{+1}_O$. The line shape for the excitation of an electron from 𝑉$^{+1}_O$ to the CBM agrees well with an absorption band at 6.4 eV, while the recombination of an electron with 𝑉$^{+2}_O$ also produces luminescence around 3.8 eV, but with a considerably broader line shape. For Al vacancies, we confirm that the most prevalent configuration in 𝛼−Al 2 ⁢O 3 is a split-vacancy configuration 𝑉 Al,s , and we calculate migration barriers for different directions in the corundum crystal. We predict absorption and emission spectra for the excitation and recombination of an electron localized at 𝑉$^{−3}_{Al}$ and 𝑉$^{−3}_{Al,s}$ sites with the CBM. The line shapes of the two 𝑉 Al configurations overlap and are considerably broader than the spectra corresponding to 𝑉 O .

36 MATERIALS SCIENCE↗

Improving accuracy of electrochemical capacitance and solvation energetics in first-principles calculations

Here, reliable first-principles calculations of electrochemical processes require accurate prediction of the interfacial capacitance, a challenge for current computationally efficient continuum solvation methodologies. We develop a model for the double layer of a metallic electrode that reproduces the features of the experimental capacitance of Ag(100) in a non-adsorbing, aqueous electrolyte, including a broad hump in the capacitance near the potential of zero charge and a dip in the capacitance under conditions of low ionic strength. Using this model, we identify the necessary characteristics of a solvation model suitable for first-principles electrochemistry of metal surfaces in non-adsorbing, aqueous electrolytes: dielectric and ionic nonlinearity, and a dielectric-only region at the interface. The dielectric nonlinearity, caused by the saturation of dipole rotational response in water, creates the capacitance hump, while ionic nonlinearity, caused by the compactness of the diffuse layer, generates the capacitance dip seen at low ionic strength. We show that none of the previously developed solvation models simultaneously meet all these criteria. We design the nonlinear electrochemical soft-sphere solvation model which both captures the capacitance features observed experimentally and serves as a general-purpose continuum solvation model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Waltzing of a helium pair in tungsten: Migration barrier and trajectory revealed from first-principles

Despite well documented first-principles theoretical determination of the low migration energy (0.06 eV) of a single He in tungsten, fully quantum mechanical calculations on the migration of a He pair still present a challenge due to the complexity of its trajectory. By identifying the six most stable configurations of the He pair in W and decomposing its motion into rotational, translational, and rotational-translational routines, we are able to determine its migration barrier and trajectory. Our density functional theory calculations demonstrate a He pair has three modes of motion: a close or open circular two-dimensional motion in (100) plane with an energy barrier of 0.30 eV, a snaking motion along [001] direction with a barrier of 0.30 eV, and a twisted-ladder motion along [010] direction with the two He swinging in the plane (100) and a barrier of 0.31 eV. The graceful associative movements of a He pair are related to the chemical-bonding-like He-He interaction being much stronger than its migration barrier in W. The excellent agreement with available experimental measurements (0.24–0.32 eV) on He migration makes our first-principles result a solid input to obtain accurate He-W interatomic potentials in molecular dynamics simulations.

36 MATERIALS SCIENCE↗

Effects of cluster expansion on the locations of phase transition boundary as a first step to quantify uncertainty in first principles statistical mechanics framework

Predicting phase diagrams from first principle calculations eliminates the need of tedious experimental trials and errors. Fully automating first principle phase diagram calculations without any sort of human intervention has been a long daunting task and troubling scientific communities for decades. This grand problem remains not fully resolved, largely due to the vastly high-dimensional parameter space associated with density functional theory, cluster expansion, lattice Monte Carlo, and the substantial uncertainty propagating through a set of complex simulations. As a first step to tackle this grand problem, we reported a first demonstration of how sensitive phase boundary locations can be to various cluster expansion fittings and input DFT training data. To the best knowledge of the authors, this study reported the first ever attempt to quantify uncertainty in first principles statistical mechanics framework. In addition, a semi-automated phase transition detection algorithm has been devised in this paper to deal with the associated statistical errors and uncertainties from Monte Carlo method and its predecessor DFT calculations and cluster expansions. This algorithm has been applied in a classical cluster expansion Mg-Cd binary alloy system to detect phase transitions at various locations of phase diagram using different cluster expansions to demonstrate its predictive power and quantify uncertainties. The results suggested that using Chebyshev basis function shifted the transition locations toward the dilute solid solution phase and expanded the phase stability of concentrated ordered phases to wider composition range. The addition of perturbed defect configurations into training data set lowered the order-disorder transition temperature to be away from true transition temperature, suggesting the transition nature is indeed configurational disorder dominated rather than defect assisted. Finally, we found that the weighting has negligible effect on the transition locations, except for the case of using Chebyshev basis function to fit all non-weighted configurations that can be susceptible to Monte Carlo sampling hysteresis.

36 MATERIALS SCIENCE↗

Importance of charge self-consistency in first-principles description of strongly correlated systems

First-principles approaches have been successful in solving many-body Hamiltonians for real materials to an extent when correlations are weak or moderate. As the electronic correlations become stronger often embedding methods based on first-principles approaches are used to better treat the correlations by solving a suitably chosen many-body Hamiltonian with a higher level theory. The success of such embedding theories, often referred to as second-principles, is commonly measured by the quality of self-energy Σ which is either a function of energy or momentum or both. However, Σ should, in principle, also modify the electronic eigenfunctions and thus change the real space charge distribution. While such practices are not prevalent, some works that use embedding techniques do take into account these effects. In such cases, choice of partitioning, of the parameters defining the correlated Hamiltonian, of double-counting corrections, and the adequacy of low-level Hamiltonian hosting the correlated subspace hinder a systematic and unambiguous understanding of such effects. Further, for a large variety of correlated systems, strong correlations are largely confined to the charge sector. Then an adequate nonlocal low-order theory is important, and the high-order local correlations embedding contributes become redundant. Here we study the impact of charge self-consistency within two example cases, TiSe 2 and CrBr 3 , and show how real space charge re-distribution due to correlation effects taken into account within a first-principles Green’s function-based many-body perturbative approach is key in driving qualitative changes to the final electronic structure of these materials.

36 MATERIALS SCIENCE↗

Magnetotransport in semiconductors and two-dimensional materials from first principles

Here we demonstrate a first-principles method to study magnetotransport in materials by solving the Boltzmann transport equation (BTE) in the presence of an external magnetic field. Our approach employs ab initio electron-phonon interactions and takes spin-orbit coupling into account. We apply our method to various semiconductors (Si and GaAs) and two-dimensional (2D) materials (graphene) as representative case studies. The magnetoresistance, Hall mobility, and Hall factor in Si and GaAs are in very good agreement with experiments. In graphene, our method predicts a large magnetoresistance, consistent with experiments. Analysis of the steady-state electron occupations in graphene shows the dominant role of optical phonon scattering and the breaking of the relaxation time approximation. Our paper provides a detailed understanding of the microscopic mechanisms governing magnetotransport coefficients, establishing the BTE in a magnetic field as a broadly applicable first-principles tool to investigate transport in semiconductors and 2D materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗