Next Article in Journal
The Mathematical Culture in Test Items of National College Entrance Examination in China from 1978 to 2021
Previous Article in Journal
Analysis of Formation Mechanism of Slightly Inclined Bedding Mudstone Landslide in Coal Mining Subsidence Area Based on Finite–Discrete Element Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Koopman Operator and Path Integral of Quantum Free-Electron Laser Model

Department of Physics, Technion, Haifa 32000, Israel
Mathematics 2022, 10(21), 3996; https://doi.org/10.3390/math10213996
Submission received: 18 September 2022 / Revised: 13 October 2022 / Accepted: 22 October 2022 / Published: 27 October 2022

Abstract

:
A quantum model of a free-electron laser (FEL) is considered. Two different approaches for the exploration of the the FEL system are considered. In the first case, the Heisenberg equations of motion are mapped on the basis of the initial wave functions, which consists of the photon coherent states and many-dimensional electron coherent states. This mapping is an exact procedure, which makes it possible to obtain an exact equation of motion for the intensity of the laser field in a closed form. The obtained equation is controlled by a Koopman operator. The analytical expression for the evolution of the FEL intensity is obtained in the framework of a perturbation theory, which is constructed for a small time scale. The second way of the consideration is based on the construction of the many-dimensional path integrals for the evolution of the wave function. This method also makes it possible to estimate the time evolution and the gain of the FEL intensity.

1. Introduction

In this paper, we discuss a quantum model of a free-electron laser (FEL). Experimental implementation and theoretical description of the FEL is a long-lasting problem that started in the seventies of the last century. This extensively studied phenomenon is well described and reviewed in [1,2,3,4], to mention a few. Contemporary studies are also reflected in recent publications and related to both classical and quantum descriptions [5,6], including application of fractional calculus to the FEL model [7,8]. In the present study, our interest to the model relates to a quantum description of the quantum FEL model with the quantum Hamiltonian of an electromagnetic field interacting with a system of N electrons [5,9,10,11]
H ^ = j = 1 N p ^ j 2 2 m + g a ^ j = 1 N e 2 i k z ^ j + a ^ j = 1 N e 2 i k z ^ j .
Here, summation relates to the positions z j and momenta p ^ j of electrons with the mass m and the wave number k. The position–momentum commutation rule is [ z ^ k , p ^ j ] = i δ k , j . The laser mode is described by the photon annihilation and creation operators a ^ and a ^ , respectively, with the commutation rule [ a ^ , a ^ ] = 1 . The interaction parameter g couples the electron dynamics to the photon laser field. The Hamiltonian (1) has been obtained from the consideration of a non-relativistic electron in an electromagnetic field [9,10]; see also ref. [5]. Note also that the starting point of inferring the Hamiltonian (1) is the classical relativistic Hamiltonian [5,10]; see also Appendix A, where some grounds of the applicability of the model are explained.
It is worth stressing that, discussing the dynamics of the system, we also set the notation for the operators considered here. Namely, for any time-dependent operator x ^ ( t ) x ^ t , its initial value is denoted by x ^ = x ^ ( t = 0 ) . The same is concerned with its average values denoting x ( t ) x t and x = x ( t = 0 ) , where x = ( a * , a , { p j , z j } ) . The only exclusion is for the Hamiltonian, since it is the integral of motion and, correspondingly, H ^ ( t ) = H ^ .
The paper relates to mathematical aspects of the quantum mechanical treatment of the system (1) in the framework of two different approaches. The first one is the Koopman operator construction by mapping the Heisenberg equation of motion on the basis of the coherent states. The second approach is the Schrödinger–Feynman consideration of the Green function in the framework of the path integral. The main issue in the task is the laser field intensity. It is shown that both methods lead to the consistent results for the laser field output.
The dynamics of the intensity of the laser field in the operator form A ^ ( t ) = a ^ ( t ) a ^ ( t ) is described by the Heisenberg equation of motion
A ^ ˙ ( t ) = 1 i [ A ^ ( t ) , H ^ ] .
In this approach, our main concern is a technical consideration of the Heisenberg Equation (2) and finding the evolution of the averaged value of the laser field intensity in the form
A ^ ( t ) = e K t A ^ ( t = 0 ) ,
where K is a so-called Koopman operator [12]; see also, e.g., [13,14]. Another approach developed in the paper is the investigation of the evolution of the initial wave function, which relates to the construction of the evolution operator by means of the path integrals [15,16]. For the latter example, the initial wave functions are the direct product of the photon and electron wave functions
| Ψ 0 = | α | z ,
where the coherent states | α are chosen for the photon wave function [17], while the electron wave function is | z = j = 1 N | z j | z 1 , , z N , where z k | z j = δ k , j δ ( z z ) .
The paper is organized as follows. In Section 2, the Koopman operator is constructed by mapping the Heisenberg equations of motion for the intensity of the photon laser field on the basis of the coherent states, both photon and electron. The obtained equation is exact and determines the dynamics of the laser field intensity. The second approach to the system is the description in the framework of the Schrödinger–Feynman consideration, which is presented in Section 3 and Section 4. The path integral is constructed to describe both the evolution of the initial wave function and the intensity of the laser field. The obtained results and the consistency of the two suggested methods are discussed in Section 5. Concluding remarks are presented in Section 6. Appendix A is devoted to some details of the construction and the applicability of the Hamiltonian (1). The path integral construction and its evaluation in the framework of the stationary-phase approximation are presented in Appendix B and Appendix C, respectively.

2. Quantum Evolution according to Koopman Operator

For the quantum mechanical analysis of the Heisenberg equations of motion, we use a technique of mapping the Heisenberg equations on a basis of the coherent states [18,19,20,21,22]. Note that the Hamiltonian is the time-independent integral of motion, H ^ ( t ) = H ^ ( t = 0 ) . Therefore, the Hamiltonian in the Heisenberg equations can be mapped on a basis of the coherent states | α and | β , constructed at the initial moment t = 0 . First, we construct the basis of the photon coherent states. That is, at the initial moment t = 0 , one introduces the coherent states vector | α as the eigenfunction of the annihilation operators a ^ = a ^ ( t = 0 ) , such that a ^ ( t = 0 ) | α a ^ | α = α | α and, correspondingly, α | a ^ = [ a ^ | α ] = α * α | . The coherent state can also be constructed from a vacuum state | 0 as follows [23]
| α = exp [ α a ^ α * a ^ ] | 0 , a ^ | 0 = 0 .
Then, defining the c-number function as an average value of the photon field operator
A ( t ) = α | A ^ ( t ) | α ,
one maps the Heisenberg equation of motion (2) on the basis of the coherent states | α as follows
i d d t A ( t ) = α | A ^ H ^ H ^ A ^ ( t ) | α .
Accounting for Equation (4) and (5), one obtains the mapping rules
α | A ^ ( t ) a ^ | α = e | α | 2 α e | α | 2 A ( t ) ,
α | a ^ A ^ ( t ) | α = e | α | 2 α * e | α | 2 A ( t ) .
As admitted, the Hamiltonian is the integral of motion H ^ a ^ ( t ) , a ^ ( t ) = H ^ ( a ^ , a ^ ) , therefore the mapping rules (7a) and (7b) yield the equation of motion (6) for A ( t ) in the form
d d t A ( t ) = L ^ 1 A ( t ) + L ^ 2 A ,
where the photon–electron operator reads
L ^ 1 = 1 i e | α | 2 H ^ ( α , α ) H ^ ( α * , α * ) e | α | 2 = i g j = 1 N e 2 i k z j α * e 2 i k z j α ,
while the commutator related to the electron part of the Hamiltonian (defined on the test function | z ) reads
L ^ 2 = 1 2 m i [ A ( t ) , p ^ 2 ] .
Note that both operators L ^ 1 and L ^ 2 are Hermitian.
The next step is the averaging procedure for the electron part of the Hilbert space. To this end, we construct a basis of the electron coherent states at the initial moment t = 0 . We admit that the electron system in the Hamiltonian (1) is considered as a system of free spinless particles, and their commutation rules corresponds to a so-called Heisenberg–Weyl group [23]. Therefore, introducing creation b ^ and annihilation b ^ operators:
b ^ j = z j + i p ^ j 2 , b ^ = z j i p ^ j 2 ,
z j = 2 ( b ^ j + b ^ j ) , p ^ j = i 2 ( b ^ j b ^ j ) ,
where [ b ^ j , b ^ j ] = δ j , j , we introduce the coherent state basis as follows:
| β = j = 1 N | β j = j = 1 N exp [ β j b ^ j β j * b ^ j ] | 0 , b ^ j | 0 = 0 ,
which belongs to the Hilbert space of electrons.
Introducing the double average
A A ^ ( t ) = A ^ ( t ) α β
and using properties (7a) and (7b) for both a ^ , a ^ and b ^ , b ^ operators, and taking into account the commutator (10), we obtain the “Koopman equation”
d d t A = K A ,
where the Hermitian Koopman operator reads
K = i 4 m j = 1 N β j * 2 β j 2 + 2 β j β j * 2 β j * β j + 2 β j β j 2 β j * β j *
+ i g j = 1 N α * e 2 i k ( β j * + β j ) e 2 i k β j * α e 2 i k ( β j * + β j ) e 2 i k β j .
The Koopman operator consists of two parts: K = K e + K e p h . The first part K e , described by Equation (15a), corresponds to the kinetic part of free electrons, while the second part describes the electron–photon interaction, as in Equation (15b).
The solution of the Koopman Equation (14) for the initial condition A ( t = 0 ) = α * α is presented in the exponential form
A = e t K α * α = n = 0 t n n ! K n α * α .
We take into account that electrons spend a very short time inside the laser size, namely t 1 , which results from the fact that the laser size is finite and electron velocity is close to the speed of light. Then, we should take into account only a few first terms in the expansion (16), and we limit ourselves by the second order of t 2 . Obviously, for t 0 , the zero-order term corresponds to the initial condition: A ( 0 ) = | α | 2 . Correspondingly, the first-order term for t 1 reads
A ( 1 ) = t K α * α = i g t j = 1 N α * e 2 i k ( β j * + β j ) e 2 i k β j * α e 2 i k ( β j * + β j ) e 2 i k β j α * α = i g t j = 1 N α e 2 i k ( β j * + β j ) α * e 2 i k ( β j * + β j ) .
At this step, the action of the Koopman operator reduces to differentiation with respect to α and α * only. The derivatives β j and β j * do not “work” here. The situation changes for the second-order term, for which we obtain
A ( 2 ) = t 2 K A ( 1 ) = K e A ( 1 ) g 2 t 2 2 × j = 1 N α * e 2 i k ( β j * + β j ) e 2 i k β j * α e 2 i k ( β j * + β j ) e 2 i k β j · j = 1 N α e 2 i k ( β j * + β j ) α * e 2 i k ( β j * + β j ) = t 2 K e A ( 1 ) + e 4 k 2 g 2 t 2 2 j , j = 1 N e 2 i k ( β j * + β j β j * + β j ) .
We restrict the expansion by the second order in Equation (18). The last term is the main contribution to the laser field intensity
A e 4 k 2 g 2 t 2 2 j , j = 1 N e 2 i k ( β j * + β j β j * + β j ) ,
and it is the first main result, which is discussed in Section 5. It is also interesting to admit that in the coherent states, the averaged energy of electrons is conserved that immediately follows from the specific form of the Koopman operator (15) and the transformations (11). Namely, K p j 2 = 0 , then p j 2 ( t ) = p j 2 . This situation is a pure quantum effect. It can be explained by coherent Cherenkov radiation, which is described just by the same Hamiltonian as in Equation (1) [10]; see also Appendix A.

