The Topmost Water Structure at a Charged Silica ... - ACS Publications

Jul 5, 2018 - structure at the negatively charged silica/aqueous interface at pH 12 using .... hereafter call this orientation “H-up orientation” ...
0 downloads 0 Views 2MB Size
Subscriber access provided by BOSTON COLLEGE

Surfaces, Interfaces, and Catalysis; Physical Properties of Nanomaterials and Materials

The Topmost Water Structure at a Charged Silica/Aqueous Interface Revealed by Heterodyne-Detected Vibrational Sum Frequency Generation Spectroscopy Shu-hei Urashima, Anton Myalitsin, Satoshi Nihonyanagi, and Tahei Tahara J. Phys. Chem. Lett., Just Accepted Manuscript • DOI: 10.1021/acs.jpclett.8b01650 • Publication Date (Web): 05 Jul 2018 Downloaded from http://pubs.acs.org on July 7, 2018

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

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 28 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

The Journal of Physical Chemistry Letters

The Topmost Water Structure at a Charged Silica/Aqueous Interface Revealed by Heterodyne-Detected Vibrational Sum Frequency Generation Spectroscopy Shu-hei Urashima, †1 Anton Myalitsin,‡1 Satoshi Nihonyanagi,1,2 Tahei Tahara1,2* 1

Molecular Spectroscopy Laboratory, RIKEN, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan

2

Ultrafast Spectroscopy Research Team, RIKEN Center for Advanced Photonics (RAP), 2-1

Hirosawa, Wako, Saitama 351-0198, Japan.

Present addresses †

Research Institute for Science & Technology, Tokyo University of Science, 2641 Yamazaki,

Noda, Chiba 278-8510, Japan. ‡

Function Analysis Department, NISSAN ARC Ltd., 1 Natsushima, Yokosuka, Kanagawa

237-0061, Japan. Corresponding Author Corresponding Author: [email protected]

ACS Paragon Plus Environment

1

The Journal of Physical Chemistry Letters 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 28

Abstract

Despite recent significant advance of interface-selective nonlinear spectroscopy, the topmost water structure at the charged silica surface is still not clearly understood. This is because, for charged interfaces, not only interfacial molecules at the topmost layer but also a large number of molecules in the electric double layer are probed even with the second-order nonlinear spectroscopy. In the present study, we studied water structure at the negatively charged silica/aqueous interface at pH 12 using heterodyne-detected vibrational sum frequency generation spectroscopy, and demonstrated that the spectral component of the topmost water can be extracted by examining the ionic strength dependence of the Imχ(2) spectrum. The obtained Imχ(2) spectrum indicates that the dominant water species in the topmost layer is hydrogen-bonded to the negatively charged silanolate at the silica surface with one OH group. There also exists minor water species that weakly interacts with the oxygen atom of a siloxane bridge or the remaining silanol at the silica surface, using one OH group. The ionic strength dependence of the Imχ(2) spectrum indicates that this water structure of the topmost layer is unchanged in a wide ionic strength range from 0.01 M to 2 M.

ACS Paragon Plus Environment

2

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

The Journal of Physical Chemistry Letters

Silica is one of the most fundamental and abundant minerals on the earth. It is known that the silica surface at aqueous interfaces is negatively charged under usual environmental conditions (pH > 2) 1-2 because of deprotonation of the surface silanol (Si-OH → Si-O- + H+). Therefore, the charged silica surface is of fundamental and industrial importance, and a number of studies have been conducted to elucidate the surface properties such as pKa of the surface deprotonation reaction, electric double layer (EDL) structure, and water structure in the EDL.1-17 However, our molecular-level understanding on the structure of the silica/aqueous interface is still insufficient, mainly because of experimental limitations. One of the reasons of the poor understanding of interfacial structure is the difficulty in probing interfacial molecules selectively. Even-order nonlinear spectroscopy, such as vibrational sum frequency generation (VSFG),12, 18 second harmonic generation (SHG),19 and fourth-order coherent Raman,20-21 can selectively monitor interfacial molecules, and they have been intensively applied to various interfaces. In particular, the heterodyne-detected (HD-) VSFG spectroscopy enables us to experimentally determine the complex second-order susceptibility  () , whose imaginary part (Im () ) directly exhibits an absorptive line shape at the resonance as IR and Raman spectra in the bulk do.22-24 The HD-VSFG technique was initially applicable only to the gas/liquid and gas/solid interfaces, but it has successfully been applied to the solid/liquid interfaces very recently.25-27 In our previous study, we measured Im () spectra of silica/aqueous interfaces by changing pH of the solution phase, and revealed the orientation and structure of interfacial water at different pHs.25 Nevertheless, when the interface is charged at pH > 2, the  () spectra include not only the contribution from the water molecules in the close vicinity of the interface but also that in a substantially far region because the charge at the interface creates the electric field and

