Search NASA⌕ Search

SEARCH · Search NASA

Results for “mixing state”

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

Redox conditions correlated with vibronic coupling modulate quantum beats in photosynthetic pigment–protein complexes

Significance Photosynthetic organisms evolved their light-harvesting antenna complexes to optimize energy transfer. It was recently shown that the redox environment can tune the mixing of electronic and vibrational states to steer energy through different pathways of a pigment–protein complex. Quantum beating signals in the spectra of pigment–protein complexes have been used to probe the excited-state dynamics within the complexes, but the microscopic dynamics that generate these signals and their role in promoting energy transfer are not fully understood. Here, we show that the redox environment that tunes energy transfer similarly tunes the quantum beating signals in the same complex. We find that the beats report on excited-state vibrations that maintain coherence through the vibronically enhanced energy transfer process.

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

Assessing the sequestration time scales of some ocean-based carbon dioxide reduction strategies

Abstract Ocean-based carbon dioxide (CO 2 ) removal (CDR) strategies are an important part of the portfolio of approaches needed to achieve negative greenhouse gas emissions. Many ocean-based CDR strategies rely on injecting CO 2 or organic carbon (that will eventually become CO 2 ) into the ocean interior, or enhancing the ocean’s biological pump. These approaches will not result in permanent sequestration, because ocean currents will eventually return the injected CO 2 back to the surface, where it will be brought into equilibrium with the atmosphere. Here, a model of steady state global ocean circulation and mixing is used to assess the time scales over which CO 2 injected in the ocean interior remains sequestered from the atmosphere. There will be a distribution of sequestration times for any single discharge location due to the infinite number of pathways connecting a location at depth with the sea surface. The resulting probability distribution is highly skewed with a long tail of very long transit times, making mean sequestration times much longer than typical time scales. Deeper discharge locations will sequester purposefully injected CO 2 much longer than shallower ones and median sequestration times are typically decades to centuries, and approach 1000 years in the deep North Pacific. Large differences in sequestration times occur both within and between the major ocean basins, with the Pacific and Indian basins generally having longer sequestration times than the Atlantic and Southern Oceans. Assessments made over a 50 year time horizon illustrates that most of the injected carbon will be retained for injection depths greater than 1000 m, with several geographic exceptions such as the Western North Atlantic. Ocean CDR strategies that increase upper ocean ecosystem productivity with the goal of exporting more carbon to depth will have mainly a short-term influence on atmospheric CO 2 levels because ∼70% will be transported back to the surface ocean within 50 years. The results presented here will help plan appropriate ocean CDR strategies that can help limit climate damage caused by fossil fuel CO 2 emissions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Structured interactions drive abrupt transitions in the spatial organization of microbial communities

Bacteria possess diverse mechanisms to regulate their motility in response to environmental and physiological signals, enabling them to navigate complex habitats and adapt their behavior. Some of these mechanisms are species specific and enable cells to modulate their movement based on the ecological identity of neighboring species. Here, we introduce a model in which bacteria interact via local signals that either enhance or suppress the motility of neighboring cells depending on species type. Through large-scale simulations and a coarse-grained stochastic model, we demonstrate the emergence of a sharp transition driven by nucleation processes: increasing the density of motility-suppressing interactions drives the system from a fully mixed, motile phase to a state characterized by large, stationary bacterial clusters. Remarkably, in systems with a large number of interacting species, this transition can be triggered solely by altering the structure of the motility-regulation interaction matrix while maintaining species and interaction densities constant. In particular, we find that heterogeneous and modular interactions promote the transition more readily than homogeneous random ones. These findings add a dimension to the theory of motility-induced phase separation and contribute to the ongoing effort to understand microbial interactions, suggesting that structured, nonrandom ones may be key to reproducing commonly observed spatial patterns in microbial communities.

bacterial communities↗
DOE OSTI2020

Spin excitations in the frustrated triangular lattice antiferromagnet NaYbO 2

