Quantum Control of Coherent π-Electron Dynamics in Aromatic Ring Molecules

Herein we review a theoretical study of unidirectional π-electron rotation in aromatic ring molecules, which originates from two quasi-degenerate electronic excited states created coherently by a linearly polarized ultraviolet/visible laser with a properly designed photon polarization direction. Analytical expressions for coherent π-electron angular momentum, ring current and ring current-induced magnetic field are derived in the quantum chemical molecular orbital (MO) theory. The time evolution of the angular momentum and the ring current are expressed using the density matrix method under Markov approximation or by solving the time-dependent Schrödinger equation. In this review we present the results of the following quantum control scenarios after a fundamental theoretical description of coherent angular momentum, ring current and magnetic field: first, two-dimensional coherent π-electron dynamics in a non-planar (P)-2,2’-biphenol molecule; second, localization of the coherent π-electron ring current to a designated benzene ring in polycyclic aromatic hydrocarbons; third, unidirectional π-electron rotations in low-symmetry aromatic ring molecules based on the dynamic Stark shift of two relevant excited states that form a degenerate state using the non-resonant ultraviolet lasers. The magnetic fields induced by the coherent π-electron ring currents are also estimated, and the position dependence of the magnetic fluxes is demonstrated.


INTRODUCTION
Recent progress in laser science and technology has facilitated the coherent control of ultrafast charge migration dynamics in molecular systems . Controlled charge migration can generate unidirectional currents and current-induced magnetic fields, which can be used as a guiding principle for the next-generation ultrafast optoelectronic devices [22,23]. Laser control of π-electron rotations in aromatic ring molecules is a typical example. Pioneering studies on the generation of coherent π-electron rotations in high-symmetry ring molecules such as Mg porphyrin having degenerate excited states have been reported by the Manz group [24][25][26]. π-Electron rotations were induced by degenerate electronic excited states, which were coherently created by circularly polarized ultraviolet (UV) lasers. The generation mechanism of the photon angular momentum involves the photon transfer from the circularly polarized UV lasers to the π-electrons, thereby the left-or right-handed circular polarization of the applied laser defines the rotational direction of the angular momentum [27]. In contrast in ring molecules with low-symmetry π-electrons cannot be rotated by using the circularly polarized lasers, due to the absence of electronic excited states degenerate, that would receive the photon angular momentum. Thus, it was commonly understood that coherent π-electron rotations could not be generated in low-symmetry aromatic ring molecules. However, Kanno et al. [28][29][30][31] have invalidated this long-established understanding by demonstrating that π-electrons in oriented chiral aromatic ring molecules can be rotated by the coherent excitation of a pair of quasi-degenerate π-electronic excited states using a linearly polarized UV/vis laser pulse with a properly designed photon polarization direction [28,29]; the polarization direction of the pulse determines the initial direction of the π-electron rotation, whether right-or left-handed one. The duration of the unidirectional rotation is inversely proportional to the energy difference between the two quasi-degenerate excited states, and inverse rotation begins after the duration because the coherent state is not a nonstationary state rather than an eigenstate. Pump and dump laser pulses with their properly designed polarization directions for these lasers are applied to eliminate the reverse rotation. The number of unidirectional rotations during the duration can be estimated from the energy difference between the quasi-degenerate excited and ground states. This only applies to the ideal case in which any dephasing processes disturbing the electronic coherence are omitted. It is expected that unidirectional ring currents produce much stronger magnetic fields than traditional static magnetic fields [32,33].
The concept of the conventional ring current has already been established previously and plays an important role in interpreting the magnetic properties and aromaticity of conjugated molecular systems [34,35]. The molecules investigated based on the ring current are in the ground state, and the current densities for evaluation of the ring current are calculated using the first-order electronic wave function in the time-independent magnetic field and in a permanent magnetic dipole. Such a conventional ring current should be called an incoherent current.
In this article, we present the most recent results of our theoretical studies on the quantum control of coherent π→electron rotations in low-symmetry aromatic ring molecules. Low-symmetry aromatic ring molecules are not rare, but rather common: high-symmetry aromatic molecules often become low-symmetry ones by substitution of a functional group, or under environment conditions. In the next section, the fundamental theory for the quantitative analysis of coherent π-electrons in a low-symmetry aromatic ring molecule is introduced. Here, the time-dependent coherent angular momentum, ring current, and ring current-induced magnetic field are analytically derived in a closed form in the quantum chemical MO theory. Time-evolutions of these quantities are derived by the density matrix method [36][37][38][39]. The temporal behavior of coherence is determined by the off-diagonal element of density matrix. Then, the Markov approximation is ultilized for considering the dephasing effects on the conherent angular momenta and ring currents. The magnitudes of the electronic angular momenta and ring currents are expressed as the summation of the expectation values of the corresponding one-electron operators in the aromatic rings. The bond current between the nearest neighbor carbon atoms, C i and C j , is defined as an electric current flowing through a half plane perpendicular to the C i − C j bond. The coherent bond current in an aromatic ring is defined as the average of all bond currents. The application of this theory to a nonplanar chiral aromatic molecule, (P)-2,2'-biphenol, is briefly described. (P)-2,2'biphenol has four patterns of coherent π-electron rotations along with the two phenol rings because of its nonplanar geometrical structure [39,40]. A sequence of the four rotational patterns can be controlled through a coherent excitation of two electronic states with two requirements: the symmetry of the two electronic states, and their relative phase. Quantum switching of coherent π-electron rotations is proposed [40]. In Section 3, the key points are summarized for the application of the quantum optimal control method for controlling coherent ring currents in polycyclic aromatic hydrocarbons (PAHs). These molecules exhibit several localization patterns of coherent π-electron rotations. Therefore, how to set up the target state for a desired localization pattern is crucial [41,42]. However, we demonstrate that the target state can easily be set up using the Lagrange multiplier method. As examples, two types of current localizations for the simplest PAH, anthracene, are adopted [42]: current localization to the designated benzene ring, and the perimeter ring current. In Section 4, a convenient scenario involving unidirectional π-electron rotations in low-symmetry aromatic ring molecules is described [43]. The basic idea behind unidirectional electron rotations is to degenerate two nondegenerate excited states by utilizing dynamic Stark shifts. A degenerate state induced by dynamic Stark shifts is called a dynamic Stark-induced degenerate electronic state (DSIDES). Two linearly polarized continuous lasers with different frequencies and phases are used to form a DSIDES, where each laser is set to selectively interact with each electronic state through non-resonant excitation. As a result, the unidirectional π-electron rotation is driven by the lasers. This scenario was applied to toluene. The resulting angular momenta can be represented by a pulse train having an angular momentum. Each angular momentum pulse represents the unidirectional π-electrons rotation, which begins with acceleration and ends with deceleration.

