Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atomic ordering”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Designing Physicochemically‐Ordered Interphases for High‐Performance Composites

To enhance the mechanical properties of carbon fiber‐reinforced polymer composites, a physicochemical scaffold is designed incorporating microscopically architected chemically reactive nanofibers that act as a multiscale bridge between the carbon fibers and the matrix. Thermally activated nanofibers leverage their morphologically driven mechanochemical properties to form covalent bonds with adjacent polymer molecules, creating a co‐continuous network that dramatically enhances fiber‐matrix load transfer. By meticulously controlling the nanofiber architecture through variable surface area, functional group availability, and polymer chain alignment effects, the extent of covalent bonding between nanofibers and the matrix is manipulated ultimately resulting in improved carbon fiber‐matrix adhesion. Further, the concept was validated using polyacrylonitrile nanofibers within an acrylonitrile butadiene styrene matrix in a discontinuous carbon fiber‐reinforced composite system. Nanomechanical studies using atomic force microscopy and low‐field nuclear magnetic resonance spectroscopy confirmed immobilized, chemically transferred, and ordered nanostructures at the interphase. The resulting composites demonstrated ≈56% and ≈175% improvements in tensile strength and toughness, respectively, compared to composites without nanofiber. Comprehensive thermal, rheological, and X‐ray scattering analyses, along side all‐atomic molecular dynamics simulations, revealed the fundamental mechanisms behind these improvements in mechanical behavior. The versatility and efficacy of the approach have the potential to address longstanding interphase challenges in the composite industry.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Catalytic Site Requirements for N2O Decomposition on Cu-, Co-, and Fe-SSZ-13 Zeolites

N2O decomposition is investigated on Cu, Co and Fe-exchanged SSZ-13 zeolite catalysts at relatively low metal loadings. The catalysts are synthesized by solution ion exchange, and subjected to X-ray diffraction (XRD), temperature-programed-reduction by H2 (H2-TPR), temperature-programed-reaction of N2O (N2O-TPR) coupled with in-situ transmission FTIR, and finally steady-state flow reaction tests. At low N2O pressures (< 0.05 kPa), all catalysts display pseudo first-order kinetics. From Arrhenius analysis, Cu and Fe-SSZ-13 display very different apparent activation energies but similar pre-exponential factors, suggesting their similar reaction mechanisms. N2O decomposition follows a dual-site mechanism, occurring on dimeric M-O-M sites in these catalysts, and O2 is formed by the combination of two O ad-atoms from two vicinal metal sites. Under low N2O pressure (0.05 kPa) and first-order kinetic regime, the reaction is limited by N-O cleavage on bare metal active sites. In comparison to Cu-SSZ-13, the much higher N2O decomposition rate over Fe-SSZ-13 is attributed to the much lower activation barriers for the N-O cleavage step. N2O decomposition occurs on isolated Co2+ ions in Co-SSZ-13. The rate-limiting step is N-O cleavage on an O-occupied Co site in the low-pressure first order kinetic regime. This single-site mechanism leads to much higher pre-exponential factors as compared to the dual-site mechanism. This beneficial factor for reaction rate enhancement, however, is compromised by the much higher activation barriers over this catalyst.

Lin, Fan↗
DOE OSTI2022

General solution to the Kohn–Luttinger nonconvergence problem

A simple, but general solution is proposed for the Kohn–Luttinger problem, i.e., the nonconvergence of the finite-temperature many-body perturbation theory with its zero-temperature counterpart as temperature is lowered to zero under some circumstances. How this nonconvergence can be avoided by altering the reference wave function is illustrated numerically by using up to the fifth order of the perturbation theory.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Epitaxial columnar growth of strain-free antiferromagnetic Weyl semimetal Mn 3 Sn on wurtzite c -plane GaN/Al 2 O 3 (0001)

