Graphene oxide comprises graphene sheets which are decorated by certain functional groups such as hydrophilic oxygen . Recently, numerous computational simulations have been performed on the ultra-filtration using nanotubes -, zeolite  and graphene  . All these theoretical investigations reveal a rapid and effective ultra-filtration using nanomaterials in comparison to conventional membrane. Reverse osmosis has found to possess the merit of low energy consumption and operation costs in comparison to other desalination methods such as multi-stage flash, multi-effect distillation and mechanical vapor compression. Although most theoretical investigations reveal a superior filtration performance of nanomaterials, none of these nanomaterials have successfully been proved for the desalination. Recently, Tanugi & Grossman  have nominated porous graphene as a mean for the seawater desalination due to its superior mechanical strength and the negligible thickness, and they claim that such graphene possesses the highest water permeability over all existing commercial RO and zeolite . Besides, several experiments have paved a way to producing and transferring singly layered graphene sheets , which make this 2D-nano structure more applicable for large scale desalination.
However, the computational approach adopted by Tanugi & Crossman  and other computational simulations - have severe temporal and spatial constraints that are infeasible to tackle this ultra-filtration problem in longer run and larger scale. Here, we adopt the continuous approximation introduced by Cox et al.  , where the pairwise interactions between ions, molecules and functional groups are coarse grained so that the time-consuming calculations could be fasten by performing certain line integrals, and single and double surface integrals. The present authors have successfully applied such approach on solving seawater desalination problems using functionalized nanotube and graphene  -.
In this paper, we investigate the possibilities and limitations of using functionalized graphene sheet for the saltwater desalination. Carbon atoms are assumed to be smeared across the surface of the graphene sheet, and the water and sodium ion are smeared across an envisaged surface of radius 1 and 1.16 Å, respectively, so that the van der Waals interactions between any two ions/molecules could be approximated by performing surface integrals. Chloride ion is neglected because its ionic size is larger than the pore size of the graphene. In addition, the functional groups are assumed to be evenly distributed around the pore boundary so that the electrostatic interactions between the intruders and the functional groups are evaluated in terms of line integrals. Least but not last the interactions between molecules/ions and solute are tabulated from MD simulations results. Such approach has been shown to speed up computational times dramatically and yield accurate results. The resulting force fields could then be incorporated into the MD algorithm to determine the time evolution of the physical system numerically.
2. Theory and MD Simulations
In this section, we introduce the theoretical and numerical basis for the rest of the paper. For simplicity, we model the van der Waals interactions between a water molecule or a ion and the functionalized graphene by the usual 6 - 12 Lennard-Jones potential , and the electrostatic interactions arising from the functional group by electrostatic potential E, which are given, respectively by
where ρij, ε, σ, Qi, ε, A and B denote the atomic distance between two typical atoms, the Lennard-Jones potential well depth, the Lennard-Jones distance between two atoms, total charge for i atom /functional group, the permittivity, the attractive constant and repulsive constant, respectively. The numerical values for several constants are summarized in  . We assume that the carbon atoms are smeared across the graphene surface, and the oxygen atom and the sodium ion are smeared across the envisaged sphere of radius 1 and 1.16 Å, respectively so that the van der Waals forces between the oxygen or the sodium ion and the graphene sheet could be approximated by the double surface integrals, while the van der Waals forces between the hydrogen atom (point particle) and the graphene are determined by the single surface integral, which are given respectively by
where η1, η2, −dV/dρ|┴, dS1 and dS2 denote the number density of carbon atoms on the graphene, the number density of the envisaged sphere, the axial force generated by the Lennard-Jones potential given in Equation (1), the surface element of the graphene and the surface element of the envisaged sphere, respectively. On the other hand, the functional group is assumed to be smeared on the pore vicinity so that the electrostatic forces between water or sodium ion and the functional group is given by
where qi, σ and dl denote the charge of the intruder, the charge density on the pore vicinity and the length element of the pore, respectively. We comment that all the above integrals reduce the multi-body interactions into several two-body problems and the introduction of the axial forces reduces the 3D problem into degenerated 1D problem. According to the prior works done by the present authors -, certain parts of these integrals, i.e. Equations (2) and (3) could be solved analytically facilitating rapid computational times. Since the orientation of water molecules is found to be crucial in this problem , we initially fix the configuration of the two hydrogens separated by an angle 2ϕ as shown in Figure 1 and different orientations of the water molecule could be generated by two independents angles θ in [0, π] and ψ in [0, 2π] so that the kinematics between the first hydrogen, i.e. upper hydrogen in Figure 1 and an arbitrary point on the graphene is written as
where a, l, R, ϕ and Φ denote the radius of the water molecule, the distance between the pore center and an arbitrary point on the graphene, the distance between the center of the water molecule and the pore center, half of the angle between two hydrogens and the polar angle on the graphene, respectively (see Figure 1 for details). The distance between the oxygen or the sodium ion and the graphene sheet or functional group can be deduced in a similar fashion. Such atomic distances could then be incorporated into Equations (2) and (3) to determine F1vdW, F2vdW and Fie. In particular, since the forces between the hydrogen atoms and the functionalized graphene are generated by different orientations, the ensemble forces acting on the two hydrogens could be deduced using Boltzmann’s statistics
where j, FHi(R)j, β and VHi(R)j denote the different orientations, forces of Hi in j configuration, the reciprocal of the usual Boltzmann’s constant times the temperature and the potential energy between Hi and the functionalized graphene, respectively. The total molecular forces between the water molecule and the functionalized graphene is therefore written as
where FivdW and Fie denote the van der Waals forces between the i atom, where I = O, H1 and H2 and the graphene and the electrostatic forces between i atom and the functional group on the pore, respectively. In addition, Φ, Fhydra, Fapp, FW denote the electric field generated by the charges on the entire graphene, the hydraulic forces in the solute, the applied force and the molecular forces comprising van der Waals and electrostatic forces between the water molecule and the graphene, respectively. Here, the external forces Fext is also coarse grained by the sum of both the hydraulic and applied forces. It is worthy to note that the hydraulic forces are interpolated from experimental data or simulation results. The total force between the sodium ion and the functionalized graphene sheet can be deduced in a similar way, i.e. Ftot = FI ? ΔΦ + Fext, where FI denotes the ionic forces between the ion and the functionalized graphene. Given these force fields Ftot, the time evolution of the system can be numerically determined using Verlet algorithm 
Figure 1. Kinematics of a water molecule intruding axially into a pore on the graphene, where ρ denotes the atomic distance between a hydrogen atom and an arbitrary point on the graphene sheet, and the functional group is assumed to be smeared on the ring of the pore entry. Both the applied and hydraulic forces are shown in terms of different arrows.
where τ, k, M and V denote the time grid, the time step, the mass of the intruder and the velocity fields, respectively. To facilitate better numerical outcomes, multistep technique  is also used in significant force regimes.
3. Numerical Results and Discussion
In this section, we carry out some numerical results as given in Section 2. Parameters adopted in this paper could be found in . By substituting them into Equations (2)-(6) and assuming the room temperature of T = 300 K, we determine the total forces consisting of both van der Waals and electrostatics forces for a water molecule and a sodium ion with the total local surface charge Q = 0e and 0.02e. 100 different water molecule’s orientations are considered and the total forces for water molecule and the sodium ion permeating into the functionalized graphene are shown in Figure 2.
Here we ignore the case of chloride ion as the ionic size of the ion is larger than that of the pore size so that chloride ions are always repelled from the graphene. From Figure 3, for two proposed Q, both the water and sodium ion experience repulsive forces but occur in different magnitudes when they approach the graphene. Without any external forces, both the water and the sodium ion could not penetrate through the graphene. While the external forces are applied into the system, the salt water will be pushed towards the sheet in order to overcome the repulsive forces induced by the graphene. Since the repulsive forces of the sodium ion are stronger than that of the water in the case of the neutral graphene, there must exist a range of the external forces, so that for a given container size R03, the water molecule will get through the graphene while the sodium ion will be rejected from the graphene. Once the water (or the sodium ion) reach the pore center, the repulsive forces generated by the permeate side of the graphene will automatically suck them out of the sheet so that the acceptance criteria for the desalination could be written as
Figure 2. Total forces of the water molecule and sodium ion entering the functionalized graphene where the positive forces indicate the repulsive forces generated by the graphene for Q = 0e (neutral) and 0.02e. The arrow indicates the growth of repulsive forces for water molecule when positive charge of 0.02e is doped on the boundary of the pore vicinity.
Figure 3. Total forces of the water and sodium ion penetrating into the functionalized graphene for different Q.
where R0, FI and FW denote the initial position of the intruder, the ionic forces of the sodium ion and the molecular forces of the water molecule, respectively. For neutral graphene, we have the following findings: There exists an energy gap between the water molecule and the sodium ion so that given the external forces that satisfy Equation (8), only water will pass through the graphene resulting in desalination. The larger the strength of the hydraulic and the size of the container, the lower the feasibility for the desalination, especially thermal fluctuation is taken into account. However, when tiny charges for example Q = 0.02e is present on the pore surface, the validity of Equation (8) is broken due to the fact that the water encounters stronger repulsive forces than that of the sodium ion (See Figure 2). In this case, the desalination fails. To demonstrate such effect, we plot the total forces of the water and the sodium ion, where Q ranges from −1e to 1e in Figure 3.
The irregular force fields of water result from the averaging process in Equation (5). For positively charged graphene, water molecule tends to experience larger repulsive forces than that of the sodium ion due to the interaction between the two outer hydrogens and the functionalized graphene. Due to the symmetry of the sodium ion, all electrostatics forces are eliminated at the pore center for all Q resulting in a meeting point of the total forces given in Figure 3. We observe that the functionalized graphene always accepts counter ion.
Now, the force fields of both the water molecule and sodium ion are inserted into the MD simulation frameworks, i.e. Equation (7) to investigate the charge effect on the desalination. Initially, we assume constant external forces Fext = 4e−4 eV/Å and container size of R03 = 10 Å3 for which two different total charge on the pore vicinity, i.e. Q = 0e and 0.02e are considered. From Equation (8), when Fext is between 3.547e-4 and 6.047e−4 eV/Å, only water could penetrate through the neutral graphene, which is confirmed by Figure 4.
However when the graphene sheet is doped by a tiny charge Q = 0.02e and under the same external force of 4e−4 eV/Å, both water and sodium ion are blocked by the positively charged graphene. The path of the sodium ion with and without functional group collides almost with each other’s which could be partially explained by the comparable repulsive force of the sodium ion as shown in Figure 2. In other words, the positively charged graphene induces more repulsive forces on water molecule than that on sodium ion. However, when Fext lies between 1.7e−3 and 12.5e−3 eV/Å, a peculiar phenomenon occurs, where only sodium ion could pass through the positively charged graphene leaving both water and chloride ions on the feed side, which is again demonstrated in Figure 4. To summarize, negatively charged graphene produces acidic solution on the permeate side; small doped positive charges on the graphene severally hamper the efficiency of the desalination; and for certain external forces, mere sodium ions could pass through the positively charged graphene.
In this paper, we adopt both the continuous approximation and MD simulations to investigate the desalination using the functionalized graphene. Small positive surface charges are shown to severely affect its desalination functionality. More intriguingly, while the negatively charged graphene produces acidic solution on the permeate side, for certain ranges of external forces, only sodium ions are sucked into the positively charged graphene. The present methodology has the merit of rapid computational speed and generates inductive results, where the total forces are determined using surface and line integrals.
Figure 4. MD simulations for the permeation of water and sodium ion into the graphene with different Fext and Q = 0e and 0.02e.
We gratefully acknowledge the constructive discussion with Dr. Yung-Tsang Chen, and the financial sup- port from small research grant (UNNC), Ningbo Natural Science Foundation (2014A610025), (2014A610172) and (2015A610281), Zhejiang Natural Science Foundation (Q15E090001) and Qianjiang Talent Scheme (QJD1402009).
 Berezhkovskii, A. and Hummer, G. (2002) Single-File Transport of Water Molecules through a Carbon Nanotube. Phys. Rev. Lett., 89, Article ID: 064503. http://dx.doi.org/10.1103/physrevlett.89.064503
 Zuo, G., Shen, R., Ma, S. and Guo, W. (2010) Transport Properties of Single-File Water Molecules inside a Carbon Nanotube Biomimicking Water Channel. ACS Nano, 4, 205-210. http://dx.doi.org/10.1021/nn901334w
 Sansom, M.S.P., Shrivastava, I.H., Ranatunga, K.M. and Smith, G.R. (2000) Simulations of Ion Channels Watching Ions and Water Move. Trends. Biochem. Sci., 25, 368-374. http://dx.doi.org/10.1016/S0968-0004(00)01613-3
 Cox, B.J., Thamwattana, N. and Hill, J.M. (2007) Mechanics of Atoms and Fullerenes in Single-Walled Carbon Nanotubes. I. Acceptance and Suction Energies. Proc. R. Soc. London, Ser. A, 463, 461. http://dx.doi.org/10.1098/rspa.2006.1771
 Cox, B.J., Thamwattana, N. and Hill, J.M. (2007) Mechanics of Atoms and Fullerenes in Single-Walled Carbon Nanotubes. II. Oscillatory Behavior. Proc. R. Soc. London, Ser. A, 463, 477. http://dx.doi.org/10.1098/rspa.2006.1772
 Chan, Y. and Hill, J.M. (2011) A Mechanical Model for Single-File Transport of Water through Carbon Nanotube Membranes. J. Membrane Sci., 372, 57-65. http://dx.doi.org/10.1016/j.memsci.2011.01.040