Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “energy density functional”

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,333 records · Page 74

Accurate prediction of oxygen vacancy concentration with disordered A-site cations in high-entropy perovskite oxides

Abstract Entropic stabilized ABO 3 perovskite oxides promise many applications, including the two-step solar thermochemical hydrogen (STCH) production. Using binary and quaternary A-site mixed {A}FeO 3 as a model system, we reveal that as more cation types, especially above four, are mixed on the A-site, the cell lattice becomes more cubic-like but the local Fe–O octahedrons are more distorted. By comparing four different Density Functional Theory-informed statistical models with experiments, we show that the oxygen vacancy formation energies ( $${E}_{V}^{f}$$ E V f ) distribution and the vacancy interactions must be considered to predict the oxygen non-stoichiometry ( δ ) accurately. For STCH applications, the $${E}_{V}^{f}$$ E V f distribution, including both the average and the spread, can be optimized jointly to improve Δ δ (difference of δ between the two-step conditions) in some hydrogen production levels. This model can be used to predict the range of water splitting that can be thermodynamically improved by mixing cations in {A}FeO 3 perovskites.

08 HYDROGEN↗

Three-dimensional nature of anomalous Hall conductivity in YMn 6 Sn 6−x Ga x , x ≈ 0.55

The unique geometry of kagome lattices leads to topological features such as flat bands and Dirac cones. When paired with ferromagnetism and a Fermi level near Dirac points, they offer a platform for realizing topological Chern magnetotransport. This prospect recently drew interest in the ferrimagnetic kagome metal TbMn 6 Sn 6 . However, density functional theory (DFT) calculations indicate that its 2D Chern gap lies well above the Fermi energy, raising questions about its role in anomalous Hall conductivity. Here, we study YMn 6 Sn 5.45 Ga 0.55 , a structurally and electronically similar material, and find that its intrinsic anomalous Hall effect is three-dimensional. This demonstrates that the Hall response in such compounds does not originate from 2D Chern gaps. Additionally, we confirm that the newly proposed empirical scaling relation for extrinsic Hall conductivity is universally governed by spin fluctuations.

