SEARCH · Engineering Papers

Results for “SERIES EXPANSION”

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 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Utilizing digitized occurrence records of Midwestern feral Cannabis sativa to develop ecological niche models

Hemp (Cannabis sativa L.) has historically played a vital role in agriculture across the globe. Feral and wild populations have served as genetic resources for breeding, conservation, and adaptation to changing environmental conditions. However, feral populations of Cannabis, specifically in the Midwestern United States, remain poorly understood. This study aims to characterize the abiotic tolerances of these populations, estimate suitable areas, identify regions at risk of abiotic suitability change, and highlight the utility of ecological niche models (ENMs) in germplasm conservation. The Maxent algorithm was used to construct a series of ENMs. Validation metrics and MOP (Mobility-oriented Parity) analysis were used to assess extrapolation risk and model performance. We also projected the final projected under current and future climate scenarios (2021–2040 and 2061–2080) to assess how abiotic suitability changes with time. Climate change scenarios indicated an expansion of suitable habitat, with priority areas for germplasm collection in Indiana, Illinois, Kansas, Missouri, and Nebraska. This study demonstrates the application of ENMs for characterizing feral Cannabis populations and highlights their value in germplasm conservation and breeding efforts. Populations of feral C. sativa in the Midwest are of high interest, and future research should focus on utilizing tools to aid the collection of materials for the characterization of genetic diversity and adaptation to a changing climate.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

The Many-Body Expansion for Metals: II. Nonadditive Terms in Clusters Composed of Metals with n s 1 , n s 2 , and n s 2 p 1 Configurations

The many-body expansion (MBE) was applied to homometallic and heterometallic trimers of metals with ns 1 , ns 2 , and ns 2 p 1 configurations to investigate its convergence, the magnitude and nature (stabilizing/destabilizing) of the individual terms and seek an understanding of their variation across the different families of clusters. In particular, we examined the series of alkali metals (Li 3 , Na 3 , K 3 , Li 2 Na, LiNa 2 ), alkali metal borides (Li 2 B and LiB 2 ), and alkaline earth metals (Be 3 , Be 2 Mg, BeMg 3 , and Mg 3 ) trimers, as well as sodium clusters Na n , n = 2–5. Here, we found that there is no uniform contribution (stabilizing or destabilizing) across the series in the different families of trimers. For instance, the 2-B term stabilizes the ground states of the Na 3 (doublet), Na 4 (singlet), and Na 5 (doublet) clusters, and the 3-B term destabilizes them; however, the opposite holds for the quartet state of the Li 3 , Li 2 Na, LiNa 2 , and Na 3 clusters (destabilizing 2-B, stabilizing 3-B). Substituting Li with B in the quartet state of Li 3 results in a significant reduction of the 3-B term amounting to 16% (Li 2 B) and 5% (LiB 3 ) of the binding energy. On the contrary, the ground states of the alkaline earth metal clusters (Be 3 , Be 2 Mg, BeMg 3 , and Mg 3 ) are stabilized by the 3-B term, while the 2-B term destabilizes them. Overall, we find that the 3-B terms significantly stabilize the high-spin multiplicity states of the ns 1 configurations and the low-spin states of the ns 2 configurations. Finally, as the size of the metal increases, the contribution of the 3-B term to the binding energy decreases due to the longer metal–metal bond distances.

Alkali metals↗

DOE OSTI2024

Strategies to Obtain Reliable Energy Landscapes from Embedded Multireference Correlated Wavefunction Methods for Surface Reactions

Embedded correlated wavefunction (ECW) theory is a powerful tool for studying ground- and excited-state reaction mechanisms and associated energetics in heterogeneous catalysis. Several factors are important to obtaining reliable ECW energies, critically the construction of consistent active spaces (ASs) along reaction pathways when using a multireference correlated wavefunction (CW) method that relies on a subset of orbital spaces in the configuration interaction expansion to account for static electron correlation, e.g., complete AS self-consistent field theory, in addition to the adequate partitioning of the system into a cluster and environment, as well as the choice of a suitable basis set and number of states included in excited-state simulations. Here, in this work, we conducted a series of systematic studies to develop best-practice guidelines for ground- and excited-state ECW theory simulations, utilizing the decomposition of NH 3 on Pd(111) as an example. We determine that ECW theory results are relatively insensitive to cluster size, the aug-cc-pVDZ basis set provides an adequate compromise between computational complexity and accuracy, and that a fixed-clean-surface approximation holds well for the derivation of the embedding potential. Additionally, we demonstrate that a merging approach, which involves generating ASs from the molecular fragments at each configuration, is preferable to a creeping approach, which utilizes ASs from adjacent structures as an initial guess, for the generation of consistent potential energy curves involving open-d-shell metal surfaces, and, finally, we show that it is essential to include bands of excited states in their entirety when simulating excited-state reaction pathways.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

