The Group-Algebraic Formalism of Quantum Probability and Its Applications in Quantum Statistical Mechanics

Abstract

We show that the theory of quantum statistical mechanics is a special model in the framework of the quantum probability theory developed by mathematicians, by extending the characteristic function in the classical probability theory to the quantum probability theory. As dynamical variables of a quantum system must respect certain commutation relations, we take the group generated by a Lie algebra constructed with these commutation relations as the bridge, so that the classical characteristic function defined in a Euclidean space is transformed to a normalized, non-negative definite function defined in this group. Indeed, on the quantum side, this group-theoretical characteristic function is equivalent to the density matrix; hence, it can be adopted to represent the state of a quantum ensemble. It is also found that this new representation may have significant advantages in applications. As two examples, we show its effectiveness and convenience in solving the quantum-optical master equation for a harmonic oscillator coupled with its thermal environment, and in simulating the quantum cat map, a paradigmatic model for quantum chaos. Other related issues are reviewed and discussed as well.

1. Introduction

Probability signifies our ignorance of the exact cause of an event. The theory of probability, as a mathematical branch, emerged in the seventeenth century. According to Laplace [1], probability can be determined as the ratio of favorable outcomes to all equally possible alternatives. But by 1900, this classical interpretation of probability had become outdated and a new statistical interpretation gradually became prevalent. In this new interpretation, it does not make sense to talk of the probability of a single event. Indeed, as exemplified by statistical mechanics, probability is defined in terms of the distribution of an ensemble in the phase space [2], or, in a statistical sense, as frequencies with reference to infinite sequences of repeatable events [3].

The modern mathematical axiomatic framework of probability theory was proposed by Kolmogorov in 1933 [4]. In Kolmogorov’s probability theory, the probability space is defined by a triple $(\Omega ,\Sigma ,\mu )$, where $\Omega$ is an arbitrary set (called a sample space), $\Sigma$ is a Boolean $\sigma$-algebra of subsets of $\Omega$ (called measurable sets), and $\mu$ is a probability measure on $\Sigma$ [5]. A random variable associated with this probability space is defined as a measurable function $\{\omega \in \Omega \to f(\omega )\in \mathbb{R}\}$ in the sample space $\Omega$ with the probability distribution function $F\left(x\right)={\int }_{f\left(\omega \right)\le x}d\mu \left(\omega \right)$.

When quantum mechanics took shape around 1925–1927, the notion of probability with its statistical interpretation was introduced into quantum mechanics by Born [6], by which the implications of a wave function of a quantum system are successfully decoded. It is well known that the orthodox mathematical formalism of quantum theory proposed by von Neumann [7] is quite different from the Kolmogorovian model. In von Neumann’s theory, states and observables are represented with density matrices and self-adjoint operators in a Hilbert space, respectively. Consequently, the probability space of a quantum system can be defined by a triple $(\mathcal{H},\mathcal{P}\left(\mathcal{H}\right),\widehat{\rho })$, where $\mathcal{H}$ is a separable Hilbert space, $\mathcal{P}\left(\mathcal{H}\right)$ is the set of all projection operators in the closed subspaces of $\mathcal{H}$, and $\widehat{\rho }$ is a density matrix in $\mathcal{H}$. A random variable associated with this probability space is defined as a self-adjoint operator $\widehat{X}$ in the Hilbert space $\mathcal{H}$ with its spectral representation $\widehat{X}=\int \lambda d\widehat{E}\left(\lambda \right)$, $\widehat{E}\left(\lambda \right)\in \mathcal{P}\left(\mathcal{H}\right)$, and the expectation value $⟨\widehat{X}⟩=\mathbf{Tr}\left[\widehat{\rho }\widehat{X}\right]=\int \lambda \mathbf{Tr}\left[\widehat{\rho }\widehat{E}\left(\lambda \right)\right]d\lambda$ (p. 649 in Ref. [8]).

To understand the transition from classical to quantum theory of the mathematical model used for describing statistical mechanics, as just mentioned, Accardi’s viewpoint [9,10] is advisable: “In analogy with geometry, one should look at probability theory not as the study of the laws of chance but of the possible, mutually inequivalent models for the laws of chance, the choice among which for the description of natural phenomena, being a purely experimental question” (For a physical understanding of quantum-to-classical transition, see Ref. [11]). A crucial point should be emphasized here: the laws of chance for studying indeterministic (statistical) phenomena are as important as the laws of dynamics for studying deterministic (causal) motion of a physical system because they are two different kinds of laws in nature. The former is intrinsically nonlocal, while the latter only describes the locally connected phenomena.

We know that, in classical mechanics, it is well founded that Euclidean spacetime, as the conventional kinetic arena of Newtonian mechanics, has to be substituted with Minkowski spacetime in relativistic mechanics. In geometry, both Newtonian spacetime and Minkowski spacetime are treated as special cases of a Riemanian manifold. In particular, the former has a positive definite metric, while the latter has an indefinite metric. By analogy, it is natural to ask whether it is possible to construct a mathematical model of the probability theory for studying statistical mechanics such that both the Kolmogorovian model and von Neumann’s model can be dealt with in a unified mathematical framework. Indeed, a generalized version of von Neumann’s framework for quantum mechanics has been proposed by Segal [12], who proposed to consider the ${\mathcal{C}}^{\ast }$-algebra generated by all bounded observables of a quantum system as the primary subject of this theory. In this abstract algebraic formalism, the notion of Hilbert space only plays a secondary role as the representation space of the algebra, and the states of a quantum statistical ensemble (i.e., density matrices in the Hilbert space) are represented by normalized positive linear functionals on the algebra so that the probability space of a quantum system can be represented with $(\mathcal{A},\omega (a\left)\right)$, where $\mathcal{A}$ is a ${\mathcal{C}}^{\ast }$-algebra and $\omega \left(a\right)$ is a normalized positive linear functional on $\mathcal{A}$. In fact, using algebraic formalism, the algebraic structure of observables provides a useful framework for analyzing quantum statistical mechanics, particularly for a system in equilibrium state [13], as well as for the study of quantum field theory [14]. Moreover, the fact that one can handle all inequivalent representations simultaneously in the algebraic approach is a manifestation of its power superior to von Neumann’s Hilbert-space approach. In addition, the algebraic approach also allows us to deal with commutative and noncommutative observable algebra within a single mathematical framework.

However, for studying quantum statistical mechanics, the mathematical framework based on the ${\mathcal{C}}^{\ast }$-algebraic approach has a serious disadvantage, i.e., the disability to involve unbounded observables ubiquitous in quantum mechanics. As such, we have to relax the topological constraints imposed on the observables and use ∗-algebra instead of ${\mathcal{C}}^{\ast }$-algebra to represent their algebraic structure, as some mathematicians have tried to develop the theory of quantum probability [15]. However, the algebraic approach of quantum probability, though very powerful mathematically, adds little to physical understanding compared with the conventional formulation of quantum theory. Thus, it is desirable to search for natural operational principles that may lead to a specific set of “preferred observables” important for a physical understanding. In this connection, the conventional algebraic (as well as von Neumann’s) approach has ignored a priori the issue of the “preferred observables”.

In this paper, we study a certain class of algebraic probability spaces denoted as $(\mathcal{A},\omega )$, with $\mathcal{A}$ being a group algebra. We find that when $\mathcal{A}$ is generated from a Lie group G with a bi-invariant Haar measure—its Lie algebra consists of the set of “preferred observables” of a quantum system—the characteristic function extensively used in the classical probability theory can be extended to the group-algebraic probability space as a normalized non-negative definite function in the group G. Furthermore, we find that this group-theoretical characteristic function (GCF) can be adopted to replace the density matrix commonly used in quantum mechanics to represent a state of the quantum ensemble, and more than that, this GCF representation has some significant advantages for applications in quantum statistical mechanics. It is worth noting that, importantly, as shown in [16], where this group-theoretical formalism was first proposed, certain properties of a quantum system are actually rooted in the specific representation of its phase group, i.e., group G associated with the group-algebraic probability space on the state vector space of the system. Our study therefore suggests that it is possible to obtain an important physical understanding of a quantum system through its algebraic structure.

The remaining sections of this paper are organized as follows. First of all, the basics of the algebraic formulation of quantum probability will be briefly introduced in Section 2. Then, in Section 3, we will discuss in detail the application of the group-algebraic formalism in statistical mechanics. Next, In Section 4, two examples will be presented to show how the GCF representation may facilitate solving the evolution of a quantum system. Finally, we will summarize our study in Section 5. The mathematical analysis of the group-algebraic formulation of the quantum probability theory, while essential to our discussion, is more technical and will therefore be detailed in Appendix A.

2. Algebraic Probability Spaces

Let $\mathcal{A}$ be a ∗-algebra, i.e., an algebra over the set of complex numbers $\mathbb{C}$ with involution ∗ and unit $1_{\mathcal{A}}$. By definition, a linear functional $\omega \left(a\right)$ of $a\in \mathcal{A}$, taking values in $\mathbb{C}$, is (i) positive if $\omega \left(a^{\ast }a\right)\ge 0$, $\forall a\in \mathcal{A}$; (ii) normalized if $\omega \left(1_{\mathcal{A}}\right)=1$; and (iii) a state if $\omega$ is positive and normalized. If $\omega$ is a state, we have

$\omega \left(a^{\ast }\right)=\overline{\omega \left(a\right)}\mathrm{and}{\left|\omega \left(a^{\ast }b\right)\right|}^2\le \omega \left(a^{\ast }a\right)\omega \left(b^{\ast }b\right),\forall a,b\in \mathcal{A}.$

We denote the space of all linear functionals of $\mathcal{A}$ by ${\mathcal{A}}^{\prime }$ and the set of all states by $\mathcal{S}$. Thus $\mathcal{S}$ is a convex subset of ${\mathcal{A}}^{\prime }$ and its extreme points are called “pure states”.

The pair $(\mathcal{A},\omega )$ for $\omega \in \mathcal{S}$ is called an algebraic probability space, where a real random variable is represented by a real element $a=a^{\ast }\in \mathcal{A}$ and $\omega \left(a\right)$ is interpreted as the expectation value of the random variable a in the state $\omega$.

It should be remarked [15,17] that the identification of $\omega \left(a\right)$ with the expectation value of a random variable $a\in \mathcal{A}$ seems quite formal so far because it is not obvious what the relevant probability distribution (if any) is, which we need to evaluate the expectation value. To show that this is not the case, let us consider a real random variable $a\in \mathcal{A}$ and its moment sequence $m_k=\omega \left(a^k\right)$, where $k=0,1,\dots ,N$. As for any set $\{c_k\in \mathbb{C}$; $k=0,1,\dots ,N\}$ and $b={\sum }_{k=0}^N{c}_k{a}^k$ we have

$\omega \left(b^{\ast }b\right)=\sum _{j,k}^N\overline{c_j}{c}_k{m}_{j+k}\ge 0,$

Then, according to Hambuger’s theorem for the classical moment problem [18], there exists at least one probability measure ${\mu }_{\omega }^a$ such that

(1)$m_k=\omega \left(a^k\right)={\int }_{-\infty }^{\infty }{x}^k{\mu }_{\omega }^a\left(dx\right).$

Furthermore, if function ${\sum }_k{m}_k{z}^k/k!$ is analytic in the neighborhood of $z=0$, the probability measure ${\mu }_{\omega }^a$ is uniquely determined.

While this approach using moments is useful for finding a suitable probability measure satisfying Equation (1) for a single random variable a, it is frustrating when we try to generalize it to the case with a set of pairwise commuting random variables for finding their joint probability measure [18]. The conventional way out of this difficulty is to use the Hilbert-space representation of the ∗-algebra $\mathcal{A}$. In this way, first, we need to construct a representation $\mathcal{A}\ni a\to \widehat{\pi }\left(a\right)$ on a Hilbert space $\mathcal{H}$ such that a real element a is mapped to a selfadjoint operator $\widehat{\pi }\left(a\right)$ and each density matrix $\widehat{\rho }$ in $\mathcal{H}$ is represented by a state ${\omega }_{\rho }\in \mathcal{S}$ satisfying ${\omega }_{\rho }\left(a\right)=\mathbf{Tr}\left[\widehat{\rho }\widehat{\pi }\left(a\right)\right]$. Then, we have [19]

