N2 Separation Membranes

*To whom correspondence should be addressed. Email: [email protected]. Page 1 of 30. ACS Paragon Plus Environment. The Journal of Physical Chemistry. 1...
0 downloads 0 Views 2MB Size
Subscriber access provided by UNIVERSITY OF THE SUNSHINE COAST

Article 2

2

Theoretical Design of Highly-Efficient CO/N Separation Membranes Based on Electric Quadrupole Distinction Yuanyuan Qu, Feng Li, and Mingwen Zhao J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.7b04921 • Publication Date (Web): 28 Jul 2017 Downloaded from http://pubs.acs.org on July 31, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry C is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Theoretical Design of Highly-Efficient CO2/N2 Separation Membranes Based on Electric Quadrupole Distinction Yuanyuan Qua,#,Feng Li a,b, #, Mingwen Zhaoa,* a

School of Physics, Shandong University, Jinan, 250100, Shandong, China

b

School of Physics and Technology, University of Jinan, Jinan 250022, Shandong, China

#

These authors contributed equally.

Corresponding Author *To whom correspondence should be addressed. Email: [email protected]

ACS Paragon Plus Environment

1

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 30

ABSTRACT

Membrane separation of CO2/N2 in fossil fuel gas is promising for the control of greenhouse gas emission, but challenging due to close kinetic diameters. Here, we propose a generalized model for the design of efficient CO2/N2 separation membranes by taking advantage of the large difference between the electric quadrupole moments of the two molecules. The interaction between the molecular electric quadrupole moment and the built-in electric field of the membrane leads to high CO2/N2 selectivity. We validate this model in five nitrogen-rich membranes: g-C3N4, g-C3N3, C2N-h2D, g-C12N8 and p-BN, and demonstrate via molecular dynamics simulations that highly efficient CO2/N2 separation can be achieved in the theoretically predicted g-C12N8 membrane with a permeance of 2.8×105 GPU. This work offers a guidance to improve the separation efficiency of molecules with distinct electric quadrupole moments.

ACS Paragon Plus Environment

2

Page 3 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Introduction Separating CO2 from fossil fuel gas and post-combustion capture of CO2 are effective strategies for CO2 emission reduction. Traditional sorbent-based methods for CO2 capture employ amine-containing molecules frameworks

4,5

and porous carbons

6-8

1

, metal-organic frameworks

2,3

, covalent-organic

, however, face several drawbacks such as low energy

efficiency, equipment corrosion, regeneration of sorbents and toxicity 9. Compared to traditional gas separation methods, the membrane-based technology has attracted considerable attentions due to its low energy cost, small footprint and environmental friendly 10,11

. Graphene is an ideal candidate for gas separation due to its one-atom thickness and

robustness. However, both theoretical and experimental approaches indicate that a perfect graphene sheet is impermeable to any gas molecules

12,13

. Therefore, nanoporous graphene was

theoretically proposed for gas separation, e.g. He, H2, N2, CH4

14-18

. In the past few years,

experiments have demonstrated the potential applications of nanoporous graphene membranes for gas separation utility

19-24

. Inspired by the success of the nanoporous grahpene, the gas

separation performance of other one-atom-thin membranes with intrinsic pores are investigated, such as graphenylene 25, polyphenylene 26,27, graphdiyne 28-32, graphitic carbon nitride 11,27,33,34. In recent years, ultrathin membranes have been proposed theoretically for CO2 separation, such as expanded porphyrins

35

, porous graphene

36,37

, as well as graphenylene

25

. These membranes

with size-specific pores serving as selective channels were employed for CO2 separation through size sieving effect. However, separating CO2 from N2, the main component of fuel gas mixture, is challenging in the size-sieving membranes due to their very close kinetic diameters (3.3-3.9 Å

ACS Paragon Plus Environment

3

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 30

for CO2 and 3.64-3.8 Å for N2) 38. Searching for new separation mechanism is therefore highly desirable. It is noteworthy that CO2 and N2 have distinct electric quadrupole moments which are 14.27×10-40 C•m2 (CO2) and -4.65×10-40 C•m2 (N2)

39

. Under a non-uniform electric field, the

contribution of quadrupole moments will significantly affect the separation process. This huge quadropole moment distinction has been demonstrated to enhance the CO2/N2 separation efficiency both by theoretical and experimental studies

6,7,40-44

. Previous works have indicated

that adding Lewis base group (e.g. N, O) or Lewis acid group (e.g. H) to the pores of the membrane can effectively tune the energy barrier for CO2 diffusion and improve the CO2/N2 selectivity

35,36,45

. However, the quadrupole-moment-driven separation mechanism has not been

