Molecular simulation of liquid water confined inside graphite channels: Thermodynamics and structural properties

Effects of interfaces on structure and dynamics of water droplets on a graphene surface: A molecular dynamics study

The structure and dynamics of water droplets on a bilayer graphene surface are investigated using molecular dynamics simulations. The effects of solid/water and air/water interfaces on the local structure of water droplets are analyzed in terms of the hydrogen bond distribution and tetrahedral order parameter. It is found that the local structure in the core region of a water droplet is similar to that in liquid water. On the other hand, the local structure of water molecules at the solid/water and air/water interfaces, referred to as the interface and surface regions, respectively, consists mainly of three-coordinated molecules that are greatly distorted from a tetrahedral structure. This study reveals that the dynamics in different regions of the water droplets affects the intermolecular vibrational density of states: It is found that in the surface and interface regions, the intensity of vibrational density of states at ∼50 cm^-1 is enhanced, whereas those at ∼200 and ∼500 cm^-1 are weakened and redshifted. These changes are attributed to the increase in the number of molecules having fewer hydrogen bonds in the interface and surface regions. Both single-molecule and collective orientation relaxations are also examined. Single-molecule orientation relaxation is found to be marginally slower than that in liquid water. On the other hand, the collective orientation relaxation of water droplets is found to be significantly faster than that of liquid water because of the destructive correlation of dipole moments in the droplets. The negative correlation between distinct dipole moments also yields a blueshifted libration peak in the absorption spectrum. It is also found that the water-graphene interaction affects the structure and dynamics of the water droplets, such as the local water structure, collective orientation relaxation, and the correlation between dipole moments. This study reveals that the water/solid and water/air interfaces strongly affect the structure and intermolecular dynamics of water droplets and suggests that the intermolecular dynamics, such as energy relaxation dynamics, in other systems with interfaces are different from those in liquid water.

Correlation between mobility and the hydrogen bonding network of water at an electrified-graphite electrode using molecular dynamics simulation

Mobility and hydrogen bonding network of water at a graphite electrode: effects of dissolved ions and applied potential.

Effect of Hydrogen Bonds on the Dielectric Properties of Interfacial Water

The dielectric constant for water is reduced under confinement. Although this phenomenon is well known, the underlying physical mechanism for the reduction is still in debate. In this work, we investigate the effect of the orientation of hydrogen bonds on the dielectric properties of confined water using molecular dynamics simulations. We find a reduced rotational diffusion coefficient for water molecules close to the solid surface. The reduced rotational diffusion arises due to the hindered rotation away from the plane parallel to the channel walls. The suppressed rotation in turn affects the orientational polarization of water, leading to a low value for the dielectric constant at the interface. We attribute the constrained out-of-plane rotation to originate from a higher density of planar hydrogen bonds formed by the interfacial water molecules.

The spatial range of protein hydration

Proteins interact with their aqueous surroundings, thereby modifying the physical properties of the solvent. The extent of this perturbation has been investigated by numerous methods in the past half-century, but a consensus has still not emerged regarding the spatial range of the perturbation. To a large extent, the disparate views found in the current literature can be traced to the lack of a rigorous definition of the perturbation range. Stating that a particular solvent property differs from its bulk value at a certain distance from the protein is not particularly helpful since such findings depend on the sensitivity and precision of the technique used to probe the system. What is needed is a well-defined decay length, an intrinsic property of the protein in a dilute aqueous solution, that specifies the length scale on which a given physical property approaches its bulk-water value. Based on molecular dynamics simulations of four small globular proteins, we present such an analysis of the structural and dynamic properties of the hydrogen-bonded solvent network. The results demonstrate unequivocally that the solvent perturbation is short-ranged, with all investigated properties having exponential decay lengths of less than one hydration shell. The short range of the perturbation is a consequence of the high energy density of bulk water, rendering this solvent highly resistant to structural perturbations. The electric field from the protein, which under certain conditions can be long-ranged, induces a weak alignment of water dipoles, which, however, is merely the linear dielectric response of bulk water and, therefore, should not be thought of as a structural perturbation. By decomposing the first hydration shell into polarity-based subsets, we find that the hydration structure of the nonpolar parts of the protein surface is similar to that of small nonpolar solutes. For all four examined proteins, the mean number of water-water hydrogen bonds in the nonpolar subset is within 1% of the value in bulk water, suggesting that the fragmentation and topography of the nonpolar protein-water interface has evolved to minimize the propensity for protein aggregation by reducing the unfavorable free energy of hydrophobic hydration.

