We present a quantum control scheme which realizes suppression of the intramolecular vibrational energy redistribution (IVR). In this scheme, we utilize effective decomposition brought by intense CW-laser fields, which enables to exclude the doorway state coupled to background manifolds. In doing so, we introduce a helper state and make it optically coupled with the doorway state through the intense CW-laser field. We have applied the present scheme to both the Bixon-Jortner model and the SCCl2 model system.

1. Introduction

In coherent laser control of molecular vibrations, the key feature is manipulation of vibrational coherence. Thus, dissipative processes such as dephasing and population decay have been considered as formidable obstacles. One of the typical dissipative processes is intramolecular vibrational energy redistribution (IVR), which is characterized as irreversible population flow from the doorway state to the background manifolds [15]. There have been several attempts to suppress the IVR by laser irradiation. For example, the suppression of IVR is theoretically confirmed under a strong resonant CW-laser excitation, in which total population is trapped on the initial and doorway states, and it rapidly oscillates between those two states [6]. It is also reported that the IVR can be restricted by the laser field, in which the frequency is swept [7, 8]. This phenomena is explained by the concept of adiabatic passages and the population is locked on the doorway state during the chirp pulse. In those approaches, only the suppression of IVR is considered, and there is no concept of controlling the population dynamics which escape from the IVR. Thus, it is necessary to apply additional control fields if one aims to drive the system to a specific desired quantum state while avoiding the IVR. However, additional laser would interfere with those existing “population locking” or “IVR-suppression” fields.

Another possible approach towards the suppression of IVR is to design the control field which actively manipulate phase relations between the eigenstates contained in the current wave packet. This could be done by optimal control theory (OCT) or local control theory (LCT) [913]. However, optimized laser fields tend to show complicated and highly sensitive nature since they try to control numerous eigenstates with various eigenenergies, which makes difficult to prepare the control laser in practice.

The main purpose of this study is to propose an alternative control scheme applicable to dissipative systems with which one can overcome the difficulties stated above. The basic idea is very simple. First, we effectively decompose the total system into target space and its complimentary space as we have shown in our former studies [14, 15]. By doing this, we can exclude the doorway state which is responsible for the dissipative dynamics, or IVR, from the system of interest, while remaining indirect optical transitions for control process. Furthermore, choosing a small space consisting of a few states as a target space, one can apply the well-established control schemes such as 𝜋-pulse [16] or stimulated Raman adiabatic passage (STIRAP) [17], which greatly simplifies the control process. The mechanism of exclusion of the doorway state is well investigated for the 3-level system, in which two lasers are applied with static detuning [18]. The decomposition mechanism we adopt in the present work is basically the same as the one in [18] except the fact that the detuning itself is introduced as the induced Stark shift. We discuss the difference of those two approaches and clarify the advantageous features of the present scheme from a viewpoint of applicability to quantum control.

As a model system, we employ the classic IVR model system, Bixon-Jortner(BJ) model [19], to confirm the effectiveness of the present method and show how the IVR process can be excluded from the molecular system dynamics. Next, as an example of more realistic system, we choose SCCl2 molecule, whose IVR is well studied, and there exists reliable model Hamiltonian [20, 21].

2. Theoretical

In the present work, we consider a model system as shown in Figure 1. Before proceeding to the formulation, we explain several terms used. Here, |𝑖 is the initial state, which is typically taken to be the lowest vibrational state on the ground electronic state. The state |𝑚 is the intermediate state, which optically connects the two states |𝑖 and |𝑓. The state |𝑓 is the control target state, that is, we aim to transfer the initial population to |𝑓 avoiding the population loss due to IVR process. Note that |𝑚 is optically accessible from |𝑖, |𝑓, and | and also coupled with {|𝐵𝑗}(𝑗=1,2,,𝑀) through the intramolecular interactions {𝑣𝑗}. The IVR, that is, the irreversible population flow to {|𝐵𝑗}, occurs through |𝑚, thus, we call this state “doorway state”, hereafter. In addition to that, we introduce the special state | which we call “helper state” in order to suppress the IVR process. This state is coupled through the intense CW-laser resonant to the transition between | and |𝑚 with interaction strength Ω.

Next, we formally split the total system into three spaces as the follows. Control target space, which we call A space, consists of two states, initial |𝑖 and final state |𝑓. We define B space as the one spanned by a doorway state |𝑚 and a helper state |, which we utilize for realizing effective decomposition. C-space is defined as the space spanned by IVR background manifolds {𝐵𝑗}.

The total Hamiltonian matrix under the rotating wave approximation (RWA) is given as𝐇=𝐇0𝚫+𝜆𝐇𝐀𝟎𝟎𝟎𝐖𝐁𝟎𝟎𝟎𝐅𝐂+𝜆𝟎𝐕𝐀𝐁𝟎𝐕𝐁𝐀𝟎𝐕𝐁𝐂𝟎𝐕𝐂𝐁𝟎.(1) Here, we introduce the parameter 𝜆 so as to clearly specify the perturbation order. The optical interactions between |𝑖 and |𝑚 (𝑉1), and |𝑚 and |𝑓 (𝑉2) are taken to be much smaller than Ω. Minor matrices in (1) are defined as𝚫𝐀=Δ00Δ,𝐖𝐁=,𝐅Ω00Ω𝐂=𝜖10000𝜖200000000𝜖𝑀,(2) where Δ denotes the detuning introduced for AB transition (see Figure 1), while {𝜖𝑖} (𝑖=1,2,,𝑀) denotes the eigenenergies of the background manifolds in C-space. Note that we prediagonalize the B space Hamiltonian so as to make 𝐖𝐁 diagonal matrix. The matrices 𝐕𝐀𝐁 (2×2) and 𝐕𝐁𝐂 (2×𝑀) correspond to the optical and static interactions, respectively.