ACS Paragon Plus Environment

3

The Journal of Physical Chemistry Letters 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 28

orients the water molecules in the EDL. This means that the experimentally observed  () ()

()

spectra,  , are the sum of the  () components arising from the topmost interface ( ) ()

()

()

and the EDL ( ).1, 28-29 ( and  can also be called the contributions from the Stern ()

layer and Gouy-Chapman layer, respectively.30) Therefore, it is desirable to separate the  ()

component from the  for clear discussion about the interfacial structure, particularly at high pH. Here, we note that the signal from the EDL (or Gouy-Chapman layer) is often represented as  () and is termed “third-order” nonlinear susceptibility1, 4, 14-15, 29 because this component arises under the existence of three electric fields (i.e., the electric fields of the two incident lights and one static field generated by the surface charge). However, sum-frequency generation in the Gouy-Chapman layer predominantly arises from second-order optical process, and the role of the static field is to make water molecules in the EDL oriented as clearly shown in recent ()

()

theoretical studies.31-32 Therefore, we describe this component as  , not  , in this paper. ()

()

For obtaining the  component, Wen et al. determined the  spectrum from subtly charged fatty acid/water interfaces and then subtracted it from the spectra of charged monolayer/water interfaces using the factor estimated from the surface charge density.28 They ()

successfully obtained  , but it is not easy to apply this method to other interfaces because it is difficult to estimate the static electric field in the EDL. A simple strategy for experimentally ()

obtaining the  spectrum is to measure  () spectra with changing the ionic strength of the solution. Because the thickness of the EDL (i.e., Debye length) is largely dependent on the ionic strength, comparison of the spectra measured at different ionic strengths enables us to separate ()

()

 and  . This strategy does not require any rigorous knowledge about the charge density at the interface, and hence it can be a very powerful method, in particular for amorphous silica

ACS Paragon Plus Environment

4

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

The Journal of Physical Chemistry Letters

surface. This is because evaluation of the surface charge density at the amorphous silica surface is model-dependent1,

7

and its surface structure is changed by pretreatment.33-34 In fact, this

strategy has been employed in previous homodyne-detected VSFG studies,13,

15

and it was

proposed that strongly hydrogen-bonded water species and bulk-like hydrogen-bonded species are present in the EDL and at the topmost layer, respectively. However, it is often difficult to correctly interpret conventional homodyne-detected VSFG spectra that correspond to  () 



()

because of the spectral distortion. For example, in the case that two components ( and ()

 ) are overlapped and coexist in the same spectral region, it is obvious that the cross term make it very difficult to separate them: ()



()

()



() 

()



() ()∗

χ  = χ + χ  = χ  + χ  + 2Reχ χ . Thus, it is necessary to obtain Imχ(2) spectra of charged silica/water interfaces using HD-VSFG ()

()

for separating the  and  components unambiguously and for elucidating the interfacial water structure based on the observed spectra. In this study, while systematically changing ionic strength of the aqueous phase, we measured Im () spectra of the charged silica/aqueous interface at pH 12 using both of pure ()

water (H2O) and isotopically diluted water (HOD in D2O) in the aqueous phase. The Im spectrum obtained by HD-VSFG measurements is significantly different from that estimated 

from the  ()  spectra obtained by conventional VSFG measurements,13, 15 and it reveals the existence of characteristic asymmetric hydrogen bonding of water molecules at the topmost layer. Figure 1a shows Im () spectra of the silica/H2O aqueous interfaces measured at pH 12 with various ionic strengths ranging from 0.01 to 5 M, which were prepared with addition of NaCl

ACS Paragon Plus Environment

5

The Journal of Physical Chemistry Letters 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 28

Figure 1. (a) Im () spectra of silica/H2O aqueous interfaces at pH 12 with ionic strengths from 0.01 to 5 M. (b) The difference spectra calculated from the spectra shown in (a). The difference spectra are area-normalized for comparison of the spectral shape.

into NaOH solution. At this pH, deprotonation of the silanol group is saturated,3 so that the surface charge density (i.e., the surface density of silanolate groups) is considered to be independent of the ionic strength of the solution. Note that the lowest ionic strength studied here (0.01 M) is due to NaOH that is added for setting pH at 12. The spectra at low ionic strengths appear positive in the entire OH stretch region, indicating that the interfacial water molecules are net-oriented with their hydrogen atoms pointing toward the silica (we hereafter call this orientation “H-up orientation” according to our experimental configuration where the silica

ACS Paragon Plus Environment

6

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

The Journal of Physical Chemistry Letters

substrate is placed on the upper side over the solution25). This observation is consistent with previous reports for low ionic strengths.25, 35 At a glance, the OH band intensity monotonically decreases with increasing the ionic strength, and the spectral change is almost saturated at the ionic strength of 2 M. Following the discussions in previous homodyne-VSFG studies, the ()

decrease in the OH band intensity is attributed to the decrease of the  contribution owing to the screening of the electric field by the counter cation (Na+). The Im () spectra above 2 M, ()

where the  contribution is sufficiently suppressed, are predominantly attributable to the ()

topmost water layer,  .14, 36-38 It is noteworthy that the spectrum shows small but noticeable changes even above 2 M. This change is ascribable to the subtle modification of the water structure in the topmost region, which we discuss in detail later. For examining the observed spectral change more quantitatively, we calculated the difference spectrum for each increase of the ionic strength, which is shown after area-normalization in Figure 1b. The difference spectra shown in this figure are well overlapped with each other in the ionic strength range of 0.01 - 2 M. This indicates that the spectral change in this range is ascribable simply to the decrease of one spectral component. Obviously, this decreasing ()

component, i.e. the difference spectra, is attributable to Im because at lower ionic strengths, e.g., 0.01 - 0.105 M, only shrinkage of EDL is expected according to the calculated field gradient as well as counter ion distribution (see Figure S1 in the supporting information (SI)). Important ()

finding here is that such simple decrease of the  component continues up to the ionic ()

strength as high as 2 M. This means that the spectral shape of  , and hence the water structure at the topmost layer, is unchanged in such a wide range of ionic strengths from 0.01 to ()

2 M. The contribution of  component is considered sufficiently suppressed at 2 M because

ACS Paragon Plus Environment

7

The Journal of Physical Chemistry Letters 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 28

the spectra in the 3000 - 3300 cm-1 region no longer change at the ionic strength higher than 2 M. ()

Thus, the  () spectrum at 2 M can be regarded as the “ spectrum”, which is commonly included in all the spectra measured at the ionic strength less than 2 M. On the other hand, we ()

refer to the difference spectrum between the “0.01 M” and “1 M” spectra as “ spectrum”. (The 1 M spectrum is chosen simply because of the low noise level especially in the HOD spectrum which is to be compared in the following section.) At the ionic strengths higher than 2 M, the Imχ(2) spectra exhibit a subtle but noticeable change, in particular the growth of the negative OH band at around 3500 cm-1 (Figure 1a). This change is more readily recognized in the “2 M - 5 M” difference spectrum shown in Figure 1b as a narrow positive OH band at around this frequency. This spectral change in the high frequency region implies that the topmost water structure slightly changes at the ionic strength higher than 2 M, presumably due to the presence of the dense ions near the interface. ()

()

The Im and Im spectra obtained by the above-described spectral decomposition ()

()

are shown in Figure 2. In this figure, we also show the Im and Im spectra measured with isotopically diluted water (HOD-D2O; H2O : HOD : D2O = 1 : 8 : 16) for showing the effect of the inter- and intramolecular vibrational couplings which substantially affect H2O spectra.25, 39-41

(The ionic strength dependent spectra for HOD-D2O solutions are shown in Figure S2 in SI.) ()

As shown in Figure 2a, the Im spectrum of the HOD-D2O solution interface exhibits the OH band peaked at around 3450 cm-1, while the corresponding H2O spectrum shows a substantially broader and red-shifted OH band. This is a typical difference between the H2O and HOD-D2O vibrational spectra of water, which arises from the vibrational coupling existing in ()

H2O. The Im spectra are similar to vibrational spectra of the corresponding aqueous

ACS Paragon Plus Environment

8

Page 9 of 28 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

The Journal of Physical Chemistry Letters

()

()

Figure 2. Im (a) and Im (b) spectra obtained with H2O (black) and HOD-D2O ()

(red) solutions. The  spectra were obtained as the difference of the spectra measured ()

with ionic strength of 0.01 M and 1 M, whereas the  spectra were as-observed spectra measured with ionic strength of 2 M.

solutions in the bulk,

42-43

indicating that the water structure in the EDL is bulk-like, as

previously reported for charged monolayer interfaces.28, 39 ()

()

In contrast to the Im spectra, the Im spectra of H2O and HOD-D2O are very similar to each other, and they exhibit a broad positive peak at around 3200 cm-1 (Figure 2b),

ACS Paragon Plus Environment

9

The Journal of Physical Chemistry Letters 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 28

()

which is markedly different from the bulk water spectra. The Im spectra provide new insights into the topmost water structure at the charged silica interface. First, the OH band is positive and substantially red-shifted. This means that the OH group giving rise to this band is H-up oriented and forms a strong hydrogen-bond with the silica surface at pH 12. Possible hydrogen-bond accepting sites of the silica surface are the oxygen of silanolate (Si-O-), siloxane bridge (Si-O-Si), and the remaining silanol (Si-OH). Among these three, only silanolate decreases its population at lower pH, and the positive peak around 3200 cm-1 was in fact absent in the spectrum measured at pH 2.25 Thus, the positive broad band around 3200 cm-1 is safely assignable to the stretch mode of the H-up oriented OH group that forms a strong hydrogen-bond ()

with a negatively charged silanolate of the silica surface. Second, the similarity between Im spectra of H2O and HOD-D2O solutions indicates that the stretching vibrations of the two OH groups in a H2O molecule are decoupled at the topmost layer. This indicates an asymmetric hydrogen-bonding environment for the topmost water, i.e., two OH groups of the H2O molecule donate their hydrogen to different types of hydrogen-bond accepting sites so that they vibrate with substantially different frequencies. It means that a H2O molecule at the topmost layer forms a very strong hydrogen bond with surface silanolate with one OH group and forms a substantially different hydrogen bond by another OH with something else, e.g., another water ()

molecule. The Im spectra also exhibit a small negative band at around 3500 cm-1 and a small positive peak around 3600 cm-1. Although these bands are very weak, they are not noise because similar negative-positive features are recognized in all the spectra measured in the ionic strength above 0.3 M (Figure 1a). The weak positive peak around 3600 cm-1 can be assigned to the H-up OH of the topmost water molecule that weakly interacts with oxygen atom of siloxane bridge (Si-O-Si) or remaining silanol (Si-OH), as previously concluded based on the

ACS Paragon Plus Environment

10

Page 11 of 28 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

The Journal of Physical Chemistry Letters

pH-dependent change of the Im () spectrum.25 Then, the negative band around 3500 cm-1 is attributable to the stretch vibration of the H-down OH (H atom is pointing toward the bulk) of the water molecule whose another H-up OH interacts with silica surface, i.e., the oxygen of silanolate (Si-O-), siloxane bridge (Si-O-Si), or the remaining silanol (Si-OH). This H-down OH forms a hydrogen bond with other water molecules or chloride anions locating in a deeper region, making the hydrogen bonds of two OHs of the topmost water molecules highly asymmetric. We note that the vibrational decoupling of water has also been observed for positively charged surfactant/aqueous halide salt interfaces, and it was attributed to the symmetry breaking induced by a hydrogen bond with chloride counterion.44 For further discussion, we quantitatively examined the ionic strength dependence of the ()

spectral intensity. According to the literature, it is expected that the amplitude of the  linearly correlates with the surface potential (see SI for the details),1, 28 i.e., ()

()

()

()

()

 =  +  =  +   , ()

where 

(1)

is the  () response of the water in the EDL per unit surface potential  . In ()

Figure 3, we plot the experimentally evaluated  amplitude (the OH band intensity integrated from 3000 to 3300 cm-1) against the surface potential  calculated by so-called modified Poisson-Boltzmann theory (the detail of the theory is given in SI).28, 45 This plot clearly ()

shows that the amplitude of  is linear to the calculated surface potentials in the ionic strength range lower than 2 M. This excellent agreement between experiment and the (modified) Poisson-Boltzmann theory shows that the classical EDL picture is valid for a surprisingly wide ionic strength range, at least for the silica interfaces.

ACS Paragon Plus Environment

11

The Journal of Physical Chemistry Letters 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 28

()

Figure 3. Intensity of the OH stretch band ( ) plotted against the surface potential ( ). ()

 was evaluated as the OH band intensity integrated from 3000 to 3300 cm-1, and  was calculated with modified Poisson-Boltzmann theory. The black dotted line is an eye guide, which is extrapolated from the points of the ionic strength between 0.01 and 0.105 M.

()

The contribution of  in the  () spectra of the charged silica/aqueous interfaces has been discussed experimentally11,

14-15, 37, 46

and theoretically.31-32 Most of the previous

experimental studies on the ionic strength dependence were performed at pH ~6 with conventional homodyne-detected VSFG, and they were essentially consistent to each other: the ionic strength dependence of VSFG intensity looks saturated at ionic strength around 0.1 M.13-14, 47

This observation had been interpreted as the EDL is sufficiently compressed by the ionic

strength of 0.1 M. Thus, the previous reports appear to contradict the result of the present study ()

that clearly shows the continuous decrease of the  component up to the ionic strength of 2 M. We think that this apparent inconsistency is likely attributable to the different pH employed in the experiments, as well as the difficulty of interpreting conventional homodyne-detected

ACS Paragon Plus Environment

12

Page 13 of 28 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

The Journal of Physical Chemistry Letters

VSFG spectra. In our previous HD-VSFG study25, we observed negative and positive peaks in the OH stretch region of the Im () spectra measured at pH < 7, revealing that “H-down” and “H-up” oriented water coexist at the topmost layer at less deprotonated or neutral silica/water ()

interface. Because the  component is positive in the entire OH stretch region as shown in ()