Weyl semimetal thin films with excellent crystalline quality are of great interest for antiferromagnetic spintronics. Mn 3 Sn is one Weyl semimetal with great properties and promise for exciting science and applications. It has proven very challenging, however, to grow Mn 3 Sn thin films with smooth surfaces, negligible strain, and excellent crystallinity. In this work, we discuss the successful preparation of epitaxial Mn 3 Sn (0001)-oriented thin films via molecular beam epitaxial growth on c -plane wurtzite GaN which was grown by MBE on Al 2 O 3 (0001). We present the reflection high energy electron diffraction analysis along with x-ray diffraction in order to demonstrate the crystalline quality of the film, and we give atomic models to explain the epitaxial orientation relationships between the crystal lattices of the substrate, GaN layer, and Mn 3 Sn layer. Importantly, we discuss the film lattice parameters as compared to expected values, demonstrating negligible strain both in-plane and out-of-plane . Atomic force microscopy reveals an epitaxial columnar growth mode characterized by flat-top-mesa islands, while scanning tunneling microscopy shows the atomically smooth surfaces of the mesa-top structures. Finally, Rutherford backscattering informs the stoichiometry of the film as well as the layer thicknesses.

Chiral antiferromagnetic material↗
DOE OSTI2020

Effects of irradiation spectrum on the microstructural and mechanical properties of bulk metallic glasses

Zr 52.5 Cu 17.9 Ni 14.6 Al 10 Ti 5 (BAM-11) and Cu60Zr 20 Hf 10 Ti 10 bulk metallic glass (BMG) specimens were ion irradiated to characterize and compare their irradiation induced microstructural and mechanical property evolution. For the ion irradiations, samples of each BMG were exposed to 9 MeV Ni 3+ and 5.5 MeV C + ions to a midrange (~1.5 μm depth) dose of 0.5 displacements per atom (dpa) at temperatures ranging from room temperature to 360 °C in order to examine potential effects associated with low (5.5 MeV C) and high (9 MeV Ni) average primary knock on energies. BAM-11 BMG samples were also irradiated by neutrons (E > 0.1 MeV) at ~70 °C to a fluence of 1.4 × 10 20 n/cm 2 (dose of 0.1 dpa). Importantly, no pronounced irradiation spectrum effects were observed for the nanoindentation mechanical properties of the two BMGs following ion irradiation. However, significant softening was observed in the BAM-11 BMG sample irradiated by neutrons, while annealing of the unirradiated samples led to a marked increase in hardening. Finally, results of the nanoindentation experiments indicate that softening is caused by irradiation induced creation of free volume defects, while the hardening is caused by their annihilation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Identification and correction of temporal and spatial distortions in scanning transmission electron microscopy

Scanning transmission electron microscopy (STEM) has become the technique of choice for quantitative characterization of atomic structure of materials, where the minute displacements of atomic columns from high-symmetry positions can be used to map strain, polarization, octahedra tilts, and other physical and chemical order parameter fields. The latter can be used as inputs into mesoscopic and atomistic models, providing insight into the correlative relationships and generative physics of materials on the atomic level. However, these quantitative applications of STEM necessitate understanding the microscope induced image distortions and developing the pathways to compensate them both as part of a rapid calibration procedure for in situ imaging, and the post-experimental data analysis stage. Here, we explore the spatiotemporal structure of the microscopic distortions in STEM using multivariate analysis of the atomic trajectories in the image stacks. Based on the behavior of principal component analysis (PCA), we develop the Gaussian process (GP)-based regression method for quantification of the distortion function. The limitations of such an approach and possible strategies for implementation as a part of in-line data acquisition in STEM are discussed. Here, the analysis workflow is summarized in a Jupyter notebook that can be used to retrace the analysis and analyze the reader's data.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Role of Pairwise Reactions on the Synthesis of Li 0.3 La 0.57 TiO 3 and the Resulting Structure–Property Correlations