well understood, which hinders the improvement of CO2/N2 separation efficiency. In this work, we propose a simple model of quadrupole-moment-driven separation process to guide the design of highly efficient CO2/N2 separation membranes. We demonstrate that the interaction between the build-in electric field around the pores and the quadrupole moments of gas molecules plays a prominent role in the penetration process, which leads to efficient separation. We propose a factor (Q) to evaluate the CO2/N2 separation ability of a membrane. This model is validated in five nitrogen-rich membranes: g-C3N4 BN

46

44

, g-C3N3

44

, C2N-h2D

45

, p-

, g-C12N8. By means of molecular dynamics simulations, we demonstrate highly efficient

CO2/N2 separation in the g-C12N8 membrane with a permeance of 2.8×105 GPU. This work also offers a guidance to increase the separation efficiency of gas molecules with distinct electric quadrupole moments.

ACS Paragon Plus Environment

4

Page 5 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Theoretical Details Our first-principles calculations were performed within the density functional theory (DFT) using the plane-wave pseudo potential approach as implemented in the Vienna Ab initio Simulation Package (VASP) 49-51. The generalized gradient approximation (GGA) in the form of the Perdew-Burke-Ernzerhof (PBE) was adopted for the exchange-correlation functional 52. The van der Waals (vdW) interactions were included explicitly by using a corrected PBE functional (DFT-D2) 53. The electron wavefunctions were expanded using the plane-waves with the energy cutoff of 500 eV. Structural optimizations were carried out using a conjugate gradient (CG) method until the remaining force on each atom was less than 0.01 eV/Å. Vacuum space larger than 30 Å was used to avoid the interaction between adjacent images. The Monkhorst-Pack meshes of 5×5×1 were used in sampling the Brillouin zone for the 2×2 supercells of all membranes

54

. For the transition state calculations, we have performed minimum energy path

profiling using the climbing image nudged elastic band method (CNEB) as implemented in the VASP transition state tools

55,56

. The structural convergence criteria were similar to that used in

the above-mentioned structural optimization. In the subsequent calculations of the potential energy profiles for gas molecules penetrating the membrane pores, the nearest-neighbor atoms surrounding the penetrated pores were fully relaxed while the z-coordinates of the other atoms were kept fixed. The molecular dynamics simulations was performed by GROMACS package

57

in the NVT

ensemble with periodic boundary conditions in all directions. The z-coordinates (perpendicular to the membrane) of the membrane are fixed. The van der Waals parameters for the membrane were taken from the CHARMM27 force field (Table S2). Models of three partial charges were adopted for CO2 and N2 molecules

36,58,59

, where the three-charge-site N2 model has two

ACS Paragon Plus Environment

5

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 30

Lennard-Jones interaction nitrogen sites and a massless point charge located in the middle of the two nitrogen, accounting for the correct quadrupole moment (Table S3). The atomic charges of the membranes were obtained from the ESP fitting according to the Merz–Kollman scheme by Gaussian 09 at the B3LYP/6-31g(d) level

60

using a fragmental cluster of each membrane (Fig.

S1 & Table S1).

Results and Discussion A theoretical model for gas separation We separate the interaction energy W between the membrane and a gas molecule with electric quadrupole moment into two parts: Van der Waals (vdw) interaction energy Wvdw and electrostatic interaction energy Wei due to the quadrupole moment. The electrostatic interaction t r energy for a quadrupole moment D in an external electric field E is

r 1 t Wei = − D : ∇E 6

(1).

Since the membrane used for gas separation usually has symmetrical atomic arrangement around the pore center in the membrane plane (xy plane), the electrostatic interaction energy at the pore center for a linear molecule (e. g. CO2 and N2) can be simplified to 1 ∂E Wei = − Dzz z 6 ∂z

(2). x = y = z =0

Here we focus on the interaction energy at the pore center mainly due to two facts: 1) the electrostatic interaction energy reaches its maximum at the pore center; 2) the transition state of a gas molecule penetrating the pore of a membrane which corresponds to the energy barrier

ACS Paragon Plus Environment

6

Page 7 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

usually locates in the pore center. Considering the spatial distribution of the electric quadrupole moment near the pore center, we introduce a parameter α to normalize the electric quadrupole moment Dzz, and (2) becomes Wei = −

α 6

D zz

∂E z ∂z

. The electrical field in the pore center can x= y = z =0

be regarded as the sum of the electrical field of the charged atoms surrounding the pore, therefore, we should have

∂E z ∂z

where k =

1 4πε 0

= k∑ x = y = z =0

i

qi ri3

(3),

, qi represents the atomic charge (positive or negative) for each atom, ri

represents the distance between each atom and the pore center (please see detailed derivation in SI). Here, we define a separation factor of a membrane as Q = ∑ i

qi . A simple illustration of the ri3

model is shown in Fig. 1, where the charged atoms are distributed symmetrically around the pore center. The electrostatic interaction energy at the pore center then can be written as Wei = −

kα Dzz Q . 6

