Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Solvation”

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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Solutes at the liquid:liquid phase boundary—Solubility and solvent conformational response alter interfacial microsolvation

A detailed understanding of solvent structure and dynamics at liquid:liquid interfaces is a necessary precursor for control and manipulation of these phase boundaries. Experimentally, amphiphilic solutes are often used to alter transport properties across water:organic interfaces; yet, a fundamental model for the mechanism of this action has not been determined. Our study compares the solvation profiles of ampiphilic solutes that traverse the phase boundary in binary water:n-hexane, and the individual microsolvation processes for interfacial water and hexane molecules therein. Microsolvation is defined as the rare event where one solvent molecule temporarily penetrates the co-solvent phases and is fully solvated therein. The solutes tri-butyl phosphate (TBP), hydrogen di-butyl phosphate, and di-hydrogen mono-butyl phosphate have been examined as they exhibit a systematic increase in aqueous solubility and selectively partition to the interfacial region at the infinite dilution limit. The relationship between adopted configurations of the solute, orientation of the solvent, and the ability of the solute to enhance microsolvation, specifically the ability of n-hexane to penetrate the aqueous phase, is observed within a 20 Å radius of TBP.

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

Composition dependent structural organization in trihexyl(tetradecyl)phosphonium chloride ionic liquid-methanol mixtures

This article reports results from the molecular dynamics simulations on the structural arrangement of the ions and molecules in the mixtures of trihexyl(tetradecyl)phosphonium chloride ([P{sub 666,14}{sup +}][Cl{sup −}]) ionic liquid (IL) and methanol (MeOH) over the entire composition range. Effects of composition on the charge and polarity orderings have been investigated via computation of X-ray scattering structure function, S(q), and by using a partitioning scheme proposed for such multi-component mixtures. Except for the neat methanol liquid, the total S(q) shows two peaks in its intermolecular region for all the mole-fractions. The lowest q peak is dominated primarily by anion-anion, cation-anion, and methanol-anion correlations. Our results signify that the methanol bulk structure, which predominantly has short-distance characteristic correlations and is governed by polar group of methanol, is retained for x{sub IL} ≤ 0.1. Then, the mixture goes through gradual structural changes from methanol-like to the IL-like for 0.1 < x{sub IL} ≤ 0.7. The dipolar interaction between methanol molecules weakens in this range, and the structural landscape of the mixture is steered by strong ion-ion, anion-methanol, and nonpolar interactions. The IL-like structural arrangement is virtually recovered for x{sub IL} > 0.7. At all the compositions studied, while the cation head groups are predominantly solvated by anions and subsequently by methanol molecules, the polar hydroxyl group of methanol is preferentially solvated by the anions. The radial distribution functions of selected pair of atomic species have also confirmed these observations.

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

Water around fullerene shape amphiphiles: A molecular dynamics simulation study of hydrophobic hydration

Fullerene C{sub 60} sub-colloidal particle with diameter ∼1 nm represents a boundary case between small and large hydrophobic solutes on the length scale of hydrophobic hydration. In the present paper, a molecular dynamics simulation is performed to investigate this complex phenomenon for bare C{sub 60} fullerene and its amphiphilic/charged derivatives, so called shape amphiphiles. Since most of the unique properties of water originate from the pattern of hydrogen bond network and its dynamics, spatial, and orientational aspects of water in solvation shells around the solute surface having hydrophilic and hydrophobic regions are analyzed. Dynamical properties such as translational-rotational mobility, reorientational correlation and occupation time correlation functions of water molecules, and diffusion coefficients are also calculated. Slower dynamics of solvent molecules—water retardation—in the vicinity of the solutes is observed. Both the topological properties of hydrogen bond pattern and the “dangling” –OH groups that represent surface defects in water network are monitored. The fraction of such defect structures is increased near the hydrophobic cap of fullerenes. Some “dry” regions of C{sub 60} are observed which can be considered as signatures of surface dewetting. In an effort to provide molecular level insight into the thermodynamics of hydration, the free energy of solvation is determined for a family of fullerene particles using thermodynamic integration technique.

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

