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]
Here, summation relates to the positions
and momenta
of electrons with the mass
m and the wave number
k. The position–momentum commutation rule is
. The laser mode is described by the photon annihilation and creation operators
and
, respectively, with the commutation rule
. 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 , its initial value is denoted by . The same is concerned with its average values denoting and , where . The only exclusion is for the Hamiltonian, since it is the integral of motion and, correspondingly, .
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
is described by the Heisenberg equation of motion
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
where
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
where the coherent states
are chosen for the photon wave function [
17], while the electron wave function is
, where
.
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,
. Therefore, the Hamiltonian in the Heisenberg equations can be mapped on a basis of the coherent states
and
, constructed at the initial moment
. First, we construct the basis of the photon coherent states. That is, at the initial moment
, one introduces the coherent states vector
as the eigenfunction of the annihilation operators
, such that
and, correspondingly,
. The coherent state can also be constructed from a vacuum state
as follows [
23]
Then, defining the
c-number function as an average value of the photon field operator
one maps the Heisenberg equation of motion (
2) on the basis of the coherent states
as follows
Accounting for Equation (
4) and (
5), one obtains the mapping rules
As admitted, the Hamiltonian is the integral of motion
, therefore the mapping rules (
7a) and (7b) yield the equation of motion (
6) for
in the form
where the photon–electron operator reads
while the commutator related to the electron part of the Hamiltonian (defined on the test function
) reads
Note that both operators
and
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
. 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
and annihilation
operators:
where
, we introduce the coherent state basis as follows:
which belongs to the Hilbert space of electrons.
Introducing the double average
and using properties (
7a) and (
7b) for both
and
operators, and taking into account the commutator (
10), we obtain the “Koopman equation”
where the Hermitian Koopman operator reads
The Koopman operator consists of two parts:
. The first part
, 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
is presented in the exponential form
We take into account that electrons spend a very short time inside the laser size, namely
, 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
. Obviously, for
, the zero-order term corresponds to the initial condition:
. Correspondingly, the first-order term for
reads
At this step, the action of the Koopman operator reduces to differentiation with respect to
and
only. The derivatives
and
do not “work” here. The situation changes for the second-order term, for which we obtain
We restrict the expansion by the second order in Equation (
18). The last term is the main contribution to the laser field intensity
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,
, then
. 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
, which determines the evolution of the initial wave function
, where
can be an arbitrary suitable function. We consider
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
and
.
We follow the standard procedure [
15,
16] partitioning the time interval
at the condition
, while
, where the boundary (initial and finite) conditions for the path integral defined at
and
are
and
. The procedure is presented in
Appendix B. Then, the path integral arises from this partition as follows (
A9):
Therefore, the wave function at time
t is defined by Equation (
20) as follows:
Here,
. The last line in Equation (
21) describes the shift of the coherent states
by means of the shift operator
, where
. It transforms one coherent state to another [
23,
24]:
In Equations (
20)–(
22), the shift operator in the photon Hilbert space is
Correspondingly, the shifted coherent state in Equation (
22) at time
t reads
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:
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:
where
. 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
, which is
To perform this calculation, we first map the photon operator
on the shifted coherent states in Equation (
26) that yields the averaged value
The term with
is estimated by the stationary phase approximation, which is
, where
The dynamics of
are controlled by a physical pendulum with the principle Hamiltonian function
of a
j-th electron and the constant energy
, determined at the initial time
, while
is a constant resulted from the Gaussian integration, which corresponds to the second variation
at the stationary point. Note that the principle Hamiltonian function
results from the first variation
. 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
, [
25], where
is defined in Equation (A23). Therefore,
.
Our main concern, however, is the path integration of the interaction term
This part of the FEL intensity is presented in the form of the
dimensional functional integration. Therefore, the main part (M.P.) of the averaged value is defined as follows:
The stationary-phase approximation in Equation (
29) is performed according to Equation (
28). It is reasonable to take the stationary solution for
at the the asymptotic times, according to
Figure A1 in
Appendix C. Therefore, we obtain that
. Eventually, we obtain
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
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
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
. The main advantage of this approach is that the gain intensity can be easily obtained by the small-time perturbation theory, (
19), which yields
If initially all coherent states are just the same, that is
, then according to Equation (11), the distributions of the initial conditions
and
are such that
, which is defined as
points on the hyper-sphere in the
-dimensional phase space. In other words, all
and
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
which is
superradiance. Note that super-radiance supposes that the intensity gain is of the order of
.
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
. 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
Therefore, in the stationary-phase approximation, the scalar product
for the stationary-phase solution
, which is valid for all physical values of
and
and time
.
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]).