Next Article in Journal / Special Issue
Critical Temperature and Critical Current Enhancement in Arrays of Josephson Junctions: A Ginzburg–Landau Perspective
Previous Article in Journal
Rectified Lorentz Force from Thermal Current Fluctuations
Previous Article in Special Issue
Adiabatically Manipulated Systems Interacting with Spin Baths beyond the Rotating Wave Approximation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
Department of Physics, Osaka Metropolitan University, Gakuen-cho 1-1, Naka-ku, Sakai 599-8531, Japan
2
Nambu Yoichiro Institute of Theoretical and Experimental Physics (NITEP), Osaka Metropolitan University, Sugimoto 3-3-138, Sumiyoshi-ku, Osaka 558-8585, Japan
*
Author to whom correspondence should be addressed.
Physics 2024, 6(2), 579-598; https://doi.org/10.3390/physics6020038
Submission received: 23 December 2023 / Revised: 20 February 2024 / Accepted: 27 February 2024 / Published: 11 April 2024

Abstract

:
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.

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
H ^ ( t ) = H ^ + W ^ ( t ) ,
H ^ = ω 0 a ^ a ^ + 1 2 + d k ω k b ^ k b ^ k + 1 2 + d k g k ( a ^ b ^ k + a ^ b ^ k ) ,
W ^ ( t ) = f ( t ) ( a ^ + a ^ ) ,
where the hat denotes the operator, a ^ ( a ^ ) 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, = 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 b ^ k ( 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 H ^ , 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 W ^ ( t ) (3).
Here, we define a set of Heisenberg operators as a vector:
ξ ^ ( t ) : = ( a ^ ( t ) , { b ^ k ( t ) } ; a ^ ( t ) , { b ^ k ( t ) } ) T ,
where { b ^ k ( 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,
i t ξ ^ ( t ) = L ^ ( t ) ξ ^ ( t ) = L ^ ξ ^ ( t ) + F ( t ) ,
with the initial condition
ξ ^ ( 0 ) = ξ ^ = ( a ^ , { b ^ k } ; a ^ , { b ^ k } ) T .
In Equation (5), L ^ ( t ) and L ^ denote the Liouvillian superoperator defined by the commutation relation with H ^ ( t ) and H ^ , respectively:
L ^ ( t ) · [ H ^ ( t ) , · ] = [ H ^ + W ^ ( t ) , · ] , L ^ · [ H ^ , · ] .
The Liouvillian superoperator L ^ is represented by a block diagonal matrix for the annihilation operator set ( a ^ , { b ^ k } ) and the creation operator set ( a ^ , { b ^ k } ) :
L = L B 0 0 L B ,
where the submatrix L B is represented by a Hermitian matrix,
L B ω 0 { g k } { g k } T diag ( { ω k } ) = L B T .
The last term of Equation (5) is due to the interaction of the molecule with the driving field represented by F ( t ) : = ( f ( t ) , { 0 k } , f ( t ) , { 0 k } ) T .
To solve the Heisenberg Equation (5), we first find invertible transformation matrices Φ ˜ and Φ to diagonalize the time-independent Liouvillian L given by Equation (8):
Φ ˜ L Φ = Z .
Since the total Liouvillian takes the block diagonal form Equation (8), we can write
Φ ˜ = Φ ˜ B 0 0 Φ ˜ B , Φ = Φ B 0 0 Φ B , Z = Z B 0 0 Z B ,
so that
Φ ˜ Φ = I
and
Z B = diag ( z a , { z k } ) ,
where ★-transformation is defined in Equations (A10) and (A11). Below, we show how to construct Φ ˜ , Φ and obtain Z 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:
i t ϕ ^ ( t ) = Z ϕ ^ ( t ) + Φ ˜ F ( t ) ,
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 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 Appendices Appendix A and Appendix C:
Φ ˜ B L B Φ B = Z B , Φ ˜ B Φ B = Φ B Φ ˜ B = I .
The transformations Φ B and Φ ˜ B correspond to the right- and left-eigenvectors of L B , respectively:
L B Φ B = Φ B Z B , Φ ˜ B L B = Z B Φ ˜ B .
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:
L B | φ j = z j | φ j , φ ˜ j | L B = z j φ ˜ j | ( j = a , { k } ) ,
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
| φ j = | a a | φ j + d k | b k b k | φ j , φ ˜ j | = φ ˜ j | a a | + d k φ ˜ j | b k b k | ,
where we introduce the basis set of | a , { | b k } for the representation of the submatrix L B :
a | : = ( 1 , 0 , ) , b k | : = ( 0 , , 1 k , ) , | a : = a | T , | b k : = b k | T .
These eigenvectors satisfy bi-orthonormality and bi-completeness:
φ ˜ a | φ a = 1 , φ ˜ k | φ k = δ ( k k ) , 1 = | φ a φ ˜ a | + d k | φ k φ ˜ k | ,
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
P : = | a a | , Q : = 1 P = d k | b k b k | .
Using these projection operators for the eigenvalue problem of Equation (18), we can reduce it to the eigenvalue problem in the molecular vector space:
L eff ( a ) ( z j ) P | φ j = z j P | φ j , φ ˜ j | P L eff ( a ) ( z j ) = z j φ ˜ j | P ,
where L eff ( z ) is the effective Liouvillian of the molecular system given by
L eff ( a ) ( z ) = P L B P + P L B Q 1 z Q L B Q Q L B P .
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 eff ( a ) ( 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 L eff ( a ) ( 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:
Σ + ( z ) = g k 2 ( z ω k ) + d k ,
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
L eff ( a ) ( z ) = ω 0 + Σ + ( z ) .
By solving the dispersion equation,
η + ( z ) : = z ω 0 Σ + ( z ) = 0 ,
one obtains the complex eigenvalues z a corresponding to the molecular resonance eigenmode. The right eigenvector of the molecular resonance mode is obtained by
| ϕ a = P | ϕ a + Q | ϕ a = 1 + 1 ( z Q L B Q ) z = z a + Q L B P P | ϕ a
= | a + d k g k ( z ω k ) z = z a + | b k a | ϕ a .
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):
φ ^ a = φ ˜ a | a a ^ + d k φ ˜ a | b k b ^ k = N a 1 / 2 a ^ + d k g k ( z ω k ) z = z a + b ^ k ,
φ ^ k = φ ˜ k | a a ^ + d k φ ˜ k | b k b ^ k = b ^ k + g k η ( ω k ) a ^ + d k g k ω k i ϵ ω k b ^ k ,
φ ^ a = a | φ a a ^ + d k b k | φ a b ^ k = N a 1 / 2 a ^ + d k g k ( z ω k ) z = z a + b ^ k ,
φ ^ k = φ ˜ k | a a ^ + d k φ ˜ k | b k b ^ k = b ^ k + g k η d + ( ω k ) a ^ + d k g k ω k + i ϵ ω k b ^ k ,
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]:
[ φ ^ a , φ ^ a ] = 1 , [ φ ^ k , φ ^ k ] = δ ( k k ) , [ φ ^ k , φ ^ k ] = [ φ ^ k , φ ^ k ] = 0 ,
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
ξ ^ ( t ) = Φ exp [ i Z t ] Φ ˜ ξ ^ + ( i ) 0 t Φ exp [ i Z ( t τ ) ] Φ ˜ F ( τ ) d τ .
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
ω k = ω B B cos ( k a ) ,
where ω B is the center frequency of the photonic band with bandwidth 2 B , and a is a lattice constant. We take the coupling constants g k = g B 2 a / π sin ( k a ) 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]
Σ + ( z ) = 2 g 2 ( z ω B ) ( z ω B ) 2 B 2 .
As for the coherent external field, we consider a monochromatic external field with frequency Ω and amplitude f 0 :
f ( t ) = f 0 cos ( Ω t ) ( t 0 ) , 0 ( t < 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 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 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
N m ( Ω , t ) : = 0 | a ^ ( t ) a ^ ( t ) | 0
= f 0 2 4 0 t d τ ( e i Ω τ + e i Ω τ ) a | φ a φ ˜ a | a e i z a ( t τ ) + d k a | φ k φ ˜ k | a e i ω k ( t τ ) 2
f 0 2 4 N a e i z a t e i Ω t z a Ω + 1 2 π i Γ d ω 1 η Γ ( ω ) e i ω t e i Ω t ω Ω 2 ,
where, in Equation (40), we used the rotating wave approximation assuming the near-resonance 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 time-dependent exponentials of the resonance pole and the branch point effect decay in the long time limit, we evaluate
lim t N m ( Ω , t ) f 0 2 4 N a 1 z a Ω + 1 2 π i Γ d ω 1 η Γ ( ω ) 1 ω Ω 2
= f 0 2 4 a | φ a φ ˜ a | a Ω z a + d k a | φ k φ ˜ k | a Ω ω k 2
= f 0 2 4 a | 1 Ω L B | a 2 = f 0 2 4 1 | η + ( Ω ) | 2 .
We also show in the inset of Figure 2 the time evolution of the polarization of the molecule defined by
p ( Ω , t ) : = 1 2 ω 0 0 | ( a ^ ( t ) + a ^ ( t ) ) | 0 .
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
N ph ( ω k , Ω , t ) : = 0 | b ^ k ( t ) b ^ k ( t ) | 0 ρ ( ω k )
= f 0 2 ρ ( ω k ) 4 0 t d τ ( e i Ω τ + e i Ω τ ) 1 2 π i C d w e i w ( t τ ) b k | 1 w L B | a 2
f 0 2 ρ ( ω k ) 4 0 t d τ e i Ω τ 1 2 π i C d w e i w ( t τ ) b k | 1 w L B | a 2
= f 0 2 g k 2 ρ ( ω k ) 4 1 2 π C d w a | 1 w L B | a 1 w ω k e i Ω t e i w t Ω w 2 = f 0 2 g k 2 ρ ( ω k ) 4 N a e i z a t e i Ω t z a Ω 1 z a ω k + e i ω k t e i Ω t ω k Ω 1 η + ( ω k )
+ 1 2 π i Γ d ω e i ω t e i Ω t ω Ω 1 η Γ ( ω ) 1 ω ω k 2 ,
where ρ ( ω ) is the density of states of the photonic band given by
ρ ( ω ) = 1 B 2 ( ω ω B ) 2 .
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 = Re z 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 ( f 0 2 g k 2 ρ ( ω k ) / 4 | η + ( ω k ) | 2 ) t 2 , 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
f 0 2 g k 2 ρ ( ω k ) 4 N a e i z a t e i Ω t z a Ω 1 z a ω k 2 ,
the continuous mode (Figure 5b), given by
f 0 2 g k 2 ρ ( ω k ) 4 e i ω k t e i Ω t ω k Ω 1 η + ( ω k ) 2 ,
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 = Re z a ω 0 with width Im z 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 = Re z a with a width of Im z 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.

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 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.

Author Contributions

Conceptualization, S.T. and K.K.; methodology, S.T. and K.K.; software, M.K.; validation, M.K., S.T. and K.K.; formal analysis, S.T. and M.K.; investigation, M.K. and S.T.; data curation, M.K.; validation, K.K.; writing—original draft preparation, M.K.; writing—review and editing, S.T. and K.K.; visualization, M.K.; supervision, S.T.; project administration, S.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by JSPS KAKENHI Grants No. JP18K03496.

Data Availability Statement

Data are contained within the article.

Acknowledgments

We would like to express our gratitude to Y. Kayanuma for valuable discussions that greatly contributed to this research. Our thanks also go to K. Noba, S. Garmon, T. Petrosky, and J. Inoue for their insightful comments and suggestions.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A. The Symplectic Structure of the System

In this Appendix, we briefly explain the symplectic structure of the present system [39,44,45]. Using the basis ξ ^ given in Equation (6), the field-free Hamiltonian H ^ (2) is expressed by
H ^ = 1 2 ξ ^ T H ξ ^ ,
where H is a representation matrix of the Hamiltonian in terms of the ξ ^ basis space, which is regarded as a basis of a symplectic space. The field-free Liouvillian, L , is represented by
L = J H ,
where
J = 0 I I 0 ,
is a metric of the symplectic space [44,45]. It can be seen from Equations (8) and (9) that L has a symplectic symmetry:
J L J = L T .
It is known that if a matrix has this symmetry, it can be diagonalized by a symplectic transformation [46]:
Φ ˜ L Φ ˜ 1 = Z , Φ 1 L Φ = Z ,
where the transformation of Φ ˜ and Φ are the symplectic matrices satisfying
Φ ˜ T J Φ ˜ = J , Φ T J Φ = J .
The complex eigenmodes are then given by
ϕ ^ = ( φ ^ , φ ^ ) T = ( { φ ^ α } , { φ ^ α } ) T = Φ ˜ ξ ^ .
This symplecticity of the transformation is essential to preserve the Hamiltonian dynamics of the system. This is because
[ φ ^ α , φ ^ α ] = δ α , α , ( α , α ) ,
where δ α , β denotes the Kronecker delta. The complex eigenmode operators still keep the commutation relation, i.e., the dynamical variables remain the canonical pairs [44].
With this transformation, the field-free Hamiltonian is written in the following diagonalized form as
H ^ = 1 2 ξ ^ T H ξ ^ = 1 2 ϕ ^ T Φ T J T Φ Φ ˜ L Φ ϕ ^ = 1 2 ϕ ^ T J T Z ϕ ^ = d α z α φ ^ α φ ^ α + 1 2 .
It follows from Equation (11), Equation (A6) that one finds the relation
Φ ˜ B = Φ B T .
One can prove, similarly,
Φ B = Φ ˜ B T .
In Appendix B, we derive the transformation of Φ ˜ by solving the complex eigenvalue problem of the Liouvillian in terms of the Brillouin–Wigner–Feshback projection method.

Appendix B. Effective Liouvillian

In this Appendix, we derive the effective Liouvilian Equation (26) with use of the Brillouin–Wigner–Feshbach projection method. Here, we consider the right eigenvalue problem of L B in Equation (18). Acting from the left on Equation (18), the projection operators P and Q give the simultaneous equations
P L B ( z ) P | φ j ) + P L B ( z ) Q | φ j ) = z j P | φ j ) ,
Q L B ( z ) P | φ j ) + Q L B ( z ) Q | φ j ) = z j Q | φ j ) .
Rewriting Equation (A13) as
Q | φ j ) = 1 z j Q L B Q Q L B P | φ j ) ,
and substituting it to Equation (A12) leads to the eigenvalue problem of the effective Liouvillian L B :
L eff ( a ) P | φ j = z j P | φ j ,
where L eff ( a ) ( z ) is given by Equation (24). Now, we derive the explicit form of L eff ( a ) ( z ) . Associating with the matrix form of L B , one can represent L B in terms of the orthonormal mode basis of ( | a , { | b k } ) and its conjugate defined by Equation (20):
L B = ω 0 | a a | + ω k | b k b k | d k + g k ( | a b k | + | b k a | ) d k .
With the use of the relations
P L B P = ω 0 | a a | , P L B Q = g k | a k | d k , Q L B P = g k | k k | d k , Q L B Q = ω k | k k | d k ,
the explicit form of the effective Liouvillian in the molecular subspace is immediately obtained as
L eff ( a ) ( z ) = L eff ( a ) ( z ) | a a | ,
where
L eff ( a ) ( z ) = ω 0 + g k 2 z ω k d k = ω 0 + Σ ( z ) .
Since the self-energy is represented by the Cauchy integral, we choose its analytic continuation from the upper half plane, which leads to Equation (25).

