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,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Stochastic and mixed density functional theory within the projector augmented wave formalism for simulation of warm dense matter

Stochastic density functional theory (DFT) and mixed stochastic-deterministic DFT are burgeoning approaches for the calculation of the equation of state and transport properties in materials under extreme conditions. In the intermediate warm dense matter regime, a state between correlated condensed matter and kinetic plasma, electrons can range from being highly localized around nuclei to delocalized over the whole simulation cell. The plane-wave basis pseudopotential approach is thus the typical tool of choice for modeling such systems at the DFT level. Unfortunately, stochastic DFT methods scale as the square of the maximum plane-wave energy in this basis. To reduce the effect of this scaling and improve the overall description of the electrons within the pseudopotential approximation, we present stochastic and mixed DFT approaches developed and implemented within the projector augmented wave formalism. In conclusion, we compare results between the different DFT approaches for both single-point and molecular dynamics trajectories and present calculations of self-diffusion coefficients of solid density carbon from 1 to 50 eV.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Low-Temperature Competing Magnetic Energy Scales in the Topological Ferrimagnet TbMn 6 Sn 6

TbMn 6 Sn 6 is a metallic ferrimagnet displaying signatures of both topological electrons and topological magnons arising from ferromagnetism and spin-orbit coupling within its Mn kagome layers. Inelastic neutron scattering measurements find strong ferromagnetic (FM) interactions within the Mn kagome layer and reveal a magnetic bandwidth of ~230 meV. The low-energy magnetic excitations are characterized by strong FM Mn-Mn and antiferromagnetic (AFM) Mn-Tb interlayer magnetic couplings. We observe weaker, competing long-range FM and AFM Mn-Mn interlayer interactions similar to those driving helical magnetism in the YMn 6 Sn 6 system. Combined with density-functional theory calculations, we find that competing Mn-Mn interlayer magnetic interactions occur in all RMn 6 Sn 6 compounds with R=Y, Gd-Lu, resulting in magnetic instabilities and tunability when Mn-R interactions are weak. In the case of TbMn 6 Sn 6 , strong AFM Mn-Tb coupling ensures a highly stable three-dimensional ferrimagnetic network.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Modeling carrier mobility in graphene as a sensitive probe of molecular magnets

Carrier mobility in graphene on a GaAs substrate and its change due to the adsorption of molecular magnets, Mn 12 and [Mn 3 ] 2 , on the surface of graphene, is calculated from first principles. Phonon limited mobility is also calculated for comparison. For Mn 12 adsorption on graphene, the mobility is compared for different organic ligands of Mn 12 (-H, -CH 3 , and - CHCl 2 ), while for [Mn 3 ] 2 dimers, the mobility is calculated for different linkers that yield the ferromagnetic and anti-ferromagnetic configurations of [Mn 3 ] 2 , as well as for different orientations of the molecule. Significant changes in carrier mobility due to the adsorption of the molecules and due to differences in linkers suggest mobility measurement as a possible sensitive probe of magnetic molecules.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Electronic structure of Humble defects in Ge and Ge 0.8⁢ Si 0.2

The group-IV diamond-structure elements are known to host a variety of planar defects, including {001} planar defects in C and {001}, {111}, and {113} planar defects in Si and Ge. Among the {001} planar defects, the Humble defect, known for some time to occur in Ge, has recently also been observed in Si-Ge alloys, but the details of its electronic structure remain poorly understood. Here, in this work, we perform first-principles density-functional calculations to study Humble defects in both Ge and Ge 0.8⁢ Si 0.2 . We also measure the Si L 2,3 -edge electron energy-loss spectra both at the defect and in a bulk-like region far from the defect and compare with theoretical calculations on corresponding Si sites in our first-principles calculations. We find that inclusion of core-hole effects in the theory is essential for reproducing the observed L 2,3 edge spectra, and that once they are included, the results provide a set of fingerprints for different types of local atomic bonding environments in Ge 0.8⁢ Si 0.2 . Our first-principles calculations reveal that the Humble defects have a tendency to enlarge the electronic band gap, which may have potential uses in band engineering. The use of hybrid functionals for an improved description of the band gap in these systems is also discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2023

Galvanic Protected Anode Materials for Molten Regolith Electrolysis

This effort aimed to explore the use of refractory alloys as long-lifetime electrode materials for use in a molten regolith electrolysis (MRE) reactor, a technology that can produce oxygen (O2) from regolith. To find a material capable of withstanding the high temperatures and corrosive nature of molten regolith, first-principles thermodynamic calculations were performed on platinum group metals in the presence of iron (Fe) and silicon (Si) to determine favorable compositions for resisting unwanted secondary phases that could form on the surface of the material. As a result of these first-principles calculations, the kinetic properties of elements such as aluminum (Al), Fe, and Si can be more easily derived, which would provide insight into the longevity of such a material when used in an MRE reactor. The effort funded by this CIF project mapped the entire phase diagram of platinum group metals, which is beneficial not only to MRE where high temperatures and molten material are present, but also aerospace applications like turbine materials, terrestrial metal processing techniques, and energy storage applications. This project narrowed a broad spectrum of refractory metals down to a small list of alloys that would be able to survive not only the high temperatures of an MRE reactor, but also the thermodynamic forces that would cause inferior alloys to form secondary phases, which would ultimately consume the electrodes of the reactor. This research needs to be followed up with experimental demonstrations of these alloys in an MRE environment.

Kevin Grossman↗
DOE OSTI2021

Efficient First-Principles Methodology for the Calculation of the All-Phonon Inelastic Scattering in Solids

Inelastic scattering experiments are key methods for mapping the full dispersion of fundamental excitations of solids in the ground as well as nonequilibrium states. A quantitative analysis of inelastic scattering in terms of phonon excitations requires identifying the role of multiphonon processes. Here, we develop an efficient first-principles methodology for calculating the all-phonon quantum mechanical structure factor of solids. We demonstrate our method by obtaining excellent agreement between measurements and calculations of the diffuse scattering patterns of black phosphorus, showing that multiphonon processes play a substantial role. The present approach constitutes a step towards the interpretation of static and time-resolved electron, x-ray, and neutron inelastic scattering data.

36 MATERIALS SCIENCE↗
DOE OSTI2022

First-Principles Phonon Quasiparticle Theory Applied to a Strongly Anharmonic Halide Perovskite

Understanding and predicting lattice dynamics in strongly anharmonic crystals is one of the long-standing challenges in condensed matter physics. Here we propose a first-principles method that gives accurate quasiparticle (QP) peaks of the phonon spectrum with strong anharmonic broadening. On top of the conventional first-order self-consistent phonon (SC1) dynamical matrix, the proposed method incorporates frequency renormalization effects by the bubble self-energy within the QP approximation. We apply the developed methodology to the strongly anharmonic α-CsPbBr 3 that displays phonon instability within the harmonic approximation in the whole Brillouin zone. While the SC1 theory significantly underestimates the cubic-to-tetragonal phase transition temperature (T c ) by more than 50%, we show that our approach yields T c = 404–423 K, in excellent agreement with the experimental value of 403 K. Furthermore, we also demonstrate that an accurate determination of QP peaks is paramount for quantitative prediction and elucidation of phonon linewidth.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Real-space Green's function approach for intrinsic losses in x-ray spectra

Intrinsic inelastic losses in x-ray spectra originate from excitations in an interacting electron system due to a suddenly created core-hole. These losses characterize the features observed in x-ray photoemission spectra (XPS), as well as many-body effects such as satellites and edge-singularities in x-ray absorption spectra (XAS). However, they are usually neglected in practical calculations. As shown by Langreth these losses can be treated within linear response in terms of a cumulant Green's function in momentum space. Here we present a complementary ab initio real-space Green's function generalization of the Langreth cumulant in terms of the dynamically screened core-hole interaction W c (ω) and the independent particle response function. Here we find that the cumulant kernel β⁡(ω) is analogous to XAS, but with the transition operator replaced by the core-hole potential with monopole selection rules. The behavior reflects the analytic structure of the loss function, with peaks near the zeros of the dielectric function, consistent with delocalized quasiboson excitations. The approach simplifies when W c (ω) is localized and spherically symmetric. In conclusion, illustrative results and comparisons are presented for the electron gas, sodium, and some early transition metal compounds.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Benchmarks and results of the two-band Hubbard model from the Gutzwiller conjugate gradient minimization theory

Ground-state properties, such as energies and double occupancies, of a one-dimensional two-band Hubbard model are calculated using a first-principles Gutzwiller conjugate gradient minimization theory. The favorable agreement with the results from the density matrix renormalization group theory demonstrates the accuracy of our method. A rotationally invariant approach is further incorporated into the method to greatly reduce the computational complexity with a speedup of approximately 50 times. Moreover, we investigate the Mott transition between a metal and a Mott insulator by evaluating the charge gap. In conclusion, with greatly reduced computational effort, our method reproduces the phase diagram in reasonable agreement with the density matrix renormalization group theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Structural and mechanical properties of monolayer amorphous carbon and boron nitride

Amorphous materials exhibit various characteristics that are not featured by crystals and can sometimes be tuned by their degree of disorder (DOD). Here, we report results on the mechanical properties of monolayer amorphous carbon (MAC) and monolayer amorphous boron nitride (maBN) with different DOD. The pertinent structures are obtained by kinetic-Monte-Carlo (kMC) simulations using machine-learning potentials (MLP) with density-functional-theory (DFT)-level accuracy. An intuitive order parameter, namely the areal fraction F x occupied by crystallites within the continuous random network, is proposed to describe the DOD. We find that F x captures the essence of the DOD: Samples with the same F x but different sizes and arrangements of crystallites, obtained using two distinct kMC procedures, have virtually identical radial distributions functions as well as bond-length and bond-angle distributions. Furthermore, by simulating the fracture process with molecular dynamics, we found that the mechanical responses of MAC and maBN before fracture are mainly determined by F x and are insensitive to the sizes and specific arrangements and to some extent the numbers and area distributions of the crystallites. The behavior of cracks in the two materials is analyzed and found to mainly propagate in meandering paths in the CRN region and to be influenced by crystallites in distinct ways that toughen the material. Furthermore, the present results reveal the relation between structure and mechanical properties in amorphous monolayers and may provide a universal toughening strategy for 2D materials.

2-dimensional systems↗
DOE OSTI2024

Electronic structure, self-doping, and superconducting instability in the alternating single-layer trilayer stacking nickelates La 3 ⁢Ni 2 ⁢O 7

Motivated by the recently proposed alternating single-layer trilayer stacking structure for the nickelate La 3 ⁢Ni 2 ⁢O 7 , we comprehensively study this system using ab initio and random-phase approximation techniques. Here, our analysis unveils similarities between this novel La 3⁢ Ni 2 ⁢O 7 structure and other Ruddlesden-Popper nickelate superconductors, such as a similar charge-transfer gap value and orbital-selective behavior of the eg orbitals. Pressure primarily increases the bandwidths of the Ni e g bands, suggesting an enhancement of the itinerant properties of those e g states. By changing the cell volume ratio V/V 0 from 0.9 to 1.10, we found that the bilayer structure in La 3 ⁢Ni 2 ⁢O 7 always has lower energy than the single-layer trilayer stacking La 3 ⁢Ni 2 ⁢O 7 . In addition, we observe a “self-doping” effect (compared to the average 1.5 electrons per eg orbital per site of the entire structure) from the trilayer to the single-layer sublattices and this effect will be enhanced by overall electron doping. Moreover, we find a leading d x 2 -y 2 -wave pairing state that is restricted to the single layer. Because the effective coupling between the single layers is very weak, due to the nonsuperconducting trilayer in-between, this suggests that the superconducting transition temperature T c in this structure should be much lower than in the bilayer structure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2011

Simulation of Relativistic Shocks and Associated Radiation from Turbulent Magnetic Fields

Using our new 3-D relativistic particle-in-cell (PIC) code, we investigated long-term particle acceleration associated with a relativistic electron-positron jet propagating in an unmagnetized ambient electron-positron plasma. The simulations were performed using a much longer simulation system than our previous simulations in order to investigate the full nonlinear stage of the Weibel instability and its particle acceleration mechanism. Cold jet electrons are thermalized and ambient electrons are accelerated in the resulting shocks. Acceleration of ambient electrons leads to a maximum ambient electron density three times larger than the original value as predicted by hydrodynamic compression. Behind the bow shock, in the jet shock, strong electromagnetic fields are generated. These fields may lead to time dependent afterglow emission. In order to go beyond the standard synchrotron model used in astrophysical objects we have used PIC simulations and calculated radiation based on first principles. We calculated radiation from electrons propagating in a uniform parallel magnetic field to verify the technique. We also used the technique to calculate emission from electrons based on simulations with a small system. We obtain spectra which are consistent with those generated from electrons propagating in turbulent magnetic fields. This turbulent magnetic field is similar to the magnetic field generated at an early nonlinear stage of the Weibel instability. A fully developed shock within a larger system may generate a jitter/synchrotron spectrum.

Nishikawa, K.-I.↗
DOE OSTI2021

Gauge invariance of light-matter interactions in first-principle tight-binding models

In this work, we study the different ways of introducing light-matter interaction in first-principle tight-binding (TB) models. The standard way of describing optical properties is the velocity gauge, defined by linear coupling to the vector potential. In finite systems a transformation to represent the electromagnetic radiation by the electric field instead is possible, albeit subtleties arise in periodic systems. The resulting dipole gauge is a multi-orbital generalization of the Peierls substitution. In this work we investigate the accuracy for both pathways, with particular emphasis on gauge invariance, for TB models constructed from maximally localized Wannier functions. This approach accurately captures the light-matter interaction close to the Fermi level. Focusing on paradigmatic two-dimensional materials, we construct first-principle models and calculate the response to electromagnetic fields in linear response and for strong excitations. Benchmarks against fully converged first-principle calculations allow for ascertaining the accuracy of the TB models. We find that the dipole gauge provides a more accurate description than the velocity gauge in all cases. The main deficiency of the velocity gauge is an imperfect cancellation of paramagnetic and diamagnetic current. Formulating a corresponding sum rule however provides a way to explicitly enforce this cancellation. This procedure corrects the TB models in the velocity gauge, yielding excellent agreement with dipole gauge and thus restoring gauge invariance.

2-dimensional systems↗
DOE OSTI2025

Engineering 2D Square Lattice Hubbard Models in 90° Twisted GeX/SnX (X =S, Se) Moiré Superlattices

Because of the large-period superlattices emerging in moiré two-dimensional (2D) materials, electronic states in such systems exhibit low energy flat bands that can be used to simulate strongly correlated physics in a highly tunable setup. While many investigations have thus far focused on moiré flat bands and emergent correlated electron physics in triangular, honeycomb, and quasi-one-dimensional lattices, tunable moiré realizations of square lattices subject to strong correlations remain elusive. Here, in this work, we propose a feasible scheme to construct moiré square lattice systems by twisting two or more layers of 2D materials in a rectangular lattice by 90°. We demonstrate the concept with twisted GeX/SnX (X =S, Se) moiré superlattices and calculate their electronic structures from first principles. We show that the lowest conduction flat band in these systems can be described by a square lattice Hubbard model with parameters which can be controlled by varying the choice of host materials, number of layers, and external electric fields. In particular, twisted double bilayer GeSe realizes a square lattice Hubbard model with strong frustration due to the next-nearest-neighbor hopping that could host unconventional superconductivity, in close analogy to the Hubbard model for copper-oxygen planes of cuprate high-temperature superconductors. The presented scheme uses 90° twisted 2D materials with rectangular unit cells as a promising platform for realizing the physical phenomena of square lattice Hubbard models, establishing a new route for studying its rich phase diagram of magnetism, charge order, and unconventional superconductivity in a highly tunable setting.

2-dimensional systems↗
NASA NTRS1969

Self-consistent many-electron theory of electron work functions and surface potential characteristics for selected metals

Electron work functions, surface potentials, and electron number density distributions and electric fields in the surface region of 26 metals were calculated from first principles within the free electron model. Calculation proceeded from an expression of the total energy as a functional of the electron number density, including exchange and correlation energies, as well as a first inhomogeneity term. The self-consistent solution was obtained via a variational procedure. Surface barriers were due principally to many-body effects; dipole barriers were small only for some alkali metals, becoming quite large for the transition metals. Surface energies were inadequately described by this model, which neglects atomistic effects. Reasonable results were obtained for electron work functions and surface potential characteristics, maximum electron densities varying by a factor of over 60.

Smith, J. R.↗