Search NASA⌕ Search

SEARCH · Search NASA

Results for “Unified Free Energy Dynamics”

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.

9 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

An interoperable implementation of collective‐variable based enhanced sampling methods in extended phase space within the OpenMM package

Collective variable (CV)-based enhanced sampling techniques are widely used today for accelerating barrier-crossing events in molecular simulations. A class of these methods, which includes temperature accelerated molecular dynamics (TAMD)/driven-adiabatic free energy dynamics (d-AFED), unified free energy dynamics (UFED), and temperature accelerated sliced sampling (TASS), uses an extended variable formalism to achieve quick exploration of conformational space. These techniques are powerful, as they enhance the sampling of a large number of CVs simultaneously compared to other techniques. Extended variables are kept at a much higher temperature than the physical temperature by ensuring adiabatic separation between the extended and physical subsystems and employing rigorous thermostatting. Here, in this work, we present a computational platform to perform extended phase space enhanced sampling simulations using the open-source molecular dynamics engine OpenMM. The implementation allows users to have interoperability of sampling techniques, as well as employ state-of-the-art thermostats and multiple time-stepping. This work also presents protocols for determining the critical parameters and procedures for reconstructing high-dimensional free energy surfaces. As a demonstration, we present simulation results on the high dimensional conformational landscapes of the alanine tripeptide in vacuo, tetra-N-methylglycine (tetra-sarcosine) peptoid in implicit solvent, and the Trp-cage mini protein in explicit water.

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

Probabilistic and maximum entropy modeling of chemical reaction systems: Characteristics and comparisons to mass action kinetic models

We demonstrate and characterize a first-principles approach to modeling the mass action dynamics of metabolism. Starting from a basic definition of entropy expressed as a multinomial probability density using Boltzmann probabilities with standard chemical potentials, we derive and compare the free energy dissipation and the entropy production rates. We express the relation between entropy production and the chemical master equation for modeling metabolism, which unifies chemical kinetics and chemical thermodynamics. Because prediction uncertainty with respect to parameter variability is frequently a concern with mass action models utilizing rate constants, we compare and contrast the maximum entropy model, which has its own set of rate parameters, to a population of standard mass action models in which the rate constants are randomly chosen. We show that a maximum entropy model is characterized by a high probability of free energy dissipation rate and likewise entropy production rate, relative to other models. We then characterize the variability of the maximum entropy model predictions with respect to uncertainties in parameters (standard free energies of formation) and with respect to ionic strengths typically found in a cell.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1996

Spray Combustion Modeling with VOF and Finite-Rate Chemistry

A spray atomization and combustion model is developed based on the volume-of-fluid (VOF) transport equation with finite-rate chemistry model. The gas-liquid interface mass, momentum and energy conservation laws are modeled by continuum surface force mechanisms. A new solution method is developed such that the present VOF model can be applied for all-speed range flows. The objectives of the present study are: (1) to develop and verify the fractional volume-of-fluid (VOF) cell partitioning approach into a predictor-corrector algorithm to deal with multiphase (gas-liquid) free surface flow problems; (2) to implement the developed unified algorithm in a general purpose computational fluid dynamics (CFD) code, Finite Difference Navier-Stokes (FDNS), with droplet dynamics and finite-rate chemistry models; and (3) to demonstrate the effectiveness of the present approach by simulating benchmark problems of jet breakup/spray atomization and combustion. Modeling multiphase fluid flows poses a significant challenge because a required boundary must be applied to a transient, irregular surface that is discontinuous, and the flow regimes considered can range from incompressible to highspeed compressible flows. The flow-process modeling is further complicated by surface tension, interfacial heat and mass transfer, spray formation and turbulence, and their interactions. The major contribution of the present method is to combine the novel feature of the Volume of Fluid (VOF) method and the Eulerian/Lagrangian method into a unified algorithm for efficient noniterative, time-accurate calculations of multiphase free surface flows valid at all speeds. The proposed method reformulated the VOF equation to strongly couple two distinct phases (liquid and gas), and tracks droplets on a Lagrangian frame when spray model is required, using a unified predictor-corrector technique to account for the non-linear linkages through the convective contributions of VOF. The discontinuities within the sharp interface will be modeled as a volume force to avoid stiffness. Formations of droplets, tracking of droplet dynamics and modeling of the droplet breakup/evaporation, are handled through the same unified predictor-corrector procedure. Thus the new algorithm is non-iterative and is flexible for general geometries with arbitrarily complex topology in free surfaces. The FDNS finite-difference Navier-Stokes code is employed as the baseline of the current development. Benchmark test cases of shear coaxial LOX/H2 liquid jet with atomization/combustion and impinging jet test cases are investigated in the present work. Preliminary data comparisons show good qualitative agreement between data and the present analysis. It is indicative from these results that the present method has great potential to become a general engineering design analysis and diagnostics tool for problems involving spray combustion.

