Search NASA⌕ Search

SEARCH · Search NASA

Results for “target configuration”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Towards self-consistent integrated modeling of the tokamak pedestal, scrape-off layer, and divertor using SOLPS-ITER and EPED

A predictive, self-consistent integrated model of the tokamak edge plasma and neutral system has been developed and validated against experimental data from DIII-D, spanning from the divertor target and scrape-off layer to the top of the pedestal. The model uses the SOLPS-ITER fluid boundary plasma and neutral code to constrain the EPED prediction of the pedestal height and width by taking the separatrix temperature and density calculated from SOLPS as inputs to EPED. An empirical proxy function based on data from a subset of these discharges is then used to infer the ratio of density at the pedestal top to that at the separatrix based on the electron temperature at the divertor target. This allows the density at the top of the pedestal to be provided to EPED as well, so that knowledge of this quantity is no longer required a priori, while SOLPS radial transport coefficients are needed instead. The gas puff rate provided to SOLPS is then a primary model input for this forward model coupling (along with radial transport coefficients, which are kept fixed through each scan in this study). Comparisons with measurements from DIII-D density scans through a range of target conditions find significant pedestal pressure degradation as the density is increased towards detachment in three different divertor configurations, consistent with the experimental data in these ballooning-limited pedestal regimes.

core-edge integration↗
DOE OSTI2021

Very Large Array imaging rules out precessing radio jets in three DES–SDSS-selected candidate periodic quasars

ABSTRACT Periodic quasars have been suggested as candidates for hosting binary supermassive black holes (SMBHs), although alternative scenarios remain possible to explain the optical light-curve periodicity. To test the alternative hypothesis of precessing radio jet, we present deep 6 GHz radio imaging conducted with NSF’s Karl G. Jansky Very Large Array (VLA) in its C configuration for the three candidate periodic quasars, DES J024703.24−010032.0, DES J024944.66−000036.8, and DES J025214.67−002813.7. Our targets were selected based on their optical variability using 20 yr long multicolour light curves from the Dark Energy Survey (DES) and the Sloan Digital Sky Survey (SDSS). The new VLA observations show that all three periodic quasars are radio-quiet with the radio loudness parameters measured to be $R\equiv f_{6\, {\rm cm}}/f_{{\rm 2500}}$ of ≲1.0–1.5 and the k-corrected luminosities νLν[6 GHz] of ≲5–21 × 1039 erg s−1. They are in stark contrast to previously known periodic quasars proposed as binary SMBH candidates such as the blazar OJ 287 and PG1302−102. Our results rule out optical emission contributed from precessing radio jets as the origin of the optical periodicity in the three DES–SDSS-selected candidate periodic quasars. Future continued optical monitoring and complementary multiwavelength observations are still needed to further test the binary SMBH hypothesis as well as other competing scenarios to explain the optical periodicity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Data, model inputs, and analysis scripts associated with a manuscript on stream intermittency controls across spatial scales in Pacific Northwest watersheds

