Contribution of Nano- to Microscale Roughness ... - ACS Publications

Anna Rasmuson†, Eddy Pazmino‡, Shoeleh Assemi§, and William P. Johnson†. † Department of Geology and Geophysics, University of Utah, Salt Lake City, ...
0 downloads 0 Views 1MB Size
Subscriber access provided by UNIV OF NEBRASKA - LINCOLN

Article

The Contribution of Nano- to Micro-scale Roughness to Heterogeneity: Closing the Gap between Unfavorable and Favorable Colloid Attachment Conditions Anna Rasmuson, Eddy Pazmino, Shoeleh Assemi, and William P. Johnson Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.6b05911 • Publication Date (Web): 30 Jan 2017 Downloaded from http://pubs.acs.org on February 5, 2017

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

Environmental Science & Technology 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 40

Environmental Science & Technology

1

The Contribution of Nano- to Micro-scale Roughness to Heterogeneity:

2

Closing the Gap between

3

Unfavorable and Favorable Colloid Attachment Conditions

4

Anna Rasmuson1, Eddy Pazmino2, Shoeleh Assemi3, William P. Johnson1,*

5 6 7

1

8

Department of Geology and Geophysics, University of Utah,

9

Salt Lake City, Utah 84112, United States 2

10

Department of Extractive Metallurgy, Escuela Politécnica Nacional, Quito, Ecuador 3

11

Department of Metallurgical Engineering, University of Utah, Salt Lake City, Utah 84112, United States

12 13 14 15 16 17 18

*

Corresponding author. Email: [email protected]; Tel: (801)585-5033; Fax: (801)5817065.Engineering, University of Utah, Salt Lake City, Utah 84112, United States 1

ACS Paragon Plus Environment

Environmental Science & Technology

19

TOC/Abstract Art

20

2 ACS Paragon Plus Environment

Page 2 of 40

Page 3 of 40

Environmental Science & Technology

21

Abstract

22

Surface roughness has been reported to both increase as well as decrease colloid retention. In

23

order to better understand the boundaries within which roughness operates, attachment of a

24

range of colloid sizes to glass with three levels of roughness was examined under both

25

favorable (energy barrier absent) and unfavorable (energy barrier present) conditions in an

26

impinging jet system. Smooth glass was found to provide the upper and lower bounds for

27

attachment under favorable and unfavorable conditions, respectively. Surface roughness

28

decreased, or even eliminated, the gap between favorable and unfavorable attachment, and

29

did so by two mechanisms: 1) under favorable conditions attachment decreased via increased

30

hydrodynamic slip length and reduced attraction ; 2) under unfavorable conditions attachment

31

increased via reduced colloid-collector repulsion (reduced radius of curvature) and increased

32

attraction (multiple points of contact, and possibly increased surface charge heterogeneity).

33

Absence of a gap where these forces most strongly operate for smaller (< 200 nm) and larger (>

34

2 µm) colloids was observed and discussed. These observations elucidate the role of roughness

35

in colloid attachment under both favorable and unfavorable conditions.

3 ACS Paragon Plus Environment

Environmental Science & Technology

36

Introduction

37

Significant progress has been made in understanding the nature of nanoscale heterogeneity

38

responsible for colloid attachment to bulk repulsive surfaces. 1-4 Colloid-collector repulsion in

39

the environment often arises from like-charged colloids and collectors wherein counter-ions

40

that shield surface charge are “squeezed” between two approaching surfaces.5,6 This repulsion

41

is reduced or eliminated by zones of charge opposite to the colloid/collector surfaces that

42

create local zones of attraction.2,4,7,8 Because the magnitude of repulsion between like-charged

43

colloids and collectors scales directly with their radii of curvature, low radii of curvature of

44

nanoscale asperities (roughness) also locally diminish or eliminate repulsion between like-

45

charged surfaces.1,3,9 Therefore, both nanoscale charge heterogeneity and roughness may

46

counteract repulsion between like-charged surfaces.

47

