SEARCH · Engineering Papers

Results for “equilibrium potential”

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 325 records · Page 18

All sources NASA NTRS DOE OSTI

NASA NTRS2019

Recent Advancements in Modeling and Simulation of Entry Systems at NASA

This paper describes recent development of modeling and simulation technologies for entry systems in support of NASA's exploration missions. Mission-tailored research and development in modeling of entry systems occurs across the Agency (e.g., within the Orion and Mars 2020 Programs), however the aim of this paper is to discuss the broad, cross-mission research conducted by NASA's Entry Systems Modeling (ESM) Project, which serves as the Agency's only concerted effort toward advancing entry systems across a range of technical disciplines. Technology development in ESM is organized and prioritized from a system-level perspective, resulting in four broad technical areas of investment: (1) Predictive material modeling, (2) Shock layer kinetics and radiation, (3) Computational and experimental aerosciences, and (4) Guidance, navigation, and control. Investments in thermal protection material modeling are geared toward high-fidelity, predictive models capable of handling complex structures, with an eye toward optimizing design performance and quantifying thermal protection system reliability. New computational tools have been developed to characterize material properties and behavior at the microstructural level, and experimental techniques (molecular beam scattering, micro-computed tomography, among others) have been developed to measure material kinetics, morphology, and other parameters needed to inform and validate detailed simulations. Advancements have also been made in macrostructural simulation capability to enable 3-D system-scale calculations of material response with complex topological features, including differential recession of tile gaps. Research and development in the area of shock layer kinetics has focused on air and CO2-based atmospheres. Capacity and capability of the NASA Ames Electric Arc Shock Tube (EAST) have been expanded in recent years and analysis of resulting data has led to several improvements in kinetic models, while simultaneously reducing uncertainties associated with radiative heat transfer predictions. First-principles calculations of fundamental kinetic, thermodynamic, and transport data, along with state-specific models for non-equilibrium flow regimes, have also yielded new insights and have the potential to vastly improve model fidelity. Aerosciences is a very broad area of interest in entry systems, yet a number of important challenges are being addressed: Coupled fluid-structure simulations of parachute inflation and dynamics; Experimental and computational studies of vehicle dynamics; Multi-phase flow with dust particles to simulate entry environments at Mars during dust storms; Studies of roughness-induced heating augmentation relevant to tiled and woven thermal protection systems; and Advanced numerical methods to optimize computational analyses for desired accuracy versus cost. Guidance and control in the context of entry systems has focused on development of methods for multi-axis control (i.e. pitch and yaw, rather than bank angle alone) of spacecraft during entry and descent. With precision landing requirements driven by Mars human exploration goals, recent efforts have yielded 6-DOF models of multi-axis control with propulsive descent of both inflatable and rigid ellipsled-like architectures.

Barnhardt, Michael D.↗

NASA NTRS2021

Recent Advancements in Modeling and Simulation of Entry Systems at NASA

This paper describes recent development of modeling and simulation technologies for entry systems in support of NASA’s exploration missions. Mission-tailored research and development in modeling of entry systems occurs across the Agency (e.g., within the Orion and Mars 2020 Programs), however the aim of this paper is to discuss the broad, cross-mission research conducted by NASA’s Entry Systems Modeling (ESM) Project, which serves as the Agency’s only concerted effort toward advancing entry systems across a range of technical disciplines. Technology development in ESM is organized and prioritized from a system-level perspective, resulting in four broad technical areas of investment: (1) Predictive material modeling, (2) Shock layer kinetics and radiation, (3) Computational and experimental aerosciences, and (4) Guidance, navigation, and control. Investments in thermal protection material modeling are geared toward high-fidelity, predictive models capable of handling complex structures, with an eye toward optimizing design performance and quantifying thermal protection system reliability. New computational tools have been developed to characterize material properties and behavior at the microstructural level, and experimental techniques (molecular beam scattering, micro-computed tomography, among others) have been developed to measure material kinetics, morphology, and other parameters needed to inform and validate detailed simulations. Advancements have also been made in macrostructural simulation capability to enable 3-D system-scale calculations of material response with complex topological features, including differential recession of tile gaps. Research and development in the area of shock layer kinetics has focused on air and CO2-based atmospheres. Capacity and capability of the NASA Ames Electric Arc Shock Tube (EAST) have been expanded in recent years and analysis of resulting data has led to several improvements in kinetic models, while simultaneously reducing uncertainties associated with radiative heat transfer predictions. First-principles calculations of fundamental kinetic, thermodynamic, and transport data, along with state-specificmodels for non-equilibrium flow regimes, have also yielded new insights and have the potential to vastly improve model fidelity. Aerosciences is a very broad area of interest in entry systems, yet a number of important challenges are being addressed: Coupled fluid-structure simulations of parachute inflation and dynamics; Experimental and computational studies of vehicle dynamics; Multi-phase flow with dust particles to simulate entry environments at Mars during dust storms; Studies of roughness-induced heating augmentation relevant to tiled and woven thermal protection systems; and Advanced numerical methods to optimize computational analyses for desired accuracy versus cost. Guidance and control in the context of entry systems has focused on development of methods for multi-axis control (i.e. pitch and yaw, rather than bank angle alone) of spacecraft during entry and descent. With precision landing requirements driven by Mars human exploration goals, recent efforts have yielded 6-DOF models of multi-axis control with propulsive descent of both inflatable and rigid ellipsled-like architectures. Results for both configurations have demonstrated the ability to land within the 50-meter precision requirement demanded by Mars human exploration missions, while also reducing propellant requirements by enabling more efficient control through entry and descent. Ongoing research in GN&C is developing mechanical specifications for the systems and establishing engineering feasibility.

aerothermodynamics↗

NASA NTRS1982

Nonlinear dynamo oscillations

The stability of steady equilibrium amplitudes of magnetic fields generated by convection flows is investigated. Two cases are considered in detail, for which steady solutions are available from the previous work of Busse (1973) and Busse (1977). In the first case instability occurs primarily because of magnetic flux expulsion at high magnetic Reynolds numbers. In the second case the change in the velocity field caused by the Lorentz force enhances dynamo action. This subcritical finite amplitude dynamo is potentially unstable. In typical cases the nonlinear dynamo oscillations that replace the steady equilibrium solutions are investigated by numerical integration.

Busse, F. H.↗

DOE OSTI2024

Out-of-equilibrium dynamics of a grid-like Fe( II ) spin crossover dimer triggered by a two-photon excitation

The application of two-photon excitation (TPE) in the study of light-responsive materials holds immense potential due to its deeper penetration and reduced photodamage. Despite these benefits, TPE has been underutilised in the investigation of the photoinduced spin crossover (SCO) phenomenon. Here, we employ TPE to delve into the out-of-equilibrium dynamics of a SCO Fe II dimer of the form [Fe II (HL) 2 ] 2 (BF 4 ) 4 ·2MeCN (HL = 3,5-bis{6-(2,2'-bipyridyl)}pyrazole). Optical transient absorption (OTA) spectroscopy in solution proves that the same dynamics take place under both one-photon excitation (OPE) and TPE. The results show the emergence of the photoinduced high spin state in less than 2 ps and with a lifetime of 147 ns. Time-resolved photocrystallography (TRXRD) reveals a single molecular reorganisation within the first 500 ps following TPE. Additionally, variable temperature single crystal X-ray diffraction (VTSCXRD) and magnetic susceptibility measurements confirm that the thermal transition is silenced by the solvent. While the results of the OTA and TRXRD utilising TPE are intriguing, the high pump fluencies required to excite enough metal centres to the high spin state may impair its practical application. Nonetheless, this study sheds light on the potential of TPE for the investigation of the out-of-equilibrium dynamics of SCO complexes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1981

Eruptive solar magnetic fields

The quasi-steady evolution of solar magnetic fields in response to gradual photospheric changes is considered, with particular attention given to the threshold of a sudden eruption in the solar atmosphere. The formal model of an evolving, force-free field dependent on two Cartesian coordinates is extended to a field which is not force free but in static equilibrium with plasma pressure and gravity. The basic physics is illustrated through the evolution of a loop-shaped electric current sheet enclosing a potential bipolar field with footpoints rooted in the photosphere. A free-boundary problem is posed and then solved for the equilibrium configuration of the current sheet in a hydrostatically supported isothermal atmosphere. As the footpoints move apart to spread a constant photospheric magnetic flux over a larger region, the equilibria available extend the field to increasing heights.

Low, B. C.↗

DOE OSTI2025

Coherent phonon flatband generated in GaAs/AlAs superlattices via layer-selective optical pumping

Flatbands, characterized by their dispersionless energy levels in electronic, magnetic, and phononic systems, hold substantial potential for advancements in electronics and quantum information processing. Most flatbands exist in thermal equilibrium and cannot be easily created or annihilated externally, limiting their flexibility as switchable knobs for use in microelectronics and quantum applications. In our work, we demonstrate the generation of a coherent phonon flatband in a GaAs/AlAs superlattice using 800 nm femtosecond laser pulses. This coherent phonon flatband does not correspond to a phonon eigenmode at equilibrium and exhibits strong coupling with two branches of coherently excited longitudinal phonon modes. With molecular dynamics simulations, we show more generally that the coherent phonon flatband can be induced by coherently and spatially modulated optical excitations of superlattice structures. Our results highlight a pathway for coherent phonon flatband creation in the time domain that can be generalized to various superlattice systems, potentially inspiring the realization of coherent flatband generation of other quasiparticles.

Condensed-matter physics↗

NASA NTRS2020

CTCVJ2056-3014: An X-Ray-Faint Intermediate Polar Harboring an Extremely Fast-Spinning White Dwarf

We report on XMM-Newton X-ray observations that reveal CTCV J2056-3014 to be an unusual accretion-powered, intermediate polar (IP) system. It is a member of the class of X-ray-faint IPs whose space density remains unconstrained but potentially very high, with L(sub χ,0.3-12 keV of 1.8 × 10(exp 31) erg s(exp -1). We discovered a coherent 29.6 s pulsation in X-rays that was also revealed in our reanalysis of published optical data, showing that the system harbors the fastest-spinning, securely known white dwarf (WD) so far. There is no substantial X-ray absorption in the system. Accretion occurs at a modest rate (∼6 × 10(exp -12)M(sub ⨀) yr(exp -1) in a tall shock above the WD, while the star seems to be spinning in equilibrium and to have low magnetic fields. Further studies of CTCV J2056-3014 potentially have broad implications on the origin of magnetic fields in WDs, on the population and evolution of magnetic cataclysmic variables, and also on the physics of matter around rapidly rotating magnetic WDs.

Cataclysmic variable stars↗

DOE OSTI2021

Computational Fermi Level Engineering and Doping-Type Conversion of Mg:Ga 2 O 3 via Three-Step Synthesis Process

Gallium oxide (Ga 2 O 3 ) is being actively explored for electronics that can operate at high power, temperature, and frequency as well as for deep-ultraviolet optoelectronics and other applications due to its ultra-wide bandgap (UWBG) and low projected fabrication cost of large-size and high-quality crystals. Efficient n-type doping of monoclinic beta-phase of Ga 2 O 3 has been achieved, but p-type doping faces fundamental obstacles due to compensation, deep acceptor levels, and the polaron transport mechanism of free holes. However, aside from the challenges of achieving p-type conductivity, plenty of opportunity exists to engineer the position of the Fermi level for improved design of Ga 2 O 3 -based devices. We use first-principles defect theory and defect equilibrium calculations to simulate a three-step growth-annealing-quench synthesis protocol for hydrogen-assisted Mg doping in ß-Ga 2 O 3 . The simulations take into account the gas phase equilibrium between H 2 , O 2 , and H 2 O, which determines the H chemical potential. We predict Ga 2 O 3 doping-type conversion to a net p-type regime after growth under reducing conditions in the presence of H2 followed by O-rich annealing, which is a similar process to Mg acceptor activation by H removal in GaN. For equilibrium annealing with re-equilibration of compensating O vacancies, there is an optimal temperature that maximizes the Ga 2 O 3 net acceptor density for a given Mg doping level; the acceptor density is further increased in the non-equilibrium annealing scenario without re-equilibration. After quenching to operating temperature, the Ga 2 O 3 Fermi level drops below mid-gap down to about 1.5 eV above the valence band maximum, creating a significant number of uncompensated neutral Mg Ga 0 acceptors. The resulting free hole concentration in Ga 2 O 3 is very low even at elevated operating temperature (~10 8 cm -3 at 400°C) due to the deep energy level of these Mg acceptors, and hole conductivity is further impeded by the polaron hopping mechanism. However, the Fermi-level reduction and suppression of free electron density in this doping-type converted (NA>ND) Ga 2 O 3 material are important for improved designs of Ga 2 O 3 electronic devices. These results illustrate the power of computational predictions not only for new materials but also for their synthesis science.

36 MATERIALS SCIENCE↗

NASA NTRS1997

The Encounter of P/Shoemaker-Levy 9 with the Jovian Plasma and Extended Sodium Cloud

The encounter of comet P/Shoemaker-Levy 9 with Jupiter during July, 1994, provided an unprecedented opportunity to observe any potential perturbations in the Jovian plasma torus and extended sodium cloud as the comet entered the planet's atmosphere. Though the most obvious affect of the encounter was the distinctive response of the visible disk to the impact of the cometary fragments, the potential disruptions to the extended Jovian atmosphere and the restoration of the system to equilibrium also provided a test for the current interpretation of the Jovian plasma torus and sodium magneto-nebula. The observations that were performed for this grant were made by a complementary group of researchers and could not have been made if the individuals worked singly. In a sense, the exciting opportunity provided by this astronomical event also provided a mechanism to test the potential of pooling limited resources from several sources to construct a state-of-the-art spectrally resolving instrument, to acquire the necessary time and resources from institutions that maintain world-class optical telescopes, to perform the observations with the assistance of students, and to analyze the data sets.

Niciejewski, R. J.↗

DOE OSTI2022

Global biomass supply modeling for long-run management of the climate system

Bioenergy is projected to have a prominent, valuable, and maybe essential, role in climate management. However, there is significant variation in projected bioenergy deployment results, as well as concerns about the potential environmental and social implications of supplying biomass. Bioenergy deployment projections are market equilibrium solutions from integrated modeling, yet little is known about the underlying modeling of the supply of biomass as a feedstock for energy use in these modeling frameworks. We undertake a novel diagnostic analysis with ten global models to elucidate, compare, and assess how biomass is supplied within the models used to inform long-run climate management. With experiments that isolate and reveal biomass supply modeling behavior and characteristics (costs, emissions, land use, market effects), we learn about biomass supply tendencies and differences. The insights provide a new level of modeling transparency and understanding of estimated global biomass supplies that informs evaluation of the potential for bioenergy in managing the climate and interpretation of integrated modeling. For each model, we characterize the potential distributions of global biomass supply across regions and feedstock types for increasing levels of quantity supplied, as well as some of the potential societal externalities of supplying biomass. We also evaluate the biomass supply implications of managing these externalities. Finally, we interpret biomass market results from integrated modeling in terms of our new understanding of biomass supply. Overall, we find little consensus between models on where biomass could be cost-effectively produced and the implications. We also reveal model specific biomass supply narratives, with results providing new insights into integrated modeling bioenergy outcomes and differences. The analysis finds that many integrated models are considering and managing emissions and land use externalities of supplying biomass and estimating that environmental and societal trade-offs in the form of land emissions, land conversion, and higher agricultural prices are cost-effective, and to some degree a reality of using biomass, to address climate change.

09 BIOMASS FUELS↗

DOE OSTI2024

Downstream Separation of Formic Acid with Anion-Exchange Resin from Electrocatalytic Carbon Dioxide (CO 2 ) Conversion: Adsorption, Kinetics, and Equilibrium Modeling

The electrocatalytic carbon dioxide reduction reaction (CO 2 RR) that produces liquid formic acid (FA) is one potential route to minimize the CO 2 concentration in the atmosphere by utilizing CO 2 as a feedstock. However, one of the main disadvantages of this route is the high energy demand required for purification through conventional distillation. This method of purification has been shown to result in a large carbon footprint in life-cycle analysis (LCA). Therefore, finding an alternative purification technique that can integrate with an electrochemical CO 2 RR process with the potential to lower the carbon footprint would be of great interest. Under the experimental conditions of the CO 2 RR, FA remains in its basic form as the formate anion in the catholyte. Anion exchange resins with their basic properties can separate the formate anion from the catholyte mixture as an alternative to energy intensive distillation. In this work, three anion exchange resins with increasing basicity order: Amberlite IRA-96 < Amberlite IRA-910 < Ambersep 900 are tested to separate FA in the presence of aqueous solution of potassium bicarbonate catholyte with varying levels of resin and FA. Kinetics and equilibrium studies data for the FA adsorption are interpreted using several kinetics and isotherm models. The kinetics data fit well with a pseudo-first-order model at high initial FA concentrations and a pseudo-second-order model at low initial FA concentrations. The experimental data can be best explained with the Freundlich isotherm model. The fitted results show that strongly basic Ambersep 900 is more effective at separating FA, with a maximum adsorption capacity of 336.7 mg/g compared to Amberlite IRA-96 (275.2 mg/g) and Amberlite IRA-910 (209.2 mg/g) in neutral water. However, weakly basic Amberlite IRA-96 (110.8 mg/g) can separate FA more efficiently from potassium bicarbonate catholyte than Amberlite IRA-910 (99.9 mg/g) and Ambersep 900 (40.0 mg/g). As a result, competitive adsorption occurs in the presence of potassium bicarbonate catholyte, but the formic acid adsorption capacity is still sufficient for adequate separation.

20 FOSSIL-FUELED POWER PLANTS↗

NASA NTRS2002

Steady State Transportation Cooling in Porous Media Under Local, Non-Thermal Equilibrium Fluid Flow

An analytical solution to the steady-state fluid temperature for 1-D (one dimensional) transpiration cooling has been derived. Transpiration cooling has potential use in the aerospace industry for protection against high heating environments for re-entry vehicles. Literature for analytical treatments of transpiration cooling has been largely confined to the assumption of thermal equilibrium between the porous matrix and fluid. In the present analysis, the fundamental fluid and matrix equations are coupled through a volumetric heat transfer coefficient and investigated in non-thermal equilibrium. The effects of varying the thermal conductivity of the solid matrix and the heat transfer coefficient are investigated. The results are also compared to existing experimental data.

Rodriquez, Alvaro Che↗

DOE OSTI2021

Assessment of Performance of Density Functionals for Predicting Potential Energy Curves in Hydrogen Storage Applications

We report the availability of accurate computational tools for modeling and simulation is vital to accelerate the discovery of materials capable of storing hydrogen (H 2 ) under given parameters of pressure swing and temperature. Previously, we compiled the H2Bind275 data set consisting of equilibrium geometries and assessed the performance of 55 density functionals over this data set (Veccham, S. P.; Head-Gordon, M. J. Chem. Theory Comput. 2020, 16, 4963-4982). As it is crucial for computational tools to accurately model the entire potential energy curve (PEC), in addition to the equilibrium geometry, we extended this data set with 389 new data points to include two compressed and three elongated geometries along 78 PECs for H 2 binding, forming the H2Bind78 × 7 data set. By assessing the performance of 55 density functionals on this significantly larger and more comprehensive H2Bind78 × 7 data set, we identified the best performing density functionals for H 2 binding applications: PBE0-DH, ωB97X-V, ωB97M-V, and DSD-PBEPBE-D3(BJ). The addition of Hartree-Fock exchange improves the performance of density functionals, albeit not uniformly throughout the PEC. We recommend the usage of ωB97X-V and ωB97M-V density functionals as they offer good performance for both geometries and energies. In addition, we also identified B97M-V and B97M-rV as the best semilocal density functionals for predicting H 2 binding energy at its equilibrium geometry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Mapping the kinetic evolution of metastable grain boundaries under non-equilibrium processing

The kinetic evolution of a multiplicity of metastable grain boundaries (GBs) under fast driving conditions are studied by atomistic modeling. Assisted with an enhanced statistical analysis, the energetic evolution of GBs over a broad metastability-temperature space is mapped out, wherein two distinct regimes—an ageing regime and a rejuvenating regime—are retrieved with high fidelity. By comparing the results under various conditions (e.g. random perturbations, isothermal annealing, and fast heating/cooling), it is shown that such ageing/rejuvenating mechanism map is universal, irrespective of the actual stimuli used to elicit the metastable GBs. Here, the ageing/rejuvenating phenomena are demonstrated to stem from the energy imbalance of uphill climbing and downhill dropping during sequential transitions in the system's potential energy landscape. Without the necessity of introducing free parameters, such model can reconcile experimentally measured hardness variation of nanocrystalline metals subjected to femto-second laser irradiation, and it therefore provides a novel perspective on achieving a plurality of interfacial states and facilitating previously inaccessible property regimes.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Part I: Predicting performance of Purolite A532E resins for remediation of comingled contaminants in groundwater

Ion exchange (IX) resins are used in pump-and-treat (P&T) facilities to remove soluble groundwater contaminants. However, natural anions present at concentrations orders of magnitude higher than contaminants can compete for IX sites and impact resin lifecycles. Here, the Hanford Site’s 200 West Area P&T facility (Washington State, USA) was selected as a case study because it currently uses two IX resins: Purolite® A532E (A532E) to remove pertechnetate (TcO 4 - ) and DOWEX 21K (DOWEX) to remove uranium from groundwater. Nitrate (NO 3 - ), sulfate (SO 4 2- ), chloride (Cl - ), and carbonate (CO 3 2- ) anions have been identified to potentially compete for A532E and DOWEX IX sites. Hanford-relevant anion groundwater concentrations were used to design a series of laboratory-scale batch experiments to evaluate the impact of competing anions on resin performance and potential kinetic effects. These data are then modeled to obtain Cl--normalized equilibrium exchange coefficients (K) needed to predict IX resin performance. The work is presented in two parts, with IX performance evaluated for A532E in Part I and DOWEX in Part II. Part I results demonstrate that TcO 4 - uptake is not impacted by NO 3 - , SO 4 2- , Cl - , CO 3 2- (as HCO 3 - ) and U(VI) carbonate anions, with K TcO4-/Cl- > 4,000, likely due to the high selectivity of A532E trihexylammonium sites for the large, weakly hydrated TcO 4 - anion. Other anion K values were K NO3-/Cl- = 20, K SO4--/Cl- = 0.2, K HCO3-/Cl- = 0.09, K U/Cl- = 370–1000. These K values provide conservative parameters for predicting A532E performance, and demonstrate that, under these test conditions, A532E will remove TcO 4 - from current and future influent streams to meet groundwater treatment objectives.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2025

STEPs-SOL, a Peptoid Force Field Parameterization to Include Solvent Effects

As peptoids (N-substituted glycines) continue to gain popularity as a class of biomimetic polymers, the importance and demand for accurate force fields in molecular simulations also grow. Building on the vacuum-optimized Systematic and Extensible Force Field for Peptoids (STEPs) force field, here we present STEPs-SOL, a novel peptoid force field parametrization that effectively incorporates solvent effects to enhance the accuracy of peptoid simulations. The development of STEPs-SOL is based on the need for precise electrostatic modeling achieved through solvent-specific partial charge optimization. Here, our systematic approach significantly improves agreement with experimental measurements, reducing the mean absolute error in cis/trans ratio predictions (ΔG c/t ) by an average of 38% across multiple peptoid residues and solvent environments. This improved parametrization addresses computational challenges associated with nonbonded energies while maintaining a workflow that relies on high-level quantum mechanical data rather than depending solely on limited experimental equilibrium properties. By evaluating the effects of conformational bias in restrained electrostatic potential (RESP) charge generation and examining their impact on peptoid conformations in various solvents, we enhance our understanding of peptoid structural dynamics while providing a more accurate modeling framework.

force field↗

DOE OSTI2020

Nonequilibrium Solvent Effects during Photodissociation in Liquids: Dynamical Energy Surfaces, Caging and Chemical Identity

In the gas phase, potential energy surfaces can be used to provide insight into the details of photochemical reaction dynamics. In solution, however, it is unclear what potential energy surfaces, if any, can be used to describe even simple chemical reactions such as the photodissociation of a diatomic solute. In this paper, we use mixed quantum/classical (MQC) molecular dynamics (MD) to study the photodissociation of $Na$ $^+_2$ in both liquid Ar and liquid tetrahydrofuran (THF). We examine both the gas-phase potential surfaces and potentials of mean force (PMF), which assume that the solvent remains at equilibrium with the solute throughout the photodissociation process and show that neither resemble a nonequilibrium dynamical energy surface that is generated by taking the time integral of work. For the photodissociation of $Na$ $^+_2$ in liquid Ar, the dynamical energy surface shows clear signatures of solvent caging, and the degree of caging is directly related to the mass of the solvent atoms. For $Na$ $^+_2$ in liquid THF, local specific interactions between the solute and solvent lead to changes in chemical identity that create a kinetic trap that effectively prevents the molecule from dissociating. Here, the results show that nonequilibrium effects play an important role even in simple solution-phase reactions, requiring the use of dynamical energy surface to understand such chemical events.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Temperature and density dependent pair potential for deuterium under shock

Large-scale classical molecular dynamics (CMD) simulations naturally include the microscopic physics necessary for atomistic modeling of shock release at the ablator-fuel interface in an inertial confinement fusion (ICF) capsule. Here, the multi-megabar shocks utilized in ICF experiments can drive the deuterium fuel from ambient to electron volt temperatures (T) and multi-fold compression. Modeling interatomic interactions over such an extreme range of conditions is challenging for empirical bond order potentials. We generate a pair potential for deuterium with explicit temperature and mass density dependence from ab initio density functional theory molecular dynamics using the iterative Boltzmann inversion method. This potential accurately reproduces the radial distribution functions and pressures from DFT in CMD equilibrium simulations across a wide range of thermodynamic conditions, yet fails to return the expected Hugoniot relations when used in direct CMD shock simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