Search NASA⌕ Search

SEARCH · Search NASA

Results for “molecular simulation”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Computational Nanomechanics of Carbon Nanotubes and Composites

Nanomechanics of individual carbon and boron-nitride nanotubes and their application as reinforcing fibers in polymer composites has been reviewed with interplay of theoretical modeling, computer simulations and experimental observations. The emphasis in this work is on elucidating the multi-length scales of the problems involved, and of different simulation techniques that are needed to address specific characteristics of individual nanotubes and nanotube polymer-matrix interfaces. Classical molecular dynamics simulations are shown to be sufficient to describe the generic behavior such as strength and stiffness modulus but are inadequate to describe elastic limit and nature of plastic buckling at large strength. Quantum molecular dynamics simulations are shown to bring out explicit atomic nature dependent behavior of these nanoscale materials objects that are not accessible either via continuum mechanics based descriptions or through classical molecular dynamics based simulations. As examples, we discus local plastic collapse of carbon nanotubes under axial compression and anisotropic plastic buckling of boron-nitride nanotubes. Dependence of the yield strain on the strain rate is addressed through temperature dependent simulations, a transition-state-theory based model of the strain as a function of strain rate and simulation temperature is presented, and in all cases extensive comparisons are made with experimental observations. Mechanical properties of nanotube-polymer composite materials are simulated with diverse nanotube-polymer interface structures (with van der Waals interaction). The atomistic mechanisms of the interface toughening for optimal load transfer through recycling, high-thermal expansion and diffusion coefficient composite formation above glass transition temperature, and enhancement of Young's modulus on addition of nanotubes to polymer are discussed and compared with experimental observations.

Srivastava, Deepak↗
NASA NTRS2001

Chemical and Flowfield Modeling for Enhanced Analysis of Contamination Experiments

This paper describes the application of a new Direct Simulation Monte Carlo (DSMC) code, the Molecular Beam Simulator (MBS), which is designed to analyze laboratory scale molecular beam-surface (and crossed-beam) experiments. The MBS is primarily intended to model experiments associated with spacecraft contamination effects, but it can also be used to simulate a variety of surface chemistry and reactive flow measurements. The MBS code is fully three-dimensional, includes a wide-range of chemical processes, and can model one or multiple pulsed (non-steady) sources. As an example application of the MBS code, a fast, pulsed, oxygen atom-surface experiment which examines the chemistry behind erosion of graphite by oxygen atoms is analyzed. Unsteady DSMC simulations show that experimental observations of excited molecular states after the pulse has hit the surface are consistent with two distinct chemical mechanisms: a direct one where the excited molecules are formed on the surface, and a two-step mechanism where ground state molecules formed on the surface are collisionally excited after they leave the surface by trailing oxygen atoms in the pulse. Further DSMC calculations suggest experiments which can distinguish between these mechanisms.

Braunstein, Matthew↗
DOE OSTI2025

Developing reliable machine learning interatomic potential for Fe–Cr–Ni austenitic alloys

Gaining atomistic understanding of mechanical behavior of heat-resistant structural materials such as Fe–Cr–Ni-based alloys requires an approach with an accuracy close to density functional theory (DFT) that considers the intrinsic properties of the bulk lattice and important defects such as stacking faults, grain boundaries, and surfaces. This work aims to develop reliable machine learning interatomic potential (MLIAP) at cross-scale for Fe–Cr–Ni ternary alloys with a focus on the face-centered-cubic (fcc) solid solution structure. Leveraging the advantages of moment tensor potentials, which typically necessitate a relatively small training dataset and enable rapid calculations using the large-scale atomic/molecular massively parallel simulator package, we ensure the stability and accuracy of the trained potentials. Important defects such as stacking faults, grain boundaries, and surfaces for wide-range compositions are investigated. Structural, thermal, elastic, and defect properties are determined from molecular dynamics simulations comprising several thousand atoms, generated via canonical Monte Carlo simulations guided by the trained potential. The trained potential allows efficient atomic simulations of structural, thermal, and mechanical properties of fcc Fe–Cr–Ni solid solution alloys as a function of composition and temperature. Therefore, the MLIAP approach represents a major advancement from DFT calculations that are limited to small simulation sizes and traditional molecular dynamics simulations using relatively low accuracy potentials. Furthermore, this work outlines a practical foundation for further investigating the structural evolution and mechanical behavior of austenitic stainless steel and nickel-based alloys in a wide array of applications in extreme environments.

