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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 523 records · Page 29

All sources NASA NTRS DOE OSTI

DOE OSTI2023

A Deep Potential model for liquid–vapor equilibrium and cavitation rates of water

Computational studies of liquid water and its phase transition into vapor have traditionally been performed using classical water models. Here, we utilize the Deep Potential methodology—a machine learning approach—to study this ubiquitous phase transition, starting from the phase diagram in the liquid–vapor coexistence regime. The machine learning model is trained on ab initio energies and forces based on the SCAN density functional, which has been previously shown to reproduce solid phases and other properties of water. Here, we compute the surface tension, saturation pressure, and enthalpy of vaporization for a range of temperatures spanning from 300 to 600 K and evaluate the Deep Potential model performance against experimental results and the semiempirical TIP4P/2005 classical model. Moreover, by employing the seeding technique, we evaluate the free energy barrier and nucleation rate at negative pressures for the isotherm of 296.4 K. Further, we find that the nucleation rates obtained from the Deep Potential model deviate from those computed for the TIP4P/2005 water model due to an underestimation in the surface tension from the Deep Potential model. From analysis of the seeding simulations, we also evaluate the Tolman length for the Deep Potential water model, which is (0.091 ± 0.008) nm at 296.4 K. Finally, we identify that water molecules display a preferential orientation in the liquid–vapor interface, in which H atoms tend to point toward the vapor phase to maximize the enthalpic gain of interfacial molecules. We find that this behavior is more pronounced for planar interfaces than for the curved interfaces in bubbles. This work represents the first application of Deep Potential models to the study of liquid–vapor coexistence and water cavitation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Machine-Learned Force Field Modeling of Metal Organic Frameworks for CO2 Direct Air Capture

To cope with legacy greenhouse gas emissions and to achieve net-zero emissions by 2050, the U.S. Department of Energy (DOE) is funding efforts to develop direct air capture (DAC), a method for removing CO2 directly from air. Metal organic frameworks (MOFs) have been studied as DAC sorbent materials because of their structural and chemical diversity. Thermodynamic calculations using classical force fields are often used to screen MOFs for their performance in separations such as CO2 capture. Machine-learned force fields (MLFFs) can use machine learning to form quantitative relationships between a material’s chemical structure and the forces and energies predicted by more accurate quantum mechanical calculations, such as dispersion-corrected density functional theory (DFT). These descriptions of forces and energies can be used to improve the accuracy of adsorption calculations. In this work, MLFF models were developed for MOFs to achieve DFT-level accuracy for the forces and energies associated with MOF flexibility and CO2 adsorption. These methods were parametrized based on thousands of DFT calculations of CO2 in flexible MOFs and used to predict MOF structural properties as well as CO2 adsorption properties.

Findley, John↗

DOE OSTI2024

Machine Learned Force Field Modeling of Metal Organic Frameworks for CO2 Direct Air Capture

Direct air capture (DAC) is a method for removing CO2 directly from air. Metal organic frameworks (MOFs) have been studied as DAC sorbent materials because of their structural and chemical diversity. Thermodynamic calculations using classical force fields are often used to evaluate MOFs for their performance in separations such as CO2 capture. Machine-learned force fields (MLFFs) can use machine learning to form quantitative relationships between a material’s chemical structure and the forces and energies predicted by more accurate quantum mechanical calculations, such as dispersion-corrected density functional theory (DFT). These descriptions of forces and energies can be used to improve the accuracy of adsorption calculations. In this work, classical models were used to pre-screen MOFs for CO2 capture. DFT calculations were then used to examine the adsorption mechanism. Next, MLFF models were developed for MOFs to achieve DFT-level accuracy for the forces and energies associated with MOF flexibility and CO2 adsorption. These methods were parametrized based on thousands of DFT calculations of CO2 in flexible MOFs and used to predict MOF structural properties as well as CO2 adsorption properties.

Findley, John↗

DOE OSTI2021

Xenon Trapping in Metal-Supported Silica Nanocages