Experimentally roughness has been shown to increase,10-13 as well as decrease14-16 attachment.

48

This apparent discrepancy may result from different influences of roughness under favorable

49

(energy barrier absent) versus unfavorable (energy barrier present) conditions. Limited studies

50

were conducted under favorable conditions15,17,19 under which the influence of roughness was:

51

1) non-monotonic15 showing a minimum collector efficiency (η = # colloids attached per #

52

colloids introduced) for intermediate roughness (~200 nm root mean square roughness); and 2)

53

monotonic, with η decreasing19 or increasing17 as roughness increased.

54

To understand the seemingly contradictory effects of roughness on particle attachment, its

55

influence needs to be understood via the following mechanisms: 1) interception of the surface,

56

and 2) arrest on the surface (attachment).20 Colloid interception of the surface via fluid drag

4 ACS Paragon Plus Environment

Page 4 of 40

Page 5 of 40

Environmental Science & Technology

57

(fluid streamlines) is enhanced by particle settling and diffusion as encapsulated in colloid

58

filtration theory (CFT).21 Additionally, roughness may diminish interception when the near-

59

surface colloid velocities are increased by the slip layer which is the zone of fluid shear that

60

exists between asperities 22,23 (Figure 1). Roughness may enhance interception18 via local

61

protrusion of the surface into the pore domain, and by decreasing hydrodynamic drag in

62

regions between these asperities.2

63 64

Colloid arrest is governed by the balance of mobilizing and arresting torques, emanating

65

primarily from fluid drag and colloid-collector interaction forces, respectively. 4,21,26,27,31

66

Roughness influences arresting torque in several ways: 1) By decreasing the magnitude of

67

colloid-collector interactions3,9,13,24,25 such that net repulsion (unfavorable conditions) and net

68

attraction (favorable conditions) are decreased since electric double layer (EDL) and van der

69

Waals (VdW) interactions scale to the local radius of curvature of the interacting surfaces;3,19 2)

70

By increasing colloid-collector contact area by establishing multiple points of contact and

71

thereby increasing the arresting torque; 3) By potentially creating additional lever arms

72

associated with contact,27 thereby increasing the arresting torque.

73

The above review demonstrates that roughness has contrasting influences on attachment

74

under unfavorable versus favorable conditions. Experiments performed under these contrasting

75

conditions should help to elucidate these effects of roughness. For this reason, we examined

76

attachment under favorable versus unfavorable conditions for variably rough surfaces.

77

Attachment was examined in pore-scale experiments (impinging jet) where the mode of

5 ACS Paragon Plus Environment

Environmental Science & Technology

78

attachment was directly observed, and motion was tracked in a single focal plane. In contrast,

79

column scale observations do not elucidate the specific modes of attachment, which include

80

immobilization on open surface, wedging in grain-to-grain contacts, and retention without

81

attachment,288 nor do porous media experiments (even direct observation micromodels) allow

82

motion tracking over significant distances because the curved surfaces cross focal planes.

83

Direct observation of attachment on the planar surface therefore allows elucidation of the

84

influence of the hydrodynamic slip layer on near-surface colloid velocities and colloid

85

attachment. The planar surface also eliminates topological complications such as ripening via

86

funneling of colloids into aggregates in grain-to-grain contacts.29 Notably, the experiments of

87

Torkzaban and Bradford19 examined colloid retention in smooth versus rough porous media

88

under high IS conditions where favorable colloid-colloid interaction promote ripening. Such

89

interactions may confound elucidating the influence of roughness on observed colloid release

90

via ionic strength (IS) reduction, since colloid-colloid interactions also respond to IS reduction.

91

The impinging jet system represents only divergent flow at the forward flow stagnation zone,

92

and does not capture Happel sphere-in-cell flow convergence at the rear-flow stagnation zone21

93

representative of porous media. However, by capturing attachment on the open collector

94

surface, the impinging jet allows calibration of the discrete heterogeneity responsible for

95

attachment4,8 and detachment30 under unfavorable conditions. Upon calibrating the

96

contribution of roughness, the resulting representative heterogeneity can be incorporated into

97

other existing collector geometries, such as the Happel sphere-in-cell and Hemisphere-in-cell

98

models.31

6 ACS Paragon Plus Environment

Page 6 of 40

Page 7 of 40

Environmental Science & Technology

99

To ascertain the contribution of roughness, attachment of colloids ranging in diameter from 20

100

nm to 4.4 µm was examined on surfaces with root mean square (RMS) roughness ranging from

101

≤1 nm to 546 nm. This range of asperity/colloid size ratios spans from 1E-5 to 5 (Figure 1),

102

where the largest colloid size (4.4 µm) represents a practical upper limit for stably-suspended

103

non-buoyant colloids in groundwater. The smallest colloid size (20 nm) represents what was

104

possible to resolve optically in our experiments. This colloid size and the largest asperities (546

105

nm RMS) define the largest asperity/colloid size ratio examined here (~5). Above this

106

asperity/colloid size ratio it is reasonable to expect that flow field geometry will be

107

characterized by additional forward and rear flow stagnation zones where impingement and

108

retention may occur as demonstrated in porous media micromodels (e.g., Ausset and Keller32).

109

Because this becomes an issue of resolving a modified flow field (fluid streamlines) rather than

110

the specific influence of roughness on colloid-surface interactions and hydrodynamic slip, we

111

therefore consider larger asperity/colloid size ratios (e.g., > 5) as complex porous media for

112

which the path toward improved prediction is development of corresponding flow fields that

113

represent the associated forward and rear flow stagnation zones.

114

In summary, our goal was to determine the range over which roughness operates (short of

115

creating new impingement surfaces) by defining upper and lower boundaries of attachment

116

which we posit correspond to smooth surfaces under favorable (upper bound attachment) and

117

unfavorable (lower bound attachment) conditions. An additional objective was to understand

118

the influence of roughness not only on attachment but also on detachment.

119

Methods

7 ACS Paragon Plus Environment

Environmental Science & Technology

120

Collector Surfaces

121

Microscope soda lime glass slides and coverslips (Fisher Scientific, Inc.) were used as the

122

impinging surface in the cell, the coverslips were used for smaller colloids that required higher

123

magnification for optical resolution. Three levels of roughness were developed for the glass

124

surfaces by examining them as: 1) untreated; 2) after NaOH-treatment (5 N for 80 min at 90ºC);

125

and 3) after trace metal grade hydrofluoric acid (HF) treatment (27.6 N for 12 hours at 21ºC).

126

Each treatment was followed by extensive rinsing with MilliQ until the rinse solution pH was

127

6.5. The HF treated surface was then heated in air at 250ºC for 2 hours to remove any residual

128

surface impurities. All surfaces were cleaned prior to experiments via the SC-1 procedure. The

129

treatment methods that were used to produce roughness may also introduce defects and

130

uncoordinated atoms on asperities that can increase charge heterogeneity.34 However, it is

131

reasonable to expect that if roughness and charge heterogeneity co-vary on chemically-etched

132

mineral surfaces in the laboratory, then they also co-vary in aquifer media that has undergone

133

natural chemical and mechanical weathering. Our observations are keyed to roughness, since

134

charge heterogeneity is inferred, whereas roughness is measurable.

135

Roughness was measured with an atomic force microscope (model N9451A Agilent

136

Technologies; Santa Clara, CA) following SC-1 cleaning using contact mode in air with silicon

137

nitride probes (type DNP-S10; Bruker Nano, Inc.) with a nominal spring constant of 0.12 N/m.

138

Roughness was evaluated using SPIP software (Image Metrology; Hørsholm Denmark), and was

139

measured as the root mean square height (Figure 2). The minimum scan size for each surface

140

was determined by measuring the RMS roughness over a large area (~50 µm) and then

8 ACS Paragon Plus Environment

Page 8 of 40

Page 9 of 40

Environmental Science & Technology

141

sampling ten random small scan areas (5 µm) within that larger area. If the average RMS

142

roughness of the small scans approximated that of the larger area, the smaller area was

143

accepted and the process was repeated. A minimum of five random locations taken randomly

144

across the observation area of the jet (~500 µm) were used to obtain average values and

145

standard deviations of roughness parameters including: average roughness, root mean square

146

(RMS) roughness, max valley depth, and max peak height. Three dimensional images were

147

generated using Gwyddion software (Czech Metrology Institute, Brno, Czech Republic). The

148

average wavelengths (peak to peak distances) of primary (larger scale) and secondary (small-

149

scale) asperities were measured from five representative profiles (Supplementary Information

150

Figure SI-1) for each surface using the ISO 4287-1997 written standard (ISO; Geneva

151

Switzerland).

152 153

Microsphere Suspensions

154

Carboxylate-modified polystyrene latex (CML) fluorescent (λex = 505, λem = 515 nm)

155

microspheres of six sizes (0.02, 0.1, 0.25, 1.1, 2.0 (Molecular Probes, Inc., Eugene, OR), and 4.4

156

(Polysciences, Washington, PA) μm diameter) were used in the experiments. Colloid

157

suspensions were prepared from stock in relevant solution with concentrations ranging from 5

158

× 106 to 2 × 107 microspheres per milliliter. The microsphere suspension concentration was

159

determined via vacuum filtration of colloid solution (volume adjusted to ensure >20 CML per

160

view area) on 0.05 or 0.1 μm polycarbonate filters (Millipore) followed by averaging counts of

161

25 random observation areas using wide-field fluorescence for colloid illumination and scaling

9 ACS Paragon Plus Environment

Environmental Science & Technology

162

this average to the area of deposition on the filter. Suspension ionic strength (IS) was adjusted

163

using NaCl. Unfavorable solutions (IS 6 mM) were buffered with 2.2 mM MOPS (3-(N-

164

morpholino) propanesulfonic acid, 4-morpholinepropanesulfonic acid; Sigma-Aldrich Corp.)

165

with pH set to 8.0 using NaOH (0.5 M). Favorable solutions were set to pH 2.0 using HCL (1.3 M)

166

and IS 50 mM. Dissolution of the collector surface was not a concern at these pH values as the

167

solubility of amorphous silica is stable and very low at pH values below pH 9.35

168

The Péclet (Pe) and Reynolds (Re) numbers for each colloid size under the relevant fluid

169

velocities are included in the Supplementary Information (Supplementary Information, Table SI-

170

1). Pe numbers all exceed 50,000, corresponding to advection–dominated conditions, and Re

171

numbers range from 1.73 to 5.94, corresponding to laminar creeping flow conditions.

172

The CML electrophoretic mobility (EPM) was measured in suspensions using ζ-potential

173

analyzer (Mobiuζ, Wyatt Technology Corp., Santa Barbara, CA). CML ζ-potentials were

174

calculated from EPM via the Smoluchowski equation for large colloids (i.e. ≥ 100 nm diameter),

175

and the Huckel approximation when the colloid diameter was smaller than the Debye length

176

(i.e. colloids < 100 nm diameter36 (Supplementary Information Table SI-2). The glass slide EPM

177

was measured in the filtrate ( 0.95) was

228

required to quantify η. This initial slope of deposition (across the area of observation, Aobs) was

229

used to calculate the collector efficiency (η) via the following equation:

230

ߟ=



#ೌ೟೟ೌ೎೓೐೏ ቁ ೟೔೘೐ ಲ



#೔೙ೕ೐೎೟೐೏ ቁ ೟೔೘೐ ಲ

ೀಳೄ

಻ಶ೅

=

#ೌ೟೟ೌ೎೓೐೏ ೟೔೘೐

(1)

஼ೀ ொ

231

where CO is the injected concentration of colloids and Q is the flow rate of the fluid that enters

232

the cell (across the area of the jet, Ajet). The product COQ is equal to the number of particles

233

injected per unit time across Ajet. This expression was developed specifically for the impinging

234

jet geometry,4 therefore comparison of results with porous media is qualitative. However, note

235

that the conversion of attachment to η as opposed to flux does not change the relative values

236

since the influence of concentration and velocity is equivalent for η (efficiency) and flux (e.g.,

237

dimensionless parameters such as the Sherwood number).

238

Near-surface velocities were determined by tracking particle displacement at the near-surface

239

as a function of time from the images captured using constant streaming (time intervals were

240

less than 1 second between images). Velocities were averaged for 0.25, 1.1 and 2.0 µm colloids

241

for a minimum of 10 particles. Near-surface colloids were distinguished from bulk colloids using

242

TIRF-M (i.e. over the ≤ 1 nm and 38 nm RMS glass surfaces). For the 546 nm RMS surface, a

243

threshold velocity of ~400 µm/min was selected to distinguish particles in the bulk solution.

244

Near-surface velocities were determined under favorable conditions in order to eliminate

245

interactions with the secondary minimum. 13 ACS Paragon Plus Environment

Environmental Science & Technology

246 247

Favorable Simulations

248

A Lagrangian particle trajectory model developed for the smooth impinging jet system was used

249

to predict η under favorable conditions.4 The model accounts for fluid drag, hydrodynamic

250

retardation, gravity, diffusion, colloid−surface interaction forces, and virtual mass. A more

251

detailed description of the force and torque balances including fluid drag, hydrodynamic

252

retardation, gravity, diffusion, steric forces, and virtual mass is described in previous

253

publications. 4,29,31 Favorable conditions were set by assigning measured colloid and collector ζ

254

potentials for the pH 2 50 mM conditions (Supplementary Information, Table SI-2) in the model

255

input.

256

The complete experimental matrix corresponding to the range of collector surfaces (3), colloid

257

sizes (6), fluid velocities (2), IS (2), and pH (2) yielded 72 experiments without replicates and 144

258

experiments with complete replication. In order to produce a tractable set of experiments and

259

simulations, some experimental conditions were not run when the observed trends could be

260

discerned from bracketing experiments, whereas in other cases it was not possible to optically

261

resolve colloids under certain conditions. For example, particles ≤ 100 nm were not examined

262

on the 546 nm RMS surface due to inability to distinguish faint particles from light scattering.

263

Approximately 75 experiments (including replicates) were conducted. Replicates were

264

performed for all experiments conducted on the untreated and NaOH treated slides under the

265

1.7E-3 ms-1 velocity condition. The complete experimental matrix with replicates indicated is

266

shown the Supplementary Information (Supplementary Information Table SI-3).

267

14 ACS Paragon Plus Environment

Page 14 of 40

Page 15 of 40

Environmental Science & Technology

268

Results

269

Roughness

270

Roughness parameters measured using AFM (Figure 2) are summarized in Supplementary

271

Information (Supplementary Information Table SI-4), where untreated, NaOH-treated, and HF-

272

treated glass slides and coverslips showed progressively increasing roughness. The glass slides

273

and coverslips showed similar RMS values for a given treatment. Three levels of roughness

274

were generated on glass slides and coverslips: a) ≤ 1 nm RMS roughness (untreated); b) 38 nm

275

RMS roughness (NaOH-treated); and c) 546 nm RMS roughness (HF-treated). The maximum

276

valley depths and peak heights were approximately a factor of five times the RMS roughness for

277

the untreated and NaOH-treated silica, and approximately a factor of three times the RMS for

278

the HF-treated silica. The primary (large scale) wavelengths were approximately equal to the

279

RMS roughness, suggesting that the asperities were evenly spaced. The untreated glass surface

280

did not exhibit secondary roughness. The secondary wavelengths were similar for both NaOH-

281

treated and HF-treated surfaces (~ 3.5 nm), which indicates that secondary roughness did not

282

change as overall roughness increased, and that roughness appears to follow a fractal

283

relationship.

284

The maximum collector (impinging surface) asperity size was on the order of size of the largest

285

colloids, whereas the smallest collector asperity size was