Water molecules in a carbon nanotube under an applied electric field at various temperatures and pressures

Water confined in carbon nanotubes (CNTs) under the influence of an electric field exhibits behavior different to that of bulk water. Such behavior is fascinating from a nanoscience point of view and has potential application in nanotechnology. Using molecular dynamics simulations, we investigate the structure of water molecules in an ( 8 , 8 ) CNT, under an electric field at various temperatures and pressures. In the absence of an electric field, water in the CNT has an ordered (solid-like) structure at temperatures of 200 K and 250 K. The solid-like structure of water at these low temperatures exhibits ferroelectric properties. At 300 K, the structure of water is solid-like or disordered (liquid-like), i.e., an unstable structure. This indicates that a melting point occurs at around these conditions. Increasing the pressure to 10 MPa does not change the structure at 300 K. At 350 K, water is completely melted and has only a disordered structure. Under an applied electric field of 1 V/nm, water forms a solid-like structure at all simulation temperatures up to 350 K. This suggests that the electric field induces a phase transition from liquid to ice-nanotube, at temperatures as high as 350 K. The structure of the ice-nanotube under an applied electric field differs from that formed in the absence of an electric field at low temperature. The electrostatic interaction within the ice-nanotube under an electric field is stronger than that in the absence of an electric field.


