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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Chemistry-mediated Ostwald ripening in carbon-rich C/O systems at extreme conditions

Abstract There is significant interest in establishing a capability for tailored synthesis of next-generation carbon-based nanomaterials due to their broad range of applications and high degree of tunability. High pressure (e.g., shockwave-driven) synthesis holds promise as an effective discovery method, but experimental challenges preclude elucidating the processes governing nanocarbon production from carbon-rich precursors that could otherwise guide efforts through the prohibitively expansive design space. Here we report findings from large scale atomistically-resolved simulations of carbon condensation from C/O mixtures subjected to extreme pressures and temperatures, made possible by machine-learned reactive interatomic potentials. We find that liquid nanocarbon formation follows classical growth kinetics driven by Ostwald ripening (i.e., growth of large clusters at the expense of shrinking small ones) and obeys dynamical scaling in a process mediated by carbon chemistry in the surrounding reactive fluid. The results provide direct insight into carbon condensation in a representative system and pave the way for its exploration in higher complexity organic materials. They also suggest that simulations using machine-learned interatomic potentials could eventually be employed as in-silico design tools for new nanomaterials.

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

Neuromorphic learning with Mott insulator NiO

Habituation and sensitization (nonassociative learning) are among the most fundamental forms of learning and memory behavior present in organisms that enable adaptation and learning in dynamic environments. Emulating such features of intelligence found in nature in the solid state can serve as inspiration for algorithmic simulations in artificial neural networks and potential use in neuromorphic computing. In this work, we demonstrate nonassociative learning with a prototypical Mott insulator, nickel oxide (NiO), under a variety of external stimuli at and above room temperature. Similar to biological species such as Aplysia, habituation and sensitization of NiO possess time-dependent plasticity relying on both strength and time interval between stimuli. A combination of experimental approaches and first-principles calculations reveals that such learning behavior of NiO results from dynamic modulation of its defect and electronic structure. An artificial neural network model inspired by such nonassociative learning is simulated to show advantages for an unsupervised clustering task in accuracy and reducing catastrophic interference, which could help mitigate the stability–plasticity dilemma. Mott insulators can therefore serve as building blocks to examine learning behavior noted in biology and inspire new learning algorithms for artificial intelligence.

42 ENGINEERING↗
NASA NTRS1994

Modeling the radial distribution of blue stragglers in M3

Recent surveys of the blue straggler (BS) population in the Galactic globular cluster M3 (NGC 5272) give the first complete characterization of the number density of BSs as a function of radius over an entire globular cluster. The BSs in M3 are overabundant at large radii and significantly underabundant at intermediate radii. Here we present the result of a simulation of the dynamical evolution of a population of BSs in a multimass model of M3. Assuming the BSs were formed in the core through binary interactions (Hut & Verbundt 1983; Leonard 1989; Sigurdsson & Phinney 1993; Hut et al. 1992b; Davies, Benz & Hills 1994), and given some very general assumptions about the recoil that occurs during stellar mergers in interacting binaries, we find an excellent fit to the observed radial distribution of BSs, suggesting strongly that most of the BSs in M3 were formed through binary collisions in the core.

Sigurdsson, Steinn↗
DOE OSTI2023

Graphic contrastive learning analyses of discontinuous molecular dynamics simulations: Study of protein folding upon adsorption

A comprehensive understanding of the interfacial behaviors of biomolecules holds great significance in the development of biomaterials and biosensing technologies. In this work, we used discontinuous molecular dynamics (DMD) simulations and graphic contrastive learning analysis to study the adsorption of ubiquitin protein on a graphene surface. Our high-throughput DMD simulations can explore the whole protein adsorption process including the protein structural evolution with sufficient accuracy. Contrastive learning was employed to train a protein contact map feature extractor aiming at generating contact map feature vectors. Subsequently, these features were grouped using the k-means clustering algorithm to identify the protein structural transition stages throughout the adsorption process. The machine learning analysis can illustrate the dynamics of protein structural changes, including the pathway and the rate-limiting step. Our study indicated that the protein–graphene surface hydrophobic interactions and the π–π stacking were crucial to the seven-stage adsorption process. Upon adsorption, the secondary structure and tertiary structure of ubiquitin disintegrated. The unfolding stages obtained by contrastive learning-based algorithm were not only consistent with the detailed analyses of protein structures but also provided more hidden information about the transition states and pathway of protein adsorption process and structural dynamics. Our combination of efficient DMD simulations and machine learning analysis could be a valuable approach to studying the interfacial behaviors of biomolecules.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A study of nuclear structure of light nuclei at the electron–ion collider