Xenon (Xe) is a valuable and scarce noble gas used in various applications, including lighting, electronics, and anesthetics, among many others. It is also a volatile byproduct of the nuclear fission of uranium. In this work, a novel material architecture consisting of silicate nanocages in contact with a metal surface and an approach for trapping single Xe atoms in these cages is presented. The trapping is done at low Xe pressures and temperatures between 400 and 600 K, and the process is monitored in situ using synchrotron-based ambient pressure X-ray photoelectron spectroscopy. Release of the Xe from the cages occurs only when heating to temperatures above 750 K. A model that explains the experimental trapping kinetics is proposed and tested using Monte Carlo methods. Density functional theory calculations show activation energies for Xe exiting the cages consistent with experiments. This work can have significant implications in various fields, including Xe production, nuclear power, nuclear waste remediation, and nonproliferation of nuclear weapons. The results are also expected to apply to argon, krypton, and radon, opening an even more comprehensive range of applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Genesis of Active Pt/CeO 2 Catalyst for Dry Reforming of Methane by Reduction and Aggregation of Isolated Platinum Atoms into Clusters

Atomically dispersed metal catalysts offer the advantages of efficient metal utilization and high selectivities for reactions of technological importance. Such catalysts have been suggested to be strong candidates for dry reforming of methane (DRM), offering prospects of high selectivity for synthesis gas without coke formation, which requires ensembles of metal sites and is a challenge to overcome in DRM catalysis. However, investigations of the structures of isolated metal sites on metal oxide supports under DRM conditions are lacking, and the catalytically active sites remain undetermined. Data characterizing the DRM reaction-driven structural evolution of a cerium oxide- supported catalyst, initially incorporating atomically dispersed platinum, and the corresponding changes in catalyst performance are reported. X-ray absorption and infrared spectra show that the reduction and agglomeration of isolated cationic platinum atoms to form small platinum clusters/nanoparticles are necessary for DRM activity. Density functional theory calculations of the energy barriers for methane dissociation on atomically dispersed platinum and on platinum clusters support these observations. Here the results emphasize the need for in-operando experiments to assess the active sites in such catalysts. The inferences about the catalytically active species are suggested to pertain to a broad class of catalytic conversions involving the rate-limiting dissociation of light alkanes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Absorptive effects and classical black hole scattering

Abstract We describe an approach to incorporating the physical effects of the absorption of energy by the event horizon of black holes in the scattering amplitudes based post-Minkowskian, point-particle effective description. Absorptive dynamics are incorporated in a model-independent way by coupling the usual point-particle description to an invisible sector of gapless internal degrees-of-freedom. The leading order dynamics of this sector are encoded in the low-energy expansion of a spectral density function obtained by matching an absorption cross section in the ultraviolet description. This information is then recycled using the scattering amplitudes based Kosower-Maybee-O’Connell in-in formalism to calculate the leading absorptive contribution to the impulse and change in rest mass of a Schwarzschild black hole scattering with a second compact body sourcing a massless scalar, electromagnetic or gravitational field. The results obtained are in complete agreement with previous worldline Schwinger-Keldysh calculations and provide an alternative on-shell scattering amplitudes approach to incorporating horizon absorption effects in the gravitational two-body problem.

Physics↗

DOE OSTI2020

Passive NO x adsorption on Pd/H-ZSM-5: Experiments and modeling

A combined experimental and modeling study of the passive NO x adsorber (PNA) is presented that advances the understanding and prediction of the effects of various operating parameters and material properties of H-ZSM-5 and Pd/H-ZSM-5. Experiments reveal the NO x uptake process to be kinetically limited. Two microkinetic schemes are developed to explain the NO x uptake and release and NO 2 generation data. Both schemes are incorporated into a 1+1 D monolith model to predict and validate NO x uptake and temperature programmed desorption (TPD) data. Scheme I involves Z – [PdOH] + , Z – Pd 2+ Z – and Z – Pd + with Z – Pd + as the strongest NO binding site. Scheme II involves reduction of PdO 2 to PdO and nitrate formation. Model tuning utilizes a combination of infrared spectroscopy measurements and density functional theory (DFT) estimates of energy barriers. Lastly, the model is validated at different conditions including feed parameters and Pd loading.

42 ENGINEERING↗

DOE OSTI2022

Measuring the maximum capacity and thermal resistances in phase-change thermal storage devices

Thermal energy storage can increase the efficiency of the electric grid by adding flexibility to thermal systems. The value of thermal storage is a function of its energy and power density, which are driven by the capacity and thermal resistances in the storage device. Measuring these properties in-situ at the device level is an important step to understanding the performance and improving the design of thermal storage systems. Here, we present methods to measure the total capacity and thermal resistances in heat exchangers with integrated phase change materials. These methods are demonstrated on two thermal storage devices - a 570-kWh ice-based storage tank and a 0.35-kWh graphite-tetradecane composite device. The results show how thermal resistances evolve with the state of charge and discharge rate in these devices and quantify the impact of applied pressure on the contact resistance in composite phase change material heat exchangers. The proposed method allows for easy comparison between different systems and provides information on the thermal bottlenecks limiting performance. Ultimately, these measurements will allow designers to make robust, high-performance thermal storage devices for next-generation thermal systems.

25 ENERGY STORAGE↗

DOE OSTI2024

Magnetic ordering in Eu 2 In and Eu 2 Sn

Eu 2 In and Eu 2 Sn crystallize in the orthorhombic Co 2 Si-type structure (oP12, Pnma, No. 62) with In and Sn atoms occupying one 4c site and the Eu atoms filling two other 4c sites. Eu 2 In has a nearly ideal first-order magnetostructural transition (FOMT) at 55 K with a hysteresis of less than 0.1 K, a large entropy change and an adiabatic temperature change of 5.0 K in a field of 2 T. The anhysteretic nature of the FOMT is likely due to there being no change in cell symmetry and relatively small changes in the lattice parameters. There is no magnetostructural transition in Eu2Sn. In this work we present the results of powder neutron diffraction, magnetization, and Eu Mössbauer spectroscopy aimed to investigate the nature of magnetic order for both Eu 2 In and Eu 2 Sn. The Eu Mössbauer spectrum of Eu 2 In at 5 K shows two equal area components, consistent with Eu occupying two equal multiplicity crystallographic sites. However, the different hyperfine fields (B) of 27 T and 17 T suggest that the magnetic environments of the Eu moments on the two 4c sites are different. Neutron diffraction data at 2.5 K show that in Eu 2 In the order is ferromagnetic, with Eu moments on both Eu sites oriented parallel to the a-axis; moment values of 6.8 μB and 6.5 μB were found. For Eu 2 Sn measurements find two antiferromagnetic transitions, which are corroborated by neutron diffraction. Analysis of density-functional theory calculations shows negligible energy difference between differing magnetic configurations, indirectly supporting stability of multiple magnetic structures observed experimentally. While the transition at T N1 = 30 K corresponds to the formation of a simple k 1 = 0 antiferromagnetic structure with Eu-moments pointing along the b-axis, at T N2 = 13 K a coexisting second magnetic order with k 2 = [0, ½, ½] appears.

36 MATERIALS SCIENCE↗

DOE OSTI2023

A conservative Galerkin solver for the quasilinear diffusion model in magnetized plasmas

We propose a conservative Galerkin scheme for the quasilinear model in three-dimensional momentum space and three-dimensional spectral space, with cylindrical symmetry. We construct an unconditionally conservative weak form and use a discretization that preserves conservation properties independent of the wave emission probability. The discrete operators, combined with a consistent quadrature rule, preserve all the conservation laws rigorously. The proposed scheme is quite general: it works for both relativistic and non-relativistic systems, for both magnetized and unmagnetized plasmas, and even for problems with time-dependent dispersion relations. We represent the particle distribution by continuous basis functions and use discontinuous basis functions for the wave spectral energy density, which enables the application of a positivity-preserving technique. We adopt the marching simplex algorithm, designed initially for computer graphics, for numerical integration on the resonance manifold. Furthermore, the numerical examples with a bump-on-tail initial configuration show how the unstable waves produce strong momentum space diffusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Thermodynamic evaluation and modeling of the CeF 3 -LaF 3 -LiF-NaF system

The thermodynamic properties of the LiF-NaF-CeF 3 -LaF 3 pseudo-quaternary system have been evaluated, and models for the relevant phases have been revised or newly developed. The solid solubility of NaF in LaF 3 and CeF 3 was described in the model based on reported phase-equilibria observations and results from density functional theory calculations. The Gibbs energy functions for the NaCeF 4 and NaLaF 4 phases have been updated with recently reported thermochemical measurements. The Gibbs energy functions and model relations for the LiF-NaF-LaF 3 and LiF-NaF-CeF 3 pseudo-ternary systems were successfully optimized based on thermochemical properties and phase equilibria observations.

CALPHAD↗

DOE OSTI2026

A Transferable Force Field for Simulating Adsorption in Metal–Organic Frameworks with Open Metal Sites Based on the 12–6–4 Lennard-Jones Potential

Metal−organic frameworks (MOFs) that contain coordinatively unsaturated open metal sites (OMSs) provide strong host− guest interactions, making them promising sorbents for low-concentration gas adsorption applications such as direct air capture and atmospheric water harvesting. However, accurately modeling host−guest interactions involving OMSs remains challenging for classical force fields (FFs) based on the 12−6 Lennard−Jones (LJ) potential, as the polarization effect of the guest molecule induced by the positively charged OMS is not considered. Here, we introduce an FF based on the 12−6−4 LJ potential, which incorporates charge−induced dipole interactions and is parametrized against a diverse set of host−guest potential energy surfaces (PESs) obtained from density functional theory (DFT). The resulting FF, trained on a generic trimetallic cluster, performs well in both host−guest binding energetics and gas adsorption isotherms across different OMS-containing MOFs, including MOF-74 series and Cu-BTC. These results highlight the excellent transferability of our approach and its potential to enhance the accuracy and robustness of high-throughput MOF discovery workflows, particularly for gas adsorption and separation in large and diverse MOF databases.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Nitromethane Decomposition via Automated Reaction Discovery and an Ab Initio Corrected Kinetic Model

In the explore the systematic construction of kinetic models from in silico reaction data for the decomposition of nitromethane. Our models are constructed in a computationally affordable manner by using reactions discovered through accelerated molecular dynamics simulations using the ReaxFF reactive force field. The reaction paths are then optimized to determine reaction rate parameters. We introduce a reaction barrier correction scheme that combines accurate thermochemical data from density functional theory with ReaxFF minimal energy paths. We validate our models across different thermodynamic regimes, showing predictions of gas phase CO and NO concentrations and high-pressure induction times that are similar to experimental data. The kinetic models are analyzed to find fundamental decomposition reactions in different thermodynamic regimes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Accurate Force Field for Carbon Dioxide–Silica Interactions Based on Density Functional Theory

Fluid–silica interfaces are ubiquitous in chemistry, occurring in both natural geochemical environments and practical applications ranging from separations to catalysis. Simulations of these interfaces have been, and continue to be, a significant avenue for understanding their behavior. A constraining factor, however, is the availability of accurate force fields. Most simulations use traditional “mixing rules” to determine nonbonded dispersion interactions, an approach that has not been critically examined. Here, in this study, we present Lennard-Jones parameters for the interaction of carbon dioxide with silica interfaces that are optimized to reproduce density functional theory (DFT)-based binding energies. The modeling is based on the recently developed silica-DDEC force field, whose atomic charges are consistent with DFT calculations. Standard mixing rules are found to predict weaker CO 2 binding to silica than that obtained from DFT, an effect corrected by the optimized parameters given here. This behavior extends to other silica force fields (Clayff and Gulmen-Thompson), and the present Lennard-Jones parameters improve their performance as well. The effects of improved Lennard-Jones parameters on the structural and dynamical properties of condensed CO 2 in silica slit pores are also examined.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Charge Transfer and Recombination Pathways through Fullerene Guests in Porphyrin-Based MOFs

Porphyrin-based metal–organic frameworks (MOFs) offer a unique platform for building porous donor–acceptor networks that exhibit long-lived charge separation and transport upon incorporation of electron acceptor guest species. Here, in this work, porphyrin-based MOFs, PCN-222(H 2 ) and PCN-222(Zn), synthesized as nanoparticle suspensions, are successfully infiltrated with fullerene acceptor molecules, C 60 or PC 61 BM, in both polar and nonpolar solvent environments. The location and relative binding strength of these guest species are evaluated through a combination of N 2 physisorption measurements, photoluminescence quenching, and UV–vis absorption titration experiments. Semiempirical tight binding calculations are used to screen potential locations of the fullerene guest within the MOF pores, and hybrid density functional theory (DFT)-computed interaction energies confirm the energetically favorable positions. The fundamental photophysics of these donor–acceptor host–guest combinations are probed using ultrafast transient absorption spectroscopy. Sub-picosecond electron transfer involving initial exciplex population is observed, with slow charge recombination lifetimes on the order of τ ~1 ns for all systems in both dimethylformamide and 1,4-dioxane. Charge recombination occurs through population of fullerene and/or framework porphyrin triplet states depending on the porphyrin metalation status. The photophysics of the fullerene-loaded MOFs are discussed in the context of relevant porphyrin–fullerene donor–acceptor molecules to highlight the unique role of the framework environment in dictating photoinduced electron transfer and decay pathways.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Smooth Dispersion Is Physically Appropriate: Assessing and Amending the D4 Dispersion Model

The addition of dispersion corrections to density functionals is essential for accurate energy and geometry predictions. Among them, the D4 scheme is popular due to its low computational cost and high accuracy. However, due to its design, the D4 correction can occasionally lead to anomalies, such as unphysical curvature and bumps in the potential energy surface. We find these anomalies are common in the D4 model, although observable consequences are rarer than in the D3 model for reasons we explain. Nevertheless, we uncover instances of unphysical local minima and stationary points with the D4 scheme and propose two solutions that yield smoother dispersion energy as a function of nuclear position. One is trivial to implement, based on a smoother reparametrization of Gaussian weighting (D4S) to find the effective coordination number. The other replaces Gaussian weighting with soft linear interpolation (D4SL). These new approaches usually remove artificial extremum points, while maintaining accuracy.

Tkachenko, Nikolay V↗

DOE OSTI2024

Ring Size Effects on the Structures of Sandwich Compounds with a Stoichiometry of C 12 H 12 M (M = Ti–Ni)

Ring size effects on geometries and electronic structures were investigated for the (C n H n )M(C m H m ) (n = 4, 5, or 6; m = 8, 7, or 6; m + n = 12; M = Ti–Ni) systems using density functional theory. The lowest-energy C 12 H 12 M structures for the early transition metals titanium, vanadium, and chromium are the experimentally known singlet (η 5 -C 5 H 5 )Ti(η 7 -C 7 H 7 ), doublet (η 5 -C 5 H 5 )V(η 7 -C 7 H 7 ), and singlet (η 6 -C 6 H 6 ) 2 Cr, respectively. The likewise experimentally known singlet (η 6 -C 6 H 6 ) 2 Ti, doublet (η 6 -C 6 H 6 ) 2 V, and singlet (η 5 -C 5 H 5 )Cr(η 7 -C 7 H 7 ) are the secondlowest- energy structures with only a small energy difference between the two vanadium structures. For the later transition metals, dibenzenemetal complexes are the lowest-energy C 12 H 12 M species with two fully bonded hexahapto benzene rings in the lowest-energy manganese and iron derivatives and one hexahapto and one dihapto benzene ring in the lowest-energy cobalt and nickel derivatives. The lowest-energy (C 5 H 5 )M(C 7 H 7 ) structures for the later transition metals iron, cobalt, and nickel have partially bonded nonplanar C 7 H 7 rings with one or two uncomplexed C=C bonds. The (C 4 H 4 )M(C 8 H 8 ) (M = Ti–Ni) structures with the metal sandwiched between four- and eight-membered rings were found to be much higher in energy than their (C 5 H 5 )M(C 7 H 7 ) and (C 6 H 6 ) 2 M isomers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Tuning Catalytic Reactivity via Wetting Control through Oxygen Vacancies: Ru Clusters on Anatase TiO 2 and CeO 2 Supports

The shape of supported metal particles regulates their catalytic reactivity and is determined by the degree of wetting between the metal particle and the support surface. Flattened particles that wet support surfaces were reported in various catalytic systems, particularly in the subnanometer size regime. Such consequential metal–support wetting phenomena are poorly understood, and methods to study them on powder catalysts under realistic conditions are lacking. Here, we investigate the size-dependent wetting behaviors of Ru particles on two reducible-oxide supports, anatase TiO 2 (TiO 2 -A) and CeO 2 , under reducing catalytic conditions. X-ray absorption spectroscopy (XAS), low-energy ion scattering (LEIS), and density functional theory (DFT) are combined to determine the shape of Ru particles. Ru particles remain three-dimensional without wetting the TiO 2 -A support within the coverage range studied (0.06–0.98 Ru nm –2 ). In contrast, at low coverages (<0.25 Ru nm –2 ), Ru wets the CeO 2 support to form flat, disordered structures. The higher wettability of CeO 2 than TiO 2 -A is attributed to oxygen vacancies in the near-surface region. The shape difference between small Ru particles or clusters on the two supports leads to drastically contrasting catalytic reactivities in polyolefin hydrogenolysis, despite similar diameters. As a result, this work highlights the implications of metal–support wetting, or cluster shape, on catalytic behaviors of small metal clusters, while establishing the foundation for future systematic studies of such a phenomenon in realistic systems, by delivering a multitechnique methodology and revealing governing fundamental principles.

42 ENGINEERING↗