Probing the Interaction Mechanism between Air ... - ACS Publications

Oct 18, 2017 - Lei Xie, Chen Shi, Xin Cui, Jun Huang, Jingyi Wang, Qi Liu, and Hongbo Zeng*. Department of Chemical and Materials Engineering, Univers...
0 downloads 0 Views 2MB Size
Subscriber access provided by READING UNIV

Article

Probing the Interaction Mechanism between Air Bubbles and Bitumen Surfaces in Aqueous Media Using Bubble Probe AFM Lei Xie, Chen Shi, Xin Cui, Jun Huang, Jingyi Wang, Qi Liu, and Hongbo Zeng Langmuir, Just Accepted Manuscript • DOI: 10.1021/acs.langmuir.7b02693 • Publication Date (Web): 18 Oct 2017 Downloaded from http://pubs.acs.org on October 23, 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.

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

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

Probing the Interaction Mechanism between Air Bubbles and Bitumen Surfaces in Aqueous Media Using Bubble Probe AFM

Lei Xie, Chen Shi, Xin Cui, Jun Huang, Jingyi Wang, Qi Liu, Hongbo Zeng* Department of Chemical and Materials Engineering, University of Alberta, Edmonton, Alberta, T6G 1H9, Canada *Email: [email protected], Phone: +1-780-492-1044, Fax: +1-780-492-2881

ACS Paragon Plus Environment

1

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 2 of 33

ABSTRACT Surface interactions involving deformable air bubbles have attracted tremendous interest in a wide range of engineering applications, such as mineral flotation and bitumen extraction. In this work, for the first time, the interaction forces between air bubbles and bitumen surfaces in complex aqueous media of varying pH, salinity and salts were directly measured using a bubble probe atomic force microscope (AFM) technique. The AFM topographic imaging reveals that bitumen surface tends to be rougher and form distinct domains at high NaCl concentration or under strongly alkaline environment. The force measurements demonstrate the critical role of ionic strength and solution pH in bubble-bitumen interaction and attachment, which could be well described by a theoretical model based on Reynolds lubrication theory and augmented Young-Laplace equation by including the effect of disjoining pressure. In 1 mM NaCl, the electrical double layer (EDL) repulsion inhibited bubble-bitumen attachment, and such repulsive effect could be further strengthened with increasing solution pH. In 500 mM NaCl, the hydrophobic attraction could lead to bubble-bitumen attachment, while a high solution pH could weaken the hydrophobic interaction. The addition of calcium ion in 500 mM NaCl could enhance the hydrophobic interaction and facilitate the bubble-bitumen attachment, most likely attributed to the bridging effect between calcium ions and the functional groups (e.g., carboxyl group) of interface-active molecules on bitumen surfaces thus leading to higher surface roughness and hydrophobic moieties/aggregates on bitumen as confirmed by AFM imaging. Our results provide quantitative information on the interaction mechanism between air bubbles and bitumen surfaces in complex aqueous solutions at the nanoscale, which has useful implications on many related interfacial interactions in industrial processes such as oil production, oil-water separation and wastewater treatment.

ACS Paragon Plus Environment

2

Page 3 of 33

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

1. INTRODUCTION The colloidal interaction involving deformable air bubbles has attracted tremendous interest from both fundamental and practical viewpoints, for its ubiquitous presence in a wide range of biological and technological applications,1-3 such as froth flotation,4 drug and gene delivery,5 and microfluidic devices.6 Air bubbles, due to the intrinsic hydrophobicity, can readily attract hydrophobic or partially hydrophobic materials in an aqueous medium and collect them at the air/water interface.7-9 Therefore, materials with different degrees of surface hydrophobicity can be selectively separated on the basis of their attachment propensities to air bubbles, which is the fundamental principle of industrial froth flotation process.4,

10-11

The bubble-surface

attachment is governed by the drainage process of the intervening water film under the combined influences of hydrodynamic pressure and disjoining pressure attributed to surface interactions.1214

Generally, repulsion arising from van der Waals (vdW) and electrical double layer (EDL)

interactions can impede thin water film drainage and prevent bubble-surface attachment, whereas strong hydrophobic attraction can lead to the rupture of thin water film and bubble-surface attachment. The Athabasca oil sands located in Alberta, Canada are the third largest oil reserves in the world, and water-based bitumen extraction is currently used to recover bitumen from mined oil sands.15 During this process (illustrated in Figure 1A), bitumen is aerated to form bitumen-air aggregates, floating to the top of slurry in the form of bitumen-rich froth. The collected bitumen froth is treated to provide clean bitumen products, in which extremely stable emulsions are commonly encountered but generally undesirable.15-17 Recent studies have revealed that bubblebitumen attachment during aeration is one of the most critical steps to determine the bitumen recovery and product quality.15 Hence, a quantitative understanding of the interaction mechanism

ACS Paragon Plus Environment

3

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 33

and thin film drainage process between air bubble and bitumen surface under various aqueous conditions, which determines the formation of bitumen-air aggregates, is of both fundamental and practical importance to elucidate the bubble-bitumen attachment mechanism and optimize the current water-based bitumen extraction process. Over the last few years, significant progresses have been made to unravel the interaction mechanism involved in the water-based bitumen extraction process.18-21 Nanomechanical tools such as surface forces apparatus (SFA) and atomic force microscope (AFM) have been widely used to quantitatively measure the interaction forces of solid surfaces.22-26 Quantitative force measurements involving bitumen and asphaltenes (an interface-active component of bitumen) showed that the measured forces were dependent on aqueous conditions (e.g., pH, salinity, and salts), organic solvent and temperature.20-21,

27-28

Nevertheless, direct measurement of surface

forces between air bubbles and bitumen surfaces at the nanoscale with theoretical analysis of the associated thin film drainage process have not been reported, probably due to the practical difficulties in precise manipulation of deformable bubbles and interpretation of the measured results.29-31 Recently, the bubble/drop probe AFM has been developed and applied to study the interaction mechanisms involving deformable gas bubbles and emulsion drops, demonstrating the combinatory effects of hydrodynamic interaction, surface forces and bubble/drop deformation on the interaction.8,

11, 32-41

The measured force results could be successfully

described by a theoretical model based on Reynolds lubrication theory and augmented YoungLaplace equation. The film drainage process during interaction was also directly visualized by simultaneously measuring the force and separation using the bubble probe AFM combined with reflection interference contrast microscopy (RICM), which has validated the above theoretical

ACS Paragon Plus Environment

4

Page 5 of 33

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

