Ultrafast dynamics in the vicinity of quantum light-induced conical intersections

Nonadiabatic effects appear due to avoided crossings or conical intersections (CIs) that are either intrinsic properties in field-free space or induced by a classical laser field in a molecule. It was demonstrated that avoided crossings in diatomics can also be created in an optical cavity. Here, the quantized radiation field mixes the nuclear and electronic degrees of freedom creating hybrid field-matter states called polaritons. In the present theoretical study we go further and create CIs in diatomics by means of a radiation field in the framework of cavity quantum electrodynamics. By treating all degrees of freedom, that is the rotational, vibrational, electronic and photonic degrees of freedom on an equal footing we can control the nonadiabatic quantum light-induced dynamics by means of CIs. First, the pronounced difference between the the quantum light-induced avoided crossing and the CI with respect to the nonadiabatic dynamics of the molecule is demonstrated. Second, we discuss the similarities and differences between the classical and the quantum field description of the light for the studied scenario.

quantum control problem. In this framework the nonadiabatic dynamics arises due to the strong coupling between the molecular degrees of freedom and the photonic mode of the radiation field which can alter the molecular levels by controlling the dynamics of basic photophysical and photochemical processes.
In most of the theoretical descriptions the molecules are treated via a reduced number of degrees of freedom or by some simplified models assuming two-level systems. Here, only one vibrational mode which is strongly coupled to the electronic and photonic degrees of freedom is taken into account resulting in a new set of 'cavity induced' or 'polariton' surfaces in the molecular Hamiltonian [21][22][23][24][25]. These polariton surfaces can form CIs under special conditions [23] but are not expected to cross each other in general forming only 'avoided crossings'. This scenario resembles a one-dimensional (1D), semi-classical treatment of the light-induced avoided crossing occurring in the nonadiabatic gas-phase molecular dynamics. Recent studies of the time-dependent state population in the NaI molecule have successfully demonstrated that the 1D results obtained by the semi-classical approach [50] are fully consistent with the results obtained for NaI and quantized light field [24].
Adding another molecular DOF, which can either be a second vibrational mode in case of poly atomics or the rotational angle between the molecular axis and the polarization axes in the cavity, one can fully describe the photon-induced quantum dynamics in the framework of quantum LICI. CIs can even be formed in diatomics due to the availability of the two independent nuclear degrees of freedom, which are essential for forming a 2D branching space. This picture can now be extended in a straightforward way from the simplified 1D model with a second molecular mode to treat the rotational, vibrational, electronic and photonic degrees of freedom on an equal footing. In the following we call the latter description the 2D description.
By including the molecular rotation in the cavity treatment, recently we have investigated the field-dressed rovibronic spectrum of diatomics in the framework of cavity quantum electrodynamics [51]. Incorporating the concept of LICIs with the quantized radiation field, similar impact on the adiabatic spectrum has been found as in the case of classical dressed situation [52]. We have demonstrated that the 'intensity borrowing' effect can cause a significant variation in the pattern of the molecular spectra irrespective of the origin of the CIs. It can be either natural or light-induced and the latter can be created by classical or quantum light.
In the present work we go further and investigate theoretically how the dynamics of the lithium fluoride molecule, which already possesses a natural avoided crossing, is affected by quantum light-induced CI. This work complements previous theoretical investigations [24,50], where only avoided crossings were created by quantum light [24], but not CIs [50]. Here we investigate the nonadiabatic quantum dynamics by incorporating the concept of LICIs with the quantized radiation field.
The aim of this letter is two-fold. First, we would like to study the quantum light-induced nonadiabatic dynamics of the LiF molecule both in 1D and 2D frameworks to demonstrate the difference between the effects of the radiation field-induced avoided crossing (AC) and CI. The underlying dynamics is mainly governed by the interplay between one of the quantum light-induced phenomena (either LIAC or LICI) and of the natural avoided crossing which is present in the field-free molecule. Second, the similarities and the differences between the classical and the cavity Fock radiation field description of the light are investigated. We discuss for both the LIAC and LICI to what extent the different physical scenarios, a molecule in a laser field and a molecule in a cavity, show similar or different results.
The general form of the Hamiltonian in the basis of the adiabatic potential energy surfaces Σ 1 and Σ 2 of LiF can be given as Here we assume the electric dipole approximation and that the molecule interacts with only a single mode of the cavity. In equation (1), the first term represents the rovibrational kinetic energy of the LiF molecule with R and θ being the vibrational and rotational degrees of freedom, respectively. M r is the reduced mass of the LiF molecule, L θ is the angular momentum operator (with m=0 fixed) and the 1 symbol represents the 2×2 unit matrix. The second term contains the field-free adiabatic potential curves S V 1 and S V 2 (see figure 1(a)). The K field-free NAC operator is approximated as , where τ(R) is first order NAC term, shown in figure 1(b). The third term in equation (1) represents the harmonic oscillator description of the photon mode with the unit-less photon displacement coordinate x. Here ω c is the cavity mode angular frequency. The last term of equation (1) describes the interaction between the molecule and the quantized field. Here χ is the cavity coupling strength, while μ i and μ ij (i,j=Σ 1 , Σ 2 ) are the permanent and transition dipoles, respectively. In the actual calculations ÿω c =3.995 eV is taken which corresponds to the resonant coupling of the Σ 1 and Σ 2 states around R=5 a.u. The coupling strength χ ranges from 0.0012 to 0.04 to simulate moderate and strong coupling strengths.
The initial wave packet at t = −10 fs is created from the product of the rovibrational ground state of Σ 1 located around R≈3 a.u. and one of the Fock states and is placed on the lower adiabatic potential Σ 1 . The Fock states are considered eigenstates of the w - Hamiltonian.
The form of the time-dependent Hamiltonian in the classical framework in the basis of the adiabatic states Σ 1 and Σ 2 of LiF reads: In equation (2), the first term represents the vibrational and rotational kinetic energy (the same as in equation (1)), while the second term contains the field-free S V 1 and S V 2 potential curves and the K NAC operator. The third term of equation (2) describes the laser-molecule interaction in the dipole approximation. Here ε 0 is the amplitude of the electric field, ω c is the angular frequency of the laser, f (t) is the envelope function which is set to unity during the whole propagation (t final =1000fs). μ i and μ ij (i, j=Σ 1 , Σ 2 ) are the permanent and transition dipoles, respectively. The actual value of the applied laser energy was set to ÿω c =3.995 eV, which resonantly couples the electronic states at R=5 a.u., and the peak laser intensity ranges from I 0 =3×10 11 to 3×10 13 W cm −2 .
In both the cavity and classical calculations a linearly polarized resonant pump pulse is applied to initiate the dissociation dynamics. The form of this pump pulse is given as e w · ( ) · ( ) t g t cos pump pump where ε pump is the peak electric field strength, ω pump is the resonance angular frequency and = is the envelope function. Applying a laser pulse with a center frequency of ÿω pump =7 eV, T p =20 fs pulse duration, and I pump = 4.8 × 10 13 W cm −2 peak intensity, 35% of the total population is excited to Σ 2 . This wave packet then starts to oscillate and gradually dissociate on Σ 1 via the AC and LICI or LIAC. The multi configurational time-dependent Hartree (MCTDH) method [53,54] has been applied to solve the time-dependent Schrödinger-equation characterized by either equation (1) or (2). The R DOF was defined on a sin-DVR (discrete variable representation) grid (N R basis elements for R=1.6-60 a.u.). The rotational DOF, θ, was described by N θ Legendre-polynomials, q ( ) P cos l m with m=0 and l=0, 1, K, N θ −1. The photon displacement coordinate x was described by N x Hermite-polynomials, H n (x) with n=0,1, K, N x −1. In the MCTDH wave function representation, these primitive basis sets (ξ) are then used to construct the single particle functions (f) whose time-dependent linear combinations form the total nuclear wave packet (ψ) , , The actual number of basis functions were N R =1169, N θ =271 and N x =100 for the vibrational, rotational and photon modes, respectively. The number of single particle functions for the three DOFs and on both the Σ 1 and Σ 2 electronic states were ranging from 10 to 50. The values of n R =n θ and n x were chosen depending on the actual value of the χ cavity coupling strength and I 0 peak laser intensity. In order to minimize unwanted reflexions and transmissions caused by the finite length of the R-grid, complex absorbing potentials (CAP) have been employed at the last 10 a.u. of the grid. The total propagation time was set to t final =1000 fs and the state populations and dissociation rates were calculated according to In equation (5), Θ is the Heaviside step function, R D =20 a.u. is the starting point of the dissociation region, and −iW is the CAP. To calculate the potential energy, the dipole moment, and the NAC curves of the LiF molecule ( figure 1(b)), the program packages Molpro [55] was used. These quantities were calculated at the MRCI/CAS(6/12)/aug-cc-pVQZ level of theory [56]. In particular, τ(R) has been computed by finite differences of the MRCI electronic wave functions. The number of active electrons and molecular orbitals in the individual irreducible representations of the C 2v point group were A 1  2/5, B 1  2/3, B 2  2/3, A 2  0/1. We start with the dynamics of the bare molecule (no cavity). Figures 2(a) and (b) show the time dependent populations of both the ground and excited states. The length of the red interval in vertical direction denotes the amount of the population in the excited states. The initial set up can be seen at t<−10 fs, when the whole population is on the ground state. Then a linearly polarized resonant pump pulse (the details of it are given in the introduction) is applied to initiate the dissociation dynamics. From now on all population will be related to the one which is transferred to the upper states during the pump process. The wave packet of the transferred population starts to oscillate on the excited state potential curve and reaches the avoided crossing at t≈80 fs. Here about 80% of the population is transferred back to the ground state due to the intrinsic NAC. The area colored in blue shows the amount of the population on the ground state, which begins to dissociate at t ≈ 140 fs. At t≈180 fs approximately 80% of the population dissociates (marked in white color) and does not show up again in the dynamics. In the bare molecule the Hamiltonians for 1D and 2D framework are effectively the same, which results in identical dynamics. The bumpy structure on the curves at t≈800 fs denotes that the wave packet has returned to the avoided crossing region.
By switching on the cavity coupling between the photonic and the molecular degrees of freedom, the quantum light-induced nonadiabaticity now competes with the intrinsic one. In the 1D calculations only the vibrational molecular DOF mixes with the photonic mode and the electronic states, creating an avoided crossing between the polariton states. In the 2D scheme both the rotational and vibrational modes are accounted for and a LICI can be formed even for a diatomic molecule. We begin by investigating the vacuum state of the cavity mode (n=0). We performed 1D and 2D calculations for several values of the cavity coupling and the results of the calculation are shown in figure 2. It can be seen that with increasing coupling strength χ the 1D and 2D results increasingly deviate from each other and the impact of the LICI becomes more prominent. For a coupling parameter of (χ=0.02) large amplitude oscillations can be observed between the electronic surfaces both in the 1D and 2D populations. This is due to the increased splitting at the LIAC or LICI for larger coupling strengths creating decoupled dressed state surfaces [21,23]. However, at the largest coupling parameter (χ=0.04), owing to the different shapes of the upper polaritonic surface in different directions, the period of this oscillation strongly depends on the molecular orientation which-after a few periods-leads to a wash out of the fingerprints of the oscillation from the total populations depicted in the figure. The dissociation itself is less suppressed in the cavity LICI picture than in the case of the quantum light-induced avoided crossing. This is clear evidence that in the cavity the bond hardening effect [57] is less efficient in 2D than in 1D. That is, even when the quantized radiation field couples to the molecular degrees of freedom, the ultrafast decay channel created by the LICI is more efficient with respect to the fast population transfer than the avoided crossing. Although the time evolution of the quantum dynamics is determined by the interplay between the radiationfield induced and intrinsic nonadiabaticity, at sufficiently strong coupling region the quantum LICI dominates Figure 2. Population of the ground and excited electronic states as a function of time computed in 1D and 2D. Panels on the left-hand side (a)-(g) show the 1D results, while those on the right-hand side (b)-(h) depict the 2D results for a given cavity coupling strength; from top to bottom χ=0, χ=0.005, χ=0.02 and χ=0.04. The blue and red areas correspond to the population of the Σ 1 and the Σ 2 states, respectively. The white regions show the ground state population being in the dissociation region (R>20 a.u.). In all the panels the Fock vacuum state (n=0) is considered as an initial state along the cavity mode. the process. Similar effects have been found for natural avoided crossings, natural CIs, and for LIACs and LICIs induced by classical light field. In this work we have demonstrated the dynamics of quantum light-induced avoided crossings and quantum light-induced CIs for the first time.
To better understand the differences between the 1D and 2D results, we compare the dissociation rates for different photon number of Fock states (n=0, 5, and 10) in figure 3. Again, the dissociation rate in the bare molecule (χ=0) is the same in 1D and 2D. However, for the vacuum Fock state (n=0) a minor difference between the two schemes can be recognized. The structure of the dissociation curves is similar but the effect, that the dissociation is more efficient in 2D, can already be seen here. The efficiency of LICI is increasing compared to the 1D model when we use Fock states with n=5 and n=10 as an initial condition. Increasing the photon number of the initial state the dissociation rate of the 2D model less supressed than that of the 1D one. The decay channel provided by the quantum LICI is more efficient for transferring the population to the lower polaritonic state, which leads to higher dissociation rates. The effect is similar to what we experience in the case of increasing coupling strengths. Namely, the increasing photon number provides an increasing coupling strength, which is enhanced by + n 1 , leading to a more efficient LICI. The rotation of the molecule starts slowly, therefore up to t=220 fs the 1D and 2D curves are practically the same for all the studied photon number states. We now discuss figure 3(b) where we see two crossings between the n=5 and n=10 dissociation curves. Because the bond hardening effect is stronger for n=10, evidently the dissociation starts slower. At t=330 fs the effect of rotation suddenly becomes apparent resulting in crossings in the 2D curves (for n=5 and 10). Between the period of t=350 fs and t=850 fs the rotation already plays an important role in amplifying the effect of the LICI. This effect is particularly pronounced for n=10.
To gain more insight into this phenomena we have analyzed the two dimensional wave packet density functions y q and found that the n=5 and n=10 Fock states differ significantly from the n=0 vacuum. Namely, a much stronger alignment and bond hardening effect can be realized in the Fock space including photons than in the case of a vacuum which then strongly suppresses the dissociation rate. At weaker coupling strengths the dissociation rates obtained by the LICI and LIAC are more or less similar, independently from the photon number of the Fock states, while in the stronger coupling regime the LICI provides a significantly larger amount of dissociation rate. In the case of quantum light the molecular degrees of freedom interact directly with the photonic degrees of freedom and can modify the state of the light field through stimulated emission. This is not possible with a classical description of light field.
In figure 4 we compare the results of the cavity mode for a photon number n=5 with the classical description as obtained from the 1D and 2D calculations. A similar comparison for the 1D model, but for a photon number n=0, has been made between the quantum light and classical results of the Na 2 molecule ground state population [50]. This is comparable to our earlier findings where independently from the applied coupling strength practically no considerable difference has been found between the different 1D simulations. In contrast, in the 2D calculations significant differences have been found in the dissociation rates when comparing the classical and quantum light model. At weaker coupling strengths the difference is not so prominent, but becomes more significant at stronger couplings. At really short time dynamics the LIAC and LICI behave similarly but differ significantly at longer time scale.
In summary, we could show that the dynamical properties of diatomic molecules can be strongly modified by quantized light in an optical cavity. By using the LiF molecule as a showcase example, we demonstrated that the LICI created by quantum light allows for a more efficient population transfer than a LIAC given a sufficiently large coupling strength. The stronger the cavity coupling, the more prominent the effect. This difference can be explained by the fact that the LICI retains a degeneracy between the dressed states even for large coupling strengths. In contrast, the dressed state curves of the 1D model become increasingly separated for larger coupling strengths, which leads to a decreased mixing between the nuclear degrees of freedom and the electron + photon degrees of freedom. In addition, for the case of LIAC a close similarity has been found between the classical and the cavity radiation field description of the light for all the studied coupling strengths and Fock states.
A significant difference has been found between the dynamics governed by the quantum Fock state and the classical light description of the LICI in the Hamiltonian. The stronger the coupling is, the larger difference becomes between the dynamics of the LICI formed by quantum and classical light. By increasing the coupling strength in the cavity description, the dissociation rate is more and more supressed compared to the classical 2D calculation. It has also been realized that even in the quantum description of the LICI the photons which are present in the Fock state (n=5 and n=10) can also play an important role comparing to the Fock vacuum situation. The detailed dissection of the dynamics and the analysis of the interplay between photonic degrees and molecular degrees of freedom will be subject of a future study.
We note that there is more potential for exploring quantum LICIs in poly atomic molecules, covering a larger class of chemically relevant molecules. Here, the cavity induced nonadiabatic molecular dynamics can be studied in the absence of molecular rotations as they already provide the necessary number of vibrations to form CIs.