Chen, Yen-Sen↗
DOE OSTI2025

Cation Disorder of ${\text{Mg}}_{\mathbf{2}}{\text{SiO}}_{\mathbf{4}}$ in Super‐Earth Mantles

Understanding the mineralogy of exoplanets is essential for unraveling their interior structures, dynamics, and evolution. For large super-Earths, the post-post spinel ${\text{Mg}}_{\mathbf{2}}{\text{SiO}}_{\mathbf{4}}$, one of the major mantle phases, may undergo the order-disorder transition (ODT) at high temperatures. However, the ODT phase boundary of ${\text{Mg}}_{\mathbf{2}}{\text{SiO}}_{\mathbf{4}}$ has not been rigorously constrained. Additionally, fundamental thermodynamic properties of the disordered ${\text{Mg}}_{\mathbf{2}}{\text{SiO}}_{\mathbf{4}}$ remain poorly investigated. Here, we develop a unified machine learning potential (MLP) for ${\text{Mg}}_{\mathbf{2}}{\text{SiO}}_{\mathbf{4}}$ of ab initio accuracy under super-Earth mantle conditions. With the efficient MLP, we extensively calculate the free energy of post-post spinel ${\text{Mg}}_{\mathbf{2}}{\text{SiO}}_{\mathbf{4}}$ via the thermodynamic integration method. The results are used to constrain the ODT phase boundary. Furthermore, we report the P-V-T equation of state and Grüneisen parameters for post-post spinel ${\text{Mg}}_{\mathbf{2}}{\text{SiO}}_{\mathbf{4}}$ across various degrees of disorder. These thermodynamic properties are further applied to update the adiabatic thermal profiles and the mass-radius relation of super-Earths.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Unravelling ultralow thermal conductivity in perovskite Cs2AgBiBr6: dominant wave-like phonon tunnelling and strong anharmonicity

Abstract Understanding the lattice dynamics and heat transport physics in the lead-free halide double perovskites remains an outstanding challenge due to their lattice dynamical instability and strong anharmonicity. In this work, we investigate the microscopic mechanisms of anharmonic lattice dynamics and thermal transport in lead-free halide double perovskite Cs 2 AgBiBr 6 from first principles. We combine self-consistent phonon calculations with bubble diagram correction and a unified theory of lattice thermal transport that considers both the particle-like phonon propagation and wave-like tunnelling of phonons. An ultra-low thermal conductivity at room temperature (~0.21 Wm −1 K −1 ) is predicted with weak temperature dependence( ~ T −0.34 ), in sharp contrast to the conventional ~T −1 dependence. Particularly, the vibrational properties of Cs 2 AgBiBr 6 are featured by strong anharmonicity and wave-like tunnelling of phonons. Anharmonic phonon renormalization from both the cubic and quartic anharmonicities are found essential in precisely predicting the phase transition temperature in Cs 2 AgBiBr 6 while the negative phonon energy shifts induced by cubic anharmonicity has a significant influence on particle-like phonon propagation. Further, the contribution of the wave-like tunnelling to the total thermal conductivity surpasses that of the particle-like propagation above around 310 K, indicating the breakdown of the phonon gas picture conventionally used in the Peierls-Boltzmann Transport Equation. Importantly, further including four-phonon scatterings is required in achieving the dominance of wave-like tunnelling, as compared to the dominant particle-like propagation channel when considering only three-phonon scatterings. Our work highlights the importance of lattice anharmonicity and wave-like tunnelling of phonons in the thermal transport in lead-free halide double perovskites.

Chemistry↗
NASA NTRS2002

Structure and Dynamics of the Quiet Solar Chromosphere

The grant supported research on the structure of the quiet, nonmagnetic chromosphere and on wave excitation and propagation in both the nonmagnetic chromosphere and the magnetic network. The work on the structure of the chromosphere culminated in the recognition that between two competing views of the solar chromosphere, older models by Avrett and collaborators (referred to as VAL) and the newer, dynamical model by Carlsson & Stein (referred to as CS), the clear decision is in favor of the older models, and this in spite of the evident lack of physics, which does not include wave motion and oscillations. The contrast between the static VAL models and the dynamical CS model can be stated most succinctly by comparing the temperature variation implied by the VAL models and the temperature fluctuations of the CS model, which are, respectively, of the order of 10% for the VAL model (at heights where hydrogen is 50% ionized) and a factor of 10 (at the upper boundary of their chromospheric model). The huge fluctuations of the CS model have never been observed, whereas the smaller temperature variations of the VAL models are consistent with ground-based and space-based observations. While it should be obvious which model describes the Sun and which one fails, the case is far from settled in the minds of solar physicists. Thus, much educational work remains to be done and, of course, more research to develop arguments that make the case more convincing. The research on waves and oscillations has been based on a unified theory of excitation of acoustic waves in the field-free atmosphere and of transverse and longitudinal waves in magnetic flux tubes located in the magnetic network by noting, first, that impulsive excitation of all these waves in gravitationally stratified media leads to oscillations at the respective cutoff frequencies and, second, that the observed oscillation frequencies in the nonmagnetic and magnetic parts of the chromosphere match corresponding cutoff frequencies in the upper solar photosphere. The dynamical simulations by Carlsson and Stein have been most instructive and of fundamental importance for understanding wave propagation in a stratified medium by their "flaws", the most important of which is an intensity excess at the H2v emission peak in the H line of Ca II, which is surprising since the observed intensity should have been an upper limit to the simulated intensity. The only plausible explanation for a predicted intensity that is higher than observed is that energy is spread horizontally in upward propagation in the Sun, but not in the plane-wave modeling of CS (as well as by almost everyone else). Investigations of the horizontal size of the region disturbed by the upward-propagating shock in the acoustic-wave propagation implies that the waves in H2v bright-point oscillations emanate from a point source with a diameter corresponding to the width of an intergranular lane, about 100 km, and reach a size of about 4000 km in the upper layers of the chromosphere. Linear, analytic modeling of waves emanating from a point source in a stratified atmosphere shows that the upward-expanding propagation channel does not have sharp boundaries and that the shape of the wave front depends on the order of the wave behind the initial pulse. Otherwise, the behavior of the linear waves resembles that of the nonlinear shock waves observed in the Sun. Research that needs to be done to firm up the conclusions reached above concerns the numerical simulations of nonlinear waves and oscillations in a three-dimensional stratified atmosphere with impulsive excitation, and observations linking directly the horizontal size of the disturbed area in upward propagation to individual waves.

Kalkofen, Wolfgang↗
DOE OSTI2025

Entanglement, trace anomaly, and confinement in QCD

We formulate confinement in quantum chromodynamics (QCD) as an entropic surface phenomenon. Quark and gluon quantum information is localized on a transverse, entangling two-sphere of radius 𝑅 𝐸⁢𝐸 ; at this radius the QCD vacuum—partitioned by a hadron into interior and exterior regions—reaches its maximal entanglement entropy. Lattice-QCD determinations of the scalar (trace) gravitational form factors fix both 𝑅 𝐸⁢𝐸 and the transverse trace-anomaly density 𝜌 ℎ ⁡(𝑅 𝐸⁢𝐸 ), yielding a parameter-free slope 𝑐 ℎ = 8⁢𝜋 2 ⁢𝑅𝑆$^2_{𝐸⁢𝐸}$𝜌 ℎ ⁡(𝑅 𝐸⁢𝐸 ) and a mechanical entropy 𝑆 𝐸⁢𝐸 ⁡(𝑦) = 𝑐 ℎ ⁢𝑦 that grows linearly with rapidity 𝑦. The entropy gradient ∂ 𝑅 𝑆 𝐸⁢𝐸 changes sign at 𝑅 𝐸⁢𝐸 : it pushes colored degrees of freedom outward for 𝑟 <𝑅 𝐸⁢𝐸 and pulls them inward for 𝑟 >𝑅 𝐸⁢𝐸 , thereby localizing them on the codimension-2 entangling two-sphere Σ ⊥ = 𝑆$^2_{𝑅_{𝐸⁢𝐸}}$ (which, in the infinite-momentum frame (IMF), projects onto the transverse plane)—the “information wall.” This provides a high-energy (large-𝑦) entropic confinement diagnostic that complements—rather than replaces—Wilson’s area-law criterion, which probes long-distance dynamics near the rest frame (𝑦 → 0). Imposing unitarity on an entropic ansatz for the amplitude yields 𝜎⁡(𝑠) ∝ 𝑦 𝛿 . World data favor 𝛿 = 2 for elastic 𝑝⁡𝑝⁡($𝑝\bar{⁡𝑝}$) scattering and heavy-quark photoproduction, whereas 𝜙 photoproduction favors a softer 𝛿 = 0.387. All extracted cross sections remain well below the Froissart-Martin bound. These results provide a confinement criterion quantified directly from nonperturbative QCD inputs, unifying the trace anomaly, entanglement entropy, and high-energy scattering within a single quantitative framework.

Color confinement↗
DOE OSTI2025

Lewis Acid Site Engineering in Chromite Spinels Orchestrated Surface Reconstruction and Surpasses RuO 2 in Oxygen Evolution