NOTE: The manuscript associated with this data package is currently in review. The data may be revised based on reviewer feedback. Upon manuscript acceptance, this data package will be updated with the final dataset and additional metadata. This data package is associated with the manuscript "Hydroclimatic Memory and Watershed Template Shape Stream Intermittency: Multi-scale Attribution Using Process-based Simulation and Explainable ML" by Niroula et al. (2026), submitted to Water Resources Research (WRR). The study investigates the dominant controls on stream intermittency across local, reach, and watershed scales using a coupled process-based simulation and explainable machine-learning framework. Long-term daily simulations from the Advanced Terrestrial Simulator (ATS) were used to generate wetness states and ponded-depth responses over river-corridor cells. These ATS outputs were then aggregated across scales and used to train XGBoost (eXtreme Gradient Boosting) models. SHAP (SHapley Additive exPlanations) was applied to quantify the relative importance of hydroclimatic forcings, watershed template attributes, and antecedent-memory effects in shaping intermittency behavior. The analysis is carried out for three contrasting Pacific Northwest watersheds: Oak Creek (OCW), American River Watershed (ARW), and H.J. Andrews (HJA). Across these testbeds, the package contains ATS-ready watershed inputs, ATS run configuration and selected output files, model-evaluation data products, intermittency-analysis datasets, machine-learning target-feature tables, SHAP outputs, and notebooks used to organize, analyze, and visualize results. At a high level, the package documents a workflow in which ATS provides the physically based simulation backbone and explainable machine learning is used as a post-processing attribution tool. The contents are intended to support interpretation of the manuscript figures and results, provide context for how intermittency metrics were generated at multiple scales, and preserve the key artifacts needed to understand and reuse the analysis workflow. The package contains a high-level directory summary file (`summary.txt`) and four main content folders (1) `evaluation_plots` contains evaluation figures and supporting evaluation datasets; (2) `intermittency_plots` contains intermittency-focused analysis notebook and prepared datasets; (3) `ml-training-and-shap_values_plots` contains ML training inputs, SHAP outputs, and figure-generation notebooks; and (4) `watershed_mesh_and_ats_input` contains ATS model setup materials, forcing inputs, geometry, and selected run files. More specifically, the `evaluation_plots` folder contains the notebook used for ATS evaluation plotting and site-specific evaluation datasets. These include evapotranspiration and water-balance products for three watersheds, as well as an Oak Creek field-measurement discharge file. The `intermittency_plots` folder contains the notebook used for intermittency analysis and the prepared datasets used to analyze intermittent and non-intermittent wetness behavior across the study watersheds. The `ml-training-and-shap_values_plots` folder contains notebooks and outputs for the machine-learning and explainability workflow. This includes the main XGBoost and SHAP notebook(s), a beeswarm plotting notebook, target-feature tables for machine-learning training, SHAP summary tables, and per-sample SHAP value archives. The `watershed_mesh_and_ats_input` folder contains ATS-related watershed inputs and supporting materials. This includes mesh and shape products, ATS-readable LAI and meteorological forcing inputs, selected ATS spinup and transient-run files, and a watershed workflow example notebook. Subdirectories are organized by watershed where applicable.All files are .cpg (codepage files), .csv (comma-separated values), .dbf (database files), .exo (Exodus mesh format), .h5 (HDF5 format), .ipynb (Jupyter notebooks), .pkl (Python pickle), .prj (projection files), .sh (shell scripts), .shp (shapefile geometry), .shx (shapefile index), .txt (text files), or .xml (markup data).

Advanced Terrestrial Simulator↗
DOE OSTI2026

Machine learning at the Spallation Neutron Source accelerator and target

We describe the ongoing efforts to apply Machine Learning techniques to improve the performance of our accelerator and target. Specially, we are looking to minimize halo beam losses in the absence of a proper physics model, automatically detect and log anomalies in the target support systems such as cooling, and detect and prevent errant beam pulses in the linac. We also describe the infrastructure we use to acquire and stream data to the GPU cluster for training, our code development cycle, and edge computing for model inference. To minimize halo beam losses, we use a Reinforcement Learning technique tested on a virtual accelerator. The target anomaly detection is trained on archived data using incomplete physics models and is made part of the existing target reporting system. The errant beam prevention analyzes beam current and beam phase waveforms as well as accelerator configuration data to predict errant pulses. We also develop continual learning to adapt to changes in the accelerator.

Accelerator Physics↗
DOE OSTI2020

COMET/ZEUS Activation Foils Report for Irradiations IER-241 and IER-261

The goal of this work was to determine the internal neutron spectrum of the COMET/ZEUS assembly using activation foils. Similar measurements were performed on the COMET assembly when it was located at the Los Alamos Critical Experiments Facility (LACEF) at the Los Alamos National Laboratory (LANL). The present measurements were performed after the COMET assembly was moved to the National Criticality Experiments Research Center (NCERC) located in the Device Assembly Facility (DAF) at the Nevada National Security Site (NNSS). The December 2012 COMET/ZEUS irradiation was the first in a series of experiments, and the COMET/ZEUS assembly was selected because of operational readiness and authorization to run specific targets. The August 2013 irradiation was the second in the series and involved only small changes in the configuration. This report describes the unfolding of the activation foils to determine the internal neutron spectra.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2025

