1. Introduction
In classical mechanics, the states of particles are described by the position
q and momentum
p, and in classical statistical mechanics, by the probability density
[
1,
2,
3]. In quantum mechanics, the states of particles are described either by complex wave functions
[
4,
5] or by the density matrices [
6,
7,
8]. Other representations of the states, such as Wigner quasidistribution functions
[
9] and Husimi quasidistribution [
10], as well as Glauber–Sudarshan functions and their generalizations,
s-quasidistributions, are used [
11,
12]. These representations were found because the attempts to establish the classical description of quantum states, using the probability distributions of the position
q and momentum
p, had difficulties associated with the existence of uncertainty relations for the position and momentum.
The evolution of the wave functions of quantum states is described by the Schrödinger equation and the evolution of the density matrices of quantum states is described by the von Neumann equation [
13]. The evolution formalism can be also based on the approach of the path integral [
14].
Recently, the probability representation of quantum states was introduced [
15,
16]; see [
17], where the states are described by the conventional probability distributions, which are connected to other mentioned representations by integral transforms [
18]. To avoid the mentioned difficulties with the Heisenberg uncertainty relations, the probability distribution employed in the approach [
15] depends on only one position. In this representation, all quantum effects can be explained, in view of standard properties of conventional probability theory presented in [
1].
The evolutions for probability distributions were discussed in [
19], and some new relations between the probability distributions connected with calculations of the Franck–Condon principle were obtained in [
20].
The aim of this work is to consider different tomographic probability distributions describing the quantum oscillator states and states of a charged particle moving in the electric field
E, if the initial state at time
is taken to be one of known normalized states, such as the coherent state [
11] or the squeezed state, as well as the energy state of the harmonic oscillator or the quantum oscillator at temperature
T. We obtain explicit expressions for the integrals of motion of a charged particle following [
21,
22,
23]. Using these integrals of motion, we obtain probability distributions describing different states of the charged particle, if the initial state of the particle is prepared in a specific state of a harmonic oscillator, such as the coherent state or excited state.
This paper is organized as follows.
In
Section 2, the approach of symplectic tomography is briefly reviewed. In
Section 3, tomograms of the oscillator states at temperature
T are considered. In
Section 4, the tomograms of the charged particle with different initial states and path integrals are discussed. A generic approach to the evolution of these states described by explicit forms of tomographic probability distributions is studied in
Section 5. The relation of the approach to the Green function of the Schrödinger equation and the path integral is considered in
Section 6. In
Section 7, the evolution of coherent and excited states of oscillator is studied. Entropies of evolving tomographic probability distributions are discussed in
Section 8. Conclusions and prospectives are presented in
Section 9.
2. Symplectic Tomography
First, we consider tomographic probability distribution describing the particle states in classical statistical mechanics. The conventional description of such states is associated with probability distributions
in the phase space, i.e., non-negative functions depending (or not depending) on time
t and on the position
q and momentum
p of the particle. This function satisfies the normalization condition
One can introduce the classical tomographic probability distribution [
24], using the Radon transform [
25],
The inverse Radon transform reads
One has the normalization condition for the non-negative function
; it is
In view of the delta-function properties, the parameter
is the particle position in the reference frames, with rotated and rescaled momentum and position axes in the phase space.
In the case of the quantum particle, we discuss the construction of tomographic probability distribution describing the state of any system, such as an oscillator, with the wave function, where the Hamiltonian reads
here, we assume that the particle mass
m and Planck’s constant
ℏ are such that
, and the oscillator frequency is
.
The Schrödinger evolution equation has the form
and its solution, in the case of
, can be expressed in terms of the Green function
as
We consider the construction of tomographic conditional probability distribution of only one random position X; this position is determined as in the reference frame in the phase space, where the position axis q and the momentum axis p are transformed according first to rescaling and and then to the following rotation: and . This means that the parameter X is the position of the charged particle measured in the reference frame in the phase space obtained via rescaling and rotation of the initial axes q and p determined by real parameters and . Namely introducing the probability distribution of only one random variable X provides the possibility of avoiding the problems connected with the Heisenberg uncertainty relation.
One can easily check that the density operator
can be mapped onto the function
using the operator
; the map is
The inverse transform reads [
15,
17,
26,
27]
In terms of an explicit form of the wave function, relation (
8) converts into
this function is non-negative and, in view of (
8), it is obvious that
These properties mean that the function
, called the symplectic tomogram, is the conditional probability distribution; the conditions are labeled by parameters
and
, determining the reference frame where the random position
X is measured. In view of relations (
8) and (
9), the distribution
(or the state tomogram) contains all information on the density operator (
9). This means that the tomographic probability distribution completely describes the quantum state. Additionally, this is an explicit construction of the probability representation of quantum state
, which was not known in the literature.
We can find the Wigner function of the state [
9]
and check that
The case of parameter
corresponds to the optical tomogram
, which replaces the symplectic tomogram in (
8) and (
9). Namely, the optical tomogram
of photon states is measured in the experiments in quantum optics [
28], and it is used to reconstruct the Wigner function [
29,
30,
31].
As became clear, the optical tomogram can also be employed to describe quantum states, and it is not only an experimental tool to measure the object determining the Wigner function, as it was considered before symplectic tomography was introduced. Comparing Equations (
7) and (
9), we obtain a new relation for tomogram, which can be expressed in terms of the Green function of the Schrödinger equation
this formula provides the possibility of mapping the initial wave function
onto the evolving tomogram
. It can be rewritten as follows:
where
This relation can be used to connect the Feynman path integral approach, which describes the Green function while considering quantum dynamics, with dynamics of the probability distributions determining the system quantum states. To achieve this, one can rewrite the Green function in Equation (
16) in the form of the Feynman path integral, which expresses the Green function in terms of the Feynman path integral [
32].
One can associate the probability distribution
of three random variables with the tomogram, employing the Bayes formula
where the function
is a probability density, e.g., in the form of normal distribution [
33],
Analogously, one can also obtain new properties of Hermite polynomials [
34].
The density operator
of the quantum state with such a function is obtained, in view of the formula
The probability distribution
is determined by the density operator
, using the relation
the function
is non-negative, and it satisfies the normalization condition
Thus, the quantum state is described by the probability distribution of three random variables:
X,
, and
, with
X, being the position, and parameters
and
being random characteristics of reference frames in the particle phase space, where the position
X is measured. We can use different functions
in (
17) and (
18).
Born’s rule for the probability
reads
If the probability distribution
is chosen to be
the probability distribution
takes the form
this means that the probability distribution of one random position reads
and it coincides with the tomographic conditional probability distribution describing the quantum state with the density operator
(
8).
If the probability distribution
is chosen as the function of two random variables
and
, the probability distribution of the position
X reads
3. Tomograms of the Oscillator States at Temperature T
Now, we consider the tomogram of the oscillator state, corresponding to the thermal state with temperature
T, introducing the parameter
and assuming the oscillator parameters and Planck’s constant to be
. The density operators of the thermal state with the Hamiltonian
i.e.,
has, in the position representation, the matrix elements of the Gaussian form
[
32]; they are
The density operator
satisfies the condition
; this means that, for the oscillator at temperature
, we have
In the Gaussian state (
27), the mean values of the position
and the momentum
are equal to zero. Additionally, the dispersion of the position
and the dispersion of the momentum
read
In addition, the correlation parameter
r of the position and the momentum is also equal to zero, due to the equality
Thus, the Schrödinger–Robertson uncertainty relation of the position and momentum [
35,
36]
corresponds to the inequality for the state (
27), which reads
Tomogram of the thermal state
has the form of conditional normal probability distribution
Then, the symplectic tomogram identified with the oscillator state at temperature
T is explicitly described as follows:
For zero temperature,
, the parameter
, and the state with density operator (
27) is the ground state of the oscillator with the tomogram
The oscillator tomographic probability distribution of three random variables
X,
, and
reads
this probability distribution describes the state of the quantum oscillator at temperature
T.
For generic probability distribution
, we have the tomogram
Tomogram of the ground state of harmonic oscillator at temperature
T, with the wave function
, given by the integral
has the form
For the function
, the probability distribution of three random variables
X,
, and
reads
this means that the probability distribution of position
X has the form
In the case of the classical oscillator at temperature
T, with the joint probability density
the symplectic tomogram is
in the explicit form, we have
Employing the relationship
in the limit
, we arrive at the probability density, which is
This means that for the classical state,
, that is permitted for the classical oscillator, since the uncertainty relations for the position and momentum take place for the quantum oscillator and not for the classical oscillator.
The classical tomographic probability distribution of three random variables
X,
, and
, at temperature
T, reads
4. Charge in the Electric Field
The integrals of motion for classical systems such as, for example, a charged particle with the mass
and charge
e moving in the constant electric field
E along the field are obtained, in view of formulas for the position
and momentum
of the charge, which read
where
is the charge acceleration. The initial position
in (
50) and momentum
in (
51) of the charge trajectory in its phase space are the constants of motion, which are functions of time and the position
q and momentum
p of the form
Thus, we can check that for the classical charge
e, the integrals of motion
and
, corresponding to the classical Hamiltonian
, are linear forms of the position
q, momentum
p, and time
t.
For the charged particle with the charge
e and Hamiltonian (
5), where
and the constant parameter
, there exist two integrals of motion
and
Since the Green function
[
22,
23] satisfies the system of equations
with the initial conditions
, it can be found in the Gaussian form. This function can be presented in the limit of integral form
where
is the matrix element of the evolution operator, which is the Green function
.
One can check that the Gaussian solution for the Green function satisfying (
54)–(
56) reads
where the function
is determined by (
56).
Following the consideration of quantum integrals of motion for the systems with the Hamiltonians, which are quadratic forms in the position and momentum, presented in [
21,
22], we can use quantized expressions in (
54)–(
56) with operators
,
,
, and
acting on the wave functions in the system with the charge Hamiltonian
.
In the case of time-dependent electric field E, the parameter a is the corresponding function, which depends on the time. If we find the quantum integrals of motion and , we obtain the system of equations for the Green function of the Schrödinger equation related to the integrals of motion.
5. Evolution of the Density Operator of Quantum System
in the Probability Representation
The state vectors
of the pure quantum states satisfy the Schrödinger Equation (
6). The density operator
of the quantum states satisfies the von Neumann equation for the systems with the Hamiltonian of the form
The density operator
can be found as a solution of this equation, if one knows the evolution operator
satisfying the following equation:
For the known density operator
, the evolution density operator
reads
If the Hamiltonian
is the Hermitian operator, i.e.,
, the evolution operator
is the unitary operator, i.e.,
.
If the Hamiltonian
of the system under consideration does not depend on time, the evolution operator
satisfying Equation (
60) reads
Any operators
of the form
are integrals of motion, i.e.,
and this means that
satisfy the condition
The density operator
satisfying the evolution Equation (
59) is the integral of motion. Additionally, the operators
where
and
are the position and momentum operators, are the integrals of motion. Then, the operators
are the position and momentum operators satisfying, in the Heisenberg representation, the following evolution equations:
For systems with time-independent Hamiltonians, the integrals of motion
(
64) and the operators
have the property
this means that, if one has an explicit form of the integrals of motion
, the corresponding operator, in the Heisenberg representation, is also given as the integral of motion for the negative time
, and vice versa. The discussed properties provide the possibility of finding the evolution of the symplectic tomogram, e.g., of the oscillator, which is determined by the density operator
; it is
In fact, the density operator is the solution of (
59); from this fact follows that, in view of (
65), we arrive at
or
6. The State Evolution of the Oscillator and Charge Probability Distributions
Employing expression (
72), we study the time dependence of tomographic probability distributions of two systems–oscillators and their charge in the electric field.
The state evolution can be presented as the time evolution of the tomographic probability distribution determining the quantum state. This evolution can be obtained using the evolution equation for the probability distribution, which follows from the evolution equation for the density operator (
59). In our study, we employ the possibility of using the evolution equation for the case of Hamiltonians, which are quadratic forms of the position and momentum operators. For such systems, the Heisenberg position and momentum are expressed as linear forms in the position and momentum, and we use this circumstance to explicitly obtain the evolution of the tomogram due to Equation (
72).
First, we consider the oscillator initially prepared in the state of thermal equilibrium at temperature
T, then at time
, it is removed from the potential
and starts to freely move. The evolution of such a state of the oscillator is described by (
72), where the initial density operator is associated with the tomographic probability distribution (
36). For free motion of quantum oscillator with a removed spring,
and one has the equality
Then, one can express the argument of Dirac delta-function as follows:
Thus, expression (
72) describes the evolution of tomographic conditional probability distribution of the oscillator state initially given for time
by (
36), providing at time
t the following result:
At zero temperature
, we have the evolution of the tomogram of the oscillator’s initial ground state, if the potential energy
is removed (say, a spring of the oscillator is suddenly removed at time
) and the probability distribution describing the evolving state of the oscillator is given by normal time-dependent distribution of the form
At the initial large temperatures, we obtain the expression corresponding to the classical behavior of the probability of the oscillator state.
Second, we consider a case in which the initial state of the charged oscillator at temperature
T and time
starts to evolve following the Hamiltonian with potential energy
. Here, equality (
72) converts into
In view of the relations
after replacing the parameters
X,
, and
in Equation (
36) with the above parameters
,
, and
, we arrive at the expression of the tomogram of the time-dependent state of the charged oscillator evolving in the potential
. As a result, tomogram
of normal distribution describing the evolution of thermal state reads
Thus, for the charged oscillator state, we obtain a tomogram in the form of normal distribution with an evolving mean value of the quadrature
X and evolving dispersion of the quadrature.
9. Conclusions
To conclude, we point out the main results of this work.
Throughout the last century, there were many attempts to formulate various notions of quantum states using different concepts of the wave function, density operator, and quasidistribution functions such as the Wigner function; thus, the possibility of describing the notion of quantum state and identifying the state with a conventional probability distribution function was recently realized [
17]. In this paper, we present a short review of this approach based on the existence of an invertible map of the density operator; for example, of the oscillator state onto the tomographic conditional probability distribution function. Such a probability distribution also exists for the classical oscillator, but the difference in the probability distributions is connected to the existence of the uncertainty relations of the position and momentum in the quantum case and the violation of the uncertainty relation in the classical case. These aspects are the reason why the construction of the probability representation of quantum mechanics took so many decades.
We point out that an analogous consideration can be extended to the case of systems with discrete variables such as spin systems. The spin states can be described by probability distributions, and these probability distributions can also be investigated [
37,
38].
The other new important result of this work is introducing entropy of the quantum state, which is conventional Shannon entropy of the introduced probability representation. We also considered quantum statistics and constructed the conventional probability distribution identified with the quantum state at temperature
T. The explicit form of the harmonic oscillator state identified by the probability distribution is given by Equation (
36). Since this distribution is normal probability distribution, we obtained all its characteristics and the results can be extended to the case of many systems.
Additionally, another important aspect of quantum mechanics is the relation of the Green function of the Schrödinger equation to the Feynmann path integral. We also found the relation of the evolution of the probability distribution of the harmonic oscillator to the path integral (the Feynmann path integral) on the example of a harmonic oscillator.
The different kinds of entropic inequalities known in the probability theory can be discussed for the introduced tomographic probability distributions, and we plan to do this in a future publication. Additionally, the probability representation of quantum mechanics provides the possibility of introducing different distributions and relating them to the irreducible representations of symmetry groups [
33] such as unitary
symmetry of a harmonic oscillator, relating the approach to construction of the group representation operators acting in Hilbert spaces of the wave functions in the oscillator space; this will be considered in the next publication, along with the generalization of the elaborated approach to study the properties of quantum field states.