In this work, we present a neutron scattering-based study of magnetic excitations and magnetic order in NaYbO 2 under the application of an external magnetic field. The crystal electric field-split J = 7=2 multiplet structure is determined, revealing a mixed |m z $\rangle$ ground state doublet and is consistent with a recent report Ding et al. Our measurements further suggest signatures of exchange effects in the crystal field spectrum, manifested by a small splitting in energy of the transition into the first excited doublet. The field-dependence of the low-energy magnetic excitations across the transition from the quantum disordered ground state into the fluctuation driven ordered regime is analyzed. Signs of a first-order phase transition into a noncollinear ordered state are revealed at the upper-field phase boundary of the ordered regime, and higher order magnon scattering, suggestive of strong magnon-magnon interactions, is resolved within the previously reported up-up-down phase. Our results reveal a complex phase diagram of field-induced order and spin excitations within NaYbO 2 and demonstrate the dominant role of quantum fluctuations cross a broad range of fields within its interlayer frustrated triangular lattice.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Lattice vibrations and energy landscape of the isoelectronic semiconductor series CuBr, ZnSe, GaAs, and Ge: The special case of CuBr and its $\mathcal{d}$-level chemistry

Here we have examined the lattice vibrations and the energy landscape of the isoelectronic diamond and zincblende semiconductor series CuBr, ZnSe, GaAs, and Ge. Vibrations are found to be an increasing function of ionicity, with the cation sublattice always vibrating more strongly than the anion sublattice. These findings are consistent with density functional theory (DFT) calculations of the energy landscape and temperature-dependent molecular dynamics simulations of the atomic-position fluctuations. For CuBr, inclusion of the Cu 3d Hubbard U term is found necessary to stabilize the zincblende structure and to bring its vibrational amplitudes into agreement with experiment. In addition, vibrations are found to strongly affect the CuBr near-edge x-ray absorption fine structure that we have successfully modeled by including displacements in our theoretical Bethe-Salpeter equation calculations. Reverse Monte Carlo structural refinements using large atomic configurations to simultaneously fit x-ray absorption and x-ray total-scattering data support these conclusions, and they reveal strong Cu-Br first-neighbor correlations and asymmetric distributions of interatomic distances in the temperature ranges of both negative and positive thermal expansion. Delineation of the CuBr valence band photoelectron spectrum into its Cu 3d and Br 4p states uniquely reveals their covalent mixing and further supports the DFT results.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Spin-Dependent $\pi$$\pi$* Gap in Graphene on a Magnetic Substrate

Here, we present a detailed analysis of the electronic properties of graphene/Eu/Ni(111). By using angle- and spin-resolved photoemission spectroscopy and ab initio calculations, we show that the intercalation of Eu in the graphene/Ni(111) interface gives rise to a gapped freestanding dispersion of the $\pi$$\pi$* Dirac cones at the K point with an additional lifting of the spin degeneracy due to the mixing of graphene and Eu states. The interaction with the magnetic substrate results in a large spin-dependent gap in the Dirac cones with a topological nature characterized by a large Berry curvature and a spin-polarized Van Hove singularity, whose closeness to the Fermi level gives rise to a polaronic band.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Throughput-Oriented and Accuracy-Aware DNN Training with BFloat16 on GPU

Deep Neural Networks (DNNs) have transformed the field of artificial intelligence and achieved extraordinary success in many areas. The training of DNNs is commonly compute and memory-intensive, which has resulted in several optimizations in the training phase. Among them, reduced precision is a typical and widely used technique to accelerate DNN training and reduce memory requirements. However, applying a widely adopted reduced precision format such as Float16 to all involved operations in DNN training is not optimal as the use of Float16 in some operations can hurt model accuracy. Meanwhile, additional optimizations including loss scaling and autocast techniques can mitigate the accuracy loss but lead to inherent overhead and inadequate use of reduced precision. In this work, we leverage another reduced precision format, BFloat16, and introduce a throughput-oriented and accuracy-aware approach to maximize the performance potential of DNN training. Since the high throughput provided by BFloat16 format is accompanied by low precision of the floating-point representation, this approach achieves high throughput by using BFloat16 on all DNN operations and avoids the accuracy loss through a customized accuracy-aware normalization. Results show that our approach outperforms the state-of-the-art mixed-precision training by 1.21x on an NVIDIA A100 GPU.

Xie, Zhen↗
DOE OSTI2021

Measurements of Density of Liquid Oxides with an Aero-Acoustic Levitator