The performance of single-ion conductors is highly sensitive to the material’s defect chemistry. Tuning these defects is limited for solid-state reactions as they occur at particle–particle interfaces, which provide a complex evolving energy landscape for atomic rearrangement and product formation. In this report, we investigate the (1) order of addition and (2) lithium precursor decomposition temperature and their effect on the synthesis and grain boundary conductivity of the perovskite lithium lanthanum titanium oxide (LLTO). We use an intimately mixed sol–gel, a solid-state reaction of Li precursor + La 2 O 3 + TiO 2 , and Li precursor + amorphous La 0.57 TiO x as different chemical routes to change the way in which the elements are brought together. The results show that the perovskite can accommodate a wide range of Li deficiencies (upward of 50%) while maintaining the tetragonal LLTO structure, indicating that X-ray diffraction (XRD) is insufficient to fully characterize the chemical nature of the product (i.e., Li-deficient LLTO may behave differently than stoichiometric LLTO). Variations in the relative intensities of different reflections in XRD suggest variations in the La ordering within the crystal structure between synthesis methods. Furthermore, the choice of the precursor and the order of addition of the reactants lower the time required to form a pure phase. Density functional theory calculations of the formation energy of possible reaction intermediates support the hypothesis that a greater thermodynamic driving force to form LLTO leads to a greater LLTO yield. Here, the retention of lithium is correlated with the thermal decomposition temperature of the Li precursor and the starting material mixing strategy. Taking the results together suggests that cations that share a site with Li should be mixed early to avoid ordering. Such cation ordering inhibits Li motion, leading to higher Li ion resistance.

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

Atomic Origin of the Autocatalytic Reduction of Monoclinic CuO in a Hydrogen Atmosphere

Reducibility is key for the use of bulk metal oxides in chemical transformation involving redox reactions but probing microscopic processes of oxide reduction is challenging. This is because the insulating nature of bulk oxides restricts ion and electron spectroscopic measurements of oxide surfaces. Herein, using a combination of environmental transmission electron microscopy and atomistic modeling, we report direct in-situ atomic-scale observations of the surface and subsurface dynamics and show that the hydrogen induced CuO reduction occurs via the surface decay of Cu-O/Cu bilayer atomic steps, formation of partially reduced CuO superstructure by the self-ordering of O vacancies in the subsurface, and collapse of Cu-O layers in the bulk. In this work, all these substeps can be traced back to the progressively increased concentration and activity of O vacancies in the surface and subsurface of the oxide, thereby leading to the self-accelerated oxide reduction. These results demonstrate the microscopic details that may have broader applicability in modulating various redox processes.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Evidence for a single-layer van der Waals multiferroic

Here we searched for signatures of type-II multiferroic order in ultrathin van der Waals materials, to determine if such a state, comprising coexisting electrical polarization and spin order with strong magnetoelectric effects, could persist down to a single-atomic layer. We reported the first observation of such a state in triangular lattice material NiI 2 . This system has a layered structure and built-in magnetic frustration that promotes a complex magnetic ground state with proper-screw spin helices. The spin helix pattern breaks inversion symmetry and induces electrical polarization via magnetoelectric coupling. We use circular dichroic Raman measurements to directly probe the magneto-chiral ground state and its electromagnon modes originating from dynamic magnetoelectric coupling. Using birefringence and second-harmonic generation measurements, we detect a highly anisotropic electronic state simultaneously breaking three-fold rotational and inversion symmetry. The evolution of the optical signatures as a function of temperature and layer number surprisingly revealed a multiferroic state that persists down to the ultrathin limit of monolayer NiI 2 . These observations establish NiI 2 and transition metal dihalides as a new platform for studying emergent multiferroic phenomena, chiral magnetic textures and ferroelectricity in the 2D limit.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Optimization of the Reax force field for the lithium–oxygen system using a high fidelity charge model

Studies using molecular dynamics (MD) have long struggled to simulate the failure modes of materials, predicting unrealistically high ductility and failing to capture brittle fracture. The primary cause of this shortcoming is an inadequate description of bond breaking. While reactive force fields such as ReaxFF show improvements compared to traditional force fields, the charge models used yield unphysical partial charges, especially during dissociation of ionic bonds. This flaw may be remedied by using the atom-condensed Kohn–Sham density functional theory (DFT) approximated to a second order (ACKS2) charge model for determining partial charges. In this work, we present a new ACKS2-enabled Reax force field for fracture simulations of lithium oxide systems, which was obtained by training against an extensive set of DFT, multireference configuration interaction (MRCI), and MRCI+Q reference data using genetic optimization techniques. This new force field significantly improves the bond breaking behavior, but still cannot fully capture the brittle fracture in MD simulations, suggesting more research is needed to improve simulation of brittle fracture.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Kinetics of short-range order formation in GeSn alloy: MBE vs CVD