Introduction
The ability of carbon nanotubes (CNTs) to transport water very fast [1][2][3][4] has led to huge recent interest in CNT-based nanofluidics. However, many important issues still need to be clarified. From a hydrodynamics aspect, the mechanism of the enhanced fluid flow through CNTs remains unclear [5,6]. Water exhibits different physical properties when confined in nanoscale space compared with those in the bulk, and this has not yet been explored in detail [7][8][9][10][11]. The ability of CNTs to deliver fluid allows them to be applied in nanofluidic devices such as separation membranes [12][13][14][15] and nano pumps [16][17][18].
Water molecules are polar, so are easily influenced by an electric field. The combined effects of nanoscale confinement and an applied electric field induce interesting properties in water [19][20][21][22][23][24][25]. An applied electric field strengthens the hydrogen bonding of water molecules between the two plates, and induces electrofreezing when the field direction is parallel to the confining wall [19,25]. In contrast, an applied electric field perpendicular to the confining wall disrupts hydrogen bonding, causing melting of the ice [22,23]. Water in nanotubes under the influence of an electric field is of interest from a nanoscience viewpoint and is also important in biophysics. Different ion concentrations across a cell membrane can induce a strong electric field in the cell channel [26][27][28]. Studying water in CNTs under an electric field could further our understanding of the behavior of water in the cell channel [29]. The effect of an electric field on water in CNTs is also relevant for nanotechnology-based applications such as nano pumping [18,[30][31][32] and desalination [33].
The dipole moments of water molecules align parallel to the direction of an applied electric field. The change in the molecular orientation of water within a CNT then causes its anomalous behavior. The applied electric field shifts the equilibrium state toward the condition in which water molecules favor filling a (6, 6) CNT [34]. A similar result has been demonstrated in larger CNTs, where the occupancy of water molecules in the CNTs increases with an electric field, implying a preference of water to fill the CNTs [35]. The preference of water to fill CNTs under an electric field can also result in a separation effect on water-alcohol solution [36,37]. The uniform molecular orientation under an electric field forms ordered structures in the CNTs [24,35,38,39]. These structures can be helical or otherwise, depending on the CNT diameter and electric field strength [35]. The electric field induces liquid-to-solid phase transition in CNTs under ambient conditions [24,35,38,39].
In the current study, we perform molecular dynamics (MD) simulations to investigate the structure of water molecules in an (8,8) CNT, in the presence and absence of an electric field at various temperatures and pressures. To mimic a membrane, the simulation consists of a CNT, graphene sheets, and reservoirs at both ends of the CNT. This configuration differs from that used in previous studies on finite length [38,39] and infinite length [24] CNTs. An applied electric field of 1 V/nm induces the formation of an ordered structure (solid-like structure) in the CNT, at all simulation temperatures and pressures. The applied electric field causes an ice-like structure to form at temperature as high as 350 K and pressure of 0.1 MPa. In the absence of an applied electric field and at temperatures of 200 K and 250 K, water in the CNT also forms an ordered structure. The structure in the absence of an electric field differs from that in the presence of the electric field. The electrostatic interaction within the ordered structure under an electric field is stronger than that in the absence of an electric field.

Method
The configuration for the MD simulation consisted of an (8,8) CNT embedded between two graphene sheets. Water reservoirs were located on the sides of sheets opposite the CNT (Figure 1). The configuration was similar to that used in our previous study [35,36]. The length of the CNT was 2.95 nm and the diameter was 1.08 nm. As the initial condition, each reservoir was filled with 4000 water molecules. All components were placed in a rectangular box that was subjected to periodic boundary conditions in all directions (x, y, and z-axes). The simulation box was 3.83 nm in width and 3.69 nm in height. After equilibration, the length of the simulation box was about 20 nm.
The SPC (simple point charge) model was used for the water molecule [40]. The Lennard-Jones (LJ) parameters for the oxygen atom of water were σ O = 0.3166 nm and O = 0.65 kJ/mol. The corresponding parameters for the carbon atom of the CNT were σ C = 0.34 nm and C = 0.3612 kJ/mol. The LJ parameters for the interactions of different atoms were calculated using the combination rule σ OC = (σ O σ C ) 1 2 and OC = ( O C ) 1 2 . A homogeneous 1 V/nm electric field (E) was applied to the axial direction of the CNT, which is in the positive direction of the z-axis. We also performed simulations with the TIP4P water model to clarify dependence of simulation results on water models. The simulation conditions were the same as those carried out with the SPC model.
The MD simulation was carried out using the Gromacs 4.5.5 software package (University of Groningen, Netherlands) [41]. The van der Waals interactions were cut off at 1.5 nm. Electrostatics were treated using the particle mesh Ewald (PME) method [42], with a real space cut off set to 1.5 nm. The SHAKE algorithm [43] was employed to fix the chemical bonds between the atoms of a water molecule. The CNT and graphene sheets were made rigid by maintaining the length and angle of chemical bonds as constant. The simulation was performed with NL x L y P z T, where the constant pressure (P) in the z-direction was maintained using the Parrinello-Rahman technique [44]. The temperature (T) was kept constant using the Nosé-Hoover coupling scheme, with a time constant of 0.4 ps [45,46]. The simulations were carried out at various temperatures at constant pressure, i.e., 200 K, 250 K, 300 K, and 350 K at 0.1 MPa, and at various pressures at constant temperature, i.e., 0.1 MPa, 1 MPa, and 10 MPa at 300 K. The time step was set to 2 fs, and the simulation time was 25 ns. The simulation reached equilibrium within 2 ns, and analysis data was taken after 3 ns.  (8,8) carbon nanotube (CNT), two graphene sheets, and water reservoirs on both sides. Water molecules fill the reservoirs, and an electric field E is applied in the z direction.

Water Structure
An electric field produces torque on water molecules because of the electric forces acting on the positively charged hydrogen atoms and negatively charged oxygen atom. As a result, the electric field changes the molecular orientation of water (see Figure 2). Figure 2a shows the probability distribution functions (PDFs) of water orientation in an (8,8) CNT at 200 K to 350 K, ambient pressure of 0.1 MPa, and under electric fields of 0 V/nm and 1 V/nm. The orientation is represented with the angle α between the vectors of the dipole moment and the electric field. At 0 V/nm and 250 K to 350 K, water molecules tend to show two peaks in the PDFs, at cos α = 0.8 and −0.8. The peaks are higher at the lower temperature of 250 K than those at the higher temperatures of 300 K and 350 K. The two peaks mean at conditions of 0 V/nm and 250 K to 350 K, water in the CNT tends to having two dipole moment orientations, i.e., on the positive z-axis and the negative z-axis, which are opposite to each other, as shown in a previous study at 300 K [47]. When the temperature decreases at 200 K, the distribution consists of a single orientation with one peak at cos α = 0.9, instead of the two peaks observed at higher temperature. Water molecules have a uniform orientation in the (8,8) CNT at 200 K, i.e., parallel to the CNT axis. Under an applied electric field of 1 V/nm, water molecules in the CNT have a uniform orientation at all simulated temperature conditions. The orientation of the dipole moment of water is restrained parallel to the direction of the electric field. Figure 2b shows the orientation of water molecules in an (8,8) CNT at various pressures of P = 0.1 MPa to 10 MPa, constant temperature of T = 300 K, under applied electric fields of 0 V/nm and 1 V/nm. The change in pressure does not affect the orientation of water molecules in the CNT under any of the electric field conditions. Simulations with the TIP4P model show the same results (see Figure S1 in supporting information). However, the major difference of the TIP4P model from the SPC model is the peak at cos α = −0.8 at the condition of 200 K and 0 V/nm. The peak at cos α = −0.8 is significantly lower than that at cos α = 0.8. This indicates that the TIP4P water shows as a ferroelectric material in the CNT at low temperatures. This result is the same as that shown with the SPC model. Polarization was determined by calculating the average z-component of the dipole moment (µ z ) of water molecules in the CNT. The average value was calculated for 0.1 ns (100 frames) from 22 ns of simulation data, and is shown in Figure 3a,b. At a temperature of 200 K and applied electric field of 0 V/nm, the polarization is µ z ≥ 1.7 Debye (D). This indicates that water in the CNT exhibits ferroelectric behavior at 200 K. This is consistent with the water orientation result shown in Figure 2a, in which one peak is observed at 200 K. The polarization of water molecules at low temperature has also been reported in previous studies [48,49]. However, one of these studies reported that polarization does not occur in a 4-gonal ice-nanotube [49]. As explained in the next paragraph, the structure of water molecules in the CNT at 200 K observed in the current study is 4-gonal. The previous study [49] was based on water in an infinite length CNT, rather than a CNT with reservoirs at each end. The dynamics of water molecules in the reservoirs may affect the polarization of water in the CNT at low temperature. At 250 K and 0 V/nm, the direction of polarization changes with time between µ z = +0.8 D and µ z = −0.8 D. A similar result was reported for a 5-gonal ice-nanotube structure in a CNT [48]. At 300 K µ z fluctuates, and at the higher temperature of 350 K, the amplitude of µ z decreases. This result is similar to that reported in the previous study [48]. Under an applied electric field, the orientation of water molecules is uniform (see Figure 2a,b), so µ z is positive, i.e., 2.1 D. Changes in pressure up to 10 MPa do not affect the polarization of water in the CNT (see Figure 3b). Simulations with the TIP4P model show similar results (see Figure S2). However, at the condition of 300 K and 0 V/nm, the µ z change with time between µ z = +0.8 D and µ z = −0.8 D, which is different from that of the SPC model. As explained in the next paragraph, the structure and dynamics of the TIP4P water is slightly different from those with the SPC model at the condition of 300 K and 0 V/nm. This suggests the melting temperature of the TIP4P model is slightly higher than that of the SPC model.
The densities of oxygen atoms in an (8,8) CNT in the radial direction at various temperatures and pressures are shown in Figure 4a,b, respectively. At E = 1 V/nm, the density distribution is narrower and the peak is higher, compared with at E = 0 V/nm. This is observed for both various temperatures ( Figure 4a) and pressures ( Figure 4b). This implies that the water structure in an (8,8) CNT under an electric field differs from that in the absence of an electric field at all simulated temperature and pressure conditions. This is discussed in more detail in the next paragraph. At lower temperatures, the distribution is narrower and consequently the peak is higher, especially at 0 V/nm ( Figure 4a). This indicates that, at E = 0 V/nm, the water structure changes with increasing temperature. The change in pressure does not change the radial density distribution, either in the presence or absence of an electric field, as shown in Figure 4b. This suggests that increasing the pressure up to 10 MPa does not change the water structure. A similar result is observed using the TIP4P model ( Figure S3). However, the distribution with the TIP4P model is slightly narrower than that with the SPC model at 300 K and 0 V/nm. It is consistent with the result of polarization in the previous paragraph. Another difference is that the distribution with the TIP4P model is slightly wider than that with the SPC model at 350 K and 1 V/nm. This suggests that, at 350 K and 1 V/nm, both water models have different structures.     Figure 5a,b, respectively. Snapshots under applied electric fields of 0 V/nm and 1 V/nm are shown in both cases. The structure of water molecules is helical under an applied electric field of 1 V/nm. Under the applied electric field, the structure is the same at all temperature and pressure conditions. This result differs to that in the absence of an electric field. To better understand the water structure, we unrolled the water structure to obtain a two-dimensional (2D) structure. The radial positions of oxygen atoms in the CNT are averaged to determine the radius of the tube. The unrolled 2D structures are shown in Figure 6a for 200 K, 0.1 MPa, and 0 V/nm, and in Figure 6b for 200 K, 0.1 MPa, and 1 V/nm. The structure is identified using the (n, m) index [50,51], as was used in our previous study [35]. If n = m, the structure is helical. If m = 0, it is an n-gonal structure [10,50,51]. As an example, water molecules in the structure of Figure 6b form an orderly arrangement, i.e., in both the column (magenta-dotted line) and row (black-dotted line) directions. The width of the structure is five column units (magenta arrow), and the number of rows corresponding to the five unit width is two units (i.e., the length of black arrow). Thus, it is a (5, 2) helical structure. In the same manner, the structure shown in Figure 6a is a (4, 0) 4-gonal structure. The difference in the two ordered structures, (4, 0) and (5,2), suggests that an electric field induces a phase transition from solid to solid. The snapshot in Figure 6a also supports our above conclusion that at 200 K and 0 V/nm, the water exhibits ferroelectric behavior. The directions of the dipole moments of all water molecules are on the positive z-axis. The water structures at all temperature conditions are summarized in Table 1. In the absence of an electric field, water has a solid-like structure of (4, 0) at a low temperature of 200 K and 250 K. At 300 K, the structure of water is (4, 0) or disordered. This indicates that the ordered structure of (4, 0) is unstable at these conditions. The solid-like structure of (4, 0) readily changes to the disordered or liquid-like structure, and vice versa [35]. This indicates that the thermodynamic conditions of the simulation are near the melting or freezing point. A similar result was reported in a previous study, in which the formation of a 5-gonal ice-nanotube in a CNT of diameter 1.15 nm reportedly occurs at 310 ± 5 K [48]. The unstable structure at 300 K and 0 V/nm causes oscillation of the dipole orientation because of rearrangement of the water molecules into a tube-like structure (see Figure 3a). At a higher temperature of 350 K, the structure is completely disordered, and the water behaves as a liquid phase. Interestingly, under an applied electric field of E = 1 V/nm, the structure of water molecules in the (8, 8) CNT (diameter 1.08 nm) is (5, 2) at all temperatures. This structure is different from the (4, 0) structure at low temperature in the absence of an electric field. Previous studies of water in CNTs of slightly larger diameters under an electric field at ambient pressure showed different results [24,38,39]. In a CNT of diameter 1.2 nm and under a 1 V/nm field, TIP4P water has a solid-like structure at T < 105 K, and a liquid-like structure at higher temperatures [38]. In a 1.31 nm diameter CNT under a 1 V/nm field, the solid-like structure of TIP4P water reportedly melts to a liquid phase at T > 230 K [39]. These differences may have arisen from the different simulation configurations. In previous studies [38,39], the CNTs are immersed in bulk water, i.e., there are water molecules both inside and outside of the CNT wall. The interaction of water inside the CNT with water outside the CNT may result in a different freezing point. Moreover, in a 1.11 nm diameter CNT of infinite length under a 1 V/nm field, the melting point of TIP4P water reportedly occurs at 310 K [24]. In our system, there is no water outside the CNT wall. The water is instead located in reservoirs at the ends of the CNT. This may result in the different freezing points of water structures within the CNTs. Snapshots of the TIP4P water in the (8,8) CNT are shown in Figures S4 and S5, and the structures at various temperatures are summarized in Table S1. In general, the structures are similar to those of the SPC. Slight differences are at the conditions of 300 K and 0 V/nm and of 350 K and 1 V/nm. At condition of 300 K and 0 V/nm, the (4, 0) structure with the TIP4P model is more stable than that with the SPC model. This indicates that the melting point of the TIP4P is slightly higher than that of the SPC. At 350 K and 1 V/nm, the structure with the TIP4P is (4, 2) instead of (5, 2) with the SPC. Hydrogen bond dynamics analysis shows that the structure of the TIP4P at this condition is weaker than that of the SPC. This suggests that, under an electric field, the melting point of the TIP4P model is slightly lower than that of the SPC model.   Table 2 summarizes the water structures in an (8,8) CNT under various pressure conditions. Under a E = 0 V/nm field, water forms a (4, 0) or a disordered structure that is unstable, as explained above. Increasing the pressure up to 10 MPa does not yield a stable (4, 0) solid-like structure. This implies that the freezing point occurs at a pressure higher than 10 MPa. Under a E = 1 V/nm field, the structure at all pressure conditions is (5,2). This structure is the same as that at all temperature conditions in Table 1. An electric field 1 V/nm causes the formation of a (5, 2) structure at all simulated temperatures and pressures. This differs from the structure in the absence of an electric field. Structures of TIP4P at various pressures are listed in Table S2. The structures are (4, 0) and (5, 2) at 0 V/nm and 1 V/nm, respectively. Pressure up to 10 MPa does not change the structures.

Hydrogen Bond Dynamics
We calculated the hydrogen bond (HB) autocorrelation function C HB (t), to investigate the dynamic behavior of the water structures [52,53]: where h(t) = 1 if a pair of water molecules is bonded at time t, and h(t) = 0 if otherwise. The denominator h 2 is for normalization and is the number of hydrogen bonds (HBs) tagged at t = 0. The HBs between water molecules in the CNT are determined by the following geometrical criteria: R OO < 0.36 nm, R OH < 0.24 nm, and angle of HO − O < 30 • [54,55]. The HB autocorrelation functions C HB (t) within the water structures at various temperatures and 0 V/nm and 1 V/nm fields are shown in Figure 7a. At T = 200 K, 250 K, and 0 V/nm, C HB (t) decreases slightly and then becomes constant with increasing t. The value of C HB (t) is still higher than 0.8 at t = 100 ps. This indicates that the lifetime of the HBs is long. Similar results are observed for HBs within water structures under a 1 V/nm field at all temperature conditions. This indicates that the solid-like structures in the CNT, i.e., both the (4, 0)-structure at 0 V/nm, and the (5, 2)-structure at 1 V/nm, are retained by strong HBs. Different results are observed at temperatures of 300 K and 350 K under a 0 V/nm field, where the C HB (t) decreases to less than 0.2 at t = 100 ps. This indicates that more than 80% of HBs within the water structure are broken. Pairs of water molecules in this structure are not maintained by strong HBs. Consequently, pairs are easily broken, and the two molecules diffuse away from each other. As seen in Table 1, the structure at these temperatures is disordered, i.e., a liquid-like structure. The C HB (t) within TIP4P structures at various temperatures with 0 V/nm and 1 V/nm are shown in Figure S6a. The dynamics of HB of TIP4P structures is similar to that of SPC except at 300 K and 0 V/nm. At this condition, the C HB (t) of the TIP4P is higher than that of the SPC. This implies that HBs within (4, 0)-structure of the TIP4P is stronger, and the structure is more stable than that of the SPC. Moreover, this denotes that the melting point of the TIP4P is slightly higher than that of the SPC. This result supports the explanation in the previous paragraph. At the condition of 350 K and 1 V/nm, the C HB (t) of the TIP4P is slightly lower than that of the SPC. The HBs within the TIP4P structure is weaker. This suggests that, under an electric field, the melting point of the TIP4P is lower than that of the SPC.
The effect of pressure on C HB (t) under 0 V/nm and 1 V/nm fields is shown in Figure 7b. At E = 0 V/nm, the value of C HB (t) decreases to around 0.2 for all pressure conditions. A pair of water molecules within the disordered (liquid-like) structure will readily diffuse away from each other. Increasing the pressure up to 10 MPa has no significant effect on the HB strength within the structure. Under a 1 V/nm field, C HB (t) > 0.8 at t = 100. This indicates that the (5, 2)-structure within the CNT is retained by strong HBs. This supports the notion that the water behaves as a solid phase. Increasing the pressure up to 10 MPa does not change the C HB (t) within the structure in the CNT under a E = 1 V/nm field. The dynamics of HBs in the structures at 200 K and 250 K under 0 V/nm and 1 V/nm fields are similar (Figure 7a). The lifetimes of HBs within the (4, 0)-structure (at 0 V/nm) and (5, 2)-structure (at 1 V/nm) are very long. The C HB (t) with the TIP4P model at various pressures are presented at Figure S6b. The pressure up to 10 MPa does not change C HB (t).

Potential Energy
To identify the main driving force for the formation of the (4, 0)-structure and (5, 2)-structure within the CNTs, we analyzed the Lennard-Jones (LJ) and Coulomb potential energies of water molecules in the reservoir and in the CNT. In addition, under an external electric field, water has a dipole potential energy U dipole = − − → µ · − → E , where − → µ is the vector of the dipole moment, and − → E is the electric field. Table 3 shows averages of the LJ potential energy, Coulomb potential energy, and dipole potential energy per water molecule in the reservoirs and in the CNT at 200 K and 250 K. ∆ U is the difference between the potential energy in the CNT (U CNT ) and that in the reservoir (U Res. ), ∆U = U CNT − U Res . Thus, ∆U indicates the effect of the CNT wall on the water molecules in the tube. There is no electrostatic interaction of water molecules with the CNT and graphene walls because the CNT and graphene have no charge, i.e., the Coulomb potential represents only water-water interactions.
In the absence of an electric field, i.e., E = 0 V/nm, the LJ potential energy in the CNT is lower than that in the reservoirs at T = 200 K and 250 K. The values of ∆U LJ are −14.86 kJ/mol and −11.26 kJ/mol at 200 K and 250 K, respectively. The negative value of ∆U LJ indicates that the van der Waals interaction facilitates water molecules to occupy the tube. On the other hand, the Coulomb potential energy in the reservoirs is lower than that in the CNT. The energy differences ∆U Coulomb are 14.86 kJ/mol and 8.53 kJ/mol at 200 K and 250 K, respectively. This indicates that the electrostatic interaction tends to move the water molecules from the CNT into reservoir. Therefore, at 200 K and 250 K under a 0 V/nm field, the van der Waals interaction is the main driving force maintaining water within the CNT. Similar results are shown by the simulation with the TIP4P water model as presented in Table S3. The values of ∆U LJ are −14.26 kJ/mol and −11.19 kJ/mol at 200 K and 250 K, respectively. Meanwhile, the ∆U Coulomb with the TIP4P are 16.08 kJ/mol and 6.70 kJ/mol at 200 K and 250 K, respectively. To clarify the contribution for the van der Waals interaction, we calculated the LJ potential energy of the water-water interaction (LJ W−W ) and water-CNT interaction (LJ W−CNT ) in the reservoirs and in the CNT (see Table 4). The potential energy LJ W−CNT does not include any contribution from the graphene sheets. Under a 0 V/nm field, LJ W−W in the reservoirs and CNT have similar positive values. In contrast, LJ W−CNT in the CNT has negative values of −13.51 and −13.33 kJ/mol at 200 K and 250 K, respectively. This indicates that, in the absence of an electric field, the van der Waals interaction of CNT-water is the main driving force for the formation of the (4, 0)-structure in the CNT, which facilitates water molecules to occupy the CNT. Simulations with TIP4P model show similar results (Table S4). The LJ W−W in the reservoirs and CNT do not have significant different values. On the other hand, the LJ W−CNT in the CNT have negative values, −12.75 kJ/mol and −13.12 kJ/mol at 200 K and 250 K, respectively.
Under an electric field of 1 V/nm, the LJ potential energy difference between the reservoir and CNT is small. The ∆U LJ values are −1.53 kJ/mol and 0.13 kJ/mol at 200 K and 250 K, respectively (see Table 3). The TIP4P water model shows similar results (Table S3). The values of ∆U LJ are −4.47 kJ/mol and −1.11 kJ/mol at 200 K and 250 K, respectively. We now consider the contribution of the van der Waals interaction. The LJ W−CNT facilitates water molecules to occupy the CNT, while the LJ W−W tends to remove water from the CNT (see Table 4). Simulations with the TIP4P model show the same results (Table S4). In contrast, the Coulomb potential energy in the CNT is significantly lower than that in the reservoirs. The ∆U Coulomb values are −10.02 kJ/mol and −15.59 kJ/mol at 200 K and 250 K, respectively (Table 3). Therefore, under an electric field of 1 V/nm, the electrostatic interaction of water molecules in the CNT is the main driving force for the formation of the (5, 2)-structure in the CNT, and facilitates water molecules to occupy the CNT. Similarly, the presence of an electric field strengthens the electrostatic interaction within the TIP4P water structure. Coulomb potential energy per TIP4P molecule in the CNT at 1 V/nm is significantly lower than that at 0 V/nm (Table S3). At 200 K, the ∆U Coulomb is −0.43 kJ/mol, which is still higher than −4.47 kJ/mol of ∆U LJ . Therefore, at a low temperature of 200 K, the LJ W−CNT is still dominant for facilitating the TIP4P water to occupy in the CNT. However, at higher temperatures of 250 K, ∆U Coulomb becomes −11.79 kJ/mol, which is the main driving energy to facilitate the TIP4P water to fill the CNT. This is due to the Coulomb potential energy of water in reservoirs increasing at a higher temperature. At 300 K and 1 V/nm, ∆U Coulomb decreases to −16.37 kJ/mol.
We now consider the mechanism of phase transition of the water structure from (4, 0) to (5, 2) within the CNT under an electric field, using the energy differences (see Table 3). Under the electric field, the LJ potentials in the CNT increase by 13.27 kJ/mol and 11.47 kJ/mol at 200 K and 250 K, respectively, compared with those in the absence of an electric field. This main contribution to this increase is the water-water interaction (LJ W−W ). In contrast, the Coulomb potentials decrease by −24.47 kJ/mol and −24.11 kJ/mol at 200 K and 250 K, respectively. These results indicate that the phase transition of the water structure from (4, 0) to (5, 2) results from the decrease in Coulomb potential and increase in LJ potential (i.e., mainly water-water interactions). This explanation also applies for the TIP4P model. With 1 V/nm, the LJ potentials in the CNT increase by 9.88 kJ/mol and 10.10 kJ/mol at 200 K and 250 K, respectively, compared with those with 0 V/nm (Table S3). Contrarily, the Coulomb potentials decrease by −15.74 kJ/mol and −18.26 kJ/mol at 200 K and 250 K, respectively.
The dipole potential energy will be at a minimum if the direction of the dipole moment is parallel to that of the electric field. However, the direction of the dipole moment is not statically parallel to the electric field, but instead oscillates around it. The presence of the graphene sheet causes polarization in the reservoirs, and induces an internal electric field in the direction opposite to the external electric field [35]. As a result, the effective electric field in the reservoir is lower than that in the CNT. The orientation of the dipole moment of water in the reservoir therefore oscillates more than that in the CNT. This causes the average dipole potential in the reservoirs to be slightly higher than that in the CNT, so ∆U dipole becomes negative (Table 3). This means that the dipole potential energy promotes water to fill the CNT. However, the contribution of the dipole energy is minor compared to that of the Coulomb potential energy. This supports the notion above that the electrostatic interaction within the water structure is the main driving force for the transition from the (4, 0)-structure to (5, 2)-structure. The ∆U dipole of the TIP4P water is slightly higher than that of the SPC. This is owing to the fact that the dipole moment of the TIP4P is lower than that of the SPC, i.e., 2.18 D and 2.27 D, respectively.
In the absence of an electric field, the LJ potential is the main driving force for the formation of an ordered structure in the CNT. However, the electrostatic interaction within water structure tends to remove water from CNT. At the condition of 300 K and 0 V/nm, ∆U Coulomb of the SPC is higher than that of the TIP4P, i.e., 12.13 kJ/mol and 3.93 kJ/mol, respectively (see Table S5). This makes the ∆U total of the SPC higher than that of the TIP4P, i.e., 1.76 kJ/mol and −5.42 kJ/mol, respectively. As the result, at 300 K and 0 V/nm the solid-like structure of the TIP4P in the CNT is more stable, and the melting point is higher than that of the SPC. Conversely, with the presence of an electric field, the electrostatic interaction within water structure is the main driving force for the formation of an ordered structure. At 350 K and 1 V/nm, ∆U Coulomb of the SPC is lower than that of the TIP4P, i.e., −19.52 kJ/mol and −14.18 kJ/mol, respectively (see Table S6). This makes ∆U total of the SPC lower than that of the TIP4P, i.e., −21.64 kJ/mol and −19.86 kJ/mol, respectively. This causes the structure of the SPC at a high temperature of 350 K under an electric field to be stronger than that of the TIP4P. This suggests that other water models having higher dipole moments, such as TIP3P and TIP4P/2005, may have higher melting point within CNTs. Table 3. Average potential energy per water molecule in the reservoirs and in the (8,8) CNT. ∆U is the difference between the potential energies in the CNT (U CNT ) and in the reservoirs (U Res ), ∆U = U CNT − U Res . The simulation conditions are T = 200 K and 250 K, P = 0.1 MPa, under applied electric fields of E = 0 V/nm and 1 V/nm.  Table 4. Average the Lennard-Jones (LJ) potential energy for the water-water interaction (LJ W−W ) and water-CNT interaction (LJ W−CNT ) per water molecule, in the reservoirs and in the (8,8) CNT. ∆U is the difference between the potential energies in the CNT (U CNT ) and in the reservoirs (U Res ), ∆U = U CNT − U Res .

Conclusions
MD simulations have been performed to clarify the structure of water molecules in an (8,8) CNT, in the presence and absence of an electric field. Simulations have been performed at various temperatures of 200 K to 350 K at constant pressure of 0.1 MPa, and at various pressures of 0.1 MPa to 10 MPa at constant temperature of 300 K. Under a 0 V/nm field at low temperatures of 200 K and 250 K, water in the CNT has an ordered solid-like structure. This indicates that water in the CNT forms an ice-nanotube. At 200 K, the ice-nanotube exhibits ferroelectric behavior. At T = 300 K, the solid-like structure remains but is unstable. The structure readily changes to a disordered structure or liquid-like structure, and vice versa. This suggests that this temperature is near the melting point. At 300 K, increasing the pressure up to 10 MPa does not change the structure. Under these conditions, water has both ordered and disordered structures, i.e., an unstable structure. At a higher temperature of 350 K, the structure is completely disordered, and water in the CNT is in the liquid phase.
Under an applied electric field of 1 V/nm, water molecules form an ordered structure in the CNT at all simulated temperature conditions. This implies that the electric field induces a phase transition from liquid to ice-nanotube in the CNT, at temperatures as high as 350 K. The structure of the ice-nanotube under the applied electric field differs from that formed in the absence of the electric field at low temperatures of 200 K and 250 K. The electrostatic interactions or hydrogen bonding within the ice-nanotube structure under the electric field are stronger than those in the absence of the electric field. The electrostatic interaction is the main driving force for the solid-solid phase transition at low temperatures. The structure of the ice-nanotube under a 1 V/nm field does not change with increasing pressure of up to 10 MPa.