Understanding the substructure of atomic nuclei, particularly the clustering of nucleons inside them, is essential for comprehending nuclear dynamics. Various cluster configurations can emerge depending on excitation energy, the number and types of core clusters, and the presence of excess neutrons. Despite the prevalence of tightly bound cluster formations in low-lying states, understanding the correlation between clusters and their formation mechanisms remains incomplete. This exploring study investigates nuclear clustering at the electron–ion collider (EIC) using simulations based on the modified BeAGLE model. By simulating collisions involving e+ 9 Be, e+ 12 C, and e+ 16 O nuclei, we find that the average energy of particles and the system size ratios of particles at forward rapidity exhibit sensitivity to alpha clustering and its various configurations. These findings offer valuable insights into the dynamics of nuclear clustering and its implications for future studies at the EIC.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Mitigating Band Tailing in Kesterite Solar Absorbers: Ab Initio Quantum Dynamics

Open-circuit voltage deficits are limiting factors in kesterite solar cells. Addressing this issue by suppressing band tailing and nonradiative charge recombination is essential for enhancing the performance. We employ ab initio nonadiabatic molecular dynamics to elucidate the origin of band tailing and charge losses and propose a mitigation strategy. The simulations show that Cu–Zn disorder, associated with antisite defect clusters [Cu Zn +Zn Cu ], is a significant source of band tailing in kesterites, as evidenced by the much larger Urbach energy in disordered than ordered kesterites. Cu–Zn disorder gives rise to new sulfur-centered coordination polyhedra, increases structural inhomogeneity, changes electrostatic potential at sulfur centers, and shifts the S(3p) orbital energy. Differences in the S(3p)/Cu(3d) and S(3p)/Sn(5s) hybridization strengths and the S(3p) orbital energy shift reduce the band gap by 0.37 eV. Furthermore, Cu–Zn disorder enhances vibrational motion of sulfur anions and surrounding cations, increasing band gap fluctuations by 15 meV. The stronger electron–phonon interactions reduce charge carrier lifetimes and limit the kesterite solar cell efficiency. Partial substitution of Zn with Cd facilitates structural ordering and significantly suppresses band tailing, particularly in disordered systems. The improvement can be attributed to the larger atomic radius and mass of Cd, which weakens bonding around the anion, suppresses S-related vibrations within the covalent tetrahedra, and reduces nonadiabatic coupling, thereby increasing charge carrier lifetimes. The reported results establish the key influence of cation disorder on band tailing and reduced charge carrier lifetimes in kesterites and highlight cation disorder engineering as a strategy to achieve high-efficiency kesterite solar cells.

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

Structural adaptation of oxygen tolerance in 4-hydroxybutyrl-CoA dehydratase, a key enzyme of archaeal carbon fixation

Abstract Autotrophic microorganisms that convert inorganic carbon into organic matter were key players in the evolution of life on Earth. As the early atmosphere became oxygenated, microorganisms needed to develop mechanisms for oxygen protection, especially those relying on enzymes containing oxygen-sensitive metal clusters (e.g., Fe-S). Here we investigated how 4-hydroxybutyryl-CoA dehydratase (4HBD) - the key enzyme of the 3-hydroxypropionate/4-hydroxybutyrate (HP/HB) cycle for CO 2 -fixation - adapted as conditions shifted from anoxic to oxic. 4HBD is found in both anaerobic bacteria and aerobic ammonia-oxidizing archaea (AOA). The oxygen-sensitive bacterial 4HBD and oxygen-tolerant archaeal 4HBD share 59 % amino acid identity. To examine the structural basis of oxygen tolerance in archaeal 4HBD, we determined the atomic resolution structure of the enzyme. Two tunnels providing access to the canonical [4Fe-4S] cluster in oxygen-sensitive bacterial 4HBD were closed with four conserved mutations found in all aerobic AOA and other archaea. Further biochemical experiments and molecular dynamics simulations support our findings that restricting access to the active site is the key to oxygen tolerance, explaining how active site evolution drove a major evolutionary transition. Significance statement Autotrophy (primary production) was the first life strategy on Earth. Before photosynthesis (using solar energy to fix carbon dioxide), life relied on chemical reactions for energy. These chemosynthetic reactions are present in all domains of life, including archaea possessing the most energy-efficient carbon fixation pathway - the 3-hydroxypropionate/4-hydroxybutyrate cycle. This efficiency results from enzyme modifications, including enhanced enzyme stability and catalysis of multiple reactions. We reveal the first structure of aerobic 4-hydroxybutyryl-CoA dehydratase (4HBD) from ammonia-oxidizing archaea. These archaea are among the most abundant organisms on the planet, and their 4HBD active site evolved oxygen tolerance to support aerobic metabolism. This modification can provide further insight into enzyme evolution on early earth, as photosynthesis developed and began oxygenating the atmosphere.

DeMirci, Hasan (ORCID:0000000291355397)↗
DOE OSTI2022

Determination of Molecular Structure and Dynamics of Molten Salts by Advanced Neutron and X-ray Scattering Measurements and Computer Modeling

The design and development of fully functional Molten-Salt Reactors (MSR) require detailed knowledge of the molten salt properties in order to understand and predict the salt’s behavior. Fundamental properties of interest include molecular structure, speciation, and dynamics (such as diffusion coefficients) of salt components and dissolved corrosion and fission products. Computer modeling is necessary to predict changes in physical and chemical properties due to irradiation, burning of dissolved fuel, and corrosion. The modeling requires experimental data, and advanced neutron and x-ray scattering and spectroscopy provide the most reliable and direct determination of the structure (Pair-Distribution Functions, PDF), and dynamics of ions in the melt. This project dealt with both fluoride and chloride salts. The PDFs have been measured by a combination of neutron and x-ray diffraction. We utilized the techniques of isotope substitutions, a very powerful tool available for neutron-scattering, to extract the details of the liquid structure. Although similar measurements have been done before, modern advanced neutron and x-ray-scattering techniques allow collecting the data at much higher resolution and in a wider range of temperatures. Importantly, we were among the first to study fluoride salts by neutron scattering. The importance of impurities and their effects on salt properties have become apparent recently and so new methods of salt purification were developed. We took advantage of these developments to produce reliable data, which have been used for computer simulations of both clean salts and those with added fission and corrosion products most relevant for MSRs. Ab initio molecular dynamics simulations have been performed to understand the multi-component liquid solution, in particular solubility of impurities and thermodynamic interactions in relation to the ionic-cluster structure of the fluid. We applied machine learning to regress from the simulation and experimental data in order to develop a fast-acting model that can handle molten salt with an arbitrary (≥ 10) number of chemical elements and be able to predict chemical potential as a function of composition and temperature. This project resulted in a number of experimental and computer-simulation publications, a patent application, and numerous conference presentations (American Physical Society, American Chemical Society, and The Electrochemical Society among others). Multiple students and postdocs participated and collaborated on aspects of this project. This project seeded new collaborations between MIT and other institutions, such as the University of Massachusetts Lowell, the University of Illinois Urbana-Champaign, the University of California Berkeley, and Oak Ridge and Los Alamos National Labs. As such, this project has had a broad and lasting impact beyond its original scientific scope.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1984

When clusters are superclusters

Recent research on substructure in clusters and superclusters of galaxies is critically surveyed, with an emphasis on the distinction (based on physical size and dynamic state) between the two classes of structures. The cosmological and large-scale structural implications of substructure in rich clusters are analyzed, using N-body-simulation results, and illustrated for the clusters A98 and A548. Consideration is also given to gravitational lensing of distant quasars by intervening clusters and the interpretation of the matter distribution on scales larger than about 5 Mpc. It is determined that A98 and A548 are superclusters.

Geller, M. J.↗
NASA NTRS1993

An X-ray method for detecting substructure in galaxy clusters - Application to Perseus, A2256, Centaurus, Coma, and Sersic 40/6

We use the moments of the X-ray surface brightness distribution to constrain the dynamical state of a galaxy cluster. Using X-ray observations from the Einstein Observatory IPC, we measure the first moment FM, the ellipsoidal orientation angle, and the axial ratio at a sequence of radii in the cluster. We argue that a significant variation in the image centroid FM as a function of radius is evidence for a nonequilibrium feature in the intracluster medium (ICM) density distribution. In simple terms, centroid shifts indicate that the center of mass of the ICM varies with radius. This variation is a tracer of continuing dynamical evolution. For each cluster, we evaluate the significance of variations in the centroid of the IPC image by computing the same statistics on an ensemble of simulated cluster images. In producing these simulated images we include X-ray point source emission, telescope vignetting, Poisson noise, and characteristics of the IPC. Application of this new method to five Abell clusters reveals that the core of each one has significant substructure. In addition, we find significant variations in the orientation angle and the axial ratio for several of the clusters.

Mohr, Joseph J.↗
DOE OSTI2020

Scalable molecular dynamics on CPU and GPU architectures with NAMD

NAMD is a molecular dynamics program designed for high-performance simulations of very large biological objects on CPU- and GPU-based architectures. NAMD offers scalable performance on petascale parallel supercomputers consisting of hundreds of thousands of cores, as well as on inexpensive commodity clusters commonly found in academic environments. It is written in C++ and leans on Charm++ parallel objects for optimal performance on low-latency architectures. NAMD is a versatile, multipurpose code that gathers state-of-the-art algorithms to carry out simulations in apt thermodynamic ensembles, using the widely popular CHARMM, AMBER, OPLS, and GROMOS biomolecular force fields. Here, we review the main features of NAMD that allow both equilibrium and enhanced-sampling molecular dynamics simulations with numerical efficiency. We describe the underlying concepts utilized by NAMD and their implementation, most notably for handling long-range electrostatics; controlling the temperature, pressure, and pH; applying external potentials on tailored grids; leveraging massively parallel resources in multiple-copy simulations; and hybrid quantum-mechanical/molecular-mechanical descriptions. We detail the variety of options offered by NAMD for enhanced-sampling simulations aimed at determining free-energy differences of either alchemical or geometrical transformations and outline their applicability to specific problems. Last, we discuss the roadmap for the development of NAMD and our current efforts toward achieving optimal performance on GPU-based architectures, for pushing back the limitations that have prevented biologically realistic billion-atom objects to be fruitfully simulated, and for making large-scale simulations less expensive and easier to set up, run, and analyze. NAMD is distributed free of charge with its source code at www.ks.uiuc.edu.

high-performance computing↗
DOE OSTI2020

Identification of preferred multimodal ligand‐binding regions on IgG1 F C using nuclear magnetic resonance and molecular dynamics simulations

Abstract In this study, the binding of multimodal chromatographic ligands to the IgG1 F C domain were studied using nuclear magnetic resonance and molecular dynamics simulations. Nuclear magnetic resonance experiments carried out with chromatographic ligands and a perdeuterated 15 N‐labeled F C domain indicated that while single‐mode ion exchange ligands interacted very weakly throughout the F C surface, multimodal ligands containing negatively charged and aromatic moieties interacted with specific clusters of residues with relatively high affinity, forming distinct binding regions on the F C . The multimodal ligand‐binding sites on the F C were concentrated in the hinge region and near the interface of the C H 2 and C H 3 domains. Furthermore, the multimodal binding sites were primarily composed of positively charged, polar, and aliphatic residues in these regions, with histidine residues exhibiting some of the strongest binding affinities with the multimodal ligand. Interestingly, comparison of protein surface property data with ligand interaction sites indicated that the patch analysis on F C corroborated molecular‐level binding information obtained from the nuclear magnetic resonance experiments. Finally, molecular dynamics simulation results were shown to be qualitatively consistent with the nuclear magnetic resonance results and to provide further insights into the binding mechanisms. An important contribution to multimodal ligand‐F C binding in these preferred regions was shown to be electrostatic interactions and π–π stacking of surface‐exposed histidines with the ligands. This combined biophysical and simulation approach has provided a deeper molecular‐level understanding of multimodal ligand–F C interactions and sets the stage for future analyses of even more complex biotherapeutics.

Gudhka, Ronak B.↗
DOE OSTI2024

Structural and electronic characteristics of amorphous Ge 8 Sb 2 Te 11

GeTe-rich phase-change materials can be utilized in rewriteable optical memory due to the large contrast in reflectivity between amorphous and crystalline phases. Here we explored the structure and electronic properties of amorphous Ge 8 Sb 2 Te 11 using ab initio molecular dynamics simulations. The results indicate that amorphous Ge 8 Sb 2 Te 11 is dominantly composed of Ge-, Sb- and Te-centered octahedrons with distortions, while 30.4 % of Ge-centered clusters are in the form of tetrahedrons which are randomly distributed. The 5-fold rings possess a large proportion, and the Gesingle bondTe and Sbsingle bondTe bonds present larger formation energies than other bonds, leading to the ABAB bonding sequence (A: Ge and Sb, B: Te). The lone pair electrons locating at the opposite direction of bonds possess a large fraction of 14.8 %, which may enhance the distortions of local clusters. Importantly, these structural properties lead to the high stability of amorphous Ge 8 Sb 2 Te 11 and thus long data retention in the information storage.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Transport Mechanisms of Hydrogen Gas and Ions in Poly(potassium acrylate) Electrolytes

Poly(potassium acrylate) (PAAK) in high concentration KOH solutions is a promising electrolyte for use in alkaline zinc batteries. Atomistic molecular dynamics simulations are used to investigate the transport mechanisms of hydrogen gas and zincate ions through both KOH solutions and PAAK solutions with comparable concentrations of KOH. Hydrogen molecules do not interact with the ions or the polymer, and simply diffuse through the system. Furthermore, their diffusion constant decreases with both increasing KOH concentration and when in the polymer, due to the increased viscosity of the solution. By contrast, zincate ions form long-lived clusters with the PAAK chains which significantly slow their diffusion.

