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,171 records · Page 65

Pseudo–Jahn–Teller Distortion in a One-Dimensional π-Conjugated Polymer

Structural distortions in low-dimensional π-conjugated systems profoundly influence their electronic properties, but the control of such behavior in laterally extended systems remains challenging. Here, in this study, we demonstrate that a one-dimensional conjugated polymer─ poly -(difluorenoheptalene-ethynylene) (PDFHE)─undergoes a pronounced out-of-plane backbone distortion, equivalent to a spontaneous symmetry breaking (SSB) of its mirror symmetry. We synthesized PDFHE on noble metal surfaces and characterized its structure and electronic states using low-temperature scanning tunneling microscopy (STM). Rather than adopting a planar, high-symmetry conformation, PDFHE relaxes into nonplanar isomers stabilized by a pseudo–Jahn–Teller (PJT) distortion having mirror-odd out-of-plane character. The distortion lowers the total energy and increases the band gap, providing a concise rationale for the observed symmetry breaking. Density functional theory calculations corroborate these findings, providing a microscopic explanation for the SSB. Our results show that even in mechanically robust extended π-systems, subtle electron–lattice coupling can spontaneously drive significant structural rearrangements, even in mechanically robust extended π-systems.

conformation↗

Role of water model on ion dissociation at ambient conditions

In this work, we study ion pair dissociation in water at ambient conditions using a combination of classical and ab initio approaches. The goal of this study is to disentangle the sources of discrepancy observed in computed potentials of mean force. In particular, we aim to understand why some models favor the stability of solvent-separated ion pairs vs contact ion pairs. We found that some observed differences can be explained by non-converged simulation parameters. However, we also unveil that for some models, small changes in the solution density can have significant effects on modifying the equilibrium balance between the two configurations. We conclude that the thermodynamic stability of contact and solvent-separated ion pairs is very sensitive to the dielectric properties of the underlying simulation model. In general, classical models are very robust in providing a similar estimation of the contact ion pair stability, while this is much more variable in density functional theory-based models. The barrier to transition from the solvent-separated to contact ion pair is fundamentally dependent on the balance between electrostatic potential energy and entropy. This reflects the importance of water intra- and inter-molecular polarizability in obtaining an accurate description of the screened ion–ion interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Analysis of Unpolarised p+p¿ Photoproduction with the GlueX Experiment

This thesis presents measurements of spin-density matrix elements in unpolarised p+p? photoproduction on a proton target. The dominant resonance contribution to the dipion system is the r(770) meson. Due to the large production cross section for this resonance, a valid comparison can be made between the obtained final results and r(770) spin-density matrix elements measured previously with other experiments. The measurement was performed over the 3:0 ? 11:6 GeV beam energy regime, which is a more extensive energy range than has ever been studied previously for the r(770). Results were obtained by analysing data from the GlueX experiment based at Jefferson Lab. Extended maximum likelihood fits were applied to extract three spin-density matrix elements using Markov chain Monte Carlo based parameter estimations. This was performed using various binning configurations to probe the energy, mass, and four-momentum transfer dependence of the determined physics observables. Spin-density matrix elements are shown to be consistent with the model of s-channel helicity conservation at low ?t. The effects of pomeron and f2 exchanges are clearly visible in the energy dependence of the measured observables. These observations provide valuable insights into the relative strengths of both processes as a function of the photon energy, and may enable theorists to disentangle the f2=P coupling ratio. Spin-density matrix elements are seen to be highly dependent on the reconstructed resonance mass. This observation is likely to be a result of non-resonant S-wave background processes, and emphasises the need for a more detailed model of the p+p? angular distribution that considers all of the competing angular momentum components that contribute to the measured final state. The statistical precision of measurements performed for this thesis surpass what was achievable in previous studies of the r(770) by several orders of magnitude. Studies of the energy and four-momentum transfer dependence, and insights into the effects of non r(770) background contributions provide valuable input for production models. This will help inform the choice of wave-sets used for partial wave analyses, supporting GlueX in its search for exotic hybrid meson states.

Fitches, James↗

Using DAPPER to extract the photon strength function of 58 Fe using the inverse Oslo and shape methods

