M.
Zimmer
a,
F.
Rupp
b,
P.
Singer
b,
F.
Walz
c,
F.
Breher
*c,
W.
Klopper
*d,
R.
Diller
*b and
M.
Gerhards
*a
aChemistry Department and Research Center Optimas, TU Kaiserslautern, Erwin-Schrödinger-Straße 52, 67663 Kaiserslautern, Germany. E-mail: gerhards@chemie.uni-kl.de; Fax: +49 631 2052750; Tel: +49 631 2052537
bPhysics Department, TU Kaiserslautern, Erwin-Schrödinger-Straße 46, 67663 Kaiserslautern, Germany. E-mail: diller@physik.uni-kl.de; Fax: +49 631 2053902; Tel: +49 631 2052323
cInstitute of Inorganic Chemistry, Karlsruhe Institute of Technology (KIT), Engesserstr. 15, 76131 Karlsruhe, Germany. E-mail: breher@kit.edu; Fax: +49 721 60847021; Tel: +49 721 60844855
dInstitute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Fritz-Haber-Weg 2, 76131 Karlsruhe, Germany. E-mail: klopper@kit.edu; Fax: +49 721 60847225; Tel: +49 721 60847263
First published on 1st May 2015
This paper presents a combined spectroscopic and theoretical analysis of a trinuclear [Pd3{Si(mtMe)3}2] complex (mtMe = methimazole) which has been demonstrated to be a potential catalyst for coupling reactions. It is a highly symmetric model system (D3 in the electronic ground state) for the investigation of electronic states and the structure of polynuclear transition metal complexes. Different time-resolved IR spectroscopic methods covering the femtosecond up to the microsecond range as well as density functional computations are performed to unravel the structure and character of this complex in the electronically excited state. These are the first time-resolved IR studies on a trinuclear Pd complex. Based on the interplay between the computational results and those from the IR studies a 3A state is identified as the lowest lying triplet state which has C2 symmetry.
In this context a very interesting trinuclear palladium complex29 (1, Fig. 1) has been investigated by our research groups with a variety of state-of-the-art photophysical and computational methods.30 These analyses gave a detailed insight into the photophysical properties and into the electronically excited states. Based on these investigations we could characterize the lowest lying triplet state from which a long-lived (τAr = 1.3–1.7 μs, τAir = 222 ± 3 ns)30 luminescence occurs. Extensive theoretical work was involved to obtain conclusions on the electronically excited states of 1. However, no conclusive investigations on the structures of the electronic states have been performed up to now, which would yield more precise information on these states and thus a better and more reliable assignment. For the investigations presented here we use structure sensitive time-resolved IR spectroscopy from the femto- to the microsecond time range in combination with further theoretical analyses that go beyond our previous work;30 this is the first application of such a strategy to a trinuclear Pd complex. Our chosen combination of methods is especially suited to analyse the structure of the complex in its ground and excited states.
To cover the nano- and microsecond range, time-resolved step-scan FTIR spectroscopy (TR-FTIR) is applied. Publications of TR-FTIR measurements on transition metal complexes exist since the early 1990s.22–25 Most of those investigations are about electronically excited states, especially with respect to metal-to-ligand charge transfer (MLCT) transitions in complexes including one or two metal centers. In these early investigations so called “reporter” ligands (carbonyl and cyanide) have been introduced which bear the advantage of a huge oscillator strength for the ν(CO) and ν(CN) stretching vibration(s) and a high sensitivity for structural changes.22 Investigated systems without a “reporter” ligand are for example Ru(bpy) complexes (bpy = 2,2′-bipyridine) including the well analyzed [Ru(bpy)3]2+.26,27 Regarding the number of metal atoms very few examples for tri- and tetrametallic investigations with TR-FTIR can be found in the literature, some of them focus on the energy transfer in rhenium complexes, e.g. [(bpy)(CO)3Re(trans-dppene)Re(phen)(CO)2(trans-dppene)Re(bpy)(CO)3]3+ (phen = phenanthroline, dppene = 1,2-bis(diphenylphosphino)ethene) and Re4(CO)12(4,4′-bpy)4Cl4.24,28
Having none of the above mentioned “reporter” groups in our system (1) our investigations are even more demanding. The transient (step-scan) FTIR spectroscopy still benefits from the advantages of the FTIR technique among them the possibility to record a complete spectral region simultaneously (multiplex advantage). Thus it is possible to capture the whole mid-IR region (and in principle also NIR and FIR). In our case the mid-IR region, especially from 1600 to 1100 cm−1, is of interest.
After a detector rise time of about 30 ns our step-scan set up enables us to perform measurements with a resolution of about 10 ns (see Section 2). Beyond this, ultrafast transient IR spectroscopy is applied to the trinuclear palladium complex to investigate the primary photoinduced vibrational dynamics in the femto and picosecond time domain, thereby covering the entire excited state lifetime from its very beginning. Femtosecond transient IR spectroscopy has been applied successfully to the investigation of divers photoinduced phenomena in metal complexes such as CO release,31–33 excited state dynamics in Cr(acac)334 and a Ru(II) catalyst35 as well as high-spin state formation in Fe(II) spin-crossover systems.21 Thereby ultrafast vibrational dynamics can be related to electronic dynamics as intersystem crossing and internal conversion.36,37
A mM solution of 1 in argon purged DMSO-d6 was pumped by a peristaltic pump (Masterflex Easy-Load II, variable flow rate) through an optical cell that is part of a closed-cycle continuous flow tubing system including a sample reservoir. The cell was composed of housing, CaF2 windows, separated by two Mylar® spacers with a variable thickness (0.006 to 0.5 mm). In our case 0.26 mm spacers and a flow rate of about 40 mL min−1 are chosen. The sample cell as well as the sample reservoir were purged by argon to assure an excited state lifetime of τAr = 1.3–1.7 μs, since it is quenched by oxygen to τAir = 222 ns.30 The 355 nm excitation laser beam was attenuated to 0.6–0.8 mJ per shot at a diameter of about 9 mm. The temporal resolution was set to 20 ns for these measurements. Five adjacent 20 ns intervals were joined and averaged, so a temporal resolution of 100 ns was obtained. The step-scan experiment was triggered 400 ns before the Q-Switch of the laser was triggered which means that the laser beam reached the sample 600 ns after the experiment starts. Hence this time was set as zero point in all spectra. A total number of 465 coadditions at each interferogram point was recorded. The spectral region was limited by undersampling to 0–2633 cm−1 with a spectral resolution of 2 cm−1 (resulting in 2665 interferogram points). For the prevention of problems by performing a Fourier transformation an IR longpass filter (transmission 4.73–9.00 μm) was used (i.e. no IR intensity outside the measured region should be observed).
Both amplified AC and DC signals of the MCT detector were fed to a 14-bit transient recorder board (Spectrum Germany, M3I4142, 400 MS/s), which in turn was connected to an external computer running Bruker OPUS software for spectrometer control and data manipulation/analysis. For data analysis the recorded DC and AC spectra were used to calculate the ΔA (change in absorbance) spectra with:
![]() | (1) |
The femtosecond pump–probe experiments were based on a Ti:sapphire regenerative amplifier system to generate ultra-short pump and probe pulses (CPA-2001, Clark MXR: rep. rate 1 kHz, pulse energy 1 mJ, central wavelength 775 nm, pulse duration 150 fs (FWHM)).
A home-build two stage non-collinear optical parametric amplifier was used to generate femtosecond pump pulses at 330 nm with energy of 250–300 nJ at the sample position. Another portion of the fundamental beam was fed in a two stage optical parametric amplifier with consecutive difference frequency generation, yielding the spectrally tunable probe pulse between 6–12 μm with a spectral width of 50–60 cm−1 as described before.38 Pump and probe pulses were superimposed in the sample with focal width of about 200 and 100 μm, respectively. The sample was rotated to exchange the excited volume between two consecutive laser shots. The transmitted mid-IR probe light was spectrally dispersed and detected by a polychromator with a 32-element MCT (mercury–cadmium–telluride) detector array (Infrared-Systems-Development) with ≈2 cm−1 spectral width per element. In all experiments the pump beam was chopped at half the repetition rate (500 Hz), allowing the calculation of pump-induced difference absorbance spectra ΔOD(λ,t) on a shot-to-shot basis as a function of probe wavelength and pump–probe delay. Time zero of the experiment was determined via a transient absorption experiment on a thin silicon wafer. The pump/probe cross correlation time was typically ≈350 fs (FWHM, Gauss).
The transient data was analyzed by a global fit routine, describing the time dependent absorbance change ΔOD(λ,t) as a sum of exponential decay functions
![]() | (2) |
For gaining insight into the observed energy lowering due to geometry relaxation, we computed a correlation diagram (cf. ESI,† Fig. S2) of the frontier orbitals as a function of a distortion coordinate x, which is defined as follows: x = 0 corresponds to an unbiased D3-symmetric geometry that is obtained by averaging the 3A1 and 3A2 geometries optimized in ref. 30; x = −1 corresponds to the C2-symmetric geometry of the 3B state and x = +1 corresponds to the C2-symmetric geometry of the 3A state. Separate linear interpolations of the structures are performed between the x = −1 and x = 0 geometries on the one hand and between the x = 0 and x = +1 geometries on the other. The LUMO is stabilized for both positive (x > 0) and negative (x < 0) distortions. This can be understood as pseudo-Jahn–Teller effect42 due to interactions with a close-lying two-fold degenerate e level (LUMO+1, 72e). At the same time, the HOMO is destabilized towards x = +1 while the HOMO−1 is destabilized towards x = −1. The latter even becomes the HOMO in the 3B geometry. The interplay between destabilization of the relevant ground-state occupied orbitals and stabilization of the LUMO appears to be responsible for the energy lowering upon distortion. Consistent with the correlation diagram, Fig. S3 (ESI†) shows the energies of the excited states of interest as a function of the distortion coordinate, as obtained from TDDFT computations at the B3LYP/TZVP level using the lowest triplet state as reference.
The Cartesian coordinates of the ground state and the two triplet states are given in the ESI† (Tables S1–S3).
![]() | ||
Fig. 3 TR-FTIR spectrum 1.0 μs after excitation; top: in comparison with calculated excitation into a 3A state; bottom: in comparison with calculated excitation into a 3B state. |
In the first 100 ns interval after excitation the bands at 1303/1289, 1385/1373 and 1458/1449 cm−1 are already observable (see inset of Fig. 2). To obtain a better S/N ratio and since no dynamics in this temporal range are visible spectra are also averaged to 500 or 1000 ns steps. The signal disappears completely after about 3.5 μs, confirming the TCSPC lifetime measurements30 and giving a hint that the argon purged sample chamber and argon purged housing of the sample reservoir prevent the sample from oxygen. (For a better illustration of the decay within microseconds, the transient IR spectra up to 3.5 μs are shown in steps of 500 ns in the ESI,†cf. Fig. S7).
In order to compare calculated harmonic frequencies with experimental values the computed vibrational frequencies have to be scaled. A typical empirical factor for the chosen theoretical ansatz (cf. Section 3) is 0.98 (cf. e.g.ref. 43). The reliability of this factor can be tested by fitting the calculated IR spectrum of the electronic ground state (1A) as well as the spectra of the electronically excited states (3A and 3B) to the experimental values. Indeed the best agreement is achieved for a factor of 0.98. (In case of the 3A state the scaling factors vary between 0.979 and 0.983 whereas for the 3B state the variation is larger, between 0.970 and 0.983, indicating a better description of the experimental frequencies for the 3A state.)
In the DFT calculations a red shift of the 3A state vibrations compared to the ground state 1A (modes) in the order of about 10 cm−1 is observable, for example 16 cm−1 for the 1285/1301 cm−1 peaks or 7 cm−1 for the 1373/1380 cm−1 peaks (Fig. S8). In contrast, the vibrational modes of the 3B state are closer to the values of the 1A state. It should be mentioned that the investigated spectral region shows small but significant differences for the calculated 3A and 3B spectra. Other parts of the spectrum are either less significant or cannot be observed with our experimental set-ups (cf. Fig. S9 of the ESI†).
In order to quantify the comparison between calculated and experimentally observed frequencies the calculated spectrum of the electronic ground state (1A) has to be subtracted from the calculated spectra obtained for the 3A and 3B states. In order to simulate the vibrational spectra, Lorentzian functions are used. By fixing the band positions (scaling factor 0.98), different full widths at half maximum (FWHM) are chosen. By testing a constant width or variable widths (between 10 and 20 cm−1) within one fit it turns out that 14 cm−1 is the best average value. The goodness of the fit can be described by the R2/χ2 parameter which is about a factor of three to four better for the 3A state compared to 3B (by using a scaling factor of 0.98 and a constant band width of 14 cm−1 or a variable band width between 10 and 20 cm−1 optimized for each transition). Fig. 3 shows the comparison of experiment and theory: Fig. 3 (top) ground state (1A) subtracted from 3A and Fig. 3 (bottom) 1A subtracted from 3B. It can be seen that the calculated band positions (relative to the experimental values) are better predicted for the 3A state, especially in the case of the transitions at 1289, 1373 and 1449 cm−1. In case of the 3B state also an additional band around 1310 cm−1 would be expected which is not observed in the experimental spectrum. Furthermore, the calculated difference spectra reveal also that only the 3A state gives an appropriate description of the experimentally observed band intensity ratios, especially for the transitions at 1373/1385 cm−1. In case of this most prominent pair of transitions both increase and depletion in the ΔA spectrum are of similar intensity for the 3A state (in agreement with the experiment) whereas for the 3B state the increase should be by a factor of 3 larger than the depletion. Taking the qualitative and quantitative aspects of calculated frequencies and intensities into account, it can be concluded that the excited state probed on the nano- and microsecond time scale is a 3A state.
Another experiment was performed with a frequency doubled Nd:YLF laser (edgewave INNOSLAB IS8 II – ET, output of about 0.4 mJ per shot at 261.5 nm) that was operated at 100 Hz. According to our previous experimental study30 it is known that the UV/Vis spectrum shows two strong and broad absorptions: the first electronic transition has its maximum around 333 nm and a second one is around 262 nm. By exciting the second electronic transition with a wavelength of 261.5 nm, the same IR transitions as for the excitation of the first electronic transitions (with 355 nm) are observed in the step-scan FTIR spectrum. These finding clearly show that after electronic excitation the same long-lived electronic state is reached.
The stationary difference spectrum at 340 ps exhibits negative peaks at 1176, 1303, 1386, 1402 and 1460 cm−1 and positive peaks at 1164, 1291, 1375 and 1450 cm−1. As in the TR-FTIR spectra the bleach peaks (neg.) are assigned to the electronic ground state and the positive vibrational bands are assigned to the long living excited state, the latter showing a rather uniform bathochromic shift of about 10 cm−1 with respect to the former. The difference spectrum at early times (0.5 ps) is very similar to that at 340 ps, however, the respective difference bands are of larger amplitude with increased positive absorbance on the low-energy side. This feature is easily explained by vibrational cooling, i.e. the dissipation of the excess energy to the solvent: Anharmonic coupling between thermally populated low-energy vibrational modes and high-energy modes (fingerprint region) causes asymmetrically broadened vibrational bands. Thus, the observed apparent blue-shift on the picosecond time scale indicates the cooling process of the vibrationally hot (unrelaxed) molecule in the excited electronic state. It should be noted that the spectral shapes (widths) of the positive bands in the transient spectra are not uniform since they are determined by band-specific parameters such as (1) spectral width (dephasing time) of the particular band in the excited electronic state and (2) the anharmonic coupling constants of the particular band to all other vibrational modes in the excited electronic state, especially to the low frequency modes.
Although the relation between vibrational temperature and time-dependency of the observed spectral changes is intrinsically non-exponential, the global fit according to eqn (2) is used in order to derive a characteristic time constant. This yields τVR = 7.8 ± 0.4 ps with the decay associated spectrum A1(λ) given in Fig. 5. It now describes the observed blue-shift of each product vibrational band (Fig. 4) in terms of monoexponentially decaying positive (low energy lobe) and negative (high energy lobe) absorbance strength as an approximation for the cooling process. A0(λ) is basically identical with the 340 ps difference spectrum, confirming the quality of the fit. Note the excellent match of A0(λ) and the microsecond-FTIR spectra (Fig. 3). No further, sub-ps time constant (as found in the transient UV/Vis data reported earlier (380 fs)30) could be derived from the data. This might be due to the highly transitory character of the Franck–Condon state, leading to spectrally extremely broad infrared bands with possibly low absorption cross section. In addition the given (limited) time resolution and S/N-ratio add to this effect. Thus the IR data together with the transient UV/Vis data reported earlier (ref. 30) indicate that ultrafast intersystem-crossing (ISC), internal conversion (IC) and thermalization via intramolecular vibrational relaxation (IVR) lead to the vibrationally hot 3A excited state within less than 400 fs. Subsequently the excess energy of about 15500 cm−1 (according to λexc = 330 nm and λFluor = 675 nm30) is dissipated into the solvent.
![]() | ||
Fig. 5 Decay associated mid-IR spectra (DAS, according to eqn (2)), describing the cooling process (A1(λ)) and the “stationary” difference absorption spectrum (A0(λ)) (cf. text for details). |
The ultrafast dynamics, i.e. ISC within a few hundred femtoseconds followed by vibrational cooling on the picosecond time scale is typical for photoexcited transition-metal complexes36 and has been frequently reported.21,35 In addition, it turns out that the femtosecond transient IR spectra, by their strongly structured spectra, allow for a much clearer identification of the cooling/relaxation process as compared to broad and unstructured UV/Vis spectra, where corresponding spectral broadening and shifts are often much weaker or not even detectable as in our earlier work30 on this complex.
Furthermore the femtosecond transient IR data reveal vibrational cooling characterized by a monoexponential time constant of 7.8 ps and no further changes up to 1 ns. Thus the vibrational excess energy dissipates on a picosecond time scale.
The presented extensive spectroscopic and theoretical studies on a highly symmetric trinuclear Pd complex allow an assignment of electronically excited states and its dynamical behaviour. A long-lived electronically excited state is formed by interaction of the three Pd atoms. The geometry of the complex changes from D3 symmetry in the electronic ground state to C2 symmetry in the 3A excited state, i.e. a pseudo-Jahn–Teller effect is observed. The trinuclear Pd complex has proved to be a rewarding object to study photophysics of a oligonuclear metal complex by the interplay of experiment and theory. The application of a variety of spectroscopic methods (UV/Vis, luminescence, molecular beam)30 including now the structure sensitive IR techniques in combination with high level DFT and TDDFT calculations provide a detailed analysis on the structural changes by going from the electronic ground to an electronically excited state. These investigations are a starting point for further examples of interacting metal centers which can be of photophysical, magnetic or catalytical relevance.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c5cp00959f |
This journal is © the Owner Societies 2015 |