Grigory
Smolentsev
*a,
Alexander A.
Guda
b,
Markus
Janousch
a,
Cristophe
Frieh
a,
Gaudenz
Jud
a,
Flavio
Zamponi
c,
Murielle
Chavarot-Kerlidou
d,
Vincent
Artero
d,
Jeroen A.
van Bokhoven
ac and
Maarten
Nachtegaal
a
aPaul Scherrer Institut, WLGA 217, 5232, Villigen, Switzerland. E-mail: grigory.smolentsev@psi.ch; Tel: +41 56 310 5173
bResearch Center for Nanoscale Structure of Matter, Southern Federal University, Sorge 5, 344090, Rostov-on-Don, Russia
cETH Zurich, Wolfgang Pauli-Str. 10, 8093, Zurich, Switzerland
dLaboratory of Chemistry and Biology of Metals, Universite Grenoble Alpes, CEA, CNRS, France
First published on 14th April 2014
In order to probe the structure of reaction intermediates of photochemical reactions a new setup for laser-initiated time-resolved X-ray absorption (XAS) measurements has been developed. With this approach the arrival time of each photon in respect to the laser pulse is measured and therefore full kinetic information is obtained. All X-rays that reach the detector are used to measure this kinetic information and therefore the detection efficiency of this method is high. The newly developed setup is optimized for time-resolved experiments in the microsecond range for samples with relatively low metal concentration (∼1mM). This setup has been applied to study a multicomponent photocatalytic system with a Co(dmgBF2)2 catalyst (dmg2− = dimethylglyoximato dianion), [Ru(bpy)3]2+ chromophore (bpy = 2,2′-bipyridine) and methyl viologen as the electron relay. On the basis of the analysis of hundreds of Co K-edge XAS spectra corresponding to different delay times after the laser excitation of the chromophore, the presence of a Co(I) intermediate is confirmed. The calculated X-ray transient signal for a model of Co(I) state with a 0.14 Å displacement of Co out of the dmg ligand plane and with the closest solvent molecule at a distance of 2.06 Å gives reasonable agreement with the experimental data.
Scheme 1 Structure of Co(dmgBF2)2 (R = methyl) and Co(dpgBF2)2 (R = phenyl). Sol is a solvent molecule (CH3CN). |
The key requirements to the techniques used to identify the photocatalytic intermediates are sensitivity to the structural and electronic changes of the catalytic center, selectivity that allows to distinguish the contribution of the catalyst from all other possible changes in the complicated chemical system and appropriate time range and resolution. Stopped flow and freeze quench methods are too slow to track the intermediates of H2 evolving catalysts, especially for the most promising systems that produce hydrogen rapidly. Furthermore, since most of the reaction steps are intermolecular, the intermediates have lifetimes in the microsecond time range. Some initial reaction steps, for example, charge transfer from the chromophore to the catalytic center in the supramolecular systems can be intramolecular and occur within picoseconds,19 but the charge-separated state has to be long lived to participate in the subsequent intermolecular reactions. Thus measurements in the microsecond range with sub-microsecond resolution are required for such applications.
Among time-resolved spectroscopic methods that are typically applied to study photocatalytic intermediates and reaction mechanism are transient optical absorption spectroscopy20,21 and time-resolved IR.22 Time-resolved X-ray absorption spectroscopy (XAS) has a few advantages in comparison with optical and IR methods. First of all XAS is element-specific and therefore the contribution from the metal center of the catalyst is separated from all other changes in the system. Second, XAS spectra contain structural information about the local arrangement of the metal (in the fine structure of the near edge X-ray absorption spectrum known as XANES) and information about the oxidation state (mainly in the position of the absorption edge). Setups for time-resolved hard X-ray absorption spectroscopy measurements with ∼100 ps resolution have been developed at a few beamlines worldwide.23–28 Experiments in the microsecond range are quite rare and can be performed using spatially separated laser and X-ray beams (pump–flow–probe mode)29,30 gating of the detector (pump–probe mode)31 or recording the arrival time of all photons (pump–sequential-probes mode).29,32–34
Pump–flow–probe experiments use non-overlapping laser excitation and X-ray probing beams focused on the sample jet. The solution is circulating in the flow system and a cylindrical jet with a liquid speed of a few meters per second is formed. A continuous wave laser focused on the jet initiates the reaction, the distance between laser and X-ray beams and the velocity of the liquid flow determine the time delay, while the focusing of both beams defines the time resolution. The pump–flow–probe method is characterized by a high detection efficiency and the resolution is limited to ∼40 μs with easily achievable flow parameters (flow speed ∼5m s−1, jet diameter ∼1mm) and moderate X-ray spot size (∼100 μm). Pump–probe experiments with a microsecond resolution can use the same data acquisition electronics as pump–probe experiments with picosecond resolution by increasing the width of the X-ray detector gate pulse. Since many X-ray bunches delivered by the synchrotron will fall within the time window defined by such a gate pulse, the time resolution is defined by the width of this pulse. The detection efficiency of pump–probe methods also depends on the gate, its shortening corresponds to higher time resolution but lower efficiency. Acquiring kinetic data with this approach requires scanning the delay between laser and gate pulses.
In the present manuscript we report a pump–sequential-probes setup to study photocatalytic intermediates that have lifetimes in the microsecond range. This setup is realized at the SuperXAS beamline of the Swiss Light Source and described in the experimental setup section. With this approach the arrival time of each photon in respect to the laser pulse is measured (method known as time-tagged photon counting) and therefore, in contrast to currently available pump–probe setups, full kinetic information is obtained without any delay scan. In the application section we demonstrate how this method is used to identify the intermediate states of a Co(dmgBF2)2 catalyst in the multicomponent system and to probe the structure of the Co(II) and Co(I) states of the catalyst in solution.
For our source, the maximal focusing of the X-ray beam is ∼100 ×100 μm2. The laser spot with diameter D has to be bigger in order to probe only the volume excited efficiently and taking into account that it is difficult to achieve a stable flow of the sample in the liquid jet with a speed v higher than ∼5m s−1, the maximal repetition rate f = v/D ∼ 30 kHz is required. In our setup we used a Xiton IDOL laser with a repetition rate up to 40 kHz and maximal output power ∼2W at 15kHz. Available wavelengths for this laser are 447 nm and 671 nm. Maximal pulse energy is 125 μJ which is more than enough to excite samples of reasonable optical density (0.2 or higher) with a metal concentration of ∼1 mM in an irradiated volume of ∼200 μm diameter and ∼1 mm thickness.
To achieve synchronization between X-rays, laser and DAQ we used the radio frequency (RF) signal from the storage ring (500 MHz), which is intrinsically synchronized with the X-ray pulses. It serves as the input for the event receiver (Micro Research, VME-EVR-230RF) that divides the frequency of the RF signal and generates square pulses corresponding to the required repetition rate of the experiment. These signals are further processed with the delay card (Micro Research, 4CHTIM-200) generating pulses of required duration (2 μs for the laser trigger, 40ns for the digital X-ray processor trigger) and with required relative delay. The signal sent to the laser allows to control its acousto-optic modulator and thus to achieve synchronization between the actual laser pulse and DAQ better than 2 ns. HDF5 files generated by XMAP are processed online during the incoming X-ray energy scan with a short delay (time of measurement of 1–2 energy points). Histograms with the number of events as a function of time for selected energies of fluorescent photons are calculated. Further processing such us calculation of transient spectra is performed on the basis of these data.
After each laser pulse intensities I0 and If are measured as a function of time. A spectrum with good statistics can be obtained by averaging X-ray data following at least 105 laser pulses. However if the time position of the laser pulse is fixed to the camshaft position, the I0 as a function of time would exhibit dips in intensity, as shown in Fig. 2 by I0 synch. Instead, in the asynchronous mode the trigger pulses follow with the frequency so that the filling pattern of the storage ring is shifted relative to the laser pulse by 2 ns (1/480 of the storage ring cycle) after each laser period. As a result measuring each point of the spectrum for at least 10 s (corresponding to 104 laser cycles) one observes a uniform intensity distribution for the incoming beam, as shown in Fig. 2 by I0 asynch. It allows to avoid “bad time points” that correspond to the measurements of the kinetic histogram during the time between multibunches of the storage ring. Additionally, we eliminated the need of incoming intensity measurements for each time point I0(t). Only the average value I0 has to be measured. In the synchronized mode fluorescence intensity has to be normalized to the incoming intensity for each time point since bunches in the storage ring have a different current. Simplification of I0 detection in the asynchronous mode has a positive impact on the signal to noise ratio of measured XAS spectra since the contribution of noise from I0 measurements becomes negligible.
[Ru(bpy)3]2+* + MV2+ → [Ru(bpy)3]3+ + MV+ |
[Ru(bpy)3]2+* + Co(II) → [Ru(bpy)3]3+ + Co(I) |
MV+ + Co(II) → Co(I) + MV2+ |
[Ru(bpy)3]3+ + Co(II) → [Ru(bpy)3]2+ + Co(III) |
According to the time-resolved optical data reported for Co(dpgBF2)221 the catalyst in the Co(III) state is formed in parallel with Co reduction at the conditions without any proton source. In most Co-based H2-evolution systems,5 Co(I) species play a crucial role in the catalytic mechanism: protonation of Co(I) yields Co-hydride intermediates which are the actual active species involved in H2 evolution, either in the Co(III)-H state or in the Co(II)-H state obtained after a subsequent reducing event.
Our first goal has been to determine how many Co intermediates formed in the system. In the bottom panel of Fig. 3 we show the transient signal corresponding to different delay times between laser pump and registered X-ray photons. The width of the time window for the data corresponding to the long delay is significantly larger than for other two spectra to compensate the reduction of statistics due to the signal reduction, one can also notice that the data measured for the (0.25, 0.75) μs delay is noisier than the spectrum corresponding to a six times wider (−0.25, 2.75) μs time window. The amplitude of the transients has been re-normalized for easier comparison. The time-dependence of the signal monitored at a fixed X-ray energy corresponding to the maximum of the transient signal is shown in the top panel of Fig. 3. Since a liquid jet with continuous flow has been used to refresh the sample between two laser pulses the slow component of the decay can be caused by the movement of the excited volume out of the probed area, while the quick decay of the signal during the first 20 μs is related to the real kinetics of the system. Kinetic measurements utilizing the flow system are thus not as precise as those that one can imagine in the experiment with the full exchange of the sample during the pause before each laser pulse (the repetition rate of such measurements will be extremely low). Nevertheless, the experimental spectra (bottom panel of Fig. 3) measured at different time delays have very similar shapes indicating the formation of only one novel species upon irradiation of the photocatalytic system. The same result has been obtained from the analysis of the shape of the transient signals as a function of time over a set of 200 experimental spectra using the method reported in ref. 37, which indicated the presence of only 2 independent spectral components in our data corresponding to the initial species together with one intermediate species.
Fig. 4 shows the experimental XANES spectrum of the photocatalytic system measured before the laser pulse together with the transient spectrum that corresponds to the 10 μs time window after the laser pulse. The position of the absorption edge of the intermediate species is shifted towards lower energies with regards to the initial Co(II) (first peak of the transient signal is positive), indicative of the formation of a more reduced species, likely the Co(I) intermediate mentioned above and expected to be the entry point into a catalytic cycle for H2 evolution.
The local structure around the metal center for molecules in solution can differ from that in the crystal phase measured using XRD and from that obtained by DFT optimization in vacuum. This is especially true for metal complexes possessing vacant coordination sites or labile ligands. Co(dmgBF2)2 belongs to this class of metal complexes and therefore we have started the simulations from the refinement of the structure of the initial Co(II) state. Both XRD38 and DFT17 structures of [Co(dmgBF2)2(CH3CN)2] have rather distant solvent molecules with bond lengths of 2.25 Å and 2.33 Å correspondingly, while some other Co(II) complexes with distorted octahedral environment and acetonitrile in the first coordination shell have significantly shorter Co–NSol bonds (2.08–2.09 Å) in the crystallized state.39,40 Therefore in addition to XRD and DFT structures we have constructed models (3) and (4) with short Co–NSol bonds (2.06 Å and 2.0 Å respectively). To take into account the expected disorder of the solvent molecules and avoid biases due to the strong multiple scattering along the Co–N–C path as found in the crystal and DFT-based structures, we have approximated the contribution of solvent to the XANES by taking into account the scattering from the nearest N atoms only.
Fig. 5 shows that the agreement between the experimental spectrum and the spectra calculated for XRD (1) and DFT (2) models is not satisfactory. The relative intensity of peaks B and C is wrong. In addition, both calculated spectra display a shoulder between minimum D and maximum E that is not so well pronounced in the experimental data. Both models with short Co–NSol bonds (3 and 4) agree reasonably well with the experiment. Variation of the Co–NSol distance within these limits influences slightly the intensity of shoulder A which is observable, but rather weak in the experimental spectrum. As an additional small adjustment we have reduced the Co–Neq distances for models (3) and (4) by 0.03 Å relative to the XRD model. It has not changed the overall shape of the spectrum but has allowed to achieve a better agreement in the position of minimum D.
Fig. 5 Theoretical Co K-edge XANES spectra of [CoII(dmgBF2)2(CH3CN)2] calculated for a model obtained using either X-ray diffraction38 (red curve (1)), DFT optimization17 (blue curve (2)), a model with disordered solvent molecules with Co–NSol distances of 2.06 Å (black curve (3)) and 2.00 Å (magenta curve (4)). Curve (5) is the experimental spectrum of the multicomponent photocatalytic system with [Co(dmgBF2)2(CH3CN)2] catalyst measured before the laser pulse. The insert shows the structure of [Co(dmgBF2)2(CH3CN)2]. Grey atoms are C, red O, blue N, green B, yellow F, magenta Co. H atoms are omitted for simplicity. |
Experimentally determined reference structures for Co(I) species with glyoxime ligands are very rare in the literature. The only available crystal structure of a Co(I) complex with 2 diphenylglyoximato ligands and acetonitrile has been reported in ref. 12. Since the atomic cluster with radius 4.9 Å has been used to calculate X-ray absorption spectra, the substitution of methyl by phenyl group does not influence significantly the shape of XANES. According to this XRD model Co(I) is five-coordinated in a distorted square-pyramidal environment with a short 1.97 Å Co–NSol bond and 0.27 Å displacement of Co out of the plane formed by the equatorial nitrogen atoms. Acetonitrile is oriented almost perpendicularly to the equatorial plane with a Co–N–C angle of 174°. Hu et al. note that “solvent molecules were severely disordered and could not be refined accurately”.12 The DFT optimized structure for Co(I) state reported by Muckerman and Fujita17 has a significantly different orientation of the solvent molecules. Two of them were included in this model and the authors have found that both of them are oriented by the CH3 group towards the Co atom with large Co–C distances of 3.7 Å and 4.8 Å, thus the Co atom is 4-coordinated.
Model (3) of the Co(I) intermediate has been constructed starting from the model of the Co(II) species with disordered solvents (curve (3) in Fig. 5) that gave a good agreement with the experiment. One of the solvent molecules has been moved away from the Co center to the distance 2.18 Å, while the second CH3CN has been fixed at 2.06 Å. Additionally we have introduced a moderate out of plane displacement of the Co atom (0.14 Å) and a small contraction of Co–Neq (0.02 Å) as it has been suggested in both XRD and DFT models. Fig. 6 compares theoretical transient spectra calculated for these three models with the experimental pump–probe data. In addition to the structural changes we have taken into account 1 eV chemical shift of the absorption edge which is due to the movement of the core level as a result of changes of the screening by valence electrons. As one can see the shape of the transient signal is completely wrong for XRD model (1). The reason is similar to those observed also for the XRD structure of the ground state Co(II) species. Strong multiple scattering from the collinear Co–N–C chain of the solvent molecule influences the XANES significantly and therefore the model that ignores the disorder of solvent molecule fixing the solvent at idealized orthogonal orientation is not realistic. The DFT model is in better agreement with the experiment. A correct overall trend of the transient spectrum is seen, but all the main features are shifted in energy. Please note that even if the absolute scale of the energy calculations is not always correct we have performed already a few eV shift of the ground state spectrum and therefore there is no additional free parameter related with the independent shift of the transient spectrum. Among the 3 considered structures, model (3) gives the best agreement with the experiment. The theoretical spectrum is more structured between maxima A and B, which can be partially due to the fact that for the ground state peak A is more visible in the theory then in the experimental spectrum. Further improvement of the theoretical model can be based on more accurate modeling of disordered solvent around the metal complex by averaging many configurations of solvent molecules, for example using the approach previously proposed in ref. 41. Nevertheless, on the basis of our data we can conclude that the model with disordered solvent molecules with the shortest Co–NSol distance ∼2.06 Å and moderate displacement of Co out of the plane formed by N atoms of dmg ligands is the most realistic while the formation of a complex with perpendicular orientation of acetonitrile molecules has not been confirmed.
Discussing the experimental results shown in Fig. 4 that are plotted with the error bars we would like to mentioned that these measurements were taken during ∼11 h and as one can see we have obtained a transient signal of very high quality. Please note that in comparison to the supramolecular systems,42,43 or simple chromophorores,44–46 the pump–probe signal for multicomponent systems is ∼10 times weaker even if we are not limited by the pulse energy of the initial photoexcitation. Nevertheless, intermediates are relatively long-lived which somehow compensates the decrease of the statistics from the reduction of the signal amplitude. The incoming intensity of X-ray radiation at the bending magnet beamline of SLS is lower than, for example, at the undulator beamlines of the APS synchrotron. Paradoxically it can be seen also as an advantage since lower intensity allows to increase the detection efficiency using the single photon counting approach and at the same time to decrease the X-ray induced damage of the sample. In our previous work on the multicomponent system with similar Co-based catalyst that has been performed at 11ID-D beamline of the APS we observed a significant damage of the sample due to X-rays (incident X-ray flux ∼6 × 1012 photons s−1) while it was negligible at the SuperXAS beamline of the SLS (flux 3 × 1011 photons s−1).29 The number of photons registered by the fluorescence detector at the APS is around 3–4 per X-ray bunch for the samples with ∼1mM concentration; therefore single photon counting is not possible and the data acquisition system works in the current mode. At SuperXAS of SLS for analogous sample we have a bit less than 1 photon μs−1 and therefore single photon counting is efficient. An additional advantage of our approach that has not been used for the present application of the setup is the possibility to use energy-resolving detectors, such as silicon drift detectors, that have a resolution of ∼120 eV and therefore allow for a further increase of the signal to noise ratio for diluted samples.
This journal is © The Royal Society of Chemistry 2014 |