model.32 Using the bubble/drop probe AFM technique, various surface forces at air/water and oil/water interfaces were precisely quantified, such as structural and hydrophobic interactions, with sub-nN resolution.14, 42-43 One of the most recent applications of this bubble/drop probe AFM technique was to measure the interaction mechanism of oil drops in water in the presence of asphaltenes, indicating that asphaltenes stabilized the oil-in-water emulsion by strengthening the repulsive EDL interaction and inducing steric repulsion.36 In this study, the bubble probe AFM technique was employed to quantitatively measure the interaction forces between air bubble and bitumen surface in complex aqueous media, and the measured force results were analyzed using the theoretical model based on Reynolds lubrication theory and augmented Young-Laplace equation. The surface energy of bitumen was measured using three-probe-liquid method, and the effects of solution conditions (e.g., ionic strength and solution pH) on the morphology of bitumen surface were also investigated. Our results for the first time provide valuable quantitative information on the interaction mechanism between air bubbles and bitumen surfaces at the nanoscale, with useful implications for related interfacial processes in heavy oil production and wastewater treatment.

ACS Paragon Plus Environment

5

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 33

Figure 1. (A) Illustration of bubble-bitumen-solid interactions in complex aqueous medium in water-based bitumen extraction process to recover bitumen from mined oil sands. (B) Schematic of force measurements between air bubble and bitumen surface in an aqueous solution using the bubble probe AFM technique. It is noted that the bubble size typically ranges from tens of µm to several mm in industrial bitumen extraction process (Figure 1A),10 and the typical bubble size in AFM force measurements can be varied from tens to hundreds of um (Figure 1B).

2. MATERIALS AND METHODS 2.1. Materials. Sodium chloride (NaCl, ACS reagent grade) and calcium chloride (CaCl2, ACS reagent grade) were purchased from Fisher Scientific and used as received without further purification. All aqueous solutions were prepared using Milli-Q water (Millipore deionized, 18.2 MΩ·cm resistivity). Hydrochloric acid (HCl, ACS reagent grade) and sodium hydroxide (NaOH, ACS reagent grade) purchased from Fisher Scientific were used to adjust solution pH. 2.2. Preparation of Bitumen Surface. Bitumen solution was prepared by adding bitumen (provided by Shell) into toluene at a concentration of 3 mg/mL and then centrifuging at 7000 rpm (or ~600 times the force of gravity) for 30 min to remove possible fine solids. After filtration with 0.2 µm PTFE membrane, the bitumen solution was spin coated on silica wafers at 2000 rpm for 40 s to ensure the full and uniform bitumen coverage on silica wafers; thereafter, the prepared bitumen surface was dried overnight to remove the remaining toluene. 2.3. Surface Energy Measurement. Three-probe-liquid method was used to determine the surface energy of bitumen. A contact angle goniometer (ramé-hart instrument Co., NJ, USA) was used to measure the contact angles of three probing liquids, including one nonpolar (i.e., diiodomethane) and two polar (i.e., water and ethylene glycol) liquids, on bitumen surface using a sessile drop method. For the same type of liquids, at least two different surfaces and three different positions on each surface were tested, and the average contact angle was reported.

ACS Paragon Plus Environment

6

Page 7 of 33

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

Surface energy γ with Lifshitz-vdW γ LW and Lewis acid-base γ AB (electron acceptor γ + and electron donor γ − ) components can be derived by a method developed by van Oss et al. in eq 1.44-45

γ = γ LW + γ AB = γ LW + 2 γ +γ −

(1)

The relation between liquid contact angle on bitumen surface θ and surface energy components of bitumen surface and liquid can be given by eq 2, where the subscript B or L represents bitumen or liquid, respectively.44-45

γ L (cosθ + 1) = 2( γ BLW γ LLW + γ B+γ L− + γ B−γ L+ )

(2)

To determine surface energy components of bitumen surface, γ BLW , γ B+ and γ B− , three different probing liquids of known surface energy components should be used for contact angle measurements, and then three equations can be obtained based on eq 2. Thus, the energy components of bitumen surface can be obtained by eq 3, where L1, L2 and L3 represent three different probe liquids.

 γ LW  B  γ B+   γ−  B

   γ LW    L1  = 2  γ LLW2       γ LW    L3

γ

− L1

γ

− L2

γ L−3

γ   γ L+2   γ L+3   + L1

−1

 γ L1 (cosθ1 + 1)   γ (cosθ + 1)  2  L2   γ (cosθ + 1)  3  L3 

(3)

2.4. Surface Characterization. Tapping mode imaging was applied to characterize the topography and roughness of bitumen surface both in air and aqueous solution with MFP-3D AFM system (Asylum Research, Santa Barbara, CA, USA). The silicon probe (RTESP-150, Bruker) and silicon nitride probe (SCANASYST-FLUID, Bruker) were respectively used for tapping mode imaging in air and liquid. A Zetasizer Nano (Malvern Instruments Ltd., United Kingdom) was used to measure the zeta potential of bitumen drop in 1 mM NaCl at different pHs.

ACS Paragon Plus Environment

7

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 8 of 33

2.5. Force Measurements. The interaction between air bubble and bitumen surface was measured using bubble probe AFM technique. Air bubbles were generated and immobilized on the glass disk of a fluid cell by purging air through an ultra-sharp glass pipette into aqueous solution. A custom-made rectangular silicon cantilever (400 × 70 × 2 µm) with a circular gold patch was used to pick up an air bubble of suitable size (typically 60-90 µm radius) to create a bubble probe. Prior to creating the bubble probe, the glass disk and AFM cantilever were hydrophobized following an established method.8, 11, 32 The spring constant of the cantilever was determined to be 0.3-0.4 N/m using the Hutter and Bechhoefer method.46 The cantileveranchored bubble was positioned over bitumen surface and then driven to approach/retract toward the surface until a fixed deflection of the cantilever was reached or bubble attachment occurred. The force measurements were conducted at a fixed driving velocity of 1 µm/s to minimize the hydrodynamic effect unless the effect of hydrodynamic condition on the bubble-bitumen interaction was considered. The movement of the cantilever with anchored air bubble and the corresponding interaction forces were recorded as a function of time by AFM software. After each force measurement, the surface tension of water was measured. The measured surface tension remained unchanged before and after the experiments, indicating the possible interfaceactive components (e.g., asphaltenes) in bitumen had negligible influence on the force measurements in this work. A schematic of typical experiment setup using the bubble probe AFM is illustrated in Figure 1B. It is noted that in common practice of froth flotation, chemical reagents such as surfactants (as frother) can be adsorbed at air/water interface, changing the air/water interfacial properties (e.g., interfacial tension) and thereby affecting the bubble-bubble and bubble-substrate interactions. In this work, we only considered a model system without the addition of chemical reagent.