MURR LEU structural and thermal hydraulics analyses: Part II – Impacts of irradiation thermo-mechanical behavior on thermal hydraulics safety analyses

A series of structural analyses have been performed to support the conversion of the University of Missouri Research Reactor (MURR) from the use of highly enriched uranium (HEU; ≥20 wt% U-235) to low-enriched uranium (LEU; <20 wt% U-235) fuel. The irradiation thermo-mechanical analysis evaluated the effects of fuel swelling, irradiation creep, thermal expansion, as well as thermal resistance from the oxide layer growth for the MURR LEU element in prototypic thermal and irradiation conditions as presented in Part I of this article. Overall, this irradiation thermo-mechanical analysis predicts smaller gap thickness reductions in previously limiting regions, and larger reductions in the middle of the outermost end channels where power density is not typically a maximum. Due to substantial differences between the channel gap reductions assumed for the previous safety analyses and those predicted by the irradiation thermo-mechanical analysis, a need to evaluate their impact on the thermal hydraulics safety analyses arose. This article presents the results from the steady-state safety analysis for normal operation as well as the two most limiting accident scenarios. The calculation models were revised in order to account for the spatial and temporal variation of the channel gap thicknesses. The results show that sufficient safety margins are still maintained for normal operation as well as during the postulated accident transients. This work provides a methodology of incorporating the irradiation thermo-mechanical behavior of plate-type fuel into the thermal hydraulics safety analyses.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2026

Viability of perturbative expansion for quantum field theories on neurons

Neural Network (NN) architectures that break statistical independence of parameters have been proposed as a new approach for simulating local quantum field theories (QFTs) [1]. In the infinite neuron number limit, single-layer NNs can exactly reproduce QFT results. This paper examines the viability of this architecture for perturbative calculations of local QFTs for finite neuron number N using scalar ϕ 4 theory in d Euclidean dimensions as an example. We find that the renormalized O(1/N ) corrections to two-and four-point correlators yield perturbative series which are sensitive to the UV cut-off and therefore have a weak convergence. We propose a modification to the architecture to improve this convergence and discuss constraints on the parameters of the theory and the scaling of N which allow us to extract accurate field theory results.

