Solvation Dynamics at the Water/Zirconia Interface: Molecular Dynamics Simulations

The sensitive aspects of modelling polymer-ceramic composite solid-state electrolytes using molecular dynamics simulations

Solid-state composite electrolytes have arisen as one of the most promising materials classes for next-generation Li-ion battery technology. These composites mix ceramic and solid-polymer ion conductors with the aim of combining the advantages of each material. The ion-transport mechanisms within such materials, however, remain elusive. This knowledge gap can to a large part be attributed to difficulties in studying processes at the ceramic-polymer interface, which are expected to play a major role in the overall ion transport through the electrolyte. Computational efforts have the potential of providing significant insight into these processes. One of the main challenges to overcome is then to understand how a sufficiently robust model can be constructed in order to provide reliable results. To this end, a series of molecular dynamics simulations are here carried out with a variation of certain structural (surface termination and polymer length) and pair potential (van der Waals parameters and partial charges) models of the Li7La3Zr2O12 (LLZO) poly(ethylene oxide) (PEO) system, in order to test how sensitive the outcome is to each variation. The study shows that the static and dynamic properties of Li-ion are significantly affected by van der Waals parameters as well as the surface terminations, while the thickness of the interfacial region - where the structure-dynamic properties are different as compared to the bulk-like regime - is the same irrespective of the simulation setup.

Deciphering the molecular mechanism of water boiling at heterogeneous interfaces

Water boiling control evolution of natural geothermal systems is widely exploited in industrial processes due to the unique non-linear thermophysical behavior. Even though the properties of water both in the liquid and gas state have been extensively studied experimentally and by numerical simulations, there is still a fundamental knowledge gap in understanding the mechanism of the heterogeneous nucleate boiling controlling evaporation and condensation. In this study, the molecular mechanism of bubble nucleation at the hydrophilic and hydrophobic solid-water interface was determined by performing unbiased molecular dynamics simulations using the transition path sampling scheme. Analyzing the liquid to vapor transition path, the initiation of small void cavities (vapor bubbles nuclei) and their subsequent merging mechanism, leading to successively growing vacuum domains (vapor phase), has been elucidated. The molecular mechanism and the boiling nucleation sites' location are strongly dependent on the solid surface hydrophobicity and hydrophilicity. Then simulations reveal the impact of the surface functionality on the adsorbed thin water molecules film structuring and the location of high probability nucleation sites. Our findings provide molecular-scale insights into the computational aided design of new novel materials for more efficient heat removal and rationalizing the damage mechanisms.

Water on graphene surfaces

In this paper, we summarize the main results obtained in our group about the behavior of water confined inside or close to different graphene surfaces by means of molecular dynamics simulations. These include the inside and outside of carbon nanotubes, and the confinement inside a slit pore or a single graphene sheet. We paid special attention to some thermodynamical (binding energies), structural (hydrogen-bond distributions) and dynamic (infrared spectra) properties, and their comparison to their bulk counterparts.

PACKMOL: a package for building initial configurations for molecular dynamics simulations

Adequate initial configurations for molecular dynamics simulations consist of arrangements of molecules distributed in space in such a way to approximately represent the system's overall structure. In order that the simulations are not disrupted by large van der Waals repulsive interactions, atoms from different molecules must keep safe pairwise distances. Obtaining such a molecular arrangement can be considered a packing problem: Each type molecule must satisfy spatial constraints related to the geometry of the system, and the distance between atoms of different molecules must be greater than some specified tolerance. We have developed a code able to pack millions of atoms, grouped in arbitrarily complex molecules, inside a variety of three-dimensional regions. The regions may be intersections of spheres, ellipses, cylinders, planes, or boxes. The user must provide only the structure of one molecule of each type and the geometrical constraints that each type of molecule must satisfy. Building complex mixtures, interfaces, solvating biomolecules in water, other solvents, or mixtures of solvents, is straightforward. In addition, different atoms belonging to the same molecule may also be restricted to different spatial regions, in such a way that more ordered molecular arrangements can be built, as micelles, lipid double-layers, etc. The packing time for state-of-the-art molecular dynamics systems varies from a few seconds to a few minutes in a personal computer. The input files are simple and currently compatible with PDB, Tinker, Molden, or Moldy coordinate files. The package is distributed as free software and can be downloaded from http://www.ime.unicamp.br/~martinez/packmol/.

Water/hydrocarbon interfaces: effect of hydrocarbon branching on single-molecule relaxation

Water/hydrocarbon interfaces are studied using molecular dynamics simulations in order to understand the effect of hydrocarbon branching on the dynamics of the system at and away from the interface. A recently proposed procedure for studying the intrinsic structure of the interface in such systems is utilized, and dynamics are probed in the usual laboratory frame as well as the intrinsic frame. The use of these two frames of reference leads to insight into the effect of capillary waves at the interface on dynamics. The systems were partitioned into zones with a width of 5 A, and a number of quantities of dynamical relevance, namely, the residence times, mean squared displacements, the velocity auto correlation functions, and orientational time correlations for molecules of both phases, were calculated in the laboratory and intrinsic frames at and away from the interface. For the aqueous phase, translational motion is found to be (a) diffusive at long times and not anomalous as in proteins or micelles, (b) faster at the interface than in the bulk, and (c) faster upon reduction of the effect of capillary waves. The rotational motion of water is (a) more anisotropic at the interface than in the bulk and (b) dependent on the orientation of the covalent O-H bond with respect to the plane of the interface. The effect of hydrocarbon branching on aqueous dynamics was found to be small, a result similar to the effect on the interfacial water structure. The hydrocarbon phase shows a larger variation for all dynamical probes, a trend consistent with their interfacial structure.

Solvation of coumarin 314 at water/air interfaces containing anionic surfactants. I. Low coverage

Through the use of molecular dynamics techniques, we analyze equilibrium and dynamical aspects of the solvation of Coumarin 314 adsorbed at water/air interfaces in the presence of sodium dodecyl sulfate surfactant molecules. Three different coverages in the submonolayer regime were considered, 500, 250, and 100 A(2)/SDS molecule. The surfactant promotes two well-differentiated solvation environments, which can be clearly distinguished in terms of their structures for the largest surfactant coverage considered. The first one is characterized by the probe lying adjacent or exterior to two-dimensional spatial domains formed by clustered surfactant molecules. A second type of solvation environment is found in which the coumarin appears embedded within compact surfactant domains. Equilibrium and dynamical aspects of the interfacial orientation of the probe are investigated. Our results show a gradual transition from parallel to perpendicular dipolar alignment of the probe with respect to the interface as the concentration of surfactant rho(s) increases. The presence of the surfactant leads to an increase in the roughness and in the characteristic width of the water/air interface. These modifications are also manifested by the decorrelation times for the probe reorientational dynamics, which become progressively slower with rho(s) in both solvation states, although much more pronounced for the embedded ones. The dynamical characteristics of the solvation responses of the charged interfaces are also analyzed, and the implications of our findings to the interpretation of available experimental measurements are discussed.

Molecular Dynamics Simulation of Liquid Water Confined inside Graphite Channels: Dielectric and Dynamical Properties

Electric and dielectric properties and microscopic dynamics of liquid water confined between graphite slabs are analyzed by means of molecular dynamics simulations for several graphite-graphite separations at ambient conditions. The electric potential across the interface shows oscillations due to water layering, and the overall potential drop is about -0.28 V. The total dielectric constant is larger than the corresponding value for the bulklike internal region of the system. This is mainly due to the preferential orientations of water nearest the graphite walls. Estimation of the capacitance of the system is reported, indicating large variations for the different adsorption layers. The main trend observed concerning water diffusion is 2-fold: on one hand, the overall diffusion of water is markedly smaller for the closest graphite-graphite separations, and on the other hand, water molecules diffuse in interfaces slightly slower than those in the bulklike internal areas. Molecular reorientational times are generally larger than those corresponding to those of unconstrained bulk water. The analysis of spectral densities revealed significant spectral shifts, compared to the bands in unconstrained water, in different frequency regions, and associated to confinement effects. These findings are important because of the scarce information available from experimental, theoretical, and computer simulation research into the dielectric and dynamical properties of confined water.

Aqueous Solvation Dynamics at Metal Oxide Surfaces

Broadband transient absorption (TA) spectroscopy, three-pulse photon echo peak shift (3PEPS), and anisotropy decay measurements were used to study the solvation dynamics in bulk water and interfacial water at ZrO(2) surfaces, using Eosin Y as a probe. The 3PEPS results show a multiexponential behavior with two subpicosecond components that are similar in bulk and interfacial water, while a third component of several picoseconds is significantly lengthened at the interface. The bandwidth correlation function from TA spectra exhibits the same behavior, and the TA spectra are well reproduced using the doorway-window picture with the time constants from PEPS. Our results suggest that interfacial water is restricted to a thickness of less than 5 A. Also the high-frequency collective dynamics of water does not seem to be affected by the interface. On the other hand, the increase of the third component may point to a slowing down of diffusional motion at the interface, although other effects, may play a role, which are discussed.

Molecular origin of time-dependent fluorescence shifts in proteins

Time-resolved fluorescence spectroscopy is used increasingly to probe molecular motions at the aqueous interfaces of biological macromolecules and membranes. By recording the time variation of the fluorescence frequency, thermal atomic fluctuations in the vicinity of the chromophore can be probed. From such fluorescence Stokes shift (FSS) experiments, it has been inferred that water motions in the hydration layer are slowed down by 1-3 orders of magnitude. To provide a more secure foundation for the interpretation of FSS data, we use molecular dynamics simulations to examine the molecular origin of the FSS from a tryptophan residue in a protein. By using linear response theory to decompose the FSS into its water and protein components, we find that the water component dominates the static FSS but decays rapidly. Thus, after a few picoseconds, the FSS essentially reflects protein dynamics, including the self-motion of the chromophore. Because of its collective nature, the FSS response is insensitive to the motion of individual water molecules. Collective water displacement by slowly fluctuating protein groups introduces a long-time tail in the water autocorrelation function, but this dynamic coupling is hardly manifested in the observed FSS. Our analysis reconciles FSS data with the picture of a highly dynamic hydration layer, derived mainly from magnetic relaxation dispersion and simulation studies, and calls for a revision of previous interpretations of FSS decays in terms of slow hydration dynamics at biomolecular and other interfaces.

Solvation Dynamics in Reverse Micelles: The Role of Headgroup−Solute Interactions

We present molecular dynamics simulation results for solvation dynamics in the water pool of anionic-surfactant reverse micelles (RMs) of varying water content, w(0). The model RMs are designed to represent water/aerosol-OT/oil systems, where aerosol-OT is the common name for sodium bis(2-ethylhexyl)sulfosuccinate. To determine the effects of chromophore-headgroup interactions on solvation dynamics, we compare the results for charge localization in model ionic diatomic chromophores that differ only in charge sign. Electronic excitation in both cases is modeled as charge localization on one of the solute sites. We find dramatic differences in the solvation responses for anionic and cationic chromophores. Solvation dynamics for the cationic chromophore are considerably slower and more strongly w(0)-dependent than those for the anionic chromophore. Further analysis indicates that the difference in the responses can be ascribed in part to the different initial locations of the two chromophores relative to the surfactant interface. In addition, slow motion of the cationic chromophore relative to the interface is the main contributor to the longer-time decay of the solvation response to charge localization in this case.