3. Evolution of Initial Wave Function

It is also instructive to obtain an analytical expression for the evolution operator U ( t ) = e i H ^ t / , which determines the evolution of the initial wave function Ψ 0 = ψ 0 ( z ) | α , where ψ 0 can be an arbitrary suitable function. We consider ψ 0 ( z ) = | z defined in Equation (3). The evolution operator can be presented in the form of the path integral related to the electron Hilbert space. In this case, the interaction term in the Hamiltonian (1) plays the role of the potential, which is however the operator valued the function of the operators a ^ and a ^ .
We follow the standard procedure [15,16] partitioning the time interval t = R Δ t at the condition Δ t 0 , while R , where the boundary (initial and finite) conditions for the path integral defined at r = 0 and r = R are z ( r = 0 ) = z and z ( R ) = z ( t ) . The procedure is presented in Appendix B. Then, the path integral arises from this partition as follows (A9):
e i t ( K ^ + V ^ ) / = z ( 0 ) = z z ( t ) = z ( t ) D [ z ( τ ) ] × exp i 0 t d τ j = 1 N m 2 ( z ˙ j ) 2 g j = 1 N e 2 i k z j a ^ + e 2 i k z j a ^
Therefore, the wave function at time t is defined by Equation (20) as follows:
| Ψ ( t , α , z ) = e i H ^ t / ψ 0 ( z ) | α = z z ( t ) D [ z ( τ ) ] × exp i 0 t d τ j = 1 N m 2 ( z ˙ j ( τ ) ) 2 × exp i g 0 t d τ j = 1 N e 2 i k z j ( τ ) a ^ + e 2 i k z j ( τ ) a ^ | α | z .
Here, ψ 0 ( z ) = | z . The last line in Equation (21) describes the shift of the coherent states | α by means of the shift operator D ^ μ ( t ) = e μ ( t ) a ^ μ * ( t ) a ^ , where μ ( t ) = i g 0 t j = 1 N e 2 i k z j ( τ ) d τ . It transforms one coherent state to another [23,24]:
D ^ μ ( t ) | α = e i I m μ ( t ) α * | α + μ ( t ) .
In Equations (20)–(22), the shift operator in the photon Hilbert space is
D ^ ( μ ( t ) ) = e μ ( t ) a ^ μ * ( t ) a ^ = exp i g 0 t d τ j = 1 N e 2 i k z j ( τ ) a ^ + e 2 i k z j ( τ ) a ^ .
Correspondingly, the shifted coherent state in Equation (22) at time t reads
| α ( t , z ) = e i I m μ ( t ) α * | α + μ ( t ) = exp i g 2 0 t d τ j = 1 N α * e 2 i k z j ( τ ) + α e 2 i k z j ( τ ) × α i g 0 t d τ j = 1 N e 2 i k z j ( τ ) .
Note that the wave function in Equation (24) reflects the electron–photon interaction. Substituting Equation (24) in Equation (21), we obtain the wave function in the formal form. Eventually, the exact expression for the evolution of the wave function reads as follows:
| Ψ ( t , α , z ) = z z ( t ) D [ z ( τ ) ] exp i 0 t d τ j = 1 N m 2 ( z ˙ j ( τ ) ) 2 × exp i g 2 0 t d τ j = 1 N α * e 2 i k z j ( τ ) + α e 2 i k z j ( τ ) × α i g 0 t d τ j = 1 N e 2 i k z j ( τ ) | z ,
where the initial wave function is defined in Equation (3). The path integration is performed in the semiclassical or the stationary-phase approximation, defined by the path integral as follows:
G z ; z ( t ) | α + μ ( t ) | z = j = 1 N z j z j ( t ) D [ z j ( τ ) ] exp i j = 1 N × 0 t d τ m 2 ( z ˙ j ( τ ) ) 2 + g | α | cos ϕ α + 2 k z j ( τ ) × α + i g 0 t d τ j = 1 N e 2 i k z j ( τ ) | z ,
where α = | α | e i ϕ α . Note that it is instructive to apply the stationary-phase approximation for the averages, namely for the intensity of the laser field.

4. Evolution of the FEL Intensity

Let us estimate the evolution of the intensity of the laser field a ^ ( t ) a ^ ( t ) , which is
a ^ ( t ) a ^ ( t ) = Ψ ( t , α , z ) | a ^ a ^ | Ψ ( t , α , z ) .
To perform this calculation, we first map the photon operator a ^ ( t ) a ^ ( t ) on the shifted coherent states in Equation (26) that yields the averaged value
[ α * + μ * ( t ) ] [ α + μ ( t ) ] α + μ ( t ) | α + μ ( t ) = [ α * + μ * ( t ) ] [ α + μ ( t ) ] e 1 2 | μ ( t ) μ ( t ) | 2 .
The term with α is estimated by the stationary phase approximation, which is α G [ z e ( τ ) ] , where
G z ; z ( t ) = j = 1 N D [ z j ( τ ) ] e i S [ z j ; z j ( t ) ] G [ z e ( τ ) ] = j = 1 N A j e i S e [ z j , e ] .
The dynamics of z e , j ( τ ) are controlled by a physical pendulum with the principle Hamiltonian function S e ( z j , e ) of a j-th electron and the constant energy E j = p j 2 2 m 2 k | α | cos ( ϕ α + 2 k z j ) , determined at the initial time τ = 0 , while A j is a constant resulted from the Gaussian integration, which corresponds to the second variation δ ( 2 ) S at the stationary point. Note that the principle Hamiltonian function S [ z e , j ( τ ) ] results from the first variation δ S [ z j ; z j ( t ) ] = 0 . Details of the calculation are shown in Appendix C.
The dynamics are described by the elliptic functions with a separatrix, which separates unlimited flight trajectories from the confined elliptic ones. Therefore, the stationary points are determined by the elliptic amplitude function am ( τ , κ ) , [25], where κ = κ ( E j ) 1 is defined in Equation (A23). Therefore, 2 k z j = am ( τ , κ ) ϕ α .
Our main concern, however, is the path integration of the interaction term
i g 0 t d t j = 1 N e 2 i k z j ( t ) .
This part of the FEL intensity is presented in the form of the 2 N dimensional functional integration. Therefore, the main part (M.P.) of the averaged value is defined as follows:
M . P . a ^ ( t ) a ^ ( t ) = g 2 n = 1 N D [ z n ( τ ) ] n = 1 N D [ z n ( τ ) ] j = 1 N j = 1 N 0 t d t 1 0 t d t 2 e 2 i k [ z j ( t 1 ) z j ( t 2 ) ] × exp i j = 1 N 0 t d τ m 2 ( z ˙ j ( τ ) ) 2 g | α | cos ϕ α + 2 k z j ( τ ) × exp i j = 1 N 0 t d τ m 2 ( z ˙ j ( τ ) ) 2 g | α | cos ϕ α + 2 k z j ( τ ) × α i g 0 t d τ j = 1 N e 2 i k z j ( τ ) α i g 0 t d τ j = 1 N e 2 i k z j ( τ ) .
The stationary-phase approximation in Equation (29) is performed according to Equation (28). It is reasonable to take the stationary solution for κ = 1 at the the asymptotic times, according to Figure A1 in Appendix C. Therefore, we obtain that 2 k z j , e ( τ ) = 2 π ϕ α . Eventually, we obtain
a ^ ( t ) a ^ ( t ) g 2 N 2 t 2 j = 1 N | A j | 2 .

5. Discussion

We collect some physical consequences of the calculations of the laser intensity and discuss them in this section. The main contribution to the intensity of the laser photon field is due to the interaction term in the Hamiltonian (1), that is
A g 2 a ^ j = 1 N e 2 i k z j + a ^ j = 1 N e 2 i k z j 2 .
Two different approaches for the evaluation of the laser intensity are suggested in the framework of the Heisenberg and the Schrödinger representations of quantum mechanics.
In the first approach, the Heisenberg equations of motion are mapped on the basis of the coherent states, which are constructed at the initial moment of time t = 0 for both the photon laser field and electrons. Since the algebras of the photon and electron operators belong to the same Heisenberg–Weyl algebras, the analytical forms of the coherent states (as the eigenstates of the annihilation operators) are the same, although the wave functions | α and | β belong to the different Hilbert spaces. The Heisenberg equations of motions for the laser photon field intensity are exactly mapped on the basis of the coherent states and the obtained “Koopman equation”. The work in (14) describes the gain process, and it is controlled by the Koopman operator K . The main advantage of this approach is that the gain intensity can be easily obtained by the small-time perturbation theory, (19), which yields
A e 4 k 2 g 2 t 2 2 j , j = 1 N e 2 i k ( β j * + β j β j * + β j ) .
If initially all coherent states are just the same, that is β j = β , j , then according to Equation (11), the distributions of the initial conditions p j and z j are such that ( p j 2 + z j 2 ) / 2 = | β | 2 j , which is defined as N 2 points on the hyper-sphere in the 2 N -dimensional phase space. In other words, all p j and z j are independent initial values defined on the hyper-sphere, which is just the direct product of N circles for the non-interacting electrons. In this case, the gained intensity is
A g 2 N 2 t 2 ,
which is superradiance. Note that super-radiance supposes that the intensity gain is of the order of N 2 .
The second approach relates to the construction of the path integral for the evolution operator in Equation (20). This makes it possible to estimate the FEL intensity according to the interaction term. This procedure is performed in the stationary-phase approximation. A detailed analysis is presented in Appendix C. The maximum gain effect is estimated for the large-time asymptotics of the stationary-phase solution for E γ t 2 1 . This result is described by Equation (30) and eventually coincides with Equation (32). This corresponds to super-radiance as well.
Note also that the scalar product of the wave function in Equation (25) is
Ψ ( t , α , z ) | Ψ ( t , α , z ) = 1 .
Therefore, in the stationary-phase approximation, the scalar product
α i g 0 t d τ j = 1 N e 2 i k z j ( τ ) α i g 0 t d τ j = 1 N e 2 i k z j ( τ ) = 1
for the stationary-phase solution 2 k z j , e ( τ ) = am ( E γ τ , κ ) ϕ α , which is valid for all physical values of κ and E γ and time τ [ 0 , t ] .

6. Conclusions

A quantum model of a free-electron laser (FEL) is considered. Two different approaches for the evaluation of the laser field intensity are suggested, and these are the Heisenberg and the Schrödinger representations of quantum mechanics. In the first case, the Heisenberg equations of motion are mapped on the many-dimensional basis, which consists of the photon coherent states and many-dimensional electron coherent states. This mapping is an exact procedure, which makes it possible to obtain an equation of motion for the laser intensity in a closed form. The obtained equation is controlled by a Koopman operator. The analytical expression for the evolution of the FEL intensity is obtained in the framework of a perturbation theory, which is constructed for a small time scale.
The second way of the solution is based on the construction of the many-dimensional path integral for the evolution of the wave function. This method also makes it possible to estimate the time evolution and the gain of the FEL intensity.
It should be admitted that the laser field outputs presented in Equations (30) and (32) obtained by these different methods are mathematically consistent. However, the description of the electron dynamics underlying the FEL is completely different. In particular, in the path integral consideration, the laser field intensity is estimated by the stationary-phase approximation, where the electron dynamics are controlled by a physical pendulum and described semiclassically by the elliptic functions. Contrary to that, the Koopman operator describes the electron dynamics exactly, and the latter is a free-particle dynamic. This phenomenon is pure quantum and can be explained by the Cherenkov radiation, which is another scenario of the FEL [10]. It should be admitted that this behavior cannot be explained by the wiggler mechanism of interaction. Citing from the Tamm and Frank’s paper [26], we note that “if one takes in account the fact that an electron moving in a medium does radiate light even if it is moving uniformly provided that its velocity is greater than the velocity of light in the medium”. (See also [27]).

Funding

This research received no external funding.

Conflicts of Interest

The author declares no conflict of interest.

Appendix A. Quantum FEL Hamiltonian

The basic Hamiltonian (1), which governs a variety of different situations has been inferred in refs. [9,10]. We also follow ref. [5]. The starting point is the classical relativistic Hamiltonian
H = [ ( c p e A ) 2 + m 0 2 c 2 ] 1 2 ,
which describes a single electron, with the rest mass m 0 and the elementary charge e interacting with two electromagnetic plane wave fields A = A L e A W . The laser field
A L = A ˜ L ϵ e i k L ( c t z ) + C . C . ,
travels in positive z-direction, while the wiggler field
A W = A ˜ W ϵ e i k W ( c t + z ) + C . C . ,
propagates in the opposite direction. Here, A ˜ L , A ˜ W are amplitudes of the vector potentials, and c is the speed of light while the circular polarization is chosen, when ϵ · ϵ * = 1 and ϵ 2 = ( ϵ * ) 2 = 0 . In a so-called Bambini–Renieri frame [28,29], the wave numbers k L and k W coincide, a.e, k L = k W = k .
The non-relativistic approximation is performed by expanding the relativistic square root and neglecting the relativistic terms that yield [5]
H ( z , p ) = p 2 2 m + e 2 m A ˜ w A ˜ L e 2 i k z + A ˜ L * e 2 i k z ,
where m = m 0 1 + a 0 2 is the shifted mass with the wiggler parameter [5,30]
a 0 = 2 e | A ˜ W | m 0 c .

Appendix A.1. Quantization

In the quantization procedure of the classical Hamiltonian (A5), both the electron dynamics and the laser field are quantized. Therefore, the electron coordinate and the momentum are operators ( z , p ) ( z ^ , p ^ ) , where their commutation relation is [ z ^ , p ^ ] = i . Then, treating the laser field amplitude A ˜ L as the quantized field, the photon annihilation a ^ and creation a ^ d a g operators are introduced, and the commutation relation is [ a ^ , a ^ ] = 1 . Then, the quantized amplitudes of the laser field are defined by the substitutions
A ˜ L A L a ^ A ˜ L * A L a ^ ,
where A L is the amplitude of the quantized laser field. However, the wiggler field is considered as an external classical parameter A ˜ W * = A ˜ W = const [5,10].
Thus, the quantized Hamiltonian reads
H ^ = p ^ 2 2 m + g a ^ e 2 i k z ^ + a ^ e 2 i k z ^ ,
where the coupling constant is
g = e 2 m A L A ˜ W ,
which relates to the wiggler scenario of the FEL. Eventually, the straightforward generalization of the Hamiltonian (A6) for the N electrons leads to the Hamiltonian (1).

Appendix A.2. Cherenkov Radiation Scenario of the FEL

A slightly different consideration has been suggested in ref. [10], which is also suitable for the Cherenkov radiation (CR) scenario of the FEL. In the case of stimulated CR, the wiggler is replaced by a medium with a refractive index n > 1 . Since A W = 0 , the interaction term is due to the transverse canonical momentum of the electron p 0 , while linear polarization ϵ = ϵ * is supposed in Equation (A2). The relativistic Hamiltonian now reads
H = [ c 2 p 2 2 c e ϵ A L + m 0 2 c 4 ] 1 2 .
As shown in ref. [10], quantization of the Hamiltonian (A8) for the N electrons systems yields
H ^ = ω a ^ a ^ + j p j 2 2 m + g j a ^ e i ( k z j + ω t ) + a ^ e i ( k z j + ω t ) ,
where ω = c k / n 2 1 is the frequency of the laser field. The coupling constant now is g = e m c ϵ · p A L , while m 2 = m 0 2 + p 2 ,10].
In the interaction representation
H ^ e i ω a ^ a ^ t H ^ e i ω a ^ a ^ t
and re-scaling k 2 k , one arrives at the Hamiltonian H ^ in Equation (1), obtained by means of the commutation rule [31],
e η b ^ c ^ e η b ^ = n = 0 ( η ) n n ! [ b ^ , [ b ^ [ b ^ , c ^ ] n ] ,
where over-brace contains n operators b ^ . Taking into account that b ^ a ^ a ^ , while c ^ = ( a ^ , a ^ ) , we obtain
e i ω a ^ a ^ t a ^ e i ω a ^ a ^ t = a ^ e i ω t ,
e i ω a ^ a ^ t a ^ d a g e i ω a ^ a ^ t = a ^ d a g e i ω t .

Appendix B. Path Integral

We follow the standard procedure [15,16] partitioning the time interval t = R Δ t at the condition Δ t 0 , while R . Then, the path integral arises from this partition as follows:
e i t ( K ^ + V ^ ) / = lim R e i Δ t ( K ^ + V ^ ) / R = lim R e i Δ t K ^ / e i Δ t V ^ / + O ( Δ t 2 ) R = lim R r = 1 R 1 d N z r r = 0 R 1 z r + 1 | e i Δ t K ^ / e i Δ t V ^ / | z r .
Here, K ^ + V ^ = H ^ is the Hamiltonian (1), such that
K ^ = j = 1 N p ^ j 2 2 m , V ^ = g a ^ j = 1 N e 2 i k z j + a ^ j = 1 N e 2 i k z j .
The operator V ^ leads to shifts in photon Hilbert space, while in the electron Hilbert space it is just the multiplication operator, which is diagonal and behaves like a potential
e i Δ t V ^ / | z r = | z r e i Δ t V ^ ( z r ) / .
Note that it can also be taken at the middle point e i Δ t V ^ z r + z r + 1 2 / that, however, yields the same result in the limit Δ t 0 [15].
To estimate the matrix elements z r + 1 e i t K ^ ^ / | z r , we insert a complete set of momentum states
I = d N p | p p | , p | z = ( 2 π ) 1 / 2 e i p z / .
From Equation (A13), this yields
j z j , r | e i Δ t K ^ ( p ^ j ) / | z j , r + 1 = j d p j z j , r | e i Δ t K ^ ( p ^ j ) | p j , r p j , r | z j , r + 1 ( 2 π ) N j d p j , r exp i Δ t 2 m p j , r 2 e i p j , r ( z j , r + 1 z j , r ) / = ( m / 2 i π Δ t ) N / 2 j exp i m Δ t 2 ( z j , r + 1 z j , r ) Δ t 2
Eventually, we arrive at the expression
e i t ( K ^ + V ^ ) / = lim R ( m / 2 i π Δ t ) R N / 2 × exp i Δ t r = 0 R 1 j = 1 N m 2 ( z j , r + 1 z j , r ) Δ t 2 g j = 1 N e 2 i k z j a ^ + e 2 i k z j a ^ ,
which defines the evolution of the initial state | z according to the Hamiltonian. The boundary condition of the integration for r = 0 is this initial state. Taking the limit R , we arrive at the path integral
e i t ( K ^ + V ^ ) / = z ( 0 ) = z z ( t ) = z ( t ) D ( z ( τ ) ) exp i 0 t d τ j = 1 N m 2 ( z ˙ j ) 2 g j = 1 N e 2 i k z j a ^ + e 2 i k z j a ^

Appendix C. Stationary-Phase Approximation

Let us estimate the integral (26) in the framework of the stationary-phase approximation.
G = j = 1 N z j z j ( t ) D ( z j ( τ ) ) exp i 0 t d τ m 2 ( z ˙ j ( τ ) ) 2 g | α | cos ϕ α + 2 k z j ( τ ) .
where α = | α | e i ϕ α . This approximation yields
G = j = 1 N A j e i S [ ( z j ; z j ( t ) ] j = 1 N A j e i S e ( z j , e ) ,
where the Hamilton principal function S e ( z j , e ) results from the first variation δ S [ ( z j ; z j ( t ) ] = 0 t L [ z ˙ j ( τ ) , z j ( τ ) ] d τ = 0 , which yields d d τ L z ˙ j L z j = 0 . Correspondingly, z e , j ( τ ) is the extreme trajectory determined by the Euler–Lagrange equation
ϕ ˙ j = 2 E + 2 γ 2 cos ϕ ,
where ϕ ϕ ( τ ) = π + ϕ α + 2 k z j ( τ ) , while γ 2 = 4 k 2 g | α | / m and the energy
E j = ( 2 k p j ) 2 / 2 m 2 + γ 2 cos ( ϕ α + 2 k z j ) = 2 ( k p j ) 2 / m 2 + γ 2 cos ( ϕ 0 )
is determined by the initial position z j and the momentum p j of electrons. In the classical limit, without restriction of generality, we can set the energy E j = E = const , j . Therefore, for all j, we obtain the solutions in the form of the elliptic integral as follows:
τ = 1 E γ ϕ 0 ϕ d ϕ 1 κ 2 sin 2 ϕ 2 , κ 2 = 2 γ 2 E + γ 2 = 4 γ 2 E γ ,
where E γ = 2 ( E + γ 2 ) . The solutions are described by the elliptic functions with a separatrix where κ = 1 , which separates unlimited flight trajectories with κ < 1 from the confined elliptic ones with κ > 1 . The physical solutions are for κ 1 , which are amplitude functions am ( E γ ) τ , κ ) , [25]. Therefore,
ϕ ( τ ) = am ( E γ τ , κ ) .
It is instructive to set κ = 1 , which yields Equation (A23) in the form of the topological solution [32] as follows:
E γ τ 2 = ϕ 0 ϕ d ϕ 2 cos ϕ 2 = ln tan π 4 + ϕ 4 ln tan π 4 + ϕ 0 4 .
Then, we have
ϕ ( τ ) = 4 arctan tan π 4 + ϕ 0 4 e E γ τ 2 π .
Figure A1. The solution ϕ ( τ ) in Equation (A26) for E γ τ [ 0 , 10 ] and ϕ 0 = 0 , 0.1 , 0.2 , 0.3 . The linear fitting plots correspond to the initial and the large time asymptotes, which are ϕ ( τ ) = 0.25 + 0.8 E γ τ and ϕ ( τ ) = π .
Figure A1. The solution ϕ ( τ ) in Equation (A26) for E γ τ [ 0 , 10 ] and ϕ 0 = 0 , 0.1 , 0.2 , 0.3 . The linear fitting plots correspond to the initial and the large time asymptotes, which are ϕ ( τ ) = 0.25 + 0.8 E γ τ and ϕ ( τ ) = π .
Mathematics 10 03996 g0a1
A graphical representation of the solution ϕ ( τ ) is presented in Figure A1. As it follows from the plots, the initial solution for E γ τ < 1 is well fitted by the linear approximation, which is z e , j ( τ ) z j + E γ τ / 2 k . In the opposite case, when E γ τ 1 , the solution is a constant value z j , e ( τ ) = π / k , j .

Gaussian Integration

The second step of the approximation is calculation of A j from the second variation δ 2 S , which leads to the Gaussian integration [16,33]. Therefore, accounting for the extremum solution z j , e ( τ ) j , we perform the variable change z j ( τ ) = z j , e ( τ ) + η j ( t a u ) , where η j ( 0 ) = η j ( t ) = 0 . Then, Equation (A20) reads
G e i S e ( z j , e ) j = 1 N e i S e ( z j , e ) 0 , 0 0 , t D η j ( τ ) exp i δ 2 S [ η j ( τ ) ]
with
δ 2 S [ η j ( τ ) ] = 0 t m η ˙ j 2 m γ 2 cos ϕ α + 2 i k z j , e ( τ ) η j 2 ( τ ) d τ .
For the extreme solution, we set the asymptotic one with z j , e ( τ ) = π / k , j . Therefore, δ 2 S [ η j ( τ ) ] = 0 t [ m ( η ˙ j ( τ ) ) 2 + m γ 2 η j 2 ( τ ) ] d τ . It is supposed that η j O for 0 [16] that determines the quantum fluctuations in the vicinity of the classical extreme trajectory. Eventually, we arrive at the dynamics in the inverted potential η j 2 , whose path integral can be easily related to the linear “oscillator” with imaginary γ i γ . In this case, the result of the path integration (for A j ) is well known [15,16] and reads
A j = m γ 2 π i sinh γ t 1 2 .

References

  1. Bonifacio, R.; Casagrande, F.; Cerchioni, G.; De Salvo Souza, L.; Pierini, P.; Piovella, N. Physics of the high-gain FEL and superradiance. Riv. Del Nuovo C. 1990, 13, 2–69. [Google Scholar] [CrossRef]
  2. Colson, W. Classical Free Electron Laser theory. In Laser Handbook; Colson, W., Pellegrini, C., Renieri, A., Eds.; North Holland: Amsterdam, The Netherlands, 1990; Volume VI. [Google Scholar]
  3. Dattoli, G.; Renieri, A.; Torre, A. Lectures on the Free Electron Laser Theory and Related Topics; World Scientific: Singapore, 1993. [Google Scholar] [CrossRef]
  4. Saldin, E.; Schneidmiller, E.; Yurkov, M. The Physics of Free Electron Laser; Springer: Heidelberg/Berlin, Germany, 2000. [Google Scholar]
  5. Kling, P.; Giese, E.; Endrich, R.; Preiss, P.; Sauerbrey, R.; Schleich, W.P. What defines the quantum regime of the free-electron laser? New J. Phys. 2015, 17, 123019. [Google Scholar] [CrossRef]
  6. Piovella, N.; Volpe, L. A review of high-gain free-electron laser theory. Atoms 2021, 9, 28. [Google Scholar] [CrossRef]
  7. Artioli, M.; Dattoli, G.; Licciardi, S.; Pagnutti, S. Fractional Derivatives, Memory Kernels and Solution of a Free Electron Laser Volterra Type Equation. Mathematics 2017, 5, 73. [Google Scholar] [CrossRef] [Green Version]
  8. Iomin, A. Fox H-Functions in Self-Consistent Description of a Free-Electron Laser. Fractal Fract. 2021, 5, 263. [Google Scholar] [CrossRef]
  9. Becker, W.; McIver, J.K. Fully quantized many-particle theory of a free-electron laser. Phys. Rev. A 1983, 27, 1030–1043. [Google Scholar] [CrossRef]
  10. Becker, W.; McIver, J. Many-particle quantum theory for a class of free-electron devices. Phys. Rep. 1987, 154, 205–245. [Google Scholar] [CrossRef]
  11. Kling, P.; Giese, E.; Carmesin, C.M.; Sauerbrey, R.; Schleich, W.P. High-gain quantum free-electron laser: Emergence and exponential gain. Phys. Rev. A 2019, 99, 053823. [Google Scholar] [CrossRef] [Green Version]
  12. Koopman, B.O. Hamiltonian Systems and Transformation in Hilbert Space. Proc. Natl. Acad. Sci. USA 1931, 17, 315–318. [Google Scholar] [CrossRef] [PubMed]
  13. Gaspard, P. Chaos, Scattering and Statistical Mechanics; Cambridge University Press: Cambridge, MA, USA, 1988. [Google Scholar]
  14. Cvitanović, P.; Artuso, R.; Mainieri, R.; Tanner, G.; Vattay, G. Chaos: Classical and Quantum; ChaosBook.org; 2022. [Google Scholar]
  15. Feynman, R.P.; Hibbs, A.R. Quantum Mechanics and Path Integrals; McGraw-Hill: New York, NY, USA, 1965. [Google Scholar]
  16. Schulman, L. Techniques and Applications of Path Integration; Wiley: New York, NY, USA, 1981. [Google Scholar]
  17. Becker, W.; Scully, M.O.; Zubairy, M.S. Generation of Squeezed Coherent States via a Free-Electron Laser. Phys. Rev. Lett. 1982, 48, 475–477. [Google Scholar] [CrossRef]
  18. Sinitsyn, Y.A.; Tsukernik, V.M. Manifestation of quantum properties in evolution of physical quantities for nonlinear Hamiltonian systems. Phys. Lett. A 1982, 90, 339–341. [Google Scholar] [CrossRef]
  19. Berman, G.; Iomin, A.; Zaslavsky, G. Method of quasiclassical approximation for c-number projection in coherent states basis. Phys. D Nonlinear Phenom. 1981, 4, 113–121. [Google Scholar] [CrossRef]
  20. Berman, G.P.; Iomin, A.M.; Kolovsky, A.R.; Tarkhanov, N.N. On dynamics of four-wave interactions in nonlinear quantum chain. Kirensky Inst. Physics. Krasn. 1986. (In Russian) [Google Scholar] [CrossRef]
  21. Iomin, A. Quantum continuous time random walk in nonlinear Schrödinger equation with disorder. Chaos Solitons Fractals 2016, 93, 64–70. [Google Scholar] [CrossRef]
  22. Iomin, A. Subdiffusion in classical and quantum nonlinear Schrödinger equations with disorder. Comput. Math. Appl. 2017, 73, 914–930. [Google Scholar] [CrossRef] [Green Version]
  23. Perelomov, A. Generalized Coherent States and Their Applications; Springer: Berlin, Germany, 1986. [Google Scholar]
  24. Carruthers, P.; Nieto, M.M. Phase and Angle Variables in Quantum Mechanics. Rev. Mod. Phys. 1968, 40, 411–440. [Google Scholar] [CrossRef] [Green Version]
  25. Abramovitz, M.; Stegun, I.A. Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables; Dover Publications: New York, NY, USA, 1972. [Google Scholar]
  26. Frank, I.; Tamm, I. Coherent radiation of fast electrons in a medium. Dokl. Akad. Nauk SSSR 1937, 14, 109–114. [Google Scholar]
  27. Tamm, I.E. Selected Papers; Bolotovskii, B.M., Frenkel, V.Y., Eds.; Springer: Berlin/Heidelberg, Germany, 1991. [Google Scholar]
  28. Bambini, A.; Renieri, A. The Free Electron Laser: A Single-Particle Classical Model. Lett. Nuovo Cimento 1978, 21, 399–404. [Google Scholar] [CrossRef]
  29. Bambini, A.; Renieri, A.; Stenholm, S. Classical theory of the free-electron laser in a moving frame. Phys. Rev. A 1979, 19, 2013–2025. [Google Scholar] [CrossRef]
  30. Schmüser, P.; Dohlus, M.; Rossbach, J. Ultraviolet and Soft X-ray Free-Electron Lasers; Springer: Heidelberg, Germany, 2008. [Google Scholar]
  31. Louisell, W.H. Quantum Statistical Properties of Radiation; Wiley: New York, NY, USA, 1973. [Google Scholar]
  32. Rajaraman, R. Solitons and Instantons; North-Holland Publishing Company: Amsterdam, The Netherlands, 1982. [Google Scholar]
  33. Gel’fand, I.M.; Yaglom, A.M. Integration in function spaces and its application to quantum physics. Uspekhi Mat. Nauk 1956, 11, 77–114. [Google Scholar]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Iomin, A. Koopman Operator and Path Integral of Quantum Free-Electron Laser Model. Mathematics 2022, 10, 3996. https://doi.org/10.3390/math10213996

AMA Style

Iomin A. Koopman Operator and Path Integral of Quantum Free-Electron Laser Model. Mathematics. 2022; 10(21):3996. https://doi.org/10.3390/math10213996

Chicago/Turabian Style

Iomin, Alexander. 2022. "Koopman Operator and Path Integral of Quantum Free-Electron Laser Model" Mathematics 10, no. 21: 3996. https://doi.org/10.3390/math10213996

APA Style

Iomin, A. (2022). Koopman Operator and Path Integral of Quantum Free-Electron Laser Model. Mathematics, 10(21), 3996. https://doi.org/10.3390/math10213996

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop