# Suppression of Intramolecular Vibrational Energy Redistribution by Intense CW-Laser Fields

Authors: M. Sugawara
Journal: Advances in Physical Chemistry
 (2011)

---

## Abstract

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.


---

## Body

## 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 [1–5]. 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) [9–13]. 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 Figure1. Before proceeding to the formulation, we explain several terms used. Here, |i〉 is the initial state, which is typically taken to be the lowest vibrational state on the ground electronic state. The state |m〉 is the intermediate state, which optically connects the two states |i〉 and |f〉. The state |f〉 is the control target state, that is, we aim to transfer the initial population to |f〉 avoiding the population loss due to IVR process. Note that |m〉 is optically accessible from |i〉, |f〉, and |h〉 and also coupled with {|Bj〉}(j=1,2,…,M) through the intramolecular interactions {vj}. The IVR, that is, the irreversible population flow to {|Bj〉}, occurs through |m〉, thus, we call this state “doorway state”, hereafter. In addition to that, we introduce the special state |h〉 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 |h〉 and |m〉 with interaction strength Ω.Figure 1
Schematic energy level diagram of a general IVR model system together with definitions of A, B and C spaces.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|i〉 and final state |f〉. We define B space as the one spanned by a doorway state |m〉 and a helper state |h〉, which we utilize for realizing effective decomposition. C-space is defined as the space spanned by IVR background manifolds {Bj}.The total Hamiltonian matrix under the rotating wave approximation (RWA) is given as(1)H=H0+λH′≡(ΔA000WB000FC)+λ(0VAB0VBA0VBC0VCB0).
Here, we introduce the parameter λ so as to clearly specify the perturbation order. The optical interactions between |i〉 and |m〉 (V1), and |m〉 and |f〉 (V2) are taken to be much smaller than Ω. Minor matrices in (1) are defined as(2)ΔA=(Δ00Δ),WB=(Ω00-Ω),FC=(ϵ10000ϵ20000⋱0000ϵM),
where Δ denotes the detuning introduced for AB transition (see Figure 1), while {ϵi} (i=1,2,…,M) denotes the eigenenergies of the background manifolds in C-space. Note that we prediagonalize the B space Hamiltonian so as to make WB diagonal matrix. The matrices VAB (2×2) and VBC (2×M) correspond to the optical and static interactions, respectively.The Schrödinger equation in the matrix representation is given as(3)H⋅U=U⋅E,
where E is a diagonal matrix with eigenvalues {Ei} (i=1,2,…,N) in its diagonal elements, and U is the unitary matrix consisting of corresponding eigenvectors ci (i=1,2,…,N) as U≡(c1,c2,…,cN). Here, we expand U and E as(4)U=U(0)+λU(1)+λ2U(2)+⋯,E=E(0)+λE(1)+λ2E(2)+⋯,
where(5)U(n)=(UAA(n)UAB(n)UAC(n)UBA(n)UBB(n)UBC(n)UCA(n)UCB(n)UCC(n)),E(n)=(EA(n)000EB(n)000EC(n)).
Here, EA(n), EB(n), and EC(n) are diagonal matrices, whose diagonal elements are nth perturbation energies. Inserting  (1), (4) into (3) and comparing λn (n=0,1,2) terms gives the followings:(6)H0⋅U(0)=U(0)⋅E(0),(7)H0⋅U(1)-U(1)⋅E(0)+H′⋅U(0)=U(0)⋅E(1),(8)H0⋅U(2)-U(2)⋅E(0)-U(1)⋅E(1)+H′⋅U(1)=U(0)⋅E(2).
From (6), one obtains the following equations for minor matrices:(9)ΔA⋅UAA(0)=UAA(0)⋅EA(0),(10)WB⋅UBB(0)=UBB(0)⋅EB(0),(11)FC⋅UCC(0)=UCC(0)⋅EC(0),
while off-diagonal blocks of U(0) are all zero matrices. As ΔA commutes with UAA(0),  (9) gives EA(0)=ΔA, which leaves UAA(0) undetermined. On the other hand, from (10), (11) one obtains that UBB(0) and UCC(0) are both determined as identity matrices of 2×2 and an M×M, respectively, whereas EB=WB and EC=FC. The first two terms of the left hand side of (7) are given as (12)H0⋅U(1)-U(1)⋅E(0)=(0ΔA⋅UAB(1)-UAB(1)⋅WBΔA⋅UBC(1)-UBC(1)⋅FCWB⋅UBA(1)-UBA(1)⋅ΔAWB⋅UBB(1)-UBB(1)⋅WBWB⋅UBC(1)-UBC(1)⋅FCFC⋅UCA(1)-UCA(1)⋅ΔAFC⋅UCB(1)-UCB(1)⋅WBFC⋅UCC(1)-UCC(1)⋅FC),whereas other terms are(13)H′⋅U(0)=(0VAB⋅UBB0VBA⋅UAA(0)0VBC⋅UCC0VCB⋅UBB0),UAA(0)⋅EAA(1)=(UAA(0)⋅EAA(1)000UBB(0)⋅EBB(1)000UCC(0)⋅ECC(1)).
Comparing the component of A space, one obtains the equation UAA(0)·EAA(1)=0, which leads to EAA(1)=0 since UAA(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(14)WB⋅UBA(1)-UBA(1)⋅ΔA+VBA⋅UAA(0)=0.
Considering UBA(1)·ΔA=ΔB·UBA(1),  (14) can be rewritten as(15)(WB-ΔB)⋅UBA(1)=-VBA⋅UAA(0),
which leads to(16)UBA(1)=-(WB-ΔB)-1⋅VBA⋅UAA(0).
Here, ΔB denotes the identity matrix of B space multiplied by Δ. The prefactor matrix (WB-ΔB)-1 in (16) is given in a simple form as(17)(WB-ΔB)-1=(1Ω-Δ001-Ω-Δ).
Note that UBA becomes zero matrix as Ω becomes large or Ω≫Δ, which leads to the suppression of AB-space mixing.Comparing CA-block of  (7), one obtains(18)FC⋅UCA(1)-UCA(1)⋅ΔA=0,
which is rewritten as(19)(FC-ΔC)⋅UCA(1)=0,
where ΔC is M×M identity matrix multiplied by Δ. It is readily seen that UCA(1)=0 since FC-ΔC≠0. Thus, there is no first-order CA-mixing.Next, we consider the second-order terms. WithE(1)=0,  (8) becomes(20)H0⋅U(2)-U(2)⋅E(0)+H′⋅U(1)=U(0)⋅E(2).
Since we focus on the dynamics of A space, we extract AA-block of  (20) as(21)ΔA⋅UAA(2)-UAA(2)⋅ΔA+VAB⋅UBA(1)=UAA(0)⋅EA(2).
Considering that the first two terms of (21) are canceled out and inserting (16), one obtains(22)VAB⋅(WB-ΔB)-1⋅VAB⋅UAA(0)=UAA(0)⋅EA(2).
Note that (22) is an eigenvalue equation, which determines the second-order energy correction EA(2) onto A space as well as corresponding eigenvectors UAA(0). Here, we define the A space effective Hamiltonian HA(eff) as(23)HA(eff)≡ΔA+VAB⋅(WB-ΔB)-1⋅VAB,
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 H′ 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 fieldVi, 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|m′〉 adjacent to the intermediate state |m〉. Typically, the level spacing Δ′ between |m〉 and |m′〉 is about 1000 cm-1. The CW-laser (Ω) applied for decomposition turns out to be off-resonant by Δ′, which is large enough to exclude the undesired transition to |m′〉.


## 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 |m〉 with uniform strength v. Shown in Figure 2(a) is the population dynamics under field-free condition while the initial state is taken to be |m〉. The coupling parameter v 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 {|Bj〉} occurs due to the dephasing process between numerous eigenstates contained in the initial state |m〉.Population dynamics of Bixon-Jortner model. The dark gray, light gray, and white areas correspond to the population of|i〉, |f〉, and |m〉, respectively. The hatched area denoted by IVR corresponds to the population distributed over background manifolds |Bj〉. (a) |m〉 is taken to be the initial state. (b) Two CW-lasers, V1 and V2, are applied corresponding to the |i〉-|m〉 and |m〉-|f〉 transitions. Both V1 and V2 are taken to be 0.1Ω. (c) Intense CW-laser is introduced to make the helper state |h〉 strongly coupled with |m〉 together with the relatively weak lasers V1 and V2.
(a)(b)(c)Next, we try a naive population control in which control target is the population transfer from|i〉 to |f〉. We introduce two resonant CW-lasers corresponding to the transitions, |i〉↔|m〉, and |m〉↔|f〉, respectively. The total Hamiltonian corresponding to (1) is given as
(24)H=(Δ0V10⋯00ΔV20⋯0V1V20vvv00vϵ1⋯0⋮⋮v⋮ϵi⋮00v0⋯ϵM),
where the vertical and horizontal dots are all zeroes. Shown in Figure 2(b) is the population dynamics with the initial condition |i〉. Both optical interaction parameters V1 and V2 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 |f〉 is only 25%. Thus, the control objective is hardly achieved with this approach.Now, we introduce a helper state|h〉 in order to suppress the population loss due to IVR. The total Hamiltonian represented by the original basis set is given as
(25)H̃=(Δ0V100⋯00ΔV200⋯0V1V20Ωvvv00Ω00⋯000v0ϵ1⋯0⋮⋮v⋮⋮ϵi⋮00v00⋯ϵM),
where Ω is taken to be the optical interaction between |m〉 and |h〉 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,
(26)H=(Δ0V′V′0⋯00ΔV′′V′′0⋯0V′V′′Ω0v′v′v′V′V′′0-Ωv′v′v′00v′v′ϵ1⋯0⋮⋮v′v′⋮ϵi⋮00v′v′0⋯ϵM),
where V′=V1/2, V′′=V2/2, and v′′=v/2. Shown in Figure 2(c) is the population dynamics under the condition V1=V2=0.01Ω. It is shown that the population oscillates between |i〉 and |f〉, and there is very little population leak into the background manifolds {|Bj〉}. 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(27)HA(eff)=(Δ00Δ)+ΔΔ2-Ω2(V12V1V2V2V1V22).
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](28)ĤM=ĤM(0)+ĤM′,
where(29)ĤM(0)=∑i=16ωi(ni+12)+∑i=16χii(ni+12)2,ĤM′=Va3-Δn.
Here Δn=∑i|ni-ni′|. 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 |n1,n2,n3,n4,n5,n6〉. 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 |i〉 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 |f〉≡|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 V1=V2=30cm−1 and Δ=30cm−1. The corresponding laser amplitudes for V1 and V2 are estimated as 9.6×109W/cm2 and 5.0×109W/cm2, respectively [20].Figure 3
Schematic potential curve with energy levels for SCCl2 molecule.Shown in Figure4 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 (|i〉 and |f〉). 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 |f〉 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)).Population dynamics of SCCl2 molecule with changing coupling parameter Ω. Corresponding laser power is as follows: (a) 0.0 W/cm2, (b) 2.48×109 W/cm2, (c)9.91×109W/cm2, (d) 3.96×1010W/cm2, (e) 8.92×1010 W/cm2. The hatched area denoted by IVR corresponds to the population distributed over background manifolds |Bj〉. The white area in (a) denotes the population of |m〉, whereas the area painted by black background with horizontal lines in (b) corresponds to the population of |h〉. Note that the population of |m〉 and |h〉 is so little and hardly seen in Figures 4(c) and 4(e).
(a)
Ω=0 cm-1(b)
Ω=50cm-1(c)
Ω=100cm-1(d)
Ω=200cm-1(e)
Ω=300cm-1Here, 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, Δ≫V1,V2 in the static detuning case whereas Ω≫Δ,  Ω≫V1,V2, 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 (|i〉 and |f〉) and B space (|m〉 and |h〉). Instead of omitting C-space (background states), we introduce a phenomenological dephasing parameter γ associated with the intermediate state |m〉, that is, we introduce an imaginary component γ for the eigenvalue of |m〉. 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 V1=V2=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 γ=30cm-1. 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 |h〉, 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.Comparison of suppression efficiency of Stark shift detuning (present method) and static detuning approaches. Shown in each figure is population dynamics of 4-level system with phenomenological dephasing parameterγ=30 cm−1 on |m〉. V1 and V2 are taken to be 30 cm−1. Other parameters are taken to be (a) Ω=0 cm−1, Δ=30 cm−1, (b) Ω=300 cm−1, Δ=30 cm−1, (c) Ω=0 cm−1, Δ=300 cm−1, and (d) Ω=0 cm−1, Δ=3000 cm−1.
(a)(b)(c)(d)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] (30)HA(eff)=(Δ00Δ)+Δ(Δ+iγ)Δ-Ω2(V12V1V2V2V1V22),
whereas that of the static detuning case is [18](31)HA(eff)=(Δ00Δ)+1Δ+iγ(V12V1V2V2V1V22).
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 V1 and V2 can be taken lower than the decomposition field Ω by one order, which denotes that there are very few side processes caused by Vi 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/cm2~109W/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.



---

*Source: 584082-2011-07-24.xml*