Dynamics of System States in the Probability Representation of Quantum Mechanics

Abstract

A short description of the notion of states of quantum systems in terms of conventional probability distribution function is presented. The notion and the structure of entangled probability distributions are clarified. The evolution of even and odd Schrödinger cat states of the inverted oscillator is obtained in the center-of-mass tomographic probability description of the two-mode oscillator. Evolution equations describing the time dependence of probability distributions identified with quantum system states are discussed. The connection with the Schrödinger equation and the von Neumann equation is clarified.

1. Introduction

The motivation of this work is to show, using a simple example of a quantum oscillator and its evolution, that the oscillator states can be described by conventional probability distributions. The dynamics of the oscillator states are described by the time dependence of the probability distribution identified with the quantum state. This description of the quantum system is valid for all the oscillator systems, including the inverted oscillator. The explicit description of this time evolution of the probability distribution, called the center-of-mass tomographic probability distribution of the two-mode inverted oscillator, is one of the goals of our work. Another objective of this paper is to study the structure of probability distribution functions, which describe the entangled states of the quantum system.

The conventional languages of quantum mechanics and classical mechanics are very different from each other. The language of classical mechanics operates with definitions, such as functions, point-wise multiplication, and probabilities. The quantum mechanical language is much more complicated. It operates with such definitions as operators, density matrices, and state vectors. At the dawn of quantum mechanics, in the time of Dirac, it was even called matrix mechanics because of its mathematical apparatus. There are some phenomena and concepts that have passed from classical mechanics to quantum mechanics. The principle of superposition in quantum mechanics appeared as a natural extension of the phenomenon of wave interference in classical mechanics, and Bohr’s quantization rule was inspired by the classical condition for maximum wave interference, which was applied to the de Broglie wave. Some phenomena of quantum mechanics have no classical analogs, such as spin phenomenon or entanglement of states. It would be interesting to enrich the classical theory with some concepts that appear in quantum mechanics, and that were not previously in the classical theory. One such concept is the entangled probability distribution, which we will consider in this paper.

The states of the quantum system in standard formulation of quantum mechanics [1] are determined either by vectors $|\psi \rangle$ in the Hilbert space [2] or by density operators $\widehat{\rho }$ acting in the Hilbert space [3]. The vectors in the Hilbert space are associated with wave functions $\psi (x)$ of pure quantum states, and the density operators are associated with pure or mixed states described by density matrices [4] or matrix elements of the density operators in some representations. Different representations of quantum states were constructed, e.g., Wigner functions $W(q,p)$, which are quasiprobability distributions [5,6] that have some properties of probabilistic distributions. In classical mechanics, the system states are described by objects that are probability distribution functions, and their properties are described by conventional probability theory [7]. The probability theory is also used to study different aspects of quantum system properties [8] as well as in connection with quantum mechanical methods applications to other areas of science [9]. Some new aspects of quantum system correlation properties, such as entanglement phenomena, were discussed in [10,11]. The entanglement phenomenon in quantum physics provides the possibility to apply this notion in classical probability theory [12]. The functions that define the states of a quantum system, and that are probability distribution functions, were introduced in [13]; they were named symplectic tomograms, and this representation was called probability representation of quantum mechanics (see also [14,15,16,17,18]). Some mathematical aspects of the probability representation of quantum and classical states were discussed in [19,20]. The possibility of finding the probability representation of quantum states is based on the existence of the invertable map of the density operators of mixed or pure quantum states, which is mapped onto conventional probability distribution functions. We will discuss this map in our paper, using the example of quantum oscillator states.

The tomograms and the entanglement phenomenon in the two-mode squeezed states and two-mode even and odd coherent states were considered in [21]. Stimulated Raman scattering and stimulated Brillouin scattering of light were considered within the frame of the symplectic tomography scheme in [22,23]. Furthermore, the entanglement phenomenon in the processes of stimulated light scattering of different types, and its connection with the probability distribution functions determining the states of photon and phonon modes, was discussed [22,23,24,25]. In [26], the evolution of different kinds of states in the Kerr medium, including maximum entangled states, were theoretically studied within the frame of the optical tomography scheme (which is a partial case of the symplectic tomography scheme). The instability of the reconstructed tomogram determining the state was considered in [27] in connection with the Radon transform properties. In [28], it was shown that, in classical mechanics, the Hermitian operators can be introduced, and the concepts of classical mechanics can be formulated in a language analogous to the quantum mechanics language. In [29], a review of classical probability representations of quantum states and observables is presented. New fundamental aspects of quantum mechanics based on the groupoid approach are investigated in [30]. In [31], the evolution of states of a system containing quantum and classical parts was studied. The cosmology features were considered within the frame of the probability representation of quantum states in [32,33]. The density matrix properties, using the symplectic representation of quantum mechanics, are given in [34]. In some tomographic methods, the quantization is based on the associative star product of the functions; applications of these approaches in different kinds of experiments were discussed in [35,36,37,38,39,40,41,42,43,44].

The idea to construct the probability representation of quantum states is based on the method of mapping operators onto functions called symbols of operators. This method is the same method that is used to construct the Wigner function [5] and other quasidistributions, such as the Husimi function [45] and the Glauber–Sudarshan P-function [46,47].

The aim of this paper is to study properties of the probability representation and to consider the probability distributions, describing the quantum states in the case of continuous variables. We will consider the dynamics of the quantum oscillator states as the dynamics of the probability distributions, including the superpositions of the wave functions and the superposition principle. Moreover, some examples of the probability distributions for continuous variables (called tomographic probability distributions) will be studied for quantum oscillator systems.

The paper is organized as follows: The notion of entangled probability distributions describing the quantum states in the probability representation of quantum mechanics is discussed in Section 2. A specific example of the entangled probability distribution for a two-mode oscillator is considered in Section 3. The time dependence of states in different representations of quantum mechanics is described in Section 4. The probability representation of quantum states is described in Section 5, using the method of quantizer–dequantizer operators, as well as the evolution equation for the probability distributions and other functions corresponding to quasiprobability representations of system states. The symplectic tomography of oscillator system states is discussed in Section 6, and the dynamics of operator symbols for the Hamiltonians, which are quadratic forms of position and momentum operators, are considered in Section 7. The center-of-mass tomography and dynamics of the Schrödinger cat states of the ordinary and inverted two-mode oscillators, including explicit expressions for time evolution of the center-of-mass tomograms, are presented in Section 8. The conclusions and prospectives of the probability representation of quantum mechanics for studying entanglement and dynamics of the states of quantum systems are presented in Section 9.

2. Entangled Probability Distributions of Random Variables

In order to discuss the entangled probability distribution notion, we will address the concept of the conditional probability notion using an example of a probability distribution of two random variables $P(X,a)$, where X and a are real continuous parameters and the non-negative function satisfies the normalization condition

$$\int P(X,a)dXda=1.$$

(1)

Then, the conditional probability distribution $P(X|a)$ is related to the function $P(X,a)$ by the Bayes formula

$$P(X|a)=\frac{P(X,a)}{\int dXP(X,a)},$$

where a the parameter describing the condition of measuring variable X. The important property of this formula is the normalization condition of the function $P(X|a)$, which reads

$$\int dXP(X|a)=1.$$

(2)

This means that for conditional probability distributions, the integration of the function as the function of the random variable X gives the result that does not depend on the condition parameter a. Analogous properties take place for conditional probability distributions of general random variables, which we will use in our construction of entangled probability distributions in an example of the function of two random variables $X_1$ and $X_2$.

Following [12], we introduce the concept of separable and entangled probability distributions using the notion of entangled states in quantum mechanics, and, as introduced in [13,16,48], the notion of probability representation of quantum states. In this representation of a quantum system, the density operators of separable states can be written as a convex sum of tensor products of the density operators of the subsystems. Using the probability representation of the density operators, we formulate the new notion in the conventional probability theory using, as an example, the probability distribution of two random variables that are obtained using the invertible map of the density operators, which are mapped onto the probability distributions. Definition: the conditional probability distribution $P(X_1,X_2|a_1,a_2)$ is deemed separable if it can be represented as the convex sum of the products of the probability distributions $P^{(k)}(X_1|a_1)$ and $P^{(k)}(X_2|a_2)$ of the form

$$P(X_1,X_2|a_1,a_2)=\sum _k{\mathcal{P}}_k{P}_1^{(k)}(X_1|a_1)P_2^{(k)}(X_2|a_2).$$

(3)

Here, $P(X_1,X_2|a_1,a_2)\ge 0$, $P_1(X_1|a_1)\ge 0$, $P_2(X_2|a_2)\ge 0$, coefficients ${\mathcal{P}}_k\ge 0,{\sum }_k{\mathcal{P}}_k=1$ and

$$\int P(X_1,X_2|a_1,a_2)d{X}_{1}d{X}_2=1.$$

(4)

The probability distribution $P(X_1,X_2|a_1,a_2)$ is called the entangled probability distribution if it cannot be presented as the convex sum of the form (3), i.e.,

$$P(X_1,X_2|a_1,a_2)\ne \sum _k{\mathcal{P}}_k{P}_1^{(k)}(X_1|a_1)P_2^{(k)}(X_2|a_2).$$

(5)

For separable probability distribution

$$\int P(X_1,X_2|a_1,a_2)d{X}_2=\sum _k{\mathcal{P}}_k{P}_1^{(k)}(X_1|a_1)$$

(6)

and

$$\int P(X_1,X_2|a_1,a_2)d{X}_1=\sum _k{\mathcal{P}}_k{P}_2^{(k)}(X_2|a_2).$$

(7)

For the entangled probability distributions, we apply the probability distribution $\mathsf{\Pi }(X_1|a_1)$ as the integral (6)

$$\int P(X_1,X_2|a_1,a_2)d{X}_2=\mathsf{\Pi }(X_1|a_1),$$

(8)

and it cannot be presented as a convex sum as in (6).

The proof of the independence of the integral (8) on the parameter $a_2$ is analogous to the proof of (1), where one uses the Bayes formula for the function of two random variables, $X_1$ and $X_2$, and two conditions described by the parameters $a_1$ and $a_2$, which can be multi-component parameters.

3. Examples of the Entangled Probability Distributions

The entangled probability distribution can be related to probability distributions realized by using the superposition principle of quantum state wave functions; for example, the superposition of Fock states resembling states of two–mode oscillators with the wave functions. We consider the very simple model of state ${\psi }_{+}(x_1,x_2)$ of the form

$${\psi }_{+}(x_1,x_2)=\frac{1}{\sqrt{2}}\left({\psi }_0(x_1){\psi }_1(x_2)+{\psi }_1(x_1){\psi }_0(x_2)\right)=\frac{x_1+x_2}{\sqrt{\pi }}exp\left(-\frac{x_1^2}{2}-\frac{x_2^2}{2}\right).$$

(9)

The function (9) is the superposition of wave functions of two-mode oscillators. The functions $|{\psi }_0(x_1)\rangle$ and $|{\psi }_0(x_2)\rangle$ are ground states of the first and second oscillators, i.e.,

$${\psi }_0(x_1)=\frac{e^{-\frac{x_1^2}{2}}}{{\pi }^{1/4}},{\psi }_0(x_2)=\frac{e^{-\frac{x_2^2}{2}}}{{\pi }^{1/4}},$$

(10)

and the function ${\psi }_1(x_1)$ is the first excited state of the first oscillator, and ${\psi }_1(x_2)$ is the first excited state of the second oscillator, i.e.,

$${\psi }_1(x_1)=\frac{\sqrt{2}{x}_1}{{\pi }^{1/4}}{e}^{-\frac{x_1^2}{2}},{\psi }_1(x_2)=\frac{\sqrt{2}{x}_2}{{\pi }^{1/4}}{e}^{-\frac{x_2^2}{2}}.$$

(11)

One can extend the construction of the superposition function (9) by following the superposition states of the oscillators studied in the literature, which are even and odd states [49]. The state (9) is the even state, i.e., ${\psi }_{+}(-x_1,-x_2)={\psi }_{+}(x_1,x_2)$ and one can consider an odd state, which reads

$${\psi }_{-}(x_1,x_2)=\frac{1}{\sqrt{2}}\left({\psi }_0(x_1){\psi }_1(x_2)-{\psi }_1(x_1){\psi }_0(x_2)\right)=\frac{x_2-x_1}{\sqrt{\pi }}exp\left(-\frac{x_1^2}{2}-\frac{x_2^2}{2}\right).$$

(12)

This state is the odd state, i.e., ${\psi }_{-}(-x_1,-x_2)=-{\psi }_{-}(x_1,x_2)$, which is the analog of odd coherent states (odd Schrödinger cat states). The more general superposition state ${\psi }_b(x_1,x_2)$ can be given by the following construction

$$\begin{array}{ccc}& & {\psi }_b(x_1,x_2)=\frac{1}{\sqrt{2}}\left({\psi }_0(x_1){\psi }_1(x_2)+e^{ib}{\psi }_1(x_1){\psi }_0(x_2)\right)\hfill \\ & & =\frac{1}{\sqrt{\pi }}\left(x_2+e^{ib}{x}_1\right)exp\left(-\frac{x_1^2}{2}-\frac{x_2^2}{2}\right).\hfill \end{array}$$

(13)

An example of a simple separable state is a state that is not a superposition state, for instance,

$${\psi }_s(x_1,x_2)={\psi }_0(x_1){\psi }_1(x_2).$$

(14)

Using the relationship between the symplectic tomogram and the wave function [50]

$$\begin{array}{ccc}& & w(X_1,X_2|{\mu }_1,{\mu }_2,{\nu }_1,{\nu }_2)=\hfill \\ & & \frac{1}{4{\pi }^2|{\nu }_1\left|\right|{\nu }_2|}{\left|\int \psi (x_1,x_2)exp\left(\frac{i{\mu }_1}{2{\nu }_1}{x}_1^2+\frac{i{\mu }_2}{2{\nu }_2}{x}_2^2-\frac{i{X}_1{x}_1}{2{\nu }_1}-\frac{i{X}_2{x}_2}{2{\nu }_2}\right)d{x}_{1}d{x}_2\right|}^2,\hfill \end{array}$$

(15)

one can obtain the explicit form of the conditional probability distribution $w_{+}(X_1,X_2|{\mu }_1,{\nu }_1,{\mu }_2,{\nu }_2)$ for the even state (9), i.e.,

$$\begin{array}{ccc}& & w_{+}(X_1,X_2|{\mu }_1,{\nu }_1,{\mu }_2,{\nu }_2)=\frac{\left({\nu }_2^2+{\mu }_2^2\right)X_1^2+2\left({\nu }_1{\nu }_2+{\mu }_1{\mu }_2\right)X_1{X}_2+\left({\nu }_1^2+{\mu }_1^2\right)X_2^2}{\pi {\left({\nu }_1^2+{\mu }_1^2\right)}^{3/2}{\left({\nu }_2^2+{\mu }_2^2\right)}^{3/2}}\hfill \\ & & \times exp\left[-\frac{X_1^2}{{\mu }_1^2+{\nu }_1^2}-\frac{X_2^2}{{\mu }_2^2+{\nu }_2^2}\right].\hfill \end{array}$$

(16)

For the particular case where ${\nu }_1={\nu }_2=1,{\mu }_1={\mu }_2=0,$ one gets

$$w_{+}(X_1,X_2|{\mu }_1=0,{\nu }_1=1,{\mu }_2=0,{\nu }_2=1)=\frac{1}{\pi }{(X_1+X_2)}^{2}exp\left(-X_1^2-X_2^2\right).$$

(17)

One can check that the function $w_{+}(X_1,X_2|{\mu }_1,{\nu }_1,{\mu }_2,{\nu }_2)$ (17) satisfies the condition

$$\int \int w_{+}(X_1,X_2,|{\mu }_1,{\nu }_1,{\mu }_2,{\nu }_2)d{X}_{1}d{X}_2=1.$$

(18)

As we know, this probability distribution function corresponding to superposition of the wave functions (9) determines the quantum state, which is the entangled state. Due to this, we call this probability distribution the entangled probability distribution. In quantum mechanics, the wave functions of two–mode oscillators, which are obtained by means of superposition of two different wave functions, are entangled pure states. In connection with this, the tomographic probability distribution is described by the probability distribution function (16), and it cannot be represented in the form of Equation (3). On the other hand, it can be represented by the integral