Crystal structure↗
DOE OSTI2026

Effect of Solvent on the Local Structure, Dynamics, and Vibrational Density of States in Sn-BEA Zeolite

Lewis acid zeolites are attractive catalysts for epoxidation and biomass valorization, as they are highly active and selective in the liquid phase and can operate at or near ambient conditions. While a rich experimental literature exists on liquid-phase Lewis acid zeolite catalysis, our understanding of the molecular organization and solvent dynamics in the vicinity of Lewis acid sites with differing metal site speciation remains limited. In this work, we investigate the molecular coordination and diffusion of two common solvents (methanol and water) around the closed and open Sn-BEA zeolite active sites using molecular dynamics simulations with a machine-learned interatomic potential trained on ab initio molecular dynamics trajectories. Molecular dynamics simulations reveal that introducing active sites significantly enhances local order in the first and second solvation shells compared to the pure silica case. For methanol, both closed and open active sites are singly coordinated, while more than two water molecules coordinate the open site. In contrast to methanol, we observed that water molecules dissociate, leading to the formation of additional Sn-OH and silanol groups away from the active site. The diffusion coefficients of water and methanol are functions of the solvent population in the pore. Here, our work provides insights into how active site speciation in Lewis acid zeolites affects solvent coordination, diffusion, and vibrational signature. This information is foundational for catalyst design and optimization of liquid-phase catalytic processes in zeolites. It also demonstrates the suitability of machine-learned interatomic potentials for modeling reactive systems, enabling sufficiently long trajectories for appropriate statistical averaging.

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

Cross-Linked Nanotube Materials with Variable Stiffness Tethers

The constitutive properties of a cross-linked single-walled carbon nanotube material are predicted with a multi-scale model. The material is modeled as a transversely isotropic solid using concepts from equivalent-continuum modeling. The elastic constants are determined using molecular dynamics simulation. Some parameters of the molecular force field are determined specifically for the cross-linker from ab initio calculations. A demonstration of how the cross-linked nanotubes may affect the properties of a nanotube/polyimide composite is included using a micromechanical analysis.

Frankland, Sarah-Jane V.↗
DOE OSTI2025

Shock-induced chemistry and high strain-rate viscoelastic behavior of a phenolic polymer

We use impact experiments and a finite element model (up to 1.2 GPa), and molecular dynamics simulations (up to 60 GPa), to examine the behavior of a phenolic polymer under shock compression, spanning both nonreactive and reactive regimes. In the nonreactive regime, relaxation following compression at strain rates of ∼105 s−1 can be explained by viscoelasticity observed at ordinary laboratory rates (≲1 s−1) by accounting for the temperature dependence of the phenolic β-transition. Reasonable agreement is found between the measured shock Hugoniot up to 1.2 GPa and molecular dynamics simulation for cross-linked structures of comparable density. We also observed a first-order mechanical transition near 0.36 GPa shock stress and estimated a spall strength of 0.102 GPa and Hugoniot elastic limit of 1–2 GPa. The shock stress is found to vary up to 24% among phenolics made with different resin and/or cure processes. Finally, molecular dynamics simulations are used to identify a reactive regime at shock pressures ≳20 GPa that is characterized by chemically driven, rate-dependent relaxation processes, including dehydrogenation and dehydration reactions that promote the formation of a dense, highly cross-linked carbonaceous solid and the release of light volatiles.

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

Solvent Isotherms and Structural Transitions in Nanoparticle Superlattice Assembly

Here, we introduce a Molecular Theory for Compressible Fluids (MOLT-CF) that enables us to compute free energies and other thermodynamic functions for nanoparticle superlattices with any solvent content, including the dry limit. Quantitative agreement is observed between MOLT-CF and united-atom molecular dynamics simulations performed to assess the reliability and precision of the theory. Among other predictions, MOLT-CF shows that the amount of solvent within the superlattice decreases approximately linearly with its vapor pressure and that in the late stages of drying, solvent-filled voids form at lattice interstitials. Applied to single-component superlattices, MOLT-CF predicts fcc-to-bcc Bain transitions for decreasing vapor pressure and for increasing ligand length, both in agreement with experimental results. We explore the stability of other single-component phases and show that the C14 Frank-Kasper phase, which has been reported in experiments, is not a global free-energy minimum. Implications for precise assembly and prediction of multicomponent nanoparticle systems are discussed.

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

