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

Radiation damage effects in beryllium for next generation neutrino beam targetry (Final Technical Report)

Current and future high-power accelerators put severe requirements on materials used for target and beam windows and target facilities have been recognized as a critical challenge in development of future particle accelerators. In accelerators, window and target materials are exposed to extreme conditions, which include bombardment with very high energy protons (1- 100 GeV) and thermomechanical shock waves. Radiation can cause direct damage in the material, and it leads to production of transmutation products (especially helium), both phenomena having a potential adverse effect on the stability and durability of the target/window material. At high enough temperatures, He can aggregate to form gas bubbles, which in turn cause significant dimensional changes (swelling), enable easy crack propagation, and eventually cause failure by fracture. On the other hand, if the temperature is too low, radiation damage accumulates in the form of internal defects (e.g., dislocations), leading to hardening and a decreased ductility of the material. In this project, we will focus on beryllium since it is considered to be one of the candidate materials for beam windows and targets in the next-generation proton accelerators, e.g., the Long Baseline Neutrino Facility (LBNF). Radiation effects in Be have been studied in the context of nuclear fusion reactor applications. However, key differences exist between reactor and accelerator conditions, including neutron vs. proton irradiation, continuous vs. pulsed beam flux, much higher energies of bombarding particles in accelerators, and higher operating temperatures for typical reactors. For example, the impact of beam pulsing on the radiation damage and the He bubble kinetics is largely unknown. While results obtained on Be from fusion research might not be directly transferrable to understanding target materials, there is an opportunity to bring state-of-the-art tools from materials research in nuclear reactors to aid design of target and beam window materials in high-power accelerators. To this end, the overarching goal of this project are to develop an experimentally-validated computational framework capable of predicting radiation damage evolution in beryllium relevant to beam window and target conditions, focusing on He bubble formation and growth as a function of irradiation temperature. Our model will be based on the cluster dynamics formalism, where size distribution of defects and He bubbles is simulated as a function of time, temperature, and radiation dose. Parameters for the model will be taken from published experiments and from high-fidelity atomistic simulations proposed in this project. In addition, we will carry out a series of targeted ex-situ and in-situ dual-beam experiments using low-energy protons to provide critical data for validation of the model on the effects of radiation on He clustering, He bubble distribution, and dislocation loop density/size in proton irradiated Be.

36 MATERIALS SCIENCE↗
DOE OSTI2021

First-Principles Grand-Canonical Simulations of Water Adsorption in Proton-Exchanged Zeolites

Water appears by design or as impurities in many important reactive systems. For those catalyzed by porous solid acids, such as widely used zeolites, experimentally quantifying the amount or elucidating the structure of adsorbed water clusters at reaction conditions is challenging, while computational studies (e.g., first-principles molecular dynamics simulations) often need to assume the loading to examine solvation effects. Furthermore, we perform first-principles grand-canonical simulations to predict water adsorption to H-ZSM-5 zeolites under specified experimental conditions. Presampling with inexpensive force fields and a pool-based parallelization algorithm are used to improve simulation efficiency, while molecular dynamics is used to sample configurations involving hydronium species. We observe that H + exchange dramatically increases the hydrophilicity of zeolite MFI and an appreciable amount of water is present at very low relative humidities. At all conditions examined, the zeolitic protons are found to dissociate readily from the surface basic sites and become mobile by participating in the hydrogen-bonded chains of adsorbed water molecules.

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

Direct observation of ultrafast cluster dynamics in supercritical carbon dioxide using X-ray Photon Correlation Spectroscopy

Supercritical fluids exhibit distinct thermodynamic and transport properties, making them of particular interest for a wide range of scientific and engineering applications. These anomalous properties emerge from structural heterogeneities due to the formation of molecular clusters at conditions above the critical point. While the static behavior of these clusters and their effects on the thermodynamic response functions have been recognized, the relation between the ultrafast cluster dynamics and transport properties remains elusive. By measuring the intermediate scattering function in carbon dioxide at conditions near the critical point with X-ray photon correlation spectroscopy, we directly capture the cross-over dynamics between 4 and 13 picoseconds, revealing the transition between ballistic and diffusive motion. Complementary analysis using large-scale molecular dynamics simulations reveals that this behavior arises from collisions between unbound molecules and clusters. This study provides direct evidence of the ultrafast momentum exchange between clusters, which has significant impact on transport properties, solvation processes, and reaction kinetics in supercritical fluids.