Densities of liquid oxide melts with melting temperatures above 2000 °C are required to establish mixing models in the liquid state for thermodynamic modeling and advanced additive manufacturing and laser welding of ceramics. Accurate measurements of molten rare earth oxide density were recently reported from experiments with an electrostatic levitator on board the International Space Station. In this work, we present an approach to terrestrial measurements of density and thermal expansion of liquid oxides from high-speed videography using an aero-acoustic levitator with laser heating and machine vision algorithms. The following density values for liquid oxides at melting temperature were obtained: Y 2 O 3 4.6 ± 0.15; Yb 2 O 3 8.4 ± 0.2; Zr 0.9 Y 0.1 O 1.95 4.7 ± 0.2; Zr 0.95 Y 0.05 O 1.975 4.9 ± 0.2; HfO 2 8.2 ± 0.3 g/cm 3 . The accuracy of density and thermal expansion measurements can be improved by employing backlight illumination, spectropyrometry and a multi-emitter acoustic levitator.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Igiugig Site Visit Report

The National Renewable Energy Laboratory (NREL) team conducted a site visit during January 23-25, 2019 with the Igiugig Village Council (IVC) and other stakeholders to assist the community of Igiugig in refining their long-term energy strategy. The agenda included a tour of the village, data collection, presentations, long-term visioning exercises, identification of energy scenarios, and a school presentation. Participants in the site visit activities spanned a range of local, regional, government, and industry participation. There were 24 attendees at the all-day community workshop on Thursday, January 24th, with participants from the IVC, NREL, Bristol Bay Native Association (BBNA), Bristol Bay Native Corporation (BBNC), Lake and Peninsula Borough, Alaska Energy Authority, Intergrid, ORPC, Deer Stone Consulting, the Southwest Alaska Municipal Conference, and University of Alaska Fairbanks - Alaska Center for Energy and Power. Igiugig is a well-organized rural Alaskan community, with a leadership team that appears to have broad community support. The use of consensus-based decision making is one tangible example of how the leadership team actively invokes and promotes community-focused thinking. The NREL team observed anecdotal evidence of how this approach seems to be strengthening the community: active participation in the visioning exercise, a student hosted fund-raiser dinner at the school, a clean and organized landfill, and well-maintained roads and buildings. Igiugig has wind, solar, and river hydrokinetic resources readily available within the community. Wind has shown to be a promising resource in the region, and has been integrated into microgrids around the state, but the community has had mixed success with wind technologies: some devices failed quickly and others continue to operate. The economics of solar energy are improving in Alaska, and economical solar projects are being installed around the state. Considering that economic activity in Igiugig peaks during the summer sport fishing season, solar could prove to be a valuable supplement to the electrical system. Igiugig has been a test site for two different river hydrokinetic devices, and they have began a third project to operate ORPC's RivGen device, which is delivering valuable device performance data, operations and maintenance experience, and design refinement information. Already, this device has made over 7-million revolutions, and delivered over 8 MWh of power to the community. Igiugig's river resource is fairly unique because it is available year-round and has the potential to provide reliable base-load power for months at a time. A preliminary investigation of this diverse resource mix suggests that Igiugig could achieve very high levels (70% or more) of annual renewable generation contributions. A critical step in pursuing this path is identifying the mix of energy assets (generation and storage) that best meets the community's budget, needs, and goals. The assets that a community installs early in their grid-modernization initiative can constrain the options that are economical at later stages, which may lead to sub-optimal solutions. This is where technical and economic analysis of potential scenarios (i.e., distinct mixes of energy assets) can be useful in identifying the most promising pathways so that a community can make informed and strategic decisions about the assets they install. These analyses are most accurate and informative when they are based on actual technology performance and cost data. As the RivGen project continues to operate and generate this data, we will be better prepared to evaluate the technology's long-term viability and to identify research areas that would improve it. Igiugig's energy projects are at the cutting edge of two intersecting technology areas: 1) deploying an operational river hydrokinetic turbine, and 2) integrating renewable energy sources to achieve very high percent renewables. If these projects are successful, the lessons learned and technologies developed could be valuable for other microgrids around the world. Igiugig's unique resource mix make it an ideal location for this work, and the community's organizational strength make it an ideal partner in this ground-breaking work. Ongoing support for technical assistance to manage and address technical challenges along the way will maximize the probability of project success.

13 HYDRO ENERGY↗
DOE OSTI2025

Modeling of surface-state induced inter-electrode isolation of n -on- p devices in mixed-field and γ -irradiation environments