Recently, short-range order (SRO) has attracted significant attention, challenging the conventional view of the atomic positions in alloys being random. Furthermore, the presence of SRO has been predicted to have profound effects on the electronic and topological properties of group-IV alloys, offering a different direction in designing group-IV materials for photoelectronic and quantum devices. However, due to the limited understanding of the formation mechanisms, developing effective methods to manipulate SRO in epitaxy is still challenging. To address this, we propose a mechanism for the GeSn alloy, revealing that surface diffusion plays a key role in SRO formation. Building on this mechanism, we show that the distinct surface conditions in MBE and CVD lead to the formation of SRO with enhanced Sn–Sn pairing in MBE-grown samples, while CVD-grown samples remain random alloys. Furthermore, our findings provide an initial understanding of the kinetic process of SRO formation, providing guidance for the design of experiments to manipulate SRO.

Alloys↗
DOE OSTI2023

World beyond the nearest neighbors

The structure beyond the nearest neighbor atoms in liquid and glass is characterized by the medium-range order (MRO). In the conventional approach, the MRO is considered to result directly from the short-range order (SRO) in the nearest neighbors. To this bottom–up approach starting with the SRO, we propose to add a top–down approach in which global collective forces drive liquid to form density waves. The two approaches are in conflict with each other, and the compromise produces the structure with the MRO. The driving force to produce density waves provides the stability and stiffness to the MRO, and controls various mechanical properties. This dual framework provides a novel perspective for description of the structure and dynamics of liquid and glass.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Precision measurement of lifetime and branching ratios of the 4⁢𝑓 13⁡ 5⁢𝑑⁢6⁢𝑠 1 ⁢[5/2] 5/2 state in Yb + ions

We report spectroscopic and time-resolved experimental observations to characterize the [Xe]⁢4⁢𝑓 13 ⁡( 2 𝐹$^𝑜_{5/2}$)⁢5⁢𝑑⁢6⁢𝑠⁡( 1 𝐷)⁢ 1 ⁢[5/2]$^𝑜_{5/2}$ state in 172 Yb + ions. We access this state from the metastable 4⁢𝑓 14 ⁡5⁢𝑑⁡( 2 𝐷 3/2,5/2 ) manifold and observe an unexpectedly long lifetime of 𝜏 = 37.9⁢(9) ⁢µ⁢s that allows visible Rabi oscillations and resolved-sideband spectroscopy. Using a combination of coherent population dynamics, high-fidelity detection and heralded state preparation, and optical pumping methods, we measure the branching ratios to the 2 𝐷 3/2 , 2 𝐷 5/2 , and 2 𝑆 1/2 states to be 0.359⁢(2), 0.639(2), and 0.0023⁢(16), respectively. The branching ratio to the 4⁢𝑓 13⁡ 6⁢𝑠 2 ⁡( 2 𝐹 7/2 ) is compatible with zero within our experimental resolution. We also report measurements of Landé 𝑔-factor of the 1 ⁢[5/2]$^𝑜_{5/2}$ state. Further, the branching ratio of the 2 𝐷 5/2 to 2 𝑆 1/2 decay in 172 Yb + was measured to be 0.188(3), improving its relative uncertainty by an order of magnitude. Our measurements provide experimental benchmarks for better understanding the atomic structure of Yb + ions, which still lacks accurate numerical descriptions, and the use of high-lying excited states for partial detection and qubit manipulation in the omg architecture.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

matsim-agents v1.0

matsim-agents is a multi-agent AI framework for atomistic materials simulation and discovery. It orchestrates large language models (LLMs), machine-learned interatomic potentials (MLIPs), and DFT codes into a single agentic loop running on laptops and DOE leadership-class supercomputers. MULTI-AGENT ORCHESTRATION A LangGraph state machine with three nodes: a Planner that converts a natural-language research objective into structured tasks; an Executor that dispatches atomistic tools and loops until the queue is empty; and an Analyst that summarizes results into a human-readable report. State is checkpointed after every step and human-in-the-loop gates can be inserted at any edge. HYPOTHESIS-DRIVEN DISCOVERY CHAT An interactive REPL (matsim-agents chat) that couples LLM dialogue with atomistic simulation. Chemical formulas are automatically detected in conversation turns and trigger a full crystal-phase exploration: structure generation → relaxation → stability scoring → result injection back into the conversation, creating a closed hypothesis-refinement loop. CRYSTAL PHASE ENUMERATION Given a composition, the phase explorer enumerates prototypes by stoichiometry: elemental (fcc/bcc/hcp/sc/diamond), binary 1:1 (rocksalt/CsCl/zincblende/ wurtzite/fluorite/rutile), ternary 1:1:3 (cubic perovskite), ternary 1:2:4 (perovskite + spinel), quaternary 1:1:2:6 (Fm-3m double perovskite). 2-D prototypes (graphene, h-BN, MoS2 2H/1T) and multilayer stacking are also supported via --include-2d and --num-layers. SUPERCELL GENERATION AND SITE DECORATION Auto-tiling to a minimum atom count (--min-atoms), explicit NxNxN tiling (--supercell), symmetry-distinct site decorations (--n-orderings), and isotropic lattice-scale sweeps (--lattice-scales) for volume bracketing. MLFF RELAXATION AND STABILITY SCORING HydraGNN (multi-headed GNN) drives structure relaxation via ASE with FIRE, BFGS, or BFGSLineSearch. Stability output: delta-E/atom ranking across phases and a max-residual-force dynamical-stability proxy. Other MLIPs (MACE, NequIP, Orb) can be plugged in through the same interface. DFT BACKENDS Quantum ESPRESSO pw.x and VASP 6.6 are first-class labellers. Both have validated GPU builds and SLURM/PBS launchers for three DOE platforms: Frontier (AMD MI250X, ROCm), Aurora (Intel PVC, oneAPI), Perlmutter (NVIDIA A100, CUDA). QE produces ~100 binaries (pw.x, ph.x, epw.x, ...). VASP supports scf, relax, vc-relax, and vc-relax-shape run types. ACTIVE-LEARNING LOOP matsim-agents al run CONFIG.yaml drives an iterative HydraGNN-DFT loop: MD generates candidates → ensemble/MC-dropout uncertainty selects the most informative → DFT labels them in parallel inside one allocation → dataset grows → HydraGNN retrains → repeat. DFT backend is a single YAML toggle (dft.backend: vasp | qe). LLM-generated seed structures are supported (no curated POSCAR library needed). Config uses ${VAR}, ${VAR:-default}, ${VAR:?msg} shell-style substitution for cross-user/cross-site portability. LLM BACKENDS Ollama (local, default), vLLM (HPC multi-GPU serving), OpenAI, Anthropic, HuggingFace Transformers+Accelerate. Selected at runtime via flag or env var with no code changes. HPC PORTABILITY Same Python entry points run on Frontier (ROCm 7.2), Aurora (oneAPI), and Perlmutter (CUDA 12). DFT and ML stacks are never co-loaded in the same shell; they couple through the scheduler and filesystem. Advanced multi-node launchers (serve, discovery-chat, single-relaxation, active-learning, QE warm-start) are provided for all three platforms. CODABENCH COMPETITION BUNDLE A self-contained benchmark: 159 atomistic test structures across 11 material classes, 5 tasks (formation energy, forces, ML relaxation, AI-DFT relaxation, phase stability ranking), public/private leaderboard split (30/70), and four ready-to-run baselines: MACE-MP-0, HydraGNN, UMA, AllScAIP.

