Nanobubble Retention in Saturated Porous Media under Repulsive

May 6, 2019 - One-dimensional column experiments were therefore conducted to investigate the transport, retention, and release of NBs in glass beads ...
0 downloads 0 Views 662KB Size
Subscriber access provided by The Library | University of Bath

Interfaces: Adsorption, Reactions, Films, Forces, Measurement Techniques, Charge Transfer, Electrochemistry, Electrocatalysis, Energy Production and Storage

Nano-Bubble Retention in Saturated Porous Media under Repulsive van der Waals and Electrostatic Conditions Shoichiro Hamamoto, Takuya Sugimoto, Takato Takemura, Taku Nishimura, and Scott A. Bradford Langmuir, Just Accepted Manuscript • DOI: 10.1021/acs.langmuir.9b00507 • Publication Date (Web): 06 May 2019 Downloaded from http://pubs.acs.org on May 14, 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 37 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

Langmuir

Nano-Bubble Retention in Saturated Porous Media under Repulsive van der Waals and Electrostatic Conditions Shoichiro Hamamoto,1,* Takuya Sugimoto,1

Takato Takemura,2 Taku Nishimura,1

Scott A. Bradford3

1Graduate

School of Agricultural and Life Sciences, The University of Tokyo, 1-1-1,

Yayoi, Bunkyoku, Tokyo, 113-8657, Japan.

2College

of Humanities and Sciences, Nihon University, 3-25-40, Sakurajousui,

Setagaya, Tokyo, 156-8550, Japan

3US

Salinity Laboratory, USDA, ARS, 450 W. Big Springs Road, Riverside, CA

92507-4617, USA

KEYWORDS : Nano-Bubble; Retention; Hydrophobic Interaction; Ion Species; Surface Roughness

1 ACS Paragon Plus Environment

Langmuir 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

ABSTRACT - An understanding of nano-bubble (NB) migration in porous media is needed for potential environmental applications. The solution chemistry is well-known to be a critical factor in determining interactions of other colloids and nanoparticles with surfaces. However, little quantitative research has examined the influence of solution chemistry on NB transport. One-dimensional column experiments were therefore conducted to investigate the transport, retention, and release of NBs in glass beads under different solution chemistry conditions. NB concentrations in the effluent were reduced with an increase in ionic strength (IS) or a decrease in pH due to a reduction in the repulsive force between the glass surface and NBs, especially when the solution contained Ca2+ as compared to Na+ and for larger NBs. This result was somewhat surprising because electrostatic and van der Waals interactions for NBs were both repulsive on a homogeneous glass bead surface. NB retention on the surface was explained by ubiquitous nanoscale roughness on the glass beads that significantly lowered the energy barrier, and localized attractive charge heterogeneity and/or hydrophobic interactions. In contrast to Na+, adsorbed Ca2+ ions produced charge heterogeneity that enhanced NB retention and inhibited release with IS reduction.

2 ACS Paragon Plus Environment

Page 2 of 37

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

Langmuir

ABSTRACT ART

3 ACS Paragon Plus Environment

Langmuir 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 37

INTRODUCTION Nano-scale bubbles (NBs) have recently gained increasing attention due to their unique physicochemical properties (large surface area, charged surface, stability, and high gas internal pressure) and many potential applications such as a use in water treatment by flotation  1,2,

and in sterilization applying ozone gas 3, and for accelerating the metabolism in shellfishes

and vegetables 4. In addition, NBs show considerable potential for soil remediation of organic chemicals in combination with bioremediation

5–8.

An understanding of NB transport

characteristics in soils is therefore required for successful soil remediation and agricultural applications using NBs. However, fundamental studies on NB behavior in porous media are very limited 9–11. The transport of solid colloids (e.g., microorganisms, clays, and nanoparticles) in porous media is well-known to depend on the solution pH, ionic strength (IS), and ion composition 12–14. The interaction of solid colloids with porous medium is usually explained by considering the sum of van der Waals and electric double layer interactions

15.

The van der

Waals interaction between solid surfaces is attractive 16, whereas it is repulsive for NBs and a porous media surface

17.

Most porous medium surfaces, including glass beads, exhibit a net

negative charge at neutral or alkaline pH conditions 18. NBs also exhibit a negative charge under these pH conditions due to the preferential adsorption of OH– ions

4 ACS Paragon Plus Environment

19,20,21.

Consequently, a

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

Langmuir

repulsive electrostatic interaction is expected between glass beads and NBs, although compression of the double layer thickness and alteration of the zeta potentials with increasing IS are expected to decrease this repulsion

22.

This indicates that consideration of only van der

Waals and double layer interactions cannot explain the observed retention of NBs in porous media. Furthermore, no systematic studies have examined the influence of solution IS and ionic composition on NB retention to date. Alternative explanations for observed NB retention in porous media therefore requires consideration of other interactions, factors, and/or processes. The air-water interface has been reported to be extremely hydrophobic 23. However, previous research examining the transport of NBs has not yet considered the potential influence of hydrophobic interactions 9. Ubiquitous nanoscale roughness and chemical heterogeneity on porous medium surface have also been demonstrated to locally reduce and/or eliminate the energy barrier to attachment of solid colloids on porous media surfaces

