Werner
Fuß
*a,
Wolfgang
Rettig
b,
Wolfram E.
Schmid
a,
Sergei A.
Trushin
ac and
Tomoyuki
Yatsuhashi
ad
aMax-Planck-Institut für Quantenoptik, D-85741 Garching, Germany
bHumboldt-Universität zu Berlin, Institut für Chemie, Brook-Taylor-Str. 2, D-12489 Berlin, Germany
cB.I. Stepanov Institute of Physics, Belarus Academy of Sciences, 220602, Minsk, Belarus
dDepartment of Chemistry, Graduate School of Science, Osaka City University, 3-3-138 Sugimoto, Sumiyoshi, 558-8585, Osaka, Japan
First published on 13th April 2004
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. Parent and fragment ion signals show components with time constants <100 fs which we attribute to ultrafast relaxation to the Lb (S1) state. From this short time we infer a conical intersection between the Lb and La surfaces. The wave packet can branch there, one part going temporarily to a strongly displaced state. Its shift was concluded from an anisotropy observed only there. The only excited state known to have a large displacement is the charge-transfer (CT) state. The positive anisotropy indicates that the CT state belongs to the 2A species for both molecules. For pyrrolylbenzonitrile, this is in contrast to previous assignments. 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.
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). The transition is weak and polarized perpendicularly to the long molecular axis. 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). 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). The excitation wavelength is near the absorption maximum of this band.
As will be further discussed in sec. 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). 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 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 Along the same coordinates, the lower Lb state rises and crosses the 2A surface in a symmetry-allowed conical intersection (CI). After initial excitation to La, a wave packet can thus rapidly reach the CI and branch there to the CT and Lb states. The temporary population in CT is then depleted towards Lbvia a path leading around the lower cone of the CI (see sec. 4). We found this scheme to apply for DMABN and PIPBN8–10 and present evidence that it is also suitable for PYRBN. Surprisingly, it seems also to apply for PBN. 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 On the other hand, another calculation26 only reported on a CT state of 2A type. In pyrrolylbenzene two CT states, one with B and the other with A symmetry, were found very close-by.27
The pyrrole ring is slightly twisted versus the benzene ring in S0 and Lb of pyrrolylbenzene15 and probably similarly in PBN (symmetry group C2). Both rings are expected to be planar. 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 This was also inferred from photoelectron spectra.16
(Such spectra were reported for PBN, too.17,29) 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 The kinetics and fluorescence of PYRBN30–36 and PBN17,18,29,37 were previously investigated in solution. 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.
The shape of the signals did not depend on the pump intensity (varied by a factor of about three). This proves that they are produced by a single pump photon. 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.
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. 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. 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. By contrast, although the PBN signals show some irregularities, an oscillation could not unambiguously be extracted from them. The time evolution was measured up to 600 ps (not shown); in this part, the signals simply decayed doubly exponentially.
![]() | ||
Fig. 1 The time-dependent PYRBN parent (mass 172 u) and fragment (171 u) ion signals (normalized to the maxima), measured at the indicated probe wavelengths at the magic angle between pump and probe polarizations. The parent is actually about 10 times stronger than the fragment. The solid lines represent simulations of the nonperiodic parts (see text). |
![]() | ||
Fig. 2 Upper panel: The PYRBN parent ion signal with probing at 1310 nm. The thin solid lines show the simulated contributions 172σ1L1![]() ![]() |
![]() | ||
Fig. 3 The time-dependent PBN parent (mass 168 u) and fragment (140 u) ion signals (normalized to the maxima), measured at two probe wavelengths at the magic angle between pump and probe polarizations. The parent is actually about 10 times stronger than the fragment. The solid lines represent simulations of the nonperiodic parts (see text). |
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. 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. This part could also be simulated by a damped oscillation of the type
exp(−t/τdeph)cos(2πt/τosc![]() ![]() | (1) |
Exponential parts of PYRBN | ||||||
---|---|---|---|---|---|---|
i | 1 | 2 | 3 | 4 | 5 | Anisotropy r |
τ i | 5![]() ![]() |
40![]() ![]() |
540![]() ![]() |
160![]() ![]() |
1100![]() ![]() |
|
172 σ i (810 nm) | 1 | 0.243 | 0.235 | 0.262 | 0.157 | 0 |
171 σ i (810 nm) | 1 | 0.188 | 0.175 | 0.225 | 0.135 | 0.2 |
172 σ i (1300 nm) | 1 | 0.40 | 0.144 | 0.128 | ||
171 σ i (1300 nm) | 0 | 1 | 0.674 | 0.550 | ||
172 σ i (2100 nm) | 1 | 0.277 | 0.025 | 0.018 |
Exponential parts of PBN | ||||||
---|---|---|---|---|---|---|
i | 1 | 2 | 3 | 4 | 5 | Anisotropy r |
τ i | 5![]() ![]() |
48![]() ![]() |
910![]() ![]() |
104![]() ![]() |
352![]() ![]() |
|
168 σ i (810 nm) | 1 | 0.680 | 0.373 | 0.389 | 0.123 | 0 |
140 σ i (810 nm) | 0 | 0.890 | 1 | 0.940 | 0.205 | 0.3 |
168 σ i (1300 nm) | 1 | 0.488 | 0.447 | 0.335 | ||
140 σ i (1300 nm) | 0 | 1 | 0.673 | 0.496 | 0.162 |
For a given molecule, the set of time constants is the same, independent of probe conditions. Only the coefficients of the exponentials depend on the probe wavelength and the ion mass observed. This is interpreted by population Li flowing consecutively from one observation window Li to the next. This flow is described by rate equations, whose solution is a sum of exponentials for each population (the same, independent of probe conditions). 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. They are simulation parameters and are also given in the table. The data are redundant for deducing the τi and the populations Li. 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. The very short τ1 is needed to reproduce the steep rise of the parent signals. τ2 causes, for instance, the delayed rise of the fragment as compared to the parent, clearly visible at 1300 nm. τ3 represents the decay of the initial spike, which is evident in particular at the longer probe wavelengths. 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. The decays in the range up to 600 ps (τ4, τ5) are not shown.
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![]() ![]() ![]() ![]() ![]() ![]() | (2) |
![]() | ||
Fig. 4 Time dependence of the anisotropy r![]() ![]() ![]() ![]() ![]() ![]() |
![]() | ||
Fig. 5 Potential energy surfaces and relaxation paths. The reaction coordinate q changes direction along the path, but probably everywhere contains the amino-group twist as an important component. The second coordinate spanning the branching space of the conical intersection (CI) must have b symmetry and probably contains the amino-group pyramidalization as a component. The τ3 equilibration follows a path around the lower cone. This diagram is used for both PYRBN and PIPBN. |
It is instructive to have a closer look to the ionic resonance. Energy differences of ion states are just equal to orbital energy differences in the neutral molecules. Orbitals and their energies are easy to understand by the help of their nodal pattern. The HOMO of PYRBN has the symmetry species b1 and is π-antibonding between donor and acceptor moieties (Fig. 6). The ground state of the parent ion (2B1) has a hole in this orbital (single occupancy). The HOMO−
1 has a node containing the amino group; it has the species a2. 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). 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. 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. There are not many coordinates with this property. 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. A strong pyramidalization of the amino group would probably be less effective, and planarization would do the opposite.
![]() | ||
Fig. 6 Bottom: π orbitals and energies of neutral PYRBN (a) and PBN (b)
versus twist angle θ as calculated by Parusel.20 Top: Ionic states resulting from a hole in the corresponding orbitals; for instance, the ground state (2B1) of the ion PYRBN+ has a hole in the HOMO (b1 orbital). Ionic excitation energies are equal to orbital energy differences. The suggested resonance with the probe laser is indicated. Symmetry species in C2v are only given for zero twist angle. For θ![]() ![]() |
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). The conclusion that the path leads through a CI can again be obtained from the shortness of τ2. The next time constant (τ3) corresponds to Lb↔
CT equilibration around the lower cone of the CI. As before, the much longer τ4 is ascribed to intersystem crossing, probably to Tn, which then can rapidly equilibrate with T1. From there, the molecule decays within τ5 to the ground state, from which ionization was too weak to be detected. (Any excited state would be detectable.) This process may be accompanied by decomposition such as H elimination, as suggested in PIPBN.9
There are not many vibrations with such a low wavenumber in this molecule. 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. 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. Hence, the observed oscillation should be assigned to the lower-frequency arching or to the twist vibrations. Arching involves the inversion motion of the amino nitrogen.
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°. 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
(As mentioned in the Introduction, a molecular-mechanics calculation predicted θ
=
15° in S0.12 This would imply also a minor barrier at θ
=
0°. It would not much change the following considerations.) 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 Therefore both types of vibration will occur between two turning points that are symmetry-equivalent. 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. A vibration (belonging to an antisymmetric species) in such a symmetric potential will therefore give rise to signal maxima twice in each period. (This is the equivalent to the spectroscopic selection rule Δv
=
2 for antisymmetric vibrations in electronic transitions. It will also apply for the 69-cm−1 structure in the PYRBN supersonic beam spectrum,11–13 if our assignment is correct.) Therefore we should divide by two to convert the observed period to wavenumbers:
ν![]() ![]() | (3) |
Also the La potential will be symmetric to θ=
0 and ω
=
0, again due to the partial double bond. Therefore these vibrations will not (or nearly not) be Franck–Condon active in the pump transition. They are thus not expected to be excited by the pump laser. Instead, they will be stimulated by acceleration on the last slope on entering the Lb well from the CI. The direction of this slope therefore contains either θ or ω as an important component. 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. 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 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). 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.
In this scheme, ionization from the CT state would leave a hole in the a2 orbital, giving an ion ground state 2A2. The positive anisotropy found (Fig. 4) indicates a resonance in the ion parallel to the molecular axis. We thus would need a 2A2→
2A2 transition in the ion or the corresponding transition between a pair of a2 orbitals. 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. On the other hand, the pair of b1 orbitals just does the required energy approach along the twist coordinate. 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). 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 We suppose that attaching a para-cyano substituent (resulting in PBN) will lower the 2A state more than the 2B state. In fact, another CASPT2 calculation reports for PBN at low energy only a 2A state with CT character.26 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. This implies that the same scheme for the surfaces and relaxation paths can be used as for PYRBN (Fig. 5).
In contrast to PYRBN, no oscillations were observed in PBN (although some irregularities in the signals might indicate that there are some). 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.
It may be worth repeating8–10 a comment on the amino-group twist as a component of the reaction coordinate. 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. A CASPT2 calculation for PBN in fact found such effects.26 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).
This journal is © The Royal Society of Chemistry 2004 |