Atomic-scale engineering of chromite spinels featuring redox-active tetrahedral A-sites and strong Cr–O covalency offers a promising route to superior platinum-group-metal-free oxygen evolution reaction (OER) catalysts. However, comprehensive studies addressing how cation substitution influences surface chemistry and governs OER activity and durability in chromite spinels remain limited. Here, in this work, a systematic investigation of the multicationic chromite series Ni x Fe y Cr 3−x−y O 4 is presented, identifying composition-dependent Lewis acidity as a descriptor of superior OER performance. It is further demonstrated that tuning surface acidity directly controls dynamic reconstruction processes and lattice-oxygen participation during spinel-based electrocatalysis. Following activation, the optimized Ni 0.8 Fe 0.3 Cr 1.9 O 4 catalyst delivers a current density of 10 mA cm −2 at an overpotential of 235 mV, surpassing RuO 2 , with excellent long-term stability. Integrating microscopic and spectroscopic analysis with operando impedance spectroscopy, it shows that activation generates an oxyhydroxide overlayer and reveals a previously unrecognized link between surface Lewis acidity and the growth kinetics and activity of these shells. Density functional theory calculations indicate that Fe incorporation at octahedral sites raises the O 2p-band center and lowers oxygen-vacancy formation energy, promoting lattice-oxygen activation and triggering reconstruction, yielding enhanced OER. This work integrates cation-driven surface-acidity modulation, acidity-governed reconstruction, and OER activity enhancement into a unified predictive framework for designing earth-abundant spinel-based catalysts.

operando impedance spectroscopy↗
DOE OSTI2025

Unified Universal Control and Coordination of Inverter-Based Resources, and Validation for a PV + Battery Hybrid Plant

As renewable energy deployment grows, hybrid power plants (HPPs) combining photovoltaic (PV) and battery systems must evolve to offer both energy and grid stability services. These systems typically include a mix of grid-following (GFL) and grid-forming (GFM) inverters, presenting unique coordination and control challenges. This Department of Energy–funded project developed and validated a Unified Universal Control and Coordination (UUCC) framework for such PV + battery hybrid plants, enabling seamless and stable operation, including ultrafast black start, autonomous synchronization, and robust frequency and voltage regulation, under different grid conditions. The project significantly advanced the understanding of inverter-based resource (IBR) control by developing and validating three complementary system-level approaches for hybrid GFL/GFM operation: 1. A combined Virtual Resistance (VR)-based GFL and Virtual Oscillator Control (VOC)-based GFM method, where each inverter type is governed by a specialized control strategy. Together, these achieve stable, fast-response coordination, eliminating inrush current and enabling smooth black start and grid synchronization across a wide range of grid strengths. 2. A Deadbeat-based UUCC strategy, which uses discrete-time, switching-cycle-level control for both GFL and GFM inverters. This approach replaces traditional PI/PLL control with a control parameter-free, high-bandwidth framework that supports stable LVRT and instantaneous synchronization under all conditions. 3. A benchmark comparison with Siemens’ commercial GFM microgrid controller, which provided a fast baseline platform. The commercial approach decoupled v & f control was implemented on a commercial microgrid controller.The baseline commercial benchmark helped highlight superior transient response and black start performance offered by the deadbeat and VOC approaches. These technical contributions offer substantial improvements over conventional inverter control schemes, which often rely on slow phase-locked loop (PLL)-based synchronization, require careful control parameters tuning, and prone to unstable in weak grids with GFL inverters and in stiff grid with GFM inverters therefore challenging for hybrid GFL+GFM under all grid conditions. The deadbeat-based UUCC framework enables simpler, faster, and more robust operation of hybrid IBR systems using wide-bandgap (WBG) devices such as SiC power semiconductors. The rapid expansion of hybrid distributed energy resources (DERs), including residential and commercial PV-BESS installations such as Tesla Powerwall, PV with vehicle-to-grid (V2G) capability, and other integrated configurations, presents complex operational challenges for medium-voltage radial distribution feeders. These networks are subject to frequent disturbances such as faults, switching operations, rapid reclosing sequences, and feeder reconfigurations, all of which introduce dynamic stress on IBRs. In addition, planned feeder segmentation and deliberate islanding for resilience will require DERs that can autonomously perform blackstart, establish voltage and frequency references, and resynchronize with the main grid. The advanced deadbeat-based UUCC control and blackstart functionalities developed in this project directly address these requirements, enabling decentralized and autonomous operation of inverter-dominated DERs in distribution systems under a wide range of fault and reconfiguration scenarios. From a public benefit perspective, these innovations enable more reliable and cost-effective integration of renewable energy into distribution networks. The ability to autonomously black start and stabilize grids under varying grid conditions support accelerates recovery from outages and support decentralized resilient energy systems. By reducing system complexity and improving performance, this project lays critical groundwork for future inverter-dominated power grids that are clean, reliable, and accessible to all.

14 SOLAR ENERGY↗