Hydrated proton and hydroxide charge transfer at the liquid/vapor interface of water

The role of the solvated excess proton and hydroxide ions in interfacial properties is an interesting scientific question with applications in a variety of aqueous behaviors. The role that charge transfer (CT) plays in interfacial behavior is also an unsettled question. Quantum calculations are carried out on clusters of water with an excess proton or a missing proton (hydroxide) to determine their CT. The quantum results are applied to analysis of multi-state empirical valence bond trajectories. The polyatomic nature of the solvated excess proton and hydroxide ion results in directionally dependent CT, depending on whether a water molecule is a hydrogen bond donor or acceptor in relation to the ion. With polyatomic molecules, CT also depends on the intramolecular bond distances in addition to intermolecular distances. The hydrated proton and hydroxide affect water’s liquid/vapor interface in a manner similar to monatomic ions, in that they induce a hydrogen-bonding imbalance at the surface, which results in charged surface waters. This hydrogen bond imbalance, and thus the charged waters at the surface, persists until the ion is at least 10 Å away from the interface.

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

Water network-mediated, electron-induced proton transfer in [C{sub 5}H{sub 5}N ⋅ (H{sub 2}O){sub n}]{sup −} clusters

The role of proton-assisted charge accommodation in electron capture by a heterocyclic electron scavenger is investigated through theoretical analysis of the vibrational spectra of cold, gas phase [Py ⋅ (H{sub 2}O){sub n=3−5}]{sup −} clusters. These radical anions are formed when an excess electron is attached to water clusters containing a single pyridine (Py) molecule in a supersonic jet ion source. Under these conditions, the cluster ion distribution starts promptly at n = 3, and the photoelectron spectra, combined with vibrational predissociation spectra of the Ar-tagged anions, establish that for n > 3, these species are best described as hydrated hydroxide ions with the neutral pyridinium radical, PyH{sup (0)}, occupying one of the primary solvation sites of the OH{sup −}. The n = 3 cluster appears to be a special case where charge localization on Py and hydroxide is nearly isoenergetic, and the nature of this species is explored with ab initio molecular dynamics calculations of the trajectories that start from metastable arrangements of the anion based on a diffuse, essentially dipole-bound electron. These calculations indicate that the reaction proceeds via a relatively slow rearrangement of the water network to create a favorable hydration configuration around the water molecule that eventually donates a proton to the Py nitrogen atom to yield the product hydroxide ion. The correlation between the degree of excess charge localization and the evolving shape of the water network revealed by this approach thus provides a microscopic picture of the “solvent coordinate” at the heart of a prototypical proton-coupled electron transfer reaction.

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

Carbon dioxide in an ionic liquid: Structural and rotational dynamics

Ionic liquids (ILs), which have widely tunable structural motifs and intermolecular interactions with solutes, have been proposed as possible carbon capture media. In order to inform the choice of an optimal ionic liquid system, it can be useful to understand the details of dynamics and interactions on fundamental time scales (femtoseconds to picoseconds) of dissolved gases, particularly carbon dioxide (CO 2 ), within the complex solvation structures present in these uniquely organized materials. The rotational and local structural fluctuation dynamics of CO 2 in the room temperature ionic liquid 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide (EmimNTf2) were investigated by using ultrafast infrared spectroscopy to interrogate the CO 2 asymmetric stretch. Polarization-selective pump probe measurements yielded the orientational correlation function of the CO 2 vibrational transition dipole. It was found that reorientation of the carbon dioxide occurs on 3 time scales: 0.91 ± 0.03, 8.3 ± 0.1, 54 ± 1 ps. The initial two are attributed to restricted wobbling motions originating from a gating of CO 2 motions by the IL cations and anions. The final (slowest) decay corresponds to complete orientational randomization. Two-dimensional infrared vibrational echo (2D IR) spectroscopy provided information on structural rearrangements, which cause spectral diffusion, through the time dependence of the 2D line shape. Analysis of the time-dependent 2D IR spectra yields the frequency-frequency correlation function (FFCF). Polarization-selective 2D IR experiments conducted on the CO 2 asymmetric stretch in the parallel- and perpendicular-pumped geometries yield significantly different FFCFs due to a phenomenon known as reorientation-induced spectral diffusion (RISD), revealing strong vector interactions with the liquid structures that evolve slowly on the (independently measured) rotation time scales. To separate the RISD contribution to the FFCF from the structural spectral diffusion contribution, the previously developed first order Stark effect RISD model is reformulated to describe the second order (quadratic) Stark effect—the first order Stark effect vanishes because CO 2 does not have a permanent dipole moment. Through this analysis, we characterize the structural fluctuations of CO 2 in the ionic liquid solvation environment, which separate into magnitude-only and combined magnitude and directional correlations of the liquid’s time dependent electric field. This new methodology will enable highly incisive comparisons between CO 2 dynamics in a variety of ionic liquid systems.

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

