Competitive Adsorption of Methanol−Acetone on Surface

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...
1 downloads 0 Views 2MB Size
Subscriber access provided by Macquarie University

Functional Nanostructured Materials (including low-D carbon)

Competitive Adsorption of Methanol-Acetone on Surface Functionalization (-COOH, -OH, -NH2 and -SO3H)#GCMC and DFT Simulations Yang Guo, Zheng Zeng, Liqing Li, Changqing Su, Ruofei Chen, Chunhao Wang, Ke Zhou, Xiang Xu, and Hailong Li ACS Appl. Mater. Interfaces, Just Accepted Manuscript • DOI: 10.1021/acsami.9b10804 • Publication Date (Web): 28 Aug 2019 Downloaded from pubs.acs.org on August 30, 2019

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 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 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.

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 29 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

Competitive

ACS Applied Materials & Interfaces

Adsorption

of

Methanol−Acetone

on

Surface

Functionalization (-COOH, -OH, -NH2 and -SO3H):GCMC and DFT Simulations Yang Guo, Zheng Zeng*, Liqing Li*, Changqing Su, Ruofei Chen, Chunhao Wang, Ke Zhou, Xiang Xu, Hailong Li School of Energy Science and Engineering, Central South University, Changsha 410083, Hunan, China

ABSTRACT: The capture and separation properties of surface functionalized activated carbons (AC-Rs, R= COOH, -OH, -NH2 and -SO3H) for the methanol-acetone mixture were investigated for the first time by grand canonical Monte Carlo (GCMC) simulation and density functional theory (DFT). The effects of surface functional groups and structural characteristics of AC-Rs on the adsorption and separation behaviors of methanol and acetone had been clarified. The surface functional group with strong electron-donating or electron-accepting capacity (i.e., -NH2, -OH, and -SO3H) was a crucial factor for the methanol-acetone capture and separation performance at the lower pressure range, and the accessible surface area was found to be another determinative factor. AC-NH2 with the relatively large accessible surface area (4497 m2/g) exhibited an efficient capture performance for the single component (15.7 mol/kg for methanol and 6.7 mol/kg for acetone) and the highest methanol/acetone selectivity (~23) at 0.02 kPa. While at high pressures, the surface functionalization and available pore volume of AC-Rs played pivotal roles in the adsorptive separation process. This study provided mechanistic insights on how the surface functional groups affected the capture and separation properties of ACs, which would further provide a rational alternative strategy in the preparation and synthesis of ACs for the effective gas mixture separation. KEYWORDS: activated carbons, functional groups, methanol/acetone adsorption, separation, GCMC, DFT 1

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 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 29

1. INTRODUCTION Methanol-acetone mixture, used in both preparative chemistry and chemical industry, is one of the most widely used organic solvents.1 The mixture of these two compounds exhibits easy diffusion as volatile organic compounds (VOCs) into the atmosphere, which may cause a negative effect on the human health and the atmosphere.2-5 Hence, establishing effective strategies in separation and capture of the methanol-acetone mixture will be significant to the recycling and concentration of the solvent. Adsorption has been considered as one of the most favorable capture and separation technologies for its low energetic consumption and high efficiency.6-8 Many investigations on different materials, i.e., activated carbon (AC),9-11 metal oxides, 12-14 and metalorganic frameworks (MOFs),15-17 for VOCs capture have been conducted. Among these adsorbents, AC materials have been regarded as the competitive candidates because of their high specific surface area, variable surface functional groups, eco-friendly property and relatively low-cost.18,19 The VOCs adsorption capacity of AC is determined by its structural characteristics and surface modification.20 Various functional groups have been added to carbon materials to enhance their VOCs capture performance.21-27 However, since it’s difficult to calculate the selectivity of VOCs mixtures by ideal adsorbed solution theory (IAST) for the large nonideality in the adsorbed phase,

28,29

the separation performances of ACs for VOCs were rarely

reported. The use of molecular simulation allows us to have an intuitive understanding of the gas-framework interaction and the adsorptive separation process at a molecular level. Molecular simulations have been widely used to study the adsorption behavior of MOFs,30-32 covalent organic frameworks (COFs),33,34 ACs,35,36 and

2

ACS Paragon Plus Environment

Page 3 of 29 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

ACS Applied Materials & Interfaces

zeolites.37,38 In terms of AC, the slit pore model has been most commonly adopted to study the influence of functional groups and pore size distribution (PSD) on its adsorption capacity. 39-41 However, this simplified model is not indicative of presenting the complex nature of real carbon materials. Recently, a random structure AC model42, consisting of several coronene-shaped graphitic basic structural unit (BSU), was reported to have a good match with the pore topology and morphology of the real carbon. This random structure model has been successfully and widely used to predict the adsorption performance of CO2/CH4 in ACs, and the calculated adsorption amounts are well consistent with the experimental results.43 More importantly, functional groups can be connected to the graphitic BSU to constitute the modified AC frameworks. For example, Lu et al.44 predicted the adsorption and separation performance of a binary CO2/CH4 mixture in edge functionalized ACs, and found that adding appropriate surface functional groups to the ACs could improve the capture and separation performance of the CO2/CH4 mixture. In this work, the random structure AC model was constructed to simulate internal structure of real carbons. The graphitic BSU was intercepted from a sheet of graphene that consists of 24 carbon atoms, and four functional groups (-COOH, -OH, -NH2 and -SO3H) were connected to the BSU to constitute the modified frameworks. Grand canonical Monte Carlo (GCMC)45,46 simulations were used to predict the thermodynamic equilibrium properties of single- and mixture-component of methanol-acetone in ACs. Density functional theory (DFT) was performed to assess the most stable structure of the graphitic BSUs, and to establish the optimal adsorption position of molecules on the surface of ACs. The intrinsic enhancement mechanism in adsorption capacity was elucidated by analyzing the pore topology and electrostatic interaction contribution. The calculated results of potential energies of adsorbate-framework and adsorbate-adsorbate interactions (electrostatic and van der Waals

3

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 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 29

(vdW)) as well as the optimal adsorption position clarified the different adsorption behaviors of methanol and acetone in ACs at both low and high pressures. The adsorption of a binary methanol-acetone mixture in ACs was preformed to investigate the separation performance of the chosen functional groups. This study provided the detailed information on the capture, competitive adsorption and separation of the binary methanol-acetone mixture on the surface functionalized ACs, which paves a new way for designing ACs for the effective gas mixture separation.

2. SIMULATION METHODOLOGY To study the capture and separation properties of surface functionalized activated carbons (AC-Rs, R= COOH, -OH, -NH2 and -SO3H) for methanol-acetone mixtures, functional groups were connected to the graphitic BSU to constitute the modified AC frameworks (see Figure 1).

Figure 1. The structure of graphitic BSU and AC frameworks. All the functionalized BSUs were optimized with the DFT geometry optimization in Gaussian 09 package.47 Mulliken charge was used as the atomic partial charge and the B3LYP functional combined with the 6-31+G(d,p) basis set was used to do the spin-unrestricted all-electron quantum-chemistry calculation.48,49 The optimized BSU 4

ACS Paragon Plus Environment

Page 5 of 29 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

ACS Applied Materials & Interfaces

structure and the atom charge were presented in Figure S1 and Table S1, and the calculated atom charge was further compared with literatures which proved to be reasonable as shown in Table S2. The AC models were then constructed from a collection of geometry optimized BSUs at a mass density of 0.542 g/cm3.42 The visualization and detailed modeling process of the unit cell of the five AC-Rs (R=None, COOH, OH, NH2, SO3H) models were shown in Figure S2-S4. The minimum energy configurations of methanol/acetone molecule adsorbed on the BSUs were calculated by density functional dispersion correction

(DFT-D) in Dmol3 program package.50 Generalized gradient

approximation (GGA) with the Perdew–Burke–Ernzerhof (PBE)51 functional was chosen to predict the adsorbatesurface interaction for its accuracy in describing the exchange–correlation interaction.52,53 The DFT semi-core pseudopotential (DSPP)52 processing and the double numerical including polarization (DNP) basis set were chosen for the atoms. GCMC simulation was carried out to estimate the adsorption isotherms for single or mixed gas components. All the GCMC calculations were performed in the simulation software RASPA.54,55 The GCMC simulations were began with an equilibration run of 1.0×105 cycles, followed by a production run of 5.0×105 cycles. The potential energy of the host-sorbate as well as the sorbate-sorbate system was calculated by the sum of the electrostatic and the vdW energy. The calculation formulas of the vdW and the electrostatic potential energy were presented in the SI. In addition, the Lennard Jones (12-6) model56 with a cutoff radius of 18 Å was established to describe the vdW energy, while the electrostatic energy was calculated by the Ewald summation. The L-J potential parameters, for the AC-Rs, were taken from Dreiding force field57. The L-J potential parameters are provided in Table S3. The L-J potential parameters and partial charges of the methanol and acetone adsorbates were obtained from

5

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 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 29

TraPPE-UA (united atom) force field.58,59 The L-J potential parameters and atomic partial charges of the gas molecule are shown in Table S4.In this force field, the methyl groups in the methanol and acetone molecules were considered to be a united-atom located on a carbon atom, which has been widely adopted in predicting the adsorption behavior in AC23 and MOFs31.

3. RESULTS AND DISCUSSION 3.1. Structural Characteristics of AC-Rs. To investigate how the functionalization affects the pore topology and morphology of AC, the physical characteristics of the five AC models are calculated (Table 1). The structural features of the AC-Rs in Table 1 were calculated in RASPA. And the relevant computational methods and definitions are elaborated in the SI. Compared to AC-None (1.339 cm3/g), the presence of -SO3H (1.457 cm3/g), -COOH (1.396 cm3/g) or -OH (1.379 cm3/g) functional groups leads to the enlargement of the total available pore volume (Vtot). The variation of the maximum pore diameter (Dm) is in accordance with that of Vtot, that is, the AC-SO3H exhibits the largest expansion in Dm values (19.13 Å) compared with AC-None (14.13 Å), then followed by AC-COOH (15.92 Å) and AC-OH (14.66 Å). The above results indicate that all the selected functional groups enlarge the Vtot and Dm of AC except AC-NH2, which has a Vtot of 1.302 cm3/g and a maximum pore diameter of 13.77 Å. Table 1. Structural characteristics of the AC-Rs. -None

-COOH

-OH

-NH2

-SO3H

40

20

40

30

20

Dimensions of the model (Å )

32.81

37.03

39.22

35.35

42.62

Porosity, Фc (%)

72.49

75.66

74.67

70.55

78.97

AC-Rs Number of BSUs 3

6

ACS Paragon Plus Environment

Page 7 of 29 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

ACS Applied Materials & Interfaces

Total available pore volume, Vtot (cm3/g)

1.339

1.396

1.379

1.302

1.457

Maximum pore diameter, Dm(Å)

14.13

15.92

14.66

13.77

19.13

Accessible surface area (m2/g)

4390

4043

4443

4497

3684

The accessible surface area of the framework is another key factor to obtain a higher adsorption capacity, since a larger accessible surface area means more physical adsorption sites.60 The AC-NH2 (4497 m2/g ) and ACOH (4443 m2/g) has the relatively large accessible surface area among the functionalized ACs, which is larger than that of carbon nanotubes (CNTs) (153-1315 m2/g)61 and mesoporous high-surface-area activated carbon (931-2191 m2/g)62. However, the accessible surface areas of other surface functionalized AC models are smaller than that of the AC-None. Hence, the -NH2 functional group has a positive effect on the accessible surface area but a negative effect on the Vtot, to which the -SO3H, -COOH and -OH functional groups are adverse. This is further verified by analyzing the pore size distributions (PSDs) shown in Figure 2. And the visualization of the pore shape and connectivity is presented in Figure S5. The introducing of -NH2 functional group moves the pore sizes of AC-None (3.80-14.13 Å) to a smaller range (3.38-13.77 Å), whereas the other three functional groups lead to an expansion in the PSDs, especially the -SO3H (4.00-19.13 Å). More detailly, the pore size of AC-NH2, with two peaks at ∼7.1 and ∼9.1 Å, is dominantly distributed under 10 Å, which means large proportion of smaller pores. Obviously, the smaller pores have a higher contribution to the accessible surface area but a lower contribution to the Vtot. Instead, when it comes to AC-SO3H, whose pores are dominantly distributed in 10.0018.0 Å, it exhibits the largest Vtot but the smallest accessible surface area among all the ACs. All the PSDs are within a pore size range of 3.38-19.13Å, dominated by supermicropores (within 7.00-20.0 Å) with less proportion of ultramicropores ( AC-SO3H > AC-OH >

8

ACS Paragon Plus Environment

Page 9 of 29 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

ACS Applied Materials & Interfaces

AC-COOH> AC-None at the low pressures. Among the five AC-Rs models, AC-NH2 has the highest uptake (15.7 and 6.7 mol/kg for methanol and acetone, respectively) at 0.02 kPa, while AC-COOH shows a relatively weak adsorption capacity of two adsorbates. Adsorption isotherms also demonstrate that the adsorption capacity for methanol on selected AC-Rs at low pressures were more exothermic than those for acetone except AC-None.

(a)

(b)

Figure 3. Absolute adsorption isotherms of methanol (a) and acetone (b) on AC-Rs. Moreover, the adsorption energies (Eads) between the adsorbent and the adsorbate at 0.02 kPa are shown in Figure 4. The Eads of the two adsorbates on ACs follows the order of AC-SO3H>AC-NH2>AC-OH>ACCOOH>AC-None. The introduction of appropriate functional groups has a positive effect on the interaction between the adsorbates and ACs at ultralow pressures. Note that the selected functional groups enhance the absolute value of Eads from ∼23 to ∼73 kJ/mol for methanol and from ∼29 to ∼63 kJ/mol for acetone. AC-SO3H exhibits the strongest interaction with the two adsorbates (-73.2 kJ/mol for methanol and -63.4 kJ/mol for acetone), but a weaker adsorption capacity than AC-NH2 at 0.02 kPa, since the accessible surface area for providing physical adsorption sites of AC-NH2 (4497 m2/g) is larger than that of AC-SO3H (3684 m2/g). The largest uptake of AC-NH2 at low pressures is a consequence of a high Eads combined with the largest accessible surface area.

9

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 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 29

Figure 4. Adsorption energies for methanol and acetone on AC-Rs. The isotherms of methanol and acetone nearly reach their saturation at 12 and 18 kPa, respectively. Moreover, the largest enhancement of methanol and acetone uptake on AC-SO3H at high pressures can be attributed to its largest total available pore volume (1.457 cm3/g) and its highest Eads. Hence, the adsorption capacity of AC-Rs under high pressures is related to its total pore volume and functionalization property. In terms of isotherm type, the absolute adsorption isotherms of methanol and acetone exhibit type-I Langmuir adsorption behavior,63 except the adsorption isotherm of methanol on the AC-None. Apparently, the adsorption of methanol on the AC-None shows a low adsorption capacity at low pressures but a rapid rise at higher pressures, which is similar to the adsorption behavior of water on the graphite surface and exhibits the type-Ⅲ according to the IUPAC. 3.3. Influences of Electrostatic and Van Der Waals Interactions. Consider that the introduction of surface

10

ACS Paragon Plus Environment

Page 11 of 29 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

ACS Applied Materials & Interfaces

functional groups can simultaneously enhance the adsorption capacity and Eads of methanol and acetone on ACs, the electrostatic interaction (EI) contribution to the adsorption capacity of the two adsorbates is investigated. The EI contribution is defined as Eq. 1. EI contribution =

|𝑄𝑤𝑖𝑡ℎ ― 𝑄𝑤𝑖𝑡ℎ𝑜𝑢𝑡| 𝑄𝑤𝑖𝑡ℎ

(1)

x100%

where 𝑄𝑤𝑖𝑡ℎ and 𝑄𝑤𝑖𝑡ℎ𝑜𝑢𝑡 represent the absolute adsorption capacity with or without the gas-adsorbent EI. The 𝑄𝑤𝑖𝑡ℎ𝑜𝑢𝑡 was calculated by setting the atomic partial charges of the adsorbent (AC-Rs) to zero. Figure 5 shows the EI contributions of gas-framework at different pressure ranges. Surface-functionalized AC-Rs exhibit a larger EI contribution than the AC-None at low pressures, indicating that the functionalization enhances the EI between the adsorbent and the adsorbate. The selected electron-donating groups (-NH2, -OH) and electron-accepting group (-SO3H, -COOH) exhibit strong electrostatic interactions for their donating/accepting electron density capacity (a)

(b)

Figure 5. The gas-framework EI contributions on the uptake of (a) methanol and (b) acetone. Figure 6 exhibits that the N atoms in BSU-NH2 show an electronegative property by gaining electrons from the direct-connected C atoms, while the H atoms in the -NH2 group and the peripheral C in the activated carbon exhibit a high electropositivity by donating electrons. Those electropositive/electronegative atoms on the 11

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 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 29

activated carbon surface become polar adsorption sites for methanol and acetone molecules, which effectively strength the electrostatic interactions between the adsorbate and framework especially at low pressures. Moreover, stronger electropositive/electronegative property of the functional group leads to more effects on the atomic partial charge of activated carbon. As shown in Table S5, the N atom in -NH2 group exhibits stronger electronegativity than the O atoms in -OH and -COOH functional groups, whereas the S atom in the -SO3H functional group shows a high electropositivity.

0.296 -0.711 0.360 0.323 -0.822 -0.302 0.317

0.291 -0.666

0.311 -0.515

0.326 -0.820 0.324

-0.053

-0.508

1.024

0.619

0.061

-0.400

0.350

H N C

0.331 0.302 0.325 -0.699 0.311 -0.720 -0.668 0.335 -0.532

0.070

0.076

0.046

0.118

-0.798 -0.623

0.302

-0.097

0.807

0.030

1.338

-0.426

1.012

-0.348

-0.305

0.327

0.352

-0.248

-0.676 -0.850 0.295 0.357 0.331 0.304

0.285 -0.668 0.322

0.324 -0.672 0.290

The direct-connected C with N atom

The peripheral C

Figure 6. Partial charges for the BSU-NH2. It can also be found that EI contributions of the gas-framework on the adsorption capacity decrease with the rising pressure. The adsorbate gas molecules firstly take up the optimal adsorption sites on the surface of the ACRs because of the strong adsorbate-framework electrostatic interactions at low pressures. However, when these 12

ACS Paragon Plus Environment

Page 13 of 29 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

ACS Applied Materials & Interfaces

physical adsorption sites reach their saturation, the monolayer adsorption is changed into multilayer one64 and the EIs between the adsorbate and adsorbent become less significant to the adsorption. Moreover, EI contributions of gas-framework on the uptake of methanol remain stable at a relatively high level before 1 kPa, and then follow a rapid drop at high pressures, which is much different from that of acetone (decreasing in a stepwise manner). It can be inferred that the methanol with a higher electrostatic affinity is more sensitive to electrostatic interactions at the low pressure. It can be verified from the adsorption isotherms (see Figure 2). The surface functionalized activated carbon exhibits a lager methanol uptake than that of acetone at the low pressures. In addition, to clarify the adsorption mechanisms of the two adsorbates on ACs at both low and high pressures, we take AC and AC-SO3H as examples, and calculate the potential energies of four interaction cases by GCMC (Figure 7): (i) the adsorbate-framework electrostatic interactions (EIs), (ii) the adsorbate-framework van der Waals interactions (VIs), (iii) the adsorbate-adsorbate EIs, and (iv) the adsorbate-adsorbate VIs.

13

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 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

(a)

(b)

(c)

(d)

Page 14 of 29

Figure 7. Potential energies of (a, c) methanol and (b, d) acetone on AC-None and AC-SO3H in four cases: the adsorbate-framework EIs, the adsorbate-framework VIs, the adsorbate-adsorbate EIs, the adsorbate-adsorbate VIs. As shown in Figure 7a, the potential energy of EI between the methanol and AC-None starts at a relatively high level at 0.02 kPa, then drops constantly with an increase of pressure (from -20.1 to -11.9 kJ/mol). When it comes to acetone (see Figure 7b), the vdW force dominates the gas-framework interaction. The VI energy between the acetone and AC-None reaches a high level (-38.6 kJ/mol) at 0.02 kPa, which is remarkably higher than that of methanol (-4.5 kJ/mol). This should be attributed to different configurations of methanol and acetone. As shown in Figure S6, the two methyl groups of acetone interact strongly with the AC-None, while the methanol with typical polar hydroxyl groups shows remarkable EI affinity (hydrogen bonding interaction).65 The varied interaction behaviors of methanol and acetone lead to their different isotherm trends (type-Ⅲ vs. type-I) on AC14

ACS Paragon Plus Environment

Page 15 of 29 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

ACS Applied Materials & Interfaces

None at low pressures. As a consequence, the AC-None with less polar adsorption sites exhibits a low methanol uptake at the low pressures. In terms of the AC-SO3H (see Figure 7c and d). The EI energy between the methanol and AC-SO3H achieves an extremely high level (-69.6 kJ/mol) at 0.02 kPa, while VI energy still remains at a low level (about -3.4 kJ/mol). Meanwhile, the dominating interaction between the acetone and framework changes from VI to EI after the introduction of -SO3H, since the strong polar adsorption sites on the AC-SO3H surface exhibit distinctly EI with the carbonyl group in acetone (see Figure S6). In addition, Figure 7 also shows that the potential energy of gas-gas interactions (both VIs and EIs) increases with the increased pressures. And the interaction between the methanol and methanol is stronger than that between the acetone and acetone. Moreover, the DFT-D calculated results also indicate that the adsorption distance between the methanol and methanol is smaller than that between the acetone and acetone (Figure S7). Obviously, smaller adsorption distance provides more electrostatic interaction contributions for polar adsorptions sites, which finally leads to a larger saturated methanol adsorption capacity on the AC-Rs. 3.4. Separation of Methanol/Acetone Mixture by AC-Rs. The adsorption isotherms of the methanol-acetone mixture (mole bulk composition 50:50) are calculated by the GCMC simulation. The selectivity (S) of methanol over acetone is computed by Eq. 2: 𝑥1

𝑦2

S = 𝑥2 × 𝑦1

(2)

where 𝑥1 and 𝑥2 represent the mole fractions of methanol and acetone in the adsorbed phase, while 𝑦1 and 𝑦2 are the corresponding mole fractions in the bulk phase. Figure 8 presents the methanol-acetone selectivity on ACs at 298 K. The selectivity of methanol over acetone follows the order of AC-NH2 >AC-OH > AC-SO3H >AC-COOH > AC-None under the entire pressure range.

15

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 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 29

Introduction of -NH2, -OH, -SO3H functional groups significantly enhances the selectivity of methanol over acetone at the ultra-low pressure for the strong electron donating/accepting ability of the functional groups (see section 3.3). Surface functionalization provides polar adsorption sites for methanol molecules at the ultra-low pressure resulted from their high electrostatic affinity. Moreover, the tendency of selectivity curves (a rapid descent followed by a steady rise) of the -NH2, -OH, and -SO3H modified AC-Rs is similar to the CO2/CH4 selectivity curves of nanoporous carbons66 and PPN-1.67

Figure 8. Selectivity of methanol over acetone in AC-Rs at 298 K To clarify the separation mechanism of the mixtures, the density distribution contours of the two adsorbates on the AC-NH2 at 0.02, 0.50, 2.00 kPa are presented in Figure 9. The results also give an explanation of the initial decrease in the selectivity curve. At the adsorption onset (0.02 kPa), the methanol molecules with stronger polar hydroxyl functional groups68 quickly take up most of the polar functionality sites by the hydrogen bond interaction, while the acetone molecules only occupy a partial ultramicropore space (see Figure 9a, d). The obvious competitive effects lead to a relatively high methanol-acetone separation rate. When the pressure rises from 0.02 to 0.50 kPa (see Figure 9b, e), the acetone molecules begin to be adsorbed at benzene ring areas due to their 16

ACS Paragon Plus Environment

Page 17 of 29 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

ACS Applied Materials & Interfaces

stronger vdW interaction with the benzene rings. The acetone uptake increases rapidly at this pressure range, which results in a decrease in the methanol-acetone selectivity. With the increase of pressure from 0.50 to 2.00 kPa, methanol tends to stuff the unoccupied pore space (e.g. supermicropores space) (see Figure 9c, f), which results in a constant increase in the methanol-acetone selectivity with a rise of pressure as shown in Figure 8. (a)

(c)

(b)

(b)

Acetone adsorbed at the benzene rings area.

Methanol occupies the

Methanol

polar functionality

unoccupied pore space.

sites.

(d)

stuffs

the

(f)

(e)

Figure 9. Density distribution contours of methanol (a, b, c) and acetone (d, e, f) on AC-NH2 at 0.02 (left), 0.50 (middle), 2.00 kPa (right), respectively. Green represents a strong adsorption site, light blue represents a medium attachment site, and dark blue represents a weak adsorption site. Moreover, to analyze the binary mixture adsorption process at the high pressure range, the Eads of four stable 17

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 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 18 of 29

adsorbate-adsorbate interaction models are established by DFT-D in RASPA (as shown in Figure 10). And the specific steps are presented in SI (Figure S8). The 𝐸𝑎𝑑𝑠 can be described as follows: 𝐸𝑎𝑑𝑠 = 𝐸𝑠𝑖𝑛𝑔𝑙𝑒 𝑎𝑑𝑠𝑜𝑟𝑏𝑎𝑡𝑒 𝑠𝑦𝑠𝑡𝑒𝑚 + 𝐸𝑎𝑑𝑠𝑜𝑟𝑏𝑎𝑡𝑒 ― 𝐸𝑡𝑤𝑜 𝑎𝑑𝑠𝑜𝑟𝑏𝑎𝑡𝑒𝑠 𝑠𝑦𝑠𝑡𝑒𝑚

(3)

where 𝐸𝑠𝑖𝑛𝑔𝑙𝑒 𝑎𝑑𝑠𝑜𝑟𝑏𝑎𝑡𝑒 𝑠𝑦𝑠𝑡𝑒𝑚 is the energy of the BSU with a single adsorbate molecule (methanol- acetone) , 𝐸𝑎𝑑𝑠𝑜𝑟𝑏𝑎𝑡𝑒 is the energy of the methanol or acetone molecule , and

𝐸𝑡𝑤𝑜 𝑎𝑑𝑠𝑜𝑟𝑏𝑎𝑡𝑒𝑠 𝑠𝑦𝑠𝑡𝑒𝑚 is the energy of the

BSU with two adsorbate molecules (the second gas molecule adsorbed on the previous adsorbed one). The 𝐸𝑎𝑑𝑠 for the adsorbed methanol molecule (adsorbed on the polar functionality sites) to catch another methanol and acetone are -47.81 and -21.88 kJ/mol respectively (see Figure 10a, b). And the Eads for the adsorbed acetone molecule (adsorbed on the benzene ring areas) to catch another methanol and acetone are -36.63 kJ/mol and -19.52 kJ/mol respectively (see Figure 10c, d). Obviously, the Eads of the single adsorbate system to capture another methanol molecule is larger than that to capture another acetone molecule, which also suggests the presence of a competitive adsorption effect at the high pressures. Moreover, the acetone molecules adsorbed at benzene ring areas become new polar functionality sites for the methanol molecule, which further enhances the methanol selectivity (see Figure 9c). In addition, we find that AC-SO3H exhibits a higher EI contribution (see Figure 5) but a weaker separation performance than AC-NH2 and AC-OH. The relatively weaker separation performance of AC-SO3H is a result of its smaller accessible surface area (3684 m2/g). Therefore, the effect of accessible surface area on separation is also investigated. Another two AC-SO3H (4133 and 2845 m2/g) models with different accessible surface areas are built and the related physical characteristics is shown in Table S6.

18

ACS Paragon Plus Environment

Page 19 of 29 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

ACS Applied Materials & Interfaces

(b)

(a)

𝐸𝑎𝑑𝑠= -21.88 kJ/mol

𝐸𝑎𝑑𝑠= -47.81 kJ/mol

(d)

(c)

𝐸𝑎𝑑𝑠= -36.63 kJ/mol

𝐸𝑎𝑑𝑠= -19.52 kJ/mol

Figure 10. Stable adsorption configurations of the two adsorbate systems: (a)methanol-methanol; (b) methanolacetone; (c) acetone-acetone; (d) acetone-methanol. The selectivity of methanol from an equimolar methanol–acetone mixture on the three kinds of AC-SO3H 19

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 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 29

at 298 K are shown in Figure 11. The selectivity follows the order of AC-SO3H (4133 m2/g)>AC-SO3H (3684 m2/g)>AC-SO3H (2845 m2/g) in the entire pressure range, which is well consistent with the accessible surface area order. Meanwhile, the adsorption isotherms for the methanol and acetone mixture in the three AC-SO3H models are shown in Figure S9. The isotherms of acetone nearly reach their saturation at 1 kPa, while that of methanol increase continuously with the increasing pressure. Obviously, the selectivity at the high pressure range is dominated by the saturated adsorption capacity of methanol and acetone. The AC-SO3H (4133 m2/g) with a larger Vtot (1.994 cm3/g) provides more space for methanol molecules, leading to a prominent separation performance at high pressures.

Figure 11. Selectivity of methanol over acetone in three AC-SO3H models. Thereby, the obtained results in this section reveal the competitive adsorption effects between the two adsorbates in the entire pressure range. The separation performance of AC-Rs is ruled by the polar functionality site and accessible surface area at low pressures, whereas the functionalization and total available pore volume turns to be more important at high pressures.

20

ACS Paragon Plus Environment

Page 21 of 29 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

4.

ACS Applied Materials & Interfaces

CONCLUSIONS In this study, the adsorption and separation performance of AC-Rs for methanol-acetone has been

investigated. AC-NH2 with the relatively large accessible surface area (4497 m2/g) and Eads (-73.2 kJ/mol for methanol and -63.4 kJ/mol for acetone) exhibits the highest uptake for methanol and acetone at low pressures. AC-SO3H shows the maximum saturated methanol and acetone adsorption capacity at high pressures for its extremely large available pore volume (1.457 cm3/g). Meanwhile, the different configurations between methanol and acetone provide us with a novel approach to separate a binary methanol-acetone mixture. Methanol has a typically polar hydroxyl group and a methyl group, while acetone has a carbonyl group and two methyl groups. The ACs materials with strong electron-donating or electron-accepting functional groups and relative larger accessible surface area (i.e., AC-NH2, AC-OH and AC-SO3H) exhibit a strong methanol adsorption capacity, and thus have a significant selectivity towards methanol. Overall, strong electron-donating or electron-accepting functional groups, large accessible surface area and total available pore volume are the key factors for improving the methanol-acetone separation performance of activated carbons, which provides a superior alternative strategy in the preparation and synthesis of activated carbon for the methanol-acetone mixture separation.

21

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 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 29

ASSOCIATED CONTENT * Supporting Information Lennard-Jones parameters and partial charges for the adsorbents and the adsorbates; Visualization of BSUs of AC-Rs; Stable adsorption configurations of methanol and acetone on BSUs; Physical Characteristics of the three AC-SO3H models; The adsorption isotherms for the equimolar methanol/acetone mixture on the three ACSO3H models.

AUTHOR INFORMATION Corresponding Author * E-mail: [email protected], [email protected] Notes The authors declare no competing financial interest

ACKNOWLEDGMENTS This work was supported by the National Natural Science Foundation of China (No. 21878338), the Key Research and Development Project of Hunan Province, China (No. 2018SK2038).

22

ACS Paragon Plus Environment

Page 23 of 29 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

ACS Applied Materials & Interfaces

REFERENCES (1) Idrissi, A.; Hantal, G.; Jedlovszky, P. Properties of the Liquid–Vapor Interface of Acetone–Methanol Mixtures, as Seen from Computer Simulation and Itim Surface Analysis. Phys. Chem. Chem. Phys. 2015, 17, 8913-8926. (2) Zhang, G.; Wang, H.; Guo, S.; Wang, J.; Liu, J. Synthesis of Cu/TiO2/Organo-Attapulgite Fiber Nanocomposite and Its Photocatalytic Activity for Degradation of Acetone in Air. Appl. Surf. Sci. 2016, 362, 257264. (3) Liu, Z.; Zhang, K.; Wu, Y.; Xi, H. New Functionalized IRMOF-10 with Strong Affinity for Methanol: A Simulation Study. Appl. Surf. Sci. 2018, 440, 351-358. (4) Liu, Z.; Zhang, K.; Wu, Y.; Xi, H. Effective Enhancement on Methanol Adsorption in Cu-BTC by Combination of Lithium-Doping and Nitrogen-Doping Functionalization. J. Mater. Sci. 2018, 53, 6080-6093. (5) Li, M. S.; Wu, S. C.; Peng, Y. H.; Shih, Y. H. Adsorption of Volatile Organic Vapors by Activated Carbon Derived from Rice Husk under Various Humidity Conditions and Its Statistical Evaluation by Linear Solvation Energy Relationships. Sep. Purif. Technol. 2016, 170, 102-108. (6) Sun, X.; Li, Y.; Xi, H.; Xia, Q. Adsorption Performance of a MIL-101 (Cr)/Graphite Oxide Composite for a Series of N-Alkanes. RSC Advances 2014, 4, 56216-56223. (7) Zhang, X.; Gao, B.; Creamer, A. E.; Cao, C.; Li, Y. Adsorption of Vocs onto Engineered Carbon Materials: A Review. J. Hazard. Mater. 2017, 338, 102-123. (8) Tang, L.; Li, L.; Chen, R.; Wang, C.; Ma, W.; Ma, X. Adsorption of Acetone and Isopropanol on Organic Acid Modified Activated Carbons. J. Environ. Chem. Eng. 2016, 4, 2045-2051. (9) Li, L.; Ma, X.; Chen, R.; Wang, C.; Lu, M. Nitrogen-Containing Functional Groups-Facilitated Acetone Adsorption by ZIF-8-Derived Porous Carbon. Materials 2018, 11, 159-164. (10) Ma, X.; Li, L.; Chen, R.; Wang, C.; Zhou, K.; Li, H. Porous Carbon Materials Based on Biomass for Acetone Adsorption: Effect of Surface Chemistry and Porous Structure. Appl. Surf. Sci. 2018, 459, 657-664. (11) Jahandar Lashaki, M.; Atkinson, J. D.; Hashisho, Z.; Phillips, J. H.; Anderson, J. E.; Nichols, M. The Role of Beaded Activated Carbon's Surface Oxygen Groups on Irreversible Adsorption of Organic Vapors. J. Hazard. Mater. 2016, 317, 284-294. 23

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 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 29

(12) Narjinary, M.; Rana, P.; Sen, A.; Pal, M. Enhanced and Selective Acetone Sensing Properties of SnO2MWCNT Nanocomposites: Promising Materials for Diabetes Sensor. Mater. Design 2017, 115, 158-164. (13) Chang, S. J.; Hsueh, T. J.; Chen, I. C.; Hsieh, S. F.; Huang, B. R. Highly Sensitive ZnO Nanowire Acetone Vapor Sensor with Au Adsorption. IEEE T. Nanotechnol. 2008, 7, 754-759. (14) Baur, G. B.; Yuranov, I.; Kiwi-Minsker, L. Activated Carbon Fibers Modified by Metal Oxide as Effective Structured Adsorbents for Acetaldehyde. Catal. Today 2015, 249, 252-258. (15) Wang, D.; Wu, G.; Zhao, Y.; Cui, L.; Shin, C. H.; Ryu, M. H.; Cai, J. Study on the Copper(II)-Doped MIL101(Cr) and Its Performance in Vocs Adsorption. Environ. Sci. Pollut. R. 2018, 25, 28109–28119. (16) Shafiei, M.; Alivand, M. S.; Rashidi, A.; Samimi, A.; Mohebbi-Kalhori, D. Synthesis and Adsorption Performance of a Modified Micro-Mesoporous MIL-101(Cr) for Vocs Removal at Ambient Conditions. Chem. Eng. J. 2018, 341, 164-174. (17) Jangodaz, E.; Alaie, E.; Safekordi, A. A.; Tasharrofi, S. Adsorption of Ethylbenzene from Air on Metal– Organic Frameworks MIL-101(Cr) and MIL-53(Fe) at Room Temperature. J. Inorg Organomet. P. 2018, 28, 2090–2099. (18) Lashaki, M. J.; Atkinson, J. D.; Hashisho, Z.; Phillips, J. H.; Anderson, J. E.; Nichols, M. The Role of Beaded Activated Carbon’s Pore Size Distribution on Heel Formation During Cyclic Adsorption/Desorption of Organic Vapors. J. Hazard. Mater. 2016, 315, 42-51. (19) Sun, X.; Li, Y.; Xi, H.; Xia, Q. Adsorption Performance of a MIL-101(Cr)/Graphite Oxide Composite for a Series of N-Alkanes. Rsc Advances 2014, 4, 56216-56223. (20) Vega, E.; Lemus, J.; Anfruns, A.; Gonzalez-Olmos, R.; Palomar, J.; Martin, M. J. Adsorption of Volatile Sulphur Compounds onto Modified Activated Carbons: Effect of Oxygen Functional Groups. J. Hazard. Mater. 2013, 258-259, 77-83. (21) Gao, J.; Li, L.; Zeng, Z.; Ma, X.; Chen, R.; Wang, C.; Zhou, K. Superior Acetone Uptake of Hierarchically N-Doped Potassium Citrate-Based Porous Carbon Prepared by One-Step Carbonization. J. Mater. Sci. 2019, 54, 6186-6198. (22) Zhu, J.; Chen, R.; Zeng, Z.; Su, C.; Zhou, K.; Mo, Y.; Guo, Y.; Zhou, F.; Gao, J.; Li, L. Acetone Adsorption Capacity of Sulfur-Doped Microporous Activated Carbons Prepared from Polythiophene. Environ. Sci. Pollut. R. 24

ACS Paragon Plus Environment

Page 25 of 29 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

ACS Applied Materials & Interfaces

2019, 26, 16166–16180. (23) Shevade, A. V.; Jiang, S.; Gubbins, K. E. Molecular Simulation Study of Water–Methanol Mixtures in Activated Carbon Pores. J. Chem. Phys. 2000, 113, 6933-6942. (24) Ma, X.; Li, L.; Zeng, Z.; Chen, R.; Wang, C.; Zhou, K.; Su, C.; Li, H. Synthesis of Nitrogen-Rich Nanoporous Carbon Materials with C3N-Type from ZIF-8 for Methanol Adsorption. Chem. Eng. J. 2019, 363, 49-56. (25) Li, L.; Sun, Z.; Li, H.; Keener, T. C. Effects of Activated Carbon Surface Properties on the Adsorption of Volatile Organic Compounds. J. Air Waste Manage 2012, 62, 1196-1202. (26) Qiu, W.; Kang, D.; Ying, Z.; Huang, H.; Chen, Y.; Lu, H. Hierarchical Pore Structure of Activated Carbon Fabricated by CO2/Microwave for Vocs Adsorption. Chin. J. Chem. Eng. 2017, 26, 81-88. (27) Greathouse, J. A.; Hart, D. B.; Ochs, M. E. Alcohol and Thiol Adsorption on (Oxy)Hydroxide and Carbon Surfaces: Molecular Dynamics Simulation and Desorption Experiments. J. Phys. Chem. C 2015, 116, 26756– 26764. (28) Myers, A. L.; Prausnitz, J. M. Thermodynamics of Mixed-Gas Adsorption. AlChE J. 2010, 11, 121-127. (29) Ushiki, I.; Ota, M.; Sato, Y.; Inomata, H. Prediction of Vocs Adsorption Equilibria on Activated Carbon in Supercritical Carbon Dioxide over a Wide Range of Temperature and Pressure by Using Pure Component Adsorption Data: Combined Approach of the Dubinin–Astakhov Equation and the Non-Ideal Adsorbed Solution Theory (NIAST). Fluid Phase Equilib. 2014, 375, 293-305. (30) Pham, T.; Forrest, K. A.; Banerjee, R.; Orcajo, G.; Eckert, J.; Space, B. Understanding the H2 Sorption Trends in the M-MOF-74 Series (M = Mg, Ni, Co, Zn). J. Phys. Chem. C 2014, 119, 1078-1090. (31) Wu, Y.; Chen, H.; Xiao, J.; Liu, D.; Liu, Z.; Qian, Y.; Xi, H. Adsorptive Separation of Methanol-Acetone on Isostructural Series of Metal-Organic Frameworks M-BTC (M = Ti, Fe, Cu, Co, Ru, Mo): A Computational Study of Adsorption Mechanisms and Metal-Substitution Impacts. ACS Appl. Mater. Interfaces 2015, 7, 2693026940. (32) Pillai, R. S.; Pinto, M. L.; Pires, J.; Jorge, M.; Gomes, J. R. B. Understanding Gas Adsorption Selectivity in IRMOF-8 Using Molecular Simulation. ACS Appl. Mater. Interfaces 2015, 7, 624-637. (33) Assfour, B.; Seifert, G. Adsorption of Hydrogen in Covalent Organic Frameworks: Comparison of 25

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 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 26 of 29

Simulations and Experiments. Microporous Mesoporous Mater. 2010, 133, 59-65. (34) Choi, Y. J.; Choi, J. H.; Choi, K. M.; Kang, J. K. Covalent Organic Frameworks for Extremely High Reversible CO2 Uptake Capacity: A Theoretical Approach. J. Mater. Chem. 2011, 21, 1073-1078. (35) Li, X.; Xue, Q.; He, D.; Lei, Z.; Teng, Z. Sulfur-Nitrogen Codoped Graphite Slit-Pore for Enhancing Selective Carbon Dioxide Adsorption: Insights from Molecular Simulations. ACS Sustainable Chem. Eng. 2017, 5,8815-8823. (36) Wang, S.; Lu, L.; Wu, D.; Lu, X.; Cao, W.; Yang, T.; Zhu, Y. Molecular Simulation Study of the Adsorption and Diffusion of a Mixture of CO2/CH4 in Activated Carbon: Effect of Textural Properties and Surface Chemistry. J. Chem. Eng. Date 2016, 61, 4139-4147. (37) Chu, X.-z.; Liu, S.-s.; Zhou, S.-y.; Zhao, Y.-j.; Xing, W.-h.; Lee, C.-H. Adsorption Behaviors of CO2 and CH4 on Zeolites JSR and NanJSR Using the GCMC Simulations. Adsorption 2016, 22, 1065-1073. (38) Hui, F.; Wang, Y.; Zhang, T.; Yang, C.; Shan, H. Adsorption and Separation Mechanism of Thiophene/Benzene in MFI Zeolite: A GCMC Study. J. Phys. Chem. C 2017, 121, 25818-25826. (39) Liang, X.; Chi, J.; Yang, Z. The Influence of the Functional Group on Activated Carbon for Acetone Adsorption Property by Molecular Simulation Study. Microporous Mesoporous Mater. 2018, 262, 77-88. (40) Li, X.; Xue, Q.; He, D.; Zhu, L.; Du, Y.; Xing, W.; Zhang, T. Sulfur–Nitrogen Codoped Graphite Slit-Pore for Enhancing Selective Carbon Dioxide Adsorption: Insights from Molecular Simulations. ACS Sustainable Chem. Eng. 2017, 5, 8815-8823. (41) Li, X.; Xue, Q.; Chang, X.; Zhu, L.; Ling, C.; Zheng, H. Effects of Sulfur Doping and Humidity on CO2 Capture by Graphite Split Pore: A Theoretical Study. ACS Appl. Mater. Interfaces 2017, 9, 8336-8343. (42) Kumar, K. V.; Müller, E. A.; Rodríguez-Reinoso, F. Effect of Pore Morphology on the Adsorption of Methane/Hydrogen Mixtures on Carbon Micropores. J. Phys. Chem. C 2012, 116, 11820-11829. (43) Biase, E. D.; Sarkisov, L. Systematic Development of Predictive Molecular Models of High Surface Area Activated Carbons for Adsorption Applications. Carbon 2013, 64, 262-280. (44) Lu, X.; Jin, D.; Wei, S.; Zhang, M.; Zhu, Q.; Shi, X.; Deng, Z.; Guo, W.; Shen, W. Competitive Adsorption of a Binary CO2-CH4 Mixture in Nanoporous Carbons: Effects of Edge-Functionalization. Nanoscale 2015, 7, 1002-1012. 26

ACS Paragon Plus Environment

Page 27 of 29 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

ACS Applied Materials & Interfaces

(45) García-Pérez, E.; Parra, J.; Ania, C.; García-Sánchez, A.; Van Baten, J.; Krishna, R.; Dubbeldam, D.; Calero, S. A. Computational Study of CO2, N2, and CH4 Adsorption in Zeolites. Adsorption 2007, 13, 469-476. (46) Alder, B. J.; Wainwright, T. E. Studies in Molecular Dynamics. Ⅱ. Behavior of a Small Number of Elastic Spheres. J. Chem. Phys. 1960, 33, 1349-1451. (47) Frisch, M.; Trucks, G.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. Gaussian 09, Revision A.02, Gaussian. Inc., Wallingford, CT 2009, 200. (48) Xu, S. C.; Irle, S.; Musaev, D. G.; Lin, M. C. Quantum Chemical Prediction of Pathways and Rate Constants for Reactions of CO and CO2 with Vacancy Defects on Graphite (0001) Surfaces. J. Phys. Chem. C 2009, 113, 18772-18777. (49) Radovic, L. R. Active Sites in Graphene and the Mechanism of CO2 Formation in Carbon Oxidation. J. Am. Chem. Soc. 2009, 131, 17166-17175. (50) Mo, Y.; Li, H.; Zhou, K.; Ma, X.; Guo, Y.; Wang, S.; Li, L. Acetone Adsorption to (BeO) 12,(MgO) 12 and (ZnO) 12 Nanoparticles and Their Graphene Composites: A Density Functional Theory (DFT) Study. Appl. Surf. Sci. 2019, 469, 962-973. (51) Zhou, K.; Li, L.; Ma, X.; Mo, Y.; Chen, R.; Li, H.; Li, H. Activated Carbons Modified by Magnesium Oxide as Highly Efficient Sorbents for Acetone. RSC Advances 2018, 8, 2922-2932. (52) Ao, Z.; Li, S.; Jiang, Q. Correlation of the Applied Electrical Field and CO Adsorption/Desorption Behavior on Al-Doped Graphene. Solid State Commun. 2010, 150, 680-683. (53) Lim, D.-H.; Negreira, A. S.; Wilcox, J. Dft Studies on the Interaction of Defective Graphene-Supported Fe and Al Nanoparticles. J. Phys. Chem. C 2011, 115, 8961-8970. (54) Dubbeldam, D.; Calero, S.; Ellis, D. E.; Snurr, R. Q. Raspa: Molecular Simulation Software for Adsorption and Diffusion in Flexible Nanoporous Materials. Mol. Simulat. 2015, 42, 81-101. (55) Dubbeldam, D.; Torres-Knoop, A.; Walton, K. S. On the Inner Workings of Monte Carlo Codes. Mol. Simulat. 2013, 39, 1253-1292. (56) Martin, M. G.; Siepmann, J. I. Novel Configurational-Bias Monte Carlo Method for Branched Molecules. Transferable Potentials for Phase Equilibria. 2. United-Atom Description of Branched Alkanes. J. Phys. Chem. B 1999, 103, 4508-4517. 27

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 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 28 of 29

(57) Mayo, S. L.; Olafson, B. D.; Goddard, W. A. Dreiding: A Generic Force Field for Molecular Simulations. J. Phys. Chem. 1990, 94, 8897-8909. (58) Maerzke, K. A.; Schultz, N. E.; Ross, R. B.; J Ilja, S. Trappe-Ua Force Field for Acrylates and Monte Carlo Simulations for Their Mixtures with Alkanes and Alcohols. J. Phys. Chem. B 2009, 113, 6415-6425. (59) Chen, B.; Potoff, J. J.; Siepmann, J. I. Monte Carlo Calculations for Alcohols and Their Mixtures with Alkanes. Transferable Potentials for Phase Equilibria. 5. United-Atom Description of Primary, Secondary, and Tertiary Alcohols. J. Phys. Chem. B 2015, 105, 3093-3104. (60) Bhatia, S. K.; Myers, A. L. Optimum Conditions for Adsorptive Storage. Langmuir 2006, 22, 1688-1700. (61) Peigney, A.; Laurent, C.; Flahaut, E.; Bacsa, R.; Rousset, A. Specific Surface Area of Carbon Nanotubes and Bundles of Carbon Nanotubes. Carbon 2001, 39, 507-514. (62) Hu, Z.; Srinivasan, M. P.; Ni, Y. Preparation of Mesoporous High‐Surface‐Area Activated Carbon. Adv. Mater. 2000, 12, 62-65. (63) Brunauer, S.; Deming, L. S.; Deming, W. E.; Teller, E. On a Theory of the Vander Waals Adsorption of Gases. J. Am. Chem. Soc. 1940, 62, 1723-1732. (64) Nguyen, V. T.; Do, D. D.; Nicholson, D.; Jagiello, J. Effects of Temperature on Adsorption of Methanol on Graphitized Thermal Carbon Black: A Computer Simulation and Experimental Study. J. Phys. Chem. C 2011, 115, 16142-16149. (65) Shevade, A. V.; Jiang, S.; Gubbins, K. E. Molecular Simulation Study of Water–Methanol Mixtures in Activated Carbon Pores. J. Chem. Phys. 2000, 113, 6933-6942. (66) Dang, Y.; Guo, W.; Zhao, L.; Zhu, H. Porous Carbon Materials Based on Graphdiyne Basis Units by the Incorporation of the Functional Groups and Li Atoms for Superior CO2 Capture and Sequestration. ACS Appl. Mater. Interfaces 2017, 9, 30002-30013. (67) Lu, W.; Yuan, D.; Zhao, D.; Schilling, C. I.; Plietzsch, O.; Muller, T.; Bräse, S.; Guenther, J.; Blümel, J.; Krishna, R. Porous Polymer Networks: Synthesis, Porosity, and Applications in Gas Storage/Separation. Chem. Mater. 2010, 22, 5964-5972. (68) Szymanowski, D. H. I. J.; Voelkel, I. A.; Binte, H. J. Polarity of Some Oxime Extractants as Determined by Gas‐Liquid‐Chromatography. Adv. Synth. Catal. 2010, 327, 353-361. 28

ACS Paragon Plus Environment

Page 29 of 29

Table of Content graphic

Methanol

Methanol-Acetone

Selectivity (Methanol/Acetone)

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

ACS Applied Materials & Interfaces

AC-Rs

Pressure (kPa)

(R= COOH, OH, NH2, SO3H) BSUs

29

ACS Paragon Plus Environment