Sen, Srimoyee [Iowa State University, Ames, IA (Un↗

DOE OSTI2024

Powered By: Transmission Planning

Presented as part of the NREL Grid Planning and Analysis Center's (GPAC) "Powered By" Series, this iteration focuses not on a specific tool but instead on a domain of expertise/capability being built out at NREL: transmission planning. Various sets of tools are integrated together in this presentation to present a combination of transmission focus areas, case studies, tools and future directions of transmission research.

capacity expansion↗

DOE OSTI2026

The high explosives & affected targets (HEAT) dataset

Artificial Intelligence (AI) surrogate models offer a computationally efficient alternative to full-physics simulations, yet no existing datasets are publicly available for training, testing, and validation of machine learning models of the dynamics of high-explosive driven shocks through multiple materials. Shock propagation through materials is a computationally challenging problem because simulations must include material-specific equations of state (EOS) along with descriptions of other physical processes such as plastic deformation, phase change, damage processes, fluid instabilities, and multi-material interactions. Shocks are typically initiated by high-velocity impacts or explosive loading. The latter case necessitates the addition of models of reactive materials to represent high-explosive (HE) detonation. Here, to address the lack of an expansive dataset for multi-material shock propagation in the AI/ML community, we present the High-Explosives and Affected Targets (HEAT) Dataset. HEAT is a physics-rich collection of two-dimensional, cylindrically symmetric, simulations generated using an Eulerian, multi-material, shock-propagation code developed at Los Alamos National Laboratory. The dataset includes two partitions: (1) the expanding shock-cylinder (CYL) simulations, Figs. 1, and (2) the Perturbed Layered Interface (PLI) simulations, Fig. 2. Entries in both partitions consist of time series of arrays of thermodynamic fields (pressure, density, and temperature), kinematic fields (position and velocity), and additional fields that depend on thermodynamic and/or kinematic fields (e.g., material stress). Materials in the CYL partition include solids (aluminium, copper, depleted uranium, stainless steel, tantalum, and a generic polymer), a liquid (water), gases (air, nitrogen), and a generic detonating material (high explosive, HE). The PLI partition spans a highly varying geometry but consists of fixed materials across entries: Copper, aluminium, stainless steel, generic polymer, and generic HE. HEAT captures critical phenomena such as momentum transfer, shock propagation, plastic deformation, and thermal effects, making HEAT a valuable benchmark for development of AI/ML emulation of multi-material shock propagation.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Locating the QCD critical point through contours of constant entropy density

We propose a new method to investigate the existence and location of the conjectured high-temperature critical point of strongly interacting matter via contours of constant entropy density. By approximating these lines as a power series in the baryon chemical potential 𝜇 𝐵 , one can extrapolate them from first-principle results at zero net-baryon density, and use them to locate the quantum chromodynamics (QCD) critical point, including the associated first-order and spinodal lines. As a proof of principle, we employ currently available continuum-extrapolated lattice data from the Wuppertal-Budapest collaboration to find a critical point at a temperature and a baryon chemical potential of 𝑇 𝑐 = 114.3 ± 6.9 MeV and 𝜇 𝐵,𝑐 = 602.1 ± 62.1 MeV, respectively, at expansion order 𝒪⁡(𝜇$^{2}_{𝐵}$). We advocate for a more precise determination of the required expansion coefficients via lattice QCD simulations as a means of pinpointing the location of the critical endpoint in the phase diagram of strongly interacting matter.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

Single nucleotide variants drive evolutionary phage-host arms race in anaerobic carbon dioxide-converting microbiome

Microbial bioconversions are shaped by environmental perturbations and the adaptation of resident microbiomes. Prokaryotes coexist with bacteriophages, yet their coevolutionary trajectories remain underexplored. Here, we investigate the effects of a cultivation vessel leak on an anaerobic consortium performing carbon dioxide reduction. Using time-series shotgun metagenomic sequencing, we reconstruct microbial and viral genomes to track community shifts. We further apply single-nucleotide variant profiling and CRISPR array analysis to monitor viral microdiversity and host defense mechanisms. After bioaugmentation restores bioconversion efficiency, the consortium undergoes pronounced restructuring, with new dominant taxa emerging from the rare biosphere. We identify patterns consistent with phage predation selectively removing certain species, while others exhibit resilience to infection. This shift aligns with a widespread viral outbreak and a transient increased frequency of single nucleotide variants in bacterial CRISPR–Cas defense genes. Expansion of CRISPR spacers further supports that CRISPR-mediated processes influence microbial resilience. Concurrently, phages infecting resilient hosts exhibited adaptive evolution, marked by high genetic heterogeneity. Selective pressure varies across their genomes, targeting infectivity genes and protospacer-adjacent motifs. These findings highlight a dynamic evolutionary arms race driven by the selection of beneficial genetic variants, providing a mechanistic framework for multi-omics investigations, and informing biotechnological applications, including phage-based microbiome manipulation.

Ghiotto, G↗

DOE OSTI2026

GenASiS: General Astrophysical Simulation System. II. Self-gravitating Baryonic Matter*

GenASiS (General Astrophysical Simulation System) is a code being developed initially and primarily, though not exclusively, for the simulation of core-collapse supernovae on the world's leading capability supercomputers. This paper---the second in a series---documents capabilities for Newtonian self-gravitating fluid dynamics, including tabulated microphysical equations of state treating nuclei and nuclear matter (`baryonic matter'). Computation of the gravitational potential of a spheroid, and simulation of the gravitational collapse of dust and of an ideal fluid, provide tests of self-gravitation against known solutions. In multidimensional computations of the adiabatic collapse, bounce, and explosion of spherically symmetric pre-supernova progenitors---which we propose become a standard benchmark for code comparisons---we find that the explosions are prompt and remain spherically symmetric (as expected), with an average shock expansion speed and total kinetic energy that are inversely correlated with the progenitor mass at the onset of collapse and the compactness parameter.

Cardall, Christian [ORNL] (ORCID:000000020086105X)↗

DOE OSTI2026

Alcock-Paczynski Blinding Scheme for the Ly-$α$ Forest Analysis

We present and validate a blinding method for the Lyman-$α$ (Ly$α$) forest analysis based on a modification of the Alcock-Paczynski test. In order to hide the background expansion history, the method employs a geometrical shift of each quasar (QSO) forest in wavelength space, once the quasar continuum has been fitted and the fluctuation field is extracted. The redshift positions for the QSO sample are also changed in a consistent manner. We show that the method remains effective when applied to real data, where contamination from metals and Lyman-$β$ is intrinsically mixed with the Lyman-$α$ forest. This limitation is primarily visible in the 1D correlation function, where other blinding strategies can mitigate the effect. To assess its effectiveness, the prescription is tested against a series of datasets of increasing complexity: from idealized low-noise mocks, to realistic DESI year one synthetic datasets, and finally to data from DESI first data release (DR1), using both the auto (Ly$α\times$Ly$α$) and cross (Ly$α\times$ QSO) correlations. We find that the method robustly shifts the BAO peak position from the 3D correlation functions to the expected value for cosmology changes of around 5% in the matter content, without altering the shape of the posteriors in the model parameters. In conclusion, this catalog-level blinding strategy is a viable method for cosmological inference with the Lyman-$α$ forest, particularly if a cross-analysis with other tracers using the same blinding strategy is pursued.

Perez-Sanchez, G. [Guanajuato U.] (ORCID:000900096↗

DOE OSTI2025

Coupled Thermo-Hydrological-Mechanical-Chemical Behavior of Anisotropic Granite for Geologic Disposal of High-Level Radioactive Waste: A Core-Scale Laboratory Investigation

The coupled thermo-hydrological-mechanical-chemical (THMC) behavior of rock within an Excavation Damaged Zone (EDZ) is critical for the safety and long-term performance of a geological repository for high-level radioactive wastes. While many laboratory experiments have been conducted to investigate EDZ rocks, the flow and deformation characteristics resulting from anisotropic rock textures and microcrack distribution under triaxial loading and elevated temperatures remain poorly understood. Particularly, cracks at various scales serve as fast paths for fluid flow and solute transport and present as focal points of mechanical weakness, which complicate the coupled THMC processes in anisotropic EDZ rocks and challenge modeling predictions. Here, in this study, a series of core-scale experiments was conducted on three granite samples under repository-relevant conditions. These rock samples were obtained from the Grimsel Underground Research Laboratory (URL), featured by anisotropic minerals (represented by bedding layers) and microcrack distributions and coarse cm-scale grain sizes. During the experiments, samples were subjected to an elevated temperature at 90 °C and different triaxial loading conditions either by radial (normal to bedding layers) or axial (parallel to bedding layers) compaction. Water was injected into the samples, and the rock permeability evolutions and effluent water chemistry were monitored closely. For intact samples, thermal expansion of minerals at 90 °C resulted in a large, 75% irreversible permeability reduction and rock strengthening under radial compaction, while thermal impact was limited to a 15% permeability reduction under axial compaction. In contrast, for a sample containing open cracks, the growth of fractures during the experiment resulted in an abrupt permeability increase and fast failure at 90 °C. The effluent water chemistry indicates much more considerable mineral dissolution from large shear sliding than that in rocks dominated by mechanical compaction. These results helped better understand the coupled THMC processes in anisotropic rocks containing cracks, evaluate the behaviors of EDZ rocks, and predict the long-term evolution of EDZ for the performance of the repositories.

Coupled THMC processes↗

DOE OSTI2025

A parallel-kinetic-perpendicular-moment model for magnetised plasmas

We describe a new model for the study of weakly collisional, magnetised plasmas derived from exploiting the separation of the dynamics parallel and perpendicular to the magnetic field. This unique system of equations retains the particle dynamics parallel to the magnetic field while approximating the perpendicular dynamics through a spectral expansion in the perpendicular degrees of freedom, analogous to moment-based fluid approaches. In so doing, a hybrid approach is obtained that is computationally efficient enough to allow for larger-scale modelling of plasma systems while eliminating a source of difficulty in deriving fluid equations applicable to magnetised plasmas. We connect this system of equations to historical asymptotic models and discuss advantages and disadvantages of this approach, including the extension of this parallel-kinetic-perpendicular moment beyond the typical region of validity of these more traditional asymptotic models. This paper forms the first of a multi-part series on this new model, covering the theory and derivation, alongside demonstration benchmarks of this approach that include shocks and magnetic reconnection.

astrophysical plasmas↗

DOE OSTI2026

Ultrafast Relaxation of MLCT Excited States in Manganese Tricarbonyl Complexes: Insights from Polarization-Resolved Femtosecond X-ray Absorption Spectroscopy at the Mn and Br K-Edges

Optical absorption features that are often described as metal-to-ligand charge transfer (MLCT) bands underlie the utility of many metal coordination complexes by harnessing the energy of light to drive otherwise inaccessible chemical reactions. Excitation of these bands triggers rapid electronic dynamics that can be challenging to understand, due to complicated potential energy landscapes and highly correlated electronic and nuclear degrees of freedom in metal-containing compounds. The lowest-energy absorption bands in Mn complexes containing alpha-diimine, carbonyl, and halide ligands are particularly interesting, due to the delocalized nature of charge transfer upon excitation. Here, in this study, we report experimental ultrafast dynamics measurements for the series of compounds Mn(CO) 3 ( R bpy)Br ( R bpy = 4,4′-disubstituted-2,2′-bipyridine; R = H, CF 3 , or NO 2 ) using polarization-resolved, femtosecond X-ray absorption spectroscopy (XAS) at both the Mn and Br K edges. The appearance of a new absorption feature in the Br pre-edge spectrum upon optical excitation reveals the instantaneous formation of an electronic hole that is partially localized on the Br atom and has a lifetime that depends on the electron withdrawing character of the R bpy ligand. For two of the complexes (R = H or CF 3 ), a large expansion of the Mn–Br distance results in a rapid redistribution of the hole and a corresponding decrease in anisotropy of the absorption feature within 50 fs, after which the absorption into the hole disappears on a time scale shorter than 300 fs. We observe a different result for the NO 2 substituted compound, for which the Mn–Br bond contracts and the absorption into the Br-centered hole persists beyond the 2 ps time scale of our measurement. The Mn pre-edge spectrum also reveals structural changes for these three complexes, but the new Mn absorption features become evident only after the nuclei respond to the initial excitation, which allows mixing of Mn 3d and 4p orbitals. The combined use of ultrafast Mn and Br K-edge spectroscopy provides unique insight into the ways in which bipyridine substitution alters the excited-state dynamics, including a very different structural response for the most strongly electron withdrawing substituent.

Otolski, Christopher J. [Univ. of Kansas, Lawrence↗

DOE OSTI2024

MultiGreen: A multiplexing architecture for GreenGate cloning

Genetic modification of plants fundamentally relies upon customized vector designs. The ever-increasing complexity of transgenic constructs has led to increased adoption of modular cloning systems for their ease of use, cost effectiveness, and rapid prototyping. GreenGate is a modular cloning system catered specifically to designing bespoke, single transcriptional unit vectors for plant transformation—which is also its greatest flaw. MultiGreen seeks to address GreenGate’s limitations while maintaining the syntax of the original GreenGate kit. The primary limitations MultiGreen addresses are 1) multiplexing in series, 2) multiplexing in parallel, and 3) repeated cycling of transcriptional unit assembly through binary intermediates. MultiGreen efficiently concatenates bespoke transcriptional units using an additional suite of level 1acceptor vectors which serve as an assembly point for individual transcriptional units prior to final, level 2, condensation of multiple transcriptional units. Assembly with MultiGreen level 1 vectors scales at a maximal rate of 2*⌈ log 6 n ⌉+3 days per assembly, where n represents the number of transcriptional units. Further, MultiGreen level 1 acceptor vectors are binary vectors and can be used directly for plant transformation to further maximize prototyping speed. MultiGreen is a 1:1 expansion of the original GreenGate architecture’s grammar and has been demonstrated to efficiently assemble plasmids with multiple transcriptional units. MultiGreen has been validated by using a truncated violacein operon from Chromobacterium violaceum in bacteria and by deconstructing the RUBY reporter for in planta functional validation. MultiGreen currently supports many of our in-house multi transcriptional unit assemblies and will be a valuable strategy for more complex cloning projects.