Figure 1. The illustration of the model. The pore is formed by positively or negatively charged atoms.

ACS Paragon Plus Environment

7

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 30

Supposing that the two molecules with similar kinetic diameters (e.g. CO2 and N2) experience similar vdw energy profiles when penetrating the pore, the electrostatic interactions arising from the quadrupole moments of gas molecules determine their energy barrier difference. The energy barrier difference at the pore center ∆W can be approximated replaced by their electrostatic energy difference ∆Wei . We have

∆W = −

αk 6

Q∆D

(4),

where ∆D = D1 - D2 is the quadrupole moment difference between these two linear molecules. The selectivity of a membrane to the two molecules can be evaluated by S = e − ∆W / RT , according to Arrhenius equation 61,62, as the pore center corresponds to the transition state. In the following parts, we take the CO2/N2 separation as a model system to validate our model, as well as to determine the parameter α.

Five membranes for CO2/ N2 separation According to the above model, highly-efficient CO2/N2 separation should have a large Q value. Graphene-like carbon nitride materials have natural porous configurations suitable for gas separation 11,27,34. The spiral distribution of the charged N and C atoms around the pores leads to build-in electric field which facilitates CO2/N2 separation. Therefore, four graphene-like carbon nitride membranes, namely g-C3N4, g-C3N3, C2N-h2D, g-C12N8, and a porous boron nitride (pBN) as shown in Fig. 2 are considered, which are expected to have large Q values due to the polar nature. The g-C3N4, g-C3N3, C2N-h2D membranes used for gas separation (He or H2) have been extensively investigated theoretically in previous works

11,34,63-65

. The optimized lattice

ACS Paragon Plus Environment

8

Page 9 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

constants are 7.13 Å for g-C3N4 lattice, 7.12 Å for g-C3N3 lattice, 8.33 Å for C2N-h2D lattice, 8.35 Å for g-C12N8 lattice and 6.88 Å for p-BN lattice respectively, consistent with previous literatures 11,34,48,64. The sizes of the pores are indicated by circles in Fig. 2. The atomic charges qi and the distance between atoms and the pore center ri for each membrane can be found in Table S1. The separation factor Q, together with the energy and the energy differences for CO2 and N2 at the pore center of each membrane are summarized in Table 1. The variation of the energy difference ∆W with the separation factor Q are plotted in Fig. 3. A linear fitting of the data give rise a slope of 10.02 VÅ3, corresponding to α=6.95, with a correlation coefficient of 0.96, indicating the feasibility and validity of our model.

ACS Paragon Plus Environment

9

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 30

Figure 2. Top view of (2×2) supercell of five membranes. The yellow circles indicate the inscribed circles of the pores. The brown balls indicate C atoms, the blue balls indicate N atoms and the green balls indicate B atoms, respectively. (A) g-C3N4 supercell; (B) p-BN supercell; (C) g-C3N3 supercell; (D) C2N-h2D supercell; (E) g-C12N8 supercell.

It is noteworthy that the built-in electric field in g-C3N3, C2N-h2D and g-C12N8 membranes leads to negative (positive) W value for CO2 (N2) and thus decreases (increases) the energy of CO2 (N2) at pore center, which increase the CO2/N2 selectivity. Although g-C4N3 has the largest ∆W value corresponding to highest CO2/N2 selectivity, the large positive WCO2 decreases the CO2 permeance. In view of this, we focus on C2N-h2D and g-C12N8 membranes with negative WCO2 and high ∆W values as suitable candidates for CO2/N2 separation for further investigation.

Table 1. The energy (eV), the energy differences (eV) for CO2 and N2 at the pore center of each membrane and the separation factor Q (eV/Å3) for each membrane.

g-C3N4

p-BN

g-C3N3

C2N-h2D

g-C12N8

WCO2

1.197

0.433

-0.186

-0.243

-0.513

WN 2

2.327

0.719

0.641

0.570

0.017

∆W

-1.13

-0.286

-0.827

-0.813

-0.530

Q

-0.104

-0.018

-0.076

-0.063

-0.051

ACS Paragon Plus Environment

10

Page 11 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 3. (A) Energy barrier difference ∆W (black squares) varies with separation factor Q. The red line is the linear fitting of ∆W verses Q.

Molecular dynamics simulation In order to estimate the CO2 permeance, we performed all-atom molecular dynamics simulation to investigate the CO2/N2 separation efficiency by the g-C12N8 and C2N-h2D membranes. In the simulation system, two sheets of each membrane were placed in the center of a simulation box of dimension of 4.176 nm×4.176 nm×50 nm (g-C12N8) or 4.165nm×3.6 nm×50 nm (C2N-h2D), with a distance of 16 nm. At the initial state, 140 CO2 and 140 N2 molecules are randomly placed in between the two sheets as shown in Fig. 4A. The simulations were carried out at 500 K which is the temperature in fuel gases.

