1
Ultrafast temporary charge transfer in pyrrolidinyl–benzonitrile and pyrrolyl–benzonitrile in the gas phase

2
4-(N-pyrrolidinyl)benzonitrile and 4-(N-pyrrolyl)benzonitrile were excited by an ultrashort pulse at 270 nm to their La (S2) states and then probed by ionization at long wavelengths.

3
Parent and fragment ion signals show components with time constants <100 fs which we attribute to ultrafast relaxation to the Lb (S1) state.

4
From this short time we infer a conical intersection between the Lb and La surfaces.

5
The wave packet can branch there, one part going temporarily to a strongly displaced state.

6
Its shift was concluded from an anisotropy observed only there.

7
The only excited state known to have a large displacement is the charge-transfer (CT) state.

8
The positive anisotropy indicates that the CT state belongs to the 2A species for both molecules.

9
For pyrrolylbenzonitrile, this is in contrast to previous assignments.

10
The anisotropy, and a coherent oscillation observed in pyrrolidinylbenzonitrile, support the idea that the amino-group twist is an important component of the CT reaction coordinate.

Introduction

11
On UV excitation in solvents of suitable polarity, many donor–acceptor substituted aromatics such as 4-dimethylaminobenzonitrile (DMABN) show dual fluorescence, originating from an Lb-like and a charge-transfer (CT) state (reviews.1–4)

12
The pathway for charge transfer is still controversial (see e.grefs. 5–7.)., although the amino-group twist is most popular as the reaction coordinate.

13
Gas-phase investigations with femtosecond time resolution have resulted in new insights, revealing the key role of a conical intersection (CI),8–10 a feature only beginning to be considered in CT reactions.

14
Furthermore, coherent oscillations observed with this method gave evidence of the relaxation direction and thus also provides information about components of the reaction coordinate.9,10

15
The method involves exciting such a molecule by a weak femtosecond UV pump pulse (in our case at 270 nm) and then probing it by ionization by an intense IR pulse, measuring the ion yields mass-selectively as a function of pump–probe delay.

16
Using various probe wavelengths (810, 1300 and 2100 nm), we investigated in this work 4-(N-pyrrolidinyl)–benzonitrile (PYRBN) and 4-(N-pyrrolyl)–benzonitrile (PBN).

17
The donor group in the former is a saturated amino group, whereas in the latter it is the aromatic pyrrole.

18
According to spectroscopy in cold molecular beams (PYRBN,11–13 PBN;14 see also the work on pyrrolylbenzene15), the longest-wavelength absorption leads to the Lb state (1B2 in the idealized symmetry C2v).

19
The transition is weak and polarized perpendicularly to the long molecular axis.

20
This band is insensitive to para-substituents (since the latter are in a node of the Lb wavefunction) and the medium (due to the small transition moment and little change of dipole moment).

21
A much stronger absorption, parallel-polarized and partially covering the weaker one in solution (solution spectra of PYRBN:,16 of PBN:17–19) leads to the La state (symmetry 2A1).

22
The excitation wavelength is near the absorption maximum of this band.

23
As will be further discussed in sec.

24
4, in the compounds with a saturated amino group the La state involves promotion of an electron from the HOMO to the LUMO, the former being π-antibonding and the latter weakly π-bonding between the amino group and the benzene ring (see e.g.refs. 20 and 21).

25
As to conclude from the fluorescence polarization, the CT state also belongs to the A species in these cases;22 the La and CT states are thus just two locations on the same 2A surface.1,4,8,21,23,24

26
This surface is lowered from La to CT by interaction with a higher π-antibonding nA state on excursion along certain coordinates, a candidate for which is the amino-group twist.8,21

27
Along the same coordinates, the lower Lb state rises and crosses the 2A surface in a symmetry-allowed conical intersection (CI).

28
After initial excitation to La, a wave packet can thus rapidly reach the CI and branch there to the CT and Lb states.

29
The temporary population in CT is then depleted towards Lbvia a path leading around the lower cone of the CI (see sec.

30
4).

31
We found this scheme to apply for DMABN and PIPBN8–10 and present evidence that it is also suitable for PYRBN.

32
Surprisingly, it seems also to apply for PBN.

33
This was not expected, because quantum chemical calculations predicted a 2B symmetry for the CT state, involving excitation from another pyrrole orbital to the LUMO.20,25

34
On the other hand, another calculation26 only reported on a CT state of 2A type.

35
In pyrrolylbenzene two CT states, one with B and the other with A symmetry, were found very close-by.27

36
The pyrrole ring is slightly twisted versus the benzene ring in S0 and Lb of pyrrolylbenzene15 and probably similarly in PBN (symmetry group C2).

37
Both rings are expected to be planar.

38
In PYRBN, however, the amino N is slightly pyramidal, and the two rings are not twisted against each other (symmetry group Cs) according to a DFT calculation.28

39
This was also inferred from photoelectron spectra.16

40
(Such spectra were reported for PBN, too.17,29)

41
A previous prediction of a slight twist and puckering of the pyrrolidino ring by molecular mechanics12 is probably less reliable than the DFT calculation.28

42
The kinetics and fluorescence of PYRBN30–36 and PBN17,18,29,37 were previously investigated in solution.

43
It was found that the Lb → CT charge-transfer reaction required a solvent of at least medium polarity (to lower the CT state) for PYRBN, whereas PBN shows this reaction already in alkanes.

Experimental

44
The preparation of PYRBN and PBN is described in refs. 17, 29 and 33.

45
The two species were investigated in the gas phase in the set-up (details described in ref. 38) heated to 130 °C to avoid adsorption; the substances themselves were at 20 °C.

46
The molecules were pumped by a frequency-tripled pulse of a Ti:sapphire laser system (270 nm, 75 fs, 109 W cm–2) and then probed by nonresonant dissociative ionization at the fundamental wavelength (810 nm, 50 fs, 1013 W cm–2); in further experiments, probing was done by the output of a parametric generator (Topas) at 1300 and 2100 nm (40 fs, 1013 W cm–2).

47
The two beams had linear polarizations with an angle of 55° (“magic angle”) between each other in order to avoid time dependences induced by molecular rotation.

48
In some experiments at 810 nm we also compared the signals at parallel and perpendicular polarizations.

49
More experimental details are described in refs. 8 and 10.

50
The shape of the signals did not depend on the pump intensity (varied by a factor of about three).

51
This proves that they are produced by a single pump photon.

52
Varying the probe intensity by a similar factor did not change the time constants, although the later parts of the signals (coming from lower molecular states) were raised at higher probe intensities.

Results

53
The UV spectra of PYRBN and PBN (not shown), measured in the gas phase at 80 °C, look very similar to that of DMABN.39

54
The strong structureless La band has a maximum at 273 (PYRBN) and 270 nm (PBN) (4.54 and 4.59 eV).

55
A weak precursor at longer wavelengths is certainly the Lb band for both molecules; according to the supersonic-beam spectra, their 00 bands are at 32 01811–13 and 34 424 cm–114 (3.97 and 4.26 eV), respectively.

56
Hence the initial excitation at 270 nm places the molecules at an energy of 0.57 and 0.33 eV above the Lb minimum.

57
The ion signals investigated with time resolution were the parent ion (M+) and the most intense fragment (about 10% of M+), which was M+ – 1 for PYRBN and M+ – 28 for PBN.

58
Time dependent signals at different probe wavelengths are shown for the first few picoseconds in Figs. 1 and 2 for PYRBN and 3 for PBN.

59
Similarly as found in PIPBN10,40 and DMABN,9 PYRBN exhibits a periodic modulation of the signals between about 0.5 and 3 ps which is extracted (see below) and Fourier-transformed in Fig. 2 from the data for the parent ion with 1.3 µm probing.

60
By contrast, although the PBN signals show some irregularities, an oscillation could not unambiguously be extracted from them.

61
The time evolution was measured up to 600 ps (not shown); in this part, the signals simply decayed doubly exponentially.

