Search NASA⌕ Search

SEARCH · Search NASA

Results for “umbrella sampling”

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 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↗
DOE OSTI2026

Molecular Mechanisms Underlying Surfactant-Based Plastics De-Inking

Surfactant-mediated binder removal is critical for de-inking pretreatment in plastics recycling, yet the molecular mechanisms governing surfactant performance remain poorly understood. We used all-atom (AA) and coarse-grained (CG) molecular dynamics (MD) simulations alongside alkaline surfactant washing experiments to investigate interactions between a series of surfactants and a polyether urethane (PEU) binder in solution and on a polyethylene (PE) surface. Experiments reveal a range of de-inking efficiencies, ranging from <25 to >95% depending upon surfactant headgroup charge and tail length. AA simulations reveal that charged surfactants reach stable levels of surfactant coverage, while nonionic surfactants aggregate on the binder. CG umbrella sampling calculations quantify the thermodynamics of binder desorption in water. In ∼0.25 M surfactant solutions, up to a 52% reduction in the free energy barrier is computed, with trends in good agreement (R 2 = 0.92, Pearson’s r = –0.96, Spearman’s ρ = –0.83) with experimental de-inking efficiencies. We find that charged surfactants are more effective than nonionic surfactants for de-inking and propose three regimes of surfactant de-inking processes: good de-inking occurs in surfactants that promote PEU desorption with a low radius of gyration (R g ); moderate de-inking occurs when surfactants stabilize PEU but increase R g ; poor de-inking occurs in surfactants that aggregate on the binder and promote extensive anchoring to the surface. Overall, these molecular-level insights have the potential to guide the design of surfactant formulations for plastics recycling applications.

ink removal↗
DOE OSTI2024

Thermodynamic Analysis of Silk Fibroin–Graphite Hybrid Materials and Their Morphology

Silk fibroin (SF) is a β-sheet-rich protein that is responsible for the remarkable tensile strength of silk. In addition to its mechanical properties, SF is biocompatible and biodegradable, making it an attractive candidate for use in biotic/abiotic hybrid materials. A pairing of particular interest is the use of SF with graphene-based nanomaterials (GBNs). The properties of this interface drive the formation of well-ordered nanostructures and can improve the electronic properties of the resulting hybrid. It was previously demonstrated that SF can form lamellar nanostructures in the presence of graphite; however, the equilibrium morphology and associated driving interactions are not fully understood. Here, in this study, we characterize these interactions between SF and SF lamellar with graphite using molecular dynamics (MD) simulations and umbrella sampling (US). We find that SF lamellar nanostructures have strong orientational and spatial preferences on graphite that are driven by the hydrophobic effect, destabilizing solvent–protein interactions and stabilizing protein–protein and protein–graphite interactions. Finally, we show how careful consideration of these underlying interactions can be applied to rationally modify the nanostructure morphology.

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

Molecular Modeling and Molecular Dynamics Simulation of a Packed and Intact Bacterial Microcompartment

Bacterial microcompartments (BMCs) are protein-bound organelles found in some bacteria which encapsulate enzymes for enhanced catalytic activity. These compartments spatially sequester enzymes within semipermeable shell proteins and are packed full of enzyme cargoes and metabolites as they fulfill their function. Coupling together recent SAXS and proteomics work, it is possible to develop molecular models for these microcompartments and interrogate enzyme and metabolite dynamics within. Our primary goal of this study is to quantify the permeability of metabolite glyceraldehyde-3-phosphate (G3P) and dihydroxyacetone phosphate (DHAP) across the BMC shell through classical molecular dynamics simulation. The Haliangium ochraceum model of BMC shell (PDB: 6MZX) was used to model an intact BMC of approximately 10 million atoms. Working at this scale presented its own challenges in managing large data sets, with multiple challenges and hardware advances discussed that facilitated this work. Over approximately 750 ns of aggregate simulation, we see multiple permeation events for these metabolites that were added at high concentration through the pores present within BMC shell tiles. When compared to independent permeability estimates for the same metabolites determined through replica exchange umbrella sampling simulations, the permeabilities varied by approximately 3 orders of magnitude. Regardless, the permeability coefficients for both G3P and DHAP are highly similar and very high, such that only very small concentration gradients can be maintained across the BMC shell between the cytosol and BMC interior. The large simulation systems also facilitated comparisons for molecular diffusivity in the crowded environment within the BMC shell. By our estimates, the viscosity within a packed BMC shell is at least 10-fold higher than it would be in neat solution and is the real driver for varying permeability estimates we obtained through simulation. These findings will be used as design inputs for future bioengineering efforts to make products from BMCs, highlighting how permeable BMC shells can be.

Diffusion↗
DOE OSTI2023

Solvent effects on extractant conformational energetics in liquid–liquid extraction: a simulation study of molecular solvents and ionic liquids

Extractant design in liquid–liquid extraction (LLE) is a research frontier of metal ion separations that typically focuses on the direct extractant–metal interactions. However, a more detailed understanding of energetic drivers of separations beyond primary metal coordination is often lacking, including the role of solvent in the extractant phase. In this work, we propose a new mechanism for enhancing metal-complexant energetics with nanostructured solvents. Using molecular dynamics simulations with umbrella sampling, we find that the organic solvent can reshape the energetics of the extractant's intramolecular conformational landscape. Here, we calculate free energy profiles of different conformations of a representative bidentate extractant, n-octyl(phenyl)-N,N-diisobutyl carbamoyl methyl phosphinoxide (CMPO), in four different solvents: dodecane, tributyl phosphate (TBP), and dry and wet ionic liquid (IL) 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([EMIM][Tf 2 N]). By promoting reorganization of the extractant molecule into its binding conformation, our findings reveal how particular solvents can ameliorate this unfavorable step of the metal separation process. In particular, the charge alternating nanodomains formed in ILs substantially reduce the free energy penalty associated with extractant reorganization. Importantly, using alchemical free energy calculations, we find that this stabilization persists even when we explicitly include the extracted cation. These findings provide insight into the energetic drivers of metal ion separations and potentially suggest a new approach to designing effective separations using a molecular-level understanding of solvent effects.

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

Investigating the Effect of Water on the Mechanical Properties of Cellulose from Multiscale Molecular Dynamics Simulations

Classical molecular dynamics (MD) simulations provide insight into the structure and physicochemical properties of materials with atomic resolution. However, the length and time scales accessible to atomistic MD are orders of magnitude smaller than many relevant processes such as the response of a bulk material to experimentally accessible strain rates, which presents challenges when comparing models to experimental measurements. Bottom-up coarse-graining provides a means for systematically mapping atomistic information to lower resolution models to increase the length and time scales achievable by simulation. Cellulose is an abundant carbohydrate biopolymer with applications to many fields of research, such as materials science and renewable energy, due to its desirable mechanical properties and viability for conversion into biofuel. The effect of moisture content on the Young's modulus of cellulose is of special interest due to its native environment often being in the hydrated secondary plant cell wall and the grinding energy requirements for biomass feedstock preprocessing. The current work investigates the effects of water solvent on the Young's modulus of cellulose calculated from coarse-grained MD mechanical stress simulations. The coarse-grained model was parametrized from atomistic MD calculations of cellulose-cellulose potentials of mean force using umbrella sampling techniques under vacuum and solvated conditions. The Young's moduli of the coarse-grained cellulose assemblies parametrized from cellulose in vacuum or solvated in water were computed via mechanical stress simulations to highlight the importance of capturing solvent interactions for modeling the mechanical behavior of cellulose.

BASIC BIOLOGICAL SCIENCES,RADIATION PROTECTION AND↗
DOE OSTI2024

SARS-CoV-2 wastewater variant surveillance: pandemic response leveraging FDA’s GenomeTrakr network

ABSTRACT Wastewater surveillance has emerged as a crucial public health tool for population-level pathogen surveillance. Supported by funding from the American Rescue Plan Act of 2021, the FDA‘s genomic epidemiology program, GenomeTrakr, was leveraged to sequence SARS-CoV-2 from wastewater sites across the United States. This initiative required the evaluation, optimization, development, and publication of new methods and analytical tools spanning sample collection through variant analyses. Version-controlled protocols for each step of the process were developed and published on protocols.io. A custom data analysis tool and a publicly accessible dashboard were built to facilitate real-time visualization of the collected data, focusing on the relative abundance of SARS-CoV-2 variants and sub-lineages across different samples and sites throughout the project. From September 2021 through June 2023, a total of 3,389 wastewater samples were collected, with 2,517 undergoing sequencing and submission to NCBI under the umbrella BioProject,PRJNA757291. Sequence data were released with explicit quality control (QC) tags on all sequence records, communicating our confidence in the quality of data. Variant analysis revealed wide circulation of Delta in the fall of 2021 and captured the sweep of Omicron and subsequent diversification of this lineage through the end of the sampling period. This project successfully achieved two important goals for the FDA’s GenomeTrakr program: first, contributing timely genomic data for the SARS-CoV-2 pandemic response, and second, establishing both capacity and best practices for culture-independent, population-level environmental surveillance for other pathogens of interest to the FDA. IMPORTANCE This paper serves two primary objectives. First, it summarizes the genomic and contextual data collected during a Covid-19 pandemic response project, which utilized the FDA’s laboratory network, traditionally employed for sequencing foodborne pathogens, for sequencing SARS-CoV-2 from wastewater samples. Second, it outlines best practices for gathering and organizing population-level next generation sequencing (NGS) data collected for culture-free, surveillance of pathogens sourced from environmental samples.

Microbiology↗
DOE OSTI2023

Realistic application of short-lived fission product delayed neutron, gamma-ray analysis for simultaneous nondestructive trace quantification of U, Pu mixtures on cellulose swipes

Detection and characterization of fissile traces are of interest to the international nuclear nonproliferation community, including the International Atomic Energy Agency. Pre-inspection check samples are analyzed by neutron activation analysis at the High Flux Isotope Reactor operated by the Oak Ridge National Laboratory under the umbrella of the IAEA Network of Analytical Laboratories. The simultaneous quantification of U and Pu mixtures was accomplished using the combined delayed neutron (DN) delayed gamma-ray (DG) method to analyze cellulose swipes with actinide loading <1ng in a blind field trial. The total fissile quantity was measured by the DN counts and the relative proportions of U, Pu were determined by calibration of the 104 Tc / 141 Ba fission product count ratio using known mixtures. The DNDG method demonstrated high accuracy in flagging the presence of 239 Pu in uranium down to <100 pg mass loading. In conclusion, peak significance tests helped to control false positive Pu flagging and simultaneous quantification of U and Pu loading was accomplished on samples that passed the significance tests.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Learnability Transitions in Monitored Quantum Dynamics via Eavesdropper’s Classical Shadows

Monitored quantum dynamics—unitary evolution interspersed with measurements—has recently emerged as a rich domain for phase structure in quantum many-body systems away from equilibrium. Here we study monitored dynamics from the point of view of an eavesdropper who has access to the classical measurement outcomes, but not to the quantum many-body system. We show that a measure of information flow from the quantum system to the classical measurement record—the —undergoes a phase transition in correspondence with the measurement-induced phase transition (MIPT). This transition determines the eavesdropper’s (in)ability to learn properties of an unknown initial quantum state of the system, given a complete classical description of the monitored dynamics and arbitrary classical computational resources. We make this learnability transition concrete by defining classical shadow protocols that the eavesdropper may apply to this problem, and show that the MIPT manifests as a transition in the sample complexity of various shadow-estimation tasks, which become harder in the low-measurement phase. We focus on three applications of interest: Pauli expectation values (where we find the MIPT appears as a point of optimal learnability for typical Pauli operators), many-body fidelity, and global charge in U ( 1 ) -symmetric dynamics. Our work unifies different manifestations of the MIPT under the umbrella of and gives this notion a general operational meaning via classical shadows. Published by the American Physical Society 2024

Physics↗