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,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Raman spectrum and polarizability of liquid water from deep neural networks

We introduce a scheme based on machine learning and deep neural networks to model the environmental dependence of the electronic polarizability in insulating materials. Application to liquid water shows that training the network with a relatively small number of molecular configurations is sufficient to predict the polarizability of arbitrary liquid configurations in close agreement with ab initio density functional theory calculations. In combination with a neural network representation of the interatomic potential energy surface, the scheme allows us to calculate the Raman spectra along 2-nanosecond classical trajectories at different temperatures for H 2 O and D 2 O. The vast gains in efficiency provided by the machine learning approach enable longer trajectories and larger system sizes relative to ab initio methods, reducing the statistical error and improving the resolution of the low-frequency Raman spectra. Decomposing the spectra into intramolecular and intermolecular contributions elucidates the mechanisms behind the temperature dependence of the low-frequency and stretch modes.

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

Structural effects of composition tuning in A-site disordered perovskite La0.5Li0.5−xMxTiO3 (M = Na, K) nanorods for fast interfacial transport for solid composite electrolyte design

Polymer-ceramic composite electrolytes (CPE) have emerged as promising replacements for currently used liquid organic electrolytes, which pose as safety hazards, in Li-ion batteries. Limits in conductivity enhancement in CPEs is often attributed to a high resistance for Li-ion transport along the interface due to the incompatibility of the polymer and ceramic phases. A clear understanding of the interfacial structure and how this impacts interfacial Li-ion transport is needed in order to efficiently design a CPE with optimal ionic conductivity. In this study, density functional theory (DFT) calculations, in conjunction with scanning transmission electron microscopy (STEM) and electron energy loss spectroscopy (EELS), are used to unveil the bulk and surface structure of La0.5Li0.5−xMxTiO3 (M = Na, K) (LMTO) nanorods, a newly-reported ceramic system with promising applications in CPEs, and the LMTO interactions with the poly(lithium sulfonyl (trifluoromethane sulfonyl)imide methacrylate) (p(MTFSILi)) polymerized ionic liquid. Substitution of lithium for sodium and potassium onto the A-site perovskite lattice is shown to increase the stability of the preferred pseudocubic perovskite phase because of their increased cation size which allows them to reduce the TiO6 octahedral rotations in the bulk and at the surface. STEM and EELS results show that LMTO nanorods with differing compositions have Ti-enriched (110)-oriented surfaces. Increased sodium and potassium compositions in LMTO is also shown using the DFT calculations to weaken the binding of the lithium atom to the MTFSI unit when LiMTFSI is adsorbed to the LMTO surface. This weakening of lithium binding to the polymer at the LMTO interface indicates that the lithium mobility is increased, correlating to an increase in interfacial Li-ion transport. Overall, this work provides insights into how to tune the interfacial interactions between the polymer (p(MTFSI)) and ceramic (LMTO) for optimal CPE design.

