Super square carbon nanotube network: a new promising water desalination membrane

Super square (SS) carbon nanotube (CNT) networks, acting as a new kind of nanoporous membrane, manifest excellent water desalination performance. Nanopores in SS CNT network can efficiently filter NaCl from water. The water desalination ability of such nanoporous membranes critically depends on the pore diameter, permitting water molecule permeatration while salt ion obstruction. On the basis of the systematical analysis on the interaction among water permeability, salt concentration limit and pressure on the membranes, an empirical formula is developed to describe the relationship between pressure and concentration limit. In the meantime, the nonlinear relationship between pressure and water permeability is examined. Hence, by controlling pressure, optimal plan can be easily made to efficiently filter the saltwater. Moreover, steered molecular dynamics (MD) method uncovers bending and local buckling of SS CNT network that leads to salt ions passing through membranes. These important mechanical behaviours are neglected in most MD simulations, which may overestimate the filtration ability. Overall, water permeability of such material is several orders of magnitude higher than the conventional reverse osmosis membranes and several times higher than nanoporous graphene membranes. SS CNT networks may act as a new kind of membrane developed for water desalination with excellent filtration ability. Simulations by researchers in China have shown that square networks of carbon nanotubes can desalinate water at exceptionally high flow rates. Jian Lu and colleagues at the City University of Hong Kong conducted a theoretical study combined with molecular dynamics simulations to understand the salt ion rejection behavior of a carbon nanotube (CNT) mesh with strictly controlled nanopore dimensions. They found that their ‘super square’ CNT network has high permeability for water, permitting flow rates several orders of magnitude higher than conventional reverse osmosis membranes while achieving 100% salt ion rejection. The superior filtration performance is attributed to the high strength and rigidity of the CNT network, which allows it to maintain optimal pore size without flex. These simulations demonstrate the potential of CNT membrane technology for reliable high-flow desalination.


INTRODUCTION
Water scarcity remains a highlighted issue worldwide and becomes more serious due to the various factors such as the growth of population, the increasing global warming and consumptions in industry and agriculture. Although 71% of the earth is covered with water, only~3% of the water all over the world is fresh and suitable for human consumption. 1 In spite of the wide availability of seawater, environmentally friendly and economic purification techniques are needed to produce fresh water from seawater for sustainability. Nowadays, the most widely adopted commercial desalination techniques are electrodialysis, 2 nanofiltration 3 and reverse osmosis (RO). 4 Although ROs, which use high pressure to force water through porous membranes, are believed to have the greatest practical potential, the water permeability is still not satisfactory, meanwhile energy and capital consumption is intensive. In contrast to classical RO membranes, nanoporous membranes can realise fast water permeation. Nanoporous graphene is one kind of nanoporous membranes. Small molecules can pass through them, whereas larger ones cannot. [5][6][7] In addition, nanoporous graphene membranes also can be used as selective molecular sieves for gas purification. 8,9 Therefore, nanoporous graphene is a newly developed membrane for filtration technology. On the other hand, molecular dynamics (MD) studies have predicted 10 and experiments have demonstrated fast transport of gases 11,12 and liquids [13][14][15][16] through carbon nanotube (CNT) membranes. The well-aligned CNTs can serve as robust pores in membranes for water desalination and decontamination applications. [17][18][19] Nanoporous membranes open a bright future for water desalination technology.
Super CNT networks can be obtained by assembling single-walled CNTs (SWCNTs) into ordered micronetworks, which is first designed for application in composites. 20,21 Theoretical and experimental studies show that CNTs can be merged covalently using controlled electron irradiation at high temperatures 22 or atomic welders during heat treatment. 23 Some kinds of CNT networks have already successfully been made and exhibit special mechanical properties. [24][25][26] Along with the production of CNT networks, theoretical models consisting of covalently bonded two-and threedimensional networks have also been developed. 27 MD simulations, atomic finite element method, continuum models and molecular structure mechanics method are developed to study super CNTs. [28][29][30][31] Although existing studies have already shown the potential application of super CNT networks in composites, the potential role of such material for water desalination remains unexplored.
In this work, super square (SS) CNT network is selected as a representative for the study on the water desalination ability. To be specific, the SS CNT networks are constructed from SWCNT arms (6,6), denoted as SS@ (6,6) in this work. The detailed introduction can be found in section 'Methods'. The computational results indicate that SS@ (6,6) can effectively separate salt from water in the desalination system. On the basis of the results obtained from classical MD simulations, optimal plan can be made to filter saltwater with specific salt concentration by controlling the pressure on SS@ (6,6). Furthermore, the deformation and failure mode is also discussed under different pressure on SS@ (6,6), which cannot be taken into consideration in other MD simulations 5-10 because the membranes are completely fixed in their systems. Our calculations demonstrate that water molecules can flow across SS@(6,6) with very high permeability (from 171 to 421 l/cm 2 /day/MPa) while still rejecting salt ions perfectly. This work will be helpful in developing a new kind of excellent nanoporous membrane based on super CNT networks.