$\widehat{\pi }\left(a\right)=\int xd{\widehat{P}}^a\left(x\right),{\mu }_{\rho }^a\left(dx\right)=\mathbf{Tr}\left[\widehat{\rho }d{\widehat{P}}^a\left(x\right)\right],\mathrm{and}{\omega }_{\rho }\left(a\right)={\int }_{-\infty }^{\infty }x{\mu }_{\rho }^a\left(dx\right),$

where $d{\widehat{P}}^a\left(x\right)$ is the projection valued measure (PVM) of the self-adjoint operator $\widehat{\pi }\left(a\right)$.

In fact, when the ∗-algebra $\mathcal{A}$ in the probability space $(\mathcal{A},\omega )$ is a ${\mathcal{C}}^{\ast }$-algebra (a normed ∗-algebra with $\parallel a^{\ast }{a\parallel =\parallel a\parallel }^2$, $\forall a\in \mathcal{A}$), we have the following GNS (after Gelfand, Neumark, and Segal) construction theorem (p. 85 in Ref. [20]):

Theorem 1.

Let $\mathcal{A}$ be a unital ${\mathcal{C}}^{\ast }$-algebra and ω be a pure state of $\mathcal{A}$. There exists a triple $({\mathcal{H}}_{\omega },|{\psi }_{\omega }⟩,{\pi }_{\omega })$, where ${\mathcal{H}}_{\omega }$ is a Hilbert space, $|{\psi }_{\omega }⟩\in {\mathcal{H}}_{\omega }$ is a unit vector, and ${\pi }_{\omega }$: $\mathcal{A}\to \mathfrak{B}\left({\mathcal{H}}_{\omega }\right)$ is a ∗-homomorphism from $a\in \mathcal{A}$ to a bounded operator ${\widehat{\pi }}_{\omega }\left(a\right)$ in ${\mathcal{H}}_{\omega }$, such that (i) $|{\psi }_{\omega }⟩$ is a cyclic vector, namely, ${\widehat{\pi }}_{\omega }\left(\mathcal{A}\right)|{\psi }_{\omega }⟩$ is dense in ${\mathcal{H}}_{\omega }$; (ii) $⟨{\psi }_{\omega }\left|{\widehat{\pi }}_{\omega }\left(a\right)\right|{\psi }_{\omega }⟩=\omega \left(a\right),\forall a\in \mathcal{A}$; and (iii) the representation $\mathcal{A}\ni a\to {\widehat{\pi }}_{\omega }\left(a\right)$ is irreducible.

We call this representation of ∗-algebra, $\mathcal{A}$, the cyclic representation associated with the state $\omega$ (or the cyclic vector $|{\psi }_{\omega }⟩$). Noting that all bounded symmetric operators are self-adjoint and any non-zero vector in ${\mathcal{H}}_{\omega }$ is cyclic with its associated representation unitarily equivalent to $\mathcal{A}\ni a\to {\widehat{\pi }}_{\omega }\left(a\right)$, therefore, when $\mathcal{A}$ is a ${\mathcal{C}}^{\ast }$-algebra, GNS construction provides a Hilbert-space representation $\mathcal{A}\ni a\to {\widehat{\pi }}_{\omega }\left(a\right)$ such that any real element $a=a^{\ast }\in \mathcal{A}$ is mapped to a self-adjoint operator and any density matrix $\widehat{\rho }$ in ${\mathcal{H}}_{\omega }$ can be represented by a state ${\omega }_{\rho }\left(a\right)=\mathbf{Tr}\left[\widehat{\rho }{\widehat{\pi }}_{\omega }\left(a\right)\right]$.

The GNS construction of a Hilbert-space representation can be generalized to the case where the ∗-algebra is not a $C^{\ast }$-algebra. Let ${\mathcal{D}}_{\omega }$ be a pre-Hilbert space (a dense subspace of a Hilbert space) and L(${\mathcal{D}}_{\omega }$) be the set of linear operators from ${\mathcal{D}}_{\omega }$ into itself which admits adjoints so that L(${\mathcal{D}}_{\omega }$) becomes a ∗-algebra. Let $|{\psi }_{\omega }⟩\in {\mathcal{D}}_{\omega }$ be a unit vector, then there exists a GNS representation ${\pi }_{\omega }$ of $(\mathcal{A},\omega )$ such that (i) ${\pi }_{\omega }$: $\mathcal{A}\to L\left({\mathcal{D}}_{\omega }\right)$ is a ∗-homomorphism; (ii) ${\widehat{\pi }}_{\omega }\left(\mathcal{A}\right)|{\psi }_{\omega }⟩={\mathcal{D}}_{\omega }$; and (iii) $⟨{\psi }_{\omega }\left|{\widehat{\pi }}_{\omega }\left(a\right)\right|{\psi }_{\omega }⟩=\omega \left(a\right),\forall a\in \mathcal{A}$. However, in contrast with the $C^{\ast }$ case, now ${\widehat{\pi }}_{\omega }\left(a\right)$ can be an unbounded operator in ${\mathcal{D}}_{\omega }$ and it is difficult to identify the precise conditions for the unbounded symmetric operator ${\widehat{\pi }}_{\omega }\left(a\right)$ in ${\mathcal{D}}_{\omega }$ being essentially self-adjoint in the sense of the Hilbert space theory (p. 37 and p. 310 in Ref. [17]).

In the following section, we will propose a special class of algebraic probability spaces, $(\mathcal{A},\omega )$, with $\mathcal{A}$ being a group algebra. Our analysis will show that, in this mathematical group-algebraic framework, the quantum generalization of the characteristic functions extensively used in the classical probability theory can be realized in a natural way. Moreover, the representation of the ∗-algebra, $\mathcal{A}$, can be generated directly from the representation of the phase group G without using the more sophisticated GNS technique. In addition, for a statistical system described with the group-algebraic probability space, its states may have some remarkable features relevant to the group representation.

3. Algebraic Probability Spaces as the Mathematical Foundation of Classical and Quantum Statistical Mechanics

The primary subject for statistical theory of classical mechanics is to address the microscopic motion of molecules in gases. The kinetic theory established to this end appeared in the middle of the nineteenth century, featuring probability that is completely extraneous to the deterministic Newtonian mechanics. It was brought to a mature stage by Clausius, Maxwell, and Boltzmann, with explicit recognition that the random motion of molecules can be characterized with probability distributions of the positions and velocities of individual molecules, and thermal energy is nothing but the kinetic energy of the microscopic motion.

Later, an important contribution to the mathematical framework of classical statistical mechanics was made by Gibbs. He introduced the concept of “statistical ensemble” in the phase space spanned by pairs of canonical variables, in contrast with probability distributions of positions and velocities. Taking Kolmogorov’s measure-theoretical approach to probability theory, Gibbs adopted the phase space as the sample space of random variables. The advantage is that the volume element of the phase space is invariant under the dynamical evolution of the system, and hence it owns the property of an invariant measure.

It was totally unanticipated until the 1930s that probability would also play a crucial role in quantum mechanics deemed to describe the microscopic world. The statistical interpretation of the wave function of a pure quantum state was first introduced by Born [6] and this interpretation was further developed by von Neumann into a mathematically rigorous theory by recognizing that states and observables are the most essential concepts of quantum mechanics. In von Neumann’s theory, states and observables of a quantum system are represented with state vectors (rays) and self-adjoint operators, respectively, in a Hilbert space, while the square of the modulus of the inner product of state vectors provides the probability introduced by Born. Furthermore, von Neumann also pointed out that for establishing the probability theory as the theory of frequencies, it is necessary to introduce the concept of quantum statistical ensemble represented by density matrices in the Hilbert space (p. 298 in Ref. [7]).

Note that although it is still a subject of debate if the quantum probability model proposed by von Neumann can be regarded as a special case of Kolmogorov’s measure-theoretical probability model, now it has been well established, with both physically verified quantum theory and rigorous mathematical criteria, that in nature, one can find some sets of data concerning mutually incompatible sets of observables of quantum systems that can be described with von Neumann’s model rather than Kolmogorov’s model [9].

As mentioned in the Introduction, a remarkable generalization and extension of von Neumann’s mathematical framework is due to Segal [12], who proposed to consider as the primary subject, the ${\mathcal{C}}^{\ast }$-algebra generated by bounded observables in the mathematical framework of quantum mechanics. However, for studying quantum statistical mechanics, the mathematical framework based on this algebraic approach suffers from the following problems: (i) Unbounded observables are common in quantum mechanics. As such, we have to relax the topological constraints and use ∗-algebra instead of ${\mathcal{C}}^{\ast }$-algebra to represent the algebraic structure of observables, so that the statistical behavior of a quantum system can be described with the probability space outlined in Section 2; (ii) For all specific quantum theories, they commonly proceed with the input of a classical theory that tells what the relevant observables are and how they are related. In this connection, the algebraic approach (as well as von Neumann’s approach) dispenses, a priori, with the choice of the “preferred observables” altogether; (iii) The Poisson structure of phase space plays a crucial role in classical statistical mechanics, but the corresponding specific structure in the algebraic approach of the mathematical framework of quantum statistical mechanics has not been explicitly explored yet.

To solve problems (ii) and (iii), more restrictive assumptions are needed to impose on the specific formalism of the algebraic probability space. There are various possible routes to this end. One was proposed by one of the authors in Ref. [21], where it was suggested to restrict the ∗-algebra of an algebraic probability space to a specific group algebra ${\mathcal{A}}_G$ (denoted by ${\mathcal{D}}^{\prime }\left(G\right)$ for compact Lie groups and ${\mathcal{E}}^{\prime }\left(G\right)$ for locally compact Lie groups, respectively, in Appendix A). The resultant formalism was referred to “the group-theoretical formalism” [16]. Here, we will review, discuss in more detail, and improve this formalism by laying it on a more solid mathematical foundation.

Let us assume that the quantum system of interest is characterized with a basic set of independent dynamical variables $\{X_1,X_2,\dots ,X_n\}$ that forms a basis of a Lie algebra $\mathfrak{g}$. We denote by G the Lie group generated by $\mathfrak{g}$ and ${\mathcal{A}}_G$ as the corresponding group algebra. Then, the statistical behavior of this quantum system can be described based on the algebraic probability space associated with the group algebra ${\mathcal{A}}_G$.

For the applications of this group-theoretical formalism in studying quantum statistical mechanics, we would like to emphasize the following notable aspects.

(a) Determining the “preferred observables” of the system.

Let X be an element of the Lie algebra $\mathfrak{g}$ and $p\left(X\right)\in {\mathcal{A}}_G$ be a real polynomial function of X. For an arbitrary unitary representation of G in a Hilbert space $\mathcal{H}$, the element $p\left(iX\right)\in {\mathcal{A}}_G$ will be mapped to a self-adjoint operator in the Hilbert space $\mathcal{H}$. In this case, it is natural to assume that the set of dynamical variables $\{i{X}_1,i{X}_2,\dots ,i{X}_n\}$ play the role of “preferred observables” of the system.

(b) Investigating the evolution of the system.

