1. 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 high-order 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 high-order 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.
2. 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
. In this paper, we set the reduced Planck constant,
. 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
(
) are the annihilation and creation boson operators for the radiation field with wavenumber
k and frequency
. 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
. 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
(3).
Here, we define a set of Heisenberg operators as a vector:
where
and
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),
and
denote the Liouvillian superoperator defined by the commutation relation with
and
, respectively:
The Liouvillian superoperator
is represented by a block diagonal matrix for the annihilation operator set
and the creation operator set
:
where the submatrix
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
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
in Appendices
Appendix A and
Appendix 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.
3. 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
given in Equation (
8).
Since the Liouvillian matrix
in Equation (
8) is in block diagonal form, we find the transformation of
and
by solving the complex eigenvalue problem of the inifinite dimensional hermitian submatrix
in Equation (
9), where we note that the matrix structure of
is the same as the Friedrichs model Hamiltonian used for spontaneous emission from a two-level atom [
37]. We diagonalize
with use of the symplectic transform
and
as shown in the Appendices
Appendix A and
Appendix C:
The transformations
and
correspond to the right- and left-eigenvectors of
, respectively:
Since the structure of
is the same as that of the Friedrichs Hamiltonian, we use the results obtained by the Brillouin–Wigner–Feshbach projection method using
instead of the Friedrichs Hamiltonian:
where
and
are the right and left eigenvectors of
with the eigenvalue
, 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
for the representation of the submatrix
:
These eigenvectors satisfy bi-orthonormality and bi-completeness:
where
is the Dirac delta function.
To reduce the eigenvalue problem of the infinite dimensional matrix
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
onto the molecular system and its complement
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
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
(
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
coincide with those of
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
corresponding to the molecular resonance eigenmode. The right eigenvector of the molecular resonance mode is obtained by
Since the submatrix
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
and
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
(
A24) is the normalization constant, and in Equation (33) the suffix
d of
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
and its ★-conjugate
[
29,
39]:
where we should emphasize that
and
are no longer hermitian conjugates of each other when
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
, 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.
4. 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
is the center frequency of the photonic band with bandwidth
, and
a is a lattice constant. We take the coupling constants
in this paper, where
g denotes a dimensionless coupling strength, which is the same as in Ref. [
39]. In this paper, we set
and
as the energy unit and length unit, respectively. The wavenumber
k takes a continuous value in
. 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
:
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 counter-rotating 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
,
, obtained as a solution of Equation (
27), for the molecular resonance eigenmode as a function of
, where we take
,
g = 0.2,
= 1. The real and imaginary parts of
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
, one has complex eigenvalues of
, 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 near-resonance situation, i.e.,
≃
. 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
, defined in Equation (
27), corresponds to the contour
in the Riemann sheets [
34].
Figure 2 illustrates the time evolution of
(solid curve) for parameters
,
,
,
, and
. In this case, the molecular oscillator is initially in the ground state, with the external field starting to act at
. One observes that
oscillates in time with a Rabi frequency (
) and eventually reaches a stationary value. We evaluate the long time limit of
. Since the time-dependent 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–frequency-resolved 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
at
, 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
: one is a luminescence whose peak appears at
independent of
, and the other is a Rayleigh scattering whose peak is energetically correlated with
, i.e., appears at
. 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 and , respectively. One can see that two variables, 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 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 .
However, caution should be made before concluding that coherent scattering always occurs at
, while incoherent luminescence appears at
. In
Figure 4, we show the time–frequency-resolved photon emission spectrum for a driving frequency
different from the molecular resonance frequency
, 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
. The broad spectrum is then separated into Rayleigh scattering at
and luminescence at
during
. 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
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
with width
as a function of the emitted photon energy
. On the other hand, Rayleigh scattering at
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
. This specific contribution is derived from the resonance characteristic of
, 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
with a width of
. 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
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.
5. 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 light–matter 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
Figure 2,
Figure 3,
Figure 4,
Figure 5 and
Figure 6. This is due to the fact that the spectral density evaluated by
[
40] approaches zero at the photonic band edges,
, a behavior illustrated in
Figure 1b. In addition, the molecular frequency
is significantly distant from these band edges. Nevertheless, the branch point effect becomes significant when
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 . 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.