Lupo Pasini, Massimiliano [Oak Ridge National Labo↗
DOE OSTI2020

Modeling Catalyzed Growth of Single-Wall Carbon Nanotubes (Final Report)

Single-walled carbon nanotubes (SWCNTs) are among the most promising nanostructures because of their special and unusual physico-chemical properties that make them strong candidates for a variety of applications. Since their discovery, chemical vapor deposition (CVD) methods have become the preferred synthesis method. CVD uses a carbon-containing precursor species, which is decomposed over a catalyst yielding C atoms that eventually evolve in carbon nanostructures. CVD temperatures are in the order of 1000K or lower, and pressures in the range of 1-5 atm, although some processes use higher pressures. Depending on the nature of the catalyst, and on the conditions of pressure, temperature, and presence of other inert compounds, the resultant nanostructures may include nanofibers, nano-onions, multiple-, or single-wall carbon nanotubes, among others. When the CVD synthesis process is appropriately tuned, however, it is possible to obtain only SWCNTs, which are geometrically characterized by their diameter and by their chiral angle. The chiral indexes (n,m) determine the degree of helical twist of the graphite lattice along the nanotube axis, where the term chirality is derived from the left and right-handed helicity of the tubes defined by their chiral indexes (n,m) and (m,n) respectively. The main goal of this research is to build fundamental understanding that will help elucidating the nature and structure of the actual catalyst in a combined catalyst/support system, which we claim is essential especially when the catalyst is of nanodimensions. Given that we are considerably advanced in the understanding of selective single-walled carbon nanotube catalyst synthesis, we take this process as a case study. Thus, we focus on the understanding of the selective growth of SWCNTs on supported nanocatalysts in systems where the catalyst synthesis takes place in the same reactor preceding the nanotube synthesis. Once we understand what type of catalyst structure would be necessary for this specific growth, we will try to elucidate the catalytic synthesis process that leads to such catalytic structure. Thus, our aim of controlling nanostructure synthesis may start to be delineated.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Towards computational polar-topotronics: Multiscale neural-network quantum molecular dynamics simulations of polar vortex states in SrTiO3/PbTiO3 nanowires

Recent discoveries of polar topological structures ( e.g ., skyrmions and merons) in ferroelectric/paraelectric heterostructures have opened a new field of polar topotronics. However, how complex interplay of photoexcitation, electric field and mechanical strain controls these topological structures remains elusive. To address this challenge, we have developed a computational approach at the nexus of machine learning and first-principles simulations. Our multiscale neural-network quantum molecular dynamics molecular mechanics approach achieves orders-of-magnitude faster computation, while maintaining quantum-mechanical accuracy for atoms within the region of interest. This approach has enabled us to investigate the dynamics of vortex states formed in PbTiO 3 nanowires embedded in SrTiO 3 . We find topological switching of these vortex states to topologically trivial, uniformly polarized states using electric field and trivial domain-wall states using shear strain. These results, along with our earlier results on optical control of polar topology, suggest an exciting new avenue toward opto-electro-mechanical control of ultrafast, ultralow-power polar topotronic devices.

Linker, Thomas↗
DOE OSTI2022

Nitrogen-sulfur-carbon nanocomposites and their application as cathode materials in lithium-sulfur batteries

The invention is directed in a first aspect to electron-conducting porous compositions comprising an organic polymer matrix doped with nitrogen atoms and having elemental sulfur dispersed therein, particularly such compositions having an ordered framework structure. The invention is also directed to composites of such S/N-doped electron-conducting porous aromatic framework (PAF) compositions, or composites of an S/N-doped mesoporous carbon composition, which includes the S/N-doped composition in admixture with a binder, and optionally, conductive carbon. The invention is further directed to cathodes for a lithium-sulfur battery in which such composites are incorporated.

25 ENERGY STORAGE↗
DOE OSTI2023

Reinforcement learning based hybrid bond-order coarse-grained interatomic potentials for exploring mesoscale aggregation in liquid–liquid mixtures

Exploring mesoscopic physical phenomena has always been a challenge for brute-force all-atom molecular dynamics simulations. Although recent advances in computing hardware have improved the accessible length scales, reaching mesoscopic timescales is still a significant bottleneck. Coarse-graining of all-atom models allows robust investigation of mesoscale physics with a reduced spatial and temporal resolution but preserves desired structural features of molecules, unlike continuum-based methods. Here, we present a hybrid bond-order coarse-grained forcefield (HyCG) for modeling mesoscale aggregation phenomena in liquid–liquid mixtures. The intuitive hybrid functional form of the potential offers interpretability to our model, unlike many machine learning based interatomic potentials. We parameterize the potential with the continuous action Monte Carlo Tree Search (cMCTS) algorithm, a reinforcement learning (RL) based global optimizing scheme, using training data from all-atom simulations. The resulting RL-HyCG correctly describes mesoscale critical fluctuations in binary liquid–liquid extraction systems. cMCTS, the RL algorithm, accurately captures the mean behavior of various geometrical properties of the molecule of interest, which were excluded from the training set. The developed potential model along with the RL-based training workflow could be applied to explore a variety of other mesoscale physical phenomena that are typically inaccessible to all-atom molecular dynamics simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