# Triple resonance for a three-level system of a chiral molecule.

I. IntroductionChirality is a fascinating phenomenon; its influence prevails in a number of fields, as described in many books and articles. (1) At an early stage of quantum mechanics, Hund already noticed a serious paradox for this new theory, namely the optical activity would vanish in quantum mechanics for any molecule in a stationary state. (2,3) In order to solve this problem, he thought of a tiny perturbation term between enantiomers, which could mix the left-handed and right-handed molecules by the tunneling motion through a potential barrier caused by the perturbation, but the rate of the motion would be so slow that the optical activity would be preserved practically forever. In 1957 a new finding of parity violation was reported. (4,5) This result means that the R- and S-enantiomers are not degenerate in energy. Both Hund's perturbation and parity violation terms are quite small for almost all the "stable" chiral molecules, in other words, the Hamiltonian of a chiral molecule is almost completely symmetric with respect to inversion, and the measurement of optical rotation has so far been nearly the only way to differentiate enantiomers. Even highly sophisticated spectroscopic methods such as the one with highly stabilized carbon dioxide lasers as sources (6) have achieved the resolution to [DELTA]v/v = [10.sup.-13], which is still at least three orders of magnitude larger than theoretical estimates for the contribution of parity violation. (7-9)

In view of the importance of enantiomer differentiation, quite a large variety of spectroscopic methods have been proposed, but most of them have proven far from successful. It might be appropriate here to cite several attempts for enantiomer differentiation using molecular complexes combined with high-resolution spectroscopic methods. Le Barbu-Debus, Zehnacker-Rentien and their collaborators (10) employed UV and UV/IR double-resonance spectroscopy to record electronic and vibrational spectra of chiral complexes generated by supersonic expansion. Howard and his collaborators, (11,12) Xu and her associates, (13-15) and Caminati's group (16) have applied Fourier-transform microwave spectroscopy to molecular complexes consisting of chiral-achiral and chiral-chiral molecules.

In these circumstances I decided to examine closely chiral molecules of [C.sub.1] symmetry, namely molecules without any symmetry operation, in order to trace the real origin of molecular chirality and to detect its consequences on spectra of chiral molecules. The other group of chiral molecules belong to [C.sub.2] group, which consists of a [C.sub.2] operation, but I shall treat this group of molecules separately, because they give only one type of rotational transitions, due to the permanent dipole moment being aligned to the [C.sub.2] axis. On the other hand, as has been well known, [C.sub.1] chiral molecules exhibit all three a-, b-, and c-type rotational transitions. This is certainly very surprising, but not many people have appreciated the significance of this important fact. It has been well established that any electric-dipole transitions take place between two levels of opposite parity, because the dipole moment is of odd parity, and certainly this rule should also apply to any pure rotational transitions of molecules. We have, however, many examples of [C.sub.1] chiral molecules in the literatures, for which all a-, b-, and c-types of rotational transitions were observed and reported. Thus the observations on these chiral molecules remained a puzzle for us for a long time. We may wonder how we can understand the rotational spectra of [C.sub.1] chiral molecules. This is the start of the present research; I suspected that this puzzle was closely related to Hund's paradox. I keep in mind mainly pure rotational spectroscopy of molecules such as microwave spectroscopy as an experimental tool of the present study, because it is straightforward and yet amenable to various types of spectroscopic experiments such as double resonance quite easily. Of course, laser spectroscopy is very similar to it in many respects, and the present treatment may readily be extended to other wavelength regions, infrared, visible, and ultraviolet, and may thus be applicable to [C.sub.2] chiral molecules as well, for which two of the three types of transitions can be vibrational, for example.

As already stressed, one of the most crucial points of the present study lies in the characteristic feature of [C.sub.1] chiral molecules that all the three a-, b-, and c-types of rotational transitions take place simultaneously. Of course, they can be observed routinely one by one by using a conventional spectroscopic method, but such observations, even with highest resolution available at the present, have not provided us with any new information on molecular chirality. Obviously a fundamentally new approach is required to clarify the real essence of molecular chirality, and I have decided to focus the attention to the dynamical behavior of a three-level system of a [C.sub.1] chiral molecule by introducing a triple-resonance method, combined with time-dependent perturbation theory.

II. Time-dependent perturbation theory

We follow the formulation described by Landau and Lifshitz in their textbook, (17) and reproduce some essential part of the theory for the sake of convenience. The Hamiltonian of a three-level system is expressed as

H = [H.sub.0] + V(t). [1]

The eigenfuction of the unperturbed time-independent Hamiltonian [H.sub.0] in a stationary state k is given by [[PSI].sub.k.sup.(0)] with the eigenvalue [E.sub.k.sup.(0)]; [[PSI].sub.k.sup.(0)] = [[psi].sub.k.sup.(0)] exp(-[iE.sub.k.sup.(0)]t/h). When a time-dependent perturbation V(t) is added, the solution of the Schrodinger equation:

ih[partial derivative][PSI]/[partial derivative]t =([H.sub.0] + V (t))[PSI] [2]

may be expressed by a sum of [[PSI].sub.k.sup.(0)] over k: [PSI] = [summation].sub.k][a.sub.k](t)[[PSI].sub.k.sup.(0)]. Inserting [PSI] in Eq. [2] leads to the following equation for the coefficient [a.sub.k](t):

[ihda.sub.m]/dt = [[summation].sub.k] [V.sub.mk](t)[a.sub.k](t), [3]

where

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. [4]

We assume the perturbation periodic with the frequency [omega], i.e. of the form:

V(t) = F exp(-i[omega]t) + G exp(i[omega]t). [5]

Because V(t) is a Hermitian operator,

[G.sub.nm] = [F.sub.mn.sup.*] [6]

holds, and Eq. [4] is given by

[V.sub.kn](t) = [F.sub.kn] exp[i([[omega].sub.kn] - [omega])t]

+ [F.sub.nk.sup.*] exp[i([[omega].sub.kn] + [omega])t]. [4a]

III. Application to two- and three-level systems

First we review the result on a two-level system worked out by Landau and Lifshitz, (17) which will allow us to make an important comparison with the results on a three-level system given below.

(1) Two-level system. The frequency of the external perturbation [omega] is assumed to be very close to [[omega].sub.mn]: [E.sub.m.sup.(0)] - [E.sub.n.sup.(0)] = ([omega] + [epsilon]), where [epsilon] represents a small frequency shift. Retaining only terms depending on the low frequency [[omega].sub.mn] - [omega] results in

[ihda.sub.m]/dt = [F.sub.mn] exp(i[epsilon]t)[a.sub.n], [7a]

[ihda.sub.n]/dt = [F.sub.mn.sup.*] exp(-i[epsilon]t)[a.sub.m]. [7b]

Because these equations involve the second time derivative of [a.sub.m], they yield two independent solutions for the coefficients, which are given by

[a.sub.n] = A exp(i[[alpha].sub.1]t),

[a.sub.m] = -A(h[[alpha].sub.1]/[F.sub.mn.sup.*]) exp(i[[alpha].sub.2]t), [8a]

[a.sub.n] = B exp(-i[[alpha].sub.2]t),

[a.sub.m] = B(h[[alpha].sub.2]/[F.sub.mn.sup.*]) exp(-i[[alpha].sub.1]t), [8b]

where

[[alpha].sub.1] = -[epsilon]/2 + [OMEGA], [[alpha].sub.2] = [epsilon]/2 + [OMEGA], [OMEGA] = [[[([epsilon]/2).sup.2] + [([absolute value of [F.sub.mn]]/h).sup.2]].sup.1/2] [9]

and A and B are determined by the normalization condition. The system at time t is, in general, expressed by a linear combination of the two solutions Eqs. [8].

[FIGURE 1 OMITTED]

Suppose that the system lies in the m-th state [[PSI].sub.m.sup.(0)] at t = 0, then after t its wavefunction changes to

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] [10]

The absolute square of the coefficient of [[PSI].sub.n.sup.(0)] in this equation is the probability for the state being in the n-th state after t, which is given by

[W.sub.n [left arrow] m] = ([[absolute value of [F.sub.mn]].sup.2]/2[h.sup.2][[OMEGA].sup.2])[1 - cos(2[OMEGA]t)], [11]

which, at the exact resonance [epsilon] = 0, reduces to

[W.sub.n [left arrow] m] = (1/2)[1 - cos(2[absolute value of [F.sub.mn]]t/h)]. [11a]

This varies between 0 and 1 as the time passes with the periodicity of [pi]h/[absolute value of [F.sub.mn]] and, when multiplied by a damping factor exp(-t/[tau]) account for the effect of intermolecular collision with the average collision time interval [tau], gives a free induction decay signal function.

(2) Three-level system. This treatment on the two-level system has been extended in the present study to a three-level system 1, 2, and 3, with the level separations given in frequency unit by [[omega].sub.a] = 2 - 1, [[omega].sub.b] = 3 - 2, and [[omega].sub.c] = 3 - 1 (the suffixes a, b, and c of the frequencies imply the types of rotational transitions). Here without loss of generality, the three levels are assumed to be arranged, as shown in Fig. 1. The equations for the coefficients, which correspond to Eqs. [7] in the case of the two-level system, are

[ihda.sub.1]/dt = [F.sub.21.sup.*] exp(-i[[epsilon].sub.a]t)[a.sub.2] + [F.sub.31.sup.*] exp(- i[epsilon].sub.c]t)[a.sub.3], [12a]

[ihda.sub.2]/dt = [F.sub.21] exp(i[[epsilon].sub.a]t)[a.sub.1] + [F.sub.32.sup.*] exp(- i[[epsilon].sub.b]t)[a.sub.3], [12b]

[ihda.sub.3]/dt = [F.sub.31] exp(i[[epsilon].sub.c]t)[a.sub.1] + [F.sub.32] exp(i[[epsilon].sub.b]t)[a.sub.2], [12c]

where

[[omega].sub.21] = [[omega].sub.a] + [[epsilon].sub.a], [13a]

[[omega].sub.31] = [[omega].sub.c] + [[epsilon].sub.c], [13b]

[[omega].sub.32] = [[omega].sub.b] + [[epsilon].sub.b], [13c]

In the following, for the sake of simplicity, the exact resonance is assumed for the three transitions. The coefficient [a.sub.1] satisfies the following equation:

[d.sup.3][a.sub.1]/[dt.sup.3] + [([F.sub.21] * [F.sub.21] + [F.sub.31] * [F.sub.31] + [F.sub.32] * [F.sub.32])/[h.sup.2]] [da.sub.1]/dt

- i[([F.sub.31] * [F.sub.32][F.sub.21] + [F.sub.21] * [F.sub.32] * [F.sub.31])/[h.sup.3]][a.sub.1] = 0. [14]

We shall assume the following trial form for [a.sub.1]:

[a.sub.1] = [A.sub.1] exp(i[lambda]t)[15]

and insert it in Eq. [14]. Then we obtain a third-order equation for [lambda] given by

[[lambda].sup.3] - [D.sub.2.sup.2] [lambda] + [D.sub.3.sup.3] = 0, [16]

where

[D.sub.2.sup.2] = ([[absolute value of [F.sub.12]].sup.2] + [[absolute value of [F.sub.13]].sup.2] + [[absolute value of [F.sub.23]].sup.2])/[h.sup.2], [16a]

[D.sub.3.sup.3] = [[F.sub.13][F.sub.32][F.sub.21] + [([F.sub.13][F.sub.32][F.sub.21]).sup.*]]/[h.sup.3]. [16b]

The left-hand side of Eq. [16] is plotted in Fig. 2 as a function of [lambda], where [D.sub.2] is arbitrarily chosen to be 3.0 and [D.sub.3] to take three values 2.0, 0.0, and -2.0.

The properties of Eq. [16] may be summarized as follows:

(1) There are three roots for [lambda]: [[lambda].sub.1], [[lambda].sub.2], and [[lambda].sub.3], which are taken to satisfy the condition:

[[lambda].sub.1] [greater than or equal to] [[lambda].sub.2] [greater than or equal to] [[lambda].sub.3]. [17]

(2) The left-hand side of Eq. [16]: f([lambda]) = [[lambda].sup.3] - [D.sub.2.sup.2][lambda] [D.sub.3.sup.3] reaches extremes at [lambda] = + [[lambda].sub.x] (minimum) and [lambda] = -[[lambda].sub.x] (maximum) with [[lambda].sub.x] denoting [D.sub.2]/[square root of 3]. The values of f([lambda]) at these extremes are

f ([+ or -] [[lambda].sub.x]) = [??] (2/3[square root of 3])[D.sub.2.sup.3] + [D.sub.3.sup.3]. [18]

(3)-1 If [D.sub.3] > 0 and (2/3[square root of 3]) [D.sub.2.sup.3] > [D.sub.3.sup.3], there are three roots: [[lambda].sub.1] and [[lambda].sub.2] > 0, [[lambda].sub.3] < 0 and (2/3[square root of 3])[D.sub.2.sup.3] = [D.sub.3.sup.3], [[lambda].sub.1] = [[lambda].sub.2] = [D.sub.2]/[square root of 3], [[lambda].sub.3] < 0 and (2/3[square root of 3])[D.sub.2.sup.3] < [D.sub.3.sup.3], only one root exists: [[lambda].sub.3] < 0.

(3)-2 If [D.sub.3] = 0, [[lambda].sub.1] = -[[lambda].sub.3] = [D.sub.2] > 0, [[lambda].sub.2] = 0.

(3)-3 If [D.sub.3] < 0 and (2/3[square root of 3]) [D.sub.2.sup.3] > [absolute value of [D.sub.3.sup.3]], there are three roots: [[lambda].sub.1] > 0, [[lambda].sub.2] and [[lambda].sub.3] < 0 and (2/3[square root of 3])[D.sub.2.sup.3] = [absolute value of [D.sub.3.sup.3]], [[lambda].sub.1] > 0, [[lambda].sub.2] = [[lambda].sub.3] = - [D.sub.1]/[square root of 3] and (2/3[square root of 3])[D.sub.2.sup.3] < [absolute value of [D.sub.3.sup.3]], only the root is [[lambda].sub.1] > 0.

[FIGURE 2 OMITTED]

Equation [15] assumed for [a.sub.1] is then inserted in Eqs. [12b] and [12c] to derive the solutions for [a.sub.2] and [a.sub.3]:

[a.sub.1] = [A.sub.1] exp(i[[lambda].sub.1]t), [19a]

[a.sub.2] = [A.sub.1][[[lambda].sub.1][F.sub.21] - [F.sub.23][F.sub.31]/h]h/[- [[lambda].sub.1.sup.2] + [[absolute value of [F.sub.23]].sup.2]/[h.sup.2] exp(i[[lambda].sub.1]t), [19b]

[a.sub.3] = [A.sub.1][[[lambda].sub.1][F.sub.31] - [F.sub.32][F.sub.21]/h]h/[- [[lambda].sub.1.sup.2] + [[absolute value of [F.sub.23]].sup.2]/[h.sup.2] exp(i[[lambda].sub.1]t), [19c]

where [lambda] = [[lambda].sub.1] is chosen. The other roots: [[lambda].sub.2] and [[lambda].sub.3], if available, are inserted in Eqs. [19] in place of [[lambda].sub.1] to derive the corresponding [a.sub.1], [a.sub.2], and [a.sub.3], respectively.

IV. Time-dependent behavior of the 2 [left arrow] 1 transition under a triple-resonance condition

Suppose that the three-level system under consideration lies in [[PSI].sub.1.sup.(0)] at t = 0. Using the results Eqs. [19], the system at t is expressed as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] [20]

where

[C.sub.21] = [[lambda].sub.1][F.sub.21] - [F.sub.32] * [F.sub.31]/h]/h/[-[[lambda].sub.1.sup.2] + [F.sub.32] * [F.sub.32]/[h.sup.2]], [21a]

[C.sub.22] = [[lambda].sub.2][F.sub.21] - [F.sub.32] * [F.sub.31]/h/h/[-[[lambda].sub.2.sup.2] + [F.sub.32] * [F.sub.32]/[h.sup.2]], [21b]

[C.sub.23] = [[lambda].sub.3][F.sub.21] - [F.sub.32] * [F.sub.31]/h]/h/[-[[lambda].sub.3.sup.2] + [F.sub.32] * [F.sub.32]/[h.sup.2]], [21c]

[C.sub.31] = [[lambda].sub.1][F.sub.31] - [F.sub.32][F.sub.21]/h]/h/[-[[lambda].sub.1.sup.2] + [F.sub.32] * [F.sub.32]/[h.sup.2]], [21d]

[C.sub.32] = [[lambda].sub.2][F.sub.31] - [F.sub.32][F.sub.21]/h]/h/[-[[lambda].sub.2.sup.2] + [F.sub.32] * [F.sub.32]/[h.sup.2]], [21e]

[C.sub.33] = [[lambda].sub.3][F.sub.31] - [F.sub.32][F.sub.21]/h]/h/[-[[lambda].sub.3.sup.2] + [F.sub.32] * [F.sub.32]/[h.sup.2]], [21f]

