The interfacial structure of water droplets in a hydrophobic liquid

Abstract

Nanoscopic and microscopic water droplets and ice crystals embedded in liquid hydrophobic surroundings are key components of aerosols, rocks, oil fields and the human body. The chemical properties of such droplets critically depend on the interfacial structure of the water droplet. Here we report the surface structure of 200 nm-sized water droplets in mixtures of hydrophobic oils and surfactants as obtained from vibrational sum frequency scattering measurements. The interface of a water droplet shows significantly stronger hydrogen bonds than the air/water or hexane/water interface and previously reported planar liquid hydrophobic/water interfaces at room temperature. The observed spectral difference is similar to that of a planar air/water surface at a temperature that is ∼50 K lower. Supercooling the droplets to 263 K does not change the surface structure. Below the homogeneous ice nucleation temperature, a single vibrational mode is present with a similar mean hydrogen-bond strength as for a planar ice/air interface.

ARTICLE

Received 19 Aug 2016 | Accepted 5 Apr 2017 | Published 24 May 2017

DOI: 10.1038/ncomms15548 OPEN

The interfacial structure of water droplets in a hydrophobic liquid

Nikolay Smolentsev1, Wilbert J. Smit2, Huib J. Bakker2 & Sylvie Roke1

Nanoscopic and microscopic water droplets and ice crystals embedded in liquid hydrophobic surroundings are key components of aerosols, rocks, oil elds and the human body. The chemical properties of such droplets critically depend on the interfacial structure of the water droplet. Here we report the surface structure of 200 nm-sized water droplets in mixtures of hydrophobic oils and surfactants as obtained from vibrational sum frequency scattering measurements. The interface of a water droplet shows signicantly stronger hydrogen bonds than the air/water or hexane/water interface and previously reported planar liquid hydro-phobic/water interfaces at room temperature. The observed spectral difference is similar to that of a planar air/water surface at a temperature that is B50 K lower. Supercooling the droplets to 263 K does not change the surface structure. Below the homogeneous ice nucleation temperature, a single vibrational mode is present with a similar mean hydrogen-bond strength as for a planar ice/air interface.

1 Laboratory for Fundamental BioPhotonics (LBP), Institute of Bioengineering (IBI), and Institute of Materials Science (IMX), School of Engineering (STI), and Lausanne Centre for Ultrafast Science (LACUS),cole Polytechnique Fdrale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland. 2 AMOLF, Science Park 104, 1098 XG Amsterdam, The Netherlands. Correspondence and requests for materials should be addressed to S.R. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 8:15548 | DOI: 10.1038/ncomms15548 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15548

Aqueous hydrophobic interfaces play a dening role on nano- and microscopic length scales. Pockets and droplets of water in a hydrophobic environment are omnipresent

in the atmosphere (as ice particles and cloud droplets1), the earth2,3 (in oil elds and inside internal pores of many geological materials), and in chemical4,5 and biological processes68 (as vehicles for medicine delivery9). For all of these processes, the nanoscopic water droplet/hydrophobic interface plays a critical role in determining the fate of the system. The reactivity of conned water is crucial for atmospheric science1, earth science3 and biology10. The molecular surface properties of water droplets are commonly inferred from data obtained from the macroscopic air/water interface1116 or from aqueous solutions of solvated hydrophobes17, systems that are very different in both their size, chemical composition and temperature dependence from actual droplet surfaces. As such, the surface structure and reactivity of water droplets and conned water remains elusive, despite its relevance.

Here we present temperature-controlled sum frequency scattering (SFS) experiments on nanoscopic water droplets in a liquid hydrophobic environment. We nd that compared to planar water/air or water/hexane interfaces the spectrum reveals a greater order corresponding to that of a planar water surface at a B50 or 40 K lower temperature. This increase in order is explained by the formation of an extended hydrophobic network, which does not exist on planar interfaces and in solution. In addition, on supercooling the droplets, the surface spectrum does not change shape. Cooling the water droplets below the homogeneous ice nucleation temperature results in a spectrum with a single symmetric peak that is similar to the one found for the basal ice/air interface.

ResultsRoom temperature droplet experiments. Nanoscopic (100 nm radius) water droplets are prepared via sonication in d34-hexadecane (Fig. 1) or in a 1:1 mixture of decane and cyclohexane (Fig. 2), both with 5 mM of a hydrophobic surfactant molecule that partially covers the interface (Span80, Fig. 1a). The waterhexadecane partitioning coefcient for these hydro-phobic liquids is 5.5 10 4 (Span80 (ref. 18)), 41.6 10 6

(decane) and 1.6 10 4 (cyclohexane)19, providing a liquid

hydrophobic/water interface in all systems. The interfacial structure of the water droplets and ice nanocrystals is measured with vibrational SFS20,21, a method that combines interface-specic sum frequency generation (SFG)11,12,14,15,22,23 with light scattering24 (illustrated in Fig. 1b). The scattered vibrational spectra report on the average orientational order and H-bond network structure of water molecules in the rst few monolayers of water at the interfacial region21.

