Search NASA⌕ Search

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.

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

At least 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Reducing module soiling with scalable and robust photocatalytic coatings

The air-glass interface at the front of a photovoltaic (PV) module reflects approximately 4% of incident light, decreasing the potential power output of the module by the same amount. Today’s modules reduce this loss by adding a low-refractive-index (1.25-1.30) SiO2 coating to the sunward side of the module glass; this antireflection coating recovers approximately 3% of the 4% light that would otherwise be lost. While such antireflection coatings work very well on clean, new modules, they do not inhibit soiling—the accumulation of soilants such as dust, pollen, soot, or other foreign material—on the module glass, and soilants reflect and scatter incident light. An improved coating would serve provide not only an antireflection effect, but also an anti-soiling effect. The goal of this project was to develop such a coating and provide a path for it to be manufactured in the U.S. The project successfully designed and fabricated coatings that provided >3% transmittance gain compared to bare glass (matching the performance of commercial antireflection coatings) and displayed anti-soiling behavior in standard laboratory soiling effects. This was achieved by using a Swift Coat proprietary coating deposition technique, aerosol impact-driven assembly (AIDA), to control the porosity and thus refractive index of coatings of photocatalytic materials—such as TiO2—that would otherwise increase (instead of decrease) reflection. These combined antireflection/anti-soiling coatings passed PV industry standard module reliability tests as well as coating-specific abrasion tests, showing that they have the durability needed for decades in the field. Swift Coat scaled the AIDA hardware and deposition process to make mini-modules that were monitored for nearly two years during field tests administered by a third party, as well as demonstrated scaling to the widths of full-sized modules. The fielded mini-modules outperformed reference modules (with commercial antireflection coatings) in two locations, providing a 1% absolute average performance boost and larger increases during periods of heavier soiling. Swift Coat’s cost analysis indicated a coating manufacturing cost below the sales price of today’s antireflection coatings. More than five module manufacturers sampled and assessed the coatings, and three provided letters of support. The coating developed in this project increases the energy output of PV modules, thereby decreasing the cost per kilowatt-hour of solar energy generated. Cheaper solar electricity benefits the public by accelerating the transition to a stable, affordable, carbon-free energy economy. In addition, for select applications in which PV modules are highly visible—such as on residential rooftops—the coating provides an aesthetic benefit because it stays cleaner than today’s modules. Finally, Swift Coat and its prospective customers are U.S. companies, and successful commercialization of this technology will provide U.S. jobs and a secure solar supply chain.

14 SOLAR ENERGY↗
DOE OSTI2025

X-ray thermal diffuse scattering as a texture-robust temperature diagnostic for dynamically compressed solids

We present a model of x-ray thermal diffuse scattering (TDS) from a cubic polycrystal with an arbitrary crystallographic texture, based on the classic approach of Warren [B. E. Warren, Acta Crystallogr. 6, 803 (1953)]. We compare the predictions of our model with femtosecond x-ray diffraction patterns gathered from ambient and dynamically compressed rolled copper foils obtained at the High Energy Density instrument of the European X-Ray Free-Electron Laser facility and find that the texture-aware TDS model yields more accurate results than does the conventional powder model owed to Warren. Nevertheless, we further show: with sufficient angular detector coverage, the TDS signal is largely unchanged by sample orientation and in all cases strongly resembles the signal from a perfectly random powder; shot-to-shot fluctuations in the TDS signal resulting from grain-sampling statistics are at the percent level, in stark contrast to the fluctuations in the Bragg-peak intensities (which are over an order of magnitude greater); and TDS is largely unchanged even following texture evolution caused by compression-induced plastic deformation. We conclude that TDS is robust against texture variation, making it a flexible temperature diagnostic applicable just as well to off-the-shelf commercial foils as to ideal powders.

Crystal lattices↗
DOE OSTI2023

Physicochemical characterization of free troposphere and marine boundary layer ice-nucleating particles collected by aircraft in the eastern North Atlantic

Abstract. Atmospheric ice nucleation impacts the hydrological cycle and climate by modifying the radiative properties of clouds. To improve our predictive understanding of ice formation, ambient ice-nucleating particles (INPs) need to be collected and characterized. Measurements of INPs at lower latitudes in a remote marine region are scarce. The Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA) campaign, in the region of the Azores islands, provided the opportunity to collect particles in the marine boundary layer (MBL) and free troposphere (FT) by aircraft during the campaign's summer and winter intensive operation period. The particle population in samples collected was examined by scanning transmission X-ray microscopy with near-edge X-ray absorption fine structure spectroscopy. The identified INPs were analyzed by scanning electron microscopy with energy-dispersive X-ray analysis. We observed differences in the particle population characteristics in terms of particle diversity, mixing state, and organic volume fraction between seasons, mostly due to dry intrusion events during winter, as well as between the sampling locations of the MBL and FT. These differences are also reflected in the temperature and humidity conditions under which water uptake, immersion freezing (IMF), and deposition ice nucleation (DIN) proceed. Identified INPs reflect typical particle types within the particle population on the samples and include sea salt, sea salt with sulfates, and mineral dust, all associated with organic matter, as well as carbonaceous particles. IMF and DIN kinetics are analyzed with respect to heterogeneous ice nucleation rate coefficients, Jhet, and ice nucleation active site density, ns, as a function of the water criterion Δaw. DIN is also analyzed in terms of contact angles following classical nucleation theory. Derived MBL IMF kinetics agree with previous ACE-ENA ground-site INP measurements. FT particle samples show greater ice nucleation propensity compared to MBL particle samples. This study emphasizes that the types of INPs can vary seasonally and with altitude depending on sampling location, thereby showing different ice nucleation propensities, which is crucial information when representing mixed-phase cloud and cirrus cloud microphysics in models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Machine learning the metastable phase diagram of covalently bonded carbon

Abstract Conventional phase diagram generation involves experimentation to provide an initial estimate of the set of thermodynamically accessible phases and their boundaries, followed by use of phenomenological models to interpolate between the available experimental data points and extrapolate to experimentally inaccessible regions. Such an approach, combined with high throughput first-principles calculations and data-mining techniques, has led to exhaustive thermodynamic databases (e.g. compatible with the CALPHAD method), albeit focused on the reduced set of phases observed at distinct thermodynamic equilibria. In contrast, materials during their synthesis, operation, or processing, may not reach their thermodynamic equilibrium state but, instead, remain trapped in a local (metastable) free energy minimum, which may exhibit desirable properties. Here, we introduce an automated workflow that integrates first-principles physics and atomistic simulations with machine learning (ML), and high-performance computing to allow rapid exploration of the metastable phases to construct “metastable” phase diagrams for materials far-from-equilibrium. Using carbon as a prototypical system, we demonstrate automated metastable phase diagram construction to map hundreds of metastable states ranging from near equilibrium to far-from-equilibrium (400 meV/atom). We incorporate the free energy calculations into a neural-network-based learning of the equations of state that allows for efficient construction of metastable phase diagrams. We use the metastable phase diagram and identify domains of relative stability and synthesizability of metastable materials. High temperature high pressure experiments using a diamond anvil cell on graphite sample coupled with high-resolution transmission electron microscopy (HRTEM) confirm our metastable phase predictions. In particular, we identify the previously ambiguous structure of n -diamond as a cubic-analog of diaphite-like lonsdaelite phase.

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

Electron irradiation effects on the optical properties of Hf- and Zn-doped β -Ga 2 O 3

Optical and electrical properties of Hf- and Zn-doped β-Ga 2 O 3 samples, which are n-type and insulating, respectively, were altered via high-energy electron irradiation at 2.5 or 0.5 MeV. The β-Ga 2 O 3 :Hf samples irradiated with 2.5 MeV electrons experienced a color change from blue to yellow and a large drop in conductivity, attributed to the creation of gallium vacancies, which compensate donors. This irradiation resulted in the absence of free carrier absorption and the presence of Cr 3+ photoluminescence (PL). PL mapping prior to irradiation revealed optically active ZnO precipitates that formed during the growth of β-Ga 2 O 3 :Zn. These precipitates have a 384 nm (3.23 eV) stacking fault emission in the core; in the outer shell of the precipitate, the PL blue-shifts to 377 nm (3.29 eV) and a broad defect band is observed. After 0.5 MeV electron irradiation, the defect band broadened and increased in intensity. The blue PL band (435 nm) of β-Ga 2 O 3 was enhanced for both Hf- and Zn-doped samples irradiated with 0.5 MeV. This enhancement is correlated with an increase in oxygen vacancies.

36 MATERIALS SCIENCE↗
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↗
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↗
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↗