Comparisons of 2D IR measured spectral diffusion in rotating frames using pulse shaping and in the stationary frame using the standard method.

Multidimensional visible spectroscopy using pulse shaping to produce pulses with stable controllable phases and delays has emerged as an elegant tool to acquire electronic spectra faster and with greatly reduced instrumental and data processing errors. Recent migration of this approach using acousto-optic modulator (AOM) pulse shaping to the mid-infrared region has proved useful for acquiring two dimensional infrared (2D IR) vibrational echo spectra. The measurement of spectral diffusion in 2D IR experiments hinges on obtaining accurate 2D line shapes. To date, pulse shaping 2D IR has not been used to study the time-dependent spectral diffusion of a vibrational chromophore. Here we compare the spectral diffusion data obtained from a standard non-collinear 2D IR spectrometer using delay lines to the data obtained from an AOM pulse shaper based 2D IR spectrometer. The pulse shaping experiments are performed in stationary, partially rotating, and fully rotating reference frames and are the first in the infrared to produce 2D spectra collected in a fully rotating frame using a phase controlled pulse sequence. Rotating frame experiments provide a dramatic reduction in the number of time points that must be measured to obtain a 2D IR spectrum, with the fully rotating frame giving the greatest reduction. Experiments were conducted on the transition metal carbonyl complex tricarbonylchloro(1,10-phenanthroline)rhenium(I) in chloroform. The time dependent data obtained from the different techniques and with different reference frames are shown to be in agreement.

Investigation of nanostructure in room temperature ionic liquids using electronic excitation transfer.

Donor-donor electronic excitation transfer among 9-phenylanthracene (9PA) chromophores was measured in the room temperature ionic liquid (RTIL) 1-methyl-3-octylimidazolium chloride using time dependent fluorescence depolarization. 9PA, which is uncharged and nonpolar, will partition into the organic regions of the RTIL. The excitation transfer rate, which is sensitive to the distribution of chromophores in the RTIL, is modeled using different spatial configurations of 9PA molecules in the RTIL solution. The models are an isotropic distribution (random distribution) and a clustered sphere model to represent hydrophobic regions of a nanostructured environment. Model calculations were performed to demonstrate the sensitivity of excitation transfer to different distributions of chromophores. When compared to the experiment, the isotropic model can adequately match the data. From a Bayesian analysis of the sensitivity of the excitation transfer to the models for the spatial distribution of chromophores, an upper limit of 6 Å radius is placed on the size of hydrophobic domains in the RTIL.

Structural dynamics of a catalytic monolayer probed by ultrafast 2D IR vibrational echoes.

Ultrafast two-dimensional infrared (2D IR) vibrational echo spectroscopy has proven broadly useful for studying molecular dynamics in solutions. Here, we extend the technique to probing the interfacial dynamics and structure of a silica surface-tethered transition metal carbonyl complex--tricarbonyl (1,10-phenanthroline)rhenium chloride--of interest as a photoreduction catalyst. We interpret the data using a theoretical framework devised to separate the roles of structural evolution and excitation transfer in inducing spectral diffusion. The structural dynamics, as reported on by a carbonyl stretch vibration of the surface-bound complex, have a characteristic time of ~150 picoseconds in the absence of solvent, decrease in duration by a factor of three upon addition of chloroform, and decrease another order of magnitude for the bulk solution. Conversely, solvent-complex interactions increase the lifetime of the probed vibration by 160% when solvent is applied to the monolayer.

Dynamics of water at the interface in reverse micelles: measurements of spectral diffusion with two-dimensional infrared vibrational echoes.