RESULTS
Different filtration stages with the variation of pressure We first investigate the pressure on SS@(6,6) (P ss ) and the distance between the two layers of SS@(6,6) (d ss ) during water desalination. The variation of d ss and P ss with the increase of time under constant pulling velocity are shown in Figure 1a. There are mainly three stages during the simulation. In stage I, the filtration velocity gradually reduces accompanied by the increase of pressure (green segment). Note that only very small amount of water molecules are filtered across SS@ (6,6) in stage I, which is called pre-filtration stage. In stage II, the P ss keeps in a stable range and the filtration velocity is almost a constant and equal to v = 0.023 Å/ps (orange segment, filtration velocity is equal to ∂d ss /∂t), which is called stable pressure filtration (SPF) stage, the corresponding stable pressure is denoted as P c . Note that the salt rejection is 100% in stage II. When the salt concentration reaches a critical value (defined as SPF concentration limit), filtration enters into stage III. The filtration velocity gradually reduces once again, accompanied by the increase of pressure (yellow segment). This stage is named as the ascending pressure filtration (APF) stage. However, the salt ions may pass through SS@ (6,6) with the increase of pressure in stage III. Therefore, some rules can be summarised: (1) filtration velocity can be kept as a constant under P c . (2) There is a SPF concentration limit under constant filtration velocity v and corresponding stable pressure P c . (3) High-concentration salt solution beyond the SPF concentration limit can be filtered by gradually increasing pressure provided that no salt ion passes through SS@(6,6). As increasing pressure could induce mechanical problems such as bending, buckling and failure of the SS CNT networks, stage II is an important regime to be paid attention to establish an optimisation criterion. For salt solution with different concentration, we can choose suitable P c to filter it under relative high filtration velocity.
Pressure, concentration limit and water permeability To investigate the relationship between the stable pressure P c and the SPF concentration limit, a series of constant pulling velocities are selected for saltwater desalination. Figure 1b shows the relationship between d ss and P ss under different constant pulling velocities. All of the cases have a relatively stable pressure (P c ) region corresponding to a series of constant filtration velocities from 0.023 to 0.082 Å/ps. Although the velocity is very fast. The filtration velocity can be considered as is equivalent to the water velocity, as the filtration velocity is almost uniform during the stable filtration stage. The water velocity~2.3 to 8.2 m/s is achievable in real life. In addition, we choose the high filtration velocities to control the pressure on SS@(6,6) in a desired range (115-170 MPa) for comparison with the applied pressure (100-200 MPa) on nanoporous graphene. 6 As shown in Figure 1b, it is obvious that there is a negative correlation between P c and SPF concentration limit. However, there are two kinds of circumstances needed to be distinguished. (1) Low SPF concentration limit induced by low P c can be further filtered in APF stage. (2) Low SPF concentration limit induced by high P c cannot be improved in APF stage. Such difference occurs as the salt ions can pass through SS@(6,6) easier under large constant filtration velocities, which can be verified in combination with the data in Table 1 and Figure 1b. For the curve of constant filtration velocity v = 0.032 Å/ps, the salt ions pass through SS@(6,6) when d ss = 1.6 nm, the corresponding P ss is 249.2 MPa, which is much higher than P c . Furthermore, if the constant filtration velocity is even smaller (v = 0.023 Å/ps), yet no salt ions flow across SS@(6,6) even when d ss = 1.44 nm. However, when the constant filtration velocity is large, salt ions can pass through SS@ (6,6) at the very beginning of APF stage. The critical pressure for ion passage is 175 MPa, much smaller than the slow filtration cases that are over 200 MPa (Table 1). So APF is not very suitable for large filtration velocity. Besides, it is discovered that Na + ions always flow across SS@(6,6) earlier than Cl − ions. As the size of Na + ions is smaller than that of Cl − ions, we believe ion size should be an important reason for earlier escape of Na + ions. Figure 2a further exhibits the relationship between stable pressures P c and SPF concentration limits. Note that the salt rejection is 100% for all the MD simulation cases under different P c in Figure 2a. Regions R1 and R2 are consistent to the low SPF concentration limits (as mentioned above) induced by low P c and high P c . If P c is in region R1, the filtration limit can be increased by APF, whereas APF hardly works for further filtration if P c is in region R2. In addition, it can be observed that there exists a critical P c crt corresponding to a critical SPF concentration limit. Using the curve fitting to the solution of the SPF concentration limit (Figure 2a), an empirical formula is established to describe the relationship between the stable pressure and the SPF concentration limit as follow: Where A p is the area of each pore on SS@ (6,6), μ is the total pore area fraction of one layer of SS@ (6,6), and k = 2.69 × 10 6 , h = − 2.78 × 10 4 and c 0 = 130.0 are the fitting parameters. The optimum pressure for water desalination can be achieved using the formula, depending on two key factors: (1) pore size on the SS CNT network; (2) total pore area fraction on the SS CNT network. It can be seen that: (1) the curve based on equation (1) well coincides with the MD simulation results (Figure 2a). (2) μ is positive correlated with the SPF concentration limit (C), whereas A p is negative correlated with the SPF concentration limit (C). It is easy to imagine that if there are more pores on unit area of the SS CNT network, the filtration should be easier with no sacrifice of salt rejection ability, as there is no change of pore size. On the other hand, if the area of each pore (A p ) is enlarged, the salt rejection ability will reduce and result in the reduction of SPF concentration limit (C). Besides, if the pore size A p is extremely small, SPF concentration limit (C) approaches infinity. This situation can be imagined as no salt ions can pass through the membrane without pores, nothing to do with the salt concentration. (3) As a new kind of RO membrane, applied pressure is needed to be exerted on the SS CNT networks to overcome osmotic pressure. Thus, if P c is close to zero, SPF concentration limit (C) converges to zero in equation (1). Furthermore, if P c approaches infinity, the SPF concentration limit (C) also converges to zero in equation (1). Therefore, equation (1) is a reasonable empirical formula to predict the relationship between stable pressures P c and SPF concentration limits. On the basis of the equation (1), the critical stable pressure P c crt is equal to 124.5 MPa and the critical SPF concentration limit is equal to 99.5 g/l for SS@ (6,6). Moreover, all the SPF concentration limits in Figure 2a exceed the salinity of seawater (35 g/l). Hence, such kind of the SS CNT networks is with the prospect of seawater purification.
To estimate the water filtration efficiency of SS@(6,6), we calculate the effective water permeability (W p ) achieved under different P c (Figure 2b). The water permeability, expressed in liters of output per square centimetre of SS@(6,6) per day and per unit of P c , ranges from 171 to 421 l/cm 2 /day/MPa. The filtration efficiency is remarkable. Besides, the permeability scales nonlinearly with P c . There exists the relation that the filtration velocity v rises accompanied by the increase of P c (Figure 1b). In addition, we can easily find that the increase of P c results in the improvement of W p (Figure 2b). Thus, the water permeability should be gradually improved along with the increase of filtration velocity. On the other hand, it is observed that the area of each pore during filtration (A 1 ) becomes larger than the size after relaxation (A 0 ). In addition, the variation of pore area ΔA (ΔA = A 1 -A 0 ) is larger and more intensively under higher P c . Therefore, the relationship between W p and P c is close related to two key factors: (1) filtration velocity v under different P c ; (2) ΔA under different P c during filtration. The combined effect of them results in the nonlinear relationship between P c and W p .
Hence, the relatively optimal filtration pressure can be selected with consideration of the relationship among pressure, salt concentration and water permeability discovered in this work as mentioned above.  Figure 2. (a) The relationship between stable pressures (P c ) and filtration concentration limits under stable pressure for SS CNT SS@(6,6). Filtration limit can be increased in region R1, whereas cannot be increased in region R2. There exists a critical stable pressure P c crt corresponding to a critical stable pressure filtration concentration limit. (b) The nonlinear relationship between stable pressures (P c ) and water permeability for SS CNT SS@ (6,6). Error bars represent the s.d. of pressures in the stable pressure filtration stage. The illustrations for nanopore size in CNT SS@(6,6) under different pressure. The pore sizes vary more intensively under higher pressure.  (Table 2). Figure 3 shows the nanostructure evolution of sample I during water desalination compared with nanoporous graphene. Figure 3a is one of the pictures obtained during water desalination at SPF stage and Figure 3c is the corresponding figure in which water molecules are invisible. Figure 3b is one of the pictures obtained during water desalination at APF stage and Figure 3d is the corresponding figure only showing salt ions. It can be easily observed that the salt ions can be strictly filtered by SS@ (6,6). Although the nanoporous graphene has the similar pore size and the same loading procedure as SS@(6,6), the salt ions cannot be impeded by it due to easier mechanical deformation of graphene under compression (Figure 3e). Moreover, the perfect heavy ions rejection ability (Cd 2+ , Mg 2+ and Zn 2+ ) of SS@(6,6) is also verified by similar filtration process (Figure 3f). Thus, SS CNT network has more excellent mechanical properties acting as a kind of nanoporous membrane compared with nanoporous graphene. On the other hand, the low SPF concentration limit can be improved in APF stage that is consistent with the R1 region in Figure 2a, as the filtration velocity in sample I is relatively slow. To improve the permeability of SS@ (6,6), the filtration velocity is increased by raising P c . However, increasing pressure could induce mechanical problems such as bending, local buckling and failure of the SS CNT networks (Figure 4). All these deformations occur in the APF stage through observation of nanostructure evolution during filtration. In addition, the faster the filtration velocity, the sooner the deformations occur in APF stage. The pore size is enlarged after bending or local buckling. Thus, the salt ions can flow across SS@(6,6) in the APF stage and SS@(6,6) lose the filtration effectiveness. Note that the salt ions can pass through SS@(6,6) only when it deforms. Hence, even better filtration capability can be achieved by improving the strength of SS CNT networks through structure design. As the deformation and failure mode of the SS CNT networks is affected by the chiralities and sizes of the SWCNTs, 32 they can be taken into consideration to design the SS CNT networks with higher strength in the future study.
Outstanding water permeability and salt rejection ability Overall, our results indicate that the SS CNT networks could act as a high-permeability desalination membrane with guarantee of 100% salt rejection and water passage. High-concentration salt solution can be purified efficiently by the SS CNT networks. Moreover, to illustrate the water permeability and salt rejection ability of the SS CNT networks, the water desalination performance of the SS@(6,6) from P c = 116.2 to 169 MPa is plotted along with the performance of RO membranes 33 and nanoporous graphenes 6 in Figure 5. Note that although the MD simulation method adopted in our study is different from most of other studies, 5-10 we only compare the data in SPF stage (Figure 2b) with other membranes. As the filtration velocity in SPF stage is almost uniform under constant stable pressure, we think the pressure on membranes can approximately be regarded as the applied hydrostatic pressure. Exactly as the performance of nanoporous graphenes, the salt rejection of SS@(6,6) maintains at 100%, both are superior than conventional RO membranes. On the other hand, the water permeability of SS@(6,6) is several orders of magnitude higher than RO membranes and several times higher than nanoporous graphene. It should be emphasised that the water permeability of SS@(6,6) can be improved by increasing pressure with no sacrifice of salt rejection (Figure 2b). Moreover, the functionalisation by hydrophilic chemical group (hydroxyl, OH-) instead of hydrophobic ones (hydrogen atoms, H-) on nanoporous graphene can roughly double the water flux. 6 We think sidewall functionalisation of SS CNT networks to improve hydrophilicity could be an alternative to further improvement the performance of water desalination since both CNTs and CNT networks are chemically hydrophobic. 34 Therefore, the SS CNT networks could be a new kind of nanoporous membrane possessing outstanding water permeability.
Regarding the question of whether the SS CNT networks could be practical for commercial-scale desalination, we think that the fabrication of well-aligned SS CNT networks with the extremely narrow pore size is a main challenge. The well-aligned SS CNT networks with the strict controlled pore size and shape are not achievable due to the limitation of the current fabrication technology. Nowadays, the experimentally available CNT networks [22][23][24][25][26] are composed of a random distribution of pores, in both size and geometric shape. The randomness of pore size will inevitably result in the existence of unfavourable large pores. As shown in Figure 3g, the super graphene CNT network is developed and tested by MD simulation. The worse salt rejection is observed due to the enlargement of the pore size. Thus, it is essential to fabricate well-aligned networks with the uniform pore size for practical application in the future. Although challenging,   6 and SS@ (6,6) is expressed by the interval based on the simulation results. The SS@(6,6) examined in this work could fully reject salt ions with a water permeability several orders of magnitude higher than commercial RO and also have better performance than functionalised nanoporous graphene membranes. (b) The amplified red region in a clearly exhibits the outstanding water permeability of SS@(6,6) compared with nanoporous graphene membranes.
we believe it can be overcome and the CNT networks could offer important advantages over existing RO technology.

DISCUSSION
Acting as a new kind of nanoporous membrane, the SS CNT networks are demonstrated possessing excellent water permeability and salt rejection ability in our simulations. The water permeability is several orders of magnitude higher than the conventional RO membranes and also several times higher than nanoporous graphene membranes. We find that the water permeability with guarantee of 100% salt rejection is both affected by the salt concentration of water and the pressure on the SS@ (6,6). To systematically describe the interaction among pressure, salt concentration and water permeability, the following two efforts are made: (1) an empirical formula is developed to describe the relationship between pressure and concentration limit; (2) the nonlinear relationship between pressure and water permeability is uncovered. Thus, optimal plan can easily be made to efficiently filter saltwater with specific salt concentration by controlling the pressure on the SS@ (6,6). On the other hand, steered MD simulation method is first adopted to analyse the mechanical behaviour of membranes during water desalination, which cannot be realised due to the limitation of methods adopted by other researchers. [5][6][7][8][9][10] On the basis of this method, we discover that the bending and buckling of SS@(6,6) under high pressure, which results in enlargement of the pore size, is the main reason for salt ions passing through SS@ (6,6). Thus, better filtration capability can be achieved by improving the strength of SS CNT networks through structure design and optimising the pore size in the SS CNT networks. Therefore, a new kind of nanoporous membrane is studied in this work for water desalination that yields significant improvement over existing technological methods. We expect that our study will promote the understanding of the excellent performance of well-designed SS CNT networks for water desalination technology.

MATERIALS AND METHODS
All the MD simulations are performed using the LAMMPS package. 35 Two layers of the SS CNT network are constructed from SWCNT arms (6,6), denoted as SS@(6,6) ( Figure 6). The X junction is built by joining four SWCNTs at an angle of 90°using pentagons and heptagons. Thus, uniformly distributed nanopores are fabricated in SS@(6,6) ( Figure 6c) with 25,000 water molecules, 250 Na + ions and 250 Cl − ions between the two SS@(6,6) layers, corresponding to the salt concentration of 37.125 g/l (similar to the salinity of seawater 35 g/l). The SS@(6,6) has dimensions of 11.62 × 10 nm (Figure 6c). The periodic boundary condition is imposed along X, Y and Z directions with the two vacuum layers along X direction. CHARMM force field 36 is adopted to model the interactions between the atoms in the system. Water is modelled by TIP3P potential 37 and interactions between the atoms are modelled using Lennard-Jones and Coulombic terms. Steered MD method 38 is adopted to simulate desalination process under constant pulling velocity. Note that the constant pulling velocity represents the moving speed of equilibrium position of the spring attached on the SS@(6,6), thus the SS@(6,6) are pulled under a traction force F = kΔx, where k is the spring constant and is set to 50 kcal/mol/Å, Δx is equal to L − L 0 , L and L 0 are the actual length and equilibrium length of the spring respectively. As the SS@(6,6) endure resistance from salt solution, the filtration velocity should be slower than the constant pulling velocity. Moreover, unlike other methods in which the rigid condition is applied on the membranes, 5-10 the SS@(6,6) atoms are not held fixed during the steered MD simulation, so the nanostructure evolution of SS@(6,6) under compression can be detected (Figure 6g). Besides, the salt concentration between SS@(6,6) changes continuously during the desalination process. We can easily know the performance of SS@(6,6) for water desalination with specific salt concentration under different pressure according to Figure 1. Thus, the desalination process affected by salt concentration can also be analysed.
Initially, the created system is subjected to energy minimisation using the conjugate gradient method, and then relaxed in the Nose/Hoover isobaric-isothermal ensemble (NPT) under both the pressure 0 bar and the desired temperature (300 K) with an initial Gaussian velocity distribution consistent with the temperature for 300 ps. As shown in Figure 6, the distance between the two layers of SS@(6,6) changes from d ss = 9.89 to 7.77 nm after relaxation. The pore size changes from 46.33 to 48.45 Å 2 . Note that there exists trivial size difference among the pores on the relaxed SS@ (6,6). Hereafter, desalination simulation is executed along X direction under a canonical ensemble (NVT). The total simulation time is 1,500-3,500 ps with a time step of 1 fs. A series of constant pulling velocity v are selected to simulate the desalination process. The average stress of all the carbon atoms is calculated during the simulation to detect the pressure on SS@(6,6) (P ss ). The distance d ss is also recorded to calculate the filtration velocity and salt concentration variation.