Position sensitive n-on-p silicon sensors will be utilized in the tracker and in the High Granularity Calorimeter (HGCAL) of the Compact Muon Solenoid (CMS) experiment at High Luminosity Large Hadron Collider (HL-LHC). The detrimental effect of the radiation-induced accumulation of positive net oxide charge on position resolution in n-on-p sensors has typically been countered by the application of isolation implants like p-stop or p-spray between n + -electrodes. In addition to the positively charged layer inside the oxide and close to the Si/SiO 2 -interface, surface damage introduced by ionizing radiation in SiO 2 -passivated silicon particle detectors includes the accumulation of trapped-oxide-charge and interface traps. A previous study of either n/γ (mixed field)- or γ-irradiated Metal-Oxide-Semiconductor (MOS) capacitors showed evidence of substantially higher introduction rates of acceptor- and donor-type deep interface traps (N it,acc/don ) in mixed-field environment. Here, in this work, an inter-pad and -strip resistance (or resistivity (ρ int )) simulation study of n-on-p sensors with and without p-stop isolation implants was conducted for both irradiation types. Higher levels of ρint showed correlation to higher densities of deep N it,acc/don , with the inter-pad isolation performance of the mixed-field irradiated sensors becoming independent of the presence of p-stop implant between the n + - electrodes up to about 100 kGy. The low introduction rates of deep N it,acc/don in γ-irradiated sensors resulted in high sensitivity of ρ int to the presence and peak doping of p-stop above the lowest dose of about 7 kGy in the study. As a consequence of the advantageous influence of radiation-induced accumulation of deep N it on the inter-electrode isolation, position sensitive n-on-p sensors without isolation implants may be considered for future HEP-experiments where the radiation is largely due to hadrons.

