1. Introduction
Quantum field theory is an important area of physics with various applications [1] [2] [3] [4] . It studies the dynamics of particles of various characteristics including spin, charge and mass.
In the present paper, we intend to give techniques for the study of finite time evolution within quantum field theory (See [5] for earlier attempts). More particularly we consider the case of finite time domain Dirac fields [6] - [11] evolution and study them via path integral methods. We begin from their standard Lagrangian and Hamiltonian densities, expand the fields in terms of anticommuting creation and annihilation operators and path integrate the corresponding Hamiltonian via standard Grassmannian variables techniques [12] . So, on considering vacuum to vacuum transitions we derive the corresponding Green function. It gives the whole dynamical information of the physical system. We use it in the study of Dirac particles’ interaction with an external field [6] [7] [8] , specifically a classical electromagnetic wave. We use perturbative methods with respect to the electromagnetic field and keep first-order terms. From the derived expressions we can study the time evolution of the distribution in space of quantities such as the probability density and therefore of the possible motion of the particles. More particularly here we derive for certain times contour plots of the probability density of a particle initially at rest which has a Gaussian spinor form. We observe that due to the interaction the particle’s probability density spreads and the particle moves away from its initial position.
The present paper proceeds as follows. In section 2 we give the present Dirac field system, expand the field in terms of anticommuting creation and annihilation operators and path integrate the Hamiltonian within a finite time domain to extract the transition amplitude between Grassmann states and more particularly between vacuum states. Then in section 3, we use that result to obtain the Green function. In section 4 we use that Green function in order to perturbatively extract the time evolution of the probability density of a Gaussian Dirac particle interacting with an electromagnetic wave and give appropriate contour plots. Finally, in section 5, we give our conclusions.
In the present paper we set
. Then
, where
is the
fine structure constant. Below we maintain the mass symbol m for clarity and finally in our calculations we set
as well. The electron rest energy is
while the Rydberg energy has value
. The atomic units of length, time
and energy are
,
and
respectively.
2. System and Path Integration
We consider a Dirac field
of mass m coupled to a current X. Its Lagrangian is
(1)
and its corresponding classical action
(2)
are Dirac matrices and
.
We introduce the current in order to use standard properties of the Green function and obtain it after integrating the whole system (see Equation (55) below).
So according to variational considerations, the free field
obeys the equation
(3)
called the Dirac equation.
The momentum canonically conjugate to the spinor field has the form
(4)
and the quantum Hamiltonian is given as
. (5)
Here we study fermions. So
and
must satisfy equal time anticommutation relations,
(6)
. (7)
We expand the fields in creation and annihilation operator components
(8)
(9)
and let the sources have the form
(10)
. (11)
We have set
(12)
(13)
(14)
(15)
and
(16)
where
.
The spinors satisfy the following orthogonality and completeness relations
(17)
(18)
(19)
where
for
and
for
.
In order
and
to satisfy the anticommutation relations (6, 7),
and
must satisfy
(20)
. (21)
Now, on using the above information the Hamiltonian (5) takes the following form
. (22)
We observe that it constitutes of the sum of single mode Hamiltonians of the form
(23)
where
and
are annihilation and creation operators for fermions which satisfy the anticommutation relations (20, 21). They are a reflection of the Pauli principle. We intend to evaluate the propagator. Initially we consider the free case with the following Hamiltonian for one fermion
(24)
We intend to construct a path integral representation of the evolution operator
. To proceed we introduce Grassmann variables
,
within a Grassmann algebra
and represent the Hamiltonian operator
as an operator acting on functions defined on that Grassmann algebra
. Then
,
. (25)
Functions of
constitute a vector space of dimension two.
What motivates to use Grassmann variables is the fact that they obey a similar algebraic structure with the anticommutation structure obeyed by the Fermionic operators
and
. So, we proceed via standard path integral techniques [1] [12] . At first, we suppose that t is small and get
. (26)
Then we use the group property
(27)
to write the evolution operator at finite time
(28)
where
(29)
(30)
and
. (31)
In the large n limit we obtain the path integral representation for the matrix elements of
(32)
with
. (33)
The above expressions immediately generalize to the case of the Hamiltonian (23). Then we obtain the path integral
(34)
and the action
(35)
where
. (36)
The path integral (34) with action (35) is Gaussian and therefore it can be evaluated exactly. The saddle point equation obtained by varying
yields
(37)
with solution
(38)
Similarly, the variation of
gives
(39)
with solution
. (40)
In the above equations we have taken into account boundary conditions.
Finally, we restore the variables appearing in Equation (22) and perform standard manipulations within Equation (35) to obtain the propagator
(41)
To derive that the normalization constant in Equation (41) is one we set
and use the semigroup property.
In order to obtain the final form of the Grassmannian amplitude we have to integrate diagonally between the initial variable
and the final one
. I.e.
(42)
For the similar case of coherent states see [13] [14] .
Now
(43)
and
(44)
So, on evaluating the Grassmannian integral that appears we obtain
(45)
We have set
(46)
Finally
(47)
is the step function.
We can obtain the finite time domain Green function of Dirac fields from the present considerations. We do that in the next section.
3. Green Function
According to the discussion of the previous section in the case of the Hamiltonian (23) the vacuum-to-vacuum amplitude in the presence of sources is
(48)
where
.
In the case of the Hamiltonian (22) we can proceed via similar steps to arrive to a superposition of expressions of the form (48). In our manipulations we have to take into account that
is negative. Therefore the sign of time is modified (see Equation (50) below).
Moreover, we transform the functional dependence on
back to the dependence on
. To achieve that we make use of the orthonormality of the w spinors
(49)
Then we obtain the generating functional
of the correlation functions of the Dirac Lagrangian in the case of a finite time domain in the following form
(50)
Now since
(51)
we can write
(52)
and
(53)
So on using the above equations, Equation (50) becomes
(54)
Therefore, since the Green function is given as
(55)
upon performing the functional derivations, we get
(56)
or
(57)
Further after appropriate geometric series expansions Equations (56, 57) can be written in the form
(58)
We should expect that after the integrations the above series converge and therefore, we can truncate them. Alternatively, we can bypass the poles in Equation (57) if we introduce a constant
which is smaller but close to one and write Equation (57) in the form
(59)
We can use the above expressions in the study of the propagation of Dirac fields in a finite time domain spacetime.
4. Application
Now we apply the above theory to the propagation of a Dirac particle interacting with an external (classical) electromagnetic field characterized by its potential
. Then the action will have the form
(60)
and according to variational considerations
satisfies the equation
(61)
while the corresponding Green function
obeys the equation
(62)
Green functions are important as functionally they are the inverses of the dynamic operators describing the evolution of the various systems. So combined with perturbation theory they constitute a powerful tool. Below we apply that mathematical structure to the present system.
In fact here we intend to consider the time evolution of a Dirac spinor and calculate its probability density as a function of time. Let initially the Dirac particle spinor be
. Then the final spinor
is
(63)
and its four current
has the form
(64)
Now we suppose that
is small enough to be treated as a perturbation so that
may be expressed as an (asymptotic) expansion in
. Then via standard methods [2] we can get the integral equation
(65)
where we have set (see Equation (56))
(66)
So, we suppose that
is the potential of a plane electromagnetic wave propagating along the wavevector
, with frequency
, polarization vector
and amplitude
Then
(67)
and
(68)
Further we set
(69)
On combining Equations (63, 65) we obtain to first order with respect to
the expression
(70)
Then upon using Equation (56) and substituting Equations (67, 68, 69) we get
(71)
We notice that
has dimension of energy (see the comments at the end of the introduction as well).
After appropriate manipulations (see Equation (59) and the comments above that equation) we can derive the following form for the
spinor
(72)
Moreover according to Equation (64) the probability density has the form
(73)
If
as we can easily conclude from Equation (69) the probability density is spherically symmetric. Moreover, if we let the polarization vector
be along the x-axis and the wavevector along the z-axis. I.e.
(74)
and
(75)
We can use the above considerations to extract the probability density
for various times T. In Figure 1 we give several contour plots. We observe that the probability density spreads and that the particle moves away from its initial position.
Figure 1. Contour plots of the probability density of a Dirac particle interacting with a plane electromagnetic wave which is linearly polarized along the x-axis and propagates along the z-axis. We consider contour plots on various planes and for various times T. For
the probability density is spherically symmetric. So, we give one of the identical contour plot and on setting
we suppose that (i, j, k) is a permutation of (1, 2, 3). For larger times we observe that the particle spreads and moves in space. Moreover we have used the values
,
,
,
,
and have applied natural units. I.e. we have set
.
5. Conclusions
In the present paper, we considered the case of a Dirac field in a finite time do-main. We expanded its Hamiltonian in terms of creation and annihilation operators and path integrated via Grassmann variables techniques to extract its Green function. We notice that the finite time domain Green function appears when the boundary conditions correspond to a finite time interval.
Moreover, we applied that Green function in the perturbative study of its interaction with a classical electromagnetic field and considered the time evolution of an initially Gaussian spinor. We kept first order terms with respect the field and as an application we gave contours plots of the probability density on various planes and for various times. From those contour plots we can observe that the probability density of the particle as a function of time spreads and that the particle moves away from its initial position due to the interaction.
In conclusion the present method—structure can be applied to various quantum systems and give interesting results concerning their finite time domain dynamics. In subsequent publications we intend to consider the case of other quantum fields, free or interacting, within finite time domains and study their dynamics.