carbon capture and storage↗
DOE OSTI2024

Large-scale atomistic model construction of subbituminous and bituminous coals for solvent extraction simulations with reactive molecular dynamics

Large-scale atomistic models for complex polycyclic aromatic hydrocarbon systems help understand the chemical properties and behaviors of complex feedstocks such as coal or petroleum. However, the development and utilization of large-scale models remain limited due to the difficulty in achieving the varied structural characteristics necessary to capture stochastic nature of these feedstocks. Here we demonstrate a systematic workflow to construct stochastic molecular systems from a broad analytical suite: high-resolution transmission electron microscopy (HRTEM), carbon-13 nuclear magnetic resonance spectroscopy ( 13 C NMR), laser desorption ionization mass spectroscopy (LDI-MS), and elemental analysis. We present a model construction and analysis utility of a new Python-based module. We selected one subbituminous and three high-volatile bituminous coals to construct large-scale models (~40,000 atoms). The constructed models were utilized to examine the affinity for solvent extraction (naphthalene or tetralin) and the effect of structural properties (e.g., aromatic cluster size, functional groups, and cross-linking) in reactive molecular dynamics simulations. Complex chemical reactions were monitored with bond order transitions, intermediates formation, and mass distributions. Reactive molecular dynamics simulations suggest a plausible chemical extraction process and products for the complex fossil feedstocks. The results indicated that radical formations with bond breaking of bridging oxygens and carbons were required at high temperatures to facilitate hydrogeneration and extraction of gas molecules from radical-free molecules. We observed that aliphatic chains of tetralin were easily decomposed and combined with radicals to form small size of molecules with aryl bonding, mainly increasing molecules in the 500–1000 Da, while naphthalene had little impact on chemical extraction process.

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

Non-dilute helium-related defect interactions in the near-surface region of plasma-exposed tungsten

We report a systematic energetic analysis of helium-related defect interactions that mediate helium (He) segregation on surfaces of plasma-exposed tungsten at different levels of He ion implantation. We focus on high He fluences that increase the He content in the plasma-exposed material well beyond the dilute limit of He concentration and employ atomic configurations generated by large-scale molecular dynamics simulations of He-implanted tungsten. We perform systematic molecular statics computations of cluster– defect interaction energetics in the highly defect-rich near-surface region of plasma-exposed tungsten for small mobile helium clusters as a function of the cluster distance from the surface. In this region, mobile clusters are also subjected to the stress fields generated by defects such as helium bubbles and other clusters, which govern cluster–defect interactions in addition to the cluster–surface interaction. Based on systematic investigation of individual cluster–defect interactions, we develop a mathematical framework to describe the interaction energy landscapes consisting of elastic interaction potential perturbations to finite-width square-well potentials, where the potential well accounts for cluster trapping by the defect at close range and subsequent coalescence and the perturbation potential is parameterized according to elastic inclusion theory. Superposition of all the relevant interaction potentials provides a comprehensive description of the interaction energy landscape that would be experienced by a small mobile cluster along its migration path toward the plasma-exposed surface at high He fluence. Such descriptions are particularly important for developing atomistically-informed, hierarchical multi-scale models of helium cluster dynamics in plasma-facing materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Understanding fluorine-free electrolytes via small-angle X-ray scattering

Fluorine-free electrolytes have attracted great attention because of its low-cost and environmental friendliness. However, so far, little is known about the solution structures of these electrolytes. Here, we compare the solvation phenomenon of sodium tetraphenylborate (NaBPh 4 ) salt dissolved in organic solvents of propylene carbonate (PC), 1,2-dimethoxyethane (DME), acetonitrile (ACN) and tetrahydrofuran (THF). Small-angle X-ray scattering (SAXS) reveals a unique two-peak structural feature in this salt concentrated PC electrolyte, while solutions using other solvents only have one scattering peak. Molecular dynamics (MD) simulations further reveal that there are anion-based clusters in addition to the short-range charge ordering in the concentrated NaBPh4/PC electrolyte. Raman spectroscopy confirms the existence of considerable contact ion pairs (CIPs). Here, this work emphasizes the importance of global and local structural analysis, which will provide valuable clues for understanding the structure performance relationship of electrolytes.

25 ENERGY STORAGE↗
DOE OSTI2022

Machine learning enabled quantification of the hydrogen bonds inside the polyelectrolyte brush layer probed using all-atom molecular dynamics simulations

The configuration of densely grafted charged polyelectrolyte (PE) brushes is strongly dictated by the properties and behavior of the counterions that screen the PE brush charges and the solvent molecules (typically water) that solvate the brush molecules and these screening counterions. Only recently, efforts have been made to study the PE brushes atomistically, thereby shedding light on the properties of brush-supported ions and water molecules. However, even for such efforts, there are limitations associated with using a generic definition to estimate certain properties of water and ions inside the brush layer. For example, water–water hydrogen bonds (HBs) will behave differently for locations outside and inside the brush layer, given the fact that the densely closely grafted PE brush molecules create a soft nanoconfinement where the water connectivity becomes highly disrupted: therefore, using the same definition to quantify the HBs inside and outside the brush layer will be unwise. In this paper, we address this limitation by employing an unsupervised machine learning (ML) approach to predict the water–water hydrogen bonding inside a cationic PE brush layer modeled using all-atom molecular dynamics (MD) simulations. Here, the ML method, which relies on a clustering approach and uses the equilibrium coordinates of the water molecules (obtained from the all-atom MD simulations) as the input, is capable of identifying the structural modification of water–water HBs (revealed through appropriate clustering of the data) inside the PE brush layer induced soft nanoconfinement. Such capabilities would not have been possible by using a generic definition of the HBs. Our calculations lead to four key findings: (1) the clusters formed inside and outside the brush layer are structurally similar; (2) the margin of the cluster is shorter inside the PE brush layer confirming the possible disruption of the HBs inside the PE brush layer; (3) the average “hydrogen–acceptor-oxygen–donor-oxygen” angle that defines the HB is reduced for the HBs formed inside the brush layer; (4) the use of the generic definition (definition usable for characterizing the HBs in brush-free bulk) leads to an overprediction of the number of HBs formed inside the PE brush layer.

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

Structural origin of disorder-induced ion conduction in NaFePO 4 cathode materials

Diffusion in NaFePO 4 can be enhanced through amorphization. Based on computations using DFT and machine learning potentials, we ascribe this phenomenon to the formation of less constrained Na-ion environments upon disordering. Most modern battery technologies depend on solid-state crystalline cathode materials. However, some of these materials are constrained by the low ionic conductivity of their most stable phases. An example of this is maricite (NaFePO 4 ). Interestingly, experiments have shown that maricite can improve its rate capability through disordering (amorphization). However, experimental characterization of amorphous cathode materials remains a major challenge, hindering a clear understanding of the structural origin of the disorder-induced improvement in sodium-ion mobility. To address this, we here employ molecular dynamics simulations by first training a machine learning potential for NaFePO 4 based on the atomic cluster expansion approach and a batch active learning potential parameterization scheme. This potential is then applied to explore the structural and dynamical properties of NaFePO 4 glasses as cathode materials. Specifically, we investigate the effect of glass structure on sodium-ion diffusion, revealing the relative influences of short-range and medium-range order features. We find significant heterogeneity in sodium-ion diffusivity in the glass, with fast-conducting ions residing in less constrained atomic environments with fewer P and Fe neighbors. These more mobile ions are also surrounded by larger ring-type structures. Overall, the results and developed approach present promising avenues for developing high-performance glassy cathodes for next-generation batteries.

Christensen, Rasmus↗
DOE OSTI2020

Observations of grain-boundary phase transformations in an elemental metal

The theory of grain boundary (the interface between crystallites, GB) structure has a long history 1 and the concept of GBs undergoing phase transformations was proposed 50 years ago 2,3 . The underlying assumption was that multiple stable and metastable states exist for different GB orientations 4,5,6 . The terminology ‘complexion’ was recently proposed to distinguish between interfacial states that differ in any equilibrium thermodynamic property 7 . Different types of complexion and transitions between complexions have been characterized, mostly in binary or multicomponent systems 8,9,10,11,12,13,14,15,16,17,18,19 . Simulations have provided insight into the phase behaviour of interfaces and shown that GB transitions can occur in many material systems 20,21,22,23,24 . However, the direct experimental observation and transformation kinetics of GBs in an elemental metal have remained elusive. In this work we demonstrate atomic-scale GB phase coexistence and transformations at symmetric and asymmetric $[11\bar{1}]$ tilt GBs in elemental copper. Atomic-resolution imaging reveals the coexistence of two different structures at Σ19b GBs (where Σ19 is the density of coincident sites and b is a GB variant), in agreement with evolutionary GB structure search and clustering analysis 21,25,26 . We also use finite-temperature molecular dynamics simulations to explore the coexistence and transformation kinetics of these GB phases. Our results demonstrate how GB phases can be kinetically trapped, enabling atomic-scale room-temperature observations. Our work paves the way for atomic-scale in situ studies of metallic GB phase transformations, which were previously detected only indirectly 9,15,27,28,29 , through their influence on abnormal grain growth, non-Arrhenius-type diffusion or liquid metal embrittlement.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ultrafast visualization of phase transitions in nonequilibrium warm dense matter

