Search NASA⌕ Search

SEARCH · Search NASA

Results for “entropic sampling methods”

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.

10 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Accelerating Ultrafast Spectroscopy with Compressive Sensing

We report ultrafast spectroscopy is an important tool for studying photoinduced dynamical processes in atoms, molecules, and nanostructures. Typically, the time to perform these experiments ranges from several minutes to hours depending on the choice of spectroscopic method. It is desirable to reduce this time overhead not only to shorten time and laboratory resources, but also to make it possible to examine fragile specimens that quickly degrade during long experiments. In this article, we motivate using compressive sensing to significantly shorten data acquisition time by reducing the total number of measurements in ultrafast spectroscopy. We apply this technique to experimental data from ultrafast transient absorption spectroscopy and ultrafast terahertz spectroscopy and show that good estimates can be obtained with as low as 15% of the total measurements, implying a sixfold reduction in the data acquisition time.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Biased degenerate ground-state sampling of small Ising models with converged quantum approximate optimization algorithm

The quantum alternating operator ansatz, a generalization of the quantum approximate optimization algorithm (QAOA), is a quantum algorithm used for approximately solving combinatorial optimization problems. QAOA typically uses the transverse field mixer as the driving Hamiltonian. One of the interesting properties of the transverse field driving Hamiltonian is that it results in nonuniform sampling of degenerate ground states of optimization problems. In this study, we numerically examine the fair sampling properties of the transverse field mixer QAOA, and Grover mixer QAOA (GM-QAOA), which provides theoretical guarantees of fair sampling of degenerate optimal solutions, up to a large enough p such that the mean expectation value converges to an optimal approximation ratio of 1. This comparison is performed with high-quality heuristically computed, but not necessarily optimal, QAOA angles, which give strictly monotonically improving solution quality as p increases. These angles are computed using the Julia based numerical simulation software JuliQAOA. Fair sampling of degenerate ground states is quantified using the Shannon entropy of the ground-state amplitudes distribution. The fair sampling properties are reported on several quantum signature Hamiltonians from previous quantum annealing fair sampling studies. Small random fully connected spin glasses are shown, which exhibit exponential suppression of some degenerate ground states with transverse field mixer QAOA. The transverse field mixer QAOA simulations show that some problem instances clearly saturate the Shannon entropy of 0 with a maximally biased distribution that occurs when the learning converges to an approximation ratio of 1 while other problem instances never deviate from a maximum Shannon entropy (uniform distribution) at any p step. Published by the American Physical Society 2025

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Ab initio molecular dynamics with enhanced sampling in heterogeneous catalysis

Ab initio molecular dynamics simulations combined with enhanced sampling techniques are becoming widespread methods to investigate chemical phenomena in catalytic systems. These techniques automatically include finite temperature effects, anharmonicity, and collective dynamics in their robust description of enthalpic and entropic contributions, which can have significant impact on reaction free energy landscapes. This contrasts with standard ab initio static approaches that are based on assessing reaction free energies from various coarse-grained descriptions of the reaction potential energy surface. Enhanced sampling ab initio molecular dynamics opens the way to first principles simulations of systems of increasing complexity like solid/liquid catalytic interfaces. Here, we aim at guiding the reader through the basis of these techniques, summarizing their fundamental theoretical and practical aspects, and reviewing the relevant literature in the field. After a brief introduction to the problem, we will illustrate the advantage of using molecular simulations to include finite temperature effects, examine the most common ab initio techniques currently in use, describe their application to solid state heterogeneous catalysts, and finally critically review the most popular enhanced sampling techniques used in computational catalysis.

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

First-principles calculation of the configurational energy density of states for a solid-state ion conductor with a variant of the Wang and Landau algorithm

