SEARCH · Search NASA

Results for “free energy 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Triton Field Trials (TFiT) underwater noise - University of New Hampshire Living Bridge turbine Processed Data

In July 2021, a commercial-off-the-shelf hydrophone was deployed in a free-drifting configuration to measure underwater acoustic emissions and characterize a 25 kW-rated tidal turbine at the University of New Hampshire's Living Bridge Project in Portsmouth, New Hampshire. Sampling methods and analysis were performed in alignment with the recently published IEC 62600-40 Technical Specification for acoustic characterization of marine energy converters. Results from this study indicate acoustic emissions from the turbine were below ambient sound levels and therefore did not have a significant impact on the underwater noise levels of the project site. As a component of Pacific Northwest National Laboratory's Triton Field Trials (TFiT) described in a paper published in a Special Issue of Journal of Marine Energy Science and Engineering, this study provides a valuable use case for the IEC 62600-40 Technical Specification framework and further recommendations for cost-effective technologies and methods for measuring underwater noise at future current energy converter project sites. The paper can be accessed in the link bellow.

16 TIDAL AND WAVE POWER↗

DOE OSTI2023

Microstructure Alignment Effects from Engineered Cooling during Additive Manufacturing of Alnico Magnets from Prealloyed Powder

Alnico is a rare-earth free; permanent magnet (PM) which maintains its saturation magnetization while at temperatures up to 550°C but lacks the coercivity and energy product needed for a wide range of demanding applications. With an optimal set of processing and manufacturing techniques, Alnico could replace Dy-free Nd-Fe-B magnets in certain applications if the microstructure is textured in a $\langle$001$\rangle$ direction. Additive manufacturing (AM) of near-net-shape magnets may provide such texture. To investigate the possibility, directed-energy deposition (DED) in conjunction with an actively cooled substrate was used to build samples with a compositionally modified, gas-atomized, Alnico 8 alloy. The resulting samples were solutionized and quenched, followed by magnetic annealing and a heat treatment. Magnetic properties and microstructures are compared to previous work.

36 MATERIALS SCIENCE↗

NASA NTRS1989

A study of the effects of Rotating Frame Turbulence (RFT) on helicopter flight mechanics

The turbulence actually experienced by a helicopter blade-element significantly differs from the space-fixed free atmospheric turbulence. The turbulence in the rotor disk requires a rotationally sampled description in a rotating frame of reference. It is referred to as the rotating frame turbulence or RFT which exhibits a striking phenomenon. The RFT spectral density versus frequency shows high peak values at 1P,2P, or 3P, frequencies. The energy increase at these peaks is balanced by an energy decrease primarily at the lower-than-1P frequency range. Particularly for low altitude flight regimes of pure helicopters, such as the nap-of-the-earth maneuvers, the conventional space-fixed description of turbulence is not a good approximation, since the turbulence scale length can have values comparable to the rotor radius. Accordingly the flight mechanics characteristics with RFT description are compared with those based on the conventional space-fixed turbulence description. The results demonstrate that the RFT qualitatively and quantitatively affects the prediction of helicopter flight mechanics characteristics in turbulence. Such comparisons should play an important role in the new development of handling qualities specifications for helicopters.

Schrage, D. P.↗

NASA NTRS1996

Investigation of Polymer Liquid Crystals

The positron annihilation lifetime spectroscopy (PALS) using a low energy flux generator may provide a reasonably accurate technique for measuring molecular weights of linear polymers and characterization of thin polyimide films in terms of their dielectric constants and hydrophobity etc. Among the tested samples are glassy poly arylene Ether Ketone films, epoxy and other polyimide films. One of the proposed techniques relates the free volume cell size (V(sub f)) with sample molecular weight (M) in a manner remarkably similar to that obtained by Mark Houwink (M-H) between the inherent viscosity (eta) and molecular wieght of polymer solution. The PALS has also demonstrated that free-volume cell size in thermoset is a versatile, useful parameter that relates directly to the polymer segmental molecular weight, the cross-link density, and the coefficient of thermal expansion. Thus, a determination of free volume cell size provides a viable basis for complete microstructural characterization of thermoset polyimides and also gives direct information about the cross-link density and coefficient of expansion of the test samples. Seven areas of the research conducted are reported here.

Han, Kwang S.↗

NASA NTRS2004

Towards Measurement of the Time-resolved Heat Release of Protein Conformation Dynamics

We present a way to observe time-resolved heat release using a laminar flow diffusional mixer coupled with a highly sensitive infrared camera which measures the temperature change of the solvent. There are significant benefits to the use of laminar flow mixers for time-resolved calorimetry: (1) The thermal signal can be made position and time- stationary to allow for signal integration; (2) Extremely small volumes (nl/s) of sample are required for a measurement; (3) The same mixing environment can be observed spectroscopically to obtain state occupation information; (4) The mixer allows one to do out of equilibrium dynamic studies. The hope is that these measurements will allow us probe the non-equilibrium thermodynamics as a protein moves along a free energy trajectory from one state to another.

Puchalla, Jason↗

DOE OSTI2026

Metal–bicarbonate ion pairing in alkaline aqueous solutions from multilevel embedded correlated wavefunction theory and molecular dynamics

In this work, we examine ion-pairing mechanisms of bicarbonates in alkaline aqueous solutions with the divalent metal ions most abundantly present in seawater, namely, Ca 2+ and Mg 2+ . We employ a rare-event enhanced sampling approach within first-principles molecular dynamics to explore regions of phase space spanning solvent-shared to contact ion pairs. Second-order Møller–Plesset perturbation theory (MP2) corrections are subsequently applied in an embedding framework (EMB) to refine the electronic structure of stationary states and associated reaction barriers along the free-energy profiles while retaining the extended solvent effects at the density functional theory (DFT) level. Ca 2+ –HCO 3 − was previously hypothesized to exist in a solvent-shared ion pair (SSHIP) by DFT studies with an endergonic contact ion pair (CIP) formation; however, our EMB-MP2 refinement of the DFT ion-pairing pathways reveals that Ca 2+ and HCO 3 − form a virtually barrier-free CIP in alkaline solutions, with even more energetic ease than the widely studied Ca–CO 3 ion pair. We find qualitative agreement between DFT and EMB-MP2 for Mg 2+ —unlike Ca 2+ , Mg 2+ refuses to shed its strong hydration shell, thereby preferring a SSHIP state with a significant activation barrier to crossover to the CIP forms—a trait reminiscent of ion pairing in Mg–CO 3 and closely related to the kinetic limitations underlying the famous subject of the dolomite problem. Our study highlights the importance of improved electronic structure descriptions of liquids, modeled as a condensed phase of matter lacking in long-range crystalline order. It also strongly suggests that Ca 2+ –HCO 3 − CIPs are likely precursors involved in prenucleation of CaCO 3 mineral formation in seawater.