Intense ultrafast laser excitation brings materials into highly nonequilibrium states with complex solid–liquid phase transitions. These ultrafast processes yield warm dense matter (WDM) conditions characterized by comparable thermal and Fermi energies, and strong ion–ion coupling. Here, we review recent studies employing mega-electron-volt ultrafast electron diffraction to resolve the structural dynamics of nonequilibrium WDM created by femtosecond laser heating of solids. Furthermore, for the study of warm dense gold, the results showed homogeneous melting that occurs within tens of picoseconds at absorbed energy densities of 0.4–1.4 MJ/kg. These results constrained the electron–ion coupling rate and revealed the melting sensitivity to nucleation seeds. For the study of radiation-damaged tungsten, the results showed a melting transition below the melting temperature. Molecular dynamics simulations suggested that this melting behavior is driven by vacancy defect clusters from radiation damage. These studies provided atomic-level insights into the melting behavior of materials under extreme conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The eROSITA Final Equatorial-Depth Survey (eFEDS)

In 2019, the eROSITA telescope on board the Russian-German satellite Spectrum-Roentgen-Gamma (SRG) began to perform a deep all-sky X-ray survey with the aim of identifying ~100 000 clusters and groups over the course of four years. As part of its performance verification phase, a ~140 deg 2 survey, called eROSITA Final Equatorial-Depth Survey (eFEDS), was performed. With a depth typical of the all-sky survey after four years, it allows tests of tools and methods as well as improved predictions for the all-sky survey. As part of this effort, a catalog of 542 X-ray selected galaxy group and cluster candidates was compiled. Here, in this paper, we present the optical follow-up, with the aim of providing redshifts and cluster confirmation for the full sample. Furthermore, we aim to provide additional information on the dynamical state, richness, and optical center of the clusters. Finally, we aim to evaluate the impact of optical cluster confirmation on the purity and completeness of the X-ray selected sample. We used optical imaging data from the Hyper Suprime-Cam Subaru Strategic Program and from the Legacy Survey to identify optical counterparts to the X-ray detected cluster candidates. We make use of the multi-component matched filter cluster confirmation tool (MCMF), as well as of the optical cluster finder CAMIRA to derive cluster redshifts and richnesses. MCMF provided the probabilities with which an optical structure would be a chance superposition with the X-ray candidate. These probabilities were used to identify the best optical counterpart as well as to confirm an X-ray candidate as a cluster. The impact of this confirmation process on catalog purity and completeness was estimated using optical to X-ray scaling relations as well as simulations. The resulting catalog was furthermore matched with public group and cluster catalogs. Optical estimators of the cluster dynamical state were constructed based on density maps of the red-sequence galaxies at the cluster redshift. By providing redshift estimates for all 542 candidates, we construct an optically confirmed sample of 477 clusters and groups with a residual contamination of 6%. Of these, 470 (98.5%) are confirmed using MCMF, and 7 systems are added through cross-matching with spectroscopic group catalogs. Using observable-to-observable scaling and the applied confirmation threshold, we predict that 8 ± 2 real systems have been excluded with the MCMF cut required to build this low-contamination sample. This number agrees well with the 7 systems found through cross-matching that were not confirmed with MCMF. The predicted redshift and mass distribution of this catalog agree well with simulations. Thus, we expect that these 477 systems include >99% of all true clusters in the candidate list. Using an MCMF-independent method, we confirm that the catalog contamination of the confirmed subsample is 6 ± 3%. Application of the same method to the full candidate list yields 17 ± 3%, consistent with estimates coming from the fraction of confirmed systems of ~17% and with expectations from simulations of ~20%. We also present a sample of merging cluster candidates based on the derived estimators of the cluster dynamical state.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Relationship of the Molecular Structure and Transport Properties of Imide-Based Lithium Salts of “Acetonitrile/Water-in-Salt” Electrolytes

