
==== Front
Chem Sci
Chem Sci
SC
CSHCBM
Chemical Science
2041-6520
2041-6539
The Royal Society of Chemistry

39184296
d4sc01876a
10.1039/d4sc01876a
Chemistry
Real-time capture of nuclear motions influencing photoinduced electron transfer†
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4sc01876a

https://orcid.org/0000-0002-1212-8760
Kim Pyosang ab
Roy Subhangi c
Valentine Andrew J. S. d
https://orcid.org/0000-0003-2157-265X
Liu Xiaolin d
https://orcid.org/0000-0002-2430-221X
Kromer Sarah c
Kim Tae Wu a
https://orcid.org/0000-0001-7341-6240
Li Xiaosong d
https://orcid.org/0000-0001-7546-8618
Castellano Felix N. c
https://orcid.org/0000-0002-8450-6687
Chen Lin X. ab
a Chemical Sciences and Engineering Division, Argonne National Laboratory Lemont IL 60439 USA kimp@anl.gov
lchen@anl.gov

b Chemistry Department, Northwestern University Evanston IL 60208 USA
c Chemistry Department, North Carolina State University Raleigh NC 27695-8204 USA fncastel@ncsu.edu

d Chemistry Department, University of Washington Seattle WA 98195 USA xsli@uw.edu

22 8 2024
18 9 2024
22 8 2024
15 36 1476614777
20 3 2024
8 8 2024
This journal is © The Royal Society of Chemistry
2024
The Royal Society of Chemistry
https://creativecommons.org/licenses/by-nc/3.0/ This article is licensed under a Creative Commons Attribution-Non Commercial 3.0 Unported Licence. You can use material from this article in other publications without requesting further permissions from the RSC, provided that the correct acknowledgement is given and it is not used for commercial purposes.
Although vibronic coupling phenomena have been recognized in the excite state dynamics of transition metal complexes, its impact on photoinduced electron transfer (PET) remains largely unexplored. This study investigates coherent wavepacket (CWP) dynamics during PET processes in a covalently linked electron donor–acceptor complex featuring a cyclometalated Pt(ii) dimer as the donor and naphthalene diimide (NDI) as the acceptors. Upon photoexciting the Pt(ii) dimer electron donor, ultrafast broadband transient absorption spectroscopy revealed direct modulation of NDI radical anion formation through certain CWP motions and correlated temporal evolutions of the amplitudes for these CWPs with the NDI radical anion formation. These results provide clear evidence that the CWP motions are the vibronic coherences coupled to the PET reaction coordinates. Normal mode analysis identified that the CWP motions originate from vibrational modes associated with the dihedral angles and bond lengths between the planes of the cyclometalating ligand and the NDI, the key modes altering their π-interaction, consequently influencing PET dynamics. The findings highlight the pivotal role of vibrations in shaping the favorable trajectories for the efficient PET processes.

This study explores vibronic wavepacket dynamics in photoinduced electron transfer (PET). The dynamic interplay of PET and coherent wavepacket motions highlights the crucial role of vibrations in influencing efficient PET process.

National Science Foundation 10.13039/100000001 CHE-1856210 CHE-2247821 CHE-2247822 U.S. Department of Energy 10.13039/100000015 Unassigned Office of Science 10.13039/100006132 Unassigned Basic Energy Sciences 10.13039/100006151 Unassigned Argonne National Laboratory 10.13039/100006224 DE-AC02-06CH11357 pubstatusPaginated Article
==== Body
pmcIntroduction

Natural photosynthetic reaction centers set examples for nearly 100% quantum efficiency in photoinduced charge separation, driving biochemical reactions responsible for storing the energy necessary to sustain life on Earth.1 The detection of coherent wavepacket (CWP) motions within the reaction centers stimulated investigations on the role of nuclear motions in achieving highly efficient charge separation.2–4 Subsequent studies have identified the CWP motion as an electronic–vibrational (vibronic) coherence that persists over the timescale of the photoinduced charge separation.5–7 It was proposed that the vibronic resonance across the exciton and charge transfer state could harness electronic coherence as a robust function in the primary electron transfer processes even in the fluctuating environmental conditions within the protein matrix.8–10 Vibronic coherence has also garnered significant attention for its potential role in photoinduced electron transfer (PET), a fundamental process in artificial solar energy conversion systems.11–13 For instance, a 100 fs PET process in a blend of poly(3-hexylthiophene) (P3HT) and fullerene coincided with the emergence of CWP motions attributed to the vibrations of fullerene, highlighting the interplay of electron and vibrational motions during the primary charge separation.14 Furthermore, it was experimentally demonstrated that the vibronic coherence directly influences the photocurrent generation in other organic heterojunctions solar cells.15 These observations suggest that the vibronic coherence between the electron donor and acceptor could play a pivot role in facilitating efficient PET dynamics in artificial systems.

The aforementioned studies demonstrating the role of vibronic coherence in PET processes has inspired the present study on an electron donor–acceptor system featuring a transition metal complex (TMC) as the donor. The light absorption process directly gives rise to a charge transfer (CT) state in the TMCs, either from metals to ligands or vice versa. The CT transition functions as a driving force, initiating the primary ET reactions in many applications for solar energy conversion, photocatalysis, and optoelectronics.16–19 Recent advances in ultrafast optical and X-ray spectroscopies have revealed that the CT excitation results in significant structural rearrangements, accompanied by CWP motions associated with the structural changes in the vibronic manifold of the CT excited states.20–32 Furthermore, the CWP motions evolve along the excited state trajectories to the final states during intersystem crossing and spin-crossover dynamics, shedding light on the possibility of the vibronic coherence between different spin states.29,30,33–36 These examples imply that the CWP motions following the CT transition in the TMCs can exert a significant impact on subsequent PET toward the electron acceptor.

To investigate the effects of CWP motion in PET reactions, we designed an electron donor–acceptor system (2) with a cyclometalated Pt(ii) dimer (1) as the donor and naphthalene diimide (NDI) as the acceptor (Fig. 1a). As shown in previous studies, a metal–metal-to-ligand charge transfer (MMLCT, dσ* → π*) transition in Pt(ii) dimer directly induces ∼0.2 Å or larger Pt–Pt distance contraction as the electron density from an antibonding orbital is depleted.29,32,37 Meanwhile, a significant amplitude of CWP motions associated with Pt–Pt stretching vibration can be launched in the MMLCT excited state, as evidenced by ultrafast optical transient absorption and X-ray solution scattering experiments.29,32,36 In this study, using ultrafast broadband transient absorption (BBTA) spectroscopy with 20 fs (FWHM) pump pulses, we track the CWP motions during the PET from Pt(ii) dimer to NDI to identify the key reaction coordinates in 2. The results unveil the clear evidence for the certain CWP motions to directly modulate the NDI radical anion formation dynamics and to influence the PET trajectories across the MMLCT and charge-separated (CS) states.

Fig. 1 (a) Schematic molecular structures of Pt(ii) dimer (1) as the donor and Pt(ii) dimer with NDIs (2) as the donor–acceptor complex. (b) Steady-state absorption spectra for 1 and 2 dissolved in THF solvent. The inset shows a zoom of MMLCT absorption bands and pump spectrum used in the BBTA measurements.

Results and discussion

Excited state and PET kinetics

Pt(ii) dimer (1) exhibits the typical electronic transitions of a cyclometalated Pt(ii) dimer as reported in previous studies, featuring (1) ligand-centered (LC) transitions below 350 nm, (2) a mixture of LC and metal-to-ligand charge transfer (MLCT) transitions between 350–450 nm and (3) the MMLCT transition above 450 nm (Fig. 1b, S5–S7 and Table S1†).37,38 The electron donor–acceptor complex (2) also displays the MMLCT band above 450 nm, while a strong π–π* transition of the NDI electron acceptors appears with a clear vibronic feature below 400 nm in the steady-state absorption spectrum (Fig. 1b, S5–S7 and Table S1†).39 We investigated the PET reaction dynamics from Pt(II) dimer to NDI upon MMLCT excitation using BBTA spectroscopy with 20 fs pump pulses centered at 550 nm (see the inset in Fig. 1b). As a reference, the BBTA signals of 1 display rise dynamics across the entire probe spectral range within the first 2–3 picosecond (ps) time delay window (Fig. 2a). After the first few ps, the TA signal amplitudes persist without significant changes up to 150 ps (Fig. S8†). Global fitting, performed with the Glotaran program40 using a sequential kinetic model, reveals two time constants for the early rise kinetics: 364 fs and 1.76 ps. The analysis of evolution-associated difference spectra (EADS) for each time component demonstrates that these rise kinetics are attributed to the decay kinetics of the stimulated emission (SE) signal (Fig. 2c). The broad SE signal centered around 650 nm first emerges, which subsequently redshifts to 700 nm with a time constant of 364 fs and then is completely quenched with a time constant of 1.76 ps. Since the SE signal stems from the singlet MMLCT state, the time constant of 1.76 ps is attributed to intersystem crossing (ISC) from the singlet to the triplet MMLCT state. After undergoing vibrational cooling within 71 ps,38 the 3MMLCT state persists longer than 150 ps time window, which is consistent with the 3MMLCT lifetime of 341 ns obtained from the previous nanosecond TA measurement on the same Pt(ii) dimer.37

Fig. 2 (a and b) 2D contour plots of BBTA spectra and kinetic profiles at specific probe wavelengths for 1 (a) and 2 (b). All measurements were carried out in THF solvent. (c and d) Evolution associated difference spectra (EADS) for 1 (c) and 2 (d), obtained by global kinetic analysis on the BBTA data scanned up to 150 ps (Fig. S8†). Black arrows highlight the temporal evolution of stimulated emission (SE) signal. The BBTA data up to 70 fs were truncated for the global fitting in order to avoid the contribution of coherent spike and cross-phase modulation signal around time zero.

The BBTA signals of 2 exhibit distinctly different early TA kinetics (Fig. 2b) under the same experimental conditions compared as those of 1. The formation of NDI radical anion dominates the temporal and spectral behaviors of 2 due to the PET from Pt(ii) dimer to NDI. Using the global fitting of the sequential kinetic model, analysis of EADS demonstrates that the NDI anion spectrum rises with two time constants: 95 fs and 475 fs (Fig. 2d). Based on the computations of time-dependent density functional theory (TD-DFT) and natural transition orbital (NTO) (Fig. S5–7 and Table S1†), we assigned the initial excited state corresponding to the 1MMLCT state featured by the black EADS. Despite cross-phase modulation in the first ∼100 fs after the excitation, the black EADS indicates the significant quenching of the SE signal around 650 nm due to the PET on the similar time scale. The black EADS evolves to the red EADS with a 95 fs time constant, and then to the blue EADS with a 475 fs time constant. Notably, the red EADS differs significantly from the blue EADS with characteristic NDI radical anion bands at 608, 700, and 775 nm.39 Instead, the red EADS shows a broader absorption band at 608 nm and two blue-shifted bands at 667 and 750 nm. This distinction indicates that the charge separation (blue EADS) occurs through an intermediate CT (InCT) state formed between the Pt(ii) dimer and NDIs (shown by the red EADS) following the MMLCT excitation. After the completion of the charge separation, the spectrum of the NDI anion decays over time without significant changes in its spectral features (Fig. 2d and S8†). The multiple time constants found in the decay kinetics (1.83, 15.1, and 220 ps) suggests that different relaxation processes could be present, such as vibrational relaxation or ISC during a charge recombination (CR) to the ground-state. Although a detailed investigation of these dynamics could provide insight into the mechanism of CR, the current study focuses on the forward PET processes (Pt(ii) dimer → NDI) which can be tracked by the CWP behaviors to deduce their roles in the PET processes.

CWP motions

The CWP motions in both 1 and 2 launched by the 20 fs pump pulses were detected as oscillatory signals across the probe range overlaid with the population evolution signals (Fig. 3a and b). Fast Fourier transformation (FFT) of the beating signals after the removing the population evolution in the TA signal uncovers the frequencies of the CWP motions for 1 and 2 below 1000 cm−1, as shown in Fig. 3c and d. (Note: the THF solvent vibration at 912 cm−1 is induced by an impulsive stimulated Raman scattering (ISRS) process41). Previous ultrafast TA studies on cyclometalated Pt(ii) dimers, combined with normal mode analysis, showed a frequency upshift in the CWP motion associated with the Pt–Pt stretching vibration from ∼110 cm−1 for the ground-state to ∼150 cm−1 for the MMLCT excited-state as the result of the Pt–Pt distance shortening or effective bond order increase.29,30,36 Particularly, fs X-ray solution scattering study directly observed a contraction of Pt–Pt distance by ∼0.2 Å in the MMLCT excited-state, accompanied by the Pt–Pt stretching CWP motion with the 150 cm−1 frequency.32 In the current investigation, the CWPs of 1 and 2 also exhibit significant FFT amplitudes at ∼110 and ∼150 cm−1 (Fig. 3c and d), which are similar to the aforementioned frequencies of Pt–Pt stretching vibrations in the ground- and MMLCT excited-states, respectively.

Fig. 3 (a and b) Wavelength-resolved residual map (top) and beating signals (bottom) at specific wavelengths for 1 (a) and 2 (b). The oscillatory residuals were extracted by a global kinetic analysis on the BBTA kinetic traces. We performed inverse Fourier filtering on the CWP modes in 0–800 cm−1 in order to visualize the temporal dynamics of the CWP motions obscured by strong THF solvent vibration at 912 cm−1. The raw data are provided in Fig. S9† (c and d) wavelength-resolved FFT map (top) and spectra (bottom) at specific probe wavelengths for 1 (c) and 2 (d). The FFT was carried out for the residuals within the time range of 0.07–2.5 ps to avoid the contribution of cross-phase modulation around time zero.

Our fs-BBTA spectroscopy provides a broad probe spectral range of 550–850 nm, allowing us to precisely map out a phase-flip position associated with the Pt–Pt stretching CWP motion and its temporal evolution (black dotted arrows in Fig. 3a and b). Since the phase flip occurs around the energy minimum of the potential energy surface (PES), its spectral position provides the assignments of the electronic state where the CWP motions take place.30,42–46 Specifically, the phase flip in the ground-state CWP would occur around the ground-state bleaching (GSB) maximum, while that in the excited-state CWP arises near the SE or excited-state absorption (ESA) maximum.

The ∼110 cm−1 CWP motions in 1 and 2 commonly display the strong FFT amplitudes exclusively around the GSB region of 550–600 nm, supporting that the CWPs are prepared in the ground-state PES via the ISRS process.47 It is highly likely that the phase-flip for this ground-state CWP is located around the absorption maximum (520 nm), which falls beyond the probe range of BBTA spectroscopy. An interference between the ∼110 and ∼150 cm−1 CWP motions is observed around the GSB region (∼600 nm) due to the overlap of these oscillations. This interference is distinct from the phase flip that arises from the ∼150 cm−1 excited state CWP observed beyond the GSB region.

In Fig. 3a, the wavelength-resolved residual map for 1 reveal that the phase flip emerges first around 650 nm and then redshifts to ∼700 nm, which is consistent with the redshift of the SE maximum from 650 to 700 nm with the time constant of 364 fs (Fig. 2c). The temporal evolution of the phase flip demonstrates that the 150 cm−1 CWP motion directly modulates the SE spectral shape and amplitude along the probe wavelength axis (also see Fig. S10†). In addition, the CWP undergoes a damping with a time constant of ∼610 fs (Fig. S11 and Table S2†). These results confirm that the 150 cm−1 Pt–Pt stretching CWP launched by the 1MMLCT excitation decays in the PES of the 1MMLCT excited state during the ISC to the 3MMLCT state. In contrast, the phase flip for the 150 cm−1 CWP motions in 2 disappears at the SE peak of 650 nm in <300 fs, while a new phase flip emerges at 800 nm in early time and blueshifts to 750 nm. Moreover, the 150 cm−1 CWP motions of 2 display an initial amplitude growth within ∼100 fs, similar to the TA rise time of ∼100 fs associated with forming the InCT state in the PET process (Fig. S11 and Table S2†). These features differentiate the origin of the 150 cm−1 CWPs for 2 from that for 1. The decay and rise of the phase flip at ∼650 and ∼800 nm demonstrate that the PET reaction quenches the SE signal of the 1MMLCT excited-state and concomitantly enhances the ESA signals of the InCT and CS states. In addition, the damping time of ∼1 ps indicates that the Pt–Pt stretching CWP sustains its phase during the PET reaction (Fig. S11 and Table S2†). Therefore, the 150 cm−1 CWP motion observed in 2 is attributed to the nuclear motion associated with the Pt–Pt stretching vibration in the InCT and CS states.

In Fig. 3c and d, FFT analysis also reveals other CWP motions in 200–700 cm−1 frequency range. The weak FFT amplitudes were observed at 230 and 293 cm−1 in 1, which were not clearly detected in previous work.29 The 230 and 293 cm−1 modes largely appear beyond the onset of the absorption spectra at ∼600 nm and thus can be reasonably assigned to the CWPs in the 1MMLCT excited-state. The 230 cm−1 CWP motion is also observed in 2, while its spectral amplitude distribution is distinct. The mode is predominantly detected around 608 and 667 nm, which are near the peaks of InCT state represented by red EADS (Fig. 2d). 2 also shows new CWP motions with the frequencies of 270, 308, 345, 540 and 678 cm−1 around 608 nm peak of NDI anion. Interestingly, the wavelength-resolved residual map displays the beating patterns of these modes (dominantly 678 cm−1 frequency) with the phase flip around 608 nm (black dotted arrow in Fig. 3b). Additionally, the amplitude of the 678 cm−1 beating signal gradually increases over a period of ∼0.8 ps, as shown in the oscillation at 615 nm (Fig. 3b).

Franck–Condon vs. Herzberg–Teller type CWP oscillations

The FFT amplitude spectra plotted as a function of probe wavelength provide further characterization of the CWP motions in 2 (Fig. 4a and S12†). The FFT amplitude of the 678 cm−1 mode is predominantly present around the NDI radical anion band at 608 nm. Similar behavior is observed in the FFT amplitudes of 270, 308, 345, and 540 cm−1 CWPs (Fig. S12†). The spectra show a distinct valley at 608 nm due to the phase flip of the CWP motion, which is consistent with the wavelength resolved residual map (black dotted arrow at 608 nm in Fig. 3b). These CWP motions are attributed to Franck–Condon type oscillations. As illustrated in Fig. 4b, the periodic propagation of the CWP motions along its vibrational coordinates leads to spectral modulations of the NDI radical anion due to the displacement between the PESs of the lowest (S1) and higher excited states (Sn) in the ESA region. Specifically, the 678 cm−1 CWP motion causes the peak of NDI radical anion band to oscillate between wavelengths larger (∼615 nm) and smaller (∼595 nm) than 608 nm with its vibrational frequency. Although we do not directly excite NDI, the phase flip at the peak of the NDI radical anion absorption strongly indicates that the CWP in 2 originates from the nuclear motions associated with the CS active vibration pertaining to the NDI core.48 However, the PET time scales of ∼100 and ∼500 fs are not sufficiently short to impulsively generate 50 fs period of 678 cm−1 CWP mode in the CS state. This result suggests a different mechanism for the formation of 678 cm−1 CWP motion in the CS state of 2.

Fig. 4 (a) FFT amplitude distribution plots for the 678, 230 and 150 cm−1 modes as a function of probe wavelength. EADSs with the time constants of 475 fs (red) and 1.83 ps (blue) are presented to compare the spectra of intermediate CT (InCT) and CS states with FFT amplitude distribution. Red arrows indicate the spectral position of phase-flip detected in the wavelength-resolved residual map (b) schematic illustration for the two types of CWP motions: Franck–Condon and Herzberg–Teller type oscillations. The thicker and thinner arrows in Herzberg–Teller type oscillation highlights the difference in oscillator strength between the left and right positions of PES.

The FFT amplitude for the 230 cm−1 CWP displays two bands with peaks around 608 and 667 nm. Importantly, these peaks match well with the characteristic bands at 608 and 667 nm of the red EADS corresponding to the InCT state of the PET reaction (Fig. 4a). This feature indicates that the 230 cm−1 CWP motion only modulates the amplitudes of both ∼608 and 667 nm bands without the spectral shifts, in contrast to the 678 cm−1 motion. Such CWP motions can be interpreted as Herzberg–Teller type oscillation, as illustrated in Fig. 4b.49–52 In this case, the CWP motion cannot modulate the TA spectral shape and amplitude because the displacement between the two PESs is negligible or zero; therefore, the transition energies on both sides of the PES are identical. When the nuclear motions, however, can differentiate the oscillator strength on the right and left sides of the PES, it can be possible to probe the intensity modulation while maintaining the spectral shape. This phenomenon occurs through a vibronic coupling that allows the intensity borrowing from one electronic state to another at a specific molecular structure compared to the opposite structure.49–53 Thus, the Herzberg–Teller type beating pattern provides the strong evidence for non-Condon effect of 230 cm−1 CWP motion. Given that the transition strengths at ∼608 and 667 nm are proportional to the population of the InCT state in the PET processes, it can be suggested that there is the vibronic coupling between the 1MMLCT and InCT states along the 230 cm−1 vibrational coordinates.

The Pt–Pt stretching vibration (150 cm−1 mode) exhibits FFT amplitude spectra characterized by a sharp peak at ∼608 nm and a broad peak at ∼700 nm (Fig. 4a). These features match the NDI anion bands in the CS state (blue EADS in Fig. 2d). Notably, the oscillation residual map reveals the absence of phase flips associated with the 150 cm−1 oscillation at ∼608 nm and ∼700 nm (Fig. 3b). The presence of a phase flip near 750 nm (black dotted arrow in Fig. 3b and red arrow in Fig. 4a) complicates the assignment of the origin of the 150 cm−1 CWP. Nevertheless, the FFT spectral features at ∼608 and 700 nm strongly suggests the non-Condon effect in the 150 cm−1 vibration, leading to alternation in the intensity of the NDI radical anion spectrum. In contrast to the 230 cm−1 CWP motion, the 150 cm−1 CWP only affects the TA spectra associated with the CS state (∼608 and 700 nm at the blue EADS). Thus, the non-Condon effect for the 150 cm−1 Pt–Pt stretching CWP motion implies that the vibronic coupling between the InCT and CS states occurs along the 150 cm−1 vibrational coordinate.

Solvent-polarity dependence of PET and CWP dynamics

To verify the correlations of the PET kinetics and the CWP evolution, we also investigated the solvent dependence of the PET and CWP dynamics (Fig. S13 and 14†). It is well-known from Marcus theory that the solvent reorganization energy in the ET processes changes with its dielectric constant, and thus the PET rate.54 Generally, the polar solvent stabilizes the charged species (radical cations and anions, i.e., CS state) and thus increases the driving force of the PET reaction. We anticipated that a polar solvent would enhance the PET rate and affect the damping or rise time of the CWP motions. This approach allows us to identify the CWP motions associated with the PET reaction coordinates.55–57 As shown in Fig. 5c, the rise kinetics of NDI anion peak at ∼608 nm become faster as the polarity increases; toluene (dielectric constant ε ≈ 2.38) < THF (ε ≈ 7.58) < acetonitrile (MeCN, ε ≈ 37.5). The global kinetic analysis further reveals the solvent polarity dependency of the PET processes (Fig. 2d, 5a and b): (1) the time scale of ∼100 fs corresponding to the formation of the InCT state (the red EADS) remain almost unchanged across different solvent polarities. (2) the transition from the InCT to the CS state (the blue EADS) becomes faster with increasing polarity, toluene (952 fs) > THF (475 fs) > MeCN (224 fs). (3) the CR processes also occur faster in higher polarity solvents.

Fig. 5 (a and b) EADS for 2 in toluene and acetonitrile (MeCN), obtained by global kinetic analysis on the BBTA data scanned up to 150 ps (Fig. S13†). (c) Early TA kinetic traces observed at the NDI anion peaks at ∼608 nm in toluene, THF and MeCN. Inset shows the kinetic traces up to 20 ps.

Using the short-time Fourier transformation (STFT) method, we analyze the temporal evolution of the CWP motions and their correlation with the PET kinetics observed in different solvents (Fig. 6, S16 and S17†). The STFT analysis captures the time evolution of the amplitude of the frequency content extracted from the oscillatory signals, and provides information about the dephasing dynamics of CWP motions (see the details in Fig. S15†). As shown in the wavelength-resolved STFT maps (Fig. 6a and b), the ∼100 cm−1 ground state wavepacket dephases only around the GSB region (550–600 nm). It is because this CWP doesn't undergo any electronic transition and vibrational relaxation in the ground state PES. Therefore, we exclude the ground-state wavepacket contribution on the excited state CWP dynamics beyond the ground state absorption.

Fig. 6 (a and b) Wavelength-resolved STFT maps for 2 in THF. (c) Time traces for STFT amplitudes for 230 cm−1 mode observed at ∼667 nm (the STFT amplitudes at each time point were integrated in the frequency range of 200–300 cm−1), 150 cm−1 mode observed at 820 nm, and 678 cm−1 mode observed at ∼615 nm. The time traces are also detected in toluene (top), THF (middle) and MeCN (bottom).

The amplitude of the 230 cm−1 mode significantly decreases during the PET from the 1MMLCT to InCT and almost completely vanishes within 0.6 ps. Additionally, its dephasing dynamics appear to be nearly independent of solvent polarity. The 150 cm−1 mode exhibits a noticeable dependence on solvent polarity. The time point at which the amplitude reaches half of its maximum value becomes smaller with increasing solvent polarity, ranging from ∼0.9 ps in toluene, ∼0.6 ps in THF to ∼0.4 ps in MeCN. This trend illustrates that the 150 cm−1 CWP undergoes faster dephasing dynamics in higher polar solvents. Furthermore, the STFT maps show that the amplitude of the ∼150 cm−1 CWP around 608 nm (blue dotted circle in Fig. 6a and b) increases during the PET process from 1MMLCT to InCT and becomes evident during the PET process from InCT to CS. This feature indicates that the non-Condon effect for the 150 cm−1 mode is delayed in time and affects the CS spectra later. This emergence suggests that the 150 cm−1 CWP induces the vibronic coupling between the InCT and CS states. Most interestingly, the 678 cm−1 CWP clearly show a distinct trend from the others, with its amplitude rise during the PET processes. As the solvent polarity increases, the STFT amplitude reaches its maximum faster at earlier delay times, ∼0.8 ps in THF to ∼0.6 ps in MeCN. The onset of the amplitude decay for the 678 cm−1 mode also varies with the solvent polarity, which can be interpreted as the consequence of the polar solvents accelerating both CS and CR processes. This finding demonstrates that the rise dynamics of this mode become faster as the polarity increases, which is consistent with the faster rise of NDI anion in MeCN solvent (Fig. 5c). Although the STFT amplitude at early time is the lowest in toluene, the time point at the maximum value remains similar to THF, followed by slower dephasing. This result indicates that the 678 cm−1 CWP in toluene begins the damping concurrently with the rise dynamics during the PET processes, distinct from those observed in other solvents. The 270, 308, 345 and 540 cm−1 CWPs do not show the polarity dependence (Fig. S17†).

To support the STFT analysis, we performed a Fourier filtering analysis using a super-Gaussian window.45 We selectively filtered out the 230, 150 and 678 cm−1 FFT bands and then conducted inverse Fourier transformation to retrieve the filtered frequency domain data back to the time domain one. As shown in Fig. 7, the 230 cm−1 CWP motion doesn't exhibit the solvent-polarity dependence. However, the damping of the 150 cm−1 CWP motion significantly accelerates in polar solvent, and the rise dynamics of the 678 cm−1 CWP also become faster in polar solvent. Although narrow filters require careful use due to potential artefacts, these results are consistent with STFT analysis. Notably, the beating maps for 230 and 150 cm−1 CWPs reveal the changes in the beating patterns around 667 and 608 nm during the PET processes (blue dotted boxes in Fig. 7a and b). These features indicate that the 230 and 150 cm−1 CWP motions propagate along the PET reaction coordinates, modulating the transition strengths of TA spectra associated with InCT and CS states, respectively. These results support the non-Condon effects of these CWP motions. Specifically, the non-Condon effect of the 230 cm−1 CWP indicates the vibronic coupling between the 1MMLCT and InCT along the 230 cm−1 vibrational coordinate. This coupling might explain why the dephasing dynamics of the 230 cm−1 CWP are independent of the solvent polarity, which is similar to the trend in the PET rate change from the 1MMLCT to InCT across different solvent polarities. For the 150 cm−1 CWP motion, the non-Condon effect suggests the vibronic coupling between the InCT and CS along this vibrational coordinate, strongly supported by the solvent-polarity dependence of its dephasing dynamics.

Fig. 7 (a–c) Fourier filtering results for the 230, 150 and 678 cm−1 CWP motions. Blue dotted boxes highlight the temporal evolution of beating pattern around 667 nm for 230 cm−1 CWP and 608 nm for 150 cm−1 CWP. Blue dotted arrow indicates the phase flip for the 678 cm−1 CWP around 608 nm. Black dotted arrows highlight the dephasing and rise dynamics of 150 and 678 cm−1 CWP motions, which are dependent on the solvent polarity.

Compared to the damping dynamics of 230 and 150 cm−1 CWPs, the 678 cm−1 CWP dynamics clearly exhibit an intensity growth during the PET processes in both STFT and Fourier filtering analyses. Since the amplitude rise dynamics depend on the PET rate varied by solvent polarity, we rule out the possibility that this CWP is a spectator mode orthogonal to the PET reaction coordinate and, thus, not active in the reaction. As mentioned in the previous section, the impulsive generation via the PET processes is also excluded because the PET time scales of ∼100 and 224–950 fs are much longer than the vibrational period of 678 cm−1 CWP (∼50 fs). Hence, the 678 cm−1 CWP transfer strongly suggests the vibronic mechanism for the entire PET processes, 1MMLCT → InCT → CS, which will be discussed below.

Key PET reaction coordinates

To gain insight into how the molecular vibrations of 150, 230, and 678 cm−1 modes contribute to the PET reactions, we conducted calculations involving optimized structures and normal modes in both ground- and CS states (Fig. 8a and S18–23†). Based on the knowledge of which the 150 cm−1 CWP motion originates from the Pt–Pt stretching vibration, our initial focus is on examining the relevant normal modes in the ground- and CS states. Similar to 1, the Pt–Pt stretching vibration in 2 exhibits an upshift in frequency from 113 cm−1 in the ground- to 148 cm−1 in CS state. It is because the one electron transferred to the NDI comes from a dσ* antibonding orbital, leading to Pt–Pt bond contraction. Notably, the nuclear motion triggers an out-of-plane (OOP) motion of 2-phenylpyridine (ppy) ligands (left in Fig. 8a and S19†), which arises due to the short distance (∼3.6 Å) between these ligands. More interestingly, the induced OOP motion concurrently leads to a dihedral angle (DA) change between the ppy ligand and the NDI moiety, resulting in ∼5° alteration (Fig. S21†). The DA alteration can potentially increase the overlap of π orbitals between the ppy and NDI, thereby accelerating the transfer of electrons from ppy to NDI. Such a torsional motion has been noted to enhance π-interactions, leading to rapid intramolecular CT in copolymers and covalent organic frameworks built from donor–bridge–acceptor configurations.58–60

Fig. 8 (a) Calculated normal modes in the ground-state, assigned to 150, 230 and 678 cm−1 CWP motions. (b) Illustration of PET reaction trajectories.

Based on the Pt–Pt stretching motion, we searched for an OOP vibration within the 200–300 cm−1 frequency range, potentially corresponding to the 230 cm−1 mode. The calculated OOP vibrational mode has a frequency around 282 cm−1 in both the ground- and CS states, displaying an alteration of the DA by ∼5° (middle in Fig. 8a, S19 and S22†). The OOP motion is also observable around a similar frequency in the normal mode analysis for 1 (Fig. S20†), consistent with the CWP motion observation at ∼230 cm−1 in 1 and 2 (Fig. 3c and d). An intriguing collective vibration involving both the ppy ligand and NDI core was calculated around 687 cm−1 in both ground and CS states (right in Fig. 8a). This nuclear motion alters ∼0.06 Å in the bond distance between the ppy and NDI constituents (Fig. S23†). Therefore, the molecular vibrations of 150, 230 and 678 cm−1 frequencies are expected to be PET active due to the alternation of DA and bond length between the ppy and NDI moieties.

Based on the analyses of CWP motions and their corresponding normal modes, the PET reaction mechanism for 2 can be proposed, as illustrated in Fig. 8b. In the initial PET process (∼100 fs), the 1MMLCT and InCT states are vibronically coupled via the 230 cm−1 OOP vibration initiated by the MMLCT excitation. This vibration leads to an enhancement in π-interaction between the ppy and NDI moieties through the modulation in the DA configuration. Even though the 230 cm−1 mode completely decoheres in the InCT state, the π-interaction keeps being boosted through the vibration of the Pt–Pt stretching mode (150 cm−1), which also alters in the DA configuration. Consequently, the vibronic couplings driven by the OOP and Pt–Pt stretching vibrations render the efficient PET pathway through the multiple PESs. More importantly, the PET pathway is accompanied by coherence transfer of the 678 cm−1 CWP to the CS state. The DA modulations enable vibronic wavefunction delocalization across the multiple PESs, leading to the transfer of the 678 cm−1 CWP along the entire PET processes. Along this line, the vibronic couplings driven by the OOP and Pt–Pt stretching vibrations may engender multiple conical intersections (CIs) through the PET trajectories (Fig. 8b). The CI is a crossing region where two adiabatic PESs are degenerate at a particular geometry. In this vicinity, strong vibronic mixing of the electronic and nuclear degrees of freedom dominates, leading to the breakdown of the Born–Oppenheimer approximation and enabling nonadiabatic transitions from the higher to lower surfaces.61–63 Such vibronic coupling significantly influences the CWP dynamics at the CI, facilitating its passage through the CI and subsequent population in the final state.64,65 In 2, the DA modulations enhance the nonadiabatic couplings at the CIs, conserving the 678 cm−1 wavepacket phase and leading to efficient coherence transfer to the CS state. In other words, the 678 cm−1 CWP in 2 passes through the two CIs in the crossing regions between the PESs of the MMLCT, InCT and CS states, followed by dephasing processes during the CR processes. In this context, the observation of non-Condon effects can imply that the 230 and 150 cm−1 vibrations are the coupling modes to form the CIs. In addition, it is conceivable that the InCT state represents a vibronic CT state delocalized between the ppy and NDI, vibronically bridging the 1MMLCT and CS states during the PET processes. The strong dependence of the TA kinetics for the InCT state on the solvent polarity also supports this interpretation. Furthermore, the calculated collective motion involving the ppy and NDI entities may signify a distinct intramolecular vibration associated with the vibronic CT.

Conclusions

In summary, we experimentally detected that MMLCT excitation of Pt(ii) dimer launches certain CWP motions (150, 230 and 678 cm−1) which directly alter the formation and evolution of the NDI radical anion during PET. The temporal correlation between the CWP and the PET dynamics revealed these modes as key reaction coordinates for the PET processes. Based on normal mode analysis, these findings strongly suggest that the CWP motions are coupled to the PET reaction coordinates by modulating the energies associated with vibronic states to be delocalized between the ligand and NDI. As a result, we also discovered the dynamic interplay of these vibronic motions enables the 678 cm−1 vibronic coherence to survive the entire PET processes, even in ∼1 ps CS in toluene. Our findings demonstrate new possibility to leverage the vibronic effect to determine the most desirable trajectory for an efficient PET reaction. The vibrational motions coupled to the CT transition of TMCs will offer an ideal avenue for harnessing the vibration-assisted coherent ET in artificial solar energy conversion systems.

Data availability

The data supporting this article have been included as part of the ESI.†

Author contributions

PK: performed the ultrafast BBTA experiments and analysis, as well as wrote the manuscript. PK, XL, FNC and LXC: conceptualized and supervised the project and reviewed and edited the manuscript. SR and SK: performed the synthesis and characterization, wrote the synthetic methodology. AJSV and XL: performed the computations, wrote the computational methodology. PK and TWK: built up the ultrafast BBTA spectroscopic technique.

Conflicts of interest

There are no conflicts to declare.

Supplementary Material

SC-015-D4SC01876A-s001

P. K., S. R., A. J. S. V, X. L., S. K., T. W. K., X. L., F. N. C and L. X. C. as well as computational work are supported in part by the Ultrafast Initiative of the U. S. Department of Energy, Office of Science, Office of Basic Energy Sciences, through Argonne National Laboratory under Contract No. DE-AC02-06CH11357. The development of computational method for simulating excited state dynamics is supported by the National Science Foundation (CHE-1856210 to X. L). The synthesis and initial spectroscopic characterization were supported by the National Science Foundation (CHE-2247821 to L. X. C and CHE-2247822 to F. N. C).
==== Refs
References

Fleming G. R. Grondelle R. v. Phys. Today 1994 47 48 55 10.1063/1.881413
Vos M. H. Rappaport F. Lambry J.-C. Breton J. Martin J.-L. Nature 1993 363 320 325 10.1038/363320a0
Ando K. Sumi H. J. Phys. Chem. B 1998 102 10991 11000 10.1021/jp982659l
Novoderezhkin V. I. Yakovlev A. G. van Grondelle R. Shuvalov V. A. J. Phys. Chem. B 2004 108 7445 7457 10.1021/jp0373346
Romero E. Augulis R. Novoderezhkin V. I. Ferretti M. Thieme J. Zigmantas D. van Grondelle R. Nat. Phys. 2014 10 676 682 26870153
Fuller F. D. Pan J. Gelzinis A. Butkus V. Senlik S. S. Wilcox D. E. Yocum C. F. Valkunas L. Abramavicius D. Ogilvie J. P. Nat. Chem. 2014 6 706 711 10.1038/nchem.2005 25054941
Ma F. Romero E. Jones M. R. Novoderezhkin V. I. van Grondelle R. J. Phys. Chem. Lett. 2018 9 1827 1832 10.1021/acs.jpclett.8b00108 29584941
Novoderezhkin V. I. Romero E. van Grondelle R. Phys. Chem. Chem. Phys. 2015 17 30828 30841 10.1039/C5CP00582E 25854607
Romero E. Novoderezhkin V. I. van Grondelle R. Nature 2017 543 355 365 10.1038/nature22012 28300093
Sahu A. Kurian J. S. Tiwari V. J. Chem. Phys. 2020 153 224114 10.1063/5.0029027 33317308
Brédas J.-L. Sargent E. H. Scholes G. D. Nat. Mater. 2017 16 35 44 10.1038/nmat4767
Scholes G. D. Fleming G. R. Chen L. X. Aspuru-Guzik A. Buchleitner A. Coker D. F. Engel G. S. van Grondelle R. Ishizaki A. Jonas D. M. Lundeen J. S. McCusker J. K. Mukamel S. Ogilvie J. P. Olaya-Castro A. Ratner M. A. Spano F. C. Whaley K. B. Zhu X. Nature 2017 543 647 656 10.1038/nature21425 28358065
Rafiq S. Scholes G. D. J. Am. Chem. Soc. 2019 141 708 722 10.1021/jacs.8b09059 30412671
Falke S. M. Rozzi C. A. Brida D. Maiuri M. Amato M. Sommer E. De Sio A. Rubio A. Cerullo G. Molinari E. Lienau C. Science 2014 344 1001 1005 10.1126/science.1249771 24876491
Bian Q. Ma F. Chen S. Wei Q. Su X. Buyanova I. A. Chen W. M. Ponseca C. S. Linares M. Karki K. J. Yartsev A. Inganäs O. Nat. Commun. 2020 11 617 10.1038/s41467-020-14476-w 32001688
Schultz D. M. Yoon T. P. Science 2014 343 1239176 10.1126/science.1239176 24578578
Wenger O. S. J. Am. Chem. Soc. 2018 140 13522 13533 10.1021/jacs.8b08822 30351136
Twilton J. Le C. Zhang P. Shaw M. H. Evans R. W. MacMillan D. W. C. Nat. Rev. Chem 2017 1 0052 10.1038/s41570-017-0052
Yam V. W.-W. Au V. K.-M. Leung S. Y.-L. Chem. Rev. 2015 115 7589 7728 10.1021/acs.chemrev.5b00074 26158432
Iwamura M. Takeuchi S. Tahara T. Acc. Chem. Res. 2015 48 782 791 10.1021/ar500353h 25646861
Katayama T. Northey T. Gawelda W. Milne C. J. Vankó G. Lima F. A. Bohinc R. Németh Z. Nozawa S. Sato T. Khakhulin D. Szlachetko J. Togashi T. Owada S. Adachi S.-i. Bressler C. Yabashi M. Penfold T. J. Nat. Commun. 2019 10 3606 10.1038/s41467-019-11499-w 31399565
Auböck G. Chergui M. Nat. Chem. 2015 7 629 633 10.1038/nchem.2305 26201738
Hainer F. Alagna N. Reddy Marri A. Penfold T. J. Gros P. C. Haacke S. Buckup T. J. Phys. Chem. Lett. 2021 12 8560 8565 10.1021/acs.jpclett.1c01580 34468159
Lemke H. T. Kjær K. S. Hartsock R. van Driel T. B. Chollet M. Glownia J. M. Song S. Zhu D. Pace E. Matar S. F. Nielsen M. M. Benfatto M. Gaffney K. J. Collet E. Cammarata M. Nat. Commun. 2017 8 15342 10.1038/ncomms15342 28537270
Kunnus K. Vacher M. Harlang T. C. B. Kjær K. S. Haldrup K. Biasin E. van Driel T. B. Pápai M. Chabera P. Liu Y. Tatsuno H. Timm C. Källman E. Delcey M. Hartsock R. W. Reinhard M. E. Koroidov S. Laursen M. G. Hansen F. B. Vester P. Christensen M. Sandberg L. Németh Z. Szemes D. S. Bajnóczi É. Alonso-Mori R. Glownia J. M. Nelson S. Sikorski M. Sokaras D. Lemke H. T. Canton S. E. Møller K. B. Nielsen M. M. Vankó G. Wärnmark K. Sundström V. Persson P. Lundberg M. Uhlig J. Gaffney K. J. Nat. Commun. 2020 11 634 10.1038/s41467-020-14468-w 32005815
Cho S. Mara M. W. Wang X. Lockard J. V. Rachford A. A. Castellano F. N. Chen L. X. J. Phys. Chem. A 2011 115 3990 3996 10.1021/jp109174f 21361265
Kim P. Kelley M. S. Chakraborty A. Wong N. L. Van Duyne R. P. Schatz G. C. Castellano F. N. Chen L. X. J. Phys. Chem. C 2018 122 14195 14204 10.1021/acs.jpcc.8b01636
Mewes L. Ingle R. A. Megow S. Böhnke H. Baranoff E. Temps F. Chergui M. Inorg. Chem. 2020 59 14643 14653 10.1021/acs.inorgchem.0c00902 32677823
Kim P. Valentine A. J. S. Roy S. Mills A. W. Chakraborty A. Castellano F. N. Li X. Chen L. X. J. Phys. Chem. Lett. 2021 12 6794 6803 10.1021/acs.jpclett.1c01289 34270259
Kim P. Valentine A. J. S. Roy S. Mills A. W. Castellano F. N. Li X. Chen L. X. Faraday Discuss. 2022 237 259 273 10.1039/D2FD00009A 35642929
Kim T. W. Kim P. Mills A. W. Chakraborty A. Kromer S. Valentine A. J. S. Castellano F. N. Li X. Chen L. X. J. Phys. Chem. C 2022 126 11487 11497 10.1021/acs.jpcc.2c02256
Leshchev D. Valentine A. J. S. Kim P. Mills A. W. Roy S. Chakraborty A. Biasin E. Haldrup K. Hsu D. J. Kirschner M. S. Rimmerman D. Chollet M. Glownia J. M. van Driel T. B. Castellano F. N. Li X. Chen L. X. Angew. Chem., Int. Ed. 2023 62 e202304615 10.1002/anie.202304615
Paulus B. C. Adelman S. L. Jamula L. L. McCusker J. K. Nature 2020 582 214 218 10.1038/s41586-020-2353-2 32528090
Gaynor J. D. Sandwisch J. Khalil M. Nat. Commun. 2019 10 5621 10.1038/s41467-019-13503-9 31819052
Monni R. Capano G. Auböck G. Gray H. B. Vlček A. Tavernelli I. Chergui M. Proc. Natl. Acad. Sci. U. S. A. 2018 115 E6396 10.1073/pnas.1719899115 29941568
Rafiq S. Weingartz N. P. Kromer S. Castellano F. N. Chen L. X. Nature 2023 620 776 781 10.1038/s41586-023-06233-y 37468632
Chakraborty A. Yarnell J. E. Sommer R. D. Roy S. Castellano F. N. Inorg. Chem. 2018 57 1298 1310 10.1021/acs.inorgchem.7b02736 29336558
Brown-Xu S. E. Kelley M. S. J. Fransted K. A. Chakraborty A. Schatz G. C. Castellano F. N. Chen L. X. J. Phys. Chem. A 2016 120 543 550 10.1021/acs.jpca.5b11233 26759897
Bhosale S. V. Jani C. H. Langford S. J. Chem. Soc. Rev. 2008 37 331 342 10.1039/B615857A 18197349
Snellenburg J. J. Laptenok S. P. Seger R. Mullen K. M. van Stokkum I. H. M. J. Stat. Software 2012 49 1 22
Dhar L. Rogers J. A. Nelson K. A. Chem. Rev. 1994 94 157 193 10.1021/cr00025a006
McClure S. D. Turner D. B. Arpin P. C. Mirkovic T. Scholes G. D. J. Phys. Chem. B 2014 118 1296 1308 10.1021/jp411924c 24428706
Rafiq S. Scholes G. D. J. Phys. Chem. A 2016 120 6792 6799 10.1021/acs.jpca.6b07796 27510098
Xian R. Corthey G. Rogers D. M. Morrison C. A. Prokhorenko V. I. Hayes S. A. Miller R. J. D. Nat. Chem. 2017 9 516 522 10.1038/nchem.2751 28537597
Fitzpatrick C. Odhner J. H. Levis R. J. J. Phys. Chem. A 2020 124 6856 6866 10.1021/acs.jpca.0c03912 32786657
Jaiswal V. K. Kabaciński P. Nogueira de Faria B. E. Gentile M. de Paula A. M. Borrego-Varillas R. Nenov A. Conti I. Cerullo G. Garavelli M. J. Am. Chem. Soc. 2022 144 12884 12892 10.1021/jacs.2c04565 35796759
Pollard W. T. Lee S. Y. Mathies R. A. J. Chem. Phys. 1990 92 4012 4029 10.1063/1.457815
Aster A. Bornhof A.-B. Sakai N. Matile S. Vauthey E. J. Phys. Chem. Lett. 2021 12 1052 1057 10.1021/acs.jpclett.0c03641 33470820
Kano H. Saito T. Kobayashi T. J. Phys. Chem. A 2002 106 3445 3453 10.1021/jp012493f
Yoneda Y. Sotome H. Mathew R. Lakshmanna Y. A. Miyasaka H. J. Phys. Chem. A 2020 124 265 271 10.1021/acs.jpca.9b09085 31867968
Kuramochi H. Aoyama G. Okajima H. Sakamoto A. Kanegawa S. Sato O. Takeuchi S. Tahara T. Angew. Chem., Int. Ed. 2020 59 15865 15869 10.1002/anie.202004583 32432809
Sohn S. H. Heo W. Lee C. Kim J. Joo T. J. Phys. Chem. A 2019 123 6904 6910 10.1021/acs.jpca.9b05613 31306027
Bizimana L. A. Carbery W. P. Gellen T. A. Turner D. B. J. Chem. Phys. 2017 146 084311 10.1063/1.4976995 28249416
Marcus R. A. Sutin N. Biochim. Biophys. Acta, Rev. Bioenerg. 1985 811 265 322 10.1016/0304-4173(85)90014-X
Rafiq S. Fu B. Kudisch B. Scholes G. D. Nat. Chem. 2021 13 70 76 10.1038/s41557-020-00607-9 33288893
Yoneda Y. Kudisch B. Rafiq S. Maiuri M. Nagasawa Y. Scholes G. D. Miyasaka H. J. Am. Chem. Soc. 2021 143 14511 14522 10.1021/jacs.1c01863 34474559
Rafiq S. Scholes G. D. Chem. Phys. Lett. 2017 683 500 506 10.1016/j.cplett.2017.03.020
Hwang I. Beaupré S. Leclerc M. Scholes G. D. Chem. Sci. 2012 3 2270 2277 10.1039/C2SC20078C
Roy P. Jha A. Yasarapudi V. B. Ram T. Puttaraju B. Patil S. Dasgupta J. Nat. Commun. 2017 8 1716 10.1038/s41467-017-01928-z 29170455
Kim T. W. Jun S. Ha Y. Yadav R. K. Kumar A. Yoo C.-Y. Oh I. Lim H.-K. Shin J. W. Ryoo R. Kim H. Kim J. Baeg J.-O. Ihee H. Nat. Commun. 2019 10 1873 10.1038/s41467-019-09872-w 31015440
Worth G. A. Cederbaum L. S. Annu. Rev. Phys. Chem. 2004 55 127 158 10.1146/annurev.physchem.55.091602.094335 15117250
Domcke W. Yarkony D. R. Annu. Rev. Phys. Chem. 2012 63 325 352 10.1146/annurev-physchem-032210-103522 22475338
Boeije Y. Olivucci M. Chem. Soc. Rev. 2023 52 2643 2687 10.1039/D2CS00719C 36970950
Kühl A. Domcke W. J. Chem. Phys. 2002 116 263 274 10.1063/1.1423326
Egorova D. Domcke W. J. Photochem. Photobiol., A 2004 166 19 31 10.1016/j.jphotochem.2004.04.034