ACS Paragon Plus Environment

8

Page 9 of 33

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

2.6. Theoretical Model. The measured force results between air bubble and bitumen surface was analyzed by a theoretical model based on Reynolds lubrication theory coupled with augmented Young-Laplace equation.

The drainage process of confined thin water film between air bubble and bitumen surface is described by the Reynolds lubrication theory.8, 11, 32-41

∂h 1 ∂  3 ∂p  =  rh  ∂t 12µ r ∂r  ∂r 

(4)

where h(r, t) is the thin film thickness, µ is the dynamic viscosity of water, and p(r, t) is the excess hydrodynamic pressure in the film relative to the bulk solution. Consistent with recent reports, immobile boundary condition was assumed at air/water and bitumen/water interfaces.8, 11, 14, 32-36, 42-43, 47

The bubble deformation under the combined effect of hydrodynamic pressure and disjoining pressure is described by the augmented Young-Laplace equation.8, 11, 32-41

γ ∂  ∂h 

2γ − p−∏ r  = r ∂r  ∂r  R

(5)

where γ is the air/water interfacial tension, R is the bubble radius, and Π(r, t) is the overall disjoining pressure contributed from surface forces such as vdW, EDL, and hydrophobic interactions. The disjoining pressure attributed to vdW, EDL and hydrophobic interactions ΠvdW, ΠEDL and ΠHB can be given by eqs 6, 7 and 8, respectively, where AA−W −B is the Hamaker constant for air-water-bitumen, κ −1 is the Debye length, ψ A and ψ B are the surface potential of air bubble

ACS Paragon Plus Environment

9

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 10 of 33

and bitumen surface, respectively, D0 is the decay length of hydrophobic interaction and C is a constant (N/m). 8, 11, 32, 48

∏vdW = −

∏ EDL =

AA−W − B 6π h3

(6)

2ε 0εκ 2 ( e+κ h + e−κ h )ψ Aψ B − (ψ A2 +ψ B2 ) 

∏ HB = −

(e

+κ h

− e −κ h )

2

C γ (1 − cos θ ) − h / D0 e − h / D0 = − e 2π D0 D0

(7)

(8)

The overall interaction between an air bubble and a bitumen surface F(t) is theoretically calculated by integrating p(r, t) and Π(r, t) based on Derjaguin approximation (eq 9).8, 11, 32-41 ∞

F (t ) = 2π ∫ ( p ( r , t ) + ∏ ( h ( r , t )))rdr

(9)

0

By fitting the measured force results with the theoretical model of Reynolds lubrication theory coupled with augmented Young-Laplace equation, parameters such as surface potential ψ and decay length D0 can be determined, and the evolution of confined thin film thickness h(r, t) with time can also be calculated.

3. RESULTS AND DISCUSSION 3.1. Surface Energy. The contact angles of three probing liquids (i.e., diiodomethane, water and ethylene glycol) on bitumen surface were determined using a sessile drop method. Based on the measured contact angle values and known surface energy components of these three probing liquids, the surface energy components of bitumen were determined. As shown in Table 1, the nonpolar liquid (e.g., diiodomethane) showed a lower contact angle on bitumen surface than the polar liquid (e.g., water), indicating the natural hydrophobicity of bitumen

ACS Paragon Plus Environment

10

Page 11 of 33

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

surface. The surface energy of bitumen was determined to be γ = 40.32 mJ/m2, comprising the Lifshitz-vdW component γ LW of 40.31 mJ/m2 and Lewis acid-base component γ AB of 0.01 mJ/m2. The extremely low γ AB value reveals that bitumen surface is almost non-polar. Table 1. Measured contact angle values (°) of probing liquids (diiodomethane, water and ethylene glycerol) on bitumen surface, literature reported surface energy components (mJ/m2) of probing liquids44-45 and calculated surface energy components (mJ/m2) of bitumen

γ+

γ LW

γ−

γ

γ AB

Measured contact angle on bitumen surface

Diiodomethane

50.80

0.00

0.00

0.00

50.80

38.6 ± 1.3

Water

21.80

25.50

25.50

51.00

72.80

93.2 ± 1.5

Ethylene glycerol

29.00

1.92

47.00

19.00

48.00

61.5 ± 0.9

Bitumen

40.31

0.00

0.90

0.01

40.32

(calculated)

The Lifshitz-vdW and Lewis acid-base interactions could be considered based on the free energy of interaction between air and bitumen in water. The free energy of Lifshitz-vdW

∆GALW−W −B and Lewis acid-base ∆GAAB−W −B interactions for air-water-bitumen is respectively calculated to be 15.69 and -92.36 mJ/m2 based on eqs 10 and 11.44-45 ∆GALW−W − B = 2 γ WLW ( γ BLW − γ WLW )

(10)

∆G AAB−W − B = 2[ γ W+ ( γ B− − γ W− ) + γ W− ( γ B+ − γ W+ )]

(11)

where the subscript A, B or W represents air, bitumen or water, respectively. The positive

∆GALW−W −B suggests that the repulsive vdW interaction inhibits the attachment between air and AB bitumen in water. The negative ∆GA−W −B reveals that the attractive interaction, including

hydrophobic interaction, facilitates the attachment between air and bitumen in water.

ACS Paragon Plus Environment

11

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 12 of 33

The Hamaker constant for two identical materials interacting in vacuum could be 2 estimated using γ ≈ A / 24π l0 , where l0 = 0.165 nm is the typical cut-off separation.48 Thus, the

Hamaker constant of bitumen ABB is calculated to be 8.3 × 10-20 J. Based on the Hamaker constant of water AWW = 3.7 × 10-20 J, the Hamaker constant for air and bitumen interacting across water AA-W-B is estimated to be -1.8 × 10-20 J.48 From thermodynamic consideration, the interaction energy per unit area of two planar half-spaces in water should approach their work of adhesion with separation h decreasing to 0. Therefore, the constant C in eq 8 is determined to be 0.41 N/m. 3.2. Surface Morphology. The AFM topographic image of bitumen surface in air is shown in Figure 2. The root-mean-square (rms) roughness is about 0.45 nm, which indicates that a very smooth bitumen surface can be obtained by the spin coating process.

Figure 2. AFM topographic image (5 × 5 µm2) of spin-coated bitumen surface in air

The morphology of bitumen surface in aqueous solutions with various ionic strengths and pH conditions was also investigated. As illustrated in Figure 3, the bitumen surface became noticeably rougher than that in air and some distinct domains were observed on bitumen surface, which coincided with the protrusion structure observed on bitumen droplet in water from freeze-

ACS Paragon Plus Environment

12

Page 13 of 33

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

fracture scanning electron microscope image.18, 49 Interestingly, the size of aggregation domains formed in 1 mM NaCl (Figures 3A, 3B and 3C) was found to be evidently smaller than that in 500 mM NaCl (Figures 3D, 3E and 3F). It is worth noting that all the bitumen substrates were imaged after being immersed in aqueous solutions for at least 15 min, after which the morphologies remained almost unchanged even after several hours. It is noted that the formation of nanobubbles on bitumen surface can be ruled out in this study based on AFM imaging due to the continuous aggregation patterns on bitumen surfaces, similar module of the aggregation domains and surrounding areas, and the increased size of aggregation domains with increasing pH and salt concentration, which are very different from the observed nanobubbles on hydrophobic substrates reported previously.35

Figure 3. The morphology (5 × 5 µm2) of bitumen surface in 1 mM NaCl at pH 4.0 (A), pH 5.8 (B) and pH 8.5 (C) and 500 mM NaCl at pH 4.0 (D), pH 5.8 (E) and pH 8.5 (F).

In 1 mM NaCl, the rms roughness of bitumen surface increased from 0.53 nm to 0.57 and 0.59 nm with solution pH increasing from 4.0 to 5.8 and 8.5. It is widely recognized that bitumen

ACS Paragon Plus Environment

13

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 14 of 33

contains large amounts of natural interface-active components (e.g., asphaltenes) with functional groups (e.g., carboxyl and phenolic groups) that can be deprotonated with the addition of alkaline.20,

50

Therefore, the deprotonation of polar functional groups and conformational

rearrangement of the molecules at bitumen/water interface at higher pH could lead to larger roughness. As solution pH increased from 4.0 to 5.8 and 8.5 in 500 mM NaCl, the rms roughness of bitumen surface raised from 1.28 nm to 1.50 and 2.07 nm. The enhanced surface roughness at high salt concentration is mainly attributed to the decreased electrostatic repulsion among the polar functional groups on bitumen molecules and/or aggregates. 3.3. Bubble-bitumen Interaction in 1 mM NaCl. Figure 4 shows the interaction force profiles measured between air bubble and bitumen surface in 1 mM NaCl at different pH conditions at v = 1 µm/s. It is worth noting that no bubble attachment could be observed during the approach-retraction cycle under all the pH conditions tested. With the cantilever-anchored bubble approaching bitumen surface, strong repulsion was registered that prevented the bubble from attaching to bitumen surface. During retraction of the cantilever, the measured repulsion gradually decreased until an attractive maximum was achieved due to the hydrodynamic suction effect.51-52 For the air-water-bitumen system, the vdW interaction is repulsive at any separation, and the EDL interaction between negatively charged air bubble and bitumen surface is also repulsive. In 1 mM NaCl, the Debye length is calculated to be 9.6 nm, and the relatively longrange EDL repulsion could inhibit the bubble attachment on bitumen surface. The surface potentials of air bubble in 1 mM NaCl at different pH conditions were obtained by conducting the bubble-bubble interaction measurements (Figure S1 in supporting information) and compared with the literature values in Table 2. The aforementioned theoretical model incorporating surface forces was applied to analyze the measured forces (open symbols) and the fitted results were

ACS Paragon Plus Environment

14

Page 15 of 33

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

shown as the red curve in Figures 4A, 4B and 4C. The fitted surface potentials of bitumen surface in 1 mM NaCl at different pH conditions are summarized in Figure 5, which are consistent with the measured zeta potentials and literature values.20, 50 As solution pH increased from 4.0 to 5.8 and 8.5, surface potential of bitumen changed from -35 ± 5 mV to -64 ± 8 and 80 ± 10 mV, thereby strengthening the EDL repulsion between air bubble and bitumen surface. The change of surface potential of bitumen is most likely due to the increased adsorption of hydroxide ion and the deprotonation of functional groups (e.g., carboxyl and phenolic groups) of interface-active components at bitumen/water interface with the addition of alkaline.

Figure 4. Interaction forces (A-C), calculated thin film profiles at maximum force load (D-F) and calculated disjoining pressure profiles (G-I) between air bubble and bitumen surface in 1 mM NaCl at pH 4.0 (A, D and G), pH 5.8 (B, E and H) and pH 8.5 (C, F and I) at v = 1 µm/s (open symbols for experiment results and solid curves for theoretical calculations). The

ACS Paragon Plus Environment

15

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 16 of 33

theoretical calculations were based on the interactions of a bubble vs. a flat substrate. The bubble radius R is 82 µm at pH 4.0, 81 µm at pH 5.6, and 81 µm at pH 8.5.

The calculated profiles of thin water film confined between air bubble and bitumen surface at maximum force load are illustrated in Figures 4D, 4E and 4F. The central portion of bubble surfaces was flattened due to the EDL disjoining pressure that balanced the Laplace pressure inside the bubbles. With solution pH increasing from 4.0 to 5.8 and 8.5, the central separation between air bubble and bitumen surface dramatically increased from 17.4 nm to 28.7 and 32.2 nm, suggesting that the enhanced EDL repulsion at higher solution pH could effectively impede thin water film drainage and prevent air bubble from getting close to bitumen surface. It is worth noting that surface roughness plays an important role in surface force measurements in fluids.53, 54 The rms roughness of bitumen surfaces under different pH conditions in Figure 4 was almost the same, thus the effect of surface roughness on the comparison of bubble-bitumen surface interactions under various solution conditions can be neglected here. It should be also noted that all the theoretical calculations were based on the interactions of a bubble vs. a flat substrate. The calculated disjoining pressure profiles between air bubble and bitumen surface are shown in Figures 4G, 4H and 4I. The EDL repulsion is the dominant interaction and the hydrophobic attraction is too weak to trigger the bubble-bitumen attachment.

Table 2. Comparison of theoretically fitted surface potential values of air bubble and the literature values in 1 mM NaCl55-57 pH

Fitted value (mV)

Literature value (mV)

4.0

-20 ± 5

-15

ACS Paragon Plus Environment

16

Page 17 of 33

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

-10 5.8

-32 ± 6

-35 -37

8.5

-45 ± 8

-50 -65

Figure 5. Comparison of theoretically fitted surface potential and measured zeta potential values of bitumen in 1 mM NaCl with literature reported zeta potential values in 1 mM KCl.20, 50

3.4. Bubble-bitumen Interaction in 500 mM NaCl. Figure 6 shows the interaction forces between air bubble and bitumen surface in 500 mM NaCl at various pH conditions at v = 1 µm/s. As shown in Figure 6A, a sudden “jump-in” behavior was observed during approach when the measured force reached ~6.5 nN, indicating that the bubble was attached on bitumen surface, which was observed from the optical microscope. In 500 mM NaCl, the EDL interaction is significantly screened and vdW interaction is repulsive, so the observed bubble attachment was induced mainly by the attractive hydrophobic interaction which was incorporated into the

ACS Paragon Plus Environment

17

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 18 of 33

aforementioned theoretical model. By fitting the measured forces (open symbols) with the theoretical model (red curve), the decay length of hydrophobic interaction was determined to be D0 = 0.95 ± 0.05 nm. The calculated profile of thin water film confined between air bubble and bitumen surface in Figure 6D shows that the critical central separation just before bubble attachment was calculated to be 9.2 nm. In Figures 6B and 6C, the “jump-in” behavior was also observed at the measured force of ~7.5 nN at pH 5.8 and ~8.2 nN at pH 8.5. The theoretical calculation shows that the decay length of hydrophobic interaction was predicted to be D0 = 0.78 ± 0.05 nm at pH 5.8 (Figure 6B) and D0 = 0.65 ± 0.07 nm at pH 8.5 (Figure 6C), and the critical central separation just before bubble attachment was calculated to be 8.3 nm at pH 5.8 (Figure 6E) and 6.9 nm at pH 8.5 (Figure 6F), both of which indicates the weakened hydrophobic attraction with the addition of alkaline. The hydrophobic interaction is believed to be closely correlated to the entropic effect originating from the disruption of the hydrogen bonding network between water molecules adjacent to the hydrophobic surface.48, 58-59 The bitumen wettability was evaluated by captive bubble method, where the solid substrate was immersed in aqueous solution with coated bitumen facing downward and a bubble was injected beneath the bitumen surface. In the captive bubble method, the bitumen surface properties would be the same as that in the AFM force measurements under the respective solution conditions. By analyzing the bubble profile, the water contact angle at pH 4 was measured to be 77.1° ± 2.0°, which decreased to 71.7° ± 1.5° at pH 5.8 and 66.8° ± 1.7° at pH 8.5. The decreased surface hydrophobicity is most likely due to the deprotonation of functional groups (e.g., carboxyl and phenolic groups) of interface-active molecules in alkaline environment that enables higher mobility of water molecules around bitumen surface. The change of bitumen wettability may alter the mobility of water molecules and water correlations at the vicinity of bitumen/water interface, and in response

ACS Paragon Plus Environment

18

Page 19 of 33

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

alter the magnitude and range of hydrophobic interaction. The calculated disjoining pressure profiles between air bubble and bitumen surface in Figures 6G, 6H and 6I shows the hydrophobic attraction is the driving interaction for the bubble attachment on the bitumen surface. It should be noted that in this work the water contact angles on bitumen surfaces measured using the captive bubble method in aqueous media and the sessile drop method in air show good agreement with the values reported previously.20, 50, 60 The difference of the water contact angles measured using these two methods was because that the bitumen surfaces were immersed in aqueous solution for the captive bubble method and the functional groups of natural interface-active components could be deprotonated and exposed at bitumen/water interface, thereby changing the surface morphology and wettability.

ACS Paragon Plus Environment

19

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 33

Figure 6. Interaction forces (A-C), calculated thin film profiles at maximum force load (D-F) and calculated disjoining pressure profiles (G-I) between air bubble and bitumen surface in 500 mM NaCl at pH 4.0 (A, D and G), pH 5.8 (B, E and H) and pH 8.5 (C, F and I) at v = 1 µm/s (open symbols for experiment results and solid curves for theoretical calculations). The theoretical calculations were based on the interactions of a bubble vs. a flat substrate. The bubble radius R is 73 µm at pH 4.0, 74 µm at pH 5.6, and 77 µm at pH 8.5.

3.4. Effect of Calcium Ion. In industrial processes such as oil production from the Athabasca oil sands, the process water or recycled water of high salinity generally contains a certain amount of divalent cations such as Ca2+. It is important to understand the influence of calcium ions on the surface properties (e.g., morphology) of bitumen and bubble-bitumen interaction. The morphologies of bitumen surface in 500 mM NaCl at pH 8.5 with the addition of 1 mM and 10 mM CaCl2 are shown in Figures 7A and 7B, respectively. In contrast to the morphology without CaCl2 addition in Figure 3F, the aggregation domains observed on bitumen surface were much larger, and the rms roughness increased from 2.07 nm (Figure 3F) to 3.41 and 4.28 nm after the addition of 1 mM and 10 mM CaCl2, respectively (Figures 7A and 7B). The larger aggregation domains formed on bitumen surface in the presence of calcium ion is most likely through the strong interaction between calcium ion and the functional groups of interfaceactive molecules.

ACS Paragon Plus Environment

20

Page 21 of 33

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. The morphologies (5 × 5 µm2) of bitumen surface in 500 mM NaCl at pH 8.5 in presence of 1 mM (A) and 10 mM (B) CaCl2.

Figures 8A and 8B show the interaction between air bubble and bitumen surface in 500 mM NaCl at pH 8.5 with the addition of 1 mM and 10 mM CaCl2, respectively. The “jump-in” behavior was observed during approach when the measured force reached ~5.8 nN and ~4.7 nN. The theoretical calculation shows that the decay length of hydrophobic interaction was fitted to be D0 = 1.07 ± 0.07 nm and D0 = 1.20 ± 0.10 nm with the addition of 1 mM and 10 mM CaCl2, respectively, and the corresponding critical central separation just before bubble attachment was calculated to be 11.1 nm and 12.2 nm as illustrated in Figures 8C and 8D. In contrast to the interaction without CaCl2 addition (Figure 6C), the critical force load required for bubblebitumen attachment in the presence of calcium ion is much weaker, and the decay length of hydrophobic interaction and calculated central separation are higher. It is likely that calcium ion may bridge the carboxyl groups of interface-active components, which increases the surface roughness (Figure 7) and leaves the hydrophobic moieties exposed to the aqueous solution and thereby strengthens the apparent hydrophobic interaction between air bubble and bitumen surface. The calculated disjoining pressure profiles between air bubble and bitumen surface in Figures 8E and 8F shows the dominant role of hydrophobic attraction.

ACS Paragon Plus Environment

21

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 33

Figure 8. Interaction forces (A and B), calculated thin film profiles at maximum force load (C and D) and calculated disjoining pressure profiles (E and F) between air bubble and bitumen surface in 500 mM NaCl at pH 8.5 with the addition of 1 mM (A, C and E) and 10 mM (B, D and F) CaCl2 at v = 1 µm/s (open symbols for experiment results and solid curves for theoretical calculations). The theoretical calculations were based on the interactions of a bubble vs. a flat substrate. The bubble radius R is 67 µm for 1 mM CaCl2 and 73 µm for 10 mM CaCl2. 3.5. Effect of hydrodynamic condition. The measured forces (open symbols) between air bubble and bitumen surface in 500 mM NaCl at pH 8.5 under higher hydrodynamic conditions (i.e., v = 5 and 20 µm/s) are shown in Figure 9. Compared to the bubble attachment observed during approach at v = 1 µm/s (Figure 6C), the bubble-bitumen attachment was observed during retraction at v = 5 µm/s (Figure 9A) and no bubble attachment was observed during the approach-retraction cycle at v = 20 µm/s (Figure 9B) under the same solution condition, indicating that the bubble attachment could be inhibited at higher approach velocity due to increased hydrodynamic force. The fitted decay length D0 = 0.65 ± 0.07 nm at v = 1 µm/s was used to predict the interaction force and bubble attachment behavior at v = 5 and 20 µm/s, and the theoretically

ACS Paragon Plus Environment

22

Page 23 of 33

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

calculated force results are shown as solid curves in Figure 9. It is obvious that theoretical calculation agrees very well with the measured force results, indicating that D0 = 0.65 ± 0.07 nm at pH 8.5 reflects the true hydrophobic interaction that is independent of hydrodynamic conditions. From the above results, it can be seen that the hydrodynamic force and surface forces (e.g., EDL and hydrophobic interactions) play critical roles in bubble-bitumen interaction and attachment. It is worth noting that the hydrophobic interaction was much less important than the hydrodynamic interaction at a higher velocity (see Figure S2 in Supporting Information).

Figure 9. Interaction forces between air bubble and bitumen surface in 500 mM NaCl at pH 8.5 at v = 5 µm/s (bubble radius R=75 µm) (A) and at v = 20 µm/s (bubble radius R=73 µm) (B) (open symbols for experiment results and solid curves for theoretical calculations).

CONCLUSIONS The colloidal interaction involving deformable air bubbles has attracted tremendous interest in a wide range of engineering applications. A bubble probe AFM technique was employed, for the first time, to directly measure the interaction forces between air bubbles and bitumen surfaces in complex aqueous media with varying solution salinity, pH and divalent electrolytes addition. The bitumen surface shows a rougher morphology with some distinct aggregation domains at higher ionic strength and/or higher solution pH. The force results reveal the critical role of ionic

ACS Paragon Plus Environment

23

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 24 of 33

strength and solution pH in bubble-bitumen interaction and attachment, which could be analyzed by a theoretical model based on Reynolds lubrication theory and augmented Young-Laplace equation by including the effect of disjoining pressure. As solution pH increased from 4.0 to 8.5, the strengthened EDL repulsion could more readily inhibit bubble attachment in 1 mM NaCl. In 500 mM NaCl, the EDL repulsion was significantly suppressed and the hydrophobic attraction between bubble and bitumen played a critical role in the bubble-bitumen attachment. While it is also noted that the hydrophobic attraction between bubble and bitumen became weaker with increasing pH from 4.0 to 8.5 in 500 mM NaCl, which was mainly due to enhanced surface hydrophilicity of bitumen under higher pH condition. It was found that the presence of calcium ions in 500 mM NaCl could strengthen the hydrophobic attraction between air bubble and bitumen surface and facilitate the bubble-bitumen attachment, which was most likely due to the bridging interaction between calcium ions and the carboxyl groups of interface-active molecules on bitumen surfaces thereby leading to higher surface roughness and hydrophobic moieties/aggregates on bitumen as confirmed by AFM imaging. This work provides valuable quantitative information on the interaction mechanism between air bubbles and bitumen surfaces at the nanoscale, with important implications in heavy oil production. The methodology can be readily extended to many other interfacial interactions in industrial processes such as oil-water separation and wastewater treatment.

SUPPORTING INFORMATION Force profiles for bubble-bubble Interaction in 1 mM NaCl, and calculated profiles of disjoining pressure and hydrodynamic pressure under different interaction velocities.

ACS Paragon Plus Environment

24

Page 25 of 33

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

ACKNOWLEDGEMENTS We gratefully acknowledge the financial support from the Natural Sciences and Engineering Research Council of Canada (NSERC), the Future Energy Systems under the Canada First Research Excellence Fund, the Canada Foundation for Innovation (CFI), the Alberta Advanced Education & Technology Small Equipment Grants Program (AET/SEGP) and the Canada Research Chairs Program (H. Zeng).

ACS Paragon Plus Environment

25

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 26 of 33

REFERENCES 1. Lindner, J. R. Microbubbles in Medical Imaging: Current Applications and Future Directions. Nat. Rev. Drug Discov. 2004, 3, 527-533. 2. Prentice, P.; Cuschieri, A.; Dholakia, K.; Prausnitz, M.; Campbell, P. Membrane Disruption by Optically Controlled Microbubble Cavitation. Nature Physics 2005, 1, 107-110. 3. Yec, C. C.; Zeng, H. C. Nanobubbles within a Microbubble: Synthesis and Self-Assembly of Hollow Manganese Silicate and Its Metal-Doped Derivatives. ACS nano 2014, 8, 6407-6416. 4. Rao, S. R. Surface Chemistry of Froth Flotation: Volume 1: Fundamentals; Springer Science & Business Media, 2013. 5. Ferrara, K.; Pollard, R.; Borden, M. Ultrasound Microbubble Contrast Agents: Fundamentals and Application to Gene and Drug Delivery. Annu. Rev. Biomed. Eng. 2007, 9, 415-447. 6. Prakash, M.; Gershenfeld, N. Microfluidic Bubble Logic. Science 2007, 315, 832-835. 7. Dickinson, E.; Ettelaie, R.; Kostakis, T.; Murray, B. S. Factors Controlling the Formation and Stability of Air Bubbles Stabilized by Partially Hydrophobic Silica Nanoparticles. Langmuir 2004, 20, 8517-8525. 8. Shi, C.; Chan, D. Y.; Liu, Q.; Zeng, H. Probing the Hydrophobic Interaction between Air Bubbles and Partially Hydrophobic Surfaces Using Atomic Force Microscopy. J. Phys. Chem. C 2014, 118, 25000-25008. 9. Safouane, M.; Langevin, D.; Binks, B. Effect of Particle Hydrophobicity on the Properties of Silica Particle Layers at the Air-Water Interface. Langmuir 2007, 23, 11546-11553. 10. Fuerstenau, M. C.; Jameson, G. J.; Yoon, R.-H. Froth Flotation: A Century of Innovation; SME, 2007.

ACS Paragon Plus Environment

26

Page 27 of 33

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

11. Xie, L.; Shi, C.; Wang, J.; Huang, J.; Lu, Q.; Liu, Q.; Zeng, H. Probing the Interaction between Air Bubble and Sphalerite Mineral Surface Using Atomic Force Microscope. Langmuir 2015, 31, 2438-2446. 12. Ralston, J.; Fornasiero, D.; Hayes, R. Bubble–Particle Attachment and Detachment in Flotation. Int. J. Miner. Process. 1999, 56, 133-164. 13. Horn, R. G.; Asadullah, M.; Connor, J. N. Thin Film Drainage: Hydrodynamic and Disjoining Pressures Determined from Experimental Measurements of the Shape of a Fluid Drop Approaching a Solid Wall. Langmuir 2006, 22, 2610-2619. 14. Tabor, R. F.; Chan, D. Y.; Grieser, F.; Dagastine, R. R. Structural Forces in Soft Matter Systems. J. Phys. Chem. Lett. 2011, 2, 434-437. 15. Masliyah, J.; Zhou, Z. J.; Xu, Z.; Czarnecki, J.; Hamza, H. Understanding Water‐Based Bitumen Extraction from Athabasca Oil Sands. Can. J. Chem. Eng. 2004, 82, 628-654. 16. Gafonova, O. V.; Yarranton, H. W. The Stabilization of Water-in-Hydrocarbon Emulsions by Asphaltenes and Resins. J. Colloid Interface Sci. 2001, 241, 469-478. 17. Yarranton, H. W.; Hussein, H.; Masliyah, J. H. Water-in-Hydrocarbon Emulsions Stabilized by Asphaltenes at Low Concentrations. J. Colloid Interface Sci. 2000, 228, 52-63. 18. Wu, X.; Czarnecki, J.; Hamza, N.; Masliyah, J. Interaction Forces between Bitumen Droplets in Water. Langmuir 1999, 15, 5244-5250. 19. Legrand, J.; Chamerois, M.; Placin, F.; Poirier, J.; Bibette, J.; Leal-Calderon, F. Solid Colloidal Particles Inducing Coalescence in Bitumen-in-Water Emulsions. Langmuir 2005, 21, 64-70. 20. Liu, J.; Xu, Z.; Masliyah, J. Studies on Bitumen-Silica Interaction in Aqueous Solutions by Atomic Force Microscopy. Langmuir 2003, 19, 3911-3920.

ACS Paragon Plus Environment

27

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 28 of 33

21. Long, J.; Zhang, L.; Xu, Z.; Masliyah, J. H. Colloidal Interactions between LangmuirBlodgett Bitumen Films and Fine Solid Particles. Langmuir 2006, 22, 8831-8839. 22. Pashazanusi, L.; Lwoya, B.; Oak, S.; Khosla, T.; Albert, J. N.; Tian, Y.; Bansal, G.; Kumar, N.; Pesika, N. S. Enhanced Adhesion of Mosquitoes to Rough Surfaces. ACS Appl. Mater. Interfaces 2017, 9, 24373–24380. 23. Zhang, X.; Kumar, A.; Scales, P. J. Effects of Solvency and Interfacial Nanobubbles on Surface Forces and Bubble Attachment at Solid Surfaces. Langmuir 2011, 27, 2484-2491. 24. Alcantar, N.; Israelachvili, J.; Boles, J. Forces and Ionic Transport between Mica Surfaces: Implications for Pressure Solution. Geochim. Cosmochim. Acta 2003, 67, 1289-1304. 25. Tian, Y.; Pesika, N.; Zeng, H.; Rosenberg, K.; Zhao, B.; McGuiggan, P.; Autumn, K.; Israelachvili, J. Adhesion and Friction in Gecko Toe Attachment and Detachment. Proc. Natl. Acad. Sci. U. S. A. 2006, 103, 19320-19325. 26. Kristiansen, K.; Stock, P.; Baimpos, T.; Raman, S.; Harada, J. K.; Israelachvili, J. N.; Valtiner, M. Influence of Molecular Dipole Orientations on Long-Range Exponential Interaction Forces at Hydrophobic Contacts in Aqueous Solutions. ACS nano 2014, 8, 10870-10877. 27. Natarajan, A.; Xie, J.; Wang, S.; Liu, Q.; Masliyah, J.; Zeng, H.; Xu, Z. Understanding Molecular Interactions of Asphaltenes in Organic Solvents Using a Surface Force Apparatus. J. Phys. Chem. C 2011, 115, 16043-16051. 28. Zhang, L.; Xie, L.; Shi, C.; Huang, J.; Liu, Q.; Zeng, H. Mechanistic Understanding of Asphaltene Surface Interactions in Aqueous Media. Energy & Fuels 2016, 31, 3348–3357. 29. Gillies, G.; Kappl, M.; Butt, H.-J. Direct Measurements of Particle–Bubble Interactions. Adv. Colloid Interface Sci. 2005, 114, 165-172.

ACS Paragon Plus Environment

28

Page 29 of 33

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

30. Johnson, D. J.; Miles, N. J.; Hilal, N. Quantification of Particle–Bubble Interactions Using Atomic Force Microscopy: A Review. Adv. Colloid Interface Sci. 2006, 127, 67-81. 31. Butt, H.-J. A Technique for Measuring the Force between a Colloidal Particle in Water and a Bubble. J. Colloid Interface Sci. 1994, 166, 109-117. 32. Shi, C.; Cui, X.; Xie, L.; Liu, Q.; Chan, D. Y.; Israelachvili, J. N.; Zeng, H. Measuring Forces and Spatiotemporal Evolution of Thin Water Films between an Air Bubble and Solid Surfaces of Different Hydrophobicity. ACS nano 2014, 9, 95-104. 33. Vakarelski, I. U.; Manica, R.; Tang, X.; O’Shea, S. J.; Stevens, G. W.; Grieser, F.; Dagastine, R. R.; Chan, D. Y. Dynamic Interactions between Microbubbles in Water. Proc. Natl. Acad. Sci. U. S. A. 2010, 107, 11177-11182. 34. Tabor, R. F.; Manica, R.; Chan, D. Y.; Grieser, F.; Dagastine, R. R. Repulsive Van Der Waals Forces in Soft Matter: Why Bubbles Do Not Stick to Walls. Phys. Rev. Lett. 2011, 106, 064501. 35. Cui, X.; Shi, C.; Xie, L.; Liu, J.; Zeng, H. Probing Interactions between Air Bubble and Hydrophobic Polymer Surface: Impact of Solution Salinity and Interfacial Nanobubbles. Langmuir 2016, 32, 11236–11244. 36. Shi, C.; Zhang, L.; Xie, L.; Lu, X.; Liu, Q.; Mantilla, C. A.; van den Berg, F. G.; Zeng, H. Interaction Mechanism of Oil-in-Water Emulsions with Asphaltenes Determined Using Droplet Probe Afm. Langmuir 2016, 32, 2302-2310. 37. Shi, C.; Yan, B.; Xie, L.; Zhang, L.; Wang, J.; Takahara, A.; Zeng, H. Long‐Range Hydrophilic Attraction between Water and Polyelectrolyte Surfaces in Oil. Angew. Chem., Int. Ed. 2016, 55, 15017-15021.

ACS Paragon Plus Environment

29

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 30 of 33

38. Shi, C.; Zhang, L.; Xie, L.; Lu, X.; Liu, Q.; He, J.; Mantilla, C. A.; Zeng, H. Surface Interaction of Water-in-Oil Emulsion Droplets with Interfacially Active Asphaltenes. Langmuir 2017, 33, 1265-1274. 39. Xie, L.; Shi, C.; Cui, X.; Zeng, H. Surface Forces and Interaction Mechanisms of Emulsion Drops and Gas Bubbles in Complex Fluids. Langmuir 2017, 33, 3911-3925. 40. Xie, L.; Wang, J.; Yuan, D.; Shi, C.; Cui, X.; Zhang, H.; Liu, Q.; Liu, Q.; Zeng, H. Interaction Mechanisms between Air Bubble and Molybdenite Surface: Impact of Solution Salinity and Polymer Adsorption. Langmuir 2017, 33, 2353-2361. 41. Cui, X.; Shi, C.; Zhang, S.; Xie, L.; Liu, J.; Jiang, D.; Zeng, H. Probing the Effect of Salinity and Ph on Surface Interactions between Air Bubbles and Hydrophobic Solids: Implications for Colloidal Assembly at Air/Water Interfaces. Chem. Asian J. 2017, 12, 1568-1577. 42. Tabor, R. F.; Lockie, H.; Mair, D.; Manica, R.; Chan, D. Y.; Grieser, F.; Dagastine, R. R. Combined Afm−Confocal Microscopy of Oil Droplets: Absolute Separations and Forces in Nanofilms. J. Phys. Chem. Lett. 2011, 2, 961-965. 43. Tabor, R. F.; Wu, C.; Grieser, F.; Dagastine, R. R.; Chan, D. Y. Measurement of the Hydrophobic Force in a Soft Matter System. J. Phys. Chem. Lett. 2013, 4, 3872-3877. 44. Van Oss, C. J. Interfacial Forces in Aqueous Media; CRC press, 2006. 45. Van Oss, C. Acid—Base Interfacial Interactions in Aqueous Media. Colloids Surf., A 1993, 78, 1-49. 46. Hutter, J. L.; Bechhoefer, J. Calibration of Atomic‐Force Microscope Tips. Rev. Sci. Instrum. 1993, 64, 1868-1873. 47. Tabor, R. F.; Chan, D. Y.; Grieser, F.; Dagastine, R. R. Anomalous Stability of Carbon Dioxide in Ph‐Controlled Bubble Coalescence. Angew. Chem., Int. Ed. 2011, 123, 3516-3518.

ACS Paragon Plus Environment

30

Page 31 of 33

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

48. Israelachvili, J. N. Intermolecular and Surface Forces: Revised Third Edition; Academic press, 2011. 49. Wu, X.; Dabros, T.; Czarnecki, J. Determining the Colloidal Forces between Bitumen Droplets in Water Using the Hydrodynamic Force Balance Technique. Langmuir 1999, 15, 8706-8713. 50. Liu, J.; Xu, Z.; Masliyah, J. Colloidal Forces between Bitumen Surfaces in Aqueous Solutions Measured with Atomic Force Microscope. Colloids Surf., A 2005, 260, 217-228. 51. Christenson, H. K.; Claesson, P. M. Cavitation and the Interaction between Macroscopic Hydrophobic Surfaces. Science 1988, 239, 390-392. 52. Tsao, Y.-H.; Evans, D. F.; Wennerstrom, H. Long-Range Attractive Force between Hydrophobic Surfaces Observed by Atomic Force Microscopy. Science 1993, 262, 547-550. 53. Valtiner, M.; Kristiansen, K.; Greene, G. W.; Israelachvili, J. N. Effect of Surface Roughness and Electrostatic Surface Potentials on Forces between Dissimilar Surfaces in Aqueous Solution. Adv. Mater. 2011, 23, 2294-2299. 54. Zeng, H. Polymer Adhesion, Friction, and Lubrication; John Wiley & Sons, 2013. 55. Yang, C.; Dabros, T.; Li, D.; Czarnecki, J.; Masliyah, J. H. Measurement of the Zeta Potential of Gas Bubbles in Aqueous Solutions by Microelectrophoresis Method. J. Colloid Interface Sci. 2001, 243, 128-135. 56. Cho, S.-H.; Kim, J.-Y.; Chun, J.-H.; Kim, J.-D. Ultrasonic Formation of Nanobubbles and Their Zeta-Potentials in Aqueous Electrolyte and Surfactant Solutions. Colloids Surf., A 2005, 269, 28-34. 57. Han, M.; Kim, M.; Shin, M. Generation of a Positively Charged Bubble and Its Possible Mechanism of Formation. J. Water Supply Res. Technol. AQUA 2006, 55, 471-478.

ACS Paragon Plus Environment

31

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 32 of 33

58. Israelachvili, J.; Pashley, R., The Hydrophobic Interaction Is Long Range, Decaying Exponentially with Distance. Nature 1982, 300, 341-342. 59. Meyer, E. E.; Rosenberg, K. J.; Israelachvili, J. Recent Progress in Understanding Hydrophobic Interactions. Proc. Natl. Acad. Sci. U. S. A. 2006, 103, 15739-15746. 60. Vargha-Butler, E. I.; Potoczny, Z. M.; Zubovits, T. K.; Budziak, C. J.; Neumann, A. W. Surface Tension of Bitumen from Contact Angle Measurements on Films of Bitumen. Energy & Fuels 1988, 2, 653-656.

ACS Paragon Plus Environment

32

Page 33 of 33

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

TOC GRAPHICS

ACS Paragon Plus Environment

33