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 1,045 records · Page 58

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↗

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↗

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↗

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↗

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↗

Comparing time-dependent density functional theory with many-body perturbation theory for semiconductors: Screened range-separated hybrids and the GW plus Bethe-Salpeter approach

We present band structure and optical absorption spectra obtained from density functional theory (DFT) and linear response time-dependent DFT (TDDFT) calculations using a screened range-separated hybrid (SRSH) functional, including spin-orbit coupling, for seven prototypical semiconductors. The results are compared to those obtained from highly converged many-body perturbation theory calculations using the GW approximation and the GW plus Bethe-Salpeter equation (GW-BSE) approaches. We use a single empirical parameter for our SRSH calculations, fit such that the SRSH band gap reproduces the GW band gap at the Γ point. We then find that ground-state generalized Kohn-Sham SRSH eigenvalues accurately reproduce the band structure obtained from GW calculations, typically to within 0.1–0.2 eV, and optical absorption spectra obtained using TDDFT with the SRSH functional agree well with those of GW-BSE, with a mean deviation of 0.03 and 0.11 eV for the location of the first and second absorption peaks, respectively, at a fraction of the computational cost.

36 MATERIALS SCIENCE↗

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↗

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↗

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↗

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↗

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↗

In Situ Electron Microscopy Investigation of Sodiation of Titanium Disulfide Nanoflakes

Two-dimensional (2D) metal sulfides show great promises for their potential applications as electrode materials of sodium ion batteries because of the weak interlayer van der Waals interactions which allow the reversible accommodation and extraction of sodium ions. The sodiation of metal sulfides can undergo a distinct process compared to that of lithiation, which is determined by their metal and structural types. Yet, the structural and morphological evolution during its electrochemical sodiation is still unclear. In this report, we studied the sodiation reaction dynamics of TiS 2 by employing in situ transmission electron microscopy and first-principles calculations. During the sodium ion intercalation process, we observed multiple intermediate phases (phase II, phase Ib, and phase Ia), different from its lithiation counterpart, with varied sodium occupation sites and interlayer stacking sequences. Further insertion of Na ions prompted a multi-step extrusion reaction which led to the phase separation of Ti metal from the Na 2 S matrix, with its 2D morphology expanded to a 3D morphology. In contrast to regular conversion electrodes, TiS 2 still maintained a compacted structure after a full sodiation. Frist-principles calculations reveal that the as-identified phases are thermodynamically preferred at corresponding intercalation/extrusion stages compared to other possible phases. The current work provides the fundamental mechanistic understanding of the sodiation process of 2D transition metal sulfides.

25 ENERGY STORAGE↗

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↗

Prediction of topological phases in metastable ferromagnetic $MPX_3$ monolayers

Density functional theory calculations are carried out to study the electronic and topological properties of MPX 3 (M = Mn, Fe, Co,Ni; X = S, Se) monolayers in the ferromagnetic (FM) metastable magnetic state. Here we find that FM MnPSe 3 monolayers host topological semimetal signatures that are gapped out when spin-orbit coupling is included. These findings are supported by explicit calculations of the Berry curvature and the Chern number. The choice of the Hubbard-U parameter to describe the d-electrons is thoroughly discussed, as well as the influence of using a hybrid-functional approach. The presence of band inversions and the associated topological features are found to be formalism-dependent. Nevertheless, routes to achieve the topological phase via the application of external biaxial strain are demonstrated. Within the hybrid-functional picture, topological band structures are recovered under a pressure of 15% (17 GPa). The present work provides a potential avenue for uncovering new topological phases in metastable ferromagnetic phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

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↗

Starting-point-independent quantum Monte Carlo calculations of iron oxide

Quantum Monte Carlo (QMC) methods are useful for studies of strongly correlated materials because they are many body in nature and use the physical Hamiltonian. Typical calculations assume as a starting point a wave function constructed from single-particle orbitals obtained from one-body methods, e.g., density functional theory. However, mean-field-derived wave functions can sometimes lead to systematic QMC biases if the mean-field result poorly describes the true ground state. In this study, we examine the accuracy and flexibility of QMC trial wave functions using variational and fixed-node diffusion QMC estimates of the total spin density and lattice distortion of antiferromagnetic iron oxide (FeO) in the ground state B1 crystal structure. We found that for relatively simple wave functions the predicted lattice distortion was controlled by the choice of single-particle orbitals used to construct the wave function, rather than by subsequent wave function optimization techniques within QMC. By optimizing the orbitals with QMC, we then demonstrate starting-point independence of the trial wave function with respect to the method by which the orbitals were constructed by demonstrating convergence of the energy, spin density, and predicted lattice distortion for two qualitatively different sets of orbitals. The results suggest that orbital optimization is a promising method for accurate many-body calculations of strongly correlated condensed phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Computational design of high-entropy rare earth aluminum garnets for advanced thermal and environmental barrier coatings

To enhance the protection of Ni-based superalloys in gas turbine engines’ high-temperature environments, it’s crucial to develop advanced thermal/environmental barrier coating (T/EBC) materials with a balanced combination of thermal and mechanical properties. This optimization is essential to safeguard against chemical and thermal challenges. Here, in this study, we harness the power of density functional theory (DFT) in conjunction with combinatorial chemistry methodologies to engineer high-performance high-entropy rare earth disilicates of the RE 3 Al 5 O 12 family (where RE denotes Y, Gd, Er, and Yb). These materials are meticulously designed to exhibit superior phase stability, a targeted coefficient of thermal expansion (CTE), low lattice thermal conductivity, and robust mechanical properties. The determination of CTE values is accomplished through phonon calculations at various volume settings within the quasi-harmonic approximation, while lattice thermal conductivities are rigorously assessed employing the Debye-Callaway model, accounting for three distinct phonon processes. Our findings highlight the remarkable attributes of the solid solution (Y 1/4 Gd 1/4 Er 1/4 Yb 1/4 ) 3 Al 5 O 12 , which displays a reduction in lattice thermal conductivity compared to its individual constituents while maintaining a favorable range of CTE values. The novel T/EBC material, distinguished by their multifaceted functionalities, are poised to use in substantial enhancements in the performance of engines.

36 MATERIALS SCIENCE↗

Impact of classical statistics on thermal conductivity predictions of BAs and diamond using machine learning molecular dynamics

Machine learning interatomic potentials (MLIPs) have greatly enhanced molecular dynamics (MD) simulations, achieving near-first-principles accuracy in thermal conductivity studies. In this work, we reveal that this accuracy, observed in BAs and diamond at sub-Debye temperatures, stems from an accidental error cancelation: classical statistics overestimates specific heat while underestimating phonon lifetimes, balancing out in thermal conductivity predictions. However, this balance is disrupted when isotopes are introduced, leading MLIP-based MD to significantly underpredict thermal conductivity compared to experiments and quantum statistics-based Boltzmann transport equation. This discrepancy arises not from classical statistics affecting phonon–isotope scattering rates but from its impact on the interplay between phonon–isotope and phonon–phonon scattering in the normal scattering-dominated BAs and diamond. In conclusion, this work underscores the limitations of MLIP-based MD for thermal conductivity studies at sub-Debye temperatures.

36 MATERIALS SCIENCE↗