Ultrafast photo-induced ligand solvolysis of cis-[Ru(bipyridine)2(nicotinamide)2]2+: experimental and theoretical insight into its photoactivation mechanism

Mechanistic insight into the photo-induced solvent substitution reaction of cis-[Ru(bipyridine)2(nicotinamide)2](2+) (1) is presented. Complex 1 is a photoactive species, designed to display high cytotoxicity following irradiation, for potential use in photodynamic therapy (photochemotherapy). In Ru(II) complexes of this type, efficient population of a dissociative triplet metal-centred ((3)MC) state is key to generating high quantum yields of a penta-coordinate intermediate (PCI) species, which in turn may form the target species: a mono-aqua photoproduct [Ru(bipyridine)2(nicotinamide)(H2O)](2+) (2). Following irradiation of 1, a thorough kinetic picture is derived from ultrafast UV/Vis transient absorption spectroscopy measurements, using a 'target analysis' approach, and provides both timescales and quantum yields for the key processes involved. We show that photoactivation of 1 to 2 occurs with a quantum yield ≥0.36, all within a timeframe of ~400 ps. Characterization of the excited states involved, particularly the nature of the PCI and how it undergoes a geometry relaxation to accommodate the water ligand, which is a keystone in the efficiency of the photoactivation of 1, is accomplished through state-of-the-art computation including complete active space self-consistent field methods and time-dependent density functional theory. Importantly, the conclusions here provide a detailed understanding of the initial stages involved in this photoactivation and the foundation required for designing more efficacious photochemotherapy drugs of this type.


Introduction
The unique photochemical properties of ruthenium-pyridyl complexes have been exploited for a plethora of technological and medicinal applications, such as: light-harvesting, 1,2 lightemitting devices, [3][4][5] fluorescence imaging 6,7 and cytotoxic action. 7,8ne of the more novel uses of these complexes is as photoactivatable prodrugs, where light induces formation of an active species from a stable, inert precursor molecule; a technique known by the term photodynamic therapy (PDT), 9 or more generally as photochemotherapy.0][21] For all these complexes, there is only a limited understanding of the photochemistry and photophysics that underpins their transformation from the inert to the active species.Perhaps the most heavily studied inorganic chromophores are polypyridyl complexes of ruthenium(II), of which, the archetypal tris-(2,2 0bipyridine)ruthenium(II), or [Ru(bpy) 3 ] 2+ , is a member. 224][25] The 'averaged' 3 MLCT state therefore involves electron transfer from a metal d-orbital to the p* orbital of one of the bipyridine (bpy) ligands and is formally This journal is © the Owner Societies 2014 termed [Ru III (bpy) 2 (bpy À )] 2+ .At a similar energy to this 3 MLCT state, lies the lowest energy triplet metal-centred ( 3 MC, d-d ligand field) state, the accessibility of which is dependent on solvent and temperature. 26The 3 MC state is known to exhibit significant distortions of the geometry relative to the ground state or that of the MLCT states, 27 and therefore an increased ligand labilization.For fully bidentate-chelated complexes such as [Ru(bpy) 3 ] 2+ , this dissociation channel is not usually a high yielding relaxation pathway 28,29 and such complexes primarily decay radiatively back to the ground state from the 3 MLCT and/or 3 MC states. 25When one of the bpy ligands is replaced, for example, by two monodentate pyridines (or related aromatic heterocycles, L) to form a complex of the class [Ru(bpy) 2 (L 0 )(L 00 )] 2+ , a solvent substitution reaction may ensue, depending on the nature of the solvent. 30As such, the replacement of an L ligand with a solvent molecule (solvolysis), following photodetachment, may provide a further excited state deactivation mode.The relative ordering of the energies of the 3 MLCT and 3 MC states becomes variable depending on the polar and steric effects of the substituent ligand, L. 26 It has been postulated that in order to create complexes with high photoproduct quantum yields (i.e.those which undergo ligand dissociation), the 3 MC state must be made the lowest energy excited state or thermally accessible from the 3 MLCT state. 262][33][34] Such investigations may provide fundamental understanding of the mode(s) of activation, highlighting inefficiencies in the photoreaction of these complexes and so aid the design of more efficacious complexes, ultimately delivering improved healthcare technologies.Turro and co-workers have studied the ultrafast dynamics of the bipyridyl ruthenium acetonitrile complex cis-[Ru(bpy) 2 (CH 3 CN) 2 ] 2+ following photoexcitation at 310 nm, in both water and acetonitrile solutions. 32For aqueous transient UV/Vis absorption experiments at early times (o10 ps), typical absorption features of the 3 MLCT excited state and an absorption, assigned to a penta-coordinated Ru(II) intermediate (PCI), i.e. [Ru(bpy) 2 (CH 3 CN)] 2+ , were observed.Due to the near simultaneous appearance of these features, and the relatively slow decay (B50 ps) of the 3 MLCT signature, the authors proposed a kinetic picture of direct formation of the PCI from the initially populated singlet states (MLCT and/or pp* ligand centred).The PCI may then recombine with its original partner ligand or, especially in the case of aqueous solutions, undergo solvolysis.Importantly, Turro and co-workers observed an absorption peak at 458 nm, appearing after 10 ps, which they assigned to the monoaqua complex, i.e., [Ru(bpy) 2 (CH 3 CN)(H 2 O)] 2+ . 32ore recently, Lamberti and co-workers have carried out a combined experimental and theoretical study, using transient X-ray and visible absorption spectroscopies, and time-dependent density functional theory calculations to unravel the 3 MLCT excited state structure of cis-[Ru(bpy) 2 (py) 2 ] 2+ (py = pyridine). 34heir results showed that there is little reorganization of the geometry from the ground state to the 3 MLCT state; only a small shortening of the Ru-N(bpy) bonds.Calculations of the 3 MC state, however, showed a significant lengthening of one of the Ru-N(py) bonds and minimal re-optimizations of the bpy ligands.The transient X-ray absorption data provided experimental verification of the 3 MLCT structure, but experimental determination of the 3 MC geometries remains a challenge.The small intensity and energy differences in the transient visible absorptions between the ground state and mono-aquated photoproduct, cis-[Ru(bpy) 2 (py)(H 2 O)] 2+ , made accurate determination of the dissociation time of the pyridine ligand difficult, with an estimated dissociation time of o3 ns after excitation.
Depending on the size of the labile ligands and the nature of the extended solvent structure, it can be expected that there will be an ensemble of PCI species with their geminate partners having various degrees of separation. 35These nascent photodetached ligands and partnering PCI complexes, surrounded by a solvation shell, can be classed as caged species and may either geminately reform the original complex or undergo (some degree) of diffusional separation, with the relative propensities for these processes depending on the exact nature of the system.In systems that undergo separation (or where the ligand is initially ejected through the solvation shell, e.g.smaller ligands), photodetached ligands may in principle still re-encounter the PCI complex, but after a longer diffusional period, to afford diffusive recombination.To date, no work has placed an emphasis on understanding the structural changes surrounding the PCI of a [Ru(bpy) 2 (L 0 )(L 00 )] 2+ type complex during ligand solvolysis and the impact of caging on the degree of solvolysis (i.e. the quantum yield for formation of the monoaquated species), which will undoubtedly contribute to its efficacy as a pro-drug.
To this end, in this work we have investigated the structural and photochemical properties of the novel ruthenium complex, cis-[Ru(bpy) 2 (NA) 2 ] 2+ (1) (bpy = 2,2 0 -bipyridine and NA = nicotinamide, pyridine-3-carboxamide), the structure of which is shown in Fig. 1 inset and is termed 1 hereon.7][38] Given the high levels of cytotoxicity shown in similar Ru(II) anticancer complexes in the recent literature, 7,21,[38][39][40][41][42] 1 is likely to provide high cytotoxicity once activated.In aqueous solution, upon irradiation with UV or blue light, 1 undergoes solvolysis to form the mono-aquated photoproduct [Ru(bpy) 2 (NA)(H 2 O)] 2+ (2).Continued irradiation (ca. 1 h) leads to the formation of the bis-aquated photoproduct [Ru(bpy) 2 (H 2 O) 2 ] 2+ (3).Here we report the first detailed investigation into the excited state dynamics of the transformation of 1 to 2 (the transformation of 2 to 3 is much slower and beyond the time-window of our experiment) using a combination of ultrafast UV/Vis transient absorption spectroscopy (TAS) in conjunction with detailed theoretical calculations using density functional theory and its time-dependent analogue, in addition to the complete active space self-consistent field method.

Synthesis and characterization
NA was purchased from Sigma Aldrich and used as received.
[Ru(bpy) 2 Cl 2 ] was prepared following the literature procedure outlined in ref. 43.1Á(PF 6 ) 2 , was synthesized in an analogous manner to other ruthenium(II) complexes of the type [Ru(bpy) 2 (L) 2 ](PF 6 ) 2 . 14[Ru(bpy) 2 Cl 2 ] (104 mg, 0.2 mmol) was added to 10 mL degassed water and the solution was heated in the dark at 353 K for 15 min.NA (122 mg, 1 mmol) was added and the solution was heated for a further 2 h.The resulting dark red solution was allowed to cool, and any insoluble material was removed by filtration.NH 4 PF 6 (163 mg, 1 mmol) was added to give an orange precipitate.The solution was cooled for 2 h (ca.277 K) and filtered under suction.

Ultrafast transient UV/Vis absorption spectroscopy
The UV/Vis TAS setup has been described previously. 44Femtosecond laser pulses are produced by a commercial Ti-Sapphire regenerative amplified laser system (Spectra-Physics, Spitfire XP), which outputs 40 fs, 3 mJ pulses centred at 800 nm with a 1 kHz repetition rate.The beam is split into three parts of equal intensity.Two of these 1 W beams are used in the transient absorption setup: one beam is split into 95% and 5% parts for harmonic generation (not used in these experiments) and white light continuum (WLC) generation, respectively; the second beam is used to pump an optical parametric amplifier (OPA) (Light Conversion, TOPAS-C (UV-VIS)).The OPA provides the pump source with tuneable wavelengths in the range 240-1160 nm and typical powers of 8-90 mJ per pulse.A neutral density filter in the pump path provides variable pump fluences.The WLC probe is generated, by focusing the weak 800 nm beam into a 1 mm thick calcium fluoride (CaF 2 ) window.A variable aperture and a neutral density filter control the focusing and power of the 800 nm beam.A piezo actuator continuously translates the CaF 2 window vertically, in order to prevent WLC instabilities and damage to the CaF 2 . 45The WLC is recollimated using a spherical mirror ( f = 50 mm) and is passed through a colour filter (Hoya CM-500) to remove the 800 nm fundamental and to reduce the intensity of the red component of the WLC.The usable WLC probe, spanning a range of 340-700 nm, is focused into the sample by a second spherical mirror ( f = 200 mm).A hollow gold retroreflector (Edmund Optics) mounted on a motorized translation stage (Physik-Instrumente, M-505.4DG) in the pre-WLC path creates a variable pump-probe time delay ranging from femtoseconds (fs) to picoseconds (ps).Use of optics on flip mounts and prepositioned magnetic bases gives extra delay lines for select longer, nanosecond (ns) pump-probe time delays.
The pump beam is focused 10 mm behind the sample by a CaF 2 lens ( f = 500 mm) to ensure that the diameter of the pump at the sample is greater than that of the probe.The pump and probe beams intersect the sample (see below) with a crossing angle of B41.Rotation of a l/2 waveplate in the pre-WLC path allows control of the polarization of the WLC; this is set at magic angle (54.71) 46 relative to the pump polarization.Probe pulses at each pump-probe time delay are recorded using a fibre-coupled spectrometer (Avantes, AvaSpec-ULS1650F).The spectrometer captures at 1 kHz with 1000 shots per frame-grab.A mechanical chopper in the pump beam blocks every-other pulse so that the detector records 'sample-pumped' and 'sample-not-pumped' spectra sequentially.Data acquisition and the calculation of difference spectra, as well as control of the translation stage, are performed using a purpose built LabVIEW program.
Sample, at room temperature (295 K), is delivered using a 1 mm path length flow through quartz cuvette (Starna UK) and a PTFE tubing peristaltic pump (Masterflex).A homebuilt shutter based around the design in ref. 47 is used to block both pump and probe beams when not recording data to reduce sample degradation.This is also controlled through the LabVIEW program.

Theoretical calculations
The ground state geometries of complexes 1, 2 and the [Ru(bpy) 2 (NA)] 2+ PCI were optimized using density functional theory (DFT).Analytical Hessian evaluation confirmed the nature of the optimized geometries as minima.Several different functionals and basis sets were compared (PBE0, 48 M06L, 49 CAM-B3LYP 50 and B3LYP 51,52 ).In particular, regarding basis sets the use of non-, quasi-, and fully relativistic effective core potentials (ECPs) was investigated in conjunction with the cc-pVTZ valence basis on the metal and non-metals, respectively. 53It was found that such basis-set effects on geometry and electronic spectroscopy (vide infra) were relatively minor compared to the choice of functional.As such, a non-relativistic This journal is © the Owner Societies 2014 SDD 28 electron (1s 2 2s 2 2p 6 3s 2 3p 6 3d 10 ) ECP was used for the Ru centre, while the cc-pVTZ basis set 54 was selected for the remaining C, N, O and H atoms.The electronic spectroscopy (singlet and triplet excitations) of 1, 2 and the PCI was subsequently modelled with time-dependent density functional theory (TD-DFT), using the same basis set.For the functional, by far the best comparison with experimental UV/Vis spectra was found to be the B3LYP functional.A polarizable continuum model (PCM) was used to model the effects of solvation in H 2 O.
To investigate the excited state reaction pathways in the PCI species, which determine the subsequent relaxation mechanism following light absorption, complete active space self-consistent field (CASSCF) calculations were performed.6][57] Alternative schemes including quasi-natural orbital approaches, such as those discussed in ref. 58 and 59,  were also performed at selected geometries to ensure the active space description was balanced and stable.Generally, these two schemes gave very similar results although the quasi-natural orbital approach was more problematic both in terms of convergence and characterization of the metal centred components, due to initial strong orbital mixing.Given the size of the systems under investigation, and the associated computational expense, the basis sets used for CASSCF calculations were a nonrelativistic SDD 28 electron (1s 2 2s 2 2p 6 3s 2 3p 6 3d 10 ) ECP basis for the ruthenium; the 4-31G basis for carbon, oxygen, and nitrogen; and the STO-3G basis for hydrogen.This basis was calibrated (at the ground state geometry of the PCI) against a larger basis consisting of the same metal centred part, plus the 6-31G(d) basis for the non-metals.Geometrical effects in using the larger basis were very small.For calculations in the singlet and triplet manifolds many-electron configuration state function (CSF) bases were used, corresponding to hS ˆzi = 0, hS ˆ2i = 0 and hS ˆzi = 1, hS ˆ2i = 2, respectively.For calculations of singlet/triplet interactions a Slater determinant many-electron basis with hS ˆzi = 0 was used.The orbital rotation derivatives were neglected in solving the coupled-perturbed multi-configurational self-consistent field (MCSCF) equations in geometry optimizations involving state-averaged orbitals.Gaussian09 60 was used for DFT and TD-DFT calculations, while Gaussian03 61 was used for CASSCF calculations.

Static UV/Vis absorption spectra
UV/Vis absorption spectra of 1 in aqueous solution (40 mM) before, and following excitation at 465 nm using a 50 mW blue LED, are shown in Fig. 1.Complex 1 has a broad visible absorption maximum, l max = 420-450 nm.Following 50 seconds of irradiation, the ground state absorption of 1 decayed and only the absorption of 2 was present with l max = 465 nm.The identity of the mono-aqua adduct 2 was confirmed by MS data shown in the ESI † (Fig. S1).Furthermore, HPLC analysis indicated that only one isomer of 2 was present after irradiation, in-line with previous literature, which reported that only the cis-2, and not trans-2, isomer will be generated. 30We discuss the lack of isomerisation (and its origins) further in Section 3.6.A further 1 h 50 min of irradiation was required to convert the entire sample to the bis-aqua complex, 3, further red-shifting the l max to 483 nm.Power dependence studies of a similar complex (cis-[Ru(bpy) 2 (CH 3 CN) 2 ] 2+ ) and the times of conversion indicated that the stepwise ligand exchange has two separate photochemical mechanisms. 32n order to benchmark our complementary TD-DFT calculations against experiments, the calculated UV/Vis absorption spectra for complexes 1 and 2, using the TD-B3LYP method, are presented in Fig. 2a and b, respectively.As alluded to in Section 2.3, the B3LYP functional together with an SDD 28 electron ECP (Ru) and a cc-pVTZ (H, C, O and N) basis set, was selected based on the best comparison with experimental findings.For complex 1, Fig. 2a shows very good agreement between the calculated spectra and the experimental spectrum.The TD-B3LYP calculations on 1 were performed both on the isolated gas phase species, as well as in the presence of an H 2 O PCM, in an attempt to capture some affects (if any) of solvation on the singlet electronic transitions.The results of both of these calculations are given in Fig. 2a, and show that for 1, negligible changes to the energies of singlet electronic excitations are observed between the gas phase and H 2 O PCM TD-B3LYP calculations, although an overall enhancement in the transition dipole moments (TDMs) for the dominant 'bright' transitions is observed (cf.oscillator strengths, f, given in Table 1), the strongest of which belong to excitations of 1 MLCT states at 428 and 408 nm.
Results of analogous calculations for the mono-aqua photoproduct complex 2 are also presented in Fig. 2b.These generate a UV/Vis absorption spectrum which is blue shifted by B30 nm relative to experiment, for both the gas phase and H 2 O PCM, although in good qualitative agreement with the experimentally observed profile.Once again, the major absorption signatures arise from bright 1 MLCT transitions, as with complex 1 (see Table 1 and ESI, † Fig. S2, for orbital transitions).Most notably, experimental findings clearly show an absorption red shift of B30 nm upon evolution from the starting complex 1 to photoproduct 2. Computationally, this red shift is only qualitatively captured upon inclusion of the H 2 O PCM, while gas phase calculations predicted that 1 and 2 absorb at similar wavelengths.However, even with the inclusion of the H 2 O PCM, the calculations only predict a red shift of B5-10 nm between 1 and 2 at this level of theory (N.B. similar shifts were observed with other functionals, e.g.PBE0).There may be several reasons for this discrepancy.(i) This may be a direct effect of solvation, and more in-depth solvation models (such as hybrid QM/MM methods 62,63 ), although beyond the scope of the present work, may be necessary to capture this behaviour fully.(ii) The TD-B3LYP calculations predict a number of singlet electronic states to lie below the bright 1 MLCT states (in both 1 and 2), which have minimal transition strengths from the equilibrium ground states of 1 and 2. Vibronic coupling however, may lend transition strength to excitations through thermal population of low frequency vibrations modes in the ground state, which are not captured in these vertical excitation calculations.][66][67] This spin-orbit coupling will lend transition strength for direct excitation to these (formally forbidden) 3 MLCT states, and may also be responsible for contributing to the observed differences between the experimental and computational UV/Vis spectra.

