Laser-Control of Ultrafast π-Electron Ring Currents in Aromatic Molecules : Roles of Molecular Symmetry and Light Polarization

Being motivated by the recent progress in attosecond laser technology, we theoretically explore the strategy of inducing ultrafast electron dynamics inherent to aromatic molecules, i.e., ring currents by means of polarized laser pulses. The main topic of discussion is how to control the direction of ring currents in an aromatic molecule of low symmetry, for which the design of an efficient control pulse cannot be achieved intuitively. We first consider a system with a single aromatic ring and show that coherent π-electron angular momentum, which oscillates with time, can be produced and controlled by a polarized laser pulse with its ellipticity and orientation properly chosen. Nonadiabatic couplings with molecular vibration gradually weaken the angular momentum, while the vibrational amplitude strongly depends on the polarization of incident light. This suggests the conversion of the polarization dependence of ring current into that of subsequent vibration, which may open a way to detect laser-driven ultrafast electron dynamics by vibrational spectroscopy. The laser-control scheme for the ring current is then extended to a molecule with two aromatic rings, which exhibits characteristic phenomena absent in that with a single ring. We demonstrate that two-dimensional switching of the direction of angular momentum is possible in such molecules. In addition, ring current can be localized at a specific ring by tailored lasers. The application of the present control method to polycyclic aromatic hydrocarbons will lead to the development of next-generation organic optical switching devices.


Introduction
Chemical reaction (i.e., chemical bond formation and cleavage) can be viewed as a drastic change in nuclear positions, accompanied by electron transfer or migration between atoms in matter.In photochemical reactions, nuclear motion is initiated by the generation of a coherent superposition of vibrational eigenstates, i.e., a quantum vibrational wave packet (WP) with light.In general, the energy spacing between vibrational eigenstates is in the infrared (IR) region.Hence, nuclei pass through a transition state of photochemical reactions within tens or hundreds of femtoseconds and the process of the passage can be directly measured by means of femtosecond laser pulses [1].Continuous progress in laser technology has been promoting this research field called femtosecond chemistry or femtochemistry [2][3][4][5], which aims at real-time observation of chemical reactions, since the 1980s.
Currently, attosecond pulse generation [6][7][8][9][10][11][12][13][14], which was achieved in the early 2000s, is opening a door to unexplored possibilities of chemistry and physics.Attosecond light pulses have been expected to be a powerful tool not only to observe but also to drive electron motion in molecules, which is much faster than nuclear motion and thus steers chemical reactions [15,16].In analogy to nuclear dynamics, the fundamental mechanism of laser-induced ultrafast electron dynamics is to create a coherent superposition of electronic eigenstates, i.e., an electronic WP.The bandwidth as well as the central frequency of an attosecond laser pulse produced by high harmonic generation [17][18][19][20][21][22] lies in the extreme ultraviolet (XUV) region, which enables simultaneous transitions from the electronic ground state to many excited states.The electronic wave function of the superposition state varies in time periodically according to the energy gaps between the eigenstates involved, giving rise to large-scale electron motion in molecules.The period of the electronic state change is only a few femtoseconds for visible-light energy gaps or even of an attosecond order for UV energy gaps, where nuclei are almost static.The instantaneous and dramatic change in molecular electronic states caused by attosecond laser irradiation may lead to the control of molecular functionalities and chemical reactions, which is one of the ultimate goals of chemistry.The development of the new discipline "attosecond chemistry" will bring innovation to various areas, such as materials science and nanotechnology.
The theoretical proposal for ultrafast charge migration was put into practice in 2014 for phenylalanine by pump-probe spectroscopy [48].A sub-300-as XUV pump pulse triggers ionization and subsequent charge migration; a 4-fs visible/near-IR probe pulse then generates an immonium dication fragment whose yield exhibits a 4-fs oscillation as a sign of the charge migration.To the best of our knowledge, this is the first experimental application of attosecond laser pulses to charge migration in molecules.Shortly thereafter, control of charge migration was demonstrated for pre-oriented iodoacetylene HCCI by using intense near-IR pulses (instead of attosecond XUV pulses) [49,50].The linear polarization direction chosen perpendicular or parallel to its molecular axis affects the superposition of cationic ground and excited states created by multiphoton ionization and thus governs the fate of hole propagation.This experiment tells us that laser control of ultrafast electron dynamics in molecular systems is becoming reality.
There are two key factors to attain a desired linear combination of electronic eigenstates: molecular symmetry and light polarization.This can be readily understood from the fact that the probability of an optical transition between two states is dominated by the inner product of the corresponding transition dipole moment, whose direction depends on molecular symmetry, and the electric field vector of incident light.For instance, it is impossible to populate an excited state whose transition moment from the ground state is perpendicular (orthogonal) to the polarization vector.Manipulating the populations of eigenstates that participate in the superposition state is not enough; the relative quantum phase between the eigenstates is crucial, that is, coherent control is required.The initial value of the relative quantum phase can in principle be tuned by laser polarization for a spatially oriented molecule and then it evolves periodically until the molecule starts to vibrate.The initial relative phase thereby dictates electron motion inside a molecule, e.g., the timing of electron localization.The choice of polarization is relatively easier and often made intuitively for highly symmetric molecules such as linear ones; the experiment of charge migration in HCCI + is a typical example (although the control was not fully coherent) [49,50].However, many systems of interest (e.g., complex organic or biologically relevant molecules) have low or no symmetry.The strategy of designing a polarized laser pulse to control electron dynamics in such molecules is unclear and needs to be established.
In the present review article, we focus on the circulation of π electrons along an aromatic ring, i.e., the so-called ring current.This phenomenon, where π electrons travel in either a clockwise or anticlockwise direction around the ring, is a model case to study the manipulation of directional electron motion in complex molecules.Numerical simulations showed that much stronger ring currents can be produced in benzene [51] and Mg porphyrin [52][53][54] by a circularly polarized laser pulse than by a static magnetic field available with present technology.For these molecules the direction of electron flow is determined by the helicity of circular polarization.The mechanism of laser-driven ring currents in the highly symmetric aromatic molecules (D 6h for benzene and D 4h for Mg porphyrin) is unambiguous as will be briefly described in Section 2.1.How can ring currents be induced and controlled in aromatic molecules of lower symmetry by polarized light?
The rest of this article offers a solid answer to the above question obtained in a series of our theoretical efforts.In Section 2, our quantum WP simulations reveal the laser-polarization effects on nonadiabatically coupled ring current and molecular vibration.We first propose a general pulse design scheme for ring currents in molecules with a single aromatic ring, which is applicable to systems of both high and low symmetry.Numerical examples are presented to illustrate the characteristic behavior of controlled ring currents in 2,5-dichloropyrazine (C 2h ), which is a lower-symmetry derivative of heterocyclic pyrazine (D 2h ).Nonadiabatic coupling between ring current and molecular vibration is significant when the current lasts as long as vibrational periods (tens of femtoseconds); however, the coupling was completely ignored in the previous studies [51][52][53][54] based on the frozen-nuclei approximation.In Section 3, the pulse design scheme is extended to molecules with plural (in particular, two) aromatic rings.The extended theory is applied to (P)-2,2 -biphenol (C 2 ), which is a nonplanar system with axial chirality.We demonstrate electron dynamics inherent to multi-ring molecules, i.e., localization of ring current at a specific ring.Finally, Section 4 concludes the article.

Laser Control of Ring Currents in a Molecule Having a Single Aromatic Ring
In this section, we present the results of theoretical analyses and numerical calculations on laser-driven ring currents in a molecule with a single aromatic ring as the simplest system.To quantify the strength and direction of ring currents, the concept of electronic angular momentum eigenstates is essential since in general angular momentum serves as a measure of circular motions of a particle.We begin this section by proposing the basic idea to generate ring currents by polarized lasers from a viewpoint of the relationship between molecular symmetry and angular momentum eigenstates.

Basic Idea of Laser-Driven Ring Currents: Molecular Symmetry and Angular Momentum Eigenstates
First of all, we briefly summarize the mechanism of ring currents in benzene excited by a circularly polarized laser pulse [51] as an example.Benzene belongs to the D 6h point group at its equilibrium geometry and its C 6 axis is chosen as the z axis.Its highest occupied and lowest unoccupied molecular orbitals (HOMO and LUMO) are e 1g and e 2u orbitals, respectively.As for multielectron states constructed from MOs, there exist doubly degenerate 1 E 1u excited states whose major components are the HOMO-LUMO excitations.By approximating the hexagonal structure of benzene with a complete cylindrical ring, its symmetry becomes D ∞h and the z component of electronic angular momentum is quantized in units of the Dirac constant h.The doubly degenerate 1 E 1u excited states can be regarded as the eigenstates of the angular momentum operator lz with the quantum number m = ±1.[55][56][57][58] We denote the 1 E 1u states with m = +1 and −1 by 1 1 E 1u+ and 1 1 E 1u− , respectively.The angular momentum eigenstates 1 1 E 1u± have a complex wave function that can be expressed as a linear combination of two real ones 1 1 E 1ux and 1 1 E 1uy : (1) Recalling that complex atomic orbitals (AOs) 2p +1 and 2p −1 are angular momentum eigenstates of an electron in a hydrogen atom, the relation in Equation ( 1) is similar to that between the complex AOs and real ones 2p x and 2p y .The transition moments 1 1 A 1g | μ|1 1 E 1ux and 1 1 A 1g | μ|1 1 E 1uy , where 1 1 A 1g is the ground state and μ is the electric dipole moment operator, have the same magnitude and are perpendicular to each other in the xy plane.The x and y axes can thus be taken parallel to 1 1 A 1g | μ|1 1 E 1ux and 1 1 A 1g | μ|1 1 E 1uy , respectively.When a circularly polarized laser pulse propagates along the z axis, the spin angular momentum of a photon selects either 1 1 E 1u+ or 1 1 E 1u− and accordingly π electrons travel clockwise or anticlockwise around a ring depending on the sign of m.This is the origin of the unique correspondence between the direction of ring currents and the helicity of circular polarization.One can immediately find it impossible to rotate π electrons in benzene by a linearly polarized laser pulse, which has no spin angular momentum.We next assume that the symmetry of benzene is lowered, e.g., by introducing functional groups and/or replacing some carbon atoms in its ring with heteroatoms.The degeneracy between 1 1 E 1u+ and 1 1 E 1u− ( 1 1 E 1ux and 1 1 E 1uy ) is then broken and they split into a pair of quasi-degenerate states.We denote the real wave functions of the states closer to 1 1 E 1ux and 1 1 E 1uy by |X and |Y , respectively.Even though there exists no excited state that is an exact eigenstate of lz in such a system, sufficiently short laser pulses can prepare a coherent superposition of the optically allowed quasi-degenerate states as mentioned in Section 1.The approximate angular momentum eigenstates can be defined in analogy to Equation (1) as follows: where the matrix elements ± lz ± are close to ±h.After a laser pulse prepares either |+ or |− , it subsequently evolves in time as a coherent nonstationary state because of the nonzero energy gap between |X and |Y : where Ĥ0 is the field-free electronic Hamiltonian, ω X (ω Y ) is the angular frequency of |X (|Y ), and 2∆ω ≡ ω Y − ω X .Note that ∆ω is positive (negative) when |X is lower (higher) in energy than |Y .The approximate angular momentum eigenstates can be transiently created within the period of the electronic state change T ≡ π/|∆ω|.This means that transient ring currents can be achieved by selective generation of an approximate angular momentum eigenstate.In Section 2.2, we will propose the strategy for generating predominantly either |+ or |− by a polarized laser field.

Laser-Polarization Effects on Ring Currents within the Frozen-Nuclei Approximation
We analyze the effects of light polarization on an optical excitation of a molecule with a single aromatic ring using a three-level model of the ground state and two excited states.Here, the frozen-nuclei approximation, which ignores all nuclear degrees of freedom, is adopted.We analytically and numerically show that the coherent superposition of the quasi-degenerate states and the resultant angular momentum can be controlled by the ellipticity and orientation of the polarization ellipse of an incident laser.

