Complex Formations between Surfactants and Polyelectrolytes of the

Jul 7, 2017 - Amirhossein Mafi†‡, Dan Hu†, and Keng C. Chou†. † Department of Chemistry, University of British Columbia, Vancouver, British ...
0 downloads 0 Views 988KB Size
Subscriber access provided by UNIVERSITY OF CONNECTICUT

Article

Complex Formations Between Surfactants and Polyelectrolytes of the Same Charge On Water Surface Amirhossein Mafi, Dan Hu, and Keng C. Chou Langmuir, Just Accepted Manuscript • DOI: 10.1021/acs.langmuir.7b01246 • Publication Date (Web): 07 Jul 2017 Downloaded from http://pubs.acs.org on July 7, 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 27

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

Complex Formations Between Surfactants and Polyelectrolytes of the Same Charge On Water Surface Amirhossein Mafi1,2, Dan Hu1, and Keng C. Chou1*

1

Department of Chemistry, University of British Columbia, Vancouver, BC V6T 1Z1, Canada

2

Department of Chemical and Biological Engineering, University of British Columbia,

Vancouver, BC V6T 1Z3, Canada

* EMAIL ADDRESS: [email protected]

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 27

Abstract

The mechanism of complex formation between surfactants and polyelectrolytes of the same charge on water surface was investigated using molecular dynamics simulations and phasesensitive sum-frequency generation vibrational spectroscopy. While the complex formation between highly-charged surfactants and polyelectrolytes of the same charge is generally expected to be prohibited by the electrostatic repulsive force, our study shows it is possible to form thermodynamically stable complexes when excess ions are present in the solution. We found that anionic partially hydrolyzed polyacrylamide (HPAM) could interact with anionic sodium dodecyl sulfate (SDS) on water surface in the presence of salts. With excess Na+ ions in the solution, the charge screening effect allows HPAM to weakly interact with SDS via hydrogen bonds. In the presence of divalent Ca2+ ions, the surfactant and the polymer is strongly coupled by forming Ca2+ ion bridges and hydrogen bonds. Our calculation shows that the presence of Ca2+ ions creates a steep binding energy of ~30 KJ/mol near the water surface. These results were qualitatively verified using phase-sensitive sum-frequency generation vibrational spectroscopy.

ACS Paragon Plus Environment

2

Page 3 of 27

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

Introduction The interaction between surfactants and polyelectrolytes has been extensively studied because of their broad applications in pharmaceuticals, personal care products, and many industrial processes.1-7 The interaction leads to the formation of surfactant-polymer complexes, which modify the interfacial properties. Investigations of surfactant-polyelectrolyte interactions are motivated in part by seeking new types of surfactant-polymer complexes with novel interfacial or dispersion properties. In the past several decades, the interactions between oppositely-charged surfactants and polyelectrolytes have been thoroughly studied.8 On the other hand, few studies have been reported on the interaction between surfactants and polyelectrolytes of the same charge. Generally, surfactants and polyelectrolytes of the same charge do not form complexes because of the long-range electrostatic repulsion force. It has been reported that amphiphilic polyelectrolytes in bulk solutions may form complexes with surfactants of the same charge via hydrophobic interactions.9-14 In these cases, the hydrophobic interaction may overcome the electrostatic repulsion and distort the polymer backbone such that the polymers form complexes with surfactants of the same charge. On a water surface, the hydrophobic interaction between surfactants and polyelectrolytes of the same charge diminishes because surfactants are known to have their hydrophilic groups pointing toward the liquid water. Therefore, complex formation via the aforementioned hydrophobic interaction becomes difficult on water surface. It has been reported that cationic polyehtylenimine (PEI) and cationic alkyltrimethylammonium bromide (CTAB) form ordered mesostructured films on water surface.15 However, the films were not thermodynamically stable as they showed a loss of structure from their neutron reflectometry profiles with time. The authors concluded that the dominant interaction between the polymer and surfactant is a

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 27

neutral/cationic interaction, where the dipole on the amine groups of PEI interacts with the charged CTAB ammonium groups. Here we study the interaction between surfactants and polymers of the same charge on water surface via a different mechanism: ion-bridged ionic/ionic interactions and hydrogen bonds. Divalent cations, such as Ca2+, may lead to a charge reversal of a charged surface.16-17 It has been reported that the addition of Ca2+ results in binding of an anionic surfactant onto an anionic surface. On the other hand, monovalent ions (Na+) do not show any evidence of such adsorption.18 The current study investigates whether excess ions, such as Na+ or Ca2+, present in the solution can neutralize the charge of surfactants and allow polymers of the same charge to approach the surface and form complex with the surfactant on water surface. We used the anionic surfactant sodium dodecyl sulfate (SDS), as shown in Figure 1a, and the anionic polymer partially hydrolyzed polyacrylamide (HPAM), as shown in Figure 1b, to explore this cationbridged anionic/anionic interactions. SDS is an anionic surfactant used in many industrial processes and consumer products. Anionic HPAM, a copolymer of poly(acrylamide) (PAM) and poly(acrylic acid) (PAA),19 is widely used for enhanced oil recovery20 and water treatments21-22. HPAM and SDS are also commonly being used as chemical flooding agents to reduce the mobility of the aqueous phase and interfacial tension between water and oil.23-24 It has been reported that the addition of SDS decreases the viscosity of HPAM solutions.25 However, the mechanism for the decrease in viscosity is not fully understood. Samanta et al.26 suggested that SDS affects the viscosity of HPAM solutions through a charge-shielding mechanism but no details of the mechanism were described. Previous studies have shown that HPAM do not form complexes with SDS.25,