Figure 1c shows the scattered SF spectrum of D2O nanodroplets in oil (296 K) and reection SFG spectra of a planar liquid D2O/air and a planar D2O/hexane interface. The spectra were normalized taking into account linear absorption effects as well as discontinuities in the normal component of the electromagnetic elds as they cross the droplet interface (see Methods section). The SFG spectra of the two planar interfaces of Fig. 1c have three main features: A peak at 2,370/2,400 cm 1, a peak at 2,500 cm 1 and a peak at 2,745/2,725 cm 1. The frequency of the OD stretching vibrations is correlated to the strength of the interfacial hydrogen (H)-bond interaction: the OD frequency decreases with increasing H-bond strength. Hence, the feature at 2,500 cm 1 reports on water molecules that are more weakly H-bonded than the ones that attribute to the feature at 2,370/2,400 cm 1. The peak at 2,745/2,725 cm 1 originates from interfacial OD bonds that are not H-bonded25. The 296 K droplet/hydrophobic

interface spectrum contains the same 2,370 and 2,500 cm 1 features (indicated by the dotted lines), but the higher frequency components are signicantly less intense than they are in the planar interface spectra. This indicates that there are relatively more strongly H-bonded water molecules at the water droplet interface than at the air/water and hexane/water interface recorded at the same temperature. The SFG spectrum of the Span80/water interface is similar to that of D2O/hexane and

D2O/air and shown in the Supplementary Fig. 1.

The peak ratio of the 2,370 and 2,500 cm 1 modes is temperature dependent. This temperature-dependent ratio is plotted in Fig. 1d for the air/water and hexane/water interfaces. Lowering the temperature increases26 the ratio, indicating that the population of stronger H-bonds increases over that of weaker H-bonds. The peak ratio obtained for the water droplets is also indicated (black marker). Extrapolating the ratios found for the air/water and hexane/water spectra (blue and red lines) the same ratio as for the water droplets would be found at a hypothetical planar water/air or water/hexane interface of B245/260 K (that is, B50/40 K colder). Thus, the water molecules of the water droplet/hydrophobic interface appear to be much more structured than the water molecules at a planar water/air or water/hexane interface. Note that this conclusion was arrived at by comparing amplitudes at particular frequencies. The spectrum of water is known to be highly complex, involving strong inhomogeneous broadening27,28, excitonic coupling effects29 and a Fermi resonance with the overtone of the bending vibration30. As such, aside from the feature that high-frequency modes generally refer to more weakly H-bonded water molecules than low-frequency modes, there are many ways to interpret the spectral information15,16,25,26,28,3135. Our analysis addresses the temperature dependence of the ratio of weakly and strongly H-bonded water molecules, and is as such quite model independent.

Very similar SFS spectra are obtained for other combinations of hydrophobic oils (mixtures of decane and cyclohexane, uorinated oil) and surfactants (Tween80, tri-block PEG 900, Supplementary Fig. 2), which indicates that the observed red-shift of the spectrum is a common phenomenon for water droplets in various hydrophobic liquids. The absence of unbonded OD groups in the droplet spectrum can be partially explained by the presence of OH groups on the Span80 molecule that can form H-bonds with interfacial water molecules that would otherwise be unbonded (see Supplementary Note 1, Supplementary Fig. 3 and Supplementary Table 1). As a result, the peak at 2,745 cm 1, corresponding to OD groups that are not H-bonded, vanishes.

From the detection limit of the SFS system36, it follows that the imbalance between OD groups that point in or out of the water droplet surface is less than one free OD group per 27 nm2. The difference with the droplet spectrum probably arises from a difference in structure, which we will discuss in more detail further on. Note that impurities are not expected to be a problem here36 as the sample volume is small (B50 ml), and the surface to volume ratio is large (B105 cm 1; which is at least three orders of magnitude larger than in standard planar reection SFG measurements).

Supercooling the droplets. Next, we repeat the experiment at reduced temperatures, supercooling the droplets (Fig. 2a) and ultimately freezing them by cooling below the homogeneous nucleation point of D2O (237 K (ref. 37), Fig. 2b). Figure 2a shows

SFS spectra recorded at different temperatures (296 K, red, 273 K, green and 263 K, blue) using a mixture of decane and cyclohexane as the main phase. The droplet spectrum from Fig. 1c is also shown (black line). It can be seen that despite the different

2 NATURE COMMUNICATIONS | 8:15548 | DOI: 10.1038/ncomms15548 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15548 ARTICLE

a Span80

Hexane

Hexadecane

c

IR frequency for H2O (cm1)

O

3,100

3,200

3,300

3,400

3,500

3,600

3,700

O

D2O droplets/oil

O

Normalized SFS/SFG intensity

2.0

1.5

1.0

0.5

OH

OH

OH

b

Sum frequency scattering and molecular structure

VIS

SFS

2,300

2,400

2,500

2,600