Water dynamics inside of reverse micelles made from the surfactant Aerosol-OT (AOT) were investigated by observing spectral diffusion, orientational relaxation, and population relaxation using two-dimensional infrared (2D IR) vibrational echo spectroscopy and pump-probe experiments. The water pool sizes of the reverse micelles studied ranged in size from 5.8 to 1.7 nm in diameter. It is found that spectral diffusion, characterized by the frequency-frequency correlation function (FFCF), significantly changes as the water pool size decreases. For the larger reverse micelles (diameter 4.6 nm and larger), the 2D IR signal is composed of two spectral components: a signal from bulk-like core water, and a signal from water at the headgroup interface. Each of these signals (core water and interfacial water) is associated with a distinct FFCF. The FFCF of the interfacial water layer can be obtained using a modified center line slope (CLS) method that has been recently developed. The interfacial FFCFs for large reverse micelles have a single exponential decay (∼1.6 ps) to an offset plus a fast homogeneous component and are nearly identical for all large sizes. The observed ∼1.6 ps interfacial decay component is approximately the same as that found for bulk water and may reflect hydrogen bond rearrangement of bulk-like water molecules hydrogen bonded to the interfacial water molecules. The long time offset arises from dynamics that are too slow to be measured on the accessible experimental time scale. The influence of the chemical nature of the interface on spectral diffusion was explored by comparing data for water inside reverse micelles (5.8 nm water pool diameter) made from the surfactants AOT and Igepal CO-520. AOT has charged, sulfonate head groups, while Igepal CO-520 has neutral, hydroxyl head groups. It is found that spectral diffusion on the observable time scales is not overly sensitive to the chemical makeup of the interface. An intermediate-sized AOT reverse micelle (water pool diameter of 3.3 nm) is analyzed as a large reverse micelle because it has distinct core and interface regions, but its core region is more constrained than bulk water. The interfacial FFCF for this intermediate-sized reverse micelle is somewhat slower than those found for the larger reverse micelles. The water nanopools in the smaller reverse micelles cannot be separated into core and interface regions. In the small reverse micelles, the FFCFs are biexponential decays to an offset plus a fast homogeneous component. Each small reverse micelle exhibits an ∼1 ps decay time, which may arise from local hydrogen bond fluctuations and a slower, ∼6-10 ps decay, which is possibly due to slow hydrogen bond rearrangement of noninterfacial water molecules or topography fluctuations at the interface.

Extracting 2D IR frequency-frequency correlation functions from two component systems.

The center line slope (CLS) method is often used to extract the frequency-frequency correlation function (FFCF) from 2D IR spectra to delineate dynamics and to identify homogeneous and inhomogeneous contributions to the absorption line shape of a system. While the CLS method is extremely efficient, quite accurate, and immune to many experimental artifacts, it has only been developed and properly applied to systems that have a single vibrational band, or to systems of two species that have spectrally resolved absorption bands. In many cases, the constituent spectra of multiple component systems overlap and cannot be distinguished from each other. This situation creates ambiguity when analyzing 2D IR spectra because dynamics for different species cannot be separated. Here a mathematical formulation is presented that extends the CLS method for a system consisting of two components (chemically distinct uncoupled oscillators). In a single component system, the CLS corresponds to the time-dependent portion of the normalized FFCF. This is not the case for a two component system, as a much more complicated expression arises. The CLS method yields a series of peak locations originating from slices taken through the 2D spectra. The slope through these peak locations yields the CLS value for the 2D spectra at a given T(w). We derive analytically that for two component systems, the peak location of the system can be decomposed into a weighted combination of the peak locations of the constituent spectra. The weighting depends upon the fractional contribution of each species at each wavelength and also on the vibrational lifetimes of both components. It is found that an unknown FFCF for one species can be determined as long as the peak locations (referred to as center line data) of one of the components are known, as well as the vibrational lifetimes, absorption spectra, and other spectral information for both components. This situation can arise when a second species is introduced into a well characterized single species system. An example is a system in which water exists in bulk form and also as water interacting with an interface. An algorithm is presented for back-calculating the unknown FFCF of the second component. The accuracy of the algorithm is tested with a variety of model cases in which all components are initially known. The algorithm successfully reproduces the FFCF for the second component within a reasonable degree of error.

Temperature dependent equilibrium native to unfolded protein dynamics and properties observed with IR absorption and 2D IR vibrational echo experiments.

Dynamic and structural properties of carbonmonoxy (CO)-coordinated cytochrome c(552) from Hydrogenobacter thermophilus (Ht-M61A) at different temperatures under thermal equilibrium conditions were studied with infrared absorption spectroscopy and ultrafast two-dimensional infrared (2D IR) vibrational echo experiments using the heme-bound CO as the vibrational probe. Depending on the temperature, the stretching mode of CO shows two distinct bands corresponding to the native and unfolded proteins. As the temperature is increased from low temperature, a new absorption band for the unfolded protein grows in and the native band decreases in amplitude. Both the temperature-dependent circular dichroism and the IR absorption area ratio R(A)(T), defined as the ratio of the area under the unfolded band to the sum of the areas of the native and unfolded bands, suggest a two-state transition from the native to the unfolded protein. However, it is found that the absorption spectrum of the unfolded protein increases its inhomogeneous line width and the center frequency shifts as the temperature is increased. The changes in line width and center frequency demonstrate that the unfolding does not follow simple two-state behavior. The temperature-dependent 2D IR vibrational echo experiments show that the fast dynamics of the native protein are virtually temperature independent. In contrast, the fast dynamics of the unfolded protein are slower than those of the native protein, and the unfolded protein fast dynamics and at least a portion of the slower dynamics of the unfolded protein change significantly, becoming faster as the temperature is raised. The temperature dependence of the absorption spectrum and the changes in dynamics measured with the 2D IR experiments confirm that the unfolded ensemble of conformers continuously changes its nature as unfolding proceeds, in contrast to the native state, which displays a temperature-independent distribution of structures.

Water dynamics in small reverse micelles in two solvents: two-dimensional infrared vibrational echoes with two-dimensional background subtraction.

Water dynamics as reflected by the spectral diffusion of the water hydroxyl stretch were measured in w(0) = 2 (1.7 nm diameter) Aerosol-OT (AOT)/water reverse micelles in carbon tetrachloride and in isooctane solvents using ultrafast 2D IR vibrational echo spectroscopy. Orientational relaxation and population relaxation are observed for w(0) = 2, 4, and 7.5 in both solvents using IR pump-probe measurements. It is found that the pump-probe observables are sensitive to w(0), but not to the solvent. However, initial analysis of the vibrational echo data from the water nanopool in the reverse micelles in the isooctane solvent seems to yield different dynamics than the CCl(4) system in spite of the fact that the spectra, vibrational lifetimes, and orientational relaxation are the same in the two systems. It is found that there are beat patterns in the interferograms with isooctane as the solvent. The beats are observed from a signal generated by the AOT/isooctane system even when there is no water in the system. A beat subtraction data processing procedure does a reasonable job of removing the distortions in the isooctane data, showing that the reverse micelle dynamics are the same within experimental error regardless of whether isooctane or carbon tetrachloride is used as the organic phase. Two time scales are observed in the vibrational echo data, ~1 and ~10 ps. The slower component contains a significant amount of the total inhomogeneous broadening. Physical arguments indicate that there is a much slower component of spectral diffusion that is too slow to observe within the experimental window, which is limited by the OD stretch vibrational lifetime.

Temperature and hydration-dependent rotational and translational dynamics of a polyether oligomer.

Temperature-dependent rotational diffusion of tetraethylene glycol dimethyl ether (TEGDE) is measured by optical heterodyne-detected optical Kerr effect (OHD-OKE) spectroscopy and compared to previous measurements of rotational diffusion as a function of water content. Both types of data, temperature-dependent and hydration-dependent, follow the Debye-Stokes-Einstein (DSE) equation and agree quantitatively with hydrodynamic calculations. Of particular importance is the result that both types of data show nearly identical dependence on the viscosity divided by the temperature (η/T). We also compare the translational diffusion constants as previously measured by pulsed field gradient spin-echo (PFG-SE) NMR as a function of both temperature and water content. The temperature-dependent data follow the Stokes-Einstein (SE) equation. Similar to the rotation, the low water content mixtures obey the SE equation and show the same proportionality to η/T as the temperature-dependent data. At higher water fractions, the data do not obey the SE equation. The principal results are that the influence of temperature on dry TEGDE orientational relaxation is the same as the influence of water content at fixed temperature, and that the influence of temperature on translational diffusion of dry TEGDE is the same as the influence of water content over a range of relatively low water concentrations. The results demonstrate that there are no large TEGDE structural changes or specific, long-lived water-polyether interactions in the solutions over the entire concentration range.

Dynamics of a myoglobin mutant enzyme: 2D IR vibrational echo experiments and simulations.

Myoglobin (Mb) double mutant T67R/S92D displays peroxidase enzymatic activity in contrast to the wild type protein. The CO adduct of T67R/S92D shows two CO absorption bands corresponding to the A(1) and A(3) substates. The equilibrium protein dynamics for the two distinct substates of the Mb double mutant are investigated by using two-dimensional infrared (2D IR) vibrational echo spectroscopy and molecular dynamics (MD) simulations. The time-dependent changes in the 2D IR vibrational echo line shapes for both of the substates are analyzed using the center line slope (CLS) method to obtain the frequency-frequency correlation function (FFCF). The results for the double mutant are compared to those from the wild type Mb. The experimentally determined FFCF is compared to the FFCF obtained from molecular dynamics simulations, thereby testing the capacity of a force field to determine the amplitudes and time scales of protein structural fluctuations on fast time scales. The results provide insights into the nature of the energy landscape around the free energy minimum of the folded protein structure.

Conformational switching between protein substates studied with 2D IR vibrational echo spectroscopy and molecular dynamics simulations.

Myoglobin is an important protein for the study of structure and dynamics. Three conformational substates have been identified for the carbonmonoxy form of myoglobin (MbCO). These are manifested as distinct peaks in the IR absorption spectrum of the CO stretching mode. Ultrafast 2D IR vibrational echo chemical exchange experiments are used to observed switching between two of these substates, A(1) and A(3), on a time scale of <100 ps for two mutants of wild-type Mb. The two mutants are a single mutation of Mb, L29I, and a double mutation, T67R/S92D. Molecular dynamics (MD) simulations are used to model the structural differences between the substates of the two MbCO mutants. The MD simulations are also employed to examine the substate switching in the two mutants as a test of the ability of MD simulations to predict protein dynamics correctly for a system in which there is a well-defined transition over a significant potential barrier between two substates. For one mutant, L29I, the simulations show that translation of the His64 backbone may differentiate the two substates. The simulations accurately reproduce the experimentally observed interconversion time for the L29I mutant. However, MD simulations exploring the same His64 backbone coordinate fail to display substate interconversion for the other mutant, T67R/S92D, thus pointing to the likely complexity of the underlying protein interactions. We anticipate that understanding conformational dynamics in MbCO via ultrafast 2D IR vibrational echo chemical exchange experiments can help to elucidate fast conformational switching processes in other proteins.

Theory of interfacial orientational relaxation spectroscopic observables.

The orientational correlation functions measured in the time-resolved second-harmonic generation (TRSHG) and time-resolved sum-frequency generation (TRSFG) experiments are derived. In the laboratory coordinate system, the Y(l) (m)(Omega(lab)(t))Y(2) (m)(Omega(lab)(0)) (l=1,3 and m=0,2) correlation functions, where the Y(l) (m) are spherical harmonics, describe the orientational relaxation observables of molecules at interfaces. A wobbling-in-a-cone model is used to evaluate the correlation functions. The theory demonstrates that the orientational relaxation diffusion constant is not directly obtained from an experimental decay time in contrast to the situation for a bulk liquid. Model calculations of the correlation functions are presented to demonstrate how the diffusion constant and cone half-angle affect the time-dependence of the signals in TRSHG and TRSFG experiments. Calculations for the TRSHG experiments on Coumarin C314 molecules at air-water and air-water-surfactant interfaces are presented and used to examine the implications of published experimental results for these systems.

Room temperature ionic liquid-lithium salt mixtures: optical Kerr effect dynamical measurements.

The addition of lithium salts to ionic liquids causes an increase in viscosity and a decrease in ionic mobility that hinders their possible application as an alternative solvent in lithium ion batteries. Optically heterodyne-detected optical Kerr effect spectroscopy was used to study the change in dynamics, principally orientational relaxation, caused by the addition of lithium bis(trifluoromethylsulfonyl)imide to the ionic liquid 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide. Over the time scales studied (1 ps-16 ns) for the pure ionic liquid, two temperature-independent power laws were observed: the intermediate power law (1 ps to approximately 1 ns), followed by the von Schweidler power law. The von Schweidler power law is followed by the final complete exponential relaxation, which is highly sensitive to temperature. The lithium salt concentration, however, was found to affect both power laws, and a discontinuity could be found in the trend observed for the intermediate power law when the concentration (mole fraction) of lithium salt is close to chi(LiTf(2)N) = 0.2. A mode coupling theory (MCT) schematic model was also used to fit the data for both the pure ionic liquid and the different salt concentration mixtures. It was found that dynamics in both types of liquids are described very well by MCT.

Hydrogen bond migration between molecular sites observed with ultrafast 2D IR chemical exchange spectroscopy.

Hydrogen-bonded complexes between phenol and phenylacetylene are studied using ultrafast two-dimensional infrared (2D IR) chemical exchange spectroscopy. Phenylacetylene has two possible pi hydrogen bonding acceptor sites (phenyl or acetylene) that compete for hydrogen bond donors in solution at room temperature. The OD stretch frequency of deuterated phenol is sensitive to which acceptor site it is bound. The appearance of off-diagonal peaks between the two vibrational frequencies in the 2D IR spectrum reports on the exchange process between the two competitive hydrogen-bonding sites of phenol-phenylacetylene complexes in the neat phenylacetylene solvent. The chemical exchange process occurs in approximately 5 ps and is assigned to direct hydrogen bond migration along the phenylacetylene molecule. Other nonmigration mechanisms are ruled out by performing 2D IR experiments on phenol dissolved in the phenylacetylene/carbon tetrachloride mixed solvent. The observation of direct hydrogen bond migration can have implications for macromolecular systems.

Dynamics in organic ionic liquids in distinct regions using charged and uncharged orientational relaxation probes.

The temperature-dependent fluorescence anisotropy decay (orientational relaxation) of perylene and sodium 8-methoxypyrene-1,3,6-sulfonate (MPTS) were measured in a series of 1-alkyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide (alkyl = ethyl, butyl, hexyl, octyl) organic room temperature ionic liquids (RTIL). The two fluorescent probe molecules display markedly different rotational dynamics when analyzed using Stokes-Einstein-Debye theory, demonstrating that they are located in distinct environments within the RTILs and have very different interactions with their surroundings. Perylene rotates with subslip behavior, becoming increasingly subslip as the length of ionic liquid alkyl chain is increased. The dynamics approach those of perylene in an organic oil. In contrast, MPTS shows superstick behavior, likely reflecting very strong coordination with the RTIL cations. These results are consistent with different elements of rotational friction within the ionic liquid structure, which are available to solutes depending on their chemical functionality.

Solvent control of the soft angular potential in hydroxyl-pi hydrogen bonds: inertial orientational dynamics.

Ultrafast polarization and wavelength selective IR pump-probe spectroscopy is used to measure the inertial and long time orientational dynamics of pi-hydrogen bonding complexes. Inertial orientational relaxation is sensitive to the angular potential associated with the hydrogen bond. The complexes studied are composed of phenol-OD (hydroxyl hydrogen replaced by deuterium) and various pi-base solvents with different electron donating or withdrawing substituents (chlorobenzene, bromobenzene, benzene, toluene, p-xylene, mesitylene, 1-pentyne). The different substituents provide experimental control of the hydrogen bond strength. The inertial orientational relaxation of the complexes, measured at the center frequency of each line, is independent of the hydrogen bond strength, demonstrating the insensitivity of the OD inertial dynamics, and therefore the H-bond angular potential, to the hydrogen bond strength. OD stretch absorption bands are inhomogeneously broadened through interactions with the solvent. The hydrogen bonding complexes all have similar wavelength dependent inertial orientational relaxation across their inhomogeneously broadened OD stretch absorption lines. The wavelength dependence of the inertial reorientation across each line arises because of a correlation between local solvent structure and the angular potential. These two results imply that local solvent structure acts as the controlling influence in determining the extent of inertial orientational relaxation, and therefore the angular potential, and that variation in the pi-hydrogen bond strength is of secondary importance.

Water dynamics at neutral and ionic interfaces.

The orientational dynamics of water at a neutral surfactant reverse micelle interface are measured with ultrafast infrared spectroscopy of the hydroxyl stretch, and the results are compared to orientational relaxation of water interacting with an ionic interface. The comparison provides insights into the influence of a neutral vs. ionic interface on hydrogen bond dynamics. Measurements are made and analyzed for large nonionic surfactant Igepal CO-520reverse micelles (water nanopool with a 9-nm diameter). The results are compared with those from a previous study of reverse micelles of the same size formed with the ionic surfactant Aerosol-OT (AOT). The results demonstrate that the orientational relaxation times for interfacial water molecules in the two types of reverse micelles are very similar (13 ps for Igepal and 18 ps for AOT) and are significantly slower than that of bulk water (2.6 ps). The comparison of water orientational relaxation at neutral and ionic interfaces shows that the presence of an interface plays the dominant role in determining the hydrogen bond dynamics, whereas the chemical nature of the interface plays a secondary role.

Proton transfer and proton concentrations in protonated Nafion fuel cell membranes.

Proton transfer in protonated Nafion fuel cell membranes is studied using several pyrene derivative photoacids. Proton transfer in the center of the Nafion nanoscopic water channels is probed with the highly charged photoacid 8-hydroxypyrene-1,3,6-trisulfonate (HPTS). At high hydration levels, both the time-integrated fluorescence spectrum and the fluorescence kinetics of HPTS permit the determination of hydronium concentration of the interior of the water pools in Nafion. The proton transfer kinetics of HPTS in protonated Nafion at maximum hydration are identical to the kinetics displayed by HPTS in a 0.5 M HCl solution. The hydronium concentration near the water interface in Nafion is estimated with rhodamine-6G to be 1.4 M. Excited state proton transfer (ESPT) is followed in the nonpolar side chain regions of Nafion with the photoacid 8-hydroxy-N,N,N',N',N'',N''-hexamethylpyrene-1,3,6-trisulfonamide (HPTA). Excited state proton transfer of HPTA is possible in protonated Nafion only at the highest hydration level due to a relatively high local pH.

Water dynamics at the interface in AOT reverse micelles.

The orientational dynamics of water molecules at the interface in large Aerosol-OT (AOT) reverse micelles are investigated using ultrafast infrared spectroscopy of the OD stretch of dilute HOD in H(2)O. In large reverse micelles ( approximately 9 nm diameter or larger), a significant amount of the nanoscopic water is sufficiently distant from the interface that it displays bulk-like characteristics. However, some water molecules interact with the interface and have vibrational absorption spectra and dynamics distinct from bulk water. The different characteristics of these interfacial waters allow their contribution to the data to be separated from the bulk. The infrared absorption spectrum of the OD stretch is analyzed to show that the interfacial water molecules have a spectrum that peaks near 2565 cm(-1) in contrast to 2509 cm(-1) in bulk water. A two-component model is developed that simultaneously describes the population relaxation and orientational dynamics of the OD stretch in the spectral region of the interfacial water. The model provides a consistent description of both observables and demonstrates that water interacting with the interface has slower vibrational relaxation and orientational dynamics. The orientational relaxation of interfacial water molecules occurs in 18 +/- 3 ps, in contrast to the bulk water value of 2.6 ps.

Geometry and nanolength scales versus interface interactions: water dynamics in AOT lamellar structures and reverse micelles.

To determine the relative importance of the confining geometry and nanoscopic length scale versus water/interface interactions, the dynamic interactions between water and interfaces are studied with ultrafast infrared spectroscopy. Aerosol OT (AOT) is a surfactant that can form two-dimensional lamellar structures with known water layer thickness as well as well-defined monodispersed spherical reverse micelles of known water nanopool diameter. Lamellar structures and reverse micelles are compared based on two criteria: surface-to-surface dimensions to study the effect of confining length scales, and water-to-surfactant ratio to study water/interface interactions. We show that the water-to-surfactant ratio is the dominant factor governing the nature of water interacting with an interface, not the characteristic nanoscopic distance. The detailed structure of the interface and the specific interactions between water and the interface also play a critical role in the fraction of water molecules influenced by the surface. A two-component model in which water is separated into bulk-like water in the center of the lamellar structure or reverse micelle and interfacial water is used to quantitatively extract the interfacial dynamics. A greater number of perturbed water molecules are present in the lamellar structures as compared to the reverse micelles due to the larger surface area per AOT molecule and the greater penetration of water molecules past the sulfonate head groups in the lamellar structures.

Ion-water hydrogen-bond switching observed with 2D IR vibrational echo chemical exchange spectroscopy.

The exchange of water hydroxyl hydrogen bonds between anions and water oxygens is observed directly with ultrafast 2D IR vibrational echo chemical exchange spectroscopy (CES). The OD hydroxyl stretch of dilute HOD in H(2)O in concentrated (5.5 M) aqueous solutions of sodium tetrafluoroborate (NaBF(4)) displays a spectrum with a broad water-like band (hydroxyl bound to water oxygen) and a resolved, blue shifted band (hydroxyl bound to BF(4)(-)). At short time (200 fs), the 2D IR vibrational echo spectrum has 4 peaks, 2 on the diagonal and 2 off-diagonal. The 2 diagonal peaks are the 0-1 transitions of the water-like band and the hydroxyl-anion band. Vibrational echo emissions at the 1-2 transition frequencies give rise to 2 off-diagonal peaks. On a picosecond time scale, additional off-diagonal peaks grow in. These new peaks arise from chemical exchange between water hydroxyls bound to anions and hydroxyls bound to water oxygens. The growth of the chemical exchange peaks yields the time dependence of anion-water hydroxyl hydrogen bond switching under thermal equilibrium conditions as T(aw) = 7 +/- 1 ps. Pump-probe measurements of the orientational relaxation rates and vibrational lifetimes are used in the CES data analysis. The pump-probe measurements are shown to have the correct functional form for a system undergoing exchange.