Bhandari, Hari [University of Notre Dame, IN (Unit↗

Spatial locality of electronic correlations in LiFeAs

In this work, we address the question of the degree of spatial nonlocality of the self-energy in the iron-based superconductors, a subject which is receiving considerable attention. Using LiFeAs as a prototypical example, we extract the self-energy from angular-resolved photoemission spectroscopy data. We use two distinct electronic structure references: density functional theory in the local density approximation and linearized quasiparticle self-consistent GW (LQSGW). We find that with the LQSGW reference, spatially local dynamical correlations provide a consistent description of the experimental data, and account for some surprising aspects of the data such as the substantial out-of-plane dispersion of the electron Fermi surface having dominant xz/yz character. Hence, correlations effects can be separated into static nonlocal contributions well described by LQSGW and dynamical local contributions. Hall effect and resistivity data are shown to be consistent with this description.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Accurate localization of Kosterlitz-Thouless-type quantum phase transitions for one-dimensional spinless fermions

We investigate the charge-density wave (CDW) transition for one-dimensional spinless fermions at half band filling with nearest-neighbor electron transfer amplitude t and interaction V. The model is equivalent to the anisotropic XXZ Heisenberg model for which the Bethe Ansatz provides an exact solution. For V>V c =2t, the CDW order parameter and the single-particle gap are finite but exponentially small, as is characteristic for a Kosterlitz-Thouless transition. It is notoriously difficult to locate such infinite-order phase transitions in the phase diagram using approximate analytical and numerical approaches. Second-order Hartree-Fock theory is qualitatively applicable for all interaction strengths, and predicts the CDW transition to occur at V$^{(2)}_{c,2}$≈1.5t. Second-order Hartree Fock theory is almost variational because the density of quasiparticle excitations is small. We apply the density-matrix renormalization group (DMRG) for periodic boundary conditions for system sizes up to 514 sites, which permits a reliable extrapolation of all physical quantities to the thermodynamic limit, apart from the critical region. We investigate the ground-state energy, the gap, the order parameter, the momentum distribution, the quasiparticle density, and the density-density correlation function to locate V c from the DMRG data. In conclusion, tracing the breakdown of the Luttinger liquid and the peak in the quasiparticle density at the band edge permits us to reproduce V c with an accuracy of one percent.

1-dimensional spin chains↗

Critical role of magnetic moments in heavy-fermion materials: Revisiting mysteries of SmB 6

Heavy-fermion family exhibits fascinating and often puzzling properties due to the presence of open-shell f ions and the complexity of the associated charge, orbital, and spin degrees of freedom. SmB 6 is a prototypical heavy-fermion compound that is electrically insulating but yet it displays quantum oscillations, which are a telltale signature of the metallic state. Adding to the enigma is the possibility that SmB 6 is a topological Kondo insulator. Here, by treating the spin degree of freedom on an equal footing with other degrees of freedom using the parameter-free strongly-constrained and appropriately-normed (SCAN) density functional, we explore the ground-state electronic structure of SmB6. A number of competing magnetic phases lying very closely in energy are found, indicating the key role of spin fluctuations in the material. The computed band structure, crystal-field splittings in the f-electron complex, the heavy effective electron mass at the Fermi energy, and the large specific heat are all in good agreement with the corresponding experimental results. In particular, our predicted FS explains the experimentally observed bulk quantum oscillations as well as the low electrical conductivity of SmB 6 . The topological Kondo state of SmB6 is shown to be robust regardless of its magnetic configuration. The excellent performance of SCAN in heavy-fermion systems is explained in terms of its ability to treat self-interaction errors and symmetry breaking within the framework of the density functional theory. Our study provides a new approach for modeling heavy-fermion materials

36 MATERIALS SCIENCE↗

Density of States and Spectral Function of a Superconductor out of a Quantum-Critical Metal

We analyze the validity of a quasiparticle description of a superconducting state above a metallic quantum-critical point (QCP). A normal state at a QCP is a non-Fermi liquid with no coherent quasiparticles. A superconducting order gaps out low-energy excitations, except for a sliver of states for non-s-wave gap symmetry, and at a first glance, restores coherent quasiparticle behavior. We argue that this does not necessarily hold as the fermionic self-energy may remain singular above the gap edge. This singularity gives rise to markedly non-BCS behavior of the density of states and to the appearance of a nondispersing mode at the gap edge in the spectral function. We analyze the set of quantum-critical models with an effective dynamical four-fermion interaction V(Ω) ∝ 1/Ω γ , where Ω is a frequency of a boson, which mediates the interaction. We show that coherent quasiparticle behavior in a superconducting state holds for γ < 1/2, but breaks down for larger γ. Here, we discuss signatures of quasiparticle breakdown and compare our results with the photoemission data for Bi2201 and Bi2212.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Structures and Free Energies of Cerium Ions in Acidic Electrolytes

The Ce 3+ /Ce 4+ redox potential changes with the electrolyte, which could be due to unequal anion complexation free energies between Ce 3+ and Ce 4+ or a change in the solvent electrostatic screening. Ce complexation with anions and solvent screening also affect the solubility of Ce and charge transfer kinetics for electrochemical reactions involving waste remediation and energy storage. Here, we report the structures and free energies of cerium complexes n seven acidic electrolytes based on Extended X-ray Absorption Fine Structure, LTV vis, and Density Functional Theory calculations. Ce 3+ coordinates with nine water molecules as [Ce(H 2 O) 9 ] 3+ in all studied electrolytes. However, Ce 4+ complexes with anions in all electrolytes except HClO 4 . Thus, our results suggest that Ce 4+ anion complexation leads to the large shifts in standard redox potential. Long range screening effects are smaller than the anion complexation energies but could be responsible for changes in the Ce solubility with acid.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Tuning Discharge Behavior of Hollandite α-MnO 2 in Hydrated Zinc Ion Battery by Transition Metal Substitution

The tunnel-type hollandite α-MnO 2 is a promising cathode material for rechargeable aqueous zinc-ion batteries (ZIBs) due to its low cost in synthesis and high energy density. However, irreversible structural degradation upon continuous cycling prevents the cathode from being utilized commercially. Herein, density functional theory (DFT) was utilized to conduct a systematic study on tuning the behavior of α-MnO 2 by substituting the Mn ions on the tunnel wall with a transition metal (V or Cr) during the H + -/Zn 2+ -intercalation in hydrated ZIB. Our study revealed that both substituents aid cyclability and capacity retention with Cr outperforming V. In term of discharge voltage, only the Cr-substitution displays clear promotion at the early stage of discharge. The superior performance of substituted Cr 4+ comes from its unique atomic and electronic structures. Upon discharge, it can be reduced to Cr 3+ more readily than Mn 4+ and thereby limits the formation of unstable Mn 3+ or Mn 2+ centers; the formed Cr 3+ is more stable than Mn 3+ and Mn 2+ from the reduction of Mn 4+ ; and Cr 3+ can also greatly stabilize the neighboring Mn ions. Here this study highlights the significant tuning effect of transition metal substitution on the electrochemical and physical performance of α-MnO 2 as a cathode in hydrated ZIBs.

25 ENERGY STORAGE↗

High-Energy, High-Power Sodium-Ion Batteries from a Layered Organic Cathode

Sodium-ion batteries (SIBs) attract significant attention due to their potential as an alternative energy storage solution, yet challenges persist due to the limited energy density of existing cathode materials. In principle, redox-active organic materials can tackle this challenge because of their high theoretical energy densities. However, electrode-level energy densities of organic electrodes are compromised due to their poor electron/ion transport and severe dissolution. Here, we report the use of a low-bandgap, conductive, and highly insoluble layered metal-free cathode material for SIBs. It exhibits a high theoretical capacity of 355 mAh g –1 per formula unit, enabled by a four-electron redox process, and achieves an electrode-level energy density of 606 Wh kg –1 electrode (90 wt % active material) along with excellent cycling stability. It allows for facile two-dimensional Na+ diffusion, which enables a high intrinsic rate capability. Growth of the active cathode material in the presence of as little as 2 wt % carboxyl-functionalized carbon nanotubes improves charge transport and charge transfer kinetics and further enhances the power performance. Altogether, these allow the construction of SIB cells built from an affordable, sustainable organic small molecule, which provide a cathode energy density of 472 Wh kg –1 electrode when charging/discharging in 90 s and a top specific power of 31.6 kW kg –1 electrode .

36 MATERIALS SCIENCE↗

Static Subspace Approximation for Random Phase Approximation Correlation Energies: Applications to Materials for Catalysis and Electrochemistry

Modeling complex materials using high-fidelity, ab initio methods at low cost is a fundamental goal for quantum chemical software packages. The GW approximation and random phase approximation (RPA) provide a unified description of both electronic structure and total energies using the same physics in a many-body perturbative approach that can be more accurate than generalized-gradient density functional theory (DFT) methods. However, GW/RPA implementations have historically been limited to either specific materials classes or application toward small chemical systems. Here, the static subspace approximation allows for reduced cost full-frequency GW/RPA calculations and has previously been benchmarked thoroughly for GW calculations. Here, we describe our approach to including partial occupations of electronic orbitals in full-frequency GW and RPA calculations for the study of electrocatalysts. We benchmarked RPA total energy calculations using the subspace approximation across a diverse test suite of materials for a variety of computational parameters. The benchmarking quantifies the impact of different extrapolation procedures for representing the static polarizability at infinite screened cutoff, and shows that using screened cutoffs above 20-25 Ryd result in diminishing accuracy returns for predicting RPA total energies. Additionally, for moderately sized electrocatalytic models, 2-3 times fewer computational resources are used to compute RPA total energies by representing the static polarizability with 20-30% of the static subspace basis, with an error of approximately 0.01 eV or better in RPA adsorption energy calculations. Finally, we show that for these electrochemical models RPA can shift DFT adsorption energy shifts by up to 0.5 eV and that GW can frequently shift DFT eigenvalues of surface and adsorbate states by approximately 0.5-1 eV.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Optimization of Magnesium‐Doped Lithium Metal Anode for High Performance Lithium Metal Batteries through Modeling and Experiment

Abstract Lithium (Li)‐magnesium (Mg) alloy with limited Mg amount, which can also be called Mg‐doped Li (Li‐Mg), has been considered as a potential alternative anode for high energy density rechargeable Li metal batteries. However, the optimum doping‐content of Mg in Li‐Mg anode and the mechanism of the improved performance are not well understood. Herein, density functional theory (DFT) calculations are used to investigate the effect of Mg amount in Li‐Mg anode. The Li‐Mg with about 5 wt. % Mg (abbreviated as Li‐Mg5) has the lowest absorption energy of Li, thus all the surface area can be “controlled” by Mg atoms, leading to the smooth and continuous deposition of Li on the surface around the Mg center. A localized high concentration electrolyte enables Li‐Mg5 to exhibit the best cycling stability in Li metal batteries with high‐loading cathode and lean electrolyte under 4.4 V high‐voltage, which is approaching the demand of practical application. This electrolyte also helps generate an inorganic‐rich solid electrolyte interphase, which leads to smooth, compact and less corrosion layer on the Li‐Mg5 surface. Both theoretical simulations and experimental results prove that Li‐Mg5 has optimum Mg content and gives best battery cycling performance.

Gao, Peiyuan↗

Grain Boundary‐Derived Cu + /Cu 0 Interfaces in CuO Nanosheets for Low Overpotential Carbon Dioxide Electroreduction to Ethylene

Abstract Electrochemical CO 2 reduction reaction can be used to produce value‐added hydrocarbon fuels and chemicals by coupling with clean electrical energy. However, highly active, selective, and energy‐efficient CO 2 conversion to multicarbon hydrocarbons, such as C 2 H 4 , remains challenging because of the lack of efficient catalysts. Herein, an ultrasonication‐assisted electrodeposition strategy to synthesize CuO nanosheets for low‐overpotential CO 2 electroreduction to C 2 H 4 is reported. A high C 2 H 4 Faradaic efficiency of 62.5% is achieved over the CuO nanosheets at a small potential of −0.52 V versus a reversible hydrogen electrode, corresponding to a record high half‐cell cathodic energy efficiency of 41%. The selectivity toward C 2 H 4 is maintained for over 60 h of continuous operation. The CuO nanosheets are prone to in situ restructuring during CO 2 reduction, forming abundant grain boundaries (GBs). Stable Cu + /Cu 0 interfaces are derived from the low‐coordinated Cu atoms in the reconstructed GB regions and act as highly active sites for CO 2 reduction at low overpotentials. In situ Raman spectroscopic analysis and density functional theory computation reveal that the Cu + /Cu 0 interfaces offer high *CO surface coverage and lower the activation energy barrier for *CO dimerization, which, in synergy, facilitates CO 2 reduction to C 2 H 4 at low overpotentials.

36 MATERIALS SCIENCE↗

Electronic structure, magnetic properties, and exchange splitting of gadolinium intermetallics

Here, we report here the electronic structure calculations of the magnetic moment, magneto-crystalline anisotropy energy (MAE), and density of states (DOS) of Gd-metal, GdAl 2 , GdAl 3 , GdCo 2 , GdCo 3 , and GdCo 5 . From the DOS, the exchange splitting energy is derived to map ferrimagnetic to paramagnetic phase transition (Curie temperature, T C ). The employed theoretical method here is based on advanced density functional theory, specifically the linearized augmented plane wave method (LAPW) within the local spin density approximation (LSDA), including spin–orbit coupling (SOC) and Hubbard model parameter (U). DOS analysis of both Gd-Al and Gd-Co systems reveals a strong correlation between exchange splitting energy and T C . The Gd-Co system exhibits an enormous exchange splitting energy among Co-3d electrons, vastly overwhelming the vanishing Gd 5d exchange splitting. Finally, these findings suggest that exchange splitting of Co-3d electrons plays a key role in the ground state magnetism and magnetic transition temperatures in these materials.

36 MATERIALS SCIENCE↗

Electronic structure of LaNiO 2 and CaCuO 2 from a self-consistent vertex-corrected GW approach

In this work, the electronic structure of one of the nickelates (LaNiO 2 ) and one of the cuprates (CaCuO 2 ) is studied with three self-consistent GW-based methods: scGW, sc(GW+vertex), and quasiparticle self-consistent GW. Low-energy features obtained in our study are in many respects similar to the features reported in previous density functional theory plus dynamical mean-field theory (DFT+DMFT) studies. Consistent with the DFT+DMFT conclusion, we find LaNiO 2 to be more correlated than CaCuO 2 . However, correlation effects included in our study change the DFT Fermi surface near the Γ point differently from that reported in DMFT studies. Features that are a few electronvolts away from the Fermi level are broader in our calculations than in the DFT+DMFT, which reflects the differences between the DFT and the GW methods. Our results are in qualitative agreement with previous G0W0 results, but the self-consistency results in quantitative differences. Generally, correlation effects are found to be sufficiently weak in both materials, which allows one to use totally ab initio diagrammatic approaches such as sc(GW+Vertex) and to avoid the methods with adjustable parameters (DFT+U or DFT+DMFT). However, the possibility of some strong correlations at low energy that cannot be captured by perturbative methods cannot be completely excluded. For instance, differences in the Fermi surface should be resolved, thus experimental studies are necessary.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Resonantly Enhanced Electromigration Forces for Adsorbates on Graphene

We investigate the electromigration forces for weakly bonded adsorbates on graphene by using density-functional based calculations. We find that the nature of electromigration forces on an adsorbate critically depends on the energy level alignment between the adsorbate state and the Fermi level of the graphene. For a resonant adsorbate, whose frontier orbitals lie close to the Fermi level, the electromigration force is dominated by the electron wind force that is strongly enhanced along the electron flow direction, irrespective of the sign of the adsorbate charge. For a nonresonant adsorbate, the electromigration force is essentially the direct force that depends on the adsorbate charge. We also show that the magnitude of electromigration forces can be continuously tunable through electrostatic gating for resonant adsorbates. Our results provide new insight for understanding and controlling how nanoscale objects behave in or on host materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Analytical harmonic vibrational frequencies with VV10-containing density functionals: Theory, efficient implementation, and benchmark assessments

VV10 is a powerful nonlocal density functional for long-range correlation that is used to include dispersion effects in many modern density functionals, such as the meta-generalized gradient approximation (mGGA), B97M-V, the hybrid GGA, ωB97X-V, and the hybrid mGGA, ωB97M-V. While energies and analytical gradients for VV10 are already widely available, this study reports the first derivation and efficient implementation of the analytical second derivatives of the VV10 energy. The additional compute cost of the VV10 contributions to analytical frequencies is shown to be small in all but the smallest basis sets for recommended grid sizes. Here, this study also reports the assessment of VV10-containing functionals for predicting harmonic frequencies using the analytical second derivative code. The contribution of VV10 to simulating harmonic frequencies is shown to be small for small molecules but important for systems where weak interactions are important, such as water clusters. In the latter cases, B97M-V, ωB97M-V, and ωB97X-V perform very well. The convergence of frequencies with respect to the grid size and atomic orbital basis set size is studied, and recommendations are reported. Finally, scaling factors to allow comparison of scaled harmonic frequencies with experimental fundamental frequencies and to predict zero-point vibrational energy are presented for some recently developed functionals (including r2SCAN, B97M-V, ωB97X-V, M06-SX, and ωB97M-V).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

From Design to Device: Challenges and Opportunities in Computational Discovery of p -Type Transparent Conductors

A high-performance p -type transparent conductor (TC) does not yet exist but could lead to advances in a wide range of optoelectronic applications and enable new architectures for, e.g., next-generation photovoltaic (PV) devices. High-throughput computational material screenings have been a promising approach to filter databases and identify new p -type TC candidates and some of these predictions have been experimentally validated. However, most of these predicted candidates do not have experimentally achieved properties on par with n -type TCs used in solar cells and therefore have not yet been used in commercial devices. Thus, there is still a significant divide between transforming predictions into results that are actually achievable in the laboratory and an even greater lag in scaling predicted materials into functional devices. In this perspective, we outline some of the major disconnects in this materials discovery process—from scaling computational predictions into synthesizable crystals and thin films in the laboratory to scaling laboratory-grown films into real-world solar devices—and share insights to inform future strategies for TC discovery and design. Published by the American Physical Society 2024

14 SOLAR ENERGY↗

Development of a Reactive Force Field for Simulating Photoinitiated Acrylate Polymerization

Light-driven and photo-curable polymer based additive manufacturing (AM) has enormous potential due to its excellent resolution and precision. Acrylated radical chain-growth polymerized resins are widely used in photopolymer AM due to their fast kinetics, and often serve as a departure point for developing other resin materials for photopolymer-based AM technologies. For successful control of the photopolymer resins, the molecular basis of the acrylate free-radical polymerization has to be understood in detail. We present an optimized reactive force field (ReaxFF) for molecular dynamics (MD) simulations of acrylate polymer resins that captures radical polymerization thermodynamics and kinetics. The force field is trained against an extensive training set including density functional theory (DFT) calculations of reaction pathways along the radical polymerization from methyl acrylate to methyl butyrate, bond dissociation energies, and structures and partial charges of several molecules and radicals. We also found that it was critical to train the force field against an incorrect, nonphysical reaction pathway observed in simulations that used parameters not optimized for acrylate polymerization. As a result, the parameterization process utilizes a parallelized search algorithm, and the resulting model can describe polymer resin formation, crosslinking density, conversion rate, and residual monomers of the complex acrylate mixtures.

36 MATERIALS SCIENCE↗