24.

Multivalent cations are reported to produce greater

amounts of solid colloid retention at the same IS than monovalent cations because of charge neutralization or reversal, and/or cation bridging

25–27.

However, no research to date has

considered the potential influence of nanoscale heterogeneities on NB retention. Straining of NB has previously been reported as one of major contributors to NB retention 10. However, they have not considered the potentially confounding effects of hydrophobic interactions and

5 ACS Paragon Plus Environment

Langmuir 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 37

nanoscale heterogeneities on retention. Consequently, there are still many unanswered questions about factors that influence NB retention. NB release has also previously been observed with a decreasing in solution IS

10,11.

The electric double layer interaction is typically more repulsive under lower IS conditions because of expansion of the double layer thickness and larger magnitudes in the zeta potentials  14,28.

Consequently, NB release with IS reduction cannot be explained by traditional interaction

energy calculations due to the increased height of energy barrier on smooth surfaces

29.

In

addition, colloid release with IS reduction is reported to be much greater when solid colloids are deposited in the presence of monovalent cations than divalent cations such as Ca2+

25–27.

Nanoscale roughness has been reported to theoretically explain the release of solid and hollow colloids with IS reduction 30, but experimental information for NBs is still lacking, especially in the presence of divalent cations. This research investigates factors that control the transport, retention, and release of air NBs in glass beads. One-dimensional column transport studies with air NBs were conducted under different solution IS, ionic composition, and pH conditions. Collected data were simulated using a mathematical model that considered advective and dispersive transport, first-order kinetic retention, and transient release. Results were interpreted with the aid of interaction energy calculations that considered electric double layer, van der Waals, and Lewis

6 ACS Paragon Plus Environment

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

Langmuir

acid-base interactions on surfaces with different amounts of nanoscale roughness. Findings from this study help to explain the critical roles of hydrophobic interactions, chemical heterogeneity, and nanoscale roughness in controlling the retention and release of NBs in porous media under different solution chemistry conditions.

MATERIALS AND METHODS Materials NBs were produced by a pressurized dissolution method (FZIN-10-I, IDEC Co., Ltd., Japan). In this method, air was allowed to dissolve in water by applying a pressure of approximately 300 kPa and water was flushed through a nozzle inside the water pool. The supersaturated condition was released due to the reduction in the pressure into the expelled water in the form of bubbles. DI (deionized) water, where distilled water was ion-exchanged and filtered (Autostill WA53, Yamamoto Scientific Co., Ltd., Japan), was used for producing the NBs. The NB generator was operated for around 60 min in 10 L of water, resulting in NB water with peak bubble size of around 170 nm and bubble concentration ranging from 1.0-2.0 x 105 numbers ml-1. The bubble-size (diameter) distributions and the concentrations of NB water were measured using a resonant mass measurement (Archimedes, Malvern Instruments Ltd.) 31. We studied the effect of freezing and thawing on NB suspension by keeping the samples at a

7 ACS Paragon Plus Environment

Langmuir 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

temperature of – 18 oC for 24 hr and then leaving to thaw at room temperature. As shown in Fig. S1, the majority of NBs disappeared after thawing without any significant agglomeration, suggesting the rupture of the gas filled bubbles 32,33. However, the small amount of nano-entities which have not disappeared after thawing likely suggests the existence of impurities whose densities are lower than water; e.g., oil contaminants. Some recent theoretical and experimental works showed that organic (hydrophobic) substances partially cover the NBs, resulting in enhancement of NB stability in water

34,35.

The total organic carbon (TOC-VCHP, Shimadzu,

Co., Ltd., Japan) of a NB suspension was measured to quantify any residual oil contamination. The measured TOC values for the NB suspension were below the lower limit of quantification. The density of NBs in this study was therefore assumed to have the same density as air for calculations of bubble diameter. Glass beads (Toshin Riko Co., Ltd., Japan) composed of silicon dioxide (70–73%), sodium oxide and potassium oxide (13–16%), calcium oxide (7–12%), magnesium oxide (1– 4%), aluminum oxide (0.5–1.5%), ferric oxide (0.1–0.15%), and residuals ( 0.85). Values of 𝑘𝑑𝑒𝑝 increased with the IS at the same pH and electrolyte solution (i.e., NaCl or CaCl2), again representing the reduced repulsive force between NBs and glass beads. In addition, the higher value of 𝑆𝑚𝑎𝑥 in the presence of Ca2+ at IS=10 mM and pH=8 or at IS= 1 mM and pH= 6 as compared to the ones in the presence of Na+ with the same IS and pH indicate more locations that were favorable for NB retention, likely produced by cation bridging or charge neutralization or reversal 42–44.

17 ACS Paragon Plus Environment

Langmuir 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

Page 18 of 37

Table 1. Fitted parameters for NBs deposition and release models. pH

