Scholarly article on topic 'Water desalination with a single-layer MoS2 nanopore'

Water desalination with a single-layer MoS2 nanopore Academic research paper on "Nano-technology"

Share paper
Academic journal
Nature Communications
OECD Field of science

Academic research paper on topic "Water desalination with a single-layer MoS2 nanopore"



Received 15 Apr 2015 | Accepted 9 Sep 2015 | Published 14 Oct 2015^^BoOl:ia1o38/ncomms9616| OPEN

Water desalination with a single-layer MoS2 nanopore

Mohammad Heiranian1, Amir Barati Farimani1 & Narayana R. Aluru1

Efficient desalination of water continues to be a problem facing the society. Advances in nanotechnology have led to the development of a variety of nanoporous membranes for water purification. Here we show, by performing molecular dynamics simulations, that a nanopore in a single-layer molybdenum disulfide can effectively reject ions and allow transport of water at a high rate. More than 88% of ions are rejected by membranes having pore areas ranging from 20 to 60 A2. Water flux is found to be two to five orders of magnitude greater than that of other known nanoporous membranes. Pore chemistry is shown to play a significant role in modulating the water flux. Pores with only molybdenum atoms on their edges lead to higher fluxes, which are ~70% greater than that of graphene nanopores. These observations are explained by permeation coefficients, energy barriers, water density and velocity distributions in the pores.

1 Department of Mechanical Science and Engineering, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA. Correspondence and requests for materials should be addressed to N.R.A. (email:

NATURE COMMUNICATIONS | 6:8616 | DOI: 10.1038/ncomms96161

Producing fresh water is currently a great challenge facing the society1-4. High capital costs and low efficiency of current desalination technology motivate the need for advances in desalination technology5'6. Approximately, half of the current desalination plants use reverse osmosis (RO) technologies2'5. RO based on traditional polymeric membranes faces several challenges including slow water transport7'8. Advances in nanotechnology open up opportunities to design energy-efficient membranes for water desalination9'10. Nanopores with diameters ranging from a few Angstroms to several nanometres can be drilled in membranes to fabricate molecular sieves11-13. As the diameter of the nanopore approaches the size of the hydrated ions, various types of ions can be rejected by nanoporous membranes promising efficient water desalination. Among nanoscale materials' graphene and carbon nanotubes were extensively studied for both water transport and desalination14-18. Graphene, a single-atom-thick membrane (0.34 nm) was demonstrated to have several orders of magnitude higher flux rates compared with conventional zeolite membranes6'11'15'16'19'20. Since water flux through a membrane scales inversely with the membrane's thickness11' graphene is attractive over most other materials due to its single-atom thickness12'16.

It has been shown that chemical functionalization of a graphene nanopore (for example' adding hydroxyl groups) can enhance its permeability19'20' but reduces desalination efficiency19. Hydroxyl groups provide hydrophilic sites at the edge of the pore' which give rise to the attraction of water molecules and enhanced flux due to denser packing of water inside the pore19. Adding precise functional groups to the edge of nanopores requires complex fabrication21; therefore' identifying a single-atom-thick membrane with hydrophilic sites can lead to further advances in water desalination technology.

Recently' a nanopore in a single-layer molybdenum disulfide (MoS2) has been investigated for DNA sequencing and has been shown to provide better results compared with graphene nanopores9'22. Compared with graphene' a MoS2 single layer has two types of atoms' that is' molybdenum (Mo) and sulfur (S). A single-layer MoS2 has a thickness of ~ 1.0 nm (ref. 23) and is a mechanically strong material with an effective Young's modulus of 270 ±100 GPa' that is comparable to that of steel24. The possibility to craft the pore edge with Mo' S or both provides flexibility to design the nanopore with desired functionality. Recently' it has been shown that a nozzle-like structure of protein channels and other nanoscale membranes enhances water permeation25. The fish-bone structure of MoS2 makes it amenable for a nozzle-like sub-nanometer pore for fast water permeation25.

Although theoretical studies of membrane efficiency are important in desalination technology' there are other aspects concerning fabrication and manufacturability of membranes such as large-area synthesis with defect-free' well-defined sealed membranes and precise pore generation that need to be addressed. Using a highly focused electron beam' and transmission electron microscope' versatile nanopores with diameters ranging from 1 to 10 nm were sculpted successfully in MoS2 membranes9. Waduge et al.26 reported that a large-area' well-sealed membrane with nanopores as tiny as 2.8 nm can be fabricated. Compared with graphene' the contamination of these membranes can be lower as carbon atoms in graphene are more susceptible to contamination during chemical vapour deposition (CVD) growth. Feng et al.27 also achieved high-quality scalable fabrication of nanopores in a single-layer MoS2 with sub-nanometre precision using electrochemical reaction. Several other studies have been performed on the synthesis of large-

area MoS2 monolayers28-37. Recently' a few groups29'34'37 have

successfully used CVD to produce highly crystalline MoS2 of

centimetre dimensions. In another study36, a refined CVD method was proposed to create high-quality monolayer MoS2 crystals in which the grain boundaries of MoS2 were faceted more strongly than that of graphene resulting in mechanically more stable MoS2 monolayers. Membrane sealing also plays an essential role in the synthesis of large-area membranes required in desalination. Waduge et al.26 showed that their CVD approach resulted in almost fully sealed MoS2 membranes. Combination of these results9'13'26-37 and the recent focus on a single-layer MoS2 fabrication is promising for the large-scale manufacturing of a single-layer MoS2.

Here we demonstrate that a single-layer MoS2 can effectively separate ions from water. Using molecular dynamics simulations, we investigate water desalination in MoS2 as a function of pore size, chemistry, geometry and applied hydrostatic pressure.


Water fluxes. A typical simulation box consists of a single-layer MoS2, a graphene sheet (acting as a rigid piston to apply the external pressure), water and ions (Fig. 1a). Here three pore edge types for MoS2 are considered to study the effect of terminating atoms and pore chemistry on the rate of water permeation and ion rejection. The first type of pore, which is labelled as mixed in this study, is a combination of molybdenum and sulfur atoms. The other two pore types are labelled as Mo only and S only, as these are terminated by molybdenum and sulfur atoms, respectively (Fig. 1b). Water fluxes through various MoS2 nanopores as a function of the applied pressure gradient are presented in Fig. 2a. Three MoS2 pore types (mixed, Mo only and S only) were studied to explore their rejection rate and flux. To investigate the relative performance of MoS2 over other two-dimensional materials, a graphene nanopore, which has been shown to be

Rigid piston Salt water MoS2 nanopore Fresh water

S only

y Mixed

Figure 1 | Simulation box and different pore architectures. (a) Schematic of the simulation box consisting of a MoS2 sheet (molybdenum in blue and sulfur in yellow), water (transparent blue), ions (in red and green) and a graphene sheet (in gray). (b) Left: Mo only pore type. Right: Sonly pore type. Bottom: mixed pore type.

50 100 150 200 250 300 350 External pressure (MPa)

-O- Mixed A=36.16 A2 + Mixed A=55.45 A2 —■— Mo only A = 18.02 A2

--Mo only A=51.03 A2

- - Mo only A=56.42 A2 <]- S only A=19.01 A2 -x- S only A=57.38 A2 ♦ Graphene A=59.67 A2

50 100 150 200 250 300 350 External pressure (MPa)

C Mo only at 250 MPa

Time (ns)

Figure 2 | Water permeation and salt rejection. (a) Water flux as a function of the applied pressure for mixed, Mo only, Sonly and graphene nanopores with similar pore areas. (b) Percentage of ion rejection by various pores as a function of the applied pressure. Pores with different edge chemistries as well as various pore areas (denoted by A) are considered. (c) Number of water molecules (#) filtered through Mo only pores as a function of simulation time for different pore areas at a fixed pressure of 250 MPa.

promising for water desalination, is also considered11,19. For the sake of comparison, the three MoS2 pores and the graphene pore have approximately equivalent accessible pore areas (mixed, A = 55.45 A2; Mo only, A = 56.42 A2; S only, A = 57.38 A2; and graphene, A = 59.67 A2). Our results indicate that the Mo only pore has the highest rate of water permeation followed by the mixed, S only and the graphene pore for all the applied pressures (Fig. 2a). Water flux through the mixed pore is intermediary between Mo only and S only nanopores. The higher water fluxes through MoS2 nanopores compared with graphene nanopores imply that for a desired water flux, a smaller applied pressure is needed with MoS2 nanopores. Later, in this paper, we will explain the physical chemistry and geometrical foundations of MoS2 pore that give rise to a higher flux.

Salt rejection efficiency. The other important aspect in water desalination is the ability of the membrane to reject ions. The percentage of total ions rejected by the MoS2 and graphene pores is plotted as a function of the applied pressure in Fig. 2b. The rejection is calculated after 1,700 water molecules have filtered through the pores for all pressures. Pore sizes ranging from 20 to 60 A2 are considered for the three types of MoS2 pores. The ion rejection decreases at higher pressures as high pressures induce higher forces on the ions giving rise to more ion translocation events. The ion rejection of small pores (for example, 18.02 A2) is found to be 100% for all types of pores. For larger pore sizes, ions escape through the pore reducing the rejection efficiency. For the pores with equivalent areas (mixed, A = 55.45 A2; Mo only, A = 56.42 A2; S only, A = 57.38 A2; and graphene, A = 59.67 A2), the general trend for ion rejection is quite similar regardless of the type of the pore (Fig. 2b). In other words, ion rejection is mainly dependent on the pore area and the type of the pore plays a less important role, for example, for the four pores considered, the difference in rejection is < 10% even at a high pressure of 350 MPa.

As shown in Fig. 2c, the water filtration rate increases sharply as the pore area increases from ~ 20 to ~ 50 À2. The sharp change in the water flow rate is due to the formation of single-file chain of water in small pores (~20À2). As shown in ref. 11, the water flow rate is considerably reduced because of the weak hydrogen bonding in single-file chains. For efficient water desalination, pore sizes should be chosen such that both the ion rejection and water filtration rate are optimized since very small pores lack high permeation rates and large pores (wider than 60 À2) fail to effectively reject ions.

As observed by Cohen-Tangui et al.19 for graphene, the polarizability of water also has a little effect on ion rejection in MoS2 nanopores. To introduce the effect of polarization, the flexible simple point charge (SPC/F) model38 was used. The ion rejection percentages associated with the flexible water model are within 2% of those modelled with the SPC/E water.

Permeation coefficient. To quantify the water permeability through various pores, we compute the permeability coefficient, p, across the pore. For dilute solutions39,

- VwDCs +

where Jw is the flux of water (# ns 1), Vw is the molar volume of water (18.91 ml mol- 1), DCs is the concentration gradient of the solute (1.0 M), Na is the Avogadro number, kB is the Boltzmann constant, T is the temperature (300 K) and DP is the applied hydrodynamic pressure (MPa). The permeability coefficients of the mixed, Mo only, S only and graphene pores were calculated to be 71.64, 83.61, 62.69 and 59.32 #ns-1, respectively. These coefficients are expected to also hold true for small applied pressures (<10 MPa), which are normally used in water desalination, since the relationship between the external

NATURE COMMUNICATIONS | 6:8616 | DOI: 10.1038/ncomms96161

pressure and the rate of water permeation is observed to be quite linear (Fig. 2a). Previous studies40'41 also show that water flux in small nanochannels is linear with respect to external pressure. The permeation rates through various pores (Mo only > mixed > S only > graphene) can also be explained by the energy barrier that a water molecule needs to overcome to enter the pore. These barriers were computed to be DEmo only = 8.50 kBT' DEmixed = 8.84 kBT' DEs only = 9.01 k^T, DEgraphene = 11.05 kBT, which are consistent with the results in Fig. 2a. The details on the energy barrier calculations are documented in Supplementary Fig. 1.

Physical chemistry and geometry of the pore. Water flux (Q) is a function of density (p) inside the pore, velocity (V) of water through the pore and the area of the pore (A), (Q = p ■ V■ A). In water desalination, increasing the area of the pore limits the salt rejection capability of the pore. As the area of the pore increases, the efficiency of rejection decreases25, leaving p and V as the control parameters to increase water flux through the pore.

As shown above, Mo only pore exhibits the highest rate of water permeation. This can be explained by the higher water density (p) and velocity (V) in the Mo only pore compared with those of the Sonly and mixed pores (Fig. 3a-c). The average density of water follows the order of Mo only > mixed > S only (1.47, 1.37 and 1.31 gcm"3, respectively). The denser packing of water molecules at the Mo only pore can be attributed to the hydrophilic nature of Mo sites42 at the edge of the nanopore, which attracts water molecules to the pore interior. It has been shown that the molybdenum surface has a water contact angle close to 0° (molybdenum is a transition metal with a large atomic diameter)42. Attraction of water molecules towards Mo sites becomes more obvious by comparing the mixed and S only pores densities (Fig. 3a). In the mixed pore, the existence of 50% Mo sites gives rise to higher density in the centre of the pore compared with that of S only pore (Fig. 3a).

Next, we explored the velocity profiles in the pore for all the three different pores. The velocities are also higher in Mo only pores compared with mixed and Sonly pores (Fig. 3c). The average velocity of water is 8.26, 7.53 and 7.51ms"1 for Mo only, mixed and S only pores, respectively.

To shed deeper insight into the physical understanding of why the velocity of Mo only pore is higher compared with mixed and S only pores, we computed velocity profiles at the sites of S and Mo for both pore types of Mo only and S only (Fig. 4a,b). This is achieved by binning both pore types at Mo and S sites and averaging velocity at each point for a large number of sets of simulations. We observed that in the Mo only pore, the velocity is

higher at Mo site compared with the S sites. Unlike Mo only pore, we did not observe the velocities to be higher in Mo site in the S only pore, (Fig. 4a,b) which implies that the arrangement of Mo and S sites matter for velocity profiles (see Supplementary Fig. 2 for more evidence on geometry dependency of the velocity in the pore).

It has been shown that conical nanopores have higher fluxes and permeation rates25,43,44. Many biological nanopores, including aquaporin25,45,46, have an hourglass shape, which facilitates rapid water permeation47. Solid-state nanopores have also been designed for conical/hourglass shape to enhance solute and DNA transport48,49. Here in Mo only pores, due to the fish-bone structure of MoS2 (ref. 9), the pore can be tailored13,27 to an hourglass shape at sub-nanometre length scale (see cartoon representation of comparison between Mo only, S only and graphene pores in Fig. 4c). Mo only pore has a contraction centre with hydrophobic S sites at the entrance and S only pore has an expanding centre (Fig. 4c). Graphene has a flat entrance and exit geometry with a single-atom-type exposure at the pore surface50. Water molecules slip on the hydrophobic edges of S and are attracted by the hydrophilic sites of Mo at the pore centre in Mo only case. This arrangement of hydrophobic and hydrophilic atoms along with the conical shape of the pore enhances the flux of water. Also, the water flux highly correlates with the energy barrier of each pore type. The computed potential of mean force for water molecules in each pore type is the reflection of pore chemistry and geometry. In Mo only pore, the potential of mean force is the lowest because of the conical/ hourglass and the hydrophobic-hydrophilic arrangement of the pore atoms (Supplementary Fig. 1). The fundamental advantage of Mo only pore architecture over other pores is the interplay of geometry and chemistry to produce a higher flux of water.


Ion rejection and water flux are the two important factors defining the effectiveness and performance of a water desalination membrane. In Fig. 4d, ion rejection and water permeation rate are plotted for various nanomembrane materials51 (MFI-type zeolite52, commercial polymeric seawater RO53, brackish RO53, nanofiltration53 and high-flux RO53) including MoS2 and graphene investigated in this work. As shown in Fig. 4d, water permeation rate is theoretically enhanced by five orders of magnitude using MoS2 compared with conventional MFI-type zeolite. Also, there is a 70% improvement in the permeation rate of MoS2 compared with graphene. In the study by Cohen-Tanugi et al.19, the permeation rate for graphene is shown to be higher than the rate we observed for graphene. This is because, in our

0 2 4 6

Distance from the center of the pore (Â)

'« 8 (m

2 -0 -

0 2 4 6

Distance from the center of the pore (Â)

Figure 3 | Water density and velocity profiles. (a) Water density distribution in the radial direction in the mixed, Mo only and S only pores with equivalent pore sizes (mixed, A = 55.45 A2 ; Mo only, A = 56.42 A2; Sonly, A = 57.38 A2 ) at a fixed pressure of 250 MPa. (b) Density map of water distribution in Mo only (i) and Sonly (ii) pores. Blue denotes a zero probability of finding a water molecule and red indicates the highest probability of observing a water molecule. (c) Axial velocity of water molecules in the radial direction for mixed, Mo only and Sonly nanopores.

Mo only

S only

Distance from the center of the pore (Â)

0 2 4 6 8

Distance from the center of the pore (Â)

Mo only MoS2

^ » » > * X%< miw

Graphene Hydrophobic edges

tion 100

| 90 of

n e rc er

1 1 1 1 1 1 • . M -

■ MFI zeolite .

• Seawater RO

▲ Brackish RO _

+ Nanofiltrartion

• High-flux RO ■

• Graphene

1 1 < Mo only MoS2 -i i i i

S only

1E-5 1E-4 1E-3 0.01 0.1 1 10 Water permeation rate (l cm-2 per day per MPa)

Figure 4 | Effect of pore type on water permeation and salt rejection. (a) Axial velocity of water molecules in the radial direction at the location of S and Mo atom layers in the Mo only nanopore of A = 56.42 A2 at 250 MPa. (b) Axial velocity of water molecules in the radial direction at the location of S and Mo atom layers in the Sonly nanopore of A = 57.38 A2 at 250 MPa. (c) Cartoon representation of the pore architecture for Mo only, Sonly and graphene nanopore. (d) Performance of various membranes in terms of their ion rejection and water permeation rate. Water permeation rate is expressed per unit area of the membrane and per unit pressure as l cm " 2 per day per MPa.

simulations, the porosity (the ratio of the pore area to the membrane area) is smaller, which decreases the permeation rate per unit area of the membrane. In this work, the comparison of MoS2 and graphene is performed by keeping all conditions identical in the simulations. Thus, MoS2 is potentially an efficient membrane for water desalination.

We have also investigated the potential performance of other transition metal dichalcogenide (MoSe2, MoTe2, WS2, WSe2 and so on) membranes. It was found that the transition metal atom plays a more important role than the chalcogen atom in desalination. More specifically, varying the Lennard-Jones (LJ) parameters of the chalcogen atom does not lead to a significant change in the ion rejection and water permeation (Supplementary Fig. 3).

In conclusion, we have shown that MoS2 membranes are promising for water purification and salt rejection. Mo only pores perform the best among all possible MoS2 pore architectures. MoS2 nanopores with water accessible pore areas ranging from 20 to 60 A2 strongly reject ions allowing <12% of the ions (depending on pore areas) to pass through the porous membranes even at theoretically high pressures of 350 MPa. The water permeation rates associated with these MoS2 porous membranes are found to be two to five orders of magnitude greater than that of currently used membrane materials (MFI-type zeolite, commercial polymeric seawater RO, brackish RO, nanofiltration and high-flux RO) and 70% better than the graphene nanopore. The fish-bone, hourglass architecture of Mo only pore with special arrangement of hydrophobic edges and hydrophilic centre within 1-nm length enhances water permeation to a large extent compared with its other counterparts.


Molecular dynamics (MD) simulations were performed using the LAMMPS package54. The graphene sheet, which acts as a rigid piston to exert external pressure on saline water, along with the MoS2 sheet, water molecules and ions were created by the Visual Molecular Dynamics55. The saline water box was placed between the graphene and MoS2 sheet and pure water was added on the other side of the MoS2 sheet as shown in Fig. 1. A nanopore was drilled in MoS2 by removing the desired atoms. The accessible pore areas considered range from 20 to 60 A2 (Supplementary Fig. 4 for details on pore area calculations). The system dimensions are 4 x 4 x 13 nm in x, y and z, respectively. The box contains B16,000 atoms and the ions (sodium and chloride) have a molarity of B1.0, which is higher than the usual salinity of seawater (0.599 M) because of the computational cost associated with low-salinity solutions.

The extended simple point charge water model was used and the SHAKE algorithm was employed to maintain the rigidity of the water molecule. For non-bonded interactions, the mixing rule was used to obtain the LJ parameters except for carbon-water interactions, which were modelled by the force-field parameters given in ref. 50. The LJ parameters are tabulated in Supplementary Table 1. The LJ cutoff distance was 12 A. The long-range electrostatic interactions were calculated by the Particle Particle Particle Mesh56. Periodic boundary conditions were applied in all the three directions.

For each simulation, first, the energy of the system was minimized for 10,000 steps. Next, the system was equilibrated in constant number of particles, pressure and temperature (NPT) ensemble for 1 ns at a pressure of 1 atm and a temperature of 300 K. Graphene and MoS2 atoms were held fixed in space during equilibration and the NPT simulations allow water to reach its equilibrium density (1 gcm _ 3). Then, an additional constant number of particles, volume and temperature (NVT) simulation was performed for 2 ns to further equilibrate the system. Temperature was maintained at 300 K using the Nose-Hoover thermostat with a time constant of 0.1 ps (refs 57,58). Finally, the production non-equilibrium simulations were carried out in NVT ensemble for 10 ns where different external pressures were applied on the rigid graphene sheet (no longer frozen in space) to characterize the water filtration through the MoS2 nanopores (Supplementary Movie 1). In the production runs, the MoS2 atoms were again held fixed in space to study solely the water transport and ion rejection properties of MoS2. To accelerate the MD simulations and gather enough statistics in the 10-ns simulations, high external pressures ranging from 50 to 350 MPa were considered in this work. Trajectories of atoms were collected every

NATURE COMMUNICATIONS | 6:8616 | DOI: 10.1038/ncomms9616 |

picosecond to obtain the results. For accurate velocity calculations, however, the trajectories were dumped every femtosecond and the data were averaged over 25 sets of simulations with different initial thermal velocity distributions.


1. Elimelech, M. & Phillip, W. A. The future of seawater desalination: energy, technology, and the environment. Science 333, 712-717 (2011).

2. Zhao, S. F., Zou, L., Tang, C. Y. Y. & Mulcahy, D. Recent developments in forward osmosis: opportunities and challenges. J. Membr. Sci. 396,1-21 (2012).

3. Shannon, M. A. et al. Science and technology for water purification in the coming decades. Nature 452, 301-310 (2008).

4. Fritzmann, C., Lowenberg, J., Wintgens, T. & Melin, T. State-of-the-art of reverse osmosis desalination. Desalination 216, 1-76 (2007).

5. Khawaji, A. D., Kutubkhanah, I. K. & Wie, J. M. Advances in seawater desalination technologies. Desalination 221, 47-69 (2008).

6. Humplik, T. et al. Nanostructured materials for water desalination. Nanotechnology 22, 292001 (2011).

7. Geise, G. M. et al. Water purification by membranes: the role of polymer science. J. Polym. Sci. B Polym. Phys. 48, 1685-1718 (2010).

8. Lee, K. P., Arnot, T. C. & Mattia, D. A review of reverse osmosis membrane materials for desalination-development to date and future potential. J. Membr. Sci. 370, 1-22 (2011).

9. Liu, K., Feng, J. D., Kis, A. & Radenovic, A. Atomically thin molybdenum disulfide nanopores with high sensitivity for DNA translocation. ACS Nano 8, 2504-2511 (2014).

10. Fornasiero, F. et al. Ion exclusion by sub-2-nm carbon nanotube pores. Proc. Natl Acad. Sci. USA 105, 17250-17255 (2008).

11. Suk, M. E. & Aluru, N. R. Water transport through ultrathin graphene. J. Phys. Chem. Lett. 1, 1590-1594 (2010).

12. O'Hern, S. C. et al. Selective molecular transport through intrinsic defects in a single layer of CVD graphene. ACS Nano 6, 10130-10138 (2012).

13. Zhao, Y. D. et al. Two-dimensional material membranes: an emerging platform for controllable mass transport applications. Small 10, 4521-4542 (2014).

14. Zhang, D. S. et al. Enhanced capacitive deionization performance of graphene/ carbon nanotube composites. J. Mater. Chem. 22, 14696-14704 (2012).

15. Mishra, A. K. & Ramaprabhu, S. Functionalized graphene sheets for arsenic removal and desalination of sea water. Desalination 282, 39-45 (2011).

16. Celebi, K. et al. Ultimate permeation across atomically thin porous graphene. Science 344, 289-292 (2014).

17. Joseph, S. & Aluru, N. R. Why are carbon nanotubes fast transporters of water? Nano Lett. 8, 452-458 (2008).

18. Suk, M. E., Raghunathan, A. V. & Aluru, N. R. Fast reverse osmosis using boron nitride and carbon nanotubes. Appl Phys. Lett. 92, 133120 (2008).

19. Cohen-Tanugi, D. & Grossman, J. C. Water desalination across nanoporous graphene. Nano Lett. 12, 3602-3608 (2012).

20. Sint, K., Wang, B. & Kral, P. Selective ion passage through functionalized graphene nanopores. J. Am. Chem. Soc. 130, 16448-16449 (2008).

21. Tang, Q., Zhou, Z. & Chen, Z. F. Graphene-related nanomaterials: tuning properties by functionalization. Nanoscale 5, 4541-4583 (2013).

22. Farimani, A. B., Min, K. & Aluru, N. R. DNA base detection using a single-layer MoS2. ACS Nano 8, 7914-7922 (2014).

23. Radisavljevic, B., Radenovic, A., Brivio, J., Giacometti, V. & Kis, A. Single-layer MoS2 transistors. Nat. Nanotechnol. 6, 147-150 (2011).

24. Bertolazzi, S., Brivio, J. & Kis, A. Stretching and breaking of ultrathin MoS2. ACS Nano 5, 9703-9709 (2011).

25. Gravelle, S. et al. Optimizing water permeability through the hourglass shape of aquaporins. Proc. Natl Acad. Sci. USA 110, 16367-16372 (2013).

26. Waduge, P. et al. Direct and scalable deposition of atomically thin low-noise MoS2 membranes on apertures. ACS Nano 9, 7352-7359 (2015).

27. Feng, J. et al. Electrochemical reaction in single layer MoS2: nanopores opened atom by atom. Nano Lett. 15, 3431-3438 (2015).

28. Lee, Y. H. et al. Synthesis of large-area MoS2 atomic layers with chemical vapor deposition. Adv. Mater. 24, 2320-2325 (2012).

29. Dumcenco, D. et al. Large-area epitaxial monolayer MoS2. ACS Nano 9, 4611-4620 (2015).

30. Coleman, J. N. et al. Two-dimensional nanosheets produced by liquid exfoliation of layered materials. Science 331, 568-571 (2011).

31. Smith, R. J. et al. Large-scale exfoliation of inorganic layered compounds in aqueous surfactant solutions. Adv. Mater. 23, 3944-3948 (2011).

32. Liu, K. K. et al. Growth of large-area and highly crystalline MoS2 thin layers on insulating substrates. Nano Lett. 12, 1538-1544 (2012).

33. Zhan, Y. J., Liu, Z., Najmaei, S., Ajayan, P. M. & Lou, J. Large-area vapor-phase growth and characterization of MoS2 atomic layers on a SiO2 substrate. Small 8, 966-971 (2012).

34. Jeon, J. et al. Layer-controlled CVD growth of large-area two-dimensional MoS2 films. Nanoscale 7, 1688-1695 (2015).

35. Mann, J. et al. 2-dimensional transition metal dichalcogenides with tunable direct band gaps: MoS2(1-x)Se2x monolayers. Adv. Mater. 26, 1399-1404 (2014).

36. van der Zande, A. M. et al. Grains and grain boundaries in highly crystalline monolayer molybdenum disulphide. Nat. Mater. 12, 554-561 (2013).

37. Yu, Y. F. et al. Controlled scalable synthesis of uniform, high-quality monolayer and few-layer MoS2 films. Sci. Rep. 3, 6 (2013).

38. Mizan, T. I., Savage, P. E. & Ziff, R. M. Molecular-dynamics of supercritical water using a flexible spc model. J. Phys. Chem. 98, 13067-13076 (1994).

39. Rosenberg, P. A. & Finkelstein, A. Water permeability of gramicidin-a-treated lipid bilayer membranes. J. Gen. Physiol. 72, 341-350 (1978).

40. Suk, M. E. & Aluru, N. R. Effect of induced electric field on single-file reverse osmosis. Phys. Chem. Chem. Phys. 11, 8614-8619 (2009).

41. Zhu, F. Q., Tajkhorshid, E. & Schulten, K. Pressure-induced water transport in membrane channels studied by molecular dynamics. Biophys. J. 83, 154-160 (2002).

42. Liu, Y. H. et al. The preparation of a strawberry-like super-hydrophilic surface on the molybdenum substrate. Colloid Surf. A 404, 52-55 (2012).

43. Gravelle, S., Joly, L., Ybert, C. & Bocquet, L. Large permeabilities of hourglass nanopores: from hydrodynamics to single file transport. J. Chem. Phys. 141, 18C526 (2014).

44. Bocquet, L. & Tabeling, P. Physics and technological aspects of nanofluidics. Lab Chip 14, 3143-3158 (2014).

45. Farimani, A. B., Aluru, N. R. & Tajkhorshid, E. Thermodynamic insight into spontaneous hydration and rapid water permeation in aquaporins. Appl. Phys. Lett. 105, 83702-83702 (2014).

46. Jung, J. S., Preston, G. M., Smith, B. L., Guggino, W. B. & Agre, P. Molecular-structure of the water channel through aquaporin chip - the hourglass model. J. Biol. Chem. 269, 14648-14654 (1994).

47. Hou, X., Guo, W. & Jiang, L. Biomimetic smart nanopores and nanochannels. Chem. Soc. Rev. 40, 2385-2401 (2011).

48. Dekker, C. Solid-state nanopores. Nat. Nanotechnol. 2, 209-215 (2007).

49. Heins, E. A., Siwy, Z. S., Baker, L. A. & Martin, C. R. Detecting single porphyrin molecules in a conically shaped synthetic nanopore. Nano Lett. 5, 1824-1829 (2005).

50. Wu, Y. B. & Aluru, N. R. Graphitic carbon-water nonbonded interaction parameters. J. Phys. Chem. B 117, 8802-8813 (2013).

51. Pendergast, M. M. & Hoek, E. M. V. A review of water treatment membrane nanotechnologies. Energy Environ. Sci. 4, 1946-1971 (2011).

52. Li, L. X., Dong, J. H., Nenoff, T. M. & Lee, R. Desalination by reverse osmosis using MFI zeolite membranes. J. Membr. Sci. 243, 401-404 (2004).

53. Guillen, G. & Hoek, E. M. V. Modeling the impacts of feed spacer geometry on reverse osmosis and nanofiltration processes. Chem. Eng. J. 149, 221-231 (2009).

54. Plimpton, S. Fast parallel algorithms for short-range molecular-dynamics. J. Comput. Phys. 117, 1-19 (1995).

55. Humphrey, W., Dalke, A. & Schulten, K. VMD: visual molecular dynamics. J. Mol. Graph. 14, 33-38 (1996).

56. Hockney, R. W. & Eastwood, J. W. Computer Simulation Using Particles (Taylor & Francis, Inc., 1988).

57. Nose, S. A unified formulation of the constant temperature molecular-dynamics methods. J. Chem. Phys. 81, 511-519 (1984).

58. Hoover, W. G. Canonical dynamics-equilibrium phase-space distributions. Phys. Rev. A 31, 1695-1697 (1985).


This work is supported by AFOSR under grant # FA9550-12-1-0464, NSF under grants 1264282, 1420882 and 1506619, and the Beckman Institute for Advanced Science and Technology at the University of Illinois. We acknowledge the use of the parallel computing resource Blue Waters provided by the University of Illinois and National Center for Supercomputing Applications.

Author contributions

M.H. and A.B.F. carried out the research and wrote the paper. N.R.A. supervised the work.

Additional information

Supplementary Information accompanies this paper at naturecommunications

Competing financial interests: The authors declare no competing financial interests.

Reprints and permission information is available online at reprintsandpermissions/

How to cite this article: Heiranian, M. et al. Water desalination with a single-layer MoS2 nanopore. Nat. Commun. 6:8616 doi: 10.1038/ncomms9616 (2015).

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit