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

Hidden Beneath the Layers: Extending the Core/Shell Model of Reverse Micelles

Reverse micelles (RMs) provide a unique and highly tunable model system to study water in confined environments. The complex properties of water within RMs arise from the disruption of extended hydrogen bond (H-bond) networks that mediate local and long-range dynamics in bulk aqueous systems. Modulating the water pool size influences its H-bond dynamics, with smaller RMs increasingly restricting the H-bond network rearrangements leading to slower dynamics; however, within small confined systems, the dynamics of the surfactants also influence the water dynamics. Using ultrafast two-dimensional infrared spectroscopy, we investigate the effects of RM size on the surfactant headgroup rotamer populations and picosecond interfacial H-bond dynamics of aerosol-OT surfactants. We find that the increased water penetration accelerates H-bond dynamics, with larger RMs showing faster dynamics. These results imply that the changes in the RM structure alter the physical structure of the RM interface and thus alter the solvation dynamics. The findings in this study can be used for developing models for structure-specific solvation dynamics that account for the surfactant packing and hydration at the interface.

Ion/Water Network Structural Dynamics in Highly Concentrated Lithium Chloride and Lithium Bromide Solutions Probed with Ultrafast Infrared Spectroscopy

The structural dynamics of highly concentrated LiCl and LiBr aqueous solutions were observed from 1-4 to 1-16 water molecules per ion pair using ultrafast polarization-selective pump-probe (PSPP) experiments on the OD stretch of dilute HOD. At these high salt concentrations, an extended ion/water network exists with complex structural dynamics. Population decays from PSPP experiments highlight two distinct water components. From the frequency-dependent amplitudes of the decays, the spectra of hydroxyls bound to halides and to water oxygens are obtained, which are not observable in the FT-IR spectra. PSPP experiments also measure frequency-dependent water orientational relaxation. At short times, wobbling dynamics within a restricted angular cone occurs. At high concentrations, the cone angles are dependent on frequency (hydrogen bond strength), but at higher water concentrations (>10 waters per ion pair), there is no frequency dependence. The average cone angle increases as the ion concentration decreases. The slow time constant for complete HOD orientational relaxation is independent of concentration but slower in LiCl than in LiBr. Comparison to structural MD simulations of LiCl from the literature indicates that the loss of the cone angle wavelength dependence and the increase in the cone angles as the concentration decreases occur as the prevalence of large ion/water clusters gives way to contact ion pairs.

Water Dynamics in Aqueous Poly-N-Isopropylacrylamide Below and Through the Lower Critical Solution Temperature

Poly-N-isopropylacrylamide (PNIPAM) is a thermo-responsive polymer that exhibits a reversible structural change from extended chains to aggregates in aqueous solution above its lower critical solution temperature (LCST). Using polarization-selective IR pump-probe spectroscopy, the water orientational dynamics in PNIPAM from below to above the LCST were examined and compared to those of its monomer solution, N-isopropylacrylamide (NIPAM), polyacrylamide, and an acrylamide monomer solution, which are not thermo-responsive. The OD stretch of dilute HOD in H₂O is used as a vibrational probe of water orientational dynamics. Below the LCST of the polymer, NIPAM and PNIPAM solutions exhibited identical water dynamics that were significantly different from those of bulk water, containing both faster and slower components due to solute-water interactions. Therefore, there is no difference in the nature of water interactions with a single NIPAM moiety and a long polymer chain. For all systems, including PNIPAM below and above the LCST, the orientational dynamics were modeled with a bulk water component and a polymer/monomer-associated component based on previous experimental and computational findings. Above the LCST, PNIPAM showed fast water orientational relaxation but much slower long-time dynamics compared to those of NIPAM. The slow component in PNIPAM, which was too slow to be accurately measured due to the limited OD vibrational lifetime, is ascribed to water confined in small voids (<2 nm in diameter) of PNIPAM globules. These results highlight important details about thermo-responsive polymers and the dynamics of their solvation water as they undergo a significant structural change.

Water Dynamics and Structure of Highly Concentrated LiCl Solutions Investigated Using Ultrafast Infrared Spectroscopy

In highly concentrated salt solutions, the water hydrogen bond (H-bond) network is completely disrupted by the presence of ions. Water is forced to restructure as dictated by the water-ion and ion-ion interactions. Using ultrafast polarization-selective pump-probe (PSPP) spectroscopy measurements of the OD stretch of dilute HOD, we demonstrate that the limited water-water H-bonding present in concentrated lithium chloride solutions (up to four waters per ion pair) is, on average, stronger than that occurring in bulk water. Furthermore, information on the orientational dynamics and the angular restriction of water H-bonded to both water oxygens and chloride anions was obtained through analysis of the frequency-dependent anisotropy decays. It was found that, when the salt concentration increased, the water showed increasing restriction and slowing at frequencies correlated with strong H-bonding. The angular restriction of the water molecules and strengthening of water-water H-bonds are due to the formation of a water-ion network not present in bulk water and dilute salt solutions. The structural evolution of the ionic medium was also observed through spectral diffusion of the OD stretch using 2D IR spectroscopy. Compared to bulk water, there is significant slowing of the biexponential spectral diffusion dynamics. The slowest component of the spectral diffusion (13 ps) is virtually identical to the time for complete reorientation of HOD measured with the PSPP experiments. This result suggests that the slowest component of the spectral diffusion reflects rearrangement of water molecules in the water-ion network.

Identical Water Dynamics in Acrylamide Hydrogels, Polymers, and Monomers in Solution: Ultrafast IR Spectroscopy and Molecular Dynamics Simulations

The dynamics and structure of water in polyacrylamide hydrogels (PAAm-HG), polyacrylamide, and acrylamide solutions are investigated using ultrafast infrared experiments on the OD stretch of dilute HOD/H₂O and molecular dynamics simulations. The amide moiety of the monomer/polymers interacts strongly with water through hydrogen bonding (H-bonding). The FT-IR spectra of the three systems indicate that the range of H-bond strengths is relatively unchanged from bulk water. Vibrational population relaxation measurements show that the amide/water H-bonds are somewhat weaker but fall within the range of water/water H-bond strengths. A previous study of water dynamics in PAAm-HG suggested that the slowing observed was due to increasing confinement with concentration. Here, for the same concentrations of the amide moiety, the experimental results demonstrate that the reorientational dynamics (infrared pump-probe experiments) and structural dynamics (two-dimensional infrared spectroscopy) are identical in the three acrylamide systems studied. Molecular dynamics simulations of the water orientational relaxation in aqueous solutions of the acrylamide monomer, trimer, and pentamer are in good agreement with the experimental results and are essentially chain length independent. The simulations show that there is a slower, low-amplitude (<7%) decay component not accessible by the experiments. The simulations examine the dynamics and structure of water H-bonded to acrylamide, in the first solvent shell, and beyond for acrylamide monomers and short chains. The experiments and simulations show that the slowing of water dynamics in PAAm-HG is not caused by confinement in the polymer network but rather by interactions with individual acrylamide moieties.

Simulations of the IR and Raman spectra of water confined in amorphous silica slit pores

Water in nano-scale confining environments is a key element in many biological, material, and geological systems. The structure and dynamics of the liquid can be dramatically modified under these conditions. Probing these changes can be challenging, but vibrational spectroscopy has emerged as a powerful tool for investigating their behavior. A critical, evolving component of this approach is a detailed understanding of the connection between spectroscopic features and molecular-level details. In this paper, this issue is addressed by using molecular dynamics simulations to simulate the linear infrared (IR) and Raman spectra for isotopically dilute HOD in D₂O confined in hydroxylated amorphous silica slit pores. The effect of slit-pore width and hydroxyl density on the silica surface on the vibrational spectra is also investigated. The primary effect of confinement is a blueshift in the frequency of OH groups donating a hydrogen bond to the silica surface. This appears as a slight shift in the total (measurable) spectra but is clearly seen in the distance-based IR and Raman spectra. Analysis indicates that these changes upon confinement are associated with the weaker hydrogen-bond accepting properties of silica oxygens compared to water molecules.

Versatility of Reverse Micelles: From Biomimetic Models to Nano (Bio)Sensor Design

This paper presents an overview of the principal structural and dynamics characteristics of reverse micelles (RMs) in order to highlight their structural flexibility and versatility, along with the possibility to modulate their parameters in a controlled manner. The multifunctionality in a large range of different scientific fields is exemplified in two distinct directions: a theoretical model for mimicry of the biological microenvironment and practical application in the field of nanotechnology and nano-based sensors. RMs represent a convenient experimental approach that limits the drawbacks of the conventionally biological studies in vitro, while the particular structure confers them the status of simplified mimics of cells by reproducing a complex supramolecular organization in an artificial system. The biological relevance of RMs is discussed in some particular cases referring to confinement and a crowded environment, as well as the molecular dynamics of water and a cell membrane structure. The use of RMs in a range of applications seems to be more promising due to their structural and compositional flexibility, high efficiency, and selectivity. Advances in nanotechnology are based on developing new methods of nanomaterial synthesis and deposition. This review highlights the advantages of using RMs in the synthesis of nanoparticles with specific properties and in nano (bio)sensor design.

Intramolecular Vibrational Energy Relaxation of CO2 in Cross-Linked Poly(ethylene glycol) Diacrylate-Based Ion Gels

Ultrafast two-dimensional infrared spectroscopy (2D-IR) and Fourier transform infrared spectroscopy (FTIR) were used to measure carbon dioxide (CO₂) in 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([emim][Tf₂N]), cross-linked low-molecular-weight poly(ethylene glycol) diacrylate (PEGDA), and an ion gel composed of a 50 vol % blend of the two. The center frequency of the antisymmetric stretch, ν₃, of CO₂ shifts monotonically to lower wavenumbers with increasing polymer content, with the largest line width in the ion gel (6 cm^-1). Increasing polymer content slows both spectral diffusion and vibrational energy relaxation (VER) rates. An unexpected excited-state absorbance peak appears in the 2D-IR of cross-linked PEGDA due to VER from the antisymmetric stretch into the bending mode, ν₂. Thirty-two response functions are necessary to describe the observed features in the 2D-IR spectra. Nonlinear least-squares fitting extracts both spectral diffusion and VER rates. In the ion gel, CO₂ exhibits spectral diffusion dynamics that lie between that of the pure compounds. The kinetics of VER reflect both fast excitation and de-excitation of the bending mode, similar to the ionic liquid (IL), and slow overall vibrational population relaxation, similar to the cross-linked polymer. The IL-like and polymer-like dynamics suggest that the CO₂ resides at the interface of the two components in the ion gel.

Bulk-like and Interfacial Water Dynamics in Nafion Fuel Cell Membranes Investigated with Ultrafast Nonlinear IR Spectroscopy

The water confined in the hydrophilic domains of Nafion fuel cell membranes is central to its primary function of ion transport. Water dynamics are intimately linked to proton transfer and are sensitive to the structural features and length scales of confinement. Here, ultrafast polarization-selective pump-probe and two-dimensional infrared vibrational echo (2D IR) experiments were performed on fully hydrated Nafion membranes with sodium counterions to explicate the water dynamics. Like aerosol-OT reverse micelles (AOT RMs), the water dynamics in Nafion are attributed to bulk-like core water in the central region of the hydrophilic domains and much slower interfacial water. Population and orientational dynamics of water in Nafion are slowed by polymer confinement. Comparison of the observed dynamics to those of AOT RMs helps identify local interactions between water and sulfonate anions at the interface and among water molecules in the core. This comparison also demonstrates that the well-known spherical cluster morphology of Nafion is not appropriate. Spectral diffusion of the interfacial water, which arises from structural dynamics, was obtained from the 2D IR experiments taking the core water to have dynamics similar to bulk water. Like the orientational dynamics, spectral diffusion was found to be much slower at the interface compared to bulk water. Together, the dynamics indicate slow reorganization of weakly hydrogen-bonded water molecules at the interface of Nafion. These results provide insights into proton transport mechanisms in fuel cell membranes, and more generally, water dynamics near the interface of confining systems.

Molecular Dynamics Simulations of Macromolecular Crystals

The structures of biological macromolecules would not be known to their present extent without X-ray crystallography. Most simulations of globular proteins in solution begin by surrounding the crystal structure of the monomer in a bath of water molecules, but the standard simulation employing periodic boundary conditions is already close to a crystal lattice environment. With simple protocols, the same software and molecular models can perform simulations of the crystal lattice, including all asymmetric units and solvent to fill the box. Throughout the history of molecular dynamics, studies of crystal lattices have served to investigate the quality of the underlying force fields, correlate the simulated ensembles to experimental structure factors, and extrapolate the behavior in lattices to behavior in solution. Powerful new computers are enabling molecular simulations with greater realism and statistical convergence. Meanwhile, the advent of exciting new methods in crystallography, including femtosecond free-electron lasers and image reconstruction for time-resolved crystallography on slurries of small crystals, is expanding the range of structures accessible to X-ray diffraction. We review past fusions of simulations and crystallography, then look ahead to the ways that simulations of crystal structures will enhance structural biology in the future.

Anomalous water dynamics at surfaces and interfaces: synergistic effects of confinement and surface interactions

In nature, water is often found in contact with surfaces that are extended on the scale of molecule size but small on a macroscopic scale. Examples include lipid bilayers and reverse micelles as well as biomolecules like proteins, DNA and zeolites, to name a few. While the presence of surfaces and interfaces interrupts the continuous hydrogen bond network of liquid water, confinement on a mesoscopic scale introduces new features. Even when extended on a molecular scale, natural and biological surfaces often have features (like charge, hydrophobicity) that vary on the scale of the molecular diameter of water. As a result, many new and exotic features, which are not seen in the bulk, appear in the dynamics of water close to the surface. These different behaviors bear the signature of both water-surface interactions and of confinement. In other words, the altered properties are the result of the synergistic effects of surface-water interactions and confinement. Ultrafast spectroscopy, theoretical modeling and computer simulations together form powerful synergistic approaches towards an understanding of the properties of confined water in such systems as nanocavities, reverse micelles (RMs), water inside and outside biomolecules like proteins and DNA, and also between two hydrophobic walls. We shall review the experimental results and place them in the context of theory and simulations. For water confined within RMs, we discuss the possible interference effects propagating from opposite surfaces. Similar interference is found to give rise to an effective attractive force between two hydrophobic surfaces immersed and kept fixed at a separation of d, with the force showing an exponential dependence on this distance. For protein and DNA hydration, we shall examine a multitude of timescales that arise from frustration effects due to the inherent heterogeneity of these surfaces. We pay particular attention to the role of orientational correlations and modification of the same due to interaction with the surfaces.

Ultrafast dynamics of ionic liquids in colloidal dispersion

Ionic liquid (IL)-surfactant complexes have significance both in applications and fundamental research, but their underlying dynamics are not well understood. We apply polarization-controlled two-dimensional infrared spectroscopy (2D-IR) to study the dynamics of [BMIM][SCN]/surfactant/solvent model systems. We examine the effect of the choice of surfactants and solvent, and the IL-to-surfactant ratio (W-value), with a detailed analysis of the orientation and structural dynamics of each system. Different surfactants create very different environments for the entrapped ILs, ranging from a semi-static micro-environment to a fluxional environment that evolves even faster than the bulk IL. The oil-phase also clearly affects the microscopic dynamics. The anisotropy decay for entrapped ILs completes within 10 ps, which is similar to free thiocyanate ion in water, while a significant reorientation-induced spectral diffusion (RISD) effect is observed. The entrapped ionic liquid are highly dynamic for all W-values, and no core-shell structure is observed. We hypothesize that, instead of an ionic liquid-reverse micelle (IL-RM), the microscopic structure of this system is small colloidal dispersions or pairs of IL and surfactants. A detailed analysis of the polarization-controlled 2D-IR spectra of AOT system reveals a potential ion-exchange mechanism.

Dynamics in the isotropic phase of nematogens using 2D IR vibrational echo measurements on natural-abundance 13CN and extended lifetime probes

The long time scale orientational relaxation of nematogens in the isotropic phase is associated with the randomization of pseudonematic domains, which have a correlation length that grows as the isotropic-to-nematic phase transition temperature is approached from above. Here we begin to address the fast dynamics of the nematogen molecules within the domains using two-dimensional infrared (2D IR) vibrational echo experiments. The problems of performing ultrafast IR experiments in pure liquids are discussed, and solutions are presented. In addition, the issue of short vibrational lifetimes, which limit the ability of 2D IR experiments to examine dynamics over a wide range of times, is addressed. The experiments were performed on the nematogen 4-cyano-4'-pentylbiphenyl (5CB), with the CN stretch initially used as the vibrational probe. Although the CN stretch has a small transition dipole, because the sample is a pure liquid it is necessary to use an exceedingly thin sample to perform the experiments. The small sample volume leads to massive heating effects that distort the results. In addition, the high concentration in the pure liquid can result in vibrational excitation transfer that interferes with the measurements of structural dynamics, and the CN vibrational lifetime is very short (3.6 ps). These problems were overcome by performing the experiments on the natural-abundance (13)CN stretch (5(13)CB), which greatly reduced the absorbance, eliminating the heating problems; also, this stretch has a longer lifetime (7.9 ps). Experiments were also performed on benzonitrile, which showed that the heating problems associated with pure liquids are not unique to 5CB. Again, the problems were eliminated by conducting measurements on the (13)CN stretch, which has an even longer lifetime (20.2 ps) compared with the (12)CN stretch (5.6 ps). Finally, to extend the range of the dynamical measurements, 4-pentyl-4'-thiocyanobiphenyl (5SCB) was synthesized and studied as a dilute solute in 5CB. The CN stretch of 5SCB has a vibrational lifetime of 103 ps, which permits dynamical measurements to 200 ps, revealing the full range of fast structural dynamics in the isotropic phase of 5CB. It is shown that the 5SCB probe reports essentially the same dynamics as 5(13)CB on the short time scale that is observable with the 5(13)CB vibrational probe.

Signatures of nanoconfinement on the linear and nonlinear vibrational spectroscopy of a model hydrogen-bonded complex dissolved in a polar solvent

The one-dimensional IR (1D-IR) absorption and IR pump-probe spectra of a hydrogen stretch in a model hydrogen-bonded complex dissolved in a polar solvent confined in spherical hydrophobic cavities of different sizes were simulated using ground-state mixed quantum-classical dynamics. Due to a thorough analysis of key properties of the complex and solvent from equilibrium trajectory data, we were able to gain insight into the microscopic details underlying the spectra. Both the 1D-IR and IR pump-probe spectra manifested the effects of confinement on the relative stabilities of the covalent and ionic forms of the complex through pronounced changes in their peak intensities and numbers. However, in contrast to the 1D-IR spectra, the time-resolved pump-probe spectra were found to be uniquely sensitive to the changes in the molecular dynamics as the cavity size is varied. In particular, it was found that the variations in the time evolutions of the peak intensities in the pump-probe spectra reflect the differences in the solvation dynamics associated with the various forms of the complex in different locations within the cavities. The ability to detect these differences underscores the advantage of using pump-probe spectroscopy for studying nanoconfined systems.

Computational study of the one- and two-dimensional infrared spectra of a proton-transfer mode in a hydrogen-bonded complex dissolved in a polar nanocluster

The signatures of nanosolvation on the one- and two-dimensional (1D and 2D) IR spectra of a proton-transfer mode in a hydrogen-bonded complex dissolved in polar solvent molecule nanoclusters of varying size are elucidated by using mixed quantum-classical molecular dynamics simulations. For this particular system, increasing the number of solvent molecules successively from N=7 to N=9 initiates the transition of the system from a cluster state to a bulk-like state. Both the 1D and 2D IR spectra reflect this transition through pronounced changes in their peak intensities and numbers, but the time-resolved 2D IR spectra also manifest spectral features that uniquely identify the onset of the cluster-to-bulk transition. In particular, it is observed that in the 1D IR spectra, the relative intensities of the peaks change such that the number of peaks decreases from three to two as the size of the cluster increases from N=7 to N=9. In the 2D IR spectra, off-diagonal peaks are observed in the N=7 and N=8 cases at zero waiting time, but not in the N=9 case. It is known that there are no off-diagonal peaks in the 2D IR spectrum of the bulk version of this system at zero waiting time, so the disappearance of these peaks is a unique signature of the onset of bulk-like behavior. Through an examination of the trajectories of various properties of the complex and solvent, it is possible to relate the emergence of these off-diagonal peaks to an interplay between the vibrations of the complex and the solvent polarization dynamics.

Site-specific measurements of lipid membrane interfacial water dynamics with multidimensional infrared spectroscopy

One route to accessing site-specific dynamical information available with ultrafast multidimensional infrared spectroscopy is the development of robust and versatile vibrational probes. Here we synthesize and characterize a vibrationally labeled cholesterol derivative, (cholesteryl benzoate) chromium tricarbonyl, to probe model lipid membranes, focusing specifically on the membrane-water interface. Utilizing FTIR and polarized-ATR spectroscopies, we determine the location of the chromium tricarbonyl motif to be situated at the water-membrane interface with an orientation of 46 ± 2° relative to the vector normal to the membrane surface. We test the dynamical sensitivity of the (cholesteryl benzoate) chromium tricarbonyl label with two different nonlinear infrared spectroscopy methods, both of which show that the probe is well-suited to the study of membrane dynamics as well as the dynamics of water at the membrane interface. The metal carbonyl vibrational probe located at the surface of a bicelle exhibits spectral diffusion dynamics induced by membrane hydration water that is roughly three times slower than observed using a nearly identical vibrational probe in bulk water.

Dynamics in the interior of AOT lamellae investigated with two-dimensional infrared spectroscopy

The dynamics inside the organic regions of aerosol-OT (AOT)/water mixtures in the lamellar mesophase, bicontinuous cubic (BC) phase, and in an analogous molecule without the charged sulfonate headgroup are investigated by observing spectral diffusion, orientational relaxation and population relaxation using ultrafast two-dimensional infrared (2D IR) vibrational echo spectroscopy and IR pump-probe experiments on the asymmetric CO stretch of a vibrational probe, tungsten hexacarbonyl (W(CO)6). The water layer thickness between the bilayer planes in the lamellar phase was varied. For comparison, the dynamics of W(CO)6 in the normal liquid bis(2-ethylhexyl) succinate (EHS), which is analogous to AOT but has no charged sulfonate headgroup, were also studied. The 2D IR experiments measure spectral diffusion, which results from the structural evolution of the system. Spectral diffusion is quantified by the frequency-frequency correlation function (FFCF). In addition to a homogeneous component, the FFCFs are biexponential decays with fast and slow time components of ∼12.5 and ∼150 ps in the lamellar phase. Both components of the FFCF are independent of the number of water molecules per headgroup for the lamellae, but they slow somewhat in the BC phase. The dynamics in the ordered phases are in sharp contrast to the dynamics in EHS, which displays fast and slow components of the FFCF of 5 and 80 ps, respectively. As the hydration level of AOT increases, vibrational lifetime decreases, suggesting some change in the local environment of W(CO)6 with water content.

Vibrational spectroscopy of water at interfaces

Understanding liquid water's behavior at the molecular level is essential to progress in fields as disparate as biology and atmospheric sciences. Moreover, the properties of water in bulk and water at interfaces can be very different, making the study of the hydrogen-bonding networks therein very important. With recent experimental advances in vibrational spectroscopy, such as ultrafast pulses and heterodyne detection, it is now possible to probe the structure and dynamics of bulk and interfacial water in unprecedented detail. We consider here three aqueous interfaces: the water liquid-vapor interface, the interface between water and the surfactant headgroups of reverse micelles, and the interface between water and the lipid headgroups of aligned multi-bilayers. In the first case, sum-frequency spectroscopy is used to probe the interface. In the second and third cases, the confined water pools are sufficiently small that techniques of bulk spectroscopy (such as FTIR, pump-probe, two-dimensional IR, and the like) can be used to probe the interfacial water. In this Account, we discuss our attempts to model these three systems and interpret the existing experiments. For the water liquid-vapor interface, we find that three-body interactions are essential for reproducing the experimental sum-frequency spectrum, and presumably for the structure of the interface as well. The observed spectrum is interpreted as arising from overlapping and canceling positive and negative contributions from molecules in different hydrogen-bonding environments. For the reverse micelles, our theoretical models confirm that the experimentally observed blue shift of the water OD stretch (for dilute HOD in H(2)O) arises from weaker hydrogen bonding to sulfonate oxygens. We interpret the observed slow-down in water rotational dynamics as arising from curvature-induced frustration. For the water confined between lipid bilayers, our theoretical models confirm that the experimentally observed red shift of the water OD stretch arises from stronger hydrogen bonding to phosphate oxygens. We develop a model for heterogeneous vibrational lifetime distributions, and we implement the model to calculate isotropic and anisotropic pump-probe decays. We then compare these results with experimental data. Clearly, recent experimental advances in vibrational spectroscopy have led to beautiful new results, providing information about the structure and dynamics of water at interfaces. These experimental and concomitant theoretical advances (particularly the unified theoretical framework of non-linear response functions) have greatly contributed to our understanding of this unique and important substance.

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.

Reduced coupling of water molecules near the surface of reverse micelles

We report on vibrational dynamics of water near the surface of AOT reverse micelles studied by narrow-band excitation, mid-IR pump-probe spectroscopy. Evidence of OH-stretch frequency splitting into the symmetric and asymmetric modes is clearly observed for the interfacial H(2)O molecules. The polarization memory of interfacial waters is preserved over an exceptionally extended >10 ps timescale which is a factor of 100 longer than in bulk water. These observations point towards negligibly small intermolecular vibrational coupling between the water molecules as well as strongly reduced water rotational mobility within the interfacial water layer.

Time-dependent fifth-order bands in nominally third-order 2D IR vibrational echo spectra

Progress in the field of 2D IR vibrational spectroscopy has been bolstered by the production of intense mid-IR laser pulses. As higher-energy pulses are employed, a concomitant increase occurs in the likelihood of fifth-order contributions to the 2D IR spectra. We report the appearance of fifth-order signals in 2D IR spectra of CO bound to the active site of the enzyme cytochrome P450(cam) with the substrate norcamphor. Two bands with novel time dependences, one on the diagonal and one off-diagonal, are not accounted for by normal third-order interactions. These bands are associated with a ν = 1-2 vibrational transition frequency. Both bands decay to 0 and then grow back in with opposite sign. The diagonal band is positive at short time, decays to 0, reappears with negative sign, before eventually decaying to 0. The off-diagonal band is negative at short time, decays to 0, reappears positive, and then decays to 0. The appearance and time dependence of these bands are characterized. Understanding these fifth-order bands is useful because they may be misidentified with time-dependent bands that arise from other processes, such as chemical exchange, vibrational coupling, or energy transfer. The presence and unusual time dependences of the fifth-order bands are reproduced with model calculations that account for the fact that vibrational relaxation from the ν = 2 to 1 level is approximately a factor of 2 faster than that from the ν = 1 to 0 level.

Quantum process tomography of excitonic dimers from two-dimensional electronic spectroscopy. I. General theory and application to homodimers

Is it possible to infer the time evolving quantum state of a multichromophoric system from a sequence of two-dimensional electronic spectra (2D-ES) as a function of waiting time? Here we provide a positive answer for a tractable model system: a coupled dimer. After exhaustively enumerating the Liouville pathways associated to each peak in the 2D-ES, we argue that by judiciously combining the information from a series of experiments varying the polarization and frequency components of the pulses, detailed information at the amplitude level about the input and output quantum states at the waiting time can be obtained. This possibility yields a quantum process tomography (QPT) of the single-exciton manifold, which completely characterizes the open quantum system dynamics through the reconstruction of the process matrix. In this manuscript, we present the general theory as well as specific and numerical results for a homodimer, for which we prove that signals stemming from coherence to population transfer and vice versa vanish upon isotropic averaging, therefore, only allowing for a partial QPT in such case. However, this fact simplifies the spectra, and it follows that only two polarization controlled experiments (and no pulse-shaping requirements) suffice to yield the elements of the process matrix, which survive under isotropic averaging. Redundancies in the 2D-ES amplitudes allow for the angle between the two site transition dipole moments to be self-consistently obtained, hence simultaneously yielding structural and dynamical information of the dimer. Model calculations are presented, as well as an error analysis in terms of the angle between the dipoles and peak amplitude extraction. In the second article accompanying this study, we numerically exemplify the theory for heterodimers and carry out a detailed error analysis for such case. This investigation reveals an exciting quantum information processing (QIP) approach to spectroscopic experiments of excitonic systems, and hence, bridges an important gap between theoretical studies on excitation energy transfer from the QIP standpoint and experimental methods to study such systems in the chemical physics community.