Electrolyte in NB water

IS

NB deposition Smax/C0

kdep

mM

m3 kg-1

NB release R2

fnr

ksw

R2

min-1

 

 

min-1

 

8

NaCl (Case 1)

10

1.50E-03

0.03

0.99

0.85

0.17

0.98

8

NaCl

50

3.45E-03

0.13

0.98

 

 

 

8

NaCl

100

3.53E-02

0.41

0.85

 

 

 

8

CaCl2

10

9.93E+01

0.48

0.87

 

 

 

8

NaCl (Case 2)

10

2.82E-03

0.15

0.97

n.d.

n.d.

 

8

CaCl2 (Case 3)

10

7.31E-01

0.69

0.92

0.92

0.08

0.90

6

NaCl

1

8.05E-03

0.09

0.98

 

 

 

6

NaCl

5

1.32E+01

0.23

0.91

 

 

 

6

NaCl

10

4.60E+01

0.42

0.85

 

 

 

6

CaCl2

1

4.00E-02

0.26

0.92

 

 

 

6

CaCl2

10

n.d.

n.d.

n.d.

 

 

 

C0 – input bubble concentration

18 ACS Paragon Plus Environment

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

Langmuir

Figure 4. Total interaction energies between NBs and glass surface at different electrolyte solutions with IS = 10 or 100 mM (a) without considering hydrophobic interaction and (b) with hydrophobic interaction. Blue lines in (b) show the total interaction energies when considering surface roughness with a roughness height of 38 nm and different roughness fraction (fr). kB is the Boltzmann constant and Tk is the absolute temperature taken as 293 K.

The calculated interaction energy for NBs in NaCl and CaCl2 electrolyte solutions of different IS (10 and 100 mM) at pH 8 conditions is shown in Fig. 4. When the hydrophobic interaction was not considered (Fig. 4a) a repulsive force was prominent, especially in the presence of 100 mM NaCl and CaCl2 (IS=10 and 100 mM). This is because the electrostatic and van der Waals interactions were both repulsive. The van der Waals interaction was repulsive due to the negative effective Hamaker constant in the system of water-NBs (i.e., air)-glass beads. These interaction energy calculations are clearly not consistent with the experimental observations of significant NB retention in Figure 3. In contrast to Figure 4a, the interaction energy profiles in Figure 4b exhibit a distinct primary minimum for all conditions when van der

19 ACS Paragon Plus Environment

Langmuir 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 37

Waals, double layer, and hydrophobic interaction are considered. The primary minimum for IS=10 mM NaCl was deepest in all conditions but a significant energy barrier (~ 18 kBTK) was observed. Shen et al.

45

reported that few colloids can transport over the energy barrier and be

deposited in the primary minimum when the energy barrier is larger than 10 kBTK. Therefore, NBs are not expected to be attached in a primary minimum on the glass surface under this condition. In contrast, our experiments showed attachment of NBs under this condition (i.e., peak relative turbidity was around 0.8, Fig. 3a). Discrepancies between theoretical interaction energy predictions and experimental observations have also been found in numerous previous studies 46. One possible explanation is due to the influence of nanoscale roughness and/or chemical heterogeneity on glass bead surfaces. Glass has been demonstrated to exhibit significant amounts of nanoscale roughness. For example, Rasmuson et al.47 reported a root mean square roughness height of 38 nm for a glass microscope slide that was cleaned with sodium hydroxide. Nanoscale roughness can have a profound influence on the interaction energy profile. The energy barrier height and the depth of the primary minimum can be reduced or eliminated depending on the local scale roughness parameters

24,30,48.

For example, Table S1 shows the predicted influence of different nanoscale roughness fractions with a height of 38 nm on interaction energy parameters when considering van der Waals, double layer, and hydrophobic interactions for the various solution chemistry conditions at pH=8. Figure 4b also

20 ACS Paragon Plus Environment

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

Langmuir

showed reduced energy barrier height at IS=10 mM NaCl when considering surface roughness. Note that nanoscale roughness can create locations that are favorable and unfavorable for NB retention, and that the influence of roughness changes with the solution chemistry30. This can partially explain the observed NB retention in Figure 3. In addition to nanoscale roughness, glass beads that were mainly composed of SiO2 also

contain various metal oxides on their surface that produces charge heterogeneity. At the solution pH maintained during experimentation (pH=8), the metal oxides will be near-neutrally charged or carry a slightly positive charge in comparison to the bulk SiO2 surface which will be highly negatively charged

49.

Hence, these surface charge heterogeneities are expected to provide

favorable deposition sites on what is otherwise an unfavorable surface for deposition

13,24,50.

Adsorbed Ca2+ can also produce cation bridging, and charge neutralization and/or reversal 42–44. This can likely explain the enhanced retention of NBs in the presence of Ca2+ in comparison to Na+ for the same IS and pH conditions. However, it should be mentioned that nanoscale roughness has been demonstrated to control the interaction energy profile when considering both nanoscale roughness and charge heterogeneity simultaneously 24. Figure 5 shows bubble size distributions of applied NBs and effluent NBs at the time corresponding to peak effluent turbidity for each experiment, except for conditions with negligible amounts of NBs in the effluent. At both pH conditions, peak bubble size in the

21 ACS Paragon Plus Environment

Langmuir 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 37

effluents was similar to the one in the applied NB suspension. However, there was a tendency for larger NB sizes to be reduced in the effluent in comparison to the influent with increasing IS of NaCl at each pH condition (i.e., NBs larger than 250 nm for IS=50 mM at pH 8, and NBs larger than 210 nm for IS=5 mM at pH 6). Similar bubble size distribution at 1 mM of CaCl2 electrolyte solution was observed with one at 5 mM of NaCl electrolyte solution. These results suggest that NB attachment to the glass beads are dependent on the bubble size and larger NBs are more favorable for attachment, especially at higher IS and CaCl2 electrolyte solution. However, filtration theory predicts that smaller colloids will be more efficiently removed from the liquid phase because of greater rates of mass transfer to the solid phase 51. Therefore, this discrepancy might be attributed to physical trapping caused by surface straining at locations with macroscopic roughness or grain-grain contacts retention has been also reported in Hamamoto et al.

52.

11

Bubble size dependency on the NBs

where higher retention of oxygen NBs

with larger size was obtained especially at lower pH.

22 ACS Paragon Plus Environment

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

Langmuir

Figure 5. Bubble size distributions of applied and effluent NB water at (a) pH 8 and (b) pH 6. The values in the legends represent ionic strength.

Figure 6a shows the breakthrough curves of NBs at IS=10 mM (pH 8) with different electrolytes types for the NB suspension and the control (resident and eluting) solution as follows: NB suspension and control solution used NaCl (Case1); NB suspension used NaCl and control solution used CaCl2 (Case 2); and NB suspension used CaCl2 and control solution used NaCl (Case 3). Figure 6b shows ion concentrations of Na+ and Ca2+ in the effluents of Cases 2 and 3. In Case 2, the peak relative turbidity was highly reduced as compared to Case 1, suggesting introduction of CaCl2 influenced retention of the subsequently applied NB suspension with NaCl. As shown in Fig. 6b, relative Ca2+ and Na+ concentrations during NB injection were not completely 0 or 1, respectively, indicating the Ca2+ in both bulk and exchangeable phases was not fully replaced by Na+. Therefore, residual Ca2+ in the phases contributed to enhanced NB retention. The Ca2+ also worked as a cation bridge (ionic bonds)

23 ACS Paragon Plus Environment

Langmuir 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

between NBs and glass surface

37,53.

In Case 3, almost negligible NB effluents were obtained.

The repulsive force between NBs and glass beads was reduced when Na+ in bulk and exchangeable phases was replaced by Ca2+ during NB water injection (Fig. 6b) and this enhanced NB retention. Even after switching to the control solution with NaCl, detachment of NBs was not observed even though the interaction between NBs and glass surface was more unfavorable. The NB transport model well captured the breakthrough curves in Cases 2 and 3 (Fig. 6a), and higher values of 𝑘𝑑𝑒𝑝 and 𝑆𝑚𝑎𝑥 were obtained in the presence of 10 mM Ca2+ (Case 3) (Table 1). The value of 𝑆𝑚𝑎𝑥 was larger in Case 2 in comparison to Case 1. This indicates an increase in the number of favorable sites for NB retention on the glass beads due to introduction

of CaCl2 prior to NB injection.

Figure 6. (a) Breakthrough curves for NBs at experiments with different combinations of control solutions and NB water at IS = 10 mM, (b) Ion concentrations of Na+ and Ca2+ in the effluents. The dotted vertical line represents pore volume when starting injection of control solution. Solid lines in Fig. 6a represent model simulations with Eqs. (1)-(3).

24 ACS Paragon Plus Environment

Page 24 of 37

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

Langmuir

Figure 7. Bubble concentrations and relative electrical conductivity after DI water injection. Solid lines represent model simulations with Eqs. (5)-(7).

Figure 7 shows the bubble concentration and relative EC after the DI water injection following the application of the NB water at IS of 10 mM (i.e., Cases 1, 2, 3 in Fig. 6a). The EC dramatically decreased after approximately one pore volume, indicating that the fluid throughout the whole column was replaced by DI water from the control solution. A decrease of EC in the effluent was accompanied by a rapid increase in the NB concentration for Cases 1 and 3, while no increase in NB concentration was observed for Case 2. The NB concentrations for Cases 1 and 3 subsequently decreased with increasing pore volumes. Thus, a decrease of IS in the aqueous phase enhanced the detachment of NBs that were trapped inside the column during the injection of NB water for Cases 1 and 3. Diffusive release of attached NBs with an IS decrease is expected to occur when the energy barrier to detachment decreases below a critical threshold of about 5.7 to 7 54. On smooth surfaces and in the presence of chemical heterogeneity the energy barrier to detachment is actually expected to increase with IS

25 ACS Paragon Plus Environment

55.

Conversely, the

Langmuir 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

energy barrier to detachment has been demonstrated to decrease with decreasing IS for certain roughness conditions

55.

Case 3 showed lower peak bubble concentration and more gradual

decrease after the peak as compared to Case 1. Similarly, no NB effluents were observed for the DI water injection experiments using NB water at IS of 10 mM CaCl2 (Fig. 3a). These observations likely reflect an increasing role of charge heterogeneity and/or cation bridging in the presence of adsorbed Ca2+, which would increase the energy barrier to detachment. Thus, existence of Ca2+ highly reduced the mobility of NBs in the porous media. The release model (Eqs. (5)-(7)) during DI water injection provided a good description of NB release for Cases 1 and 3. The values of 𝑓𝑛𝑟 were similar for Case 1 and 3, while a lower value of 𝑘𝑠𝑤 was obtained for Cased 3. This also reflects the influence of adsorbed Ca2+ on the reduction of NB release due to

the presences of localized favorable sites.

CONCLUSIONS In this study, effects of solution chemistry on the transport, retention, and release of NBs in glass beads were investigated. Overall, the results from this study revealed that with increasing IS and decreasing pH the mobilities of NBs were decreased due to reduced repulsive force between NBs and glass surface. The NB retention was explained by attractive hydrophobic interaction and reduced energy barrier caused by nanoscale surface roughness and surface charge heterogeneities. Divalent ion, Ca2+ in the electrolyte solution and adsorbed Ca2+ 26 ACS Paragon Plus Environment

Page 26 of 37

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

Langmuir

onto glass beads, highly reduced the mobility of NBs in the porous media because of induced charge neutralization or reversal, and/or cation bridging.

ASSOCIATED CONTENT

Supporting Information

Bubble size distribution before freezing and after thawing of NB suspensions, schematic illustration on the experimental setup, and details of DLVO interaction energy calculations are provided as supporting information.

AUTHOR INFORMATION

Corresponding Author *[email protected]

Author Contributions The manuscript was written through contributions of all authors. All authors have given approval to the final version of the manuscript.

ACKNOWLEDGEMENTS This work was supported by Grant-in-Aid for Scientific Research of Japan Society for the

27 ACS Paragon Plus Environment

Langmuir 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

Promotion of Science (JSPS) (Leading Initiative for Excellent Young Researchers, and 16H04411), and Grant-in-Aid for Exploratory Research for Research Institute for Sustainable Humanosphere, Kyoto University. We would like to especially thank Prof. Ohshita and Assoc. Prof. Nihei from the Graduate School of Agricultural and Life Sciences, The University of Tokyo for their assistance to zeta potential measurements and NBs water production.

ABBREVIATIONS NB, Nano-Bubble; IS, Ionic Strength; DI, Deionized Water

REFERENCES (1)

Burns, S. E.; Yiacoumi, S.; Tsouris, C. Microbubble Generation for Environmental and Industrial Separations. Sep. Purif. Technol. 1997, 11 (3), 221–232. https://doi.org/10.1016/S1383-5866(97)00024-5.

(2)

Khuntia, S.; Majumder, S. K.; Ghosh, P. Microbubble-Aided Water and Wastewater Purification: A Review. Rev. Chem. Eng. 2012, 28 (4–6), 191–221. https://doi.org/10.1515/revce-2012-0007.

(3)

Li, P.; Tsuge, H. Ozone Transfer in a New Gas-Induced Contactor with Microbubbles. J. Chem. Eng. Japan 2006, 39 (11), 1213–1220. https://doi.org/10.1252/jcej.39.1213.

(4)

Park, J.; Kurata, K. Application of Microbubbles to Hydroponic Solution Promotes Lettuce Growth. HortThecnology 2009, 19 (December), 212–215. 28 ACS Paragon Plus Environment

Page 28 of 37

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

Langmuir

(5)

Choi, Y. J.; Park, J. Y.; Kim, Y. J.; Nam, K. Flow Characteristics of Microbubble Suspensions in Porous Media as an Oxygen Carrier. Clean - Soil, Air, Water 2008, 36 (1), 59–65. https://doi.org/10.1002/clen.200700146.

(6)

Jauregi, P.; Varley, J. Colloidal Gas Aphrons: Potential Applications in Biotechnology. Trends Biotechnol. 1999, 17 (10), 389–395. https://doi.org/10.1016/S0167-7799(99)01363-3.

(7)

Jenkins, K. B.; Michelsen, D. L.; Novak, J. T. Application of Oxygen Mirobubbles for in Situ Biodegradation of P-Xylene-Contaminated Groundwater in a Soil Column. Biotechnol. Prog. 1993, 9 (4), 394–400. https://doi.org/10.1021/bp00022a006.

(8)

Hu, L.; Xia, Z. Application of Ozone Micro-Nano-Bubbles to Groundwater Remediation. J. Hazard. Mater. 2018, 342, 446–453. https://doi.org/10.1016/j.jhazmat.2017.08.030.

(9)

Wan, J.; Veerapaneni, S.; Gadelle, F.; Tokunag. Generation of Stable Microbubbles and Their Transport through Porous Media. Water Resour. Res. 2001, 37 (5), 1173–1182.

(10)

Hamamoto, S.; Nihei, N.; Ueda, Y.; Moldrup, P.; Nishimura, T. Effects of Flow Rate and Gas Species on Microbubble and Nanobubble Transport in Porous Media. J. Environ. Eng 2017, 143 (7), 1–7. https://doi.org/10.1061/(ASCE)EE.1943-7870.0001216.

(11)

Hamamoto, S.; Takemura, T.; Suzuki, K.; Nishimura, T. Effects of PH on Nano-Bubble Stability and Transport in Saturated Porous Media. J. Contam. Hydrol. 2018.

29 ACS Paragon Plus Environment

Langmuir 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

https://doi.org/10.1016/j.jconhyd.2017.12.001. (12)

Ryan, J. N.; Elimelech, M.; Ard, R. A.; Harvey, R. W.; Johnson, P. R. Bacteriophage PRD1 and Silica Colloid Transport and Recovery in an Iron Oxide-Coated Sand Aquifer. Environ. Sci. Technol. 1999, 33 (1), 63–73. https://doi.org/10.1021/es980350+.

(13)

Tufenkji, N.; Elimelech, M. Deviation from Colloid Filtration Theory in the Presence of Repulsive Electrostatic Interactions: Implications to Microbial Transport. Langmuir 2004, No. 20, 10818–10828.

(14)

Zhou, D.; Wang, D.; Cang, L.; Hao, X.; Chu, L. Transport and Re-Entrainment of Soil Colloids in Saturated Packed Column: Effects of PH and Ionic Strength. J. Soils Sediments 2011, 11 (3), 491–503. https://doi.org/10.1007/s11368-010-0331-2.

(15)

Verwey, E. J. W. Theory of the Stability of Lyophobic Colloids. J. Phys. Colloid Chem. 1947, 51 (3), 631–636. https://doi.org/10.1021/j150453a001.

(16)

Israelachvili, J. Intermolecular and Surface Forces; Acadmic Press, London, 2011. https://doi.org/10.1016/C2009-0-21560-1.

(17)

Khaled Abdella Ahmed, A.; Sun, C.; Hua, L.; Zhang, Z.; Zhang, Y.; Marhaba, T.; Zhang, W. Colloidal Properties of Air, Oxygen, and Nitrogen Nanobubbles in Water: Effects of Ionic Strength, Natural Organic Matters, and Surfactants. Environ. Eng. Sci. 2017, 35 (7), ees.2017.0377. https://doi.org/10.1089/ees.2017.0377.

30 ACS Paragon Plus Environment

Page 30 of 37

Page 31 of 37 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

Langmuir

(18)

Litton, G. M.; Olson, T. M. Colloid Deposition Rates on Silica Bed Media and Artifacts Related to Collector Surface Preparation Methods. Environ. Sci. Technol. 1993, 27 (1), 185–193. https://doi.org/10.1021/es00038a022.

(19)

Calgaroto, S.; Wilberg, K. Q.; Rubio, J. On the Nanobubbles Interfacial Properties and Future Applications in Flotation. Miner. Eng. 2014, 60, 33–40. https://doi.org/10.1016/j.mineng.2014.02.002.

(20)

Takahashi, M. ζ Potential of Microbubbles in Aqueous Solutions: Electrical Properties of the Gas−Water Interface. J. Phys. Chem. B 2005, 109 (46), 21858–21864. https://doi.org/10.1021/jp0445270.

(21)

Ishida, N. Direct Measurement of Hydrophobic Particle-Bubble Interactions in Aqueous Solutions by Atomic Force Microscopy: Effect of Particle Hydrophobicity. Colloids Surfaces A Physicochem. Eng. Asp. 2007, 300 (3 SPEC. ISS.), 293–299. https://doi.org/10.1016/j.colsurfa.2007.02.003.

(22)

Sakai, M.; Murata, T.; Kamio, K.; Mukae, K.; Yamauchi, A.; Moroi, Y.; Sugihara, G.; Norde, W. Investigation of Gas/Liquid Interface of Small Bubbles Formed in Solutions of Different Alkylammonium Chlorides. Colloids Surfaces A Physicochem. Eng. Asp. 2010, 359 (1–3), 6–12. https://doi.org/10.1016/j.colsurfa.2009.12.028.

(23)

Flury, M.; Aramrak, S. Role of Air-Water Interfaces in Colloid Transport in Porous

31 ACS Paragon Plus Environment

Langmuir 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

Media: A Review. Water Resour. Res. 2017, 53, 5247–5275. (24)

Bradford, S. A.; Torkzaban, S. Colloid Interaction Energies for Physically and Chemically Heterogeneous Porous Media. Langmuir 2013, 29 (11), 3668–3676. https://doi.org/10.1021/la400229f.