Equations of Motion for Coherent π-Electron Rotations Induced by Ultrashort UV/Vis Lasers
The expectation values of the coherent π−electron angular momentum and the ring current operators in an aromatic molecule are generally expressed as Here, O(r) is a single-electron operator for the angular momentum or current. Ψ(t) is the wave function of the π−electrons in laser field F(t) at time t, n denotes the number of electrons, and r i express the ith electron coordinates. Since the optically-allowed electronic excited states of the conjugated aromatic rings are of our interest, the electronic wave function can be expanded by the two electronic configurations, i.e., the ground configuration Φ 0 , and the singly excited ones Φ α as where Φ 0 is defined by a single Slater determinant as Φ 0 (r 1 , /, r n ) ϕ 1 /ϕ a /ϕ b /ϕ n with ϕ n ≡ ϕ n (r n ). Here, ϕ a and ϕ b are the occupied orbitals, Φ α is the electronic wave function for a single electron excited configuration α : a → a′, i.e., the single electron transition from the occupied molecular orbital (MO) a to the unoccupied MO a′, when Φ α (r 1 , /, r n ) ϕ 1 /ϕ a' /ϕ b /ϕ n . The coherent and incoherent temporal behaviors of the electrons induced by a laser field F(t) can be obtained directly by solving the coupled electronic equations of motion expressed by the density matrix elements ρ αβ (t) under the initial population conditions ρ 00 (0) 1 and ρ αα (0) 0 for α ≠ 0 and ρ αβ (0) 0 for α ≠ β, that is, there is no electronic coherence at the initial time such that Here, ρ αβ (t) ≡ c α (t)c p β (t) is the density matrix element, and V αc (t) is the coupling matrix element between states α and c via the molecule-laser interaction V(t) −μ · F(t), where μ denotes the transition dipole moment operator, c αβ are the dephasing constants in the Markov approximation [38,39] and ω βα is the frequency difference between the two electronic states α and β. Here, c α (c β ) is the non-radiative transition rate constant of the electronic state α(β), and c (d) αβ are the pure dephasing constants induced by the elastic interaction between the heat bath and molecule of interest.
Because we are interested in the coherent behaviors of dipoleallowed quasi-degenerate π-electronic excited states in the visible or UV region of an aromatic ring molecule, Eq. 1 can be rewritten in terms of singly excited configurations {Φ κ } if κ ≠ 0 as where In Eq. 4, the coherence between singly excited state configurations is considered, and the coherence between the ground state and the excited state is neglected because this coherence is much shorter compared to that between the singly excited state configurations. In this case, the coherence time is proportional to the inverse of the energy gap between the two electronic states.
The coherent π-electron angular momentum is formulated in the quantum chemical MO theory, and the π-orbital ϕ k associated with the optical transition is expanded in terms of a linear combination of atomic orbitals χ i as where χ i denotes the atomic orbital, and c k,i indicates the molecular orbital coefficient. Equation 4 can be expressed by using the Eq. 5 as Here, the temporal behavior of the expectation value 〈 O(t)〉 can be expressed using the off-diagonal density matrix element ρ βα (t). The suffixes (a, a') ((b, b')) corresponds to the electronic configurations α(β), respectively.

Coherent π-Electron Angular Momentum
Consider the spatially fixed aromatic ring molecule with N aromatic rings. The electron angular momentum operator is defined as Here l Z,K is the electronic angular momentum operator of the component perpendicular to ring K. Coordinates x K and y K are defined on the aromatic ring K, and n ⊥,K is the unit vector perpendicular to the aromatic ring. The expectation value of Kth angular momentum operator is given in terms of a 2p z carbon atomic orbitals (AOs) as Note that in Eq. 7 the dependence of the electron angular momentum on the laser intensity is reflected in the imaginary part of the off-diagonal density matrices Im ρ βα (t).

Coherent π-Electron Ring Current
From Eq. 6, the perpendicular component of the time-dependent electric current flowing through a surface S (See Figure 1) is generally defined as Here is the nabla operating on the atomic orbital on the right-hand side (left-hand side).
Equation 8 can be expressed as where Im ρ βα (t) δ ab c p a′i c b′j − c p b′i c a′j + δ a'b' c p ai c bj − c p bi c aj J ij , where n ⊥,S is a unit vector perpendicular to a surface S, which is given as n ⊥,S r j −r i |rj−ri| , and the surface S is set at the center of the carbon bond C i − C j . The surface integrations in Eqs. 8, 10b are carried out over the half-plane S (see Figure 1). Using Slater type AOs for {χ i }, J ij in Eq. 10b can be expressed in analytical form [40]. We introduce the bridge bond current J B,pq (t) which is a specific case of an interatomic bond current, and is defined as the current bridging two aromatic rings P and Q, passing from the nearest neighbor carbon atom at site p to q. Here, each carbon atom C p , or C q belongs to a different aromatic ring P or Q which are in the neighbors. (see Figure 1).
The bridge bond current is given in terms of the interatomic current, J pq (t), as Here, θ d is the dihedral angle between the two rings P and Q [39]. We now define an averaged ring current along ring K, J K (t), which is defined by taking the average of all bond currents as where n K is the number of bonds in ring K. respectively. (C) Linearly polarized unit vector of laser e + (e − ) used to generate the in-phase superposition (a+b 1 ) (out-of-phase superposition (a-b 1 )).
Frontiers in Physics | www.frontiersin.org July 2021 | Volume 9 | Article 675134 2.4 Application to Nonplanar (P)-2,2'-biphenol (P)-2,2'-biphenol is one of the typical nonplanar chiral aromatic ring molecules, which has two aromatic rings connected through the C-C bridge bond. For convenience, hereafter, we denote L (R) for the left-(right-) hand side phenol ring of (P)-2,2'-biphenol. Figure 2A exhibits the geometrical structure of (P)-2,2'-biphenol with the transition dipole moments vectors between the ground and excited states. (P)-2,2'-biphenol was assumed to be fixed on a surface by a non-conjugated chemical bond or was oriented in the space by the molecular orientation techniques by laser [44][45][46]. The laboratory-fixed Y-axis was set parallel to the single chemical bond bridging two phenol molecules, and the rotational axis of the molecule belonging to point group C 2 was set along the laboratory-fixed Z-axis which was parallel to the surface normal. The ground state geometry of (P)-2,2'biphenol was optimized using the DFT-B3LYP level theory with the 6-31G+(d,p) basis set in the GAUSSIAN09 code [47]. The dihedral angle θ d between the two phenol rings, was found to be 108.9 + based on the density functional theory (DFT) calculations. The calculated geometrical structures are provided in Refs. [39,40]. The dephasing constants were set as c b1b2 c ab1 γ ab2 0.01 eV (∼1/50 fs −1 ). The arrows above the C -C bonds indicate the initial directions of the currents. Note that the bridge bond current J 1,7 0 for the (b 1 +b 2 ) electronic coherence, while J 1,7 ≠ 0 for the (a +b 1 ) and (a +b 2 ) electronic coherences. Reprinted with permission from Ref. [39] Copyright (2013) American Institute of Physics.
Frontiers in Physics | www.frontiersin.org July 2021 | Volume 9 | Article 675134 To create the coherent angular momenta and ring currents in (P)-2,2'-biphenol, we focused on the three dipole-allowed electronic excited states (a, b 1 , and b 2 ) as shown in Figure 2B. The transition energies from the ground (g) to the a, b 1 , and b 2 states, which were calculated under the optimized ground state geometry using the TD-DFT B3LYP level of theory [39,40] were 6.67, 6.78, and 6.84 eV, respectively.