Ultrafast transient UV/Vis absorption spectra
Ultrafast UV/Vis TAS was used to follow the dynamics of an 890 mM aqueous solution of 1 flowing through a flow cell, following excitation with a 340 nm, 650 mW pump pulse.Fig. 3 shows TAS at pump-probe time delays ranging from t = À1 ps to 1.2 ns.We begin by focusing on the spectra recorded at early pump-probe delays shown in Fig. 3a (t o 25 ps).Inspection of these spectra shows that several distinct regions can be identified immediately following excitation, related to the spectral location.A strong bleach signal is observed, centred on 420 nm (feature ii), which closely matches the spectral profile of the steady state absorption of 1 (Fig. 1, orange line).As such, feature ii is assigned to the ground state bleach (GSB) signal of 1. Relative to the absorption profile of 1 in Fig. 1, a narrowing of the GSB feature ii is also observed at the 'wings' of the bleach (B400 and B475 nm).This is due to overlap with positive DOD signal either side of this feature, where DOD denotes the change in the optical density (absorbance).The strongest of these adjacent features, labelled i in Fig. 3a, is centred at B370 nm.UV/Vis spectroelectrochemistry measurements of [Ru(bpy) 3 ] 2+ 68 and related complexes with functionalized bpy derivatives 67 indicate that feature i can be assigned to an excited Table 1 Calculated vertical excitation energies (D E), transition wavelengths and associated oscillator strengths (f) for significant triplet and singlet electronic excitation transitions in complexes 1 and 2, using TD-B3LYP with an SDD 28 electron ECP (Ru) and cc-pVTZ (H, C, O and N) basis set.Major molecular orbital transitions and their relative contributions (%) to a given state are also provided (L = LUMO and H = HOMO) together with vertical energy differences (dE) between significant states in 1 (see Fig. 7 for orbital transitions in 1).Central wavelengths for the experimentally observed absorption bands of 1 and 2, assigned predominantly to the 'bright' 1 MLCT state(s), are also given (cf.Fig. 1) state absorption (ESA) of the 3 MLCT excited state, and specifically corresponds to an absorption from the bpy anion (bpy À ) present within the formally charge-separated character of the 3 MLCT state (i.e.[Ru III (bpy)(bpy À )(NA) 2 ] 2+ ). 69Further comparison with earlier literature also indicates that there should be a second (less strong) ESA absorption signal associated with the 3 MLCT state within the region B450-550 nm. 70Between 450-500 nm any such signal will be strongly convoluted with the GSB feature ii.At l 4 500 nm however, a broad plateau of transient absorption signal is observed.We therefore suggest that the absorption between 500-550 nm is primarily due to the ESA of the 3 MLCT state.At much longer wavelengths (l 4 550 nm), another absorbing species must be contributing to this plateau of transient absorption signal, labelled feature iii.Once again, comparison with earlier studies can prove instructive here.Very recently, further UV/Vis TAS measurements on [Ru(bpy) 3 ] 2+ and its methylated derivatives by Hauser and co-workers, 71 have suggested that absorption in this region of the TAS is due to ESA of the 3 MC state.In these bidentate-ligand chelated complexes, at larger Ru-N(bpy) bond distances the 3 MC surface will exhibit notable bound character, as complete Ru-N(bpy) bond fission is necessarily aborted (see Fig. 4 in ref. 71), enabling population to become trapped transiently in the 3 MC state, prior to population transfer back to the ground state.This behaviour will extend the lifetime of population in the 3 MC state and make it spectroscopically detectable in the TAS measurements.However, the analogous 3 MC state in 1 will be purely repulsive with respect to bond fission, meaning that any population transferred to the 3 MC state will undergo rapid (ultrafast) dissociation into a PCI complex and NA ligand, rather than being transiently trapped on this surface.Calculated profiles of the 3 MC state along the Ru-N(py) dissociation coordinate in the related cis-[Ru(bpy) 2 (py) 2 ] 2+ complex support this picture. 31Such a scenario would make it unlikely that any fleeting population in the 3 MC state is observable in 1.In earlier TAS work on cis-[Ru(bpy) 2 (CH 3 CN) 2 ] 2+ by Turro and co-workers, a similar signal was instead assigned to absorption of a PCI complex, although limited justification was provided for this assignment. 32To add weight to this postulate here, we have also performed TD-B3LYP calculations in an attempt to predict the absorption profile of the PCI complex, the results of which are shown in Fig. 2c.Unlike complex 1 and the mono-aquated photoproduct 2 (Fig. 2a and b, respectively), the calculated absorption profiles for the PCI (in the gas phase and H 2 O PCM) indeed qualitatively suggest that: (i) its absorption cross-section is weaker than 1 and 2; and (ii) it should absorb weakly into the 'red end' of our experimental detection window in the TAS (l 4 550 nm).Given this, we are therefore inclined to assign feature iii to absorbance of the PCI complex, [Ru(bpy) 2 (NA)] 2+ , in concord to earlier work by Turro and co-workers. 32The presence of feature iii in the TAS at the earliest pump-probe delays (t = 0.5 ps) also confirms that any Ru-NA bond fission occurs on a sub-picosecond timeframe.
Having assigned the main features in Fig. 3a, we now consider any observed temporal evolution at early time-delays (t o 25 ps).In general, all three labelled features in Fig. 3a begin to deplete within the first 25 ps.Perhaps most significant though is the evolution of the 3 MLCT ESA feature i, which both begins to deplete, spectrally narrow and blue-shift over this timeframe.Such behaviour can be broadly attributed to vibrational cooling of the nascent 3 MLCT state occurring over the timescale of a few picoseconds and we return to consider this behaviour in greater detail below.Similar depletion of the signal associated with the PCI feature iii may also indicate rapid cooling of any nascent PCI species, or alternatively may be due to a structural rearrangement of the PCI after it is formed -we examine this further in Section 3.6.Finally, the apparent depletion of the GSB feature ii is most likely attributable to the vibrational cooling dynamics in the spectrally overlapped 3 MLCT ESA signal, rather than any rapid repopulation of the ground state of 1.
Fig. 3b follows the dynamics of these features to longer time delays (25 ps o t o 1.2 ns).The GSB feature ii recovers whilst the 3 MLCT ESA feature i and the PCI feature iii concomitantly decay over the following few hundred picoseconds and a new absorption feature centred at 475 nm (labelled feature iv) emerges at t 4 300 ps, reaching a maximum intensity after B500 ps.The growth of feature iv results in a quasi-isosbestic point at 490 nm as this new feature rises and the 3 MLCT ESA and PCI absorption features at l 4 500 nm decay over this time window.We attribute this absorption to mono-aquation of the PCI and the formation of photoproduct 2, as the location of feature iv is concordant with the static UV/Vis absorption spectrum of 2 in Fig. 1 (cyan line).The formation of a photoproduct is in-line with the fact that the GSB feature ii never fully recovers back to a baseline signal of DOD = 0, confirming that some portion of initially photo-excited population in 1 never returns to the ground state.Further confirmation of the assignment of feature iv to the formation of 2 can be seen when TAS measurements on 1 are recorded from an acetone solution (ESI, † Fig. S3); acetone is a weakly coordinating ligand compared to H 2 O.In this comparative data set (i) no obvious transient absorption signal of a ligand substituted photo-product is present within our probe window of 340-700 nm and (ii) the GSB feature fully recovers to a baseline of DOD = 0, suggesting an absence of any significant photoproduct yield (even if its absorption signature were to lie outside our probe window).
In earlier ultrafast transient absorption studies of related Ru(II) complexes, a more quantitative insight into the timescales for 3 MLCT state decay (and vibrational cooling), GSB recovery, decay of the PCI species and the formation of the mono-aquated photoproduct was attempted by integrating 'slices' through different spectral regions of the recorded TAS and then fitting the decay/growth of these traces to either mono-or bi-exponential functions. 32However, given that many of the spectral signatures for each of these individual species are heavily convoluted within the recorded TAS (vide supra), such an analysis can return misleading timescales and quantum yields (f) for the dissociation dynamics (e.g.Fig. 1 clearly shows that the spectral signatures for 1 and 2 are strongly overlapped).In an attempt to circumvent this issue, we elect instead to perform a 'target analysis' of the TAS by fitting each of the time-resolved spectra with appropriate basis functions (associated with each of the absorption/bleach features) using the spectral analysis package KOALA (see ESI, † Section S5 for further details). 72Integration of each of these fitted basis functions as a function of time then returns kinetic traces that solely correlate to the population dynamics of individual species involved in the dissociation of 1.
Representative fit examples, obtained for time-resolved spectra recorded at t = 1 ps and 1 ns, are given in Fig. 4a and b, respectively.The static UV/Vis absorption spectra recorded for 1 and 2 in Fig. 1 are used as basis functions for modelling the evolution of the GSB of 1 and the formation of photoproduct 2 (orange and cyan lines, respectively).][69][70] Finally, the spectral profile of the PCI species present at long wavelengths is modelled with a broad Gaussian function (green line).In the case of the PCI, we acknowledge that, like 1 and 2, TD-B3LYP calculations predict this species should also absorb in the blue end of our detection window (l o 450 nm), albeit less strongly (see Fig. 2c).However, attempts to include this behaviour make the target analysis far more cumbersome, due to the multiple absorbing species in this spectral region.To make the analysis feasible, we therefore elect to only model the PCI absorption at l 4 550 nm.The representative fits in Fig. 4 serve to highlight the convoluted nature of the final TAS.This is particularly clear at longer time-delays (Fig. 4b, t = 1 ns) where the remaining GSB signal for 1 and the photoproduct 2 are directly overlapped at l o 450 nm, and the sum of the basis function profiles generated in the fit (grey line) replicates the experimental spectrum (black dashed line) exceptionally well.

Kinetic model for dissociation of cis-[Ru(bpy) 2 (NA) 2 ] 2+
Integration of the fitted basis functions from the target analysis for the 3 MLCT ESA, GSB recovery of 1, the PCI and the growth of photoproduct 2 are presented in Fig. 5a-d.Before proceeding with any fitting of these traces, we first establish a full kinetic model for the dissociation of 1 by considering the dynamical pathways that give rise to Ru-NA bond fission, as well as any (dominant) competing relaxation processes.We herein discuss these with reference to the scheme in Fig. 6a.
The predominant models for photo-induced ligand exchange/ dissociation in similar Ru(II) complexes in the literature, propose extremely efficient traversal of excited state flux from the vertical Franck-Condon excitation region of the 1 MLCT state(s) to the PCI via 3 MLCT state(s) and a dissociative 3 MC state.Following ISC and geometry rearrangements, the PCI can then accommodate the addition of an H 2 O molecule to yield 2. To help construct such a model for 1, and to gain some insight into the relative energies of the 1 MLCT, 3 MLCT and 3 MC states postulated to drive ligand dissociation, we once again turn to our complementary TD-DFT calculations.The results of these TD-B3LYP calculations are summarized in Table 1 and show an optically bright 1 MLCT centred around 2.90 eV (428 nm), in both the gas phase and with an H 2 O PCM.This state is the central component of the 'triply degenerate' state (split by only B0.15 eV) that one would expect in a d 6  This journal is © the Owner Societies 2014 complex with metal-p interactions.Analysis of the (dominant) orbital transition associated with this state confirms it to be of MLCT type, with the relevant orbital transition shown in Fig. 7a (LUMO ' HOMO À 2).A higher energy, but more weakly absorbing, 1 MLCT component is identified to lie at 3.04 eV (408 nm), with the associated orbital transition for this state shown in Fig. 7b (LUMO + 1 ' HOMO À 2); for both of these 1 MLCT states the associated orbital transitions further highlight their formally charge separated [Ru III (bpy)(bpy À )(NA) 2 ] 2+ character (vide supra).Based on the oscillator strengths, f, returned by the TD-B3LYP calculations, the absorption of photons with energy around 2.90 eV or higher, viz.l excitation = 340 nm (4.76 eV), will most likely result in the initial population of highly electronically excited 1 MLCT states, rather than 1 MC states (see Table 1), followed by ultrafast internal conversion to the lower energy 1 MLCT states.We note for completeness that the TAS pump wavelength of 340 nm is selected since it lies outside our probe window and thus avoids contaminating part of our TAS spectra.

metal
The associated 3 MLCT state (of analogous LUMO ' HOMO À 2 character) is calculated to lie at B2.70 eV (B458 nm); B0.2 eV lower in energy than 1 MLCT (see Table 1).Similarly, the 3 MLCT state dominated by the LUMO + 1 ' HOMO À 2 transition lies lower in energy again at B2.65 eV (B468 nm).The high-density of states will promote very efficient vibronic coupling via multiple Jahn-Teller (JT) and pseudo-JT type couplings between these 1 MLCT and 3 MLCT potential energy surfaces. 57For the archetypal complex, [Ru(bpy) 3 ] 2+ , 1 MLCT -3 MLCT ISC is reported to: (i) possess a (near) unity quantum yield (f B 1); 73 and (ii) be ultrafast (B100 fs 24,[64][65][66][67] ), leading to a non-thermally equilibrated (vibrationally excited/hot) 3 MLCT state, henceforth termed 3 MLCT v40 in Fig. 6a. 74McGarvey and co-workers 75 found there to be remarkably little difference in 1 MLCT -3 MLCT v40 coupling within the [Ru(bpy) 3 ] 2+ complex and other similar complexes, where similar sub-picosecond population of 3 MLCT v40 state(s) is also observed.Such behaviour is also consistent with observations from our TAS spectra of 1 in Fig. 3, which indicate that population is present in the 3 MLCT state at our earliest recorded time-delay (t = 0.5 ps), and Fig. 5 (a-d) Kinetic traces for the time-dependent evolution of the 3 MLCT state population (blue), PCI (green), photoproduct 2 (cyan) and GSB recovery of 1 (orange), obtained by integration of the basis functions (see Fig. 4) used in the 'target analysis' of the TAS in Fig. 3. Kinetic fits to these traces, in accordance with the scheme in Fig. 6a, are shown by the solid lines -see main text for details.(e) Kinetic trace reflecting the vibrational cooling of the 3 MLCT state at t o 25 ps (grey squares), obtained by integration of the 'red edge' of the 3 MLCT feature i in Fig. 3a over a range l = 375-380 nm.The solid grey line represents a fit to the trace with a bi-exponential decay function.Once populated, intramolecular and/or intermolecular (solute-solvent mediated) vibrational relaxation (collectively termed vibrational energy transfer (VET) 76 ) of 3 MLCT v40 may subsequently occur with a lifetime t VET , leading to the formation of a vibrationally-cold 3 MLCT state, termed 3 MLCT v=0 , as indicated in the schematic in Fig. 6a.Once again, this is consistent with our observations for band narrowing of the 3 MLCT ESA feature i at t o 25 ps in Fig. 3a.Integrating the signal on the 'red edge' (l = 375-380 nm) of the 3 MLCT ESA feature i in Fig. 3a yields the kinetic trace in Fig. 5e, which when fitted to a bi-exponential decay returns a time-constant of 3.6 ps for t VET (see Table 2).Once cooled, the population located in 3 MLCT v=0 may then undergo relaxation (via either ISC or phosphorescence) back to the singlet electronic ground state of 1 (ground state recovery, GSR) with a time constant t GSR , which based on observations in Fig. 3b and the kinetic trace for 3 MLCT in Fig. 5a, must occur over a longer timeframe of hundreds of picoseconds.
In addition to the VET and recovery of the ground state of 1 from the 3 MLCT state, population in this state must also couple onto a neighbouring dissociative 3 MC surface to drive Ru-NA bond fission. 30Our TD-B3LYP calculations in Table 1 predict the lowest energy of these 3 MC states to lie between B2.9-3.0 eV (425-413 nm), B0.25 eV higher in energy than the lower lying 3 MLCT (LUMO ' HOMO À 2) state in the vertical Franck-Condon region.The dominant orbital transition associated with this 3 MC state (LUMO + 10 ' HOMO, see Fig. 7c) populates an antibonding MC orbital which has strong 4d x 2 Ày 2 character and is dissociative with respect to both of the Ru-NA dative bonds (as well as the equatorial Ru-N(bpy) bonds).For completeness, the analogous1 MC state is calculated to lie at B3.6 eV (B345 nm).We recall that, as with the 3 MLCT ESA feature i in Fig. 3, a transient absorption signature for the nascent PCI species, formed as a result of Ru-NA dissociation, is also observed at the earliest time-delay of t = 0.5 ps in the TAS, indicating that Ru-NA bond fission must also occur on a sub-500 fs timeframe, with a time-constant t d .Based on simple kinetics grounds, if coupling onto the 3 MC surface were barrierless (i.e.accessible from both the 3 MLCT v40 and 3 MLCT v=0 levels), the observed lifetime of the 3 MLCT state should in turn be solely defined by the rapid dissociation lifetime t d .However, this picture is clearly not consistent with general observations from Fig. 3 and the kinetic trace for 3 MLCT in Fig. 5a, which decays over hundreds of picoseconds, rather than a sub-500 fs timeframe.Given this, we invoke a simple model where a barrier to dissociation (coupling onto the 3 MC state) can only be surmounted from a vibrationally hot 3 MLCT v40 state; the presence of a barrier to coupling between 3 MLCT and 3 MC states is generally supported by the fact that our TD-B3LYP calculations predict the 3 MC to lie above the 3 MLCT in the vertical Franck-Condon region. 25Within this model, VET out of the 3 MLCT v40 level over t VET = 3.6 ps acts as a 'time dependent switch', preventing population transfer between 3 MLCT and 3 MC states as the process becomes energetically inaccessible (a barriered process) after population is funnelled into 3 MLCT v=0 .Fig. 6a therefore shows that population subsequently trapped in 3 MLCT v=0 following VET may then only undergo relaxation back to the ground state of 1. 27 Immediately after dissociation, the size of both the PCI complex and NA ligand will mean that they undergo strong caging within the solvation shell, and can be viewed as 'loosely associated', as illustrated pictorially in Fig. 6b.These nascent caged photoproducts, henceforth labelled PCIÁNA (Fig. 6), may then geminately reform the initial complex 1 with a timeconstant t GR , or alternatively undergo some degree of diffusional separation on a timescale t s .The later of these processes will yield the formation of a PCI complex that may then accept a free H 2 O solvent molecule into its now open coordination site, simply termed 'separated' PCI in Fig. 6, and undergo monoaquation to form the final photoproduct 2 on a timescale t MA .
The above kinetic picture of dissociation is broadly consistent with phenomenological models used to treat VET, bond fission, caging and diffusional separation/escape dynamics observed in ultrafast photodissociation studies of smaller organic species. 77,78However, before proceeding further we first briefly justify some of the inherent assumptions surrounding the aforementioned picture, for the dissociation of 1 in H 2 O. (i) After dissociation the 'caged' PCIÁNA pair are in close proximity and as a result we make the assumption that direct formation of 2 from this species will be strongly hindered on steric grounds -hence this pathway is absent in Fig. 6. (ii) After some degree of diffusional separation, the 'separated' PCI and NA species may in principle re-encounter one another, leading to a 'diffusive recombination' pathway for the reformation of 1, rather than the production of 2. However, given the vast excess of H 2 O solvent molecules, and their relative lability compared to the bulkier PCI and NA, we propose that t MA will strongly out-compete any timescale for diffusive recombination, t DR (i.e.t DR c t MA ).As such, t DR is neglected in our kinetic model for the dissociation of 1 in H 2 O in Fig. 6a, and is represented as a kinetically inactive channel for the post-dissociation processes illustrated in Fig. 6b, which we justify further below.

Kinetic analysis
Using this kinetic model for the photodissociation of 1 in H 2 O, full kinetic fits to the traces were performed.These traces, presented in Fig. 5a-d, were simultaneously fit to the analytical solutions for the time-dependent populations of each of these species within our model in Fig. 6a; these first-order differential kinetic equations are detailed in the ESI † online (Section S5).We note that during this fitting procedure, the t VET timescale of 3.6 ps, extracted earlier from Fig. 5e, is held fixed.Given that our kinetic model in Fig. 6a displays branched kinetic pathways, fixing t VET helps to make the fitting procedure more reliable (and less cumbersome).The results of these fits are presented as solid lines through the data sets in Fig. 5a-d, and the time-constants returned are collated in Table 2.
For the case of the 3 MLCT feature in Fig. 5a, this trace is fitted to both the kinetic equations for 3 MLCT v40 and 3 MLCT v=0 , given that they will both contribute to this feature.
Similarly, the PCI feature in Fig. 4b is fitted with the equations for both the 'caged' PCIÁNA and 'separated' PCI species in our model, as they will both possess the same spectral signature.In this case though, the amplitude of the 'separated' PCI component required to fit this trace is found to be negligible.This is in-line with the time-constants returned by the kinetic fits, which indicate that the steady-state population of the 'separated' PCI complex will be minimal, as it is rapidly consumed to generate the monoaquated photoproduct 2, with a timescale of t MA o 1 ps.
More generally the results of this analysis confirm that Ru-NA bond fission, via coupling from 3 MLCT v40 onto the dissociative 3 MC state, does indeed occur on a sub-picosecond timeframe (t d = 0.4 ps).After dissociation into 'caged' PCIÁNA species, geminate reformation of 1 or diffusive separation, forming 'separated' PCI complexes, occur.Both are strongly competitive processes, taking place on the order of a few hundred picoseconds (t GR = 263 ps and t s = 377 ps).For the fraction of PCIÁNA that do go on to form 'separated' PCI complexes, the PCI may then rapidly accommodate an H 2 O solvent molecule into the now open coordination site (t MA o 1 ps), and our analysis indicates that it is the diffusional separation/ reorganization process from 'caged' PCIÁNA to 'separated' PCI that acts as the rate limiting step to the formation of photoproduct 2 (i.e.t MA { t s ) -hence the observed build-up of 2 over hundreds of picoseconds in the TAS.Moreover, the hundreds of picoseconds timeframe extracted for diffusive separation, further serves to confirm our assertion from Section 3.3 that diffusive recombination, t DR , of the separated PCI and NA in aqueous solution (which is also a diffusion limited process) must also occur over a similar timeframe post-separation, and as such can be negated from our kinetic model in Fig. 6, as t MA { t DR .In tandem with these dissociative pathways, a fraction of the population initially present in 3 MLCT v40 will be vibrationally quenched into a cold 3 MLCT v=0 state (t VET = 3.6 ps), and from here undergo relaxation back to the ground state of 1 with a time constant of t GSR = 180 ps.

Quantum yields
In the above kinetic analysis, we have assumed that all population initially imparted to the optically bright 1 MLCT state by the excitation pulse is transferred through ISC to the 3 MLCT surface. 73n assuming this initial 1 MLCT -3 MLCT v40 ISC quantum yield to be unity (f = 1), subsequently we can extract values for f for all of the remaining processes within our kinetic model in Fig. 6, which are also collated in Table 2.
Based on our model in Fig. 6a, the quantum yield for dissociation (f d ) can be determined by the relative magnitudes of the rate constants for VET (k VET ) and dissociation (k d ), where the rate constant k is related to the lifetime t through k = 1/t.Formally, f d = k d /(k d + k VET ) and returns a value of f d = 0.89 for the photodissociation of 1 in H 2 O.Meanwhile, the remaining fraction of population (f3 MLCT v=0 = 1 À f d ) will be funnelled into vibrationally cold 3 MLCT v=0 .In deriving these values for f d and f3 MLCT v=0 , an assumption is made that the timescale for VET remains constant for the duration of the vibrational cooling process in 3 MLCT.However, this is an over simplification and, in reality, the initial VET rate around the hot 3 MLCT v40 levels is likely to be much more rapid than the value for t VET reported in Table 2, due to the comparatively high density of vibrational levels and the smaller energy spacing (cf.closer to 3 MLCT v=0 ).As vibrational energy is quenched away from 3 MLCT v40 , relaxation will populate lower vibrational levels where the density of states is reduced and the inter-state spacing is increased, leading to a reduction in the VET rate as 3 MLCT v=0 is approached.VET out of these lower vibrational levels will therefore act as a rate-limiting step in vibrational cooling, [76][77][78] and our value for t VET in H 2 O most likely reflects these later stages of this process.Taking this into consideration, the value we extract for f d from this 'zero-order' model should be considered as an absolute upper limit for the fraction of population evolving towards dissociation.Even so, the value for f d derived here is still instructive and leads to a general conclusion that f d Z f3 MLCT v=0 for 1 -such a conclusion is broadly supported by previous work on other related Ru(II) systems, [31][32][33][34] which also suggest f d will be large.Given this, it is striking that the relative intensity of the 3 MLCT ESA features in the TAS are comparable to, or greater than, the spectral signature for the PCI complex formed through Ru-NA dissociation (cf.feature i vs. feature iii in Fig. 3a).This generally suggests that the ESA crosssection from the 3 MLCT state is significantly greater than the absorption cross-section of the PCI species.Similar findings can also be seen in the TAS for 1 in acetone (ESI, † Fig. S3), and are qualitatively supported by the predicted weak absorption intensity of the PCI complex at l 4 550 nm by the TD-B3LYP calculations in Fig. 2c.After dissociation, some fraction of the 'caged' PCIÁNA pair will geminately reform 1 (f GR ), the quantum yield for which is determined according to f GR = [k GR /(k GR + k s )] Â f d .This returns a value of f GR = 0.53.The remaining fraction of the PCIÁNA pair will also undergo a degree of diffusional rearrangement into 'separated' PCI and NA.In H 2 O, the portion of PCIÁ NA that diffusively separates will rapidly form 2, with a quantum yield of f 2 = 0.36 (defined as f 2 = f d À f GR ).Knowing this, and based on our above discussion that f d must be taken as an absolute upper limit, we can state more definitively that 0.36 r f d r 0.89 must be true.Consistent with f 2 , the quantum yield for the ground state recovery of 1 (f 1 ) will be defined as f 1 = 1 À f 2 = f GR + f3 MLCT v=0 = 0.64.These values for f 1 and f 2 are in excellent agreement with the percentage recovery of the GSB determined through our target analysis (see Fig. 5d), which recovers by B65%.We draw attention to the fact that had f 1 and f 2 simply been estimated through inspection of the percentage recovery of GSB feature ii in the TAS in Fig. 3b, the returned quantum yields would have been incorrect (cf.f 2 B 0.1), due to the strong spectral overlap of the signatures for 1 and 2 in the TAS.This emphasizes that, where possible, target analysis is a vital tool for extracting reliable kinetic information and quantum yields from UV/Vis TAS.

Formation and relaxation dynamics of the intermediate [Ru(bpy) 2 (NA)] 2+
The methodologies used here are not sensitive to tracing the precise pathways for the formation of the 'caged' PCIÁNA pair from the 3 MLCT state of 1, after flux leaves the vertical Franck-Condon region.After population leaves the 3 MLCT state, not only do we lose our optical handle on any population which may be transiently present in the 3 MC surface (recall the discussion in Section 3.2), but identifying precise geometries where strong electronic state couplings occur along the Ru-NA dissociation coordinate with theoretical calculations, such as CASSCF methods, becomes problematic due to active space impurities incurred by mixing with the large number of neardegenerate states.However, the schematic potential energy profiles in Fig. 8a depict a qualitative picture of the Ru-NA bond fission in 1, and the involvement of the MLCT and MC Fig. 8 (a) Schematic potential energy profiles of the 1 MLCT, 3 MLCT, 3 MC, 1 MC and singlet ground state (GS) along the Ru-NA bond dissociation coordinate in 1, including the major relaxation channels after excitation (hn) to the 'bright' 1 MLCT state.(b) Schematic potential energy surfaces for the significant electronic states in the PCI complex.From the top left, downwards: excited state population from the 3 MC of 1 (black) may land on the 3 MC PCI (purple) following NA ligand elimination; displacement through the 3 MC/ 1 MC ISC seam minimum (Fig. 9a) to the 1 MC PCI (red); and then through the 1 MC/GS CI minimum (Fig. 9b) leading to the GS of the PCI (green).On the GS, two possible isomers of PCI exist, labelled pro-cis and pro-trans (see ESI, † Table S6, for enlarged structures), into which a solvent molecule may be accommodated to complete the formation of 2.
This journal is © the Owner Societies 2014 states (based on earlier calculated potentials in related systems 64,73 ).Coupling between the 3 MLCT state and the 3 MC state most likely involves motion along anti-symmetric Ru-N stretch vibrations, which is then energetically stabilized by the distorted 3 MC state geometry, 27 as previously described for related complexes. 34oreover, our above kinetic analysis and quantum yields indicate very efficient mixing of the 3 MLCT v40 and 3 MC states (vide supra), and find that crossing between these states must take place on a timescale of o400 fs (cf.t d in Table 2).
As with simpler d 6 metal systems, the 3 MC state, which involves a singly occupied 4d x 2 Ày 2 orbital (see Fig. 7c), will become strongly stabilised upon Ru-NA dissociation and at infinite separation will transform into the lowest energy excited 3 MC state of the PCI complex -see ref. 79 and 80 for a more expanded discussion of this behaviour.0][81][82] Thus, in order to gain a greater understanding of how these electronically excited 3 MC PCI species relax to their electronic ground state (GS) before generating 2, we close by discussing the results of complementary CASSCF calculations aimed at identifying regions of strong electronic state coupling (either via internal conversion or ISC) within the excited state potential energy landscape of the PCI species.The findings from these calculations are summarized by the schematic potential energy surfaces for the PCI complex in Fig. 8b.
These calculations locate a region in nuclear co-ordinate space for strong ISC between the electronically excited 3 MC and 1 MC states (of the same character) in the PCI.The geometry presented in Fig. 9a is a minimum energy crossing point (MECP) on the 3N À 7 (where N is the number of atoms) 83 dimensional seam of electronic state degeneracy 84 between the 3 MC and 1 MC surfaces ( 3 MC/ 1 MC ISC seam), involving predominantly population of the 4d x 2 Ày 2 orbital.Due to the presence of a seam around the ISC minimum, 3 MC -1 MC population transfer via ISC is expected to be extremely efficient in the PCI.This MECP geometry has a combination of both trigonal bipyramid (TBP) and square pyramidal (SP) character around the Ru(II) metal centre, and also shows the vibrational motion responsible for lifting the degeneracy between the 3 MC and 1 MC surfaces (labelled g in Fig. 9a).This vibration corresponds to a contraction/expansion of all Ru-N bonds in the PCI, with the most pronounced motion lying along the axial Ru-N(bpy) coordinate.
Further to the 3 MC/ 1 MC ISC seam, a conical intersection (CI) 84 geometry, of quasi-TBP nature, subsequently connects the 1 MC surface with the electronic GS of the PCI, the calculated MECP of which is shown in Fig. 9b ( 1 MC/GS CI).Minimal distortion away from the (higher energy) MECP geometry of 3 MC/ 1 MC ISC seam is subsequently required to access the MECP geometry of the 1 MC/GS CI, which is pseudo-JT like due to the degeneracy imposed by partial occupancy of both the 4d xy and 4d x 2 Ày 2 orbitals. 85As can be seen in Fig. 9b, the MECP structure is 'skewed' from a pure TBP in that the axial-equatorial angles between bpy ligands differs from 901 with one angle at approximately 771, and the other at approximately 1001 -hence the structural designation of quasi-TBP.The equatorial-axial angle to the NA ligand is around 901, and the equatorialequatorial angles are all close to 1201.Further analysis of the CASSCF calculations indicate that the topography of this 1 MC/GS CI is of 'peaked' type in the terminology of Ruedenberg and coworkers 86 (as one would anticipate for a pseudo-JT type CI 84 ), and is calculated to lie around 0.5 eV below the MECP of the 3 MC/ 1 MC ISC seam discussed above (see Fig. 8b).The two nuclear motions that remove the degeneracy between the 1 MC state and GS, termed the gradient difference (g) and derivative coupling (h), are also presented in Fig. 9b, where the latter is typically responsible for enabling population transfer through the CI. 84Both of these motions lie in the equatorial plane, and evolution along a linear combination of these vibrational modes will lead to the creation of a PCI structure with an open coordination site on the electronic GS (i.e.quasi-SP).We note that the 1 MC/GS CI and associated connected minima identified in the PCI are very similar to those identified in related species (see, for example, ref. 55, 58 and 59), for both open-and closed-shell transition metal complexes, that undergo ligand photodissociation and subsequent non-adiabatic relaxation in their associated PCI complexes.
Depending on how population evolves through the 1 MC/GS CI, two different minima may be accessed on the electronic GS of the PCI, corresponding to two different structural isomers.The lowest energy of these GS isomers is calculated to be around 2.1 eV below the MECP of the1 MC/GS CI.In Fig. 8b this minimum is labelled pro-cis-PCI, together with its calculated structure, given that insertion of an H 2 O solvent molecule into its free coordinate site will give rise to the cis isomer of 2. Alternatively, rearrangement of the bpy ligands upon evolution through the CI could yield the pro-trans-PCI to enable subsequent formation of trans-2 mono-aquated photoproducts; the pro-trans-PCI is calculated to lie B0.4 eV above the procis-PCI species using B3LYP/cc-pVTZ-SDD.Note that, isomerisation as a result of excited state PCI relaxation has been reported previously for other metal complexes. 58,59Here though, we only observe experimental evidence for the formation of the cis isomer of 2, suggesting that the pro-trans-PCI structure is never accessed.Unlike the pro-cis-PCI, the pro-trans-PCI clearly exhibits much greater deviation from a true SP geometry due to steric distortion of the two bpy ligands induced by the adjacent NA (see structures in Fig. 8b), hence why the pro-cis-PCI lies lower in energy.This may be a contributing factor in the sole formation of pro-cis-PCI.Perhaps more importantly though, the electronic state relaxation of the PCI will occur within the 'caged' PCIÁNA pair, and the close proximity of the neighbouring NA partner will aid to prevent isomerisation to the pro-trans-PCI structure on steric grounds as flux relaxes through the 1 MC/GS CI.
The electronic state relaxation scheme presented for the PCI complex from these CASSCF calculations implies that during this process, a structural evolution from quasi-SPquasi-TBPquasi-SP must occur.Following dissociation in simpler systems, this rearrangement process was observed to be ultrafast (sub-200 fs). 80,82his is also likely to be the case here, meaning that any spectral signature for this process will be too fast to observe in our current TAS measurements.However, after relaxation to the GS, the nascent pro-cis-PCI will necessarily be formed vibrationally hot.Experimentally, we noted in Section 3.2 that the early time TAS in Fig. 3a (t o 25 ps) display some rapid depletion of the PCI feature iii, on the order of B10 ps.These earlier time dynamics may, in part, be due to relaxation of a vibrationally hot pro-cis-PCI in the GS, undergoing vibrational cooling through VET to the solvent. 76Therefore, our combined experimental and computational findings lead us to broadly conclude that after Ru-NA bond fission in 1, electronically excited PCI complexes in the nascent PCIÁNA pair can relax via an efficient sequential 3 MC -1 MC -GS coupling process, mediated by appropriate 3 MC/ 1 MC ISC and 1 MC/GS CI seams, on a sub-500 fs timescale, followed by vibrational relaxation of the pro-cis-PCI on the order of r10 ps, after which diffusional separation can occur and enable the formation of photoproduct 2.

Conclusions
One of the primary aims of the study presented here has been to increase our understanding of the key mechanistic pathways of photoactivation of a potential prodrug.Specifically, a synergy between experiment and theory has revealed that cis-[Ru(bpy) 2 (NA) 2 ] 2+ (1) has several efficient excited state pathways that are pivotal in the control of the branched kinetics that ultimately dictate the quantum yield for forming the reactive photoproduct, [Ru(bpy) 2 (NA)(H 2 O)] 2+ (2).The ability of target analysis to decipher accurate timescales and quantum yields from congested transient absorption spectra is critical to these investigations.In identifying these important 'crossroads', we can propose the following arguments, which afford the relatively high reported quantum yield of 2, f 2 = 0.36: (i) complete conversion of the initial excited 1 MLCT state population to the vibrationally excited 3 MLCT with a lifetime of o100 fs; (ii) efficient traversal of flux from the vibrationally excited 3 MLCT state to the dissociative 3 MC state, and subsequent dissociation with the quantum yield 0.36 r f d r 0.89, occurs on a o400 fs timescale; (iii) the nascent photodetached ligand may be viewed as part of a 'caged' product pair within the solvation shell, PCIÁNA; and (iv) the nascent PCIÁNA pair may either geminately reform 1 or undergo a degree of diffusional separation, both of which are kinetically competitive processes, where the later subsequently enables ultrafast (sub-1 ps) solvent substitution to finally form 2.
The results presented here provide an excellent basis for understanding the initial stages involved in photoactivation of potential metal complex prodrugs and how we can tailor molecular properties to improve their performance.One clear notion is that bulkier ligands will suffer severe caging in the immediate wake of dissociation, which will, in turn, significantly dictate the formation of the target species 2. The quantum yield for forming the target product will ultimately be governed by a combination of (i) the quantum yield for dissociation and (ii) the relative rates of diffusive separation versus geminate recombination for the nascent caged product pair.This work serves to highlight the significance of structure-dynamics-function relationships, with particular emphasis here, on the design of more efficacious photo chemotherapeutic agents.

Fig. 2
Fig. 2 Calculated static UV/Vis absorption spectra for (a) 1, (b) 2 and (c) PCI, both in the gas phase (dotted lines) and with an H 2 O PCM (solid lines).Calculated ground state structures shown inset.Spectra were calculated at the TD-B3LYP level of theory using a 28 electron SDD ECP for Ru and a cc-pVTZ basis set for C, O, N and H atoms. Calculated spectra were generated by convolution of the calculated squared transition dipole moments (TDM 2 ) for singlet excitations with Gaussian functions (50 nm FWHM).Experimentally measured spectra for 1 and 2 from Fig. 1 are shown in panels (a) and (b), respectively, for comparison (dashed lines).

Fig. 3
Fig. 3 UV/Vis TAS of 890 mM aqueous solution of 1 over pump-probe delay ranges (a) t = À1 to 25 ps and (b) t = À1 ps to 1.2 ns, following photoexcitation at 340 nm.Growth/decay features i-iv highlighted indicate the spectral location assigned to the 3 MLCT ESA, GSB of 1, absorption of the PCI and, finally, photoproduct 2, respectively.See main text for further details.

Fig. 4
Fig. 4 Representative fits (grey) to transient absorption spectra (black dashed) at (a) t = 1 ps and (b) t = 1 ns, showing the individual basis functions used in the 'target analysis' of the TAS in Fig. 2. Basis functions have been assigned to the 3 MLCT ESA (blue), PCI complex (green), GSB recovery of 1 (orange) and formation of the photoproduct 2 (cyan).Static UV/Vis absorption spectra in Fig. 1 are used as basis functions for 1 and 2 -see text for further details.

Fig. 6
Fig. 6 (a) Kinetic scheme illustrating the major pathways and associated time constants (t), used to model the photodissociation dynamics of 1 in H 2 O.The individual steps involved in this kinetic scheme are discussed further in the main text (Section 3.3).(b) Pictorial schematic of the post-dissociation processes involving the PCI and formation of the photoproduct 2.

Fig. 9
Fig.9(a) MECP of the CASSCF calculated3 MC/ 1 MC ISC seam.Equatorial and axial angles of this geometry as well as the gradient difference motion (g) are presented.(b) MECP of the CASSCF calculated 1 MC/GS CI, of quasi-JT type, and branching space.Equatorial and axial angles of this geometry as well as the gradient difference (g) and derivative coupling motions (h) are presented.

Table 2
Time constants (t) and quantum yields (f) extracted for the dissociation of 1 in both H 2 O. Values for t are extracted from kinetic fits to the traces in Fig. 5, in accordance with the kinetic model presented in Fig. 6.See main text for further details Open Access Article.Published on 17 July 2014.Downloaded on 31/03/2017 21:36:01.This article is licensed under a Creative Commons Attribution 3.0 Unported Licence.