Science & Technology - Other Topics↗

DOE OSTI2024

Techno-economic and life-cycle analysis of strategies for improving operability and biomass quality in catalytic fast pyrolysis of forest residues

Many of the challenges faced by the first commercial biorefineries were associated with feedstock handling, quality, and cost. Strategies are needed to enable further expansion of biorefineries and meet the growing demand for bio-based fuels and products. Here, we examine 2 key feedstock challenges and mitigation strategies in the context of a catalytic fast pyrolysis (CFP) biorefinery: (1) the operability of the feed system, which may be improved by modifying the minimum particle size fed to the reactor, and (2) the quality of the biomass, which may be improved by employing air classification to remove undesirable material and increase fuel yields. We conduct techno-economic analysis (TEA) and life-cycle analysis for these strategies, employing a discrete event simulation model for biomass preprocessing combined with a series of correlations developed from literature data and a rigorous CFP conversion model. Our results highlight the importance of balancing increased cost and material losses from preprocessing against improved operability and fuel yields. Economics and sustainability were optimized when operating at the lowest minimum particle size, emphasizing the importance of minimizing material losses while maintaining the operability of the process. Economically, additional costs and material losses from air classification could be acceptable due to improved biomass conversion, and an optimum air classification speed was identified; however, the fuel GHG emissions were minimized when air classification was not used. Valorizing material removed during preprocessing as a coproduct could improve economics and sustainability, decreasing the burden of material losses.

09 - BIOMASS FUELS↗

DOE OSTI2024

R-matrix calculations for opacities: I. Methodology and computations

Abstract An extended version of the R -matrix methodology is presented for calculation of radiative parameters for improved plasma opacities. Contrast and comparisons with existing methods primarily relying on the distorted wave approximation are discussed to verify accuracy and resolve outstanding issues, particularly with reference to the opacity project (OP). Among the improvements incorporated are: (i) large-scale Breit–Pauli R -matrix calculations for complex atomic systems including fine structure, (ii) convergent close coupling wave function expansions for the ( e + ion) system to compute oscillator strengths and photoionization cross sections, (iii) open and closed shell iron ions of interest in astrophysics and experiments, (iv) a treatment for plasma broadening of autoionizing resonances as function of energy-temperature-density dependent cross sections, (v) a ‘top-up’ procedure to compare convergence with R -matrix calculations for highly excited levels, and (vi) spectroscopic identification of resonances and bound ( e + ion) levels. The present R -matrix monochromatic opacity spectra are fundamentally different from OP and lead to enhanced Rosseland and Planck mean opacities. An outline of the work reported in other papers in this series and those in progress is presented. Based on the present re-examination of the OP work, opacities of heavy elements might require revisions in high temperature-density plasma sources.

Pradhan, A. K. (ORCID:0000000187753643)↗

DOE OSTI2024

Gold–Thiolate Nanocluster Dynamics and Intercluster Reactions Enabled by a Machine Learned Interatomic Potential

Monolayer protected metal clusters comprise a rich class of molecular systems and are promising candidate materials for a variety of applications. While a growing number of protected nanoclusters have been synthesized and characterized in crystalline forms, their dynamical behavior in solution, including prenucleation cluster formation, is not well understood due to limitations both in characterization and first-principles modeling techniques. Recent advancements in machine-learned interatomic potentials are rapidly enabling the study of complex interactions such as dynamical behavior and reactivity on the nanoscale. Here, we develop an Au-S-C-H atomic cluster expansion (ACE) interatomic potential for efficient and accurate molecular dynamics simulations of thiolate-protected gold nanoclusters (Au n (SCH 3 ) m ). Trained on more than 30,000 density functional theory calculations of gold nanoclusters, the interatomic potential exhibits ab initio level accuracy in energies and forces and replicates nanocluster dynamics including thermal vibration and chiral inversion. Long dynamics simulations (up to 0.1 μs time scale) reveal a mechanism explaining the thermal instability of neutral Au 25 (SR) 18 clusters. Specifically, we observe multiple stages of isomerization of the Au 25 (SR) 18 cluster, including a chiral isomer. Additionally, we simulate coalescence of two Au 25 (SR) 18 clusters and observe series of clusters where the formation mechanisms are critically mediated by ligand exchange in the form of [Au-S] n rings.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