11

Only anionic polymers with a very pronounced hydrophobic nature,

such as poly-(1-decene-co-maleic acid) and poly(l-octadecene-co-maleic acid), are able to form

ACS Paragon Plus Environment

4

Page 5 of 27

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

complexes with SDS in the bulk solution.9

Figure 1. The molecular structure of SDS (a) and HPAM (b) in aqueous solutions.

To gain a molecular-level understanding of how ions affect the interaction between the anionic surfactant and the anionic polymer, we carried out a combined study using molecular dynamics (MD) simulation and phase-sensitive sum-frequency generation (SFG) vibrational spectroscopy. As expected, the MD simulation shows that the presence of SDS on the water surface pushes HPAM into the bulk liquid. However, in the presence of excess Na+ ions, the charge screening effect allows HPAM to approach the water surface and weakly interact with the SDS via hydrogen bonding. The addition of divalent cation Ca2+ induces a much stronger interaction between SDS and HPAM by forming ion-bridges and hydrogen bonds. The results of the MD simulation were confirmed by the SFG spectra of the water surface. We observed that introducing excess Na+ ions produced only a small change in the SFG spectrum of the water surface, indicating a minor change in the ordering of surface water molecules. On the other hand, Ca2+ ions significantly decreased the SFG intensity and disturbed the ordering of surface water molecules.

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 27

Material and Methods MD Simulations. We employed full atomistic MD simulations to study the interaction between SDS and HPAM in the presence of NaCl and CaCl2. The all-atom optimized potentials for liquid simulation (OPLS-AA)27 were used for SDS and HPAM. The parameters reported by Li et al.28 were used to describe the intramolecular and intermolecular interactions of SDS molecules. HPAM molecules consisting of 20 monomers with an amine-to-carboxylate functional groups ratio of 4:1 was used in the simulation. The HPAM was parameterized using the OPLS-AA force field,27, 29 which was provided in the GROMACS 5.1.230-32 database. The parameters provided by Aqvist33 were used for NaCl and CaCl2. The behavior of water molecules was described by the extended single point charge (SPC/E) model34. Although the OPLS-AA was originally designed to work with the TIP3P water model, the combination of SPC/E water model with the OPLS-AA force field has been applied successfully to study aqueous interfaces.28,

35-36

Because SDS on the water surface orients interfacial water

molecules,37 a minimum thickness of 5 nm is required to obtain a bulk-like water region. Therefore, a simulation box of 5 nm × 5 nm × 10 nm with two water surfaces was created. Two oligomers of HPAM were initially placed within 2 nm from the water surface on each side. The simulation box was filled with 8067 water molecules. The surface coverage of SDS was fixed at its critical micelle concentration with a surface coverage of 0.44 nm2/molecule,13 which was equivalent to 57 SDS molecules on each side of the water surface. The SDS were randomly distributed on water surface using the PACKMOL package38 such that the hydrophilic head groups were pointing toward the liquid water. To maintain the electrical neutrality of the system, 122 water molecules were randomly substituted by 122 Na+ ions (114 from SDS and 8 from

ACS Paragon Plus Environment

6

Page 7 of 27

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

HPAM). Two vapor regions with a thickness of 25 nm were created on each side of the simulation box. To examine the effect of excess salts, 60 and 90 water molecules were randomly replaced by 30 NaCl and 30 CaCl2, respectively. The steep-descent energy minimization (10000 steps) was used to correct the position of each atom before the simulation. All simulations were performed using GROMACS 5.1.2 GPU computation algorithm30-32 in the canonical ensemble. The temperature of water and the non-water groups were maintained at 293 K independently using the V-rescale39 thermostat with a temperature constant of 0.1 ps. The OH bonds of water were constrained by the SETTLE algorithm40 to enable a simulation time step of 2 fs. The rest of bonds were constrained using the P-LINCS algorithm30 with the LINCS order of 4. The Lennard-Jones interactions were truncated with a cutoff radius of 2.1 nm. Unlikeatom interactions were computed using the geometric combination rule, which the OPLS-AA was designed for. Periodic boundary conditions were applied to all three directions. The cutoff radius for the Coulomb potential was 2.1 nm. The long-range Coulomb interaction was treated with the particle mesh Ewald (PME)41 algorithm with a grid spacing of 0.16 nm. The simulations were carried out for 40 ns, but 30 ns was needed to reach equilibrium. Therefore, the trajectories of the final 10 ns were used for the analysis presented in the current study. The visualization was prepared using VMD 1.9.2.42 Enhanced Sampling Methods. Typical MD simulations suffer from the lack of sampling all possible metastable states, which a system may have. The transition between metastable states often requires overcoming high energy barriers which rarely happens in a typical nanosecond or even microsecond simulation. To overcome this limitation, we carried out well-tempered metadynamics simulations43 in which an external Gaussian potential favorably biases the collective variables (CV) to release the system from the local minima and explore all possible

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 27

states during the simulations. The CV considered in our simulations is the distance between the z-component of the center of mass of the four acrylic acid groups in HPAM to the z-component of the center of mass of all S atoms of SDS. This distance was biased separately for each interface. The biased potential V ( S , t ) is the sum of the Gaussians deposited on the CVs space: d

τf

V (S, t) = ∑W0τ e

(

−V (S ,T )

kB ∆T

)e



∑( i=1

Si ( R)−Si (R(t )))

2

2σ i2

(1)

τ =1

where W0 refers to the initial Gaussian amplitude, kB is the Boltzmann constant, and τ is the time step between bias deposition. The first exponential factor on the right-hand side of eq. (1) scales the amplitude to prevent overfilling the energy landscape. The value of ∆T controls the decay rate of the amplitude in the CVs space and is scaled by the bias factor γ = ( T + ∆T ) T , where T is the system temperature. In the second exponential term, σ i is the width of the Gaussian for the i-th CV. Si(R) represents the i-th CV function. A bias factor of 20 was used. The bias was deposited with a Gaussian width of 0.05 nm, an initial Gaussian amplitude of 2.0 KJ/mol, and a deposition period of 0.4 ps. The well-tempered metadynamics simulations were carried out using the PLUMED 244. To expedite the metadynamics convergence, we put two harmonic potentials with the following form for the distance between the HPAM and the SDS to limit the region of the phase space accessible during the run:

Vκ = κ (z − a)2

(2)

where κ is the wall energy constant and 1500 KJ.mol-1.nm-1 was used, z is the z-component of Cartesian coordinates, and a is the position of the wall. In this study, we set the wall at the distance of 6 nm between the center of mass of the four acrylic acid groups in HPAM to the zcomponent of the center of mass of all S atoms of SDS for each surface. To start the well-

ACS Paragon Plus Environment

8

Page 9 of 27

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

tempered metadynamics, the last trajectory of the typical MD simulation was used as the initial configuration for each system. The temperature of water and non-water groups were maintained at 293 K independently using the V-rescale39 thermostat with a temperature constant of 0.1 ps. The OH bonds of water were constrained by the SETTLE algorithm40 to enable a simulation time step of 2 fs. The rest of bonds were constrained using the P-LINCS algorithm30 with the LINCS order of 4. The Lennard-Jones cutoff radius was 1.0 nm where the interaction was smoothly reduced to 0 after 0.9 nm. The unlike-atom interactions were computed using the geometric combination rule. Periodic boundary conditions were applied to all three directions. The cutoff radius for the Coulomb potential was 1.0 nm. The long-range Coulomb interaction was treated with the particle mesh Ewald (PME)41 algorithm with a grid spacing of 0.16 nm.41 The simulations were carried out for 100 ns, where the convergence was attained. Finally, to calculate the binding free energy between the HPAM and the SDS, the reweighting algorithm45 developed by Tiwary and Parrinello was used. Materials and Sample Preparation. HPAM (average molecular weight ~520,000 g/mol with 80 wt % acrylamide), NaCl (> 99%), CaCl2 (> 99%), and SDS (> 99%) were purchased from Sigma-Aldrich. Water with resistivity > 18.2 MΩ·cm was obtained from a Millipore system. A 0.008M SDS solution was made by dissolving SDS into pure water. A solution with 0.008 M SDS and 10-8 M HPAM was made by adding 0.016 M SDS solution to an equal volume of 2×10-8 M HPAM solution under vigorous stirring. Mixtures of SDS, HPAM, and the salts were made by mixing an equal volume of the SDS/salt solution and the HPAM solution. All solutions were freshly prepared right before spectroscopic measurement. All experiments were performed at 20°C and 1 atm.

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 27

Phase-sensitive SFG setup. A femtosecond Ti-sapphire laser (120 fs, 800 nm, 1 kHz, and 1 mJ/pulse) was used to pump an optical parametric amplifier for generating a femtosecond IR beam. The broad-band IR beam and a narrow-band picosecond 800 nm beam were aligned collinearly.46-48 The incident angle was 60°. A reference SFG was obtained by focusing the IR and picosecond 800 nm beams into a 50 µm thick quartz crystal. The IR, 800 nm, and reference SFG beams were then focused again on the sample. The reference SFG and the SFG generated at the sample went through a time-delay, a polarizer, a band pass filter, a lens, and a monochromator, and then the interference pattern was recorded by a camera. The polarization combination used in this study was SSP (s-polarized SFG, s-polarized 800 nm and p-polarized IR). The energy of the 800 nm and IR beam were ~10 µJ/pulse and ~3 µJ/pulse, respectively. Spectra presented in the paper were acquired over a period of 20 min.

Results and Discussion The snapshots of the MD simulation without excess ion, with excess Na+, and with excess Ca2+ are shown in Figure 2. Without excess ions (Figure 2a), HPAM are not able to overcome the repulsive forces from the surfactants and stay in the bulk water. Excess Na+ ions in the solution led to a weak association between SDS and one of the HPAM polymers (Figure 2b). In the presence of Ca2+ (Figure 2c), both HPAM polymers approach the surface and interact with the SDS layer on water surface.

ACS Paragon Plus Environment

10

Page 11 of 27

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

(a)

(b)

(c)

Figure 2. Simulation box with dimensions of 5.0 × 5.0 × 10.0 nm3 filled with ~8000 water molecules, 114 SDS, 2 HPAM, and 122 Na+. (a) No excess ions. (b) 30 excess NaCl. (c) 30 excess CaCl2. Color codes for atom types: red (oxygen), green (carbon), blue (nitrogen), yellow (sulfur), white (hydrogen), magenta (Na+), orange (Cl-), black (Ca2+), and cyan (water).

In the MD simulations, the weak association between the HPAM and SDS in the presence of excess Na+ often results in a nonsymmetrical structure between the upper and lower water surfaces in Figure 2b as the binding between the SDS and HAPM may break and form repeatedly. To address this issue and to test whether the complex is thermodynamically stable, we calculated the binding free energy between the HPAM and the SDS using the well-tempered metadynamics43 as described in the Enhanced Sampling Methods section. The binding free energy profile for the salt-free system (Figure 3a) shows that the repulsive electrostatic force dominates the interaction between SDS and HPAM. In this case, the HPAM needs to overcome an energy barrier of ~20 KJ/mol to reach the averaged position of the S atoms in SDS, at which the x-axis is defined as 0 in Figure 3. Therefore, both HPAM polymers prefer to stay in the bulk as indicated in Figure 2a. Figure 3b shows that the presence of excess Na+ significantly reduces

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 27

the repulsive force and creates a small and gradual attractive potential for HPAM. The relatively weak potential explains why it is rare to have both HPAM polymers attached to the surfaces. In contrast, Figure 3c shows that the presence of Ca2+ creates a steep potential well of ~ 30 KJ/mol, which is higher than typical hydrogen bonds. Therefore, the Ca2+ bridged HPAM-SDS complex, shown in Figure 2c, is thermodynamically stable.

(a)

(b)

(c)

ACS Paragon Plus Environment

12

Page 13 of 27

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 3. Calculated binding free as function of distance between the z-component of the center of mass of the four acrylic acid groups in HPAM to the z-component of the center of mass of the S atoms in all SDS for the systems with no excess ions (a), 30 excess NaCl (b), and 30 excess CaCl2 (c). The red and blue curves represent results from each water surface in the simulation box. The negative distance refers to states that the position of center of mass of the acrylic acid groups in HPAM is higher than the center of mass of the S atoms.

To better understand the nature of the surfactant-polymer interaction, Figure 4 shows the time-averaged number density profiles of the –SO4 groups in SDS, the H atoms in the NH2 groups of HPAM (HP), the O atoms in the CO2 groups of HPAM (OP), and the intrinsic Na+ ions, which come with the SDS and HPAM. Even without excess ions (Figure 4a), a significant number of the intrinsic Na+ ions accumulate near the surface, but they do not allow HPAM to overcome the strong electrostatic repulsion force. This result is consistent with the previous experimental observation that HPAM and SDS do not associate with each other in pure water.25 Figure 4b shows that the addition of NaCl in the solution results in more Na+ ions accumulating at the surface. The charge screening produced by the additional Na+ ions allows the polymers to approach the surface. In the presence of excess Ca2+ (Figure 4c), nearly all Ca2+ ions are located near the surface, and the distributions of the surfactant and the polymers become overlapped near the surface.

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 27

Figure 4. Time-averaged number density profiles of the –SO4- groups in SDS, the H atoms in the NH2 groups of HPAM (HP), the O atoms in the CO2 groups of HPAM (OP), Na+, Cl-, and Ca2+. (a) no excess ions. (b) 30 excess NaCl. (c) 30 excess CaCl2. The values of HP, OP, Na+, Cl-, and Ca2+ are magnified by 5, 20, 5, 10, and 5 times, respectively for a better visibility.

A detailed analysis of the simulation data reveals that the complex formation between the surfactants and the polymers in the presence of Ca2+ ions is dominated by the Ca2+ bridged

ACS Paragon Plus Environment

14

Page 15 of 27

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

interaction between the OP in HPAM and the ionic O atoms (Oi) of SDS. Figure 5a shows the radial distribution function of Oi with respect to OP in the presence of excess Na+ and Ca2+ ions. It is clear that Ca2+ ions induce a distinguishable layered structure between OP and Oi, (red curve in Figure 5a) while Na+ ions only allow OP and Oi to approach each other (blue curve in Figure 4a). Figure 5b shows that the ester O atoms of SDS (Oe) are not directly involved in the interaction with the polymer as Oe mostly resides in the second layer.

Figure 5. Radial distribution functions between the O atoms in the CO2 groups of HPAM (OP) and the O atoms in SDS: (a) ionic O atoms (Oi); (b) ester O atoms (Oe). OP is the reference point of the radial distribution functions.

Our simulation indicates that the hydrogen bond formation plays an important role in the interaction between HPAM and SDS. Hydrogen bonds may form between the O atoms of SDS

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 27

(both Oi and Oe) and the H atoms in the NH2 group of HPAM (HP). Figure 6 illustrates the radial distribution functions of the Oi and Oe with respect to HP in the presence of Ca2+ and Na+. Overall, Ca2+ promotes the hydrogen bond formation more than Na+ does. In Figure 6, narrower peaks within 1 nm originate from the neighboring H atoms of the polymer (NH2 groups), and the broad peak above 1 nm originate from the O atoms in the neighboring SDS.

Figure 6. Radial distribution functions between the H atoms in the NH2 group of HPAM (HP) and (a) the ionic O atoms (Oi) and (b) the ester O atoms (Oe) of SDS. HP is the reference point of the radial distribution functions.

To quantitatively study the mechanism for the complex formation between SDS and HPAM, the number of ion-bridges were calculated. It was assumed that an ion bridge (OP-ion-Oi

ACS Paragon Plus Environment

16

Page 17 of 27

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

or OP-ion-Oe) is established if an ion (either Na+ or Ca2+) resides within a cutoff distance from the C atom in the CO2 groups of HPAM ( Rcutoff (ion−C ) ) and from the S atom of SDS head groups (

Rcutoff (ion−S ) ). The cutoff distances were determined by the first minimum in the radial distribution function:

Rcutoff (Ca − C ) = 0.45nm ,

Rcutoff ( Na − C ) = 0.40nm ,

Rcutoff (Ca − S ) = 0.44nm ,

Rcutoff ( Na − C ) = 0.40nm , and Rcutoff ( Na − S ) = 0.43nm . The numbers of various types of ionbridges are presented in Table 1. Without salts in the solution, the intrinsic Na+ ions from SDS and HPAM do not produce any ion bridges between SDS and HPAM. The addition of NaCl to the system allows a small amount of Na+ ions to participate in the formation of ion bridges between SDS and HPAM. On the other hand, the presence of Ca2+ significantly facilitates the complex formations between SDS and HPAM. When Ca2+ ion bridges are formed between SDS and HPAM, Na+ ions further stabilized the SDS-HPAM complex as the number of Na+ ion bridges increases significantly. Table1: The number of ion-bridges and hydrogen bond for different systems in terms of the excess ion.

Mediate Ions

Ion bridges (per ns)

Hydrogen bonds (per ns)

Na+ (no salts)

0

0

Na+ (with NaCl)

7

2568

Na+ (with CaCl2)

364

Ca2+ (with CaCl2)

193

3239

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 27

To quantitatively study the hydrogen bond formation between SDS and HPAM, the number of hydrogen bonds was calculated using geometric criteria. A hydrogen bond is formed between HP and Oi (or Oe) if the distance between the N atom of NH2 groups of HPAM and Oi (or Oe) is less than the cutoff distance ( Rcutoff (N − Oi ) = 0.46nm ; Rcutoff (N − Oe ) = 0.60nm ), which was derived from the first minimum in the radial distribution function. Additionally, the angle

∠H...N...O should be less than 30°, which is commonly used to identify hydrogen bonding.49 The number of hydrogen bonds between SDS and HPAM is presented in Table 1. With the addition of NaCl, a significant number of hydrogen bonds are formed although very few Na+ bridges are present indicating that hydrogen bonds are the dominating interaction between SDS and HPAM with the addition of NaCl. While a greater number of hydrogen bonds are formed in the presence of Ca2+, the stronger electrostatic interaction of Ca2+ ions play a more important role in linking SDS and HAPM. Phase-sensitive SFG vibrational spectra were collected at air/water interfaces to verify the MD simulation results. SFG has been known for its high surface sensitivity because under the electric-dipole approximation the second-order optical process is forbidden in a centrosymmetric medium, such as bulk water. While the traditional SFG vibrational spectroscopy measures only the amplitude of the 2nd-order nonlinear optical susceptibility |χ(2)|, phase-sensitive SFG also measures the phase of χ(2).20, 50-52 For the OH symmetric stretch of water surface, the imaginary part of χ(2), Im(χ(2)), can be positive or negative, depending on the sign of the OH projection with respect to the surface normal: a positive peak indicates water molecules with the hydrogen pointing toward the air (up), and a negative peak indicates the OHs pointing toward the bulk (down).53-54 Figure 7 illustrates the Im(χ(2)) of water surface with (a) SDS only, (b) SDS + HPAM, (c) SDS + HPAM + excess Na+, (d) SDS + HPAM + Ca2+, and (e) SDS + Ca2+. As

ACS Paragon Plus Environment

18

Page 19 of 27

shown in Figure 7, because of the negative charge of SDS, SDS (magenta curve) produces an ordered water structure with the water’s OHs pointing up (a positive SFG peak). This observation is in agreement with previous studies revealing that anionic surfactants lead to flipping of water molecules at the surface.47, 53 Adding HPAM with or without excess Na+ does not significantly alter the SFG spectrum. However, the presence of Ca2+ dramatically changes the SFG spectrum of water. The decrease in SFG intensity suggests that the ordering of water is significantly disturbed in the presence of Ca2+. Although Ca2+ ions alone perturb the hydrogen network of surface water molecule (black curve), addition of HPAM further disrupts the ordering of water molecules, which is consistent with the MD simulations suggesting that SDS and HPAM form stable complexes at the surface via ion bridges and hydrogen bonds.

