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 1,117 records · Page 62

Alternative first-principles calculation of entropy for liquids

Here, w present an alternative method for interpreting the velocity autocorrelation function (VACF) of a fluid with application to extracting the entropy in a manner similar to the methods developed by Lin et al. [J. Chem. Phys. 119, 11792 (2003)] and improved upon by Desjarlais [Phys. Rev. E 88, 062145 (2013)]. The liquid VACF is decomposed into two components, one gas and one solid, and each contribution's entropic portion is calculated. But, we fit both the gas and solid portions of the VACF in the time domain. This approach is applied to a single-component liquid (a two-phase model of liquid Al at the melt line) and two different two-component systems: a superionic-to-superionic (bcc to fcc) phase transition in H 2 O at high temperatures and pressures and a metastable liquid state of MgO. Finally, for all three examples, comparisons to existing results in the literature demonstrate the validity of our alternative.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

First-principles engineering of charged defects for two-dimensional quantum technologies

Charged defects in two-dimensional (2D) materials have emerging applications in quantum technologies such as quantum emitters and quantum computation. The advancement of these technologies requires a rational design of ideal defect centers, demanding reliable computation methods for the quantitatively accurate prediction of defect properties. Here, we present an accurate, parameter-free, and efficient procedure to evaluate the quasiparticle defect states and thermodynamic charge transition levels of defects in 2D materials. Importantly, we solve critical issues that stem from the strongly anisotropic screening in 2D materials, that have so far precluded the accurate prediction of charge transition levels in these materials. Using this procedure, we investigate various defects in monolayer hexagonal boron nitride (h-BN) for their charge transition levels, stable spin states, and optical excitations. We identify C B V N (nitrogen vacancy adjacent to carbon substitution of boron) to be the most promising defect candidate for scalable quantum bit and emitter applications.

2D materials↗

First-principles study of thermal conductivities of uranium aluminides

Here, we investigate the lattice and electronic thermal conductivities of uranium aluminide compounds, i.e., UAl 2 , UAl 3 and UAl 4 , by density functional theory calculations. The lattice and electronic contributions to the total thermal conductivity of uranium aluminides are compared, and the contributions of phonons with different mean free times and modes to the lattice thermal conductivity of UAl 3 are analyzed from the calculated cumulative and spectral thermal conductivities. The effect of Mo alloying and U vacancy on the thermal conductivity of UAl 3 and UAl 4 is studied by considering the elastic phonon scattering due to mass difference. Both Mo alloying and U vacancy are found to significantly reduce the thermal conductivities of ideal UAl 3 and UAl 4 even with very low concentrations. The currently predicted thermal conductivities of uranium aluminide compounds are expected to be useful to the evaluation of the effective thermal conductivity of the interaction layer formed in U-Mo/Al dispersion fuel.

36 MATERIALS SCIENCE↗

First-principles study of carbon capture and storage properties of porous MnO 2 octahedral molecular sieve OMS-5

Based on the experimentally determined framework structure of porous MnO 2 octahedral molecular sieve (OMS)-5, we used density functional theory-based calculations to evaluate the effect of Na + cation on pore dimensionality and structural stability, and the interaction between CO 2 and OMS-5. Furthermore, we quantified the formation energy of one CO 2 /unit tunnel and two CO 2 /unit tunnel, and projected the electronic density of states on the OMS-5 framework, CO 2 molecules, and Na+ cations to reveal their individual contributions and bonding nature. Partial charge densities were also calculated to investigate CO 2 adsorption behavior in the OMS-5. Our studies predict the initial stage and driving force for the adsorption of CO 2 in the OMS-5, guiding the OMS material design for carbon capture and storage applications.

36 MATERIALS SCIENCE↗

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↗

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↗

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↗

First-principles study of coupled effect of ripplocations and S-vacancies in MoS 2

Recent experiments have revealed ripplocations, atomic-scale ripplelike defects on samples of MoS 2 flakes. We use quantum mechanical calculations based on density functional theory to study the effect of ripplocations on the structural and electronic properties of single-layer MoS 2 , and, in particular, the coupling between these extended defects and the most common defects in this material, S-vacancies. We find that the formation of neutral S-vacancies is energetically more favorable in the ripplocation. In addition, we demonstrate that ripplocations alone do not introduce electronic states into the intrinsic bandgap, in contrast to S-vacancies. We study the dependence of the induced gap states on the position of the defects in the ripplocation, which has implications for the experimental characterization of MoS 2 flakes and the engineering of quantum emitters in this material. Our specific findings collectively aim to provide insights into the electronic structure of experimentally relevant defects in MoS 2 and to establish structure-property relationships for the design of MoS 2 -based quantum devices.

36 MATERIALS SCIENCE↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

First-principles search for alloying elements that increase corrosion resistance of Mg with second-phase particles of transition metal impurities

A significant challenge for applications of Mg alloys is their poor corrosion resistance and hence Mg alloys designs with built-in corrosion resistance are of significant interest. Corrosion can result from the coupling of anodic dissolution of Mg and cathodic reduction of water on impurities such as Iron (Fe)-rich second-phase particles. Experiments have shown that small quantities of Arsenic (As) or Germanium (Ge) can inhibit Mg corrosion, possibly slowing the hydrogen evolution reaction (HER) as the cathodic reaction on Fe surfaces. Since a broader experimental search across the periodic table for other Mg corrosion inhibiting elements is unavailable, we designed thermodynamic and HER criteria, and used high-throughput computations to search a pool of 68 elements including As and Ge that can segregate from bulk Mg to surfaces of Fe particles and impede the HER there. Our computational procedure predicts that six p-block elements meet these criteria, and they rank according to their ability to reduce H adsorption energies and the HER rate as follows: As > Ge > Si > Ga > P ≈ Al. Results for As, the most effective corrosion-inhibiting element, and Ge are in qualitative accord with recent experiments. Further, while none of the 68 elements was found to enhance H adsorption, the six p-block elements reduce H adsorption via strong orbital overlap (Pauli repulsion) between their outer-shell orbitals and the s orbitals of H adsorbates. These p-block elements are also found to have the potential to reduce HER on surfaces of Ni second-phase particles in Mg according to the same criteria, but not on surfaces of Cu second-phase particles.

36 MATERIALS SCIENCE↗

Quantification of uncertainties in thermoelectric properties of materials from a first-principles prediction method: An approach based on Gaussian process regression

We present the electron-phonon averaged via Gaussian process regression (EPA-GPR) method, in which the electron-phonon coupling matrix is represented as a function of two energies and is in turn modeled as a Gaussian process. The EPA-GPR method can be used as an efficient method to estimate thermoelectric properties of materials for fast-screening applications, comparable to the original electron-phonon averaged (EPA) method and the electron-phonon averaged via moving-least-squares (EPA-MLS) method. The EPA-GPR method does not require specification of any open parameter, unlike the other EPA-related methods, since all the hyperparameters in the model can be unambiguously estimated within the type II maximum likelihood (ML-II) approximation. Thus, the EPA-GPR method is a parameter-free estimation method. Additionally, the concept of Gaussian processes in the EPA-GPR method allows us to quantify the uncertainty in estimated properties of thermoelectric materials. One can randomly realize the electron-phonon coupling coefficients from the identified Gaussian process, and those realized samples can be further analyzed in the solution process of the semiclassical Boltzmann transport equation for charge carriers. The results of the semiclassical Boltzmann transport equation provide the statistical properties of the thermoelectric properties of interest. The means, standard deviations, histograms, and confidence intervals of the Seebeck coefficient, the electrical conductivity, and the power factor can be constructed and analyzed. Furthermore the proposed EPA-GPR method is applied to a p-type half-Heusler compound, i.e., HfCoSb, as a case example, the results of which clearly present the advantages of the method.

30 DIRECT ENERGY CONVERSION↗