Formation of Water Layers on Graphene Surfaces

Although graphitic materials were thought to be hydrophobic, recent experimental results based on contact angle measurements show that the hydrophobicity of graphitic surfaces stems from airborne contamination of hydrocarbons. This leads us to question whether a pristine graphitic surface is indeed hydrophobic. To investigate the water wettability of graphitic surfaces, we use molecular dynamics simulations of water molecules on the surface of a single graphene layer at room temperature. The results indicate that a water droplet spreads over the entire surface and that a double-layer structure of water molecules forms on the surface, which means that wetting of graphitic surfaces is possible, but only by two layers of water molecules. No further water layers can cohere to the double-layer structure, but the formation of three-dimensional clusters of liquid water is confirmed. The surface of the double-layer structure acts as a hydrophobic surface. Such peculiar behavior of water molecules can be reasonably explained by the formation of hydrogen bonds: The hydrogen bonds of the interfacial water molecules form between the first two layers and also within each layer. This hydrogen-bond network is confined within the double layer, which means that no "dangling hydrogen bonds" appear on the surface of the double-layer structure. This formation of hydrogen bonds stabilizes the double-layer structure and makes its surface hydrophobic. Thus, the numerical simulations indicate that a graphene surface is perfectly wettable on the atomic scale and becomes hydrophobic once it is covered by this double layer of water molecules.

The Structural Properties and Diffusion of a Three-Dimensional Isotropic Core-Softened Model Fluid in Disordered Porous Media. Molecular Dynamics Simulation

The microscopic structure and dynamic properties of an isotropic three-dimensional core-softened model fluid in disordered matrices of Lennard-Jones particles have been studied. Molecular dynamics computer simulations in Grand Canonical ensemble were used as the methodological tools. It was shown that the microscopic structure of the fluid is characterized by anomalies similar to those found in a bulk model, but that it is affected by the fluid-matrix interactions. The dynamic properties also exhibit anomalous dependence on fluid density, but the magnitude of these anomalies is suppressed in comparison to the bulk fluid model. The anomalous behaviour of the diffusion coefficient is attributed to structural changes in the first coordination shell of a given fluid particle. It seems that the anomalies can only be suppressed at matrix densities which are higher than those studied in the present work.

Dielectric constant of water as a function of separation in a slab geometry: A molecular dynamics study

Water in confining geometries shows various anomalous properties related to its structure and dynamics compared with bulk water. Here, the dielectric constant of water as a function of separation in a graphite slab geometry was studied using molecular dynamics simulations. The dielectric constants of water were calculated from the orientational polarization of water molecules when an external electric field was applied parallel and normal to the slabs. The reduction of the dielectric constant of water compared with bulk water can be explained by investigating the structure and dynamics of water in slab geometries. We found a preferred orientation of water molecules in the layer closest to the graphite surface. The self-diffusion coefficient distribution of water molecules along the direction normal to the slabs was also computed. Highly mobile water molecules in the intermediate region were generated by the weak hydrogen bonding produced by the preferred orientation of water molecules in the layer. We concluded that the dielectric constant of water in the slab geometry is lower than that of bulk water because of the reduction of the polarization of water and the highly mobile water molecules in the intermediate region arising from the preferred orientation of water molecules.

Probing the effects of 2D confinement on hydrogen dynamics in water and ice adsorbed in graphene oxide sponges

We studied the single particle dynamics of water and ice adsorbed in graphene oxide (GO) sponges at T = 293 K and T = 20 K. We used Deep Inelastic Neutron Scattering (DINS) at the ISIS neutron and muon spallation source to derive the hydrogen mean kinetic energy, 〈EK〉, and momentum distribution, n(p). The goal of this work was to study the hydrogen dynamics under 2D confinement and the potential energy surface, fingerprinting the hydrogen interaction with the layered structure of the GO sponge. The observed scattering is interpreted within the framework of the impulse approximation. Samples of both water and ice adsorbed in GO show n(p) functions with almost harmonic and anisotropic line shapes and 〈EK〉 values in excess of the values found at the corresponding temperatures in the bulk. The hydrogen dynamics are discussed in the context of the interaction between the interfacial water and ice and the confining hydrophilic surface of the GO sponge.

Computational chemistry for graphene-based energy applications: progress and challenges

Research in graphene-based energy materials is a rapidly growing area. Many graphene-based energy applications involve interfacial processes. To enable advances in the design of these energy materials, such that their operation, economy, efficiency and durability is at least comparable with fossil-fuel based alternatives, connections between the molecular-scale structure and function of these interfaces are needed. While it is experimentally challenging to resolve this interfacial structure, molecular simulation and computational chemistry can help bridge these gaps. In this Review, we summarise recent progress in the application of computational chemistry to graphene-based materials for fuel cells, batteries, photovoltaics and supercapacitors. We also outline both the bright prospects and emerging challenges these techniques face for application to graphene-based energy materials in future.

Size effects on water adsorbed on hydrophobic probes at the nanometric scale

Molecular dynamics simulations of liquid water at ambient conditions, adsorbed at the external walls of (n,n) single-walled armchair carbon nanotubes have been performed for n = 5, 9, 12. The comparison with the case of water adsorbed on graphene has also been included. The analysis of Helmholtz free energies reveals qualitatively different ranges of thermodynamical stability, eventually starting at a given threshold surface density. We observed that, in the framework of the force field considered here, water does not wet graphene nor (12,12) tubes, but it can coat thinner tubes such as (9,9) and (5,5), which indicates that the width of the carbon nanotube plays a role on wetting. On the other hand, density profiles, orientational distributions of water, and hydrogen-bond populations indicate significant changes of structure of water for the different surfaces. Further, we computed self-diffusion of water and spectral densities of water and carbon molecules, which again revealed different qualitative behavior of interfacial water depending on the size of the nanotube. The crossover size corresponds to tube diameters of around 1 nm.

Squeezing water clusters between graphene sheets: energetics, structure, and intermolecular interactions

The behavior of water confined at the nanoscale between graphene sheets has attracted much theoretical and experimental attention recently.

Tunable giant anisotropic diffusion of water sub-monolayers between graphene layers

We investigate the in-plane confinement effect of two graphene layers on the diffusion behaviour of water sub-monolayers using molecular dynamics simulations. An unexpected fast diffusion state with giant anisotropy is observed when the two confining graphene walls have certain shifts applied to their relative positions. The phenomenon is mainly attributed to the smooth one-dimensional potential channels produced by the composition effect of the potential energy landscapes of the two graphene walls, and the concerted motion of water molecules due to hydrogen bonding. Unique duality in the diffusion mechanism is observed in the fast diffusion state, as is ballistic motion along the potential channels and Fickian diffusion across such channels. The smooth potential channels can be created in certain directions simply by shifting the confining walls, which provides a novel measure to manipulate the motion of confined molecules in real-time.

The interfacial-organized monolayer water film (MWF) induced "two-step" aggregation of nanographene: both in stacking and sliding assembly pathways

A computational investigation was carried out to understand the aggregation of nanoscale graphene with two typical pathways of stacking assembly and sliding assembly in water. The interfacial-organized monolayer water film (MWF) induced "two-step" aggregation of nanographene in both stacking and sliding assembly pathways was reported for the first time. By means of potential mean forces (PMFs) calculation, no energy barrier was observed during the sliding assembly of two graphene nanosheets, while the PMF profiles could be impacted by the contact forms of nanographene and the MWF within the interplate of two graphene nanosheets. To explore the potential physical basis of the "hindering role" of self-organized interfacial water, the dynamical and structural properties as well as the status of hydrogen bonds (H-bonds) for interfacial water were investigated. We found that the compact, ordered structure and abundant H-bonds of the MWF could be taken as the fundamental aspects of the "hindering role" of interfacial water for the hydrophobic assembly of nanographene. These findings are displaying a potential to further understand the hydrophobic assembly which mostly dominate the behaviors of nanomaterials, proteins etc. in aqueous solutions.

Confinement of anomalous liquids in nanoporous matrices

Using molecular dynamics simulations, we investigate the effects of different nanoconfinements on complex liquids-e.g., colloids or protein solutions-with density anomalies and a liquid-liquid phase transition (LLPT). In all the confinements, we find a strong depletion effect with a large increase in liquid density near the confining surface. If the nanoconfinement is modeled by an ordered matrix of nanoparticles, we find that the anomalies are preserved. On the contrary, if the confinement is modeled by a disordered matrix of nanoparticles, we find a drastically different phase diagram: the LLPT shifts to lower pressures and temperatures, and the anomalies become weaker, as the disorder increases. We find that the density heterogeneities induced by the disordered matrix are responsible for the weakening of the LLPT and the disappearance of the anomalies.

Specific ion effects in aqueous eletrolyte solutions confined within graphene sheets at the nanometric scale

The underlying mechanisms of specific ion effects on structure and dynamics of aqueous solutions have been long debated. On the other hand, the role of polarization at hydrophobic interfaces when aqueous electrolytes are present is of great importance, as it has been observed at the air-vapor interface. In this work, we have explored influence of ionic species on microscopical properties of aqueous sodium halide solutions constrained inside a double layer graphene channel, as a model for a realistic hydrophobic interface. Our systems have been simulated by molecular dynamics techniques, explicitly including polarization in water molecules and ions. Water and ionic density profiles showed the tendency of ionic species to occupy the whole space available, in good agreement with spectroscopic experimental data. The exception to this general behavior was fluoride, which preferred to stay away from interfaces. Two main regions were defined: interfaces and the central part of the slab, the bulklike region. Ionic hydration numbers at interfaces were lower than those at the bulklike area by about one to two units. We have also analyzed water-ion orientations and polarization distributions and obtained a marked dependence on ionic concentration. Residence time of anions suffered important fluctuations and tended to be largest at interfaces. Large variations of the static permittivity between interfacial and bulklike regions were observed. Ionic diffusion was found to be between 10(-5) and 10(-6) cm(2) s(-1) and showed to be mainly dependent on the concentration, whereas the type of anion considered and the polarizability had significantly less relevance. Conductivities were found to be dependent on ionic concentrations and the polarizabilities of anions, as well as on the spatial direction considered.

Wetting and prewetting of water on top of a single sheet of hexagonal boron nitride

Wetting of a single hexagonal boron nitride sheet by liquid water has been investigated by molecular dynamics simulations within a temperature range between 278 and 373 K. The wetting temperature was found to be ~310 K, while the onset of prewetting happens around the much higher temperature of 354 K. The static (hydrogen-bond populations, density profiles, energy per molecule) and dynamic (diffusion coefficients) properties of water in the stable phases in this temperature range were also studied and compared to those of water on graphene. The results indicate that hydrophobicity of boron nitride is milder than that of graphene.

Single-particle and collective dynamics of methanol confined in carbon nanotubes: a computer simulation study

We present the results of computer simulations of methanol confined in carbon nanotubes. Different levels of confinement were identified as a function of the nanotube radius and characterized using a pair-distribution function adapted to the cylindrical geometry of these systems. Dynamical properties of methanol were also analysed as a function of the nanotube size, both at the level of single-particle and collective properties. We found that confinement in narrow carbon nanotubes strongly affects the dynamical properties of methanol with respect to the bulk phase, due to the strong interaction with the carbon nanotube. In the other cases, confined methanol shows properties quite similar to those of the bulk phase. These phenomena are related to the peculiar hydrogen bonded network of methanol and are compared to the behaviour of water confined in similar conditions. The effect of nanotube flexibility on the dynamical properties of confined methanol is also discussed.

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.

Molecular simulation of water confined in nanoporous silica

This paper reports on a molecular simulation study of the thermodynamics, structure and dynamics of water confined at ambient temperature in hydroxylated silica nanopores of a width H = 10 and 20 Å. The adsorption isotherms for water in these nanopores resemble those observed for experimental samples; the adsorbed amount increases continuously in the multilayer adsorption regime until a jump occurs due to capillary condensation of the fluid within the pore. Strong layering of water in the vicinity of the silica surfaces is observed as marked density oscillations are observed up to 8 Å from the surface in the density profiles for confined water. Our results indicate that water molecules within the first adsorbed layer tend to adopt a H-down orientation with respect to the silica substrate. For all pore sizes and adsorbed amounts, the self-diffusivity of confined water is lower than the bulk, due to the hydrophilic interaction between the water molecules and the hydroxylated silica surface. Our results also suggest that the self-diffusivity of confined water is sensitive to the adsorbed amount.

Collective properties of water confined in carbon nanotubes: A computer simulation study

The collective properties of water confined in the (10,10), (8,8) and (6,6) carbon nanotubes are studied by analysing the longitudinal-current autocorrelation function, calculated from computer-simulated trajectories. The corresponding spectra clearly show the presence of two excitations, but their behaviour is quite different from that observed in the case of bulk water. Instead of the strong positive dispersion of the hydrodynamic sound mode characteristic of bulk water (the fast-sound phenomenon), the sound dispersion relation of confined water is observed to flatten into a non-propagating mode, while a second excitation appears at a higher frequency. This behaviour is analysed in terms of the localized oscillation modes of the hydrogen-bond network.

Description of ferrocenylalkylthiol SAMs on gold by molecular dynamics simulations

Molecular dynamics simulations of mixed monolayers consisting of Fc(CH2)12S-/C10S-Au SAMs are carried out to calculate structural (density profiles, angular distributions, positions of atoms) and energetic properties. The purpose of this paper is to explore the possible inhomogeneity of the neutral ferrocene moieties within the monolayer. Five systems have been studied using different grafting densities for the ferrocenylalkylthiolates. The angular distributions are described in terms of the relative contributions from isolated and clustered ferrocene moieties in the binary SAMs. It is shown that the energetic contributions strongly depend on the state of the ferrocene. The ability of molecular dynamics simulations to enable better understanding the SAM structure is illustrated in this work.

Hydration structure on crystalline silica substrates

The structure of interfacial water at the silica solid surfaces was investigated using molecular dynamics simulations. Different degrees of surface hydroxylation were employed to assess the effect of the surface chemistry on the structure of interfacial water. Density profiles, in-plane radial distribution functions, in-plane density distribution, and hydrogen-bond profiles were calculated. Our results show that the surface hydroxylation affects the structure, orientation, and hydrogen-bond network of interfacial water molecules. Data analysis suggests that the degree of hydroxylation controls the amount of water molecules in the first interfacial layer as well as the distance between the first adsorbed layer and the substrate. Well-organized and uniform structures of interfacial water appear on the homogeneously hydroxylated surface, while a heterogeneous interfacial structure, characterized by extensive water-water hydrogen bonds, forms on the partially hydroxylated surface. We demonstrate that both the local surface chemistry and water-water hydrogen bonds are the dominant factors that determine the structural properties of interfacial water.

Molecular dynamics simulations of supercritical water confined within a carbon-slit pore

We report the results of a series of molecular dynamics simulations of water inside a carbon-slit pore at supercritical conditions. A range of densities corresponding from liquid (0.66gcm;{-3}) to gas environments (0.08gcm;{-3}) at the supercritical temperature of 673K were considered. Our findings are compared with previous studies of liquid water confined in graphene nanochannels at ambient and high temperatures, and indicate that the microscopic structure of water evolves from hydrogen bond networks characteristic of hot dense liquids to looser arrangements where the dominant units are water monomers and dimers. Water permittivity was found to be very small at low densities, with a tendency to grow with density and to reach typical values of unconfined supercritical water at 0.66gcm;{-3}) . In supercritical conditions, the residence time of water at interfaces is roughly similar to that of water in the central regions of the slabs, if the size of the considered region is taken into account. That time span is long enough to compute dynamical properties such as diffusion or spectral densities. Water diffusion in supercritical states is much faster at low densities, and it is produced in such a way that, at interfaces, translational diffusion is mainly produced along planes parallel to the carbon walls. Spectral frequency shifts depend on several factors, being temperature and density effects the most relevant. However, we can observe corrections due to confinement, important both at the graphene interface and in the central region of the water slab.

Liquid water confined in carbon nanochannels at high temperatures

Structure, hydrogen bonding, electrostatics, dielectric, and dynamical properties of liquid water confined in flat graphene nanochannels are investigated by molecular dynamics simulations. A wide range of temperatures (between 20 and 360 degrees C) have been considered. Molecular structure suffers substantial changes when the system is heated, with a significant loss of structure and hydrogen bonding. In such case, the interface between adsorbed and bulk-like water has a marked tendency to disappear, and the two preferential orientations of water nearby the graphite layers at room temperature are essentially merging above the boiling point. The general trend for the static dielectric constant is its reduction at high temperature states, as compared to ambient conditions. Similarly, residence times of water molecules in adsorbed and bulk-like regions are significantly influenced by temperature, as well. Finally, we observed relevant changes in water diffusion and spectroscopy along the range of temperatures analyzed.

Effects of Confinement on Small Water Clusters Structure and Proton Transport

Analyses of the structure of two to four water molecule clusters confined between two benzene and between two naphthalene molecules have been performed using ab initio methods. The water clusters tend to maximize the number of hydrogen bonds via formation of a cyclic network. The oxygen atoms locate approximately in the middle of the confined geometry, and the dipole vectors arrange either parallel or pointing to the surfaces. Energy barriers for proton transfer calculated for H3O+-(H2O) complexes in the same confined geometries suggest that there is a specific range of confinement that helps to lower the energy barriers of the proton transfer. When the walls are too close to each other, at a separation of 4 A, the energy barriers are extremely high. Confinement does not lower the barrier energies of proton transfer when the H3O+-(H2O) complexes are located further from each of the surfaces by more than approximately 8 A.