Appendix C. The Complex Eigenmodes

In order to obtain Φ ˜ given in Equation (11), we find the transformation of Φ ˜ B and Φ B which diagonalize the subspace Liouvillian L B .
Φ ˜ B L B Φ B = Z B ,
which is rewritten by the right-eigenvalue problem,
L B Φ B = Φ B Z B ,
and the left-eigenvalue problem,
Φ ˜ B L B = Z B Φ ˜ B .
Since the subspace Liouvillian L B has the same structure as the Friedrichs Hamiltonian, one obtains the complex eigenstates of L B . The right eigenstates for the resonance state was obtained [25] by
| φ a = N a 1 / 2 | a + d k g k ( z ω k ) z = z a + | b k ,
where N a is the normalization constant defined by
N a : = 1 d d z Σ + ( z ) | z = z a 1 .
The complex eigenvalues of the resonance state | ϕ a are obtained as the solution of the dispersion equation,
η + ( z ) : = z ω 0 Σ + ( z ) = 0 ,
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:
φ ˜ a | = N a 1 / 2 a | + d k g k ( z ω k ) z = z a + b k | .
One also obtains the continuous eigenstates
| φ k = | b k + g k η d + ( ω k ) | a + d k g k ω k + i ϵ ω k | b k ,
φ ˜ k | = b k | + g k η ( ω k ) a | + d k g k ω k i ϵ ω k b k | ,
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
( Φ B ) a , a = a | φ a , ( Φ B ) a , k = a | φ k ,
( Φ B ) k , a = b k | φ a , ( Φ B ) k , k = b k | φ k ,
and the left-transformation matrix Φ ˜ B as
( Φ ˜ B ) a , a = φ ˜ a | a , ( Φ ˜ B ) a , k = φ ˜ a | b k ,
( Φ ˜ B ) k , a = φ ˜ k | a , ( Φ ˜ B ) k , k = φ ˜ k | b k .
It follows from Equations (11), (A7), (A10), (A11) and (A29)–(A32) that
φ ^ a = φ ˜ a | a a ^ + d k φ ˜ a | b k b ^ k = N a 1 / 2 a ^ + d k g k ( z ω k ) z = z a + b ^ k ,
φ ^ k = φ ˜ k | a a ^ + d k φ ˜ k | b k b ^ k = b ^ k + g k η ( ω k ) a ^ + d k g k ω k i ϵ ω k b ^ k ,
φ ^ a = a | φ a a ^ + d k b k | φ a b ^ k = N a 1 / 2 a ^ + d k g k ( z ω k ) z = z a + b ^ k ,
φ ^ k = a | φ k a ^ + d k b k | φ k b ^ k = b ^ k + g k η d + ( ω k ) a ^ + d k g k ω k + i ϵ ω k b ^ k .

Appendix D. Solution of the Heisenberg Equation and the Spectrum

Since the Heisenberg equation in terms of the complex eigenmodes is decoupled, it is immediately solved, and the solutions are
φ ^ a ( t ) = e i z a t φ ^ a + ( i ) 0 t d τ e i z a ( t τ ) φ ˜ a | a f ( τ ) ,
φ ^ k ( t ) = e i ω k t φ ^ k + ( i ) 0 t d τ e i ω k ( t τ ) φ ˜ k | a f ( τ ) ,
φ ^ a ( t ) = e i z a t φ ^ a + i 0 t d τ e i z a ( t τ ) a | φ a f ( τ ) ,
φ ^ k ( t ) = e i ω k t φ ^ k + i 0 t d τ e i ω k ( t τ ) a | φ k f ( τ ) .
Multiplying Φ from the left-hand side of the solutions (A37)–(A40), one obtains the solutions of the Heisenberg Equation (5), ξ ^ ( t ) = Φ ϕ ^ ( t ) , which read
a ^ ( t ) = a | φ a e i z a t φ ^ a + d k a | φ k e i ω k t φ ^ k
+ ( i ) 0 t d τ f ( τ ) a | φ a φ ˜ a | a e i z a ( t τ ) + d k a | φ k φ ˜ k | a e i ω k ( t τ ) , b ^ k ( t ) = b k | φ a e i z a t φ ^ a + d k b k | φ k e i ω k t φ ^ k
+ ( i ) 0 t d τ f ( τ ) b k | φ a φ ˜ a | a e i z a ( t τ ) + d k b k | φ k φ ˜ k | a e i ω k ( t τ ) , a ^ ( t ) = φ ˜ a | a e i z a t φ ^ a + d k φ ˜ k | a e i ω k t φ ^ k
+ i 0 t d τ f ( τ ) a | φ a φ ˜ a | a e i z a ( t τ ) + d k a | φ k φ ˜ k | a e i ω k ( t τ ) , b ^ k ( t ) = φ ˜ a | b k e i z a t φ ^ a + d k φ ˜ k | b k e i ω k t φ ^ k
+ i 0 t d τ f ( τ ) a | φ a φ ˜ a | b k e i z a ( t τ ) + d k a | φ k φ ˜ k | b k e i ω k ( t τ ) .
The molecular excitation is calculated by the solutions (A41)–(A44):
N m ( t ) = 0 | a ^ ( t ) a ^ ( t ) | 0 = ( i ) 0 t d τ f ( τ ) a | φ a φ ˜ a | a e i z a ( t τ ) + d k a | φ k φ ˜ k | a e i ω k ( t τ ) 2 .
For the curly bracket, we calculate
{ } = a | e i L B ( t τ ) | φ a φ ˜ a | + d k | φ k φ ˜ k | | a = a | e i L B ( t τ ) | a = 1 2 π i C d w e i w ( t τ ) a | 1 w L B | a = 1 2 π i C d w e i w ( t τ ) 1 η + ( w ) .
By deforming the contour of the integral denoted by C, we extract the resonance pole and the branch point contributions as shown in Figure A1a,b:
{ } = e i z a ( t τ ) a | φ a φ ˜ a | a + 1 2 π i Γ d w e i w ( t τ ) 1 η Γ ( w ) ,
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 ) = 0 t d τ f ( τ ) e i z a ( t τ ) a | φ a φ ˜ a | a + 1 2 π i Γ d w e i w ( t τ ) 1 η Γ ( w ) 2 .
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.
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.
Physics 06 00038 g0a1
Similarly, one can calculate the time–frequency-resolved photon emission spectrum:
N k ( t ) : = 0 | b ^ k ( t ) b ^ k ( t ) | 0 = ( i ) 0 t d τ f ( τ ) b k | φ a φ ˜ a | a e i z a ( t τ ) + d k b k | φ k φ ˜ k | a e i ω k ( t τ ) 2 = 0 t d τ f ( τ ) b k | e i L B ( t τ ) | a 2 = 0 t d τ f ( τ ) 1 2 π i C d w e i w ( t τ ) b k | 1 w L B | a 2
The Green’s function is similarly evaluated by the contour deformation shown in Figure A1c,d. As a result, one obtaibs
N k ( t ) = g k 2 0 t d τ f ( τ ) e i z a ( t τ ) a | φ a φ ˜ a | a z a ω k + e i ω k ( t τ ) 1 η + ( ω k ) + 1 2 π i Γ d w e i w ( t τ ) 1 w ω k 1 η Γ ( w ) 2 .
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:
A ph ( ω k , Ω , t ) : = f 0 2 g k 2 4 π C d w a | 1 w L B | a 1 w ω k e i Ω t e i w t Ω w ,
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 i w t 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
lim t A ph ( ω k ( Ω ) , Ω , t ) = i f 0 g k 2 e i Ω t 1 Ω ω k 1 η + ( Ω ) 1 Ω ω k i ϵ 1 η + ( ω k ) .
At ω k = Ω , the amplitude is given by
lim t A ph ( ω k ( = Ω ) , Ω , t ) = π f 0 g k 2 e i Ω t η + ( Ω ) δ ( ω k Ω ) ,
where ϵ is a positive infinitesimal. Combining Equations (A52) and (A53), one obtains spectral amplitude of the resonant photon emission in the long time limit as
lim t A ph ( ω k , Ω , t ) = i f 0 g k 2 e i Ω t 1 Ω ω k + i ϵ 1 η + ( Ω ) 1 Ω ω k i ϵ 1 η + ( ω k ) .

References

  1. Weisskopf, V.; Wigner, E. Berechnung der natürlichen Linienbreite auf Grund der Diracschen Lichttheorie. Z. Phys. 1930, 63, 54–73. [Google Scholar] [CrossRef]
  2. Milonni, P.W. The Quantum Vacuum: An Introduction to Quantum Electrodynamics; Academic Press, Inc.: San Diego, CA, USA, 1994. [Google Scholar] [CrossRef]
  3. Toyozawa, Y. Resonance and relaxation in light scattering. J. Phys. Soc. Jpn. 1976, 41, 400–411. [Google Scholar] [CrossRef]
  4. Kotani, A.; Toyozawa, Y. Resonant light scattering and luminescence in weakly coupled localized electron–phonon system. J. Phys. Soc. Jpn. 1976, 41, 1699–1706. [Google Scholar] [CrossRef]
  5. Kayanuma, Y. Resonant secondary radiation in strongly coupled localized electron–phonon system. J. Phys. Soc. Jpn. 1988, 57, 292–301. [Google Scholar] [CrossRef]
  6. Protopapas, M.; Keitel, C.H.; Knight, P.L. Atomic physics with super-high intensity lasers. Rep. Prog. Phys. 1997, 60, 389–486. [Google Scholar] [CrossRef]
  7. Wehner, M.U.; Ulm, M.H.; Chemla, D.S.; Wegener, M. Coherent control of electron–LO–phonon scattering in bulk GaAs. Phys. Rev. Lett. 1998, 80, 1992–1995. [Google Scholar] [CrossRef]
  8. Chemla, D.S.; Shah, J. Many-body and correlation effects in semiconductors. Nature 2001, 411, 549–557. [Google Scholar] [CrossRef]
  9. Corkum, P.B.; Krausz, F. Attosecond science. Nat. Phys. 2007, 3, 381–387. [Google Scholar] [CrossRef]
  10. Misochko, O.V.; Lebedev, M.V. Coherent lattice dynamics in opaque crystals: Testing the adequacy of two-tensor model. Phys. Rev. B 2016, 94, 184307. [Google Scholar] [CrossRef]
  11. Corkum, P.B. Plasma perspective on strong field multiphoton ionization. Phys. Rev. Lett. 1993, 71, 1994–1997. [Google Scholar] [CrossRef]
  12. Brabec, T.; Krausz, F. Intense few-cycle laser fields: Frontiers of nonlinear optics. Rev. Mod. Phys. 2000, 72, 545–591. [Google Scholar] [CrossRef]
  13. Ghimire, S.; Dichiara, A.D.; Sistrunk, E.; Agostini, P.; Dimauro, L.F.; Reis, D.A. Observation of high-order harmonic generation in a bulk crystal. Nat. Phys. 2011, 7, 138–141. [Google Scholar] [CrossRef]
  14. Schubert, O.; Hohenleutner, M.; Langer, F.; Urbanek, B.; Lange, C.; Huttner, U.; Golde, D.; Meier, T.; Kira, M.; Koch, S.W.; et al. Sub-cycle control of terahertz high-harmonic generation by dynamical Bloch oscillations. Nat. Photon. 2014, 8, 119–123. [Google Scholar] [CrossRef]
  15. Vampa, G.; Hammond, T.J.; Thiré, N.; Schmidt, B.E.; Légaré, F.; McDonald, C.R.; Brabec, T.; Corkum, P.B. Linking high harmonics from gases and solids. Nature 2015, 522, 462–464. [Google Scholar] [CrossRef] [PubMed]
  16. Yoshikawa, N.; Tamaya, T.; Tanaka, K. High-harmonic generation in graphene enhanced by elliptically polarized light excitation. Science 2017, 356, 736–738. [Google Scholar] [CrossRef] [PubMed]
  17. Uchida, K.; Otobe, T.; Mochizuki, T.; Kim, C.; Yoshita, M.; Tanaka, K.; Akiyama, H.; Pfeiffer, L.N.; West, K.W.; Hirori, H. Coherent detection of THz-induced sideband emission from excitons in the nonperturbative regime. Phys. Rev. B 2018, 97, 165122. [Google Scholar] [CrossRef]
  18. Vampa, G.; McDonald, C.R.; Orlando, G.; Klug, D.D.; Corkum, P.B.; Brabec, T. Theoretical analysis of high-harmonic generation in solids. Phys. Rev. Lett. 2014, 113, 073901. [Google Scholar] [CrossRef] [PubMed]
  19. Wu, M.; Ghimire, S.; Reis, D.A.; Schafer, K.J.; Gaarde, M.B. High-harmonic generation from Bloch electrons in solids. Phys. Rev. A 2015, 91, 043839. [Google Scholar] [CrossRef]
  20. Luu, T.T.; Wörner, H.J. High-order harmonic generation in solids: A unifying approach. Phys. Rev. B 2016, 94, 115164. [Google Scholar] [CrossRef]
  21. Lappas, D.G.; Fedorov, M.V.; Eberly, J.H. Spectrum of light scattered by a strongly driven atom. Phys. Rev. A 1993, 47, 1327–1335. [Google Scholar] [CrossRef]
  22. Bogatskaya, A.V.; Volkova, E.A.; Popov, A.M. Spontaneous emission of atoms in a strong laser field. J. Exp. Theor. Phys. 2017, 125, 587–596. [Google Scholar] [CrossRef]
  23. Gorlach, A.; Neufeld, O.; Rivera, N.; Cohen, O.; Kaminer, I. The quantum-optical nature of high harmonic generation. Nat. Commun. 2020, 11, 4598. [Google Scholar] [CrossRef] [PubMed]
  24. Bohm, A.R.; Scurek, R.; Wikramasekara, S. Resonances, Gamow vectors and time asymmetric quantum theory. arXiv 1999, arXiv:nucl-th/9902076. [Google Scholar] [CrossRef]
  25. Petrosky, T.; Prigogine, I.; Tasaki, S. Quantum theory of non-integrable systems. Phys. A Stat. Mech. Its Appl. 1991, 173, 175–242. [Google Scholar] [CrossRef]
  26. Hatano, N.; Nelson, D.R. Localization Transitions in non-Hermitian quantum mechanics. Phys. Rev. Lett. 1996, 77, 570–573. [Google Scholar] [CrossRef] [PubMed]
  27. Bender, C.M.; Boettcher, S. Real spectra in non-Hermitian Hamiltonians having PT symmetry. Phys. Rev. Lett. 1998, 80, 5243–5246. [Google Scholar] [CrossRef]
  28. Petrosky, T.; Ordonez, G.; Prigogine, I. Quantum transitions and nonlocality. Phys. Rev. A 2000, 62, 042106. [Google Scholar] [CrossRef]
  29. Antoniou, I.E.; Gadella, M.; Karpov, E.; Prigogine, I.; Pronko, G. Gamov algebras. Chaos Solitons Fractals 2001, 12, 2757–2775. [Google Scholar] [CrossRef]
  30. De la Madrid, R. The role of the rigged Hilbert space in quantum mechanics. Eur. J. Phys. 2005, 26, 287–312. [Google Scholar] [CrossRef]
  31. Petrosky, T.; Prigogine, I. The Liouville space Eetension of quantum mechanics. In Advances in Chemical Physics. Resonances, Instability, and Irreversibility; Prigogine, I., Rice, S.A., Eds.; John Wiley & Sons, Inc.: New York, NY, USA, 1977; pp. 1–120. [Google Scholar] [CrossRef]
  32. Böhm, A. The Rigged Hilbert Space and Quantum Mechanics. Springer: Berlin/Heidelberg, Germany, 1978. [Google Scholar] [CrossRef]
  33. Gel’fand, I.M.; Vilenkin, N.Y. Generalized Functions. Volume 4: Applications of Harmonic Analysis; Academic Press, Inc.: New York, NY, USA, 1964; Available online: https://archive.org/details/gelfand-vilenkin-generalized-functions-vol-4-applications-of-harmonic-analysis/ (accessed on 22 February 2024).
  34. Tanaka, S.; Garmon, S.; Petrosky, T. Nonanalytic enhancement of the charge transfer from adatom to one-dimensional semiconductor superlattice and optical absorption spectrum. Phys. Rev. B 2006, 73, 115340. [Google Scholar] [CrossRef]
  35. Fukuta, T.; Garmon, S.; Kanki, K.; Noba, K.i.; Tanaka, S. Fano absorption spectrum with the complex spectral analysis. Phys. Rev. A 2017, 96, 052511. [Google Scholar] [CrossRef]
  36. Yamane, H.; Tanaka, S. Ultrafast dynamics of high-harmonic generation in terms of complex Floquet spectral analysis. Symmetry 2018, 10, 313. [Google Scholar] [CrossRef]
  37. Friedrichs, K.O. On the perturbation of continuous spectra. Commun. Pure Appl. Math. 1948, 1, 361–406. [Google Scholar] [CrossRef]
  38. Chruściński, D. Quantum mechanics of damped systems. II. Damping and parabolic potential barrier. J. Math. Phys. 2004, 45, 841–854. [Google Scholar] [CrossRef]
  39. Tanaka, S.; Kanki, K. Dissipative dynamical Casimir effect in terms of complex spectral analysis in the symplectic Floquet space. Prog. Theor. Exp. Phys. 2020, 2020, 12A107. [Google Scholar] [CrossRef]
  40. Weiss, U. Quantum Dissipative Systems; World Scientific Publishing Co. Ltd.: Singapore, 2012. [Google Scholar] [CrossRef]
  41. Khalfin, L.A. Contribution to the decay theory of a quasi-stationary state. Sov. Phys. JETP 1958, 6, 1053–1063. Available online: http://jetp.ras.ru/cgi-bin/e/index/e/6/6/p1053?a=list (accessed on 22 February 2024).
  42. Zhang, W.-M.; Lo, P.-Y.; Xiong, H.-N.; Tu, M.W.-Y.; Nori, F. General non-markovian dynamics of open quantum systems. Phys. Rev. Lett. 2012, 109, 170402. [Google Scholar] [CrossRef]
  43. Ali, M.M.; Zhang, W.-M. Nonequilibrium transient dynamics of photon statistics. Phys. Rev. A 2017, 95, 033830. [Google Scholar] [CrossRef]
  44. Arvind Dutta, B.; Mukunda, N.; Simon, R. The real symplectic groups in quantum mechanics and optics. Pramana 1995, 45, 471–497. [Google Scholar] [CrossRef]
  45. Horoshko, D.B.; La Volpe, L.; Arzani, F.; Treps, N.; Fabre, C.; Kolobov, M.I. Bloch-Messiah reduction for twin beams of light. Phys. Rev. A 2019, 100, 013837. [Google Scholar] [CrossRef]
  46. Meyer, K.R.; Hall, G.R. Introduction to Hamiltonian Dynamical Systems and the N-Body Problem; Springer Science+Business Media: New York, NY, USA, 2013. [Google Scholar] [CrossRef]
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 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.
Physics 06 00038 g001
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 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.
Physics 06 00038 g002
Figure 3. Photon emission spectrum N ph ( ω k , Ω , t ) at a long time t = 100 . The parameters used are as in Figure 2.
Figure 3. Photon emission spectrum N ph ( ω k , Ω , t ) at a long time t = 100 . The parameters used are as in Figure 2.
Physics 06 00038 g003
Figure 4. Time–frequency-resolved photon emission spectrum for Ω = 10.0 . Other parameters used are as in Figure 2.
Figure 4. Time–frequency-resolved photon emission spectrum for Ω = 10.0 . Other parameters used are as in Figure 2.
Physics 06 00038 g004
Figure 5. Spectral components of time–frequency-resolved photon emission spectrum of Figure 4: (a) resonant mode contribution, (b) continuous mode contribution, and (c) interference component. The parameters used are as in Figure 4.
Figure 5. Spectral components of time–frequency-resolved photon emission spectrum of Figure 4: (a) resonant mode contribution, (b) continuous mode contribution, and (c) interference component. The parameters used are as in Figure 4.
Physics 06 00038 g005
Figure 6. Resonance (a,b,c) and continuous (d,e,f) spectral components of time–frequency-resolved photon emission for t = 5 and Ω = 0 (a,d), t = 14 and Ω = 10.0 (b,d), and t = 14 and Ω = 10.7 (c,f). Other parameters used are as in Figure 4.
Figure 6. Resonance (a,b,c) and continuous (d,e,f) spectral components of time–frequency-resolved photon emission for t = 5 and Ω = 0 (a,d), t = 14 and Ω = 10.0 (b,d), and t = 14 and Ω = 10.7 (c,f). Other parameters used are as in Figure 4.
Physics 06 00038 g006
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Katayama, M.; Tanaka, S.; Kanki, K. Ultrafast Resonant Photon Emission from a Molecule Driven by a Strong Coherent Field in Terms of Complex Spectral Analysis. Physics 2024, 6, 579-598. https://doi.org/10.3390/physics6020038

AMA Style

Katayama M, Tanaka S, Kanki K. Ultrafast Resonant Photon Emission from a Molecule Driven by a Strong Coherent Field in Terms of Complex Spectral Analysis. Physics. 2024; 6(2):579-598. https://doi.org/10.3390/physics6020038

Chicago/Turabian Style

Katayama, Maito, Satoshi Tanaka, and Kazuki Kanki. 2024. "Ultrafast Resonant Photon Emission from a Molecule Driven by a Strong Coherent Field in Terms of Complex Spectral Analysis" Physics 6, no. 2: 579-598. https://doi.org/10.3390/physics6020038

APA Style

Katayama, M., Tanaka, S., & Kanki, K. (2024). Ultrafast Resonant Photon Emission from a Molecule Driven by a Strong Coherent Field in Terms of Complex Spectral Analysis. Physics, 6(2), 579-598. https://doi.org/10.3390/physics6020038

Article Metrics

Back to TopTop