Structural, dynamical, and transport properties of the hydrated halides: How do At{sup −} bulk properties compare with those of the other halides, from F{sup −} to I{sup −}?

The properties of halides from the lightest, fluoride (F{sup −}), to the heaviest, astatide (At{sup −}), have been studied in water using a polarizable force-field approach based on molecular dynamics (MD) simulations at the 10 ns scale. The selected force-field explicitly treats the cooperativity within the halide-water hydrogen bond networks. The force-field parameters have been adjusted to ab initio data on anion/water clusters computed at the relativistic Möller-Plesset second-order perturbation theory level of theory. The anion static polarizabilities of the two heaviest halides, I{sup −} and At{sup −}, were computed in the gas phase using large and diffuse atomic basis sets, and taking into account both electron correlation and spin-orbit coupling within a four-component framework. Our MD simulation results show the solvation properties of I{sup −} and At{sup −} in aqueous phase to be very close. For instance, their first hydration shells are structured and encompass 9.2 and 9.1 water molecules at about 3.70 ± 0.05 Å, respectively. These values have to be compared to the F{sup −}, Cl{sup −}, and Br{sup −} ones, i.e., 6.3, 8.4, and 9.0 water molecules at 2.74, 3.38, and 3.55 Å, respectively. Moreover our computations predict the solvation free energy of At{sup −} in liquid water at ambient conditions to be 68 kcal mol{sup −1}, a value also close the I{sup −} one, about 70 kcal mol{sup −1}. In all, our simulation results for I{sup −} are in excellent agreement with the latest neutron- and X-ray diffraction studies. Those for the At{sup −} ion are predictive, as no theoretical or experimental data are available to date.

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

Water anomalous thermodynamics, attraction, repulsion, and hydrophobic hydration

A model composed of van der Waals-like and hydrogen bonding contributions that simulates the low-temperature anomalous thermodynamics of pure water while exhibiting a second, liquid-liquid critical point [P. H. Poole et al., Phys. Rev. Lett. 73, 1632 (1994)] is extended to dilute solutions of nonionic species. Critical lines emanating from such second critical point are calculated. While one infers that the smallness of the water molecule may be a relevant factor for those critical lines to move towards experimentally accessible regions, attention is mainly focused on the picture our model draws for the hydration thermodynamics of purely hydrophobic and amphiphilic non-electrolyte solutes. We first focus on differentiating solvation at constant volume from the corresponding isobaric process. Both processes provide the same viewpoint for the low solubility of hydrophobic solutes: it originates from the combination of weak solute-solvent attractive interactions and the specific excluded-volume effects associated with the small molecular size of water. However, a sharp distinction is found when exploring the temperature dependence of hydration phenomena since, in contrast to the situation for the constant-V process, the properties of pure water play a crucial role at isobaric conditions. Specifically, the solubility minimum as well as enthalpy and entropy convergence phenomena, exclusively ascribed to isobaric solvation, are closely related to water’s density maximum. Furthermore, the behavior of the partial molecular volume and the partial molecular isobaric heat capacity highlights the interplay between water anomalies, attraction, and repulsion. The overall picture presented here is supported by experimental observations, simulations, and previous theoretical results.

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

Anharmonic exciton dynamics and energy dissipation in liquid water from two-dimensional infrared spectroscopy