2,700

2,800

IR frequency for D2O (cm1)

d

Peak ratio

250

300

350

Temperature (K)

Figure 1 | Surface structure of water droplets. (a) Chemical structure of hexadecane, hexane and Span80. (b) Illustration of a vibrational SFS experiment. A femtosecond infrared (IR) and a picosecond visible (VIS) laser beam overlap in space and time in the sample and generate scattered sum frequency light from the water droplets. The inset illustrates the surface hydrogen-bond structure of a water droplet embedded in a hydrophobic liquid. In this gure, we sketch a possible molecular arrangement of the water molecules at the surface of a droplet surrounding an alkane tail that protrudes into the water phase (looking from the inside of the droplet). Please note that the surfactant is not present in the illustration. (c) SFS spectra of D2O droplets in d34-hexadecane with 5 mM Span80 (black) and reection SFG spectra of the planar D2O/air interface (blue) and planar D2O/hexane interface (red). The SFS/SFG spectra are collected with horizontally (P) polarized infrared and vertically polarized (S) VIS and SF beams. The normalization procedure takes into account the infrared pulse shape in the sample as well as discontinuities of the electromagnetic elds at the interface (see Methods section). The top axis shows the corresponding frequencies axis for H2O ( 1.3515). (d) The ratio between the low- and high-frequency bands of the SFG spectrum of the planar water/air

and water/hexane interface26 (blue and red markers) and an extrapolation with a quadratic polynomial t to a lower temperature range (blue and red lines). The ratio for the room temperature water nanodroplet spectrum is shown as a black marker.

hydrophobic phase the water spectrum does not show a signicant difference. In addition, reducing the liquid temperature below 273 K does not result in any spectral changes. Figure 2b shows the SF spectrum of D2O nanocrystals (233 K) and the ice/

air interface for different temperatures (170, 200 and 230 K) in the frequency range around the strongly H-bonded mode at B2,350 (3,170) cm 1. The resonance in the SF spectrum of the ice nanocrystals has a similar frequency as is observed for the planar ice surface of the same temperature. This indicates that the difference in surface structure observed for the liquid interface is not present in the solid phase. The absence of the effect for the frozen droplets can be understood from the fact that when the droplets freeze, the mutual waterwater interactions become relatively more important for the resulting water SFG spectrum (and the detected H-bond strength) than the interactions with the hydrophobic oil molecules. For crystalline planar ice, the observed SF ice/air spectrum originates in part from OH groups connecting the different layers in the ice crystal38. These deeper layers of water molecules are much less inuenced by the hydrophobic groups at the interface. Hence, for planar ice and frozen water droplets of the same temperature the H-bonds of the probed OH groups are of similar strength.

DiscussionThermodynamic39,40 and Raman multivariate curve resolution (MCR) measurements on solvated alkanols17 show water molecules hydrating a hydrophobic solute are more tetrahedrally ordered and thus more constrained in congurational space than water molecules in bulk liquid water (as sketched in the inset of Fig. 1b). These solutes are very small and full the criterion of Chandler4 that states that hydrophobic solutes/surfaces lead to an enhancement of the tetrahedral ordering of water if the radius of the solute/surface is o1 nm (ref. 4). In this perspective, the present observations of an enhanced ordering of water molecules at the surface of water nanodroplets embedded in oil is perhaps surprising, as the droplet radius (100 nm) is much larger than 1 nm. However, we can consider the oil surrounding the water droplet surface as an extended curved array of hydrophobic perturbations, that is, hydrophobic groups protruding into the surface of the water droplet. The water molecules near the hydrophobic protrusions will fold their H-bond network around these (o1 nm)

protrusions. We speculate that a high density of hydrophobic protrusions in combination with a modest overall curvature of the water surface (corresponding to a diameter of a few hundred

NATURE COMMUNICATIONS | 8:15548 | DOI: 10.1038/ncomms15548 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15548

a

IR frequency for H2O (cm1)

3,000

3,250

3,500

3,750

Normalized SFS intensity

b

Normalized SFS/SFG intensity

2,200

2,300

2,400

2,500

2,600

2,700

2,800

IR frequency for D2O (cm1)

Figure 2 | Supercooled water droplets and ice nanocrystals. (a) SFS spectra of D2O droplets at room temperature (296 K, in a mixture of decane and cyclohexane (red) and in d34-hexadecane (black)) and under supercooled conditions (green, 273 K and blue, 263 K, both in a mixture of decane and cyclohexane). (b) SFS spectra of D2O droplets at 233 K, below the homogeneous ice nucleation temperature (black) and SFG spectra of the basal

H2O/air interface at 170 K (blue), 200 K (green) and 230 K (red). All droplets samples contain 5 mM Span80. Spectral data from H2O and D2O are interchangeable. Frequency conversion between H2O and D2O was obtained by multiplying the D2O frequencies by 1.3515. The SFG spectra are collected with horizontally (P) polarized infrared and vertically (S) polarized VIS and SF beams.

nanometres), leads to a strong structuring of the water hydrogen-bond network, resulting in an extended highly structured corrugated droplet surface. For smaller droplets (for example, within micelles/solutes), the surface is too strongly curved to display such a favourable interplay between short-range and long-range interactions. For larger droplets, and at surfaces, the curvature is too weak to get such a favourable interplay, leading to local disruptions of the H-bond network and an enhanced inhomogeneity. This picture of the enhanced ordering of curved water surfaces will require support from high-quality molecular dynamics simulations.