$$w_{+}(X_1|{\mu }_1,{\nu }_1)=\int w_{+}(X_1,X_2|{\mu }_1,{\nu }_1,{\mu }_2,{\nu }_2)d{X}_2=\frac{exp\left(-\frac{X_1^2}{{\mu }_1^2+{\nu }_1^2}\right)}{\sqrt{\pi \left({\mu }_1^2+{\nu }_1^2\right)}}\left[\frac{1}{2}+\frac{X_1^2}{{\mu }_1^2+{\nu }_1^2}\right].$$

(19)

One can check that

$$\int w_{+}(X_1|{\mu }_1,{\nu }_1)d{X}_1=1.$$

(20)

The function $w_{+}(X_1|{\mu }_1,{\nu }_1)$ (19) is a marginal conditional probability distribution of position $X_1$, which is the position of the first oscillator, and the conditions are labeled by the real parameters ${\mu }_1$ and ${\nu }_1$. Moreover, if we repeat analogous calculations for the second oscillator, we get

$$w_{+}(X_2|{\mu }_2,{\nu }_2)=\int w_{+}(X_1,X_2|{\mu }_1,{\nu }_1,{\mu }_2,{\nu }_2)d{X}_1=\frac{exp\left(-\frac{X_2^2}{{\mu }_2^2+{\nu }_2^2}\right)}{\sqrt{\pi \left({\mu }_2^2+{\nu }_2^2\right)}}\left[\frac{1}{2}+\frac{X_2^2}{{\mu }_2^2+{\nu }_2^2}\right].$$

(21)

One can check that

$$\int w_{+}(X_2|{\mu }_2,{\nu }_2)d{X}_2=1.$$

(22)

The function $w_{+}(X_2|{\mu }_2,{\nu }_2)$ is a marginal conditional probability distribution of position $X_2$, which is the position of the second oscillator, and the conditions are labeled by the real parameters ${\mu }_2$ and ${\nu }_2$.

The function (16) is the probability distribution function; it has the form of the sum of three functions that contain products of Gaussian functions and different terms of position products of $X_1$ and $X_2$. The two terms are the probability distribution functions. The third term that is obtained from the integral (15) is not a probability distribution function, but when added to the two terms mentioned above it gives the probability distribution function (tomographic probability distribution).

An analogous procedure can be followed in the case of the odd states, and one gets the explicit form of the conditional probability distribution $w_{-}(X_1,X_2|{\mu }_1,{\nu }_1,{\mu }_2,{\nu }_2)$ for the odd state (12), i.e.,

$$\begin{array}{ccc}& & w_{-}(X_1,X_2|{\mu }_1,{\nu }_1,{\mu }_2,{\nu }_2)=\frac{\left({\nu }_2^2+{\mu }_2^2\right)X_1^2-2\left({\nu }_1{\nu }_2+{\mu }_1{\mu }_2\right)X_1{X}_2+\left({\nu }_1^2+{\mu }_1^2\right)X_2^2}{\pi {\left({\nu }_1^2+{\mu }_1^2\right)}^{3/2}{\left({\nu }_2^2+{\mu }_2^2\right)}^{3/2}}\hfill \\ & & \times exp\left[-\frac{X_1^2}{{\mu }_1^2+{\nu }_1^2}-\frac{X_2^2}{{\mu }_2^2+{\nu }_2^2}\right].\hfill \end{array}$$

(23)

Additionally, in the case of state with the phase (13), we obtain the explicit form of the conditional probability distribution

$$\begin{array}{ccc}& & w_b(X_1,X_2|{\mu }_1,{\nu }_1,{\mu }_2,{\nu }_2)=\hfill \\ & & \frac{\left({\nu }_2^2+{\mu }_2^2\right)X_1^2+2X_1{X}_{2}cos\left[b-arccos\left(({\mu }_1^2+{\nu }_1^2)({\mu }_2^2+{\nu }_2^2)\right)\right]+\left({\nu }_1^2+{\mu }_1^2\right)X_2^2}{\pi {\left({\nu }_1^2+{\mu }_1^2\right)}^{3/2}{\left({\nu }_2^2+{\mu }_2^2\right)}^{3/2}}\hfill \\ & & \times exp\left[-\frac{X_1^2}{{\mu }_1^2+{\nu }_1^2}-\frac{X_2^2}{{\mu }_2^2+{\nu }_2^2}\right].\hfill \end{array}$$

(24)

For the separable state (14), the tomographic probability distribution given by the Formula (15) reads

$$w_{0,1}(X_1,X_2|{\mu }_1,{\nu }_1,{\mu }_2,{\nu }_2)=\frac{2X_2^2}{\sqrt{{\pi }^2({\mu }_1^2+{\nu }_1^2)}{({\mu }_2^2+{\nu }_2^2)}^{3/2}}exp\left(-\frac{X_1^2}{{\mu }_1^2+{\nu }_1^2}-\frac{X_2^2}{{\mu }_2^2+{\nu }_2^2}\right).$$

(25)

The probability distribution (25) has the form of the product of the probability distributions of each mode state $w_0(X_1|{\mu }_1,{\nu }_1)$ and $w_1(X_2|{\mu }_2,{\nu }_2)$.

Our assertion is as follows: The structure of an entangled probability distribution in the general case always has the form of the sum of two terms, namely, the first term is the convex sum of the products of the probability distributions, and the second term is the sum of the products of the terms that are not probabilities. However, when summed, these two terms lead to an entangled probability distribution. This entangled probability distribution has all the properties of the conventional probability distribution.

4. Evolution of States in Different Representations

Let us recount the description of quantum state dynamics in the Hilbert space $\mathcal{H}$, where the pure quantum state is associated with the state vector $|\psi \rangle$[2], and the other states, including the pure states, are also described by the density operators $\widehat{\rho }$ [3,4] acting on the vectors in the Hilbert space $\mathcal{H}$. The dynamics of the states are described by the Schrödinger equation ($\hslash =1$)

$$i\frac{\partial |\psi (t)\rangle }{\partial t}=\widehat{H}|\psi (t)\rangle ,$$

(26)

where $\widehat{H}$ is the system Hermitian Hamiltonian $(\widehat{H}={\widehat{H}}^{†})$.

For the time-independent Hamiltonian, the state vector $|\psi \rangle$ evolves by means of the evolution operator $\widehat{u}(t)=exp(-i\widehat{H}t)$ of the form

$$|\psi (t)\rangle =\widehat{u}(t)|\psi (0)\rangle ,\widehat{u}(0)=\widehat{1}.$$

(27)

For the pure state with the state vector $|\psi (t)\rangle$, the density operator is given by the formula $\widehat{\rho }(t)=|\psi (t)\rangle \langle \psi (t)|$, and the Schrödinger Equation (27) provides the equation for the density operator of the form (the von Neumann equation)

$$i\frac{\partial (|\psi (t)\rangle \langle \psi (t)|)}{\partial t}=\widehat{H}|\psi (t)\rangle \langle \psi (t)|-|\psi (t)\rangle \langle \psi (t)|\widehat{H}.$$

(28)

This equation is also valid for mixed states with the Hermitian density operator $\widehat{\rho }(t)={\sum }_k{\lambda }_k|{\psi }_k(t)\rangle \langle {\psi }_k(t)|$. Here, the parameters ${\lambda }_k$ are probabilities describing mixed states. The equation can be given in the following form

$$\frac{\partial \widehat{\rho }(t)}{\partial t}+i[\widehat{H},\rho (t)]=0.$$

(29)

The solution of this equation, corresponding to the solution of Equation (27), reads

$$\widehat{\rho }(t)=\widehat{u}(t)\widehat{\rho }(0){\widehat{u}}^{†}(t).$$

(30)

The operators such as position $\widehat{q}$ and momentum $\widehat{p}$ operators in the Heisenberg representation, namely, ${\widehat{q}}_H(t)$ and ${\widehat{p}}_H(t)$, are given as follows

$${\widehat{q}}_H(t)={\widehat{u}}^{†}(t)\widehat{q}\widehat{u}(t),{\widehat{p}}_H(t)={\widehat{u}}^{†}(t)\widehat{p}\widehat{u}(t).$$

(31)

The integrals of motion ${\widehat{q}}_0(t)$ and ${\widehat{p}}_0(t)$, which have the initial values ${\widehat{q}}_0(t=0)=\widehat{q}$ and ${\widehat{p}}_0(t=0)=\widehat{p}$ and satisfy Equation (29), are connected with the Heisenberg position and momentum operators for time-independent Hamiltonian by the relationship

$${\widehat{q}}_0(-t)={\widehat{q}}_H(t),{\widehat{p}}_0(-t)={\widehat{p}}_H(t).$$

(32)

The stationary states of a system $|{\psi }_E(t)\rangle$ satisfying the Schrödinger equation (26) have the form

$$|{\psi }_E(t)\rangle =\widehat{u}(t)|{\psi }_E(0)\rangle =exp(-iEt)|{\psi }_E(0)\rangle ,$$

(33)

where the vector $|{\psi }_E(0)\rangle$ is the eigenvector of the Hamiltonian operator, i.e.,

$$\widehat{H}|{\psi }_E(0)\rangle =E|{\psi }_E(0)\rangle .$$

(34)

The eigenvalue parameter E describes the energy level of the system. The superposition principle of quantum states means that the vector $|\psi (t)\rangle$ of the form

$$|\psi (t)\rangle =\sum _k{C}_k|{\psi }_{E_k}(t)\rangle ,$$

(35)

where $C_k$ are complex numbers, is the solution of the Schrödinger equation (26). Furthermore, it means that the density operator $\widehat{\rho }(t)$ of the form

$$\widehat{\rho }(t)=\sum _k\sum _{k^{\prime }}{C}_k{C}_{k^{\prime }}^{*}|{\psi }_{E_k}(t)\rangle \langle {\psi }_{E_{k^{\prime }}}(t)|$$

(36)

is the solution of the von Neumann equation (29). Moreover, it means that due to Equations (33) and (34), we have

$$\widehat{\rho }(t)=\sum _k\sum _{k^{\prime }}{C}_k{C}_{k^{\prime }}^{*}exp\left(i(E_k-E_{k^{\prime }})t\right)|{\psi }_{E_k}(t)\rangle \langle {\psi }_{E_{k^{\prime }}}(t)|,$$

(37)

or

$$\widehat{\rho }(t)=\sum _k|C_k{|}^2|{\psi }_{E_k}(0)\rangle \langle {\psi }_{E_k}(0)|+\sum _k\sum _{k^{\prime }\ne k}{C}_k{C}_{k^{\prime }}^{*}exp\left(i(E_k-E_{k^{\prime }})t\right)|{\psi }_{E_k}(0)\rangle \langle {\psi }_{E_{k^{\prime }}}(0)|.$$

(38)

The dynamics of the state density operator are determined for all the states, which can be represented as superpositions of energy level states by Formula (38) since the vectors $|{\psi }_{E_k}(t)\rangle$ form the complete system of vectors in the Hilbert space $\mathcal{H}$.

5. Quantum States Definition

Now we consider different representations of quantum states using the formalism of quantizer–dequantizer operators $\widehat{D}(\overrightarrow{x})$ and $\widehat{U}(\overrightarrow{x})$[51], where $\overrightarrow{x}$ is a set of parameters $(x_1,x_2,\dots ,x_n)$ such that the density operators $\widehat{\rho }$ can be mapped onto the set of functions $f_{\rho }(\overrightarrow{x})$, which are named symbols of operators, i.e.,

$$f_{\rho }(\overrightarrow{x})=\mathrm{Tr}\widehat{\rho }\widehat{U}(\overrightarrow{x}).$$

(39)

The operator $\widehat{U}(\overrightarrow{x})$ is a dequantizer operator. It maps the operator on its symbol. The density operator can be reconstructed from the symbol of the density operator with the help of the inverse transform

$$\widehat{\rho }=\int f_{\rho }(\overrightarrow{x})\widehat{D}(\overrightarrow{x})d\overrightarrow{x}.$$

(40)

The operator $\widehat{D}(\overrightarrow{x})$ is a quantizer operator. All the state representations, including the Wigner function [5], Husimi function [45], and the Glauber–Sudarshan function [46,47], and corresponding symbols of other operators are formulated using corresponding quantizer–dequantizer operators. Thus, the quantum mechanics can be formulated using the formalism of operators acting in the Hilbert space or their symbols that contain the same information on quantum states. One can transform the quantum mechanics formalism and obtain equations (differential or integral) for the density operator symbols. An important novelty is that the possibility of describing quantum states by conventional probability distributions exists [48,52].

All known functions that are quasiprobability distributions and describe the states of quantum systems can be obtained using various pairs of the quantizer operator $\widehat{D}(\overrightarrow{x})$ and the dequantizer operator $\widehat{U}(\overrightarrow{x})$, where $\overrightarrow{x}=x_1,x_2,\dots x_n$. These operators make it possible to map any operator $\widehat{A}$ acting in a Hilbert space, where the position operator $\widehat{q}$ and momentum operator $\widehat{p}$ act on the function $f_A(\overrightarrow{x})$, which is called the symbol of the operator $\widehat{A}$, using the following general mapping of the operators $\widehat{A}\to f_A(\overrightarrow{x})$ to functions, i.e.,

$$f_A(\overrightarrow{x})=\mathrm{Tr}\left(\widehat{A}\widehat{U}(\overrightarrow{x})\right).$$

(41)

The operator $\widehat{A}$ can be reconstructed from its symbol $f_A(\overrightarrow{x})$ using the inverse transform $f_A(\overrightarrow{x})\to \widehat{A}$, i.e.,

$$\widehat{A}=\int f_A(\overrightarrow{x})\widehat{D}(\overrightarrow{x})d\overrightarrow{x}.$$

(42)

So, the operator $\widehat{A}$ can be reconstructed from its symbol if the quantizer operator $\widehat{D}(\overrightarrow{x})$ and the symbol of operator $f_A(\overrightarrow{x})$ are known. The map given by Equations (41) and (42) provides the possibility to introduce the star product of functions $f_A(\overrightarrow{x})$ and $f_B(\overrightarrow{x})$, which are symbols of operators $\widehat{A}$ and $\widehat{B}$. The symbol of operator $\widehat{A}\widehat{B}$, which is product of operators $\widehat{A}$ and $\widehat{B}$, is

$$f_{AB}(\overrightarrow{x})=\mathrm{Tr}\left(\widehat{A}\widehat{B}\widehat{U}(\overrightarrow{x})\right).$$

(43)

We can present the star product of the functions $f_A(\overrightarrow{x})$ and $f_B(\overrightarrow{x})$

$$\left(f_A\star f_B\right)(\overrightarrow{x})=f_{AB}(\overrightarrow{x})$$

(44)

using the relationships (41)–(43) in the integral form

$$\left(f_A\star f_B\right)(\overrightarrow{x})=\int f_A({\overrightarrow{x}}_1)f_B({\overrightarrow{x}}_2)K({\overrightarrow{x}}_1,{\overrightarrow{x}}_2,\overrightarrow{x})d{\overrightarrow{x}}_{1}d{\overrightarrow{x}}_2.$$

(45)

The kernel $K({\overrightarrow{x}}_1,{\overrightarrow{x}}_2,\overrightarrow{x})$ is expressed in terms of the quantizer and dequantizer operators

$$K({\overrightarrow{x}}_1,{\overrightarrow{x}}_2,\overrightarrow{x})=\mathrm{Tr}\left(\widehat{D}({\overrightarrow{x}}_1)\widehat{D}({\overrightarrow{x}}_2)\widehat{U}(\overrightarrow{x})\right).$$

(46)

The associativity condition for the product of operators, i.e., $((\widehat{A}\widehat{B})\widehat{C})=(\widehat{A}(\widehat{B}\widehat{C}))$ causes the star product of operator symbols to be associative as well. The formalism of quantizer–dequantizer operators can be used to determine the evolution equation for the symbols of the density operators. The von Neumann equation for the oscillator density operator $\widehat{\rho }(t)$ is written as (we use $m=\omega =\hslash =1$)

$$\frac{\partial \widehat{\rho }}{\partial t}+i\left(\widehat{H}(t)\widehat{\rho }(t)-\widehat{\rho }(t)\widehat{H}(t)\right)=0.$$

(47)

If we introduce the symbol $f_{\rho }(\overrightarrow{x},t)$ of the density operator $\widehat{\rho }(t)$ and the symbol $f_H(\overrightarrow{x},t)$ of the Hamiltonian operator $\widehat{H}(t)$, using an arbitrary pair of quantizer and dequantizer operators, then Equation (47) becomes

$$\frac{\partial f_{\rho }(\overrightarrow{x},t)}{\partial t}+i\left(f_H\star f_{\rho }-f_{\rho }\star f_H\right)(\overrightarrow{x},t)=0.$$

(48)

The equation for the evolution of the density operator symbol for the given Hamiltonian $\widehat{H}(t)$ has the general form of an integral equation

$$\frac{\partial f_{\rho }(\overrightarrow{x},t)}{\partial t}+i\int \left(f_H({\overrightarrow{x}}_1,t)f_{\rho }({\overrightarrow{x}}_2,t)-f_{\rho }({\overrightarrow{x}}_1,t)f_H({\overrightarrow{x}}_2,t)\right)K({\overrightarrow{x}}_1,{\overrightarrow{x}}_2,\overrightarrow{x})d{\overrightarrow{x}}_{1}d{\overrightarrow{x}}_2=0.$$

(49)

Here, the symbol of the Hamiltonian $f_H({\overrightarrow{x}}_1,t)=\mathrm{Tr}(\widehat{H}(t)\widehat{U}({\overrightarrow{x}}_1))$ and the symbol of density operator $f_{\rho }({\overrightarrow{x}}_2,t)=\mathrm{Tr}(\widehat{\rho }(t)\widehat{U}({\overrightarrow{x}}_2)).$ Using (39), (46), and (49), one obtains

$$\begin{array}{ccc}& & \frac{\partial f_{\rho }(\overrightarrow{x},t)}{\partial t}+i\int \left[\mathrm{Tr}(\widehat{H}(t)\widehat{U}({\overrightarrow{x}}_1))\mathrm{Tr}(\widehat{\rho }(t)\widehat{U}({\overrightarrow{x}}_2))-\mathrm{Tr}(\widehat{\rho }(t)\widehat{U}({\overrightarrow{x}}_1))\mathrm{Tr}(\widehat{H}(t)\widehat{U}({\overrightarrow{x}}_2))\right]\hfill \\ & & \times \mathrm{Tr}(\widehat{D}({\overrightarrow{x}}_1)\widehat{D}({\overrightarrow{x}}_2)\widehat{U}(\overrightarrow{x}))d{\overrightarrow{x}}_{1}d{\overrightarrow{x}}_2=0.\hfill \end{array}$$

(50)

Equation (50) can be written in the form of a kinetic equation for a probability distribution function

$$\frac{\partial f_{\rho }(\overrightarrow{x},t)}{\partial t}+i\int f_{\rho }({\overrightarrow{x}}_2,t)\mathcal{K}(\overrightarrow{x},{\overrightarrow{x}}_2,t)d{\overrightarrow{x}}_2=0,$$

(51)

where

$$\mathcal{K}(\overrightarrow{x},{\overrightarrow{x}}_2,t)=\int \left(K({\overrightarrow{x}}_1,{\overrightarrow{x}}_2,x,t)-K({\overrightarrow{x}}_2,{\overrightarrow{x}}_1,x,t)\right)f_H({\overrightarrow{x}}_1,t)d{\overrightarrow{x}}_1,$$

(52)

and the symbol of density operator $f_{\rho }(\overrightarrow{x},t)$ is a probability distribution. For the symplectic tomogram, the inverse quantum Radon transform reads [53,54]

$$\widehat{\rho }=\frac{1}{2\pi }\int w(X|\mu ,\nu )exp\left(i(X\widehat{1}-\mu \widehat{q}-\nu \widehat{p})\right)dXd\mu d\nu .$$

(53)

This means that the quantizer operator for the symplectic tomography method has the form

$$\widehat{D}(X|\mu ,\nu )=\frac{1}{2\pi }exp\left(i(X\widehat{1}-\mu \widehat{q}-\nu \widehat{p})\right).$$

(54)

Thus, we have $\overrightarrow{x}=X,\mu ,\nu$, and the dequantizer reads

$$\widehat{U}(X|\mu ,\nu )=\delta \left(i(X\widehat{1}-\mu \widehat{q}-\nu \widehat{p})\right).$$

(55)

The existence of (55) provides the possibility to map the density operator $\widehat{\rho }$ (applying Formula (39)) onto the function (tomogram of the state). The existence of quantizer (54) provides the possibility to reconstruct the density operator $\widehat{\rho }$ using the tomogram of the quantum state. Such pairs of quantizer–dequantizer operators exist for all the other quantum systems, including the two-mode usual and inverted oscillators discussed in this paper.

The kernel describing the star product of the operators in the symplectic tomography is expressed as follows

$$\begin{array}{ccc}& & K(X_1,{\mu }_1,{\nu }_1,X_2,{\mu }_2,{\nu }_2,X,\mu ,\nu )=\frac{1}{4{\pi }^2}\mathrm{Tr}\left[exp\left(i(X_1\widehat{1}-{\mu }_1\widehat{q}-{\nu }_1\widehat{p})\right)\right.\hfill \\ & & \times \left.exp\left(i(X_2\widehat{1}-{\mu }_2\widehat{q}-{\nu }_2\widehat{p})\right)\delta \left(i(X\widehat{1}-\mu \widehat{q}-\nu \widehat{p})\right)\right].\hfill \end{array}$$

(56)

In an explicit form it reads

$$\begin{array}{ccc}& & K(X_1,{\mu }_1,{\nu }_1,X_2,{\mu }_2,{\nu }_2,X,\mu ,\nu )=\frac{1}{4{\pi }^2}\delta \left(\mu ({\nu }_1+{\nu }_2)-\nu ({\mu }_1+{\mu }_2)\right)\hfill \\ & & \times exp\left(\frac{i}{2}({\nu }_1{\mu }_2-{\nu }_2{\mu }_1+2X_1+2X_2-2\frac{{\nu }_1+{\nu }_2}{\nu }X)\right).\hfill \end{array}$$

(57)

In the case of a harmonic oscillator in the tomographic probability representation, the symbol of density operator $\widehat{\rho }(t)$ is given by the probability distribution function, $(\overrightarrow{x}=X,\mu ,\nu ),$

$$w_{\rho }(X|\mu ,\nu ,t)=f_{\rho }(\overrightarrow{x},t)=\mathrm{Tr}\widehat{\rho }(t)\delta \left(X\widehat{1}-\mu \widehat{q}-\nu \widehat{p}\right),$$

(58)

The Hamiltonian $\widehat{H}$ can be mapped onto its symbol

$$f_{\widehat{H}}(X,\mu ,\nu )=\mathrm{Tr}\left(\widehat{H}\delta (X\widehat{1}-\mu \widehat{q}-\nu \widehat{p})\right).$$

(59)

The symplectic tomogram (58) is the symbol of the density operator $\widehat{\rho }$, and it is the probability distribution of position X [13] depending on extra parameters determining the reference frame in the phase space where the position X is measured. For symplectic tomography, the integral linear Equation (51) has the form

$$\frac{\partial w_{\rho }(X|\mu ,\nu ,t)}{\partial t}+i\int w_{\rho }(X_2|{\mu }_2,{\nu }_2,t)\mathcal{K}(X,\mu ,\nu ,X_2,{\mu }_2,{\nu }_2,t)d{X}_{2}d{\mu }_{2}d{\nu }_2=0.$$

(60)

Here,

$$\begin{array}{ccc}& & \mathcal{K}(X,\mu ,\nu ,X_2,{\mu }_2,{\nu }_2,t)=\hfill \\ & & \int \left[K(X_1,{\mu }_1,{\nu }_1,X_2,{\mu }_2,{\nu }_2,t)-K(X_2,{\mu }_2,{\nu }_2,X_1,{\mu }_1,{\nu }_1,t)\right]f_H(X_1,{\mu }_1,{\nu }_1,t)d{X}_{1}d{\mu }_{1}d{\nu }_1.\hfill \end{array}$$

(61)

The product of operators $\widehat{A}·\widehat{B}$ is mapped onto the star product of their symbols

$$(\widehat{A}\widehat{B})\leftrightarrow (A\star B)(X,\mu ,\nu )=\mathrm{Tr}\left(\widehat{A}\widehat{B}\delta (X\widehat{1}-\mu \widehat{q}-\nu \widehat{p})\right)$$

(62)

with the kernel of the star product defined by means of the expression

$$\begin{array}{ccc}& & (A\star B)(X,\mu ,\nu )=\hfill \\ & & \int A(X_1,{\mu }_1{\nu }_1)B(X_2,{\mu }_2,{\nu }_2)K(X_1,{\mu }_1{\nu }_1,X_2,{\mu }_2,{\nu }_2,X,\mu ,\nu )d{X}_{1}d{X}_{2}d{\mu }_{1}d{\mu }_{2}d{\nu }_{1}d{\nu }_2.\hfill \end{array}$$

(63)

This formula is the application of general Formula (46) for the kernel of the star product of the symbols. This formula can be used to study the entanglement phenomena of states that are superpositions of two-mode oscillator states.

6. Symplectic Tomography of Oscillators

One can calculate the tomographic probability distribution $w(X|\mu ,\nu )$, called the symplectic tomogram of the state with density operator ${\widehat{\rho }}_{|\psi \rangle }=|\psi \rangle \langle \psi |$, using the formula analogous to (15) expressed in terms of wave function $\psi (y)$ of the pure state in position representation, which reads [50]

$$w_{|\psi \rangle }(X|\mu ,\nu )=\frac{1}{2\pi \left|\nu \right|}{\left|\int \psi (y)exp\left(\frac{i\mu }{2\nu }{y}^2-\frac{iXy}{\nu }\right)dy\right|}^2.$$

(64)

The function is non-negative and satisfies the normalization condition

$$\int w_{|\psi \rangle }(X|\mu ,\nu )dX=1.$$

(65)

The physical meaning of the real parameters $\mu$ and $\nu$ is that they, due to using $\delta (X\widehat{1}-\mu \widehat{q}-\nu \widehat{p})$ to determine the dequantizer $\widehat{U}(x)$ as a delta function, describe the axes of reference frames in the phase space of position $\widehat{q}$ and momentum $\widehat{p}$, where the position $X\widehat{1}=\mu \widehat{q}+\nu \widehat{p}$ is measured. Thus, the tomogram $w(X|\mu ,\nu )$ is the conditional probability distribution determining the density operator for the state. If $\mu =1,\nu =0$, it is the density matrix diagonal elements $\rho (qq)$, and for $\mu =0,\nu =1$, the tomogram is the diagonal matrix element $\rho (pp)$. This means that if one knows the probability distributions of position and momentum in all the reference frames in the phase space, the state (state density operator) is known.

In the case of two-mode oscillators, the relationship between the symplectic tomogram and wave function of the state is obtained by Equation (15). Using (15), one can obtain the symplectic tomogram of the ground state of two-mode oscillators in the explicit form

$$w_0(X_1,X_2|{\mu }_1,{\mu }_2,{\nu }_1,{\nu }_2)=\frac{1}{\pi \sqrt{{\mu }_1^2+{\nu }_1^2}\sqrt{{\mu }_2^2+{\nu }_2^2}}exp\left(-\frac{X_1^2}{{\mu }_1^2+{\nu }_1^2}-\frac{X_2^2}{{\mu }_2^2+{\nu }_2^2}\right),$$

(66)

and the tomogram of coherent state of two-mode oscillators in the Gaussian form

$$w_{\alpha }(X_1,X_2|{\mu }_1,{\mu }_2,{\nu }_1,{\nu }_2)=\frac{1}{\pi \sqrt{{\mu }_1^2+{\nu }_1^2}\sqrt{{\mu }_2^2+{\nu }_2^2}}exp\left(-\frac{{(X_1-{\overline{X}}_1)}^2}{{\mu }_1^2+{\nu }_1^2}-\frac{{(X_2-{\overline{X}}_2)}^2}{{\mu }_2^2+{\nu }_2^2}\right),$$

(67)

where ${\overline{X}}_1=\sqrt{2}{\mu }_1\mathrm{Re}\alpha +\sqrt{2}{\nu }_1\mathrm{Im}\alpha$, ${\overline{X}}_2=\sqrt{2}{\mu }_2\mathrm{Re}\alpha +\sqrt{2}{\nu }_2\mathrm{Im}\alpha$, and $\alpha$ is a complex number.

7. Dynamics of Operator Symbols for Quadratic Hamiltonians in Position and Momentum

Let us discuss the problem of finding the tomographic probability distribution evolution for the systems with Hamiltonians, which are quadratic forms in position and momentum operators. Such systems have integrals of motion that are linear in position and momentum operators. Furthermore, the position and momentum operators ${\widehat{q}}_H(t)$ and ${\widehat{p}}_H(t)$ are linear forms of the position $\widehat{q}$ and momentum $\widehat{p}$ operators with time-dependent coefficients [55]. Due to this, we can explicitly obtain the time dependence of the tomographic probability distributions describing the quantum states and corresponding to solutions of the Schrödinger equation for wave functions and the von Neumann equation for the density operators. The idea to obtain the solutions of these equations was formulated in [48,52]. It is based on the following observation: Since the system state tomogram is given by the symbol of the density operator $\widehat{\rho }(t)$, i.e., (58), where the density operator evolution for the von Neumann equation is described by the evolution operator $\widehat{u}(t)$, i.e., $\widehat{\rho }(t)=\widehat{u}(t)\widehat{\rho }(0){\widehat{u}}^{†}(t)$ the symbol of the density operator can be rewritten in the form $\mathrm{Tr}(\widehat{\rho }(0)\delta (X-\mu {\widehat{u}}^{†}(t)\widehat{q}\widehat{u}(t)-\nu {\widehat{u}}^{†}(t)\widehat{p}\widehat{u}(t)).$ Here, ${\widehat{u}}^{†}(t)\widehat{q}\widehat{u}(t)={\widehat{q}}_H(t)$ and ${\widehat{u}}^{†}(t)\widehat{p}\widehat{u}(t)={\widehat{p}}_H(t)$ are the Heisenberg position and momentum operators. Such properties also occur in multi-mode systems with Hamiltonians, which are any quadratic forms in position and momentum operators; for example, for two-dimensional oscillators, both ordinary,

$${\widehat{H}}^{(1)}=\frac{{\widehat{p}}_1^2}{2}+\frac{{\widehat{p}}_2^2}{2}+\frac{{\widehat{q}}_2^2}{2}+\frac{{\widehat{q}}_1^2}{2},$$

(68)

and for two-dimensional oscillators, both inverted,

$${\widehat{H}}^{(2)}=\frac{{\widehat{p}}_1^2}{2}+\frac{{\widehat{p}}_2^2}{2}-\frac{{\widehat{q}}_2^2}{2}-\frac{{\widehat{q}}_1^2}{2}.$$

(69)

The Hamiltonian ${\widehat{H}}_2$ corresponds to the motion of the inverted oscillator. For such Hamiltonians, one has time-dependent Heisenberg operators of position and momentum of the following forms: for the ordinary oscillator with the Hamiltonian (68)

$${\widehat{q}}_{H^{(1)};1}(t)=cost·{\widehat{q}}_1+sint·{\widehat{p}}_1,{\widehat{q}}_{H^{(1)};2}(t)=cost·{\widehat{q}}_2+sint·{\widehat{p}}_2,$$

(70)

$${\widehat{p}}_{H^{(1)};1}(t)=-sint·{\widehat{q}}_1+cost·{\widehat{p}}_1,{\widehat{p}}_{H^{(1)};2}(t)=-sint·{\widehat{q}}_2+cost·{\widehat{p}}_2;$$

(71)

and for the inverted oscillator with the Hamiltonian (69)

$${\widehat{q}}_{H^{(2)};1}(t)=cosht·{\widehat{q}}_1+sinht·{\widehat{p}}_1,{\widehat{q}}_{H^{(2)};2}(t)=cosht·{\widehat{q}}_2+sinht·{\widehat{p}}_2,$$

(72)

$${\widehat{p}}_{H^{(2)};1}(t)=sinht·{\widehat{q}}_1+cosht·{\widehat{p}}_1,{\widehat{p}}_{H^{(2)};2}(t)=sinht·{\widehat{q}}_2+cosht·{\widehat{p}}_2.$$

(73)

Developed formalism provides the possibility to obtain the description of time evolution for all the multi-mode systems with time-dependent quadratic Hamiltonians. For such systems, the Heisenberg position and momentum operators are linear forms with time-dependent coefficients of usual positions and momenta.

8. Center-of-Mass Tomography

Let us introduce a dequantizer operator for the two-mode oscillator $\widehat{U}(X_1,X_2,{\mu }_1,{\nu }_1,{\mu }_1,{\nu }_2)$. Then, the symplectic tomogram reads

$$w(X_1,X_2|{\mu }_1,{\nu }_1,{\mu }_2,{\nu }_2)=\mathrm{Tr}\left(\widehat{\rho }\delta (X_1\widehat{1}-{\mu }_1{\widehat{q}}_1-{\nu }_1{\widehat{p}}_1)\delta (X_2\widehat{1}-{\mu }_2{\widehat{q}}_2-{\nu }_2{\widehat{p}}_2)\right).$$

(74)

The dequantizer operator $\widehat{U}(\overrightarrow{x})$ in the case of the symplectic tomogram is $\delta (X_1\widehat{1}-{\mu }_1{\widehat{q}}_1-{\nu }_1{\widehat{p}}_1)\delta (X_2\widehat{1}-{\mu }_2{\widehat{q}}_2-{\nu }_2{\widehat{p}}_2)$. The density operator can be reconstructed from the symplectic tomogram with the help of the quantizer operator $\widehat{D}(X_1,X_2,{\mu }_1,{\mu }_2,{\nu }_1,{\nu }_2)=\frac{1}{4{\pi }^2}exp(i{X}_1\widehat{1}-{\mu }_1{\widehat{q}}_1-{\nu }_1{\widehat{p}}_1)exp(i{X}_2\widehat{1}-{\mu }_2{\widehat{q}}_2-{\nu }_2{\widehat{p}}_2)$. Then, the symplectic tomogram of the first mode of oscillator is related to (74) as

$$w(X_1|{\mu }_1,{\nu }_1)=\int w(X_1,X_2|{\mu }_1,{\nu }_1,{\mu }_2,{\nu }_2)d{X}_2.$$

(75)

There is another type of tomography, named center-of-mass tomography. It was introduced in [56] and developed in [57,58]. In center-of-mass tomography, the state is determined by the center-of-mass tomogram. The center-of-mass tomogram is a symbol of the density operator

$$w_{cm}(X|{\mu }_1,{\nu }_1,{\mu }_2,{\nu }_2)=\mathrm{Tr}\left(\widehat{\rho }\delta (X\widehat{1}-{\mu }_1{\widehat{q}}_1-{\nu }_1{\widehat{p}}_1-{\mu }_2{\widehat{q}}_2-{\nu }_2{\widehat{p}}_2)\right).$$

(76)

The random variable X, which we named the center-of-mass coordinate, is measured in phase space in rotated and scaled reference frames, which are determined by parameters ${\mu }_1,{\nu }_1,{\mu }_2,{\nu }_2$. The dequantizer operator in center-of-mass tomography is

$$\widehat{U}(X,{\mu }_1,{\nu }_1,{\mu }_2,{\nu }_2)=\delta (X\widehat{1}-{\mu }_1{\widehat{q}}_1-{\nu }_1{\widehat{p}}_1-{\mu }_2{\widehat{q}}_2-{\nu }_2{\widehat{p}}_2).$$

(77)

The density operator can be reconstructed from the center-of-mass tomogram with the help of the quantizer operator $\widehat{D}(X,{\mu }_1,{\nu }_1,{\mu }_1,{\nu }_2)$, i.e.,

$$\widehat{\rho }=\frac{1}{4{\pi }^2}\int w_{cm}(X|{\mu }_1,{\nu }_1,{\mu }_2.{\nu }_2)exp\left(i\left(X\widehat{1}-{\mu }_1\widehat{q}-{\nu }_1{\widehat{p}}_1-{\mu }_2{\widehat{q}}_2-{\nu }_2{\widehat{p}}_2\right)\right)dXd{\mu }_{1}d{\nu }_{1}d{\mu }_{2}d{\nu }_2.$$

(78)

The center-of-mass tomogram of odd and even coherent states is of the form

$$\begin{array}{ccc}& & w_{cm,\alpha }(X|{\mu }_1,{\nu }_1,{\mu }_2,{\nu }_2)=\hfill \\ & & \frac{1}{\sqrt{\pi }\sqrt{\sigma }{N}_{\pm }^2(\alpha )}\left[exp\left(-{(X-\sqrt{2}\mathrm{Re}{\alpha }_1{\mu }_1-\sqrt{2}\mathrm{Re}{\alpha }_2{\mu }_2-\sqrt{2}\mathrm{Im}{\alpha }_1{\nu }_1-\sqrt{2}\mathrm{Im}{\alpha }_2{\nu }_2)}^2/\sigma \right)\right.\hfill \\ & & \pm exp\left(-2|{\alpha }_1|-2|{\alpha }_2|-{(X-i\sqrt{2}\mathrm{Im}{\alpha }_1{\mu }_1-i\sqrt{2}\mathrm{Im}{\alpha }_2{\mu }_2+i\sqrt{2}\mathrm{Re}{\alpha }_1{\nu }_1+i\sqrt{2}\mathrm{Re}{\alpha }_2{\nu }_2)}^2/\sigma \right)\hfill \\ & & \pm exp\left(-2|{\alpha }_1|-2|{\alpha }_2|-{(X+i\sqrt{2}\mathrm{Im}{\alpha }_1{\mu }_1+i\sqrt{2}\mathrm{Im}{\alpha }_2{\mu }_2-i\sqrt{2}\mathrm{Re}{\alpha }_1{\nu }_1-i\sqrt{2}\mathrm{Re}{\alpha }_2{\nu }_2)}^2/\sigma \right)\hfill \\ & & +\left.exp\left(-{(X+\sqrt{2}\mathrm{Re}{\alpha }_1{\mu }_1+\sqrt{2}\mathrm{Re}{\alpha }_2{\mu }_2+\sqrt{2}\mathrm{Im}{\alpha }_1{\nu }_1+\sqrt{2}\mathrm{Im}{\alpha }_2{\nu }_2)}^2/\sigma \right)\right],\hfill \end{array}$$

(79)

where $\sigma ={\mu }_1^2+{\mu }_2^2+{\nu }_1^2+{\nu }_2^2$ and $N_{\pm }^2(\alpha )=2\left(1\pm exp(-2|{\alpha }_1{|}^2-2|{\alpha }_2{|}^2)\right)$. These tomograms (79) determine the nonclassical even and odd coherent states in the probability representation of quantum mechanics.

Following the method described in Section 7, we obtain the time-dependent center-of-mass tomogram of Schrödinger cat states. This means that in Formula (79), we have to replace ${\mu }_1,{\nu }_1,{\mu }_2,{\nu }_2$ with time-dependent Heisenberg parameters in the case of evolution with the Hamiltonian of ordinary oscillator (68) of the form

$$\begin{array}{ccc}& & {\mu }_{H^{(1)};1}={\mu }_{1}cost-{\nu }_{1}sint,{\mu }_{H^{(1)};2}={\mu }_{2}cost-{\nu }_{2}sint,\hfill \\ & & {\nu }_{H^{(1)};1}={\mu }_{1}sint+{\nu }_{1}cost,{\nu }_{H^{(1)};2}={\mu }_{2}sint+{\nu }_{2}cost.\hfill \end{array}$$

(80)

So, for the initial center-of-mass tomogram of odd and even states given by (79) after evolution with the Hamiltonian (68), one obtains the explicit expression

$$\begin{array}{ccc}& & w_{cm,\alpha }(X|{\mu }_1,{\nu }_1,{\mu }_2,{\nu }_2,t)=\hfill \\ & & \frac{1}{\sqrt{\pi }\sqrt{\sigma }{N}_{\pm }^2(\alpha )}\left[exp\left(-(X-\sqrt{2}\mathrm{Re}{\alpha }_1({\mu }_{1}cost-{\nu }_{1}sint)-\sqrt{2}\mathrm{Re}{\alpha }_2({\mu }_{2}cost-{\nu }_{2}sint)\right.\right.\hfill \\ & & \left.-\sqrt{2}\mathrm{Im}{\alpha }_1({\mu }_{1}sint+{\nu }_{1}cost)-\sqrt{2}\mathrm{Im}{\alpha }_2({\mu }_{2}sint+{\nu }_{2}cost){)}^2/\sigma \right)\hfill \\ & & \pm exp\left(-2|{\alpha }_1|-2|{\alpha }_2|-(X-i\sqrt{2}\mathrm{Im}{\alpha }_1({\mu }_{1}cost-{\nu }_{1}sint)\right.\hfill \\ & & -i\sqrt{2}\mathrm{Im}{\alpha }_2({\mu }_{2}cost-{\nu }_{2}sint)+i\sqrt{2}\mathrm{Re}{\alpha }_1({\mu }_{1}sint+{\nu }_{1}cost)\hfill \\ & & \left.+i\sqrt{2}\mathrm{Re}{\alpha }_2({\mu }_{2}sint+{\nu }_{2}cost){)}^2/\sigma \right)\hfill \\ & & \pm exp\left(-2|{\alpha }_1|-2|{\alpha }_2|-(X+i\sqrt{2}\mathrm{Im}{\alpha }_1({\mu }_{1}cost-{\nu }_{1}sint)\right.\hfill \\ & & +i\sqrt{2}\mathrm{Im}{\alpha }_2({\mu }_{2}cost-{\nu }_{2}sint)-i\sqrt{2}\mathrm{Re}{\alpha }_1({\mu }_{1}sint+{\nu }_{1}cost)\hfill \\ & & \left.-i\sqrt{2}\mathrm{Re}{\alpha }_2({\mu }_{2}sint+{\nu }_{2}cost){)}^2/\sigma \right)\hfill \\ & & +exp\left(-(X+\sqrt{2}\mathrm{Re}{\alpha }_1({\mu }_{1}cost-{\nu }_{1}sint)+\sqrt{2}\mathrm{Re}{\alpha }_2({\mu }_{2}cost-{\nu }_{2}sint)\right.\hfill \\ & & \left.\left.+\sqrt{2}\mathrm{Im}{\alpha }_1({\mu }_{1}sint+{\nu }_{1}cost)+\sqrt{2}\mathrm{Im}{\alpha }_2({\mu }_{2}sint+{\nu }_{2}cost){)}^2/\sigma \right)\right],\hfill \end{array}$$

(81)

where $\sigma ={\mu }_1^2+{\mu }_2^2+{\nu }_1^2+{\nu }_2^2$.

For inverted oscillators with Hamiltonian (69), the initial center-of-mass tomogram given by (79) takes the form of conditional probability distribution of one random variable X

$$\begin{array}{ccc}& & w_{cm,\alpha }(X|{\mu }_1,{\nu }_1,{\mu }_2,{\nu }_2,t)=\hfill \\ & & \frac{1}{\sqrt{\pi }\sqrt{\sigma }{N}_{\pm }^2(\alpha )}\left[exp\left(-(X-\sqrt{2}\mathrm{Re}{\alpha }_1({\mu }_{1}cosht+{\nu }_{1}sinht)-\sqrt{2}\mathrm{Re}{\alpha }_2({\mu }_{2}cosht+{\nu }_{2}sinht)\right.\right.\hfill \\ & & \left.-\sqrt{2}\mathrm{Im}{\alpha }_1({\mu }_{1}sinht+{\nu }_{1}cosht)-\sqrt{2}\mathrm{Im}{\alpha }_2({\mu }_{2}sinht+{\nu }_{2}cosht){)}^2/\sigma \right)\hfill \\ & & \pm exp\left(-2|{\alpha }_1|-2|{\alpha }_2|-(X-i\sqrt{2}\mathrm{Im}{\alpha }_1({\mu }_{1}cosht+{\nu }_{1}sinht)\right.\hfill \\ & & -i\sqrt{2}\mathrm{Im}{\alpha }_2({\mu }_{2}cosht+{\nu }_{2}sinht)+i\sqrt{2}\mathrm{Re}{\alpha }_1({\mu }_{1}sinht+{\nu }_{1}cosht)\hfill \\ & & \left.+i\sqrt{2}\mathrm{Re}{\alpha }_2({\mu }_{2}sinht+{\nu }_{2}cosht){)}^2/\sigma \right)\hfill \\ & & \pm exp\left(-2|{\alpha }_1|-2|{\alpha }_2|-(X+i\sqrt{2}\mathrm{Im}{\alpha }_1({\mu }_{1}cosht+{\nu }_{1}sinht)\right.\hfill \\ & & +i\sqrt{2}\mathrm{Im}{\alpha }_2({\mu }_{2}cosht+{\nu }_{2}sinht)-i\sqrt{2}\mathrm{Re}{\alpha }_1({\mu }_{1}sinht+{\nu }_{1}cosht)\hfill \\ & & \left.-i\sqrt{2}\mathrm{Re}{\alpha }_2({\mu }_{2}sinht+{\nu }_{2}cosht){)}^2/\sigma \right)\hfill \\ & & +exp\left(-(X+\sqrt{2}\mathrm{Re}{\alpha }_1({\mu }_{1}cosht+{\nu }_{1}sinht)+\sqrt{2}\mathrm{Re}{\alpha }_2({\mu }_{2}cosht+{\nu }_{2}sinht)\right.\hfill \\ & & \left.\left.+\sqrt{2}\mathrm{Im}{\alpha }_1({\mu }_{1}sinht+{\nu }_{1}cosht)+\sqrt{2}\mathrm{Im}{\alpha }_2({\mu }_{2}sinht+{\nu }_{2}cosht){)}^2/\sigma \right)\right],\hfill \end{array}$$

(82)

where $\sigma =cosh2t({\mu }_1^2+{\mu }_2^2+{\nu }_1^2+{\nu }_2^2)+2sinh2t({\mu }_1{\nu }_1+{\mu }_2{\nu }_2)$.

9. Conclusions

For almost a hundred years since the discovery and development of quantum mechanics, there has been a problem with understanding its foundations and mysteries. One of the main mysteries is the following problem. Why a classical object, such as the Moon orbiting the Earth, exhibits behavior that we understand. Newton’s law of classical elliptical motion ensures that the Moon rotates around the Globe with constant energy. This is intuitive and fits our daily experience. However, since the beginning of the quantum age, we find that the electron in the hydrogen atom revolves around the proton and has a spectrum of energy levels that is completely contrary to our classical intuition. To understand this quantum behavior, a wave function dynamics formalism was developed that obeys the Schrödinger equation and completely contradicts our classical intuition. The novelty of this article lies in the fact that the results obtained within it open up the possibility of extending them to all other quantum systems, and thus enhancing our understanding of quantum phenomena, closer to the classical intuition. As an alternative to wave functions and density operators describing the states of quantum objects, in the last century, a probability distribution was proposed to describe the states of electrons.

Of course, quantum objects, such as electrons, move around protons according to the equations of quantum evolution, such as the Schrödinger equation. They have energy corresponding to quantized energy levels, and their trajectories do not coincide with the trajectories of the Moon’s rotation around the Globe. In this sense, there is no direct similarity between the quantum dynamics of electrons and the classical dynamics of the Moon. However, the similarity between understanding quantum motion and classical motion lies in the fact that both the motion of the Moon and the motion of electrons can be described by conventional probability distribution functions.

This description is closer to our intuition and experience of the statistical description of random motions. In science, this description has been the subject of lengthy research, and over the past two decades it has been achieved. The time evolution of oscillators, both conventional and inverted, is also studied in the present paper. This consideration opens up the possibility of extending the probabilistic approach describing quantum states to all other systems in order to apply these results to the development of new quantum technologies, such as quantum computing, quantum information technology, etc.

To conclude, we summarize the main results of our paper. We developed the probability representation of quantum states in which the system states are described by standard probability distribution functions. These functions determine the density operators of the states. For this, we considered two different schemes of such construction, namely, symplectic tomography probability distributions [13] and center-of-mass tomographic probability distributions [56]. In our work, we considered time evolution of the tomographic probability distributions, using an example of the Schrödinger cat states of the two-mode oscillator. The main aim was to determine the time evolution and the explicit expressions of center-of-mass tomographic probability distributions for even and odd coherent states of two-mode oscillators in ordinary (81) and inverted (82) forms. The main result is that the obtained probability distributions describe the entangled states of two-mode oscillator and its evolution. In an example of such a state, we constructed the entangled probability distributions and their dynamics. Furthermore, we studied the structure of generic entangled probability distributions. The entangled probability distributions are the new forms of standard probability distributions [12]. The possibility to construct such new probability distributions can be studied considering multi-mode oscillators with time-dependent parameters. The entangled probability distributions are new forms of distributions introduced using quantum mechanics. There are other new aspects of the classical probability theory that can be found and formulated in view of the existence of the quantum formalism of the Hilbert spaces and operators acting in the Hilbert spaces, such as Bell inequalities, which can be considered as consequences of the entangled probability distributions as well as several entropic inequalities that are obvious in quantum mechanics; however, these relations are poorly clarified, and were not even discussed in the classical probability theory. We will consider these problems and entropic properties of such probability distributions in future publications.