For the case where the relevant Lie group G admits a bi-invariant Haar measure $dg$, let $G\ni g\to \widehat{U}\left(g\right)$ be an irreducible unitary representation in a Hilbert space $\mathcal{H}$ and $\widehat{\rho }$ a density matrix in $\mathcal{H}$. Then, for each $\widehat{\rho }$ in $\mathcal{H}$, there exists a group-theoretical characteristic function (GCF) ${\phi }_{\rho }\left(g\right)=\mathrm{Tr}\left[\widehat{\rho }\widehat{U}\left(g\right)\right]$ with the well-defined inverse transform $\widehat{\rho }={\int }_G{\phi }_{\rho }\left(g\right)\widehat{U^{†}}\left(g\right)dg$ (see Appendix A). In other words, after giving an irreducible unitary representation of G in $\mathcal{H}$, the state of the quantum system, conventionally represented with the density matrix $\widehat{\rho }$ in $\mathcal{H}$, can be represented with the corresponding GCF, ${\phi }_{\rho }$. Since ${\phi }_{\rho }\left(g\right)$ is a c-number (infinitely differentiable function on G), it is advisable to represent the dynamical evolution equation of the system as differential equations of ${\phi }_{\rho }\left(g\right)$ instead, as it is generally easier to manipulate than the corresponding equation in terms of the density matrix, especially when the dimension of the representation space $\mathcal{H}$ is large. In addition, by using the GCF representation, we can handle all inequivalent representations of the group algebra ${\mathcal{A}}_G$ simultaneously in the same mathematical approach.

(c) Dealing with composite quantum systems.

Let $S_1$ and $S_2$ be two quantum systems described with probability spaces associated with group algebras, ${\mathcal{A}}_{G_1}$ and ${\mathcal{A}}_{G_2}$, respectively. Consider the composite system $S=S_1+S_2$ described with the probability space $({\mathcal{A}}_G,\phi \left(g\right))$, where $G=G_1⨂G_2$ and $\phi \left(g\right)$ is a state on ${\mathcal{A}}_G$. The reduced state $\phi \left(g\right)$ for subsystems $S_1$ and $S_2$ is now simply given by ${\phi }_1\left(g_1\right)=\phi (g_1⨂e_2)$ and ${\phi }_2\left(g_2\right)=\phi (e_1⨂g_2)$, respectively, where $g_j\in G_j$ and $e_j$ is the identity of $G_j$ for $j=1,2$. Furthermore, in group-theoretical formalism, the state of the composite system $\phi \left(g\right)$ is separable with respect to subsystems $S_1$ and $S_2$ if and only if $\phi \left(g\right)={\sum }_{k=1}^n{p}_k{\kappa }_k\left(g_1\right){\omega }_k\left(g_2\right)$, where $p_k>0$, ${\sum }_{k=1}^n{p}_k=1$, and $\{{\kappa }_k\left(g_1\right),{\omega }_k\left(g_2\right),k=1,\dots ,n\}$ are states on ${\mathcal{A}}_{G_1}$ and ${\mathcal{A}}_{G_2}$, respectively [22].

(d) Revealing the quantum-to-classical transition.

Let $\mathfrak{g}$ be the Lie algebra of a connected Lie group G. We take its dual vector space ${\mathfrak{g}}^{\ast }$ as the phase space of a classical mechanical system. Let $\mathcal{E}\left({\mathfrak{g}}^{\ast }\right)$ be the space of all $C^{\infty }$ functions on ${\mathfrak{g}}^{\ast }$. It is known that $\mathcal{E}\left({\mathfrak{g}}^{\ast }\right)$ admits a Poisson structure defined by linear Poisson brackets [23]:

(2)${\{f,k\}}_{\xi }=\sum _{\alpha ,\beta }(\xi ,[X_{\alpha },X_{\beta }]){\displaystyle \frac{\partial f}{\partial {\xi }_{\alpha }}}{\displaystyle \frac{\partial k}{\partial {\xi }_{\beta }}}=\sum _{\alpha ,\beta ,\gamma }{C}_{\alpha ,\beta }^{\gamma }{\xi }_{\gamma }{\displaystyle \frac{\partial f}{\partial {\xi }_{\alpha }}}{\displaystyle \frac{\partial k}{\partial {\xi }_{\beta }}},f,k\in \mathcal{E}\left({\mathfrak{g}}^{\ast }\right),$

where $\xi \in {\mathfrak{g}}^{\ast },\{X_{\alpha },\alpha =1,\dots ,n\}$ is a basis of $\mathfrak{g}$, and ${\xi }_{\alpha }=(\xi ,X_{\alpha })$ and $C_{\alpha ,\beta }^{\gamma }$ are the structure constants of $\mathfrak{g}$. The Lie–Poisson phase space ${\mathfrak{g}}^{\ast }$ can be foliated into a family of leaves, i.e., coadjoint orbits of G (p. 226 in Ref. [24]), where each leaf is a G-invariant symplectic manifold (p. 20 in Ref. [25]). These symplectic leaves are called superselection sectors of the phase space ${\mathfrak{g}}^{\ast }$.

To discuss the applications of the group-algebraic formalism, in the following, we will focus on the quantum systems of canonical variables. First, we show that the statistical behavior of these systems can be described using the algebraic probability spaces associated with the Heisenberg–Weyl group [26]. To this end, let us consider a quantum system with canonical variables $(q,p)$ whose commutation relation is $[\widehat{q},\widehat{p}]=i\hslash$. This commutation relation can be regarded as the first of the following three pairs of Lie brackets of a three-dimensional Lie algebra $\mathfrak{g}$ with basis $\{X_1,X_2,X_3\}$

$[X_1,X_2]=X_3,[X_1,X_3]=0,\mathrm{and}[X_2,X_3]=0,$

by identifying $X_1$, $X_2$, and $X_3$ with $\widehat{q}$, $\widehat{p}$, and $i\hslash$, respectively. Exponentiating the Lie algebra $\mathfrak{g}$, we obtain the Heisenberg–Weyl group $H\left(1\right)$, of which each element $g\in H\left(1\right)$ can be written as $exp(\mathbf{a}·X)$ with $\mathbf{a}\in {\mathbb{R}}^3$. Moreover, the multiplication law of group elements is given by $exp(\mathbf{a}·X)exp(\mathbf{b}·X)=exp(\mathbf{c}·X)$, with

(3)$\mathbf{c}=(a_1+b_1,a_2+b_2,a_3+b_3-{\displaystyle \frac{1}{2}}(a_1{b}_2-a_2{b}_1)).$

Now, let us turn to the classical systems associated with the Heisenberg–Weyl group $H\left(1\right)$. Let ${\mathfrak{g}}^{\ast }$ be the dual space of $\mathfrak{g}$, $\xi \in {\mathfrak{g}}^{\ast }$, and ${\xi }_j=(\xi ,X_j)$. Noting that the coadjoint action of $g=exp(\mathbf{a}·X)\in H\left(1\right)$ on the vector $\xi \in {\mathfrak{g}}^{\ast }$ is infinitesimally generated by

$A{d}_{\mathbf{a}·X}^{\ast }\xi =(\xi ,[X,\mathbf{a}·X])=(a_2{\xi }_3,-a_1{\xi }_3,0),$

we can conclude that there are two classes of coadjoint orbits in ${\mathfrak{g}}^{\ast }$: (i) the two-dimensional planes perpendicular to axis ${\xi }_3$ when ${\xi }_3\ne 0$ and (ii) all points on the plane of ${\xi }_3=0$. According to Equation (2), the Poisson brackets on the two-dimensional coadjoint orbits can be written as

${\{f,k\}}_{\xi }=({\displaystyle \frac{\partial f}{\partial {\xi }_1}}{\displaystyle \frac{\partial k}{\partial {\xi }_2}}-{\displaystyle \frac{\partial f}{\partial {\xi }_2}}{\displaystyle \frac{\partial k}{\partial {\xi }_1}}){\xi }_3,f,k\in \mathcal{E}\left({\mathfrak{g}}^{\ast }\right).$

Substituting ${\xi }_1$ and ${\xi }_2$ with q and p, respectively, and considering the special symplectic leaf ${\xi }_3=1$, we immediately obtain the standard expressions of Poisson brackets, the crucial structure of classical mechanics.

(e) Quantization of classical systems of compact phase spaces by constructing the corresponding algebraic probability space.

The compactness of the phase space implies that, for the corresponding quantized system, the number of phase cells in the phase space, or equivalently, the dimension of the state vector space $\mathcal{H}$, is a finite integer. As an illustrating example, let us consider a system whose degree of freedom is one and the phase space is a two-dimensional torus. Without loss of generality, we suppose that the torus has a unit area with a unit period in both dimensions. As such, the dimension of $\mathcal{H}$ is $N=1/\left(2\pi \hslash \right)$. Following the scheme suggested in Ref. [27], the first step of quantization is to construct a pair of conjugate orthonormal bases in the vector space $\mathcal{H}$, $\left\{\right|j⟩\}$ and $\left\{\right|\overline{l}⟩\}$, with $j,l=0,\dots ,N-1$, such that

$|\overline{l}⟩={\displaystyle \frac{1}{\sqrt{N}}}\sum _{j=0}^{N-1}{e}^{i2\pi jl/N}|j⟩,|j⟩={\displaystyle \frac{1}{\sqrt{N}}}\sum _{l=0}^{N-1}{e}^{-i2\pi jl/N}|\overline{l}⟩.$

Next, we can define a pair of canonical variables $(\widehat{q},\widehat{p})$ in $\mathcal{H}$, by assuming that $\left\{\right|j⟩\}$ and $\left\{\right|\overline{l}⟩\}$ are their eigenstates, respectively, i.e., $\widehat{q}|j⟩={\displaystyle \frac{j}{N}}|j⟩,j=0,\dots ,N-1$ and $\widehat{p}|\overline{l}⟩={\displaystyle \frac{l}{N}}|\overline{l}⟩,l=0,\dots ,N-1$. It is worth noting that the the operator pair $(\widehat{q},\widehat{p})$ thus defined does not obey the conventional Heisenberg commutation relation for the flat two-dimensional phase space; rather, it follows the Weyl commutation rule:

$e^{i2\pi n\widehat{p}}{e}^{i2\pi n\widehat{q}}=e^{i2\pi nm/N}{e}^{i2\pi n\widehat{q}}{e}^{i2\pi n\widehat{p}},(m,n)\in {\mathbb{Z}}^2.$

Now, let us consider a compact Lie group G, any element of which $g\in G$ is represented with three parameters, $\{m,n,z\}$, where $(m,n)\in {\mathbb{Z}}^2$, $0\le m,n<N$, and $0\le z<2\pi$. The multiplication law is

$\{m_1,n_1,z_1\}\{m_2,n_2,z_2\}=\{m_1+n_1,m_2+n_2,z_1+z_2+{\displaystyle \frac{\pi }{N}}(m_2{n}_1-m_1{n}_2)\}.$

Comparing this multiplication law with that of the Heisenberg–Weyl group $H\left(1\right)$ (see Equation (3)), the group G can be regarded as a discretized version of $H\left(1\right)$. Noting that the Weyl commutation relation gives a projective representation of G in the state vector space $\mathcal{H}$, we can therefore lift this projective representation to an irreducible unitary representation $G\ni g\to \widehat{U}\left(g\right)=e^{i(z+\pi nm/N)}{e}^{i2\pi m\widehat{q}}{e}^{i2\pi n\widehat{p}}$. Thus, using the quantization procedure mentioned above, the behavior of the quantized system can be described with the algebraic probability space associated with the group algebra ${\mathcal{A}}_G$.

4. Application Examples

As discussed in the previous section, the group-theoretical formalism provides not only a more self-consistent platform conceptually to address the problems of quantum statistical mechanics, but also one more optional tool for solving the problems technically. Indeed, in some cases, it may greatly facilitate calculations. In the following part of this section, we will illustrate its superiority with two examples. One is a harmonic oscillator coupled with a thermal environment; its relaxation process to the equilibrium state will be discussed. Another is the cat map, a well-known paradigm for chaos. We will discuss its relaxation process characterized by the coarse-grained entropy.

The group-theoretical formalism can be advantageous in other situations as well. For example, it has been found to be much more convenient than the conventional methods for solving the problem of quantum Brownian motion. See Ref. [28] for a detailed study based on the Caldeira–Leggett master equation and its extension to the Lindblad type.

4.1. The Harmonic Oscillator in a Thermal Environment

The evolution of a quantum system that interacts with its environment is of fundamental importance in the theory of open quantum systems [29,30]. A harmonic oscillator moving irreversibly in its thermostatic surroundings is presumably the simplest model for probing this issue. The reduced density matrix $\rho$ of the oscillator as a function of time is assumed to follow the quantum-optical master equation

(4)${\displaystyle \frac{d\rho }{dt}}=-i\omega [a^{†}a,\rho ]+\Gamma (N_{\beta }+1)(2a\rho a^{†}-\{a^{†}a,\rho \})+\Gamma N_{\beta }(2a^{†}\rho a-\{a{a}^{†},\rho \}).$

Here, $\Gamma$ is the damping constant, $\beta =1/\left(k_{B}T\right)$ with $k_B$ and T being the Boltzmann constant and the environmental temperature, respectively, $N_{\beta }=1/(e^{\hslash \beta \omega }-1)$ is the expected number of thermal quanta, and

$a={\displaystyle \frac{1}{\sqrt{2\hslash M\omega }}}(M\omega q+ip),a^{†}={\displaystyle \frac{1}{\sqrt{2\hslash M\omega }}}(M\omega q-ip).$

This master equation was first derived by Louisell [31] for addressing a damped mode of the radiation field in a cavity and was treated by Haake [32] as an example of the generalized master equation constructed by Nakajima and Zwanzig. It plays a central role not only in the theoretical interpretations of many quantum-optical experiments [33,34], but also in revealing the environment-induced decoherence of open quantum systems [29].

In the following, we will show how to solve Equation (4) with the GCF representation and scrutinize the relaxation process of a pure quantum state towards a mixed state. For the pair of canonical variables $(q,p)$, as discussed in Section 3 (d), their commutation relation $[q,p]=i\hslash$ can be seen as a pair of Lie brackets of a three-dimensional Heisenberg-Weyl Lie algebra. Denote the nilpotent Lie group generated by this Lie algebra as $H\left(1\right)$ and let

$g\to \mathbf{U}\left(g\right)=e^{i(xp+yq+z\hslash )},g\in H\left(1\right)$

be an irreducible unitary representation of $H\left(1\right)$ on the state vector space $\mathfrak{H}$ of the system, where $\{(x,y)\in {\mathbb{R}}^2$, $0\le z<2\pi /\hslash$} represents a canonical system of coordinates on the group $H\left(1\right)$, the GCF of the density matrix $\rho \left(t\right)$ of the system at time t is (see Section 3 (b) and Appendix A)

$\phi (t,g)={\mathrm{Tr}}_{\mathfrak{H}}\left[\rho \left(t\right)\mathbf{U}\left(g\right)\right]=e^{iz\hslash +v(x,y,t)},$

where $v(x,y,t)$ is a complex-number function. Consequently, the inverse mapping has the form

$\rho \left(t\right)={\int }_{H\left(1\right)}\phi (t,g){\mathbf{U}}^{†}\left(g\right)dg,$

where $dg={(\hslash /2\pi )}^{2}dxdydz$ is the bi-invariant Haar measure on $H\left(1\right)$. If we use complex variables ($a,a^{†}$) instead of $(q,p)$, the unitary representation of $H\left(1\right)$ in $\mathfrak{H}$ has the form $\mathbf{U}\left(g\right)=e^{u{a}^{†}-\overline{u}a+iz\hslash }$ with $u=u_1+i{u}_2=-\sqrt{{\displaystyle \frac{\hslash M\omega }{2}}}x+i\sqrt{{\displaystyle \frac{\hslash }{2M\omega }}}y$, which gives an alternative expression for the GCF

$\phi (t,g)={\mathrm{Tr}}_{\mathfrak{H}}\left[\rho \left(t\right)\mathbf{U}\left(g\right)\right]=e^{iz\hslash +v(u,t)}.$

Using formulae $a\mathbf{U}\left(g\right)=(-{\displaystyle \frac{\partial }{\partial \overline{u}}}+{\displaystyle \frac{u}{2}})\mathbf{U}\left(g\right)$ and $a^{†}\mathbf{U}\left(g\right)=({\displaystyle \frac{\partial }{\partial u}}+{\displaystyle \frac{\overline{u}}{2}})\mathbf{U}\left(g\right)$, the quantum-optical master Equation (4) can be converted into the following equation for $v(u,t)$:

(5)${\displaystyle \frac{\partial v(u,t)}{\partial t}}=-\{[(\Gamma u_1+\omega u_2){\displaystyle \frac{\partial }{\partial u_1}}+(\Gamma u_2-\omega u_1){\displaystyle \frac{\partial }{\partial u_2}}]v\left(u\right)+(1+2N_{\beta }){\Gamma \left|u\right|}^2\},$

or, alternatively, for $v(x,y,t)$

(6)${\displaystyle \frac{\partial v(x,y,t)}{\partial t}}=[({\displaystyle \frac{y}{M}}-\Gamma x){\displaystyle \frac{\partial }{\partial x}}-(M{\omega }^{2}x+\Gamma y){\displaystyle \frac{\partial }{\partial y}}]v(x,y)-{\displaystyle \frac{\hslash }{2}}(1+2N_{\beta })\Gamma (M\omega x^2+{\displaystyle \frac{y^2}{M\omega }}).$

In the following, we discuss the solution to these two equations with two different initial conditions, respectively.

(a) Case 1: An initial Gaussian state

When the initial state of the system is a Gaussian state centered at $(\overline{q},\overline{p})$, i.e.,

$⟨q|{\psi }_G⟩={\left({\displaystyle \frac{1}{2\pi {\sigma }^2}}\right)}^{{\displaystyle \frac{1}{4}}}{e}^{-{\displaystyle \frac{{(q-\overline{q})}^2}{4{\sigma }^2}}+{\displaystyle \frac{i}{\hslash }}\overline{p}(q-\overline{q})},$

the GCF of the initial density matrix $\rho \left(0\right)$ reads $\phi (0,g)=⟨{\psi }_G\left|\mathbf{U}\left(g\right)\right|{\psi }_G⟩=e^{iz\hslash +v(x,y,0)}$, where

$v(x,y,0)=i(x\overline{p}+y\overline{q})-({\displaystyle \frac{{\hslash }^2{x}^2}{8{\sigma }^2}}+{\displaystyle \frac{{\sigma }^2{y}^2}{2}}).$

Then, it is straightforward to write down the solution of the master Equation (6) as

$v(x,y,t)=A\left(t\right)x+B\left(t\right)y+a\left(t\right)x^2+b\left(t\right)y^2+c\left(t\right)xy$

with

