Robert
Berger
Chemistry Department, Technical University Berlin, Strasse des 17. Juni 135, D-10623 Berlin, Germany. E-mail: Robert.Berger@www.chem.tu-berlin.de
First published on 25th November 2002
One of the greatest challenges in the field of molecular physics is the still outstanding measurement of the parity violating energy difference (ΔEpv) which electroweak theory predicts to exist between the enantiomers of chiral molecules. Apparently, previous attempts to detect signals of parity violation in molecules almost exclusively focused on electronic ground state systems. But as this study reveals, there are good reasons to believe that compounds which are chiral in an electronically excited state are particularly suited for a possible proof of molecular parity violation. Here, first calculations of the parity violating energy difference in electronically excited states are presented and routes to its measurement are discussed.
In the past, several experimental approaches have been proposed and various, but in the outcome unsuccessful or at least doubtful, attempts have been made to measure parity violating effects in chiral molecules. These approaches can generally be classified as to detect certain frequency shifts for example in the nuclear magnetic resonance, in the microwave, infrared or even Mössbauer-spectra of the enantiomers6–14 or to measure the parity violating coupling in a special class of molecules with intermediate tunneling15 (see also ref. 16) or directly the spectroscopic energy difference between the left- and the right-handed compound17 (see also ref. 5). While compared to a direct measurement of the energy difference between enantiomers the detection of line splittings appears to be in principle less difficult, the interpretation of those splittings in terms of ΔEpv is complicated.
The intensified efforts over the past few years to detect parity violating effects in molecules have largely been motivated by the fact that opposed to early computations (see for instance refs. 18–21) more recent calculations4,22–27 predict the absolute value of the parity violating energy difference often to be larger by an order of magnitude, which has profound consequences for possible experiments.
Surprisingly, it seems that all experimental schemes proposed so far aim at a detection of molecular parity violation in the electronic ground state and that parity violating potentials have only been calculated for this state, with the following two exceptions: In a figure caption of ref. 28 it has been noted that the role of the ground and excited state could be reversed in the experiment proposed in ref. 17 (see below) and Ivanov and Letokhov29 performed explicit calculations of molecular parity violating effects in electronically excited states, namely on splittings of electronically excited states in optically inactive molecules.
As will be detailed later in this article, a different class of systems, that is molecules which are optically active in electronically excited states, is particularly attractive for a measurement of ΔEpv for the following reasons: Firstly, parity violating potentials in electronically excited states might be larger in magnitude than those in the electronic ground state. Secondly, direct measurements of parity violating energy differences in these states can in principle be performed with a simplified scheme (see also caption to Fig. 1 in ref. 28). Thirdly, a large pool of molecular systems becomes accessible for a possible measurement of ΔEpv.
![]() | ||
Fig. 1 One-dimensional cut through a multi-dimensional Born-Oppenheimer potential energy surface V(![]() |
Thus, first calculations of parity violating energy differences in electronically excited states are presented in this paper and an experimental approach to measure them is discussed.
In the non-relativistic approximation the main contribution to the parity violating potential is expected to be due to the following one-electron operator (in atomic units)
![]() | (1) |
In order to obtain non-vanishing parity violating potentials, coupling between spin and coordinate space is essential. The most pronounced coupling effect is expected to be due to spin–orbit coupling, for which the corresponding operator reads in the Breit–Pauli form (in atomic units)
![]() | (2) |
Within the multi-configuration linear response (MCLR) approach, the parity violating potential Vpv is then given by:
Vpv![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() | (3) |
In traditional second order perturbation theory, the parity violating potential reads as
![]() | (4) |
Molecules of interest for a calculation of parity violating potentials in electronically excited states and for the type of experiment to be outlined below are for instance those carbonyl compounds XYCO which are planar in the electronic ground state and pyramidal in the lowest excited singlet state. For several of these compounds the n→
π* state is non-planar and, depending on the substituents, the barrier for stereomutation can be large. This in turn leads to a small tunneling splitting ΔE±, which is of advantage for an experimental detection of molecular parity violation (see below). Other suitable molecular systems are for instance formaldimine derivatives or 1,2-disubstituted ethylene with chiral equilibrium geometries in an electronically excited state.
In early spectroscopic studies of formyl fluoride,31–33 one of the simplest compounds of the type XYCO, it has been concluded that the first electronically excited singlet states of this molecule has a (chiral) pyramidal equilibrium structure as opposed to a planar nuclear arrangement in the electronic ground state. The corresponding double minimum potential in the first electronically excited singlet state was investigated in more detail in later studies34,35 and tunneling splittings in various vibronic states were reported. Stanton and Gauss36 in a computational study used the equation-of-motion coupled-cluster singles and doubles (EOM-CCSD) method to calculate the equilibrium structure of HFCO in the S1 state.
In the present study the complete active space-self consistent field (CASSCF) approach has been employed with all 18 valence electrons distributed in a set of 13 molecular orbitals described with the (atomic) augmented correlation consistent polarized valence triple-ζ basis set (usually referred to as aug-cc-pVTZ). Equilibrium structures and harmonic vibration wavenumbers obtained for the ground and the first electronically excited singlet state are shown in Table 1. Within the multi-configuration linear response approach to electroweak quantum chemistry the parity violating potential at the chiral equilibrium geometry of the lowest excited singlet state of formyl fluoride is computed to be −139×
10−20Eh for the R-enantiomer. The parity violating potential in the electronic ground state at the same geometry is +37
×
10−20Eh and thus significantly smaller in magnitude. These findings are in accordance with what might be expected from sum-over-states arguments. SOC values between the first excited singlet state and the lowest two triplet states are on the order of 2 cm−1 and 46 cm−1, while the reference state is energetically sufficiently well-removed from the triplet states (corresponding (de)excitation wavenumbers
>
2600 cm−1). Results for ethylene and 1,2-difluoroethylene, which will be presented in a forthcoming publication, indicate the same trend of an increase of the parity violating potential in electronic singlet states with close lying triplet states. Here, however, the effect for the first electronically excited state is less pronounced since 90°-twisted ethylene is a perfect biradical with a large energy gap between the S1 and the S0 state whereas the splitting between the S0 and the lowest triplet state is very small.
Parameter | S0 | S1 | ||
---|---|---|---|---|
CASSCF | CCSDa | CASSCF | EOM-CCSDa | |
aug-cc-pVTZ | cc-pVDZ | aug-cc-pVTZ | cc-pVDZ | |
a Ref. 36. | ||||
r(CO) | 118.75 | 118.74 | 135.51 | 134.62 |
r(CF) | 133.75 | 134.03 | 134.47 | 134.58 |
r(CH) | 108.44 | 109.97 | 107.88 | 109.77 |
∠OCF | 122.27 | 123.13 | 109.18 | 109.82 |
∠OCH | 127.64 | 127.35 | 114.70 | 116.10 |
τ | 0.0 | 0.0 | 46.51 | 46.3 |
ω 1 | 3206.85 | 3173.51 | 3254.24 | 3147.33 |
ω 2 | 1860.77 | 1923.85 | 1226.95 | 1217.57 |
ω 3 | 1457.77 | 1418.33 | 1421.64 | 1384.70 |
ω 4 | 1114.77 | 1126.65 | 1162.66 | 1149.10 |
ω 5 | 682.87 | 678.41 | 509.44 | 462.35 |
ω 6 | 1080.84 | 1046.72 | 1107.84 | 996.45 |
V pv(S0) | 0 | 37 | ||
V pv(S1) | 0 | −139 |
In the following section the principle idea of a direct measurement of ΔEpv is outlined and a scheme for an observation of ΔEpv in molecules that are chiral in electronically excited states is discussed.
In the absence of parity violation, the molecular Hamiltonian is parity conserving and thus commutes with the parity operator
([Ĥpc,
]
=
0). States |L〉 and |R〉 localized on the left- or the right-hand side respectively (henceforth referred to as chiral or handed states) and which are related by the parity operation (|L〉
=
|R〉) then have identical energy expectation values (Fig. 1, left hand side). However, they are not eigenstates of the parity-conserving Hamiltonian Ĥpc since those must as well be eigenstates of the parity operator and thus be either symmetric or antisymmetric with respect to inversion at the origin. The parity eigenstates |+〉 and |−〉, which are proportional to the plus- and the minus-combination of the chiral states, will in contrast to the handed states differ in energy as a result of the chiral state's tunneling through the barrier. In the high barrier limit, however, this energy difference is negligibly small.
Due to parity violation, the states of well-defined parity will be perturbed by the parity violating part Ĥpv of the total electroweak Hamiltonian Ĥ=
Ĥpc
+
Ĥpv. If we restrict ourselves to a two level system only and work in the basis of the parity eigenstates, the strength of the mixing depends on the tunneling splitting E± and on the parity violating coupling element Vpv
=
〈+|Ĥpv|−〉. Assuming real states for simplicity, we arrive at the following eigenvalue problem (see also refs. 15,17,5):
![]() | (5) |
![]() | (6) |
If the absolute value of the tunneling splitting ΔE±=
E+
−
E− is much smaller than the absolute value of the parity violating matrix element Vpv, we obtain (assuming Vpv
<
0)
D
≈
|2Vpv| as well as 1
±
ΔE±/D
≈
1. Consequently, the true eigenstates are the handed states (|1〉
≈
{|+〉
+
|−〉}/√2, |2〉
≈
{|+〉
−
|−〉}/√2) with E1
−
E2
≈
2Vpv. In the limit of an infinitely high barrier, E1
−
E2 corresponds to the parity violating energy difference ΔEpv. In the opposite case (assuming ΔE±, Vpv
<
0) we obtain D
≈
|ΔE±| and 1
−
ΔE±/D
≈
2, 1
+
ΔE±/D
≈
0. The true eigenstates are then the states of well-defined parity (|1〉
≈
|+〉, |2〉
≈
−|−〉) with E1
−
E2
≈
ΔE±.
Harris and Stodolski15 noted that when |2Vpv|≪
|ΔE±|, a system initially prepared in a chiral state |L〉 or |R〉 would show time-dependent oscillations of the optical activity θ with a period τ of h/ΔE±. The time-averaged optical rotation 〈θ〉t however would vanish. When |Vpv| is no longer negligibly small compared to |ΔE±|, the optical activity would oscillate with a reduced time period around a non-vanishing average value. The time-dependence of the quantum beats gives experimental information about the parity violating coupling and thus indirectly also on ΔEpv by reference to the tunneling splitting according to 〈θ〉t/θmax
=
V2pv/[V2pv
+
(ΔE±/2)2].
Quack17 suggested considering the reverse situation: When |ΔE±|≪
|2Vpv|, which should be the case for typical chiral molecules, ΔEpv becomes a directly measurable quantity. A system initially prepared in a state of well-defined parity would show oscillations between the parity eigenstates. The time-period τ for these oscillations is then determined by the parity violating energy difference, according to τ
=
h/ΔEpv. It has been proposed to create this state of well-defined parity in the electronic ground state of a typical chiral molecule with a large barrier for stereomutation via optical excitation to and deexcitation from an electronically excited state with an achiral equilibrium geometry where the parity of the rovibrational eigenstates is well-defined, since |ΔE±|
≫
|2Vpv|. Due to the strong parity selection rules in optical dipole transitions, this excitation–deexcitation scheme can be used to create a state of well-defined parity in an initially unpopulated rovibrational level of the electronic ground state of the chiral molecule. Influenced by the parity violating energy difference in this rovibrational state, a transformation to the state of opposing parity will take place which can be probed with a third transition to a state of well-defined parity in the electronically excited state of the molecule.
There are several difficulties connected with this approach (see below and ref. 17) which one, however, could overcome under favorable conditions. Nevertheless, although this proposal is 15 years old and although various aspects of this approach can now be considered to be under control, the preparation step still seems to be a problematic aspect, to a lesser extent in view of technical difficulties, but rather because promising molecular candidates that are ideally suited for this scheme are apparently still to be discovered. For several chiral molecules the parity violating energy difference has been calculated (see for instance refs. 19,20,23,37–47) and for S2Cl247 it has been explicitly shown that the parity violating energy difference exceeds the tunneling splitting in the lowest vibrational levels of the electronic ground state by orders of magnitude, but for most of those molecules we do not expect that electronically excited states with an achiral equilibrium geometry exist for which a rotationally resolvable and analysable spectrum may be obtained, a necessary prerequisite to the experiment.
Therefore, I propose to search for molecular parity violation in electronically excited states. Since on the one hand we significantly enlarge our pool of systems applicable in our quest for molecular parity violation and on the other hand might even profit from a gain in order of magnitude for the quantity we are aiming at, electronically excited states appear to be promising for a successful measurement of |ΔEpv|: As illustrated in Fig. 2, for a molecule with an achiral equilibrium geometry in the electronic ground state (more precisely: for which |ΔE±|≫
|2Vpv| holds in certain rovibrational levels of the electronic ground state) and with a chiral equilibrium geometry in an electronically excited state we can due to the strong parity selection rules directly prepare a state of well-defined parity in the electronically excited state via an optical dipole transition from the initial vibronic ground state. If |ΔE±| in this excited vibronic state is much smaller than |Vpv|, the time for parity change is almost exclusively determined by the parity violating energy difference. The time-dependent increase of the population p of the initially forbidden parity state, which is in the absence of decay processes given by p
=
sin2[ΔEpvt/(2ħ)], can then selectively be probed by a second transition to a level of well-defined parity, which may for instance be a rovibrational level of an achiral electronic state (not necessarily the electronic ground state) or a level of the electronically excited state in an energy region, where the tunneling splitting is large compared to |ΔEpv|. This scheme is in principle also applicable to states with different spin symmetry. Then, however, other selection rules apply. Note, that in this type of experiment it is not necessary to resolve the parity violating splitting since the strong parity selection rules for electric dipole transitions allow for a selective preparation and detection of one of the two parity eigenstates.
![]() | ||
Fig. 2 Schematic representation of a measurement of the parity violating energy difference for systems which are chiral in an electronically excited state (|ΔE±|![]() ![]() ![]() ![]() |
Compared to the experiment proposed in ref. 17 the preparation of a state of well-defined parity is greatly simplified, which is of particular importance, since up to now no definite molecular example has been presented that is according to explicit calculations of the tunneling splitting, the parity violating potential and either computations or experimental knowledge of the parity conserving potential of the intermediate electronic state qualified for the measurement suggested in ref. 17. While already in this reference it has been suggested, among other possibilities, to use allenes with achiral excited states for a spectroscopic experiment with this scheme, only very recently it has been calculated quantitatively in a different context that 1,3-dimethylallene has an achiral electronically excited state.48 With the approach outlined in the present paper (see also caption to Fig. 1 in ref. 28) we significantly increase the number of classes of molecular systems suitable for a parity violation experiment and can benefit from the presumably larger parity violating energy differences in electronically excited states.
These merits, however, come at the expense of further disadvantages specific for an experiment performed in electronically excited states, while other difficulties associated with the experiment outlined in ref. 17 (collisions; magnetic dipole, electric quadrupole transitions; electric fields etc.) which have already been examined more closely within that letter, continue to exist for the approach discussed here. The main complicating factors will be addressed briefly in the following:
(i) For an experiment conducted in the electronically excited state the most critical step will be the evolution phase. As has been discussed before, in the absence of decay processes the signal for parity change will evolve as p=
sin2[ΔEpvt/(2ħ)]. Given that the value of Vpv computed for formyl fluoride in its first electronically excited singlet state translates to a time for parity change of about 30 s, it appears to be extremely difficult to record an entire period of the parity oscillation. However, this is not necessary. Since sin2(x) behaves approximately as x2 for small x, it has been noted for the electronic ground state17 that a measurement of the quadratically increasing signal in the initial time evolution (typically in the milli- or microsecond regime46) would be sufficient to determine ΔEpv. This also holds for electronically excited states although here additional complicating factors arise. For an isolated molecule, the prepared metastable state will decay due to internal conversion, intersystem crossing, predissociation as well as fluorescence and these processes will compete with the dynamics of parity change induced by Vpv. In order to ensure a significant population of the electronically excited state over this time span (milli- to microseconds) and thus a reasonably intense signal, the radiative lifetime should be at least of the same order of magnitude. Accordingly, the natural fluorescence (or phosphorescence) lifetime has to be long enough which requires a sufficiently small oscillator strength for a radiative transition between the electronically excited state and the ground state (see, however, below). The criterion for long radiative lifetimes (larger than microseconds) may for instance be met by various systems in their lowest electronically excited state with the same spin symmetry as the electronic ground state, by chiral high Rydberg states and in particular by electronically excited states with different spin symmetry showing phosphorescence lifetimes typically in the range of milliseconds to hundreds of seconds.
Interactions between the chiral molecular system and its environment especially due to collisions, however, will modify the dynamics of parity change considerably. An abundant literature on this particular point exists (see for instance refs. 5,15,49–52), since the hypothesis has been raised that environmental influences stabilize handed states of ordinary chiral molecules with respect to their states of well-defined parity. This necessitates working under single molecule conditions in the experiment outlined above, which can to a large extent be realized within a molecular beam or possibly within specially designed molecular traps. Severe consequences on the dynamics are due to parity changing collisions that result in a loss of phase coherence (so-called pure dephasing). Following the estimate in ref. 17 we would obtain after 1 ms time of flight in the molecular beam a relative signal of about 3×
10−9
(here ΔEpv
=
2Vpv and Vpv
=
−139
×
10−20Eh as calculated for HFCO in its first electronically excited singlet state have been used), while for the traveling time of 100 ms used in ref. 17 the relative signal increases to 3
×
10−5. The probability for a collision within the 100 ms time of flight has been estimated as 2.5
×
10−7, which translates to 2.5
×
10−9 during the 1 ms flight. Even if these collisions have a high probability for parity change, which would lead to a significant background signal in the experiment, the signal caused by the incoherent collision events increases only linearly with time while the signal for coherent parity change displays quadratic dependence, so that the different contributions can be identified with additional measurements downstream the beam.17
Further dynamical information on the various decay channels including collisional deactivation of the electronically excited states can be obtained from control experiments which record the time evolution of the two states of well-defined parity after their direct preparation.
As has been noted in ref. 17, if a superselection rule exists for the handed states of chiral molecules (see ref. 53), this type of experiment can provide a sensitive check of this hypothesis.
(ii) The preparation phase has been significantly simplified compared to the experiment in ref. 17, since the state of well-defined parity is now prepared with a single electric dipole transition. However, excitations with other parity selection rules may occur due to accompanying magnetic dipole and electric quadrupole transitions. While these transitions are typically five or eight orders of magnitude weaker than electric dipole transitions, they nonetheless impose lower limits on the strength of the electric dipole transitions. A smaller electric dipole transition moment is on the one hand favorable as it provides long radiative lifetimes but on the other hand a larger oscillator strength is required to reduce the background as a result of transitions with other parity selection rules. A compromise between these two requirements appears to be well possible, since the population of the state of unwanted parity caused by accompanying transitions can be removed by “cleaning” the prepared state. This can be performed for instance via probing the initially forbidden component with a strongly allowed electric dipole transition immediately after the preparation step which may remove the states with contributions of the unwanted parity. This procedure may also be applied to avoid complications due to rovibronic coupling in these states with relatively high energy.
(iii) Calculations of parity violating potentials of chiral molecules in electronically excited states are in principle more difficult than for electronic ground state systems and this will also affect the interpretation of future experiments. The approach used here for the calculation of ΔEpv in electronically excited states marks the first step in the direction of reliable theoretical estimates.
As is evident from the discussion above, formyl fluoride itself is not well suited for this type of experiment since the tunneling splitting in the lowest vibrational level of the first electronically excited state is, at approximately 0.2 cm−1 (ref. 35), much larger than ΔEpv. Moreover the fluorescence lifetime ranges from 130 ns to <10 ns (ref. 35) and is therefore too small for this type of experiment. For the deuterated compound DFCO the fluorescence lifetime falls in the range of 400 ns to 800 ns (ref. 35) and we may expect a further increase for TFCO so that the microsecond regime might finally be reached. Nonetheless, according to preliminary calculations, the tunneling splitting in DFCO and TFCO are on the order of 10−2 cm−1 and consequently still orders of magnitude larger than the parity violating energy difference. A more favorable situation may however be realized in other compounds of the general type XYCO.
I shall note here that electronically excited molecular states are also of interest for other approaches to detect molecular parity violation, for instance in a measurement of the dispersion of optical rotation or of the electronic circular dichroism. One might also consider to search for parity violating signals in electronic transitions, where for both the initial and the final rovibronic state the parity violating potential prevails against the tunneling splitting. However, one would have to beat the natural linewidth in order to resolve possible linesplittings in electronic transitions of the enantiomers and even if one succeeds, frequency shifts due to molecular parity violation have to be distinguished from splittings of various other sources, a difficult task even in the infrared region as recent experiments indicate.54
This journal is © the Owner Societies 2003 |