()

this work, destructive interference between the  and  components may give rise to the 

ionic strength range where the  ()  spectrum looks rather insensitive to the change of ionic strength although the suppression of the EDL (and hence the spectral change of Im () ) is still continuing. In fact, our preliminary experiments support that this is the case (not shown). We also note that some previous homodyne-detected VSFG studies proposed that water in the topmost layer and those in EDL show OH stretch bands at around 3400 cm-1 and 3200 cm-1, respectively, at neutral pH13 as well as a wide pH range of 2-12.15 However, direct measurements of Im () spectra in the present study show that this is not the case: the OH band of the topmost water is peaked around 3200 cm-1 whereas the OH stretch band of water in EDL is peaked at a higher frequency, i.e., 3200-3500 cm-1 for H2O and 3450 cm-1 for HOD-D2O. This highlights the importance of the direct measurement of Im () spectra, which is realized with the heterodyne-detection in VSFG measurements. Finally, we mention two observations that provide information about the counterion distribution in the topmost water layer at the silica/water interface. First, as already described, the obtained Im () spectra show that the topmost water structure is unchanged in the ionic strength range of 0.01-2 M. This implies that one of the followings is the case: The topmost water structure is insensitive to the presence of Na+, or the interfacial density of Na+ is saturated at 0.01 M and unchanged at higher ionic strength. The modified Poisson-Boltzmann theory indicates that the latter is more likely because the local concentration of the hydrated cation at

ACS Paragon Plus Environment

13

The Journal of Physical Chemistry Letters 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 28

the topmost layer is calculated to be nearly saturated even at the ionic strength as low as 0.01 M (Figure S1 in SI). Second, the Im () signal does not vanish even at a very high ionic strength such as 5 M. This indicates that Na+ does not displace the topmost water molecules even at the ionic strength of 5 M, as suggested previously.36 This means that the contact adsorption of Na+ does not increase with ionic strength above 0.01 M, i.e., it is already saturated below 0.01 M even if the contact adsorption exists. In fact, a good correlation is seen between the OH band intensity and the surface potential calculated by the theory in the ionic strength range of 0.01-2 M as shown in Figure 3, which is consistent with the absence of the noticeable change in the contact adsorption in this range. In summary, we measured the Im () spectra of the charged silica/aqueous interface at a wide range of the ionic strength with HD-VSFG, and demonstrated that the spectra can be experimentally decomposed into two components, i.e., the Im () spectrum of the topmost layer and that of the EDL, on the basis of the ionic strength dependence of the spectra. The isotopic dilution experiment revealed that two OH stretch vibrations of the topmost water molecule are decoupled, indicating that two OH bonds of the water form significantly different hydrogen bonds. The concluded water structure at the charged silica/aqueous interface at pH 12 is sketched in Figure 4: The dominant water species in the topmost layer is hydrogen-bonded to the negatively charged silanolate at the silica surface with one OH group which is H-up oriented (H atom is pointing toward the silica) (A). There also exists a minor water species that weakly interacts with the oxygen atom of a siloxane bridge or the remaining silanol using its one OH group with H-up orientation (B). The other OH group of these water molecules in the topmost layer is H-down oriented (H atom is pointing toward the bulk aqueous phase), forming hydrogen bonds with the surrounding water or Cl- located in slightly deeper aqueous phase (C). This water

ACS Paragon Plus Environment

14

Page 15 of 28 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

The Journal of Physical Chemistry Letters

Figure 4. Proposed water structures at the silica/aqueous interface at bulk pH of 12 at the ionic strength of 0.01-2 M. Only the water molecules observable by HD-VSFG spectroscopy are depicted. The red dotted lines represent hydrogen-bonding. The letters A, B, and C represent the OH groups which interact with silanolate (A), siloxane bridge or silanol (B), and chloride anion or another water molecule (C).

structure in the topmost layer is unchanged in the ionic strength range between 0.01 M and 2 M. We did not find any experimental evidence suggesting contact adsorption of Na+ onto the charged silica surface in the ionic strength region examined in the present study (0.01 - 5 M). On the other hand, the water molecules in the EDL are net “H-up” orientated due to the negative charge at the silica surface, forming bulk-like hydrogen bonds with surrounding water molecules. The intensity of the Im () signal from the water in the EDL is almost perfectly correlated with the surface potential (which is calculated by modified Poisson-Boltzmann theory) up to 2 M.

ACS Paragon Plus Environment

15

The Journal of Physical Chemistry Letters 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 28

Lastly, we wish to clearly mention the significant improvement of our molecular-level picture on the water structure at the charged silica/water interface. In our previous study on pH dependence of water structure at the silica/water interface, we could only conclude that interfacial water is H-up oriented at pH 12.25 In the present study, on the other hand, the asymmetric hydrogen-bonding structure of the topmost water molecules was clarified for the first ()

()

time, thanks to the separation of  and  . The experimental approach employed in this study opens a new way for systematic elucidation of interfacial water structures at mineral interfaces using interface-selective nonlinear spectroscopy.

Supporting Information. (1) Experimental setup, (2) brief description for the derivation of eq. 1, (3) brief descriptions of the modified Poisson-Boltzmann theory and the EDL properties calculated by the theory, (4) Im () spectra of the silica / HOD-D2O solution interfaces, (5) details of the Fresnel analysis, 

and (6)  ()  spectra reconstructed from the complex  () spectra.

AUTHOR INFORMATION Corresponding author *Email: [email protected]; Tel: +81-48-467-4592. Present addresses S.U.: Research Institute for Science & Technology, Tokyo University of Science, 2641 Yamazaki, Noda, Chiba 278-8510, Japan.

ACS Paragon Plus Environment

16

Page 17 of 28 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

The Journal of Physical Chemistry Letters

A. M.: Function Analysis Department, NISSAN ARC Ltd., 1 Natsushima, Yokosuka, Kanagawa 237-0061, Japan. Notes The authors declare no competing financial interests.

ACKNOWLEDGMENT This work was supported by JSPS KAKENHI Grant Numbers JP25104005, JP18H03905, and 18H05265.

References (1)

Ong, S.; Zhao, X.; Eisenthal, K. B., Polarization of Water Molecules at a Charged

Interface: Second Harmonic Studies of the Silica/Water Interface. Chem. Phys. Lett. 1992, 191, 327-335. (2)

Darlington, A. M.; Gibbs-Davis, J. M., Bimodal or Trimodal? The Influence of Starting

pH on Site Identity and Distribution at the Low Salt Aqueous/Silica Interface. J. Phys. Chem. C 2015, 119, 16560-16567. (3)

Azam, M. S.; Darlington, A.; Gibbs-Davis, J. M., The Influence of Concentration on

Specific Ion Effects at the Silica/Water Interface. J. Phys. Condens. Matter. 2014, 26, 244107. (4)

Azam, M. S.; Weeraman, C. N.; Gibbs-Davis, J. M., Specific Cation Effects on the

Bimodal Acid-Base Behavior of the Silica/Water Interface. J. Phys. Chem. Lett. 2012, 3, 1269-1274.

ACS Paragon Plus Environment

17

The Journal of Physical Chemistry Letters 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

(5)

Page 18 of 28

Scales, P. J.; Grieser, F.; Healy, T. W.; White, L. R.; Chan, D. Y. C., Electrokinetics of

the Silica-Solution Interface: a Flat Plate Streaming Potential Study. Langmuir 1992, 8, 965-974. (6)

Brown, M. A.; Abbas, Z.; Kleibert, A.; Green, R. G.; Goel, A.; May, S.; Squires, T. M.,

Determination of Surface Potential and Electrical Double-Layer Structure at the Aqueous Electrolyte-Nanoparticle Interface. Phys. Rev. X 2016, 6, 011007. (7)

Brown, M. A.; Goel, A.; Abbas, Z., Effect of Electrolyte Concentration on the Stern

Layer Thickness at a Charged Interface. Angew. Chem. Int. Ed. 2016, 55, 3790-3794. (8)

Totland, C.; Nerdal, W., Experimental Determination of Water Molecular Orientation

near a Silica Surface Using NMR Spectroscopy. J. Phys. Chem. C 2016, 120, 5052-5058. (9)

Morag, J.; Dishon, M.; Sivan, U., The Governing Role of Surface Hydration in Ion

Specific Adsorption to Silica: an AFM-Based Account of the Hofmeister Universality and Its Reversal. Langmuir 2013, 29, 6317-6322. (10)

Li, R.; Todd, B. A., Multivalent Ion Screening of Charged Glass Surface Studied by

Streaming Potential Measurements. J. Chem. Phys. 2013, 139, 194704. (11)

Du, Q.; Freysz, E.; Shen, Y. R., Vibrational Spectra of Water Molecules at

Quartz/Water Interfaces. Phys. Rev. Lett. 1994, 72, 238-241. (12)

Shen, Y. R.; Ostroverkhov, V., Sum-Frequency Vibrational Spectroscopy on Water

Interfaces: Polar Orientation of Water Molecules at Interfaces. Chem. Rev. 2006, 106, 1140-1154. (13)

Jena, K. C.; Hore, D. K., Variation of Ionic Strength Reveals the Interfacial Water

Structure at a Charged Mineral Surface. J. Phys. Chem. C 2009, 113, 15364-15372.

ACS Paragon Plus Environment

18

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

The Journal of Physical Chemistry Letters

(14)

Jena, K. C.; Covert, P. A.; Hore, D. K., The Effect of Salt on the Water Structure at a

Charged Solid Surface: Differentiating Second- and Third-order Nonlinear Contributions. J. Phys. Chem. Lett. 2011, 2, 1056-1061. (15)

Darlington, A. M.; Jarisz, T. A.; DeWalt-Kerian, E. L.; Roy, S.; Kim, S.; Azam, M. S.;

Hore, D. K.; Gibbs, J. M., Separating the pH-Dependent Behavior of Water in the Stern and Diffuse Layers with Varying Salt Concentration. J. Phys. Chem. C 2017, 121, 20229-20241. (16)

Flores, S. C.; Kherb, J.; Konelick, N.; Chen, X.; Cremer, P. S., The Effects of

Hofmeister Cations at Negatively Charged Hydrophilic Surfaces. J. Phys. Chem. C 2012, 116, 5730-5734. (17)

Dewan, S.; Carnevale, V.; Bankura, A.; Eftekhari-Bafrooei, A.; Fiorin, G.; Klein, M. L.;

Borguet, E., Structure of Water at Charged Interfaces: a Molecular Dynamics Study. Langmuir 2014, 30, 8056-8065. (18)

Tian, C. S.; Shen, Y. R., Recent Progress on Sum-Frequency Spectroscopy. Surf. Sci.

Rep. 2014, 69, 105-131. (19)

Eisenthal, K. B., Liquid Interfaces Probed by Second-Harmonic and Sum-Frequency

Spectroscopy. Chem. Rev. 1996, 96, 1343-1360. (20)

Yamaguchi, S.; Tahara, T., Interface-Specific χ(4) Coherent Raman Spectroscopy in the

Frequency Domain. J. Phys. Chem. B 2005, 109, 24211-24214. (21)

Yamaguchi, S.; Tahara, T., χ(4) Raman Spectroscopy for Buried Water Interfaces.

Angew. Chem. Int. Ed. 2007, 46, 7609-7612.

ACS Paragon Plus Environment

19

The Journal of Physical Chemistry Letters 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

(22)

Page 20 of 28

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. (23)

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. (24)

Shen, Y. R., Phase-Sensitive Sum-Frequency Spectroscopy. Annu. Rev. Phys. Chem.

2013, 64, 129-150. (25)

Myalitsin, A.; Urashima, S.; Nihonyanagi, S.; Yamaguchi, S.; Tahara, T., Water

Structure at the Buried Silica/Aqueous Interface Studied by Heterodyne-Detected Vibrational Sum-Frequency Generation. J. Phys. Chem. C 2016, 120, 9357-9363. (26)

Vanselous, H.; Petersen, P. B., Extending the Capabilities of Heterodyne-Detected

Sum-Frequency Generation Spectroscopy: Probing Any Interface in Any Polarization Combination. J. Phys. Chem. C 2016, 120, 8175-8184. (27)

Ge, A.; Peng, Q.; Qiao, L.; Yepuri, N. R.; Darwish, T. A.; Matsusaki, M.; Akashi, M.;

Ye, S., Molecular Orientation of Organic Thin Films on Dielectric Solid Substrates: a Phase-Sensitive Vibrational SFG Study. Phys. Chem. Chem. Phys. 2015, 17, 18072-18078. (28)

Wen, Y.-C.; Zha, S.; Liu, X.; Yang, S.; Guo, P.; Shi, G.; Fang, H.; Shen, Y. R.; Tian, C.,

Unveiling Microscopic Structures of Charged Water Interfaces by Surface-Specific Vibrational Spectroscopy. Phys. Rev. Lett. 2016, 116, 016101.

ACS Paragon Plus Environment

20

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

The Journal of Physical Chemistry Letters

(29)

Ohno, P. E.; Wang, H. F.; Geiger, F. M., Second-Order Spectral Lineshapes from

Charged Interfaces. Nat. Commun. 2017, 8, 1032. (30)

Bard, A. J.; Faulkner, L. R., Electrochemical Methods: Fundamentals and Applications,

2nd ed. JOHN WILEY & SONS, INC.: 2001. (31)

Joutsuka, T.; Hirano, T.; Sprik, M.; Morita, A., Effects of Third-Order Susceptibility in

Sum Frequency Generation Spectra: a Molecular Dynamics Study in Liquid Water. Phys. Chem. Chem. Phys. 2018, 20, 3040-3053. (32)

Pezzotti, S.; Galimberti, D. R.; Shen, Y. R.; Gaigeot, M. P., Structural Definition of the

BIL and DL: a New Universal Methodology to Rationalize Non-Linear χ(2)(ω) SFG Signals at Charged Interfaces, Including χ(3)(ω) Contributions. Phys. Chem. Chem. Phys. 2018, 20, 5190-5199. (33)

Kasuya, M.; Hino, M.; Yamada, H.; Mizukami, M.; Mori, H.; Kajita, S.; Ohmori, T.;

Suzuki, A.; Kurihara, K., Characterization of Water Confined Between Silica Surfaces Using the Resonance Shear Measurement. J. Phys. Chem. C 2013, 117, 13540-13546. (34)

Dalstein, L.; Potapova, E.; Tyrode, E., The Elusive Silica/Water Interface: Isolated

Silanols Under Water as Revealed by Vibrational Sum Frequency Spectroscopy. Phys. Chem. Chem. Phys. 2017, 19, 10343-10349. (35)

Ostroverkhov, V.; Waychunas, G. A.; Shen, Y. R., New Information on Water

Interfacial Structure Revealed by Phase-Sensitive Surface Spectroscopy. Phys. Rev. Lett. 2005, 94, 046102.

ACS Paragon Plus Environment

21

The Journal of Physical Chemistry Letters 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

(36)

Page 22 of 28

Lovering, K. A.; Bertram, A. K.; Chou, K. C., New Information on the

Ion-Identity-Dependent Structure of Stern Layer Revealed by Sum Frequency Generation Vibrational Spectroscopy. J. Phys. Chem. C 2016, 120, 18099-18104. (37)

Schaefer, J.; Gonella, G.; Bonn, M.; Backus, E. H. G., Surface-Specific Vibrational

Spectroscopy of the Water/Silica Interface: Screening and Interference. Phys. Chem. Chem. Phys. 2017, 19, 16875-16880. (38)

Covert, P. A.; Jena, K. C.; Hore, D. K., Throwing Salt into the Mix: Altering Interfacial

Water Structure by Electrolyte Addition. J. Phys. Chem. Lett. 2014, 5, 143-148. (39)

Nihonyanagi, S.; Yamaguchi, S.; Tahara, T., Water Hydrogen Bond Structure near

Highly Charged Interfaces Is Not like Ice. J. Am. Chem. Soc. 2010, 132, 6867-6869. (40)

Kusaka, R.; Ishiyama, T.; Nihonyanagi, S.; Morita, A.; Tahara, T., Structure at the

Air/Water Interface in the Presence of Phenol: a Study Using Heterodyne-Detected Vibrational Sum Frequency Generation and Molecular Dynamics Simulation. Phys. Chem. Chem. Phys. 2017, 20, 3002-3009. (41)

Sovago, M.; Kramer Campen, R.; Bakker, H. J.; Bonn, M., Hydrogen Bonding Strength

of Interfacial Water Determined with Surface Sum-Frequency Generation. Chem. Phys. Lett. 2009, 470, 7-12. (42)

Rull, F.; de Saja, J. A., Effect of Electrolyte Concentration on the Raman Spectra of

Water in Aqueous Solutions. J. Raman. Spectrosc. 1986, 17, 167-172. (43)

Wyss, H. R.; Falk, M., Infrared Spectrum of HDO in Water and in NaCl Solution. Can.

J. Chem. 1970, 48, 607-614.

ACS Paragon Plus Environment

22

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

The Journal of Physical Chemistry Letters

(44)

Singh, P. C.; Nihonyanagi, S.; Yamaguchi, S.; Tahara, T., Interfacial Water in the

Vicinity of a Positively Charged Interface Studied by Steady-State and Time-Resolved Heterodyne-Detected Vibrational Sum Frequency Generation. J. Chem. Phys. 2014, 141, 18C527. (45)

Kilic, M. S.; Bazant, M. Z.; Ajdari, A., Steric Effects in the Dynamics of Electrolytes at

Large Applied Voltages. I. Double-Layer Charging. Phys. Rev. E 2007, 75, 021502. (46)

Eftekhari-Bafrooei, A.; Borguet, E., Effect of Electric Fields on the Ultrafast Vibrational

Relaxation of Water at a Charged Solid–Liquid Interface as Probed by Vibrational Sum Frequency Generation. J. Phys. Chem. Lett. 2011, 2, 1353-1358. (47)

Yang, Z.; Li, Q.; Chou, K. C., Structures of Water Molecules at the Interfaces of

Aqueous Salt Solutions and Silica: Cation Effects. J. Phys. Chem. C 2009, 113, 8201-8205.

ACS Paragon Plus Environment

23

The Journal of Physical Chemistry Letters 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

284x284mm (96 x 96 DPI)

ACS Paragon Plus Environment

Page 24 of 28

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

The Journal of Physical Chemistry Letters

215x317mm (96 x 96 DPI)

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 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

98x151mm (190 x 190 DPI)

ACS Paragon Plus Environment

Page 26 of 28

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

The Journal of Physical Chemistry Letters

104x76mm (190 x 190 DPI)

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 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

262x218mm (96 x 96 DPI)

ACS Paragon Plus Environment

Page 28 of 28