arxiv: v1 [cond-mat.soft] 26 Nov 2018

Transcription

1 Water diffusion in rough carbon nanotubes Bruno Henrique da Silva e Mendonça a,, Patricia Ternes a, Evy Salcedo b, Alan B. de Oliveira c, Marcia C. Barbosa a a Instituto de Física, Universidade Federal do Rio Grande do Sul, Porto Alegre, RS , Brazil b Coordenadoria Especial de Física, Química e Matemática, Universidade Federal de Santa Catarina, Rua Pedro João Pereira, 150, , Araranguá, SC, Brazil c Departamento de Física, Universidade Federal de Ouro Preto, Ouro Preto, MG, , Brazil arxiv: v1 [cond-mat.soft] 26 Nov 2018 Abstract Keywords: 1. Introduction Confined Water; Mobility; Carbon Nanotube; Diffusion. Some physical and chemical properties observed in nanomaterials differ from the properties presented by the same material in micro or macroscopic scale.[1] When the material is hollow, it becomes crucial to know the behavior of this material in the presence of a solvent. In this context, the behavior of water confined in carbon nanotubes (CNTs) has been the subject of a large number of works [2, 3, 4, 5, 6]. It is widely established that the confined water in carbon nanotubes presents a different behavior from that observed in bulk systems, including results for diffusion coefficient.[7, 8] The bulk water diffusion constant, at 298 K and 1 atm determined experimentally is cm 2 s 1, [3] while the bulk diffusion constant, determined through computational study, vary from 1.4 [9] up to 5.06(9), [10] in units of 10 5 cm 2 s 1, depending on the model of liquid water.[11, 12] For water confined in carbon nanotubes, the diffusion coefficient vary with the nanotube radius. Farimani et al. (2011) [13] studied the average axial diffusion constant for SPC/E water, at 300 K, in different size of armchair nanotubes. They showed that for diameters, d, above 6 nm the diffusion is close to the bulk value. As the diameter decreases the diffusion constant increases until reaching a maximum value for d = 2.7 nm. For diameters from 0.95 up to 2.2 nm the diffusion constant is less than the bulk value and presents a minimum for d = 1.22 nm. This minimum occurs due to the ice-like structure assumed by water in the (9,9) CNT, at 300 K, leading an anomalous diffusion behavior.[13] Several works reports on diffusion of water confined in CNT. Different water models, carbon nanotubes, carbon-oxygen interaction, and different methodologies were applied. [14, 15, 16] Generally, these studies are performed with perfect nanotubes. However, in nature, Corresponding author. address: brunnohennrique13@gmail.com (Bruno Henrique da Silva e Mendonça) Preprint submitted to Physica A November 27, 2018

2 the probability of occurrence of perfect CNTs is very small since the nanotubes are subject to compression, and molecules can be adsorbed and deposited in the nanotubes. Recently, Mendonça et al. reported on the study of TIP4P/2005 water diffusion confined in kneaded nanotubes. They systematically analyzed the effect of the degree of eccentricity of the nanotube on the diffusion constant, for different armchair CNTs radius. They found that large deformations (eccentricity) favor the formation of hydrogen bonds, which often drastically reduces diffusion, but for the (9,9) CNT the behavior is quite different. The high eccentricity promotes the breakdown of hydrogen bonds, that melts the anomalous ice-like structure observed in perfect (9,9) nanotube, leading to a non-null diffusion. Although the system studied by Mendonça et al. presents a radial deformation, it preserves an axial symmetry, which is not observed experimentally.[17, 18] In this work we calculated the TIP4P/2005 water diffusion constant confined in carbon nanotubes with non uniform wrinkles, creating a system closer to reality in the diffusive process of water. The effect of pore size on water diffusion was analyzed. In addition, to understanding the effect of nanotube topology to the diffusion behavior, different nanotube helicity, and different carbon-oxygen interaction, were considered. 2. The model, simulation and methods 2.1. Nanotube models We study the water behavior confined in different carbon nanotubes. Each perfect nanotube was submitted to two different deformation processes. First the entire nanotube was uniformly kneaded in the y direction until the nanotube reaches an elliptical cylinder shape with eccentricity of 0.8. The second process was to randomly kneaded the nanotube in the y direction until the nanotube has non-uniform wrinkles. The perfect nanotube will be identify by P (n, m), the uniformly kneaded nanotube will be identify by K(n, m), and the randomly wrinkled nanotube will be identify by W (n, m), were (n, m) is related with the nanotube helicity and the nanotube diameter, d. Figure 1 show the perfect, the uniformly kneaded and the randomly wrinkled nanotubes for the armchair (9,9) case. Figure 1: Snapshot of typical nanotubes used. 2

3 Two armchair nanotubes (n = m) and two zigzag nanotubes (m = 0) were used. The details of each nanotube are listed in the table 1. It s interesting to note that the CNT(9,9) has almost the same diameter that the CNT(16,0), and the CNT(12,12) has almost the same diameter that the CNT(21,0). Table 1: Parameters for the simulated perfect carbon nanotubes. Density ρ A for the usual and hydrofilic nanotubes and the density ρ B of hydrophobic nanotubes. CNT d (nm) L z (nm) Number of H 2 O ρ A (g/cm 3 ) ρ B (g/cm 3 ) (9,9) (12,12) (16,0) (21,0) Simulation details and methods Molecular dynamics (MD) simulations at constant number of particles, volume and temperature are employed to analyzed the diffusion of water when confined in carbon nanotubes with different radius, helicity and defects. The water model used was the TIP4P/2005,[12] that presents a good agreement with many experimental results including results for diffusion coefficient.[19] The carbon-water interaction was modeled through the Lennard-Jones potential [20] between the carbon and oxygens atoms. Several different values for the carbonoxygen energy, ɛ CO, and effective carbon-oxygen diameter, σ CO, were found in previous works. Here we test the effect of three different set of LJ parameters: ɛ CO = kcal/mol and σ CO = 3.28 Å (set to usual), ɛ CO = kcal/mol and σ CO = 3.28 Å (set to hydrophilic) and ɛ CO = kcal/mol and σ CO = 3.28 Å (set to hydrophobic).[21, 22] Periodic boundary conditions in the z direction and a cutoff radius of 12 Å were used. The density was determined considering the volume excluded due the LJ interaction between carbon and oxygen atoms. Therefore the effective density will be ρ = 4M/[π(d σ CO ) 2 L z ], where M is the total water mass into the pore, and L z is the nanotube length.[23] The simulations were carried out with LAMMPS MD package.[24] The PPPM method was used to compute the long-range Coulomb interactions and the structure of water molecules were constraint through the SHAKE algorithm with a precision tolerance of The temperature was set at 300 K and controlled by the Nosé-Hoover thermostat with a damping time of 100 fs and a time step of 1 fs. The system was equilibrated during 5 ns, then the properties were stored every 0.01 ns during 5 ns, totaling simulations of 10 ns. The dynamics of confined water is characterized by calculation of the diffusion coefficient. The diffusion is minimal in the radial direction, therefore only the axial diffusion is considered, and is determined through the one dimensional Einstein relation 1 d D z = lim z 2 (τ), (1) τ 2 dτ 3

4 where z 2 (τ) = [z(τ 0 τ) z(τ 0 )] 2 is the mean square displacement (MSD) in the axial direction. To characterize the water structure we calculated the number of hydrogen bonds and several colormaps in the xy, xz and yz directions. The hydrogen bonds were calculated considering geometrical definition, that is α 30 r OO 3.50 Å, where α is the OH O angle and r OO is the distance between two oxygens. The xy colormaps are obtained dividing the xy plane in bins with length of 0.1 Å and cont the number of oxygen in each bin. The xz and yz colormaps are obtained following the same method. 3. Results and discussion We analyzed the (9,9), (12,12), (16,0), and (21,0) nanotubes under different degrees of deformation. Some set of LJ parameters for the carbon-oxygen interaction were used, but no significant difference was observed. In the next results only the usual set of parameters is showed. Figure 2 (a) shows the axial diffusion coefficient with the deformation of CNT. The results for the armchair perfect and kneaded nanotubes are in accordance with the previously reported.[13] Water confined on the P (12, 12) CNT present a diffusion higher than the water confined in the P (9, 9) CNT, since in the P (9, 9) nanotube the water forms an ice-like structure. This structure is due to the compatibility of the nanotube radius with the distance of the water-water hydrogen bond, favoring the formation of an ice-like network.[13] (a) (b) Figure 2: Axial diffusion constant (a), and average number of hydrogen bond per water molecule (b), as a function of deformation degree for different nanotubes radius and helicity. 4

5 The diffusion in the K(12, 12) is less than the value in the P (12, 12). The decrease of diffusion in the kneaded nanotubes when compared with the respective perfect nanotube is expected for most armchair nanotubes, with the exception of K(9, 9) CNT, where the radial deformation melts the ice-like structure, leading a breakdown of hydrogen bonds, increasing the diffusion. The average number of hydrogen bond showed in figure 2 (b) confirm that the kneaded (12,12) nanotube present a higher number of HB that the perfect (12,12) nanotube, while the K(9, 9) present a lower number of hydrogen bond when compared with P (9, 9). The water confined in the zigzag nanotubes present a diffusion greater than water confined in armchair nanotubes for almost all CNTs analyzed, with exception for the K(16, 0) that present a diffusion value close to the value for K(9, 9) and K(12, 12). The presence of radial defects induces an increase in the hydrogen bond number, leading a decrease of diffusion constante for the uniform kneaded zigzag nanotube. Generally, the presence of defects in the axial and radial direction causes the diffusion to be smaller in relation to the perfect nanotube, and greater in relation to the CNT with defect only in the radial direction. This is not true only in the (9, 9) CNT where an anomalous diffusion behavior is observed in the perfect case. When both radial and axial defects are taken into account, the diffusion is not affected by pore size, but depends on the topology of nanotubes. The wrinkled zigzag nanotubes have a greater diffusion than the wrinkled armchair CNTs. Although the diffusion value is practically the same for W (16, 0) and W (21, 0), the average number of hydrogen bonds is different in each situation. The same behavior is observed in the armchair case. Water confined in W (9, 9) and W (12, 12) CNTs presents the same diffusion, but different number of HB. These results indicate that in the presence of radial and axial defects, the diffusion mechanism depends not only on the amount of hydrogen bonds that the water can form. For a better understanding of the diffusion mechanism under different deformation degrees, we analyzed the structure assumed by the confined water. Figure 3 shows some informations about the structure assumed by confined water in the perfect nanotubes. The colormaps (A), (B), (C), and (D), represent a radial distribution of water molecules over all nanotube length for P (9, 9), P (16, 0), P (21, 0), and P (12, 12), respectively. The figure 3 (E) shows the excluded distance as function of nanotube radius. 5

6 Figure 3: (Color online) (A) Colormaps of xy direction for P (9, 9), P (16, 0) (B), P (21, 0) (C), and P (12, 12) (D). Dark blue regions have low probability to find water molecules, while red regions have high probability to find water molecules. (E) Excluded distance as a function of perfect nanotube radius. Although the nanotubes (9,9) and (16,0) have practically the same radius, R, the difference of nanotube topology promotes more or less carbon-water repulsion. For the armchair (9,9) the carbon-water repulsion is greater, and only one layer of water is formed, and the anomalous diffusion is observed. For the P (21, 0), and P (12, 12) (figures 3 (C) and (D)) the water presented two layers. Although the number of layers is equal, the diffusion for the case P (21.0) is almost double the value observed for the case P (12, 12). This increased diffusion is due to a greater carbon-water repulsion for P (21, 0) than for P (12, 12). A similar analyzes were made to the kneaded nanotubes. The colormaps (A), (B), (C), and (D), represents the radial distribution of water molecules over all nanotube length for K(9, 9), K(16, 0), K(21, 0), and K(12, 12), respectively. The figure 4 (E) shows the excluded distance as function of nanotube radius only over the semi-minor axis as indicate in scheme of Fig. 4 (E). 6

7 Figure 4: (Color online) (A) Colormaps of xy direction for K(9, 9), K(16, 0) (B), K(21, 0) (C), and K(12, 12) (D). Dark blue regions have low probability to find water molecules, while red regions have high probability to find water molecules. (E) Excluded distance as a function of nanotube semi-minor axis. For the kneaded nanotubes, the major excluded radius occurs for the K(21, 0) CNT, that is the system that present the higher diffusion value. The water in kneaded (16,0) nanotube form single-file structures and have the minor excluded radius, leading to greater diffusion decrease when compared to its respective perfect nanotube. The radius excluded depends on nanotube topology and pore size, that is, by varying these two parameters, the water-water interaction becomes more, or less, relevant in relation to the water-wall interaction. For the wrinkled nanotube is not possible to make a similar analyzes of the excluded radius, because there is no symmetry along the axial direction. Therefore, we analyzed several xy colormaps of different nanotube regions, with length of 25 nm each, in order to extract information from the radial structuring. Additionally, xz and yz colormaps are analyzed. Figure 5 shows these colormaps for the W (9, 9) case. 7

8 Figure 5: Colormaps of xz, yz and xy directions for the (9, 9) randomly wrinkled nanotube. The xy colormaps are dependent of the nanotube region, so each xy colormap is related to one nanotube region with 25nm in z direction. Dark blue regions have low probability to find water molecules, while red regions have high probability to find water molecules. The xy maps show several different radial structures. To water transition from one radial structure to another, the radial hydrogen bonds are broken, and the transition between regions with different degrees of defects is abrupt. The xz and yz colormaps show that water 8

9 tends to maintain the axial hydrogen bonds, performing a smooth transition of structure between regions with different degrees of defects. The same behavior was found for W (16, 0) nanotube. Figure 6 shows the colormaps for the W (21, 0) case. Again, the xy colormaps show several distinct radial structures, but different from that observed in figure 5, the difference between the structures is small. The transition between the radial structures is smooth and the water is able to make radial hydrogen bonds. The xz and yz colormaps confirm that water tends to maintain the axial hydrogen bonds. The same behavior was found for W (12, 12) nanotube. 9

10 Figure 6: Colormaps of xz, yz and xy directions for the (21, 0) randomly wrinkled nanotube. The xy colormaps are dependent of the nanotube region, so each xy colormap is related to one nanotube region with 25nm in z direction. Dark blue regions have low probability to find water molecules, while red regions have high probability to find water molecules Theses results explain the values obtained for the hydrogen bonds shown in figure 2 (b). Axially, each water molecule makes two hydrogen bonds. For nanotubes with smaller pore sizes, where the water has many radial transitions, the average number of hydrogen bonds is close to 2. 10

11 4. Conclusions In this work we analyzed the water diffusion constant under confinement with different degrees of defects. Different pore sizes and diferentes pore topology were considered. In general, perfect nanotubes lead to a larger diffusion than defective nanotubes, with the exception of P (9, 9) CNT, where the water has an anomalous diffusion. Therefore, defects in the nanotube suppress the anomalous diffusion enhancement of water inside CNTs. For all CNTs analyzed, the diffusion for the nanotubes with defects only in the radial direction was smaller than the diffusion for the nanotubes with defects in the radial and axial directions. Analyzing the structure along these defective nanotubes, we showed that the axial defect promotes to a breakdown of hydrogen bond in the radial direction, leading to an increased diffusion, if we compared with the kneaded results. The diffusion observed in the wrinkled nanotubes presented a strong relation with the topology of the nanotube. The diffusion is larger for zigzag CNTs and smaller for armchair CNTs. A very interesting result is that the combination of axial and radial defects causes the diffusion to be independent of the nanotube radius. The value of the diffusion constant depends in a non trivial way on a balance between the structure assumed by the water, the water-water interaction, and the water-carbon interaction. Acknowledgments This work is partially supported by Brazilian agencies CNPq and CAPES, Universidade Federal de Ouro Preto and Universidade Federal do Rio Grande do Sul. References References [1] M. A. Ratner, D. Ratner, Prentice Hall Professional 1 (2003) 133. [2] A. Alexiadis, S. Kassinos, Chemical Reviews 108 (2008) [3] K. R. Harris, L. A. Woolf, Journal of the Chemical Society, Faraday Transactions 1: Physical Chemistry in Condensed Phases 76 (1980) [4] J. Wang, Y. Zhu, J. Zhou, X.-H. Lu, Physical Chemistry Chemical Physics 6 (2004) [5] Y.-C. Liu, J.-W. Shen, K. E. Gubbins, J. D. Moore, T. Wu, Q. Wang, Physical Review B 77 (2008) [6] Y.-g. Zheng, H.-f. Ye, Z.-q. Zhang, H.-w. Zhang, Physical Chemistry Chemical Physics 14 (2012) [7] A. I. Kolesnikov, J.-M. Zanotti, C.-K. Loong, P. Thiyagarajan, A. P. Moravsky, R. O. Loutfy, C. J. Burnham, Physical Review Letters 93 (2004) [8] N. Naguib, H. Ye, Y. Gogotsi, A. G. Yazicioglu, C. M. Megaridis, M. Yoshimura, Nano Letters 4 (2004) [9] G. S. Fanourgakis, S. S. Xantheas, The Journal of Physical Chemistry A 110 (2006) [10] W. L. Jorgensen, J. Chandrasekhar, J. D. Madura, R. W. Impey, M. L. Klein, The Journal of chemical physics 79 (1983) [11] M. W. Mahoney, W. L. Jorgensen, The Journal of Chemical Physics 114 (2001) [12] J. L. F. Abascal, C. Vega, The Journal of Chemical Physics 123 (2005) [13] A. Barati Farimani, N. Aluru, The Journal of Physical Chemistry B 115 (2011)

12 [14] B. Xu, Y. Li, T. Park, X. Chen, The Journal of chemical physics 135 (2011) [15] G. Arora, N. J. Wagner, S. I. Sandler, Langmuir 20 (2004) [16] H. Ye, Y. Zheng, Z. Zhang, Z. Chen, H. Zhang, Scientific reports 5 (2015) [17] S. Iijima, nature 354 (1991) 56. [18] T. I. S. Iijima, Nature 363 (1993) 603. [19] C. Vega, J. L. F. Abascal, Physical Chemistry Chemical Physics 13 (2011) [20] J. E. Lennard-Jones, Proceedings of the Physical Society 43 (1931) 461. [21] G. Hummer, J. C. Rasaiah, J. P. Noworyta, Nature 414 (2001) 188. [22] I. Moskowitz, M. A. Snyder, J. Mittal, The Journal of chemical physics 141 (2014) 18C532. [23] J. R. Bordin, A. B. de Oliveira, A. Diehl, M. C. Barbosa, The Journal of chemical physics 137 (2012) [24] S. Plimpton, J. Comput. Phys. 117 (1995) 1, 12