The initial conditions at t = 0, which we assumed, correspond to

[A.sub.1] + [A.sub.2] + [A.sub.3] = 1, [22a]

[A.sub.1][C.sub.21] + [A.sub.2][C.sub.22] + [A.sub.3][C.sub.23] = 0, [22b]

[A.sub.1][C.sub.31] + [A.sub.2][C.sub.32] + [A.sub.3][C.sub.33] = 0. [22c]

Equations [22] are solved for [A.sub.i] as follows:

[A.sub.1] = [l.sub.1] ([[lambda].sub.2] - [[lambda].sub.3])/[LAMBDA], [23a]

[A.sub.2] = [l.sub.2] ([[lambda].sub.3] - [[lambda].sub.1])/[LAMBDA], [23b]

[A.sub.3] = [l.sub.3] ([[lambda].sub.1] - [[lambda].sub.2])/[LAMBDA], [23c]

where

[LAMBDA] = [[lambda].sub.1][[lambda].sub.2](- [[lambda].sub.1] + [[lambda].sub.2]) + [[lambda].sub.2] [[lambda].sub.3] (-[[lambda].sub.2] + [[lambda].sub.3]) + [[lambda].sub.3][[lambda].sub.1](- [[lambda].sub.3] + [[lambda].sub.1]), [24]

[l.sub.1] = -[[lambda].sub.1.sup.2] + [[lambda].sub.0.sup.2], [25a]

[l.sub.2] = -[[lambda].sub.2.sup.2] + [[lambda].sub.0.sup.2], [25b]

[l.sub.3] = -[[lambda].sub.3.sup.2] + [[lambda].sub.0.sup.2], [25c]

[[lambda].sub.0.sup.2] = [F.sub.32] * [F.sub.32]/[h.sup.2]. [26]

The probability of the state being in the state 2 at t is given by the absolute square of the coefficient of [[PSI].sub.2.sup.(0)](t) in Eq. [20], i.e.

[W.sub.2-1] = [[absolute value of [A.sub.1] [C.sub.21] exp (i [[lambda].sub.1]t) + [A.sub.2] [C.sub.22] exp (i[[lambda].sub.2]t) + [A.sub.3] [C.sub.33] exp (i[[lambda].sub.3]t)].sup.2], [27]

which, when Eqs. [21a], [21b], and [21c] are inserted, reduces to

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] [28a]

with the new definitions:

[P.sub.i] + [iQ.sub.i] = [[[lambda].sub.i][F.sub.21] - [F.sub.32] * [F.sub.31]/h]/h (i = 1, 2, 3). [28b]

We thus finally obtain

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] [29]

V. Implications for enantio-differentiation

We shall focus attention to the rotational motion of a chiral molecule. The perturbation Eq. [4] is then ascribed to the interaction of the molecule with the external microwave electric field through the dipole moment of the molecule, namely due to the Stark effect. Because the dipole moment is a vector, it is reversed in sign when one enantiomer is transformed to the other, which means that all the matrix elements [F.sub.ij] change sign. As shown in a previous section III, because [D.sub.3] as defined by Eq. [16b] changes its sign, the three characteristic frequencies [[lambda].sub.1], [[lambda].sub.2], and [[lambda].sub.3] will be transformed accordingly as follows:

[lambda].sub.1] [right arrow] - [[lambda].sub.3], [30a]

[lambda].sub.2] [right arrow] - [[lambda].sub.2], [30b]

[lambda].sub.3] [right arrow] - [[lambda].sub.2], [30c]

and [P.sub.i] + [iQ.sub.i], as defined by Eq. [28b], varies as follows:

[P.sub.1] + [iQ.sub.1] [right arrow] [P.sub.3] + [iQ.sub.3], [31a]

[P.sub.2] + [iQ.sub.2] [right arrow] [P.sub.2] + [iQ.sub.2], [31b]

[P.sub.3] + [iQ.sub.3] [right arrow] [P.sub.1] + [iQ.sub.1], [31c]

and hence

[P.sub.1.sup.2] + [Q.sub.1.sup.2] [right arrow] [P.sub.3.sup.2] + [Q.sub.3.sup.2], [31d]

[P.sub.2.sup.2] + [Q.sub.2.sup.2] [right arrow] [P.sub.2.sup.2] + [Q.sub.2.sup.2], [31e]

[P.sub.3.sup.2] + [Q.sub.3.sup.2] [right arrow] [P.sub.1.sup.2] + [Q.sub.1.sup.2], [31f]

[P.sub.1][P.sub.2] + [Q.sub.1][Q.sub.2] [right arrow] [P.sub.2][P.sub.3] + [Q.sub.2][Q.sub.3], [31g]

([P.sub.1][Q.sub.2] - [P.sub.2][Q.sub.1]) [right arrow] - ([P.sub.2][Q.sub.3] - [P.sub.3][Q.sub.2]), [31h]

[P.sub.2][P.sub.3] + [Q.sub.2][Q.sub.3] [right arrow] - [P.sub.1][P.sub.2] + [Q.sub.1][Q.sub.2], [31i]

([P.sub.2][Q.sub.3] - [P.sub.3][Q.sub.2] [right arrow] - ([P.sub.1][Q.sub.2] - [P.sub.2][Q.sub.1]), [31j]

[P.sub.3][P.sub.1] + [Q.sub.3][Q.sub.1] [right arrow] - [P.sub.1][P.sub.3] - [Q.sub.1][Q.sub.3], [31k]

([P.sub.3][Q.sub.1] + [P.sub.1][Q.sub.3]) [right arrow] - ([P.sub.3][Q.sub.1] - [P.sub.1][Q.sub.3]), [31l]

while the frequencies and their trigonometric functions in Eq. [29] are changed as follows:

[([[lambda].sub.2] - [[lambda].sub.3]).sup.2] [right arrow] [([[lambda].sub.1] - [[lambda].sub.2]).sup.2], [32a]

[([[lambda].sub.3] - [[lambda].sub.1]).sup.2] [right arrow] [([[lambda].sub.3] - [[lambda].sub.1]).sup.2], [32b]

[([[lambda].sub.1] - [[lambda].sub.2]).sup.2] [right arrow] [([[lambda].sub.2] - [[lambda].sub.3]).sup.2], [32c]

([[lambda].sub.2] - [[lambda].sub.3])([[lambda].sub.3] - [[lambda].sub.1]) [right arrow] ([[lambda].sub.3] - [[lambda].sub.1])([[lambda].sub.1] - [[lambda].sub.2]), [32d]

([[lambda].sub.3] - [[lambda].sub.1])([[lambda].sub.1] - [[lambda].sub.2]) [right arrow] ([[lambda].sub.2] - [[lambda].sub.3])([[lambda].sub.3] - [[lambda].sub.1]), [32e]

([[lambda].sub.1] - [[lambda].sub.2])([[lambda].sub.2] - [[lambda].sub.3]) [right arrow] ([[lambda].sub.1] - [[lambda].sub.2])([[lambda].sub.2] - [[lambda].sub.3]), [32f]

cos([[lambda].sub.1] - [[lambda].sub.2])t [right arrow] cos([[lambda].sub.2] - [[lambda].sub.3])t, [33a]

sin([[lambda].sub.1] - [[lambda].sub.2])t [right arrow] sin([[lambda].sub.2] - [[lambda].sub.3])t, [33b]

cos([[lambda].sub.2] - [[lambda].sub.3])t [right arrow] cos([[lambda].sub.1] - [[lambda].sub.2])t, [33c]

sin([[lambda].sub.2] - [[lambda].sub.3])t [right arrow] sin([[lambda].sub.1] - [[lambda].sub.2])t, [33d]

cos([[lambda].sub.3] - [[lambda].sub.1])t [right arrow] cos([[lambda].sub.3] - [[lambda].sub.1)t, [33e]

sin([[lambda].sub.3] - [[lambda].sub.1])t [right arrow] sin([[lambda].sub.3] - [[lambda].sub.1)t. [33f]

As a result, the three sine terms in Eq. [29] will change their signs from one enantiomer to the other. Therefore, we may conclude that we might have a chance to differentiate two enantiomers through the observation of their time-dependent rotational spectra under a triple resonance condition.

VI. Discussion and conclusion

VI-1. The [D.sub.3] constant. One of the most important results of the present study lies in the discovery of the [D.sub.3] constant, which, as defined by Eq. [16b], consists of products of three Stark matrix elements, each corresponding to a different-type transition between two of the three levels of the system, as illustrated in Fig. 1. This constant may thus be regarded to bear an essential feature of molecular chirality.

It should be pointed out that this constant is not compatible with conventional theory of molecular rotation (i.e. the rotation of a rigid body). The Stark effect for the molecular rotation is expressed by products of the components of the permanent dipole moment and the direction cosines between the dipole moment components and the space-fixed axis along which the external electric field is applied. For the pure rotation, the dipole moment is a constant, and the matrix elements of the three direction cosines are normally chosen such that two of them real and one imaginary. Therefore, the rotation theory concludes [D.sub.3] to be zero. We may generalize this treatment such that the three Stark matrices result in a real [D.sub.3], to be consistent with its definition, namely by choosing the Stark matrices to be either (i) all real or (ii) one real and two imaginary. Then we would arrive at an insignificant conclusion that [P.sub.i] and [Q.sub.i] defined by Eq. [28b] satisfy either [P.sub.i] [not equal to] 0 and [Q.sub.i] = or [P.sub.i] = 0 and [Q.sub.i] [not equal to] 0, making all the sine terms in Eq. [29] vanish. The present study indicates that the rotational levels of a [C.sub.1] chiral molecule do not belong to any definite parity, rather to a mixed parity, and hence they do not follow the theory of molecular rotation.

In deriving Eq. [29], we implicitly assumed that there are three solutions of [lambda]. This assumption, which is indispensable for chirality, looks plausible, namely the condition [D.sub.3] > 0 and (2/3[square root of 3])[D.sub.2.sup.3] > [D.sub.3.sup.3] or [D.sub.3] < 0 and (2/3[square root of 3] [D.sub.2.sup.3] hold in most cases.

VI-2. The mixed parity. Next we need clarify the meaning of parity mixing for the rotational levels of a [C.sub.1] chiral molecule. As mentioned in the Introduction, there are two terms known at the present time that give rise to molecular chirality or to parity mixing of rotational levels: Hund's perturbation and parity violation. Here I shall derive a formal expression for the mixed parity, although the two terms originate from backgrounds quite different in physical nature (see Ref. 3, pp. 207-213 for the physical basis of the parity violation). We designate the two terms (in the energy unit) by [E.sub.F] and [E.sub.PV], respectively, and write down the wavefunctions for the enantiomers: |R> and |S> in terms of the two basis sets of definite parity: |+> and |-> as follows:

|R> = cos [phi]|+> + sin [phi]|->, [34a]

|S> = sin [phi]|+> - cos [phi]|->, [34b]

where

cos [phi] = {([E.sub.T] + [E.sub.PV])/[(2[E.sub.T])}.sup.1/2], [35a]

sin [phi] = {([E.sub.T] - [E.sub.PV])/[(2[E.sub.T])}.sup.1/2], [35b]

[E.sub.T] = [[[E.sub.PV.sup.2] + [[E.sub.F.sup.2]].sup.1/2]. [36]

If we assume [E.sub.PV] susceptible to the inversion operation P as follows:

P [E.sub.PV] = -[E.sub.PV], [37]

we have

P cos [phi] = sin [phi], [38a]

P sin [phi] = cos [phi], [38b]

P|[+ or -]> = [+ or -]|[+ or -]>, [38c]

and arrive at

P|R> = sin [phi]|+> - cos [phi]|-> = |S), [39a]

P|S> = cos [phi]|+> + sin [phi]|-> = |R>, [39b]

as we expected. Note that the assignment of |R> and |S> is immaterial here, in other words, they may be interchanged.

IV-3. The parity mixing of rotational levels of a chiral molecule. The parity mixing expressed by Eqs. [34] is of great significance for enantiomer differentiation as represented by the sine terms in Eq. [29]. In order to make clear this point further, we examine the parity of pure rotational states in some detail.

In order to include the rotational motion in the discussion of chirality, we rewrite the inversion operation as follows:

P = [[sigma].sub.g'g"] [C.sub.2.sup.g], [40]

where g, g', and g" denote the three different principal inertial axes, a, b, and c. The first factor [[sigma].sub.g'g"] of Eq. [40] denotes the reflection on the g'g" plane and the second one the two-fold rotation about the g axis. We designate a rovibronic state as |+, r> or |-, r), where the vibronic part is of definite parity either D or -, as shown above, and r denotes the rotational part, which is one of the four group members for an asymmetric rotor: [K.sub.a][K.sub.c] = ee, eo, oo, and oe. (18) The transition probability for the [r.sub.1]-[r.sub.2] transition is the absolute square of the matrix element of the dipole moment component along the space-fixed Z axis [[mu].sub.Z].

Let us consider two simple examples of chiral molecules, in order to illustrate the role of the rotational motion. The first example is the isopropanol molecule: [(C[H.sub.3]).sub.2] CHOH, which has been confirmed to exist in two rotational isomers, trans and gauche, with respect to the dihedral angle between the C-O-H and the H-C-O planes. (19-23) There are two equivalent gauche forms, which may be regarded as a chiral pair, while we simply ignore trans, because it is well isolated from gauche. Hirota (24) and Hirota and Kawashima (25) found that the rotational spectra observed for gauche were well reproduced by assuming the a inertial axis antisymmetric, whereas the b and c axes symmetric, and established that the two gauche forms were split into a symmetric (+) and an antisymmetric (-) state with the energy separation of 46798.901(60) and 4431.4613(17) MHz, respectively, for [(C[H.sub.3]).sub.2] CHOH and [(C[H.sub.3]).sup.2] CHOD. Although the gauche forms of isopropanol are by no means optically stable in a conventional sense, as indicated by the large energy separations, they are of great use in discussing the role of rotational motion in chirality.

We may thus interpret the matrix elements of the dipole moment Z component of the gauche chiral pair of isopropanol using Eq. [40] as follows:

<[+ or 0], [r.sub.1] [[absolute value of [[mu].sub.Z]] [+ or -], [r.sub.2]> = [[mu].sub.b] <[r.sub.1] [[absolute value of [[PHI].sub.bZ]] [r.sub.2]'> + [[mu].sub.c] <[r.sub.1] [[absolute value of [[PHI].sub.cZ]] [r.sub.2]">, [41a]

<[+ or -], [r.sub.1] [[absolute value of [[mu].sub.Z]] [??], [r.sub.2]> = [[mu].sub.a] <[r.sub.1] [absolute value of [[PHI].sub.aZ]] [r.sub.2]'">, [41b]

where the ([r.sub.1], [r.sub.2]) pair in the ([r.sub.1] [absolute value of [[PHI].sub.bZ]] [r.sub.2]'> block, the matrix element of the direction cosine [[PHI].sub.bZ], is either (ee, oo) or (eo, oe), that in the <[r.sub.1] [absolute value of [[PHI].sub.bZ]]] [r.sub.2]"> block is either (ee, oe) or (oo, eo), and that in the <[r.sub.1] [[absolute value of [[PHI].sub.aZ]] [r.sub.2]"'> block is either (ee, eo) or (oo, oe). Eqs. [41] exactly provide a basis for the theoretical framework, which Hirota employed to analyze the rotational spectra of gauche isopropanol.

We then derive the matrix element of the dipole moment in the enantiomer basis, provided that the enantiomer is optically stable, i.e. localized in a limited area of phase space, by employing Eqs. [34]:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] [42a]

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] [42b]

<R, [r.sub.1] [[absolute value of [[mu].sub.Z]] S, [r.sub.2]> = ([E.sub.PV]/[E.sub.T]) [[mu].sub.a] <[r.sub.1] [[absolute value of [[PHI].sub.aZ]] [r.sub.2]"'>. [42c]

Although the relative magnitudes of [E.sub.F] and [E.sub.PV] are different for different molecules, we may estimate [E.sub.F]/[E.sub.T] to be of the order of one, in view of the fact that the parity violation term is as small as [10.sup.-13] [cm.sup.-1]. We may thus conclude that the rotational spectra of an enantiomer consist of all three a-, b-, and c-types transitions, with one component, a-type in the present example, opposite in sign for the two enantiomers, in agreement with what we found for [D.sub.3] and others. It is interesting to note that, when the parity violation is sizable compared with the Hund term, we may observe transitions between the two different enantiomers, as Eq. [42c] indicates.

The second example of the chiral molecule is propylene oxide, C[H.sub.3]CH(O)C[H.sub.2]. An optically pure sample of this molecule is available commercially, which means that the enantiomer of this molecule is optically stable, and in fact, all three types of rotational spectra were observed and analyzed in a conventional way; one trivial example of three transitions, which form a loop similar to that depicted in Fig. 1, is the c-type [1.sub.11]-[1.sub.01], the a-type [2.sub.12]-[1.sub.11] and the b-type [2.sub.12]-[1.sub.01] transitions. (26-28) The antisymmetric inertial axis of this molecule is the c axis, and thus all the results derived for isopropanol hold for propylene oxide as well, provided that the c axis plays the role in place of the a axis.

IV-4. Comments on enantio-differentiation experiments. It should be noted that the observation of time-dependent microwave spectra based upon Eq. [29], if successful, tells us only the difference between the enantiomers; the assignment of the observed spectra to either R-or S-enantiomers individually needs additional information. It should also be noted that the level scheme depicted in Fig. 1 and the enantiomer-dependent signal might be correlated with each other.

The prediction of enantio-differentiation based upon Eq. [29] looks reasonable and an experiment based on this prediction is certainly worth being carried out. It requires three microwave radiations of different frequencies, and care must be taken to estimate an optimum power for these radiations. In recent years we employ molecular beams in a routine way, which eliminate the burden of pressure broadening to a considerable extent. Furthermore, the advances in electronics allow us to use microwave sources with a wide variety of power and frequency, and thus no serious limitations exist in this respect in conducting triple resonance experiments. There are two factors to which we should pay attention: the off-resonance terms [[epsilon].sub.a], [[epsilon].sub.b], and [[epsilon].sub.c] in Eqs. [13] and the relative phases of the three radiations. The former is deliberately ignored in the present study primarily to simplify the theoretical development, but would play significant roles when we are interested in the line shape of the triple resonance signals, as are the case for double resonance, (29) where the intermolecular collisions will compete with the effects of multiple resonances. The phases of the three external radiations might affect triple resonance conditions to a considerable extent, but can be exploited to extract some further useful information from the triple resonance signals.

IV-5. The nature of "tunneling" in Hund's paradox. The mechanism leading to the tunneling in Hund's paradox has not all been clarified in detail. The chiral molecule is normally "defined" in such a way that one must break at least one chemical bond of this chiral molecule to transform one enantiomer to the other, and thus the transformation and/or the racemization will resemble a sort of chemical reaction, although the reaction paths have not been identified clearly for most chiral molecules. The gauche forms of isopropanol briefly discussed in IV-2 are not considered as chiral molecules in light of this conventional definition of molecular chirality. However, as I demonstrate, the gauche forms certainly behave like chiral molecules, but the tunneling splitting is quite large and thus enantiomers are almost impossible to isolate one from the other. The tunneling motions in these cases are thought to approximate some known internal motions. In the example of isopropanol, the internal rotation about the CH-OH axis has been suspected to be the main "reaction" coordinate. However, it should be pointed out that the energy differences between the two gauche forms determined for the OH and OD species of isopropanol do not satisfy the relation expected from a simple internal-rotation model, (25) indicating close examination of the "reaction coordinate" indispensable to clarify the real mechanism of tunneling in this molecule.

Presumably because the tunneling in Hund' s paradox remains ambiguous in nature, several investigations have recently been carried out, and here I cite a study of Trost and Hornberger (30) as such an example. They treated a specific molecule [D.sub.2][S.sub.2] by assuming the dominant collisional decoherence mechanism to stabilize the configuration states, i.e. the enantiomer states, of the molecule. They stated that the chirally sensitive parts of the London dispersion interaction in their treatment consist of (ED-EQ)(ED-ED) van der Waals interaction, where ED denotes the electric dipole and EQ the electric quadrupole transitions, and that the associated dispersion term is proportional to [r.sup.-7], rather than of the ordinary [r.sup.-6] dependence on the distance r between the interacting atoms and/or molecules.

The mechanism proposed by Trost and Hornberger can play an important role in the tunneling problem. It should be pointed out, however, that most asymmetric chemical syntheses so far employed are based on chiral catalyzers, and we certainly have to carry out much more extensive researches in order to clarify how significant the collisional decoherence is in enormous number of chiral catalytic materials.

doi: 10.2183/pjab.88.120

(Received Aug. 30, 2011; accepted Jan. 30, 2012)

References

(1) See, for example, Gardner, M. (1990) The New Ambidextrous Universe. Symmetry and Asymmetry from Mirror Reflections to Superstrings, 3rd revised ed., W. H. Freeman and Co., New York; Japanese translation (1992) translated by Tsuboi, C., Fujii, A. and Kojima, H., Kinokuniya, Tokyo.

(2) Hund, F. (1927) Zur Deutung der Molekelspektren. III. Bemerkungen uber das Schwingungs- und Rotationsspektrum bei Molekeln mit mehr als zwei Kernen. Z. Phys. 43, 805-826.

(3) Barron, L.D. (2004) Molecular Light Scattering and Optical Activity. 2nd ed., revised and enlarged, Cambridge University Press, Cambridge.

(4) Lee, T.D. and Yang, C.N. (1956) Question of parity conservation in weak interactions. Phys. Rev. 104, 254-258.

(5) Wu, C.S., Ambler, E., Hayward, R.W., Hoppes, D.D. and Hudson, R.P. (1957) Experimental test of parity conservation in beta decay. Phys. Rev. 105, 1413-1415.

(6) Daussy, Ch., Marrel, T., Amy-Klein, A., Nguyen, C.T., Borde, Ch.J. and Chardonner, Ch. (1999) Limit on the parity nonconserving energy difference between the enantiomers of a chiral molecule by laser spectroscopy. Phys. Rev. Lett. 83, 1554-1557.

(7) Quack, M. (2002) How important is parity violation for molecular and biomolecular chirality? Angew. Chem. Int. Ed. 41, 4618-4630.

(8) Berger, R. and Stuber, J.L. (2007) Electroweak interactions in chiral molecules: two-component density functional theory study of vibrational frequency shifts in polyhalomethanes. Mol. Phys. 105, 41-49.

(9) Figgen, D. and Schwerdtfeger, P. (2009) Structures, inversion barriers, and parity violation effects in chiral SeOXY molecules (X, Y = H, F, Cl, Br, or I). J. Chem. Phys. 130, 054306 and earlier references cited therein.

(10) A most recent publication is Le Barbu-Debus, K., Broquier, M., Mahjoub, A. and Zehnacker-Rentien, A. (2008) Chiral recognition between [alpha]- Hydroxylesters: A double-resonance IR/UV study of the complexes of methyl mandelate with methyl glycolate and methyl lactate. J. Phys. Chem. A 112, 9731-9741 and earlier references cited therein.

(11) King, A.K. and Howard, B.J. (2001) A microwave study of the hetero-chiral dimer of butan-2-ol. Chem. Phys. Lett. 348, 343-349.

(12) Hearn, J.P.I., Cobley, R.V. and Howard, B.J. (2005) High-resolution spectroscopy of induced chiral dimers: A study of the dimers of ethanol by Fourier transform microwave spectroscopy. J. Chem. Phys. 123, 134324.

(13) Su, Z., Borho, N. and Xu, Y. (2006) Chiral self-recognition: Direct spectroscopic detection of the homochiral and heterochiral dimers of propylene oxide in the gas phase. J. Am. Chem. Soc. 128, 17126-17131.

(14) Borho, N. and Xu, Y. (2007) Molecular recognition in 1:1 hydrogen-bonded complexes of oxirane and trans-2,3-dimethyloxirane with ethanol: a rotational spectroscopic and ab initio study. Phys. Chem. Chem. Phys. 9, 4514-4520.

(15) Borho, N. and Xu, Y. (2007) Lock-and-key principle on a microscopic scale: The case of the propylene oxide ... ethanol complex. Angew. Chem. Int. Ed. 46, 2276-2279.

(16) Maris, A., Giuliano, B.M., Bonazzi, D. and Caminati, W. (2008) Molecular recognition of chiral conformers: A rotational study of the dimers of glycidol. J. Am. Chem. Soc. 130, 13860-13861.

(17) Landau, L.D. and Lifshitz, E.M. (1967) Quantum Mechanics; Japanese translation (translated by Sasaki, K. and Koumura, S.) Tokyo Tosho Co., Tokyo, Chap. 6.

(18) See, for example, Gordy, W. and Cook, R.L. (1984) Microwave Molecular Spectra. 3rd Ed., John Wiley & Sons, New York, Chap. 7.

(19) Tanaka, C. (1962) Vibrational spectra of isopropyl alcohol and its deuterated species [(C[H.sub.3]).sub.2] CHOD, [(C[D.sub.3]).sub.2]CHOH, [(C[H.sub.3]).sub.2] CHOD. Nippon Kagaku Zasshi 83, 521-528.

(20) Tanaka, C. (1962) Normal vibrations of isopropyl alcohol. Nippon Kagaku Zasshi 83, 657-660.

(21) Tanaka, C. (1962) Infrared absorption spectra of deuterated isopropyl alcohol [(C[H.sub.3]).sub.2] CDOH, [(C[H.sub.3]).sub.2] CDOD. Nippon Kagaku Zasshi 83, 661-667.

(22) Kondo, S. and Hirota, E. (1970) Microwave spectrum and internal rotation of isopropyl alcohol. J. Mol. Spectrosc. 34, 97-107.

23) Inagaki, F., Harada, I. and Shimanouchi, T. (1973) Far-infrared spectra and barriers to internal rotation of isopropyl alcohol and alkyl mercaptans. J. Mol. Spectrosc. 46, 381-396.

(24) Hirota, E. (1979) Internal rotation in isopropyl alcohol studied by microwave spectroscopy. J. Phys. Chem. 83, 1457-1465.

25) Hirota, E. and Kawashima, Y. (2001) Internal rotation of the hydroxyl group in isopropanol and the chirality of the gauche form: Fourier transform microwave spectroscopy of [(C[H.sub.3]).sub.2] CHOD. J. Mol. Spectrosc. 207, 243-253.

(26) Swalen, J.D. and Herschbach, D.R. (1957) Internal barrier of propylene oxide from the microwave spectrum. I. J. Chem. Phys. 27, 100-108.

(27) Herschbach, D.R. and Swalen, J.D. (1958) Internal barrier of propylene oxide from the microwave spectrum. II. J. Chem. Phys. 29, 761-776.

(28) Hirota, E. and Kawashima, Y. unpublished.

(29) See, for example, Field, R.W., Hirota, E., Maier, J.P. and Tsuchiya, S. (eds.) (1998) Nonlinear Spectroscopy for Molecular Structure Determination. Blackwell Science, Oxford, UK.

(30) Trost, J. and Hornberger, K. (2009) Hund's paradox and the collisional stabilization of chiral molecules. Phys. Rev. Lett. 103, 023202.

By Eizi HIROTA * [1], ([dagger])

(Communicated by Toshimitsu YAMAZAKI, M.J.A.)

* [1] The Graduate University for Advanced Studies, Hayama, Kanagawa, Japan.

([dagger]) Correspondence should be addressed: E. Hirota, The Graduate University for Advanced Studies, Hayama, Kanagawa 240-0193, Japan (e-mail: ehirota@triton.ocn.ne.jp).

Printer friendly Cite/link Email Feedback | |

Author: | Hirota, Eizi |
---|---|

Publication: | Japan Academy Proceedings Series B: Physical and Biological Sciences |

Article Type: | Report |

Geographic Code: | 9JAPA |

Date: | Mar 1, 2012 |

Words: | 7105 |

Previous Article: | Functional cDNA expression cloning: pushing it to the limit. |

Next Article: | Two distinct pathways for essential metabolic precursors for isoprenoid biosynthesis. |

Topics: |