“Water-in-salt” (WIS) electrolytes exhibit excellent safety and electrochemical performance. However, they possess high concentrations, relatively low diffusion coefficient, and high viscosity. “Acetonitrile/water in salt” (AWIS) electrolytes can overcome the disadvantages of WIS electrolytes. Under relatively low concentrations, AWIS electrolytes show good electrochemical performance comparable to WIS electrolytes and low conductivity. Herein, we investigate the relationship between the solvation structures and the transport properties using small-angle X-ray scattering and molecular dynamics simulation. We observed two solvation behaviors of AWIS: anions dissolved in acetonitrile forming small acetonitrile/anion clusters and additional water further dissolving the small acetonitrile/anion clusters. The introduction of acetonitrile weakens the water–solute interaction and enhances the cation–anion interaction, which results in an enhanced dynamical slowdown as the concentration increases. In conclusion, this work provides molecular-level understanding of the connection between two-stage solvation structures and transport properties for imide-based lithium salt solutions.

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

Effective Storage of Electrons in Water by the Formation of Highly Reduced Polyoxometalate Clusters

Aqueous solutions of polyoxometalates (POMs) have been shown to have potential as high-capacity energy storage materials due to their potential for multi-electron redox processes, yet the mechanism of reduction and practical limits are currently unknown. Herein, we explore the mechanism of multi-electron redox processes that allow the highly reduced POM clusters of the form {MO 3 } y to absorb y electrons in aqueous solution, focusing mechanistically on the Wells-Dawson structure X 6 [P 2 W 18 O 62 ], which comprises 18 metal centers and can uptake up to 18 electrons reversibly (y = 18) per cluster in aqueous solution when the countercations are lithium. This unconventional redox activity is rationalized by density functional theory, molecular dynamics simulations, UV-vis, electron paramagnetic resonance spectroscopy, and small-angle X-ray scattering spectra. These data point to a new phenomenon showing that cluster protonation and aggregation allow the formation of highly electron-rich meta-stable systems in aqueous solution, which produce H 2 when the solution is diluted. Finally, we show that this understanding is transferrable to other salts of [P 5 W 30 O 110 ] 15- and [P 8 W 48 O 184 ] 40- anions, which can be charged to 23 and 27 electrons per cluster, respectively.

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

DONKEY: A Flexible and Accurate Algorithm for Clustering

We propose an accurate clustering algorithm suitable for the varied and multidimensional data sets that correspond to temporal snapshots from on-the-fly nonadiabatic trajectory-based simulations of photoexcited dynamics. The algorithm approximates the underlying probability density function using variable kernel density estimation, with local maxima corresponding to cluster centers. Each data point is then assigned to one of the maxima by employing a maximization procedure. Finally, clusters artificially separated by minor fluctuations in the probability density are merged. The algorithm does not require parameter tuning, which ensures flexibility and reduces the risk of bias. It is tested on several synthetic data sets, where it consistently outperforms conventional clustering algorithms. As a final example, the algorithm is applied to the excited dynamics of the norbornadiene ⇌ quadricyclane (C 7 H 8 ) molecular photoswitch, demonstrating how distinct reaction pathways can be identified.

algorithms↗
DOE OSTI2024

From ionic clusters dynamics to network constraints in ionic polymer solutions

Physical networks formed by ionizable polymers with ionic clusters as crosslinks are controlled by coupled dynamics that transcend from ionic clusters through chain motion to macroscopic response. Here, the coupled dynamics, across length scales, from the ionic clusters to the networks in toluene swollen polystyrene sulfonate networks, were directly correlated, as the electrostatic environment of the physical crosslinks was altered. The multiscale insight is attained by coupling neutron spin echo measurements with molecular dynamics simulations, carried out to times typical of relaxation of polymers in solutions. The experimental dynamic structure factor is in outstanding agreement with the one calculated from computer simulations, as the networks are perturbed by elevating the temperature and changing the electrostatic environment. In toluene, the long-lived clusters remain stable over hundreds of ns across a broad temperature range, while the polymer network remains dynamic. In conclusion, though the size of the clusters changes as the dielectric constant of the solvent is modified through the addition of ethanol, they remain stable but morph, enhancing the polymer chain dynamics.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Accurate Fe–He machine learning potential for studying He effects in BCC-Fe

