SEARCH · Engineering Papers

Results for “coordinated energy system”

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 757 records · Page 42

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Electrochemical Random-Access Memory: Progress, Perspectives, and Opportunities

Non-von Neumann computing using neuromorphic systems based on analogue synaptic and neuronal elements has emerged as a potential solution to tackle the growing need for more efficient data processing, but progress toward practical systems has been stymied due to a lack of materials and devices with the appropriate attributes. Recently, solid state electrochemical ion-insertion, also known as electrochemical random access memory (ECRAM) has emerged as a promising approach to realize the needed device characteristics. ECRAM is a three terminal device that operates by tuning electronic conductance in functional materials through solid-state electrochemical redox reactions. This mechanism can be considered as a gate-controlled bulk modulation of dopants and/or phases in the channel. Early work demonstrating that ECRAM can achieve nearly ideal analogue synaptic characteristics has sparked tremendous interest in this approach. More recently, the realization that electrochemical ion insertion can be used to tune the electronic properties of many types of materials including transition metal oxides, layered two-dimensional materials, organic and coordination polymers, and that the changes in conductance can span orders of magnitude has further attracted interest in ECRAM as the basis for analogue synaptic elements for inference accelerators as well as for dynamical devices that can emulate a wide range of neuronal characteristics for implementation in analogue spiking neural networks. At its core, ECRAM shares many fundamental aspects with rechargeable batteries, where ion insertion materials are used extensively for their ability to reversibly store charge and energy. Computing applications, however, present drastically different requirements: systems will require many millions of devices, scaled down to tens of nanometers, all while achieving reliable electronic-state tuning at scaled-up rates and endurances, and with minimal energy dissipation and noise. Further, in this review, we discuss the history, basic concepts, recent progress, as well as the challenges and opportunities for different types of ECRAM, broadly grouped by their primary mobile ionic charge carrier, including Li, protons, and oxygen vacancies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Synopsis of NREL's Automated Mobility District (AMD) Research Program and Associated Publications

An automated mobility district (AMD) envisions a system of integrated mobility options that serves major activity centers such as campuses, central business districts, and large medical facilities. The National Renewable Energy Laboratory (NREL) has been investigating the implementation prospects for fully automated passenger transport systems that are deployed to operate within dense urban settings. This document provides a synopsis of findings revealed over the last three phases of work, which have yielded insights into the creation and management of AMDs anticipated to use automated vehicle (AV) technology over the next decade. Phase I and Phase II tracked the deployment and lessons learned from 10 early-stage demonstrations of automated shuttle deployments, and their associated insights into the challenges for automated driving systems to achieve safe operations within district-scale deployments. Phase III began in-depth investigations of critical subsystem components, as automation, electrification, and on-demand service continue to converge within initial AMD operations. These directed studies focus on elements of electrification, curbfront/station management, the role of infrastructure sensing, and overall integration of AMD safety management in central, simultaneous coordination of multiple AMD fleets. Future research in AMDs includes systems engineering methodology (more frequently referred to as "digital twins") for planning, design, testing, and ongoing operation of AMDs; location (or co-location) of management functions; and human supervision and passenger communications for safety and security in unattended vehicles. The synopsis references the foundational research products (papers and presentations) that have been published through conference proceedings, journal articles, and NREL reports.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2023

Skyrme–Hartree–Fock–Bogoliubov mass models on a 3D mesh: IIb. Fission properties of BSkG2

Large-scale models of nuclear structure are currently the only way to provide consistent datasets for the many properties of thousands of exotic nuclei that are required by nucleosynthesis simulations. In [W. Ryssens et al., Eur. Phys. J. A 58, 246 (2022)], we recently presented the new BSkG2 model based on an energy density functional of the Skyrme type. Relying on a flexible three-dimensional coordinate representation of the nucleus, the model takes into account both triaxial deformation and time-reversal symmetry breaking. BSkG2 achieves a state-of-the-art global description of nuclear ground state (g.s.) properties and reproduces in particular the known masses with a root-mean-square (rms) deviation of 678 keV. Moving beyond g.s. properties, the model also reproduces all empirical values for the primary and secondary barriers as well as isomer excitation energies of actinide nuclei with rms deviations below 500 keV, i.e. with unprecedented accuracy. Here we discuss in detail the extension of our framework to the calculation of the fission barriers of 45 actinide nuclei, including odd-mass and odd-odd systems. We focus in particular on the impact of symmetry breaking which is key to the accuracy of the model: we allow systematically for axial, reflection and time-reversal symmetry breaking. The effect of the latter on the fission properties of odd-mass and odd-odd nuclei is small, but we find that allowing for shapes with triaxial or octupole deformation, as well as shapes with both, is crucial to achieving this accuracy. Here, the numerical accuracy of our coordinate space approach, the variety of nuclear configurations explored and the simultaneous successful description of fission properties and known masses makes BSkG2 the tool of choice for the large-scale study of nuclear structure.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Electrification Futures Study: Power Systems Operation with Newly Electrified and Flexible Loads

The Electrification Futures Study (EFS) is a multi-year study designed to analyze the impacts of widespread electrification in the U.S. energy system. In this webinar we present the detailed power systems operational analysis from the final report in the EFS series. For this analysis, we simulate multiple 2050 power systems for the conterminous United States to assess how variations in the magnitude and shape of electricity demand driven by electrification, and the extent of load participation to more-actively provide grid services, might impact the hourly operation, operational costs, and emissions of various power systems in 2050. The impacts of electrification and demand flexibility are overlayed across systems with significantly greater penetrations of variable renewable energy than today. Overall, we find that the high electrification scenarios envisioned in the EFS with significant VRE penetration (66% of annual national generation) can be operated to meet future increased levels of electrified demand. We also find that demand-side flexibility can enhance operational efficiency and reduce overall annual production costs by $5-$10 billion. The complementary relationship between flexible electric vehicle charging loads and solar generation is particularly pronounced, but in the absence of demand-side flexibility electrification can lead to increased wind and solar curtailment. By helping to reduce renewable curtailment, flexibility can also reduce power sector CO2 emissions. The analysis highlights the value of increased integration and coordination of demand- and supply-side resources in future electric system planning and operations—particularly under high electrification futures.

43 PARTICLE ACCELERATORS↗

DOE OSTI2023

Electrolytes with Solvating Inner Sheath Engineering for Practical Na–S Batteries

Sodium–sulfur (Na–S) batteries with durable Na-metal stability, shuttle-free cyclability, and long lifespan are promising to large-scale energy storages. However, meeting these stringent requirements poses huge challenges with the existing electrolytes. Herein, a localized saturated electrolyte (LSE) is proposed with 2-methyltetrahydrofuran (MeTHF) as an inner sheath solvent, which represents a new category of electrolyte for Na–S system. Unlike the traditional high concentration electrolytes, the LSE is realized with a low salt-to-solvent ratio and low diluent-to-solvent ratio, which pushes the limit of localized high concentration electrolyte (LHCE). The appropriate molecular structure and solvation ability of MeTHF regulate a saturated inner sheath, which features a reinforced coordination of Na + to anions, enlarged Na + -solvent distance, and weakened anion-diluent interaction. Such electrolyte configuration is found to be the key to build a sustainable interphase and a quasi-solid–solid sulfur redox process, making a dendrite-inhibited and shuttle-free Na–S battery possible. With this electrolyte, pouch cells with decent cycling performance under rather demanding conditions are demonstrated.

25 ENERGY STORAGE↗

DOE OSTI2021

Bacterial xylan utilization regulons: systems for coupling depolymerization of methylglucuronoxylans with assimilation and metabolism

Abstract Bioconversion of lignocellulosic resources offers an economically promising path to renewable energy. Technological challenges to achieving bioconversion include the development of cost-effective processes that render the cellulose and hemicellulose components of these resources to fermentable hexoses and pentoses. Natural bioprocessing of the hemicellulose fraction of lignocellulosic biomass requires depolymerization of methylglucuronoxylans. This requires secretion of endoxylanases that release xylooligosaccharides and aldouronates. Physiological, biochemical, and genetic studies with selected bacteria support a process in which a cell-anchored multimodular GH10 endoxylanase catalyzes release of the hydrolysis products, aldotetrauronate, xylotriose, and xylobiose, which are directly assimilated and metabolized. Gene clusters encoding intracellular enzymes, including α-glucuronidase, endoxylanase, β-xylosidase, ABC transporter proteins, and transcriptional regulators, are coordinately responsive to substrate induction or repression. The rapid rates of glucuronoxylan utilization and microbial growth, along with the absence of detectable products of depolymerization in the medium, indicate that assimilation and depolymerization are coupled processes. Genomic comparisons provide evidence that such systems occur in xylanolytic species in several genera, including Clostridium, Geobacillus, Paenibacillus, and Thermotoga. These systems offer promise, either in their native configurations or through gene transfer to other organisms, to develop biocatalysts for efficient production of fuels and chemicals from the hemicellulose fractions of lignocellulosic resources.

Biotechnology & Applied Microbiology↗

DOE OSTI2020

Nonadiabatic Dynamics of Photoexcited cis-Stilbene Using Ab Initio Multiple Spawning

The photochemistry of cis-stilbene proceeds through two pathways: cis-trans isomerization and ring closure to 4a,4b-dihydrophenanthrene (DHP). Despite serving for many decades as a model system for photoisomerization, the photodynamics of cis-stilbene is still not fully understood. In this work, we use ab initio multiple spawning on a SA-2-CASSCF(2,2) potential energy surface to simulate the nonadiabatic dynamics of isolated cis-stilbene. We find the cyclization (to DHP and cis-stilbene) and isomerization (to trans- and cis-stilbene) reaction coordinates to be orthogonal; branching between the two pathways is determined on the S1 excited state within 150 fs of photoexcitation. Furthermore ,trajectory basis functions (TBFs) undergoing cyclization decay rapidly to the ground state in 250 fs, while TBFs moving along the isomerization coordinate remain on the excited state longer, with the majority decaying between 300 and 500 fs. We observe three avoided crossing regions in the dynamics: two along the isomerization coordinate (displaying pyramidalization and migration of an ethylenic hydrogen or phenyl group), and one DHP-like conical intersection along the cyclization coordinate. The isomeric form of the vibrationally hot photoproducts (as determined by measurement 2 ps after photoexcitation) is determined within less than 50 fs of decay to the ground state mediated by passage through a conical intersection. Excess vibrational energy of ground state cis- and trans-stilbene is channelled into phenyl torsions (with mostly opposing directionality). Finally, our simulations are validated by direct comparison to experiment for the absorption spectrum, branching ratio of the three photoproducts (44:52:4 cis-stilbene:trans-stilbene:DHP), and excited state lifetime (520 ± 40 fs).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Circularly Polarized Attosecond Pulses Enabled by an Azimuthal Phase and Polarization Grating

High-harmonic generation (HHG) is an extreme nonlinear optical process that can map the properties of an infrared driving laser beam onto short wavelength attosecond pulse trains. However, current techniques for generating circularly polarized high harmonics for probing magnetic materials and chiral systems have limitations: two-color collinear counter-rotating driving lasers result in a low cutoff photon energy, while single-color non-collinear counter-rotating schemes suffer from low conversion efficiency. In this work, we generate circularly polarized attosecond pulse trains by using a structured laser driver which has a rotating polarization and phase grating along the azimuthal coordinate. Furthermore, our experimental and numerical results demonstrate the production of left and right circularly polarized harmonics, which naturally separate upon propagation. Our approach uses a single laser color in a collinear geometry, that can be scaled for high efficiency. Simulations show this scheme can extend into the soft x-ray region when driven by mid-infrared driving lasers, while preserving the same high phase-matching cutoff photon energy as for linearly-polarized high harmonics.

attosecond science↗

DOE OSTI2021

Robust Distributed State Estimator for Interconnected Transmission and Distribution Networks (Final Report RPPR-1)

This project’s objective is to develop a combined transmission and distribution state estimator which accounts for very large system size and model complexity (by way of distributing the computations) and large number of solar PV units connected to the distribution system on multiple feeders. The project not only provides a robust formulation and solution to this problem but also tests the solution by implementing it on a well-established large utility system. It introduces several improvements with respect to the state of the art in existing state estimation software: (a) The developed state estimator (SE) allows robust and accurate monitoring of bidirectional flows in distribution systems which result due to the distributed energy sources which are not observable and thus not incorporated in generation dispatch; (b) Large utility systems with tens of thousands of transmission buses and hundreds of thousands of distribution nodes are difficult to model as a single integrated system. This shortcoming is addressed by developing a “scalable distributed computational framework” which allows splitting the ultra large system models into several small subsystems and coordinating their solution by a robust and practical state estimation formulation; (c) Measurement errors irrespective of their locations are detected and removed by the developed state estimator. Historically, transmission and distribution systems were analyzed and operated as two independent systems. Given the non-transposed short feeder sections, unevenly loaded phases, strictly radial configuration and unidirectional power flows in the absence of remote generation, distribution system analysis was customized to account for these characteristics. However, some of these assumptions are no longer valid (non-radial configuration, bidirectional power flows) and thus distribution system analysis should be revisited. Furthermore, in the past, the interaction between the transmission and distribution systems was quite passive, where distribution substations were modeled as lumped loads in the transmission system model. With substantial generation injected by renewable generation located in the distribution systems, such modeling will no longer be accurate. The developed state estimator facilitates proper monitoring of the interactions between the transmission and distribution systems and enables smart dispatch of these units which are made observable by the state estimator.

24 POWER TRANSMISSION AND DISTRIBUTION↗

NASA NTRS1981

Simulation of the planetary boundary layer with the UCLA general circulation model

A planetary boundary layer (PBL) model is presented which employs a mixed layer entrainment formulation to describe the mass exchange between the mixed layer with the upper, laminar atmosphere. A modified coordinate system couples the mixed layer model with large scale and sub-grid scale processes of a general circulation model. The vertical coordinate is configured as a sigma coordinate with the lower boundary, the top of the PBL, and the prescribed pressure level near the tropopause expressed as coordinate surfaces. The entrainment mass flux is parameterized by assuming the dissipation rate of turbulent kinetic energy to be proportional to the positive part of the generation by convection or mechanical production. The results of a simulation of July are presented for the entire globe.

Suarez, M. J.↗

NASA NTRS1983

Determination of neutron flux distribution by using ANISN, a one-dimensional discrete S sub n ordinates transport code with anisotropic scattering

The purpose of this project was to use a one-dimensional discrete coordinates transport code called ANISN in order to determine the energy-angle-spatial distribution of neutrons in a 6-feet cube rock box which houses a D-T neutron generator at its center. The project was two-fold. The first phase of the project involved adaptation of the ANISN code written for an IBM 360/75/91 computer to the UNIVAC system at JSC. The second phase of the project was to use the code with proper geometry, source function and rock material composition in order to determine the neutron flux distribution around the rock box when a 14.1 MeV neutron generator placed at its center is activated.

Ghorai, S. K.↗

DOE OSTI2021

Multi-Sectoral Urban Interactions: Fundamental Science Needs to Inform Pathways to More Resilient Communities in a Changing Climate Multi-Sectoral Urban Interactions: Fundamental Science Needs to Inform Pathways to More Resilient Communities in a Changing Climate

Urban areas and the supply networks that support their resource use are inherently multisectoral systems composed of infrastructural, environmental, and socio-institutional components. These systems are vulnerable to accelerating and interacting stresses from climate change, population growth, resource scarcity, and land-use pressure at the same time as they have a major influence on regional and global systems. For instance, a majority of the world’s greenhouse gas emissions, food consumption, and economic activity can be attributed to urban areas. Urban areas are highly heterogenous, both across and within cities in terms of their socio-demographic, environmental, and infrastructural characteristics. This heterogeneity shapes how urban systems interact and co-evolve and contributes to different economic, environmental, and health outcomes for communities within urban areas. Urban heterogeneity also generates differential vulnerabilities to stressors and differential capacity for adapting to change. The evolution of urban space is thus critical in shaping how human societies respond to global change as they seek to improve resilience to stressors, support prosperous and equitable communities, and use natural resources in a sustainable manner. Developing fundamental scientific understanding of urban heterogeneity and system interactions across sectors and scales is critical for mapping the resilience, sustainability, and equity implications of alternative future pathways. On July 2021 the MSD Urban Systems Working Group organized this workshop to provide a venue for coordination and identification of shared objectives, research themes, and major knowledge gaps, as well as developing a shared strategy for addressing these gaps. This report summarizes the outcomes of this workshop.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2025

Crystalline Order Yet Glass-Like Heat Transport Driven by Hidden Local Distortions as the Structural Origin of Ultralow Thermal Conductivity in AgErTe 2

Given their rich chemical diversity and the interplay among the p-, d-, and f-orbitals of chalcogens, transition metals, and lanthanides, respectively, rare-earth transition-metal chalcogenides exhibit a wide variety of structural, magnetic, and transport phenomena. As a result, they form a particularly appealing platform for investigating structure–property relationships, emergent electronic and magnetic behaviors, and thermal transport. Here we investigate AgErTe 2 as a model system to understand phonon-glass behavior in ordered crystalline solids, which establishes the design principles for thermal barrier coatings and next-generation thermoelectrics. The local bonding asymmetry and lattice softness suppress the inherently low lattice thermal conductivity, resembling the characteristics of amorphous materials. This suppression is significantly influenced by local off-centering of Ag atoms, which breaks lattice periodicity while maintaining global crystallinity. The presence of antibonding states just below the Fermi level, arising from Ag 4d and Te 5p orbital interactions, leads to lattice softening and destabilizes ideal tetrahedral coordination, resulting in a pseudo Jahn–Teller distortion. Furthermore, the coexistence of weaker, more polarizable Ag–Te bonds and stronger Er–Te bonds creates a complex vibrational landscape enriched with low-frequency modes and enhanced phonon scattering. A pronounced disparity in interatomic force constants gives rise to highly localized, low-energy optical phonons linked to Ag rattling. These flat vibrational modes exhibit strong coupling with transverse acoustic phonons, resulting in ultrashort phonon lifetimes and mean free paths approaching interatomic distances. These features collectively enhance phonon scattering across a broad range of length and energy scales. This work offers a framework for engineering suppressed thermal conductivity in crystalline systems without the introduction of alloying elements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Exploring Cr and molten salt interfacial interactions for molten salt applications

Molten salts play an important role in various energy-related applications such as high-temperature heat transfer fluids and reaction media. However, the extreme molten salt environment causes the degradation of materials, raising safety and sustainability challenges. A fundamental understanding of material–molten salt interfacial evolution is needed. This work studies the transformation of metallic Cr in molten 50/50 mol% KCl–MgCl 2 via multi-modal in situ synchrotron X-ray nano-tomography, diffraction and spectroscopy combined with density functional theory (DFT) and ab initio molecular dynamics (AIMD) simulations. Notably, in addition to the dissolution of Cr in the molten salt to form porous structures, a δ-A15 Cr phase was found to gradually form as a result of the metal–salt interaction. This phase change of Cr is associated with a change in the coordination environment of Cr at the interface. DFT and AIMD simulations provide a basis for understanding the enhanced stability of δ-A15 Cr vs. bcc Cr, by revealing their competitive phase thermodynamics at elevated temperatures and probing the interfacial behavior of the molten salt at relevant facets. This study provides critical insights into the morphological and chemical evolution of metal–molten salt interfaces. Finally, the combination of multimodal synchrotron analysis and atomic simulation also offers an opportunity to explore a broader range of systems critical to energy applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Molecular Dynamics and Free Energy Calculations of Dicyclohexano-18-crown-6 Diastereoisomers with Sm 2+ , Eu 2+ , Dy 2+ , Yb 2+ , Cf 2+ , and Three Halide Salts in Tetrahydrofuran and Acetonitrile Using the AMOEBA Force Field

With the continual development of lanthanides (Ln) in current technological devices, an efficient separation process is needed that can recover greater amounts of these rare elements. Dicyclohexano-18-crown-6(DCH18C6) is a crown ether that may be a promising candidate for Ln separation, but additional research is required. As such, molecular dynamics(MD) simulations have been performed on four divalent lanthanide halide salts(Sm 2+ , Eu 2+ , Dy 2+ , and Yb 2+ ) and one divalent actinide halide salt (Cf 2+ ) bound to three diastereoisomers of DCH18C6. Dy 2+ , Yb 2+ , Cf 2+ , DCH18C6, and tetrahydrofuran (THF) solvent were parameterized for the AMOEBA polarizable force field for the first time, whereas existing parameters for Sm 2+ and Eu 2+ were utilized from our previous efforts. A coordination number (CN) of six for Ln 2+ /An 2+ –O solvated in THF indicated that the cations interacted almost entirely with the oxygens of the polyether ring. A CN of one for Ln 2+ /An 2+ -N solvated in acetonitrile for systems containing iodide suggested that theN atom of acetonitrile was competitive with I – for cation interactions. Fluctuation between five and six CNs for Dy 2+ and Yb 2+ suggested that although the cations remained in the polyether ring, the size of the ring may not be an ideal fit as these cations possess comparatively smaller ionic radii. Gibbs binding free energies of Sm 2+ in all DCH18C6 diastereoisomers solvated in THF were calculated. The binding free energy of the cis-syn-cis diastereoisomer was the most favorable, followed by cis-anti-cis, and then trans-anti-trans. Lastly, two major types of conformation were observed for each diastereoisomer that were related to the electrostatic interactions and charge density of the cations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1995

Korea's developmental program for superconductivity

Superconductivity research in Korea was firstly carried out in the late 70's by a research group in Seoul National University (SNU), who fabricated a small scale superconducting magnetic energy storage system under the financial support from Korea Electric Power Company (KEPCO). But a few researchers were involved in superconductivity research until the oxide high Tc superconductor was discovered by Bednorz and Mueller. After the discovery of YBaCuO superconductor operating above the boiling point of liquid nitrogen (77 K)(exp 2), Korean Ministry of Science and Technology (MOST) sponsored a special fund for the high Tc superconductivity research to universities and national research institutes by recognizing its importance. Scientists engaged in this project organized 'High Temperature Superconductivity Research Association (HITSRA)' for effective conducting of research. Its major functions are to coordinate research activities on high Tc superconductivity and organize the workshop for active exchange of information. During last seven years the major superconductivity research has been carried out through the coordination of HITSRA. The major parts of the Korea's superconductivity research program were related to high temperature superconductor and only a few groups were carrying out research on conventional superconductor technology, and Korea Atomic Energy Research Institute (KAERI) and Korea Electrotechnology Research Institute (KERI) have led this research. In this talk, the current status and future plans of superconductivity research in Korea will be reviewed based on the results presented in interim meeting of HITSRA, April 1-2, 1994. Taejeon, as well as the research activity of KAERI.

Hong, Gye-Won↗

DOE OSTI2025

Coordinated Thermal Safety Attack and Defense on EV Battery Management Systems

Battery temperature sensor and battery current sensor data which are key sensing inputs to the Battery Management Controllers in electric vehicles, are vulnerable to possible cyber/ physical manipulation due to known vulnerabilities inherited from CAN bus technology that is used for in-vehicle communications between electronic control units that transfer sensing and control data. In this paper, we first create a simulation that enables us to evaluate impact of cyber physical attacks on electric vehicle battery management system in a controlled environment that violates thermal safety. Specifically, we emulate a Level 3 - DC fast charging system with SAE J1772/CCS, integrated with standard charging controls and thermal safety controls on EVs, and various sensing data flows. Second, we propose a coordinated current and battery temperature attack that has crippling economic, and safety impacts. Third, we quantify the usability, economic and safety impacts of such attacks as a function of the extent of data manipulation. Finally, we propose a physics model driven detection technique to detect presence of such attacks.

25 ENERGY STORAGE↗

DOE OSTI2024

Simulating Crystallization in a Colloidal System Using State Predictive Information Bottleneck Based Enhanced Sampling

Here, we investigate crystal nucleation in supersaturated colloid suspensions using enhanced molecular dynamics simulations augmented with machine learning techniques. The simulations reveal that crystallization in the model colloidal system studied here, with particles interacting through a repulsive screened Coulomb Yukawa potential, proceeds from vapor to dense liquid droplet to crystalline phases across multiple high barriers. Employing a one-dimensional reaction coordinate derived from the State Predictive Information Bottleneck framework, our simulations capture back-and-forth phase transitions across multiple barriers effectively in biased metadynamics simulations. We obtain relative free energy differences between different phases and also quantify the roles of different molecular level features in driving the phase changes.

Chemistry↗