Generation of Two-state Electronic Coherence Using Linearly Polarized UV Pulses
For a creation of the coherent angular momentum and ring current using linearly polarized laser pulses, it is essential to prepare for an electronic coherence with a fixed relative phase between two electronic states, a and b 1 , i.e., in-phase (a + b 1 ) or out-of-phase (ab 1 ). The principle behind the preparation of the electronic coherence (a + b 1 ) ((a − b 1 )) using a linearly polarized laser with polarization unit vector e + (e -) is schematically demonstrated in Figure 2C [39,40].
In the frame work of the three-excited state model, there are three types of two-electronic coherent states represented as b 1 ± b 2 , a ± b 1 and a ± b 2 , respectively. At the initial time each electronic coherence can be generated by applying a linearly polarized UV laser with a properly selected laser polarized direction. Because of the non-planar geometry of the molecule, the angular momentum is two-dimensional, and the two ring currents flow are on the two different planes, the direction of total angular momentum is dependent on the symmetry of the coherent state created by the laser: The Z-directional angular momentum (ring current) is generated for the created coherent state with the A-irreducible representation of the C 2 point group, while the X-directional angular momentum (ring current) is generated for the coherent state with B.
It is remarkable that even though the third excited state is located between the two excited states in the three-excited states system, the coherent electronic state can still be created if the applied linearly polarized lasers satisfy the following conditions. For example, for the (a ± b 2 ) electronic coherence, the conditions for the linearly polarized lasers with polarization vectors e ± are given as μ ga · e ± ± μ gb2 · e ± and μ gb1 · e ± 0, or equivalently, Thus, if the laser overlaps the electronic states a, b 1 , and b 2 , the (a + b 2 ) or (ab 2 ) electronic coherent state can be created selectively. Figure 3 exhibits the temporal behaviors of the angular momenta calculated for three types of electronic coherences, (b 1 +b 2 ), (a + b 1 ) and (a + b 2 ), each of which is created by a linearly polarized UV laser pulse with a properly selected polarization direction of laser. For the (b 1 +b 2 ) electronic coherence, the total angular momentum parallel to the Z-axis, l Z , is created together with two ring components l L l R , whereas for both the (a + b 1 ) and (a + b 2 ) electronic coherences, the total angular momenta parallel to the X-axis, l X , are generated with l L − l R , with a π-phase shift. Similarly, for out-of-phase electronic coherences,

Coherent Angular Momentum Quantum Beats and Bond Currents
, the angular momentum can be given by the corresponding in-phase electronic coherence with a π-phase shift. The simple sinusoidal temporal behavior originates from the coherence of two electronic states with the oscillation period corresponding to a frequency difference between the twoelectronic states. This is referred to as the angular momentum quantum beats, which are similar to the fluorescence quantum beats which originate from the vibronic coherence [48]. We note that π-electrons rotate a few times in a unidirectional manner within a half cycle of the oscillation. This unidirectional π-electron rotation can produce a unidirectional ring current and corresponding current-induced magnetic flux. In principle, this enables the design of ultrafast switching devices which consist of organic aromatic ring molecules.
The angular momenta l Z (t) and l X (t), shown in Figure 3, are obtained by the summation of the angular momenta created in the L and R aromatic rings l L (t) and l R (t) using the following relationship l Z (t) 2l L (t) · e Z sin θ d 2 with l L (t) · e Z l R (t) · e Z , and l X (t) 2l L (t) · e X cos θ d 2 (where l L (t) · e X −l R (t) · e X ). Table 1 lists the angular momenta of (P)-2,2'-biphenol with the dihedral angle θ d 108.9 + where the electronic coherence is maximum at time t tp, i.e., when the magnitude of the imaginary part of the density matrix element is maximized (Im ρ αβ (tp) − 1 2 ) and dephasing effects are neglected. This results in l X (tp) > 0 or l X (tp) < 0 for the (a + b 1 ) and (a + b 2 ) electronic coherences, respectively, and in l Z (tp) < 0 for the The magnitudes of the bond current, J ij calculated at the maximum coherence time, are presented in detail in [40]. The magnitudes of the averaged ring current over the C-C bonds at the maximum coherence time, J are on the order of tens of μA, i.e., J 161, 86.5 and 63.4 μA for (b 1 + b 2 ), (a + b 1 ), and (a + b 2 ) coherences, respectively.
Effects of dephasings on coherent π-electron angular momentum and ring currents were treated in the Markov approximation, and time-independent dephasing constants were used under the assumption of instantaneous interactions between the system and phonon baths. In a system such as condensed phases, the Markov approximation is broken down. Non-Markov 1 | Angular momenta of the two phenol rings l L ≡ l L (tp) · n L and l R ≡ l R (tp) · n R , and the resulting angular momenta l X and l Z at the maximum coherence time t t p1) .
Frontiers in Physics | www.frontiersin.org July 2021 | Volume 9 | Article 675134 response of coherent should be essential. Time evolution of coherent ring currents were calculated in a hierarchical master equation approach beyond the Markov approximation has been treated [49].

Ultrafast Quantum Switching of Angular Momentum
Consider the quantum control of π-electron rotations for twodimensional angular momentum switching based on the results shown in Table 1. Here, two-dimensional quantum switching is defined as a sequential pulse of the electronic angular momentum with its constant sign (positive or negative) along the Z-or X-axis. Note that any quantum switching step should be completed before the reverse rotation of the π-electrons begins, because it may otherwise disturb the signal. Consider a sequential four-step control, which is expressed as l Z (−)→l X (+)→l Z (+)→l X (−). This indicates the switching of rotational patterns in the order CC → AC → AA → CA, where the symbol C (A) means clockwise (anti-clockwise) direction, and for example, CA rotation means clockwise and anti-clockwise rotations along the phenol rings L and R, respectively. Here, l X (+) (l Z (−)) means the π-electron angular momentum along the X-(Z-) axis with a positive (negative) sign, i.e., anti-clockwise (clockwise) rotation of π-electrons around the corresponding axis. Figure 4A presents a 3D plot of the angular momentum as it switches based on the sequential four-step control scheme. From Figure 4A, we can see that the π-electron rotations are successfully manipulated by the pulses depicted in Figure 4B. That is, both the rotational axis parallel to the Z-or X-axis and the rotational directions around those axes, clockwise or anticlockwise, are manipulated by the sequential four-step process. In Figure 4B, each switching step of control was performed using pump and dump pulses with specific polarization directions and phases. The laser pulse with an amplitude of 1.2 GV/m ( 1.9 × 10 11 W/cm 2 ) was used in the second and fourth steps. The dynamic Stark shifts between electronic states a, b 1 and b 2 were on the order of 0.01 eV [39], indicating that the Stark effects could be omitted in Figure 4A.
The pulses shown in Figure 4B have two features. The first feature is that the pump (dump) pulse for each step has polarizations, e (−) αβ (e (+) αβ ) or e (+) αβ (e (−) αβ ). Each pulse has an energy width that is large enough to coherently excite two quasi-degenerate electronic excited states, as shown in Figure 2. The second feature is that the pump and dump laser pulses partially overlap. In the creation of the CC rotation, for example, the electric field of the pump pulse was Here, T b1b2 ( 60.9 fs) indicates the oscillation period between b 1 and b 2 states, E 0 b1b2 is the amplitude of the laser pulse, ω c,b1b2 is the central frequency between the two excited states b 1 and b 2 ; and t pd b1b2 (the time interval between the pump and dump laser pulses) was set to T b1b2 /2 (Supporting information in Ref. [39]). The angle between the two polarization vectors, e (+) b1b2 and e (−) b1b2 , was 113.5°. With respect to ultrafast quantum switching, it is crucial to create the overlap of the pump and dump pulses, where the resulting electric field rotates as an elliptically polarized electric field in the overlapped time domain, and the dump laser pulse reverse the rotation that occurs during this region. As a result, the angular momentum of the π-electrons is nullified.

Coherent Ring Current-Induced Magnetic Field
There have been interesting reports on the evaluation of the magnetic fields of atoms and oriented heteronuclear diatomic molecules, AlCl and BeO [32,50,51]. Strong and unidirectional magnetic fields are generated from the degenerate electronic states of these atoms and molecules excited by circularly polarized intense laser pulses. We estimated the magnetic fields (magnetic field flux density) generated by the coherent ring currents of (P)-2,2'-biphenol. The results may provide fundamental information for designing ultrafast switching devices controlled by current-induced magnetic fields as well as coherent ring currents [52][53][54].
As an example, consider the magnetic field induced by the ring current for the (b 1 + b 2 ) electronic coherence. In Figure 5, B K (tp, h) represents the current-induced magnetic field along the central axis perpendicular to ring K (L or R) as a function of the height h above the Z-axis at t tp under the maximum coherence condition (when Im ρ b1b2 (tp) −1/2). An expression for the magnetic field B K (tp, h) was derived by taking into account the interatomic bond currents with the 2p z Slater AOs [40]. Note that B L (tp, h) B R (tp, h) for the (b 1 + b 2 ) electronic coherence ( Table 1). It is interesting to compare the behaviors of B K (tp, h) with those calculated using which was derived using a simple ring loop (SRL) model. Here, μ 0 4π · 10 −7 [W b /(A · m)] is the magnetic constant, r ( r ij ) is the ring radius, and η ( sin −1 (r/ r 2 + h 2 √ )) is defined as the angle between the Z-axis and a straight line depicted from the point on ring K that crosses the Z-axis. B K (t, 0) at η π/2 is the magnetic field measured at the center of ring K. The magnitude of B K (t p , 0) induced by J K (t p ) 100 μA is 448 mT at the center of the ring K with r 0.14 nm. It can be found from Figure 5 that the magnitudes of B SRL K (t p , h) are overestimated near the aromatic ring plane 0 ≤ h < 1 Å, while the magnitudes are reasonable for h > 1 Å, although slightly different results can be observed between the two magnetic fields for large values of h. At the center of the aromatic ring, B K (t p , h 0) 0.66B SRL K (t p , h 0). This can be understood from the result that the π-electron current density is densely distributed over the aromatic ring. A subtle difference between B SRL K (tp, h)and B K (tp, h) is observed for large h, although both magnetic fields B SRL K (tp, h → ∞) 0 and B K (tp, h → ∞) 0 approach zero. This deviation originates from the approximation that the aromatic ring is not a considered to be a perfect ring.
It is necessary to check whether the value of B K (t, 0) is larger than the one corresponding to the magnetic field B Laser induced by the applied laser field F when the induced magnetic fields are created during an ultrasfast laser pulse excitation at the early time regime. The magnitude of magnetic field B Laser can be estimated from a simple formula |B Laser | |F|/c with c 3.0 × 10 8 ms −1 . The calculated magnitude of magnetic field B Laser with |F| 1.0 GV/ m is approximately 7.4B K (tp, h), which is on the same order as the field induced by the ring current J 100 μA. This implies that we need a careful examination to observe the current-induced magnetic flux, or to use the electro-magnetic device as a switching control tool.
Thus far, we have taken into account the ring current-induced magnetic fields of the low-symmetry aromatic ring molecule (P)-2,2'-biphenol, in which nondegenerate two electronic excited states are coherently excited by the linearly polarized UV lasers. Here we briefly discuss the ring current-induced magnetic field for the degenerate electronic excited states of an aromatic molecule induced by the intense circularly polarized UV laser pulse. Yuan and Bandrauk [2,33] have numerically demonstrated that the circularly polarized ultrashort pulses are generated from the molecular high-order harmonic generation using the intense linearly and circularly polarized laser pulses. This indicates a possibility for the ultrashort circularly polarized UV laser pulses to create the ring current-induced magnetic fields in high symmetric aromatic ring molecules. We can estimate the magnitude of the ring current-induced magnetic field of benzene within the SRL model using Eq. 15. Here, J(t) is the electric ring current of the degenerated electronic state. The electronic spectrum of benzene is characterized by the dipole-allowed transition from the ground state to the third singlet electronic excited state ( 1 E 1u ). For an equal population between the two states at t tp, B SRL (tp, 0) 874 mT is obtained using the maximum value of the coherent ring current J(tp) 195 μA evaluated under the π-electron approximation [40,55]. For comparison, we obtained B SRL K (tp, 0) 874 mT and B K (tp, 0) 579 mTwhen J(tp) 195 μAfor (P)-2,2'-biphenol. Note that the same magnitudes of the induced magnetic field B SRL K (tp, 0) for (P)-2,2'-biphenol and B SRL (tp, 0) 874 mT for benzene were obtained in the SRL model because the two parameters in the SRL model, the radius of the ring r, and the ring current J(t p ) have the same values for both types of ring molecules. The same tendency in the magnitudes between B SRL K (t p , 0) and B K (t p , 0) is also observed for benzene.

CURRENT LOCALIZATION IN POLYCYCLIC AROMATIC HYDROCARBONS
In this section we consider a localization of coherent ring current in polyatomic aromatic hydrocarbons (PAH) [41,42]. There exist Frontiers in Physics | www.frontiersin.org July 2021 | Volume 9 | Article 675134 various current localization patterns in PAH simply because of their geometrical structures consisting of many benzene units. Quantum optimal control method is a general and reliable one to choose a desired current localization pattern from the various patterns using control lasers. The quantum control method has successfully been applied to manipulation of molecules such as coherent control of chemical reactions [56,57]. After a brief introduction of quantum optimal control method [41,42,58], we demonstrate how the target state is set up for the designated ring current using the Lagrange multiplier method. Here, the π-electron ring current is expressed in terms of the interatomic bond currents between two adjoining C-C atoms. At final, the target states are derived for the two types of current localization: the localized ring current, which indicates that the ring current is localized to the designated aromatic ring in linear PAH, and the perimeter ring current in linear PAH [41,42].
where H 0 is a Hamiltonian in absence of field, and Ψ(t) is the time-dependent wave function in the electric field F(t). Here, the The penalty factor α 0 is introduced to suppress the magnitude of the electric field F(t). ξ(t)is the time-dependent Lagrange multiplier. The third term in Eq. 16 is added to decouple the boundary conditions of the equations for Ψ(t) and ξ(t) as indicated in Eq. 17. Taking the variational condition δJ[F] 0, the following coupled equations are obtained, where Here, Ψ(t) satisfies the initial condition Ψ(0) Φ 0 , and ξ(T) satisfies the condition, ξ(T) O T Ψ(T) at final time T. Note that by solving Eq. 17, we obtain the optimal solution Ψ(T), which is equal to Ψ T .

