Hydrophobic interaction and hydrogen-bond network for a methane pair in liquid water

Segregation of Chemical Groups at PMMA/H2O Interface Leads to Different Local Hydrophobicity

The local hydrophilicity of polymer surfaces is essential for many applications, such as coatings and biocompatibility, but revealing its structural origin is challenging. Here, we used poly(methyl methacrylate) (PMMA) film as a model and excavated several SFG spectral features, generated by femtosecond sum-frequency generation vibrational spectroscopy (SFG-VS), to elucidate the nature of microscopic hydrophilicity of a polymer/H2O interface. For the first time, we successfully probed the SFG spectra of the bend-libration combination band of interfacial water, which exhibits high sensitivity to solvent-water interactions. Two local hydrophilic domains are observed at the PMMA/H2O interface. The segregation of -OCH3 and -CH3 groups to the PMMA/H2O interface results in the formation of a local hydrophobic domain, where weak solvent-water interactions, slow vibrational dynamics of OH stretching, and no ice-like interfacial water are detected. In contrast, when both C═O and -OCH3 groups segregate to the PMMA/water interface, a local hydrophilic domain is formed, leading to strong solvent-water interactions, fast vibrational dynamics of OH stretching, and the presence of ice-like interfacial water. The water molecules around the hydrophobic domains of the PMMA surface are mainly liquid-like water rather than ice-like water. This work contributes to a molecular-level understanding of the local hydrophilicity of polymer surfaces.

Three-body interaction of gold nanoparticles: the role of solvent density and ligand shell orientation

Molecular dynamics simulations are used to study the effective interactions of alkanethiol passivated gold nanoparticles in supercritical ethane at two- and three-particle levels with different solvent densities. Effective interaction is calculated as the potential of mean force (PMF) between two nanoparticles, and the three-body effect is estimated as the difference in PMFs calculated at the two- and three-particle levels. The variation in the three-body effect is examined as a function of solvent density. It is found that effective interaction, which is completely repulsive at very high solvent concentrations, progressively turns attractive as solvent density declines. On the other hand, the three-body effect turns out to be repulsive and increases exponentially with decreasing solvent density. Further, the structure of the ligand shell is analyzed as a function of nanoparticle separation, and its relationship with the three-body effect is investigated. It is observed that the three-body effect arises when the ligand shell begins to deform due to van der Waals repulsion between ligand shells. The study provides a deep insight into good understanding of the solvent evaporation-assisted nanoparticle self-assembly and can aid in experiments.

Early Steps in the O2 Scavenger Process in the Aqueous Phase: Hydrazine vs DEHA

We investigate the first direct proton abstraction reactions from reducing agents (RAHs) hydrazine and diethyl hydroxylamine (DEHA), toward dioxygen (O2) in the aqueous phase, spanning ambient to high-temperature conditions. Quantum chemistry methods and molecular dynamics simulations are employed in this study. Quantum chemistry methods are used to analyze the quasi-equilibrium between a reactive conformation and a transition state in the [RAH,O2] cluster. On the other hand, molecular dynamics simulations estimate the probability of observing a reactive conformation of the [RAH,O2] cluster in the solution. In this study, we assume that the energy barrier of the quasi-equilibrium is sufficiently high for the RAH/O2 association process to be at equilibrium. Our findings indicate that the first proton abstraction process from a reactive conformation cluster by DEHA is energetically favored compared to hydrazine. Conversely, the association process of hydrazine and O2 in solution is more favorable than that of DEHA. Consequently, the rate constant for the first proton abstraction process is similar for both hydrazine and DEHA, particularly at high temperatures, with activation energies of approximately 21.5 ± 1.5 kcal mol-1 for both compounds. These results align with recent experiments investigating the complete O2 scavenger process in liquid water with hydrazine and DEHA. Therefore, our findings support the assumption that first proton abstraction reactions are the rate-determining steps in O2 scavenger processes in the aqueous phase.

Speciation of La3+-Cl- Complexes in Hydrothermal Fluids from Deep Potential Molecular Dynamics

The stability of rare earth element (REE) complexes plays a crucial role in quantitatively assessing their hydrothermal migration and transformation. However, reliable data are lacking under high-temperature hydrothermal conditions, which hampers our understanding of the association behavior of REE. Here a deep learning potential model for the LaCl3-H2O system in hydrothermal fluids is developed based on the first-principles density functional theory calculations. The model accurately predicts the radial distribution functions compared to ab initio molecular dynamics (AIMD) simulations. Furthermore, species of La-Cl complexes, the dissociation pathway of the La-Cl complexes dissociation process, and the potential of mean forces and corresponding association constants (logK) for LaCln3-n (n = 1-4) are extensively investigated under a wide range of temperatures and pressures. Empirical density models for logK calculation are fitted with these data and can accurately predict logK data from both experimental results and AIMD simulations. The distribution of La-Cl species is also evaluated across a wide range of temperatures, pressures, and initial chloride concentration conditions. The results show that La-Cl complexes are prone to forming in a low-density solution, and the number of bonded Cl- ions increases with rising temperature. In contrast, in a high-density solution, La3+ dominates and becomes the more prevalent species.

Water model for hydrophobic cavities: structure and energy from quantum-chemical calculations

This ab initio study aims to design a series of large water clusters having a hollow clathrate-like cage able to host hydrophobic solutes of various sizes. Starting from the (H₂O)_n (n = 18, 20, 24 and 28) hollow cages, water layers have been added in a stepwise manner in order to model the configuration of water molecules beyond the primary shell. The large (H₂O)₁₀₀, (H₂O)₁₂₀ and (H₂O)₁₄₀ clusters complete the hydrogen bonding network of the cage with optimal and regular tiling of the do-, tetra-decahedron and hexa-decahedron, respectively. This study is corroborated by an investigation of dense water clusters up to the (H₂O)₁₂₃ one, being highly consistent with experimental data on ice concerning the electronic and zero-point energies for aggregate formation at 0 K and enthalpy and entropy at 273 K. The cavity creation profoundly alters the orientation of water molecules compared with those found in dense clusters. Nevertheless, such a large reorganization is necessary to maximize the water-water attraction by making it similar to the one found in dense clusters. The cage formation is an endothermic process; however, the computed values are large compared with previous reports for hydrocarbon aqueous solutions. Larger clusters are required for a more fruitful comparison.

On the Trail of Molecular Hydrophilicity and Hydrophobicity at Aqueous Interfaces

Uncovering microscopic hydrophilicity and hydrophobicity at heterogeneous aqueous interfaces is essential as it dictates physico/chemical properties such as wetting, the electrical double layer, and reactivity. Several molecular and spectroscopic descriptors were proposed, but a major limitation is the lack of connections between them. Here, we combine density functional theory-based MD simulations (DFT-MD) and SFG spectroscopy to explore how interfacial water responds in contact with self-assembled monolayers (SAM) of tunable hydrophilicity. We introduce a microscopic metric to track the transition from hydrophobic to hydrophilic interfaces. This metric combines the H/V descriptor, a structural descriptor based on the preferential orientation within the water network in the topmost binding interfacial layer (BIL) and spectroscopic fingerprints of H-bonded and dangling OH groups of water carried by BIL-resolved SFG spectra. This metric builds a bridge between molecular descriptors of hydrophilicity/hydrophobicity and spectroscopically measured quantities and provides a recipe to quantitatively or qualitatively interpret experimental SFG signals.

The effect of a mixture of an ionic liquid and organic solvent on oxygen reduction reaction kinetics

Li-O₂ batteries attract great attention due to their promising theoretical energy density. One of the main obstacles on the way to achieving high energy density and good cyclability is positive electrode passivation by the Li₂O₂ discharge product as well as the presence of parasitic reactions that degrade electrode and electrolyte materials. To overcome these issues new electrolytes are being extensively searched for to ensure the bulk-mediated mechanism of the oxygen reduction reaction and inhibition of parasitic reactions. Different additives to organic solvents can significantly change the properties of electrolytes. This work is devoted to the effect of ionic liquids (ILs), which are proposed as an additive to the solvent due to their excellent solvation properties, high stability, low volatility and flammability. Using molecular dynamics simulations we investigate mixtures of the Pyr₁₄TFSI ionic liquid and dimethoxyethane (DME) with different volume fractions of the IL. Our calculations show that the presence of the ionic liquid in the electrolyte stabilises solvation shells around the ions, both involved in the oxygen reduction and parasitic reactions, slowing down the kinetics of Li⁺ and O₂^- association. This makes the usage of such mixtures promising for electrolyte design for Li-O₂ batteries.

The behavior of methane-water mixtures under elevated pressures from simulations using many-body potentials

Non-polarizable empirical potentials have been proven to be incapable of capturing the mixing of methane-water mixtures at elevated pressures. Although density functional theory-based ab initio simulations may circumvent this discrepancy, they are limited in terms of the relevant time and length scales associated with mixing phenomena. Here, we show that the many-body MB-nrg potential, designed to reproduce methane-water interactions with coupled cluster accuracy, successfully captures this phenomenon up to 3 GPa and 500 K with varying methane concentrations. Two-phase simulations and long time scales that are required to fully capture the mixing, affordable due to the speed and accuracy of the MBX software, are assessed. Constructing the methane-water equation of state across the phase diagram shows that the stable mixtures are denser than the sum of their parts at a given pressure and temperature. We find that many-body polarization plays a central role, enhancing the induced dipole moments of methane by 0.20 D during mixing under pressure. Overall, the mixed system adopts a denser state, which involves a significant enthalpic driving force as elucidated by a systematic many-body energy decomposition analysis.

Three-body aggregation of guest molecules as a key step in methane hydrate nucleation and growth

Gas hydrates have an important role in environmental and astrochemistry, as well as in energy materials research. Although it is widely accepted that gas accumulation is an important and necessary process during hydrate nucleation, how guest molecules aggregate remains largely unknown. Here, we have performed molecular dynamics simulations to clarify the nucleation path of methane hydrate. We demonstrated that methane gather with a three-body aggregate pattern corresponding to the free energy minimum of three-methane hydrophobic interaction. Methane molecules fluctuate around one methane which later becomes the central gas molecule, and when several methanes move into the region within 0.8 nm of the potential central methane, they act as directional methane molecules. Two neighbor directional methanes and the potential central methane form a three-body aggregate as a regular triangle with a distance of ~6.7 Å which is well within the range of typical methane-methane distances in hydrates or in solution. We further showed that hydrate nucleation and growth is inextricably linked to three-body aggregates. By forming one, two, and three three-body aggregates, the possibility of hydrate nucleation at the aggregate increases from 3/6, 5/6 to 6/6. The results show three-body aggregation of guest molecules is a key step in gas hydrate formation.

On the size, shape and energetics of the hydration shell around alkanes

A large number of clathrate-like cages have been proposed as the very first hydration shell of alkanes. The cages include canonical structures commonly found in clathrate hydrates and many others, not previously reported, derived from the carbon fullerene cavities. These structures have a rich and variegated form, which can adapt to the shape and conformation of the solute. They avoid "wasting" hydrogen bonds, while minimizing the volume cage and maximizing the solute-solvent van der Waals interactions. DFT/M06-2X and MP2 ab initio calculations give comparable structural and energetic results although the latter predicts slightly larger cages for a given solute. It is shown that the van der Waals interactions are substantial and the large exoenergetic values found for isobutane and cyclopentane provide an explanation for the surprising high melting points of related hydrates at room pressure. The encaging enthalpy for various hydrocarbons is similar to the enthalpy of solution measured at a temperature just above the melting point of aqueous hydrocarbon solutions, thus indicating that water molecules should not deviate too much from the configuration with O-H bonds tangentially oriented with respect to the solute surface. The computed trend differs from the enthalpy of solution measured at room temperature, thus the very first hydration shell departs, up to a certain degree, from the clathrate-like structures.

Atomic Force Microscopy Study of Non-DLVO Interactions between Drops and Bubbles

The heterointeraction between liquid drops and air bubbles dispersed in another immiscible liquid is studied with the application of the atomic force microscopy (AFM) probe techniques. The tetradecane drops and air bubbles readily coalescence to form a lens-like structure in 100 mM sodium chloride aqueous solution, demonstrating strong hydrophobic (HB) attraction. The interaction range and strength of this hydrophobic attraction between oil drops and air bubbles is investigated by fine control of electrical double layer thicknesses related to specific electrolyte concentrations, and a midrange term in combination with a short-range term is found to present a proper characterization of this hydrophobic attraction. A further step is taken by introducing a triblock copolymer (Pluronic F68) into the aqueous solution, with results indicating that a relatively long-range steric hindrance (SH) furnished by a polymer "brush" surmounts the hydrophobic attraction. Finally, the interaction between a water drop and an air bubble in tetradecane is also measured as a comparison. The repelling action between a hydrophobic body (air bubble) and water drop indicates a strong repulsion. The present results show an interesting understanding of hydrophobic interactions between drops and bubbles, which is of potential application in controlling dispersion stability.

Role of water model on ion dissociation at ambient conditions

We study ion pair dissociation in water at ambient conditions using a combination of classical and ab initio approaches. The goal of this study is to disentangle the sources of discrepancy observed in computed potentials of mean force. In particular, we aim to understand why some models favor the stability of solvent-separated ion pairs vs contact ion pairs. We found that some observed differences can be explained by non-converged simulation parameters. However, we also unveil that for some models, small changes in the solution density can have significant effects on modifying the equilibrium balance between the two configurations. We conclude that the thermodynamic stability of contact and solvent-separated ion pairs is very sensitive to the dielectric properties of the underlying simulation model. In general, classical models are very robust in providing a similar estimation of the contact ion pair stability, while this is much more variable in density functional theory-based models. The barrier to transition from the solvent-separated to contact ion pair is fundamentally dependent on the balance between electrostatic potential energy and entropy. This reflects the importance of water intra- and inter-molecular polarizability in obtaining an accurate description of the screened ion-ion interactions.

Origin of the Surprising Mechanical Stability of Kinesin's Neck Coiled Coil

Kinesin-1 is a motor protein moving along a microtubule with its two identical motor heads dimerized by two neck linkers and a coiled-coil stalk. When both motor heads bind the microtubule, an internal strain is built up between the two heads, which is indispensable to ensure proper coordination of the two motor heads during kinesin-1's mechanochemical cycle. The internal strain forms a tensile force along the neck linker that tends to unwind the neck coiled coil (NCC). Experiments showed that the kinesin-1's NCC has a high antiunwinding ability compared with conventional coiled coils, which was mainly attributed to the enhanced hydrophobic pressure arising from the unconventional sequence of kinesin-1's NCC. However, hydrophobic pressure cannot provide the shearing force which is needed to balance the tensile force on the interface between two helices. To find out the true origin of the mechanical stability of kinesin-1's NCC, we perform a novel and detailed mechanical analysis for the system based on molecular dynamics simulation at an atomic level. We find that the needed shearing force is provided by a buckle structure formed by two tyrosines which form effective steric hindrance in the presence of tensile forces. The tensile force is balanced by the tensile direction component of the contact force between the two tyrosines which forms the shearing force. The hydrophobic pressure balances the other component of the contact force perpendicular to the tensile direction. The antiunwinding strength of NCC is defined by the maximum shearing force, which is finally determined by the hydrophobic pressure. Kinesin-1 uses residues with plane side chains, tryptophans and tyrosines, to form the hydrophobic center and to shorten the interhelix distance so that a high antiunwinding strength is obtained. The special design of NCC ensures exquisite cooperation of steric hindrance and hydrophobic pressure that results in the surprising mechanical stability of NCC.

Synergistic and Antagonistic Effects of Aromatics on the Agglomeration of Gas Hydrates

Surfactants are often used to stabilize aqueous dispersions. For example, surfactants can be used to prevent hydrate particles from forming large plugs that can clog, and sometimes rupture pipelines. Changes in oil composition, however dramatically affect the performance of said surfactants. In this work we demonstrate that aromatic compounds, dissolved in the hydrocarbon phase, can have both synergistic and antagonistic effects, depending on their molecular structure, with respect to surfactants developed to prevent hydrate agglomerations. While monocyclic aromatics such as benzene were found to disrupt the structure of surfactant films at low surfactant density, they are expelled from the interfacial film at high surfactant density. On the other hand, polycyclic aromatics, in particular pyrene, are found to induce order and stabilize the surfactant films both at low and high surfactant density. Based on our simulation results, polycyclic aromatics could behave as natural anti-agglomerants and enhance the performance of the specific surfactants considered here, while monocyclic aromatics could, in some cases, negatively affect performance. Although limited to the conditions chosen for the present simulations, the results, explained in terms of molecular features, could be valuable for better understanding synergistic and antagonistic effects relevant for stabilizing aqueous dispersions used in diverse applications, ranging from foodstuff to processing of nanomaterials and advanced manufacturing.

Shaft Function of Kinesin-1's α4 Helix in the Processive Movement

INTRODUCTION

Kinesin-1 motor is a molecular walking machine constructed with amino acids. The understanding of how those structural elements play their mechanical roles is the key to the understanding of kinesin-1 mechanism.

METHODS

Using molecular dynamics simulations, we investigate the role of a helix structure, α4 (also called switch-II helix), of kinesin-1's motor domain in its processive movement along microtubule.

RESULTS

Through the analysis of the structure and the interactions between α4 and the surrounding residues in different nucleotide-binding states, we find that, mechanically, this helix functions as a shaft for kinesin-1's motor-domain rotation and, structurally, it is an amphipathic helix ensuring its shaft functioning. The hydrophobic side of α4 consists strictly of hydrophobic residues, making it behave like a lubricated surface in contact with the core β-sheet of kinesin-1's motor domain. The opposite hydrophilic side of α4 leans firmly against microtubule with charged residues locating at both ends to facilitate its positioning onto the intra-tubulin groove.

CONCLUSIONS

The special structural feature of α4 makes for an effective reduction of the conformational work in kinesin-1's force generation process.

Does a pair of methane molecules aggregate in water?

Molecular dynamics (MD) simulations of methane-water mixtures were performed using ab initio force fields for the CH₄-H₂O, H₂O-H₂O, and CH₄-CH₄ interactions. Both methane and water molecules were polarizable. From these calculations, the potential of mean force (PMF) between two methane molecules was extracted. Our results are compared with PMFs from a density-functional-theory (DFT) based Born-Oppenheimer type MD (BOMD) simulation, from a Monte Carlo (MC) simulation with ab initio-based force fields, and from MD simulations with empirical force fields. Our PMF is qualitatively similar to that obtained from the simulations with empirical force fields but differs significantly from those resulting from the DFT-BOMD and MC simulations. The depth of the PMF global minimum obtained in the present work is in a much better agreement with the experimental estimate than the result of the DFT-BOMD simulation, possibly due to the inability of DFT to describe the dispersion interactions and the lack of extensive sampling in the BOMD simulations. Our work indicates that, for a pair of methane molecules, there are configurations where the solvent increases the attraction between the solutes, but there are also conformations in which the solvent causes a weak net repulsion. On average, the methane molecules are more likely to be in the configuration where they are separated by a water molecule than in the one in which they are in contact even though the minimum of the PMF at the latter configuration is deeper than that at the former. Finally, we found that the water structure around methane solutes does not show a greater tetrahedral ordering than in neat bulk water.

Prediction of the near-IR spectra of ices by ab initio molecular dynamics

Use of ab initio molecular dynamics to predict the near-IR spectra of ices and application to astronomical models.

Influence of glycerol on the cooling effect of pair hydrophobicity in water: relevance to proteins’ stabilization at low temperature

Glycerol reduces the cooling effect of pair hydrophobicity (reduction of hydrophobicity with decreasing temperature) in water.

Computational biochemical investigation of the binding energy interactions between an estrogen receptor and its agonists

We present the energy profiles of estrogen receptor–agonist ligand interactions in atomic detail using a quantum biochemical approach.

Energy benchmarks for methane-water systems from quantum Monte Carlo and second-order Møller-Plesset calculations

The quantum Monte Carlo (QMC) technique is used to generate accurate energy benchmarks for methane-water clusters containing a single methane monomer and up to 20 water monomers. The benchmarks for each type of cluster are computed for a set of geometries drawn from molecular dynamics simulations. The accuracy of QMC is expected to be comparable with that of coupled-cluster calculations, and this is confirmed by comparisons for the CH4-H2O dimer. The benchmarks are used to assess the accuracy of the second-order Møller-Plesset (MP2) approximation close to the complete basis-set limit. A recently developed embedded many-body technique is shown to give an efficient procedure for computing basis-set converged MP2 energies for the large clusters. It is found that MP2 values for the methane binding energies and the cohesive energies of the water clusters without methane are in close agreement with the QMC benchmarks, but the agreement is aided by partial cancelation between 2-body and beyond-2-body errors of MP2. The embedding approach allows MP2 to be applied without loss of accuracy to the methane hydrate crystal, and it is shown that the resulting methane binding energy and the cohesive energy of the water lattice agree almost exactly with recently reported QMC values.

Free energetics of carbon nanotube association in aqueous inorganic NaI salt solutions: Temperature effects using all-atom molecular dynamics simulations

In this study, we examine the temperature dependence of free energetics of nanotube association using graphical processing unit-enabled all-atom molecular dynamics simulations (FEN ZI) with two (10,10) single-walled carbon nanotubes in 3 m NaI aqueous salt solution. Results suggest that the free energy, enthalpy and entropy changes for the association process are all reduced at the high temperature, in agreement with previous investigations using other hydrophobes. Via the decomposition of free energy into individual components, we found that solvent contribution (including water, anion, and cation contributions) is correlated with the spatial distribution of the corresponding species and is influenced distinctly by the temperature. We studied the spatial distribution and the structure of the solvent in different regions: intertube, intratube and the bulk solvent. By calculating the fluctuation of coarse-grained tube-solvent surfaces, we found that tube-water interfacial fluctuation exhibits the strongest temperature dependence. By taking ions to be a solvent-like medium in the absence of water, tube-anion interfacial fluctuation shows similar but weaker dependence on temperature, while tube-cation interfacial fluctuation shows no dependence in general. These characteristics are discussed via the malleability of their corresponding solvation shells relative to the nanotube surface. Hydrogen bonding profiles and tetrahedrality of water arrangement are also computed to compare the structure of solvent in the solvent bulk and intertube region. The hydrophobic confinement induces a relatively lower concentration environment in the intertube region, therefore causing different intertube solvent structures which depend on the tube separation. This study is relevant in the continuing discourse on hydrophobic interactions (as they impact generally a broad class of phenomena in biology, biochemistry, and materials science and soft condensed matter research), and interpretations of hydrophobicity in terms of alternative but parallel signatures such as interfacial fluctuations, dewetting transitions, and enhanced fluctuation probabilities at interfaces.

Hydrophobic Ambivalence: Teetering on the Edge of Randomness

Processes ranging from oil-water phase separation to the formation of solid clathrate hydrates send mixed messages regarding whether oil molecules hate or love to be surrounded by water. Recent experimental and theoretical results help decipher these mixed messages by illuminating the conditions under which the stability of a hydrophobic contact is expected to exceed thermal energy fluctuations - thus facilitating hydrophobic self-assembly and the emergence of structure from randomness. Important open questions remain regarding the dependence of hydrophobic interactions on molecular size and temperature, as well as the balance of direct and water-mediated interactions.

Contacts Between Alcohols in Water Are Random Rather than Hydrophobic

Given the importance of water-mediated hydrophobic interactions in a wide range of biological and synthetic self-assembly processes, it is remarkable that both the sign and the magnitude of the hydrophobic interactions between simple amphiphiles, such as alcohols, remain unresolved. To address this question, we have performed Raman hydration-shell vibrational spectroscopy and polarization-resolved femtosecond infrared experiments, as well as random mixing and molecular dynamics simulations. Our results indicate that there are no more hydrophobic contacts in aqueous solutions of alcohols ranging from methanol to tertiary butyl alcohol than in random mixtures of the same concentration. This implies that the interaction between small hydrophobic groups is weaker than thermal energy fluctuations. Thus, the corresponding water-mediated hydrophobic interaction must be repulsive, with a magnitude sufficient to negate the attractive direct van der Waals interaction between the hydrophobic groups.

Ab initio Kinetic Monte Carlo simulations of dissolution at the NaCl-water interface

We have used ab initio molecular dynamics (AIMD) simulations to study the interaction of water with the NaCl surface. As expected, we find that water forms several ordered hydration layers, with the first hydration layer having water molecules aligned so that oxygen atoms are on average situated above Na sites. In an attempt to understand the dissolution of NaCl in water, we have then combined AIMD with constrained barrier searches, to calculate the dissolution energetics of Na(+) and Cl(-) ions from terraces, steps, corners and kinks of the (100) surface. We find that the barrier heights show a systematic reduction from the most stable flat terrace sites, through steps to the smallest barriers for corner and kink sites. Generally, the barriers for removal of Na(+) ions are slightly lower than for Cl(-) ions. Finally, we use our calculated barriers in a Kinetic Monte Carlo as a first order model of the dissolution process.

Electrostatic contribution from solvent in modulating single-walled carbon nanotube association

We perform all-atom molecular dynamics simulations to compute the potential of mean force (PMF) between two (10,10) single-walled carbon nanotubes solvated in pure nonpolarizable SPC/E and polarizable TIP4P-FQ water, at various temperatures. In general, the reversible work required to bring two nanotubes from a dissociated state (free energy reference) to contact state (free energy minimum) is more favorable and less temperature-dependent in TIP4P-FQ than in SPC/E water models. In contrast, molecular properties and behavior of water such as the spatially-resolved water number density (intertube, intratube, or outer regions), for TIP4P-FQ are more sensitive to temperature than SPC/E. Decomposition of the solvent-induced PMF into different spatial regions suggests that TIP4P-FQ has stronger temperature dependence; the opposing destabilizing/stabilizing contributions from intertube water and more distal water balance each other and suppress the temperature dependence of total association free energy. Further investigation of hydrogen bonding network in intertube water reveals that TIP4P-FQ retains fewer hydrogen bonds than SPC/E, which correlates with the lower water number density in this region. This reduction of hydrogen bonds affects the intertube water dipoles. As the intertube volume decreases, TIP4P-FQ dipole moment approaches the gas phase value; the distribution of dipole magnitude also becomes narrower due to less average polarization/perturbation from other water molecules. Our results imply that the reduction of water under confinement may seem trivial, but underlying effects to structure and free energetics are non-negligible.

Characterization of the Methane-Graphene Hydrophobic Interaction in Aqueous Solution from Ab Initio Simulations

In this article, the interaction between a methane molecule and a graphene plane in liquid water has been characterized employing DFT-based free energy Molecular Dynamics calculations. This system represents a good model to understand the generic interaction between a small hydrophobic solute (methane molecule) and an extense hydrophobic surface (graphene plane). The structural and dynamical properties of graphene and methane hydration water are analyzed and found to be closely related to the main features of the potential of mean force. The results could be used in coarse-grained models to take into account the effect of the hydrophobic interaction in realistic systems relevant to experiment.

Soft collective fluctuations governing hydrophobic association

The interaction between two associating hydrophobic particles has traditionally been explained in terms of the release of entropically frustrated hydration shell water molecules. However, this picture cannot account for the kinetics of hydrophobic association and is therefore not capable of providing a microscopic description of the hydrophobic interaction (HI). Here, Monte Carlo simulations of a pair of molecular-scale apolar solutes in aqueous solution reveal the critical role of collective fluctuations in the hydrogen bond (HB) network for the microscopic picture of the HI. The main contribution to the HI is the relaxation of solute-water translational correlations. The existence of a heat capacity maximum at the desolvation barrier is shown to arise from softening of non-HB water fluctuations and the relaxation of many-body correlations in the labile HB network. The microscopic event governing the kinetics of hydrophobic association has turned out to be a relatively large critical collective fluctuation in hydration water displacing a substantial fraction of HB clusters from the inner to the outer region of the first hydration shell.