(25)

Wang, Y.; Li, Y.; Pennel, D, K. Influence of Electrolyte Species and Concentration on the Aggregation and Transport of Flluerene Nanoparticles in Quartz Sands. Environ. Toxicol. Chem. 2008, 27 (9), 1860–1867.

(26)

Jaisi, D.; Elimelech, M. Single-Walled Carbon Nanotubes Exhibit Limited Transport in Soil Columns. Env. Sci Technol 2009, 43 (24), 9161–9166.

(27)

Shen, C.; Lazouskaya, V.; Jin, Y.; Li, B.; Ma, Z.; Zheng, W.; Huang, Y. Coupled Factors Influencing Detachment of Nano- and Micro-Sized Particles from Primary Minima. J. Contam. Hydrol. 2012, 134–135, 1–11. https://doi.org/10.1016/j.jconhyd.2012.04.003.

(28)

Bradford, S. A.; Torkzaban, S.; Kim, H.; Simunek, J. Modeling Colloid and Microorganism Transport and Release with Transients in Solution Ionic Strength. Water Resour. Res. 2012, 48 (9), 1–10. https://doi.org/10.1029/2012WR012468.

(29)

Xu, S.; Qi, J.; Chen, X.; Lazouskaya, V.; Zhuang, J.; Jin, Y. Coupled Effect of Extended DLVO and Capillary Interactions on the Retention and Transport of Colloids through Unsaturated Porous Media. Sci. Total Environ. 2016, 573, 564–572.

32 ACS Paragon Plus Environment

Page 32 of 37

Page 33 of 37 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

Langmuir

https://doi.org/10.1016/j.scitotenv.2016.08.112. (30)

Bradford, S. A.; Kim, H.; Shen, C.; Sasidharan, S.; Shang, J. Contributions of Nanoscale Roughness to Anomalous Colloid Retention and Stability Behavior. Langmuir 2017, 33 (38), 10094–10105. https://doi.org/10.1021/acs.langmuir.7b02445.

(31)

Kobayashi, H.; Maeda, S.; Kashiwa, M.; Fujita, T. Measurement and Identification of Ultrafine Bubbles by Resonant Mass Measurement Method. Proc. SPIE, Aya, N.; Iki, N.; Shimura, T.; Shirai, T., Eds 2014, 9232, 92320S.

(32)

Jin, J.; Feng, Z.; Yang, F.; Gu, N. Bulk Nanobubbles Fabricated by Repeated Compression of Microbubbles. Langmuir 2019, 35 (12), 4238–4245. https://doi.org/10.1021/acs.langmuir.8b04314.

(33)

Nirmalkar, N.; Pacek, A. W.; Barigou, M. On the Existence and Stability of Bulk Nanobubbles. Langmuir 2018, 34, 10964–10973. https://doi.org/10.1021/acs.langmuir.8b01163.

(34)

Yasui, K.; Tuziuti, T.; Kanematsu, W.; Kato, K. Dynamic Equilibrium Model for a Bulk Nanobubble and a Microbubble Partly Covered with Hydrophobic Material. Langmuir 2016, 32 (43), 11101–11110. https://doi.org/10.1021/acs.langmuir.5b04703.

(35)

Sugano, K.; Miyoshi, Y.; Inazato, S. Study of Ultrafine Bubble Stabilization by Organic Material Adhesion. Japanese J. Multiph. Flow 2017, 31 (3), 299–306.

33 ACS Paragon Plus Environment

Langmuir 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

https://doi.org/10.3811/jjmf.31.299. (36)

Bradford, S. A.; Torkzaban, S.; Shapiro, A. A Theoretical Analysis of Colloid Attachment and Straining in Chemically Heterogeneous Porous Media. Langmuir 2013, 29 (23), 6944–6952. https://doi.org/10.1021/la4011357.

(37)

Torkzaban, S.; Bradford, S. A.; Wan, J.; Tokunaga, T.; Masoudih, A. Release of Quantum Dot Nanoparticles in Porous Media: Role of Cation Exchange and Aging Time. Environ. Sci. Technol. 2013, 47 (20), 11528–11536. https://doi.org/10.1021/es402075f.

(38)

Bradford, S. A.; Torkzaban, S.; Leij, F.; Simunek, J. Equilibrium and Kinetic Models for Colloid Release under Transient Solution Chemistry Conditions. J. Contam. Hydrol. 2015, 181, 141–152. https://doi.org/10.1016/j.jconhyd.2015.04.003.

(39)

Wang, Y.; Li, Y.; Pennell, K. D. Influence of Electrolyte Species and Concentration on the Aggregation and Transport of Fullerene Nanoparticles in Quartz Sands. Environ. Toxicol. Chem. 2008, 27 (9), 1860–1867. https://doi.org/10.1897/08-039.1.

(40)

Elimelech, M.; O’Melia, C. R. Kinetics of Deposition of Colloidal Particles in Porous Media. Environ. Sci. Technol. 1990, 24 (10), 1528–1536. https://doi.org/10.1021/es00080a012.