Nanostructured Ferritic Alloys (NFAs) containing oxide nanoparticles are prospective advanced structural materials in future fusion reactors. Lacking fusion neutron sources with an adequate flux, modeling is critical to understand radiation damage in the NFAs, including helium effects. Machine learning interatomic potentials (MLPs) have been demonstrated to be capable of describing atomic interactions in chemically complex systems with comparable accuracy to density functional theory (DFT) but with much less computational costs. Hence, MLPs are chosen to study He bubble accumulation in chemically complex NFAs. As a preliminary step, we develop a Fe–He potential, based on ~10,000 atomic configurations. The developed MLP accurately predicts the bulk properties of BCC-Fe compared to DFT, including the lattice constant (2.832 Å), elastic constants (c 11 = 271 GPa, c 12 = 141, c 44 = 93), and phonon frequencies (maximum error < 3.5%). The MLP predicts He binding energies in He n V bubbles and He n clusters with a mean absolute error (MAE) of only ~70 meV, which is ~3 to ~5 times smaller than the MAE of binding energies estimated using empirical potentials. The MLP also accurately predicts the migration barrier of interstitial He (64 meV). Furthermore, the MLP correctly predicts the maximum number of He atoms (critical size) at which a He n V bubble (He n cluster) can kick-out a self-interstitial atom (SIA), forming a bigger He n V 2 bubble (He n V bubble). Subsequently, we use the MLP to perform molecular dynamics simulations to investigate the effect of temperature on the critical size of He n V bubbles and He n clusters. We find that the critical size decreases with increasing temperature — indicating that smaller bubbles and clusters can kick-out an SIA at higher temperatures, thereby escalating the He bubble growth. We report the current Fe–He MLP can be further developed to include all chemical interactions in NFAs.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Elucidation of Agonist and Antagonist Dynamic Binding Patterns in ER-α by Integration of Molecular Docking, Molecular Dynamics Simulations and Quantum Mechanical Calculations

Estrogen receptor alpha (ERα) is a ligand-dependent transcriptional factor in the nuclear receptor superfamily. Many structures of ERα bound with agonists and antagonists have been determined. However, the dynamic binding patterns of agonists and antagonists in the binding site of ERα remains unclear. Therefore, we performed molecular docking, molecular dynamics (MD) simulations, and quantum mechanical calculations to elucidate agonist and antagonist dynamic binding patterns in ERα. 17β-estradiol (E2) and 4-hydroxytamoxifen (OHT) were docked in the ligand binding pockets of the agonist and antagonist bound ERα. The best complex conformations from molecular docking were subjected to 100 nanosecond MD simulations. Hierarchical clustering was conducted to group the structures in the trajectory from MD simulations. The representative structure from each cluster was selected to calculate the binding interaction energy value for elucidation of the dynamic binding patterns of agonists and antagonists in the binding site of ERα. The binding interaction energy analysis revealed that OHT binds ERα more tightly in the antagonist conformer, while E2 prefers the agonist conformer. The results may help identify ERα antagonists as drug candidates and facilitate risk assessment of chemicals through ER-mediated responses.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Ab Initio Molecular Dynamics Study of Pd Nucleation on a γ-Al 2 O 3 Surface

We perform ab initio molecular dynamics (AIMD) free energy simulations to investigate the stability and dynamics of dispersed Pd atoms and subnanometer clusters on γ-Al 2 O 3 . We consider nucleation for varying Pd loadings on dry γ-Al 2 O 3 (100), dry γ-Al 2 O 3 (110) and wet γ-Al 2 O 3 (110) and examine temperature, entropy and water coverage effects on the stability and dynamics of Pd clusters. At low Pd loading (1.07 Pd atoms/nm 2 on the (100) and 0.75 Pd/nm 2 on the (110)), association of the Pd atoms is thermodynamically unfavorable on the dry surfaces and almost thermoneutral on the wet (110) owing to water blocking of the AlIII defect sites. In agreement with past studies, the Pd atoms are quite mobile but the diffusion free energy barriers do not increase with the water coverage, as previously suggested, but rather remain similar in magnitude owing to translational entropy at finite temperatures. At these Pd loadings, oxidized Pd(II) atoms on the wet (110) surface are dispersed and immobile. At higher Pd loading (1.12 and 1.5 Pd atoms/nm 2 ) on the (110) surface, Pd 3 and Pd 4 clusters are thermodynamically very stable and form fast. At the 1.12 Pd atoms/nm 2 loading, the formation of the Pd 3 cluster is two orders of magnitude slower when the surface is hydroxylated as the system becomes somewhat trapped before the dinuclear cluster Pd 2 collides with a Pd atom.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