In summary, vibrational SFS spectra from 100 nm radius water droplets embedded in a hydrophobic liquid at different temperatures show that the H-bond network of the interfacial water possesses a much enhanced (in-plane) tetrahedral structure compared to a planar air/water interface. This greater order manifests itself as a red-shifted spectrum that has a 100% increased peak ratio of the 2,370 and 2,500 cm 1 modes.

This increase is much larger than observed in previous studies of planar liquid hydrophobic/water interfaces. The observed

structure of the water droplet interface corresponds to that of an air/water interface that is B50 K colder. This increase in order is explained by the formation of an extended network of hydrophobic protrusions into the water droplet surface, which does not exist on planar interfaces and in solution. On supercooling the droplets, the surface spectrum does not change shape. Cooling the sample below the homogeneous ice nucleation temperature results in a spectrum with a single symmetric peak at a comparable frequency as that of a planar basal ice/air interface of the same temperature.

The presented experiments demonstrate the possibility of quantifying the interfacial structure of water droplets and illustrate the effect of nanoscopic hydrophobic surfaces on the H-bond structure of water. With a reduction of the effective surface temperature we expect that the reactivity of water droplets embedded in a hydrophobic environment41, such as in cloud droplets, atmospheric ice particles1, rocks3 or medicine carriers9 is lower than what one might think based on the actual temperature. Future studies on such nanoscopic systems may reveal important information on the role and reactivity of interfacial water in aerosol formation, protein folding, pore functioning4,6,7,42 and the charging and stabilization of hydrophobic interfaces.

Methods

Chemicals. Before use all glassware was cleaned with a 1:3 H2O2:H2SO4 solution, after which it was thoroughly rinsed with ultra-pure water (H2O, Milli-Q UF plus, Millipore, Inc., electrical resistance of 18.2 MO cm). Hexadecane (Fluka, 99.8%), decane (Fluka, 99.8%), cyclohexane (Sigma, 499.7%), Span80 (Sigma, GC quality), Tween80 (Sigma, GC quality), sodium dodecyl sulfate, SDS (99% BioMol), d34-hexadecane (C16D34, 98% d, Cambridge Isotope Laboratories), uorinatedoil (Novec HFE7500, 3-Ethoxy-dodecauor-2-triuormethyl-hexan), sulfuricacid (9597%, ISO, Merck) and H2O2 (30%, Reactolab SA) were used as received. The purity of alkanes was veried with a Zisman test43,44. All samples for SFS measurements were prepared using D2O (99.8% Armar, 42 MO cm).

The D2O for the water/air SFG experiment was obtained from Cambridge Isotope Laboratories (99.9%).

Water nanodroplets. Water nanodroplets were prepared using a sonication procedure with 1 v.v.% D2O in oil concentration and 5 mM Span80 (or

Span80:Tween80 mixture), according to the procedure described in ref. 45. For some samples (shown in Supplementary Fig. 2) 10 mM SDS was added to the water phase before the emulsication. The droplet-size distribution was measured with dynamic light scattering (DLS, Zetasizer Nano ZS, Malvern). For all samples, the droplets have a mean diameter of B200 nm with a polydispersity index of B0.2.

Infrared spectra were recorded with a Bruker Vertex 70 instrument equipped with silver mirrors.

Monocrystalline ice. Monocrystalline ice was grown following the manual of Basu et al.46. Five pieces of PVC pipe with small notches at the bottom were placed in an aluminium pan, which was chilled from underneath with a cooling liquid. Millipore water (18.2 MO cm) was degassed and added until the bottoms of the moulds were covered. The temperature in the pan was set to 0.5 C and after

equilibration ice nucleation was initiated by adding a small piece of ice to the centre of the pan. The small notches at the bottom of each PVC mould allow only one ice crystal to propagate through, resulting in a single ice crystal in each mould. For 3 subsequent days, 1 cm water was added to each mould and the temperature was lowered by 0.3 C. At day 4, the ice was pushed outside of the pipe pieces and stored in a freezer, where all further handling was done. The crystal orientation was determined using a Rigsby-type universal stage47 and cut to the basal plane with a band saw. The ice surface was polished by a microtome and a 4 mm thick slice was cut off for the experiment.

