Ultrafast Resonant Photon Emission from a Molecule Driven by a Strong Coherent Field in Terms of Complex Spectral Analysis

: In this study, we investigate the time–frequency-resolved resonant photon emission from a molecular vibrational oscillator driven by a monochromatic coherent external field. Using the complex spectral analysis of the Liouvillian, which integrates irreversible dissipative phenomena into quantum theory, we elucidate the fundamental processes of photon emission. Indeed, our analytical approach successfully decomposes the emission spectrum into two intrinsic contributions: one from a resonance eigenmode and another from continuous eigenmodes. These components are responsible for incoherent luminescence and coherent scattering photon emission processes, respectively. Our results show that while spontaneous emission dominates in the early stages of the emission process, coherent scattering gradually becomes more pronounced with time. Furthermore, destructive quantum interference between the two components plays a key role in determining the overall shape of the emission spectrum.


Introduction
Spontaneous emission from an excited atom is a fundamental process in light-matter interactions.It is an irreversible dissipative process due to quantum fluctuations of the electromagnetic field and energetic resonance between the atomic excited states and the radiation field.The Weisskopf-Wigner theory explains the Lorentzian-shaped spontaneous emission spectrum of an excited atom with use of the solutions of the Heisenberg equation, where the dissipation effect of an atomic excitation is taken into account phenomenologically by introducing a decay constant [1,2].This theory assumes that the atom is in the excited state as the initial state and does not address how the excitation process affects the photon emission process.In contrast, Yutaka Toyozawa and collaborators studied a correlation of excitation and photon emission processes in a resonant light emission, treating it as a sequential coherent quantum process.They revealed that two different types of photon emission processes coexist: coherent light scattering and incoherent luminescence [3][4][5].While in the former the emitted photon energy is correlated with the incident light energy, in the latter, the correlation disappears due to the dissipation effect.
Recent progress in the development of strong laser sources allows us to explore various new types of light emission from light-driven materials [6][7][8][9][10], such as high harmonic generation and high-order sideband generation [11][12][13][14][15][16][17].These are often considered as highorder nonlinear optical phenomena, where coherently driven electronic systems give rise to substantial electronic currents or polarizations that subsequently act as sources of classical radiation [15,[18][19][20].This leads to fundamental questions: how is quantum mechanical spontaneous emission in a weak excitation related to coherent classical radiation in a strong coherent driving situation, and how can coherent scattering and incoherent luminescence coexist?Answering these questions requires a quantum mechanical treatment of light emission processes, considering the whole system of light and matter.Although several studies have analyzed radiation processes quantum mechanically, they have been largely limited to the lower-order perturbation regime [21][22][23].
Even with a higher-order perturbation method, a fundamental challenge remains: how to describe dissipative spontaneous emission in a way that is consistent with quantum theory.Traditional approaches yield only time-reversible unitary evolution from hermitian Hamiltonian dynamics.To fill this gap, several extensions of quantum theory have been proposed to include irreversible dissipation [24][25][26][27][28][29][30].Among these, the contributions by Ilya Prigogine and collaborators [25,28,31] and by Arno Bohm and collaborators [24,30,32] are notable; they extend the vector space beyond the conventional Hilbert space to the rigged Hilbert space, or the Gel'fand triplet [33].In this extended framework, called complex spectral analysis, the time evolution generators-specifically, Hamiltonian and Liouvillian-can have complex eigenvalues, thus allowing irreversible processes.
Using complex spectral analysis, we have studied optical dissipative processes to elucidate the role of resonance states characterized by complex eigenvalues [34,35].Our recent investigations used this analysis in conjunction with Floquet theory to study highorder sideband generation in a coherently driven two-level system [36], revealing both transient and stationary spectral profiles.Through this analysis, we distinguished the resonant state contributions in the photon emission spectrum from the dressed continuous state contributions and identified instances of destructive interference between the two.These studies assumed an initial excited state for the two-level system and a modification of its energy levels by the driving field.However, this initial condition precludes the clarification of the energy correlation between incident and emitted light and, thus, hinders the study of the interplay between scattering and luminescence under strong driving fields.
In this paper, we present a theoretical formulation for the resonant photon emission of a harmonic oscillator driven by a strong coherent field.Our approach uniquely considers the quantized radiation field as an integral part of the whole system.Unlike in the previous studies [15,[18][19][20], we focus here on the irreversible spontaneous photon emission process under a strong driving field, which requires an extension beyond conventional quantum theory restricted to Hilbert space.We address this by solving the complex eigenvalue problem of the Liouvillian, the generator of the Heisenberg equation, to elucidate the resonance mode responsible for the irreversible spontaneous emission that manifests as incoherent luminescence.
In terms of complex spectral analysis, our results reveal the physical origins of resonance photon emission under strong coherent drive.The solution we provide is exact, eliminating ambiguities in the physical interpretation of the spectrum.We present the time-frequency-resolved photon emission spectrum and show that incoherent luminescence contributions can exist even under strong coherent drive.Furthermore, our results emphasize the important role of coherent photon emission at the molecular resonance.In particular, destructive interference plays a key role during the initial stages of photon emission, profoundly influencing the build-up of the spectral structure.
The structure of this paper is organized as follows.In Section 2, we introduce our model system and describe the Liouvillian, which serves as the generator for the time evolution of the quantum dynamical variables.Section 3 details the derivation of complex eigenmodes and presents the solution of the Heisenberg equation.The calculated results, including molecular excitation and the resonant photon emission spectrum, are discussed in Section 4. We conclude in Section 5 with a summary of our results and further discussions.We introduce the symplectic structure of the present system in Appendix A, derive the effective Liouvillian in Appendix B, and obtain the complex eigenmode operators in Appendix C. The solution of the Heisenberg equations and the spectral expressions are given in Appendix D.

Model
In this paper, we study the photon emission from a polarizable diatomic molecule driven by a time-dependent strong coherent driving electric field, where the stretching vibrational mode is described by a harmonic oscillator interacting with the time-dependent external field.The Hamiltonian operator for this system is then given by Ĥ where the hat denotes the operator, â ( â † ) is the annihilation (creation) boson operator for the harmonic oscillator representing the stretching vibration of a diatomic molecule with natural frequency ω 0 .In this paper, we set the reduced Planck constant, h = 1.In what follows, we refer to this oscillator as the "molecular oscillator" or just the "molecule".The emitted radiation field is treated quantum mechanically, as described by the second term in Equation ( 2), where bk ( b † k ) are the annihilation and creation boson operators for the radiation field with wavenumber k and frequency ω k .Specific wavenumber dependencies of these parameters are used in the numerical calculation described in Section 4. The third term in Equation ( 2) represents the interaction between the molecular oscillator and the radiation field with the real coupling constant g k .The Hamiltonian Ĥ, known as the Friedrichs model, has been used extensively to describe quantum mechanically an irreversible spontaneous emission process from an excited atom [2,25,37].The interaction of the molecular oscillator with an external time-dependent driving field is represented by Ŵ(t) (3).
Here, we define a set of Heisenberg operators as a vector: where { bk (t)} and { b † k (t)} are infinite dimensional row vectors of annihilation and creation operators, respectively, for continuous k variables, and T is the transpose.We solve the Heisenberg equation, with the initial condition In Equation (5), L(t) and L denote the Liouvillian superoperator defined by the commutation relation with Ĥ(t) and Ĥ, respectively: The Liouvillian superoperator L is represented by a block diagonal matrix for the annihilation operator set ( â, { bk }) and the creation operator set where the submatrix L B is represented by a Hermitian matrix, The last term of Equation ( 5) is due to the interaction of the molecule with the driving field represented by To solve the Heisenberg equation ( 5), we first find invertible transformation matrices Φ and Φ to diagonalize the time-independent Liouvillian L given by Equation ( 8): Since the total Liouvillian takes the block diagonal form Equation ( 8), we can write so that and where ⋆-transformation is defined in Equations (A10) and (A11).Below, we show how to construct Φ, Φ and obtain Z in Appendixes A and C in terms of the Brillouin-Wigner-Feshbach projection method.With use of the transformation Φ, we define the complex eigenmode of the total system as Equation (A7): where ξ is given in Equation ( 6).
Multiplying Φ from the left of the Heisenberg Equation ( 5) and using Equations ( 10) and ( 14) yields a decoupled differential equation in terms of the eigenmodes: Therefore, the identification of the complex eigenmodes should be useful not only for the calculation, but also for the interpretation of the physical origin of the spectrum.Our method differs from traditional approaches in that conventional methods typically solve the time-dependent Schrödinger equation first and then treat the irreversible photon emission process using first-order perturbation theory [21][22][23].This standard approach does not correctly describe an irreversible dissipation process from the molecule to the free radiation field.

Complex Eigenmodes and the Spectrum
The spontaneous emission process is one of the fundamental irreversible processes due to the combination of the quantum vacuum fluctuation and the resonance singularity.In order to correctly treat the irreversible process starting from the hermitian Hamiltonian, we resort to quantum dynamics interpreted in terms of the complex eigenmodes based on the solutions of the complex eigenvalue problems of the Liouvillian [24,25,32,38].In this Section, we obtain the complex eigenmodes of the field-free Liouvillian L given in Equation (8).
Since the Liouvillian matrix L in Equation ( 8) is in block diagonal form, we find the transformation of Φ B and ΦB by solving the complex eigenvalue problem of the inifinite dimensional hermitian submatrix L B in Equation ( 9), where we note that the matrix structure of L B is the same as the Friedrichs model Hamiltonian used for spontaneous emission from a two-level atom [37].We diagonalize L B with use of the symplectic transform Φ B and ΦB as shown in the Appendixes A and C: The transformations Φ B and ΦB correspond to the right-and left-eigenvectors of L B , respectively: Since the structure of L B is the same as that of the Friedrichs Hamiltonian, we use the results obtained by the Brillouin-Wigner-Feshbach projection method using L B instead of the Friedrichs Hamiltonian: where |φ j ⟩ and ⟨ φj | are the right and left eigenvectors of L B with the eigenvalue z j , belonging to the rigged Hilbert space [29,30,33].The eigenvectors are represented by the unperturbed basis vectors as where we introduce the basis set of (|a⟩, {|b k ⟩}) for the representation of the submatrix L B : These eigenvectors satisfy bi-orthonormality and bi-completeness: where δ(•) is the Dirac delta function.
To reduce the eigenvalue problem of the infinite dimensional matrix L B to a finite dimensional problem, we use the Brillouin-Wigner-Feshbah projection method.First, we derive the effective Liouvillian for the molecular oscillator with use of the projection operator P onto the molecular system and its complement Q defined by Using these projection operators for the eigenvalue problem of Equation (18), we can reduce it to the eigenvalue problem in the molecular vector space: where L eff (z) is the effective Liouvillian of the molecular system given by In Equation (24), the first term is the unperturbed part and the second term represents the self-energy of the molecular system.Note that the effective Liouvillian depends on its own eigenvalues, so the eigenvalue problem of L (a) eff (z) (23) becomes nonlinear, indicating that the nonperturbative effect of the interaction of the molecular oscillator with the free radiation field has been included.It is this nonlinear feature that makes the eigenvalues of eff (z) coincide with those of L B of the total system.
The denominator of the self-energy exhibits the resonance singularity which is the cause of the dissipation of the molecular excitation energy into the continuous radiation field.By projecting onto the molecular vector space, we obtain the self-energy for the molecular oscillator as a Cauchy integral: where the upper plus indicates the direction of the analytic continuation from the upper half complex plane.As shown in Appendix B, the effective Liouvillian of the molecule is given by By solving the dispersion equation, one obtains the complex eigenvalues z a corresponding to the molecular resonance eigenmode.The right eigenvector of the molecular resonance mode is obtained by Since the submatrix L B is the same as the Friedrichs Hamiltonian [25], the eigenvectors are the same as those eigenvectors belonging to the rigged Hilbert space.Similarly, we have obtained all the other eigenvectors given by Equations (A26)-(A28) in Appendix C. The transformation matrices of Φ B and ΦB are given by right-and left-eigenvectors, respectively, shown in Equations (A29)-(A32).Then, the complex eigenmodes are obtained by the transformation Φ defined in Equation (11): where N a (A24) is the normalization constant, and in Equation ( 33) the suffix d of η + d (z) denotes the delayed analytic continuation such that the contour path is taken so as to avoid the resonance pole in the second Riemann sheet [25].Of essential importance is the direction of the analytic continuations in Equations ( 30)- (33), which ensures the canonical commutation relation between φa and its ⋆-conjugate φ⋆ a [29,39]: where we should emphasize that φj and φ⋆ j are no longer hermitian conjugates of each other when z a is complex.
Since the Heisenberg equation is decoupled in terms of complex eigenmodes Equation (15), it is easily solved, and the solutions are given in Equations (A37)-(A40).
Multiplying Φ from the left of φ(t), we rigorously obtain the solutions of the Heisenberg Equation (5) as In Section 4, we calculate the time-frequency-resolved photon emission under strong coherent driving of a molecular oscillator and clarify the relationship between the scattering and luminescence components.

Time-Frequency-Resolved Photon Emission Spectrum
For a specific calculation, in this paper, we consider the one-dimensional photonic band for the emitted radiation field, which is represented by a semi-infinite tight-binding model.The dispersion of the photonic band is given by where ω B is the center frequency of the photonic band with bandwidth 2B, and a is a lattice constant.We take the coupling constants g k = gB √ 2a/π sin(ka) in this paper, where g denotes a dimensionless coupling strength, which is the same as in Ref. [39].In this paper, we set B = 1 and a = 1 as the energy unit and length unit, respectively.The wavenumber k takes a continuous value in 0 < k < π.The self-energy for the molecular oscillator is obtained analytically as [35] Σ As for the coherent external field, we consider a monochromatic external field with frequency Ω and amplitude f 0 : In the present paper, we show the results for the parameters where the rotating wave approximation is appropriate for the interaction of the molecule with the driving field, although the formulation in Section 3 is valid for the general cases including the counterrotating wave terms.Now, we apply the method presented in the previous Sections to this specific case.In Figure 1, we show the complex eigenvalues of L B , z a , obtained as a solution of Equation ( 27), for the molecular resonance eigenmode as a function of ω 0 , where we take ω B = 10.0, g = 0.2, f 0 = 1.The real and imaginary parts of z a are shown in Figure 1a and Figure 1b, respectively where the solutions on the first and second Riemann sheets are represented by the red and blue curves, respectively.The level scheme for these parameters is shown in Figure 1c.When the molecular oscillator is in resonance with the photonic band for ω B − B ≲ ω 0 ≲ ω B + B, one has complex eigenvalues of z a , indicating that the excited molecular levels decay by emitting a photon into the photonic band.The branch points where the imaginary parts appear are known as the exceptional points.
Before showing the results of the photon emission spectrum, let us study the time evolution of the molecular excitation when the initial state is the unperturbed ground state.We derive the general expression of the molecular excitation with the use of the decomposition by the complex eigenmodes in Equation (A48).Substituting Equation (38) to Equation (A48), one has where, in Equation ( 40), we used the rotating wave approximation assuming the nearresonance situation, i.e., Ω≃ω 0 , ω B .The detail of the derivation is explained in Appendix D.
In Equation ( 40), the first term corresponds to the resonance mode.The second term corresponds to the so-called branch point effect [25,28,34], where Γ denotes the contour of the integral deformed to extract the branch point effect, as shown in Figure A1.The contour deformation to extract the branch point effect is given in Equation (A47) in Appendix D.
The analytic continuation of η Γ (z), defined in Equation ( 27), corresponds to the contour Γ in the Riemann sheets [34].Figure 2 illustrates the time evolution of N m (Ω, t) (solid curve) for parameters ω 0 = 10.5, Ω = 10.0, ω B = 10.0, g = 0.2, and f 0 = 1.In this case, the molecular oscillator is initially in the ground state, with the external field starting to act at t = 0.One observes that N m (Ω, t) oscillates in time with a Rabi frequency (Re[z a ] − Ω) and eventually reaches a stationary value.We evaluate the long time limit of N m (Ω, t).Since the timedependent exponentials of the resonance pole and the branch point effect decay in the long time limit, we evaluate We also show in the inset of Figure 2 the time evolution of the polarization of the molecule defined by This indicates that although the excitation population of the molecule remains constant in the steady state, the polarization oscillates rapidly with the external field frequency Ω.Our method effectively captures the time evolution of the entire system from the quantum vacuum to this nonequilibrium steady state under coherent driving.In particular, the contribution of the resonance mode is found to be dominant over the branch point effect, as long as the frequency of the molecular oscillator is far from the band edge.We now shift our focus to the photon emission spectrum.The time-frequencyresolved photon emission spectrum is obtained as follows: Substituting Equation ( 38) to Equation (A50), one has where ρ(ω) is the density of states of the photonic band given by In Equation ( 46), we use the rotating wave approximation.As shown in Appendix D, the contour of the integral C is deformed to extract the contributions of the resonance pole, the real pole, and the branch point effect, which are represented by the first, second, and third terms in Equation (48), respectively.One can see from Equation (A50) that the resonance pole contribution is attributed to the resonance eigenmode, while the real pole and the branch point contributions are attributed to the continuous eigenmode.This decomposition in terms of the complex eigenmodes gives a new perspective for understanding the photon emission spectrum.
Figure 3 shows the photon emission spectrum N ph (ω k , Ω, t) at t = 100, where the system is assumed to have reached the nonequilibrium steady state, as seen from Figure 2. One can see that there are two types of distinct spectral structures as a function of the driving field and the emitted photon energies, i.e., Ω and ω k : one is a luminescence whose peak appears at ω k ≃ ω 0 = 10.5 independent of Ω, and the other is a Rayleigh scattering whose peak is energetically correlated with Ω, i.e., appears at ω k = Ω.Our complex spectral analysis clarifies the different physical origins of the two peak structures.We give the explicit expression of the stationary photon emission spectrum in Equation (A54).The luminescence is attributed to the resonance mode described by the first term of Equation ( 48), which has a Lorentzian spectral shape with peak position and width determined by Re[z a ] and Im[z a ], respectively.One can see that two variables, ω k and Ω, are separated in the spectral function of the resonance mode component, so that the peak due to the resonance mode component appears just at ω k = Rez a ≃ ω 0 independently of Ω.In other words, there is no energetic correlation between the driving field and luminescence.On the other hand, the coherent scattering is due to the real pole term of Equation ( 48), which makes the variables inseparable.The denominator of the first factor of the real pole term is mainly responsible for the Rayleigh scattering peak at ω k = Ω.
However, caution should be made before concluding that coherent scattering always occurs at ω k = Ω, while incoherent luminescence appears at ω k ≃ ω 0 .In Figure 4, we show the time-frequency-resolved photon emission spectrum for a driving frequency Ω = 10.0 different from the molecular resonance frequency ω 0 = 10.5, so that the Rayleigh peak and the luminescence can be seen separately, all other parameters being the same as in Figure 3.We observe that initially a broad spectrum is gradually built up with time, covering almost the entire photonic band (ω B − B = 9.0 < ω k < ω B + B = 11.0).The broad spectrum is then separated into Rayleigh scattering at ω k = Ω and luminescence at ω k ≃ ω 0 during t ≃ 10 ∼ 20.Thereafter, the Rayleigh scattering peak grows proportional to ( , while the luminescence shows an oscillatory behavior.In order to clarify the origin of these spectral structures, we show in Figure 5 the time evolution of the spectral components of Figure 4, where the contributions of the resonance mode (Figure 5a), given by the continuous mode (Figure 5b), given by and the interference (Figure 5c), given by the cross-term between the two components, are shown separately.The resonance mode and continuous mode components, represented by the first and second terms in Equation (47), contribute to the spectrum, while the branch point effect represented by the third term in Equation ( 48), is negligible since ω 0 is far from the photonic band edges.As can be seen in Figure 5a, the resonance mode is responsible for the incoherent luminescence of an excited molecule, which peaks at ω k = Rez a ≃ ω 0 with width Imz a as a function of the emitted photon energy ω k .On the other hand, Rayleigh scattering at ω k = Ω is solely due to the continuous mode component, as shown in Figure 5b.However, it is important to note that the continuous mode also plays a significant role in the spectral structure around ω k ≃ ω 0 .This specific contribution is derived from the resonance characteristic of η + (ω k ), as indicated in the second term of Equation (48).
In order to reveal the characteristic difference between coherent scattering and incoherent luminescence, we show in Figure 6 the spectral functions of the resonance mode Equation (50) and continuous mode Equation (51) for different time t and driving field frequency Ω: Figure 6a-c for the resonance mode and Figure 6d-f for the continuous mode.It is seen that the spectral profile of the resonance mode component remains the same irrespective of variations in t and Ω, manifesting a Lorentzian spectral shape that peaks at ω k = Rez a with a width of Imz a .The resonance mode component is therefore interpreted as incoherent luminescence from the excited molecule.In contrast, the spectral profile of the continuous mode with a peak around ω k ≃ ω 0 varies with t and Ω.These results indicate that the photon emission associated with the continuous mode is coherent scattering, as the emitted photons are correlated with the incoming driving field.This finding challenges the prevailing assumption that photon emission at the molecular excitation energy represents incoherent luminescence from an excited molecule.While the resonance component shown in Figure 5a converges to a constant in a manner consistent with the molecular excitation population shown in Figure 2, the spectral structure of the continuous mode component at the molecular resonance frequency shows an oscillation, as shown in Figure 5b.This oscillation effect contributes to the spectral structure at the molecular resonance.In addition, Figure 5c shows the presence of destructive interference between the resonance mode and the continuous mode components.This interference is responsible for the gradual build-up and broadness of the spectral structure observed during the early stages of the spectral evolution in Figure 4.

Concluding Remarks
We have presented a formulation for resonant photon emission from a molecular oscillator driven by a strong coherent laser field using complex spectral analysis of Liouvillian.In this paper, we considered the periodic shift of the equilibrium position of the oscillator induced by the strong driving field.In our theory, we treat the emitted radiation field quantum mechanically, including the radiation field as a part of the whole system as well as the molecular oscillator.Moreover, the theory extends the vector space to the rigged Hilbert space, which allows us to interpret the irreversible spontaneous emission process in the framework of quantum dynamics.We obtained the analytical solution of the Heisenberg equation in terms of the complex spectral analysis, with which the molecular excitation and the resonant photon emission spectrum were calculated.Thus, our theory succeeds in describing spontaneous emission, which requires both quantum vacuum fluctuation and irreversibility, in a unified way.
In contrast to the earlier studies [21][22][23], which solved a time-dependent Schrödinger equation for the matter prior to taking into account the interaction with the free radiation field, our theory first identifies the complex eigenmodes of the Liouvillian for the lightmatter interacting system which are crucial for the spontaneous emission.Once we obtained the complex eigenmodes, the exact solution of the Heisenberg equation including the external field was immediately obtained without any approximation.Since the calculation is exact, our method continuously covers the weak coupling up to the strong coupling cases.
Our method clarifies the physical origin of the spectral structure by spectral decomposition in terms of eigenmodes.We confirm that the incoherent luminescence and coherent scattering are mainly due to the resonant mode and continuous components, respectively.However, our detailed spectral decomposition of the time-frequency-resolved photon emission shows that simplistic interpretations of the physical origin based on energetic correlations between the incident driving field and the emitted photon frequencies may be misleading.We found that the spectral structure around the molecular excitation is not only due to the resonance mode, but is also significantly influenced by the continuous mode, which differs from the time evolution of the molecular excitation.Furthermore, our study shows that destructive interference between these modes suppresses the spectral intensity in the early stages, leading to the formation of a broadband spectral structure.
In our study, we observed that branch point effects are insignificant in the results of Figures 2-6.This is due to the fact that the spectral density evaluated by ImΣ(ω k ) [40] approaches zero at the photonic band edges, ω k = ω B ± B, a behavior illustrated in Figure 1b.In addition, the molecular frequency ω 0 is significantly distant from these band edges.Nevertheless, the branch point effect becomes significant when ω 0 is close to the band edges, especially when the spectral density has large values at these points.This can happen when the molecular dipole interacts with an infinite one-dimensional photonic band, where the spectral density diverges at the band edges by the Van Hove singularities [34].Since the branch point effect has been identified as a key factor in the non-Markovian relaxation processes of excited states [41], the exaggeration of the non-Markovian effect due to the Van Hove singularities is worth further investigation [42,43].
In the present study, we showed the results for a long pulse case, where the envelope function is constant f 0 .Our method can be applied to arbitrary envelope functions.Specifically, we found that the relative contribution of the resonance mode to that of the continuous modes is increased by shortening the pulse duration, making the incoherent luminescence prominent compared to the Rayleigh scattering.It is expected that the predictions of our theory can be tested experimentally using state-of-the-art pulse manipulation techniques.
Since the molecular oscillator in our study is a linear system, nonlinear phenomena such as high harmonic generation are not observed.Investigating the quantum properties of photon emission in high-harmonic generation through complex spectral analysis is an interesting topic for future research.
where N a is the normalization constant defined by The complex eigenvalues of the resonance state |ϕ a ⟩ are obtained as the solution of the dispersion equation, where the self-energy is defined in Equation ( 25), and its explicit form of the present paper is given by Equation (37).One also has the left-resonance eigenstate as follows: One also obtains the continuous eigenstates where η(z) is defined by Equation (A25) and the suffix d of η + d (z) denotes the delayed analytic continuation such that the contour path is taken so as to avoid the resonance pole in the second Riemann sheet [25].
With the use of these eigenvectors, one can form the right-transformation matrix Φ B as where the second term represents the branch point effect obtained by the difference of the contour integrals extending to infinity on the complex lower half-plane between the first and second Riemann sheets.As a result, we obtain the molecular excitation as N m (t) = Similarly, one can calculate the time-frequency-resolved photon emission spectrum: (A50) Lastly, we give the explicit expression of the resonance photon spectrum in the long time limit.It follows from Equation (47) that the amplitude of the resonant photon emission reads: where the contour of the integral C is shown in Figure A1c.We deform the contour to decompose the resonance pole, real pole, and the branch point contributions, as shown in Figure A1d, where the factors e −iwt of the resonance pole and the branch point contributions decay in the long time limit.Except for the singular point at ω k = Ω of the Rayleigh peak position, the amplitude in the long time limit is given by

Figure 1 .
Figure 1.The real (a) and imaginary (b) parts of the complex eigenvalues z a of L as a function of ω 0 for ω B = 10.0 and g = 0.2.The red and blue curves represent the solutions on the first and second Riemann sheets, respectively.(c): The level schemes for the parameters used.See text for details.

Figure 2 .
Figure 2. Time evolution of N m (Ω, t) for ω 0 = 10.5, Ω = 10.0, ω B = 10.0, g = 0.2, and f 0 = 1.The contribution of the resonance mode is shown by the dashed curve, which is almost overlapped by the total spectrum shown by the solid curve.The inset shows the time evolution of the polarization, p(Ω, t), of the molecule.

Figure 3 .
Figure 3. Photon emission spectrum N ph (ω k , Ω, t) at a long time t = 100.The parameters used are as in Figure 2.

Figure A1 .
Figure A1.Contour deformation of the integrals of Equation (A46) (a,b) and of Equation (A49) (c,d).The black color indicates the contours and the poles that lie in the first Riemann sheet, while the red color indicates those in the second Riemann sheet.