General Form of Polarized Laser Pulses
Before entering into the analysis of optical excitation processes in aromatic molecules, let us describe the mathematical form of a classical laser field of arbitrary polarization.It is given in the dipole approximation by: where t is time, E 0 is the peak field strength, ω is the central frequency, ϕ is the optical phase, and e is the complex polarization unit vector.The envelope function f (t) varies between zero and unity for [0, t f ] with t f being the pulse duration; otherwise f (t) = 0.The polarization vector e for an arbitrarily polarized light is: e = e −iδ (cos β)e +1 + e iδ (sin where e +1 and e −1 are the spherical unit vectors corresponding to positive and negative helicities (spins), respectively.When the laser field propagates in the z direction and thus oscillates in the xy plane, the spherical unit vectors can be defined as follows: where e x and e y are the unit vectors in the x and y directions, respectively (in some references the signs ∓ are placed before 2 −1/2 ).Following the convention in optics, we refer to e +1 and e −1 as left and right circular polarizations, respectively.In Equation ( 5), δ denotes the orientation angle of the major axis of the polarization ellipse with respect to the x axis and β is termed the ellipticity angle in this article (Figure 1).The range of β is [0, π/2] and the minor-to-major axial ratio of the polarization ellipse is |tan (π/4 − β)|; for example, β = 0, π/4, and π/2 designate left circular, linear, and right circular polarizations, respectively.The general form of E(t) in Equation ( 4) will be utilized in the model analysis presented in Section 2.2.2.
where t is time, E 0 is the peak field strength, ω is the central frequency, φ is the optical phase, and e is the complex polarization unit vector.The envelope function f(t) varies between zero and unity for [0, tf] with tf being the pulse duration; otherwise f(t) = 0.The polarization vector e for an arbitrarily polarized light is: where e+1 and e−1 are the spherical unit vectors corresponding to positive and negative helicities (spins), respectively.When the laser field propagates in the z direction and thus oscillates in the xy plane, the spherical unit vectors can be defined as follows: where ex and ey are the unit vectors in the x and y directions, respectively (in some references the signs ∓ are placed before 2 −1/2 ).Following the convention in optics, we refer to e+1 and e−1 as left and right circular polarizations, respectively.In Equation ( 5), δ denotes the orientation angle of the major axis of the polarization ellipse with respect to the x axis and β is termed the ellipticity angle in this article (Figure 1).The range of β is [0, π/2] and the minor-to-major axial ratio of the polarization ellipse is |tan (π/4 − β)|; for example, β = 0, π/4, and π/2 designate left circular, linear, and right circular polarizations, respectively.The general form of E(t) in Equation ( 4) will be utilized in the model analysis presented in Section 2.2.2.Modified from [56] with permission of the American Chemical Society.

Three-Level Model Analysis of Optical Excitation
We formulate the polarization effects on the optical excitation in aromatic molecules of low symmetry.The time-dependent Schrödinger equation for a molecule in a laser field E(t) is given in the length gauge by: Modified from [56] with permission of the American Chemical Society.

Three-Level Model Analysis of Optical Excitation
We formulate the polarization effects on the optical excitation in aromatic molecules of low symmetry.The time-dependent Schrödinger equation for a molecule in a laser field E(t) is given in the length gauge by: The system is initially in its ground state |G .We assume that only the ground and quasi-degenerate excited states are involved in the excitation process.The energy of |G is defined to be zero and the central frequency ω of the laser field is set on resonance with the average energy of |X and |Y : ω = ω X + ∆ω = ω Y − ∆ω.The state vector |Ψ(t) is then expanded as: All diagonal elements of μ are assumed to be zero for simplicity, although this assumption is not necessary for our model analysis to be valid.For off-diagonal elements, we denote µ n ≡ G| μ|n = n| μ|G = 0 (n = X and Y), while X| μ|Y = Y| μ|X = 0.As defined in Section 2.1, the z axis is perpendicular to the ring plane, which coincides with the xy plane.Because the excitations to |X and |Y correspond to ππ * transitions, both µ X and µ Y lie in the xy plane.Substituting Equations ( 4) and (8) into Equation ( 7) yields the equation of motion for the expansion coefficient vector where and X and Y) are the complex Rabi frequencies.The initial condition is C(0) = (1, 0, 0) T .To simplify Equation (10a,b), we apply the rotating-wave approximation (RWA) [59]: Near resonance (i.e., ω >> |∆ω|) is assumed so that the contribution of the first terms in Equation (10a,b), which are rapidly oscillating exponentials, to the time evolution of C(t) is averaged out and removed.Under the RWA, Equation ( 9) is transformed into: In order to efficiently prepare an approximate angular momentum eigenstate |+ or |− defined in Equation (2), the quasi-degenerate states |X and |Y should be populated equally.The laser field is therefore assumed to satisfy where θ is the relative phase between can be fulfilled with the ellipticity angle β and orientation angle δ determined by the absolute value and argument of the right-hand side of Equation ( 12), respectively.To render Equation ( 11) more compact, we introduce an alternative pair of superposition states: where One can easily find that they are normalized and orthogonal to each other.We rewrite Equation (8) in terms of these superposition states as: where Then, the time evolution of the new coefficient vector D(t) ≡ (c G (t), c + (t), c − (t)) T follows: with the initial condition D(0) = (1, 0, 0) T .

Case of Degenerate Excited States
Before proceeding to the solution of Equation ( 16) for a quasi-degenerate system, let us consider the degenerate case of ∆ω = 0 for highly symmetric aromatic molecules such as benzene and Mg porphyrin.In this case, Equation ( 16) becomes simple, namely: The system practically consists of the two states |G and |U + .It is easy to integrate Equation (17): where F(t) ≡ t 0 dt f (t ).This indicates that the so-called π pulse [59], which satisfies ΩF t f = π, achieves complete population inversion from |G to |U + .Inserting Equation (18) into Equation (14)  yields the following: The relative quantum phase between |X and |Y coincides with that between the Rabi frequencies Ω X and Ω Y , i.e., θ, which is independent of time and controllable by the ellipticity angle β and orientation angle δ.A desired superposition of the degenerate states can be created by tuning the polarization of incident light properly.The populations of |+ and |− , which are the exact eigenstates of lz in this case, are given by: and the expectation value of lz is: The sign of l z (t) (in other words, the direction of ring currents) depends on the relative quantum phase θ and remains unchanged throughout the time evolution: Ring current flows in a sole direction even after a laser pulse vanishes at t = t f .
We revert to Equation (12), which connects the relative quantum phase θ to the ellipticity angle β and orientation angle δ.As stated for benzene, where |G = 1 1 A 1g , |X = 1 E 1ux , and |Y = 1 E 1uy , in Section 2.1, µ X and µ Y of highly symmetric aromatic molecules have the same magnitude and are perpendicular to each other: ||µ X || = ||µ Y || and µ X • µ Y = 0.The x and y axes can thus be chosen to be parallel to µ X and µ Y , respectively.As a result, Equation ( 12) is reduced to: which involves neither transition moments nor spherical unit vectors and is purely imaginary.Given a value of β, there exist two solutions of Equation ( 22): The major axis of the polarization ellipse with δ = π/4 bisects the angle between µ X and µ Y ; that with δ = − π/4 is perpendicular to it.In both cases, the value of l z (t f ) for the π-pulse excitation is hsinθ = hcos2β.We present illustrative examples and the case of linear polarization (β = π/4) is the first one; Equation ( 5) is then rewritten as: e = (cos δ)e x + (sin δ)e y (24) and δ determines the polarization direction.The Rabi frequencies Ω X and Ω Y are thus real-valued and their relative phase θ takes either zero or π.Henceforth, we denote the linear polarization vectors for θ = 0 (e iθ = +1) and π (e iθ = −1) by e (+) and e (−) , respectively.Do not confuse them with the spherical unit vectors e +1 and e −1 .The values of δ for e (+) and e (−) are π/4 and −π/4, respectively, from Equation (23).As in Equation (19), a linearly polarized laser pulse with e (+) [e (−) ] creates an in-phase (out-of-phase) superposition |X + |Y (|X − |Y ), which is an equal mixture of |+ and |− : P + (t) = P − (t) and thus l z (t) = 0 for all t.It is impossible to trigger ring currents in a degenerate system by a linearly polarized laser pulse.
Next, for circular polarization, β = 0 (π/2) provides e = e −iδ e +1 (e iδ e −1 ) and δ acts as an additional optical phase, which can be set to any value since Equation (19) tells us that the optical phase does not influence the superposition of the degenerate states.From the definition of θ as an argument of Ω Y /Ω X , we obtain θ = ± (η Y − η X ) for e ±1 , where η X (η Y ) is the angle of µ X (µ Y ) with respect to the x axis; θ = ±π/2 (e iθ = ±i) in a degenerate system.The superposition states |U ± are therefore identical to the exact angular momentum eigenstates |± (|∓ ) for left (right) circular polarization, while optical excitation occurs only to |U + .Ring current flows in the direction inherent to the eigenstate generated as noted in Section 2.1.
The phase factor e iθ yielded by elliptical polarization is neither real nor purely imaginary regardless of whether δ = π/4 or −π/4.The resultant behavior of π electrons is intermediate between those for linear and circular polarizations: The direction of ring currents is determined by the more populated of |+ or |− , while the magnitude of l z (t f ) is less than h even for the π-pulse excitation.

Case of Quasi-Degenerate Excited States
Now, we move on to the quasi-degenerate case of ∆ω 0 for aromatic molecules of low symmetry.It is difficult to analytically integrate Equation ( 16) with an arbitrary envelope function f (t).For this reason, f (t) is restricted to a rectangular function: f (t) = 1 for [0, t f ] and otherwise zero.During irradiation we have: and its solution is derived as: The generalized Rabi frequency Ω ≡ Ω 2 + (2∆ω) 2 1/2 is the root mean square of the Rabi frequencies and the detuning frequency 2∆ω.The pulse area of a laser field for a quasi-degenerate system is defined as ΩF(t f ), which is equal to Ωt f in the rectangular-envelope case.The ground state |G can be fully emptied by a rectangular-envelope pulse that meets 2|∆ω| ≤ Ω and has a pulse area of 2arccos − 2∆ω/Ω 2 .The pulse area of this laser field (termed the full-excitation pulse in this article) is larger than that of the π pulse, that is, 2arccos − 2∆ω/Ω 2 > π (but not larger than 2π); however, in fact, the full-excitation pulse is very close to the π pulse for ∆ω 0. The electronic WP of a quasi-degenerate system is eventually written down as: where and ϑ(t) ≡ 2argα(t).The relative quantum phase between |X and |Y is θ − ϑ(t) and its temporal behavior is the main subject of discussion in Section 2.2.The populations P ± (t) are derived as: and the angular-momentum expectation value l z (t) is where λ ≡ + lz + = − − lz − (> 0).In sharp contrast to Equation ( 21), the time-dependent phase ϑ(t) may change the direction of ring currents during irradiation.
For aromatic molecules of low symmetry, where ||µ 22) and ( 23) do not strictly hold.In this case, the linear polarization vectors e (+) (θ = 0) and e (−) (θ = π) are determined by [55,57,[60][61][62]: which fulfills ±Ω X = Ω Y .Once µ X and µ Y are given, the values of δ for e (+) and e (−) can be readily evaluated from Equation (31).It is straightforward to find a pair of β and δ for the other values of θ from Equation (12).Circular polarizations do not exactly satisfy Equation (12) unless the magnitudes of µ X and µ Y happen to be equal.If ∆ω 0, the relationship among β, δ, and θ should be close to that for degenerate systems; for example, e (+) and e (−) should be almost perpendicular to each other, and we approximately have ±iΩ X Ω Y (θ ± π/2) for e ±1 .
The electronic WP propagates freely after a laser pulse is turned off at t = t f .Similarly to Equation (3), the relative quantum phase in the free propagation is θ − ϑ(t f ) − 2∆ω(t − t f ) and thereby the coherent superposition of the quasi-degenerate states oscillates with the period T. The approximate angular momentum eigenstates |+ and |− are alternately created as predicted in Section 2.1.As a consequence, the direction of ring currents switches between clockwise and anticlockwise in marked distinction from the degenerate case, in which the current direction is unchanged and l z (t) is constant after the end of the laser pulse.
It should be emphasized that the initial relative phase θ of the superposed quasi-degenerate states (i.e., the initial direction of ring currents) can be controlled by the ellipticity angle β and orientation angle δ of the applied laser.Generation of polarization-shaped femtosecond UV laser pulses has become experimentally realizable [63][64][65][66].The efficient scheme to produce circularly polarized attosecond XUV laser pulses has been theoretically proposed [67].In Section 2.2.5, we will compare this model analysis in the rectangular-envelope case with numerical electronic WP simulations for pulse excitation with smooth rise and decay.

Numerical Demonstration for 2,5-Dichloropyrazine
The model system with a single aromatic ring adopted for numerical simulations is 2,5-dichloropyrazine (hereafter abbreviated as DCP), whose molecular formula is illustrated in Figure 2a.The molecule is assumed to be pre-oriented, e.g., by the nonadiabatic optical alignment technique [68,69].We executed ab initio electronic structure calculations for DCP with the 6-31G* Gaussian basis set [70] by using the quantum chemistry program MOLPRO [71].The geometry of DCP was optimized in the ground state at the level of the second-order Møller-Plesset perturbation theory (MP2) [72] and the optimized geometry was of C 2h symmetry.The excited-state properties at this geometry were obtained by the single-point calculation at the state-averaged complete-active-space self-consistent field (CASSCF) [73,74] level of theory.The active space was composed of ten electrons distributed among eight (four a u and four b g ) orbitals [abbreviated as CASSCF (10,8)].At the optimized geometry in the ground state |G = 1 1 A g , there exists a pair of optically allowed quasi-degenerate excited states |X = 3 1 B u and |Y = 4 1 B u with the energy gap 2h∆ω = 0.44 eV.The difference between their average energy and the ground-state energy is 9.62 eV.The magnitudes of the transition moments µ X and µ Y are 1.66ea 0 and 1.58ea 0 , respectively, where e is the elementary charge and a 0 is the Bohr radius.The directions of µ X and µ Y are illustrated in Figure 2b and the angle between them is η Y − η X = 0.35π.The linear combinations of |X and |Y defined by Equation (2) are the approximate angular momentum eigenstates |+ and |− with λ = 0.98h.Positive (negative) angular momentum designates π electrons traveling anticlockwise (clockwise) around the ring in Figure 2a.We performed electronic WP simulations for linear (β = π/4) and circular (β = 0 and π/2) polarizations.The directions of the linear polarization vectors e (+) and e (−) evaluated from Equation (31) are depicted in Figure 2b.For circular polarizations (e+1 and e−1), we set δ = 0. Instead of a rectangular envelope, we assumed a sin 2 function: for [0, tf] and otherwise zero.
The common laser parameters for the above four types of excitations were tf = 7.26 fs, ω = 9.62 eV/ħ [vacuum UV (VUV) light with a wavelength of 129 nm], and φ = 0.In analogy to the full-excitation pulse in the rectangular-envelope case, the peak field strength E 0 was chosen so that the pulse area ΩF(tf) (equal to Ωtf/2 for a sin 2 envelope) was . The values of E 0 so determined for e (+) , e (−) , e+1, and e−1 were 5.99, 9.80, 7.22, and 7.22 GVm −1 , respectively.Since ||µX||  ~ ||µY|| for DCP, the condition  approximately holds for circular polarizations.For the four types of excitations, we numerically integrated Equation ( 9) without the RWA by the fourth-order Runge-Kutta method [75].The time increment for the numerical integration was 0.024 fs (2.4 as).
Figure 3 displays the temporal behavior of lz(t) for the four types of excitations.The timing of the initial increase in the magnitude of lz(t) and its sign [i.e., the phase in the oscillation of lz(t)] manifestly depend on laser polarization.In Figure 3a, the temporal changes in lz(t) for e (+) and e (−) are identical with each other except for opposite signs: For the former (latter), π electrons start to circulate clockwise (anticlockwise), while in both cases the amplitude of lz(t) reaches almost ħ at t = 5.4 fs designated by the vertical broken line in Figure 3a.This is consistent with the difference in the initial relative phase θ by π between e (+) and e (−) predicted in the model analysis.The magnitudes of lz(t) for circular polarizations in Figure 3b grow earlier than for linear ones.Besides, a peak of lz(t) for e+1 (e−1) indicated by the vertical broken lines in Figure 3b appears at t = 7.0 fs (4.4 fs), when a laser pulse is almost off (near the maximum intensity).The peak positions of lz(t) in Figure 3b, which oscillate with different phases from those for e (+) and e (−) , are fairly consistent with θ = ± (ηY − ηX) = ± 0.35π for e±1; the 1.6-fs difference for e+1 from the peak positions at t = 5.4 fs in Figure 3a agrees perfectly with that in the field-free condition, 0.35π/(2∆ω) = 1.6 fs, while the difference for e−1 is shorter (1.0 fs) because We performed electronic WP simulations for linear (β = π/4) and circular (β = 0 and π/2) polarizations.The directions of the linear polarization vectors e (+) and e (−) evaluated from Equation (31) are depicted in Figure 2b.For circular polarizations (e +1 and e −1 ), we set δ = 0. Instead of a rectangular envelope, we assumed a sin 2 function: f (t) = sin 2 πt/t f for [0, t f ] and otherwise zero.The common laser parameters for the above four types of excitations were t f = 7.26 fs, ω = 9.62 eV/h [vacuum UV (VUV) light with a wavelength of 129 nm], and ϕ = 0.In analogy to the full-excitation pulse in the rectangular-envelope case, the peak field strength E 0 was chosen so that the pulse area ΩF(t f ) (equal to Ωt f /2 for a sin 2 envelope) was 2arccos − 2∆ω/Ω 2 under the constraint 2|∆ω| ≤ Ω.The values of E 0 so determined for e (+) , e (−) , e +1 , and e −1 were 5.99, 9.80, 7.22, and 7.22 GVm −1 , respectively.Since ||µ X || ||µ Y || for DCP, the condition |Ω X | |Ω Y | approximately holds for circular polarizations.
For the four types of excitations, we numerically integrated Equation ( 9) without the RWA by the fourth-order Runge-Kutta method [75].The time increment for the numerical integration was 0.024 fs (2.4 as).
Figure 3 displays the temporal behavior of l z (t) for the four types of excitations.The timing of the initial increase in the magnitude of l z (t) and its sign [i.e., the phase in the oscillation of l z (t)] manifestly depend on laser polarization.In Figure 3a, the temporal changes in l z (t) for e (+) and e (−) are identical with each other except for opposite signs: For the former (latter), π electrons start to circulate clockwise (anticlockwise), while in both cases the amplitude of l z (t) reaches almost h at t = 5.4 fs designated by the vertical broken line in Figure 3a.This is consistent with the difference in the initial relative phase θ by π between e (+) and e (−) predicted in the model analysis.The magnitudes of l z (t) for circular polarizations in Figure 3b grow earlier than for linear ones.Besides, a peak of l z (t) for e +1 (e −1 ) indicated by the vertical broken lines in Figure 3b appears at t = 7.0 fs (4.4 fs), when a laser pulse is almost off (near the maximum intensity).The peak positions of l z (t) in Figure 3b, which oscillate with different phases from those for e (+) and e (−) , are fairly consistent with θ = ± (η Y − η X ) = ± 0.35π for e ±1 ; the 1.6-fs difference for e +1 from the peak positions at t = 5.4 fs in Figure 3a agrees perfectly with that in the field-free condition, 0.35π/(2∆ω) = 1.6 fs, while the difference for e −1 is shorter (1.0 fs) because of the strong laser intensity.These agreements with the model analysis confirm the controllability of the initial relative phase θ of the superposed quasi-degenerate states by light polarization.As anticipated from Equation (30), the magnitude of l z (t) turns to decrease or even the direction of ring currents is reversed before the laser pulse ceases at t = t f = 7.26 fs.Finally, for all the four types of excitations, the total population in |X and |Y , P X (t f ) + P Y (t f ), reaches 0.995 and thus the amplitude of l z (t) achieves 0.97h.This indicates that a smoothly rising and decaying pulse with a pulse area of 2arccos − 2∆ω/Ω 2 is able to excite almost all the population in |G .At t > 7.26 fs, l z (t) oscillates with the period of T = π/∆ω = 9.4 fs as expected.The oscillation of l z (t) can be prevented to achieve unidirectional ring currents by a pump-dump control method [55,57,60].The application of this method to a molecule with two aromatic rings will be presented in Section 3.2.
anticipated from Equation ( 30), the magnitude of lz(t) turns to decrease or even the direction of ring currents is reversed before the laser pulse ceases at t = tf = 7.26 fs.Finally, for all the four types of excitations, the total population in X and Y , PX(tf) + PY(tf), reaches 0.995 and thus the amplitude of lz(t) achieves 0.97ħ.This indicates that a smoothly rising and decaying pulse with a pulse area of is able to excite almost all the population in G .At t > 7.26 fs, lz(t) oscillates with the period of T = π/∆ω = 9.4 fs as expected.The oscillation of lz(t) can be prevented to achieve unidirectional ring currents by a pump-dump control method [55,57,60].The application of this method to a molecule with two aromatic rings will be presented in Section 3.2.Such a strong VUV laser pulse as used in Figure 3 may induce multiphoton excitations to higher excited states or ionizations but the contribution of these additional processes is not large [55,57,60] The conclusions of this article are unaffected by the use of weaker laser pulses, although a smaller population is transferred to the quasi-degenerate excited states.

Beyond the Frozen-Nuclei Approximation: Nonadiabatically Coupled Vibronic Dynamics
Without nuclear motions, ring current continues even after the end of a laser pulse until a spontaneous emission relaxes the molecule; however, nonadiabatic coupling with molecular vibration must be taken into account to clarify the actual duration of the current.It is also very Such a strong VUV laser pulse as used in Figure 3 may induce multiphoton excitations to higher excited states or ionizations but the contribution of these additional processes is not large [55,57,60] The conclusions of this article are unaffected by the use of weaker laser pulses, although a smaller population is transferred to the quasi-degenerate excited states.

Beyond the Frozen-Nuclei Approximation: Nonadiabatically Coupled Vibronic Dynamics
Without nuclear motions, ring current continues even after the end of a laser pulse until a spontaneous emission relaxes the molecule; however, nonadiabatic coupling with molecular vibration must be taken into account to clarify the actual duration of the current.It is also very interesting how laser polarization affects vibrational motions of the molecule through ring currents.For this purpose, we next carried out nuclear WP simulations on effective potential energy surfaces (PESs) with respect to selected vibrational modes of DCP.

Selection of Effective Vibrational Modes and Method of Nuclear Wave Packet Propagation
To pick out effective vibrational degrees of freedom, excited-state geometry optimization was executed for |X and |Y at the CASSCF(10,8) level of theory.The optimized geometries of |X and |Y are of C 2h symmetry as that in |G .Hence, the displacements from the optimized geometry of |G to that of |X and |Y are totally symmetric.What is more, vibrational modes that nonadiabatically couple two 1 B u states |X and |Y are totally symmetric A g modes as well.We therefore selected two A g normal modes with large potential displacements and calculated the relevant nonadiabatic coupling matrix element between |X and |Y ; they are the breathing and distortion modes, whose vibrational vectors are illustrated in Figure 2c,d, respectively.The harmonic wave numbers of the breathing and distortion modes in |G are 1160 and 1570 cm −1 , respectively.We omitted nonadiabatic couplings between the ground and two excited states because |X and |Y lie far above |G and there exists no potential crossing between them near the Franck-Condon region.The adiabatic PESs of |X and |Y with respect to the normal coordinates Q of the breathing and distortion modes were evaluated at the CASSCF(10,8) level of theory.As will be shown in Section 2.3.4,there exists an avoided crossing (not a conical intersection) between the PESs, although whether it is an avoided crossing or a conical intersection does not affect the qualitative discussion below.At the crossing point, the gap between the two PESs is about 190 cm −1 .We also calculated the PESs at the level of the second-order CAS perturbation theory (CASPT2) [76,77] and confirmed that the geometry at the crossing point remains unchanged when dynamical electron correlation is included, while the PESs are lowered by about 3 eV.
The initial nuclear WP is the vibrational ground-state wave function of |G .The system is then electronically excited by a laser pulse E(t) defined by Equation ( 4).The subsequent propagation of the WP was carried out in the diabatic basis.Rigorous adiabatic-to-diabatic transformation requires the line integral of the derivative coupling matrix, which is in general path dependent [78,79].We instead used the quasi-diabatization scheme implemented in MOLPRO, which is based on an analysis of configuration interaction vectors [80].The total wave function of DCP was then expanded in terms of the three diabatic states, each of which is a linear combination of the adiabatic states |G , |X and |Y .The diabatic WPs ψ D n (Q, t) can be propagated by solving the coupled equations [81]: where ∇ 2 is the Laplace operator with respect to Q. V D nn (Q) are the diabatic potentials (n = n ) and couplings (n = n ) and µ D nn (Q) are the transition moments between the two diabatic states.The coupled equations were numerically integrated by the split-operator method [82] with the fast Fourier transform algorithm.The resultant diabatic WPs ψ D n (Q, t) were converted to adiabatic WPs ψ n (Q, t).For each of the breathing and distortion modes, the domain [−1.60 u 1/2 a 0 , 1.55 u 1/2 a 0 ] with u being the unified atomic mass unit was divided into 64 grid points at intervals of 0.05 u 1/2 a 0 so as to represent nuclear WPs.The time increment for the WP propagation was the same as used in Figure 3.

Population Transfer Dynamics
We employed the same laser pulses as used in Figure 3.The temporal change in P X (t) and P Y (t) for e (+) , e (−) , e +1 , and e −1 are plotted in Figure 4.When the laser pulse fully decays at t = t f = 7.26 fs, P X (t f ) + P Y (t f ) reaches 0.84, 0.93, 0.89, and 0.91 for e (+) , e (−) , e +1 , and e −1 , respectively.The polarization-dependent reductions in P X (t f ) + P Y (t f ) from those under the frozen-nuclei approximation (0.995) are attributed to the coordinate dependence of µ X (Q) and µ Y (Q) because the laser pulses used were designed with the transition moments at the optimized geometry of |G .
We employed the same laser pulses as used in Figure 3.The temporal change in PX(t) and PY(t) for e (+) , e (−) , e+1, and e−1 are plotted in Figure 4.When the laser pulse fully decays at t = tf = 7.26 fs, PX(tf) + PY(tf) reaches 0.84, 0.93, 0.89, and 0.91 for e (+) , e (−) , e+1, and e−1, respectively.The polarizationdependent reductions in PX(tf) + PY(tf) from those under the frozen-nuclei approximation (0.995) are attributed to the coordinate dependence of μX(Q) and μY(Q) because the laser pulses used were designed with the transition moments at the optimized geometry of G .The subsequent behaviors of PX(t) and PY(t) are totally different among the four types of excitations.For e (+) (Figure 4a), a small portion of the population shifts from X to Y by t ~ 10 fs and then a downward population transfer takes place around t ~ 10-14 fs by nonadiabatic transition.
In the case of e (−) (Figure 4b), a substantial amount of the population is transferred from Y to X and consequently   X Pt is more than seven times as large as   Y Pt at t ~ 10 fs.Afterwards, the direction of population transfer flips periodically with the rather small fraction of the population transferred.In contrast, for e+1 (Figure 4c), a part of the population is continuously exchanged between X and Y and, thus, electronic relaxation is completed at t > 30 fs.The behaviors of PX(t) and PY(t) for e−1 (Figure 4d) are intermediate between those for e (+) and e (−) : A considerable population transfer from Y to X is observed around t ~ 8-13 fs.The very different evolutions of PX(t) and PY(t) will be explained in Section 2.3.4.These results manifest the notable effects of laser polarization on the nonadiabatic transition between X and Y , which develops mainly after irradiation.

Angular Momentum Expectation Value and Vibrational Spectrum
Figure 5 shows lz(t) for the four types of excitations.Its behavior is analogous to that in Figure 3 while PX(t)  ~ PY(t) (t < 4 fs) except for its smaller amplitude due to the less efficiency of excitation and the coordinate dependence of λ(Q).The amplitude of lz(t), however, diminishes gradually, which is a characteristic feature absent in Figure 3.This is ascribed to two factors: decrease of the overlap In each panel, the populations of |X and |Y are denoted by the solid and dotted lines, respectively.The laser pulses fully decay at t = 7.26 fs.Reprinted from [56] with permission of the American Chemical Society.
The subsequent behaviors of P X (t) and P Y (t) are totally different among the four types of excitations.For e (+) (Figure 4a), a small portion of the population shifts from |X to |Y by t ~10 fs and then a downward population transfer takes place around t ~10-14 fs by nonadiabatic transition.In the case of e (−) (Figure 4b), a substantial amount of the population is transferred from |Y to |X and consequently P X (t) is more than seven times as large as P Y (t) at t ~10 fs.Afterwards, the direction of population transfer flips periodically with the rather small fraction of the population transferred.In contrast, for e +1 (Figure 4c), a part of the population is continuously exchanged between |X and |Y and, thus, electronic relaxation is completed at t > 30 fs.The behaviors of P X (t) and P Y (t) for e −1 (Figure 4d) are intermediate between those for e (+) and e (−) : A considerable population transfer from |Y to |X is observed around t ~8-13 fs.The very different evolutions of P X (t) and P Y (t) will be explained in Section 2.3.4.These results manifest the notable effects of laser polarization on the nonadiabatic transition between |X and |Y , which develops mainly after irradiation.

Angular Momentum Expectation Value and Vibrational Spectrum
Figure 5 shows l z (t) for the four types of excitations.Its behavior is analogous to that in Figure 3 while P X (t) P Y (t) (t < 4 fs) except for its smaller amplitude due to the less efficiency of excitation and the coordinate dependence of λ(Q).The amplitude of l z (t), however, diminishes gradually, which is a characteristic feature absent in Figure 3.This is ascribed to two factors: decrease of the overlap between the WPs moving on the two excited-state PESs, which occurs even without nonadiabatic couplings, and population transfer due to nonadiabatic couplings shown in Figure 4.Both of the two factors cause the loss of a superposition of |X and |Y ; in particular, the latter contributes to the loss dominantly.Indeed, we performed nuclear WP simulations within the Born-Oppenheimer approximation (BOA) [83], where the WPs simply propagated on the individual PESs, and found that the amplitude of l z (t) exhibits periodic reduction and recovery in turn owing to the temporal change in the WP overlap but never decays monotonically [56][57][58].The curves of l z (t) for e (+) and e +1 can be approximated by a sinusoidal exponential decay.The lifetime of the decay is ~6 fs for e (+) and ~18 fs for e +1 ; this difference stems from the different rates of nonadiabatic transition in Figure 4a,c.Meanwhile, the amplitudes of l z (t) for e (−) and e −1 do not undergo a monotonic decrease but make a small transient recovery [in the time ranges of t ~14-20 fs for e (−) and t ~18-24 fs for e −1 ].This arises from the regeneration of the superposition of |X and |Y due to the upward population transfer in the respective time ranges in Figure 4b,d.We also noticed that l z (t) oscillates with slightly different periods among the four types of excitations.The oscillation period of l z (t) is determined by the energy gap between the two PESs in the region where the WPs run; the region depends on the applied pulse as will be shown in Section 2.3.4.The results in Figure 5 demonstrate that nonadiabatic couplings weaken laser-driven ring currents within tens of femtoseconds.
factors cause the loss of a superposition of X and Y ; in particular, the latter contributes to the loss dominantly.Indeed, we performed nuclear WP simulations within the Born-Oppenheimer approximation (BOA) [83], where the WPs simply propagated on the individual PESs, and found that the amplitude of lz(t) exhibits periodic reduction and recovery in turn owing to the temporal change in the WP overlap but never decays monotonically [56][57][58].The curves of lz(t) for e (+) and e+1 can be approximated by a sinusoidal exponential decay.The lifetime of the decay is ~ 6 fs for e (+) and ~ 18 fs for e+1; this difference stems from the different rates of nonadiabatic transition in Figure 4a,c.Meanwhile, the amplitudes of lz(t) for e (−) and e−1 do not undergo a monotonic decrease but make a small transient recovery [in the time ranges of t ~ 14-20 fs for e (−) and t ~ 18-24 fs for e−1].This arises from the regeneration of the superposition of X and Y due to the upward population transfer in the respective time ranges in Figure 4b,d.We also noticed that lz(t) oscillates with slightly different periods among the four types of excitations.The oscillation period of lz(t) is determined by the energy gap between the two PESs in the region where the WPs run; the region depends on the applied pulse as will be shown in Section 2.3.4.The results in Figure 5 demonstrate that nonadiabatic couplings weaken laser-driven ring currents within tens of femtoseconds.
Q being the operator of Q.As seen in Figure 6a, the behaviors of Q(t) for linear polarizations are strongly dependent on the polarization direction: The vibrational amplitude for e (−) is more than twotimes as large as that for e (+) .In contrast, in Figure 6b, the trajectories of Q(t) for circular polarizations (e+1 and e−1) differ only a little and are located between those for linear ones.The behavior of Q(t) being the operator of Q.As seen in Figure 6a, the behaviors of Q(t) for linear polarizations are strongly dependent on the polarization direction: The vibrational amplitude for e (−) is more than two-times as large as that for e (+) .In contrast, in Figure 6b, the trajectories of Q(t) for circular polarizations (e +1 and e −1 ) differ only a little and are located between those for linear ones.The behavior of Q(t) hardly depends on laser polarization under the BOA [56][57][58], which implies the importance of nonadiabatic couplings.The frequency spectrum of ψ X (Q, t) after the nonadiabatic transition from |Y to |X can be obtained by the Fourier transform of its autocorrelation function [84]: The parameter τ was introduced to smooth the spectra and set to 39.6 fs, which is longer than the vibrational periods of the breathing and distortion modes (28.8 and 21.2 fs).The values of t i for e (+) , e (−) , e +1 , and e −1 were 14.0, 10.0, 34.0, and 13.0 fs, respectively, and u = 99.1 fs for all of them.The zero of ω was set to the minimum of the lower PES.
X Y to X can be obtained by the Fourier transform of its autocorrelation function [84]: The parameter τ was introduced to smooth the spectra and set to 39.6 fs, which is longer than the vibrational periods of the breathing and distortion modes (28.8 and 21.2 fs).The values of ti for e (+) , e (−) , e+1, and e−1 were 14.0, 10.0, 34.0, and 13.0 fs, respectively, and u = 99.1 fs for all of them.The zero of ω was set to the minimum of the lower PES.Reprinted from [56] with permission of the American Chemical Society.
The frequency spectra for the four types of excitations are depicted in Figure 7.For e (+) , the maximum of  

X
 is identified at  ~ 1400 cm −1 and another peak emerges at  ~ 2500 cm −1 in Figure 7a; for e (−) , the strongest peak of   X  is located at  ~ 2500 cm −1 and a couple of strong peaks are also found at  > 3000 cm −1 .The wave numbers of 1400, 2500, and 3000 cm −1 are almost equal to those of the lowest three vibrational states in G owing to the analogy between G and X in the PES around its minimum.The frequency spectra for linear polarizations in Figure 7a corroborate that at t > ti,    Reprinted from [56] with permission of the American Chemical Society.
The frequency spectra for the four types of excitations are depicted in Figure 7.For e (+) , the maximum of σ X (ω) is identified at ν ~1400 cm −1 and another peak emerges at ν ~2500 cm −1 in Figure 7a; for e (−) , the strongest peak of σ X (ω) is located at ν ~2500 cm −1 and a couple of strong peaks are also found at ν > 3000 cm −1 .The wave numbers of 1400, 2500, and 3000 cm −1 are almost equal to those of the lowest three vibrational states in |G owing to the analogy between |G and |X in the PES around its minimum.The frequency spectra for linear polarizations in Figure 7a corroborate that at t > t i , ψ X (Q, t) mainly consists of low (high) vibrational states for e (+) [e (−) ].The spectral features for circular ones in Figure 7b are very similar to each other: The primary peaks of σ X (ω) for e +1 and e −1 are both at ν ~1400 cm −1 , while the intensities of the other peaks are a little stronger in the latter case.This indicates that ψ X (Q, t) is composed of the same frequency components for e +1 and e −1 in the time domain after the nonadiabatic transition.  +) and e (−) are denoted by the solid and dotted lines, respectively; in panel (b), those for e+1 and e−1 are denoted by the solid and dotted lines, respectively.In each case, the values of σX(ω) were scaled so that the maximum value is unity.Reprinted from [56] with permission of the American Chemical Society.
The remarkable difference in the vibrational amplitude between e (+) and e (−) suggests that ring currents in aromatic molecules such as DCP irradiated by linearly polarized laser pulses can be distinguished experimentally by vibrational spectroscopy in the frequency domain without timeresolved measurements.This suggestion may be useful for rapid identification of molecular chirality.For instance, if two enantiomers of a polar chiral molecule are oriented so that their permanent dipoles are parallel to each other, the linear polarization vector e (+) or e (−) designed for one of the enantiomers does not necessarily coincide with that for the other enantiomer.Therefore, the strength of the induced ring current should be different between the enantiomers and the current direction can be opposite.This difference will be reflected in vibrational motions of the enantiomers.  +) and e (−) are denoted by the solid and dotted lines, respectively; in panel (b), those for e +1 and e −1 are denoted by the solid and dotted lines, respectively.In each case, the values of σ X (ω) were scaled so that the maximum value is unity.Reprinted from [56] with permission of the American Chemical Society.