Water’s extended hydrogen-bond network results in rich and complex dynamics on the sub-picosecond time scale. In this paper, we present a comprehensive analysis of the two-dimensional infrared (2D IR) spectrum of O–H stretching vibrations in liquid H{sub 2}O and their interactions with bending and intermolecular vibrations. By exploring the dependence of the spectrum on waiting time, temperature, and laser polarization, we refine our molecular picture of water’s complex ultrafast dynamics. The spectral evolution following excitation of the O–H stretching resonance reveals vibrational dynamics on the 50–300 fs time scale that are dominated by intermolecular delocalization. These O–H stretch excitons are a result of the anharmonicity of the nuclear potential energy surface that arises from the hydrogen-bonding interaction. The extent of O–H stretching excitons is characterized through 2D depolarization measurements that show spectrally dependent delocalization in agreement with theoretical predictions. Furthermore, we show that these dynamics are insensitive to temperature, indicating that the exciton dynamics alone set the important time scales in the system. Finally, we study the evolution of the O–H stretching mode, which shows highly non-adiabatic dynamics suggestive of vibrational conical intersections. We argue that the so-called heating, commonly observed within ∼1 ps in nonlinear IR spectroscopy of water, is a nonequilibrium state better described by a kinetic temperature rather than a Boltzmann distribution. Our conclusions imply that the collective nature of water vibrations should be considered in describing aqueous solvation.

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

IR spectral assignments for the hydrated excess proton in liquid water

The local environmental sensitivity of infrared (IR) spectroscopy to a hydrogen-bonding structure makes it a powerful tool for investigating the structure and dynamics of excess protons in water. Although of significant interest, the line broadening that results from the ultrafast evolution of different solvated proton-water structures makes the assignment of liquid-phase IR spectra a challenging task. In this paper, we apply a normal mode analysis using density functional theory of thousands of proton-water clusters taken from reactive molecular dynamics trajectories of the latest generation multistate empirical valence bond proton model (MS-EVB 3.2). These calculations are used to obtain a vibrational density of states and IR spectral density, which are decomposed on the basis of solvated proton structure and the frequency dependent mode character. Decompositions are presented on the basis of the proton sharing parameter δ, often used to distinguish Eigen and Zundel species, the stretch and bend character of the modes, the mode delocalization, and the vibrational mode symmetry. We find there is a wide distribution of vibrational frequencies spanning 1200-3000 cm -1 for every local proton configuration, with the region 2000-2600 cm -1 being mostly governed by the distorted Eigen-like configuration. Finally, we find a continuous red shift of the special-pair O∙∙∙H + ∙∙∙O stretching frequency, and an increase in the flanking water bending intensity with decreasing δ. Also, we find that the flanking water stretch mode of the Zundel-like species is strongly mixed with the flanking water bend, and the special pair proton oscillation band is strongly coupled with the bend modes of the central H 5 O 2 + moiety.

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

On combining the conductor-like screening model and optimally tuned range-separated hybrid density functionals

Range-separated hybrid functionals whose range-separation parameter γ has been nonempirically tuned to a particular molecule have been shown to yield frontier orbital energies and other properties in very good agreement with experiments. Yet, many cases, such as organic optoelectronic devices, require the description of molecules embedded in an environment. This can be done by combining the γ-tuning procedure with polarizable continuum models in general and the very versatile conductor-like screening model in particular. There are at least two different ways of performing this combination. The partially vertical γ-tuning employs equilibrium solvation throughout. The strictly vertical γ-tuning, on the other hand, employs nonequilibrium solvation to obtain ionization energies. Here, we compare ground-state and excited-state properties of several different molecules relevant to organic optoelectronics that were obtained using both of the two different tuning procedures. While there are significant differences in the ground-state properties, we see virtually no difference in the excited-state properties. Given these conclusions, we find that both tuning procedures have to be used in conjunction for the correct description of both ground-state and excited-state properties.

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

Molecular insights into ionic liquid/aqueous interface of phosphonium based phase-separable ionic liquids

Phosphonium based phase-separable ionic liquids (PSILs) are promising green solvents for dissolution of cellulose and lignin, a necessary step for conversion of biomass to fuels and chemicals. The knowledge of interfacial behavior of ionic liquid/solvent systems is critical for designing efficient separation processes. Here, molecular dynamics simulations are carried out for aqueous interface of tetraalkylphosphonium ionic liquids with chloride and acetate as anions to investigate IL miscibility with water. The transition zone from miscible to immiscible behavior was observed for alkyl chain lengths of 6 to 8. Emulsion phase was observed for [P8888]+ ion and multiple IL/water interface was observed for [P12121212]+ phosphonium cation. IL/water interface is observed to be enriched with solvated anions with phosphorous atom oriented towards the aqueous phase. The potential of mean force calculations suggest that as the alkyl chain of the phosphonium cations increase, so does the driving force for aggregation of cations in the aqueous phase. Due to the difference in the molecular shape and size, anions are better solvated by the water molecules, and these ions diffuse into the aqueous phase. Furthermore, the analysis of selected pair interactions provides insights into the nature of intermolecular forces and the role of the alkyl side chains on the interfacial properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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

Influence of non-adiabatic effects on linear absorption spectra in the condensed phase: Methylene blue

Modeling linear absorption spectra of solvated chromophores is highly challenging as contributions are present both from coupling of the electronic states to nuclear vibrations and from solute–solvent interactions. In systems where excited states intersect in the Condon region, significant non-adiabatic contributions to absorption line shapes can also be observed. Here, we introduce a robust approach to model linear absorption spectra accounting for both environmental and non-adiabatic effects from first principles. This model parameterizes a linear vibronic coupling (LVC) Hamiltonian directly from energy gap fluctuations calculated along molecular dynamics (MD) trajectories of the chromophore in solution, accounting for both anharmonicity in the potential and direct solute–solvent interactions. The resulting system dynamics described by the LVC Hamiltonian are solved exactly using the thermalized time-evolving density operator with orthogonal polynomials algorithm (T-TEDOPA). The approach is applied to the linear absorption spectrum of methylene blue in water. We show that the strong shoulder in the experimental spectrum is caused by vibrationally driven population transfer between the bright S1 and the dark S2 states. The treatment of the solvent environment is one of many factors that strongly influence the population transfer and line shape; accurate modeling can only be achieved through the use of explicit quantum mechanical solvation. The efficiency of T-TEDOPA, combined with LVC Hamiltonian parameterizations from MD, leads to an attractive method for describing a large variety of systems in complex environments from first principles.

Dunnett, Angus J. (ORCID:0000000285797826)↗
DOE OSTI2023

Multiscale modeling of solute diffusion in triblock copolymer membranes

We develop a multiscale simulation model for diffusion of solutes through porous triblock copolymer membranes. The approach combines two techniques: self-consistent field theory (SCFT) to predict the structure of the self-assembled, solvated membrane and on-lattice kinetic Monte Carlo (kMC) simulations to model diffusion of solutes. Solvation is simulated in SCFT by constraining the glassy membrane matrix while relaxing the brush-like membrane pore coating against the solvent. The kMC simulations capture the resulting solute spatial distribution and concentration-dependent local diffusivity in the polymer-coated pores; we parameterize the latter using particle-based simulations. We apply our approach to simulate solute diffusion through nonequilibrium morphologies of a model triblock copolymer, and we correlate diffusivity with structural descriptors of the morphologies. We also compare the model’s predictions to alternative approaches based on simple lattice random walks and find our multiscale model to be more robust and systematic to parameterize. Furthermore, our multiscale modeling approach is general and can be readily extended in the future to other chemistries, morphologies, and models for the local solute diffusivity and interactions with the membrane.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Dynamic and relativistic effects on Pt–Pt indirect spin–spin coupling in aqueous solution studied by ab initio molecular dynamics and two- vs four-component density functional NMR calculations

Treating 195Pt nuclear magnetic resonance parameters in solution remains a considerable challenge from a quantum chemistry point of view, requiring a high level of theory that simultaneously takes into account the relativistic effects, the dynamic treatment of the solvent–solute system, and the dynamic electron correlation. A combination of Car-Parrinello molecular dynamics (CPMD) and relativistic calculations based on two-component zeroth order regular approximation spin–orbit Kohn–Sham (2c-ZKS) and four-component Dirac–Kohn–Sham (4c-DKS) Hamiltonians is performed to address the solvent effect (water) on the conformational changes and JPtPt1 coupling. A series of bridged PtIII dinuclear complexes [L1–Pt2(NH3)4(Am)2–L2]n+ (Am = α–pyrrolidonate and pivalamidate; L = H2O, Cl−, and Br−) are studied. The computed Pt–Pt coupling is strongly dependent on the conformational dynamics of the complexes, which, in turn, is correlated with the trans influence among axial ligands and with the angle N–C–O from the bridging ligands. The J-coupling is decomposed in terms of dynamic contributions. The decomposition reveals that the vibrational and explicit solvation contributions reduce JPtPt1 of diaquo complexes (L1 = L2 = H2O) in comparison to the static gas-phase magnitude, whereas the implicit solvation and bulk contributions correspond to an increase in JPtPt1 in dihalo (L1 = L2 = X−) and aquahalo (L1 = H2O; L2 = X−) complexes. Relativistic treatment combined with CPMD shows that the 2c-ZKS Hamiltonian performs as well as 4c-DKS for the JPtPt1 coupling.

Chemistry↗
DOE OSTI2025

Key Factors in Semi-Generic Coarse-Grained Modeling of Solid Polymer Electrolytes

Given the long length and time scales of interest and strong interactions present in solid polymer electrolytes, coarse-grained molecular models can be useful in understanding the molecular basis of their structure and transport properties. Rather than representing atoms individually, coarse-grained models use beads representing groups of atoms, making the system simpler and more efficient. Generic coarse-grained models, that are not built based on matching a particular atomistic system, can provide insight into the important physical considerations that apply across different chemical systems, though it can be unclear how to best match to experimental systems and capture relevant experimental behaviors with as simple of a model as possible. Here, we build on generic bead-spring models, but use stiff angle potentials, set different bead properties for different polymer types, and include additional ion parameters, creating a semi-generic coarse-grained model that can map more specifically to polymers and copolymers with different chain architectures, component glass transition temperatures, and ion solvation behaviors. We set most parameters with basic homopolymer data such as glass transition temperatures, Kuhn length, density, and dielectric constant, with further adjustment based on data from the polymer electrolyte system. We specifically model polystyrene-block-poly(oligo-oxyethylene methyl ether methacrylate) (PS-b-POEM) with lithium triflate salt, considering the POEM to be made of a poly(methylmethacrylate) backbone with poly(ethylene oxide) side chains. Solvation of ions is accounted for by additional polymer-ion interactions of the form -S/r4, plus additional lithium-polymer Lennard-Jones interactions. We find these potentials have different effects and discuss strategies for setting these parameters.

Zhang, Yuanhao↗
DOE OSTI2025

Structure, dynamics, and electrochemistry of choline chloride/ethylene glycol eutectic solvents at an electrode surface explored by molecular dynamics simulations

Choline chloride and ethylene glycol mixtures with 1:2, 1:4, and 1:6 molar ratios on the surfaces of graphite and gold electrodes were studied using classical molecular dynamics simulations. Both neutral and charged electrodes were considered. The liquid composition, solvation structure, molecular orientation, and dynamics at the electrode surface are significantly different from those of the bulk liquid. These properties strongly depend on the electrode material and charge density, whereas they are less sensitive to the overall solvent composition. The effect of the electrode on the composition, structure, and orientation of the liquid fades beyond ∼10 Å from the surface of the electrode. This distance corresponds to about two layers of the structured solvent, despite the fact that the layered structure extends to at least five layers or about 25 Å. However, the electrode influences solvent dynamics over a longer distance. The electrochemical properties of the eutectic solvent at both electrode surfaces were also studied. The simulations captured the experimental differential capacitance shapes for both electrode systems, although the magnitudes and exact shapes differ. The simulations further revealed that the solvent in the first solvation layer can both overscreen and underscreen the electrode charges depending on the electrode material and electrode potential.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