Water diffusion in carbon nanotubes: Interplay between confinement, surface deformation, and temperature.

In this article, we investigate, through molecular dynamics simulations, the diffusion behavior of the TIP4P/2005 water confined in pristine and deformed carbon nanotubes (armchair and zigzag). To analyze different diffusive mechanisms, the water temperature was varied as 210 ≤ T ≤ 380 K. The results of our simulations reveal that water presents a non-Arrhenius to Arrhenius diffusion crossover. The confinement shifts the diffusion transition to higher temperatures when compared with the bulk system. In addition, for narrower nanotubes, water diffuses in a single line, which leads to its mobility independent of the activation energy.

Water diffusion in rough carbon nanotubes.

We use molecular dynamics simulations to study the diffusion of water inside deformed carbon nanotubes with different degrees of deformation at 300 K. We found that the number of hydrogen bonds that water forms depends on nanotube topology, leading to enhancement or suppression of water diffusion. The simulation results reveal that more realistic nanotubes should be considered to understand the confined water diffusion behavior, at least for the narrowest nanotubes, when the interaction between water molecules and carbon atoms is relevant.

Relation between boundary slip mechanisms and waterlike fluid behavior.

The slip of a fluid layer in contact with a solid confining surface is investigated for different temperatures and densities using molecular dynamic simulations. We show that for an anomalous waterlike fluid the slip goes as follows: for low levels of shear, defect slip appears and is related to the particle exchange between the fluid layers; at high levels of shear, global slip occurs and is related to the homogeneous distribution of the fluid in the confining surfaces. The oscillations in the transition velocity from defect to global slip are shown to be associated with changes in the layering distribution in the anomalous fluid.

Single-file mobility of water-like fluid in a generalized Frenkel-Kontorova model.

In this work, we used a generalized Frenkel-Kontorova model to study the mobility of water molecules inside carbon nanotubes with small radius at low temperatures. Our simulations show that the mobility of confined water decreases monotonically increasing the amplitude of the substrate potential at fixed commensurations. On the other hand, the mobility of the water molecules shows a non-monotonic behavior when varying the commensuration. This result indicates that the mobility of the confined fluid presents different behavior regimes depending on the amplitude of the water-nanotube interaction. In order to qualitatively understand these results, we study analytically the driven Frenkel-Kontorova model at finite temperatures. This analysis allows us to obtain the curves of the mobility versus commensurations, at fixed substrate potentials. Such curves show the existence of three regimes of mobility behavior as a function of the commensuration ratio. Additionally, our study indicates a nontrivial and strong dependence of the mobility with a quantity that can be interpreted as an effective amplitude of the substrate potential, depending on the bare amplitude of the substrate potential, the commensuration ratio, and temperature.

Relation between occupation in the first coordination shells and Widom line in core-softened potentials.

Three core-softened families of potentials are checked for the presence of density and diffusion anomalies. These potentials exhibit a repulsive core with a softening region and at larger distances an attractive well. We found that the region in the pressure-temperature phase diagram in which the anomalies are present increases if the slope between the core-softened scale and the attractive part of the potential decreases. The anomalous region also increases if the range of the core-softened or of the attractive part of the potential decreases. We also show that the presence of the density anomaly is consistent with the non-monotonic changes of the radial distribution function at each one of the two scales when temperature and density are varied. Then, using this anomalous behavior of the structure we show that the pressure and the temperature at which the radial distribution function of one of the two length scales equals the radial distribution function of the other length scales identify the Widom line.

Thermodynamic, dynamic, structural, and excess entropy anomalies for core-softened potentials.

Using molecular dynamic simulations, we study three families of continuous core-softened potentials consisting of two length scales: a shoulder scale and an attractive scale. All the families have the same slope between the two length scales but exhibit different potential energy gap between them. For each family three shoulder depths are analyzed. We show that all these systems exhibit a liquid-liquid phase transition between a high density liquid phase and a low density liquid phase ending at a critical point. The critical temperature is the same for all cases suggesting that the critical temperature is only dependent on the slope between the two scales. The critical pressure decreases with the decrease of the potential energy gap between the two scales suggesting that the pressure is responsible for forming the high density liquid. We also show, using the radial distribution function and the excess entropy analysis, that the density, the diffusion, and the structural anomalies are present if particles move from the attractive scale to the shoulder scale with the increase of the temperature indicating that the anomalous behavior depends only in what happens up to the second coordination shell.

Core-softened fluids, water-like anomalies, and the liquid-liquid critical points.

Molecular dynamics simulations are used to examine the relationship between water-like anomalies and the liquid-liquid critical point in a family of model fluids with multi-Gaussian, core-softened pair interactions. The core-softened pair interactions have two length scales, such that the longer length scale associated with a shallow, attractive well is kept constant while the shorter length scale associated with the repulsive shoulder is varied from an inflection point to a minimum of progressively increasing depth. The maximum depth of the shoulder well is chosen so that the resulting potential reproduces the oxygen-oxygen radial distribution function of the ST4 model of water. As the shoulder well depth increases, the pressure required to form the high density liquid decreases and the temperature up to which the high-density liquid is stable increases, resulting in the shift of the liquid-liquid critical point to much lower pressures and higher temperatures. To understand the entropic effects associated with the changes in the interaction potential, the pair correlation entropy is computed to show that the excess entropy anomaly diminishes when the shoulder well depth increases. Excess entropy scaling of diffusivity in this class of fluids is demonstrated, showing that decreasing strength of the excess entropy anomaly with increasing shoulder depth results in the progressive loss of water-like thermodynamic, structural and transport anomalies. Instantaneous normal mode analysis was used to index the overall curvature distribution of the fluid and the fraction of imaginary frequency modes was shown to correlate well with the anomalous behavior of the diffusivity and the pair correlation entropy. The results suggest in the case of core-softened potentials, in addition to the presence of two length scales, energetic, and entropic effects associated with local minima and curvatures of the pair interaction play an important role in determining the presence of water-like anomalies and the liquid-liquid phase transition.

Entropy, diffusivity and the energy landscape of a waterlike fluid.

Molecular dynamics simulations and instantaneous normal mode (INM) analysis of a fluid with core-softened pair interactions and waterlike liquid-state anomalies are performed to obtain an understanding of the relationship between thermodynamics, transport properties, and the potential energy landscape. Rosenfeld scaling of diffusivities with the thermodynamic excess and pair correlation entropy is demonstrated for this model. The INM spectra are shown to carry information about the dynamical consequences of the interplay between length scales characteristic of anomalous fluids, such as bimodality of the real and imaginary branches of the frequency distribution. The INM spectral information is used to partition the liquid entropy into two contributions associated with the real and imaginary frequency modes; only the entropy contribution from the imaginary branch captures the nonmonotonic behavior of the excess entropy and diffusivity in the anomalous regime of the fluid.

Effects of the attractive interactions in the thermodynamic, dynamic, and structural anomalies of a two length scale potential.

Using molecular dynamic simulations, we study a system of particles interacting through a continuous core-softened potentials consisting of a hard core, a shoulder at closest distances, and an attractive well at further distance. We obtain the pressure-temperature phase diagram of this system for various depths of the tunable attractive well. Since this is a two length scale potential, density, diffusion, and structural anomalies are expected. We show that the effect of increasing the attractive interaction between the molecules is to shrink the region in pressure in which the density and the diffusion anomalies are present. If the attractive forces are too strong, particle will be predominantly in one of the two length scales and no density of diffusion anomaly is observed. The structural anomalous region is present for all the cases.

Thermodynamic, dynamic, and structural anomalies for shoulderlike potentials.

Using molecular dynamic simulations we study a family of continuous core-softened potentials consisting of a hard core, a shoulder at closest distances, and an attractive well at further distance. The repulsive shoulder and the well distances represent two length scales. We show that if the first scale, the shoulder, is repulsive or has a small well, the potential has a region in the pressure-temperature phase diagram with density, diffusion, and structural anomalies. However, if the closest scale becomes a deep well, the regions in the pressure-temperature phase diagram where the three anomalies are present shrink and disappear. This result helps in defining two length scales potentials that exhibit anomalies.