Laser Control of Wave Packet Interference
The remarkable difference in the vibrational amplitude between e (+) and e (−) suggests that ring currents in aromatic molecules such as DCP irradiated by linearly polarized laser pulses can be distinguished experimentally by vibrational spectroscopy in the frequency domain without time-resolved measurements.This suggestion may be useful for rapid identification of molecular chirality.For instance, if two enantiomers of a polar chiral molecule are oriented so that their permanent dipoles are parallel to each other, the linear polarization vector e (+) or e (−) designed for one of the enantiomers does not necessarily coincide with that for the other enantiomer.Therefore, the strength of the induced ring current should be different between the enantiomers and the current direction can be opposite.This difference will be reflected in vibrational motions of the enantiomers.

Laser Control of Wave Packet Interference
Figure 8 shows the propagations of the adiabatic WPs on the two excited-state PESs.The probability densities |ψ X (Q, t)| 2 and |ψ Y (Q, t)| 2 at t ~5 fs resemble that of the initial WP in |G and then the WPs start to move along the gradient of each PES.As they approach the avoided crossing, the nonadiabatic effects on vibronic dynamics appear differently for the four types of excitations.For e (+) (Figure 8a), the WP in |Y is diminished by nonadiabatic transition at t ~12 fs and the contour map of |ψ X (Q, t)| 2 clearly displays the node arising from the interference; thereafter, the WPs on the two PESs are distorted significantly.In contrast, for e (−) (Figure 8b), the WPs maintain a Gaussian-like form even after the nonadiabatic transition already in progress at t ~8 fs.For circular polarizations, the WPs behave as expected from the behavior of the populations in Figure 4c,d.In the case of e +1 (Figure 8c), only a small fraction of ψ Y (Q, t) has been transferred to |X until t ~12 fs and no interference patterns are observed in the WPs.For e −1 (Figure 8d), the WPs exhibit intermediate features between those for e (+) and e (−) : No clear interference is observed in the WPs but the shape of |ψ X (Q, t)| 2 is deformed from a Gaussian especially at t ~12 fs.As a result of nonadiabatic couplings, the center of ψ X (Q, t), which contributes dominantly to Q(t) in Figure 6, runs in a low-energy (high-energy) region on the PES for e (+) [e (−) ], whereas the motions of the center for circular polarizations are similar to each other. Q , which contributes dominantly to Q(t) in Figure 6, runs in a low-energy (high- energy) region on the PES for e (+) [e (−) ], whereas the motions of the center for circular polarizations are similar to each other.The polarization effects on the expectation values lz(t) and Q(t) as well as on the populations PX(t) and PY(t) can be explained in terms of WP interference.As defined in Section 2.2, a linearly polarized laser pulse with e (−) initially produces a superposition of   phase (θ = π).Their relative quantum phase evolves during irradiation [by −ϑ(t) in the rectangularenvelope cases] and also gains a dynamical phase associated with the shape of each PES as the WPs propagate.In nonadiabatic transition, an additional phase shift is further imposed on the WP created by nonadiabatic couplings, which interferes with that on the other PES.We do not quantify the additional phase in DCP, but the downward population transfer around t ~ 5-10 fs in Figure 4b implies opposite interferences on the two PESs: The WPs are almost in phase and interfere constructively on the lower PES, while those on the higher PES are out of phase with destructive interference.The constructive interference on the lower PES works particularly on high vibrational states in   , X t  Q .The direction of the population transfer flips as the relative quantum phase grows.
For e (+) , the WPs prepared on the two PESs are in phase (θ = 0).The effects of WP interference are therefore reversed from those for e (−) : For DCP, the interference is destructive on the lower PES but constructive on the higher one around t ~ 5-10 fs.The amount of the WP created by the  The origin of the PESs is the optimized geometry of |G and the avoided crossing is designated by a circle.The bold contours represent the probability densities |ψ X (Q, t)| 2 and |ψ Y (Q, t)| 2 .The arrows indicate the motion of the center of the WPs.Modified from [56] with permission of the American Chemical Society.
The polarization effects on the expectation values l z (t) and Q(t) as well as on the populations P X (t) and P Y (t) can be explained in terms of WP interference.As defined in Section 2.2, a linearly polarized laser pulse with e (−) initially produces a superposition of ψ X (Q, t) and ψ Y (Q, t) out of phase (θ = π).Their relative quantum phase evolves during irradiation [by −ϑ(t) in the rectangular-envelope cases] and also gains a dynamical phase associated with the shape of each PES as the WPs propagate.In nonadiabatic transition, an additional phase shift is further imposed on the WP created by nonadiabatic couplings, which interferes with that on the other PES.We do not quantify the additional phase in DCP, but the downward population transfer around t ~5-10 fs in Figure 4b implies opposite interferences on the two PESs: The WPs are almost in phase and interfere constructively on the lower PES, while those on the higher PES are out of phase with destructive interference.The constructive interference on the lower PES works particularly on high vibrational states in ψ X (Q, t).The direction of the population transfer flips as the relative quantum phase grows.
For e (+) , the WPs prepared on the two PESs are in phase (θ = 0).The effects of WP interference are therefore reversed from those for e (−) : For DCP, the interference is destructive on the lower PES but constructive on the higher one around t ~5-10 fs.The amount of the WP created by the nonadiabatic transition from |X to |Y , which interferes with ψ Y (Q, t), is smaller than that for the transition from |Y to |X , which interferes with ψ X (Q, t); this is the reason why the resultant upward population transfer is small in Figure 4a.The WPs on the two PESs then reach the avoided crossing and the reverse population transfer occurs around t ~10-14 fs.The interference enhances low vibrational states in ψ X (Q, t).
In the cases of circular polarizations, ψ X (Q, t) and ψ Y (Q, t) are initially neither in phase nor out of phase because θ = ± (η Y − η X ) for e ±1 ; for DCP, θ = ± 0.35π.Hence, no fully constructive or destructive interference takes place on either PES at t ~5-10 fs, resulting in the vibrational amplitudes in Figure 6b, which are intermediate between those for e (+) and e (−) in Figure 6a.As the relative quantum phase of DCP evolves in the negative direction, for e −1 (θ = − 0.35π), it takes more time than for e (−) (θ = π) but less time than for e (+) (θ = 0) to experience constructive and destructive interferences on the lower and higher PESs, respectively.As a consequence, the downward population transfer for e −1 appears at t ~8-13 fs in Figure 4d.For e +1 (θ = 0.35π), the WPs pass through the avoided crossing before the requirement of interference for population transfer is met; the nonadiabatic transition is therefore not completed until they come closer again to the avoided crossing at t ~30 fs as in Figure 4c.
The present results indicate that the initial relative phase θ governs not only ring currents but also the subsequent molecular vibration through nonadiabatic couplings.Furthermore, it is possible to control the time-dependent phase due to irradiation, i.e., ϑ(t) by adjusting the other laser parameters such as the peak field strength E 0 and pulse duration t f .This suggests that a superposition of the quasi-degenerate states with an arbitrary relative phase can be created by applying a laser pulse of arbitrary polarization.The interference between nuclear WPs in nonadiabatic transition can be controlled as desired by attosecond/several-femtosecond laser pulses, leading to sophisticated control of molecular vibrations.

Coherent π-Electron Angular Momentum and Current of a Molecule Having Two Aromatic Rings
So far, we have considered the coherent quantum dynamics of π electrons in aromatic molecules with a single aromatic ring [55][56][57][58][60][61][62]85].In this section, we present the results of a theoretical study on coherent quantum dynamics of π electrons in a molecule with two aromatic rings.For such an aromatic molecule, ultrafast coherent π-electron dynamics that cannot be induced in a molecule with a single aromatic ring appears.Typical examples of the coherent π-electron dynamics are localization of angular momentum (ring current) at a specific aromatic ring and transfer of localized angular momentum (ring current) from one of the aromatic rings to the other one.Localization of π-electron rotation driven by pulsed lasers is expected to play an essential role in the development of organic optoelectronic switching devices for next generation as well as femtosecond and attosecond chemistry [86].

Two Ring Components of π-Electron Angular Momentum and Current
In this subsection, first, we derive an expression for the coherent π-electron angular momentum of aromatic ring molecules and that for the ring current by using the density matrix method [87].Next, we apply the theoretical expressions to (P)-2,2 -biphenol as a real chiral molecule with two benzene rings.This is a typical nonplanar aromatic ring molecule with axial chirality.

Decomposition of Angular Momentum and Current into Atomic Orbital Components
The expectation values of the coherent electronic angular momentum and current operators of an aromatic ring molecule excited by visible or UV light are expressed as [87] Here, Ô → r 1 is a one-electron operator for the angular momentum or current, Ψ(t) is an electronic wave function at time t, n is the number of electrons, and → r i denotes the electron coordinates of the i-th electron.
Ψ(t) is expanded in terms of electronic configurations, the ground electronic Φ 0 and singly excited electronic ones Φ α , as Here, Φ 0 is given as . φ a and φ b are among the n occupied orbitals, and Φ α is given as (LCAO-MO) approximation, a π orbital φ k is expanded as φ k = ∑ i c ki χ i (k = a, a , b, b ), where χ i is i-th AO and c ki is the MO coefficient.Equation ( 35) can be expressed within the MO approximations as follows: where the electronic-state density matrix element ρ αβ (t) is defined as ρ αβ (t) ≡ C α (t)C * β (t).Let us now express the expectation value of the coherent angular momentum operator and that of ring current operator of a space-fixed aromatic ring molecule with two aromatic rings connected with a single bond.

Coherent Angular Momentum
The total electron angular momentum operator of an aromatic ring molecule can be expressed as a sum of the angular momentum operator of each aromatic ring, which is defined as is the electronic angular momentum operator of the z-component of ring K (= L or R).Here, L and R denote the ring on the left-hand side and that on the right-hand side, respectively.Coordinates x K and y K operate on AOs belonging to ring K, and → n K is the unit vector perpendicular to the ring.The expectation value of the angular momentum operator is given, for example, in terms of a 2p z carbon AOs as [87]: where Here, a 2 is a constant related to the orbital exponent of the 2p z AO and r ij is the distance between i and j carbon atom sites.

Coherent Ring Current
The ring current passing through a surface S at time t is defined as follows: Here, → n ⊥ is a unit vector perpendicular to surface S, and Ĵ → r , t ≡ Ψ(t) Ĵ → r Ψ(t) , which can be expressed as: where Ĵ → r = (eh/2m e i)  39) is carried out over a planar plane S perpendicular to a C-C bond of an aromatic ring.
Since the L and R rings are not round but consist of nonequivalent C-C bonds, we introduce the bond current J(t) ij from the nearest neighbor atoms at sites j to i, which is given in terms of inter-atomic bond current (IABC) J I ABC ij as where Here, surface S is set to be perpendicular to the C i and C j bond at the center.Equation ( 42) is given in the 2p carbon AO basis set, {χ i }, as [87]: Here, Θ = π for the bond current of the chemical bond belonging to one of the two aromatic rings, L or R, while Θ = Θ d for the bridge bond current of bond C i − C j with a dihedral angle between the two rings Θ d .C i (C j ) refers to the bridge carbon belonging to the R (L) ring.We now define an effective ring current along ring K, J(t) K , by taking the average over all of the bond currents as [87]: where N K is the number of bonds of ring K. Time-dependent behaviors of electronic dynamics induced by a pulsed laser field, E(t), can be directly obtained by solving the coupled equations of motion of ρ αβ (t) with the initial conditions ρ 00 (0) = 1 and ρ αα (0) = 0 for α = 0; ρ αβ (0) = 0 for α = β (hereafter Greek subscripts refer to electronic states): Here, V αγ (t) denotes the coupling between electronic states α and γ through the molecule-field interaction, V(t) = − μ • E(t), ω βα is the angular frequency difference between two electronic states α and β, and αβ ] is the dephasing constant in Markov approximation [88].Here, γ α (γ β ) is the nonradiative transition rate constant of state α (β), and γ αβ is the pure dephasing constant that is induced by elastic interactions between the molecule of interest and heat baths including intramolecular vibrations.

Application to (P)-2,2 -Biphenol
For generation of coherent angular momentum and ring current of (P)-2,2 -biphenol, we focus on the three optically allowed excited states (a, b 1 , and b 2 ) whose energies range from 6.67 to 6.84 eV as shown in Figure 9a [87].The geometry of (P)-2,2 -biphenol in the ground state was optimized by using the density functional theory (DFT) at the B3LYP level of theory in the GAUSSIAN09 program [89].(P)-2,2 -biphenol was assumed to be fixed on a surface by a non-conjugated chemical bond or in a space by laser molecular orientation techniques.The laboratory-fixed Y-axis is set to be parallel to the single bond bridging the two phenol groups, and the rotation axis of point group C 2 is placed along the laboratory-fixed Z-axis parallel to the surface normal.The dihedral angle between the two phenol groups, Θ d , was found to be 108.9degrees from our DFT calculations [86].
fixed Y-axis is set to be parallel to the single bond bridging the two phenol groups, and the rotation axis of point group C2 is placed along the laboratory-fixed Z-axis parallel to the surface normal.The dihedral angle between the two phenol groups, Θd, was found to be 108.9degrees from our DFT calculations [86].The transition energies to the a (A), b1 (B), and b2 (B) states, which were calculated at the optimized ground state geometry using the time-dependent DFT at the B3LYP level of theory, are 6.67, 6.78, and 6.84 eV, respectively [86]. Figure 9b shows the directions of transition dipole moments.The three optically allowed excited states associated with the present laser-induced electron dynamics have one-electronic excitation configurations since they basically consist of optically allowed ππ* excited states of the two phenol rings.The transition dipole moment from the ground state to excited state a (A), μga [= (0, 0, −0.77) in atomic units], is parallel to the Z-axis, and the others, μgb1 [= (0.08, 1.93, 0)] and μgb2 [= (1.24, −0.34, 0)], are nearly orthogonal to each other in the XY-plane as shown in Figure 9.

Generation of Two-State Electronic Coherence by Linearly Polarized UV Pulses
We can generate three types of two-state electronic coherence in the present three excited state model, which can be represented as (b1 ± b2), (a ± b1), and (a ± b2).The phase of each electronic coherence can be set by applying a linearly polarized UV laser pulse with a properly selected polarization direction.Because of the nonplanar geometrical structure, the angular momentum is twodimensional and the ring current is on two different surfaces: Z-directional angular momentum is generated for the created coherent state with A symmetry and corresponding ring current is generated on the Z-directional angular momentum, while X-directional angular momentum is generated for the coherent state with B symmetry and ring current is generated on the X-directional angular momentum.
Figure 10 shows time-dependent angular momenta calculated for three types of electronic coherence, (b1 + b2), (a + b1), and (a + b2), each of which is created by a linearly polarized UV pulse with a properly selected polarization direction.The envelope function of the pulse for each coherence is shown in the right-hand side.The pulse amplitudes for (b1 + b2), (a + b1), and (a + b2) coherences are 0.19, 0.83 and 3.32 TW/cm 2 , respectively.Angular momentum along the Z-axis, lZ, is generated for (b1 + b2) electronic coherence, while those along the X-axis, lX, are generated with π-phase shift for (a + b1) electronic coherence as well as (a + b2) electronic coherence.Similarly, for out-of-phase electronic coherences (b1 − b2), (a − b1), and (a − b2), each angular momentum can be expressed by that for the corresponding in-phase coherence with π-phase shift.The time-dependent angular momenta in the The transition energies to the a (A), b 1 (B), and b 2 (B) states, which were calculated at the optimized ground state geometry using the time-dependent DFT at the B3LYP level of theory, are 6.67, 6.78, and 6.84 eV, respectively [86]. Figure 9b shows the directions of transition dipole moments.The three optically allowed excited states associated with the present laser-induced electron dynamics have one-electronic excitation configurations since they basically consist of optically allowed ππ* excited states of the two phenol rings.The transition dipole moment from the ground state to excited state a (A), µ ga [= (0, 0, −0.77) in atomic units], is parallel to the Z-axis, and the others, µ gb1 [= (0.08, 1.93, 0)] and µ gb2 [= (1.24, −0.34, 0)], are nearly orthogonal to each other in the XY-plane as shown in Figure 9.

Generation of Two-State Electronic Coherence by Linearly Polarized UV Pulses
We can generate three types of two-state electronic coherence in the present three excited state model, which can be represented as (b 1 ± b 2 ), (a ± b 1 ), and (a ± b 2 ).The phase of each electronic coherence can be set by applying a linearly polarized UV laser pulse with a properly selected polarization direction.Because of the nonplanar geometrical structure, the angular momentum is two-dimensional and the ring current is on two different surfaces: Z-directional angular momentum is generated for the created coherent state with A symmetry and corresponding ring current is generated on the Z-directional angular momentum, while X-directional angular momentum is generated for the coherent state with B symmetry and ring current is generated on the X-directional angular momentum.
Figure 10 shows time-dependent angular momenta calculated for three types of electronic coherence, (b 1 + b 2 ), (a + b 1 ), and (a + b 2 ), each of which is created by a linearly polarized UV pulse with a properly selected polarization direction.The envelope function of the pulse for each coherence is shown in the right-hand side.The pulse amplitudes for (b 1 + b 2 ), (a + b 1 ), and (a + b 2 ) coherences are 0.19, 0.83 and 3.32 TW/cm 2 , respectively.Angular momentum along the Z-axis, l Z , is generated for (b 1 + b 2 ) electronic coherence, while those along the X-axis, l X , are generated with π-phase shift for (a + b 1 ) electronic coherence as well as (a + b 2 ) electronic coherence.Similarly, for out-of-phase electronic coherences (b 1 − b 2 ), (a − b 1 ), and (a − b 2 ), each angular momentum can be expressed by that for the corresponding in-phase coherence with π-phase shift.The time-dependent angular momenta in the Z and X directions are, respectively, expressed in terms of l L (t) as l Z (t) = 2l L (t) sin(Θ d /2) and l X (t) = 2l L (t) cos(Θ d /2).l Z (t) was calculated using the relation l L (t) = l R (t) and l X (t) was calculated using l L (t) = −l R (t) as shown in Figure 10.The simple sinusoidal time-dependence with a damping originates from creation of two-state electronic coherence, the oscillation period of which corresponds to the frequency difference between the two states.The damping behaviors of the oscillation amplitudes are due to dephasing between the two excited states.Here, dephasing constants γ ab1 = γ ab2 = γ b1b2 = 0.01 eV were adopted [87].The oscillations in the angular momenta should be called angular momentum quantum beats, similar to quantum beats originating from vibronic coherence in fluorescence [90].It should be noted that in angular momentum quantum beats, π electrons rotate in a unidirectional way a few times within the half cycle of the oscillation.The unidirectional π-electron rotation can produce unidirectional ring current and current-induced magnetic flux.This is the basic principle to design ultrafast switching devices made from organic ring molecules.
Z and X directions are, respectively, expressed in terms of

 
Z lt was calculated using the relation lL(t) = lR(t) and   X lt was calculated using lL(t) = −lR(t) as shown in Figure 10.The simple sinusoidal time-dependence with a damping originates from creation of two-state electronic coherence, the oscillation period of which corresponds to the frequency difference between the two states.The damping behaviors of the oscillation amplitudes are due to dephasing between the two excited states.Here, dephasing constants γab1 = γab2 = γb1b2 = 0.01 eV were adopted [87].The oscillations in the angular momenta should be called angular momentum quantum beats, similar to quantum beats originating from vibronic coherence in fluorescence [90].It should be noted that in angular momentum quantum beats, π electrons rotate in a unidirectional way a few times within the half cycle of the oscillation.The unidirectional π-electron rotation can produce unidirectional ring current and current-induced magnetic flux.This is the basic principle to design ultrafast switching devices made from organic ring molecules.
Figure 10.Time-dependent coherent angular momentum of (P)-2,2′-biphenol for three types of electronic coherence.Modified from [87] with permission of the American Institute of Physics.The pulse form for (a + b2) electronic coherence in [87] should be replaced by that in the figure.[ℏ] Figure 10.Time-dependent coherent angular momentum of (P)-2,2 -biphenol for three types of electronic coherence.Modified from [87] with permission of the American Institute of Physics.The pulse form for (a + b 2 ) electronic coherence in [87] should be replaced by that in the figure.

Bond Currents and Time Evolution of Coherent Ring Current
Figure 11 shows the magnitudes of bond current, J ij , which were calculated using Equation (42) at the maximum coherence time together with the directions of bond currents for the three types of electronic coherence in the first period of oscillation.The arrow along each bond of (P)-2,2 -biphenol denotes the initial direction of the bond current.Magnitudes of the averaged ring current at the maximum coherence time are of an order of tens of µA: J = 161, 86.5 and 63.4 µA for (b 1 + b 2 ), (a + b 1 ), and (a + b 2 ) coherences, respectively [87].If the phenol ring is assumed to be a ring with radius r and the magnitude of the angular momentum l = rp with linear momentum of electrons, p, is given, the ring current J can be estimated by using the relation 2πrJ = ep/m e .The magnitudes of the ring currents for the (b 1 + b 2 ), (a + b 1 ), and (a + b 2 ) coherences are 216, 163.5 and 25.5 µA, respectively.The magnitude of each ring current is large compared with that of the corresponding averaged ring current except for (a + b 2 ) coherence.The difference in the magnitudes between the averaged bond current at the maximum coherence time and the ring current calculated within a perfect ring model mainly originates from the fact that phenol ring has a functional OH group.Barth et al. reported an electronic current of 84.5 µA for Mg porphyrin excited by a circularly polarized UV intense laser pulse to yield population inversion [51].It is interesting that (P)-2,2 -biphenol and Mg porphyrin have the same order of magnitude in ring currents, though they have different numbers of electrons and different radii of their rotation.
Figure 11 shows the magnitudes of bond current, Jij, which were calculated using Equation (42) at the maximum coherence time together with the directions of bond currents for the three types of electronic coherence in the first period of oscillation.The arrow along each bond of (P)-2,2′-biphenol denotes the initial direction of the bond current.Magnitudes of the averaged ring current at the maximum coherence time are of an order of tens of μA: J = 161, 86.5 and 63.4 μA for (b1 + b2), (a + b1), and (a + b2) coherences, respectively [87].If the phenol ring is assumed to be a ring with radius r and the magnitude of the angular momentum l = rp with linear momentum of electrons, p, is given, the ring current J can be estimated by using the relation 2πrJ = ep/me.The magnitudes of the ring currents for the (b1 + b2), (a + b1), and (a + b2) coherences are 216, 163.5 and 25.5 μA, respectively.The magnitude of each ring current is large compared with that of the corresponding averaged ring current except for (a + b2) coherence.The difference in the magnitudes between the averaged bond current at the maximum coherence time and the ring current calculated within a perfect ring model mainly originates from the fact that phenol ring has a functional OH group.Barth et al. reported an electronic current of 84.5 μA for Mg porphyrin excited by a circularly polarized UV intense laser pulse to yield population inversion [51].It is interesting that (P)-2,2′-biphenol and Mg porphyrin have the same order of magnitude in ring currents, though they have different numbers of electrons and different radii of their rotation.Bond currents for the three types of the electronic coherence at a time when the maximum electronic coherence is created.Modified from [87] with permission of the American Institute of Physics.

Quantum Switching of Coherent π-Electron Rotations in a Nonplanar Chiral Aromatic Molecule
In general, π electrons in a planar aromatic molecule generate electronic angular momentum perpendicular to their aromatic ring.In a nonplanar aromatic ring molecule, on the other hand, it is possible for π electron rotations to generate multi-dimensional angular momenta, each component of which has a particular angular momentum vector with a definite direction and generates the corresponding angular momentum.Therefore, nonplanar aromatic ring molecules are promising as a sophisticated switching device with several controllable variables [86].

Control Scheme for Angular Momentum Switching
As shown in Figure 12a, four rotational patterns (CC, AA, CA, and AC), in which C (A) means clockwise (anticlockwise) rotation along the L or R phenol ring, are possible for (P)-2,2′-biphenol, depending on two factors: one is the polarization direction of the applied linearly polarized laser pulse and the other is the symmetry of the created coherent state.Let us first clarify how these patterns can be prepared, and then how ultrafast switching between these rotational patterns can be realized.As already shown by Equation (31), the polarization direction is expressed in terms of a unit vector,    Bond currents for the three types of the electronic coherence at a time when the maximum electronic coherence is created.Modified from [87] with permission of the American Institute of Physics.

Quantum Switching of Coherent π-Electron Rotations in a Nonplanar Chiral Aromatic Molecule
In general, π electrons in a planar aromatic molecule generate electronic angular momentum perpendicular to their aromatic ring.In a nonplanar aromatic ring molecule, on the other hand, it is possible for π electron rotations to generate multi-dimensional angular momenta, each component of which has a particular angular momentum vector with a definite direction and generates the corresponding angular momentum.Therefore, nonplanar aromatic ring molecules are promising as a sophisticated switching device with several controllable variables [86].

Control Scheme for Angular Momentum Switching
As shown in Figure 12a, four rotational patterns (CC, AA, CA, and AC), in which C (A) means clockwise (anticlockwise) rotation along the L or R phenol ring, are possible for (P)-2,2 -biphenol, depending on two factors: one is the polarization direction of the applied linearly polarized laser pulse and the other is the symmetry of the created coherent state.Let us first clarify how these patterns can be prepared, and then how ultrafast switching between these rotational patterns can be realized.As already shown by Equation (31), the polarization direction is expressed in terms of a unit vector, e   Modified from [86] with permission of the American Chemical Society.For example, CA rotation means clockwise and anticlockwise rotations along the phenol rings L and R, respectively.(b) Resultant angular momenta l Z and l X for the four rotational patterns.l L (l R ) denotes the electronic angular momentum of the L (R) ring.Modified from [86] with permission of the American Chemical Society.
Figure 12b shows the resultant angular momentum l = l L + l R .Here, the π-electron angular momentum l L (l R ) is perpendicular to the L (R) phenol ring, and the resultant total angular momentum l Z (l X ) is parallel to the Z-(X-) axis on the ZX-plane.The directions of l Z indicate that clockwise (anticlockwise) rotation around the Z-axis is created by e excitations induce CA and AC rotations, respectively.Let us now consider two-dimensional current switching on the basis of the results shown in Figure 12.Here, two-dimensional switching means a sequential variation of electronic angular momentum with its definite sign (plus or minus) along the Zor X-axis.The switching should be carried out before reverse rotation of π electrons begins since the prepared coherent states are not eigenstates.For this purpose, we take a sequential four-step control as an example, indicated as l Z (−) → l X (+) → l Z (+) → l X (−), which is generated by the switching of rotational patterns CC → AC → AA → CA.Here, l Z (−) [l X (+)] denotes the π-electron angular momentum along the Z-(X-) axis with a negative (positive) sign, i.e., clockwise (anticlockwise) rotation of π electrons around the corresponding axis.

Simulation Results of Angular Momentum Switching
Figure 13a shows a 3D plot of the resultant angular momentum switching based on the sequential four-step scheme, which demonstrates that π-electron rotations are successfully controlled by the pulses depicted in Figure 13b.Control at each switching step was carried out by using pump and dump pulses with specific polarization directions.The pulsed laser with amplitude of 1.2 GV/m (= 1.9 × 10 11 W/cm 2 ) [14] was used at the second and fourth steps [86].
irreducible representations.π Electrons start to rotate along the two phenol rings in opposite directions to each other (CA or AC rotation).Therefore, the total angular momentum parallel to the Z-axis is zero and the nonzero resultant angular momentum lX is generated along the X-axis on the ZX-plane.If the two phenol rings are in the same plane without twist, both lX and lZ vanish.
excitations induce CA and AC rotations, respectively.Let us now consider two-dimensional current switching on the basis of the results shown in Figure 12.Here, twodimensional switching means a sequential variation of electronic angular momentum with its definite sign (plus or minus) along the Z-or X-axis.The switching should be carried out before reverse rotation of π electrons begins since the prepared coherent states are not eigenstates.For this purpose, we take a sequential four-step control as an example, indicated as lZ (−) → lX (+) → lZ (+) → lX (−), which is generated by the switching of rotational patterns CC → AC → AA → CA.Here, lZ (−) [lX (+)] denotes the π-electron angular momentum along the Z-(X-) axis with a negative (positive) sign, i.e., clockwise (anticlockwise) rotation of π electrons around the corresponding axis.

Simulation Results of Angular Momentum Switching
Figure 13a shows a 3D plot of the resultant angular momentum switching based on the sequential four-step scheme, which demonstrates that π-electron rotations are successfully controlled by the pulses depicted in Figure 13b.Control at each switching step was carried out by using pump and dump pulses with specific polarization directions.The pulsed laser with amplitude of 1.2 GV/m (= 1.9 × 10 11 W/cm 2 ) [14] was used at the second and fourth steps [86].αβ ] polarization.Each pulse has a bandwidth that is sufficient to coherently excite two quasi-degenerate electronic excited states.Secondly, the pump and dump pulses partially overlap.For the first step, i.e., creation of CC rotation, for example, the electric field of the pump pulse was given as e (+) b1b2 (t) = E (+) b1b2 E 0 b1b2 sin 2 (πt/T b1b2 ) sin(ω c,b1b2 t), while that of the dump pulse was E [86].Here, E 0 b1b2 is the amplitude of the pulse, T b1b2 (= 60.9 fs) is the oscillation period between the two excited states, b 1 and b 2 states, with the final condition ξ t Here, ξ(t) is time-dependent Lagrange multiplier, which satisfies a time-dependent equation: with the final condition . In order to carry out quantum control simulation for localization of the π-electron angular momentum of (P)-2,2′-biphenol, we need a set of angular momentum eigenstates defined in each aromatic ring.The localized angular momentum eigenstates ( A K W and C K W with K = L, R) can be expressed in terms of the three excited states (Φa, Φb1, Φb2) as [91]: The transformation matrices in the above equations were calculated by using the electronic angular momentum operator.The target states for the optimal control are, for example, expressed as follows: ˆAA  Here, O means no rotation.Reprinted from [91] with permission of the American Chemical Society.

Simulation Results of Angular Momentum Localization
In order to carry out quantum control simulation for localization of the π-electron angular momentum of (P)-2,2 -biphenol, we need a set of angular momentum eigenstates defined in each aromatic ring.The localized angular momentum eigenstates (W A K and W C K with K = L, R) can be expressed in terms of the three excited states (Φ a , Φ b1 , Φ b2 ) as [91]: The transformation matrices in the above equations were calculated by using the electronic angular momentum operator.The target states for the optimal control are, for example, expressed as follows: Here, O means no rotation on the L or R ring.

Simulation Results of Angular Momentum Localization
Figure 15 shows the quantum control results of the first step for the localization of angular momentum.The control time was set to t f = 30 fs, which is the inverse of the energy differences between two electronic excited states.In Figure 15a, the expectation value at t f = 30 fs, l L (t f ) [l R (t f )], reaches the maximum value for the [AO] ([OA]) target, while that of l R (t f ) [l L (t f )] is suppressed to zero.This demonstrates that the π-electron localization on one benzene ring is successfully achieved.Figure 15b represents the time-dependent electronic coherences, i.e., imaginary parts of ρ ab1 (t), ρ b1b2 (t), and ρ ab2 (t).The imaginary part of ρ b1b2 (t) for the [AO] and that for [OA] targets provide the same time dependence, while the imaginary parts of ρ ab1 (t) and ρ ab2 (t) for the [AO] target are shifted by π compared to those for the [OA] target.These phase relations originate from that of the target operator in Equation (51). Figure 15c shows Re E(ω), real parts of the optimal electric fields in the frequency domain, for the [AO] target (solid line) and the [OA] target (dashed line).The vertical dotted lines refer to the energies of the three eigenstates.The peaks of E X (ω), E Y (ω), and E Z (ω) appeared at the resonant positions of b 2 , b 1 , and a states, respectively.The phases of the z-component are shifted by π between [AO] and [OA] targets, while phase shifts for the other two components, E X (ω) and E Y (ω), are zero.This indicates that the π-phase shift of the electric field E Z (ω) is the key to control π electron rotations for the [AO] and [OA] localization.It should be noted that E Z (ω) excites only the a state.The difference in the target coefficients between W A L and W A R is that of Φ a .This is the reason why the π-phase shift appeared only in the z-component.The phase shift of E Z (ω) results in π-phase shift in coefficients of Φ a , Φ b1 , and Φ b2 , and finally the [AO] and [OA] localizations are realized at t = 30 fs.The maximum amplitude of the optimal field is about 1 GVm −1 (0.13 TWcm −2 ), which is much lower than the so-called threshold of multiphoton processes [91].[91] with permission of the American Chemical Society.

Simulation Results of Maintaining Angular Momentum Localization
Figure 16 shows the quantum control results of the second step for the localization, i.e., maintaining the unidirectional (anticlockwise or clockwise) rotation of π electrons around the L-or R-axis.After the localization has been completed, the localized angular momentum undergoes free propagation if there is no control field.Figure 16a shows the time dependence of the expectation values of the π-electron angular momentum under the field-free condition after the [AO] localization  [91] with permission of the American Chemical Society.

Simulation Results of Maintaining Angular Momentum Localization
Figure 16 shows the quantum control results of the second step for the localization, i.e., maintaining the unidirectional (anticlockwise or clockwise) rotation of π electrons around the Lor R-axis.After the localization has been completed, the localized angular momentum undergoes free propagation if there is no control field.Figure 16a shows the time dependence of the expectation values of the π-electron angular momentum under the field-free condition after the [AO] localization is achieved.Here, the initial time t = 0 was set to be the time when the first step of localization has been achieved.The expectation value of the π-electron angular momentum around the L-axis, l L (t), starts to oscillate, and that around the R-axis, l R (t), deviates from zero and oscillates in the opposite way. Figure 16b shows the results of the second step control to the initial localization demonstrated in Figure 15a.The control time was again set at t f = 30 fs. Figure 16b indicates that l L (t) keeps positive values under the optimal field.It is possible to keep the localization at any longer times t f > 30 fs.Although there exists small deviation, l R (t) is kept at around zero.The time-averaged expectation value of the target states defined as J L = t f 0 dt Ψ(t) Ô Ψ(t) /t f was calculated to be 0.85, while that of the free propagation was 0.22.Figure 16c shows the absolute values of the optimal electric fields, E X , E Y , and E Z in the frequency domain.The main peaks of all components appear at around 6.9 eV just above the highest eigenstates b 2 .The optimal field E(t) in the time domain is also shown in the inserted figure to indicate the phase relations.It can be seen that compared to E X (t), E Z (t) is shifted by π, and E Y (t) is shifted by π/2, and these relations are almost kept during the control time.This means that the resultant optimal field has an elliptically polarized character.We now briefly mention how to observe the localization of π-electron rotations around one of the benzene rings.This may be carried out by separately measuring the ring currents and magnetic fields generated on the two phenol groups since the two benzene rings are not planar but tilted with the dihedral angle of 108.8  was calculated to be 0.85, while that of the free propagation was 0.22.Figure 16c shows the absolute values of the optimal electric fields, EX, EY, and EZ in the frequency domain.The main peaks of all components appear at around 6.9 eV just above the highest eigenstates b2.The optimal field E(t) in the time domain is also shown in the inserted figure to indicate the phase relations.It can be seen that compared to EX(t), EZ(t) is shifted by π, and EY(t) is shifted by π/2, and these relations are almost kept during the control time.This means that the resultant optimal field has an elliptically polarized character.We now briefly mention how to observe the localization of π-electron rotations around one of the benzene rings.This may be carried out by separately measuring the ring currents and magnetic fields generated on the two phenol groups since the two benzene rings are not planar but tilted with the dihedral angle of 108.8°.The simulations in this section were performed using a simplified model omitting both the effects of vibrational motions and those of heat bath modes.In real aromatic ring molecules these effects influence on π-electron rotations.Mineo et al. recently found that intramolecular vibrational effects make significant contributions to dephasing of π-electron rotations under finite temperature conditions [95] as well as bath mode effects in condensed phases beyond the Markov approximation [96].It has been also shown in Section 2.3 that nonadiabatic transitions strongly influence π-electron rotations [97].Application of optimal control theory to such open systems is one of the important subsequent issues.[91] with permission of the American Chemical Society.

Summary and Outlook
In this article, we have focused on several fundamental issues of angular momentum (ring current) of π electrons in a low-symmetry aromatic ring molecule.One of the issues is on characteristics of the angular momentum generated by ultrashort polarized laser pulse(s).The The simulations in this section were performed using a simplified model omitting both the effects of vibrational motions and those of heat bath modes.In real aromatic ring molecules these effects influence on π-electron rotations.Mineo et al. recently found that intramolecular vibrational effects make significant contributions to dephasing of π-electron rotations under finite temperature conditions [95] as well as bath mode effects in condensed phases beyond the Markov approximation [96].It has been also shown in Section 2.3 that nonadiabatic transitions strongly influence π-electron rotations [97].Application of optimal control theory to such open systems is one of the important subsequent issues.

Summary and Outlook
In this article, we have focused on several fundamental issues of angular momentum (ring current) of π electrons in a low-symmetry aromatic ring molecule.One of the issues is on characteristics of the angular momentum generated by ultrashort polarized laser pulse(s).The generated angular momentum, called coherent angular momentum, is time-dependent, i.e., its direction changes in time, the duration of which is inversely proportional to the energy difference between quasi-degenerate, two electronic eigenstates.The phase in the oscillatory behavior of angular momentum, which is directly related to the relative quantum phase of the superposed quasi-degenerate states, can be controlled by the ellipticity and orientation of the incident laser.This is contrasted with the angular momentum generated in a highly symmetric aromatic molecule, which is time-independent, i.e., unidirectional.That is, stationary angular momentum is generated by an excitation of a degenerate excited state by a circularly polarized laser.
Another issue is the nonadiabatic coupling between laser-driven π-electron ring currents and molecular vibration.The coherent angular momentum is gradually attenuated by the nonadiabatic coupling.On the other hand, the amplitude of induced molecular vibration is prominently dependent on the direction of linear polarization vectors but insensitive to the helicity of circular polarization.The characteristic feature in vibrational amplitudes is ascribed to the WP interference effects in nonadiabatic transition.This offers a new possibility of attosecond/several-femtosecond polarized laser pulses as a promising tool to manipulate molecular vibration (and chemical reactions) through the WP interference in nonadiabatic transition.
Recently, there has been new progress on the issue of coherent angular momentum.Yamaki et al. have performed optimal control simulations of coherent π-electron rotations in (P)-2,2 -biphenol by taking into account two types of control targets: one is generation of the maximum angular momentum of the x component, which corresponds to CC or AA rotation of π electrons shown in Figure 12, and the other is the maintaining of the generated unidirectional angular momentum during a setting time duration [98].The optimal control pulse for each target was designed.The analysis of the simulation results shows that the effective maintaining of unidirectional angular momentum can be realized by applying the 2π pulse to one of the two excited states forming a coherent electronic state.The 2π pulse prevents the reverse rotation of π electrons by dumping the excited state population to the ground state and subsequently by pumping the population back to the excited state.In addition, Yamaki et al. have found an effectively stationary angular momentum in a low-symmetry molecule [99].This is realized by applying two linearly polarized lasers with the same frequency and with an appropriate relative phase to quasi-degenerate two excited states.This is due to creation of a dressed state with equal populations for the two excited eigenstates.Mineo et al. have proposed another laser-control scenario of generation of unidirectional angular momentum in a low-symmetry aromatic ring molecule [100].This is intuitively explained in terms of dynamic Stark shifts of two relevant excited states, which are induced by using two linearly polarized stationary lasers: Each laser is set to selectively interact with one of the two excited states.The lower and higher excited states are shifted up and down with the same rate, respectively, and the two excited states become degenerate at their midpoint.The control scenario with the ground state as the initial condition was applied to toluene.The derived time-dependent angular momentum consists of a train of unidirectional angular momentum.Further theoretical investigations on coherent π-electron rotation dynamics of aromatic ring molecules having no degenerate excited states are required by taking into account vibrational modes.
So far, we took simple low-symmetry molecules of a single aromatic ring or two aromatic rings for demonstration of fundamental features of coherent π-electron dynamics.From the viewpoint of practical application, organic materials having a large number of aromatic rings such as polycyclic aromatic hydrocarbons (PAHs) are much more interesting.Attention has already been paid to PAHs as organic molecular electronic materials, particularly for field-effect transistor devices [101].It is important to explore the possibility of such an optical response in PAHs.Mineo and Fujimura have theoretically demonstrated the control of π electrons in PAHs by lasers to create ring currents and current-induced magnetic fields [102,103].This is expected to serve as a fundamental research for next-generation organic optical switching devices.
The theoretical approach proposed in this article is useful to investigate not only the polarization dependence of intramolecular electron dynamics caused by single-photon absorption but also that of multiphoton electronic excitation.Hertel et al. have experimentally discovered that the multiphoton ionization probabilities of a xenon atom and C 60 fullerene irradiated by an intense femtosecond near-IR laser greatly changes with the ellipticity of the laser [104,105].The application of the present approach to the analysis of multiphoton excitation processes induced by polarized light has been reported elsewhere [106].

Figure 1 .
Figure 1.Orientation of the polarization ellipse.The laser field propagates in the z direction and oscillates in the xy plane.The orientation angle of the major axis of the polarization ellipse with respect to the x axis is represented by δ.The minor-to-major axial ratio is |tan (π/4 − β)|.Modified from[56] with permission of the American Chemical Society.

Figure 1 .
Figure 1.Orientation of the polarization ellipse.The laser field propagates in the z direction and oscillates in the xy plane.The orientation angle of the major axis of the polarization ellipse with respect to the x axis is represented by δ.The minor-to-major axial ratio is |tan (π/4 − β)|.Modified from[56] with permission of the American Chemical Society.
charge and a0 is the Bohr radius.The directions of µX and µY are illustrated in Figure2band the angle between them is ηY − ηX = 0.35π.The linear combinations of X and Y defined by Equation (2) are the approximate angular momentum eigenstates  and  with λ = 0.98ħ.Positive (negative) angular momentum designates π electrons traveling anticlockwise (clockwise) around the ring in Figure2a.

Figure 2 .
Figure 2. (a) Molecular formula of 2,5-dichloropyrazine (DCP).(b) Directions of transition moments µX and µY at the ground-state optimized geometry of DCP as well as those of linear polarization vectors e (+) and e (−) defined by Equation (31).The angle between e (+) and e (−) is 0.48π.Normal mode vectors of the (c) breathing and (d) distortion modes are represented by thick arrows.Modified from [56] with permission of the American Chemical Society.

Figure 2 .
Figure 2. (a) Molecular formula of 2,5-dichloropyrazine (DCP).(b) Directions of transition moments µ X and µ Y at the ground-state optimized geometry of DCP as well as those of linear polarization vectors e (+) and e (−) defined by Equation (31).The angle between e (+) and e (−) is 0.48π.Normal mode vectors of the (c) breathing and (d) distortion modes are represented by thick arrows.Modified from [56] with permission of the American Chemical Society.