(41)

Ryan, J. N.; Gschwend, P. M. Effects of Ionic Strength and Flow Rate on Colloid Release: Relating Kinetics to Intersurface Potential Energy. Journal of Colloid And

34 ACS Paragon Plus Environment

Page 34 of 37

Page 35 of 37 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

Langmuir

Interface Science. 1994, pp 21–34. https://doi.org/10.1006/jcis.1994.1139. (42)

Torkzaban, S.; Wan, J.; Tokunaga, T. K.; Bradford, S. A. Impacts of Bridging Complexation on the Transport of Surface-Modified Nanoparticles in Saturated Sand. J. Contam. Hydrol. 2012, 136–137, 86–95. https://doi.org/10.1016/j.jconhyd.2012.05.004.

(43)

Kloster, N.; Brigante, M.; Zanini, G.; Avena, M. Aggregation Kinetics of Humic Acids in the Presence of Calcium Ions. Colloids Surfaces A Physicochem. Eng. Asp. 2013, 427, 76–82. https://doi.org/10.1016/j.colsurfa.2013.03.030.

(44)

Majzik, A.; Tombácz, E. Interaction between Humic Acid and Montmorillonite in the Presence of Calcium Ions II. Colloidal Interactions: Charge State, Dispersing and/or Aggregation of Particles in Suspension. Org. Geochem. 2007, 38 (8), 1330–1340. https://doi.org/10.1016/j.orggeochem.2007.04.002.

(45)

Shen, C.; Li, B.; Huang, Y.; Jin, Y. Kinetics of Coupled Primary- and Secondary-Minimum Deposition of Colloids under Unfavorable Chemical Conditions. Environ. Sci. Technol. 2007, 41 (20), 6976–6982. https://doi.org/10.1021/es070210c.

(46)

Ryan, J. N.; Elimelech, M. Colloid Mobilization and Transport in Groundwater. Colloids Surfaces A Physicochem. Eng. Asp. 1996, 107 (95), 1–56. https://doi.org/10.1016/0927-7757(95)03384-X.

(47)

Rasmuson, A.; Pazmino, E.; Assemi, S.; Johnson, W. P. Contribution of Nano- to

35 ACS Paragon Plus Environment

Langmuir 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

Microscale Roughness to Heterogeneity: Closing the Gap between Unfavorable and Favorable Colloid Attachment Conditions. Environ. Sci. Technol. 2017, 51 (4), 2151– 2160. https://doi.org/10.1021/acs.est.6b05911. (48)

Shen, C.; Wang, F.; Li, B.; Jin, Y.; Wang, L. P.; Huang, Y. Application of DLVO Energy Map to Evaluate Interactions between Spherical Colloids and Rough Surfaces. Langmuir 2012, 28 (41), 14681–14692. https://doi.org/10.1021/la303163c.

(49)

Torkzaban, S.; Bradford, S. A. Critical Role of Surface Roughness on Colloid Retention and Release in Porous Media. Water Res. 2016, 88, 274–284. https://doi.org/10.1016/j.watres.2015.10.022.

(50)

Song, L.; Johnson, P. R.; Elimelech, M. Kinetics of Colloid Deposition onto Heterogeneously Charged Surfaces in Porous Media. Environ. Sci. Technol. 1994, 28 (6), 1164–1171. https://doi.org/10.1021/es00055a030.

(51)

Yao, K.-M.; Habibian, M. T.; O’Melia, C. R. Water and Waste Water Filtration-Concepts and Applications. Environ. Sci. Technol. 1971, 5 (11), 1105–1112.

(52)

Bradford, S. A.; Torkzaban, S. Colloid Transport and Retention in Unsaturated Porous Media: A Review of Interface-, Collector-, and Pore-Scale Processes and Models. Vadose Zo. J. 2008, 7 (2), 667. https://doi.org/10.2136/vzj2007.0092.

(53)

Grosberg, A. Y.; Nguyen, T. T.; Shklovskii, B. I. Low Temperature Physics at Room

36 ACS Paragon Plus Environment

Page 36 of 37

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

Langmuir

Temperature in Water: Charge Inversion in Chemical and Biological Systems. Int. J. High Speed Electron. Syst. 2002, 12 (2), 235–265. https://doi.org/10.1142/S012915640200123X. (54)

Bradford, S. A.; Torkzaban, S. Determining Parameters and Mechanisms of Colloid Retention and Release in Porous Media. Langmuir 2015, 31 (44), 12096–12105. https://doi.org/10.1021/acs.langmuir.5b03080.

(55)

Shen, C.; Bradford, S. A.; Li, T.; Li, B.; Huang, Y. Can Nanoscale Surface Charge Heterogeneity Really Explain Colloid Detachment from Primary Minima upon Reduction of Solution Ionic Strength? J. Nanoparticle Res. 2018, 20 (6). https://doi.org/10.1007/s11051-018-4265-8.

37 ACS Paragon Plus Environment