Search NASA⌕ Search

SEARCH · Search NASA

Results for “molecular cluster structure”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Revealing the Structural Evolution of Green Rust Synthesized in Ionic Liquids by In Situ Molecular Imaging

Abstract Ionic liquids are green solvents that have wide applications in material synthesis, catalysis, and separation. A model switchable ionic liquid (SWIL) consisting of 1,8‐diazabicycloundec‐7‐ene (DBU) and 1‐hexanol with carbon dioxide (CO 2 ) as the trigger gas is chosen to synthesize nanocrystalline green rust. Under anoxic conditions, by adding iron (II) acetate (Fe(C 2 H 3 O 2 ) 2 ) and methanol to the degassed SWIL, a nanoparticulate green rust with carbonate (nano‐GR) is synthesized. More importantly, the molecular structural change of the ionic liquid resulting in green rust crystallines is first characterized using in situ liquid using time‐of‐flight secondary ion mass spectrometry (ToF‐SIMS) coupled with the system for analysis at the liquid vacuum interface (SALVI), a vacuum compatible microfluidic reactor in this work. Dimers and ionic clusters consisting ferric and ferrous ions are identified in liquid ToF‐SIMS spectra. Spectral principal component analysis is used to confirm that these ion pairs including dimers and cluster ions are different from nonionic liquids. The results show that liquid ToF‐SIMS can be a useful tool to study complex liquids at the molecular level providing insights into predicative synthesis of nanomaterials using environmentally friendly green solvents and giving unique visualization of the evolving material interface due to nucleation.

Shen, Yanjie↗
DOE OSTI2025

Ionic Associations and Hydration in the Electrical Double Layer of Water-in-Salt Electrolytes

Water-in-Salt-Electrolytes (WiSEs) are an exciting class of concentrated electrolytes finding applications in energy storage devices because of their expanded electrochemical stability window, good conductivity and cation transference number, and fire-extinguishing properties. These distinct properties are thought to originate from the presence of an anion-dominated ionic network and interpenetrating water channels for cation transport, which indicates that associations in WiSEs are crucial to understanding their properties. Currently, associations have mainly been investigated in the bulk, while little attention has been given to the electrolyte structure near electrified interfaces. Here, we develop a theory for the electrical double layer (EDL) of WiSEs, where we consistently account for the thermoreversible associations of species into Cayley tree aggregates. The theory predicts an asymmetric structure of the EDL. At negative voltages, hydrated Li + dominates, and cluster aggregation is initially slightly enhanced before disintegration at larger voltages. At positive voltages, when compared to the bulk, clusters are strictly diminished. Performing atomistic molecular dynamics (MD) simulations of the EDL of WiSE provides EDL data for validation and bulk data for parametrization of our theory. Validating the predictions of our theory against MD showed good qualitative agreement. Furthermore, we performed electrochemical impedance measurements to determine the differential capacitance of the studied LiTFSI WiSE and also found reasonable agreement with our theory. Overall, the developed approach can be used to investigate ionic aggregation and solvation effects in the EDL, which, among other properties, can be used to understand the precursors for solid-electrolyte interphase formation.

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

Bespoke Liquid/Liquid Interfaces (Final Technical Report)

The goal of DE-SC0001815 was to advance the basic science of liquid:liquid interface formation, to develop a deeper understanding of the mechanisms of phase separation and the essential relationships between solution composition, organization and dynamics that underlie the kinetic regime of solvent extraction. This included learning how interfacial organization and dynamics alters the properties of the primary coordination sphere of ions and the free energy of transport of ions complexes across a phase boundary. We relied primarily upon classical molecular dynamics studies to determine the equilibrium ensembles of these complex systems, but also utilized ab-initio MD and cluster-based density functional theory (DFT) calculations when more detailed investigation of the electronic structure was needed. We continued development of graph-theory based analyses to elucidate hierarchical correlations and expanded into geometric topology methods to quantify the collectively organized structures that can organize at a liquid/liquid interface during solute transport. One of the main conclusions was from the observation of two distinct mechanisms for solute transport - those that derive from amplifications of interfacial heterogeneity and surface roughness, and those wherein surface roughness has been dampened and instead collectively organized macrostructures work to bring solutes into the organic phase. It was our aim to create a concrete chemical model of the underlying driving forces behind interfacial primary and secondary structure formation and to map out the energetic features of solute transport so that tailored liquid/liquid can be developed that have characteristic kinetic features associated with mass transport.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

A High-Fidelity Molecular Model of the Cu(111) Repeating Unit

Dynamic processes at surfaces are central to heterogeneous catalysis, but their atomistic mechanism(s) can prove difficult to elucidate due to variations in material structure and the corresponding impact on reactivity. Moreover, disparities between reaction conditions and those employed for spectroscopic characterization at surfaces can inhibit detailed understanding of catalysis-relevant chemistries. Herein, we substantiate the so-called “cluster-surface” analogy by leveraging a low-valent tricopper architecture ( 1 ) as a model system for small molecule activation at Cu(111). Two reaction classes are explored: the adsorption of carbon monoxide (CO) and the dissociative adsorption of dihydrogen (H 2 ). These processes serve as an ideal testbed to compare the reactivity of a molecular cluster ( 1 ) to that of a heterogeneous surface, as both reactions have empirical data from measurements performed on crystalline Cu(111). Cluster 1 reversibly binds CO. Variable temperature NMR analysis with 13 CO reveals a favorable enthalpy but large negative entropy (−5.1 kcal × mol –1 and −22.9 cal × mol –1 × K –1 , respectively) for CO binding, affording a process that is marginally endergonic at room temperature (ΔG ads (298.15 K) = 1.7 ± 0.5 kcal × mol –1 ). Similarly, analogous to a Cu(111) surface, 1 is shown to oxidatively add (chemisorb) H 2 . Kinetic parameters were determined for this process and the activation enthalpy (8.4 ± 0.5 kcal × mol –1 ) closely mirrors that established for H 2 binding at the Cu(111) facet (6.0 to 12.4 kcal × mol –1 ). Together, these results showcase that a trinuclear cluster can reproduce the small molecule binding and activation energetics of a bulk crystalline surface, setting the stage for studying less-defined surface processes in an atomically precise molecular setting.

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

Nucleation Rate Theory for Coordination Number: Elucidating Water-Mediated Formation of a Zigzag Na 2 SO 4 Morphology

Predicting and controlling nanostructure formation during nucleation can pave the way to synthesizing novel energy materials via crystallization. However, such control over nucleation and crystallization remains challenging due to an inadequate understanding of critical factors that govern evolving atomistic structures and dynamics. Here, in this paper, we utilize coordination number as a reaction coordinate and rate theory to investigate how sodium sulfate, commonly known as a phase-change energy material, nucleates in a supersaturated aqueous solution. In conjunction with ab initio and force field-based molecular dynamics simulation, the rate theoretical analysis reveals that sodium sulfate from an initially dissolved metastable state transits to a heterogeneous mixture of prenucleated clusters and finally to a large cylindrical zigzag morphology. Measurements of Raman spectra and their ab initio modeling confirm that this nucleated morphology contains a few waters for every sulfate. Rate processes such as solvent exchange and desolvation exhibit high sensitivity to the evolving prenucleation/nucleation structures, providing a means to distinguish between critical nucleation precursors. Desolvation and forming the first-shell interionic coordination structure via monomer-by-monomer addition around sulfates are found to explain the formation of large nuclei. Thus, a detailed understanding of the step-by-step structure formation across scales has been achieved. This can be leveraged to predict nucleation-related structures and dynamics and potentially control the synthesis of novel phase-change materials for energy applications.

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

Simulation of polymerization induced phase separation in model thermosets

Polymerization induced phase separation (PIPS) in a three component thermoset is studied using molecular dynamics simulations of a new coarse-grained thermoset model. The system includes two crosslinker molecules, which differ in their glass transition temperatures (T g ) and chain length and thus have the potential for phase separation. Here, one crosslinker has a high T g corresponding to a rubbery behavior, and simulations were performed for a short length (4 beads) and a long length (33 beads). The resin and other crosslinker have low T g . A coarse-grained model is developed with these features and with interaction parameters determined so that for either rubbery crosslinker length, the system is in the liquid state at the cure temperature. For sufficiently slow reaction rates, the long rubbery molecule exhibits PIPS into a bicontinuous array of nanoscale domains, but the short one does not, reproducing recent experimental results. The simulations demonstrate that the reaction rates must be slow enough to allow diffusion to yield phase separation. Particularly, the reaction rate corresponding to the secondary amine must be very slow, else the structure of crosslinked clusters and the substantially increased diffusion time will prevent PIPS.

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

Introducing metal–sulfur active sites in metal–organic frameworks via post-synthetic modification for hydrogenation catalysis

Metal–sulfur active sites play a central role in catalytic processes such as hydrogenation and dehydrogenation, yet the majority of active sites in these compounds reside on the surfaces and edges of catalyst particles, limiting overall efficiency. Here we present a strategy to embed metal–sulfur active sites into metal–organic frameworks (MOFs) by converting bridging or terminal chloride ligands into hydroxide and subsequently into sulfide groups through post-synthetic modification. We apply this method to two representative MOF families: one featuring one-dimensional metal–chloride chains and another containing discrete multinuclear metal clusters. Crystallographic and spectroscopic analyses confirm structural integrity and sulfide incorporation, and the transformation is monitored by in situ total scattering methods. The sulfided MOFs display enhanced catalytic activity in the selective hydrogenation of nitroarenes using molecular hydrogen. Density functional theory calculations indicate that sulfur incorporation promotes homolytic metal–ligand bond cleavage and facilitates H 2 activation. This work establishes an approach to construct MOFs featuring accessible metal–sulfide sites.

Catalyst synthesis↗
NASA NTRS1985

Simulation studies for surfaces and materials strength

From intermolecular force studies, it is now known that the overall non-additive contribution to the lattice enegy is positive so that analysis based on only pairwise additivity suggests a shallower intermolecular potential than the true value. Two body contributions alone are also known to be categorically unable to even qualitatively describe some configurations of molecular clusters in the gas phase or the general relaxation and reconstruction of fcc crystal surfaces. In addition, the many-body contribution was shown to play a key role in the stability of certain crystal structures. In these recent analyses, a relatively simple potential energy function (PEF), comprising only a two-body Mie-type potential plus a three-body Axilrod-Teller-type potential, was found to be extremely effective. This same parametric PEF is applied to describe the bulk stability and surface energy for the diamond cubic structure. To test the stability condition, the FCC, BCC, diamond cubic, graphite and beta-tin structures were considered.

Halicioglu, T.↗
DOE OSTI2022

Reversible motions and disordered structure of soft particles in amorphous solids

In amorphous solids, soft vibrational modes derived from normal mode analysis can be utilized to identify the soft particles that are prone to irreversible rearrangements. However, the normal mode analysis approach cannot explain why the spatial distributions of clustered soft particles do not change over time. We define a softness parameter based on the vibrational density of states calculated directly from molecular dynamics simulations with both the harmonic vibrations and anharmonic relaxations being properly captured at finite temperatures. This parameter spontaneously correlates with the real space atomic motions and the dynamics heterogeneity. Using the softness parameter, we show that the softest particles are confined within rigid cages. These particles keep rearranging reversibly without long-range diffusion. The moderately soft particles rearrange irreversibly, and the hard particles mainly participate in vibrations without rearrangement. Here, we also show that the soft particles form locally disordered structures, while the hard particles present strong ordering. These findings confirm the defective nature of soft particles, and provide insights on the nature of softness as the ability to rearrange, but not necessarily irreversibly.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Using Density-Corrected DFT to Understand Density-Driven and Functional-Dependent Errors in Ab Initio Simulations of the Hydrated Electron

The hydrated electron, an excess electron in liquid water, plays a crucial role in a plethora of chemical processes, motivating extensive research efforts to characterize its structure, dynamics, and reactivity in solution. Recent theoretical approaches to understanding this intriguing object have involved ab initio simulations based on density functional theory (DFT). Although DFT allows for the study of hydrated electron reactivity and quantum mechanical behavior, it is well-known that anionic systems can suffer from significant density-driven errors (DDEs). Density-corrected DFT (DC-DFT) provides a framework to mitigate such errors; the method reduces DDEs by replacing the self-consistent (SC) density associated with a given density functional with the Hartree–Fock (HF) density. Since HF densities tend to be more localized than DFT SC densities, the DC-DFT scheme significantly improves errors in calculations where the SC density is spuriously delocalized. Here, we investigate how the use of density correction affects the calculated properties of the DFT-simulated (PBEh) hydrated electron, a particularly challenging diffuse anionic system to simulate. First, we analyze charge delocalization in a system consisting of a model octahedral hydrated electron water cluster (the so-called Kevan structure) along with a spatially separated sulfur atom. We show that the use of density correction indeed reduces DDEs in comparison to a standard DFT global hybrid functional. We then propagate molecular dynamics trajectories of the hydrated electron using DC-DFT, where we find that DC further localizes electron density in the cavity region, a signature of reduced charge delocalization. Unfortunately, the decreased radius of gyration of the spin density and corresponding tightening of the local solvation structure from density correction causes predicted observables to deviate further from experimental measurements than when density correction is not employed. Here, we argue that DC’s worse agreement with experiment results from the removal of a fortuitous cancellation of errors that is intrinsic to the PBEh functional. This indicates that the difficulties with DFT to simulate hydrated electrons are primarily due to the inherent approximations in DFT rather than to density-driven errors.