Sharma, Vidushi [Princeton Plasma Physics Laborato↗

DOE OSTI2026

The Role of Subcloud Mesoscale Convergence in Sculpting Convective Updraft Width and Depth

The initiation of deep moist convection is governed in part by the horizontal width of updrafts near cloud base, which limits the deleterious effects of entrainment-driven dilution on buoyant thermals ascending through the free troposphere. However, the factors controlling cloud-base updraft width, which in turn dictates cloud depth, are not well understood. We track the evolving three-dimensional structure of the mesoscale subcloud forcing for vertical motion and near-cloud thermodynamic ingredients within a high-resolution ensemble of simulations of seven realistic daytime orographic convection initiation events to determine their relative roles in controlling cloud width and depth. Statistical analysis of approximately 5000 cloudy updraft samples indicates that the most important contributors to the width of cloudy updrafts across the ensemble are the depth and magnitude of the subcloud mesoscale ascent. However, the depth achieved by clouds is more consistently predicted by the near-cloud ambient relative humidity within the lower to middle free troposphere and convective available potential energy. Therefore, although the width of cloudy updrafts may be partly set at low levels by the mesoscale vertical mass and moisture flux, the likelihood of deep moist convection is governed by the generation of positive buoyancy within cumulus thermals and entrainment-driven dilution that reduces it. The persistence of the low-level mesoscale vertical forcing locally consolidates and vertically transports boundary layer moisture, helping to reduce updraft dilution. However, these factors vary in relative impacts on cloudy updrafts across individual cases, indicating multiple pathways for deep convection initiation.

Convective storms↗

DOE OSTI2022

A Study on Surface Texture and Wettability of Femtosecond Laser Treated Aluminum Alloys

This research aims to study the effects of femtosecond laser treatment on surface texture and wettability properties. The first portion of this research studies 2024-T3 aluminum alloy concerning surface texture (e.g., roughness), wettability, and surface free energy. SEM images help aid the researchers in detecting minuscule changes in surface texture. Profilometry extracts topographical data from the surface. A contact angle goniometer measures contact angle, which follows the sessile drop method. Owens, Wendt, Rabel, and Kaelble's method calculates the surface free energy in conjunction with the contact angle. Femtosecond laser modification produces relatively uniform surface morphology at the macroscale. When looking at the microscale, surface morphology differs depending on the interplay between material melting and self-assembly, where the surface becomes rougher after laser treatment. Contact angle results show that all but the control sample and sample 7 exhibit hydrophobic behavior, i.e., contact angle >90°, with Sample 1 having the most significant angle of 115°. The effects of laser parameters on surface energy is apparent. More work is needed to establish correlations among these three quantities: contact angle, surface energy, and laser parameters. The second portion of this research studies aluminum alloy 6061 as the subject material with the primary goal of creating a superhydrophobic surface. A central composite design creates the laser surface texturing design of experiments. SEM images help aid the researchers in detecting minuscule changes in surface texture. A 3D optical profiler extracts topographical data from the surface. The sessile drop method incorporates a contact angle goniometer to measure the contact angle. The response surface method builds a second-order polynomial model for the contact angle and obtains optimized parameters to maximize the contact angle. This research shows that laser surface texturing can generate a wide range of surface profiles and roughness values with geometric features ranging from hundreds of μm to submicron. All three laser parameters (pulse energy, pulse duration, repetition rate) affect surface roughness and contact angle to some degree. To quantify the relationship between the contact angle, pulse energy, and pulse repetition rate, a response surface model for Al 6061 is identified and used to find the optimal conditions of E=214 μJ, tp=10 ps, and fp=2427 Hz with a predicted maximum contact angle of 161°. A confirmation experiment produces a contact angle of 168°, in good agreement with the predicted value.

Angell, Dakota↗

NASA NTRS2003

Self-Propagating Frontal Polymerization in Water at Ambient Pressure

Advances in polymer chemistry have led to the development of monomers and initiation agents that enable propagating free-radical polymerization fronts to exist. These fronts are driven by the exothermicity of the polymerization reaction and the transport of heat from the polymerized product to the reactant monomer/solvent/initiator solution. The thermal energy transported to the reactant solution causes the initiator to decompose, yielding free radicals, which start the free radical polymerization process as discussed in recent reviews. The use of polymerization processes based on propagating fronts has numerous applications. Perhaps the most important of these is that it enables rapid curing of polymers without external heating since the polymerization process itself provides the high temperatures necessary to initiate and sustain polymerization. This process also enables more uniform curing of arbitrarily thick samples since it does not rely on heat transfer from an external source, which will necessarily cause the temperature history of the sample to vary with distance from the surface according to a diffusion-like process. Frontal polymerization also enables filling and sealing of structures having cavities of arbitrary shape without having to externally heat the structure. Water at atmospheric pressure is most convenient solvent to employ and the most important for practical applications (because of the cost and environmental issues associated with DMSO and other solvents). Nevertheless, to our knowledge, steady, self-propagating polymerization fronts have not been reported in water at atmospheric pressure. Currently, polymerization fronts require a high boiling point solvent (either water at high pressures or an alternative solvent such as dimethyl sulfoxide (DMSO) (boiling point 189 C at atmospheric pressure.) Early work on frontal polymerization, employed pressures up to 5000 atm in order to avoid boiling of the monomer/solvent/initiator solution. High boiling point solutions are needed because in order to produce a propagating front, a high front temperature is needed to produce sufficiently rapid decomposition of the free radical initiator and subsequent free radical polymerization and heat release at a rate faster than heat losses remove thermal energy from the system. (While the conduction heat loss rate increases linearly with temperature, the free radical initiator decomposition is a high activation energy process whose rate increases much more rapidly than linearly with temperature, thus as the temperature decreases, the ratio of heat loss to heat generation increases, eventually leading to extinction of the front if the temperature is too low.) In order to obtain atmospheric pressure frontal polymerization in water, it is necessary to identify a monomer/initiator combination that is water soluble and will not extinguish even when the peak temperature (T*) is less than 100 C. In this work acrylic acid (AA) was chosen as the monomer because is it one of the most reactive monomers and can polymerize readily at low temperatures even without initiators. Ammonium persulfate (AP) was chosen as the initiator because it decomposes readily at low temperatures, produces relatively few bubbles and is commercially available. The propagation rates and extinction conditions of the fronts are studied for a range of AA and AP concentrations. Small amounts of fumed silica powder (Cab-o-sil, Cabot Corporation) were added to the solutions to inhibit buoyancy induced convection in the solutions; future studies will investigate the effects of buoyant convection within the solutions.

Olten, Nesrin↗

DOE OSTI2021

Probing remdesivir nucleotide analogue insertion to SARS-CoV-2 RNA dependent RNA polymerase in viral replication

Remdesivir (RDV) prodrug can be metabolized into a triphosphate form nucleotide analogue (RDV-TP) to bind and insert into the active site of viral RNA dependent RNA polymerase (RdRp) to further interfere with viral genome replication. Here in this work, we computationally studied how RDV-TP binds and inserts to the SARS-CoV-2 RdRp active site, in comparison with natural nucleotide substrate adenosine triphosphate (ATP). To do that, we first constructed atomic structural models of an initial binding complex (active site open) and a substrate insertion complex (active site closed), based on high-resolution cryo-EM structures determined recently for SARS-CoV-2 RdRp or non-structural protein (nsp) 12, in complex with accessory protein factors nsp7 and nsp8. By conducting all-atom molecular dynamics simulation with umbrella sampling strategies on the nucleotide insertion between the open and closed state RdRp complexes, our studies show that RDV-TP can initially bind in a comparatively stabilized state to the viral RdRp active site, as it primarily forms base stacking with the template uracil nucleotide (nt +1), which under freely fluctuations supports a low free energy barrier of the RDV-TP insertion (~1.5 kcal mol -1 ). In comparison, the corresponding natural substrate ATP binds initially to the RdRp active site in Watson–Crick base pairing with the template nt, and inserts into the active site with a medium low free energy barrier (~2.6 kcal mol -1 ), when the fluctuations of the template nt are well quenched. The simulations also show that the initial base stacking of RDV-TP with the template can be specifically stabilized by motif C-S759, S682 (near motif B) with the base, and motif G-K500 with the template backbone. Although the RDV-TP insertion can be hindered by motif F-R555/R553 interaction with the triphosphate, the ATP insertion seems to be facilitated by such interactions. The inserted RDV-TP and ATP can be further distinguished by specific sugar interaction with motif B-T687 and motif A-D623, respectively.

36 MATERIALS SCIENCE↗

DOE OSTI2021

The Role of Hydrophobic Nodes in the Dynamics of Class A β-Lactamases

Class A β-lactamases are known for being able to rapidly gain broad spectrum catalytic efficiency against most β-lactamase inhibitor combinations as a result of elusively minor point mutations. The evolution in class A β-lactamases occurs through optimisation of their dynamic phenotypes at different timescales. At long-timescales, certain conformations are more catalytically permissive than others while at the short timescales, fine-grained optimisation of free energy barriers can improve efficiency in ligand processing by the active site. Free energy barriers, which define all coordinated movements, depend on the flexibility of the secondary structural elements. The most highly conserved residues in class A β-lactamases are hydrophobic nodes that stabilize the core. To assess how the stable hydrophobic core is linked to the structural dynamics of the active site, we carried out adaptively sampled molecular dynamics (MD) simulations in four representative class A β-lactamases (KPC-2, SME-1, TEM-1, and SHV-1). Using Markov State Models (MSM) and unsupervised deep learning, we show that the dynamics of the hydrophobic nodes is used as a metastable relay of kinetic information within the core and is coupled with the catalytically permissive conformation of the active site environment. Finally, our results collectively demonstrate that the class A enzymes described here, share several important dynamic similarities and the hydrophobic nodes comprise of an informative set of dynamic variables in representative class A β-lactamases.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Molecular-Resolution Electron Imaging of Defects and Dynamics at the Ice-Water Interface

Water crystallization into hexagonal ice (type I h ) is one of the most critical processes relevant to the Earth’s environment and human activities. However, despite recent breakthroughs in imaging non-equilibrium condensed ice structures, the ice-water interface has never been imaged at a molecular resolution. This is primarily due to the low stability of the hydrogen bonds in ice under high-resolution microscopy conditions and a lack of methods to prepare compatible samples. Here, this presentation describes the first molecular-resolution imaging of ice crystallized from liquid water and the ice-water interface using high-resolution transmission electron microscopy (HRTEM). By encapsulating deionized (DI) water between two amorphous carbon (a-C) TEM grids and subsequently freezing it with liquid N 2 on a cryo sample stage, we generated two types of ice: non-equilibrium, condensed ice from the atmosphere and encapsulated ice from the DI water (Fig. 1A). Condensed ice usually shows irregular, spherulitic shapes (Fig. 1B). Selected area electron diffraction (SAED) shows that they are a mixture of cubic and hexagonal crystals (Fig. 1C). On the contrary, encapsulated ice forms thin films that contain large-area single-crystalline regions of hexagonal ice oriented along the [0001] zone axis (Fig. 1D). Differential electron energy-loss spectroscopy (EELS) confirmed the high purity of the encapsulated ice free from organic contaminations that are common in other encapsulation methods for HRTEM such as graphene liquid cells. These single-crystalline areas are robust under the electron beam up to ~100 e/Å 2 s. Aberration-corrected HRTEM imaging in these areas achieved a line resolution of ~1.3 Å (Fig. 1E and F). This platform allows us to study near-equilibrium ice structures and dynamics at an unprecedented spatial resolution (Fig. 2). For example, we discovered subdomain-rich regions near the defective crystal edges despite the structure appearing single-crystalline according to diffraction criteria. These subdomains connect via low-angle grain boundaries with flat energy landscapes as a function of tilt angles (according to simulations), showing the high tolerance of ice to defect structures. When we tuned the sample temperature and electron flux rate, we observed radiolysis-controlled bubble generation and dissolution in ice single crystals near a steady state of bubble dynamics. Furthermore, rich beam-induced melting and recrystallization dynamics were observed at the ice-water interface with lattice resolution. These data represent the first observation of the ice-water phase transformation at the sub-nanometer level. In summary, the methods developed in this work enabled molecular-resolution observations of ice and the ice-water interface and shed light on the microstructures and phase transformation pathways. Finer control on the temperature, electron irradiation profile, and imaging detector could eventually lead to real-time observation of ice nucleation in water and address long-standing questions in the nucleation pathways.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2024

A Generation-Storage Coordination Dispatch Strategy for Power System Based on Causal Reinforcement Learning

In the backdrop of global energy transformation, power systems integrating high proportions of renewable energy sources are facing unprecedented challenges in operational stability and dispatch efficiency. To address these challenges, this study introduces a generation-storage coordination real-time dispatch strategy based on Causal Power System Dynamic Reinforcement Learning (CPSDRL). Diverging from traditional reinforcement learning approaches, CPSDRL innovatively incorporates causal inference within the state prediction model - the crux of model-based reinforcement learning - thereby establishing the Power Causal Dynamic Model (PCDM). Assisted by the prior knowledge of power systems, the model significantly enhances prediction accuracy and reliability through a two-stage training process. Utilizing PCDM, this study further applies a direct policy search algorithm to optimize the real-time dispatch strategy. Experimental results indicate that the proposed method improves the stability of generation-storage coordination real-time dispatch and exhibits competitive advantages in sample efficiency and computational speed, compared to traditional model-based and model-free reinforcement learning algorithms. This method is expected to enhance the practicality and adaptability of causal reinforcement learning techniques in power system scheduling and control.

causal reinforcement learning↗

DOE OSTI2021

Dark Energy Survey Year 3 Results: Galaxy mock catalogs for BAO analysis

The calibration and validation of scientific analysis in simulations is a fundamental tool to ensure unbiased and robust results in observational cosmology. In particular, mock galaxy catalogs are a crucial resource to achieve these goals in the measurement of baryon acoustic oscillation (BAO) in the clustering of galaxies. Here we present a set of 1952 galaxy mock catalogs designed to mimic the Dark Energy Survey Year 3 BAO sample over its full photometric redshift range 0.6 < zphoto < 1.1. The mocks are based upon 488 ICE-COLA fast N-body simulations of full-sky light cones and were created by populating halos with galaxies, using a hybrid halo occupation distribution – halo abundance matching model. This model has ten free parameters, which were determined, for the first time, using an automatic likelihood minimization procedure. We also introduced a novel technique to assign photometric redshift for simulated galaxies, following a two-dimensional probability distribution with VIMOS Public Extragalactic Redshift Survey data. The calibration was designed to match the observed abundance of galaxies as a function of photometric redshift, the distribution of photometric redshift errors, and the clustering amplitude on scales smaller than those used for BAO measurements. An exhaustive analysis was done to ensure that the mocks reproduce the input properties. Finally, mocks were tested by comparing the angular correlation function w(θ), angular power spectrum Cℓ, and projected clustering ξp(r⊥) to theoretical predictions and data. The impact of volume replication in the estimate of the covariance is also investigated. The success in accurately reproducing the photometric redshift uncertainties and the galaxy clustering as a function of redshift render this mock creation pipeline as a benchmark for future analyses of photometric galaxy surveys.Key words: catalogs / large-scale structure of Universe / galaxies: distances and redshifts / Galaxy: halo / methods: numerical

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Active learning of ternary alloy structures and energies

Abstract Machine learning models with uncertainty quantification have recently emerged as attractive tools to accelerate the navigation of catalyst design spaces in a data-efficient manner. Here, we combine active learning with a dropout graph convolutional network (dGCN) as a surrogate model to explore the complex materials space of high-entropy alloys (HEAs). We train the dGCN on the formation energies of disordered binary alloy structures in the Pd-Pt-Sn ternary alloy system and improve predictions on ternary structures by performing reduced optimization of the formation free energy, the target property that determines HEA stability, over ensembles of ternary structures constructed based on two coordinate systems: (a) a physics-informed ternary composition space, and (b) data-driven coordinates discovered by the Diffusion Maps manifold learning scheme. Both reduced optimization techniques improve predictions of the formation free energy in the ternary alloy space with a significantly reduced number of DFT calculations compared to a high-fidelity model. The physics-based scheme converges to the target property in a manner akin to a depth-first strategy, whereas the data-driven scheme appears more akin to a breadth-first approach. Both sampling schemes, coupled with our acquisition function, successfully exploit a database of DFT-calculated binary alloy structures and energies, augmented with a relatively small number of ternary alloy calculations, to identify stable ternary HEA compositions and structures. This generalized framework can be extended to incorporate more complex bulk and surface structural motifs, and the results demonstrate that significant dimensionality reduction is possible in thermodynamic sampling problems when suitable active learning schemes are employed.

Chemistry↗

DOE OSTI2022

Supervised learning and the finite-temperature string method for computing committor functions and reaction rates

A central object in the computational studies of rare events is the committor function. Though costly to compute, the committor function encodes complete mechanistic information of the processes involving rare events, including reaction rates and transition-state ensembles. Under the framework of transition path theory, Rotskoff et al. [Proceedings of the 2nd Mathematical and Scientific Machine Learning Conference, Proceedings of Machine Learning Research (PLMR, 2022), Vol. 145, pp. 757–780] proposes an algorithm where a feedback loop couples a neural network that models the committor function with importance sampling, mainly umbrella sampling, which collects data needed for adaptive training. Here, in this work, we show additional modifications are needed to improve the accuracy of the algorithm. The first modification adds elements of supervised learning, which allows the neural network to improve its prediction by fitting to sample-mean estimates of committor values obtained from short molecular dynamics trajectories. The second modification replaces the committor-based umbrella sampling with the finite-temperature string (FTS) method, which enables homogeneous sampling in regions where transition pathways are located. We test our modifications on low-dimensional systems with non-convex potential energy where reference solutions can be found via analytical or finite element methods, and show how combining supervised learning and the FTS method yields accurate computation of committor functions and reaction rates. We also provide an error analysis for algorithms that use the FTS method, using which reaction rates can be accurately estimated during training with a small number of samples. The methods are then applied to a molecular system in which no reference solution is known, where accurate computations of committor functions and reaction rates can still be obtained.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Solvent reaction coordinate for an S N 2 reaction

We study the prototypical S N 2 reaction Cl – + CH 3 Cl → CH 3 Cl + Cl – in water using QM/MM computer simulations with transition path sampling and inertial likelihood maximization. We have identified a new solvent coordinate that improves upon the original atom-exchange coordinate used in the classic analysis by Chandrasekhar, Smith, and Jorgensen. The new solvent coordinate quantifies instantaneous solvent induced polarization relative to equilibrium charge density at each point along the reaction pathway. In this paper, we will present free energy surfaces, transmission coefficients, and committor tests associated with this collective coordinate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Kekulé valence bond order in the honeycomb lattice optical Su-Schrieffer-Heeger model and its relevance to graphene

We perform sign-problem-free determinant quantum Monte Carlo simulations of the optical Su- Schrieffer-Heeger model on a half-filled honeycomb lattice. In particular, we investigate the model’s semi-metal (SM) to Kekulé Valence Bond Solid (KVBS) phase transition at zero and finite temper- atures as a function of phonon energy and interaction strength. Using hybrid Monte Carlo sampling methods we can simulate the model near the adiabatic regime, allowing us to access regions of parameter space relevant to graphene. Our simulations suggest that the SM-KVBS transition is weakly first-order at all temperatures, with graphene situated close to the phase boundary in the SM region of the phase diagram. Furthermore, our results highlight the important role bond-stretching phonon modes play in the formation of KVBS order in strained graphene-derived systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