Figure 3 .
Figure 3. Expectation value of the electronic angular momentum lz(t) in DCP irradiated by (a) linearly and (b) circularly polarized vacuum UV (VUV) laser pulses calculated from electronic WP simulations at the optimized geometry of G .In panel (a), the expectation values for e (+) and e (−) are denoted by the solid and dotted lines, respectively.The peaks at t = 5.4 fs are indicated by the vertical broken line.In panel (b), the expectation values for e+1 and e−1 are denoted by the solid and dotted lines, respectively.The peaks at t = 4.5 and 7.0 fs are indicated by the vertical broken lines.The laser pulses fully decay at t = 7.26 fs.Modified from [58] with permission of the Physical Society of the Republic of China.

Figure 3 .
Figure 3. Expectation value of the electronic angular momentum l z (t) in DCP irradiated by (a) linearly and (b) circularly polarized vacuum UV (VUV) laser pulses calculated from electronic WP simulations at the optimized geometry of |G .In panel (a), the expectation values for e (+) and e (−) are denoted by the solid and dotted lines, respectively.The peaks at t = 5.4 fs are indicated by the vertical broken line.In panel (b), the expectation values for e +1 and e −1 are denoted by the solid and dotted lines, respectively.The peaks at t = 4.5 and 7.0 fs are indicated by the vertical broken lines.The laser pulses fully decay at t = 7.26 fs.Modified from [58] with permission of the Physical Society of the Republic of China.

Figure 4 .
Figure 4. Temporal behavior in the populations of the quasi-degenerate excited states for (a) e (+) , (b) e (−) , (c) e+1, and (d) e−1 calculated from nuclear WP simulations.In each panel, the populations of Xand Y are denoted by the solid and dotted lines, respectively.The laser pulses fully decay at t = 7.26 fs.Reprinted from[56] with permission of the American Chemical Society.

Figure 4 .
Figure 4. Temporal behavior in the populations of the quasi-degenerate excited states for (a) e (+) , (b) e (−) , (c) e +1 , and (d) e −1 calculated from nuclear WP simulations.In each panel, the populations of |X and |Y are denoted by the solid and dotted lines, respectively.The laser pulses fully decay at t = 7.26 fs.Reprinted from[56] with permission of the American Chemical Society.

Figure 5 .
Figure 5. Expectation value of the electronic angular momentum lz(t) in DCP irradiated by (a) linearly and (b) circularly polarized VUV laser pulses calculated from nuclear WP simulations.In panel (a), the expectation values for e(+) and e(−) are denoted by the solid and dotted lines, respectively; in panel (b), those for e+1 and e−1 are denoted by the solid and dotted lines, respectively.The laser pulses fully decay at t = 7.26 fs.Reprinted from[56] with permission of the American Chemical Society.

Figure 6
Figure 6 displays the expectation value of the normal coordinates      

Figure 5 .
Figure 5. Expectation value of the electronic angular momentum l z (t) in DCP irradiated by (a) linearly and (b) circularly polarized VUV laser pulses calculated from nuclear WP simulations.In panel (a), the expectation values for e(+) and e(−) are denoted by the solid and dotted lines, respectively; in panel (b), those for e +1 and e −1 are denoted by the solid and dotted lines, respectively.The laser pulses fully decay at t = 7.26 fs.Reprinted from[56] with permission of the American Chemical Society.

Figure 6
Figure 6 displays the expectation value of the normal coordinates Q(t) ≡ Ψ(t) Q Ψ(t) with Q

Figure 6 .
Figure 6.Expectation value of the normal coordinates Q(t) of the breathing and distortion modes in DCP irradiated by (a) linearly and (b) circularly polarized VUV laser pulses.In panel (a), the expectation values for e (+) and e (−) are denoted by the solid and dotted lines, respectively; in panel (b), those for e+1 and e−1 are denoted by the solid and dotted lines, respectively.The laser pulses fully decay at t = 7.26 fs.Reprinted from[56] with permission of the American Chemical Society.


Q mainly consists of low (high) vibrational states for e (+) [e (−) ].The spectral features for circular ones in Figure 7b are very similar to each other: The primary peaks of   X  for e+1 and e−1 are both at  ~ 1400 cm −1 , while the intensities of the other peaks are a little stronger in the latter case.This indicates that the same frequency components for e+1 and e−1 in the time domain after the nonadiabatic transition.

Figure 6 .
Figure 6.Expectation value of the normal coordinates Q(t) of the breathing and distortion modes in DCP irradiated by (a) linearly and (b) circularly polarized VUV laser pulses.In panel (a), the expectation values for e (+) and e (−) are denoted by the solid and dotted lines, respectively; in panel (b), those for e +1 and e −1 are denoted by the solid and dotted lines, respectively.The laser pulses fully decay at t = 7.26 fs.Reprinted from[56] with permission of the American Chemical Society.

Figure 7 .
Figure 7.The frequency spectra of ψX(Q, t), σX(ω), defined by Equation (34) for DCP irradiated by (a) linearly and (b) circularly polarized VUV laser pulses.In panel (a), the spectra for e(+) and e (−) are denoted by the solid and dotted lines, respectively; in panel (b), those for e+1 and e−1 are denoted by the solid and dotted lines, respectively.In each case, the values of σX(ω) were scaled so that the maximum value is unity.Reprinted from[56] with permission of the American Chemical Society.

Figure 8
Figure 8 shows the propagations of the adiabatic WPs on the two excited-state PESs.The probability densities

Figure 7 .
Figure7.The frequency spectra of ψ X (Q, t), σ X (ω), defined by Equation(34) for DCP irradiated by (a) linearly and (b) circularly polarized VUV laser pulses.In panel (a), the spectra for e(+) and e (−) are denoted by the solid and dotted lines, respectively; in panel (b), those for e +1 and e −1 are denoted by the solid and dotted lines, respectively.In each case, the values of σ X (ω) were scaled so that the maximum value is unity.Reprinted from[56] with permission of the American Chemical Society.
a Gaussian especially at t ~ 12 fs.As a result of nonadiabatic couplings,

Figure 8 ..
Figure 8. Propagation of the adiabatic WPs on the two-dimensional potential energy surfaces (PESs) of X and Y for (a) e (+) , (b) e (−) , (c) e+1, and (d) e−1.The origin of the PESs is the optimized geometry of G and the avoided crossing is designated by a circle.The bold contours represent the probability densities nonadiabatic transition from X to Y , which interferes with   , Y t  Q , is smaller than that for the transition from Y to X , which interferes with   , X t  Q ; this is the reason why the resultant upward population transfer is small in Figure 4a.The WPs on the two PESs then reach the avoided crossing and the reverse population transfer occurs around t ~ 10-14 fs.The interference enhances low vibrational states in

Figure 8 .
Figure 8. Propagation of the adiabatic WPs on the two-dimensional potential energy surfaces (PESs) of |X and |Y for (a) e (+) , (b) e (−) , (c) e +1 , and (d) e −1 .The origin of the PESs is the optimized geometry of |G and the avoided crossing is designated by a circle.The bold contours represent the probability densities |ψ X (Q, t)| 2 and |ψ Y (Q, t)| 2 .The arrows indicate the motion of the center of the WPs.Modified from[56] with permission of the American Chemical Society.
the nabra operating on the AO on the right-hand (left-hand) side.Surface integration in Equation (

Figure 9 .
Figure 9. (a) Three electronic excited states and optical excitations to create three types of the electronic coherence of (P)-2,2′-biphenol (b) The geometrical structure and the electronic transition moment vectors.Modified from[87] with permission of the American Institute of Physics.

Figure 9 .
Figure 9. (a) Three electronic excited states and optical excitations to create three types of the electronic coherence of (P)-2,2 -biphenol (b) The geometrical structure and the electronic transition moment vectors.Modified from[87] with permission of the American Institute of Physics.

Figure 11 .
Figure11.Bond currents for the three types of the electronic coherence at a time when the maximum electronic coherence is created.Modified from[87] with permission of the American Institute of Physics.

Figure 11 .
Figure11.Bond currents for the three types of the electronic coherence at a time when the maximum electronic coherence is created.Modified from[87] with permission of the American Institute of Physics.
αβ with α, β = a, b 1 , b 2 .Here, e αβ • µ gα = e (−) αβ • µ gβ , respectively.The excitation by laser pulses with e αβ ] creates a linear combination of two electronic states α and β with their in-phase (out-of-phase) relation, which is called e (+) αβ [e (−) αβ ] excitation hereafter.The schematic drawings on the left-hand side of Figure 12a represent the initial rotational directions of π electrons after the preparation of a superposition of the two excited states, b 1 and b 2 .The upper (lower) drawing shows the rotational direction just after coherent e (+) b1b2 [e (−) b1b2 ] excitation, which indicates CC (AA) rotation, i.e., π electrons start to rotate along the two phenol rings in the same directions.

e
rotational directions of π electrons after the preparation of a superposition of the two excited states, b1 and b2.The upper (lower) drawing shows the rotational direction just after coherent ] excitation, which indicates CC (AA) rotation, i.e., π electrons start to rotate along the two phenol rings in the same directions.

Figure 12 .
Figure 12.(a) Four patterns of π-electron rotations induced by selective coherent electronic excitations.For example, explained in this figure caption rotation means clockwise and anticlockwise rotations along the phenol rings L and R, respectively.(b) Resultant angular momenta lZ and lX for the four rotational patterns.lL (lR) denotes the electronic angular momentum of the L (R) ring.Modified from[86] with permission of the American Chemical Society.

Figure 12b shows theFigure 12 .
Figure 12b shows the resultant angular momentum l = lL + lR.Here, the π-electron angular momentum lL (lR) is perpendicular to the L (R) phenol ring, and the resultant total angular momentum lZ (lX) is parallel to the Z-(X-) axis on the ZX-plane.The directions of lZ indicate that clockwise (anticlockwise) rotation around the Z-axis is created by     e [     e ] excitation.The schematic drawings on the right-hand side of Figure 12a (12b) represent π-electron rotations (the resultant angular αβ ] excitation.The schematic drawings on the right-hand side of Figure12a(12b) represent π-electron rotations (the resultant angular momenta) after preparation of a superposition of the two excited states, a and b 1 , with different irreducible representations.π Electrons start to rotate along the two phenol rings in opposite directions to each other (CA or AC rotation).Therefore, the total angular momentum parallel to the Z-axis is zero and the nonzero resultant angular momentum l X is generated along the X-axis on the ZX-plane.If the two phenol rings are in the same plane without twist, both l X and l Z vanish.Furthermore, e (+) ab1 and e (−) ab1

Figure 14 .
Figure 14.Superposition of coherent π-electron rotation patterns.The arrow with positive (negative) sign specifies an in-phase (out-of-phase) superposition of [AA] and [AC] rotation patterns.The inphase superposition creates [AO] rotation pattern, which means anticlockwise rotation localized on the L ring with no rotation on the R ring, while the out-of-phase superposition creates [OA] rotation pattern.Here, O means no rotation.Reprinted from [91] with permission of the American Chemical Society.
OA] localization.Here, O means no rotation on the L or R ring.

Figure 14 .
Figure 14.Superposition of coherent π-electron rotation patterns.The arrow with positive (negative) sign specifies an in-phase (out-of-phase) superposition of [AA] and [AC] rotation patterns.The in-phase superposition creates [AO] rotation pattern, which means anticlockwise rotation localized on the L ring with no rotation on the R ring, while the out-of-phase superposition creates [OA] rotation pattern.Here, O means no rotation.Reprinted from[91] with permission of the American Chemical Society.

Figure 15 .
Figure 15.Results of quantum localization of π-electron rotations for the target [AO] and for the target [OA].(a) Calculated time-dependent π-electron angular momenta lL(t) (blue solid line) and lR(t) (purple dashed line).(b) Imaginary parts of electronic coherences, ρab1(t), ρb1b2(t), and ρab2(t).(c) Real parts of the optimal electric field in the frequency domain.Results of the [AO] ([OA]) target are shown in solid (dashed) lines.Reprinted from[91] with permission of the American Chemical Society.

Figure 15 .
Figure 15.Results of quantum localization of π-electron rotations for the target [AO] and for the target [OA].(a) Calculated time-dependent π-electron angular momenta l L (t) (blue solid line) and l R (t) (purple dashed line).(b) Imaginary parts of electronic coherences, ρ ab1 (t), ρ b1b2 (t), and ρ ab2 (t).(c) Real parts of the optimal electric field in the frequency domain.Results of the [AO] ([OA]) target are shown in solid (dashed) lines.Reprinted from[91] with permission of the American Chemical Society.

Figure 16 .
Figure 16.(a) Expectation values of π-electron angular momentum for the free rotations after the [AO] localization.(b) Results of the optimal control for maintaining the localization.The control time is tf = 30 fs.(c) Absolute values of the optimal electric field.The inset figure is the field in the restricted time domain.Reprinted from[91] with permission of the American Chemical Society.

Figure 16 .
Figure 16.(a) Expectation values of π-electron angular momentum for the free rotations after the [AO] localization.(b) Results of the optimal control for maintaining the localization.The control time is t f = 30 fs.(c) Absolute values of the optimal electric field.The inset figure is the field in the restricted time domain.Reprinted from[91] with permission of the American Chemical Society.
General Form of Polarized Laser PulsesBefore entering into the analysis of optical excitation processes in aromatic molecules, let us describe the mathematical form of a classical laser field of arbitrary polarization.It is given in the dipole approximation by: