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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Machine Learned Hückel Theory: Interfacing Physics and Deep Neural Networks

The Hückel Hamiltonian is an incredibly simple tight-binding model known for its ability to capture qualitative physics phenomena arising from electron interactions in molecules and materials. Part of its simplicity arises from using only two types of empirically fit physics-motivated parameters: the first describes the orbital energies on each atom and the second describes electronic interactions and bonding between atoms. By replacing these empirical parameters with machine-learned dynamic values, we vastly increase the accuracy of the extended Hückel model. The dynamic values are generated with a deep neural network, which is trained to reproduce orbital energies and densities derived from density functional theory. The resulting model retains interpretability, while the deep neural network parameterization is smooth and accurate and reproduces insightful features of the original empirical parameterization. Altogether, this work shows the promise of utilizing machine learning to formulate simple, accurate, and dynamically parameterized physics models.

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

How well do one-electron self-interaction-correction methods perform for systems with fractional electrons?

Recently developed locally scaled self-interaction correction (LSIC) is a one-electron SIC method that, when used with a ratio of kinetic energy densities (z σ ) as iso-orbital indicator, performs remarkably well for both thermochemical properties as well as for barrier heights overcoming the paradoxical behavior of the well-known Perdew–Zunger self-interaction correction (PZSIC) method. In this work, we examine how well the LSIC method performs for the delocalization error. Our results show that both LSIC and PZSIC methods correctly describe the dissociation of $H$$^{+}_{2}$ and $H$$^{+}_{2}$ but LSIC is overall more accurate than the PZSIC method. Likewise, in the case of the vertical ionization energy of an ensemble of isolated He atoms, the LSIC and PZSIC methods do not exhibit delocalization errors. For the fractional charges, both LSIC and PZSIC significantly reduce the deviation from linearity in the energy vs number of electrons curve, with PZSIC performing superior for C, Ne, and Ar atoms while for Kr they perform similarly. The LSIC performs well at the endpoints (integer occupations) while substantially reducing the deviation. The dissociation of LiF shows both LSIC and PZSIC dissociate into neutral Li and F but only LSIC exhibits charge transfer from Li + to F – at the expected distance from the experimental data and accurate ab initio data. Overall, both the PZSIC and LSIC methods reduce the delocalization errors substantially.

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

3D reconstruction of an inertial-confinement fusion implosion with neural networks using multiple heterogeneous data sources

3D asymmetries are major degradation mechanisms in inertial-confinement fusion implosions at the National Ignition Facility (NIF). These asymmetries can be diagnosed and reconstructed with the neutron imaging system (NIS) on three lines of sight around the NIF target chamber. Conventional tomographic reconstructions are used to reconstruct the 3D morphology of the implosion using NIS [Volegov et al., J. Appl. Phys. 127, 083301 (2020)], but the problem is ill-posed with only three imaging lines of sight. Asymmetries can also be diagnosed with the real-time neutron activation diagnostics (RTNAD) and the neutron time-of-flight (nToF) suite. Since the NIS, RTNAD, and nToF each sample a different part of the implosion using different physical principles, we propose that it is possible to overcome the limitations of too few imaging lines of sight by performing 3D reconstructions that combine information from all three heterogeneous data sources. This work presents a new machine learning-based reconstruction technique to do just this. By using a simple physics model and group of neural networks to map 3D morphologies to data, this technique can easily account for data of multiple different types. A simple proof-of-principle is presented, demonstrating that this technique can accurately reconstruct a hot-spot shape using synthetic primary neutron images and a hot-spot velocity vector. In particular, the hot-spot’s asymmetry, quantified as spherical harmonic coefficients, is reconstructed to within ±4% of the radius in 90% of test cases. In the future, this technique will be applied to actual NIS, RTNAD, and nToF data to better understand 3D asymmetries at the NIF.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

R-matrix calculations for opacities: I. Methodology and computations

Abstract An extended version of the R -matrix methodology is presented for calculation of radiative parameters for improved plasma opacities. Contrast and comparisons with existing methods primarily relying on the distorted wave approximation are discussed to verify accuracy and resolve outstanding issues, particularly with reference to the opacity project (OP). Among the improvements incorporated are: (i) large-scale Breit–Pauli R -matrix calculations for complex atomic systems including fine structure, (ii) convergent close coupling wave function expansions for the ( e + ion) system to compute oscillator strengths and photoionization cross sections, (iii) open and closed shell iron ions of interest in astrophysics and experiments, (iv) a treatment for plasma broadening of autoionizing resonances as function of energy-temperature-density dependent cross sections, (v) a ‘top-up’ procedure to compare convergence with R -matrix calculations for highly excited levels, and (vi) spectroscopic identification of resonances and bound ( e + ion) levels. The present R -matrix monochromatic opacity spectra are fundamentally different from OP and lead to enhanced Rosseland and Planck mean opacities. An outline of the work reported in other papers in this series and those in progress is presented. Based on the present re-examination of the OP work, opacities of heavy elements might require revisions in high temperature-density plasma sources.

Pradhan, A. K. (ORCID:0000000187753643)↗
DOE OSTI2025

On-lattice kinetic Monte Carlo approaches for modeling molecular anisotropy in resveratrol crystallization

Stilbenes are a class of organic compounds with broad-ranging pharmaceutical and agricultural applications, which are typically isolated and purified through recrystallization. We are motivated by reducing experimental waste and optimizing yield via developing predictive simulations for processing-dependent crystal morphologies. Using resveratrol as a model stilbene system, we have developed an approach for simulating crystallization with molecular resolution using on-lattice kinetic Monte Carlo. In this work, we highlight modifications to the Stochastic Parallel PARticle Kinetic Simulator (SPPARKS) software package, which were essential to this application. Key enhancements include the incorporation of non-orthogonal cell shapes and monomer anisotropy approximations using bound hard spheres. This new SPPARKS application has been applied to resveratrol with attachment energy libraries obtained from density functional theory, resulting in excellent agreement with experimental morphology prediction.

crystallization↗
DOE OSTI2026

Machine-learned quantum molecular dynamics calculations of warm dense equation of state and ionic transport coefficients of deuterated water

White dwarf models require accurate equations of state and ionic transport coefficients in the warm dense matter regime, where kinetic theory models and tabulated equations of state are often inaccurate. In this work, spectral-partitioned density functional theory and machine-learned interatomic potentials are combined to perform large-scale, first-principles quantum molecular dynamics simulations of deuterated water (D 2 O) near the principal Hugoniot. This approach retains Kohn-Sham accuracy while achieving orders-of-magnitude speedup, yielding converged equation of state and transport properties over a broad pressure and temperature range. The results reveal the thermodynamic conditions under which ionic transport models for interdiffusivity and shear viscosity converge and identify those in closest agreement with density functional theory benchmarks at temperatures in the warm dense matter regime. The present framework extends first-principles transport calculations to higher temperatures than previously achieved, and provides an efficient, scalable, and general approach for studying transport properties in complex multicomponent mixtures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Thermoelectric Figure of Merit of a Superatomic Crystal Re 6 Se 8 I 2 Monolayer

Superatomic materials are newly emerging candidates for high-performance thermoelectric (TE) devices due to their intrinsic ultralow thermal conductivities. However, the low TE power factor becomes a huge obstacle to reaching the required the dimensionless figure of merit (ZT) values for practical applications. Here, motivated by the recently synthesized superatomic Re 6 Se 8 I 2 monolayer [J. Am. Chem. Soc. 144, 74 (2022)], we study its superior thermoelectric properties by using density functional theory combined with phonon Boltzmann transport theory and deep potential molecular dynamics. We show that the large mass and anharmonic Re-I bonds introduce strong phonon scattering and result in low lattice thermal conductivity of 1.20 W m -1 K -1 at 300 K, while the strong and harmonic Re-Se network ensures the high TE power factor of 4344 μW m -1 K -2 in b direction for n-type doping. This is an order of magnitude higher than that of other cluster-based materials. Due to the low thermal conductivity and high TE power factor, Re 6 Se 8 I 2 exhibits high ZT values of 1.20 (500 K) and 1.43 (900 K) with n-type doping along the b direction among all the cluster-based thermoelectric materials reported so far.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Phonons, Q-dependent Kondo spin fluctuations, and 4f phonon resonance in YbAl 3

4f intermediate valence (IV) compounds are canonical hosts of correlated electron physics and can contribute to our understanding of the larger class of correlated electron materials. In this work, we study the prototype IV compound YbAl 3 which exhibits a nonintegral valence with a moderately heavy fermion ground state and a large Kondo temperature (T K ~500–600K). To better characterize the correlated physics of YbAl 3 , we have measured the phonon and the magnetic excitation spectra on single crystals of this material by time-of-flight inelastic neutron scattering and inelastic x-ray scattering. We have also performed theoretical calculations of the phonon spectra. We present three findings of these measurements. First, we observe that the measured phonon spectra can be described adequately by a calculation based on standard DFT+U density functional theory. The calculated energies, however, are 10% too low compared to the measured energies. This discrepancy may reflect a hardening of the phonons due to dynamic 4f correlations. Second, the low-temperature spin fluctuations on the Kondo energy scale k B T K have a momentum (Q) dependence similar to that seen recently in the IV compound CePd 3 . For that system, the Q dependence has been attributed to particle-hole excitations in a coherent itinerant 4f correlated ground state. We suggest a similar origin for the momentum dependence seen in YbAl 3 . This Q dependence disappears as the temperature is raised towards room temperature and the 4f electron band states become increasingly incoherent. Such a coherent/incoherent crossover is expected to be generic for correlated electron systems. Third, a low-temperature magnetic peak observed in the neutron scattering near 30 meV shows dispersion identical to a particular optic-phonon branch. This 4f/phonon resonance disappears for T≥150K. The phonon spectrum appears to be unaffected by the resonance. We discuss several possibilities for the origin of this unusual excitation, which may be unique to YbAl 3 . We suggest that the excitation may arise from the large amplitude beating of the light Al atoms against the heavy Yb atoms, resulting in a dynamic 4f/3p hybridization.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Electronic states dressed by an out-of-plane supermodulation in the quasi-two-dimensional kagome superconductor CsV 3 Sb 5

CsV 3 Sb 5 has attracted much recent attention as the first quasi-two-dimensional (2D) kagome superconductor. While the kagome layers are 2D in nature, increasing evidence has pointed to the importance of out-of-plane correlation in this material. However, it remains unclear whether such correlation can change the fundamental electronic structure of the quasi-2D system. Here, we reveal this missing piece of information, using angle-resolved photoemission spectroscopy, complemented by scanning tunneling microscope measurements. The three-dimensional electronic structures in the high-temperature state are revealed, which agree well with density-functional theory calculations. Electron energy bands are observed in the low-temperature state that exhibit additional periodicities along the out-of-plane momentum. Furthermore, these results reveal a direct response to the out-of-plane electronic supermodulation in the single-particle spectral function of CsV 3 Sb 5 , thus establishing an electronic platform to examine emergent phenomena beyond 2D limit in kagome superconductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Excitonic Condensate in Flat Valence and Conduction Bands of Opposite Chirality

Excitonic Bose-Einstein condensation (EBEC) has drawn increasing attention recently with the emergence of 2D materials. A general criterion for EBEC, as expected in an excitonic insulator (EI) state, is to have negative exciton formation energies in a semiconductor. Here, using exact diagonalization of a multiexciton Hamiltonian modeled in a diatomic kagome lattice, we demonstrate that the negative exciton formation energies are only a prerequisite but insufficient condition for realizing an EI. By a comparative study between the cases of both conduction and valence flat bands (FBs) versus that of a parabolic conduction band, we further show that the presence and increased FB contribution to exciton formation provide an attractive avenue to stabilize the excitonic condensate, as confirmed by calculations and analyses of multiexciton energies, wave functions, and reduced density matrices. In conclusion, our results warrant a similar many-exciton analysis for other known and/or new candidates of EIs and demonstrate the FBs of opposite parity as a unique platform for studying exciton physics, paving the way to material realization of spinor BEC and spin superfluidity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Parallel Implicit Hydrodynamics with Material Strength for High Explosive Burn Calculations

High explosives are almost always evolving in some form of metal containment vessel. This fact requires that a materials model for the evolution of the metal containment vessel be part of any simulation of the HE. Since the actual form of the strength model to be used is an open question, we incorporate the material strength properties in a relatively agnostic fashion, which will accommodate many strength models. The only restriction is that the stress components be a function of density, specific internal energy and velocities. In addition, the timescales and rates of the chemical reactions in the HE vary between thousands of seconds and nanoseconds. A hydrodynamics capable of operating in anticipation of the eventual violent release of energy via these chemical reactions is a necessity. The central issue in creating a hydrodynamics capable of spanning these timescales is to forgo the use of a Courant time step control necessary in an explicit hydrodynamics. The Courant time-step arises because of characteristic velocities associated with the material, such as the sound speed or, in materials with strength, characteristic longitudinal and transverse speeds associated with compression and shear. The method used to circumvent the need for a Courant time-step limit is to develop an implicit calculation of the advanced particle pressure and shear modulus derived from the conservation laws of mass, momentum and specific internal energy. This process creates the implicit hydrodynamics needed for the HE calculations of interest to this study.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Theoretical studies of chemical reactions related to the formation and growth of polycyclic aromatic hydrocarbons (PAH) and molecular properties of their key intermediates (Final Progress Report)

The formation mechanisms of polycyclic aromatic hydrocarbons, (PAHs) – organic molecules carrying fused benzene rings – are of great interest to scientists and engineers due to their importance in combustion chemistry and astrochemistry. On Earth, PAHs are largely produced in incomplete combustion of fossil fuel and are considered as critical precursors to unwanted soot particles leading to combustion inefficiency and causing air pollution along with detrimental health effects. Simple PAH molecules initially formed in the gas phase, are further involved in a build-up process in combustion flames leading to larger PAH, bowl-shaped nanostructures, fullerenes, and solid-phase species including carbonaceous dust, graphene particles, and soot. In deep space, PAH and their derivatives are potential key intermediates and nucleation sites leading eventually to carbonaceous nanoparticles (“interstellar grains”). Therefore, the understanding of the key processes in the synthesis of PAHs along with their precursors and their degradation mechanisms in combustion systems and in interstellar, circumstellar, and planetary atmospheric environments will provide critical insights into how complex aromatic structures, carbonaceous nanoparticles, and fullerenes are formed and destroyed. Achieving this understanding is an important step in the development of the efficient combustion processes and of the ecofriendly devices with reduced environmental pollution as well as technological strategies for the production of hydrogen and solid carbon through thermal or plasma-assisted pyrolysis of natural gas and biomass. Also, the understanding of the key processes of PAH and soot growth will help in our comprehension of chemical evolution in the universe. Detailed information on the mechanisms and reliable rate constants of the key elementary chemical reactions involved in PAH formation and destruction processes and in inception of soot particles is often missing, with the main deficiencies being the absence of temperature- and pressure-dependent rate constants for the broad range of conditions occurring in various terrestrial and interstellar processes and the lack of data on the reaction products and their branching ratios. Complementary to experimental studies, these gaps in knowledge can be filled by using quantum chemical calculations of reaction potential energy surfaces providing us with accurate energies of reaction products, intermediates, and transition states, revealing the reaction mechanism, and giving the molecular properties required to compute rate constants for relevant reaction steps and product branching ratios using the RRKM-Master Equation (ME) method. Molecular dynamics (MD) simulations can be used in cases when a reaction rate cannot be properly described by statistical theories. During the terminal renewal project period we employed these ab initio/RRKM-ME and MD approaches to complete our studies on several key reactions relevant to the formation/growth of PAH and inception of soot particles including (1) the reaction mechanism and kinetics of the resonance stabilized fulvenallenyl radical with propargyl and C 3 H 4 isomers; (2) the reaction mechanism and kinetics for the C + indene and C 2 + styrene reactions producing naphthyl or azulenyl radicals in low-temperature environments; (3) the MD study of non-equilibrium dimerization of acepyrene and coronene and its radical. The information derived from our theoretical calculations contributed to a better fundamental understanding of the reaction mechanisms and provide missing critical kinetic data to improve combustion models of hydrocarbon fuels and astrochemical models of the growth of carbonaceous molecules and particles in cold molecular clouds, circumstellar envelopes, and planetary atmospheres.

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

Initializing BSQ with Open-Source ICCING

While it is well known that there is a significant amount of conserved charges in the initial state of nuclear collisions, the production of these due to gluon splitting has yet to be thoroughly investigated. The ICCING (Initial Conserved Charges in Nuclear Geometry) algorithm reconstructs these quark distributions, providing conserved strange, baryon, and electric charges, by sampling a given model for the g → qq¯ splitting function over the initial energy density, which is valid at top collider energies, even when µB = 0. The ICCING algorithm includes fluctuations in the gluon longitudinal momenta, a structure that supports the implementation of dynamical processes, and the c++ version is now open-source. A full analysis of parameter choices on the model has been done to quantify the effect these have on the underlying physics. We find there is a sustained difference across the different charges that indicates sensitivity to hot spot geometry.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Atomic structure of different surface terminations of polycrystalline ZnPd

The intermetallic compound ZnPd has been found to have desirable characteristics as a catalyst for the steam reforming of methanol. The understanding of the surface structure of ZnPd is important to optimize its catalytic behavior. However, due to the lack of bulk single-crystal samples and the complexity of characterizing surface properties in the available polycrystalline samples using common experimental techniques, all previous surface science studies of this compound have been performed on surface alloy samples formed through thin-film deposition. In this study, we present findings on the chemical and atomic structure of the surfaces of bulk polycrystalline ZnPd studied by a variety of complementary experimental techniques, including scanning tunneling microscopy (STM), x-ray photoelectron spectroscopy (XPS), low energy electron microscopy (LEEM), photoemission electron microscopy (PEEM), and microspot low-energy electron diffraction ( μ -LEED). These experimental techniques, combined with density functional theory (DFT)-based thermodynamic calculations of surface free energy and detachment kinetics at the step edges, confirm that surfaces terminated by atomic layers composed of both Zn and Pd atoms are more stable than those terminated by only Zn or Pd layers. DFT calculations also demonstrate that the primary contribution to the tunneling current arises from Pd atoms, in agreement with the STM results. The formation of intermetallics at surfaces may contribute to the superior catalyst properties of ZnPd over Zn or Pd elemental counterparts. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗
DOE OSTI2022

Xenon mobility in γ-uranium and uranium–molybdenum alloys

Diffusion in bcc uranium and U–Mo alloys is of great interest because fission gas and other fission products impact the performance of nuclear fuels. We investigate the mobility of xenon and molybdenum in bcc uranium (⁠$γ$-U) and metallic U–Mo alloys by calculating the migration energies of xenon and molybdenum for various local compositions using density functional theory. We also calculate the solute–vacancy binding energies of different solutes to vacancies in bcc uranium. We find that the solute–vacancy binding energy in bcc uranium is significantly higher than it is in other bcc metals (e.g., Fe and W). We also find that the migration energy of molybdenum is substantially higher than the migration energy of xenon, indicating that xenon is much more mobile than molybdenum in bcc uranium. The presence of molybdenum in the nearest-neighbor shell around a xenon atom typically increases the migration energy of xenon, which indicates a reduction of xenon mobility in U–Mo alloys compared to pure bcc uranium.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Polarization consistent dielectric screening in polarizable continuum model calculations of solvation energies

A polarization consistent framework, where dielectric screening is affected consistently in polarizable continuum model (PCM) calculations, is employed for the study of solvation energies. The computational framework combines a screened range-separated-hybrid functional (SRSH) with PCM calculations, SRSH-PCM, where dielectric screening is imposed in both PCM self-consistent reaction field (SCRF) iterations and the electronic structure Hamiltonian. We begin by demonstrating the impact of modifying the Hamiltonian to include such dielectric screening in SCRF iterations by considering the solutions of electrostatically embedded Hartree–Fock (HF) exact exchange equations. Long-range screened HF-PCM calculations are shown to capture properly the linear dependence of gap energy of frontier orbitals on the inverse of the dielectric constant, whereas unscreened HF-PCM orbital energies are fallaciously semi-constant with respect to the dielectric constant and, therefore, inconsistent with the ionization energy gaps. Similar trends affect density functional theory (DFT) calculations that aim to achieve predictive quality. Importantly, the dielectric screened calculations are shown to significantly affect DFT- and HF PCM-based solvation energies, where screened solvation energies are smaller compared to the unscreened values. Importantly, SRSH-PCM, therefore, appears to reduce the tendency of DFT-PCM to overestimate solvation energies, where we find the effect to increase with the dielectric constant and the polarity of the molecular solute, trends that enhance the quality of DFT-PCM calculations of solvation energy. Understanding the relationship of dielectric screening in the Hamiltonian and DFT-PCM calculations can ultimately benefit on-going efforts for the design of predictive and parameter free descriptions of solvation energies.

Chemistry↗