The Schrödinger equation in the matrix representation is given as𝐇𝐔=𝐔𝐄,(3) where 𝐄 is a diagonal matrix with eigenvalues {𝐸𝑖} (𝑖=1,2,,𝑁) in its diagonal elements, and 𝐔 is the unitary matrix consisting of corresponding eigenvectors 𝑐𝑖 (𝑖=1,2,,𝑁) as 𝐔(𝐜1,𝐜2,,𝐜𝑁). Here, we expand 𝐔 and 𝐄 as𝐔=𝐔(0)+𝜆𝐔(1)+𝜆2𝐔(2)+,𝐄=𝐄(0)+𝜆𝐄(1)+𝜆2𝐄(2)+,(4) where𝐔(𝑛)=𝐔(𝑛)𝐀𝐀𝐔(𝑛)𝐀𝐁𝐔(𝑛)𝐀𝐂𝐔(𝑛)𝐁𝐀𝐔(𝑛)𝐁𝐁𝐔(𝑛)𝐁𝐂𝐔(𝑛)𝐂𝐀𝐔(𝑛)𝐂𝐁𝐔(𝑛)𝐂𝐂,𝐄(𝑛)=𝐄𝐀(𝑛)𝟎𝟎𝟎𝐄𝐁(𝑛)𝟎𝟎𝟎𝐄𝐂(𝑛).(5) Here, 𝐄𝐀(𝑛), 𝐄𝐁(𝑛), and 𝐄𝐂(𝑛) are diagonal matrices, whose diagonal elements are 𝑛th perturbation energies. Inserting (1), (4) into (3) and comparing 𝜆𝑛 (𝑛=0,1,2) terms gives the followings:𝐇0𝐔(0)=𝐔(0)𝐄(0),𝐇(6)0𝐔(1)𝐔(1)𝐄(0)+𝐇𝐔(0)=𝐔(0)𝐄(1)𝐇,(7)0𝐔(2)𝐔(2)𝐄(0)𝐔(1)𝐄(1)+𝐇𝐔(1)=𝐔(0)𝐄(2).(8) From (6), one obtains the following equations for minor matrices:𝚫A𝐔(0)𝐀𝐀=𝐔(0)𝐀𝐀𝐄𝐀(0),𝐖(9)𝐁𝐔(0)𝐁𝐁=𝐔(0)𝐁𝐁𝐄𝐁(0)𝐅,(10)𝐂𝐔(0)𝐂𝐂=𝐔(0)𝐂𝐂𝐄𝐂(0),(11) while off-diagonal blocks of 𝐔(0) are all zero matrices. As Δ𝐀 commutes with 𝐔(0)𝐀𝐀, (9) gives 𝐄𝐀(0)=Δ𝐀, which leaves 𝐔(0)𝐀𝐀 undetermined. On the other hand, from (10), (11) one obtains that 𝐔(0)𝐁𝐁 and 𝐔(0)𝐂𝐂 are both determined as identity matrices of 2×2 and an 𝑀×𝑀, respectively, whereas 𝐄𝐁=𝐖𝐁 and 𝐄𝐂=𝐅𝐂. The first two terms of the left hand side of (7) are given as 𝐇0𝐔(1)𝐔(1)𝐄(0)=𝟎𝚫𝐀𝐔(1)𝐀𝐁𝐔(1)𝐀𝐁𝐖𝐁𝚫𝐀𝐔(1)𝐁𝐂𝐔(1)𝐁𝐂𝐅𝐂𝐖𝐁𝐔(1)𝐁𝐀𝐔(1)𝐁𝐀𝚫𝐀𝐖𝐁𝐔(1)𝐁𝐁𝐔(1)𝐁𝐁𝐖𝐁𝐖𝐁𝐔(1)𝐁𝐂𝐔(1)𝐁𝐂𝐅𝐂𝐅𝐂𝐔(1)𝐂𝐀𝐔(1)𝐂𝐀𝚫𝐀𝐅𝐂𝐔(1)𝐂𝐁𝐔(1)𝐂𝐁𝐖𝐁𝐅𝐂𝐔(1)𝐂𝐂𝐔(1)𝐂𝐂𝐅𝐂,(12)whereas other terms are𝐇𝐔(0)=𝟎𝐕𝐀𝐁𝐔𝐁𝐁𝟎𝐕𝐁𝐀𝐔(0)𝐀𝐀𝟎𝐕𝐁𝐂𝐔𝐂𝐂𝟎𝐕𝐂𝐁𝐔𝐁𝐁𝟎,𝐔(0)𝐀𝐀𝐄(1)𝐀𝐀=𝐔(0)𝐀𝐀𝐄(1)𝐀𝐀𝟎𝟎𝟎𝐔(0)𝐁𝐁𝐄(1)𝐁𝐁𝟎𝟎𝟎𝐔(0)𝐂𝐂𝐄(1)𝐂𝐂.(13) Comparing the component of A space, one obtains the equation 𝐔(0)𝐀𝐀𝐄(1)𝐀𝐀=𝟎, which leads to 𝐄(1)𝐀𝐀=𝟎 since 𝐔(0)𝐀𝐀 is a unitary matrix. Thus, there is no first-order correction onto the A space eigenenergies.

Next, comparing the BA-block of both sides of (7) leads to the equation𝐖𝐁𝐔(1)𝐁𝐀𝐔(1)𝐁𝐀𝚫𝐀+𝐕𝐁𝐀𝐔(0)𝐀𝐀=𝟎.(14) Considering 𝐔(1)𝐁𝐀Δ𝐀=Δ𝐁𝐔(1)𝐁𝐀, (14) can be rewritten as𝐖𝐁𝚫𝐁𝐔(1)𝐁𝐀=𝐕𝐁𝐀𝐔(0)𝐀𝐀,(15) which leads to𝐔(1)𝐁𝐀𝐖=𝐁𝚫𝐁1𝐕𝐁𝐀𝐔(0)𝐀𝐀.(16) Here, Δ𝐁 denotes the identity matrix of B space multiplied by Δ. The prefactor matrix (𝐖𝐁Δ𝐁)1 in (16) is given in a simple form as𝐖𝐁𝚫𝐁1=1001ΩΔΩΔ.(17) Note that 𝐔𝐁𝐀 becomes zero matrix as Ω becomes large or ΩΔ, which leads to the suppression of AB-space mixing.

Comparing CA-block of (7), one obtains𝐅𝐂𝐔(1)𝐂𝐀𝐔(1)𝐂𝐀𝚫𝐀=𝟎,(18) which is rewritten as𝐅𝐂𝚫𝐂𝐔(1)𝐂𝐀=𝟎,(19) where Δ𝐂 is 𝑀×𝑀 identity matrix multiplied by Δ. It is readily seen that 𝐔(1)𝐂𝐀=0 since 𝐅𝐂Δ𝐂𝟎. Thus, there is no first-order CA-mixing.

Next, we consider the second-order terms. With 𝐄(1)=𝟎, (8) becomes𝐇0𝐔(2)𝐔(2)𝐄(0)+𝐇𝐔(1)=𝐔(0)𝐄(2).(20) Since we focus on the dynamics of A space, we extract AA-block of (20) as𝚫𝐀𝐔(2)𝐀𝐀𝐔(2)𝐀𝐀𝚫𝐀+𝐕𝐀𝐁𝐔(1)𝐁𝐀=𝐔(0)𝐀𝐀𝐄𝐀(2).(21) Considering that the first two terms of (21) are canceled out and inserting (16), one obtains𝐕𝐀𝐁𝐖𝐁𝚫𝐁1𝐕𝐀𝐁𝐔(0)𝐀𝐀=𝐔(0)𝐀𝐀𝐄𝐀(2).(22) Note that (22) is an eigenvalue equation, which determines the second-order energy correction 𝐄𝐀(2) onto A space as well as corresponding eigenvectors 𝐔(0)𝐀𝐀. Here, we define the A space effective Hamiltonian 𝐇(e)𝐀 as𝐇(e)𝐀𝚫𝐀+𝐕𝐀𝐁𝐖𝐁𝚫𝐁1𝐕𝐀𝐁,(23) which determines the dynamics of A space within the second-order. Note that there is no interaction terms relevant to the C-space. This implies that as long as 𝐇 is relatively small enough and the second-order perturbation theory is applicable, the interactions between the doorway state and the background manifolds can be excluded from the system dynamics effectively.

Note that the perturbation expansion shown above is applied for the weak field 𝑉𝑖, which is responsible for interspace optical transitions. On the other hand, the strong interaction through CW-laser, which cannot be dealt with perturbation method, is precisely taken into account by prediagonalization of the B space.

Finally, we discuss possible side effect brought by optically allowed state |𝑚 adjacent to the intermediate state |𝑚. Typically, the level spacing Δ between |𝑚 and |𝑚 is about 1000 cm1. The CW-laser (Ω) applied for decomposition turns out to be off-resonant by Δ, which is large enough to exclude the undesired transition to |𝑚.

3. Results

To confirm the applicability of the present theory, we first apply it to well-known Bixon-Jortner model system [19]. In this model, all the background levels are equally spaced and coupled to a single bright doorway state |𝑚 with uniform strength 𝑣. Shown in Figure 2(a) is the population dynamics under field-free condition while the initial state is taken to be |𝑚. The coupling parameter 𝑣 is taken to be 0.1Ω where as Ω denotes the unit frequency, which determines the timescale (Ω1) of the dynamics. It is shown that the irreversible population flow to the background manifolds {|𝐵𝑗} occurs due to the dephasing process between numerous eigenstates contained in the initial state |𝑚.

Next, we try a naive population control in which control target is the population transfer from |𝑖 to |𝑓. We introduce two resonant CW-lasers corresponding to the transitions, |𝑖|𝑚, and |𝑚|𝑓, respectively. The total Hamiltonian corresponding to (1) is given as 𝐇=Δ0𝑉1000Δ𝑉2𝑉001𝑉20𝑣𝑣𝑣00𝑣𝜖10𝑣𝜖𝑖00𝑣0𝜖𝑀,(24) where the vertical and horizontal dots are all zeroes. Shown in Figure 2(b) is the population dynamics with the initial condition |𝑖. Both optical interaction parameters 𝑉1 and 𝑉2 are taken to be 1.0Ω. As shown in Figure 2(b), significant amount of total population leaks to the background manifolds, and its nature is basically irreversible though coherent oscillations are observed in early-time stage. Note that such oscillatory nature of the population dynamics is rapidly damped, and the system finally settles into stationary state in which the population of |𝑓 is only 25%. Thus, the control objective is hardly achieved with this approach.

Now, we introduce a helper state | in order to suppress the population loss due to IVR. The total Hamiltonian represented by the original basis set is given as 𝐇=Δ0𝑉10000Δ𝑉2𝑉0001𝑉20Ω𝑣𝑣𝑣00Ω00000𝑣0𝜖10𝑣𝜖𝑖00𝑣00𝜖𝑀,(25) where Ω is taken to be the optical interaction between |𝑚 and | in this case. We reconstruct the basis set for B space so as to express the Hamiltonian matrix in the form of (2), that is, 𝐇=Δ0𝑉𝑉000Δ𝑉𝑉𝑉00𝑉Ω0𝑣𝑣𝑣𝑉𝑉0Ω𝑣𝑣𝑣00𝑣𝑣𝜖10𝑣𝑣𝜖𝑖00𝑣𝑣0𝜖𝑀,(26) where 𝑉=𝑉1/2, 𝑉=𝑉2/2, and 𝑣=𝑣/2. Shown in Figure 2(c) is the population dynamics under the condition 𝑉1=𝑉2=0.01Ω. It is shown that the population oscillates between |𝑖 and |𝑓, and there is very little population leak into the background manifolds {|𝐵𝑗}. This feature denotes that the control target space, that is, A space, is effectively isolated as a 2-level system without dissipative dynamics. The effective Hamiltonian for the A space can be obtained from (23) as𝐇(e)𝐀=+ΔΔ00ΔΔ2Ω2𝑉21𝑉1𝑉2𝑉2𝑉1𝑉22.(27) It is clearly shown that the Rabi oscillation is driven by the second term which originates from the second-order interspace optical interaction between A and B spaces. Note also that this term is proportional to the detuning parameter Δ. This implies that common detuning is required for driving population dynamics in the A space.

As for a more realistic molecular system, we consider Thiophosgen (SCCl2), which has six vibrational modes, and there exists IVR process among highly excited states. It is known that only a limited set of vibrational resonances are related to the dissipative dynamics, and an effective Hamiltonian consisting of six modes for describing the IVR process is proposed as [20]𝐻𝑀=𝐻𝑀(0)+𝐻𝑀,(28) where𝐻𝑀(0)=6𝑖=1𝜔𝑖𝑛𝑖+12+6𝑖=1𝜒𝑖𝑖𝑛𝑖+122,𝐻𝑀=𝑉𝑎3Δ𝑛.(29) Here Δ𝑛=𝑖|𝑛𝑖𝑛𝑖|. We used the same parameters listed in [20]. As for the basis set, we take the direct product of the harmonic oscillator eigenfunctions with respect to each mode, which is denoted as |𝑛1,𝑛2,𝑛3,𝑛4,𝑛5,𝑛6. With RWA, the structure of the Hamiltonian matrix is essentially identical to (25). We consider the mode 1 (corresponding to CS-stretching motion) is optically allowed, and the bright doorway state is defined as |8,0,0,0,0,0 as shown in Figure 3. The initial state |𝑖 is taken to be the ground state |0,0,0,0,0,0, and we aim to transfer the population onto the target final state |𝑓|12,0,0,0,0,0, which is one of the highly excited zero-order states of the CS-stretching mode with relatively long lifetime. The laser parameters are taken to be 𝑉1=𝑉2=30cm−1 and Δ=30cm−1. The corresponding laser amplitudes for 𝑉1 and 𝑉2 are estimated as 9.6×109W/cm2 and 5.0×109W/cm2, respectively [20].

Shown in Figure 4 are the population dynamics changing the CW-laser intensity Ω from 0 cm−1 to 300 cm−1. As shown in Figure 4(a), the irreversible population flow to the background modes 2–6, that is, IVR, occurs as we see in the Bixon-Jortner model case. As the laser power increases, the dissipative dynamics is gradually suppressed (see Figures 4(b), 4(c), and 4(d)) due to the decomposition effect, which isolates A space (|𝑖 and |𝑓). Finally, as shown in Figure 4(e), a neat Rabi oscillation appears in the population dynamics of A space, which implies that the target subspace can be effectively treated as a simple 2-level system. Then, one can apply conventional 𝜋-pulse scheme for population transfer onto |𝑓 avoiding the population loss due to the IVR. Note also that the timescale of the Rabi oscillation in Figure 4(e) can be controlled by changing the detuning parameter Δ (see (27)).

Here, we should clearly state the difference in the decomposition conditions between the present method and the static detuning approach in [18]. In the present scheme, the detuning itself is introduced as the induced Stark shift, or a large energy splitting between the zero-order eigenvalues in B space, that is, ±Ω in (2). Thus, the required condition for effective decomposition is different accordingly, that is, Δ𝑉1,𝑉2 in the static detuning case whereas ΩΔ, Ω𝑉1,𝑉2, in the present method.

One of the important differences between those two schemes is the suppression efficiency with respect to the detuning value, Δ or Ω. To see this, we consider a simple 4-level system which consists of only A space (|𝑖 and |𝑓) and B space (|𝑚 and |). Instead of omitting C-space (background states), we introduce a phenomenological dephasing parameter 𝛾 associated with the intermediate state |𝑚, that is, we introduce an imaginary component 𝛾 for the eigenvalue of |𝑚. Here, we take 𝛾=30 cm−1 so as to reproduce the IVR timescale of SCCl2. Shown in Figure 5(a) is the population dynamics with Ω=0 cm−1 and 𝑉1=𝑉2=30cm−1, which corresponds to the condition of Figure 4(a). One can see that the dissipative IVR dynamics is qualitatively reproduced with the phenomenological parameter 𝛾=30cm1. The shown Figure 5(b) is the population dynamics with Ω=300 cm−1 corresponding to Figure 4(e), and the suppression effect works fine as expected. On the other hand, shown in Figure 5(c) is the result for static detuning case without using the helper state |, that is, Ω=0 cm−1 and Δ=300 cm−1. As shown in the result, the suppression effect is insufficient compared to Figure 5(b), although the detuning is taken to be 300 cm−1, in both cases. In order to achieve the suppression level of Figure 5(b), one needs to increase the detuning up to Δ=3000 cm−1 as shown in Figure 5(d). Thus, the static detuning approach requires ten times larger than that in the present case. In the practical application to SCCl2, it is difficult to take required detuning Δ=3000 cm−1 because adjacent level exists within 1000 cm−1.

Further investigation on the effective Hamiltonian clarifies the difference in suppression efficiency shown above. The effective Hamiltonian for the CW-laser-induced Stark shift detuning case is given as [15] 𝐇(e)𝐀=+ΔΔ00Δ(Δ+𝑖𝛾)ΔΩ2𝑉21𝑉1𝑉2𝑉2𝑉1𝑉22,(30) whereas that of the static detuning case is [18]𝐇(e)𝐀=+1Δ00Δ𝑉Δ+𝑖𝛾21𝑉1𝑉2𝑉2𝑉1𝑉22.(31) One can see that Δ𝛾 or 1(𝛾/Δ) is required to neglect 𝛾 for static detuning case. On the other hand, the denominator of (30) shows that 1(𝛾Δ/Ω2) is required. Note that the condition 1(𝛾Δ/Ω2) can be satisfied with smaller (Stark shift) detuning because of its Ω2 dependency. This is why ten times larger value (Δ=3000 cm−1) is needed in the static detuning approach in Figure 5(d). Note also that choosing small Δ helps to satisfy the exclusion of dissipation in the present method. Another difference is lying in the laser parameter dependence of the Rabi frequency driven by the effective Hamiltonian. In the static detuning approach, the control parameter Δ affects both decomposition condition and the time scale of the isolated system. On the other hand, in the present method, the decomposition condition is mainly determined by Ω, while the timescale of Rabi oscillation can be adjusted by changing Δ as far as ΩΔ stands. From a viewpoint of quantum control, the present approach may offer versatile control scheme because it has more handling parameters compared to the static detuning approach.

Apparent drawback of the present method is possible side effects induced by the intense CW-lasers, such as multi-photon transitions which eventually leads to ionization, or unexpected near-resonant transition. We, here, consider the limitation of the present method when it is applied to the realistic molecular system as shown above. Note that some results shown above are rather extreme cases. In practice,the Ω=300 cm−1 condition in Figure 4(e) corresponds to CW-laser intensity of ~1011 W/cm2, which may cause undesired resonant multiphoton excitation. However, milder condition, such as Figure 4(d) corresponding to ~1010 W/cm2 still works fairly well. Note also that the intensities of weak lasers for 𝑉1 and 𝑉2 can be taken lower than the decomposition field Ω by one order, which denotes that there are very few side processes caused by 𝑉𝑖 interactions. Since resonant optical processes mainly occur under this condition, it could be possible to suppress undesired transitions by choosing helper level and laser frequency avoiding unexpected (near)resonances with care. Taking as large Δ value as possible also makes it possible to shorten the irradiation time of CW Laser, since we can expect faster Rabi-oscillation dynamics for control.

It is known that general IVR time scale varies from subpico to several hundreds ps [22]. For the IVR system with 𝛾100 cm−1~10−1 cm−1, applicability of the present method becomes more plausible, because required CW-laser intensity turns out to be around 108 W/cm2109W/cm2.

Finally, we should mention applicability of RWA for the present study. As the laser power increases, validity of RWA becomes doubtful. In order to verify this, we have carried out numerical calculation in which we treat the laser fields semiclassically. Overall behaviors of population dynamics agree with the ones obtained with RWA. The only difference is additional high-frequency components with small amplitudes onto the Rabi oscillation in the semiclassical treatments. Thus, we have concluded that the RWA are appropriate approximation for the present work.

4. Summary

We have presented that the dissipative dynamics due to the dephasing, such as IVR, can be suppressed by utilizing the effective decomposition brought by the intense CW lasers. One of the unique features of the present scheme is that one can treat it within a time-independent picture with the help of RWA, which makes the problem easy to handle. Consequently, one can avoid the situation that the control laser fields becomes too complicated. Another intriguing feature is that one can control the time scale of the population dynamics by changing the detuning parameters and it can be taken even longer than that of IVR itself. This is quite advantageous feature from a viewpoint of control problem, because one can design the control field regardless of the time constant of the dissipation.


This paper was supported, in part, by the Japan Society for the Promotion of Science (JSPS), Grant-in-Aid for Scientific Research (C), KAKENHI (20550021).