Fouling behavior of zwitterionic membranes compared to polyamide membranes

Membrane fouling remains a critical bottleneck for reverse osmosis (RO) desalination, driving energy consumption and reducing membrane lifetime. Here, we employ all-atom molecular dynamics simulations to investigate the antifouling behavior of random zwitterionic amphiphilic copolymer (r-ZAC) membranes composed of sulfobetaine methacrylate (SBMA) and allyl methacrylate (AMA), benchmarked against conventional polyamide (PA) RO membranes. Structural and dynamical analyses—including radial distribution functions, coordination numbers, tetrahedral order parameters, vector orientation, and residence-time correlation functions—reveal that r-ZAC surfaces sustain tightly bound, long-lived hydration layers with preserved tetrahedrality and anisotropic water orientation, in sharp contrast to the weak and disordered hydration of PA. Steered molecular dynamics simulations demonstrate that r-ZAC membranes impose substantial free-energy barriers to foulant approach (alginate ≈ 90 kcal/mol, sucrose ≈ 35 kcal/mol, humic acid ≈ 15 kcal/mol), whereas PA membranes exhibit negligible barriers (< 1 kcal/mol) and thermodynamically favorable adsorption. Detailed foulant–surface interaction analyses show that zwitterionic hydration and electrostatic heterogeneity in r-ZAC suppress adhesion, except in the case of amphiphilic humic acid, which exploits multiple binding modes. Together, these results establish molecular-level design principles for antifouling membranes: the combination of zwitterionic hydration, structured interfacial water, and controlled amphiphilic balance in r-ZAC membranes provides superior resistance to organic fouling relative to PA.

Cross-linked polyamide↗
DOE OSTI2025

Accurate and efficient parameterization of an atomic cluster expansion (ACE) potential for ammonia under extreme conditions

We present a machine learning interatomic potential for ammonia designed to capture its complex multiphase behavior, including both molecular and superionic phases. The potential is based on the atomic cluster expansion (ACE) formulation and has been parameterized to facilitate high-fidelity molecular dynamics simulations of ammonia under extreme conditions, for pressures up to 100 GPa and for temperatures above 500 K and up to 6000 K. A diverse range of configurations was generated through high-quality ab initio molecular dynamics simulations, covering insulating and superionic ice phases, liquid ammonia, molecular nitrogen (N 2 ) and hydrogen (H 2 ), and metastable compounds that form upon dissociation, including $NH^{+}_{4}$, $H^{+}_{3}$, N 2 H 4 , and N 3 H. We demonstrate that the ammonia ACE potential accurately reproduces experimental and density functional theory predicted isotherms and Hugoniots. Crucially, the potential is able to capture the intricate phase behavior of ammonia, including the transition from insulating molecular fluid to the superionic phase. This work provides a robust interatomic potential that can be used for large-scale, accurate simulations of ammonia under extreme thermodynamic conditions, offering a powerful tool for investigating its behavior in various phases and applications.

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

Ranges of applicability for the continuum beam model in the mechanics of carbon nanotubes and nanorods

Limitations in the validity of the continuum beam model for carbon nanotubes (NTs) and nanorods are examined. Applicability of all assumptions used in the model is restricted by the two criteria for geometric parameters that characterize the structure of NTs. The key non-dimensional parameters that control the NT buckling behavior are derived via dimensional analysis of the nanomechanical problem. A mechanical law of geometric similitude for NT buckling is extended from continuum mechanics for different molecular structures. A model applicability map, where two classes of beam-like NTs are identified, is constructed for distinct ranges of non-dimensional parameters. Expressions for the critical buckling loads and strains are tailored for two classes of NTs and compared with the data provided by the molecular dynamics simulations. copyright 2001 Elsevier Science Ltd. All rights reserved.

Continuum beam model↗
DOE OSTI2024

Melting point of iron at high pressure: An assessment of uncertainties and effect of electronic temperature

An accurate calculation of the melting point of iron at various pressures in the Earth's core is important for understanding the core structure, geodynamo, and the Earth's history. Previous studies have assessed the melt line of iron at these extreme conditions using various experimental measurement techniques as well as both ab initio and classic molecular dynamics simulations. However, experimental measurements have uncertainties up to several hundred Kelvin, and inconsistencies remain among simulation results. Here in this work, we propose an iterative framework that couples density functional theory (DFT) calculations and molecular dynamics simulations performed using an ensemble of interatomic potentials to assess the effect of electronic temperature on the melting point. We systematically validate the potentials by comparing lattice constants and phonon dispersion curves at 0 K and enthalpy differences between liquid and HCP, FCC, BCC phases of iron close to the melt line at 300 GPa with DFT. Our results show that HCP iron melts at 6144 K (at 300 GPa), BCC phase is thermodynamically unstable, and FCC is metastable at this temperature. The melting points of FCC and BCC phases at 300 GPa are 5858 and 5647 K, respectively.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1999

Bone Blood Flow During Simulated Microgravity: Physiological and Molecular Mechanisms

Blood flow to bone has been shown to affect bone mass and presumably bone strength. Preliminary data indicate that blood flow to the rat femur decreases after 14 days of simulated microgravity, using hindlimb suspension (HLS). If adult rats subjected to HLS are given dobutamine, a synthetic catecholamine which can cause peripheral vasodilation and increased blood flow, the loss of cortical bone area usually observed is prevented. Further, mechanisms exist at the molecular level to link changes in bone blood flow to changes in bone cell activity, particularly for vasoactive agents like nitric oxide (NO). The decreases in fluid shear stress created by fluid flow associated with the shifts of plasma volume during microgravity may result in alterations in expression of vasoactive agents such as NO, producing important functional effects on bone cells. The primary aim of this project is to characterize changes in 1) bone blood flow, 2) indices of bone mass, geometry, and strength, and 3) changes in gene expression for modulators of nitric oxide activity (e.g., nitric oxide synthase) and other candidate genes involved in signal transduction of mechanical loading after 3, 7, 14, 21, and 28 days of HLS in the adult rat. Using a rat of at least 5 months of age avoids inadvertently studying effects of simulated microgravity on growing, rather than adult, bone. Utilizing the results of these studies, we will then define how altered blood flow contributes to changes in bone with simulated microgravity by administering a vasodilatory agent (which increases blood flow to tissues) during hindlimb suspension. In all studies, responses in the unloaded hindlimb bones (tibial shaft, femoral neck) will be compared with those in the weightbearing humeral shaft and the non-weightbearing calvarium (skull) from the same animal. Bone volumetric mineral density and geometry will be quantified by peripheral quantitative CT; structural and material properties of the long bones will be determined by 3-point bending (tibia, humerus) or compression (femoral neck) testing to failure. A unique aspect of these studies will be defining the time course of changes in gene expression in bone cell populations with unloading, accomplished with Northern blots, in situ hybridization, and immunohistochemistry. These studies have high relevance for concurrent protocols being proposed by investigators on NSBRI Cardiovascular and Muscle teams, with blood flow data available on a number of tissues other than bone. Further, dobutamine and other Beta-agonists have been tested as countermeasures for altered muscle and cardiovascular function. Results of the intervention tested in our studies have potential relevance for a number of systemic changes seen with prolonged spaceflight.

Bloomfield, Susan A.↗
DOE OSTI2026

Custom-trained Machine-learning Interatomic Potentials: ZnCl2 Aqueous Solution

This dataset was generated using an iterative active-learning strategy implemented in the ArcaNN software package (https://github.com/arcann-chem/arcann_training) to train machine-learning interatomic potentials for aqueous ZnCl2 solutions. Each active-learning cycle consisted of three stages: training, exploration, and labeling. The initial training set combined configurations generated in this work from enhanced-sampling ab initio molecular dynamics simulations with configurations from a previously reported neural-network-potential study of aqueous ZnCl2. The enhanced-sampling ab initio molecular dynamics simulations involved Zn–Cl separation and the chloride coordination number around Zn²? as collective variables. These configurations served as the seed dataset. Subsequent active-learning cycles expanded the training set by identifying and labeling configurations that were poorly represented by the current models, thereby improving coverage of ion-association states and changes in local coordination and charge-state environments relevant to the solution free-energy landscape. For all selected configurations, single-point calculations of the total energies and atomic forces were performed within density functional theory using the CP2K Quickstep module. Reference calculations employed the revPBE-D3 and r2SCAN exchange-correlation functionals. Motivated by recent work on aqueous Zn²?, the main revPBE calculations omitted D3 dispersion contributions involving Zn²?, while retaining the D3 correction for water and chloride. For comparison, fully dispersion-corrected revPBE-D3 reference calculations were also performed, with D3 applied to all species, including Zn²?. Valence electrons were treated explicitly, while core electrons were represented using norm-conserving Goedecker–Teter–Hutter pseudopotentials. The wave functions were expanded using the mixed Gaussian-and-plane-wave scheme with TZV2P-MOLOPT basis sets for all elements and a 600 Ry auxiliary plane-wave cutoff for the electron density. Self-consistent-field convergence was accelerated using the orbital-transformation and Direct Inversion in the Iterative Subspace algorithms, with a convergence threshold of 10?6. All single-point calculations were performed in periodic orthorhombic cells. The CELL_REF keyword in CP2K was used to define a fixed reference cell with a box length of 25 Å. This treatment ensured a consistent reference for configurations extracted from NpT trajectories with fluctuating cell dimensions. The resulting DFT energies and atomic forces constitute the ground-truth labels used to train the MLIPs. The resulting MLIP was trained for aqueous ZnCl2 solutions spanning concentrations from 0 to 30 molal and a broad pH range, from strongly acidic to strongly basic conditions. Representative examples of configurations included in the MLIP training dataset are provided below. These include 1) Representative configurations from the dataset labeled at the revPBE-D3 level, with D3 dispersion interactions involving Zn2+ excluded (revPBE-wo-D3). 2) Representative configurations from the dataset labeled at the fully dispersion-corrected revPBE-D3 level, with D3 interactions applied to all species, including Zn2+ (revPBE-D3). 3) Representative configurations from the dataset labeled at the r2SCAN level of theory (r2SCAN).

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗
DOE OSTI2024

Beyond nothingness in the formation and functional relevance of voids in polymer films

Abstract Voids—the nothingness—broadly exist within nanomaterials and impact properties ranging from catalysis to mechanical response. However, understanding nanovoids is challenging due to lack of imaging methods with the needed penetration depth and spatial resolution. Here, we integrate electron tomography, morphometry, graph theory and coarse-grained molecular dynamics simulation to study the formation of interconnected nanovoids in polymer films and their impacts on permeance and nanomechanical behaviour. Using polyamide membranes for molecular separation as a representative system, three-dimensional electron tomography at nanometre resolution reveals nanovoid formation from coalescence of oligomers, supported by coarse-grained molecular dynamics simulations. Void analysis provides otherwise inaccessible inputs for accurate fittings of methanol permeance for polyamide membranes. Three-dimensional structural graphs accounting for the tortuous nanovoids within, measure higher apparent moduli with polyamide membranes of higher graph rigidity. Our study elucidates the significance of nanovoids beyond the nothingness, impacting the synthesis‒morphology‒function relationships of complex nanomaterials.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Ionic Liquid-Enhanced Interfaces to Boost Reactive C O2 Capture

The addition of ionic liquids (ILs) to a mixture containing a molecular solvent and other ionic species can induce the heterogeneous redistribution of cations and anions at the gas–liquid interface. This nonuniform redistribution of cations and anions driven by the differences in the solvophilicity of ions can improve the thermophysical and interfacial properties of such mixtures, creating a local chemical environment that is conducive to some reactions. In this work, ILs are added to a mixture of potassium hydroxide (KOH) and ethylene glycol (EG), used as a reactive absorbent and electrolyte in the migration-assisted moisture-gradient (MAMG) process for CO 2 capture. Molecular dynamics (MD) simulations are employed to probe into the effects of complex ion–ion and ion–solvent interactions and to examine the chemical composition at the gas–liquid interface. A total of 12 systems are investigated using molecular simulations to identify trends in the performance of IL additives based on the choice of cation, anion, and IL concentration. The cation effects are studied using IL additives based on 1-ethyl-3-methylimidazolium ([EMIM] + ) and 1-butyl-3-methylimidazolium ([BMIM] + ), while the impact of anions is examined using additives based on dicyanamide [DCA] − , triflate [TfO] − , bistriflimide [NTf 2 ] − , and hexafluorophosphate [PF 6 ] − anions, respectively. The influence of the IL concentration is also evaluated at molar concentrations between 1% and 4%. The simulation results indicate that the use of IL additives can affect the physical CO 2 solubility, surface tension, and the localization of CO 2 around the [OH] − ions at the gas–liquid interface. It is also evident that the choice of cations, anions, and IL concentration determines the extent to which the IL additives impact the local physicochemical properties. Physical dissolution, diffusive transport, and interaction with [OH] − are critical intermediate steps toward reactive CO 2 capture using a liquid absorbent. Hence, the improvement in one or more of these properties, aided by IL additives, is expected to improve the overall CO 2 capture performance. Experiments reaffirmed the impact of IL additives on CO 2 capture performance and the sensitivity to the choice of the cation, anion, and concentration of the IL additive.

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

Release dynamics of nanodiamonds created by laser-driven shock-compression of polyethylene terephthalate

Laser-driven dynamic compression experiments of plastic materials have found surprisingly fast formation of nanodiamonds (ND) via X-ray probing. This mechanism is relevant for planetary models, but could also open efficient synthesis routes for tailored NDs. We investigate the release mechanics of compressed NDs by molecular dynamics simulation of the isotropic expansion of finite size diamond from different P-T states. Analysing the structural integrity along different release paths via molecular dynamic simulations, we found substantial disintegration rates upon shock release, increasing with the on-Hugnoiot shock temperature. We also find that recrystallization can occur after the expansion and hence during the release, depending on subsequent cooling mechanisms. Our study suggests higher ND recovery rates from off-Hugoniot states, e.g., via double-shocks, due to faster cooling. Laser-driven shock compression experiments of polyethylene terephthalate (PET) samples with in situ X-ray probing at the simulated conditions found diamond signal that persists up to 11 ns after breakout. In the diffraction pattern, we observed peak shifts, which we attribute to thermal expansion of the NDs and thus a total release of pressure, which indicates the stability of the released NDs.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Structure and Dynamics of CO 2 at the Air–Water Interface from Classical and Neural Network Potentials

The accurate description of the structure and dynamics of CO 2 at the instantaneous air–water interface, along with the effects of surface fluctuations on the CO 2 -transport processes, is essential for the development of negative emission technologies aimed at minimizing climate change. In this study, we performed molecular dynamics simulations of CO 2 at the air–water interface using neural network potentials (NNPs) trained on ab initio data generated through density-functional-theory-based molecular dynamics simulations. We compared these results with classical force fields to assess their performance in modeling interfacial CO 2 behavior. Our findings revealed that the asymmetric interactions, coupled with thermal surface fluctuations at the air–water interface, significantly influence CO 2 transport into the aqueous phase. The simulations demonstrate that classical force fields underestimate both the free energy of CO 2 transport and the strength of its interactions at the interface compared with the neural network potentials. In conclusion, the free energy and the interfacial dynamics of CO 2 are primarily influenced by the distribution of water within the instantaneous interfacial water layer, responsible for creating an asymmetric intermolecular interaction environment within the interfacial region.

Ab initio molecular dynamics↗
DOE OSTI2025

Long timescale solvation dynamics and confinement: The case of non-ionic deep eutectic solvents of lauric acid and N-methylacetamide

Microscopic segregation and molecular heterogeneities in complex liquids are the result of the interplay between different intermolecular forces, all of which contribute to the energy landscape of the system. A consequence of the intricate energy landscape is the nontrivial effect on the solvation dynamics. Here, the impact of molecular heterogeneities on the solvation dynamics is studied using infrared spectroscopies and molecular dynamics simulations. In particular, this study focuses on the dynamical effect of nanoscopic heterogeneities present in deep eutectic solvents (DESs) composed of lauric acid (LA) and N-methylacetamide (NMA). To this end, a molecular probe containing a carbon triple bond is used as an infrared reporter. The results show that the vibrational probe is likely to be located in the NMA polar domains. Furthermore, the probe solvation dynamics derived from the 2DIR spectra presents a slowdown of its timescale with increasing LA concentration in the DES. Kubo modeling of the probe solvation dynamics shows a correlation between the amplitude of its long time component and the presence of molecular heterogeneities in the sample. Semiclassical modeling of the vibrational line shape of the triple bond stretch demonstrates that the heterogeneities affect the whole solvation dynamics of the system through the amplitudes of the frequency fluctuations. Molecular dynamics simulations confirm the experimental results and their interpretation by showing a slowdown of the solvation dynamics when the LA heterogeneities are present. Overall, the study presents a molecular framework to explain the effect of confinement created by nanoscopic LA heterogeneities on the solvation dynamics of the system.

Chemical compounds↗