$\left\{\begin{array}{c}A\left(t\right)=i{e}^{-\Gamma t}(\overline{p}cos\omega t-M\overline{q}\omega sin\omega t),\hfill \\ B\left(t\right)=i{e}^{-\Gamma t}(\overline{q}cos\omega t+{\displaystyle \frac{\overline{p}}{M\omega }}sin\omega t),\hfill \\ a\left(t\right)=-{\displaystyle \frac{\hslash M\omega }{4}}(2N_{\beta }+1)(1-e^{-2\Gamma t})-{\displaystyle \frac{1}{2}}{e}^{-2\Gamma t}({\displaystyle \frac{{\hslash }^2}{4{\sigma }^2}}{cos}^2\omega t+M^2{\sigma }^2{\omega }^2{sin}^2\omega t),\hfill \\ b\left(t\right)=-{\displaystyle \frac{\hslash }{4M\omega }}(2N_{\beta }+1)(1-e^{-2\Gamma t})-{\displaystyle \frac{1}{2}}{e}^{-2\Gamma t}({\sigma }^2{cos}^2\omega t+{\displaystyle \frac{{\hslash }^2}{4M^2{\sigma }^2{\omega }^2}}{sin}^2\omega t),\hfill \\ c\left(t\right)={\displaystyle \frac{e^{-2\Gamma t}}{8M{\sigma }^2\omega }}(4M^2{\sigma }^4{\omega }^2-{\hslash }^2)sin2\omega t.\hfill \end{array}\right.$

In the GCF representation, the purity of the system has the form [16]

$P\left(t\right)=\mathrm{Tr}\left[\rho {\left(t\right)}^2\right]={\int }_{H\left(1\right)}{\left|\phi (t,g)\right|}^{2}dg={\displaystyle \frac{\hslash }{2\pi }}\int dx\int dy{e}^{2\mathbf{Re}\left[v\right(x,y,t\left)\right]},$

and thus the linear entropy [35] of the system is

(7)$\begin{array}{c}\hfill S\left(t\right)=-lnP\left(t\right)=ln{\displaystyle \frac{2}{\hslash }}\sqrt{4a\left(t\right)b\left(t\right)-c{\left(t\right)}^2}\\ \hfill ={\displaystyle \frac{1}{2}}ln[e^{-4\Gamma t}+{(1-e^{-2\Gamma t})}^2{(2N_{\beta }+1)}^2+2(1-e^{-2\Gamma t})e^{-2\Gamma t}(2N_{\beta }+1)cosh2r],\end{array}$

where the squeeze factor r is defined by $e^{-r}={\displaystyle \frac{\sigma }{{\sigma }_c}}$ with ${\sigma }_c=\sqrt{{\displaystyle \frac{\hslash }{2M\omega }}}$. As shown elsewhere [36], for a classical harmonic oscillator with an initial localized phase-space distribution, the evolution of its entropy depends on neither its environmental temperature nor its initial position in the phase space. Therefore, it is somehow perplexing to comprehend such differences between both quantum and classical dynamical behavior in the context of non-equilibrium statistical mechanics.

From Equation (7), we can see that if the initial state is a coherent state with $r=0$,

(8)$S\left(t\right)=ln[1+2N_{\beta }(1-e^{-2\Gamma t})],$

suggesting that the entropy increases monotonically from $S\left(0\right)=0$ to its equilibrium value $S_{eq}=ln(1+2N_{\beta })$, as expected. However, for an initial squeezed state with $r>0$, the evolution of entropy becomes more complicated. There exists a certain time

$t_m={\displaystyle \frac{1}{2\Gamma }}ln\left[{\displaystyle \frac{2+4N_{\beta }+4N_{\beta }^2-2(1+2N_{\beta })cosh2r}{(1+2N_{\beta })(1+2N_{\beta }-cosh2r)}}\right]$

before which, $S\left(t\right)$ increases in time but after which $S\left(t\right)$ decreases and converges to $S_{eq}$ (see Figure 1). Since $t_m$ increases from ${\displaystyle \frac{ln2}{2\Gamma }}$ to ∞ as $N_{\beta }$ increases from 0 to $N_c={sinh}^{2}r$, such an interesting non-monotonic “hump” characteristic of the time curve of $S\left(t\right)$ appears if and only if $N_{\beta }<N_c$.

(b) Case 2: An initial Fock state

When the initial state is a Fock state, i.e., an energy eigenstate of the isolated harmonic oscillator, denoted as $|n⟩$, it is convenient to take the complex coordinate system on $H\left(1\right)$ instead and consequently, the GCF can be written as

$\phi (0,g)=⟨n\left|\mathbf{U}\left(g\right)\right|n⟩=e^{iz\hslash +{\displaystyle \frac{{\left|u\right|}^2}{2}}}⟨n|e^{-\overline{u}\mathbf{a}}{e}^{u{\mathbf{a}}^{†}}|n⟩=e^{iz\hslash +v_n\left(u\right)}.$

Noting that for the coherent state $|\alpha ⟩$, we have $⟨n|\alpha ⟩={\displaystyle \frac{{\alpha }^n}{\sqrt{n!}}}{e}^{-{\displaystyle \frac{{\left|\alpha \right|}^2}{2}}}$ and as a result,

$⟨n|e^{-\overline{u}a}{e}^{u{a}^{†}}|n⟩={\displaystyle \frac{1}{\pi n!}}\int d^2\alpha {\left|\alpha \right|}^{2n}{e}^{-{\left|\alpha \right|}^2+u\overline{\alpha }-\overline{u}\alpha }$

$={\displaystyle \frac{e^{-{\left|u\right|}^2}}{\pi n!}}\sum _{m=0}^n\left({\displaystyle \genfrac{}{}{0pt}{}{n}{m}}\right)\int d{\alpha }_1{\alpha }_1^{2m}{e}^{-{({\alpha }_1-i{u}_2)}^2}\int d{\alpha }_2{\alpha }_2^{2(n-m)}{e}^{-{({\alpha }_2+i{u}_1)}^2},$

where $\alpha ={\alpha }_1+i{\alpha }_2$. Using the formula ${\int }_{-\infty }^{\infty }dx{x}^n{e}^{-{(x-y)}^2}={\displaystyle \frac{\sqrt{\pi }}{{\left(2i\right)}^n}}{H}_n\left(iy\right)$, we have

(9)$v_n\left(u\right)=v(u,0)=ln\left[{\displaystyle \frac{{(-1)}^n}{n!2^{2n}}}{e}^{-{\displaystyle \frac{1}{2}}{\left|u\right|}^2}\sum _{m=0}^n\left({\displaystyle \genfrac{}{}{0pt}{}{n}{m}}\right)H_{2m}\left(u_2\right)H_{2(n-m)}\left(u_1\right)\right],$

where $H_n\left(x\right)$ is the Hermite polynomial of degree n.

The solution of the master Equation (5) with the initial condition (9) is therefore

(10)$v(u,t)=v_n\left(e^{-\Gamma t}u\right)+{\displaystyle \frac{1}{2}}(1+2N_{\beta })(e^{-2\Gamma t}-1){\left|u\right|}^2$

and the linear entropy of the system reads

$S\left(t\right)=-ln[{\displaystyle \frac{1}{\pi }}\int d{u}_1\int d{u}_2{e}^{2\mathbf{Re}\left[v\right(u,t\left)\right]}].$

Making use of Equation (10) and letting $a^2=1+(1+2N_{\beta })(e^{2\Gamma t}-1)$, the entropy of the system (with initial state $\rho \left(0\right)=|n⟩⟨n|$) can be written explicitly as

(11)$S_n\left(t\right)=-ln\left[{\displaystyle \frac{e^{2\Gamma t}}{{(n!2^{2n})}^2\pi }}\sum _{m_1,m_2}^n\left({\displaystyle \genfrac{}{}{0pt}{}{n}{m_1}}\right)\left({\displaystyle \genfrac{}{}{0pt}{}{n}{m_2}}\right)IP(m_1,m_2)IP(n-m_1,n-m_2)\right],$

where [37]

(12)$\begin{array}{c}\hfill IP(m_1,m_2)=\int dx{e}^{-a^{2}x}{H}_{2m_1}\left(x\right)H_{2m_2}\left(x\right)\\ \hfill =2^{2(m_1+m_2)}\Gamma (m_1+m_2+{\displaystyle \frac{1}{2}}){\displaystyle \frac{{(1-a^2)}^{m_1+m_2}}{a^{2(m_1+m_2)+1}}}{}_{2}F_1(-2m_1,-2m_2,{\displaystyle \frac{1}{2}}-m_1-m_2,{\displaystyle \frac{a^2}{2(a^2-1)}}).\end{array}$

From Equations (11) and (12), we know that

$S_n\left(0\right)=0,S_n{\left(t\right)}_{t\to \infty }=S_{eq}=ln(1+2N_{\beta }).$

We also note that the expression of $S_0\left(t\right)$ given by Equation (11) is the same as that given by Equation (8), as the ground state of a harmonic oscillator itself is a coherent state.

On the other hand, the expressions of $S_n\left(t\right)$ for $n>0$ become more and more complicated as n grows. Taking the case of $n=1$ as an example, we have

$S_1\left(t\right)=ln\left[{\displaystyle \frac{{[e^{-2\Gamma t}+(1-e^{-2\Gamma t})(1+2N_{\beta })]}^3}{e^{-4\Gamma t}+{(1-e^{-2\Gamma t})}^2{(1+2N_{\beta })}^2}}\right].$

To appreciate the role played by the parameter $N_{\beta }$, we consider the extreme points of curve $S_1\left(t\right)$. They are the roots of $d{S}_1\left(t\right)/dt=0$, which are given by

$t_{\pm }={\displaystyle \frac{1}{2\Gamma }}ln\left[{\displaystyle \frac{1+2N_{\beta }-2N_{\beta }^2-4N_{\beta }^3\pm \sqrt{1+6N_{\beta }+10N_{\beta }^2-8N_{\beta }^4}}{(N_{\beta }-1){(1+2N_{\beta })}^2}}\right].$

So a necessary condition for $S_1\left(t\right)$ having extreme points is $N_{\beta }<1+d_1$ with $d_1={\displaystyle \frac{1}{2}}(\sqrt{3}-1)$∼0.366. Thus, $S_1\left(t\right)$ is non-decreasing only when $N_{\beta }$ is greater than the critical value $N_c\left(1\right)=1+d_1$. For $1<N_{\beta }<N_c\left(1\right)$, $S_1\left(t\right)$ has two extreme points located at $t=t_{-}$ and $t=t_{+}$, respectively, satisfying $S\left(t_{-}\right)>S_{eq}>S\left(t_{+}\right)$. As $N_{\beta }$ decreases from $N_c\left(1\right)$ to 1, $t_{+}$ increases to infinity, so that $S_1\left(t\right)$ exhibits only one hump when $N_{\beta }<1$. In Figure 2, $S_1\left(t\right)$, for a set of different values of $N_{\beta }$, is shown for comparison. Note that in contrast with the case where the initial state is Gaussian, here, we have an additional phase featuring two extremes or two turning points. In this phase, $S\left(t\right)$ first increases and then decreases at $t=t_{-}$ until $t=t_{+}$ when it turns up and eventually saturates.

Finally, an interesting observation is that, although the analytical expressions of $S_n\left(t\right)$ for larger n are more complicated, qualitatively, $S_n\left(t\right)$ is still characterized by the three-phase feature that is the same as $S_1\left(t\right)$, i.e., a “hump” phase, a phase of two turning points, and a monotonically increasing phase.

To summarize, for the two initial conditions investigated, thanks to the GCF, we can solve the motion equations explicitly, based on which the relaxation process of the system can be analyzed in detail. See Ref. [36] for more and detailed discussions.

4.2. The Cat Map

The cat map describes the evolution of a periodically kicked particle with a compact phase space. The Hamiltonian of the system is [38]

$H={\displaystyle \frac{1}{2}}{p}^2+{\displaystyle \frac{K}{2}}{q}^2{\delta }_1\left(t\right),$

where, without loss of generality, both the particle mass and the kicking period have been assumed to be unity. The kicks are exerted instantaneously and accordingly; ${\delta }_1\left(t\right)$ represents a sequence of $\delta$ functions spaced uniformly in time. In addition, the periodic boundary conditions of a unit period are imposed on both q and p, making the phase space a unit torus. Integrating the equations of motion over a unit time from just before the kth kick to just before the $(k+1)$th kick, we have the cat map

$\left(\begin{array}{c}{q}_{k+1}\\ p_{k+1}\end{array}\right)=\left(\begin{array}{cc}1-K\hfill & 1\\ -K\hfill & 1\end{array}\right)\left(\begin{array}{c}{q}_k\\ p_k\end{array}\right)\left(\mathrm{mod}1\right).$

Note that K is the only parameter of the map. If $|K-2|>2$, the map is an Anosov diffeomorphism on the torus and the motion it generates is strongly chaotic (in particular, the motion is mixing and ergodic). The Arnold cat map discussed in the following corresponds to $K=-1$ [39].

Based on the quantization scheme outlined in Section 3 (e) and the quantum counterpart of the classical Hamiltonian given above, the unitary time evolution operator corresponding to the classical cat map can be integrated out exactly, which is

$\widehat{M}=e^{-i\pi N{\widehat{p}}^2}{e}^{-i\pi NK{\widehat{q}}^2},$

with the second and the first factor on the r.h.s. being responsible for the kick and the free motion between two successive kicks, respectively. In terms of density operator, the evolution of the quantum system can be expressed as

(13)${\widehat{\rho }}^{k+1}=\widehat{M}{\widehat{\rho }}^k{\widehat{M}}^{†}.$

Obviously, based on this quantum map directly, it is difficult for both analytical and numerical study of quantum evolution, especially when the dimension of the Hilbert space $N=1/\left(2\pi \hslash \right)$ is large.

However, if we take advantage of the group-theoretical formalism, the quantum map can be greatly simplified [40,41,42]. Setting $\widehat{U}(m,n)=e^{i\pi mn/N}{e}^{i2\pi m\widehat{q}}{e}^{i2\pi n\widehat{p}}$, $(m,n)\in {\mathcal{Z}}^2$, the GCF for a given density operator $\widehat{\rho }$ is defined as

$\phi (m,n)=\mathrm{Tr}\left[\widehat{\rho }\widehat{U}(m,n)\right]$

and the inverse transformation is

$\widehat{\rho }={\displaystyle \frac{1}{N}}\sum _{m=m_0}^{m_0+N-1}\sum _{n=n_0}^{n_0+N-1}\phi (m,n){\widehat{U}}^{†}(m,n).$

It suggests that a quantum ensemble can be assigned equivalently by the GCF $\phi (m,n)$ over any $N\times N$ square lattice (note that $m_0$ and $n_0$ are two arbitrarily chosen integers).

By substituting the GCF into Equation (13), the corresponding quantum evolution can be rewritten as

${\phi }^{k+1}(m_{k+1},n_{k+1})={\phi }^k(m_k,n_k)$

with

$\left(\begin{array}{c}{m}_{k+1}\\ n_{k+1}\end{array}\right)=\left(\begin{array}{cc}1\hfill & K\\ -1\hfill & 1-K\end{array}\right)\left(\begin{array}{c}{m}_k\\ n_k\end{array}\right),$

implying that $\phi (m,n)$ is simply permuted on a finite $N\times N$ lattice. As a consequence, even for a large N, the quantum evolution problem can be dealt with conveniently.

For the cat map, it is straightforward to show that both its purity and the linear entropy do not vary in time. In order to appreciate how its classical mixing and ergodic dynamics affects the quantum motion, we therefore turn to its coarse-grained entropy instead:

$S_{ϵ}^k=-ln\left(\mathrm{Tr}\left[{\left({\widehat{\rho }}_{ϵ}^k\right)}^2\right]\right)=-ln\left(\sum _{m=m_0}^{m_0+N-1}\sum _{n=n_0}^{n_0+N-1}\mid {\phi }_{ϵ}^k(m,n){\mid }^2\right).$

If a Gaussian coarse-graining factor with width $ϵ$ is taken, then for $Nϵ\gg 1$, we have [41]

${\phi }_{ϵ}^k(m,n)\approx e^{-{\pi }^2{ϵ}^2(m^2+n^2)}{\phi }^k(m,n).$

From this expression, we can tell that the role played by the coarse-graining factor is a wave filter such that ${\phi }_{ϵ}^k(m,n)$ can be negligible when $(m,n)$ is far from the origin $(0,0)$.

For an initial Gaussian wavepacket of width a that is centered at $(q_c,p_c)$, we have

${\phi }^0(m,n)=e^{i2\pi (m{q}_c+n{p}_c)}{e}^{-a^2{\pi }^2(m^2+n^2)}.$

The system evolves in terms of the GCF, and its coarse-grained entropy as a function of time is shown in Figure 3. We can see that the coarse-grained entropy increases linearly at the initial stage, signaling the underlying mixing dynamics. After this relaxation stage, the coarse-grained entropy stops increasing but shows some fluctuation patterns due to the quantum coherent effect. As $N\to \infty$, the quantum coherent effect tends to decay and the quantum coarse-grained entropy converges to its classical counterpart as expected. For more detailed discussions and interesting related research on the Wigner distribution and the quantum scar and anti-scar effects, see Ref. [42].

Finally, we would like to mention that, for the numerical computations presented in Figure 3, the cost is the same for any N value. It does not increase as N increases, and this is a remarkable advantage of the GCF representation in this case. In contrast, if the simulations were carried out based on the density matrix, the cost would increase as ∼$N^3$ instead, which would be prohibitively expensive with the computing resources available today.

5. Summary

We have briefly reviewed the mathematical framework of quantum statistical mechanics and the quantum probability theory developed by mathematicians in recent years, in an attempt to address the connection between them. It is worth noting that the latter is a broad field encompassing directions developed with various mathematical assumptions. Here, we have mainly focused on the theories based on ${\mathcal{C}}^{\ast }$- and ∗-algebraic representations of the probability space, particularly analyzing their pros and cons from the perspective of physics. We therefore propose to restrict the ∗-algebra to a specific group algebra where the commutation relations of the dynamical variables of the system are integrated. In particular, we have established the quantum characteristic function in a more mathematically rigorous way than was previously studied in Ref. [16]. The properties of this group-algebraic representation have been discussed in detail, suggesting that it is not only more self-consistent conceptually but also more convenient in certain applications.

To illustrate how to use this suggested group-theoretic formalism in analytical and numerical studies and to appreciate its effectiveness, two application examples have been presented. Obviously, searching for more applicable situations is desirable and should be interesting for future studies.

Author Contributions

Conceptualization, Y.G.; writing—original draft preparation, Y.G. and J.W.; funding acquisition, J.W. All authors have read and agreed to the published version of the manuscript.

Conflicts of Interest

The authors declare no conflicts of interest.

Footnotes

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Figures

Figure 1. The entropy as a function of time for an initial Gaussian state with [Forumla omitted. See PDF.]. The four curves, from top to bottom, are for [Forumla omitted. See PDF.] (green), [Forumla omitted. See PDF.] (red), 1 (black dotted), and 0 (blue), respectively. The black dotted line is for the critical case that [Forumla omitted. See PDF.] and the blue line is for the coherent state with [Forumla omitted. See PDF.].

Figure 2. The entropy [Forumla omitted. See PDF.] for the initial state [Forumla omitted. See PDF.] with [Forumla omitted. See PDF.]. The five curves, from top to bottom, are for [Forumla omitted. See PDF.] (green), [Forumla omitted. See PDF.] (black dashed), [Forumla omitted. See PDF.] (red), [Forumla omitted. See PDF.] (black dotted), and [Forumla omitted. See PDF.] (blue), respectively. Note that the three phases [Forumla omitted. See PDF.] may undergo are separated by the two critical curves, the black dashed and the black dotted line, respectively.

Figure 3. The coarse-grained entropy for the quantum cat map. The initial Gaussian distribution is centered at (0, 0) with [Forumla omitted. See PDF.]. The Gaussian coarse-graining parameter is [Forumla omitted. See PDF.]. The curves, from bottom to top, are for [Forumla omitted. See PDF.] (dashed, magenta), [Forumla omitted. See PDF.] (dot-dashed, red), [Forumla omitted. See PDF.] (short-dashed, blue), and [Forumla omitted. See PDF.] (dotted, black), respectively. As a comparison, the gray solid line is for the corresponding coarse-grained entropy for the classical cat map of the same settings.

Appendix A. Algebraic Probability Spaces Associated with Group Algebras

Appendix A.1. GroupAlgebras of Compact Lie Groups

Let G be a compact Lie group and $dg$ a Haar measure on G, normalized by the condition ${\int }_{G}dg=1$. We denote the set of all infinitely differentiable complex-valued functions on G by $\mathcal{D}\left(G\right)$ and its dual space by ${\mathcal{D}}^{\prime }\left(G\right)$. The functional value of the distribution $f\in {\mathcal{D}}^{\prime }\left(G\right)$ on $\phi \in \mathcal{D}\left(G\right)$ is given by$$(f,\phi )={\int }_{G}f\left(g\right)\phi \left(g\right)dg.$$ Define the involution and multiplication of distributions in ${\mathcal{D}}^{\prime }\left(G\right)$ as(A1)$$f^{\ast }\left(g\right)=\overline{f}\left(g^{-1}\right)\mathrm{and}f\circ k\left(g\right)={\int }_{G}f\left(h\right)k\left(h^{-1}g\right)dh,\forall f,k\in {\mathcal{D}}^{\prime }\left(G\right),$$ respectively; then, ${\mathcal{D}}^{\prime }\left(G\right)$ forms a ∗-algebra with the multiplicative unit $\delta \left(g\right)$.

By definition, the function $\phi \in \mathcal{D}\left(G\right)$ is (i) non-negative definite if used for arbitrary sets ${\alpha }_1,\dots ,{\alpha }_m\in \mathbb{C}$ and $g_1,\dots ,g_m\in G$, $m\in \mathbb{N}$,$$\sum _{i,j=1}^m{\alpha }_i\overline{{\alpha }_j}\phi \left(g_i{g}_j^{-1}\right)\ge 0$$ and (ii) normalized if $\phi \left(e\right)=1$, where e is the identity in G. A normalized non-negative definite function $\phi \in \mathcal{D}\left(G\right)$ is termed a group-theoretical characteristic function (GCF) in the following.

We denote the space of all these characteristic functions by $\mathcal{S}\left(G\right)$, which is a convex subspace of $\mathcal{D}\left(G\right)$.

Let $\widehat{U}\left(g\right)$ be an irreducible unitary representation of G in a Hilbert space $\mathcal{H}$ of dimension N and $\left\{\right|{\mu }_j⟩$, $j=1,2,\dots ,N\}$ be a normalized orthogonal basis in $\mathcal{H}$, such that $⟨{\mu }_i\left|\widehat{U}\left(g\right)\right|{\mu }_j⟩\in \mathcal{D}\left(G\right)$, $\forall 1\le i,j\le N$. Then, there exists a ∗-representation of ${\mathcal{D}}^{\prime }\left(G\right)$ in $\mathcal{H}$ (p. 511 in Ref. [43]),$${\mathcal{D}}^{\prime }\left(G\right)\ni f\to {\widehat{\pi }}_f={\int }_{G}f\left(g\right)\widehat{U}\left(g\right)dg.$$ It is straightforward to prove that ${\widehat{\pi }}_{f^{\ast }}={\widehat{\pi }}_f^{†}$ and ${\widehat{\pi }}_{f\circ k}={\widehat{\pi }}_f{\widehat{\pi }}_k$.

On the other hand, for any density matrix $\widehat{\rho }$ in $\mathcal{H}$, we can define a GCF on G as(A2)$${\phi }_{\rho }\left(g\right)=\mathrm{Tr}\left[\widehat{\rho }\widehat{U}\left(g\right)\right]\in \mathcal{S}\left(G\right).$$ Using the orthogonality relations for matrix elements of an irreducible unitary representation of G (p. 170 in Ref. [8]), i.e.,$${\int }_G⟨\mu |\widehat{U}\left(g\right)|\nu ⟩⟨{\mu }^{\prime }|\widehat{U^{†}}\left(g\right)|{\nu }^{\prime }⟩dg={\delta }_{\mu {\nu }^{\prime }}{\delta }_{\nu {\mu }^{\prime }}/N,$$ we can derive from Equation (A2) that$$\widehat{\rho }=N{\int }_G{\phi }_{\rho }\left(g\right){\widehat{U}}^{†}\left(g\right)dg,$$ as well as the Parseval equality$${\int }_G{\left|{\phi }_{\rho }\left(g\right)\right|}^{2}dg={\displaystyle \frac{1}{N}}\mathrm{Tr}\left[{\widehat{\rho }}^2\right].$$ Furthermore, the expectation value of random variable f with respect to ${\phi }_{\rho }$ is$$(f,{\phi }_{\rho })={\int }_{G}f\left(g\right){\phi }_{\rho }\left(g\right)dg=\mathrm{Tr}\left[\widehat{\rho }{\int }_{G}f\left(g\right)\widehat{U}\left(g\right)dg\right]=\mathrm{Tr}\left[\widehat{\rho }{\widehat{\pi }}_f\right].$$ Since $(f\circ f^{\ast },{\phi }_{\rho })=\mathrm{Tr}\left[\widehat{\rho }{\widehat{\pi }}_f{\widehat{\pi }}_{f^{\ast }}\right]\ge 0$, $\forall f\in {\mathcal{D}}^{\prime }\left(G\right)$ and $\mathcal{D}\left(G\right)={\mathcal{D}}^{\prime \prime }\left(G\right)$, the set of all ${\phi }_{\rho }$ are states of the ∗-algebra ${\mathcal{D}}^{\prime }\left(G\right)$ and $({\mathcal{D}}^{\prime }\left(G\right),{\phi }_{\rho })$ is an algebraic space associated with group algebra ${\mathcal{D}}^{\prime }\left(G\right)$. Thus, after introducing an irreducible unitary representation $G\ni g\to \widehat{U}\left(g\right)$, all real elements $f=f^{\ast }\in {\mathcal{D}}^{\prime }\left(G\right)$ are mapped to self-adjoint (hermitian) matrices ${\widehat{\pi }}_f$ and any density matrix $\widehat{\rho }$ in $\mathcal{H}$ can be represented by a GCF ${\phi }_{\rho }\left(g\right)$. We denote the set of ${\phi }_{\rho }\left(g\right)$ for all density matrices in $\mathcal{H}$ by ${\mathcal{S}}_U\left(G\right)$, which forms a convex subspace of $\mathcal{S}\left(G\right)$.

Now, we generalize the mapping (A2) between a GCF $\phi \left(g\right)\in \mathcal{S}\left(G\right)$ and the corresponding density matrix to the case where $\widehat{U}\left(g\right)$ is a reducible unitary representation of the compact group G. Let $\widehat{G}=\left\{U_{\lambda \in \Lambda }\right\}$ be an arbitrary set of irreducible non-equivalent unitary representations of G and $\Lambda =\{{\lambda }_1,{\lambda }_2,\dots ,{\lambda }_n;i\ne j\Rightarrow {\lambda }_i\ne {\lambda }_j\}$. Considering a reducible unitary representation $\widehat{U}\left(g\right)={⨁}_{j=1}^n{\widehat{U}}_{{\lambda }_j}\left(g\right)$, where ${\widehat{U}}_{{\lambda }_j}$ is an irreducible unitary representation of G in a Hilbert space ${\mathcal{H}}_j$ of dimension $d_j$, and letting ${\widehat{\rho }}_j$ be a density matrix in ${\mathcal{H}}_j$ and ${\phi }_j\left(g\right)={\mathrm{Tr}}_{{\mathcal{H}}_j}\left[{\widehat{\rho }}_j{\widehat{U}}_{{\lambda }_j}\left(g\right)\right]$, we then have$${\widehat{\rho }}_j=d_j{\int }_G{\phi }_j\left(g\right){\widehat{U}}_{{\lambda }_j}^{†}\left(g\right)dg.$$

Note that, between the special subset of density matrices in the Hilbert space $\mathcal{H}={⨁}_{j=1}^n{\mathcal{H}}_j$, i.e.,(A3)$$\widehat{\rho }=\underset{j=1}{\overset{n}{⨁}}{c}_j{\widehat{\rho }}_j,c_j>0\mathrm{and}\sum _{j=1}^n{c}_j=1,$$ and the set of states $\phi \left(g\right)\in {\mathcal{S}}_U\left(G\right)$ for the reducible representations $\widehat{U}\left(g\right)={⨁}_{j=1}^n{\widehat{U}}_{{\lambda }_j}\left(g\right)$, there exists a one-to-one correspondence given by the mapping(A4)$${\phi }_{\rho }\left(g\right)={\mathrm{Tr}}_{\mathcal{H}}\left[\widehat{\rho }\widehat{U}\left(g\right)\right]=\sum _{j=1}^n{c}_j{\phi }_j\left(g\right).$$ In fact, for a given density matrix $\widehat{\rho }$ in the Hilbert space $\mathcal{H}$, the GCF given by Equation (A4) eliminates all the interference terms of the matrix $\widehat{\rho }$ between different coherent subspaces $\left\{{\mathcal{H}}_i\right\}$ in a straightforward way. Thus, if the states of a quantum system can be described with the group-algebraic formalism of quantum probability, they should be represented by the density matrices satisfying Equation (A3) as well. In other words, for such systems, the superposition of the pure states in different coherent subspaces ${\mathcal{H}}_i$ and ${\mathcal{H}}_j$ is not allowed. In the physical literature, such limitations on the superposition principle are known as the “super-selection rules” for a quantum system.

Example: $G=SU\left(2\right)$

The generators of $SU\left(2\right)$ form a three-dimensional real Lie algebra with the Lie brackets$$[X_i,X_j]=-{ϵ}_{ijk}{X}_k.$$ In the neighbourhood of the identity, we have $G\ni g=e^{\overrightarrow{\alpha }·\overrightarrow{X}}$, where $\{\overrightarrow{\alpha }\in {\mathbb{R}}^3,\alpha =|\overrightarrow{\alpha }|\le 2\pi \}$ is a canonical coordinate in G. Since $|\overrightarrow{\alpha }|=2\pi$ represents only one element, the manifold of G is isomorphic to $S^3$. In the following, we will take a non-canonical coordinate $g\to (x_0,\overrightarrow{x})$, where $x_0=cos{\displaystyle \frac{\alpha }{2}}$, $\overrightarrow{x}={\displaystyle \frac{\overrightarrow{\alpha }}{\alpha }}sin{\displaystyle \frac{\alpha }{2}}$, and $x_0^2+{\left|\overrightarrow{x}\right|}^2=1$. In this coordinate, the normalized Haar measure of G is given by$$dg={\displaystyle \frac{d\overrightarrow{x}}{2{\pi }^2\left|x_0\right|}}{|}_{x_0=\pm \sqrt{1-|\overrightarrow{x}{|}^2}}.$$

Let $g\to {\widehat{U}}_J\left(g\right)$ with $J=0,{\displaystyle \frac{1}{2}},1,{\displaystyle \frac{3}{2}},\dots$ be the set of irreducible unitary representations of G in the complex Euclidean space ${\mathbb{C}}^{2J+1}$.

(1) Two-level systems

For $J={\displaystyle \frac{1}{2}}$, we have ${\widehat{U}}_{{\displaystyle \frac{1}{2}}}\left(g\right)=x_0+\overrightarrow{x}·\widehat{\overrightarrow{\sigma }}$, where $\widehat{\overrightarrow{\sigma }}$ are Pauli matrices. The general form of the density matrix in ${\mathbb{C}}^2$ can be written as$$\widehat{\rho }={\displaystyle \frac{1}{2}}(1+\overrightarrow{P}·\widehat{\overrightarrow{\sigma }}),$$ where $\overrightarrow{P}$ is the polarization vector of the system, satisfying $|\overrightarrow{P}|\le 1$. The state represented by $\widehat{\rho }$ is a pure state if and only if $|\overrightarrow{P}|=1$. It is straightforward to verify that all the mixed state ($|\overrightarrow{P}|<1$) can be decomposed into a sum of two orthogonal pure states and this decomposition is unique except for the mixed state with $|\overrightarrow{P}|=0$. We note that, for a two-level system, the GCF of the density matrix $\widehat{\rho }$ is simply given by a monomial of the variables $\{x_{\mu },\mu =0,1,2,3\}$, i.e.,$$\phi \left(g\right)=\mathrm{Tr}\left[\widehat{\rho }{\widehat{U}}_{{\displaystyle \frac{1}{2}}}\left(g\right)\right]=x_0+i\overrightarrow{P}·\widehat{\overrightarrow{x}}.$$

(2) Three-level systems

For a three-level system, pure states are represented by vectors in the representation space ${\mathbb{C}}^3$ of ${\widehat{U}}_1\left(g\right)$. Let $a=x_0+i{x}_3$ and $b=x_2+i{x}_1$, the unitary matrix ${\widehat{U}}_1\left(g\right)$ can be written as(A5)$${\widehat{U}}_1\left(g\right)=\left(\begin{array}{ccc}{a}^2& \sqrt{2}ab& b^2\\ -\sqrt{2}a\overline{b}& {\left|a\right|}^2-{\left|b\right|}^2& \sqrt{2}\overline{a}b\\ {\overline{b}}^2& -\sqrt{2}\overline{a}\overline{b}& {\overline{a}}^2\end{array}\right).$$ The mixed states of the system are represented by density matrices in ${\mathbb{C}}^3$. Here, we are especially interested in those density matrices that are diagonalizable with SU(2) transformations given by Equation (A5). Introducing variables $u=p_0+i{p}_3$ and $v=p_2+i{p}_1$ satisfying ${\left|u\right|}^2+{\left|v\right|}^2=p_0^2+p_1^2+p_2^2+p_3^2=1$ and defining three orthonormal pure states of the system by density matrices$${\widehat{\rho }}_1=\left(\begin{array}{ccc}{\left|u\right|}^4& \sqrt{2}{\overline{u}}^{2}uv& {\overline{u}}^2{v}^2\\ \sqrt{2}{\left|u\right|}^{2}u\overline{v}& {2\left|u\right|}^2{\left|v\right|}^2& \sqrt{2}\overline{u}v{\left|v\right|}^2\\ u^2{\overline{v}}^2& \sqrt{2}u\overline{v}{\left|v\right|}^2& {\left|v\right|}^4\end{array}\right),$$$${\widehat{\rho }}_2=\left(\begin{array}{ccc}{2\left|u\right|}^2{\left|v\right|}^2& -\sqrt{2}{\left(\right|u|}^2-{\left|v\right|}^2)\overline{u}v& -2{\overline{u}}^2{v}^2\\ -\sqrt{2}{\left(\right|u|}^2-{\left|v\right|}^2)u\overline{v}& {\left(\right|u|}^2-{\left|v\right|}^2{)}^2& \sqrt{2}\overline{u}{v\left(\right|u|}^2-{\left|v\right|}^2)\\ -2u^2{\overline{v}}^2& \sqrt{2}u\overline{v}{\left(\right|u|}^2-{\left|v\right|}^2)& {2\left|u\right|}^2{\left|v\right|}^2\end{array}\right),$$$${\widehat{\rho }}_3=\left(\begin{array}{ccc}{\left|v\right|}^4& -\sqrt{2}{\left|v\right|}^2\overline{u}v& {\overline{u}}^2{v}^2\\ -\sqrt{2}{\left|v\right|}^{2}u\overline{v}& {2\left|u\right|}^2{\left|v\right|}^2& -\sqrt{2}uv{\left|u\right|}^2\\ u^2{\overline{v}}^2& -\sqrt{2}u\overline{v}{\left|u\right|}^2& {\left|u\right|}^4\end{array}\right)$$ with $\mathrm{Tr}\left[{\widehat{\rho }}_i{\widehat{\rho }}_j\right]={\delta }_{ij}$, $i,j=1,2,3$, all the density matrices that are diagonalizable with SU(2) can be expressed as $\widehat{\rho }=c_1{\widehat{\rho }}_1+c_2{\widehat{\rho }}_2+c_3{\widehat{\rho }}_3$, where $c_1+c_2+c_3=1$ and $c_j\ge 0$.

The GCF of this orthonormal set of density matrices are quadratic forms of the variables $\{x_{\mu },\mu =0,1,2,3\}$:$${\phi }_1\left(g\right)=\mathrm{Tr}\left[{\widehat{\rho }}_1{\widehat{U}}_1\left(g\right)\right]=(x_0+i(\overline{u}v+u\overline{v})x_1+(u\overline{v}-\overline{u}v)x_2+{i\left(\right|u|}^2-{\left|v\right|}^2{\left)x_3\right)}^2,$$$${\phi }_2\left(g\right)=\mathrm{Tr}\left[{\widehat{\rho }}_2{\widehat{U}}_1\left(g\right)\right]=x_0^2+a_{11}{x}_1^2+a_{22}{x}_2^2+a_{33}{x}_3^2+a_{12}{x}_1{x}_2+a_{13}{x}_1{x}_3+a_{23}{x}_2{x}_3,$$$${\phi }_3\left(g\right)=\mathrm{Tr}\left[{\widehat{\rho }}_3{\widehat{U}}_1\left(g\right)\right]=(x_0-i(\overline{u}v+u\overline{v})x_1-(u\overline{v}-\overline{u}v)x_2-{i\left(\right|u|}^2-{\left|v\right|}^2{\left)x_3\right)}^2,$$ where$$a_{11}={2\left|u\right|}^2{\left|v\right|}^2-u^2({\overline{u}}^2-2{\overline{v}}^2)+v^2(2{\overline{u}}^2-{\overline{v}}^2),$$$$a_{12}=4i({\overline{u}}^2{v}^2-u^2{\overline{v}}^2),$$$$a_{13}={4\left(\right|u|}^2-{\left|v\right|}^2)(u\overline{v}+\overline{u}v),$$$$a_{22}={2\left|u\right|}^2{\left|v\right|}^2-u^2({\overline{u}}^2+2{\overline{v}}^2)-v^2(2{\overline{u}}^2+{\overline{v}}^2),$$$$a_{23}=-{4i\left(\right|u|}^2-{\left|v\right|}^2)(u\overline{v}-\overline{u}v),$$$$a_{33}={\left|u\right|}^4+{\left|v\right|}^4-6{\left|u\right|}^2{v}^2.$$ They satisfy the following orthogonal relations:$${\int }_G{\phi }_i\left(g\right){\overline{\phi }}_j\left(g\right)dg={\displaystyle \frac{1}{3}}{\delta }_{ij},i,j=1,2,3.$$

Appendix A.2. Group Algebras of Locally Compact Lie Groups

For a unimodular locally compact Lie group G associated with a fixed bi-invariant Haar measure $dg$, let $\mathcal{E}\left(G\right)$ be the space of all infinitely differentiable functions on G with the topology of uniform convergence on compact sets. The dual space ${\mathcal{E}}^{\prime }\left(G\right)$ of $\mathcal{E}\left(G\right)$ consists of all distributions with compact supports on G so that both $\mathcal{E}\left(G\right)$ and ${\mathcal{E}}^{\prime }\left(G\right)$ are countably normed topological spaces. Just the same as the case for compact Lie groups, we define the involution and multiplication of distributions in ${\mathcal{E}}^{\prime }\left(G\right)$ by Equation (A1), so that ${\mathcal{E}}^{\prime }\left(G\right)$ forms a ∗-algebra. We denote the set of all normalized non-negative definite functions (GCF) $\phi \in \mathcal{E}\left(G\right)$ by $\mathcal{S}\left(G\right)$. Then $({\mathcal{E}}^{\prime }\left(G\right),\phi )$ is an algebraic probability space when $\phi \in \mathcal{S}\left(G\right)$ is a state of ${\mathcal{E}}^{\prime }\left(G\right)$.

Let $\widehat{U}\left(g\right)$ be an irreducible unitary representation of G in a Hilbert space $\mathcal{H}$ and $|\mu ⟩$ a non-zero vector in $\mathcal{H}$ such that $⟨\mu |\widehat{U}\left(g\right)|\mu ⟩\in \mathcal{E}\left(G\right)\bigcap L^2(G,dg)$, the vector set $\left\{\widehat{U}\left(g\right)\right|\mu ⟩,g\in G\}$ forms an invariant dense subspace in $\mathcal{H}$, denoted by ${\mathcal{L}}_U$, so that all vectors in ${\mathcal{L}}_U$ are both infinitely differentiable and square integrable with respect to G. Let $\left\{\right|{\mu }_j⟩,j=1,2,\dots \}$ be a normalized orthogonal basis in ${\mathcal{L}}_U$. Using this basis, for every density matrix $\widehat{\rho }$ in $\mathcal{H}$, we can define a GCF on G as(A6)$${\phi }_{\rho }\left(g\right)=\mathrm{Tr}\left[\widehat{\rho }\widehat{U}\left(g\right)\right]\in \mathcal{S}\left(G\right).$$ Note that, for square integrable irreducible unitary representations of a locally compact group, the orthogonality relations for matrix elements can be written as (p. 138 in Ref. [24])$${\int }_G⟨\mu |\widehat{U}\left(g\right)|\nu ⟩⟨{\mu }^{\prime }|\widehat{U^{†}}\left(g\right)|{\nu }^{\prime }⟩dg={\delta }_{\mu {\nu }^{\prime }}{\delta }_{\nu {\mu }^{\prime }}/d_U,$$ where $d_U$ is a positive constant termed as the formal dimension of the representation $\widehat{U}\left(g\right)$. In the following, we will fix the Haar measure $dg$ such that $d_U=1$. Thus, from Equation (A6), we have$$\widehat{\rho }={\int }_G{\phi }_{\rho }\left(g\right)\widehat{U^{†}}\left(g\right)dg$$ and the Parseval equality now has the form$${\int }_G{\left|{\phi }_{\rho }\left(g\right)\right|}^{2}dg=\mathrm{Tr}\left[{\widehat{\rho }}^2\right].$$ On the other hand, based on this irreducible unitary representation in the Hilbert space $\mathcal{H}$, we can define a ∗-representation of the ∗-algebra ${\mathcal{E}}^{\prime }\left(G\right)$ in the pre-Hilbert space ${\mathcal{L}}_U$ as$${\mathcal{E}}^{\prime }\left(G\right)\ni f\to {\widehat{\pi }}_f={\int }_{G}f\left(g\right)\widehat{U}\left(g\right)dg.$$ The expectation value of the random variable f with respect to ${\phi }_{\rho }$ is given by$$(f,{\phi }_{\rho })={\int }_{G}f\left(g\right){\phi }_{\rho }\left(g\right)dg=\mathrm{Tr}\left[\widehat{\rho }{\widehat{\pi }}_f\right]$$ which ensures $(f\circ f^{\ast },{\phi }_{\rho })=\mathrm{Tr}\left[\widehat{\rho }{\widehat{\pi }}_f{\widehat{\pi }}_{f^{\ast }}\right]\ge 0$, $\forall f\in {\mathcal{E}}^{\prime }\left(G\right)$. That is, ${\phi }_{\rho }$ is a state of ${\mathcal{E}}^{\prime }\left(G\right)$.

For a real element $f=f^{\ast }\in {\mathcal{E}}^{\prime }\left(G\right)$, ${\widehat{\pi }}_f$ is a well-defined symmetric operator on ${\mathcal{L}}_U$. A symmetric operator on a pre-Hilbert space is said to be essentially self-adjoint if it admits a unique self-adjoint extension on the Hilbert space. It is widely believed among physicists that, in order to ensure a proper interpretation of measurements, these real random variables should be represented by at least essentially self-adjoint operators. However, in our opinion, to provide a mathematical framework of quantum probability to quantum statistical mechanics, it is satisfying to ensure only a certain set of dynamical variables of physical importance for describing the system to be represented by essentially self-adjoint operators. In fact, the dynamical evolution of a quantum system is described by its state specified by the density matrix. In spite of some dynamical variables that are relevant to the master equation of the system, the solubility of the equation does not depend on whether they are represented by essentially self-adjoint operators or not.

For further elucidation, let us consider a subalgebra $R(G,K)$ of ${\mathcal{E}}^{\prime }\left(G\right)$ that consists of all distributions with compact supports contained in a subgroup $K\subset G$, in particular the subalgebra $R(G,\{e\left\}\right)$, which is isomorphic to the enveloped algebra $\mathcal{U}\left(\mathfrak{g}\right)$ of the Lie algebra $\mathfrak{g}$ of group G and its symmetric elements are usually supposed to be observables. Let X be an arbitrary element of the Lie algebra $\mathfrak{g}$ and $p\left(X\right)\in \mathcal{U}\left(\mathfrak{g}\right)$ be any real polynomial function of X; it is known that the unitary representation of G admits $p\left(iX\right)$ to be represented by an essentially self-adjoint operator on ${\mathcal{L}}_U$. For the expressions of observables as distributions $f\left(g\right)\in {\mathcal{E}}^{\prime }\left(G\right)$, noting that in the neighborhood of the identity of G, the unitary representation can be expressed as $\widehat{U}\left(g\right)=e^{{\sum }_j{x}_j\widehat{X_j}}$, the corresponding random variable for the real element $iX\in \mathcal{U}\left(\mathfrak{g}\right)$ in the group-algebraic probability formalism is represented by the distribution $f\left(g\right)=-i{\displaystyle \frac{\partial }{\partial x}}\delta \left(g\right)\in R(G,\left\{e\right\})$. However, it should be emphasized that, in a unitary representation of a Lie group, the symmetric elements of the enveloped algebra $\mathcal{U}\left(\mathfrak{g}\right)$ are not necessarily represented by essentially self-adjoint operators [8].

References

1. Laplace, P.S. Théorie Analytique des Probabilités; Mme ve Courcier: Paris, France, 1812.

2. Gibbs, J.W. Elementary Principles in Statistical Mechanics; Yale University Press: New Haven, CT, USA, 1902.

3. von Mieses, R. Grundlagen der Wahrscheinlichkeitrechnung. Math. Z.; 1919; 5, pp. 52-99. [DOI: https://dx.doi.org/10.1007/BF01203155]

4. Kolmogorov, A.N. Foundations of the Probability Theory; Chelsea: New York, NY, USA, 1956.

5. Halmos, P.R. Measure Theory; Springer: New York, NY, USA, 1950.

6. Born, M. Zur Quantenmechanik der Stoßvorgänge. Z. Physik; 1926; 37, pp. 863-867. [DOI: https://dx.doi.org/10.1007/BF01397477]

7. von Neumann, J. Mathematical Foundations of Quantum Mechanics; Princeton University Press: Princeton, NJ, USA, 1955.

8. Barut, A.O.; Raczka, R. Theory of Group Representations and Applications; Polish Scientific Publishers: Warsaw, Poland, 1980.

9. Accardi, L. Some Trends and Problems in Quantum Probability. Lecture Notes in Mathematics; Springer: Berlin/Heidelberg, Germany, 1984; Volume 1055.

10. Accardi, L. Quantum Probability: New Perspectives for the Laws of Chance. Milan J. Math.; 2010; 78, 481. [DOI: https://dx.doi.org/10.1007/s00032-010-0134-3]

11. Schlosshauer, M. Decoherence and the Quantum-to-Classical Transition; Springer: Berlin/Heidelberg, Germany, 2007.

12. Segal, I.E. Postulates for general quantum mechanics. Ann. Math.; 1947; 48, 930. [DOI: https://dx.doi.org/10.2307/1969387]

13. Bratteli, O.; Robibson, D.W. Operator Algebras and Quantum Statistical Mechanics; Springer: Berlin/Heidelberg, Germany, 1979; Volumes I–II.

14. Haag, R. Local Quantum Physics; Springer: Berlin/Heidelberg, Germany, 1992.

15. Meyer, P.A. Quantum Probability for Probabilists; Springer: Berlin/Heidelberg, Germany, 1995.

16. Gu, Y. Group-theoretical formalism of quantum mechanics based on quantum generalization of characteristic functions. Phys. Rev. A; 1985; 32, 1310. [DOI: https://dx.doi.org/10.1103/PhysRevA.32.1310] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/9896209]

17. Moretti, V. Fundamental Mathematical Structures of Quantum Theory; Springer: Berlin/Heidelberg, Germany, 2019.

18. Berg, C.; Christensen, J.F.R.; Ressel, P. Harmonic Analysis on Semigroups; Springer: Berlin/Heidelberg, Germany, 1984.

19. Parthasarathy, K.R. An Introduction to Quantum Stochastic Calculus; Springer: Basel, Switzerland, 1992.

20. Pedersen, G.K. C^∗-Algebras and Their Automorphism Groups; Academic Press: Cambridge, MA, USA, 1979.

21. Gu, Y. Group-theoretical formalism of quantum mechanics and classical-quantum correspondence. Sci. China A; 1992; 35, 200.

22. Korbiez, J.K.; Lewenstein, M. Group-theoretical approach to entanglement. Phys. Rev. A; 2006; 74, 022318. [DOI: https://dx.doi.org/10.1103/PhysRevA.74.022318]

23. Marsden, J.; Weinstein, A. Coadjoint orbits, vortices, and Clebsch variables for incompressible fluids. Phys. D Nonlinear Phenom.; 1983; 7, pp. 305-323. [DOI: https://dx.doi.org/10.1016/0167-2789(83)90134-3]

24. Kirillov, A.A. Elements of The Theory of Representations; Springer: Berlin/Heidelberg, Germany, 1976.

25. Dufour, J.P.; Nguyen, T.Z. Poisson Structures and Their Normal Forms; Birkhauser Verlag: Basel, Switzerland, 2005.

26. Wolf, K.B. The Heisenberg-Weyl Ring in Quantum Mechanics. Group Theory and Its Applications; Loebl, E.M. Academic Press: New York, NY, USA, 1975; Volume III, 189.

27. Balaz, N.L.; Voros, A. The quantized Baker’s transformation. Ann. Phys.; 1989; 190, pp. 1-31. [DOI: https://dx.doi.org/10.1016/0003-4916(89)90259-5]

28. Gu, Y. Steady states of quantum Brownian motion and decomposition of quantum states into ensembles of Gaussian packets having a uniform position variance. Phys. Scr.; 2019; 94, 115205. [DOI: https://dx.doi.org/10.1088/1402-4896/ab2cc1]

29. Breuer, H.P.; Petruccione, F. The Theory of Open Quantum Systems; Oxford University Press: Oxford, UK, 2002.

30. Weiss, D. Quantum Dissipative Systems; World Scientific: Singapore, 2000.

31. Louisell, W.H. Quantum Statistical Properties of Radiation; John Wiley & Sons: New York, NY, USA, 1973.

32. Haake, F. Statistical Treatment of Open Systems by Generalized Master Equations; Springer Tracts in Modern Physics Springer: Berlin/Heidelberg, Germany, 1973; Volume 66.

33. Walls, D.F.; Milburn, G.J. Quantum Optics; Springer: Berlin/Heidelberg, Germany, 2008.

34. Miguel, O. Quantum Optics; Springer: Berlin/Heidelberg, Germany, 2007.

35. Prigogine, I. Time, Irreversibility and Structure. The Physicist’s Conception of Nature; Mehra, J. Reidel: Dordrecht, The Netherlands, 1973; 579.

36. Gu, Y.; Wang, J. Relaxation of A Thermally Bathed Harmonic Oscillator: A Study Based on the Group-theoretical Formalism. arXiv; 2024; arXiv: 2012.16981

37. Gradshteyn, L.S.; Ryzhik, L.M. Table of Integrals, Series and Products; Academic Press: New York, NY, USA, 1980; 837.

38. Ford, J.; Mantica, G.; Ristow, G.H. The Arnol’d cat: Failure of the correspondence principle. Phys. D Nonlinear Phenom.; 1991; 50, pp. 493-520. [DOI: https://dx.doi.org/10.1016/0167-2789(91)90012-X]

39. Arnold, V.I.; Avez, A. Ergodic Problems of Classical Mechanics; Benjamin: New York, NY, USA, 1968.

40. Gu, Y. Evidences of classical and quantum chaos in the time evolution of nonequilibrium ensembles. Phys. Lett. A; 1990; 149, pp. 95-100. [DOI: https://dx.doi.org/10.1016/0375-9601(90)90532-S]

41. Gu, Y.; Wang, J. Time evolution of coarse-grained entropy in classical and quantum motions of strongly chaotic systems. Phys. Lett. A; 1997; 229, 208. [DOI: https://dx.doi.org/10.1016/S0375-9601(97)00194-1]

42. Wang, J.; Lai, C.-H.; Gu, Y. Ergodicity and scars of the quantum cat map in the semiclassical regime. Phys. Rev. E; 2001; 63, 056208. [DOI: https://dx.doi.org/10.1103/PhysRevE.63.056208] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/11414991]

43. Naimark, M.A.; Stern, A.I. Theory of Group Representations; Springer: Berlin/Heidelberg, Germany, 1980.