0.8

SDS SDS+HPAM SDS+HPAM+Na

(2)

Im(χ ) (arb. u.)

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+

0.6

SDS+HPAM+Ca 2+

SDS+Ca

0.4

0.2

0.0 3000

3100

3200

3300

3400

3500

-1

Wavenumber (cm ) Figure 7 Im(χ (2)) spectra of air/water interfaces in the OH regions for various aqueous solutions of 0.008M SDS (magenta), 0.008M SDS + 1× 10-8 M HPAM (green), 0.008M SDS + 1× 10-8 M HPAM + 0.1 M excess NaCl (blue), 0.008 M SDS + 1× 10-8 M HPAM + 0.1 M excess CaCl2 (red), and 0.008M SDS + 0.1 M excess CaCl2 (black),.

Conclusion

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 27

We investigated the interaction between surfactants and polymers of the same charge at air/water interfaces using MD simulation and SFG. The results indicated that excess Na+ in the system allowed the polymer to approach the surface and interact with the surfactant via hydrogen bonding. In the presence of Ca2+ ions, HPAM and SDS interact via forming Ca2+ ion bridges and hydrogen bonds. These results are qualitatively consistent with the observed SFG spectra.

Acknowledgments This work was financially supported by the Natural Sciences and Engineering Research Council of Canada. A.M was partially supported by the Four Year Doctoral Fellowship Program at the University of British Columbia.

Reference 1.

Errol Desmond Goddard, K. P. A., Interactions of Surfactants with Polymers and Proteins; CRC Press 1992.

2.

Guzman, E.; Llamas, S.; Maestro, A.; Fernandez-Pena, L.; Akanno, A.; Miller, R.; Ortega, F.; Rubio, R. G., Polymer-Surfactant Systems in Bulk and at Fluid Interfaces. Adv. Colloid Interface Sci. 2016, 233, 38-64.

3.

Pegg, J. C.; Eastoe, J., Solid Mesostructured Polymer-Surfactant Films at the Air-Liquid Interface. Adv. Colloid Interface Sci. 2015, 222, 564-572.

4.

Bureiko, A.; Trybala, A.; Kovalchuk, N.; Starov, V., Current Applications of Foams Formed from Mixed Surfactant-Polymer Solutions. Adv. Colloid Interface Sci. 2015, 222, 670-677.

ACS Paragon Plus Environment

20

Page 21 of 27

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.

Rai, K.; Johns, R. T.; Delshad, M.; Lake, L. W.; Goudarzi, A., Oil-Recovery Predictions for Surfactant Polymer Flooding. J. Pet. Sci. Eng. 2013, 112, 341-350.

6.

Noskov, B. A., Dilational Surface Rheology of Polymer and Polymer/Surfactant Solutions. Curr. Opin. Colloid Interface Sci. 2010, 15, 229-236.

7.

Taylor, D. J. F.; Thomas, R. K.; Penfold, J., Polymer/Surfactant Interactions at the Air/Water Interface. Adv. Colloid Interface Sci. 2007, 132, 69-110.

8.

Miyake, M., Recent Progress of the Characterization of Oppositely Charged Polymer/Surfactant Complex in Dilution Deposition System. Adv. Colloid Interface Sci. 2017, 239, 146-157.

9.

Mcglade, M. J.; Randall, F. J.; Tcheurekdjian, N., Fluorescence Probe Studies of AqueousSolution Interaction between Sodium Dodecyl-Sulfate and Anionic Polyelectrolytes. Macromolecules 1987, 20, 1782-1786.

10.

Mcglade, M. J.; Olufs, J. L., Adsorption of Sodium Dodecyl-Sulfate onto Alpha-Olefin Maleic-Acid Copolymers in Aqueous-Solutions. Macromolecules 1988, 21, 2346-2349.

11.

Iliopoulos, I.; Wang, T. K.; Audebert, R., Viscometric Evidence of Interactions between Hydrophobically Modified Poly(Sodium Acrylate) and Sodium Dodecyl-Sulfate. Langmuir 1991, 7, 617-619.

12.

Zana, R.; Kaplun, A.; Talmon, Y., Microstructural Aspects of Polysoap Sodium DodecylSulfate Interactions. Langmuir 1993, 9, 1948-1950.

13.

Bakeev, K. N.; Ponomarenko, E. A.; Shishanova, T. V.; Tirrell, D. A.; Zezin, A. B.; Kabanov, V. A., Complexation of Amphiphilic Polyelectrolytes with Surfactants of the Same Charge in Water Solutions. Macromolecules 1995, 28, 2886-2892.

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

14.

Page 22 of 27

Bai, G.; Catita, J. A. M.; Nichifor, M.; Bastos, M., Microcalorimetric Evidence of Hydrophobic Interactions between Hydrophobically Modified Cationic Polysaccharides and Surfactants of the Same Charge. J. Phys. Chem. B 2007, 111, 11453-11462.

15.

O'Driscoll, B. M. D.; Milsom, E.; Fernandez-Martin, C.; White, L.; Roser, S. J.; Edler, K. J., Thin Films of Polyethylenimine and Alkyltrimethylammonium Bromides at the Air/Water Interface. Macromolecules 2005, 38, 8785-8794.

16.

Kekicheff, P.; Marcelja, S.; Senden, T. J.; Shubin, V. E., Charge Reversal Seen in Electrical Double-Layer Interaction of Surfaces Immersed in 2-1 Calcium Electrolyte. J. Chem. Phys. 1993, 99, 6098-6113.

17.

Jianjun Liu , Z. X., Jacob Masliyah, Studies on Bitumen−Silica Interaction in Aqueous Solutions by Atomic Force Microscopy. Langmuir 2003, 19, 3911–3920.

18.

Wang, X. F.; Lee, S. Y.; Miller, K.; Welbourn, R.; Stocker, I.; Clarke, S.; Casford, M.; Gutfreund, P.; Skoda, M. W. A., Cation Bridging Studied by Specular Neutron Reflection. Langmuir 2013, 29, 5520-5527.

19.

Morgan, S. E.; Mccormick, C. L., Water-Soluble Polymers in Enhanced Oil-Recovery. Prog. Polym. Sci. 1990, 15, 103-145.

20.

Stiopkin, I. V.; Weeraman, C.; Pieniazek, P. A.; Shalhout, F. Y.; Skinner, J. L.; Benderskii, A. V., Hydrogen Bonding at the Water Surface Revealed by Isotopic Dilution Spectroscopy. Nature 2011, 474, 192-195.

21.

Zhao, P.; Gao, B. Y.; Yue, Q. Y.; Kong, J. J.; Shon, H. K.; Liu, P.; Gao, Y., Explore the Forward Osmosis Performance Using Hydrolyzed Polyacrylamide as Draw Solute for Dye Wastewater Reclamation in the Long-Term Process. Chem. Eng. J. 2015, 273, 316-324.

ACS Paragon Plus Environment

22

Page 23 of 27

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

22.

Sen, G.; Ghosh, S.; Jha, U.; Pal, S., Hydrolyzed Polyacrylamide Grafted Carboxymethylstarch (Hyd. Cms-G-Pam): An Efficient Flocculant for the Treatment of Textile Industry Wastewater. Chem. Eng. J. 2011, 171, 495-501.

23.

Wever, D. A. Z.; Picchioni, F.; Broekhuis, A. A., Polymers for Enhanced Oil Recovery: A Paradigm for Structure-Property Relationship in Aqueous Solution. Prog. Polym. Sci. 2011, 36, 1558-1628.

24.

Cheraghian, G., An Experimental Study of Surfactant Polymer for Enhanced Heavy Oil Recovery Using a Glass Micromodel by Adding Nanoclay. Petrol. Sci. Technol. 2015, 33, 1410-1417.

25.

Methemitis, C.; Morcellet, M.; Sabbadin, J.; Francois, J., Interactions between Partially Hydrolyzed Polyacrylamide and Ionic Surfactants. Eur. Polym. J. 1986, 22, 619-627.

26.

Samanta, A.; Bera, A.; Ojha, K.; Mandal, A., Effects of Alkali, Salts, and Surfactant on Rheological Behavior of Partially Hydrolyzed Polyacrylamide Solutions. J. Chem. Eng. Data 2010, 55, 4315-4322.

27.

Jorgensen, W. L.; Tiradorives, J., The Opls Potential Functions for Proteins - Energy Minimizations for Crystals of Cyclic-Peptides and Crambin. J. Am. Chem. Soc. 1988, 110, 1657-1666.

28.

Li, Z. F.; Van Dyk, A. K.; Fitzwater, S. J.; Fichthorn, K. A.; Milner, S. T., Atomistic Molecular Dynamics Simulations of Charged Latex Particle Surfaces in Aqueous Solution. Langmuir 2016, 32, 428-441.

29.

Jorgensen, W. L.; Maxwell, D. S.; TiradoRives, J., Development and Testing of the Opls All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids. J. Am. Chem. Soc. 1996, 118, 11225-11236.

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

30.

Page 24 of 27

Hess, B.; Kutzner, C.; van der Spoel, D.; Lindahl, E., Gromacs 4: Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. J. Chem. Theory Comput. 2008, 4, 435-447.

31.

Pronk, S., et al., Gromacs 4.5: A High-Throughput and Highly Parallel Open Source Molecular Simulation Toolkit. Bioinformatics 2013, 29, 845-854.

32.

Abraham, M. J.; Murtola, T.; Schulz, R.; Páll, S.; Smith, J. C.; Hess, B.; Lindahl, E., Gromacs: High Performance Molecular Simulations through Multi-Level Parallelism from Laptops to Supercomputers. SoftwareX 2015, 1–2, 19-25.

33.

Aqvist, J., Ion Water Interaction Potentials Derived from Free-Energy Perturbation Simulations. J. Phys. Chem-Us 1990, 94, 8021-8024.

34.

Berendsen, H. J. C.; Grigera, J. R.; Straatsma, T. P., The Missing Term in Effective Pair Potentials. J. Phys. Chem-Us 1987, 91, 6269-6271.

35.

Stephenson, B. C.; Goldsipe, A.; Beers, K. J.; Blankschtein, D., Quantifying the Hydrophobic Effect. 1. A Computer Simulation-Molecular-Thermodynamic Model for the Self-Assembly of Hydrophobic and Amphiphilic Solutes in Aqueous Solution. J. Phys. Chem. B 2007, 111, 1025-1044.

36.

Yan, H.; Guo, X. L.; Yuan, S. L.; Liu, C. B., Molecular Dynamics Study of the Effect of Calcium Ions on the Monolayer of Sdc and Sdsn Surfactants at the Vapor/Liquid Interface. Langmuir 2011, 27, 5762-5771.

37.

Hu, D.; Mafi, A.; Chou, K. C., Revisiting the Thermodynamics of Water Surfaces and the Effects of Surfactant Head Group. J. Phys. Chem. B 2016, 120, 2257-2261.

ACS Paragon Plus Environment

24

Page 25 of 27

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

38.

Martinez, L.; Andrade, R.; Birgin, E. G.; Martinez, J. M., Packmol: A Package for Building Initial Configurations for Molecular Dynamics Simulations. J. Comput. Chem. 2009, 30, 2157-2164.

39.

Bussi, G.; Donadio, D.; Parrinello, M., Canonical Sampling through Velocity Rescaling. J. Chem. Phys. 2007, 126, 01401.

40.

Miyamoto, S.; Kollman, P. A., Settle - an Analytical Version of the Shake and Rattle Algorithm for Rigid Water Models. J. Comput. Chem. 1992, 13, 952-962.

41.

Essmann, U.; Perera, L.; Berkowitz, M. L.; Darden, T.; Lee, H.; Pedersen, L. G., A Smooth Particle Mesh Ewald Method. J. Chem. Phys. 1995, 103, 8577-8593.

42.

Humphrey, W.; Dalke, A.; Schulten, K., Vmd: Visual Molecular Dynamics. J. Mol. Graph. Model. 1996, 14, 33-38.

43.

Barducci, A.; Bussi, G.; Parrinello, M., Well-Tempered Metadynamics: A Smoothly Converging and Tunable Free-Energy Method. Phys. Rev. Lett. 2008, 100, 020603.

44.

Tribello, G. A.; Bonomi, M.; Branduardi, D.; Camilloni, C.; Bussi, G., Plumed 2: New Feathers for an Old Bird. Comput. Phys. Commun. 2014, 185, 604-613.

45.

Tiwary, P.; Parrinello, M., A Time-Independent Free Energy Estimator for Metadynamics. J. Phys. Chem. B 2015, 119, 736-742.

46.

Ji, N.; Ostroverkhov, V.; Chen, C. Y.; Shen, Y. R., Phase-Sensitive Sum-Frequency Vibrational Spectroscopy and Its Application to Studies of Interfacial Alkyl Chains. J. Am. Chem. Soc. 2007, 129, 10056-10057.

47.

Dai, C.; Zhao, J.; Yan, L.; Zhao, M., Adsorption Behavior of Cocamidopropyl Betaine under Conditions of High Temperature and High Salinity. J. Appl. Polym. Sci. 2014, 131, 40424.

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

48.

Page 26 of 27

Chanda, J.; Bandyopadhyay, S., Molecular Dynamics Study of a Surfactant Monolayer Adsorbed at the Air/Water Interface. J. Chem. Theory Comput. 2005, 1, 963-971.

49.

Durrant, J. D.; McCammon, J. A., Hbonanza: A Computer Algorithm for MolecularDynamics-Trajectory Hydrogen-Bond Analysis. J. Mol. Graph. Model. 2011, 31, 5-9.

50.

Nihonyanagi, S.; Mondal, J. A.; Yamaguchi, S.; Tahara, T., Structure and Dynamics of Interfacial Water Studied by Heterodyne-Detected Vibrational Sum-Frequency Generation. Annu. Rev. Phys. Chem. 2013, 64, 579-603.

51.

Jiao, J., Polyoxyethylated Nonionic Surfactants and Their Applications in Topical Ocular Drug Delivery. Adv. Drug Deliv. Rev. 2008, 60, 1663-1673.

52.

Shen, Y. R., Basic Theory of Surface Sum-Frequency Generation (Vol 116, Pg 15505, 2012). J. Phys. Chem. C 2013, 117, 11884-11884.

53.

Nihonyanagi, S.; Yamaguchi, S.; Tahara, T., Direct Evidence for Orientational Flip-Flop of Water Molecules at Charged Interfaces: A Heterodyne-Detected Vibrational Sum Frequency Generation Study. J. Chem. Phys. 2009, 130, 204704.

54.

Hu, D.; Yang, Z.; Chou, K. C., Interactions of Polyelectrolytes with Water and Ions at Air/Water Interfaces Studied by Phase-Sensitive Sum Frequency Generation Vibrational Spectroscopy. J. Phys. Chem. C 2013, 117, 15698-15703.

ACS Paragon Plus Environment

26

Page 27 of 27

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 Graphic

ACS Paragon Plus Environment

27