Search NASA⌕ Search

SEARCH · Search NASA

Results for “Cluster Dynamics Simulations”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Sulfur-doped graphene anchoring of ultrafine Au 25 nanoclusters for electrocatalysis

The biggest challenge of exploring the catalytic properties of under-coordinated nanoclusters is the issue of stability. We demonstrate herein that chemical dopants on sulfur-doped graphene (S-G) can be utilized to stabilize ultrafine (sub-2 nm) Au 25 (PET)18 clusters to enable stable nitrogen reduction reaction (NRR) without significant structural degradation. The Au 25 @S-G exhibits an ammonia yield rate of 27.5μgNH 3 ∙ mg Au -1 ∙ h -1 at -0.5 V with faradic efficiency of 2.3%. More importantly, the anchored clusters preserve ~ 80% NRR activity after four days of continuous operation, a significant improvement over the 15% remaining ammonia production rate for clusters loaded on undoped graphene tested under the same conditions. Isotope labeling experiments confirmed the ammonia was a direct reaction product of N 2 feeding gas instead of other chemical contaminations. Ex-situ X-ray photoelectron spectroscopy and X-ray absorption near-edge spectroscopy of post-reaction catalysts reveal that the sulfur dopant plays a critical role in stabilizing the chemical state and coordination environment of Au atoms in clusters. Further ReaxFF molecular dynamics (RMD) simulation confirmed the strong interaction between Au nanoclusters (NCs) and S-G. Furthermore, this substrate-anchoring process could serve as an effective strategy to study ultrafine nanoclusters’ electrocatalytic behavior while minimizing the destruction of the under-coordinated surface motif under harsh electrochemical reaction conditions.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Bayesian prior construction for uncertainty quantification in first-principles statistical mechanics

First-principles statistical mechanics enables the prediction of thermodynamic and kinetic properties of materials, but is computationally expensive. Many approaches require surrogate models to calculate energies within Monte Carlo or molecular dynamics simulations. Inexpensive surrogates such as cluster expansions enable otherwise intractable calculations by interpolating data from higher accuracy methods, such as Density Functional Theory (DFT). Surrogate models introduce uncertainty into downstream calculations, in addition to any uncertainty inherent to DFT calculations. Bayesian frameworks address this by quantifying uncertainty and incorporating expert knowledge through priors. However, constructing effective priors remains challenging. This work introduces and describes practical strategies for building Bayesian cluster expansions, focusing on basis truncation, hyperparameter selection, and ground state replication. We analyze multiple basis truncation schemes, compare cross-validation to the evidence-approximation for hyperparameter optimization, and provide methods to find and enforce ground-state-preserving models through priors. Additionally, we compare the uncertainties between different approximations to DFT (LDA, PBE, SCAN) against the uncertainty introduced with the use of cluster expansion surrogate models. These approaches are demonstrated on the BCC Li x Mg 1-x and Li x Al 1-x alloys, which are both of interest for solid-state Li batteries. Our results provide guidelines for constructing and utilizing Bayesian cluster expansions, thereby improving the transparency of materials modeling. Furthermore, the approaches and insights developed in this work can be transferred to a wide range of cluster expansion surrogate models, including the atomic cluster expansion and related machine-learned interatomic potential architectures.

Alloy theory↗
DOE OSTI2023

Structural mapping and tuning of mixed halide ions in amorphous sulfides for fast Li-ion conduction and high deformability

Amorphous sulfides are among the most promising candidates for solid electrolytes (SEs) owing to their excellent deformability and acceptable Li-ion conductivity (σion) at room temperature. However, the complex atomic structure of these amorphous materials without long-range ordering results in a lack of structural understanding and difficulties in tuning material properties. Here, in this study, we performed structural mapping of glassy sulfide SEs composed of various anion clusters using a combinatorial atomic level analysis of synchrotron X-ray-based pair distribution function (PDF) and reverse Monte-Carlo (RMC) methods, demonstrating the potential to tune σ ion in glass SEs. First, we prepared a new glassy sulfide with a mixed anion framework of two halogens (Br and I) as atomic anions, with PS 4 3- molecular anions in a Li 2 S–P 2 S 5 -based glass SE. At a specific Br content, [(Li 2 S) 0.658 (LiI 0.9 LiBr 0.1 ) 0.342 ] 0.825 [P 2 S 5 ] 0.175 recorded a σ ion of 2.27 mS cm -1 , the highest value for any glass sulfides reported to date; however, its elastic modulus was still suppressed to 14.48 GPa. PDF and RMC calculations successfully provided structural mapping of anion clusters, including two halogens. Molecular dynamics simulations of each composition confirmed that flexible coordination caused by the rattling of small polarizable Br ions in the mixed halogens of glassy SEs contributed to the superior σ ion . Our results may provide new insights into the design of superior glassy SEs that play key roles in all-solid-state batteries requiring fast Li-ion conduction and high deformability.

36 MATERIALS SCIENCE↗
NASA NTRS1994

Star cluster evolution with primordial binaries. 3: Effect of the Galactic tidal field

We present the results of N-body simulations of tidally limited star clusters with an initial population of 0%-20% binaries. We find that (1) if enough binaries are initially present, the binary fraction may fall to a minimum value, then increase at late times; (2) the cluster evaporation timescale is quite insensitive to the details of the initial binary distribution; (3) the cluster core radius stabilizes at a few percent of the half-mass radius when binaries are present, just as in the case of isolated clusters; and (4) there may be a marked difference between the spatial distribution of low-energy and high-energy binaries as the cluster evolves. Specifically, the spatial distribution of the lower energy systems is often substantially more extended than that of the more tightly bound pairs. At no time are our simulated clusters well described by simple dynamical models that neglect the close coupling between the binding energies and the center-of-mass energies of the binaries they contain.

Mcmillan, Steve↗
DOE OSTI2020

Neon Cluster Formation and Phase Separation during White Dwarf Cooling

Recent observations of Galactic white dwarfs (WDs) with Gaia suggest there is a population of massive crystallizing WDs exhibiting anomalous cooling—the Q branch. While single-particle 22 Ne sedimentation has long been considered a possible heat source, recent work suggests that 22 Ne must separate into clusters, enhancing diffusion, in order for sedimentation to provide heating on the observed timescale. We show definitively that 22 Ne cannot separate to form clusters in C/O WDs using molecular dynamics simulations, and we further present a general C/O/Ne phase diagram showing that strong 22Ne enrichment is not achievable for 22 Ne abundance ≲30%. Here, we conclude that the anomalous heating cannot be due to 22Ne cluster sedimentation and that Q branch WDs may have an unusual composition, possibly rich with heavier elements.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Multi-site reaction dynamics through multi-fragment density matrix embedding

The practical description of disordered chemical reactions, where the reactions involve multiple species at multiple sites, is presently a challenge using correlated electronic structure methods due to their high computational cost and steep scaling. Here, we describe the gradient theory of multi-fragment density matrix embedding theory, which potentially provides a minimal computational framework to model such processes at the correlated electron level. We present the derivation and implementation of the gradient theory, its validation on model systems and chemical reactions using density matrix embedding, and its application to a molecular dynamics simulation of proton transport in a small water cluster, a simple example of multi-site reaction dynamics.

Chemistry↗
DOE OSTI2021

27 Al NMR diffusometry of Al 13 Keggin nanoclusters

Although nanometer-sized aluminum hydroxide clusters (i.e., ϵ-Al 13 , [Al 13 O 4 (OH) 24 (H 2 O) 12 ] 7+ ) command a central role in aluminum ion speciation and transformations between minerals, measurement of their translational diffusion is often limited to indirect methods. Here, in this work, 27 Al pulsed field gradient stimulated echo nuclear magnetic resonance (PFGSTE NMR) spectroscopy has been applied to the AlO 4 core of the ϵ-Al 13 cluster with complementary theoretical simulations of the diffusion coefficient and corresponding hydrodynamic radii from a boundary element-based calculation. The tetrahedral AlO 4 center of the ϵ-Al 13 cluster is symmetric and exhibits only weak quadrupolar coupling, which results in favorable T 1 and T 2 27 Al NMR relaxation coefficients for 27 Al PFGSTE NMR studies. Stokes–Einstein relationship was used to relate the 27 Al diffusion coefficient of the ϵ-Al 13 cluster to the hydrodynamic radius for comparison with theoretical simulations, dynamic light scattering from literature, and previously published 1 H PFGSTE NMR studies of chelated Keggin clusters. This first-of-its-kind observation proves that 27 Al PFGSTE NMR diffusometry can probe symmetric Al environments in polynuclear clusters of greater molecular weight than previously considered.

27Al↗
DOE OSTI2023

Time-dependent equation-of-motion coupled-cluster simulations with a defective Hamiltonian

Here simulations of laser-induced electron dynamics in a molecular system are performed using time-dependent (TD) equation-of-motion (EOM) coupled-cluster (CC) theory. The target system has been chosen to highlight potential shortcomings of truncated TD-EOM-CC methods [represented in this work by TD-EOM-CC with single and double excitations (TD-EOM-CCSD)], where unphysical spectroscopic features can emerge. Specifically, we explore driven resonant electronic excitations in magnesium fluoride in the proximity of an avoided crossing. Near the avoided crossing, the CCSD similarity-transformed Hamiltonian is defective, meaning that it has complex eigenvalues, and oscillator strengths may take on negative values. When an external field is applied to drive transitions to states exhibiting these traits, unphysical dynamics are observed. For example, the stationary states that make up the time-dependent state acquire populations that can be negative, exceed one, or even complex-valued.

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

Efficient machine learning interatomic potentials robust for liquid and multiple solid polymorphs of NaF and KF

Achieving atomic-level understanding of crystallization of molten salts is of importance to a wide range of technological applications. Recent work [Fan et al., Proc. Natl. Acad. Sci. USA 122, e2425702122 (2025)] revealed that crystal nucleation in molten LiF salt is a multistage process according to the molecular-dynamics (MD) simulations based on an atomic cluster expansion (ACE) machine-learning interatomic potential (MLIP). In order to understand the influence of increasing cation size on nucleation pathways and nucleation rates of molten fluoride salts, here we develop two new ACE MLIPs for NaF and KF. The two ACE MLIPs feature DFT-SCAN-level accuracy for liquid and multiple solid polymorphs over a wide temperature (0–2000 K) and pressure (0–100 GPa) range, and also reproduce well a number of experimental data for solid and liquid equilibrium properties. The efficiency of the two ACE MLIPs enable million-atom-scale or microsecond-scale MD simulations. The two general-purpose ACE MLIPs are expected to be useful for atomistic simulations for different purposes, in addition to studying crystallization of molten NaF and KF salts.

Crystal melting↗
DOE OSTI2024

Suppressed Cosmic-Ray Energy Densities in Molecular Clouds from Streaming Instability-regulated Transport

Abstract Cosmic rays (CRs) are the primary driver of ionization in star-forming molecular clouds (MCs). Despite their potential impacts on gas dynamics and chemistry, no simulations of star cluster formation following the creation of individual stars have included explicit cosmic-ray transport (CRT) to date. We conduct the first numerical simulations following the collapse of a 2000 M ⊙ MC and the subsequent star formation including CRT using the STAR FORmation in Gaseous Environments framework implemented in the GIZMO code. We show that when CRT is streaming-dominated, the CR energy in the cloud is strongly attenuated due to energy losses from the streaming instability. Consequently, in a Milky Way–like environment the median CR ionization rate in the cloud is low ( ζ ≲ 2 × 10 −19 s −1 ) during the main star-forming epoch of the calculation and the impact of CRs on the star formation in the cloud is limited. However, in high-CR environments, the CR distribution in the cloud is elevated ( ζ ≲ 6 × 10 −18 ), and the relatively higher CR pressure outside the cloud causes slightly earlier cloud collapse and increases the star formation efficiency by 50% to ∼13%. The initial mass function is similar in all cases except with possible variations in a high-CR environment. Further studies are needed to explain the range of ionization rates observed in MCs and explore star formation in extreme CR environments.

Astronomy & Astrophysics↗
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↗
DOE OSTI2025

Molecular Origin of Nanoscale Anion Ordering of LiTFSI Electrolytes Revealed through SAXS/WAXS and Molecular Dynamics Simulations

Recent developments in “water-in-salt” electrolytes have precipitated a renewed effort to study imide-based electrolytes. While previous small-/wide-angle X-ray scattering (SAXS/WAXS) studies have attributed the emergence of a low-Q peak in the SAXS profile of aqueous LiTFSI electrolytes to nanometer-scale anion clustering, a molecular-level understanding of the root of these clusters remains unclear. Here, in this study, we combined molecular dynamics simulations and SAXS/WAXS to study the solvation structures of LiTFSI in acetonitrile, methanol, and water. We concluded that hydrogen bonding in water and MeOH stabilizes anion clusters, while nonpolar methyl groups on methanol and acetonitrile interrupt the nanoscale ordering of TFSI anions. This causes LiTFSI in water and MeOH electrolytes to exhibit two low-Q SAXS profile peaks while LiTFSI in acetonitrile exhibits only a single peak below Q = 1 Å –1 . These findings shed light on the underlying molecular origins of nanoscale anion clusters, which may help in the design of the next generation of electrolyte chemistries.

Trojanowski, Lucas [Univ. of Michigan, Ann Arbor, ↗
DOE OSTI2021

Physics-based, neural network force fields for reactive molecular dynamics: Investigation of carbene formation from [EMIM + ][OAc - ]

Reactive molecular dynamics simulations enable detailed understanding of solvent effects on chemical reaction mechanisms and reaction rates. While classical molecular dynamics using reactive force fields allows significantly longer simulation time scales and larger system sizes compared with ab initio molecular dynamics, constructing reactive force fields is a difficult and complex task. In this work, we describe a general approach following the Empirical Valence Bond (EVB) framework for constructing ab initio reactive force fields for condensed phase simulations by combining physics-based methods with neural networks (PB/NN). The physics-based terms ensure correct asymptotic behavior of electrostatic, polarization, and dispersion interactions, and are compatible with existing solvent force fields. Neural networks are utilized for versatile description of short-range orbital interactions within the transition state region, and accurate rendering of vibrational motion of the reacting complex. Herein, we demonstrate our methodology for a simple deprotonation reaction of the 1-ethyl-3-methylimidazolium (EMIM+) cation with acetate to form 1-ethyl- 3-methylimidazol-2-ylidene and acetic acid. Our PB/NN force field exhibits ~ 1 kJ/mol MAE accuracy within the transition state region for the gas-phase complex. To characterize solvent modulation of the reaction profile, we compute potentials of mean force (PMFs) for the gas-phase reaction as well as the reaction within a four ion cluster, and benchmark against ab initio molecular dynamics simulations. We find that the surrounding ionic environment significantly destabilizes formation of the carbene product, and we show that this effect is accurately captured by the reactive force field. By construction, the PB/NN potential may be directly employed for simulations of other solvents/chemical environments without additional parameterization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2000

Construction and Utilization of a Beowulf Computing Cluster: A User's Perspective

Lockheed Martin Space Operations - Stennis Programs (LMSO) at the John C Stennis Space Center (NASA/SSC) has designed and built a Beowulf computer cluster which is owned by NASA/SSC and operated by LMSO. The design and construction of the cluster are detailed in this paper. The cluster is currently used for Computational Fluid Dynamics (CFD) simulations. The CFD codes in use and their applications are discussed. Examples of some of the work are also presented. Performance benchmark studies have been conducted for the CFD codes being run on the cluster. The results of two of the studies are presented and discussed. The cluster is not currently being utilized to its full potential; therefore, plans are underway to add more capabilities. These include the addition of structural, thermal, fluid, and acoustic Finite Element Analysis codes as well as real-time data acquisition and processing during test operations at NASA/SSC. These plans are discussed as well.

Woods, Judy L.↗
DOE OSTI2022

Adaptive Conformer Sampling for Property Prediction Using the Conductor-like Screening Model for Real Solvents

The valorization of lignocellulose-derived bioproducts requires effective separation from excessive water. Liquid–liquid extraction is a promising low-energy separation technology, but effective extraction requires solvent selection based on the thermodynamic properties of the bioproduct and solvent components. We propose a computational framework for predicting such properties by developing an adaptive conformer selection approach for use with COSMO-RS (conductor-like screening model for real solvents) calculations. In this framework, molecular dynamics simulations are used to generate many molecular structures (conformers) at representative temperatures in varying solvent environments. Conformers are then clustered based on structural metrics in a low-dimensional space and selected using a mixed-integer quadratic programming problem to iteratively insert a sampled conformer. At each iteration, we determine bioproduct properties using COSMO-RS. Here, we demonstrate the capability of the proposed framework on representative bioproducts to show convergence of the adaptive sampling toward experimentally measured properties with fewer calculations than required by random conformer sampling, enabling the improved screening of solvent systems for liquid-phase separation.

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

Influence of nanoscale interfaces on the dynamic deformation and spall failure of Cu–Fe alloy microstructures

Additively manufactured immiscible Cu–Fe alloys can exhibit a distribution of nanoscale interfaces due to the distribution of nanoscale clusters in the equilibrium and metastable phases. Molecular dynamics simulations investigate the role of such interfaces on the phase stability and transformation behavior during shock compression, as well as the mechanisms of damage nucleation during spall failure. The model multiphase Cu–Fe systems studied here comprise a distribution of Fe clusters in an FCC Cu matrix, as well as Cu clusters in a BCC Fe matrix. The length scales of the nanoscale clusters determine the energetics of the interfaces that can result in FCC/BCC (equilibrium) or FCC/FCC (metastable) interfaces in the Cu system matrix, and BCC/FCC (equilibrium) or BCC/BCC (metastable) interfaces in the Fe matrix. The MD simulations demonstrate that nanoscale metastable interface microstructures can induce plastic deformation and also stabilize phases or suppress phase transformations in metastable phase clusters. In contrast, equilibrium interfaces can influence phase transformation thresholds and serve as additional void nucleation sites during failure.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Functional Design of Peptide Materials Based on Supramolecular Cohesion

Peptide materials offer a broad platform to design biomimetic soft matter, and filamentous networks that emulate those in extracellular matrices and the cytoskeleton are among the important targets. Given the vast sequence space, a combination of computational approaches and readily accessible experimental techniques is required to design peptide materials efficiently. Here, we report here on a strategy that utilizes this combination to predict supramolecular cohesion within filaments of peptide amphiphiles, a property recently linked to supramolecular dynamics and consequently bioactivity. Using established coarse-grained simulations on 10,000 randomly generated peptide sequences, we identified 3500 likely to self-assemble in water into nanoscale filaments. Atomistic simulations of small clusters were used to further analyze this subset of sequences and identify mathematical descriptors that are predictive of intermolecular cohesion, which was the main purpose of this work. We arbitrarily selected a small cohort of these sequences for chemical synthesis and verified their fiber morphology. With further characterization, we were able to link the latent heat associated with fiber to micelle transitions, an indicator of cohesion and potential supramolecular dynamicity within the filaments, to calculated hydrogen bond densities in the simulation clusters. Based on validation from in situ synchrotron X-ray scattering and differential scanning calorimetry, we conclude that the phase transitions can be easily observed by very simple polarized light microscopy experiments. We are encouraged by the methodology explored here as a relatively low-cost and fast way to design potential functions of peptide materials.

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

Unique Conductivity Behavior in Water-In-Salt Electrolytes Driven by Ion Clusters

Understanding and predicting ion transport in aqueous electrolytes are crucial for advanced energy storage and biophysics, and many emergent technologies yet remain elusive. Herein, we introduce a unified framework to quantitatively describe and predict electrolyte conductivity that shifts from conventional molar concentration-based metrics to a volume fraction-based approach. Through analyzing a variety of electrolyte solutions via this perspective, we observe a universal conductivity peak at a 37% volume fraction. Small-angle X-ray scattering (SAXS) and molecular dynamics (MD) simulations reveal that nanometer-scale ion clusters drive this general behavior. Moreover, key geometric features of the ion transport pathwayssuch as pore size, tortuosity, and connectivityfollow a consistent dependence with respect to the volume fraction, reinforcing the argument for the universal conductivity trend. This paradigm shift opens new avenues for designing high-performance electrolytes and provides transformative insights for advancing studies in many fields, wherein molecular aggregates dictate transport properties.

Nguyen, Huong T. D.↗