ACS Paragon Plus Environment

11

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 30

Figure 4. (A) The snapshot of the initial configuration where the gas molecules are sandwiched between two g-C12N8 membranes. (B) Pure CO2 permeation through the g-C12N8 membrane at 500 K, where different CO2 curves indicate different trajectories. (C) A snapshot for a CO2 molecule passing through the pore of g-C12N8 membrane.

For the g-C12N8 membrane, a total of 10 ns simulation was performed for each trajectory and the numbers of permeate CO2 and N2 were recorded as shown in Fig. 4B. There is no N2 permeation in all three trajectories (red dashed line in Fig. 4B). There are about 70 CO2

ACS Paragon Plus Environment

12

Page 13 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

molecules passing through the g-C12N8 membrane within 10 ns on average of three trajectories, corresponding to a permeance of 2.8×105 GPU (1GPU=3.35×10-10 mol m-2 s-1 Pa-1). Fig. 4C shows a snapshot of a CO2 molecule passing through the g-C12N8 membrane. In contrary to the case of g-C12N8 membrane, none of the gas molecules (both CO2 and N2) was found to pass through the C2N-h2D membrane within 20 ns for three trajectories, suggesting the low CO2 permeance.

To reveal the origins of the permeance difference of these two membranes, we calculated the energy profiles for CO2 and N2 passing through the g-C12N8 and C2N-h2D membranes respectively, as plotted in Fig. 5A&B, using the CNEB method based on first-principles calculations. The interaction between the built-in electric field of membranes and the quadrupole moments of gas molecules leads to attractive states of CO2 molecule at the pore center for both membranes. However, for the C2N-h2D membrane, there is an energy barrier of ~0.25 eV at 4 Å away from the pore center, which greatly reduces the possibility of CO2 penetration through the pores. The lowest energy state of CO2 in the pore center has the molecular axis perpendicular to the membrane. Therefore, CO2 molecules should align perpendicular to the membrane in the collision processes to increase the penetration possibility. Meanwhile, CO2 molecules should have enough energy to overcome this energy barrier via thermal fluctuation. These two requirements lead to low CO2 permeance through the C2N-h2D membrane.

ACS Paragon Plus Environment

13

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 30

Figure 5. Energy profiles for CO2 and N2 passing through the C2N-h2D membrane (A) and gC12N8 membrane (B) based on CNEB method. Potential of mean force (PMF) of CO2 and N2 passing through the C2N-h2D membrane (C) and g-C12N8 membrane (D) by umbrella sampling. The position of the membrane corresponds to 0 Å.

To simulate the effects of molecule alignment and thermal fluctuation, we calculated the free energy profiles of CO2 and N2 molecules passing through the pores of these two membranes via potential of mean force (PMF) analysis. The PMF curves obtained by umbrella sampling method are plotted in Fig. 5C&D. Obviously, the PMF profiles differ significantly from the potential

ACS Paragon Plus Environment

14

Page 15 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

energy profiles calculated by the CNEB method. Due to the randomness of the molecule alignments in the collision process, the free energy barriers for molecules passing through both membranes are much increased, resulting in energy barriers of ~4.39 kBT ( for g-C12N8) and ~14.95 kBT (for C2N-h2D) at 500 K. It should be noted that an energy barrier of ~ 5 kBT is regarded low enough for a molecule to pass through by thermal fluctuation

66

. Therefore, CO2

can easily pass through the g-C12N8 membrane with high permeance. The C2N-h2D membrane, however, has a high free energy barrier, and thus lower CO2 permeability. Compared with CO2, N2 molecules experience much higher free energy barriers, ~12.75 kBT (for g-C12N8) and ~29.90 kBT (for C2N-h2D), indicating the high CO2/N2 selectivity of the two membranes. According to Arrhenius equation 61,62, the CO2/ N2 selectivity of the g-C12N8 membrane can be estimated to be 4.3×103.

Comparison with previous studies and applications We compared the CO2 separation performance (permeance and selectivity) of g-C12N8 membrane with that in previous works on CO2 separation by a few porous two-dimensional membranes. The results are summarized in Table 2 below. With similar permeance, the selectivity of the g-C12N8 membrane is superior to that of ame-4O, PG-4N4H and PG-ES1 membranes. Although the selectivity of the g-C12N8 membrane is lower than that of ame-2O membrane, the permeance, however, are much higher that that of ame-2O membrane. Moreover, it has been reported that CO2 permeance over 4.0×103 GPU with a CO2/N2 selectivity of 40 would reduce the CO2 capture cost below $15/ton 67. Therefore, the g-C12N8 membrane proposed in this study is a competitive membrane for CO2 separation and holds promising for CO2 capture applications.

ACS Paragon Plus Environment

15

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 30

Table 2. Comparison of CO2 permeance (GPU) and CO2/N2 ideal selectivity between difference ultrathin membranes: g-C12N8 (this study), ame-4O (ref. 35), ame-2O (ref. 35), PG-4N4H (ref. 36) and PG-ES1 (ref. 37). g-C12N8

ame-4O

ame-2O

PG-4N4H

PG-ES1

permeance

2.8×105

7.2×104

8.0×103

2.9×105

3×105

selectivity

4.3×103

2.2×102

1.4×106

3×102

6×101

Conclusions To conclude, we propose a generalized model for design of highly-efficient CO2/N2 separation membranes by taking advantage of the large difference of electric quadrupole moments between these two molecules. We define a Q factor which is related to the spatial distribution of charged atoms around the pores in the membranes to predict the separation efficiency. The interaction between the molecular electric quadrupole moment and the built-in electric field leads to high CO2/N2 selectivity. Using this model, we design five two-dimensional nitrogen-rich porous structures as candidates and demonstrate high CO2/N2 separation efficiency in the g-C12N8 membrane through molecular dynamics simulations. This model is not restricted to CO2/N2 separation, but also applicable for other linear gas molecules with distinct quadrupole moments.

ACS Paragon Plus Environment

16

Page 17 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

ASSOCIATED CONTENT

Supporting Information. Supporting Information included Supplementary Figures , Supplementary Tables and Supplementary Method.

AUTHOR INFORMATION The authors declare no competing financial interests.

ACKNOWLEDGMENT The following financial support is acknowledged: National Natural Science Foundation of China (No. 21433006, 11504204), the 111 project (No. B13029), the Fundamental Research Funds of Shandong University (Grant No. 2015HW012), the Technological Development Program in Shandong Province Education Department (Grant No. J14LJ03). The calculations were performed at the National Super Computing Centre in Jinan.

REFERENCES 1.

Dumée, L.; Scholes, C.; Stevens, G.; Kentish, S. Purification of aqueous amine solvents

used in post combustion CO2 capture: a review. Int. J. Greenh. Gas Con. 2012, 10, 443-455. 2.

Sabouni, R.; Kazemian, H.; Rohani, S. Carbon Dioxide Capturing Technologies: A

review focusing on metal organic framework materials (MOFs). Environ. Sci. Pollut. R. 2014, 21, 5427-5449.

ACS Paragon Plus Environment

17

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

3.

Page 18 of 30

Pimentel, B. R.; Parulkar, A.; Zhou, E. K.; Brunelli, N. A.; Lively, R. P. Zeolitic

imidazolate frameworks: next-generation materials for energy-efficient gas separations. ChemSusChem 2014, 7, 3202-3240.

4.

Furukawa, H.; Yaghi, O. M. Storage of hydrogen, methane, and carbon dioxide in highly

porous covalent organic frameworks for clean energy applications. J. Am. Chem. Soc. 2009, 131, 8875-8883. 5.

Dawson, R.; Cooper, A. I.; Adams, D. J. Nanoporous organic polymer networks. Prog.

Polym. Sci. 2012, 37, 530-563.

6.

Sevilla, M.; Valle-Vigón, P.; Fuertes, A. B. N-doped polypyrrole-based porous carbons

for CO2 capture. Adv. Funct. Mater. 2011, 21, 2781-2787. 7.

Hao, G. P.; Li, W. C.; Qian, D.; Lu, A. H. Rapid synthesis of nitrogen-doped porous

carbon monolith for CO2 capture. Adv. Mater. 2010, 22, 853-857. 8.

Dutta, S.; Bhaumik, A.; Wu, K. C. W. Hierarchically porous carbon derived from

polymers and biomass: effect of interconnected pores on energy applications. Energ. Environ. Sci. 2014, 7, 3574-3592.

9.

Ferey, G.; Serre, C.; Devic, T.; Maurin, G.; Jobic, H.; Llewellyn, P. L.; De Weireld, G.;

Vimont, A.; Daturi, M.; Chang, J. S. Why hybrid porous solids capture greenhouse gases? Chem. Soc. Rev. 2011, 40, 550-562.

10. Peng, Y.; Li, Y.; Ban, Y.; Jin, H.; Jiao, W.; Liu, X.; Yang, W. Metal-organic framework nanosheets as building blocks for molecular sieving membranes. Science 2014, 346, 1356-1359. 11. Qu, Y.; Li, F.; Zhou, H.; Zhao, M. Highly efficient quantum sieving in porous graphenelike carbon nitride for light isotopes separation. Sci. Rep. 2016, 6, 19952.

ACS Paragon Plus Environment

18

Page 19 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

12. Bunch, J. S.; Verbridge, S. S.; Alden, J. S.; van der Zande, A. M.; Parpia, J. M.; Craighead, H. G.; McEuen, P. L. Impermeable atomic membranes from graphene sheets. Nano Lett. 2008, 8, 2458-2462. 13. Leenaerts, O.; Partoens, B.; Peeters, F. M. Graphene: a perfect nanoballoon. Appl. Phys. Lett. 2008, 93, 193107.

14. Jiang, D. E.; Cooper, V. R.; Dai, S. Porous graphene as the ultimate membrane for gas separation. Nano Lett. 2009, 9, 4019-4024. 15. Du, H.; Li, J.; Zhang, J.; Su, G.; Li, X.; Zhao, Y. Separation of hydrogen and nitrogen gases with porous graphene membrane. J. Phys. Chem. C 2011, 115, 23261-23266. 16. Hauser, A. W.; Schwerdtfeger, P. Methane-selective nanoporous graphene membranes for gas purification. Phys. Chem. Chem. Phys. 2012, 14, 13292-13298. 17. Jungthawan, S., Reunchan, P.; Limpijumnong, S. Theoretical study of strained porous graphene structures and their gas separation properties. Carbon 2013, 54, 359-364. 18. Schrier, J. Helium separation using porous graphene membranes. J. Phys. Chem. Lett.

2010, 1, 2284-2287. 19. Koenig, S. P.; Wang, L.; Pellegrino, J.; Bunch, J. S. Selective molecular sieving through porous graphene. Nat Nanotechnol. 2012, 7, 728-732. 20. Kim, H. W.; Yoon, H. W.; Yoon, S. M.; Yoo, B. M.; Ahn, B. K.; Cho, Y. H.; Shin, H. J.; Yang, H.; Paik, U.; Kwon, S.; Choi J. Y.; Park H. B. Selective gas transport through few-layered graphene and graphene oxide membranes. Science 2013, 342, 91-95. 21. Li, H.; Song, Z.; Zhang, X.; Huang, Y.; Li, S.; Mao, Y.; Ploehn, H. J.; Bao, Y.; Yu, M. Ultrathin, molecular-sieving graphene oxide membranes for selective hydrogen separation. Science 2013, 342, 95-98.

ACS Paragon Plus Environment

19

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 30

22. O'Hern, S. C.; Boutilier, M. S.; Idrobo, J. C.; Song, Y.; Kong, J.; Laoui, T.; Atieh, M.; Karnik, R. Selective ionic transport through tunable subnanometer pores in single-layer graphene membranes. Nano Lett. 2014, 14, 1234-1241. 23. Celebi, K.; Buchheim, J.; Wyss, R. M.; Droudian, A.; Gasser, P.; Shorubalko, I.; Kye, J. I.; Lee, C.; Park, H. G. Ultimate permeation across atomically thin porous graphene. Science

2014, 344, 289-292. 24. Hu, S.; Lozada-Hidalgo, M.; Wang, F. C.; Mishchenko, A.; Schedin, F.; Nair, R. R.; Hill, E. W.; Boukhvalov, D. W.; Katsnelson, M. I.; Dryfe, R. A.; Grigorieva, I. V.; Wu, H. A.; Geim, A. K. Proton transport through one-atom-thick crystals. Nature 2014, 516, 227-230. 25. Zhu, L.; Jin, Y.; Xue, Q.; Li, X.; Zheng, H.; Wu, T.; Ling, C. Theoretical study of a tunable and strain-controlled nanoporous graphenylene membrane for multifunctional gas separation. J. Mater. Chem. A 2016, 4, 15015-15021. 26. Lu, R.; Meng, Z.; Rao, D.; Wang, Y.; Shi, Q.; Zhang, Y.; Kan, E.; Xiao, C.; Deng, K. A promising monolayer membrane for oxygen separation from harmful gases: nitrogen-substituted polyphenylene. Nanoscale 2014, 6, 9960-9964. 27. Li, Y.; Zhou, Z.; Shen, P.; Chen, Z. Two-dimensional polyphenylene: experimentally available porous graphene as a hydrogen purification membrane. Chem. Commun. 2010, 46, 3672-3674. 28. Bartolomei, M.; Carmona-Novillo, E.; Hernández, M. I.; Campos-Martínez, J.; Pirani, F.; Giorgi, G. Graphdiyne pores: 'ad hoc' openings for helium separation applications. J. Phys. Chem. C 2014, 118, 29966-29972.

29. Cranford, S. W.; Buehler, M. J. Selective hydrogen purification through graphdiyne under ambient temperature and pressure. Nanoscale 2012, 4, 4587-4593.

ACS Paragon Plus Environment

20

Page 21 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

30. Jiao, Y.; Du, A.; Hankel, M.; Zhu, Z.; Rudolph, V.; Smith, S. C. Graphdiyne: a versatile nanomaterial for electronics and hydrogen purification. Chem. Commun. 2011, 47, 11843-11845. 31. Zhang, H.; Zhao, X.; Zhang, M.; Luo, Y.; Li, G.; Zhao, M. Three-dimensional diffusion of molecular hydrogen in graphdiyne: a first-principles study. J. Physi. D Appl. Phys. 2013, 46, 5307. 32. Zhang, H.; He, X.; Zhao, M.; Zhang, M.; Zhao, L.; Feng, X.; Luo, Y. Tunable hydrogen separation in sp–sp2 hybridized carbon membranes: a first-principles prediction. J. Phys. Chem. C 2012, 116, 16634-16638.

33. Cao, S.; Yu, J. g-C3N4-based photocatalysts for hydrogen generation. J. Phys. Chem. Lett.

2014, 5, 2101-2107. 34. Li, F.; Qu, Y.; Zhao, M. Efficient helium separation of graphitic carbon nitride membrane. Carbon 2015, 95, 51-57.

35. Tian, Z.; Dai, S.; Jiang, D. E. Expanded porphyrins as two-dimensional porous membranes for CO2 separation. ACS Appl. Mater. Inter. 2015, 7, 13073-13079. 36. Liu, H.; Dai, S.; Jiang, D. E. Insights into CO2/N2 separation through nanoporous graphene from molecular dynamics. Nanoscale 2013, 5, 9984-9987. 37. Schrier, J. Carbon dioxide separation with a two-dimensional polymer membrane. ACS Appl. Mater. Inter. 2012, 4, 3745-3752.

38. Mason, J. A.; McDonald, T. M.; Bae, T. H.; Bachman, J. E.; Sumida, K.; Dutton, J. J.; Kaye, S. S.; Long, J. R. Application of a high-throughput analyzer in evaluating solid adsorbents for post-combustion carbon capture via multicomponent adsorption of CO2, N2, and H2O. J. Am. Chem. Soc. 2015, 137, 4787-4803.

ACS Paragon Plus Environment

21

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 30

39. Graham, C.; Imrie, D. A.; Raab, R. E. Measurement of the electric quadrupole moments of CO2, CO, N2, Cl2 and BF3. Mol. Phys. 1998, 93, 49-56. 40. Babarao, R.; Dai, S.; Jiang, D. E. Nitrogen-doped mesoporous carbon for carbon capture – a molecular simulation study. J. Phys. Chem. C 2012, 116, 7106-7110. 41. Wu, T.; Xue, Q.; Ling, C.; Shan, M.; Liu, Z.; Tao, Y.; Li, X. Fluorine-modified porous graphene as membrane for CO2/N2 separation: molecular dynamic and first-principles simulations. J. Phys. Chem. C 2012, 118, 7369–7376. 42. Babarao, R.; Dai, S.; Jiang, D. E. Functionalizing porous aromatic frameworks with polar organic groups for high-capacity and selective CO2 separation: a molecular simulation study. Langmuir 2011, 27, 3451-3460.

43. Babarao, R.; Custelcean, R.; Hay, B. P.; Jiang, D. E. Computer-aided design of interpenetrated tetrahydrofuran-functionalized 3D covalent organic frameworks for CO2 capture. Cryst. Growth Des. 2012, 12, 5349-5356.

44. Fu, J.; Wu, J.; Custelcean, R.; Jiang, D. E. Nitrogen-doped porous aromatic frameworks for enhanced CO2 adsorption. J. Colloid Interf. Sci. 2015, 438, 191-195. 45. Wang, Y.; Yang, Q.; Zhong, C.; Li, J. Graphene-like poly (triazine imide) as N2-selective ultrathin membrane for postcombustion CO2 capture. J. Phys. Chem. C 2016, 120, 28782-28788. 46. Guo, Q.; Yang, Q.; Yi, C.; Zhu, L.; Xie, Y. Synthesis of carbon nitrides with graphitelike or onion-like lamellar structures via a solvent-free route at low temperatures. Carbon 2005, 43, 1386-1391.

47. Mahmood, J.; Lee, E. K.; Jung, M.; Shin, D.; Jeon, I. Y.; Jung, S. M.; Choi, H. J.; Seo, J. M.; Bae, S. Y.; Sohn, S. D. Nitrogenated holey two-dimensional structures. Nat. Commun. 2015, 6, 6486.

ACS Paragon Plus Environment

22

Page 23 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

48. Zhang, H.; Tong, C. J.; Zhang, Y.; Zhang, Y. N.; Liu, L. M. Porous BN for hydrogen generation and storage. J. Mater. Chem. A 2015, 3, 9632-9637. 49. Kohn, W.; Sham, L. J. Self-consistent equations including exchange and correlation effects. Physical Review 1965, 140, 1133-1138. 50. Kresse, G.; Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. J. Phys. Rev. B 1996, 54, 11169-11186. 51. Hohenberg, P.; Kohn, W. Inhomogeneous electron gas. Physical Review 1964, 136, B864-B871. 52. Perdew, J.; Burke, K.; Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 1996, 77, 3865-3868.

53. Grimme, S. Semiempirical GGA-type density functional constructed with a long-range dispersion correction. J. Comput. Chem. 2006, 27, 1787-1799. 54. Chadi, D. J. Special points for brillouin-zone integrations. Phys. Rev. B 1977, 16, 17461747. 55. Henkelman, G.; Uberuaga, B. P.; Jónsson, H. A climbing image nudged elastic band method for finding saddle points and minimum energy paths. J. Chem. Phys. 2000, 113, 99019904. 56. Henkelman, G.; Jónsson, H. Improved tangent estimate in the nudged elastic band method for finding minimum energy paths and saddle points. J. Chem. Phys. 2000, 113, 99789985. 57. Abraham, M. J.; Murtola, T.; Schulz, R.; Páll, S.; Smith, J. C.; Hess, B.; Lindahl, E. GROMACS: high performance molecular simulations through multi-level parallelism from laptops to supercomputers. Softwarex, 2015, 1–2, 19-25.

ACS Paragon Plus Environment

23

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 30

58. Harris, J. G.; Yung, K. H. Carbon dioxides liquid-vapor coexistence curve and critical properties as predicted by a simple molecular-model. J. Phys. Chem. 1995, 99, 12021-12024. 59. Murthy, C. S.; Singer, K.; Klein, M. L.; McDonald, I. R. Pairwise additive effective potentials for nitrogen. Mol. Phys. 1980, 41, 1387-1399. 60. Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Petersson, G. A.; Nakatsuji, H, et al. Gaussian 09, Revision A. 02; Gaussian. Inc.: Wallingford, CT. 2009. 61. Arrhenius, S. Über die dissociationswärme und den einfluss der temperatur auf den dissociationsgrad der elektrolyte. Zeitschrift für physikalische Chemie 1889, 4, 226-248 62. Arrhenius, S. Über die reaktionsgeschwindigkeit bei der inversion von rohrzucker durch säuren. Zeitschrift für physikalische Chemie 1889, 4, 96-116. 63. de Silva, S. W.; Du, A.; Senadeera, W.; Gu, Y. Strained graphitic carbon nitride for hydrogen purification. J. Membrane Sci. 2017, 528, 201-205. 64. Ma, Z.; Zhao, X.; Tang, Q.; Zhou, Z. Computational prediction of experimentally possible g-C3N3 monolayer as hydrogen purification membrane. Int. J. Hydrogen Energ. 2014, 39, 5037-5042.

65. Xu, B.; Xiang, H.; Wei, Q.; Liu, J. Q.; Xia, Y. D.; Yin, J.; Liu, Z. G. Two-dimensional graphene-like C2N: an experimentally zvailable porous membrane for hydrogen purification. Phys. Chem. Chem. Phys. 2015, 17, 15115-15118.

66. Won, C. Y.; Aluru, N. R. Water permeation through a subnanometer boron nitride nanotube. J. Am. Chem. Soc. 2007, 129, 2748-2749. 67. Merkel, T. C.; Lin, H.; Wei, X.; Baker, R. Power plant post-combustion carbon dioxide capture: An opportunity for membranes. J. Membrane Sci. 2010, 359, 126-139.

ACS Paragon Plus Environment

24

Page 25 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

TOC GRAPHIC

ACS Paragon Plus Environment

25

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 1. The illustration of the model. The pore is formed by positively or negatively charged atoms. 68x35mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 26 of 30

Page 27 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 2. Top view of (2×2) supercell of five membranes. The yellow circles indicate the inscribed circles of the pores. The brown balls indicate C atoms, the blue balls indicate N atoms and the green balls indicate B atoms, respectively. (A) g-C3N4 supercell; (B) p-BN supercell; (C) g-C3N3 supercell; (D) C2N-h2D supercell; (E) g-C12N8 supercell. 177x184mm (300 x 300 DPI)

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 3. (A) Energy barrier difference ∆W (black squares) varies with separation factor Q. The red line is the linear fitting of ∆W verses Q. 81x64mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 28 of 30

Page 29 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 4. (A) The snapshot of the initial configuration where the gas molecules are sandwiched between two g-C12N8 membranes. (B) Pure CO2 permeation through the g-C12N8 membrane at 500 K, where different CO2 curves indicate different trajectories. (C) A snapshot for a CO2 molecule passing through the pore of g-C12N8 membrane. 177x164mm (300 x 300 DPI)

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 5. Energy profiles for CO2 and N2 passing through the C2N-h2D membrane (A) and g-C12N8 membrane (B) based on CNEB method. Potential of mean force (PMF) of CO2 and N2 passing through the C2N-h2D membrane (C) and g-C12N8 membrane (D) by umbrella sampling. The position of the membrane corresponds to 0 Å. 177x142mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 30 of 30