Setup of Target Operators
Consider the ring current localization to a designated ring χ in a PAH, as shown in Figure 6. From Eqs. 10a, 12, the ring current on ring κ at the target time T is expressed as (18b) Here, n is the number of electronic excited states. Hereafter, we write c α (T) as c α for simplicity.
The target state Ψ T α c α Φ α can be determined by applying the Lagrange multiplier method to the ring currents at the target time T in Eq. 18a, which provides the coupled equations in terms of the configuration interaction coefficients with {Rec α , Imc α }. The target operator is given as O T Ψ T > < Ψ T . The coupled equations for the ring current localization to aromatic ring χ can be expressed as A brief derivation of Eq. 19 is summarized in Appendix A.
The coefficients {Rec α , Imc α }with λ l included are numerically determined by applying the Newton-Raphson method to the coupled equations in Eq. 19. The target operator for the localization is obtained as O T α c α Φ α 〉〈 β c β Φ β . Now consider the bridge bond currents J Bl ′ shown in Figure 6, and the perimeter ring current of a PAH with m aromatic rings, J l , which is defined as the average of the π-electron ring currents at each aromatic ring site. For the perimeter ring current, the bridge bond currents J Bl' (1 ≤ l ′ ≤ m − 1) flowing among the nearest neighbor aromatic rings, should be zero at the target time. The target operators for the bridge bond currents, J Bl' ( 1 ≤ l ′ ≤ m − 1), and perimeter ring current, J P , can be derived from Eq. 19 by replacing J κ,βα with J Bl',βα and J P,βα , respectively, because J Bl' and J P can be written in terms of {Rec α , Imc α } as and respectively. Similar to Eq. 19, the coupled equations for the perimeter ring current are derived in Appendix B.

Application of Current Localization Control to Anthracene
We applied the quantum optimal control (QOC) method described in the preceding section to anthracene, one of the smallest PAHs [42]. Anthracene (D 2h ) was assumed to be fixed on a surface (the XY plane), or oriented spatially by lasers The molecular geometry was optimized in the MP2/6-311+g(d,p) level theory using the GAUSSIAN09 code [47]. As demonstrated in Figure 7A, anthracene consists of three aromatic rings, which are called L-, M-, and R-rings respectively.
Consider the two types of π-electron localized ring currents of anthracene: the perimeter current flowing along the outside chemical bonds of anthracene, and the middle ring current localized to the M-ring. Both types of localized ring currents belong to the irreducible representation B 1g of the D 2h point group, which are symmetry-adapted. The excited states contributing to the two ring currents need to belong to the B 3u and B 2u representations, because each corresponding electronic coherence created by the two excited states belongs to the B 1g representation. The excited state for the localized ring  Frontiers in Physics | www.frontiersin.org July 2021 | Volume 9 | Article 675134 current is S 3 with B 3u representation, and the other two excited states S 5 and S 6 with B 2u representation as shown in Figure 7B. This defines a "symmetry-adapted ring current" [41,42]. The excited state energies in Figure 7B were calculated to be E 3 5.10 eV, E 5 5.71 eV, and E 6 5.78 eV. Here, S 4 (B 3u ) with E 4 5.13 eV was excluded because the oscillator strength was negligibly small (f 0.0013) in our numerical simulation. The nonzero transition dipole moments relevant to the coherent π-electron ring current control were calculated using the time-dependent density functional theory (TDDFT) method as, μ S0,S3 (−3.98, 0, 0) , μ S0,S5 (0, −0.64, 0), and μ S0,S6 (0, −0.44, 0) (in a.u.) for the ground and excited states. The vapor absorption spectrum of anthracene shows the strongest absorption peak at 5.30 eV [62], which corresponds to E 3 5.10 eV in the TDDFT calculation. Figure 8 shows the QOC results for the generation of an anticlockwise perimeter ring current in anthracene. The target state is expressed from the results obtained using the Lagrange multiplier method as Here, the subscripts specifying the target state, for example, ACA for Ψ ACA , indicates the anti-clockwise ring currents in two aromatic rings, L and R, and a clockwise ring current in aromatic ring M. It should be noted that the target state for a generation of any coherent ring current is expressed in terms of its complex form. The control target time was set to T 60 fs. The matrix elements of π-electron ring currents J Χ,αβ (X L, M, R) and the bond current J Bi,αβ (i 1, 2), for two excited states α, β are presented in Table 2. Here the π-electron ring current flowing in an anticlockwise direction is defined as positive, whereas the bond current flowing toward the Y direction is defined as positive (See Figure 6). Figure 8A shows the temporal evolutions of the coherent ring currents J χ for the three aromatic rings of anthracene, χ L, M, R, and those for B1 and B2. It can be observed from Figure 8A that the two bond currents J B1,αβ and J B2,αβ vanish at the target control time of T 60 fs, and all ring currents in three aromatic rings exhibit positive. This indicates that the coherent π-electrons rotate clockwise along the outside (perimeter) of the aromatic ring. The laser-controlled π-electrons includes a ring current averaged over the three aromatic rings, J p ≡ (J L + J M + J R )/3 89.0 μA, which is the perimeter ring current. Figure 8B shows the temporal behavior of the ground state S 0 , and three excited states S 3 (B 3u ), S 5 (B 2u ), and S 6 (B 2u ) populations during the QOC process, which are induced by the two control lasers. Figure 8C shows the X and Y-components of the electric field F(t) generated by the control lasers, and Figure 8D shows the Fourier transformed spectra F X (ω) and F Y (ω) in the ultraviolet (UV) frequency domain. By analyzing the temporal behavior of the electric field F(t) of the control lasers in Figure 8C and the Fourier transformed spectra of the control laser fields in Figure 8D, the mechanisms of the laser-controlled ring currents in anthracene can be clarified. It is evident that from the modulation in Figure 8C the two linearly (Y-) polarized lasers with a relative phase zero induce the electronic coherence between the two excited states S 5 (B 2u ) and S 6 (B 2u ), and that the linearly polarized laser pulse parallel to the X-axis creates the perimeter ring current on the molecular plane.
Thus far, we have considered the QOC procedure for generating an anti-clockwise perimeter current. The QOC procedure for generating a clockwise perimeter current can be carried out in the same manner as described above (See Eq. 21) by considering the target state, We carry out the QOC procedure to generate the anticlockwise ring current localized to the middle ring of anthracene. The target state is expressed as The target state for the middle ring current with the clockwise flow is given by the relationship Ψ 0C0 Ψ p 0A0 . Figure 9 exhibits the QOC results, where Figure 9A displays the temporal evolutions of ring current localization control to the three aromatic rings, indicating that J M is a ring current with 64.4 μA at the target control time of 60 fs, but on the other hand the ring currents J L and J R of the other two aromatic rings vanish. This indicates that the ring current localized to the middle aromatic ring is created by the control laser pulses presented in Figures 9C,D. The temporal behaviors of the populations shown in Figure 9B are nearly the same as those of the perimeter ring currents as shown in Figure 8B. As presented in Figure 9C, the two linearly (Y-) polarized lasers induce the electronic coherence between the two excited states S 5 and S 6 with a definite relative phase π, in contrast to the results presented in Figure 8C. In the mechanism of generation between the perimeter and the middle ring currents there is no difference except the phase difference between the two excited states S 5 and S 6, because the temporal behavior in both Figure 8B and Figure 9B and the temporal behavior in the X-polarized laser pulse in both Figure 8C and Figure 9C are similar to each other.
Main difference between the temporal behaviors of the Y-polarized lasers in Figure 8C and Figure 9C can be explained briefly in the following. The two coherent states created by the Y-polarized lasers can be expressed as and For simplicity, the amplitudes of the two coherent states are assumed to have the same value, i.e., F S5 F S6 ≡ F Y . The above expressions can then be simplified to and Here, the frequency difference |ω S5 − ω S6 |/2 yields an oscillation beating period of 120 fs, which is nearly the same as that of the quantum beat frequency observed in Figure 8C and Figure 9C.
Having clarified the control mechanism of the coherent ring current generation, we can semi-quantitatively reproduce the above results in Figures 8, 9 by using an analytical method [42]. This indicates, in principle, that two types of ring current localizations, the perimeter ring current and the middle ring current in linear PAHs, can be generated in experiment using three coherent UV lasers without a sophisticated QOC device.
For QOC numerical simulations of the ring current localization control in anthracene, we have only considered the perimeter and the middle ring currents generations, which are symmetry-adapted, while we did not consider a ring current localization to the L-or R-ring of anthracene, which are created by a symmetry-broken procedure, as demonstrated for naphthalene in Ref. [41]. That is, other excited state(s) with gerade symmetry must be considered in addition to the excited FIGURE 9 | Quantum optimal control simulations for generation of the middle ring current in anthracene: (A) temporal behavior of the ring currents J L (t) J R (t) (solid line) and J M (t) (broken line); (B) temporal behavior of the S 0 population and those of the three excited state S 3 (B 2u ), S 5 (B 2u ) and S 6 (B 2u ) populations; (C) optimized X-and Y-polarized laser pulse fields F X (t) (left-hand side) and F Y (t) (right-hand side); (D) Fourier transformed spectra of the laser pulse fields F X (ω) and F Y (ω). Reprinted with permission from Ref. [42] Copyright (2017) American Institute of Physics.
Frontiers in Physics | www.frontiersin.org July 2021 | Volume 9 | Article 675134 states with ungerade symmetry as described in the previous subsection.
It is important to consider how to observe the ultrafast coherent ring currents in PAHs. We proposed a method to detect the direction of the atto-second coherent ring current by tracking the femtosecond molecular vibrational motions that can induce the nonadiabatic couplings [30]. Rodriguez and Mukamel [63] have proposed measuring the circular dichroism (CD) of the ring current using the pump-probe method. Recently, two methods have been proposed for the detection of the ultrafast coherent ring currents. One method, recently proposed by Yuan et al., Bandrauk's group [64], is utilizing the atto-second detection method of the coherent electronic dynamics in molecules with the temporal and spatial resolutions using the circularly polarized ultrashort UV pump and X-ray probe laser pulses. The other method is to utilize the time-resolved scanning microscopy (STM) and the magnetic force microscopy (MFM) to detect the ring current-induced magnetic fields during the ultrashort time [65][66][67].
It is also expected that to apply our methods to planar or nonplanar extended molecular systems, such as graphene sheets. In principle, our method can be extended to these large systems, by considering the symmetry of the molecular system and the constraints on the π-electron ring currents.
Using laser control, it is also essential to maintain the created ring current at the target region (path) at least during one vibrational period of the PAH [68]. Otherwise, the created ring current would dissipate quickly because of the vibronic interactions and/or the nonadiabatic couplings between the two electronic excited states [30,[69][70][71].

DYNAMIC STARK-INDUCED π-ELECTRON ROTATIONS IN LOW-SYMMETRY AROMATIC RING MOLECULES
In this section, we present a convenient method for inducing unidirectional π-electron rotations in aromatic ring molecules with low symmetry [43]. The basic idea behind the induction of unidirectional electron rotations is to degenerate two nondegenerate excited states by utilizing dynamic Stark shifts, as demonstrated in Figure 10. We refer to this as the dynamic Stark-induced degenerate electronic state (DSIDES) [43]. Two linearly polarized continuous lasers operating at different frequencies and phases are used to form DSIDES: Each laser is set to selectively interact with each electronic state through non-resonant excitation. The lower and higher excited states are shifted up and down with the same population, respectively, and the DSIDES is formed at the center between them. As a result, unidirectional π-electron rotation is driven by two lasers. In the laser control scenario, only one input parameter out of the four possible parameters (frequency and intensity for each laser), is required to induce the DSIDES formation.
First, the DSIDES formation is described in a three-electronic state model under the fixed nuclei condition, and the timedependent expectation value of the angular momentum operator of π-electrons is derived and analytically expressed. Next, to demonstrate the applicability of the control scenario, the results of the DSIDES treatment for toluene, which is a typical aromatic ring molecule of low symmetry belonging to Cs point group, are presented.

Formation of Dynamic Stark-Induced Degenerate Electronic State DSIDES and the Resulting Angular Momentum
Consider the coherent π-electron angular momentum in an aromatic ring molecule with low symmetry excited by nonresonant, stationary UV/visible lasers. The molecule of our interest is one oriented in a space or attached to a surface, as mentioned in the preceding sections. We adopt a three-electronic state model in the frozen nuclei approximation. The three electronic states including the ground state specified by the energy ε 0 (≡ Zω 0 ), and two excited states specified by ε 1 (≡ Zω 1 ) and ε 2 (≡ Zω 2 ), as shown in Figure 10. Here, DSIDES can be formed by two stationary linearly polarized lasers with detuning frequencies Δ 1 ω a − ω 10 < 0 and Δ 2 ω b − ω 20 > 0. The frequency difference between the two excited states is expressed as ω 21 ω 2 − ω 1 . The dynamic Stark shifts are denoted by the Rabi frequencies Ω 1 and Ω 2 , as presented in Figure 10. The wave function of the total system in the stationary lasers is defined through laser-molecular interactions as Here, the normalization condition for the coefficients c 0 (t), c 1 (t), and c 2 (t) are as follows, FIGURE 10 | Dynamic Stark-induced degenerate state with the Rabi frequencies (Ω 1 and Ω 2 ) using two stationary linearly-polarized lasers with the frequencies, ω a and ω b , and detuning frequencies Δ 1 andΔ 2 respectively. Frontiers in Physics | www.frontiersin.org July 2021 | Volume 9 | Article 675134 The system Hamiltonian interacting with the electric fields is given as where H 0 satisfies H 0 ϕ i ε i ϕ i for i 0, 1, and 2. In Eq. 29, the interaction Hamiltonian V(t) between the system and the two electric fields is written as Here μ −er is the electronic dipole moment operator, where r means the electron coordinate. F α F α e α (α a, b) is the electric field with amplitude F α and photo-polarization vector e α , ζ α is the initial phase, and ω α is the central frequency. In Eq. 30, the two electric fields denoted by a and b induce non-resonant transitions from the ground state to the two excited states. The selective transition conditions are set by the choice of the laser polarization vectors (e a and e b ) satisfying μ 02 ⊥e a and μ 01 ⊥e b , respectively.
The time-dependent Schrödinger equation can be written as The coefficients must satisfy the coupled differential equation where the interaction Hamiltonian, where the transition magnitudes are V a 01 ≡ − μ 01 · F a /(2Z), and In Eq. 34c, the two dressed states are assumed to have equal energies, i.e., Ω ≡ Ω 1 Ω 2 ( ), which is called the Rabi frequency [43]. The three conditions lead to a reduction of the input parameters of the two lasers, such that the two amplitudes (F a and F b ), and two central frequencies (ω a and ω b ) are reduced to one input parameter. We take F a ≡ F hereafter. Analytical expressions for time-dependent coefficients {c i (t)} are given in Appendix C.
The time-dependent angular momentum defined as an expectation value of an angular momentum operator L Z −iZ z zφ Z l Z , l Z −i z zφ , can be expressed, Here, l Z,12 < ϕ 1 l Z ϕ 2 > − < ϕ 2 l Z ϕ 1 > −l Z,21 and l Z,12 iIml Z,12 .

Unidirectional π-Electron Rotations in Toluene
We calculate L Z (t) (derived in the preceding subsection) in a real molecule, toluene. The simplest three-electronic state model is applied for toluene because the quasi-degenerate states S 3 (A") and S 4 (A′) in toluene (C s ) correspond to the doubly degenerate state S 3 (E 1u ) in benzene (D 6h ): Note that the S 3 (E 1u ) state is a dipole-allowed excited state in benzene (D 6h ), whereas the other two lower excited states, S 2 (B 1u ) and S 1 (B 2u ), are dipoleforbidden. The geometry optimization of toluene was carried out with the MP2 level of theory [43]. The C S symmetry of toluene indicates that one of the hydrogen atoms belonging to the methyl group is perpendicular to an aromatic ring plane ( Figure 11A). Figure 12A shows the calculated time-dependent angular momentum expectation values L Z (t) with respect to the relative phase ζ −π/2 (ζ +π/2) between the two lasers for the left (right) panel. These were calculated using Eq. 35 combined with Eq. 34. Here, the amplitude of the electric field F a is adopted as the input parameter F (≡ F a ). Other parameters are shown in Table 3A, while Table 3B shows the timedependent populations in the three electronic states.
The time-dependent angular momenta plotted in Figure 12A are comprised of angular momentum pulse trains of the same shape for each value of F. Each angular momentum pulse corresponds to the unidirectional π-electron rotation, which begins with acceleration and ends with deceleration. The direction of the π-electron rotations is determined by the relative phase ζ between the two lasers: clockwise rotation for ζ +π/2, and anti-clockwise rotation for ζ −π/2. Here we discuss how the unidirectional π-electron rotations are created. We first note that L Z (t) can be rewritten under the two conditions, Eq. 34a and Eq. 34b, for Ω≈Δ as expresses no unidirectional π-electron rotations. By introducing the third condition, Eq. 34c, which provides for the formation of the doubly degenerate dressed states with ζ ±π/2, Eq. 36 can be expressed as This indicates a unidirectional π-electron rotation whose direction is determined by the relative phase ζ, that is, L Z (t) > 0 for ζ −π/2 and L Z (t) < 0 for ζ +π/2, as demonstrated in Figure 12. The dotted line shown in Figure 12A represents L Z (t) for F 1.5 GV/m as calculated by Eq. 37, which well reproduces L Z (t) curve calculated using Eq. 35 without the approximation of Ω≈Δ approximation. Thus, the unidirectional π-electron rotation originates from the DSIDES formed by the two non-resonant lasers with a definite relative phase of ±π/2.
As demonstrated in Figure 12A, the oscillation periods of the angular momentum pulses become shorter (higher Rabi frequency Ω), and those amplitudes decrease as F increases. These two behaviors result from the third restriction condition for degeneracy (Eq. 35c) used in the derivation of L Z (t). This degeneracy condition should be satisfied for the two dressed states with equal energies to maintain their energies located at the center of the two excited states, even though the laser intensities of two lasers increase. This results in an increase in the detuning parameter Δ of the two lasers, that is, a decrease (increase) in ω a (ω b ), as shown in Table 3A. Figure 12A exhibits that the maximum magnitude of the angular momentum occurs in the vicinity of F 1.5 GV/m. This can be explained by noting that the constant in Eq. 37, V 2 Δ 2 /Ω 4 , with Rabi frequency of Ω Δ 2 + 2V 2 √ , reaches its maximum at Δ 2 2V 2 , which in turn provides the maximum angular momentum. Using the relationship Δ 2 2V 2 for the dressed states with equal energies in Eq. 34c, the electric field F is estimated as F 2 √ ZΔ/(μ 10 ) X 1.52 GV/m with Δ 0.118 eV, which is close to the F 1.5 GV/m value shown in Figure 12A. Figure 12B shows the F-dependence of the c S 0 (t) 2 , c S 3 (t) 2 and c S 4 (t) 2 populations in toluene, which were calculated using Eqs. (B2). The relationship between the excited state populations, c S 3 (t) 2 c S 4 (t) 2 , is maintained because the conditions represented by Eqs. 34a, 34b, were used for a derivation of L Z (t). Inversion of populations between the ground and excited states occurs for F 1.5 GV/m and 2.0 GV/m, as demonstrated in Figure 12B. By comparing Figures 12A,B, it can be observed that major components of the unidirectional angular momentum pulse are created during when the population inversion occurs.
We now estimate the physical constants associated with π-electron rotations, which are listed in Table 3B: one cycle count of π-electron rotation N r , and the magnitude of ring current J in one cycle of Rabi oscillation with periodicity T ≡ 2π/Ω. In particular for aromatic molecules the ring current is an important physical quantity because it is directly related to the magnetic field induced by the ring current. For simplicity, we consider the classical model of π-electron rotation, in which N r is defined as N r ≡ T/τ. Here τ 2πm e r/p 2πm e r 2 / L Z is the rotational period of a π-electron having an angular momentum averaged over one cycle period of Rabi oscillation (L Z ). m e is the electron mass, and r indicates the radius of the aromatic ring. L was calculated using Eq. 37 substituting for L Z (t) as The ring current can be approximately expressed as J ≈ L Z J 0 /Z [39][40][41]43]. Here on a single aromatic ring, J 0 200 μA was a typical value of the ring current estimated with a unit angular momentum h. Table 3B exhibits the calculated ring current J, cycle counts of π-electron rotations along the aromatic ring N r and three values of F adopted in Figure 12, together with the parameters used in this model calculation. In section 3, we presented a pump-dump control procedure for the generation of sequential unidirectional ring currents in a 2,2'-biphenol  a Parameters V, Δ, Ω, ω a , and ω b , were calculated for three values of F using Eq. 34. These parameters used for the calculation of the following physical quantities: T as a period of Rabi oscillation; LZ is the angular momentum averaged over one cycle period of Rabi oscillation, τ 2πmer/p 2πmer 2 / LZ as the rotational period of a π-electron, J as a ring current, and N r ≡T/τ as one cycle count of π-electron rotations.
Frontiers in Physics | www.frontiersin.org July 2021 | Volume 9 | Article 675134 molecule having non-degenerate excited states. The same order of magnitudes of the pump-and dump-pulse lasers were applied in the vicinity of ∼1 GV/m, considering that the number of aromatic rings differs between biphenol and toluene. This can be explained by the point that the magnitudes of the ring current J or the angular momentum L Z are proportional to that of the coherence between the two excited states, but not proportional to the number of aromatic rings. It is remarkable that in Table 3B J is not proportional to N r , it is actually proportional to N r /T. This is originated from the different behaviors between J and N r related to the F-dependence shown in Table 3B.

SUMMARY AND PERSPECTIVE
This paper briefly provides an overview of the theoretical study of quantum laser control of coherent π-electron dynamics in low-symmetry aromatic ring molecules, which we have undertaken in recent several years [39][40][41][42][43][68][69][70][71][72][73]. The essential principles to generate the π-electron angular momentum and ring current in a low symmetry aromatic ring molecule are first to create a superposition of two electronic excited (coherent) states using two linearly polarized lasers, and second to select the clockwise or anticlockwise rotational component from the non-stationary time evolution of the coherent state using pump and dump lasers with properly designed polarization directions. Here, the direct product of the irreducible representations of each excited states in the molecular point group belongs to that of the molecular rotation, (Rx, Ry or Rz), and the relative phase between the two linearly polarized lasers are fixed in order to determine the initial direction of the angular momentum and ring current. In additon to the laser conditions, aromatic ring molecules also need to satisfy the follwing conditions to generate the π-electron angular momentum and ring current. First two electronic excited states are nearly degenerate, and second, those states can be excited from the ground states using UV lasers, i.e., transitions between the electronic excited states and ground states are optically allowed. On the basis of the above principles, analytical expressions for coherent π-electron angular momentum, ring current, and ring current-induced magnetic field are derived in the quantum chemical MO theory. A coherent π-electron ring current is defined as the average of the bond currents. Temporal behaviors of these coherent quantities are calculated using the density matrix method within the Markov approximation. Quantum simulations of the coherent quantuities have been performed for (P)-2,2'-biphenol. Based on the simulation results a new quantum control method for two-dimensional ultrafast angular momentum switching was proposed. The essential point of the quantum control is to use a sequence of the overlapped pump and dump pulses with a selected relative phase and selected polarization directions between the two lasers. The results for the (P)-2,2'-biphenol indicate that this new control scheme can potentially be used for the design and realization of ultrafast multi-dimensional electronic switching devices, or fieldeffect transistor devices.
The familiar quantum optimal control procedure was applied to an optiaml control of the coherent π-electron ring currents in planar PAHs with D 2h symmetry. For this purpose, we have to set up the target state for a desired ring current pattern. The final target wave functions are determined via the Lagrange multiplier method by solving the coupled equations of motion under the constraints that the ring currents must satisfy. The method was applied to anthracene. The creation of the perimeter current and the middle ring current in anthracene were successfully realized. The control mechanisms were elucidated by analyzing the timedependent behaviors of the control laser puses and the populations of the relevant electronic states. Concerning the types of molecules used in the above laser control studies, polycyclic planar aromatic molecules like PAHs, coronene and benzoic acids dimers are applicable for a laser control of perimeter or localized ring currents.
Another laser control method for the π-electrons unidirectional rotation in low symmetric aromatic rings was presented. The basic idea of the control method is to degenerate two nondegenerate excited states by using dynamic Stark shifts. Doubly degenerate excited state was created by the Stark shift using the two linearly polarized non-resonant lasers with a definite relative phase of ± π 2 , which determines the angular momentum direction (clockwise or anti-clockwise π-electron rotation). Applying the RWA, the coherent π-electron angular momentum in a three-electronic state model was analytically expressed in a closed form, and the proposed control theory have been applied to toluene (Cs). The numerical simulations indicated that the resulting angular momentum comprised sequential angular momentum pulses with a positive or negative value depending on the relative phases between the two lasers.
Some theoretical studies on laser control of the coherent π-electron rotations that were not included in this review are noteworthy. Probing molecular chirality, right-handed or lefthanded chiral molecules, is a central issue in natural science. E.x., Phenylalanine is an essential amino acid, and L-enantiomer is found in natural system. On the other hand, D-enantiomer is synthesized artificially, and racemate phenylalanine is used for dietary supplement. However, in conventional methods such as CD spectroscopy and optical rotatory dispersion (ORD) spectroscopy [74,75], the signals are very week because these processes involve second-order evaluations. As a possible of laser control scheme of the unidirectional π-electron rotations in low symmetric aromatic ring molecules, we previously proposed an efficient enantiomer-probing scenario for chiral aromatic molecule [71,72]. A pair of nondegenerate excited states becomes degenerate by applying the dynamic Stark shift in the presence of two non-resonant UV lasers, producing an enantiomer-specific angular momentum. In that study, phenylalanine was adopted as an example for the numerical simulation. The resulting enantiomer-specific magnetic fluxes were on the order of a few Tesla, with periods of several tens of femtoseconds.
The contribution of vibrational motion to coherent π-electron rotations is also an important research topic that was not discussed in this review. In our previous research, vibrational effects on the coherent π-electron rotations in (P)-2,2'-biphenol were theoretically studied in the adiabatic approximation [73]. It was found that the low-frequency torsion mode around the bridge causes modulations in the beat of the ring current [73]. The vibrational effects on dynamic Stark-induced π-electron rotations in aromatic ring molecules with low symmetry were also studied using the displaced harmonic oscillator (DHO) model in the adiabatic approximation [72]. A pair of the lowest vibronic state in the two electronic excited states was degenerated using the two linearly polarized UV lasers. The two potential displacements between the ground state and the two electronic excited states were used as the parameters. The shapes of the sequential angular momentum pulses were affected by the potential displacements.
From one perspective, theoretical treatments beyond the adiabatic approximation [28,31,[74][75][76][77] are expected to be incorporated such that other electronic excited states interact with the two relevant excited states through nonadiabatic couplings. It is crucial to clarify the effect of nonadiabatic couplings to the unidirectional π-electron rotations and determine how to maintain π-electron rotations by way of quantum optimal control techniques. In large PAHs, once π-electron rotation is realized at a localized ring site and once a site-selective coherent ring current is generated, the localized ring current is transferred from site to site using the laser pulses, i.e., transferred ring currents, which allow for an ultrafast switching function at the selected local site. Furthermore, the site-selective coherent ring current and transferred ring current can create the induced magnetic fields. These electromagnetic fields are expected to provide ultrafast probing of local sites in large molecular systems, biomolecules with chiral aromatic ring molecules, and PAHs. Further development of theoretical treatments involving photon-induced electronic coherence in molecular systems would be promising in the near future.

AUTHOR CONTRIBUTIONS
The authors made a substantial, direct, and intellectual contribution to the work and approved it for publication. Original draft preparation, software, formal analysis, and funding acquisition, HM; Conceptualization, supervision, project administration, and validation, YF; Writing-original draft preparation, HM and YF; Investigation, writing-review and editing, HM, N-LP, and YF.