1887
banner image
No data available.
Please log in to see this content.
You have no subscription access to this content.
No metrics data to plot.
The attempt to load metrics for this article has failed.
The attempt to plot a graph for these metrics has failed.
Nonequilibrium molecular dynamics simulation of pressure-driven water transport through modified CNT membranes
Rent:
Rent this article for
USD
10.1063/1.4794685
/content/aip/journal/jcp/138/12/10.1063/1.4794685
http://aip.metastore.ingenta.com/content/aip/journal/jcp/138/12/10.1063/1.4794685

Figures

Image of FIG. 1.
FIG. 1.

Perspective snapshots of the simulation systems produced by VMD package at the beginning status of MD simulations (adapted from the figure in Ref. 12 ). Shown is the beginning empty CNT membrane model connecting two liquid filled reservoirs, and the two graphene sheets acted as the movable walls where the force f t or f b is applied on each carbon atom of the top or bottom wall. Carbon atoms in green, hydrogen atoms in white, and oxygen atoms in red.

Image of FIG. 2.
FIG. 2.

Bar chart of simulation and theoretical flow rates (molecules/ns) of unmodified and modified CNT membranes. The results are based on membrane simulations using the (12, 12) CNT, L z : 6.0 nm; ΔP: 8.0 MPa; T: 300 K. The statistical error for each flow rate data is less than 0.5%. H-P equation refers to the Hagen-Poiseuille equation (Eq. (5) ), and the values of the well-depth parameters (ɛ (CC), ɛ (HH), and ɛ (OO)) are shown in Table I .

Image of FIG. 3.
FIG. 3.

Simulation flow rate (molecules/ns) as a function of the square root of the well-depth parameter. The data points are based on membrane simulations using the (12, 12) CNT, L z : 6.0 nm; ΔP: 8.0 MPa; T: 300 K. The statistical error for each data point is less than 0.5%. The solid line is the linear trend based on the simulation results.

Image of FIG. 4.
FIG. 4.

Effect of van der Waals interactions between water and membrane, represented by the well-depth parameter of the Lennard-Jones potential (interaction increasing from red ♦, to black ■, to green ▲), on the density distributions along the radial direction: density values are averaged over thin annular sections in the pore at the indicated r values over the NEMD simulation at steady state water flow. The data points are based on membrane simulations using the (12, 12) CNT, L z : 6.0 nm; ΔP: 8.0 MPa; T: 300 K. The statistical error for each data point is less than 1%. The values of the well-depth parameters (ɛ (CC), ɛ (HH), and ɛ (OO)) are shown in Table I . The solid lines are the trend curves based on the simulation results, and the dotted lines are the pore surface represented by the effective pore radius (0.643 nm).

Image of FIG. 5.
FIG. 5.

Effect of van der Waals interactions between water and membrane, represented by the well-depth parameter of the Lennard-Jones potential (interaction increasing from red ♦, to black ■, to green ▲), on the density distributions along the flow direction: density values are averaged over thin sections about the indicated z values (cylindrical sections in the pore and cuboid sections in the water reservoirs) over the NEMD simulation at steady state water flow. The data points are based on membrane simulations using the (12, 12) CNT, L z : 6.0 nm; ΔP: 8.0 MPa; T: 300 K. The statistical error for each data point is less than 1%. The values of the well-depth parameters (ɛ (CC), ɛ (HH), and ɛ (OO)) are shown in Table I . The solid lines are the trend curves based on the simulation results. The dashed lines are the membrane boundaries along z direction separating the system in three parts from left to right: the water reservoir on the high pressure side, the pore of the CNT membrane, and the water reservoir on the low pressure side.

Image of FIG. 6.
FIG. 6.

Effect of van der Waals interactions between water and membrane, represented by the well-depth parameter of the Lennard-Jones potential (interaction increasing from red ♦, to black ■, to green ▲), on the velocity distributions (in z direction) as a function of radial position: velocity (in z direction) values are averaged over thin annular sections in the pore at the indicated r values over the NEMD simulation at steady state water flow. The data points are based on membrane simulations using the (12, 12) CNT, L z : 6.0 nm; ΔP: 8.0 MPa; T: 300 K. The statistical errors for most data points are less than 10%. N-S equation refers to the Navier-Stokes equation (Eq. (6) ), and the values of the well-depth parameters (ɛ (CC), ɛ (HH), and ɛ (OO)) are shown in Table I . The solid lines are the trend curves based on the simulation results, the dashed lines are the trends based on the Navier-Stokes equation, and the dotted lines are the pore surface represented by the effective pore radius (0.643 nm).

Image of FIG. 7.
FIG. 7.

Effect of van der Waals interactions between water and membrane, represented by the well-depth parameter of the Lennard-Jones potential (interaction increasing from red ♦, to black ■, to green ▲), on the velocity distributions (in z direction) along the flow direction: velocity (in z direction) values are averaged over thin sections about the indicated z values (cylindrical sections in the pore and cuboid sections in the water reservoirs) over the NEMD simulation at steady state water flow. The data points are based on membrane simulations using the (12, 12) CNT, L z : 6.0 nm; ΔP: 8.0 MPa; T: 300 K. The statistical error for each data point is less than 10%. The values of the well-depth parameters (ɛ (CC), ɛ (HH), and ɛ (OO)) are shown in Table I . The solid lines are the trend curves based on the simulation results. The dashed lines are the membrane boundaries along the z direction separating the system in three parts from left to right: the water reservoir at high pressure side, the pore of the CNT membrane, and the water reservoir at low pressure side.

Image of FIG. 8.
FIG. 8.

Charge patterns of the unmodified and modified CNTs: (a) unmodified CNT, (b) unwrapped CNT of the ring polarized model, and (c) unwrapped CNT of the band polarized model. This figure illustrates the relative locations of the neutral atoms in green, the positive atoms in orange, and the negative atoms in blue. Transport is from left to right as indicated by the arrows on the z-axes.

Image of FIG. 9.
FIG. 9.

Bar chart of simulation and theoretical flow rates (molecules/ns) of unmodified and modified polarized CNT membranes. The results are based on membrane simulations using the (12, 12) CNT, L z : 6.0 nm; ΔP: 8.0 MPa; T: 300 K. The statistical error for each flow rate data is less than 0.5%. H-P equation refers to the Hagen-Poiseuille equation (Eq. (5) ), and the charge patterns of the modified polarized CNTs are shown in Figs. 8(b) and 8(c) .

Image of FIG. 10.
FIG. 10.

Effect of electrostatic interactions between water and membrane, caused by the charge patterns within the CNT, on the density distributions along the radial direction: density values are averaged over thin annular sections in the pore at the indicated r values over the NEMD simulation at steady state water flow. The data points are based on membrane simulations using the (12, 12) CNT, L z : 6.0 nm; ΔP: 8.0 MPa; T: 300 K. The statistical error for each data point is less than 1%. The charge patterns of the modified polarized CNTs are shown in Figs. 8(b) and 8(c) . The solid lines are the trend curves based on the simulation results, and the dotted lines are the pore surface represented by the effective pore radius (0.643 nm).

Image of FIG. 11.
FIG. 11.

Effect of electrostatic interactions between water and membrane, caused by the charge patterns within the CNT, on the density distributions along the flow direction: density values are averaged over thin sections about the indicated z values (cylindrical sections in the pore and cuboid sections in the water reservoirs) over the NEMD simulation at steady state water flow. The data points are based on membrane simulations using the (12, 12) CNT, L z : 6.0 nm; ΔP: 8.0 MPa; T: 300 K. The statistical error for each data point is less than 1%. The charge patterns of the modified polarized CNTs are shown in Figs. 8(b) and 8(c) . The solid lines are the trend curves based on the simulation results. The dashed lines are the membrane boundaries along the z direction separating the system in three parts from left to right: the water reservoir at high pressure side, the pore of the CNT membrane, and the water reservoir at low pressure side.

Image of FIG. 12.
FIG. 12.

Effect of electrostatic interactions between water and membrane, caused by the charge patterns of the CNT, on the velocity (in z direction) distributions along the radial direction: velocity (in z direction) values are averaged over thin annular sections in the pore at the indicated r values over the NEMD simulation at steady state water flow. The data points are based on membrane simulations using the (12, 12) CNT, L z : 6.0 nm; ΔP: 8.0 MPa; T: 300 K. The statistical errors for most data points are less than 10%. N-S equation refers to the Navier-Stokes equation (Eq. (6) ), and the charge patterns of the modified polarized CNTs are shown in Figs. 8(b) and 8(c) . The solid lines are the trend curves based on the simulation results, the dashed lines are the trends based on the Navier-Stokes equation, and the dotted lines are the pore surface represented by the effective pore radius (0.643 nm).

Image of FIG. 13.
FIG. 13.

Effect of electrostatic interactions between water and membrane, caused by the charge patterns of the CNT, on the velocity (in z direction) distributions along the flow direction: velocity (in z direction) values are averaged over thin sections about the indicated z values (cylindrical sections in the pore and cuboid sections in the water reservoirs) over the NEMD simulation at steady state water flow. The data points are based on membrane simulations using the (12, 12) CNT, L z : 6.0 nm; ΔP: 8.0 MPa; T: 300 K. The statistical error for each data point is less than 10%. The charge patterns of the modified polarized CNTs are shown in Figs. 8(b) and 8(c) . The solid lines are the trend curves based on the simulation results. The dashed lines are the membrane boundaries along the z direction separating the system in three parts from left to right: the water reservoir at high pressure side, the pore of the CNT membrane, and the water reservoir at low pressure side.

Image of FIG. 14.
FIG. 14.

Orthographic snapshots of water in the (a) unmodified CNT, (b) ring polarized CNT, and (c) band polarized CNT. Each snapshot is a cross-sectional view in the center section (facing z = −1.0 nm, from z = −1.0 nm to 1.0 nm; see Fig. 8 for the coordinates). Neutral carbon atoms in green, positively charged carbon atoms in orange, negatively charged carbon atoms in blue, hydrogen atoms in white, and oxygen atoms in red. For the unmodified (neutral) wall in (a) the water orientation near the pore surface is random. However, for the polarized cases (b) and (c) there is preferential orientation of water in the CNTs.

Tables

Generic image for table
Table I.

The Lennard-Jones parameters of the unmodified and modified CNT membranes.

Loading

Article metrics loading...

/content/aip/journal/jcp/138/12/10.1063/1.4794685
2013-03-22
2014-04-16
Loading

Full text loading...

This is a required field
Please enter a valid email address
752b84549af89a08dbdd7fdb8b9568b5 journal.articlezxybnytfddd
Scitation: Nonequilibrium molecular dynamics simulation of pressure-driven water transport through modified CNT membranes
http://aip.metastore.ingenta.com/content/aip/journal/jcp/138/12/10.1063/1.4794685
10.1063/1.4794685
SEARCH_EXPAND_ITEM