Online and Offline Data Quality Monitoring for the Mu2e Calorimeter

This thesis presents the design, implementation, and validation of a calorimeter Data Quality Monitoring (DQM) toolchain for the Mu2e experiment at Fermilab. Mu2e searches for charged lepton flavor violation via coherent muon-to-electron conversion in the field of an aluminum nucleus, $\mu^- Al \rightarrow e^-Al$, a process whose observation would constitute clear evidence of physics beyond the Standard Model. Achieving target sensitivity requires stringent control of detector performance and data integrity during acquisition, as subtle issues in readout configuration, data formatting, or electronics behavior can compromise reconstruction and bias downstream analyzes. To address these challenges, this work develops a multi-layer DQM approach spanning both raw data validation and reconstructed digi-level diagnostics. At the low level, a fragment analysis component performs word- and bit-field decoding of calorimeter readout blocks, enabling sanity checks of the expected structure and producing detailed error and integrity statistics useful for commissioning and troubleshooting. At the digi level, the CaloDigiDQM analyzer is implemented within the art framework and transforms each CaloDigiCollection into a structured hierarchy of ROOT histograms designed for fast drill-down diagnostics. The module generates coherent monitoring views at global, disk, board, and channel granularity, including occupancy, waveform-derived features (baseline, RMS, peak amplitude and position), and left-right sensor consistency metrics. Detector-aware channel-to-electronics mapping is performed through the conditions system (CaloDAQMap), ensuring that diagnostics remain aligned with hardware identifiers used in operations. For end-to-end testing without reliance on live DAQ data, a synthetic CaloDigi producer is developed to generate realistic waveforms with controlled noise and pulse shapes. The resulting system supports both offline ROOT-file production and online operation, including optional histogram streaming through otsdaq via ots::HistoSender. This toolchain provides a practical and scalable foundation for calorimeter commissioning and stable data collection, enabling early detection of anomalies and reducing operational risk for Mu2e.

Vakulenko, Mark [Drew U.] (ORCID:0009000276197818)↗
DOE OSTI2024

Liquid lithium divertor analysis using coupled plasma material interaction model

A liquid lithium divertor can improve performance of future fusion devices by creating efficient power exhaust and improving the energy confinement via pumping of the hydrogen isotopes. In addition, significantly higher heat fluxes can be handled if controlled vapor shielding is used to redistribute the divertor heat flux over a wider area. Design and optimization of such a system calls for an analysis model which includes a strong two-way coupling between the plasma and divertor material. The incoming plasma heat and particle flux will affect the divertor surface temperature, which is a defining factor of the lithium evaporative and sputtered flux going into the plasma. Results of the coupled model based on the plasma edge code SOLPS-ITER and the computational fluid dynamics (CFD) code ANSYS-CFX will be presented for different configurations. An analytical slab flow model is used as a heat transfer boundary condition for SOLPS, defining particle flux from the wall via calculation of the surface temperature. At the final step, results of the SOLPS analysis are verified using a 3D CFD magnetohydrodynamics (MHD) analysis which uses heat and particle flux from SOLPS as a boundary condition. In addition to plasma heat flux, both analytical and CFD temperature models include several plasma material interaction effects, such as lithium evaporation, condensation and sputtering based on deuterium target flux. New adatom sputtering model based on the available experimental data is presented. Analytical model is expanded to include free surface axisymmetric configurations. Results of parametric studies of the divertor configurations with different lithium inlet temperature and velocity will be presented leading to the optimal design resulting in the lowest possible lithium contamination in the core.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Multi-scale structural analysis of swift heavy ion-irradiated ZrO 2 and HfO 2

The radiation-induced monoclinic-to-tetragonal phase transition in ZrO 2 and HfO 2 has been the subject of many investigations, but the transformation pathways and underlying structural mechanisms are still not well understood. In this study, microcrystalline powder samples of ZrO 2 and HfO 2 were irradiated with 946 MeV and 1470 MeV Au ions to a wide fluence range up to 3 × 10 13 ions/cm 2 . To characterize beam-induced structural modifications across all spatial length scales, complementary experimental techniques such as synchrotron X-ray diffraction and spallation neutron total scattering were used. The phase evolution of the tetragonal polymorph with increasing ion fluence is accurately described by a heterogeneous track-overlap model that incorporates both direct- and double-impact processes. These damage accumulation processes are an expression of a core-shell ion track morphology that depends on irradiation conditions and target material. Neutron pair distribution function analysis revealed that ion-beam-induced tetragonal ZrO 2 is merely a configurational average of short-range orthorhombic (Pbcn) domains stabilized by a dense network of domain walls. Furthermore, this knowledge is critical for a better understanding of how crystalline-to-crystalline phase transformations proceed at the atomic scale under extreme conditions.

Monoclinic-to-tetragonal phase transformation↗
DOE OSTI2026

Alchemy: A Model-Based Approach for 2D to 3D Autonomous Nuclear System Design

Engineering design of nuclear power plant (NPP) piping and equipment systems frequently bypasses crucial 2D system planning, instead moving straight to 3D modeling. This often leads to designs that exceed building envelope constraints, forcing expensive and time-consuming redesigns. When 2D modeling is employed, it typically involves labor-intensive manual workflows that convert 2D drawings into 3D models, resulting in inefficiencies and errors across design iterations. These workflows further suffer from poor software interoperability and dependence on proprietary software ecosystems, thereby contributing to schedule delays and cost overruns. This paper presents Alchemy, an autonomous framework that transforms 2D system definitions into Industry Foundation Classes (IFC)-compliant 3D building information models (BIMs) for expediting nuclear facility design at the conceptual preliminary phase. Using a model-based approach, the framework treats the 2D system diagram as the central reference model employed to automatically generate all subsequent outputs, ensuring consistency between the system definition and the resulting physical design. A web-based interface enables engineers to define hierarchical system topologies including associated equipment, geometric properties, and connectivity requirements. A two-phase equipment layout optimization algorithm automatically computes collision-free spatial configurations within predefined building envelopes. An artificial intelligence (AI)-assisted pipe routing module then generates orthogonal, collision-free routing paths, allowing the user to select either an A* search-based method or an Ant Colony Optimization (ACO)-based method. All outputs are authored natively in IFC format, relying on open-source technologies and standardized formats in order to ensure extensibility and eliminate proprietary software dependencies. The proposed framework is validated on two representative pressurized-water reactor (PWR)-based case studies, for which it autonomously generates IFC-compliant 3D models in minutes, drastically reducing workflows that typically require hours of manual effort. The generated model demonstrates topologically correct equipment placement, physically plausible spatial relationships, and collision-free pipe routing consistent with known PWR loop configurations. This work represents a foundational step toward digital engineering for nuclear facility preliminary design, with future ongoing development targeting design code compliance and expanded system complexity.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2023

ATF-2 Sensitivity M&C 2023 (Slides)

Irradiation experiments in the Idaho National Laboratory's Advanced Test Reactor are typically assumed to have little effect on one another. This assumption does not hold true for certain experiments in close proximity. To evaluate the impacts on safety and programmatic parameters of experiments in the center flux trap, the contents of the adjacent H and inner-A positions are modeled with a range of possible irradiation targets. First, neutron flux maps with experiments in those positions are compared to those of a water-filled configuration. Then, several safety and programmatic parameters for a generic accident-tolerant fuel test train are calculated. It is shown that these parameters can exhibit considerable sensitivity to the contents of the H and A experiment positions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Organic Rankine Cycle Integration and Optimization for High Efficiency CHP Genset Systems (Final Technical Report)

This project successfully advanced the integration of Organic Rankine Cycle (ORC) technology with reciprocating engine–based combined heat and power (CHP) systems to improve electrical efficiency, total CHP efficiency, and grid-responsive operation. Over three budget periods, the work progressed from high-temperature ORC component development and thermodynamic model validation to next-generation system design, working fluid transition, and techno-economic analysis. Key technical accomplishments include development and validation of a thermodynamic model capable of accurately predicting ORC performance across an expanded temperature and pressure envelope; successful identification and validation of low-global-warming-potential (GWP) working fluids—most notably R1233zd(E)—as viable replacements for R245fa; and demonstration of scalable ORC architectures suitable for integration with 1–20 MW class reciprocating engines. These advances enable flexible CHP configurations that can increase electrical output while maintaining high overall utilization of available thermal energy. The project also produced a clean-sheet design for a next-generation ORC system targeting substantially higher power output per unit, supported by detailed component selection, heat exchanger evaluation, and system-level modeling. Techno-economic analyses indicate that ORC-enabled flexible CHP systems can meet or exceed Department of Energy (DOE) efficiency targets while providing value to both facility operators and the electric grid.. Late-stage testing of the largest next-generation ORC prototype identified limitations related to pump net positive suction head (NPSH) requirements and condenser flooding under certain operating conditions. Although these issues constrained full validation of that configuration within the project timeframe, they provided clear and actionable design guidance for future system refinements. Importantly, validated modeling, smaller-scale testing, and working fluid evaluations confirmed the technical viability of the overall approach. In aggregate, this project met its core objectives by establishing validated design tools, de-risking key ORC technologies for CHP applications, and defining a credible pathway toward commercialization of flexible, high-efficiency CHP systems. The results form a strong foundation for continued development and deployment beyond the conclusion of the DOE-funded effort.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2026

Efficient Optimization of Plasma Radiation Detector Configurations using Imperfect Inference Models

The configurations of instruments fielded on an experiment affect the amount of information captured and the quality of subsequent inference. Here, we investigate the problem of optimizing plasma x-ray radiation detectors in a magneto-inertial fusion experiment at Sandia National Laboratories. It is impossible to directly measure properties such as the temperature of the thermonuclear fusion plasma produced in these experiments because of the extreme environment and destructive nature of the experiment. Among other diagnostics, several detectors are placed with significant standoff from the fusion target to capture the x-rays emitted by the fusion plasma, which can be used to infer some of its properties. To optimize the configuration of these detectors, a high-fidelity model (HFM) is used for simulating outputs and a low-fidelity model (LFM) is used for inference. We develop methods based on A- and L-optimality criteria that are efficient to compute while explicitly accounting for the discrepancy between the HFM and the LFM. The method allows us to find detector configurations that perform similarly to or better than the configuration obtained using an existing sampling-based optimization method while decreasing computational time by a factor of 50. Supplementary materials for this article are available online, including a standardized description of the materials available for reproducing the work.

Bayesian optimization↗
DOE OSTI2024

Software defined networking moving target defense honeypot

Disclosed herein are systems, methods, and storage media for thwarting cyber-attacks and data theft. A computing system receives packets and compares with a configuration resource. The computing system determines that the packet does not match the configuration resource and transmits a packet to a decoy environment via an SDN switch. The decoy environment is configured to generate time-out, service unavailable, or restricted access messages. In some embodiments, the computing system determines the packet does match the configuration and transmits the packet to a production network via an SDN switch. The SDN switch is communicatively connected to the decoy environment through a first channel and communicatively connected to a production environment through a separate second channel. The computing system is further configured to create and transmit to the SDN switch, rules to manage transmitting the packet to either the decoy environment or the production environment.

Lyle, Joshua A.↗
DOE OSTI2022

Unveiling the complex configurational landscape of the intralayer cavities in a crystalline carbon nitride

The in-depth understanding of the reported photoelectrochemical properties of the layered carbon nitride, poly(triazine imide)/LiCl (PTI/LiCl), has been limited by the apparent disorder of the Li/H atoms within its framework. To understand and resolve the current structural ambiguities, an optimized one-step flux synthesis (470 °C, 36 h, LiCl/KCl flux) was used to prepare PTI/LiCl and deuterated-PTI/LiCl in high purity. Its structure was characterized by a combination of neutron/X-ray diffraction and transmission electron microscopy. The range of possible Li/H atomic configurations was enumerated for the first time and, combined with total energy calculations, reveals a more complex energetic landscape than previously considered. Experimental data were fitted against all possible structural models, exhibiting the most consistency with a new orthorhombic model (Sp. Grp. Ama2) that also has the lowest total energy. In addition, a new Cu(I)-containing PTI (PTI/CuCl) was prepared with the more strongly scattering Cu(I) cations in place of Li, and most closely matching with the partially-disorder structure in Cmc2 1 . Thus, a complex configurational landscape of PTI is revealed to consist of a number of ordered crystalline structures that are new potential synthetic targets, such as with the use of metal-exchange reactions.

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

Development of a hardened THz energy meter for use on the kilojoule-scale, short-pulse OMEGA EP laser

A highly adaptable and robust terahertz (THz) energy meter is designed and implemented to detect energetic THz pulses from high-intensity (>10 18 W/cm 2 ) laser–plasma interactions on the OMEGA EP. THz radiation from the laser driven target is detected by a shielded pyrometer. A second identical pyrometer is used for background subtraction. The detector can be configured to detect THz pulses in the 1 mm to 30 μm (0.3- to 10-THz) range and pulse energies from joules to microjoules via changes in filtration, aperture size, and position. Additional polarization selective filtration can also be used to determine the THz pulse polarization. The design incorporates significant radiation and electromagnetic pulse shielding to survive and operate within the OMEGA EP radiation environment. We describe the design, operational principle, calibration, and testing of the THz energy meter. The pyrometers were calibrated using a benchtop laser and show linear sensitivity to up to 1000 nJ of absorbed energy. Furthermore, the initial results from four OMEGA EP THz experiments detected up to ~15 μJ at the detector, which can correspond to hundreds of mJ depending on THz emission and reflection models.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Effect of magnetic field ripple on parallel electron transport during microwave plasma heating in the Proto-MPEX linear plasma device

The Prototype Material Plasma Exposure eXperiment (Proto-MPEX) is a linear plasma device located at Oak Ridge National Laboratory to develop the plasma source concept for the Materials Plasma Exposure eXperiment (MPEX). Recent experiments have demonstrated the heating of electrons in a high-density deuterium plasma (n e ~ 5 × 10 19 m -3 ) via electron Bernstein waves and source plasma production via helicon waves in this linear configuration. Moreover, experimental observations suggest that the magnetic ripple adversely affects the parallel transport of the heated electrons toward the target where material samples are to be exposed to the plasma. In order to understand the transport process during microwave application, a test-particle Monte-Carlo (MC) code has been developed that incorporates the effects of Coulomb collisions, magnetic mirror adiabatic trapping and electron cyclotron interaction via a quasilinear RF heating operator. MC calculations indicate that the absence of 2nd turning points along the trajectory of cyclotron heated electrons significantly reduces adiabatic trapping and maximizes power transport to the Target surface. Thermalization of cyclotron heating electrons is analyzed with the MC code and it is found that for conditions relevant to ProtoMPEX, a significant fraction of the absorbed heating power is coupled to the Target surface via ballistic fast electrons; however, under certain conditions, up to 60% of the absorbed RF power can be dissipated via thermalization of the fast electrons on the background plasma. Results are compared with experimental measurements. Finally, the implications of the results are discussed in the context of the upcoming MPEX device.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

First Observation of the β3αp Decay of 13 O via β-Delayed Charged-Particle Spectroscopy

The β-delayed proton decay of 13 O has previously been studied, but the direct observation of β-delayed 3⁢α⁢p decay has not been reported. Rare 3⁢α⁢p events from the decay of excited states in 13 N* provide a sensitive probe of cluster configurations in 13 N*. To measure the low-energy products following β-delayed 3⁢αp decay, the Texas Active Target (TexAT) time projection chamber was employed using the one-at-a-time β-delayed charged-particle spectroscopy technique at the Cyclotron Institute, Texas A&M University. A total of 1.9 × 10 5 13 O implantations were made inside the TexAT time projection chamber. Furthermore, a total of 149 3⁢αp events were observed, yielding a β-delayed 3⁢αp branching ratio of 0.078(6)%. Four previously unknown α-decaying excited states were observed in 13 N at 11.3, 12.4, 13.1, and 13.7 MeV decaying via the 3⁢α + p channel.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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↗