62
In a first step, we simulated the nonperiodic part of the signals by a sum of exponentials with time constants τi (Table 1), convoluted with pump and probe pulses.

63
Subtracting this part from the PYRBN data left the oscillating part (Fig. 2), whose Fourier transformation (Fig. 2 bottom) gave the wavenumber and period τosc also reported in Table 1.

64
This part could also be simulated by a damped oscillation of the type exp(–t/τdeph)cos(2πt/τosc + ϕ)Using τosc from the Fourier transform, we deduced the dephasing time and phase from this simulation.

65
The simulation procedure is described in detail in refs. 10 and 38.

66
For a given molecule, the set of time constants is the same, independent of probe conditions.

67
Only the coefficients of the exponentials depend on the probe wavelength and the ion mass observed.

68
This is interpreted by population Li flowing consecutively from one observation window Li to the next.

69
This flow is described by rate equations, whose solution is a sum of exponentials for each population (the same, independent of probe conditions).

70
Each Li has a mass spectrum of its own, and their dependence on the probing conditions help to distinguish the signals from the different Li; the spectra are represented by the relative ionization cross sections mσi, i.e. the probabilities to produce an ion of mass m from window Li.

71
They are simulation parameters and are also given in the table.

72
The data are redundant for deducing the τi and the populations Li.

73
Fig. 2 shows by the thin lines the simulation components (mσi × Li), whose sums for all i add up to the signals for mass m.

74
The very short τ1 is needed to reproduce the steep rise of the parent signals.

75
τ2 causes, for instance, the delayed rise of the fragment as compared to the parent, clearly visible at 1300 nm.

76
τ3 represents the decay of the initial spike, which is evident in particular at the longer probe wavelengths.

77
The flat part of the signals extending from <1 ps to beyond the range shown is due to contributions from L3 and L4, the former decaying in the same time τ3 as the latter rises.

78
The decays in the range up to 600 ps (τ4, τ5) are not shown.

79
Fig. 4 shows the anisotropies r calculated for the two molecules from the signals at 810 nm (other wavelengths not investigated) with the probe polarization parallel and perpendicular to that of the pump: r = (S – S)/(S + 2S).Obviously, for both molecules r ≈ 0 for the parent ion, whereas for the fragment r begins with positive values and decays in the picosecond range.

80
The asymptotic value, which is difficult to determine, is set to 0 in the figure.

Discussion

Relaxation path for PYRBN

81
The PYRBN signals are very similar to those for DMABN and PIPBN.8,10

82
Because also the molecules are very similar, we can directly take over the assignments.

83
Consequently we also use the same scheme of potential surfaces and relaxation paths (Fig. 5).

84
Similarly as in DMABN, an evidence for transient population of the CT state is based on the anisotropy.

85
This can be seen as follows: The ionization itself is isotropic, which is confirmed by the zero-anisotropy observed for the parent ion.

86
Photodissociation of the parent ion by absorption of an additional probe photon can, however, give rise to an anisotropy of the fragment ions, which will be positive if there is a resonant ionic transition at 810 nm (1.53 eV) polarized parallel to the pump laser, that is parallel to the long axis of the molecule.

87
After vertical ionization from a state with near-equilibrium geometry, there is no such resonance, as to judge from the energy difference of photoelectron bands.16

88
Therefore a substantial distortion is necessary in the neutral molecule to reach a geometry with such a resonance in the ion; the only low-lying excited state with strong displacement is the CT state.

89
This assignment implies that the anisotropy decay of Fig. 4 contains a contribution from this internal relaxation (time constant τ3) and is not purely due to rotation; evidence for such a double origin was presented for DMABN in .ref. 9

90
As already discussed previously,8 the temporary population of the CT state implies that low-quantum yield CT emission should be detectable from the gas phase and in nonpolar solvents; the latter has indeed been observed in DMABN.41,42

91
It is interesting in this context that a recent work on spectroscopy of similar molecules in cold supersonic beams also suggested population of the CT state, if the excess energy is sufficient.43

92
It is instructive to have a closer look to the ionic resonance.

93
Energy differences of ion states are just equal to orbital energy differences in the neutral molecules.

94
Orbitals and their energies are easy to understand by the help of their nodal pattern.

95
The HOMO of PYRBN has the symmetry species b1 and is π-antibonding between donor and acceptor moieties (Fig. 6).

96
The ground state of the parent ion (2B1) has a hole in this orbital (single occupancy).

97
The HOMO – 1 has a node containing the amino group; it has the species a2.

98
An ionic 2A2 → 2B1 transition would be perpendicularly polarized, which would disagree with our parallel anisotropy in observation window L3 (or the strongly displaced state claimed to correspond to CT).

99
HOMO – 2 provides the right symmetry (b1), since an ionic 2B1 → 2B1 transition would be parallel; it is π-bonding between amino donor and the benzonitrile parts.

100
A displacement coordinate drastically weakening this π interaction may approach the two b1 orbitals sufficiently, so that a resonant transition in the ion will be possible, as indicated in Fig. 6.

101
There are not many coordinates with this property.

102
A twist of the amino group was predicted to do it,16,20 and Figs. 5 and 6 were therefore drawn assuming that the twist is an important (but not the only, see below) component of the reaction coordinate.

103
A strong pyramidalization of the amino group would probably be less effective, and planarization would do the opposite.

104
So, as in PIPBN10 and DMABN,8,9 we assign for PYRBN the times below 100 fs (τ1, τ2) to departure from the Franck–Condon region and further relaxation through a conical intersection (CI) to both the Lb and CT minima (Fig. 5).

105
The conclusion that the path leads through a CI can again be obtained from the shortness of τ2.

106
The next time constant (τ3) corresponds to Lb ↔ CT equilibration around the lower cone of the CI.

107
As before, the much longer τ4 is ascribed to intersystem crossing, probably to Tn, which then can rapidly equilibrate with T1.

108
From there, the molecule decays within τ5 to the ground state, from which ionization was too weak to be detected.

109
(Any excited state would be detectable).

110
This process may be accompanied by decomposition such as H elimination, as suggested in PIPBN9.

The oscillation in PYRBN

111
The observed oscillation period is 530 fs.

112
Its damping (dephasing) time is 1.4 ps.

113
The only excited state with sufficiently long lifetime is Lb (or the observation window L4, representing the Lb-CT equilibrium.

114
This equilibrium is expected to lie on the Lb side, because the CT state is higher in energy).

115
So the oscillation must be a vibration in the Lb well.

116
In fact, converting the period τosc to wavenumbers ν by ν = 1/osc, one receives 64 cm–1, very close to a value (69 cm–1) observed in supersonic-beam spectra of the Lb state of PYRBN.11–13

117
There are not many vibrations with such a low wavenumber in this molecule.

118
According to semiempirical PM3 calculations for S0, there are only three below 100 cm–1: (1) arching (42 cm–1), where the CN and pyrrolidine groups bend in phase out of plane versus the benzene ring, (2) amino group twist (54 cm–1) and (3) another arching (85 cm–1), where the CN and pyrrolidino-N bend in phase out of plane, whereas the C2H4 end group moves out of phase.

119
Because the electronic interaction of the benzene with its substituents is only minor in the Lb state, these wavenumbers are not expected to be very different in this excited state.

120
Hence, the observed oscillation should be assigned to the lower-frequency arching or to the twist vibrations.

121
Arching involves the inversion motion of the amino nitrogen.

122
However, both vibrations deserve a little more scrutiny, taking the symmetry of PYRBN into account: According to DFT calculations,28 in the ground state, the twist angle θ = 0 and the wagging angle (angle between the benzene and pyrrolidine planes) is only ω = 8°.

123
Any barrier at ω = 0 will be smaller than in DMABN (where ω = 15° in S044), where it has recently been reevaluated and found to be <100 cm–1 in S0 and <50 cm–1 in Lb.45

124
(As mentioned in the Introduction, a molecular-mechanics calculation predicted θ = 15° in S0.12

125
This would imply also a minor barrier at θ = 0°.

126
It would not much change the following considerations).

127
The Lb potential must be very similar, since the 00 transition is by far stronger than the adjacent bands in the vibrational structure.11–13

128
Therefore both types of vibration will occur between two turning points that are symmetry-equivalent.

129
Also the ion is expected, due to the partial double bond between the amino group and the benzene ring, to have its potential minimum at θ = 0 and ω ≈ 0, not displaced versus that of Lb.

130
A vibration (belonging to an antisymmetric species) in such a symmetric potential will therefore give rise to signal maxima twice in each period.

131
(This is the equivalent to the spectroscopic selection rule Δv = 2 for antisymmetric vibrations in electronic transitions.

132
It will also apply for the 69-cm–1 structure in the PYRBN supersonic beam spectrum,11–13 if our assignment is correct).

133
Therefore we should divide by two to convert the observed period to wavenumbers: ν = (2osc)–1We thus calculate ν = 32 cm–1 for the vibration.

134
At present we cannot decide between the two possible assignments to twist or arching.

135
In DMABN we observed both types of oscillations,9 whereas in PIPBN we attributed a strong modulation to the arching vibration.9

136
Its assignment was based on a substituent mass effect observed for it.

137
In PIPBN, the neutral and ion potentials are displaced against each other (due to the nonplanar piperidine ring), so that the factor of two in (3) does not apply, although it does apply in PYRBN for symmetry reasons: The potential is symmetric to ω = 0 in this molecule.

138
(We assume that the vibration takes place above any barrier)..

139
The wavenumber thus calculated for PIPBN was 31.4 cm–1, very close to that calculated by (3) for PYRBN.

140
Therefore we tend to the assignment as arching (involving amino group inversion).

141
The other vibration might be excited, too, even if it is not observed because possibly the Lb and ion potentials are parallel in this direction, so that this vibration would not be Franck–Condon active and would not give rise to a signal modulation.

142
The wavenumber observed by us (0.5 × 64 cm–1) deviates from that in the cold supersonic jet (0.5 × 69 cm–1) by much more than the experimental error (0.5 cm–1, see Table 1).

143
A plausible interpretation assumes that levels with v ≫ 1 are excited in our experiment, and the wavenumber difference is due to an anharmonic shift.

144
Also the La potential will be symmetric to θ = 0 and ω = 0, again due to the partial double bond.

145
Therefore these vibrations will not (or nearly not) be Franck–Condon active in the pump transition.

146
They are thus not expected to be excited by the pump laser.

147
Instead, they will be stimulated by acceleration on the last slope on entering the Lb well from the CI.

148
The direction of this slope therefore contains either θ or ω as an important component.

149
Both can be justified (and would be supported by an assignment): It was argued above (see also Fig. 5) and previously8–10 that a CI is probably displaced from the Lb minimum towards a direction with θ as a component.

150
On the other hand, a b-type vibration such as the inversion (ω) could be stimulated by the 2A → 1B relaxation in the CI, which needs a symmetry-breaking distortion.9

151
It is remarkable that such oscillations can be understood as described, by considering a pathway with intermediate excursion of the wave packet via a conical intersection (as also recommended in .ref. 46)

152
A previous model47 would only consider electronic matrix elements and Franck–Condon factors between La and Lb; the latter factors are ≈ 0, because the potentials are not shifted against each other in the θ and ω directions.

Relaxation path for pyrrolylbenzonitrile

153
The kinetic behavior of PBN is very similar to that of PYRBN: There seems to be the same number of steps or observation windows, and the time constants differ by not more than a factor of 2.

154
This similarity is not expected, because PBN has 4 more π electrons than PYRBN, so that the level scheme could be different.

155
In fact, the lowest-energy CT state was predicted to belong to the 1B species20,25 (in contrast to the 2A CT state of DMABN or PYRBN; see also the reviews2,4Ref. 17. found evidence that two different CT states play a role in PBN, depending on the solvent).

156
This state involves a transition from the HOMO of a2 symmetry (localized on the carbons of pyrrole, node containing the nitrogen) to the b1 LUMO of the benzonitrile moiety.

157
The orbital energies thus calculated20 are shown as functions of the twist angle θ in Fig. 6.

158
In this scheme, ionization from the CT state would leave a hole in the a2 orbital, giving an ion ground state 2A2.

159
The positive anisotropy found (Fig. 4) indicates a resonance in the ion parallel to the molecular axis.

160
We thus would need a 2A2 → 2A2 transition in the ion or the corresponding transition between a pair of a2 orbitals.

161
It is obvious from Fig. 6 that the available orbital pair does not approach each other along the twist, and it seems hard to imagine any other coordinate, relevant to a charge transfer reaction, that could do so: Both orbitals have a node containing the substituents.

162
On the other hand, the pair of b1 orbitals just does the required energy approach along the twist coordinate.

163
If we postulate that the CT state observed by us involves excitation of a b1 electron (HOMO – 1 in Fig. 6b), certainly to the LUMO (b1) so that a 2A state would result, then the positive anisotropy would be understood by an ionic 2B1 → 2B1 transition in much the same way as for PYRBN (Fig. 6a).

164
This postulate is very much encouraged by a CASPT2 calculation which predicts practically degenerate CT states of 2B and 2A symmetry in pyrrolylbenzene at perpendicular geometry.27

165
We suppose that attaching a para-cyano substituent (resulting in PBN) will lower the 2A state more than the 2B state.

166
In fact, another CASPT2 calculation reports for PBN at low energy only a 2A state with CT character.26

167
We therefore interpret our anisotropy data as a support that the lowest CT state of PBN belongs to the 2A species and is thus on the same potential energy surface as the initially excited La state.

168
This implies that the same scheme for the surfaces and relaxation paths can be used as for PYRBN (Fig. 5).

169
In contrast to PYRBN, no oscillations were observed in PBN (although some irregularities in the signals might indicate that there are some).

170
It may be that the Lb potential and that of the ion are better parallel in PBN than in PYRBN, so an existing vibration would cause less signal modulation in PBN, because it would not be Franck–Condon active in ionization.

Concluding remarks

171
It was surprising that the same potential energy scheme and relaxation paths (Fig. 5) could be used for PYRBN and PBN.

172
Apparently the states arising from the additional π electrons in PBN are higher in energy, including also a CT state of B symmetry.

173
It is also remarkable how little the early times τ1τ3 differ between the molecules DMABN,8 PIPBN,10 PYRBN and PBN investigated by our method.

174
The largest variation was found in τ3, which was 1.0, 1.9, 0.53 and 0.91 ps for the four compounds, respectively.

175
τ3–1 is the sum of the forward (Lb → CT) and backward (CT → Lb) rates in the Lb/CT equilibration, with the rate for the exothermic backward reaction dominating.

176
The CT → Lb reaction is less exothermic in the gas phase for PBN than for PYRBN.

177
This can be concluded from the fact that CT fluorescence is already observed in nonpolar environments for PBN, whereas at least medium polar solvents are needed for PYRBN.

178
Hence, the longer τ3 for the former is not surprising.

179
The same can be said about the pair PIPBN and DMABN.

180
It seems also natural that reverse differences of rates were found in solution for the forward (Lb → CT) reactions, whose reaction enthalpy is the negative of that of the CT → Lb reaction.

181
It may be worth repeating8–10 a comment on the amino-group twist as a component of the reaction coordinate.

182
Whereas this twist gives rise to observable consequences such as the anisotropy (in this work) or a coherent oscillation (this work and refs. 9 and 10) and many others,1–4 it will not be the only component: It seems natural that the electron transferred in the CT state to the benzonitrile moiety causes CC bond lengthening there.

183
A CASPT2 calculation for PBN in fact found such effects.26

184
Because bond stretching is energetically more effective than twisting, such components can give rise to steep slopes, which may explain our very fast initial relaxations (τ1,τ2).