Shepard, Lauren B [Pennsylvania State University, ↗
DOE OSTI2020

Liquid to crystal Si growth simulation using machine learning force field

Machine learning force field (ML-FF) has emerged as a potential promising approach to simulate various material phenomena for large systems with ab initio accuracy. However, most ML-FFs have been used to study the phenomena relatively close to the equilibrium ground states. In this work, we have studied a far from equilibrium system of liquid to crystal Si growth using ML-FF. Here, we found that our ML-FF based on ab initio decomposed atomic energy can reproduce all the aspects of ab initio simulated growth, from local energy fluctuations to transition temperatures, to diffusion constant, and growth rates. We have also compared the growth simulation with the Stillinger-Weber classical force field and found significant differences. A procedure is also provided to correct a systematic fitting bias in the ML-FF training process, which exists in all training models, otherwise critical results like transition temperature will be wrong.

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

Carbon nanotube thermoelectric devices by direct printing: Toward wearable energy converters

Thermoelectric devices convert thermal energy to electrical energy and are particularly well-suited for energy harvesting from waste heat. Even as the number of electronic devices used in daily life proliferates, technical advances diminish the average power such devices require to perform a given function. Here, localized thermal gradients that abound in our living environments, despite having modest energy densities, are therefore becoming increasingly viable and attractive to power such devices. With this motivation, we report the design, fabrication, and characterization of single-wall carbon nanotube thermoelectric devices (CNT-TDs) on flexible polyimide substrates as a basis for wearable energy converters. Our aqueous-solution-based film fabrication process could enable readily scalable, low-cost TDs; here, we demonstrate CNT-hydroxypropyl cellulose (HPC) composite thermoelectric films by aerosol jet printing. The electrical conductivity of the composite films is controlled through the number of CNT/HPC layers printed in combination with control of the annealing conditions.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2021

Influence of solution phase environmental heterogeneity and fluctuations on vibronic spectra: Perylene diimide molecular chromophore complexes in solution

Ensembles of ab initio parameterized Frenkel-exciton model Hamiltonians for different perylene diimide dimer systems are used, together with various dissipative quantum dynamics approaches, to study the influence of the solvation environment and fluctuations in chromophore relative orientation and packing on the vibronic spectra of two different dimer systems: a π-stacked dimer in aqueous solution in which the relative chromophore geometry is strongly confined by a phosphate bridge and a side-by-side dimer in dichloromethane involving a more flexible alkyne bridge that allows quasi-free rotation of the chromophores relative to one another. These entirely first-principles calculations are found to accurately reproduce the main features of the experimental absorption spectra, providing a detailed mechanistic understanding of how the structural fluctuations and environmental interactions influence the vibronic dynamics and spectroscopy of solutions of these multi-chromophore complexes.

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

Screening of bimetallic electrocatalysts for water purification with machine learning

Electrocatalysis provides a potential solution to NO 3 - pollution in wastewater by converting it to innocuous N 2 gas. However, materials with excellent catalytic activity are typically limited to expensive precious metals, hindering their commercial viability. Here, in response to this challenge, we have conducted the most extensive computational search to date for electrocatalysts that can facilitate NO 3 - reduction reaction, starting with 59 390 candidate bimetallic alloys from the Materials Project and Automatic-Flow databases. Using a joint machine learning- and computation-based screening strategy, we evaluated our candidates based on corrosion resistance, catalytic activity, N 2 selectivity, cost, and the ability to synthesize. We found that only 20 materials will satisfy all criteria in our screening strategy, all of which contain varying amounts of Cu. Our proposed list of candidates is consistent with previous materials investigated in the literature, with the exception of Cu–Co and Cu–Ag based compounds that merit further investigation.

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

SQMBox: Interfacing a semiempirical integral library to modular ab initio electronic structure enables new semiempirical methods

Ab initio and semiempirical electronic structure methods are usually implemented in separate software packages or use entirely different code paths. As a result, it can be time-consuming to transfer an established ab initio electronic structure scheme to a semiempirical Hamiltonian. Here we present an approach to unify ab initio and semiempirical electronic structure code paths based on a separation of the wavefunction ansatz and the needed matrix representations of operators. With this separation, the Hamiltonian can refer to either an ab initio or semiempirical treatment of the resulting integrals. We built a semiempirical integral library and interfaced it to the GPU-accelerated electronic structure code TeraChem. Equivalency between ab initio and semiempirical tight-binding Hamiltonian terms is assigned according to their dependence on the one-electron density matrix. The new library provides semiempirical equivalents of the Hamiltonian matrix and gradient intermediates, corresponding to those provided by the ab initio integral library. This enables the straightforward combination of semiempirical Hamiltonians with the full pre-existing ground and excited state functionality of the ab initio electronic structure code. We demonstrate the capability of this approach by combining the extended tight-binding method GFN1-xTB with both spin-restricted ensemble-referenced Kohn–Sham and complete active space methods. We also present a highly efficient GPU implementation of the semiempirical Mulliken-approximated Fock exchange. The additional computational cost for this term becomes negligible even on consumer-grade GPUs, enabling Mulliken-approximated exchange in tight-binding methods for essentially no additional cost.

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

Solvation of furfural at metal–water interfaces: Implications for aqueous phase hydrogenation reactions

Metal–water interfaces are central to understanding aqueous-phase heterogeneous catalytic processes. However, the explicit modeling of the interface is still challenging as it necessitates extensive sampling of the interfaces’ degrees of freedom. Herein, we use ab initio molecular dynamics (AIMD) simulations to study the adsorption of furfural, a platform biomass chemical on several catalytically relevant metal–water interfaces (Pt, Rh, Pd, Cu, and Au) at low coverages. We find that furfural adsorption is destabilized on all the metal–water interfaces compared to the metal–gas interfaces considered in this work. This destabilization is a result of the energetic penalty associated with the displacement of water molecules near the surface upon adsorption of furfural, further evidenced by a linear correlation between solvation energy and the change in surface water coverage. To predict solvation energies without the need for computationally expensive AIMD simulations, we demonstrate OH binding energy as a good descriptor to estimate the solvation energies of furfural. Using microkinetic modeling, we further explain the origin of the activity for furfural hydrogenation on intrinsically strong-binding metals under aqueous conditions, i.e., the endothermic solvation energies for furfural adsorption prevent surface poisoning. In conclusion, our work sheds light on the development of active aqueous-phase catalytic systems via rationally tuning the solvation energies of reaction intermediates.

aqueous-phase hydrogenation↗
DOE OSTI2023

Research on the mechanism of interactions between Li/Na/K atoms and electrode materials

As the demand for high-capacity battery materials continues to grow, the exploration of the interaction mechanisms between Li/Na/K atoms and electrode materials has gained tremendous attention. In this study, we introduce a descriptor (ɛɑ) aimed at efficiently identifying high-performance electrode materials, specifically from two-dimensional transition metal chalcogenides and C and N compounds (MXenes). Our density functional theory calculations disclose a robust linear relationship between this descriptor and the adsorption energy (Ead), represented as Ead = ɛɑ + b. The fitting parameters, ɛɑ and b, are individually linked to the inherent properties of the substrate and adsorbate. Explicitly, the descriptor ɛɑ represents the capacity of the electrode material to accommodate an extra electron, while the intercept b is determined by the ionization energy (EIP) of alkali atoms and the coupling energy (Ecp) between the cation to the negatively charged substrate. This finding serves as a practical guide for designing high-capacity battery materials and paves the way for future advancements in this field.

Physics↗
DOE OSTI2025

Embedded random phase approximation for magnetic systems: H 2 dissociative adsorption on Fe(110)

The random phase approximation (RPA), a method for treating electron correlation, has been shown to be superior to standard density functional theory (DFT) approximations in numerous cases. However, the RPA’s computational cost is substantially higher than that of DFT, particularly restricting its application to extended surfaces. The recently introduced embedded RPA (emb-RPA) approach [Wei et al., J. Chem. Phys. 159(19), 194108 (2023)] reduces this computational cost by approximately two orders of magnitude. While previous applications of emb-RPA focused on non-spin-polarized systems, here we extend the approach to ferromagnetic ones. Unlike other embedded correlated wavefunction methods, such as embedded complete active space self-consistent field theory, emb-RPA is advantageous for spin-polarized systems because the RPA is compatible with unrestricted DFT solutions, which are eigenfunctions of the spin angular momentum operator S z but not the total spin-squared operator S 2 . By applying emb-RPA with specific magnetization constraints, we achieved a speedup of two to three orders of magnitude (one order when accounting for the one-time embedding potential optimization cost) with only small errors (∼50 meV) compared to full periodic RPA. Moreover, emb-RPA significantly reduces the over-binding errors of DFT approximations. In conclusion, we anticipate that the acceleration enabled by the spin-polarized emb-RPA approach will broaden the applicability of RPA to magnetic materials.

Density functional theory↗
DOE OSTI2020

Selective reduction of CO to acetaldehyde with CuAg electrocatalysts

Electrochemical CO reduction can serve as a sequential step in the transformation of CO 2 into multicarbon fuels and chemicals. In this study, we provide insights on how to steer selectivity for CO reduction almost exclusively toward a single multicarbon oxygenate by carefully controlling the catalyst composition and its surrounding reaction conditions. Under alkaline reaction conditions, we demonstrate that planar CuAg electrodes can reduce CO to acetaldehyde with over 50% Faradaic efficiency and over 90% selectivity on a carbon basis at a modest electrode potential of −0.536 V vs. the reversible hydrogen electrode. The Faradaic efficiency to acetaldehyde was further enhanced to 70% by increasing the roughness factor of the CuAg electrode. Density functional theory calculations indicate that Ag ad-atoms on Cu weaken the binding energy of the reduced acetaldehyde intermediate and inhibit its further reduction to ethanol, demonstrating that the improved selectivity to acetaldehyde is due to the electronic effect from Ag incorporation. These findings will aid in the design of catalysts that are able to guide complex reaction networks and achieve high selectivity for the desired product.

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

Erosion and cathodic arc plasma of Nb–Al cathodes: composite versus intermetallic

Many properties of cathodic arcs from single-element cathodes show a correlation to the cohesive energy of the cathode material. For example, the burning voltage, the erosion rate, or, to a lesser extent, plasma properties like electron temperatures or average ion energy and charge states. For multi-element cathodes, various phases with different cohesive energies can initially be present in the cathode, or form due to arc exposure, complicating the evaluation of such correlations. To test the influence of morphology and phase composition of multi-element cathodes on cathodic arc properties, a Nb–Al cathode model system was used that includes: pure Nb and Al cathodes; intermetallic Nb 3 Al, Nb 2 Al and NbAl 3 cathodes; and three composite Nb–Al cathodes with atomic ratios corresponding to the stoichiometric ratios of the intermetallic phases. Pulsed cathodic arc plasmas from these cathodes were examined using a mass-per-charge and energy-per-charge analyzer, showing that charge-state-resolved ion energy distributions of plasmas from the intermetallic and corresponding composite cathodes are nearly identical. An examination of converted layers of eroded cathodes using x-ray diffraction and scanning electron microscopy indicates the formation of a surface layer with similar phase composition for intermetallic and their corresponding composite cathode types. The average arc voltages do not follow the trend of cohesive energies of Nb, Al and intermetallic Nb–Al phases, which have been calculated using density functional theory. Possible reasons for this effect are discussed based on the current knowledge of multi-element arc cathodes and their arc plasma available in literature.

36 MATERIALS SCIENCE↗
DOE OSTI2026

From Spin to Pseudospin Symmetry: The Origin of Magic Numbers in Nuclear Structure

Magic numbers lie at the heart of nuclear structure, reflecting enhanced stability in nuclei with closed shells. While the emergence of magic numbers beyond 20 is commonly attributed to strong spin-orbit coupling, the microscopic origin of the spin-orbit potential remains elusive, owing to its dependence on the resolution scale and renormalization scheme of nuclear forces. Here, we investigate the evolution of nuclear shell structure with varying momentum resolution in nuclear interactions derived from chiral effective field theory, using the similarity renormalization group to link different scales. We uncover a novel transition from spin symmetry to pseudospin symmetry as the resolution scale decreases, during which magic numbers emerge naturally. A similar pattern is found in calculations using relativistic one-boson-exchange potentials, underscoring the robustness of the phenomenon. This establishes a direct connection between realistic nuclear forces with a high resolution scale and effective nuclear forces at coarse-grained scales, offering a first-principles explanation for the origin of magic numbers and pseudospin symmetry in nuclear shell structure and new insights into the structure of exotic nuclei far from stability

Energy levels↗
DOE OSTI2025

Robustness of Vacancy-Bound Non-Abelian Anyons in the Kitaev Model in a Magnetic Field

Non-Abelian anyons in quantum spin liquids (QSLs) provide a promising route to fault-tolerant topological quantum computation. In the exactly solvable Kitaev honeycomb model, such anyons of the QSL state can be bound to nonmagnetic spin vacancies and endowed with non-Abelian statistics by an infinitesimal magnetic field. Here, we investigate how this approach for stabilizing non-Abelian anyons extends to a finite magnetic field represented by a proper Zeeman term. Through large-scale density-matrix renormalization group simulations, we compute the vacancy-anyon binding energy as a function of magnetic field for both the ferromagnetic and antiferromagnetic Kitaev models. Here, we find that anyon binding remains robust within the entire QSL phase for the ferromagnetic Kitaev model but breaks down already inside this phase for the antiferromagnetic Kitaev model. To compute a binding energy several orders of magnitude below the magnetic energy scale, we introduce both a refined definition and an extrapolation scheme based on carefully tailored perturbations.

Xiao, Bo [Oak Ridge National Laboratory (ORNL), Oa↗
DOE OSTI2020

Effects of filling, strain, and electric field on the Néel vector in antiferromagnetic CrSb

CrSb is a layered antiferromagnet (AFM) with perpendicular magnetic anisotropy, a high Néel temperature, and large spin-orbit coupling (SOC), which makes it interesting for AFM spintronic applications. To elucidate the various mechanisms of Néel vector control, the effects of strain, band filling, and electric field on the magnetic anisotropy energy (MAE) of bulk and thin-film CrSb are determined and analyzed using density functional theory. The MAE of the bulk crystal is large (1.2 meV per unit cell). Due to the significant ionic nature of the Cr-Sb bond, finite slabs are strongly affected by end termination. Truncation of the bulk crystal to a thin film with one surface terminated with Cr and the other surface terminated with Sb breaks inversion symmetry, creates a large charge dipole and average electric field across the film, and breaks spin degeneracy, such that the thin film becomes a ferrimagnet. The MAE is reduced such that its sign can be switched with realistic strain, and the large SOC gives rise to an intrinsic voltage controlled magnetic anisotropy. A slab terminated on both faces with Cr remains a compensated AFM, but with the compensation occurring nonlocally between mirror-symmetric Cr pairs. In-plane alignment of the moments is preferred, the magnitude of the MAE remains large, similar to that of the bulk, and it is relatively insensitive to filling.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Fathoming the anisotropic magnetoelasticity and magnetocaloric effect in GdNi

Intermetallic GdNi adopts a CrB type of crystal structure (space group Cmcm), and it orders ferromagnetically via a second-order phase transition at 70 K, exhibiting unusually strong spontaneous striction along the three independent crystallographic axes in the ferromagnetically ordered state. Here we introduce a microscopic model to describe anisotropic changes of lattice parameters and elastic contribution to magnetocaloric effect of GdNi. In the model, results of density functional theory (DFT) calculations are used as inputs into a Hamiltonian that includes elastic energy of an anisotropic crystal lattice, exchange interactions, and Zeeman effect. The magnetic and elastic Hamiltonians are coupled through an anisotropic Bean-Rodbell model of magnetoelastic interactions. This coupling gives rise to anisotropic changes in the lattice parameters observed experimentally, and the model reveals good to reasonable agreements between the current theoretical results and earlier experimental data, thus validating the model within the limits of assumptions made. We also show that DFT calculations with 4f electrons of Gd treated as core electrons lead to a more adequate estimate of elastic constants of GdNi in comparison with the LDA + U method where 4f electrons are treated as valence electrons.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Evolution of structure and spectroscopic properties of a new 1,3-diacetylpyrene polymorph with temperature and pressure

A new polymorph of 1,3-diacetylpyrene has been obtained from its melt and thoroughly characterized using single-crystal X-ray diffraction, steady-state UV–Vis spectroscopy and periodic density functional theory calculations. Experimental studies covered the temperature range from 90 to 390 K and the pressure range from atmospheric to 4.08 GPa. Optimal sample placement in a diamond anvil cell according to our previously presented methodology ensured over 80% data coverage up to 0.8 Å for a monoclinic sample. Unrestrained Hirshfeld atom refinement of the high-pressure crystal structures was successful and anharmonic behavior of carbonyl oxygen atoms was observed. Unlike the previously characterized polymorph, the structure of 2°AP-β is based on infinite π-stacks of antiparallel 2°AP molecules. 2°AP-β displays piezochromism and piezofluorochromism which are directly related to the variation in interplanar distances within the π-stacking. The importance of weak intermolecular interactions is reflected in the substantial negative thermal expansion coefficient of -55.8 (57) MK -1 in the direction of C—H∙∙∙O interactions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Rational design of Lewis base molecules for stable and efficient inverted perovskite solar cells

Lewis base molecules that bind undercoordinated lead atoms at interfaces and grain boundaries (GBs) are known to enhance the durability of metal halide perovskite solar cells (PSCs). Using density functional theory calculations, we found that phosphine-containing molecules have the strongest binding energy among members of a library of Lewis base molecules studied herein. Experimentally, we found that the best inverted PSC treated with 1,3-bis(diphenylphosphino)propane (DPPP), a diphosphine Lewis base that passivates, binds, and bridges interfaces and GBs, retained a power conversion efficiency (PCE) slightly higher than its initial PCE of ~23% after continuous operation under simulated AM1.5 illumination at the maximum power point and at ~40°C for >3500 hours. DPPP-treated devices showed a similar increase in PCE after being kept under open-circuit conditions at 85°C for >1500 hours.

14 SOLAR ENERGY↗