Akchurin, N. [Texas Tech Univ., Lubbock, TX (Unite↗
DOE OSTI2023

Simulation of Methane Oxycombustion in Supercritical Carbon Dioxide

The Allam-Fetvedt Cycle (AFC) offers the potential for electricity generation with minimal or even negative carbon emissions by using oxycombustion in supercritical carbon dioxide (sCO2). Because of the extreme high temperatures and pressures encountered in these systems, simulation of oxycombustion requires use of a complex equation of state such as the modified Soave-Redlich-Kwong (SRK) equation, coupled with equally complex correlations for transport properties. The current work simulates a model combustor using both an ideal gas equation of state and the SRK equation of state. The model combustor features a jet of CH4 mixing with and reacting with O2 in the sCO2 ambient. Two jet configurations are considered: one of pure CH4, and one of CH4 mixed with sCO2. Comparison of results shows that key output properties such as combustor temperature and flow distribution are impacted by the choice of equations of state. Inspection of simulation results shows that mixing and reaction rates, as well as the uniformity of the flow, are impacted by the choice of equation of state. These results show that even though use of a complex equation of state greatly increases the computational costs of these simulations, it is critical to accurately capturing the physics of these systems.

computational fluid dynamics↗
DOE OSTI2025

The Two Arctic Wintertime Boundary Layer States: Disentangling the Role of Cloud and Wind Regimes in Reanalysis and Observations During MOSAiC

The wintertime central Arctic atmosphere comprises a radiatively clear and a radiatively opaque state, which are linked to synoptic forcing and mixed-phase clouds. Weather and climate models often lack process representations surrounding these states, but prior work mostly treated the problem as an aggregate of synoptic conditions, resulting in partially overlapping biases. Here, we disaggregate the Arctic states and confront ERA5 reanalysis with observations from the MOSAiC campaign over the central Arctic sea ice during winter 2019/2020. Low-level winds and liquid water path (LWP) are combined to derive different synoptic classes. Results show that the clear state is primarily formed by weak/moderate winds and the absence of liquid-bearing clouds, while strong winds and enhanced LWP primarily form the radiatively opaque state. ERA5 struggles to reproduce these basic statistics, shows too weak sensitivity of thermal radiation to synoptic forcing, and overestimates thermal radiation for similar LWP amounts. The latter is caused by a warm bias, which has a pronounced inversion structure and is largest in clear and calm conditions. Under strong synoptic forcing, the warm bias is constant with height and discrepancies in mixed-phase cloud altitude appear. Separating synoptic conditions is regarded as useful for process-oriented evaluation of the Arctic troposphere in models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Operando Characterization of Fe in Doped Ni x (Fe 1– x )O y H z Catalysts for Electrochemical Oxygen Evolution

Iron-doped nickel oxyhydroxides, Ni x (Fe 1– x )O y H z , are among the most promising oxygen evolution reaction (OER) electrocatalysts in alkaline environments. Although iron (Fe) significantly enhances the catalytic activity, there is still no clear consensus on whether Fe directly participates in the reaction or merely acts as a promoter. To elucidate the Fe’s role, we performed operando X-ray spectroscopy studies supported by DFT on Ni x (Fe 1– x )O y H z electrocatalysts. We probed the reversible changes in the structure and electronic character of Ni x (Fe 1– x )O y H z as the electrode potential is cycled between the resting (here at 1.10 V RHE ) and operational states (1.66 V RHE ). DFT calculations and XAS simulations on a library of Fe structures in various NiO y H z environments are in favor of a distorted local octahedral Fe(III)O 3 (OH) 3 configuration at the resting state with the NiO y H z scaffold going from α-Ni(OH) 2 to γ-NiOOH as the potential is increased. Under catalytic conditions, EXAFS and HERFD spectra reveal changes in p-d mixing (covalency) relative to the resting state between O/OH ligands and Fe leading to a shift from octahedral to square pyramidal coordination at the Fe site. XES measurements and theoretical simulations further support that the Fe equilibrium structure remains in a formal Fe(III) state under both resting and operational conditions. These spectral changes are attributed to potential dependent structural rearrangements around Fe. The results suggest that ligand dissociation leads to the C 4v symmetry as the most stable intermediate of the Fe during OER. This implies that Fe has a weakly coordinated or easily dissociable ligand that could serve to coordinate the O–O bond formation and, tentatively, play an active role in the Ni x (Fe 1– x )O y H z electrocatalyst.

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

Twisting the Hubbard model into the momentum-mixing Hatsugai–Kohmoto model

The Hubbard model is a standard theoretical tool for studying materials with strong electron–electron interactions, such as cuprate superconductors. Unfortunately, interaction-driven phenomena, such as a transition into the strongly correlated Mott insulator phase, are difficult to treat with established theoretical techniques. However, the exactly solvable Hatsugai–Kohmoto model displays similar Mott physics. In this work we show how the Hatsugai–Kohmoto model can be deformed continuously into the Hubbard model. The trick is to systematically reintroduce all the momentum mixing that the original Hatsugai–Kohmoto model omits. This can be accomplished by grouping n momenta into a cell and hybridizing them, resulting in the momentum-mixing Hatsugai–Kohmoto model. We recover the Bethe ansatz ground-state energy of the one-dimensional Hubbard model to within 1% from only ten mixed momenta. Overall, the convergence scales as 1/n2 as opposed to the inverse linear behaviour of standard finite-cluster techniques. Our results for a square lattice reproduce all the known features from state-of-the-art simulations also with only a few mixed momenta. Consequently, we believe that the momentum-mixing Hatsugai–Kohmoto model offers an alternative tool for strongly correlated quantum matter.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Upcycling Mixed Spent Ni-Lean Cathodes into Ni-Rich Polycrystalline Cathodes

Sustainable battery recycling is vital for conserving resources and reducing environmental impacts. Current open- and closed-loop recycling strategies often focus on recovering individual components, making the reuse of mixed cathode materials a complex challenge. Meanwhile, the research on upcycling has been limited to using pristine cathode feedstocks and virgin materials for synthesis. Here, to address this issue, we present an upcycling approach for spent Ni-lean mixed cathode materials that integrate an upcycling hydrometallurgical recycling process with traditional hydrometallurgical methods. This strategy achieves a utilization of 92.31 mol % of recycled materials, enabling the regeneration of Ni-rich cathode materials while significantly reducing the reliance on virgin resources. The regenerated 83Ni cathode materials demonstrate physical properties comparable to those produced from virgin materials. Electrochemical evaluations using single-layer pouch cells show that both recycled and virgin cathodes exhibit initial specific capacity close to 201.1 mAh/g and maintain approximately 88 % capacity retention after 500 cycles. Additionally, 2Ah cells confirmed these findings, delivering 85 % capacity retention after about 900 cycles. Techno-economic analysis demonstrates notable environmental benefits, including reductions in greenhouse gas emissions and energy consumption, achieving 232.75 MJ/kg of product, which is 8.6 % lower than traditional methods and comparable to direct upcycling. Furthermore, the upcycling hydrometallurgical recycling process generates the highest profit, proving its economic viability. This scalable and versatile process is adaptable to varying transition metal compositions, facilitating a closed-loop recycling system that bridges mixed spent cathodes with next-generation cathode materials, and offers a sustainable solution for managing waste battery materials.

Hydrometallurgical recycling↗
DOE OSTI2021

Role of a 3D Structure in Energy Transfer in Mixed-Ligand Metal–Organic Frameworks

We present a detailed investigation of the photophysical properties of mixed-ligand pyrene- and porphyrin-based metal–organic frameworks (MOFs) as a function of their 3D structure. Solvothermal reactions between metal salts (InCl 3 , Zr(acac) 4 , and ZrCl 4 ) and suitable ratios of 1,3,6,8-tetrakis(p-benzoic acid)pyrene (TBAPy) and meso-tetrakis- (4-carboxyphenyl)porphyrin (TCPP) were performed to prepare a series of mixed-ligand ROD-7, NU-901, and NU-1000 MOFs. Time-resolved and steady-state fluorescence measurements were conducted on the mixed-ligand MOFs to study their photophysics. Based on the results, we concluded that upon excitation of TBAPy linkers in the MOFs, singlet excitation energy migrates across TBAPy linkers until it finds a TCPP unit. TCPP acts as an energy trap and quenches the excitation. The efficiency of TBAPy-to-TBAPy energy transfer was found to be sensitive to the structural parameters of MOFs. Analysis of steady-state and time-resolved fluorescence data revealed that excitation energy transfer (EET) is most efficient in ROD-7, followed by NU-901 and NU-1000. We propose that topology that invokes the shorter interchromophoric distances between TBAPy linkers in ROD-7 is responsible for its higher EET efficiency. The distance dependence of the EET rate constant (kEET) was investigated to gain insight into the mechanistic aspects of energy transfer in MOFs. This study revealed that (a) energy transfer in MOFs deviates from the classical Förster model and (b) the geometrical arrangement of linkers influences the mechanism of EET in MOFs. A theoretical investigation was also performed to determine energy-transfer rate constants along different directions and assess the directionality of energy transfer in these MOFs. The magnitude of rate constants indicated that energy transfer in ROD-7 should be highly anisotropic along the stacking direction. Furthermore, these findings suggest that ROD-7 is a promising candidate to play the role of the light-harvesting and energy-transfer component in solar energy conversion devices, where directional energy transfer is required.

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

Computing equilibrium free energies through a nonequilibrium quench

Many methods to accelerate sampling of molecular configurations are based on the idea that temperature can be used to accelerate rare transitions. These methods typically compute equilibrium properties at a target temperature using reweighting or through Monte Carlo exchanges between replicas at higher temperatures. A recent paper [G. M. Rotskoff and E. Vanden-Eijnden, Phys. Rev. Lett. 122, 150602 (2019)] demonstrated that accurate equilibrium densities of states can also be computed through a nonequilibrium “quench” process, where sampling is performed at a higher temperature to encourage rapid mixing and then quenched to lower energy states with dissipative dynamics. Here, we provide an implementation of the quench dynamics in LAMMPS and evaluate a new formulation of nonequilibrium estimators for the computation of partition functions or free energy surfaces (FESs) of molecular systems. We show that the method is exact for a minimal model of N-independent harmonic springs and use these analytical results to develop heuristics for the amount of quenching required to obtain accurate sampling. We then test the quench approach on alanine dipeptide, where we show that it gives an FES that is accurate near the most stable configurations using the quench approach but disagrees with a reference umbrella sampling calculation in high FE regions. We then show that combining quenching with umbrella sampling allows the efficient calculation of the free energy in all regions. Moreover, by using this combined scheme, we obtain the FES across a range of temperatures at no additional cost, making it much more efficient than standard umbrella sampling if this information is required. Finally, we discuss how this approach can be extended to solute tempering and demonstrate that it is highly accurate for the case of solvated alanine dipeptide without any additional modifications.

Chemistry↗