Supporting Information Retention and Release of Graphene Oxide in Structured Heterogeneous Porous Media under Saturated and Unsaturated Conditions Shunan Dong 1, Xiaoqing Shi 1, Bin Gao 3, Jianfeng Wu 1, Yuanyuan Sun 1 *, Hongyan Guo 2, Hongxia Xu 1, Jichun Wu 1 * 1 State Key Laboratory of Pollution Control and Resource Reuse, Key Laboratory of Surficial Geochemisty, Ministry of Education, School of Earth Sciences and Engineering, Hydrosciences Department, Nanjing University, Nanjing 210023, China 2 State Key Laboratory of Pollution Control and Resource Reuse, School of Environment, Nanjing University, Nanjing 210023, China, 3 Department of Agricultural and Biological Engineering, University of Florida, Gainesville, FL 32611 Contents of this file: 11 pages, 4 figures, and 4 tables. * Corresponding authors. Tel.: +86 25 89680835; fax: +86 25 8368 6016. E-mail address: sunyy@nju.edu.cn (Y. Sun), jcwu@nju.edu.cn (J. Wu). S1
S1 XDLVO Theory The extended Derjaguin-Landau-Verwey-Overbeek (XDLVO) interactions between GO particles and porous media, GO and GO, and GO and air-water interface can be modeled as the integration of the Lifshitz-van der Waals (LW) attraction, the electrical double layer (EDL) repulsion and the Lewis acid-base (AB) interactions between two parallel plates. 1-3 The Lifshitz-van der Waals attraction per unit area (Φ LW ), the electrostatic double-layer repulsion per unit area (Φ EDL ), the Lewis acid-base (AB) interactions per unit area (Φ AB ), and the total interaction energy per unit area (Φ tot ) can be calculated as follows: 1-3 A Φ LW = - 12πh Φ Φ EDL 132 2 kt 2 = 32ε0εrκ12 ( ) exp( - κh) (2) ve h0 h AB = ΔG exp( ) (3) - AB h0 λ (1) tot EDL LW AB Φ =Φ +Φ +Φ (4) where A 132 represents the Hamaker constant for substances l and 2 in presence of medium 3, which can be determined from the Hamaker constant of each material, h is the separation distance between the GO particles and sand surface, 0 is the dielectric permittivity of vacuum (8.85*10-12 CV -1 m -1 ), r is the dielectric constant of water (78.5), is the Debye length, k is the Boltzmann s constant, T is the absolute temperature, v is the valence of electrolyte, e is the electron charge (1.60 10-19 C), (=0.6 nm).the Hamaker constants of S2
between GO and quartz sand and between GO and air-water interface are calculated as 9.8 10-21 and -1.81 10-20 J, respectively. 1, 2 The values of and were calculated using the following equation: 1 0 r kt 2N Ie A 2 (5) 1 2 veξ1 = tanh( ) 4kT veξ 2 = tanh( ) 4kT (6) (7) where 0 and r were previous defined, N A is the Avogadro constant, I is the ionic strength, 1 and 2 were the zeta potential of GO particle and sand respectively. To obtain the values of AB Gh 0, the contact angles (θ) of three probing liquids (water, glycerol and diiodomethane) need to be acquired first. The method was similar to the study of Xu and Ong. 2, 4 According to the followed equations, the LW ( LW GO ), electron-accepting ( + GO ) and electron-donating ( ) interfacial tension values for GO could be determined through the ḠO following equations: L LW LW + - - + i ( 1+ cos θ) = 2 i + 2 i + 2 i (8) where i represents water ( L LW + w = 72.8, w = 21. 8 and w = LW + 64.0, = 34. 0 and = g and + d = - d = 0 mj m -2 ), respectively. 5 The values of g - w 3.92 and = 57.4 mj m -2 ) or diiodomethane ( L ḡ d = 25.5 mj m -2 ), glycerol ( L g = AB Gh 0 were then calculated from the interfacial tension parameters: LW = 50.8, = 50. 8 d - - - - - =2[ ( - ) ( - ) - - ] (9) G 0 AB h w GO j w w GO j w GO j GO j S3
where the subscript j represents sand ( + s = 1.4 mj m -2, - s = 47.8 mj m -2 ), and air ( + a = = 0 mj ā m -2 ). 6, 7 S2 Determination of Pore Water Velocity, Dispersion Coefficient, and Moisture Content S2.1 Saturated columns For all the saturated column experiments, coarse sand area was the fast flow domain () and fine sand area was the slow flow domain (), the pore water velocities (v and v ) corresponding to the two areas were directly calculated as: v = [ ωk uk + (1 - ω) K ][ θ ω + θ (1 - ω)] (10) v = [ ωk uk + (1 - ω) K ][ θ ω + θ (1 - ω)] (11) where the subscripts and refer to the fast flow and slow flow domains in the column, respectively; u was the Darcy velocity total applied in the saturated columns (L T -1 ); ω is the ratio of the volumes of the coarse sand and the total sand system (-); K is the saturated hydraulic conductivity (L T -1 ), which were measured using the constant head method based on the Darcy s Law; θ is the water content (-), which is the porosity of the sand under saturated conditions. The dispersion coefficient (D) of the two areas was determined as: D = δ v (12) S4
D = δ v (13) where δ is the saturated longitudinal dispersivity (L), which were determined from the tracer experiments of the corresponding saturated homogeneous columns. S2.1 Unsaturated columns For unsaturated columns, fine sand area was the fast flow domain () and coarse sand area was the slow flow domain (). Flow calculations of and were carried out using the TOUGH2 EOS9 module based on the dual-permeability conceptual model (implemented using Multiple Interacting Continua method in TOUGH2), 8 which simulated the whole draining process (during the preparation of the unsaturated columns) and the unsaturated flow according to Richards equation. 9 Capillary pressure and relative permeability functions of saturation for the dual domains are based on the van Genuchten formulation. 10 The material properties that were used for the dual domains were referred to the van Genuchten formulation 10 and modified based on the corresponding saturated permeability values obtained from experimental measurements. The equivalent pore water velocities (v) and moisture contents (θ) of and were calculated from the numerical simulator TOUGH2 and used in the GO retention and transport model (equation 1 and 2 in the manuscript). S5
The unsaturated longitudinal dispersivity for was also obtained from the tracer experiments of unsaturated homogeneous columns (with corresponding sand and moisture content), which was then used to calculate the dispersion coefficient (D). Table S1. Surface Properties of the GO and Sands and XDLVO Results Sand (0.70-0.85 mm) Sand (0.45-0.50 mm) Sand (0.35-0.45 mm) Sand (0.25-0.35 mm) Sand (0.15-0.20 mm) Sand (0.11-0.15 mm) GO IS (mm) Zeta-potential (mv) (kt μm -2 ) 1-50.5 -- 20-37.6-161 1-44.5 -- 20-34.4-165 1-41.4 -- 20-32.3-168 1-43.2 -- 20-36.8-162 1-42.0 -- 20-35.9-163 1-43.6 -- 20-33.0-167 1-42.1 -- 20-28.7-152 Table S2. Summary of the β, Peclet and Damkohler Number for All the Experiments Peclet number β Damkohler number Damkohler Damkohler number Column (cm -1 ) for mass transfer number for 1 mm for 20 mm S1 1.54 10 2 2.09 10 2 2.71 3.15 0.19 0.57 0.0277 0.103 0.80 5.57 S2 1.54 10 2 3.91 10 2 2.62 3.17 0.09 0.57 0.0241 0.288 0.72 11.4 S3 1.54 10 2 5.17 10 2 2.66 3.04 0.07 0.84 0.0219 0.581 0.63 25.4 S4 3.30 10 2 3.91 10 2 2.66 3.10 0.32 0.62 0.0601 0.115 1.91 4.45 S5 3.30 10 2 5.17 10 2 2.62 3.11 0.13 0.51 0.0507 0.205 1.62 9.31 S6 3.30 10 2 8.06 10 2 2.48 3.16 0.05 0.58 0.0465 0.625 1.45 35.0 U1 1.57 10 2 -- 3.65 7.49 0.55 54.2 0.0247 9.84 0.83 52.7 U2 2.98 10 2 -- 3.15 13.3 0.30 83.0 0.0296 15.6 1.00 84.0 U3 4.15 10 2 -- 3.31 9.83 0.41 54.7 0.0285 10.6 1.42 74.7 Where and represent fast flow domain and slow flow domain respectively. β represent geometry coefficient S ( β =, where S is the contact area between coarse and fine grains packed in the column, V is the volume of each Vθ S6
vh domain, and θ is the moisture content). Peclet number ( Pe =, where v is the velocity of pore water, H is the height of D the column, D is the dispersion coefficient) compares the advection and dispersion processes. Damkohler number βαh (Damkohler number for mass transfer: Da MT =, where α is the mass transfer coefficient; Damkohler numbers for v kh GO deposition: Da CD =, where k is the retention rate constant) compares the time scale of the two processes; when v its value is >>1 or <<1, one of the process dominants. Table S3. Summary of Recovery Rate for All Columns According to GO Transport Experiment and Model IS 1 mm IS 20 mm Column MR ** MR ** MR Total ** MR Total MR ** MR ** MR Total ** MR Total S1 71.2 23.7 94.9 93.3 43.0 3.3 46.3 46.0 S2 84.0 10.1 94.1 93.5 53.5 0.4 53.9 51.5 S3 89.0 4.9 93.9 95.7 58.7 0.2 58.9 57.2 S4 61.7 30.0 91.7 90.4 25.4 5.1 30.5 31.0 S5 75.1 16.1 91.2 90.7 34.3 0.6 34.9 35.1 S6 87.4 5.1 92.5 90.0 44.7 0.1 44.8 46.1 U1 87.7 1.8 89.5 89.4 55.2 1.0 56.2 56.2 U2 89.1 1.4 90.5 89.7 55.8 0.7 56.5 55.1 U3 88.7 1.8 90.5 90.6 46.3 0.7 47.0 46.6 Where MR represent mass recovery rate from effluent according to GO transport process; and represent fast flow domain and slow flow domain; superscript ** represent parameters determined from the GO transport model. Table S4. Summary of Recovery Rate for All Columns According to GO Release Process Where MRR represent mass recovery rate from effluent according to GO release process; and represent fast flow domain and slow flow domain. Column MRR MRR MRR Total S1 15.8 14.0 29.8 S2 14.9 6.2 21.1 S3 14.1 5.1 19.2 S4 17.3 14.7 32.0 S5 17.8 10.9 28.7 S6 14.2 9.3 23.5 U1 -- -- 12.1 S7
U2 -- -- 13.1 U3 -- -- 22.8 S8
(a) (b) (c) (d) (e) (f) Figure S1. XDLVO energy between GO and sand grain of different size (a: 0.70-0.85 mm; b: 0.45-0.50 mm; c: 0.35-0.45 mm; d: 0.25-0.35 mm; e: 0.15-0.20 mm; f: 0.11-0.15 mm) under different NaCl solution (1 mm and 20 mm NaCl). Figure S2. XDLVO energy between GO and GO under different NaCl solution (1 mm and 20 mm NaCl). S9
Figure S3. XDLVO energy between GO and air-water interface under different NaCl solution (1 mm and 20 mm NaCl). (a) (b) (c) (d) Figure S4. SEM micrographs of cleaned quartz sand (a) and the quartz sand excavated from the inlet of the packed column (b), and the corresponding EDX spectrum of cleaned quartz sand (c) and the quartz sand excavated from the inlet of the packed column (d). S10
REFERENCES (1) Liu, L.; Gao, B.; Wu, L.; Morales, V. L.; Yang, L.; Zhou, Z.; Wang, H. Deposition and transport of graphene oxide in saturated and unsaturated porous media. Chem. Eng. J. 2013, 229, 444-449. (2) Feriancikova, L.; Xu, S. Deposition and remobilization of graphene oxide within saturated sand packs. J. Hazard. Mater. 2012, 235, 194-200. (3) Wu, L.; Liu, L.; Gao, B.; Muñoz-Carpena, R.; Zhang, M.; Chen, H.; Zhou, Z.; Wang, H. Aggregation kinetics of graphene oxides in aqueous solutions: Experiments, mechanisms, and modeling. Langmuir 2013, 29 (49), 15174-15181. (4) Ong, Y.; Razatos, A.; Georgiou, G.; Sharma, M. Adhesion forces between E. coli bacteria and biomaterial surfaces. Langmuir 1999, 15 (8), 2719-2725. (5) Vanoss, C. J. Acid-base interfacial interactions in aqueous media. Colloids Surf., A 1993, 78, 1-49. (6) Vanoss, C. J. Interfacial forces in aqueous media; CRC Press: Boca Raton, U.S.A., 2006. (7) Morrow, J. B.; Stratton, R.; Yang, H. H.; Smets, B. F.; Grasso, D. Macro- and manoscale observations of adhesive behavior for several E-coli strains (O157: H7 and environmental isolates) on mineral surfaces. Environ. Sci. Technol. 2005, 39 (17), 6395-6404. (8) Pruess, K.; Oldenburg, C.; Moridis, G. TOUGH2 user's guide version 2; Lawrence Berkeley National Laboratory: Berkeley, U.S.A., 1999. S11
(9) Richards, L. A. Capillary conduction of liquids through porous mediums. J. Appl. Phys. 1931, 1 (5), 318-333. (10) van Genuchten, M. T. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci. Soc. Am. J. 1980, 44 (5), 892-898. S12