Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “first principles density functional theory”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Ab initio polaritonic potential-energy surfaces for excited-state nanophotonics and polaritonic chemistry

Advances in nanophotonics, quantum optics, and low-dimensional materials have enabled precise control of light–matter interactions down to the nanoscale. Combining concepts from each of these fields, there is now an opportunity to create and manipulate photonic matter via strong coupling of molecules to the electromagnetic field. Toward this goal, here we demonstrate a first principles framework to calculate polaritonic excited-state potential-energy surfaces, transition dipole moments, and transition densities for strongly coupled light–matter systems. In particular, we demonstrate the applicability of our methodology by calculating the polaritonic excited-state manifold of a formaldehyde molecule strongly coupled to an optical cavity. This proof-of-concept calculation shows how strong coupling can be exploited to alter photochemical reaction pathways by influencing avoided crossings with tuning of the cavity frequency and coupling strength. Therefore, by introducing an ab initio method to calculate excited-state potential-energy surfaces, our work opens a new avenue for the field of polaritonic chemistry.

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

Topological surface states on the nonpolar (110) and (111) surfaces of SmB 6

In order to clarify the controversial issue of the topological nature in a mixed-valent Kondo system, SmB 6 , we have explored the surface states on the nonpolar (110) surface of SmB 6 , employing both angle-resolved photoemission spectroscopy (ARPES) experiment and ab initio density-functional theory (DFT) band calculations. Based on ARPES spectroscopic fingerprints and the DFT surface band structures, we ascribe the observed spectral weights at $\bar{X}$ and $\bar{Y}$ on the (110) surface Brillouin zone to topological surface states (TSSs) of “topological insulator (TI)” nature and of “topological crystalline insulator (TCI)” nature, respectively. With varying the chemical potential, the double Dirac cones of the TCI nature exhibit a Lifshitz transition of Fermi surfaces with intriguing spin textures. Further, we have also examined the TSSs on the nearly nonpolar (111) surface of SmB 6 in connection with a recently reported ARPES result and proposed a way to probe the Dirac points that are buried in the bulk-projected bands.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Quantum Monte Carlo calculations of electron scattering from 12 C in the short-time approximation

The short-time approximation is a method introduced to evaluate electroweak nuclear response for systems with A ≥ 12, extending the reach of first-principle many-body quantum Monte Carlo calculations. Using realistic two- and three-body nuclear interactions and consistent one- and two-body electromagnetic currents, we calculate longitudinal and transverse response densities and response functions of 12 C. Here, we compare the resulting cross sections with experimental data for electron-nucleus scattering, finding good agreement.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Anion Association Strength as a Unifying Descriptor for the Reversibility of Divalent Metal Deposition in Nonaqueous Electrolytes

Developing next-generation battery chemistries that move beyond traditional Li-ion systems is critical to enabling transformative advances in electrified transportation and grid-level energy storage. Here, we provide the first evidence for common descriptors for improved reversibility of metal plating/stripping in nonaqueous electrolytes for multivalent ion batteries. Focusing first on the specific role of chloride (Cl – ) in promoting electrochemical reversibility in multivalent systems, rotating disk (RDE) and ring-disk electrode (RRDE) investigations were performed utilizing a variety of divalent cations (Mg 2+ , Zn 2+ , and Cu 2+ ) and the bis-(trifluoromethane sulfonyl) imide (TFSI – ) anion. By introducing varying concentrations of Cl – , a cooperative effect is observed between TFSI – and Cl – that yields the more reversible behavior of mixed electrolytes relative to electrolytes containing only TFSI – . This effect is shown to be general for Mg, Zn, and Cu electrodeposition, and mechanistic understanding of the role of Cl – in improving reversibility of TFSI-based electrolytes is obtained through the combination of R(R)DE experimental results and density functional theory (DFT) evaluation of the redox activity and thermodynamic stability of various TFSI- and Cl-based solution complexes of metal ions. The cooperative anion effect is further generalized to other mixed-anion systems, where systematic variations in anion association strength predicted from DFT (i.e., Cl – > OTf – ≈ TFSI – > BF 4 – > PF 6 – ) yield corresponding trends in redox potentials and improvements of ≥200 mV in the reversibility of metal deposition/dissolution. These results identify anion association strength as a common descriptor for the reversibility of divalent metal anodes and suggest a set of general design principles for developing new electrolytes with improved activity and stability.

25 ENERGY STORAGE↗
DOE OSTI2020

Reduction of electron-phonon coupling in warm dense iron

The electron-ion relaxation dynamics in warm dense iron were investigated by time-resolved x-ray absorption near-edge spectroscopy (XANES). A novel analysis, combining ab initio density functional theory (DFT) and two temperature model (TTM) simulations, was developed to calculate the x-ray absorption spectra as a function of delay time. Here we present experimental evidence of changes at the XAS L3 edge of iron that are consistent with the reduction of the electron-phonon coupling under warm dense matter conditions. The experimental results are in agreement with the model that takes both the electron ($T_e$) and the ion temperature ($T_i$) dependence of the thermophysical properties into consideration, while models where either constant electron-phonon coupling factor (G) or only $T_e$- dependent G are used do not agree with the observed relaxation dynamics of iron.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Epitaxial growth, magnetoresistance, and electronic band structure of GdSb magnetic semimetal films

Motivated by observations of extreme magnetoresistance (XMR) in bulk crystals of rare-earth monopnictide (RE-V) compounds and emerging applications in novel spintronic and plasmonic devices based on thin-film semimetals, we have investigated the electronic band structure and transport behavior of epitaxial GdSb thin films grown on III-V semiconductor surfaces. The Gd 3 + ion in GdSb has a high spin S=7/2 and no orbital angular momentum, serving as a model system for studying the effects of antiferromagnetic order and strong exchange coupling on the resulting Fermi surface and magnetotransport properties of RE-Vs. Here, we present a surface and structural characterization study mapping the optimal synthesis window of thin epitaxial GdSb films grown on III-V lattice-matched buffer layers via molecular-beam epitaxy. To determine the factors limiting XMR in RE-V thin films and provide a benchmark for band-structure predictions of topological phases of RE-Vs, the electronic band structure of GdSb thin films is studied, comparing carrier densities extracted from magnetotransport, angle-resolved photoemission spectroscopy (ARPES), and density-functional theory (DFT) calculations. ARPES shows a hole-carrier rich, topologically trivial, semimetallic band structure close to complete electron-hole compensation, with quantum confinement effects in the thin films observed through the presence of quantum-well states. DFT-predicted Fermi wave vectors are in excellent agreement with values obtained from quantum oscillations observed in magnetic field-dependent resistivity measurements. An electron-rich Hall coefficient is measured despite the higher hole-carrier density, attributed to the higher electron Hall mobility. The carrier mobilities are limited by surface and interface scattering, resulting in lower magnetoresistance than that measured for bulk crystals.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Quantum Monte Carlo Approaches to Na Intercalation on Bilayer Graphene

We have performed Quantum Monte Carlo (QMC) simulations on Na-intercalated bilayer graphene to study the evolution of electronic and optical properties upon Na intercalation into hard carbon layers. The objective was to model the optimal configuration of Na intercalation into a hard carbon matrix containing graphene regions. Our study showed that Na intercalation can be energetically stabilized at large interlayer distances (over 6 Å) in both AA- and AB-stacked bilayer graphene. In the QMC results, we found a significant band gap opening at the equilibrium interlayer distance of Na-intercalated bilayer graphene, while corresponding density functional theory (DFT) results showed no gap. This difference between DFT and QMC results indicates that the gap opening induced by Na intercalation into a hard carbon is underestimated within the DFT framework. In addition, a zigzag configuration of Na atoms was found to be energetically stable at interlayer distances up to 10 Å, leading us to predict the existence of a local minimum of Na intercalation at large interlayer distance. These computation and modeling results can provide guidance on how to synthesize and optimize hard carbon with bilayer graphene regions that permit a zigzag intercalation configuration that will maximize and stabilize sodium hosting.

Binding energy↗
DOE OSTI2021

Electronic structure and topology across T c in the magnetic Weyl semimetal Co 3 Sn 2 S 2

Co 3 Sn 2 S 2 is a magnetic Weyl semimetal, in which ferromagnetic ordering at 177 K is predicted to stabilize Weyl points. We perform temperature and spatial dependent angle--resolved photoemission spectroscopy measurements through the Curie temperature (T c ), which show large band shifts and renormalization concomitant with the onset of magnetism. We argue that Co 3 Sn 2 S 2 evolves from a Mott ferromagnet below T c to a correlated metallic state above T c . To understand the magnetism, we derive a tight-binding model of Co-3d x 2 -y 2 orbitals on the kagome lattice. At the filling obtained by first-principles calculations, this model reproduces the ferromagnetic ground state, and results in the reduction of Coulomb interactions due to cluster effects. Using a disordered local moment simulation, we show how this reduced Hubbard U leads to a collapse of the bands across the magnetic transition, resulting in a correlated state which carries associated characteristic photoemission signatures that are distinct from those of a simple lifting of exchange splitting. Finally, the behavior of topology across T c is discussed in the context of this description of the magnetism.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Phase Diagram of a Deep Potential Water Model

Using the Deep Potential methodology, we construct a model that reproduces accurately the potential energy surface of the SCAN approximation of density functional theory for water, from low temperature and pressure to about 2400 K and 50 GPa, excluding the vapor stability region. The computational efficiency of the model makes it possible to predict its phase diagram using molecular dynamics. Satisfactory overall agreement with experimental results is obtained. Here, the fluid phases, molecular and ionic, and all the stable ice polymorphs, ordered and disordered, are predicted correctly, with the exception of ice III and XV that are stable in experiments, but metastable in the model. The evolution of the atomic dynamics upon heating, as ice VII transforms first into ice VII" and then into an ionic fluid, reveals that molecular dissociation and breaking of the ice rules coexist with strong covalent fluctuations, explaining why only partial ionization was inferred in experiments.

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

Band structures and Z 2 invariants of two-dimensional transition metal dichalcogenide monolayers from fully relativistic Dirac-Kohn-Sham theory using Gaussian-type orbitals

Two-dimensional (2D) materials exhibit a wide range of remarkable phenomena, many of which owe their existence to the relativistic spin-orbit coupling (SOC) effects. To understand and predict properties of materials containing heavy elements, such as the transition-metal dichalcogenides (TMDs), relativistic effects must be taken into account in first-principles calculations. We present an all-electron method based on the four-component Dirac Hamiltonian and Gaussian-type orbitals (GTOs) that overcomes complications associated with linear dependencies and ill-conditioned matrices that arise when diffuse functions are included in the basis. Until now, there has been no systematic study of the convergence of GTO basis sets for periodic solids either at the nonrelativistic or the relativistic level. Here, in this study, we provide such a study of relativistic band structures of the 2D TMDs in the hexagonal (2H), tetragonal (1T), and distorted tetragonal (1T') structures, along with a discussion of their SOC-driven properties (Rashba splitting and Z 2 topological invariants). We demonstrate the viability of our approach even when large basis sets with multiple basis functions involving various valence orbitals (denoted triple- and quadruple-ζ) are used in the relativistic regime. Our method does not require the use of pseudopotentials and provides access to all electronic states within the same framework. Our study paves the way for direct studies of material properties, such as the parameters in spin Hamiltonians, that depend heavily on the electron density near atomic nuclei where relativistic and SOC effects are the strongest.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Bending effects and optical properties of WSe 2 nanoribbons of topological phase

A WSe 2 monolayer of 1T' phase is a large band gap quantum spin Hall insulator, supporting dissipationless charge and spin transports through the topologically protected edge states. Here, in this work, we explore the nanoribbon forms of 1T' phase WSe 2 by first-principles density functional calculations and the many-body perturbation GW and Bethe-Salpeter equation method. We found that the 1T' WSe 2 nanoribbon can show topological edge states with a ribbon width of ~4-6 nm. Those edge bands show crossing through the Fermi level an odd number of times, with one kind of spin-polarization connecting the valence band continuum and conduction band continuum. The topological features of the edge bands hold even under small and medium bending in the nanoribbon, while large bending induces large band splitting, resulting in a topological switch-off in the edge bands. The semiconducting 1T' WSe2 nanoribbon shows a large tunability with bending in optical absorption spectra and exciton states. The lowest-energy exciton is changed from optically dark in the flat nanoribbon to bright in the bent nanoribbons. These properties in the 1T' WSe 2 nanoribbons suggest potential applications in controllable quantum electronics and exciton-based quantum information processes.

36 MATERIALS SCIENCE↗
DOE OSTI2020

NWChem: Past, present, and future

Specialized computational chemistry packages have permanently reshaped the landscape of chemical sciences by providing tools to support and guide the experimental effort and for prediction of chemical and materials properties. In this regard, a special role has been played by electronic structure packages where complex chemical and materials processes can be modeled using first-principle-driven methodologies. Over the last few decades, the rapid development of computing technologies and a tremendous increase in computational power has offered a unique chance to study complex chemical transformations using sophisticated and predictive many-body techniques to describe correlated behavior of electrons in molecular and condensed phase systems at different levels of theory. In enabling these simulations, a critical role has been played by novel parallel algorithms capable of taking advantage of computational resources to address polynomial scaling of electronic structure methods. NWChem was among the first electronic structure codes that focused on delivering scalable performance for electronic structure simulations. Herein, we briefly review the NWChem suite of computational codes including its history, design principles, parallel tools, current capabilities, outreach and outlook.

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

Atomistic modeling of lanthanide diffusion in refractory body-centered cubic molybdenum

Lanthanide fission products can strongly interact with candidate cladding alloys, but their transport properties in refractory metals remain poorly understood. Here, in this work, we investigate the atomic-scale diffusion behavior of La, Ce, Pr, and Nd in body-centered cubic (bcc) molybdenum, a potential candidate for advanced nuclear cladding. Self-consistent mean-field transport modeling is performed to evaluate the fission product transport and vacancy mobility, informed by first-principles and nudged elastic band calculations of vacancy formation energies, migration barriers, and solute–vacancy binding characteristics. Compared with bcc Fe, lanthanide solutes in bcc Mo exhibit slower tracer diffusion due to higher vacancy formation and migration energies. Furthermore, the calculations reveal that the influence of fission products on migration barriers in bcc Mo are not as extensive in range compared to bcc Fe. Among the studied lanthanides, La exhibits the strongest vacancy binding while also being the fastest diffuser in Mo. These findings highlight how refractory bcc alloys can reduce fission product infiltration, offering valuable insight into the development of durable cladding systems for advanced reactors.

36 - MATERIALS SCIENCE↗
DOE OSTI2022

Calculated electron paramagnetic resonance g-tensor and hyperfine parameters for zinc vacancy and N related defects in ZnO

Various defects in ZnO, focused on substitutional N O and N 2 in various sites, O-site, interstitial and Zn-site are studied using first-principles calculations with the goal of understanding the electron paramagnetic resonance (EPR) center reported for N 2 in ZnO and substitutional N on the O-site. The g tensors are calculated using the gauge including projector augmented wave (GIPAW) method and compared with experiments. The g-tensor of the free $N^{+}_{2}$ and $N^{–}_{2}$ radicals and their various contributions within the GIPAW theory are analyzed first to provide a baseline reference for the accuracy of the method and for understanding the N 2 behavior in ZnO. Previous controversies on the site location of N 2 in ZnO for this EPR center and on the shallow or deep nature and donor or acceptor nature of this center are resolved. Here, we find that the N 2 on the Zn site is mostly zinc-vacancy like in its spin density and g-tensor, while for the O-site, a model with the N 2 axis lying in-the basal plane and the singly occupied π g -orbital along the c axis provides good agreement with experiment. For the interstitial location, if the N 2 is not strongly interacting with the surroundings, no levels in the gap are found and hence also no possible EPR center. The calculated g-tensors for N O and V Zn are also found to be in good agreement with experiment. The effects of different functionals affecting the localization of the spin density are shown to affect the g-tensor values.

36 MATERIALS SCIENCE↗
DOE OSTI2021

First principles study of nearly strain-free Ni/WSe 2 and Ni/MoS 2 interfaces

Metal/transition metal dichalcogenide interfaces are the subject of active research, in part because they provide various possibilities for interplay of electronic and magnetic properties with potential device applications. Here, we present results of our first principles calculations of nearly strain-free Ni/WSe 2 and Ni/MoS 2 interfaces in thin-film geometry. It is shown that while both the WSe 2 and MoS 2 layers adjacent to Ni undergo metallic transition, the layers farther from the interface remain semiconducting. In addition, a moderate value of spin-polarization is induced on interfacial WSe 2 and MoS 2 layers. At the same time, the electronic and magnetic properties of Ni are nearly unaffected by the presence of WSe 2 and MoS 2 , except a small reduction of magnetic moment at the interfacial Ni atoms. Furthermore, these results can be used as a reference for experimental efforts on epitaxial metal/transition metal dichalcogenide heterostructures, with potential application in modern magnetic storage devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Hybridization effect on the x-ray absorption spectra for actinide materials: Application to PuB 4

Studying the local moment and 5 f -electron occupations sheds insight into the electronic behavior in actinide materials. X-ray absorption spectroscopy (XAS) has been a powerful tool to reveal the valence electronic structure when assisted with theoretical calculations. However, the analysis currently taken in the community on the branching ratio of the XAS spectra generally does not account for the hybridization effects between local f orbitals and conduction states. In this paper, we discuss an approach which employs the density functional theory plus Gutzwiller rotationally invariant slave boson method to obtain a local Hamiltonian for the single-impurity Anderson model, and calculates the XAS spectra by the exact diagonalization (ED) method. A customized numerical routine was implemented for the ED XAS part of the calculation. By applying this technique to the recently discovered 5 f -electron topological Kondo insulator Pu B 4 , we determined the signature of 5 f -electronic correlation effects in the theoretical x-ray spectra. Furthermore, we found that the Pu 5 f - 6 d hybridization effect provides an extra channel to mix the j = 5 / 2 and 7 / 2 orbitals in the 5 f valence. As a consequence, the resulting electron occupation number and spin-orbit coupling strength deviate from the intermediate-coupling regime.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Boosting CO 2 reduction on Fe-N-C with sulfur incorporation: Synergistic electronic and structural engineering

Developing earth-abundant efficient catalysts for CO 2 reduction reaction (CO 2 RR) is of paramount importance for electrochemical conversion of CO 2 into value-added products. Despite numerous studies on iron and nitrogen codoped carbon (Fe-N-C) catalysts, grand challenges exist due to limited performance and understanding of catalytic mechanisms. This study reports a general strategy to boost electrocatalytic CO 2 RR activity of Fe-N-C with the incorporation of S atoms to engineer carbon support structure and electronic properties of active Fe-N sites simultaneously via a copolymer-assisted synthetic approach. The employment of N,S comonomers significantly increases the numbers of micropores and surface area, enabling dense atomic Fe-N and enhanced utilization efficiency. We report the first-principles calculations reveal that S modulation upraises the Fermi energy of Fe 3d and increases charge density on Fe atoms of Fe-N 4 , thereby enhancing intrinsic catalytic reactivity and selectivity for CO 2 reduction by strengthening the binding interaction between the Fe site and key COOH* intermediate. These integrated structural and electronic merits endow Fe-NS-C with outstanding activity (e.g., CO Faradaic efficiency of 98% at an overpotential of 490 mV) and stability (without deactivation in 30 h), ranking it one of the most active Fe-N-C reported to date. The finding offers an innovative design strategy to enable the design of advanced catalysts for CO 2 conversion.

25 ENERGY STORAGE↗
DOE OSTI2021

High-pressure phase diagram of beryllium from ab initio free-energy calculations

In this report we use first-principles molecular dynamics simulations coupled with the thermodynamic integration method to study the hexagonal close-packed (hcp) to body-centered cubic (bcc) transition and melting of beryllium up to a pressure of 1600 GPa. We derive the melting line by equating solid and liquid Gibbs free energies and represent it by a Simon-Glatzel fit T m = 1564 K [1 + P/(15.6032 GPa)] 0.383 , which is in good agreement with previous two-phase simulations <6000 K. We also derive the hcp-bcc solid-solid phase boundary and show that the quasiharmonic approximation underestimates the stability of the hcp structure, predicting lower transition pressures between hcp and bcc phases. Our results are consistent with the stability regime predicted by the phonon quasiparticle method. We also predict that the hcp-bcc-liquid triple point is located at 164.7 GPa and 4314 K. In addition, we compute the shock Hugoniot curve and show that it is in good agreement with experiments, intersecting our derived melting curve at ~235 GPa and 4900 K. Finally, we make predictions for future ramp compression experiments. Starting with an isentropic compression of the liquid, we predict the path to intersect the melting line at low pressure and temperature, then to continue along the melting line over a large temperature interval of 7000 K as the sample remains in the mixed solid-liquid state before it enters the solid phase.

36 MATERIALS SCIENCE↗