Batteries↗
DOE OSTI2024

The physical basis for solvent flow in organic solvent nanofiltration

Organic solvent nanofiltration (OSN) is an emerging membrane technology that could revolutionize chemical separations in numerous vital industries. Despite its significance, there remains a lack of fundamental understanding of solvent transport mechanisms in OSN membranes. Here, we use an extended Flory-Rehner theory, nonequilibrium molecular dynamic simulations, and organic solvent transport experiments to demonstrate that solvent flow in OSN membranes is driven by a pressure gradient. We show that solvent molecules migrate as clusters through interconnected pathways within the membrane pore structure, challenging the widely accepted diffusion-based view of solvent transport in OSN. We further reveal that solvent permeance is dependent on solvent affinity to the OSN membrane, which, in turn, controls the membrane pore structure. Our fundamental insights lay the scientific groundwork for the development of next-generation OSN membranes.

42 ENGINEERING↗
DOE OSTI2021

Understanding Solid Electrolyte Interphase Nucleation and Growth on Lithium Metal Surfaces

Experiments and theory are needed to decode the exact structure and distribution of components of a passivation layer formed at the anode surface of Li metal batteries, known as the Solid Electrolyte Interphase (SEI). Due to the inherent dynamic behavior as well as the lithium reactivity, the SEI structure and its growth mechanisms are still unclear. This study uses molecular simulation and computational chemistry tools to investigate the initial nucleation and growth dynamics of LiOH and Li2O that provide us with thermodynamics and structural information about the nucleating clusters of each species. Following the most favorable pathways for the addition of each of the components to a given nascent SEI cluster reveals their preferential nucleation mechanisms and illustrates different degrees of crystallinity and electron density distribution that are useful to understand ionic transport through SEI blocks.

25 ENERGY STORAGE↗
DOE OSTI2019

Quantum simulations of excited states with active-space downfolded Hamiltonians

Many-body techniques based on the double unitary coupled cluster ansatz (DUCC) can be used to downfold electronic Hamiltonians into low-dimensional active spaces. It can be shown that the resulting dimensionality reduced Hamiltonians are amenable for quantum computing. Recent studies performed for several benchmark systems using quantum phase estimation (QPE) algorithms demonstrated that these formulations can recover a significant portion of ground-state dynamical correlation effects that stem from the electron excitations outside of the active space. These results have also been confirmed in studies of ground-state potential energy surfaces using quantum simulators. In this letter, we study the effectiveness of the DUCC formalism in describing excited states. We also emphasize the role of the QPE formalism and its stochastic nature in discovering/identifying excited states or excited-state processes in situations when the knowledge about the true configurational structure of a sought after excited state is limited or postulated (due to the specific physics driving excited-state processes of interest). In this context, we can view the QPE algorithm as an engine for verifying various hypotheses for excited-state processes and providing statistically meaningful results that correspond to the electronic state(s) with the largest overlap with a postulated configurational structure. We illustrate these ideas on examples of strongly cor- related molecular systems, characterized by small energy gaps and high density of quasi-degenerate states around the Fermi level.

quantum chemistry, quantum computing, computationa↗
DOE OSTI2025

Simulations of Attosecond Metallization in Quartz and Diamond Probed with Inner-Shell Transient Absorption Spectroscopy

When dielectrics are hit with intense infrared (IR) laser pulses, transient metalization can occur. The initial attosecond dynamics behind this metallization are not entirely understood. Therefore, simulations are needed to understand this process and to help interpret experimental observations of it, such as with attosecond transient absorption (ATA). In this paper, we present first-principles simulations of ATA based on bulk-mimicking clusters and real-time time-dependent density functional theory (RT-TDDFT), with Koopmans-tuned range-separated hybrid functionals and Gaussian basis sets. Our method gives good agreement with the experiment for the breakdown threshold in silica and diamond. This breakdown voltage corresponds to a Keldysh parameter of approximately one and thus involves a transition to a regime where the dynamics are driven by tunneling. Pumping at an amplitude just below this value causes a mixture of multiphoton and tunneling excitations across the band gap to occur. The computed extreme ultraviolet and X-ray attosecond transient spectra also agree well with the experiment and show a decrease in optical density due to the transient population of the conduction band from the IR field. First-principles approaches such as this are valuable for interpreting the complicated modulations in a spectrum and for guiding future attosecond experiments on solids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