SEARCH · Search NASA
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.
Schematic illustration of the hybrid method model. ρ is the net charge density calculated from cDFT. v is the velocity calculated from non-equilibrium MD simulation. b and α W are the slipping length and the surface properties, respectively.
MA'AT conformational models of the phi torsion angles of O -glycosidic linkages differ from those obtained from MD simulation.
Molecular dynamics (MD) simulations to understand the thermodynamic, dynamic, and structural changes in supercritical water across the Frenkel line and the melting line have been performed. The two-phase thermodynamic model [J. Phys. Chem. B(2010), 114(24), 8191-8198] and the velocity autocorrelation functions are used to locate the Frenkel line and to calculate the thermodynamic and dynamic properties. The Frenkel lines obtained from the two-phase thermodynamic model and the velocity autocorrelation criterion do not agree with each other. Structural characteristics and the translational diffusion dynamics of water suggest that this inconsistency could arise from the two oscillatory modes in water, which are associated with the bending of hydrogen bonds and intermolecular collisions inside the first coordination shell. Finally, the overall results lead us to conclude that the universality of the Frenkel line as a dynamic crossover line from rigid to nonrigid fluids is preserved in water.
An electric field applied across the interface has been shown to enable transitions from the Cassie to the Wenzel state on superhydrophobic surfaces with miniature corrugations. Molecular dynamics (MD) simulations manifest the possibility of reversible cycling between the two states when narrow surface wells support spontaneous expulsion of water in the absence of the field. With approximately 1 nm sized wells between the surface asperities, the response times to changes in the electric field are of O(0.1) ns, allowing up to GHz frequency of the cycle. Because of the orientation preferences of interfacial water in contact with the solid, the phenomenon depends on the polarity of the field normal to the interface. The threshold field strength for the Cassie-to-Wenzel transition is significantly lower for the field pointing from the aqueous phase to the surface; however, once in the Wenzel state, the opposite field direction secures tighter filling of the wells. Considerable hysteresis revealed by the delayed water retraction at decreasing field strength indicates the presence of moderate kinetic barriers to expulsion. Known to scale approximately with the square of the length scale of the corrugations, these barriers preclude the use of increased corrugation sizes while the reduction of the well diameter necessitates stronger electric fields. Field-controlled Cassie-to-Wenzel transitions are therefore optimized by using superhydrophobic surfaces with nanosized corrugations. Abrupt changes indicate a high degree of cooperativity reflecting the correlations between the wetting states of interconnected wells on the textured surface.
Ion conducting capability is often imparted to polymeric materials through short polyether side-chains, and yet the impact of this graft polymer architecture on ion solvation and conduction has not been fully explored. Here, we use a combination of impedance spectroscopy, vibrational spectroscopy, and atomistic molecular dynamics (MD) to compare the conductivity, ionic interactions, and polymer dynamics in a series of graft polyether electrolytes. We find that in poly[(oligo ethylene oxide)methyl ether methacrylate] (POEM), a widely used graft polymer electrolyte, the ionic conductivity drops more than an order of magnitude as the side-chain length is decreased from nine ethylene oxide (EO) units to three. This difference in conductivity is unexplained by differences in the calorimetric glass transition temperature (T g ), which varies only slightly with side-chain length. Furthermore, through vibrational spectroscopy and MD simulations we demonstrate that both linear and graft polyethers solvate Li⁺ ions effectively and dissociate them from large counterions, irrespective of side-chain length. Li⁺ ions do, however, show preferential solvation by EO units far from the methacrylate backbone. Similarly, EO units far from the backbone show enhanced segmental dynamics, while those near the immobile methacrylate group move substantially more slowly, as quantified by bond vector autocorrelation relaxation times. This heterogeneity in both ion solvation and local segmental relaxation explains variation in ion conductivity where material-averaged properties such as T g and number of free ions fail to do so. Importantly, the ionic conductivity is dictated primarily by the segmental mobility of the EO units which form effective solvation sites, rather than system-wide dynamics.
We report a systematic diffraction study of two “water-in-salt” electrolytes and a “water-in-bisalt” electrolyte combining high-energy X-ray diffraction (HEXRD) with polarized and unpolarized neutron diffraction (ND) on both H 2 O and D 2 O solutions. The measurements provide three independent combinations of correlations between the different pairs of atom types that reveal the short- and intermediate-range order in considerable detail. The ND interference functions show pronounced peaks around a scattering vector Q ~ 0.5 Å -1 that change dramatically with composition, indicating significant rearrangements of the water network on a length scale around 12 Å. The experimental results are compared with two sets of Molecular Dynamics (MD) simulations, one including polarization effects and the other based on a non-polarizable force field. The two simulations reproduce the general shapes of the experimental structure factors and their changes with concentration, but differ in many detailed respects, suggesting ways in which their force fields might be modified to better represent the actual systems.
Tetrahydropyran (THP) is a five-carbon heterocyclic ether that is non-carcinogenic, non-peroxide forming, biodegradable, and economically competitive with tetrahydrofuran (THF) as a solvent. In this work, THP has been synthesized from renewable biomass at >99.8% selectivity and 98% yield via hydrogenation of furfural-derived 3,4-dihydropyran (DHP) over Ni/SiO 2 in a continuous flow reactor at 150–200 °C. The apparent activation energy of THP formation is 31 kJ mol -1 , and the reaction orders with respect to the partial pressures of H 2 , DHP, and THP are: 2, 1, and -0.3. The kinetic data has been fitted to a Hougen–Watson model where the rate limiting step is the hydrogenation of adsorbed DHP. Ni/SiO 2 is shown to have a low deactivation rate constant of 0.012 h -1 over 100 h time on stream and can be regenerated in situ. As a performance advantage, THP is shown to be resistant to ring opening polymerization under strongly acidic conditions that THF is not, revealing it to be a superior solvent. Further, the minimum selling price of THP is competitive with the market price of THF ($$900 – 1400 per ton ) at a DHP feedstock cost of $1000 per ton. Conductor-like Screening Model for Real Solvents (COSMO-RS) and molecular dynamics (MD) simulations with 1008 solvents and 8 common plastics have demonstrated that THP can serve as an alternative solvent to THF, 2-methyltetrahydrofuran (MeTHF), and cyclopentyl methyl ether (CPME) for plastic dissolution, especially low-density polyethylene (LDPE), polypropylene (PP), polystyrene (PS) and polyvinyl chloride (PVC). This work establishes THP as a green solvent with excellent thermal, chemical and peroxidative stability that can be used for numerous applications, including waste plastic recycling.
Electrochemical measurements, nonlinear spectroscopy, and MD simulations are used to explore how the sizes and identities of ions influence their non-uniform distribution in acetonitrile that is highly organized due to proximity to a polar interface.
Elucidating the influence of nanoscale confinement on the freezing behavior of salt solutions is of fundamental interest to environmental security and materials science and engineering. Specific structural information such as the coordination environment of ions in confined aqueous medium, effects on the extended hexagonal network of water molecules and their mutual non-bonding interactions in confinement are sparse in the literature along with the change in the dynamical characteristics of water in the presence of ions. To address these knowledge gaps, the current study is focused on investigating the influence of reduced dimensionality arising from nanoscale confinement on the structural evolution of salt solutions on freezing and the associated fluid–surface interactions. In this regard, operando wide-angle X-ray scattering (WAXS) measurements and classical molecular dynamics (MD) simulations are conducted with water and 0.5 M CaCl 2 , MgCl 2 and KCl solutions confined in 4 nm sized SBA-15 silica pores upon cooling from 300 K to 200 K. The freezing point of the salt solutions is depressed to 235 K which is about 10 K lower than that of confined water. The translational dynamics of confined salt solutions indicates shifts in the fragile-to-strong dynamical crossover to a lower temperature compared to pure water. The strong electrostatic attraction between the cations and the surrounding water molecules contributes to the freezing point depression in confined salt solutions. These insights unlock the molecular-scale basis and mechanisms underlying the freezing behavior of confined salt solutions.
Lanthanide ion solvation chemistry in nonaqueous phases is key to understanding and developing effective separation processes for these critical materials. Due to the complexity and inherent disorder of the solution phase, a comprehensive picture of the solvated metal ion is often difficult to generate solely from conventional spectroscopic approaches and electronic structure calculations, particularly in the extractant phase. In this work, we use classical molecular dynamics (MD) simulation with an advanced sampling technique, metadynamics, supplemented by experimental spectroscopy and speciation analysis, to measure lanthanide solvation free energy landscapes. We define coordination-based collective variables to probe the entire range of solvation configurations in the organic phase of lanthanum (La), europium (Eu), and lutetium (Lu) nitrate salts bound with a commonly used extractant, N,N′-dimethyl, N,N′-dioctylhexylethoxymalonamide (DMDOHEMA). The known lanthanide extraction trend of La ≈ Eu > Lu is readily explained by the measured free energy surfaces, which show consistent DMDOHEMA coordination from La to Eu, followed by loss of DMDOHEMA coordination from Eu to Lu. These simulations suggest how ligand crowding at the metal center can control selectivity, in this case resulting in the opposite extraction trend as observed with other conventional extractants, where the enthalpic contribution from increasing lanthanide charge density across the series dominates the extraction energetics. We also find that the presence of inner-sphere water, verified by time-resolved fluorescence, diversifies the accessible solvation structures. As a result, understanding solvation requires consideration of an entire thermodynamic ensemble, rather than the single dominant lowest-energy structure, as is often considered out of necessity in interpretation of spectroscopic data or in electronic structure-based ligand design approaches. In general, we demonstrate how metadynamics uniquely enables investigation of complex, multidimensional solvation energetic landscapes, and how it can explain selectivity trends where extraction is controlled by more complex mechanisms than simple charge density-based selectivity.
Polydimethylsiloxane is an important polymeric material with a wide range of applications. However, environmental effects like low temperature can induce crystallization in this material with resulting changes in its structural and dynamic properties. The incorporation of phenyl-siloxane components, e.g., as in a poly(dimethyl-co-diphenyl)siloxane random copolymer, is known to suppress such crystallization. Molecular dynamics (MD) simulations can be a powerful tool to understand such effects in atomistic detail. Unfortunately, all-atomistic molecular dynamics (AAMD) is limited in both spatial dimensions and simulation times it can probe. Here, to overcome such constraints and to extend to more useful length- and time-scales, we systematically develop a coarse-grained molecular dynamics (CGMD) model for the poly(dimethyl-co-diphenyl)siloxane system with bonded and non-bonded interactions determined from all-atomistic simulations by the iterative Boltzmann inversion (IBI) method. Additionally, we propose a lever rule that can be useful to generate non-bonded potentials for such systems without reference to the all-atomistic ground truth. Our model captures the structural and dynamic properties of the copolymer material with quantitative accuracy and is useful to study long-time dynamics of highly-entangled systems, sequence-dependent properties, phase behaviour, etc.
Nature routinely exploits liquid–liquid phase separation (LLPS) of proteins to control the assembly and mineralization of hybrid materials. Here, we show that fusion of the Car9 silica-binding peptide to an elastin-like polypeptide (ELP) yields temperature- and sequence-programmable soft matter templates for the synthesis of silicified architectures ranging in size from nanometers to micrometers. Specifically, we demonstrate unprecedented control over the diameter of silica nanoparticles (SiNP) in the 30–60 nm range with 4 nm precision, show that a single arginine residue (R4) in the Car9 sequence underpins the transition from micelles to proteinosomes, and find that substitutions in other basic residues modulate electrostatic repulsion and solvation to enable access to kinetically trapped species. These structures, which include interconnected micelles, small (∼200 nm) and large (>5 µm) vesicles, are readily visualized by SEM imaging following silicification. Molecular dynamics (MD) simulations and AlphaFold predictions reveal that mutations in positively charged residues alter interfacial packing, hydration, and conformational freedom of the silica-binding segments. Overall, our results establish sequence and thermal energy as synergistic levers for morphological control across length scales using solid-binding ELPs and establish mineralization as a powerful tool to visualize the structure of dynamic soft matter assemblies.
We report a drastic difference in stacking behavior of oleic acid-stabilized 4-monolayer (4 ML) CdSe nanoplatelets (NPLs) in toluene and methylcyclohexane (MCH), two nonpolar solvents that differ in the conformational flexibility of their molecules. Using liquid cell transmission electron microscopy (TEM) and small angle scattering (SAXS) techniques, we show that NPLs form microns-long ribbons consisting of 4 ML CdSe NPLs in toluene, the solvent widely used to form stable colloidal solutions of a broad range of quasi-spherical nanoparticles. In contrast, 4 ML CdSe NPLs are well dispersed in MCH. The difference in stacking behavior of NPLs in toluene and MCH suggests that the conformational flexibility of the solvent molecules, such as the ability to adopt multiple chair conformations, modulates nanoplatelet interactions. Molecular dynamics (MD) simulations reveal that solvent molecules subtly alter the structure of the organic ligand shell. These solvent-dependent changes propagate to the inorganic core, modulating the degree of CdSe nanoplatelet (NPL) twisting and, consequently, the properties of the nanoparticles. We show that toluene better solvates oleate ligands while MCH induces a bimodal oleate span distribution, which can lead to increased solubility of CdSe NPLs. In addition, the solvent can also influence the inorganic core, which, in turn, can modify the nanoparticle properties. We demonstrate that destabilization of toluene solution containing ribbons of 4 ML CdSe NPLs without CdS shells results in the formation of NPL assemblies with amplified spontaneous emission (ASE) with a low threshold of 14 µJ cm−2 that is comparable with that of CdSe/CdS core/shell NPLs. Our results emphasize that the solvent plays a major role in mediating interactions between NPLs and hence their processability for fabrication of functional structures.
Thermal neutron scattering cross-section data for light water available in the major nuclear data libraries observes significant differences especially at reactor operating temperatures. During the past few years there has been a renewed interest in reviewing the existing thermal scattering models and generating more accurate and reliable thermal scattering cross sections using existing experimental data and in some cases based on Molecular Dynamics (MD) simulations. There is a need for performing new time-of-flight experiments at high temperatures and pressures, to have a better understanding of the physics involved in the scattering process that could help improve the existing TSL data.Lack of experimental thermal scattering data for light water at high temperatures led to a new measurement campaign within the INSIDER project at the Institut de radioprotection et de sûreté nucléaire (IRSN). Double differential scattering cross section for light water have been measured at the SEQUOIA spectrometer based at the Spallation Neutron Source (SNS), Oak Ridge National Laboratory, United States. Several measurements have been carried out at different temperatures and pressures corresponding to liquid light water. Measurements at five different incident neutron energies Ei (8, 60, 160, 280 and 800 meV) have been carried out to help exploring different regions of the frequency spectrum. This paper presents the analysis of the dynamic structure factor and the derived frequency spectrum of light water. The analysis of the experimental data would provide one with better confidence, the behavior of thermal scattering cross sections for light water at high temperatures, knowledge of which is very important for the design of novel reactors as well as existing pressurized water reactors.
Here, the effect of concentration of the alloying element silicon (Si) on the shock and spall response of aluminum (Al) is studied using classical molecular dynamics (MD) simulations. Silicon is distributed as individual grains and as grain boundaries and the dynamic response under shock is studied for an impact velocity of 1000 m/s. A general trend observed here is that the spallation is almost exclusively intergranular and limited to Al-Si interfaces. The spall strengths are observed to a) decrease with an increase in alloying of Si as grains and b) show little increase for up to 50% grain boundary segregation. An increase in grain size of the whole system showed negligible differences in spall strength, highlighting the significance of concentration of the alloying species rather than its grain size. These results provide an atomistic understanding of light-weight metallic systems under shock compression and pave the way for designing multiphase metal matrix alloys and composites for defense/armor applications.
Molecular dynamics (MD) simulations are performed to investigate the failure behavior of Al/Al, Al/Fe, and Al/AlFe bi-layers under high strain-rate shock loading. Additional simulations are performed to invesigate solely the tensile response of these systems and to isolate the effects of loading history on damage and failure. Specifically, the compression stage of the shock is excluded from the tensile simulations. This study shows that local structure variation, including the introduction of serrations, plays different roles in controlling damage nucleation depending on the specific loading condition. Under shock loading, the stress for void nucleation and the resulting void distribution is insensitive to an interface structure, whereas under pure tensile loading, the opposite is true, even though the strain rate is comparable. The above difference can be explained based on the corresponding change in the total deformation prior to damage nucleation under shock loading, which is missing from the pure tensile loading. Therefore, it is concluded that whether local interface structure matters under high strain-rate loading depends on the specific loading history and the stress state: the local structure variation matters only when it alters the deformation behavior, and subsequently, damage nucleation and evolution.
Controlling electrochemical reactivity requires a detailed understanding of the charging behavior and thermodynamics of the electrochemical interface. Experiments can independently probe the overall charge response of the electrochemical double layer by capacitance measurements and the thermodynamics of the inner layer with potential of maximum entropy measurements. Relating these properties by computational modeling of the electrochemical interface has so far been challenging due to the low accuracy of classical molecular dynamics (MD) for capacitance and the limited time and length scales of ab initio MD. Here, we combine large ensembles of long-time-scale classical MD simulations with charge response from electronic density functional theory to predict the potential-dependent capacitance of a family of ideal aqueous electrochemical interfaces with different peak capacitances. We show that while the potential of maximum capacitance varies, this entire family exhibits an electrode charge of maximum capacitance (CMC) between -2.9 and -2.2 μC/cm 2 , regardless of the details in the electronic response. Simulated heating of the same interfaces reveals that the entropy peaks at a charge of maximum entropy (CME) of -5.1 ± 0.6 μC/cm 2 , in agreement with experimental findings for metallic electrodes. The CME and CMC both indicate asymmetric response of interfacial water that is stronger for negatively charged electrodes, while the difference between CME and CMC illustrates the richness in behavior of even the ideal electrochemical interface.
Supported metallic nanoparticles play a central role in catalysis. However, predictive modeling is particularly challenging due to the structural and dynamic complexity of the nanoparticle and its interface with the support, given that the sizes of interest are often well beyond those accessible via traditional ab initio methods. With recent advances in machine learning, it is now feasible to perform MD simulations with potentials retaining near-density-functional theory (DFT) accuracy, which can elucidate the growth and relaxation of supported metal nanoparticles, as well as reactions on those catalysts, at temperatures and time scales approaching those relevant to experiments. Furthermore, the surfaces of the support materials can also be modeled realistically through simulated annealing to include effects such as defects and amorphous structures. We study the adsorption of fluorine atoms on ceria and silica supported palladium nanoparticles using machine learning potential trained by DFT data using the DeePMD framework. We show defects on ceria and Pd/ceria interfaces are crucial for the initial adsorption of fluorine, while the interplay between Pd and ceria and the reverse oxygen migration from ceria to Pd control spillover of fluorine from Pd to ceria at later stages. In contrast, silica supports do not induce fluorine spillover from Pd particles.