In this work, a variant of the Wang and Landau algorithm for calculation of the configurational energy density of states is proposed. The algorithm was developed for the purpose of using first-principles simulations, such as density functional theory, to calculate the partition function of disordered sublattices in crystal materials. The expensive calculations of first-principles methods make a parallel algorithm necessary for a practical computation of the configurational energy density of states within a supercell approximation of a solid-state material. The algorithm developed in this work is tested with the two-dimensional (2d) Ising model to bench mark the algorithm and to help provide insight for implementation to a materials science application. Tests with the 2d Ising model revealed that the algorithm has good performance compared to the original Wang and Landau algorithm and the 1/$\textit{t}$ algorithm, in particular the short iteration performance. Further, a proof of convergence is presented within an adiabatic assumption, and the analysis is able to correctly predict the time dependence of the modification factor to the density of states. The algorithm was then applied to the lithium and lanthanum sublattice of the solid-state lithium ion conductor Li 0.5 La 0.5 TiO 3 . This was done to help understand the disordered nature of the lithium and lanthanum. The results find, overall, that the algorithm performs very well for the 2d Ising model and that the results for Li 0.5 La 0.5 TiO 3 are consistent with experiment while providing additional insight into the lithium and lanthanum ordering in the material. The primary result is that the lithium and lanthanum become more mixed between layers along the c axis for increasing temperature. In part, the simulation of the disordered Li 0.5 La 0.5 TiO 3 system serves as a benchmark for what size systems are currently and in the near future practical to calculate with density functional theory methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Using Metadynamics to Reveal Extractant Conformational Free Energy Landscapes

Understanding the impact of extractant functionalization on metal-binding energetics in liquid-liquid extraction is essential to guide the development of better separation processes. Traditionally, computational extractant design uses electronic structure calculations on metal-ligand clusters to determine the metal-binding energy of the lowest energy state. Although highly accurate, this approach does not account for all of the relevant physics encountered under experimental conditions. Such methodologies often neglect entropic contributions such as temperature effects and ligand flexibility, in addition to approximating solvent-extractant interactions with implicit solvent models. In this study, we use classical molecular dynamics simulations with an advanced sampling method, metadynamics, to map out extractant molecule conformational free energies in the condensed phase. Here we generate the complete conformational landscape in solution for a family of bidentate malonamide-based extractants with different functionalizations of the headgroup and the side chains. In particular, we show how such alkyl functionalization reshapes the free energy landscape, affecting the free energy penalty of organizing the extractant into the cis-like metal-binding conformation from the trans-like conformation of the free extractant in solution. Specifically, functionalizing alkyl tails to the center of the headgroup has a greater influence on increasing molecular rigidity and disfavoring the binding conformation than functionalizing side chains. These findings are consistent with trends in metal-binding energetics based on experimentally reported distribution ratios. We also consider a different bidentate extractant molecule, carbamoylmethylphosphine oxide, and show how the choice of solvent can further reshape the conformational energetic landscape. This study demonstrates the feasibility of using molecular dynamics simulations with advanced sampling techniques to investigate extractant conformational energetics in solution, which, more broadly, will enable extractant design that accounts for entropic effects and explicit solvation.

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

Thermodynamic Basis for the Stabilization of Helical Peptoids by Chiral Sidechains

Peptoids are a class of highly customizable biomimetic foldamers that retain properties from both proteins and polymers. It has been shown that peptoids can adopt peptide-like secondary structures through the careful selection of sidechain chemistries, but the underlying conformational landscapes that drive these assemblies at the molecular level remain poorly understood. Given the high flexibility of the peptoid backbone, it is essential that methods applied to study peptoid secondary structure formation possess the requisite sensitivity to discriminate between structurally similar yet energetically distinct microstates. In this work, a generalizable simulation scheme is used to robustly sample the complex folding landscape of various 12mer polypeptoids, resulting in a predictive model that links sidechain chemistry with preferential assembly into one of 12 accessible backbone motifs. Using a variant of the metadynamics sampling method, four peptoid dodecamers are simulated in water: sarcosine, N-(1-phenylmethyl)glycine (Npm), (S)-N-(1-phenylethyl)glycine (Nspe), and (R)-N-(1-phenylethyl)glycine (Nrpe)–to determine the underlying entropic and energetic impacts of hydrophobic and chiral peptoid sidechains on secondary structure formation. Our results indicate that the driving forces to assemble Nrpe and Nspe sequences into polyproline type-I helices in water are found to be enthalpically driven, with small benefits from an entropic gain for isomerization and steric strain due to the presence of the chiral center. The minor entropic gains from bulky chiral sidechains in Nrpe- and Nspe-containing peptoids can be explained through increased configurational entropy in the cis state. However, overall assembly into a helix is found to be overall entropically unfavorable. Furthermore, these results highlight the importance of considering the many various competing interactions in the rational design of peptoid secondary structure building blocks.

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

Quantifying Energetic and Entropic Pathways in Molecular Systems

When examining dynamics occurring at nonzero temperatures, both energy and entropy must be taken into account to describe activated barrier crossing events. Furthermore, good reaction coordinates need to be constructed to describe different metastable states and the transition mechanisms between them. Here we use a physics-based machine learning method called state predictive information bottleneck (SPIB) to find nonlinear reaction coordinates for three systems of varying complexity. SPIB is able to correctly predict an entropic bottleneck for an analytical flat-energy double-well system and identify the entropy- and energy-dominated pathways for an analytical four-well system. Lastly, for a simulation of benzoic acid permeation through a lipid bilayer, SPIB is able to discover the entropic and energetic barriers to the permeation process. Given these results, we thus establish that SPIB is a reasonable and robust method for finding the important entropy, energy, and enthalpy barriers in physical systems, which can then be used to enhance the understanding and sampling of different activated mechanisms.

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

Effect of Collective Dynamics and Anharmonicity on Entropy in Heterogenous Catalysis: Building the Case for Advanced Molecular Simulations

In this paper, we present a perspective on the computational determination of entropy and its effects and consequences on heterogeneous catalysis. Special attention is paid to the role of anharmonicity and the deviations from the standard harmonic oscillator approximations which can fail to provide a reliable assessment of entropy. To address these challenges, advanced methodologies are needed that can explicitly account for these thermodynamic drivers through the appropriate statistical sampling of reactive free energy surfaces. We discuss where anharmonicity should be expected, where it has been observed from a theoretical perspective, and the methods currently employed to address it. We concentrate on three types of systems where we have observed major, nonnegligible anharmonic effects: (1) supported nanoparticles, where the migration of metal atoms, complexes, and entire clusters exhibit anharmonic behavior in their dynamic motion; (2) porous solids where confinement effects distort potential energy surfaces and hinder molecular motions resulting in large entropic terms; and (3) solid/liquid interfaces, where interactions between solvent molecules and adsorbed species can result in large solvent organization free energy and unique reactivity.

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

A Computational Information Criterion for Particle-Tracking with Sparse or Noisy Data

Traditional probabilistic methods for the simulation of advection-diffusion equations (ADEs) often overlook the entropic contribution of the discretization, e.g., the number of particles, within associated numerical methods. Many times, the gain in accuracy of a highly discretized numerical model is outweighed by its associated computational costs or the noise within the data. Herein, we address the question of how many particles are needed in a simulation to best approximate and estimate parameters in one-dimensional advective-diffusive transport. To do so, we use the well-known Akaike Information Criterion (AIC) and a recently-developed correction called the Computational Information Criterion (COMIC) to guide the model selection process. Random-walk and mass-transfer particle tracking methods are employed to solve the model equations at various levels of discretization. Numerical results demonstrate that the COMIC provides an optimal number of particles that can describe a more efficient model in terms of parameter estimation and model prediction compared to the model selected by the AIC even when the data is sparse or noisy, the sampling volume is not uniform throughout the physical domain, or the error distribution of the data is non-IID Gaussian.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1999

Engineering of Novel Biocolloid Suspensions

Colloidal suspensions are materials with a variety of uses from cleaners and lubricants to food, cosmetics, and coatings. In addition, they can be used as a tool for testing the fundamental tenets of statistical physics. Colloidal suspensions can be synthesized from a wide variety of materials, and in the form of monodisperse particles, which can self-assemble into highly ordered colloidal crystal structures. As such they can also be used as templates for the construction of highly ordered materials. Materials design of colloids has, to date, relied on entropic self-assembly, where crystals form as result of lower free energy due to a transition to order. Here, our goal is to develop a completely new method for materials fabrication using colloidal precursors, in which the self-assembly of the ordered colloidal structures is driven by a highly controllable, attractive interaction. This will greatly increase the range of potential structures that can be fabricated with colloidal particles. In this work, we demonstrate that colloidal suspensions can be crosslinked through highly specific biological crosslinking reactions. In particular, the molecules we use are protein-carbohydrate interactions derived from the immune system. This different driving force for self-assembly will yield different and novel suspensions structures. Because the biological interactions are heterotypic (A binding to B), this chemical system can be used to make binary alloys in which the two colloid subpopulations vary in some property - size, density, volume fraction, magnetic susceptibility, etc. An additional feature of these molecules which is unique - even within the realm of biological recognition - is that the molecules bind reversibly on reasonable time-scales, which will enable the suspension to sample different configurations, and allow us to manipulate and measure the size of the suspension dynamically. Because of the wide variety of structures that can be made from these novel colloids, and because the suspension structure can be altered dynamically, we believe this biocolloid system will yield a novel set of materials with many technological applications, including sensors (both biological and non-biological), optical filters and separation media.

Hammer, D. A.↗