Density functional theory↗
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 OSTI2025

MBX V1.2: Accelerating Data-Driven Many-Body Molecular Dynamics Simulations

The MBX software provides an advanced platform for molecular dynamics simulations, leveraging state-of-the-art MB-pol and MB-nrg data-driven many-body potential energy functions. Developed over the past decade, these potential energy functions integrate physics-based and machine-learned many-body terms trained on electronic structure data calculated at the "gold standard" coupled-cluster level of theory. Recent advancements in MBX have focused on optimizing its performance, resulting in the release of MBX v1.2. While the inherently many-body nature of MB-pol and MB-nrg ensures high accuracy, it poses computational challenges. MBX v1.2 addresses these challenges with significant performance improvements, including enhanced parallelism that fully harnesses the power of modern multicore CPUs. In conclusion, these advancements enable simulations on nanosecond time scales for condensed-phase systems, significantly expanding the scope of high-accuracy, predictive simulations of complex molecular systems powered by data-driven many-body potential energy functions.

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

Metallization and spin fluctuations in Cu-doped lead apatite

An electronic structure and magnetic properties analysis of the recently proposed Cu-doped lead apatite is performed. Here, we show that electronic structures of differently Cu-substituted structures are characterized by localized molecular Cu-O bands at or near the Fermi level. The Cu substitutions can happen at both Pb1 and Pb2 sites, leading to metallic and semiconducting systems differently. The electronic systems in these bands are highly unstable magnetically and form clusters of rigidly ferromagnetically coupled magnetic moments on Cu and neighboring oxygen atoms with a total moment of about 1µ B . The ground state of uniformly Cu-doped lead apatite appears to be magnetic and semiconducting. The nonuniform distribution of two Cu atoms at the nearest Pb2 sites leads to an antiferromagnetic semiconducting state with formation energy close to uniformly distributed Cu configurations. The inclusion of quantum spin fluctuations confirms the stability of magnetic Cu-O clusters. Our calculations revealed the absence of the long-range magnetic order between uniformly distributed Cu-O clusters, creating the spin glass type of system.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Unravelling the mechanism of pressure induced polyamorphic transition in an inorganic molecular glass

The atomic structure of a germanium doped phosphorous selenide glass of composition Ge 2.8 P 57.7 Se 39.5 is determined as a function of pressure from ambient to 24 GPa using Monte-Carlo simulations constrained by high energy x-ray scattering data. The ambient pressure structure consists primarily of P 4 Se 3 molecules and planar edge shared phosphorus rings, reminiscent of those found in red phosphorous as well as a small fraction of locally clustered corner-sharing GeSe 4 tetrahedra. This low-density amorphous phase transforms into a high-density amorphous phase at ~6.3 GPa. The high-pressure phase is characterized by an extended network structure. The polyamorphic transformation between these two phases involves opening of the P 3 ring at the base of the P 4 Se 3 molecules and subsequent reaction with red phosphorus type moieties to produce a cross linked structure. The compression mechanism of the low-density phase involves increased molecular packing, whereas that of the high pressure phase involves an increase in the nearest-neighbor coordination number while the bond angle distributions broaden and shift to smaller angles. The entropy and volume changes associated with this polyamorphic transformation are positive and negative, respectively, and consequently the corresponding Clapeyron slope for this transition would be negative. This result has far reaching implications in our current understanding of the thermodynamics of polyamorphic transitions in glasses and glass-forming liquids.

42 ENGINEERING↗
DOE OSTI2020

Effects of surface charge and cluster size on the electrochemical dissolution of platinum nanoparticles using COMB3 and continuum electrolyte models

We study the site-dependent dissolution of platinum nanoparticles under electrochemical conditions to assess their thermodynamic stability as a function of shape and size using empirical molecular dynamics and electronic-structure models. The third-generation charge optimized many-body potential is employed to determine the validity of uniform spherical representations of the nanoparticles in predicting dissolution potentials (the Kelvin model). To understand the early stages of catalyst dissolution, implicit solvation techniques based on the self-consistent continuum solvation method are applied. It is demonstrated that interfacial charge and polarization can shift the dissolution energies by amounts on the order of 0.74 eV depending on the surface site and nanoparticle shape, leading to the unexpected preferential removal of platinum cations from highly coordinated sites in some cases.

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

Short- and medium-range orders in Al90Tb10 glass and their relation to the structures of competing crystalline phases

Molecular dynamics simulations using an interatomic potential developed by artificial neural network deep machine learning are performed to study the local structural order in Al90Tb10 metallic glass. We show that more than 80% of the Tb-centered clusters in Al90Tb10 glass have short-range order (SRO) with their 17 first coordination shell atoms stacked in a ‘3661’ or ‘15551’ sequence. Medium-range order (MRO) in Bergman-type packing extended out to the second and third coordination shells is also clearly observed. Analysis of the network formed by the ‘3661’ and ‘15551’ clusters show that ~82% of such SRO units share their faces or vertexes, while only ~6% of neighboring SRO pairs are interpenetrating. Such a network topology is consistent with the Bergman-type MRO around the Tb-centers. Moreover, crystal structure searches using genetic algorithm and the neural network interatomic potential reveal several low-energy metastable crystalline structures in the composition range close to Al 90 Tb 10 . Some of these crystalline structures have the ‘3661’ SRO while others have the ‘15551’ SRO. While the crystalline structures with the ‘3661’ SRO also exhibit the MRO very similar to that observed in the glass, the ones with the ‘15551’ SRO have very different atomic packing in the second and third shells around the Tb centers from that of the Bergman-type MRO observed in the glassy phase.

36 MATERIALS SCIENCE↗
NASA NTRS1992

Cool stars, stellar systems, and the sun; Proceedings of the 7th Cambridge Workshop, Tucson, AZ, Oct. 9-12, 1991

Consideration is given to HST observations of late-type stars, molecular absorption in the UV spectrum of Alpha Ori, EUV emission from late-type stars, Rosat observations of the Pleiades cluster, a deep ROSAT observation of the Hyades cluster, optical spectroscopy detected by EXOSAT, stellar photospheric convection, a structure of the solar X-ray corona, magnetic surface images of the BY Dra Star HD 82558, a Zebra interpretatin of BY Dra stars, optical flares on II Peg, a low-resolution spectroscopic survey of post-T tauri candidates, millimeter and sub-millimeter emission from flare stars, and activity in tidally interacting binaries. Attention is also given to modeling stellar angular momentum evolution, extended 60-micron emission from nearby Mira variables, the PANDORA atmosphere program, the global properties of active regions, oscillations in a stratified atmosphere, lithium abundances in northern RS CVn binaries, a new catalog of cool dwarf stars, the Far UV Spectrograph Explorer, and development of reflecting coronagraphs.

Giampapa, Mark S.↗
DOE OSTI2021

CASTELO: clustered atom subtypes aided lead optimization—a combined machine learning and molecular modeling method

Background: Drug discovery is a multi-stage process that comprises two costly major steps: pre-clinical research and clinical trials. Among its stages, lead optimization easily consumes more than half of the pre-clinical budget. We propose a combined machine learning and molecular modeling approach that partially automates lead optimization workflow in silico, providing suggestions for modification hot spots. Results: The initial data collection is achieved with physics-based molecular dynamics simulation. Contact matrices are calculated as the preliminary features extracted from the simulations. To take advantage of the temporal information from the simulations, we enhanced contact matrices data with temporal dynamism representation, which are then modeled with unsupervised convolutional variational autoencoder (CVAE). Finally, conventional and CVAE-based clustering methods are compared with metrics to rank the submolecular structures and propose potential candidates for lead optimization. Conclusion: With no need for extensive structure-activity data, our method provides new hints for drug modification hotspots which can be used to improve drug potency and reduce the lead optimization time. It can potentially become a valuable tool for medicinal chemists.

59 BASIC BIOLOGICAL SCIENCES↗