Vibrational sum frequency reection mode spectra. The SFG and VIS beams were polarized in the vertical (S) direction, and the infrared beam was polarized in the horizontal plane (P), leading to the SSP (SFG, VIS, IR) polarization combination using the experimental set-up described in ref. 26. The set-up was purged with N2 gas. The VIS beam with a pulse energy of 14 mJ was centred at 798.6 nm with a full-width at half-maximum (FWHM) bandwidth of 13 cm 1. The D2O/air measurements were performed with the 10 mJ 300 cm 1 FWHM infrared beam centred at a wavelength of 4 mm. The procedure for measuring hexane/water is given in ref. 26. To measure the water/Span80 interface, 0.16 mM Span80 was dissolved in 0.1 ml hexane and added to liquid D2O. Subsequently, the hexane was

4 NATURE COMMUNICATIONS | 8:15548 | DOI: 10.1038/ncomms15548 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15548 ARTICLE

evaporated off. The diameter of the teon tray for this experiment was 3.5 cm. The basal ice/air measurements were performed with a 3 mJ 550 cm 1 FWHM IR beam centred at 3 mm. The angles of incidence with respect to the surface normal were 35 (VIS) and 40 (infrared). The D2O/air and basal ice spectra were recorded with acquisition times of 30 s and 300 s, respectively. The spectral intensities were baseline subtracted and normalized with the infrared and VIS pulse energies, and the (normalized) SFG spectrum recorded from z-cut quartz.

The ice measurements were performed using a temperature cell that was cooled with liquid N2. It allows transmittance of the laser beams through a CaF2 window (see Supplementary Fig. 4). The temperature was monitored by a thermocouple welded on the edge of the ice surface with a drop of water and the desired temperature was set by a heating foil resistance, covered by a copper plate.

Vibrational SFS spectra. Vibrational SFS spectra were recorded using the set-up for SFG experiments described in refs 4850. An 800 nm regeneratively amplied Ti:sapphire system (Spitre Pro, Spectra physics) seeded with an 80 MHz 800 nm oscillator (Integral 50, Femtolasers) was operated at a 1 kHz repetition rate to pump a commercial OPG/OPA/DPG system (HE-TOPAS-C, Light Conversion), which was used to generate infrared pulses. The visible beam was split off directly from the amplier, and spectrally shaped with a home-built pulse shaper. The angle between the 10 mJ visible (VIS) beam (800 nm, FWHM 15 cm 1) and the 6 mJ infrared beam (34.5 mm, FWHM 160 cm 1) was 20 (as measured in air). The focused laser beams were overlapped in a sample cuvette with a path length of 200 mm. At a scattering angle (y, measured in air) of 55, the scattered sum frequency (SF) light was collimated using a plano-convex lens (f 15 mm, Thorlabs LA1540-B)

and passed through two short-wave pass lters (3rd Millenium, 3RD770SP). The SF light was spectrally dispersed with a monochromator (Acton, SpectraPro 2300i) and detected with an intensied CCD camera (Princeton Instruments, PI-Max3) using a gate width of 10 ns. The acquisition time for a single spectrum was 300 s. A Glan-Taylor prism (Thorlabs, GT15-B), a half-wave plate (EKSMA, 460-4215) and a polarizing beam splitter cube (CVI, PBS-800-050) and two BaF2 wire grid polarizers (Thorlabs, WP25H-B) were used to control the polarization of the SFG, VIS and infrared beams, respectively. All measurements were performed in the SSP polarization combination. For normalization purposes, a reection beam path is accessible by removable optics so that scattering and reection can be done without a realignment procedure.

Temperature-controlled SFS measurements. Temperature-controlled SFS measurements were performed using a custom-made sample cell (Quantum Northwest, Supplementary Fig. 5). The sample cuvette was placed in a metallic holder cooled with a single Peltier element. The sample chamber was closed and lled with N2 gas to avoid condensation of air. The windows of the temperature cell were made from CaF2 (for the incoming beams) and quartz (for the scattered SF beam). The temperature cell itself was thermally isolated and also placed inside a box purged with N2 gas. The hot side of the Peltier element was cooled with a ow of cold liquid ethanol (223 K, using a peristaltic pump, Ismatec ISM1200). The cold ethanol was obtained by cooling ethanol inside a copper coil that was placed in a bath containing a mixture of dry ice and ethanol. The infrared beam path outside the laser source was purged with N2 gas.

Generation of broadband infrared pulses. To probe the whole frequency range of the D2O stretch mode region for the SFS spectral recordings, we stepped the infrared frequency trough the spectral window between 2,000 and 3,400 cm 1 using 50 nm steps (leading to a total acquisition time of 5,400 s per spectrum).

Supplementary Fig. 6a shows SF spectra for the different infrared centre frequencies recorded in reection mode from a z-cut quartz crystal. The summed spectrum for one of the runs is shown in Supplementary Fig. 6b.

Spectral shape and infrared absorption. The scattered SF spectra recorded from the water droplets were baseline subtracted and summed to obtain the spectral intensity for the entire spectral region of the OD stretching mode. The resulting spectrum was then divided by the energies of the infrared and VIS beams and normalized to the summed SFG spectrum of a z-cut quartz crystal (shown for example in Supplementary Fig. 6b). Since SF scattering experiments are performed in transmission mode, we have to take into account that the oil phase is absorbing part of the infrared spectrum, leading to a transmission spectrum mFTIR. To correct

for absorption effects, we divide the obtained SFS spectrum by a Fourier transform infrared transmittance spectrum measured from the same sample in the same cuvette. Supplementary Fig. 7 shows the results of this correction as the blue spectra. It can be seen that this procedure results in signicant spectral changes in the red spectral side since here the SFG spectrum overlaps with the CD modes of the deuterated main phase (d34-hexadecane). The CH modes in the SFS spectrum originate from the Span80 molecules that are present on the oil side of the interface (Span80 is insoluble in water18).

Discontinuity in the interfacial electric elds. Frequency-dependent refraction modies the incident electromagnetic eld at the interface, which can lead to signicant frequency-dependent distortions in the obtained water spectra51,52. The

sum frequency reection mode spectra are therefore divided by the Fresnel factors14, whereby we use the refractive indices of D2O and ice from refs 5355.

The effective refractive index of the interfacial layers is calculated using a simple slab model56. For linear light scattering experiments such effects are captured by linear Mie theory57 and thus for SFS experiments by nonlinear Mie theory58. Using nonlinear Mie theory, we calculated a frequency-dependent correction factor for the SFS spectra using the input parameters from Supplementary Table 2. The refractive index values for oil and D2O are taken from refs 59,60. Supplementary

Fig. 8 shows the effects of the correction on the spectral data reported in Figs 1 and 2, where it can be seen that the inuence of refraction effects is minimal. The reason for the small difference is that the droplets are small compared to the wavelength. For simplicity, in the main text we refer to the corrected SFS intensity ISFS/IIRmFTIR as normalized SFS intensity.

Data availability. The data that support the ndings of this study are available from the corresponding author on reasonable request.

References

1. de Leeuw, G. et al. Production ux of sea spray aerosol. Rev. Geophys. 49, 139 (2011).

2. Putnis, A. Why mineral interfaces matter. Science 343, 14411442 (2014).3. Zachara, J. et al. Internal domains of natural porous media revealed: critical locations for transport, storage, and chemical reaction. Environ. Sci. Technol. 50, 28112829 (2016).

4. Chandler, D. Interfaces and the driving force of hydrophobic assembly. Nature 437, 640647 (2005).

5. Zhang, L., Chen, K., Chen, B., White, J. L. & Resasco, D. E. Factors that determine zeolite stability in hot liquid water. J. Am. Chem. Soc. 137, 1181011819 (2015).

6. Dyson, H. J., Wright, P. E. & Scheraga, H. A. The role of hydrophobic interactions in initiation and propagation of protein folding. Proc. Natl Acad. Sci. 103, 1305713061 (2006).

7. Ball, P. Water as an active constituent in cell biology. Chem. Rev. 108, 74108 (2008).

8. Mhlmann, D. T. F. Are nanometric lms of liquid undercooled interfacial water bio-relevant? Cryobiology 58, 256261 (2009).

9. Amstad, E., Kim, S.-H. & Weitz, D. A. Photo- and thermoresponsive polymersomes for triggered release. Angew. Chem. Int. Ed. 51, 1249912503 (2012).

10. Grzybowski, B. A. & Huck, W. T. S. The nanotechnology of life-inspired systems. Nat. Nanotechnol. 11, 585592 (2016).

11. Johnson, C. M. & Baldelli, S. Vibrational sum frequency spectroscopy studies of the inuence of solutes and phospholipids at vapor/water interfaces relevant to biological and environmental systems. Chem. Rev. 114, 84168446 (2014).

12. McFearin, C. L., Beaman, D. K., Moore, F. G. & Richmond, G. L. From Franklin to today: toward a molecular level understanding of bonding and adsorption at the oilwater interface. J. Phys. Chem. C 113, 11711188 (2009).

13. Gan, W., Wu, D., Zhang, Z., Feng, R.-R. & Wang, H.-F. Polarization and experimental conguration analyses of sum frequency generation vibrational spectra, structure, and orientational motion of the air/water interface. J. Chem. Phys. 124, 114705 (2006).

14. Wang, H.-F., Gan, W., Lu, R., Rao, Y. & Wu, B.-H. Quantitative spectral and orientational analysis in surface sum frequency generation vibrational spectroscopy (SFG-VS). Int. Rev. Phys. Chem. 24, 191256 (2005).

15. Bonn, M., Nagata, Y. & Backus, E. H. G. Molecular structure and dynamics of water at the waterair interface studied with surface-specic vibrational spectroscopy. Angew. Chem. Int. Ed. 54, 55605576 (2015).

16. Shen, Y. R. & Ostroverkhov, V. Sum-frequency vibrational spectroscopy on water interfaces: polar orientation of water molecules at interfaces. Chem. Rev. 106, 11401154 (2006).

17. Davis, J. G., Gierszal, K. P., Wang, P. & Ben-Amotz, D. Water structural transformation at molecular hydrophobic interfaces. Nature 491, 582585 (2012).

18. Peltonen, L. J. & Yliruusi, J. Surface pressure, hysteresis, interfacial tension, and CMC of four sorbitan monoesters at waterair, waterhexane, and hexaneair interfaces. J. Colloid Interface Sci. 227, 16 (2000).

19. Abraham, M. H., Whiting, G. S., Fuchs, R. & Chambers, E. J. Thermodynamics of solute transfer from water to hexadecane. J. Chem. Soc. Perkin Trans. 2, 291300 (1990).

20. Roke, S. et al. Vibrational sum frequency scattering from a submicron suspension. Phys. Rev. Lett. 91, 258302 (2003).

21. Roke, S. & Gonella, G. Nonlinear light scattering and spectroscopy of particles and droplets in liquids. Annu. Rev. Phys. Chem. 63, 353378 (2012).

22. Hunt, J. H., Guyot-Sionnest, P. & Shen, Y. R. Observation of C-H stretch vibrations of monolayers of molecules optical sum-frequency generation. Chem. Phys. Lett. 133, 189192 (1987).

NATURE COMMUNICATIONS | 8:15548 | DOI: 10.1038/ncomms15548 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15548

23. Harris, A. L., Chidsey, C. E. D., Levinos, N. J. & Loiacono, D. N. Monolayer vibrational spectroscopy by infrared-visible sum generation at metal and semiconductor surfaces. Chem. Phys. Lett. 141, 350356 (1987).

24. Bohren, C. F. & Huffman, D. R. Absorption and Scattering of Light by Small Particles (Wiley, 2007).

25. Du, Q., Superne, R., Freysz, E. & Shen, Y. R. Vibrational spectroscopy of water at the vapor/water interface. Phys. Rev. Lett. 70, 23132316 (1993).

26. Strazdaite, S., Versluis, J., Backus, E. H. G. & Bakker, H. J. Enhanced ordering of water at hydrophobic surfaces. J. Chem. Phys. 140, 054711 (2014).

27. Auer, B. M. & Skinner, J. L. Vibrational sum-frequency spectroscopy of the water liquid/vapor interface. J. Phys. Chem. B 113, 41254130 (2009).

28. Ni, Y. & Skinner, J. L. Communication: vibrational sum-frequency spectrum of the air-water interface, revisited. J. Chem. Phys. 145, 031103 (2016).

29. Buch, V. et al. Sum frequency generation surface spectra of ice, water, and acid solution investigated by an exciton model. J. Chem. Phys. 127, 204710 (2007).

30. Sovago, M. et al. Vibrational response of hydrogen-bonded interfacial water is dominated by intramolecular coupling. Phys. Rev. Lett. 100, 173901 (2008).

31. Scatena, L. F., Brown, M. G. & Richmond, G. L. Water at hydrophobic surfaces: weak hydrogen bonding and strong orientation effects. Science 292, 908912 (2001).

32. Stiopkin, I. V. et al. Hydrogen bonding at the water surface revealed by isotopic dilution spectroscopy. Nature 474, 192195 (2011).

33. Shultz, M. J., Schnitzer, C., Simonelli, D. & Baldelli, S. Sum frequency generation spectroscopy of the aqueous interface: ionic and soluble molecular solutions. Int. Rev. Phys. Chem. 19, 123153 (2000).

34. Nihonyanagi, S. et al. Accurate determination of complex w(2) spectrum of the air/water interface. J. Chem. Phys. 143, 124707 (2015).

35. Zhang, Z., Piatkowski, L., Bakker, H. J. & Bonn, M. Ultrafast vibrational energy transfer at the water/air interface revealed by two-dimensional surface vibrational spectroscopy. Nat. Chem. 3, 888893 (2011).

36. Jena, K. C., Scheu, R. & Roke, S. Surface impurities are not responsible for the charge on the oil/water interface: a comment. Angew. Chem. Int. Ed. 51, 1293812940 (2012).

37. Stckel, P., Weidinger, I. M., Baumgartel, H. & Leisner, T. Rates of homogeneous ice nucleation in levitated H2O and D2O droplets. J. Phys. Chem. A 109, 25402546 (2005).

38. Li Barnett, I., Groenzin, H. & Shultz, M. J. Hydrogen bonding in the hexagonal ice surface. J. Phys. Chem. A 115, 60396045 (2011).

39. Frank, H. S. & Evans, M. W. Free volume and entropy in condensed systems iii. Entropy in binary liquid mixtures; partial molal entropy in dilute solutions; structure and thermodynamics in aqueous electrolytes. J. Chem. Phys. 13, 507532 (1945).

40. Glew, D. N. Aqueous solubility and the gas-hydrates. The methane-water system. J. Phys. Chem. 66, 605609 (1962).

41. Jung, Y. & Marcus, R. A. On the theory of organic catalysis on water.J. Am. Chem. Soc. 129, 54925502 (2007).42. Aryal, P., Abd-Wahab, F., Bucci, G., Sansom, M. S. P. & Tucker, S. J. A hydrophobic barrier deep within the inner pore of the TWIK-1 K2P potassium channel. Nat. Commun. 5, 4377 (2014).

43. Day, J. P. R. & Bain, C. D. Ellipsometric study of depletion at oil-water interfaces. Phys. Rev. E 76, 041601 (2007).

44. Bigelow, W. C., Pickett, D. L. & Zisman, W. A. Oleophobic monolayers.J. Colloid Sci. 1, 513538 (1946).45. Chen, Y., Jena, K. C., Ltgebaucks, C., Okur, H. I. & Roke, S. Three dimensional nano Langmuir trough for lipid studies. Nano Lett. 15, 55585563 (2015).

46. Basu, K. et al. Determining the ice-binding planes of antifreeze proteins by uorescence-based ice plane afnity. J. Vis. Exp. 83, e51185 (2014).

47. Groenzin, H. Sum-Frequency Studies Of Single Crystalline Ice Ih. PhD thesis, Tufts University (2007).

48. de Aguiar, H. B., Samson, J.-S. & Roke, S. Probing nanoscopic droplet interfaces in aqueous solution with vibrational sum-frequency scattering: a study of the effects of path length, droplet density and pulse energy. Chem. Phys. Lett. 512, 7680 (2011).

49. de Aguiar, H. B., Scheu, R., Jena, K. C., de Beer, A. G. F. & Roke, S. Comparison of scattering and reection SFG: a question of phase-matching. Phys. Chem. Chem. Phys. 14, 68266832 (2012).

50. Scheu, R., Chen, Y., Subinya, M. & Roke, S. Stern layer formation induced by hydrophobic interactions: a molecular level study. J. Am. Chem. Soc. 135, 1933019335 (2013).

51. Pool, R. E., Versluis, J., Backus, E. H. G. & Bonn, M. Comparative study of direct and phase-specic vibrational sum-frequency generation spectroscopy: advantages and limitations. J. Phys. Chem. B 115, 1536215369 (2011).

52. Feng, R.-R., Guo, Y., L, R., Velarde, L. & Wang, H.-F. Consistency in the sum frequency generation intensity and phase vibrational spectra of the air/neat water interface. J. Phys. Chem. A 115, 60156027 (2011).

53. Warren, S. G. & Brandt, R. E. Optical constants of ice from the ultraviolet to the microwave: a revised compilation. J. Geophys. Res. 113, D14220 (2008).

54. Max, J.-J. & Chapados, C. Isotope effects in liquid water by infrared spectroscopy. III. H2O and D2O spectra from 6000 to 0cm 1. J. Chem. Phys.

131, 184505 (2009).55. Kedenburg, S., Vieweg, M., Gissibl, T. & Giessen, H. Linear refractive index and absorption measurements of nonlinear optical liquids in the visible and near-infrared spectral region. Opt. Mater. Express 2, 15881611 (2012).

56. Zhuang, X., Miranda, P. B., Kim, D. & Shen, Y. R. Mapping molecular orientation and conformation at interfaces by surface nonlinear optics. Phys. Rev. B 59, 1263212640 (1999).

57. Jackson, J. D. Classical Electrodynamics (Wiley, 1999).58. de Beer, A. G. F. & Roke, S. Nonlinear Mie theory for second-harmonic and sum-frequency scattering. Phys. Rev. B 79, 155420 (2009).

59. Bertie, J. E., Ahmed, M. K. & Eysel, H. H. Infrared intensities of liquids. 5. Optical and dielectric constants, integrated intensities, and dipole moment derivatives of water and water-d2 at 22 degree C. J. Phys. Chem. 93, 22102218 (1989).

60. Tuntomo, A., Tien, C. L. & Park, S. H. Optical constants of liquid hydrocarbon fuels. Combust. Sci. Technol. 84, 133140 (1992).

Acknowledgements

This work is supported by the Julia Jacobi Foundation, and the European Research Council (Grant number 616305). We thank Prof. Esther Amstad for providing the sample of water droplets in uorinated oil. Dion Ursem and Hinco Schoenmaker are acknowledged for technical support. Part of this work belongs to the research programme of the Netherlands Organisation for Scientic Research (NWO) and was performed at the research institute AMOLF.

Author contributions

N.S. and W.J.S. performed experiments, N.S., W.J.S., H.J.B. and S.R. interpreted the data and wrote the manuscript. S.R. conceived and supervised the work.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/ http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Smolentsev, N. et al. The interfacial structure of water droplets in a hydrophobic liquid. Nat. Commun. 8, 15548 doi: 10.1038/ncomms15548 (2017).

Publishers note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional afliations.

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecommons.org/licenses/by/4.0/

r The Author(s) 2017

6 NATURE COMMUNICATIONS | 8:15548 | DOI: 10.1038/ncomms15548 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications