Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “first principles calculation”

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,549 records · Page 86
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Heavy boron doping in superconducting carbon materials

Here we examine physical properties of heavily boron-doped sp 3 -hybridized carbon allotropes: cubic diamond, hexagonal diamond, and body centered tetragonal C 4 . The structural similarity between cubic diamond and hexagonal diamond leads to similar responses to substitutional boron doping. On the other hand, body centered tetragonal C 4 exhibits distinct structural and electronic properties because of its characteristic structure. Our study also shows that the superconducting transition temperatures up to 60 K are possible in heavily doped carbon materials despite their various atomistic structures.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Machine learning study of magnetism in uranium-based compounds

Actinide and lanthanide-based materials display exotic properties that originate from the presence of itinerant or localized f electrons and include unconventional superconductivity and magnetism, hidden order, and heavy-fermion behavior. Due to the strongly correlated nature of the 5f electrons, magnetic properties of these compounds depend sensitively on applied magnetic field and pressure, as well as on chemical doping. However, precise connection between the structure and magnetism in actinide-based materials is currently unclear. In this investigation, we established such structure-property links by assembling and mining two datasets that aggregate, respectively, the results of high-throughput density functional theory simulations and experimental measurements for the families of uranium- and neptunium-based binary compounds. Various regression algorithms were utilized to identify correlations among accessible attributes (features or descriptors) of the material systems and predict their cation magnetic moments and general forms of magnetic ordering. Descriptors representing compound structural parameters and cation f-subshell occupation numbers were identified as most important for accurate predictions. The best machine learning model developed employs the random forest regression algorithm. It can predict both spin and orbit moment size with root-mean-square error of 0.17 μ B and 0.19 μ B , respectively. Lastly, the random forest classification algorithm is used to predict the ordering (paramagnetic, ferromagnetic, and antiferromagnetic) of such systems with 76% accuracy.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Localized excitons in defective monolayer germanium selenide

Germanium selenide (GeSe) is a van der Waals–bonded layered material with promising optoelectronic properties, which has been experimentally synthesized for two-dimensional (2D) semiconductor applications. In the monolayer, due to reduced dimensionality and, thus, screening environment, perturbations such as the presence of defects have a significant impact on its properties. We apply density functional theory and many-body perturbation theory to understand the electronic and optical properties of GeSe containing a single selenium vacancy in the –2 charge state. We predict that the vacancy results in midgap “trap states” that strongly localize the electron and hole density and lead to sharp, low-energy optical absorption peaks below the predicted pristine optical gap. Analysis of the exciton wave function reveals that the 2D Wannier-Mott exciton of the pristine material monolayer is strongly modified by the presence of the defect. The lowest-energy exciton wave function is highly localized around the defect, with the Bohr radius reduced by a factor of 4 when compared to the pristine material and a dipole moment along the out-of-plane axis due to the defect-induced symmetry breaking. Altogether, these results suggest that the vacancy significantly modifies the excitonic properties of the system, demonstrating the importance of considering defects in the context of material design.

2-dimensional systems↗
DOE OSTI2021

Structural motifs and bonding in two families of boron structures predicted at megabar pressures

The complex crystal chemistry of elemental boron has led to numerous proposed structures with distinctive motifs as well as contradictory findings. Herein, evolutionary structure searches performed at 100 GPa have uncovered a series of metastable phases of boron, and bonding analyses were carried out to elucidate their electronic structure. These polymorphs, dynamically stable at 100 GPa, were grouped into two families. Here, the first was derived from the thermodynamic minimum at these conditions, α-Ga, whereas channels comprised the second. Two additional intergrowth structures were uncovered, and it was shown they could be constructed by stacking layers of α-Ga-like and channel-like allotropes on top of each other. A detailed bonding analysis revealed networks of four-center σ-bonding functions linked by two-center B-B bonds in the α-Ga based structures, and networks that were largely composed of three-center σ-bonding functions in the channel-based structures. Seven of these high-pressure phases were found to be metastable at atmospheric conditions, and their Vickers hardnesses were estimated to ≈ 36 GPa.

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

Antiphase boundary migration as a diffusion mechanism in a P3 sodium layered oxide

The electrode materials of Na-ion batteries must exhibit fast ion transport, a property that is often negatively impacted when mobile cations order. Many layered Na (and K) intercalation compounds adopt the P3 crystal structure, which tends to stabilize cation-vacancy orderings over a wide range of Na concentrations. The ordered phases have been predicted to consist of domains that are periodically separated by antiphase boundaries (APBs) to accommodate variations in Na concentration. The mechanisms and rates with which Na diffuses within such ordered phases of P3 are not understood. We have conducted a first-principles study of Na diffusion in the P3 structure, using Na x CoO 2 as a model system. We find that most simple nearest-neighbor hops in Na-vacancy ordered phases in P3 result in unstable configurations and are therefore forbidden. We identify an alternative mechanism of Na transport that is mediated by APB migration via the formation and expansion of kinks along domain boundaries. The limiting kinetic barriers for this mechanism are quite low (0.03 to 0.30 eV). Vacancy defects are not required, and furthermore have formation energies that are higher than those of the APB kink defects. Furthermore, our results suggest that APB migration is a fundamental diffusion mechanism in the P3 structure.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Dependence of the electronic structure of the EuS/InAs interface on the bonding configuration

The EuS/InAs interface has recently attracted attention for the possibility of inducing magnetic exchange correlations in a strong spin-orbit semiconductor, which could be useful for topological quantum devices. We use density functional theory with a machine-learned Hubbard U correction [Yu et al., npj Comput. Mater. 6, 180 (2020)] to elucidate the effect of the bonding configuration at the interface on the electronic structure. For all interface configurations considered here, we find that the EuS valence band maximum (VBM) lies below the InAs VBM. In addition, dispersed states emerge at the top of the InAs VBM at the interface, which do not exist in either material separately. These states are contributed mainly by the InAs layer adjacent to the interface. They are localized at the interface and may be attributed to charge transfer from the EuS to the InAs. The interface configuration affects the position of the EuS VBM with respect to the InAs VBM, as well as the dispersion of the interface states. For all interface configurations studied here, the induced magnetic moment in the InAs is small. Our results suggest that this interface, in its coherent form studied here, may not be promising for inducing equilibrium magnetic properties in InAs.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Intrinsic anharmonicity and thermal properties of ultralow thermal conductivity Ba 6 Sn 6 Se 13

Ultralow thermal conductivity materials continue to be of great interest for technologically important applications such as thermal insulators and thermoelectrics. Often, such materials possess constituents with extensive dynamic disorder, or “rattlers,” or structural disorder such as metallic glasses. Nevertheless, targeted crystalline bonding order and anharmonicity can provide effective means for the scattering of the phonons in certain materials. Our combined experimental and theoretical investigation of Ba 6 Sn 6 Se 13 reveals such a material, with ultralow thermal conductivity measured over a large temperature range. Optic phonon modes hybridize with acoustic modes at relatively low energies resulting in strong acoustic-optic scattering and limited phonon lifetimes. Moreover, our investigations reveal various avoided crossings at low energies that contribute to phonon scattering. These are derived from chiral phase symmetries and may result in nontrivial topological behaviors for phonon band crossings in this chiral material. This work contributes to the ongoing research on low thermal conductivity materials and the underlying mechanisms that affect their thermal behaviors. Manipulation of these effects may provide pathways for enhancement of material properties for targeted thermal applications using similar materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Optical and magneto-optical properties of ferromagnetic monolayer CrBr 3 : A first-principles G W and G W plus Bethe-Salpeter equation study

The discovery of atomically thin two-dimensional (2D) magnetic semiconductors has triggered enormous research interest recently. In this paper, we use first-principles many-body perturbation theory to study a prototypical 2D ferromagnetic semiconductor, monolayer chromium tribromide (CrBr 3 ). With broken time-reversal symmetry, spin-orbit coupling, and excitonic effects included through the full-spinor GW and GW plus Bethe-Salpeter equation (GW-BSE) methods, we compute the frequency-dependent layer polarizability tensor and dielectric function tensor that govern the optical and magneto-optical (MO) properties. In addition, we provide a detailed theoretical formalism for simulating magnetic circular dichroism, MO Kerr effect, and Faraday effect, demonstrating the approach with monolayer CrBr 3 . Due to reduced dielectric screening in 2D and the localized nature of the Cr d orbitals, we find strong self-energy effects on the quasiparticle band structure of monolayer CrBr 3 that give a 3.8 eV indirect bandgap. Also, excitonic effects dominate the low-energy optical and MO responses in monolayer CrBr 3 where a large exciton binding energy of 2.3 eV is found for the lowest bright exciton state with excitation energy at 1.5 eV. We further find that the MO signals demonstrate strong dependence on the excitation frequency and substrate refractive index. Furthermore, our theoretical framework for modeling optical and MO effects could serve as a powerful theoretical tool for future study of optoelectronic and spintronics devices consisting of van der Waals 2D magnets.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Bending as a control knob for the electronic and optical properties of phosphorene nanoribbons

We have assessed mechanical bending as a powerful controlling tool for the electronic structure and optical properties of phosphorene nanoribbons. We use state-of-the-art density functional approximations in our work. The overall performance of the recently developed meta-generalized gradient approximation (meta-GGA) mTASK [Phys. Rev. Materials 5, 063803 (2021)] functional establishes the method as a useful alternative to the screened hybrid HSE06 for better band gaps of phosphorene nanoribbons. We present a detailed analysis to interpret the optical absorption of bent phosphorene nanoribbons using the GW-Bethe-Salpeter approximation. Finally, we demonstrate the important role of the unoccupied in-gap state and conclude that this in-gap state in armchair nanoribbons introduced by bending can significantly affect the properties of low-energy excitons and add some useful opportunities for applications in optoelectronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2022

First-principles study of magnetic states and the anomalous Hall conductivity of M Nb 3 S 6 ( M = Co, Fe, Mn, and Ni)

Here, inspired by the observation of the extremely large anomalous Hall effect in the absence of applied magnetic fields or uniform magnetization in CoNb 3 S 6 [Nat. Commun. 9, 3280 (2018); Phys. Rev. Research 2, 023051 (2020)], we perform a first-principles study of this and related compounds of the MNb 3 S 6 type with different transition metal M ions to determine their magnetic orders and the anomalous Hall conductivity (AHC). We find that noncoplanar antiferromagnetic ordering is favored relative to collinear or coplanar order in the case of M = Co, Fe, and Ni, while ferromagnetic ordering is favored in MnNb 3 S 6 at low temperatures. The AHC in these materials with noncoplanar spin ordering can reach about e 2 /h per crystalline layer, while being negligible for coplanar and collinear cases. We also find that the AHC depends sensitively on doping and reaches a maximum for intermediate values of the local spin exchange potential between 0.3 and 0.8 eV. Our AHC results are consistent with the reported Hall measurements in CoNb 3 S 6 and suggest a possibility of similarly large anomalous Hall effects in related compounds.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Local positional and spin symmetry breaking as a source of magnetism and insulation in paramagnetic EuTi O 3

Here we consider theoretically the paramagnetic phases of EuTiO 3 that represent configurations created by two sets of microscopic degrees of freedom (m-DOFs): positional symmetry breaking due to octahedral rotations and magnetic symmetry breaking due to spin disorder. The effect of these sets of m-DOFs on the electronic structure and properties of the para phases is assessed by considering sufficiently large (super) cells with the required nominal global average symmetry, allowing, however, the $\mathcal{local}$ positional and magnetic symmetries to be lowered. We find that tendencies for local symmetry breaking can be monitored by following total energy lowering in mean-fieldlike density-functional theory, without recourse for strong correlation effects. While most nominally cubic ABO 3 perovskites are known for their symmetry breaking due to the B-atom sublattice, the case of f-electron magnetism in EuTiO 3 is associated with A-sublattice symmetry breaking and its coupling to structural distortions. We find that (i) paramagnetic $\mathcal{cubic}$ EuTiO 3 has an intrinsic tendency for both magnetic and positional symmetry breaking, while paramagnetic $\mathcal{tetragonal}$ EuTiO 3 has only magnetic symmetry lowering and no noticeable positional symmetry lowering with respect to low-temperature antiferromagnetic tetragonal phase. (ii) Properly modeled paramagnetic tetragonal and cubic EuTiO 3 have a nonzero local magnetic moment on each Eu ion, consistent with the experimental observations of local magnetism in the para phases of EuTiO 3 significantly above the Néel temperature. Interestingly, (iii) the local positional distortion modes in the short-range ordered para phases are inherited from the long-range ordered low-temperature antiferromagnetic ground-state phase.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Mixed excitonic nature in water-oxidized BiVO 4 surfaces with defects

BiVO 4 is a promising photocatalyst for efficient water oxidation, with surface reactivity determined by the structure of active catalytic sites. Surface oxidation in the presence of oxygen vacancies induces electron localization, suggesting an atomistic route to improve the charge transfer efficiency within the catalytic cycle. In this paper, we study the effect of oxygen vacancies on the electronic and optical properties at BiVO 4 surfaces upon water oxidation. We use density functional theory and many-body perturbation theory to explore the change in the electronic and quasiparticle energy levels and to evaluate the electron-hole coupling as a function of the underlying structure. Further we show that while the presence of defects alters the atomic structure and largely modifies the wave-function nature, leading to defect-localized states at the quasiparticle gap region, the optical excitations remain largely unchanged due to the substantial hybridization of defect and nondefect electron-hole transitions. Our findings suggest that defect-induced surface oxidation supports improved electron transport, both through bound and tunable electronic states and via a mixed nature of the optical transitions, expected to reduce electron-hole defect trapping.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Gaussian approximation potential for amorphous Si : H

Hydrogenation of amorphous silicon (a–Si : H) is critical for reducing defect densities, passivating midgap states and surfaces, and improving photoconductivity in silicon-based electro-optical devices. Modeling the atomic-scale structure of this material is critical to understanding these processes, which in turn is needed to describe c–Si/a–Si : H heterojunctions that are at the heart of modern solar cells with world-record efficiency. Density functional theory (DFT) studies achieve the required high accuracy but are limited to moderate system sizes of 100 atoms or so by their high computational cost. Simulations of amorphous materials have been hindered by this high cost because large structural models are required to capture the medium-range order that is characteristic of such materials. Empirical potential models are much faster, but their accuracy is not sufficient to correctly describe the frustrated local structure. Data-driven, machine-learned interatomic potentials have broken this impasse and have been highly successful in describing a variety of amorphous materials in their elemental phase. Here, we extend the Gaussian approximation potential (GAP) for silicon by incorporating the interaction with hydrogen, thereby significantly improving the degree of realism with which amorphous silicon can be modeled. We show that our Si : H GAP enables the simulation of hydrogenated silicon with an accuracy very close to DFT but with computational expense and run times reduced by several orders of magnitude for large structures. Here, we demonstrate the capabilities of the Si : H GAP by creating models of hydrogenated liquid and amorphous silicon and showing that their energies, forces, and stresses are in excellent agreement with DFT results, and their structure as captured by bond and angle distributions are in agreement with both DFT and experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Laplacian-level meta-generalized gradient approximation for solid and liquid metals

In this work, we derive and motivate a Laplacian-level, orbital-free meta-generalized-gradient approximation (LL-MGGA) for the exchange-correlation energy, targeting accurate ground-state properties of sp and sd metallic condensed matter, in which the density functional for the exchange-correlation energy is only weakly nonlocal due to perfect long-range screening. Our model for the orbital-free kinetic energy density restores the fourth-order gradient expansion for exchange to the r 2 SCAN meta-GGA [Furness et al., J. Phys. Chem. Lett. 11, 8208 (2020)], yielding a LL-MGGA we call OFR2. OFR2 matches the accuracy of SCAN for prediction of common lattice constants and improves the equilibrium properties of alkali metals, transition metals, and intermetallics that were degraded relative to the PBE GGA values by both SCAN and r 2 SCAN. We compare OFR2 to the r 2 SCAN -L LL-MGGA [Mejia-Rodriguez and Trickey, Phys. Rev. B 102, 121109(R) (2020)] and show that OFR2 tends to outperform r 2 SCAN -L for the equilibrium properties of solids, but r 2 SCAN -L much better describes the atomization energies of molecules than OFR2 does. For best accuracy in molecules and nonmetallic condensed matter, we continue to recommend SCAN and r 2 SCAN. Numerical performance is discussed in detail, and our paper provides an outlook to machine learning.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Identifying topological superconductivity in two-dimensional transition-metal dichalcogenides

Here we study the superconducting pairing instabilities and gap functions for prototypical two-dimensional (2D) transition-metal dichalcogenides (TMDCs) WS 2 , MoTe 2 , and MoS 2 in the 2H phase under both hole and electron doping at 10 K. Our first-principles quantum many-body Green's function approach allows us to treat the full d and p manifold of orbitals with strong spin-orbit coupling, yielding pairing predictions with material-specific detail. The resulting gap functions exhibit a variety of mixed-parity superconducting states, including s, p, d, f, d ± id, and p ± ip pairing modes. In particular, we predict 3% and 4% hole-doped WS 2 to be a chiral p ± ip topological superconductor. For 1% hole-doped MoS 2 , we find a competition between three doubly degenerate chiral and nonchiral instabilities. Overall, the relative pairing strengths are found to follow the Fermi-surface topology, due to nesting between the Fermi-surface sheets. Finally, we discuss our predictions in relation to available experimental data and classify the topology of the predicted superconducting pairing symmetries.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

First-principles comparative study of Cr migration in O3 and O/P hybrid-phased NaCrO 2

In layered Na transition-metal (TM) oxides, TM migration usually occurs at highly charged states and severely deteriorates the capacity and reversibility. Meanwhile, the formation of hybrid phases with the intergrowth of octahedral (O-type) and prismatic (P-type) Na layers also takes place at highly charged states. These hybrid phases are often more stable than simple O3 or P3 stackings. However, little is known about the mechanism and impact of TM migration in these hybrid phases. In this work, a comparative first-principles study is performed to understand the connections between structural changes and Cr migration in layered O3 and hybrid-phased NaCrO 2 . After Cr migration, the hybrid-phased NaCrO 2 suffers from greater layer shrinkage than the O3 phase. Three factors affect the Cr migration energy E mig : the Na concentration, local 3D configurations, and 2D in-plane geometries. Low Na concentrations and certain local 3D configurations facilitate the Cr migration. The Cr migration barriers in both O3 and hybrid-phased NaCrO 2 are positively correlated with the Cr E mig . The Cr migration in 17 doped O3 and hybrid-phased NaCrO 2 is surveyed. In these doped NaCrO 2 , a more uniform distribution of the Cr–O bond lengths usually suggests suppressed Cr migration. Optimal dopants for suppressing Cr migration are identified by considering E mig for both Cr and the dopant. Finally, our comparative study on Cr migration in O3 and hybrid-phased NaCrO 2 reveals the significant role of hybrid-phased structures in the development of layered cathode materials.

36 MATERIALS SCIENCE↗