Probing the Structure and Dynamics of Interfacial Water with Scanning Tunneling Microscopy and Spectroscopy

Summary

Here, we present a protocol to investigate the structure and dynamics of interfacial water at the atomic scale, in terms of submolecular resolution imaging, molecular manipulation, and single-bond vibrational spectroscopy.

Abstract

Water/solid interfaces are ubiquitous and play a key role in many environmental, biophysical, and technological processes. Resolving the internal structure and probing the hydrogen-bond (H-bond) dynamics of the water molecules adsorbed on solid surfaces are fundamental issues of water science, which remains a great challenge owing to the light mass and small size of hydrogen. Scanning tunneling microscopy (STM) is a promising tool for attacking these problems, thanks to its capabilities of sub-Ångström spatial resolution, single-bond vibrational sensitivity, and atomic/molecular manipulation. The designed experimental system consists of a Cl-terminated tip and a sample fabricated by dosing water molecules in situ onto the Au(111)-supported NaCl(001) surfaces. The insulating NaCl films electronically decouple the water from the metal substrates, so the intrinsic frontier orbitals of water molecules are preserved. The Cl-tip facilitates the manipulation of the single water molecules, as well as gating the orbitals of water to the proximity of Fermi level (E_F) via tip-water coupling. This paper outlines the detailed methods of submolecular resolution imaging, molecular/atomic manipulation, and single-bond vibrational spectroscopy of interfacial water. These studies open up a new route for investigating the H-bonded systems at the atomic scale.

Introduction

The interactions of water with the surfaces of solid materials are involved in various surface reaction processes, such as heterogeneous catalysis, photoconversion, electrochemistry, corrosion and lubrication et al.^1,2,3 In general, to investigate interfacial water, spectroscopic and diffraction techniques are commonly used, such as infrared and Raman spectroscopy, sum-frequency generation (SFG), X-ray diffraction (XRD), nuclear magnetic resonance (NMR), neutron scattering^4,5,6,7,8. However, these methods suffer from the limitation of spatial resolution, spectral broadening, and averaging effects.

STM is a promising technique to overcome these limitations, which combines the sub-Ångström spatial resolution, atomic manipulation, and single-bond vibrational sensitivity^{9,10,11,12,13,14}. Since the beginning of this century, STM has been extensively applied to investigate the structure and dynamics of water on solid surfaces^{3,15,16,17,18,19,20}. Additionally, vibrational spectroscopy based on STM could be obtained from the second-derivative differential tunneling conductance (d²I/dV²), also known as inelastic electron tunneling spectroscopy (IETS). However, resolving the internal structure, i.e. the H-bond directionality, and acquiring reliable vibrational spectroscopy of water are still challenging. The main difficulty lies in that water is a close shell molecule, whose frontier orbitals are far away from the E_F, thus the electrons from the STM tip can hardly tunnel into the molecular resonance states of water, leading to the poor signal-to-noise ratio of molecular imaging and vibrational spectroscopy.

Water adsorbed on the Au-supported NaCl(001) films provides an ideal system for atomic-scale investigation by STM with a Cl-terminated tip (Figure 1a), which is performed at 5 K in the ultrahigh-vacuum (UHV) environment with a base pressure better than 8×10^-11 mbar. On one hand, the insulating NaCl films decouple water molecules electronically from the Au substrate so the native frontier orbitals of water are preserved and the lifetime of the electrons residing in the molecular resonant state is prolonged. On the other hand, the STM tip could effectively tune the frontier orbital of water toward the E_F via tip-water coupling, especially when the tip is functionalized with a Cl atom. These key steps enable high-resolution orbital imaging and vibrational spectroscopy of water monomers and clusters. In addition, water molecules could be manipulated in a well-controlled manner, due to the strong electrostatic interaction between the negatively charged Cl-tip and water.

In this report, the preparation procedures of the sample and the Cl-terminated tip for STM investigation are outlined in detail in section 1 and 2, respectively. In section 3, we describe the orbital imaging technique, by which the O-H directionality of water monomer and tetramer are resolved. The tip-enhanced IETS is introduced in section 4, which allows the detection of vibrational modes of water molecules at single-bond limit, and determination of the H-bonding strength with high accuracy from the red shift in the oxygen-hydrogen stretching frequency of water. In section 5, we show how the water tetramer can be constructed and switched by controlled tip manipulation. Based on the orbital imaging, spectroscopy, and manipulation techniques, isotopic substitution experiments can be performed to probe the quantum nature of protons in interfacial water, such as quantum tunneling and zero-point motion.

Protocol

NOTE: The experiments are performed on water molecules adsorbed on the Au-supported NaCl(001) film (Figure 1a) at 5 K with an ultrahigh-vacuum (UHV) cryogenic STM equipped with Nanonis electronic controller.

1. Fabrication of Experimental Sample

1. Clean the Au(111) single crystal
   1. Pump the gas line to the pressure of ~10^-7 mbar and then flush the gas line with Ar gas. Put through the pump/flush cycle for three times.
      NOTE: Each pump/flush cycle takes about 30 min.
   2. Fill the gas line with Ar gas to the pressure of 2 bar, thus prohibiting the atmosphere to permeate through the gas line.
   3. Put the Au(111) crystal on the heater stage, which is mounted on the manipulator in the UHV chamber (base pressure of 1.4×10^-10 mbar).
   4. Clean the Au(111) single crystal by cycles of Ar⁺ ion sputtering for 15 min (p(Ar) = 5×10⁻⁵ mbar, 1.0 kV, 6 µA) and subsequent annealing at about 900 K for 5 min.
      NOTE: The annealing temperature should be decreased slowly, otherwise a high density of step edges will form on the Au surface. 3 - 5 sputtering/annealing cycles are usually used.
   5. Transfer the Au(111) sample to the STM scanning stage, and check the cleanliness with STM (Inset of Figure 1b).
2. Deposition of NaCl on the Au(111) substrate
   1. Degas the NaCl source. Slowly increase the current applied on the Knudsen cell until the temperature of the source reaches 670 K. Degas the NaCl source several times until the pressure of the chamber is below 4×10^-9 mbar.
      NOTE: The increasing rate of the current depends on the outgassing rate of the NaCl source to maintain the pressure of the chamber below 1×10^-8 mbar.
   2. Put the Au(111) sample on the manipulator and adjust the position of the Au sample to make the sample face the shutter of the Knudsen cell.
      NOTE: The temperature of the Au(111) substrate could be decreased below room temperature (77-300 K) by cooling down the manipulator head with a continuous flow of liquid nitrogen
   3. Increase the current applied on the Knudsen cell until the temperature of the source reaches 640 K, and let the evaporation flux stabilize for 5 min before opening the shutter.
   4. Open the shutter and deposit the NaCl onto the Au(111) sample held at 290 K for 2 min.
   5. Transfer the Au-supported NaCl sample to the STM scanning stage. Check the coverage and size of the bilayer NaCl(001) islands on the Au(111) substrate with STM (Figure 1b).
3. Purify the water under vacuum by freeze-pump-thaw cycles²¹ to remove remaining impurities.
   1. Prepare three sealed-off glass-UHV adapters. Place water H₂O, D₂O, and D₂O:HOD:H₂O isotopic mixture solutions (2 mL) into three adapters separately, and mount the adapters on the gas line (Figure 2).
      NOTE: The D₂O:HOD:H₂O isotopic mixtures can be obtained by mixing the ultrapure H₂O and D₂O with equal amounts under ultrasonic oscillation for 10 min.
   2. Freeze the liquid water with liquid nitrogen. Make sure that the gas line is pumped to the pressure of ~10^-7 mbar before freezing.
   3. Open the diaphragm-sealed valve and pump off the atmosphere for 15 min. Then close the diaphragm-sealed valve and thaw the solution.
      NOTE: Gas bubbles evolve from the solution when it is thawing.
      CAUTION: Let the frozen solvent thaw by itself. Thawing the solution with water bath may cause the glass vessel to break. To freeze and thaw the solution quickly, replace glass-UHV adapters with metal-UHV adapters, though the solution in the metal vessel is invisible.
   4. Repeat steps 1.3.2-1.3.3 until no gas bubbles evolve from the solution as the solution thaws. Put through the freeze-pump-thaw cycle at least three times.
   5. Close the bellows-sealed valve and leave the gas line in vacuum. Then open the diaphragm-sealed valve, and let the water vapor fill in the gas line.
4. Dose water molecules in situ onto the sample surface
   1. Decrease the temperature of the sample to 5 K. Open the leak valve slowly to make the pressure of the STM UHV chamber increase to 2×10^-10 mbar.
      NOTE: Water molecules flow into the UHV chamber through the dosing tube, which points to the shutter of the shield. The distance between the shutter and the sample (in the shield) is about 6 cm. The base pressure of STM chamber are better than 7×10^-11 mbar. The deposition rate is about 0.01 bilayer min^-1.
   2. Open the shutter. Dose the water molecules onto the Au-supported NaCl surface for 1 min. Then close the shutter and the leak valve.
   3. Check the coverage of water molecules on the Au-supported NaCl(001) surface with STM. Isolated water monomers form on the sample surface (Figure 1c).

2. Preparation of the Cl-Terminated Tip

1. Fabricate an electrochemically etched tungsten (W) tip.
   1. Place 0.3 mm W wire into a 3 Mol/L NaOH etching solution with an immersion length of about 2 mm.
   2. Apply a 5 V dc potential to the W wire with respect to a platinum ring electrode inserted into the NaOH solution.
   3. Stop the etching process when the suspended W wire fell off. Clean the etched W tip with distilled water and ethanol. Then transfer the W tip into the scanner.
      NOTE: The electrochemically etched W tip can be used for one year before exchange.
2. Apply voltage pulses (2-10 V) and controlled crashing procedures (0.25-0.4 nm) on the STM tip until the atomic Cl atoms of the NaCl surface are resolved.
   NOTE: The STM tip is poked on a clean region of the Au(111) surface.
3. Position the STM tip over the center of one Cl atom (Figure 3a). Bring the bare STM tip close to the NaCl surface in proximity with the set point of V = 5 mV and I = 5 n A (Figure 3b).
4. Retract the tip to original set point (Figure 3c) and scan the same area. Check the obtainment of the Cl-tip by improved resolution and a missing Cl atom in the STM image of NaCl (Figure 3d-e).
   NOTE: Unsuccessful cases may occasionally occur, for example, when the Cl atom doesn't transfer to the STM or multiple Cl/Na atoms adsorb on the tip. If this occurs, repeat the steps 2.2-2.5.

3. Orbital Imaging of Water Monomer

1. Shape the STM tip with voltage pulse (2-10 V) and controlled crashing (0.25-0.4 nm) procedures.
2. Scan the water molecules adsorbed on the NaCl(001) surface with 10 nm by 10 nm frame at 5 K.
3. Focus on one individual water monomer and zoom in. Scan the water monomer in a systematic way as functions of the bias (-400-400 mV) and the tunneling current (50-300 pA).
   NOTE: With a bare STM tip, the highest occupied (HOMO) and lowest unoccupied (LUMO) molecular orbitals of water appear at positive and negative bias, respectively²². Once the tip is Cl-terminated, only HOMO emerges (Figure 4a), and the LUMO feature is not observed throughout the accessible bias range (from -400 mV to 400 mV). Even under larger bias voltage, the water molecules will be unstable due to vibrational excitation.

4. Single-Molecule Vibrational Spectroscopy

1. Setup of the digital lock-in and bias spectroscopy module (Nanonis electronic controller)
   1. Setup of bias spectroscopy module: Select the current, LIX1 (dI/dV spectra signal), and LIX2 (d²I/dV² spectra signal) channels. Set the setting time as 50 ms, and integration time as 300 ms. Increase the integration time and sweep times to obtain smooth spectra. Tune the Z offset to take the bias spectroscopy at different tip heights. Make sure the Z-controller set holds and lock-ins runs during measurement
      NOTE: Setting time is defined as: the time to wait after changing the bias to the next level and before starting to acquire data to avoid transient effect induced by the bias change. Integration time is defined as: the time during which the data are acquired and averaged.
   2. Setup of lock-in module
      NOTE: The scanning tunneling spectroscopy, dI/dV and d²I/dV² spectra, are acquired simultaneously using a lock-in amplifier by demodulating the first and second harmonics of the tunneling current, respectively.
      1. Modulate the bias and demodulate the current. Set the modulation frequency as a few hundred Hz and modulation amplitude as 5-7 mV. Make sure there is no mechanical and electronic noise at the set point frequency and the corresponding second harmonic frequency.
      2. Set the first harmonic phase: Switch to the Z-controller module. Set the tip lift to 10 nm and turn off the feedback. Switch to the lock-in module and turn on the lock-in button (green). Click the first harmonic auto phase and record the phase. Repeat the auto phase at least five times and take the average. Then subtract 90 degrees from the averaged phase to get the phase of the junction.
      3. Set the second harmonic phase: Position the STM tip on the Au(111) substrate and start the bias spectroscopy sweep from -1.5 V to 1.5 V. Select the channel LIX 1 and function dY/dX, which together show the derivative of dI/dV spectrum. Find a prominent peak feature in the spectrum and set the corresponding energy as the bias. Turn on the lock-in and keep the STM system in tunneling. Click the second harmonic auto phase at least five times and take the average.
         NOTE: Since the second harmonic signals are usually very weak, the phase can fluctuate wildly. When decreasing the tip height to increase the intensity of the signal, the fluctuation of the phase will be much smaller (a few degrees) and the second harmonic phase will be more accurate.
2. Tip-enhanced IETS of a D₂O monomer
   1. Scan a water monomer with Cl-tip at the set point of V=100 mV and I=50 pA.
   2. Position the Cl-tip on the NaCl surface and take the bias spectroscopy as the background signal. Then position the Cl-tip on the water monomer and start the bias spectroscopy sweep.
   3. If the dI/dV and d²I/dV² spectra of water are featureless, simply follow the background NaCl surface (blue curves of Figure 4c-d). Decrease the tip height by tuning the Z offset until the vibrational features emerge in the spectra (red curves of Figure 4c-d).
      ​NOTE: For the IETS measurement, a long integration time (~1s) and multiple sweeps are needed. For a D₂O water monomer, set the bias range from -360 mV to 360 mV. For H₂O/HOD water monomers, sweep the bias from -475 mV to 475 mV. Comparing with D₂O, H₂O, and HOD, water monomers are more easily disturbed and even swept away during IETS measurement.
3. H-bond strength
   1. Repeat steps 4.2.2-4.2.3 and tune the sweep bias range to focus on the stretching mode of water monomers. The IETS of water D₂O, H₂O, and HOD are presented and discussed in ²³.
   2. Obtain the H-bonding energy by converting from the redshift of the H-bonded OH stretching frequency (relative to the free OH stretching energy) using this empirical formula:
      ΔH = 1.3 × √Δv    (1)
      NOTE: ΔH is the H-bonding energy, in kJ/mol; Δv is the redshift of the OH stretching mode, in cm^-1. Convert the unit of H-bonding strength to meV by: 1kJ/mol=10.4 meV/atom. To apply Eq. 1 to the OD stretching mode, the quantity Δv should be multiplied by a factor: v(OH) / v(OD) = 1.3612, where v(OH) and v(OD) are the OH and OD stretching frequencies of the free HOD molecule, respectively.

5. Molecular Manipulation

1. Construction of a water tetramer (Figure 5a)
   1. Scan an area containing four water monomers. Position the Cl-tip on top of a monomer at the set point of V=100 mV and I=50 pA. Decrease the height to the set point of V=10 mV and I=150 pA to enhance the tip-water interaction.
   2. Move the Cl-tip along the predesigned trajectories. Then retract the tip to the initial set point (V=100 mV, I=50 pA), and rescan the same area to check that the water dimer is formed.
   3. Repeat the steps 5.1.1-5.1.2 to form the water trimer and tetramer.
      NOTE: The above manipulation process could be realized by the Nanonis controller (Scan control-Follow me module). Setup of the Scan control-Follow me module:
      Bias: 10 mV
      Speed: 500 pm/s
      Z-ctrl Setpoint: 150 pA
      Switch on/off Z-Ctrl: green
      Time to wait: 1s
      Current Gain: LN 10^9
      Path: Click RECORD button and draw up the designed trajectories on the image, then click the STOP button.
      Click the EXECUTE button and the STM tip will move along the predesigned trajectories with the setpoint in the Follow me module. If the water monomer doesn't move, decrease the tip height (smaller bias and bigger current) during manipulation.
2. Chirality switching of a water tetramer (Figure 6)
   1. Scan a water tetramer with the Cl-tip. Change the set point to V=5 mV, I=5 pA, and position the tip slightly off the center of the water tetramer.
   2. In the Z-controller module, define a distance to lift the tip when the Z-controller is switched off (e.g., Tip Lift: -230 pm). Turn off the Z-controller feedback. Bring the tip close to the water tetramer (~230 pm).
   3. Record the current trace, which shows two different levels, indicating that tetramer has undergone reversible interconversion between two H-bonding chirality.
   4. Leave the current at the high level and switch on the Z-controller feedback. Retract the tip to the original set point (V=5 mV, I=5 pA). Then scan the water tetramer with the set point of V=10 mV and I=100 pA to check the chiral state of the water tetramer.
   5. Repeat steps 5.2.1-5.2.4 at least 10 times to confirm the corresponding chiral state of water tetramer at the high current level.
   6. Repeat steps 5.2.1- 5.2.4 at least 10 times, but leave the current at the low level to check the corresponding chiral state of water tetramer.
   7. Record the tunneling trace for 20 min, which contains a few hundred switching events
   8. Plot the distribution of the times the tetramer spent in the low and high level of the current trace, respectively.
   9. Fit the distribution to an exponential decay (Figure 7). Then get the fitted time constant. Inverse the time constant to yield the switching rate.

Results

Figure 1a illustrates the schematic of the STM experimental setup. First, Au(111) substrate is cleaned by sputtering and annealing cycles in the UHV chamber. The clean Au(111) sample shows 22×√3 reconstructed surface, where the atoms of the surface layer occupy both the hcp and the fcc sites forming herringbone structures (Inset of Figure 1b). The NaCl is evaporated on the Au(111) subs...

Discussion

To probe the internal structure, dynamics, and vibrational spectroscopy of water molecules adsorbed on the solid surfaces, paying particular attention to the degrees of freedom of hydrogen, some experimental steps are of crucial importance, which will be discussed in the following paragraphs.

The orbital imaging of water molecules is achieved based on two key steps. First, the insulating NaCl films decouple the water electronically from the Au substrate, second the orbital gating effect of the...

Materials

Name	Company	Catalog Number	Comments
Au(111) single crystal	MaTeck	NA
NaCl	Sigma Aldrich	450006
Water, deuterium-depleted	Sigma Aldrich	195294
Deuterium oxide	Sigma Aldrich	364312
Sealed-off glass-UHV adapters	MDC vacuum products	46300
Diaphragm-sealed valve	any	NA
Bellows-sealed valve	any	NA
Leak valve	Kurt J. Lesker	NA
Scanning tunneling microscopy	CreaTec	NA
Electronic controller.	Nanonis	NA
Tungsten wire	any		diameter:0.3 mm; purity: 99.95%