Efficient Intramolecular Vibrational Excitonic ... - ACS Publications

Dec 24, 2017 - Due to structural symmetry (D3h point group), 12 carbonyl (C O) groups in the Ru3(CO)12 complex contribute to mainly three excitonic ca...
0 downloads 9 Views 2MB Size
Article Cite This: J. Phys. Chem. B 2018, 122, 1296−1305

pubs.acs.org/JPCB

Efficient Intramolecular Vibrational Excitonic Energy Transfer in Ru3(CO)12 Cluster Revealed by Two-Dimensional Infrared Spectroscopy Xueqian Dong,†,‡ Fan Yang,*,† Juan Zhao,†,‡ and Jianping Wang*,†,‡ †

Beijing National Laboratory for Molecular Sciences; Molecular Reaction Dynamics Laboratory, CAS Research/Education Center for Excellence in Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, P. R. China ‡ University of Chinese Academy of Sciences, Beijing 100049, P. R. China S Supporting Information *

ABSTRACT: Trinuclear transition-metal carbonyl complex dodecacarbonyl triruthenium (Ru3(CO)12) is considered as one of the paradigms in cluster chemistry, which plays an important role in photocatalysis, photoenergy conversion, and synthetic chemistry. Due to structural symmetry (D3h point group), 12 carbonyl (CO) groups in the Ru3(CO)12 complex contribute to mainly three excitonic carbonyl stretching modes: E′ (radial), A2″ (axial), and E′ (axial). In this work, efficient intramolecular vibrational energy redistribution (IVR) processes among the three modes in this Ru− CO complex were observed to occur on the time scale of tens of picoseconds. The IVR processes were characterized in detail using a kinetic model and fitting to the waiting-timedependent diagonal and off-diagonal signals of ultrafast two-dimensional infrared spectroscopy. In addition, the diagonal anharmonicities of the three CO stretching modes were determined to be quite close to one another, and the couplinginduced cross peaks were invariant because this Ru3(CO)12 cluster does not show picosecond fluxionality and hence their contributions were neglected in modeling the IVR processes. Our results provide a benchmark for understanding the excitonic nature of the vibrational excited states of the carbonyl vibrators and the associated efficient vibrational energy-flow pathways, in such multicentered transition-metal complexes, which are of key importance to their functions. fluxional behavior, and photolysis of this metal carbonyl cluster. In addition, theoretical insights have been obtained on Ru3(CO)12 using several density functionals, whose geometrical parameters and energetics have been predicted.21 Intuitively speaking, under a broad-band ultrafast infrared laser excitation, because there are multiple CO stretching modes involved in Ru3(CO)12, the intramolecular vibrational relaxation (or intramolecular vibrational energy redistribution, IVR) among these CO modes could have complicated pathways and/or varied time scales, which, however, still remains unknown for this multicentered Ru−carbonyl complex. In particular, a recent IR pump−probe spectroscopic study of Os3(CO)12, which is a similar transition-metal carbonyl complex to Ru3(CO)12, showed that no vibrational energy transfer occurred among the multiple CO stretching modes in such trinuclear carbonyl cluster.22 An effective and comprehensive way to address these energy transfer and associated vibrational dynamics questions is to use two-

I. INTRODUCTION The cluster chemistry involving transition-metal carbonyls has been developed since the discovery of metal carbonyls, such as Co2(CO)8 and Fe2(CO)9.1,2 Transition-metal carbonyl clusters containing two or more metal atoms surrounded by carbonyl ligands have been widely used as solution-phase homogeneous catalysts3,4 and also synthetic precursors5,6 to other carbonyl complexes. Dodecacarbonyl triruthenium, Ru3(CO)12, is one of the representative transition-metal carbonyl clusters, which is considered as a model of trinuclear metal complex. Moreover, this trinuclear metal carbonyl cluster has been known to play an important role in photocatalysis,7 synthetic chemistry,8 as well as photochemistry.9−11 In condensed phases, the Ru3(CO)12 complex has 12 nonbridged CO groups and a D3h symmetry configuration corresponding to the most stable trinuclear core structure, which is quite typical in cluster chemistry. The overall molecular geometry and the ligand rearrangement mechanism of the Ru3(CO)12 cluster have been studied extensively. A variety of methods, such as NMR spectroscopy,12 single-crystal X-ray diffraction,9,13−15 and infrared and Raman spectroscopies,16−20 have been used to examine the molecular structure, © 2017 American Chemical Society

Received: October 11, 2017 Revised: December 1, 2017 Published: December 24, 2017 1296

DOI: 10.1021/acs.jpcb.7b10067 J. Phys. Chem. B 2018, 122, 1296−1305

Article

The Journal of Physical Chemistry B

Figure 1. Vibrational CO stretching modes of Ru3(CO)12 cluster, with their symmetry coordinates arranged from low frequency to high frequency (from left to right). Here, “axial” and “radial” represent the symmetry coordinates, which dominate in the vibrations.

a 35 fs, 1 mJ pulse energy, 800 nm laser pulse was generated by a regenerative Ti:sapphire laser amplifier. The produced nearIR pulse first went into an optical parametric amplifier and then into a difference frequency mixer to further generate a mid-IR pulse with its frequency centered at 2020 cm−1 with a spectral width of ca. 240 cm−1 (full width at half-maximum, FWHM). The mid-IR pulse was split into three excitation beams (ca. 400 nJ each) and focused together on the liquid sample to generate 2D IR signal, and a fourth weak mid-IR beam was used as a local oscillator for heterodyning amplification of the 2D IR signal. The delay time between the first and second excitation pulses (k1 and k2) is the so-called coherence time (τ), whereas that between the second and third excitation pulses (k2 and k3) is called waiting time (or dynamical time, Tw). The generated 2D IR signal is directed into an IR monochromator that is equipped with a 64-element mercury−cadmium−telluride (MCT) array detector for data acquisition along the detection time (t) using the so-called spectral interferometry. In this way, the time-domain 2D IR data are half Fourier-transformed by the MCT array detector (along the detection frequency, ωt). Rephasing (−k1 excites sample first) and nonrephasing (k2 excites sample first) 2D IR data were recorded as a function of τ at a step of 5 fs for a given Tw. A series of Tw-dependent 2D IR spectra were collected to characterize the IVR dynamics. The 2D IR measurement was carried out at room temperature (22 °C) under dark condition. II.III. Density Functional Theory (DFT) Computations. Molecular geometry optimization and harmonic normal-mode vibrational frequency analysis for the Ru3(CO)12 complex in gas phase were carried out using density functional theory (DFT) at the B3LYP level. The lanl2dz pseudopotential basis was used for the transition-metal atom (Ru), and the 6-31+G* basis set was used for other atoms (C, O). The computations were performed using Gaussian.38

dimensional infrared (2D IR) spectroscopy, as also pointed by these authors. In a specific 2D IR spectrum, a set of normal modes can be excited at the same time, including the fundamental, first overtone, and first combination transitions. These vibrational transitions are shown as diagonal and offdiagonal 2D IR signals. Waiting time (or dynamical time) between a series of 2D IR spectra is the time window to monitor the 2D IR spectral change, in both line shape, which reflects spectral diffusion that is mostly related to the evolution of solvent and neighboring chemical composition, and signal intensity, which reflects either vibrational relaxation (diagonal signals) or vibrational energy transfer among different states (off-diagonal signals). Thus, by analyzing the shape and intensity changes in a series of waiting-time-dependent 2D IR spectra, one obtains information about molecular structural dynamics (for both solute and solvent),23−31 and vibrational energy transfer (occurred intramolecularly and/or intermolecularly).32−35 We focus on the latter aspect in the study of the Ru3(CO)12 complex in CHCl3 solvent. This work is the first report on the IVR process among different carbonyl stretching modes of trinuclear metal carbonyl cluster. There are three IR-active CO stretching normal modes with large transition dipole in the Ru3(CO)12 complex, indicating their excitonic nature. Directly monitoring the time evolution of the populations of energy donor and acceptor via these exciton-like vibrational modes, as the intramolecular energy-transfer process of Ru3(CO)12 taking place, can yield a set of kinetical (reaction rate) and dynamical (structural) insights into the vibrational energy-transfer processes of these multicentered transition-metal clusters. This article is organized as follows. Section II describes experimental and computational methods used in this work. Section III contains detailed analysis on the vibrational properties of structurally optimized Ru3(CO)12 complex, and more importantly, on the timedependent 2D IR experimental results of this transition-metal complex, based on which, the IVR processes and dynamics in this molecule are modeled and discussed.

III. RESULTS AND DISCUSSION III.I. Carbonyl Stretching Modes. To assist peak assignment of the CO stretching modes in the IR spectroscopies of this work described below, we first performed DFT calculations on the Ru3(CO)12 cluster. Four IR-active carbonyl stretching modes of the optimized Ru3(CO)12 cluster are shown in Figure 1. The detailed computed vibrational parameters are also listed in Table 1. The vibrational CO stretching modes can be grouped into two types through the linear combination of the 12 local modes: one is the axial carbonyl vibration and the other is the radial carbonyl vibration. The carbonyl groups in the axial mode are opposite to each other spatially through the Ru atom, whereas those in the radial mode are opposite along the Ru−Ru bond. The difference of the intensity between the radial and axial modes is due to the π-bonding effect in the cluster.17,39,40 The highest band (with a vibrational frequency of 2122.9 cm−1)

II. MATERIAL AND METHODS II.I. Sample Preparation and Linear IR Experiment. Ru3(CO)12 was purchased from Sigma-Aldrich and used without purification. For linear and nonlinear IR measurements, the Ru 3 (CO) 12 cluster was dissolved in CHCl 3 at a concentration of 4.6 mM. The sample solution was loaded into a homemade IR sample cell, which consisted of two 2 mm thick CaF2 optical windows separated by 50 μm thick Teflon spacer. Fourier transform infrared (FTIR) spectrum (linear IR) of Ru3(CO)12 cluster was collected using a Nicolet 6700 FTIR spectrometer (Thermo Electron) with 1 cm−1 resolution and 32-scan average at room temperature (22 °C). II.II. Two-Dimensional IR Experiment. The 2D IR experimental setup has been described elsewhere.24,36,37 Briefly, 1297

DOI: 10.1021/acs.jpcb.7b10067 J. Phys. Chem. B 2018, 122, 1296−1305

Article

The Journal of Physical Chemistry B

reasonable agreement with previous studies of Ru3(CO)12 cluster in solution phases and also in solid state.18−20,48,49 In addition, it can be seen that the weak peak of the E′ (radial) mode overlaps to some extent with that of the A2″ (axial) mode in the spectrum. The assignment of the E′ (radial, B2) component in the FTIR spectrum is not reliable and hence is considered together in the E′ (radial) mode. The FTIR spectrum was fitted by Voigt line shape, and the fitting results are listed in Table 2. During the fitting, peak

Table 1. Calculated IR-Active Vibrational Frequency (ω) and Intensity (I) of Ru3(CO)12 Cluster in the CO Stretching Region ω/cm−1

assignment

I/(kM·mol−1)

2078.3 2080.1 2103.4 2122.9

E′ (radial, B2) E′ (radial, A1) A2″ (axial) E′ (axial)

22.0 764.0 2939.9 4763.2

is assigned to an axially dominated E′ mode, whose high intensity is due to the transition dipole strength borrowing from the local CO vibrators.17,18,40 The second highfrequency band (2103.4 cm−1) represents mostly a pure axial A2″ mode,39 the carbonyl groups of which vibrate in phase showing lower intensity than those vibrate out of phase (such as in the E′ axial mode).39,41 Two weak-intensity vibrations located at 2080.1 and 2078.3 cm−1 are assigned to the E′ (radial, A1) and E′ (radial, B2) modes, respectively. These two modes refer to their local symmetry coordinates, where an isolated Ru(CO)4 is considered as C2v symmetry.18,39 The two bands as the consequence of factor-group splitting caused by static field effects appear to overlap,16,18,42 and the E′ (radial, B2) mode is included together with the E′ (radial, A1) in this work because of its extremely low intensity and its overlapping frequency position with the latter. Moreover, the vibrational coupling between these CO groups was calculated using the wave function demixing method,43−46 and the results are shown in Table S1 of the Supporting Information (SI). How these local CO groups (labeled in Figure S1 of the SI) contribute to the three normal modes is further demonstrated in the pairwise vibrational coupling analysis (Table S1) and also in the potential energy distribution (PED) analysis (Table S2 of the SI).47 These results clearly indicate the vibrational excitonic nature of the limited number of IR-active carbonyl stretching modes in this system. III.II. Linear IR Spectrum. The infrared absorption spectrum and molecular structure of the Ru3(CO)12 cluster solvated in CHCl3 are shown in Figure 2. Three absorption bands are observed around the 2000−2060 cm−1 region. The peaks from high- to low-frequency side are assigned to the E′ (axial), A2″ (axial), and E′ (radial) modes, respectively. The assignment for these vibrational modes (frequency positions and intensities) is based on the DFT computations, which is in

Table 2. Fitting Parameters of Linear IR Spectrum of Ru3(CO)12 Cluster in the CO Stretching Frequency Region mode

ω/cm−1

FWHM/cm−1

area/arb unit

E′ (radial, B2) + E′ (radial, A1) A2″ (axial) E′ (axial)

2010.6 2029.1 2061.0

18.6 12.3 8.0

2.9 5.3 5.6

positions of these CO stretching modes obtained from the projection of 2D IR spectrum (see below, Section III.IV) were used. The peak frequencies of these modes in the FTIR spectra are somewhat different from the DFT results because the influence of solvent was not considered in our calculation and also because of the inaccuracy of the DFT method in predicting the vibrational frequency of molecules. The frequency-scaling factor usually applied for various density functionals was not considered either in our computation. But the computation does agree with the experimental results in relative peak positions and peak intensities of the three more major IR bands shown in Figure 2. In addition, one sees that the highestfrequency peak (E′, axial) is the narrowest among the three IR bands (listed in Table 2), suggesting more homogeneous broadening in this peak. However, the spectral line width contribution is not the major concern of this work and hence will not be discussed further. III.III. Two-Dimensional IR Spectra of the Ru3(CO)12 Cluster. The real parts of the purely absorptive 2D IR spectra of Ru3(CO)12 cluster dissolved in CHCl3 at different waiting times are shown in Figure 3. In these 2D IR graphs, both diagonal and off-diagonal peaks can be clearly seen, and their origins can be assigned without any ambiguity. The 2D IR signals always appear in pairs, including the υ0→1 and υ1→2 transitions. The positive signal in red is denoted by the υ0→1 transition, whereas the negative one in blue represents the υ1→2 transition. Along the diagonal, three positive−negative peak pairs are observed, which correspond to the peaks in the linear IR spectrum of Ru3(CO)12 cluster. However, the diagonal peak of the E′ (radial) mode shows only the negative signal because its positive signal overlaps with the negative one of peak 2 (in panel D). In the meantime, the negative signal of peak 2 appears to be smaller than its positive counterpart. The intensity of the E′ (axial, E′a, labeled in panel D) mode is the strongest, whereas that of the E′ (radial, E′r) mode is the weakest. Moreover, it is observed that all of the υ1→2 transitions decrease in frequency along the ωt axis with respect to their υ0→1 transitions as a consequence of considerable vibrational anharmonicity of the CO stretching normal mode.50,51 Therefore, the actual frequency differences between υ0→1 and υ1→2 transitions represent the diagonal anharmonicities (Δ), and more details are given in Section III.V. Six off-diagonal peak pairs are clearly shown in the purely absorptive 2D IR spectra. Peaks 1 and 2 are the off-diagonal

Figure 2. Experimental linear IR spectrum of Ru3(CO)12 cluster dissolved in CHCl3 in the CO stretching region (solid line) and its fitting using Voigt line shape (dashed line). The inset shows the molecular structure of Ru3(CO)12. The fitting residual (dashed line) is given in the bottom panel. 1298

DOI: 10.1021/acs.jpcb.7b10067 J. Phys. Chem. B 2018, 122, 1296−1305

Article

The Journal of Physical Chemistry B

Figure 3. Representative waiting-time-dependent purely absorptive 2D IR spectra of Ru3(CO)12 cluster dissolved in CHCl3 in the CO stretching region. The waiting times are given in each panel. Six representative off-diagonal peak pairs are labeled with dashed squares in (D) by numbers 1−6. In (B) and (C), the red dashed lines indicate 2D IR slices used to determine the anharmonicities (see Figures 5 and S2 of the SI). Each 2D IR spectrum is accompanied with a color bar to show its signal intensity scale.

structure on the picosecond time scale (no fluxional dynamics yet), the coupling-induced cross peaks should have an unchanged intensity for each cross peak (with respect to the intensities of any pair of couple modes). Because the observed off-diagonal 2D IR signals change as a function of waiting time (as shown clearly in Figure 3), the coupling contribution is not dominant and can be neglected on the basis of its nearly invariant nature on the picosecond time scale. The third possibility is the intramolecular vibrational energy redistribution, which fits the evolving feature of a network of cross peaks in the waiting-time-dependent 2D IR spectra shown in Figure 3. On the basis of these considerations, the change of these off-diagonal peak amplitudes in purely absorptive 2D IR spectra in Figure 3 is believed to be mainly due to the intramolecular vibrational energy transfer. Recent studies have showed that picosecond IVR processes among the CO stretching vibrational energy levels also occur in other similar transition-metal carbonyl complexes.33,37,55 However, it has to be pointed out that the contribution of anharmonic coupling to the off-diagonal signals cannot be completely separated from the IVR contributions, which will be further discussed in Section III.V. III.IV. Measuring Accurate Frequency Positions. The peak of the weak E′ (radial) mode overlaps that of the A2″ (axial) mode in the FTIR spectrum, so the accurate peak position of the E′ (radial) mode cannot be read out directly from the FTIR spectrum. Two-dimensional IR spectroscopy can be used to study the underlying vibrational transitions in the case of multivibrational modes. This is because the signal intensity in a 2D IR spectrum is proportional to the fourth power of associated transition dipole strength, whereas that in an FTIR spectrum is only proportional to the square of the transition dipole strength. An early study of Cp2Fe2(CO)4 and Cp2Ru2(CO)4 has revealed the precise transition frequency for

peaks between the A2″ (axial) and E′ (radial) modes, whereas peaks 3 and 4 are generated from the E′ (axial) and E′ (radial) modes. Similarly, peaks 5 and 6 between the E′ (axial) and A2″ (axial) modes are also present in Figure 3. At early Tw (2 ps), strong off-diagonal peak pairs (peaks 5 and 6) and weak offdiagonal peak pairs (peaks 3 and 4) appear. The off-diagonal peak 3 shows only its negative signal, and in the meantime, the negative signal of peak 5 is weaker than its positive counterpart. However, the diagonal peak of the E′ (radial) mode and the associated off-diagonal peaks (peaks 1 and 2) cannot be identified completely. As the waiting time increases (Tw = 10 ps), the off-diagonal peaks (peaks 1 and 2) become visible, but still only the negative signal of peak 1 is observable. As the waiting time further evolves, one sees that the intensities of diagonal and off-diagonal peaks generally decrease due to the vibrational relaxation. Although the off-diagonal peak signals seem to become distinct, the amplitudes generally decrease, which is particularly true for peak 6 (see Section III.VI). The appearance of off-diagonal peaks in such metal carbonyl complexes usually has three origins. The first possibility is picosecond-time-scale interchange of axial−equatorial CO groups due to the fluxional property of metal carbonyl complex.52 An early study53 indeed suggested fast fluxional dynamics for the Ru3(CO)12 cluster in solution. Later studies have shown that the activation energy for the ligand exchange process in the Ru3(CO)12 cluster is less than 5.0 kcal·mol−1.12 However, it was further shown that the ligand rearrangement occurs on the time scale of microseconds and the exchange rates are temperature-dependent.54 Thus, the fluxional origin of these picosecond-time-scale cross peaks observed in Figure 3 is unlikely. The second possibility for the off-diagonal cross peaks in Figure 3 is anharmonic vibrational coupling, which occurs inevitably for coupled modes, such as a collection of carbonyl stretches studied in this case. However, for a more or less stable 1299

DOI: 10.1021/acs.jpcb.7b10067 J. Phys. Chem. B 2018, 122, 1296−1305

Article

The Journal of Physical Chemistry B the individual isomer from off-diagonal peak regions of 2D IR spectra.33 In our case, the υ1→2 transitions instead of υ0→1 transitions of off-diagonal peaks (peaks 1 and 3) are visible in the purely absorptive 2D IR spectra of Ru3(CO)12 cluster, so it is difficult to extract the frequency position of the E′ (radial) mode from the off-diagonal peaks in a straightforward way. Thus, 2D IR spectral projection along the diagonal on the ωt axis was used to determine the peak position of the E′ (radial) mode. The diagonal projection and the FTIR spectrum are plotted together in Figure 4. As can be seen, the strong peaks of

Figure 5. Representative slice of the 2D IR real-part spectrum (empty circle) at Tw = 10 ps at ωτ = 2060 cm−1 and along the ωt-axis. The υ0→1 transition is in red, and the υ1→2 transition is in blue. Peak fitting is also shown by the black solid line.

Table 3. Fitting Parameters and Determining the Diagonal and Off-Diagonal Anharmonicities of the IR-Active Carbonyl Stretching Modes at ωτ = 2060 cm−1 a

Figure 4. Projection of the 2D IR real-part spectrum at Tw = 10 ps on the ωt-axis along the diagonal (solid line) and FTIR spectrum (dot line) of Ru3(CO)12 cluster. The peak positions of these CO stretching modes are shown by dashed lines.

mode

ωt/cm−1

transition

FWHM/cm−1

areab/arb unit

Δ/cm−1

E′ (radial)

1996.1 2009.3 2015.8 2028.8 2048.9 2061.1

υ1→2 υ0→1 υ1→2 υ0→1 υ1→2 υ0→1

16.7 17.8 12.7 9.9 10.6 10.4

−2.4 1.7 −3.6 3.9 −15.2 15.8

13.2

A2″ (axial) E′ (axial) a

the E′ (axial) and A2″ (axial) modes in the diagonal projection match their FTIR bands quite well. The E′ (radial) mode also forms an independent peak in the diagonal projection, even though the peak is quite small compared to the other two highfrequency components. On the basis of this small peak, the peak-frequency position of the E′ (radial) mode can be directly determined, which is 2010 cm−1. In addition, one has to note that the intensities of the three peaks in the 2D IR diagonal projection are not proportional to those in the FTIR spectrum because for one, the projected peaks are taken from the vibrational relaxed 2D IR diagonal signals, which may have different lifetimes, whereas in the FTIR spectrum, the intensities are measured under equilibrium condition, and for two, these signal intensities are proportional to the transition dipole strength differently, as mentioned above. III.V. Extracting Diagonal and Off-Diagonal Anharmonicities of the Carbonyl Stretches. Two-dimensional IR spectroscopy is usually sensitive to three lowest vibrational transition energy levels for a given set of anharmonic oscillators (υ = 0, 1, and 2), which can hence provide a direct measure of the vibrational anharmonicity. Previous studies have extracted the anharmonicity from the slice of purely absorptive 2D IR spectra along the probe axis at a typical pump frequency.56,57 In our case, the representative 2D IR slice (Figure 5) was taken from the 2D IR real-part spectrum at Tw = 10 ps along the ωt axis at ωτ = 2060 cm−1, which is indicated by the red dashed line in Figure 3B. This 2D IR slice was fitted by the Voigt line shape, where the positive peaks in red represent the υ0→1 transitions and the negative ones in blue indicate the υ1→2 transitions. The anharmonicity for a given mode is defined as a positive value, i.e., Δ = υ0→1 − υ1→2. During the fitting, the υ0→1 transitions of these CO stretching modes were taken from the results obtained in Figure 4, and the fitting results are listed in Table 3. Because the off-diagonal peaks (peaks 3 and 5) are mainly due to the intramolecular energy transfer, the

13.0 12.2

See Figure 5. bVoigt line shapes are used during fitting.

anharmonicity values of these off-diagonal peaks are equal to the diagonal anharmonicities of the E′ (radial) and A2″ (axial) modes, respectively. The obtained diagonal anharmonicity of the radial E′ mode (13.2 cm−1) is only slightly larger than those of the axial A2″ (13.0 cm−1) and axial E′ modes (12.2 cm−1). These results suggest the overall excitonic behavior of the three major IR bands observed in our linear and nonlinear IR spectroscopies. The diagonal anharmonicities of carbonyl modes are also in agreement with previous calculated and experimentally determined results for similar transition-metal carbonyl stretching modes.51,58 Moreover, 2D IR slices at other waiting times were also fitted, and marginal variation in the anharmonicities was observed among those diagonal and the associated off-diagonal signals that were vertically aligned (along the ωτ-axis). The results are given in Figure S2 with fitting parameters listed in Tables S3−S5 in the SI for comparison. These results also indicate insignificant anharmonicity change for a given 2D IR signal and further supporting the argument presented above for this Ru−carbonyl complex on the picosecond time scale. On the other hand, in an excitonic picture, it is known that the magnitude of a mixed-mode off-diagonal anharmonicity for a pair of vibrational modes will be non-negligible if the two modes are coupled.59 And the diagonal and off-diagonal anharmonicities are molecular-structure-dependent. This is a generally known feature of the vibrational excitonic states consisting, for example, of a group of carbonyl stretching vibrators. However, the fitting results mentioned above do not yield any explicit information about the mixed-mode offdiagonal anharmonicities among these modes. Because such offdiagonal anharmonicity is closely related to the anharmonic shift of the combination band, it gives rise to (couplinginduced) cross peaks, as mentioned in Section III.III. 1300

DOI: 10.1021/acs.jpcb.7b10067 J. Phys. Chem. B 2018, 122, 1296−1305

Article

The Journal of Physical Chemistry B

illustrated in Scheme 1. Selective absolute parts of 2D IR nonrephasing spectra at typical waiting times are given in

It should be pointed out that even at early waiting times, it is difficult to separate coupling-induced cross peak from IVRinduced cross peak if both exist. As one example, we examined early-waiting-time 2D IR spectra. The real parts of the nonrephasing 2D IR spectra at two early times (100 and 200 fs, Figure S3 in the SI) indeed show the presence of cross peaks, however, whose origin could still be partially vibrational coupling (that is structure-dependent and hence likely to be time-independent) and partially IVR (that is time-dependent). These two types of cross peaks will also contribute to the transient broad-band pump−broad-band probe infrared spectra. The latter expectation has actually been demonstrated in a previous work22 of a similar trinuclear transition-metal carbonyl complex, namely, Os3(CO)12. In addition, 2D IR spectra at Tw = 0 cannot be used to separate the coupling-induced cross peak from the IVR-induced cross peak because of the complication of spectral contribution when four IR pulses overlap at zero time in the 2D IR experiment. There is one way to measure coupling-induced cross peak, that is, at a fairly long waiting time when the IVR processes finish but the vibrational relaxation does not and the signal-tonoise ratio of 2D IR signal is still reasonably good. This is beyond the scope of this work and will not be discussed further. III.VI. Obtaining Intramolecular Vibrational Energy Redistribution Rate. Two-dimensional IR spectroscopy can provide a measure for intramolecular vibrational energy-transfer dynamics between two different modes. The peak amplitudes in 2D IR spectra as a function of waiting time can be used to obtain the IVR rate. To this end, the 2D IR cross-peak intensities can be taken from their positive and negative extremes.32 Previous studies have shown that the intensity ratio of cross peak to diagonal peak should be used to correct the vibrational relaxation when evaluating the IVR processes.33 In our recent study, the IVR time constants in Mn(CO)5Br were obtained by fitting the waiting-time evolution of the 2D IR diagonal and off-diagonal amplitudes (or simply called intensities) simultaneously.37 To fit the IVR rate, we used a series of Tw-dependent absolute nonrephasing spectra. It should be pointed out that the data collected in conventional four-wave mixing method require proper phasing of the rephasing and nonrephasing 2D IR spectra to obtain the purely absorptive 2D IR spectra, as shown in Figure 3. It is well known that such purely absorptive 2D IR spectra have many advantages (for example, clearly showing negative and positive diagonal and off-diagonal peaks). However, one has to know that a careful phasing of the timedomain rephasing and nonrephasing signals is very critical. Thus, it can be a very tedious work if a series of waiting-timedependent 2D IR spectra are needed. This problem can be circumvented using the absolute part of the 2D IR signal, if one is only interested in evaluating the signal intensity of either a diagonal or an off-diagonal signal. In this case, the intensities of both the positive and negative parts of a 2D IR signal can be taken into account, if a 2D frequency area is properly chosen. On the other hand, if this procedure is taken, then one only has to collect either the rephasing or nonrephasing part of 2D IR spectra because the absolute-magnitude 2D IR spectra for the two types are the same. Also, because usually the coherence time for the nonrephasing 2D IR spectra is shorter than the corresponding rephasing part, it is experimentally more efficient to only collect the waiting-time-dependent nonrephasing 2D IR spectra. Thus, in this work, we used a series of absolute-valued nonrephasing 2D IR spectra in modeling the IVR processes

Scheme 1. IVR Kinetic Model for a Three-Mode System in Ru3(CO)12 Clustera

a Here, E′r, A2″, and E′a represent E′ (radial), A2″ (axial), and E′ (axial) vibrational states, respectively.

Figure S4 in the SI, in which the 2D frequency regions (2 cm−1× 2 cm−1, thus it is equivalent to peak “intensity” or proportional to the “amplitude” of the positive to negative absorptive 2D IR signal) used to evaluate the diagonal and offdiagonal peak amplitudes are shown in dashed squares. The kinetic rate equations of the IVR processes shown in Scheme 1 and the corresponding solutions are given in Schemes S1 and S2 in the SI, and parameters used in the equations are listed in Table S6. Fitting details are given in Scheme S2. For the Ru3(CO)12 cluster, the IVR process involves three major CO stretching normal modes, as shown in Figure 3, and we assumed an equal vibrational relaxation rate for the involved modes during the fitting. This is a reasonable assumption to begin with because the CO stretching modes of a similar transition-metal carbonyl complex may have similar vibrational relaxation rates.22,50 In addition, considering the excitonic nature of these CO normal modes that are separated by tens of wavenumbers, their lifetimes should be quite similar. On the other hand, a very complicated analytical solution would be needed to solve the kinetic model shown in Scheme 1, if one has to use different energy relaxation time constants for these three involved vibrational modes. Furthermore, our recent study showed that varying the vibrational lifetime (T1) or the relaxation rate constant (kV = 1/T1) by 10−20% does not significant vary the IVR time.37 In addition, because of the presence of the IVR process, traditional broad-band IR pump−probe experiment under the so-called magic-angle polarization condition cannot reliably obtain the T1 values for these modes. This can be easily seen in Figure 3: the time-dependent intensity of a vibrational mode (for example, E′a) probed in the detection frequency position near ωt = 2061.0 cm−1 will contain the contribution of the two offdiagonal peaks (peaks 4 and 6) because a transient pump− probe signal at a specific pump−probe delay time is the same as the vertical summation of a 2D IR spectrum at the same waiting time. Thus, a simplified three-mode kinetic model (Scheme 1) with the same kV parameter was used to fit the diagonal and offdiagonal 2D IR peak amplitudes. In this kinetic model, E′r, A2″, and E′a represent three CO vibrational states, and there are seven variables, including both forward and backward IVR rate constants. The forward and backward rate constants are related to each other through the equilibrium constant, Keq = kbackward/ kforward. Moreover, Keq can be obtained from the transition-state theory, Keq = exp(−hΔω/kBT), where h is Planck’s constant, Δω is the frequency splitting between the two modes involved, 1301

DOI: 10.1021/acs.jpcb.7b10067 J. Phys. Chem. B 2018, 122, 1296−1305

Article

The Journal of Physical Chemistry B kB is the Boltzmann constant, and T is the temperature. The frequency splitting values between the E′ (radial) and A2″ (axial) modes, the A2″ (axial) and E′ (axial) modes, and the E′ (radial) and E′ (axial) modes are 19.5, 32.3, and 51.8 cm−1, respectively, which are obtained from the projection of the 2D IR spectrum in Figure 4. These integrated signals as a function of Tw are displayed in Figure 6.

comparable in intensity to the strong diagonal signal (E′a) at 50 ps waiting time. The IVR times obtained are all on the order of tens of picoseconds, indicating an indeed quite rapid and efficient vibrational energy exchange in this multicarbonyl system. Such efficient IVR is due to several reasons. First, the energy separations among the three normal modes are quite small (the largest gap is 51.8 cm−1), which allow efficient downward and upward energy transfers. Second, these three normal modes are anharmonically coupled; in other words, they are of excitonic nature, in a sense that they are linear combination of 12 local CO stretching motions, as pointed out in Section III.I. This is shown in Table S1 for vibrational coupling and Table S2 for potential energy distributions.47 As can be seen clearly in these results, it is nearly impossible to assign the observed IR peaks to any of the local CO stretching mode. The efficient IVR processes among the three normal modes suggest that once one of the local CO stretching mode is vibrationally excited, its neighboring CO stretching mode will instantaneously sense the excitation too. Under such circumstances, each CO group in this multicentered transition-metal complex is “equally” IR active, which is an important property of such complex when functioning as photocatalysts or in the processes of photoenergy conversion and chemical synthesis.

Figure 6. Waiting-time-dependent amplitudes from the nonrephasing 2D IR spectra of Ru3(CO)12 cluster. (A) Off-diagonal signal between the A2″ (axial) mode and E′ (radial), i.e., peak 1 in Figure 3D; (B) diagonal signal of the A2″ axial mode; (C) off-diagonal signal between the A2″ (axial) mode and E′ (axial), i.e., peak 6 in Figure 3D. Lines are the fits with functions obtained from the kinetic model shown in Scheme 1. Dashed lines are fitting residuals.

IV. CONCLUSIONS AND OUTLOOK In summary, using 2D IR spectroscopy, we directly observe intramolecular vibrational energy-transfer process among different CO stretching modes in the Ru3(CO)12 cluster. This cluster with 12 CO groups has a nonbridged D 3h configuration, forming a quite stable core structure with such a large number of ligands. The linear IR experimental measurements and DFT calculations both revealed three major CO stretching vibrational components, which are the E′ (axial) mode, A2″ (axial) mode, and E′ (radial) mode. Because the latter is weak in its linear IR intensity, 2D IR projection was used to better determine its frequency position. Two major conclusions can be drawn from this work. First, these observed CO stretching modes are of vibrational excitonic nature. This can be seen from their unbalanced absorption intensities (due to anharmonic coupling) and frequency separations (due to both the coupling and local structural differences). Spectral fitting of the selected 2D IR slices leads to a set of diagonal and off-diagonal anharmonicities as a function of waiting time. The finding is that these CO modes in the Ru3(CO)12 cluster have similar diagonal anharmonicities (13.1 cm−1 on average). The average offdiagonal anharmonicity is on the order of 12.9 cm−1, which is very close to their corresponding diagonal anharmonicities. The observed strong and intensity-increasing cross peaks, with respect to the corresponding diagonal peaks, are due to the intramolecular energy transfer, which is the second most important conclusion of this work. Here, intramolecular vibrational excitonic energy redistribution occurs efficiently

Figure 6B shows the amplitude of the diagonal peak (the A2″ (axial) mode), whereas two off-diagonal peaks (peaks 1 and 6) are given in Figure 6A,C, respectively. As is seen, a fast decrease dominates the amplitude, and the amplitudes of these CO stretching modes show similar dynamical traces. The fitting curves and residuals by applying the kinetic model shown in Scheme 1 are also shown in Figure 6. The fitting parameters and the obtained IVR rate constants among the CO vibrational states are listed in Table 4. A shorter IVR time of 11.1 ps is obtained between the E′ (radial) and A2″ (axial) modes, whereas a slightly longer IVR time of 12.7 ps exists between the E′ (radial) and E′ (axial) modes. These IVR processes occur on the time scale of picosecond, as has also been shown in previous studies of other transition-metal carbonyl complexes. 33,37,55 In addition, the vibrational relaxation rate (kV) was determined to be 0.0016 ps−1 (1/ 609.6 ps−1), which is in the time range similar to the reported relaxation time constant (400−600 ps) determined by narrowband pump−broad-band probe infrared method for a similarly structured trinuclear transition-metal carbonyl cluster, i.e., Os3(CO)12.22 The observed long vibrational lifetime allows a broad range of time scale for the intramolecular vibrational energy transfer to occur, which turned to be quite efficient in this system. As already shown in the purely absorptive 2D IR spectra (Figure 3), for example, the off-diagonal signal at peak 6 becomes

Table 4. Peak Splitting (Δω), Equilibrium Constant (Keq), and the Forward and Backward IVR Rate Constants (k) of the CO Vibrational States IVR pathway

Δω/cm−1

Keq

kforward/ps−1

kbackward/ps−1

A2″ (axial) → E′ (radial) E′ (axial) → A2″ (axial) E′ (axial) → E′ (radial)

19.5 32.3 51.8

0.910 0.856 0.779

0.090 (1/11.1) 0.086 (1/11.6) 0.079 (1/12.7)

0.082 (1/12.2) 0.074 (1/13.5) 0.062 (1/16.1)

1302

DOI: 10.1021/acs.jpcb.7b10067 J. Phys. Chem. B 2018, 122, 1296−1305

Article

The Journal of Physical Chemistry B

(3) Adams, R. D.; Barnard, T. S. Catalytic Hydrosilylation of Diarylalkynes by Layer-Segregated Platinum-Ruthenium Cluster Complexes Pt3Ru6(CO)20(μ3-RC2R)(μ3-H)(μ-H). Organometallics 1998, 17, 2567−2573. (4) Li, C.; Widjaja, E.; Garland, M. The Rh4(CO)12-Catalyzed Hydroformylation of 3,3-Dimethylbut-1-Ene Promoted with HMn(CO)5. Bimetallic Catalytic Binuclear Elimination as an Origin for Synergism in Homogeneous Catalysis. J. Am. Chem. Soc. 2003, 125, 5540−5548. (5) Schipper, D. E.; Young, B. E.; Whitmire, K. H. Transformations in Transition-Metal Carbonyls Containing Arsenic: Exploring the Chemistry of [Et4N]2[HAs{Fe(CO)4}3] in the Search for SingleSource Precursors for Advanced Metal Pnictide Materials. Organometallics 2016, 35, 471−483. (6) Bauer, S.; Hunger, C.; Bodensteiner, M.; Ojo, W.-S.; CrosGagneux, A.; Chaudret, B.; Nayral, C.; Delpech, F.; Scheer, M. Transition-Metal Complexes Containing Parent Phosphine or Phosphinyl Ligands and Their Use as Precursors for Phosphide Nanoparticles. Inorg. Chem. 2014, 53, 11438−11446. (7) Gafney, H. D.; Xu, S. P. Photocatalyzed Isomerization of 1Pentene by Ru3(CO)12 Adsorbed onto Porous Vycor Glass. Inorg. Chim. Acta 1995, 240, 645−651. (8) Peschiulli, A.; Smout, V.; Storr, T. E.; Mitchell, E. A.; Elias, Z.; Herrebout, W.; Berthelot, D.; Meerpoel, L.; Maes, B. U. W. Ruthenium-Catalyzed α-(Hetero)Arylation of Saturated Cyclic Amines: Reaction Scope and Mechanism. Chem. - Eur. J. 2013, 19, 10378−10387. (9) Kong, Q.; Lee, J. H.; Kim, K. H.; Kim, J.; Wulff, M.; Ihee, H.; Koch, M. H. J. Ultrafast X-Ray Solution Scattering Reveals Different Reaction Pathways in the Photolysis of Triruthenium Dodecacarbonyl (Ru3(CO)12) after Ultraviolet and Visible Excitation. J. Am. Chem. Soc. 2010, 132, 2600−2607. (10) Glascoe, E. A.; Kling, M. F.; Shanoski, J. E.; Harris, C. B. Nature and Role of Bridged Carbonyl Intermediates in the Ultrafast Photoinduced Rearrangement of Ru3(CO)12. Organometallics 2006, 25, 775−784. (11) Leadbeater, N. E.; Lewis, J.; Raithby, P. R.; Ward, G. N. Photochemistry of [Ru3(CO)12] with Nitrogen Heterocycles. J. Chem. Soc., Dalton Trans. 1997, 2511−2516. (12) Aime, S.; Dastru, W.; Gobetto, R.; Krause, J.; Milone, L. Evaluation of the Energy Barrier for Carbonyl Exchange in the Highly Fluxional Ru3(CO)12 System. Organometallics 1995, 14, 4435−4438. (13) Slebodnick, C.; Zhao, J.; Angel, R.; Hanson, B. E.; Song, Y.; Liu, Z. X.; Hemley, R. J. High Pressure Study of Ru3(CO)12 by X-Ray Diffraction, Raman, and Infrared Spectroscopy. Inorg. Chem. 2004, 43, 5245−5252. (14) Kong, Q.; Lee, J. H.; Plech, A.; Wulff, M.; Ihee, H.; Koch, M. H. J. Ultrafast X-Ray Solution Scattering Reveals an Unknown Reaction Intermediate in the Photolysis of [Ru3(CO)12]. Angew. Chem., Int. Ed. 2008, 47, 5550−5553. (15) Braga, D.; Grepioni, F.; Tedesco, E.; Dyson, P. J.; Martin, C. M.; Johnson, B. F. G. A Variable Temperature Study of Ru3(CO)12 in the Solid State and the Generation of Alternative Crystal Structures. Transition Met. Chem. 1995, 20, 615−624. (16) Anson, C. E.; Jayasooriya, U. A. Vibrational Spectroscopic Investigation of Metal Cluster Prototypes in the Solid State: A Novel Approach to the Carbonyl Stretching Region in Osmium Carbonyl ([Os3(CO)12]) and Ruthenium Carbonyl ([Ru3(CO)12]). Spectrochim. Acta, Part A 1990, 46, 967−974. (17) Battiston, G. A.; Sbrignadello, G.; Bor, G. Infrared Spectroscopic Studies on Metal Carbonyl Compounds. 23. A Simple Quantitative Treatment of the Infrared Band Intensity and the Induced Metal-Metal Dipole Contribution to It in Polynuclear Metal Carbonyls. An Application to the Spectrum of Dodecacarbonyltriruthenium and Dodecacarbonyltriosmium in the Carbon-Oxygen Stretching Region. Inorg. Chem. 1980, 19, 1973−1977. (18) Battiston, G. A.; Bor, G.; Dietler, U. K.; Kettle, S. F. A.; Rossetti, R.; Sbrignadello, G.; Stanghellini, P. L. Comparative Infrared and Raman Spectroscopic ν(CO) Study of Dodecacarbonyltriruthenium,

among these modes. This is manifested by the waiting-timedependent off-diagonal 2D IR peaks and further characterized by a kinetic model, from which the downhill and uphill vibrational energy transfer rates were obtained. The time scales of these energy transfers were found to be around 11−16 ps, whereas the lifetimes of these vibrational excited states are on the order of 400−600 ps.22 Such efficient IVR processes confirm the excitonic nature of the three linear IR absorption bands, in a sense that when some of the CO modes are vibrationally excited, the excitation energy will migrate to other CO modes very quickly. Such an excitonic nature appears to be very crucial to the chemical activity of the CO groups as they function in various ways. We have to point out that in fitting the IVR dynamics, the coupling-induced cross peaks were neglected because this Ru3(CO)12 cluster only shows fluxional property on the time scale of microseconds.54 In addition, we have to point out that in the present work the structural aspect of such trinuclear transition-metal complexes was not explored in detail, including equilibrium structures and their distributions, and solute−solvent interactions, which is also believed to be very important to their functions.



ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jpcb.7b10067. Molecular structure of the Ru3(CO)12 complex (Figure S1); wave function demixing analysis (Table S1); PED analysis (Table S2); representative slice of the 2D IR real-part spectra (Figure S2); determination of diagonal and off-diagonal anharmonicities of the IR-active carbonyl stretching modes at various frequencies (Tables S3− S5); real part of the nonrephasing 2D IR spectrum of Ru3(CO)12 cluster (Figure S3); representative waitingtime-dependent absolute nonrephasing 2D IR spectra of Ru3(CO)12 cluster (Figure S4); kinetic rate equations of the IVR process (Scheme S1); values of peak splitting and equilibrium constant (Table S6); analytical solutions to the kinetic rate equations (Scheme S2) (PDF)



AUTHOR INFORMATION

Corresponding Authors

*E-mail: [email protected]. Tel: (+86)-010-62563565. Fax: (+86)-010-62563167 (F.Y.). *E-mail: [email protected]. Tel: (+86)-010-62656806. Fax: (+86)-010-62563167 (J.W.). ORCID

Jianping Wang: 0000-0001-7127-869X Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS This work was supported by the National Natural Science Foundation of China (21473212 to F.Y., 21603238 to J.Z., and 21573243 to J.W.).



REFERENCES

(1) Hieber, W. Metal Carbonyls, Forty Years of Research. Adv. Organomet. Chem. 1970, 8, 1−28. (2) Powell, H. M.; Ewens, R. V. G. The Crystal Structure of Iron Enneacarbonyl. J. Chem. Soc. 1939, 286−292. 1303

DOI: 10.1021/acs.jpcb.7b10067 J. Phys. Chem. B 2018, 122, 1296−1305

Article

The Journal of Physical Chemistry B Dodecacarbonyltriosmium, Their Mixed Crystals, and the Mixed Triangulo Cluster Carbonyls Dodecacarbonylosmiumdiruthenium and Dodecacarbonyldiosmiumruthenium. Inorg. Chem. 1980, 19, 1961− 1973. (19) Gilson, T. R.; Evans, J. Vibrational Studies of the Cluster Carbonyls of Ruthenium and Osmium. Part 2. Single-Crystal Raman and Infrared Data for [Ru3(CO)12] in the CO Stretching Region. J. Chem. Soc., Dalton Trans. 1984, 155−159. (20) Butler, I. S.; Xu, Z. H.; Darensbourg, D. J.; Pala, M. FT-IR, Photoacoustic and Micro-Raman Spectra of the Dodecacarbonyltriruthenium(0) Complexes Ru3(13CO)12 and Ru3(CO)12. J. Raman Spectrosc. 1987, 18, 357−363. (21) Hunstock, E.; Mealli, C.; Calhorda, M. J.; Reinhold, J. Molecular Structures of M2(CO)9 and M3(CO)12 (M = Fe, Ru, Os): New Theoretical Insights. Inorg. Chem. 1999, 38, 5053−5060. (22) Yan, S.; Seidel, M. T.; Zhang, Z.; Leong, W. K.; Tan, H.-S. Ultrafast Vibrational Relaxation Dynamics of Carbonyl Stretching Modes in Os3(CO)12. J. Chem. Phys. 2011, 135, No. 024501. (23) Sanstead, P. J.; Stevenson, P.; Tokmakoff, A. SequenceDependent Mechanism of DNA Oligonucleotide Dehybridization Resolved through Infrared Spectroscopy. J. Am. Chem. Soc. 2016, 138, 11792−11801. (24) Yang, F.; Yu, P.; Zhao, J.; Shi, J.; Wang, J. Ultrafast Vibrational and Structural Dynamics of Dimeric Cyclopentadienyliron Dicarbonyl Examined by Infrared Spectroscopy. Phys. Chem. Chem. Phys. 2015, 17, 14542−14550. (25) Hithell, G.; Gonzalez-Jimenez, M.; Greetham, G. M.; Donaldson, P. M.; Towrie, M.; Parker, A. W.; Burley, G. A.; Wynne, K.; Hunt, N. T. Ultrafast 2D-IR and Optical Kerr Effect Spectroscopy Reveal the Impact of Duplex Melting on the Structural Dynamics of DNA. Phys. Chem. Chem. Phys. 2017, 19, 10333−10342. (26) Shin, J. Y.; Yamada, S. A.; Fayer, M. D. Dynamics of a Room Temperature Ionic Liquid in Supported Ionic Liquid Membranes vs the Bulk Liquid: 2D IR and Polarized IR Pump-Probe Experiments. J. Am. Chem. Soc. 2017, 139, 311−323. (27) Eckert, P. A.; Kubarych, K. J. Dynamic Flexibility of Hydrogenase Active Site Models Studied with 2D-IR Spectroscopy. J. Phys. Chem. A 2017, 121, 608−615. (28) Ghosh, A.; Ostrander, J. S.; Zanni, M. T. Watching Proteins Wiggle: Mapping Structures with Two Dimensional Infrared Spectroscopy. Chem. Rev. 2017, 117, 10726−10759. (29) Wang, J. Ultrafast Two-Dimensional Infrared Spectroscopy for Molecular Structures and Dynamics with Expanding Wavelength Range and Increasing Sensitivities: From Experimental and Computational Perspectives. Int. Rev. Phys. Chem. 2017, 36, 377−431. (30) Elsaesser, T. Two-Dimensional Infrared Spectroscopy of Intermolecular Hydrogen Bonds in the Condensed Phase. Acc. Chem. Res. 2009, 42, 1220−1228. (31) Perakis, F.; Hamm, P. Two-Dimensional Infrared Spectroscopy of Supercooled Water. J. Phys. Chem. B 2011, 115, 5289−5293. (32) Panman, M. R.; Newton, A. C.; Vos, J.; van den Bosch, B.; Bocokic, V.; Reek, J. N.; Woutersen, S. Ultrafast Dynamics in Iron Tetracarbonyl Olefin Complexes Investigated with Two-Dimensional Vibrational Spectroscopy. Phys. Chem. Chem. Phys. 2013, 15, 1115− 1122. (33) Anna, J. M.; King, J. T.; Kubarych, K. J. Multiple Structures and Dynamics of [CpRu(CO)2]2 and [CpFe(CO)2]2 in Solution Revealed with Two-Dimensional Infrared Spectroscopy. Inorg. Chem. 2011, 50, 9273−9283. (34) Chuntonov, L. 2D-IR Spectroscopy of Hydrogen-BondMediated Vibrational Excitation Transfer. Phys. Chem. Chem. Phys. 2016, 18, 13852−13860. (35) Hithell, G.; Shaw, D. J.; Donaldson, P. M.; Greetham, G. M.; Towrie, M.; Burley, G. A.; Parker, A. W.; Hunt, N. T. Long-Range Vibrational Dynamics Are Directed by Watson-Crick Base Pairing in Duplex DNA. J. Phys. Chem. B 2016, 120, 4009−4018. (36) Yu, P.; Yang, F.; Zhao, J.; Wang, J. Hydration Dynamics of Cyanoferrate Anions Examined by Ultrafast Infrared Spectroscopy. J. Phys. Chem. B 2014, 118, 3104−3114.

(37) Yang, F.; Zhao, J.; Wang, J. Two-Dimensional Infrared Study of C-Natural Abundant Vibrational Transition Reveals Intramolecular Vibrational Redistribution Rather than Fluxional Exchange in Mn(CO)5Br. J. Phys. Chem. B 2016, 120, 1304−1311. (38) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A.,et al. Gaussian 09, revision A.2; Gaussian, Inc.: Wallingford, CT, 2009. (39) Huggins, D. K.; Flitcroft, N.; Kaesz, H. D. Infrared Spectrum of Osmium Tetracarbonyl Trimer, Os3(CO)12; Assignment of CO and MC Stretching Absorptions. Inorg. Chem. 1965, 4, 166−169. (40) Bor, G.; Sbrignadello, G.; Marcati, F. Infrared Spectroscopic Studies on Metal Carbonyl Compounds. XV. Assignments of the C-O Stretching Infrared Spectra and Force Constants of Some (Arene)Co4(CO)9 Complexes. J. Organomet. Chem. 1972, 46, 357−368. (41) El-Sayed, M. A.; Kaesz, H. D. Assignment of the CO Stretching Absorptions in C4v Metal Pentacarbonyl Derivatives. The Method of Local Oscillating Dipoles. J. Mol. Spectrosc. 1962, 9, 310− 316. (42) Anson, C. E.; Jayasooriya, U. A. Vibrational Spectroscopic Investigation of Metal Cluster Prototypes in the Solid State: [M3(CO)12], where M = Os, Ru; Metal-Ligand Stretching and MCO Deformation Region. Spectrochim. Acta, Part A 1990, 46, 861− 869. (43) Zhao, J.; Wang, J. Understanding the Amide-II Vibrations in βPeptides. J. Phys. Chem. B 2015, 119, 14831−14839. (44) Zhao, J.; Wang, J. Dissecting Amide-I Vibration in β-Peptide Helices. J. Phys. Chem. B 2015, 119, 3387−3397. (45) Wang, J.; Hochstrasser, R. M. Anharmonicity of Amide Modes. J. Phys. Chem. B 2006, 110, 3798−3807. (46) Wang, J. Assessment of the Amide-I Local Modes in γ- and βTurns of Peptides. Phys. Chem. Chem. Phys. 2009, 11, 5310−5322. (47) Krimm, S.; Bandekar, J. Vibrational Spectroscopy and Conformation of Peptides, Polypeptides, and Proteins. Adv. Protein Chem. 1986, 38, 181−364. (48) Vergeer, F. W.; Hartl, F.; Matousek, P.; Stufkens, D. J.; Towrie, M. First Direct Observation of a CO-Bridged Primary Photoproduct of [Ru3(CO)12] by Picosecond Time-Resolved IR Spectroscopy. Chem. Commun. 2002, 1220−1221. (49) Adams, D. M.; Taylor, I. D. Solid-State Metal Carbonyls. Part 4. An Infrared and Far-Infrared Study of Triruthenium and Triosmium Dodecacarbonyls. J. Chem. Soc., Faraday Trans. 2 1982, 78, 1561− 1571. (50) Feng, M.; Yang, F.; Wang, J. Vibrational and Structural Dynamics of Mn(CO)5Br and Re(CO)5Br Examined Using Nonlinear Infrared Spectroscopy. Chin. J. Chem. Phys. 2016, 29, 81−86. (51) Baiz, C. R.; McRobbie, P. L.; Anna, J. M.; Geva, E.; Kubarych, K. J. Two-Dimensional Infrared Spectroscopy of Metal Carbonyls. Acc. Chem. Res. 2009, 42, 1395−1404. (52) Cahoon, J. F.; Sawyer, K. R.; Schlegel, J. P.; Harris, C. B. Determining Transition-State Geometries in Liquids Using 2D-IR. Science 2008, 319, 1820−1823. (53) Forster, A.; Johnson, B. F. G.; Lewis, J.; Matheson, T. W.; Robinson, B. H.; Jackson, W. G. 13C Nuclear Magnetic-Resonance Study of Trinuclear Metal-Carbonyls M3(CO)12 (M = Fe, Ru and Os) and Some Related Compounds. J. Chem. Soc., Chem. Commun. 1974, 1042−1044. (54) Aime, S.; Dastru, W.; Gobetto, R.; Krause, J.; Violano, L. Solution Structures and Dynamics of M3(CO)12‑x(NCCH3)x (M = Ru, x = 1, 2, 3; M = Os, x = 1, 2). Inorg. Chim. Acta 1995, 235, 357−366. (55) King, J. T.; Anna, J. M.; Kubarych, K. J. Solvent-Hindered Intramolecular Vibrational Redistribution. Phys. Chem. Chem. Phys. 2011, 13, 5579−5583. (56) Wang, J.; Yang, F.; Zhao, J. Selectively Probing the Structures and Dynamics of β-Peptide Aggregates Using the Amide-A Vibrational Marker. J. Phys. Chem. B 2015, 119, 15451−15459. (57) Cyran, J. D.; Nite, J. M.; Krummel, A. T. Characterizing Anharmonic Vibrational Modes of Quinones with Two-Dimensional Infrared Spectroscopy. J. Phys. Chem. B 2015, 119, 8917−8925. 13

1304

DOI: 10.1021/acs.jpcb.7b10067 J. Phys. Chem. B 2018, 122, 1296−1305

Article

The Journal of Physical Chemistry B (58) Baiz, C. R.; Kubarych, K. J.; Geva, E.; Sibert, E. L. Local-Mode Approach to Modeling Multidimensional Infrared Spectra of Metal Carbonyls. J. Phys. Chem. A 2011, 115, 5354−5363. (59) Wang, J. Conformational Dependence of Anharmonic Vibrations in Peptides: Amide-I Modes in Model Dipeptide. J. Phys. Chem. B 2008, 112, 4790−4800.

1305

DOI: 10.1021/acs.jpcb.7b10067 J. Phys. Chem. B 2018, 122, 1296−1305