The photon strength function of 58 Fe has been extracted using both the Oslo and Shape methods from particle–γ coincidence data measured using the Detector Array for Photons, Protons, and Exotic Residues, which probes nuclei utilizing (d,p) reactions in inverse kinematics. Four particle–γ coincidence matrices, each constructed with different treatments of the γ–ray energies, are explored in order to observe the impact on the resulting nuclear level density and photon strength. The final photon strength function reported is found to agree well with previous Oslo measurements of other iron isotopes. Systematic uncertainties are included, using different model parameters and their reported errors to perform the Oslo method normalization. The model-independent Shape method is explored and the functional form of the photon strength function obtained is in agreement with the Oslo method results. A low-energy enhancement is not reported for 58 Fe in this work given possible subtraction issues originating from strongly populated states.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Energy conservation in real-time nuclear–electronic orbital Ehrenfest dynamics

Real-time nuclear–electronic orbital Ehrenfest (RT-NEO-Ehrenfest) dynamics methods provide a first-principles approach for describing nonadiabatic molecular processes with nuclear quantum effects. For an efficient description of proton transfer within RT-NEO-Ehrenfest dynamics, the basis function center associated with the quantum proton can be allowed to move classically. Here, this traveling proton basis (TPB) approach effectively captures proton quantum dynamics, although its energy conservation behavior is not yet fully satisfactory. Two recently proposed TPB approaches, in principle, conserve the extended energy, which includes both the system energy and the kinetic energy associated with the proton basis function center. Herein, a thermostatted TPB approach is proposed to improve the conservation of the system energy, excluding the kinetic energy associated with the proton basis function center. In this approach, the quantum proton dynamics are modulated by dynamically rescaling the proton momentum operator to maintain the system energy conservation. With the excited-state intramolecular proton transfer of o-hydroxybenzaldehyde as an example, this approach is shown to significantly improve the system energy conservation while preserving the accuracy of the quantum proton dynamics as achieved in the original TPB approach.

Ehrenfest dynamics↗

Bridging the gap between molecules and materials in quantum chemistry with localized active spaces

The number of materials that “bridge the gap” between single molecules and extended solids, such as metal-organic frameworks and organic semiconductors, has been increasing. Consequently, there is a growing need for modeling approaches that effectively integrate the real-space molecular perspective employed by computational chemists and the reciprocal-space dispersive perspective employed by computational physicists. Here, we propose the localized active space (LAS) approach as a promising method to successfully bridge this gap. The LAS approach extends the active space concept from multiconfigurational methods such as complete active space self-consistent field theory to multiple molecular fragments via a product-form wave function ansatz. Here, we apply this method to solid state phenomena by treating each unit cell as a fragment with different sets of local quantum numbers (e.g., charge and excitation number). State interaction between these LAS states (LASSI) thus provides a comprehensive basis for the study of charge and energy transfer, meeting and surpassing the capabilities of single-reference fragmentation approaches such as constrained density functional theory (cDFT). Most centrally, we show how combining this LASSI approach with multiconfigurational pair-density functional theory (MC-PDFT) provides an elegant and efficient method to compute band structures that capture multiconfigurational character. We apply the LASSI band structure approach to the computation of band gaps in stretched hydrogen chain, polyacetylene, and bulk nickel oxide (NiO), finding good or excellent quantitative agreement with reference values in all cases. Additionally, we use the LAS basis in one-dimensional model systems to demonstrate its ability to treat difficult solid-state phenomena such as exciton transfer and excitation at p-n junctions.

method development↗

Density Functional Theory Solute-Dislocation Interactions

DFT-optimized structures of dislocations with a solute at different sites. The dataset include the VASP input files (except POTCAR) and output files OSZICAR and OUTCAR. There are also text files containing the DFT-computed interaction energies. The README file contains expressions that can be used to compute the interaction energies and chemical misfit energies using the volumetric size misfit for any solute species in Mg.

36 MATERIALS SCIENCE↗

Nonequivalent Atomic Vibrations at Interfaces in a Polar Superlattice

In heterostructures made from polar materials, e.g., AlN–GaN–AlN, the nonequivalence of the two interfaces is long recognized as a critical aspect of their electronic properties; in that, they host different 2D carrier gases. Interfaces play an important role in the vibrational properties of materials, where interface states enhance thermal conductivity and can generate unique infrared-optical activity. The nonequivalence of the corresponding interface atomic vibrations, however, is not investigated so far due to a lack of experimental techniques with both high spatial and high spectral resolution. Herein, the nonequivalence of AlN–(Al 0.65 Ga 0.35 )N and (Al 0.65 Ga 0.35 )N–AlN interface vibrations is experimentally demonstrated using monochromated electron energy-loss spectroscopy in the scanning transmission electron microscope (STEM-EELS) and density-functional-theory (DFT) calculations are employed to gain insights in the physical origins of observations. It is demonstrated that STEM-EELS possesses sensitivity to the displacement vector of the vibrational modes as well as the frequency, which is as critical to understanding vibrations as polarization in optical spectroscopies. The combination enables direct mapping of the nonequivalent interface phonons between materials with different phonon polarizations. Furthermore, the results demonstrate the capacity to carefully assess the vibrational properties of complex heterostructures where interface states dominate the functional properties.

36 MATERIALS SCIENCE↗

Ab Initio Molecular Dynamics Study of the Reduction of Acetone by the Hydrated Electron

We have investigated the reaction dynamics of the reduction of acetone by the solvated electron in water using ab initio molecular dynamics simulations at 298 and 373 K. The rate constants derived from the simulations are consistent with experimental observations that the reaction has a low activation energy. Detailed analyses carried out to shed light on the mechanism of the reaction show that solvent reorganization plays a key role as a reaction coordinate, as expected from the Marcus electron transfer theory. Furthermore, constrained density functional theory calculations indicate that the electronic coupling is large, placing the reaction in the adiabatic limit. Indeed, the activation energies and rate constants of the simulations are in accordance with the predictions of adiabatic Marcus theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Equation of state of hot, dense magnesium derived with first-principles computer simulations

Using two first-principles computer simulation techniques, path integral Monte Carlo and density functional theory molecular dynamics, we derive the equation of state of magnesium in the regime of warm dense matter, with densities ranging from 0.43 to 86.11 g cm –3 and temperatures from 20,000 K to 5×10 8 K. These conditions are relevant for the interiors of giant planets and stars as well as for shock compression measurements and inertial confinement fusion experiments. Here, we study ionization mechanisms and the electronic structure of magnesium as a function of density and temperature. We show that the L shell electrons, 2s and 2p energy bands, merge at high densities. This results in gradual ionization of the L-shell with increasing density and temperature. In this regard, Mg differs from MgO, which is also reflected in the shape of its principal shock Hugoniot curve. For Mg, we predict a single broad pressure-temperature region, where the shock compression ratio is approximately 4.9. Mg thus differs from Si and Al plasmas that exhibit two well-separated compression maxima on the Hugoniot curve for L and K shell ionizations. Finally, we study multiple shocks and effects of preheat and precompression.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Upward band gap bowing and negative mixing enthalpy in multi-component cubic halide perovskite alloys

Physical properties intermediate between constituents of alloys can be achieved as downward convex positive bowing, upward concave negative bowing, or zero bowing. Such bowing effects are essential for band gap engineering in semiconductor alloys. Upward band gap bowing effects are rather rare, hindering the exploration on half of the available physical property space of alloys. Part of this being a rare event is related to the need to stabilize an alloy with low mixing enthalpy, so it does not phase separate. Here, in this paper, we find via density functional theory that one can satisfy the simultaneous conditions of negative mixing enthalpy and upward band gap bowing in four-component AB⁢X 3 halide perovskite alloys in the cubic perovskite structure. Such perovskite alloys have the B-site occupied by a mixture of group IVB and IIB elements that have the IVB-𝑠 and IIB-s states in the valence bands and conduction bands, respectively, leaving the delocalized 𝑠 states to strongly repel each other. This 𝑠−𝑠 repulsion leads to the upward band gap bowing and negative mixing enthalpies simultaneously. Remarkably, we identify a perovskite alloy that has a band gap much larger than all its components. Analogous trends of upward band gap bowing and negative mixing enthalpy also appear in the corresponding three-component and two-component AB⁢X 3 halide perovskite alloys. These observations of upward band gap bowing and negative mixing enthalpy will significantly accelerate the design of stable upward band gap bowing alloys in a broad range of material families.

14 SOLAR ENERGY↗

Kohn-Sham density functional theory of Abelian anyons

In this work, we develop a density functional treatment of noninteracting Abelian anyons, which is capable, in principle, of dealing with a system of a large number of anyons in an external potential. Comparison with exact results for a few particles shows that the model captures the behavior qualitatively and semiquantitatively, especially in the vicinity of the fermionic statistics. We then study anyons with statistics parameter 1+1/n, which are thought to condense into a superconducting state. An indication of the superconducting behavior is the mean-field result that, for uniform density systems, the ground-state energy increases under the application of an external magnetic field independent of its direction. Our density functional theory based analysis does not find that to be the case for finite systems of anyons, which can accommodate a weak external magnetic field through density transfer between the bulk and the boundary rather than through transitions across effective Landau levels, but the “Meissner repulsion” of the external magnetic field is recovered in the thermodynamic limit as the effect of the boundary becomes negligible. We also consider the quantum Hall effect of anyons, and show that its topological properties, such as the charge and statistics of the excitations and the quantized Hall conductance, arise in a self-consistent fashion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Quantum solver for single-impurity Anderson models with particle-hole symmetry

Quantum embedding methods, such as dynamical mean-field theory (DMFT), provide a powerful framework for investigating strongly correlated materials. A central computational bottleneck in DMFT is in solving the Anderson impurity model (AIM), whose exact solution is classically intractable for large bath sizes. In this work, we benchmark a quantum-classical hybrid solver tailored for particle-hole symmetric AIMs, using the variational quantum eigensolver to prepare the ground state of the model with shallow quantum circuits. The solver uses shallow quantum ansätze and one set of variational parameters to prepare the ground state and its particle and hole excitations, enabling the construction of the impurity Green’s function through a continued-fraction expansion. We evaluate the performance of this approach across a few bath sizes and interaction strengths under noisy, shot-limited conditions. We compare three optimization routines (COBYLA, Adam, and L-BFGS-B) in terms of convergence and fidelity, assess the benefits of estimating a quantum-computed moment correction to the variational energies, and benchmark the approach by comparing the density of states computed from the impurity Green’s function against that obtained using a classical pipeline. Our results demonstrate the feasibility of Green’s function construction on near-term devices and establish practical benchmarks for quantum impurity solvers embedded within self-consistent DMFT loops.

Karabin, Mariia [ORNL]↗

Development of Li–S Battery Cells with High Energy Density and Long Cycling Life

Lithium–sulfur (Li–S) batteries are considered one of the most promising next-generation energy storage technologies for electric vehicles, owing to their high energy density (up to three times greater than current lithium-ion batteries) and the low cost of sulfur (approximately 100 times cheaper than conventional cobalt oxide). However, several challenges hinder their commercial viability, including the electronically and ionically insulating nature of sulfur, polysulfide dissolution, and lithium dendrite growth. These issues significantly limit battery energy density and cycle life and must be addressed to enable successful commercialization. Current strategies involving porous carbon structures, additives, and electrocatalyst engineering have not fully resolved the polysulfide dissolution problem. While polysulfide-free cathode materials offer a promising solution, they often suffer from low sulfur content, leading to reduced discharge capacity and sluggish reaction kinetics. To overcome these limitations, we have developed novel soluble-polysulfide-free sulfur cathode active materials with a high sulfur content (>50 wt%). By integrating these cathodes with functional binders, optimized electrolyte formulations, and refined electrode fabrication techniques, we demonstrate Li–S batteries with both high energy density and long cycle life. This advancement brings the realization of low-cost, high-energy Li–S batteries significantly closer to practical application.

25 ENERGY STORAGE↗

Effect of alloying elements on the $γ$’ antiphase boundary energy in Ni-base superalloys

Mechanical and fatigue performance of γ-γ' Ni-base superalloys are strongly affected by the antiphase boundary energy (APBE) of the γ' precipitates which, in turn, is dictated by the alloy's composition. Due to the multicomponent character of these alloys, establishing composition-APBE relationships are challenging, even though the qualitative effect of individual solutes on the APBE may be known. This work attempts to utilize density functional theory-based cluster expansion calculations to systematically assess the effect of composition on the APBE of the γ' phases in Ni-base superalloys. We aim to elucidate the influence of not only one single element but also multiple coexisting alloying elements on the γ' APBE. By explicit consideration of configurational disorder via Monte Carlo sampling, the effect of temperature on the APBE has also been analyzed. Overall, this work reveals that (1) effects of individual solute element M on the APBE energy obtained in an isolated, ternary condition (i.e. in Ni 3-x Al 1-y M x+y ) does not directly translate to a multi-solute case and that (2) the mutual synergistic interactions among different solute elements are not negligible. Based on the present results, an empirical master equation that predicts the APBE based on the composition of the γ' phases has been obtained.

36 MATERIALS SCIENCE↗

Physical, numerical, and computational challenges of modeling neutrino transport in core-collapse supernovae

The proposal that core collapse supernovae are neutrino driven is still the subject of active investigation more than 50 years after the seminal paper by Colgate and White. The modern version of this paradigm, which we owe to Wilson, proposes that the supernova shock wave is powered by neutrino heating, mediated by the absorption of electron-flavor neutrinos and antineutrinos emanating from the proto-neutron star surface, or neutrinosphere. Neutrino weak interactions with the stellar core fluid, the theory of which is still evolving, are flavor and energy dependent. The associated neutrino mean free paths extend over many orders of magnitude and are never always small relative to the stellar core radius. Thus, neutrinos are never always fluid like. Instead, a kinetic description of them in terms of distribution functions that determine the number density of neutrinos in the six-dimensional phase space of position, direction, and energy, for both neutrinos and antineutrinos of each flavor, or in terms of angular moments of these neutrino distributions that instead provide neutrino number densities in the four-dimensional phase-space subspace of position and energy, is needed. In turn, the computational challenge is twofold: (i) to map the kinetic equations governing the evolution of these distributions or moments onto discrete representations that are stable, accurate, and, perhaps most important, respect physical laws such as conservation of lepton number and energy and the Fermi–Dirac nature of neutrinos and (ii) to develop efficient, supercomputer-architecture-aware solution methods for the resultant nonlinear algebraic equations. In this review, we present the current state of the art in attempts to meet this challenge.

79 ASTRONOMY AND ASTROPHYSICS↗

Confirmation of the PPLB Derivative Discontinuity: Exact Chemical Potential at Finite Temperatures of a Model System

The landmark 1982 paper of Perdew, Parr, Levy, and Balduz (often called PPLB) laid the foundation for our modern understanding of the role of the derivative discontinuity in density functional theory, which drives much development to account for its effects. A simple model for the chemical potential at vanishing temperature played a crucial role in their argument. We investigate the validity of this model in the simplest non-trivial system to which it can be applied and which can be easily solved exactly, the Hubbard dimer. We find exact agreement in the crucial zero-temperature limit, and show the model remains accurate for a significant range of temperatures. We identify how this range depends on the strength of correlations. We extend the model to approximate free energies accounting for the derivative discontinuity, a feature missing in standard semilocal approximations. Here, we provide a correction to this approximation to yield even more accurate free energies. We discuss the relevance of these results for warm dense matter.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A General Model for the Ductility of Intermetallics Applied to Fe-Co Alloys

The mechanical properties, specifically ductility, of high performing soft magnets such as Fe-Co alloys are a limiting factor to their broader use in a number of systems. The understanding of the mechanical robustness in these materials is currently insufficient to be able to support the growing interest in applications such as magnetic shielding or electric motors. Fe-Co alloys provide the highest commercially available magnetic saturation and high magnetic permeability but have very poor ductility. The addition of vanadium to these alloys has allowed for significant commercialization and some ductility improvements, but the fundamental reasons for the observed improvements are not well understood. In most published work on mechanical properties in these alloys, the precise chemistry of the alloys investigated, often a critical aspect of intermetallics, is not reported or controlled and thermal history is unclear. This work creates a ductility model that is sensitive to changes in chemistry and can predict relative strain to failure as well as brittle fracture mode for intermetallics and is applied to Fe-Co alloys. Through the application of density functional theory (DFT), this model identifies defect stabilities, anti-phase boundary (APB) energies, the energy necessary to cross-slip, and cleavage energies and combines them through energetic competition to determine a relative failure strain. The model correctly predicts ductility improvements with the addition of vanadium as well as the transition from intergranular to transgranular cleavage, though more precise experiments are necessary to appropriately validate the various improvements observed.

36 MATERIALS SCIENCE↗