Next Article in Journal
Instability of Meissner Differential Equation and Its Relation with Photon Excitations and Entanglement in a System of Coupled Quantum Oscillators
Previous Article in Journal
Socio-Economic Sciences: Beyond Quantum Math-like Formalisms
Previous Article in Special Issue
How Does the Planck Scale Affect Qubits?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

From Geometry to Coherent Dissipative Dynamics in Quantum Mechanics

by
Hans Cruz-Prado
1,
Alessandro Bravetti
2 and
Angel Garcia-Chung
3,4,*
1
Instituto de Ciencias Nucleares, Universidad Nacional Autónoma de México, A. P. 70543, Ciudad de México 04510, Mexico
2
Instituto de Investigaciones en Matemáticas Aplicadas y en Sistemas, Universidad Nacional Autónoma de México, A. P. 70543, Ciudad de México 04510, Mexico
3
Departamento de Física, Universidad Autónoma Metropolitana—Iztapalapa, San Rafael Atlixco 186, Ciudad de México 09340, Mexico
4
Tecnológico de Monterrey, Escuela de Ingeniería y Ciencias, Carr. al Lago de Guadalupe Km. 3.5, Ciudad López Mateos 52926, Mexico
*
Author to whom correspondence should be addressed.
Quantum Rep. 2021, 3(4), 664-683; https://doi.org/10.3390/quantum3040042
Submission received: 30 July 2021 / Revised: 25 September 2021 / Accepted: 28 September 2021 / Published: 12 October 2021
(This article belongs to the Special Issue Relevance of Information Geometry in Quantum Information Science)

Abstract

:
Starting from the geometric description of quantum systems, we propose a novel approach to time-independent dissipative quantum processes according to which energy is dissipated but the coherence of the states is preserved. Our proposal consists of extending the standard symplectic picture of quantum mechanics to a contact manifold and then obtaining dissipation by using appropriate contact Hamiltonian dynamics. We work out the case of finite-level systems for which it is shown, by means of the corresponding contact master equation, that the resulting dynamics constitute a viable alternative candidate for the description of this subclass of dissipative quantum systems. As a concrete application, motivated by recent experimental observations, we describe quantum decays in a 2-level system as coherent and continuous processes.

1. Motivation and Previous Works

Dissipative quantum phenomena have been the subject of intense investigation since the early days of quantum mechanics [1]. The most widely adopted description of dissipative quantum processes is given by the Gorini–Kossakowski–Sudarshan–Lindblad (GKLS) equation [2,3,4]. This is because it has been proven that the GKLS equation describes the most general form of a non-unitary, linear, completely positive and trace-preserving dynamics; thus, one may ensure that the probabilistic character of quantum mechanics is preserved at all times. However, in the GKLS equation, in order to obtain a completely positive evolution, it is necessary to add the so-called jump term, which changes the rank of the state. As a consequence, the GKLS dynamics fundamentally couples energy dissipation and rank disipation.
On the other hand, it is possible to model energy dissipation without change in the rank of ρ by means of a time-dependent Hamiltonian, where the time dependence describes, in an effective manner, the interaction of the system with the environment. In this case, the evolution is unitary; therefore, the rank of the state is preserved. For instance, recent impressive improvements in experimental settings have allowed the observation that the decays induced by a Rabi oscillator in a 2-level system are coherent and continuous processes [5] (see also [6] for theoretical analysis).
In this work, we propose an alternative method to describe a class of time-independent dissipative quantum phenomena by non-unitary evolution-preserving statistical mixing. This is based on the symplectic formulation of quantum mechanics and due to the analogy with the geometric description of classical dissipative systems: In classical mechanics, one can describe a wide class of dissipative systems by referring to the contactification of the symplectic phase space and then using contact Hamiltonian systems to define the dynamics. It has been shown that this approach, when applicable, provides several positive features, such as relying on canonical variables and producing a generalization of canonical transformations [7], enabling an extension of both Liouville and Noether’s theorems to the dissipative case [8,9,10,11,12,13], and providing a description in terms of variational principles [14,15,16,17,18,19,20] together with a natural route to field theories with dissipation [21].
Analogously, in the geometric description of n-level quantum systems, the pure states of the corresponding Hilbert space H are points on the complex projective space CP ( H 0 ) , which is a symplectic manifold, and the Schrödinger dynamics on H is projected onto Hamiltonian dynamics on such a manifold [22,23,24,25]. This is because pure quantum states in H are actually rays, and the Schrödinger equation is invariant under dilations and multiplications by a global phase factor; that is, in the standard Schrödinger picture, one deals with redundant information in order to obtain a linear description. Interestingly, contact manifolds also appear naturally in this picture. In fact, by quotienting only over dilations, one obtains the manifold of normalized vectors in H identified as the ( 2 n 1 ) -dimensional unit sphere S 2 n 1 ( H ) with the standard contact structure. An approach to describe dissipation in the quantum case by means of contact Hamiltonian flows on S 2 n 1 ( H ) has been put forward already in [26].
Contrary to [26], we do not work on S 2 n 1 ( H ) , but we proceed by analogy with the classical case, i.e., we perform a contactification of CP ( H 0 ) and investigate the dissipative dynamics by means of contact Hamiltonian systems defined on the extended space CP ( H 0 ) × R . We take a different route that is more directly connected with the procedure of contactification used in classical mechanics.
The structure of this work is as follows: after a brief review of the standard geometric description of n-level quantum systems in Section 2, we introduce our approach in Section 3, and we show that, by choosing the contact Hamiltonian appropriately, the dynamics on CP ( H 0 ) × R is projectable onto the proper dynamics on CP ( H 0 ) , thus preserving the purity of states while at the same time dissipating the expected value of the energy of the reference system. In this manner, we obtain a dissipative dynamics on the manifold of physical quantum states. Then, in Section 4, we consider the corresponding dynamics for the density operators, which we call the contact master equation, and we prove that the resulting map is both positive and trace preserving after a brief comparison with other approaches, thus agreeing with all the prescriptions of standard quantum mechanics. As a consequence, the proposed contact dynamics is a viable candidate for describing dissipative phenomena in quantum systems when the purity of the states is preserved. As an illustration of our formalism, we consider in Section 5 the important case of radiative decay in qubit systems, finding that our description can effectively model quantum decays (or excitations) as coherent and continuous processes. Finally, in Section 6, we summarize our results and highlight future directions.

2. Geometry and Dynamics of Conservative n-Level Quantum Systems

Let us start our study by recalling some aspects of the geometric description of finite-level quantum systems. For complete reviews, we refer to [22,23,25,27,28], while further results can be found in [24,26,29,30].

2.1. Kinematics: From H 0 to CP ( H 0 )

It is well-known that, in the Hilbert space, there is a natural action of the Abelian group C 0 = C { 0 } given by the following.
| ψ λ | ψ = ϱ e i θ | ψ with ϱ > 0 ,
According to the statistical interpretation of the wave function, it turns out that pure states are better described by equivalence classes of vectors in H with respect to this action, i.e., rays in the Hilbert space H .
Let us show how, for finite-dimensional systems, the set of such equivalence classes can be given by the structure of a manifold: consider an n-level quantum system with Hilbert space H and define H 0 = H { 0 } . By selecting an orthonormal basis { | e k } k = 1 , , n in H 0 , one may introduce a Cartesian coordinate system { x k , y k } k = 1 , , n on H 0 such that for any element | ψ H 0 one has the following:
| ψ = ψ k | e k = ( x k + i y k ) | e k ,
where Einstein’s summation convention over repeated indices is assumed here and in the following equations. The action defining the equivalence classes of pure quantum states may be described infinitesimally by means of two commuting linear vector fields given in Cartesian coordinates as follows:
Δ = x k x k + y k y k and Γ = y k x k x k y k ,
where Δ is the infinitesimal generator of dilations, while Γ is the infinitesimal generator of the multiplication by a global phase factor. Now, dilations define a distribution in which the integral curves foliate the Hilbert space. Let Φ Δ denote such foliation, then the quotient space H 0 / Φ Δ is the unit sphere.
S 2 n 1 ( H ) : = | ψ H 0 | ψ | ψ = 1 .
In the following, an element of S 2 n 1 is denoted as | ψ ) , whereas | ψ is a vector in H 0 . Now, let Φ Δ , Γ be the foliation corresponding to the distribution generated by Δ and Γ . Then, the quotient space H 0 / Φ Δ , Γ is identified with the complex projective space CP ( H 0 ) , defined as follows.
CP ( H 0 ) : = { λ | ψ | λ C 0 } .
Furthermore, there is a one-to-one correspondence between elements of CP ( H 0 ) and rank-one projectors given by the following.
[ ψ ] ρ ψ : = | ψ ψ | ψ | ψ .
On CP ( H 0 ) , it is convenient to work with complex homogeneous coordinates. To introduce such coordinates, let U j CP ( H 0 ) denote the open subset where ψ j 0 . Then, on U j , one may introduce the following coordinates.
ϕ j : U j C n 1 : [ ψ 1 , , ψ n ] z 1 , , z j 1 , z j + 1 , , z n , with z k = ψ k ψ j .
The set of ( U j , ϕ j ) , with j = 1 , , n , constitutes an atlas for CP ( H 0 ) . In such coordinates, the projection from H 0 onto CP ( H 0 ) can be explicitly expressed as the following:
π : H 0 CP ( H 0 ) : | ψ | ψ ] = 1 1 + | z | 2 z 1 ,
where we used the chart ϕ n for the sake of simplicity, and we introduce the following notations.
z = z 1 z 2 z n 1 and | z | 2 = z z .
To conclude the kinematical analysis of n-level systems, let us recall that CP ( H 0 ) is a Kähler manifold [25,28,31]; in fact, considering the homogeneous coordinates in Equation (8), one can introduce the 1-form of the following.
θ FS = i 1 1 + | z | 2 ( z ¯ k , 1 ) 1 1 + | z | 2 d z k 0 + z k 1 d 1 1 + | z | 2 = 2 i z ¯ k d z k z k d z ¯ k 1 + | z | 2 .
Then, the symplectic form ω FS on CP ( H 0 ) is given by the exterior derivative ω FS = d θ FS , which reads as follows.
ω FS = i ( 1 + | z | 2 ) 2 ( 1 + | z | 2 ) d z ¯ k d z k 1 2 ( z ¯ l d z l + z l d z ¯ l ) ( z ¯ k d z k z k d z ¯ k ) .
Moreover, as proven by Wootters in [32], a natural notion of distance between quantum states is given by the Fubini–Study metric on CP ( H 0 ) , which has the following form [25,27,28,31]:
g FS = ( 1 + | z | 2 ) 2 ( 1 + | z | 2 ) d z ¯ k S d z k + 1 2 ( z ¯ k d z k z k d z ¯ k ) ( z ¯ l d z l z l d z ¯ l ) 1 2 ( z ¯ k d z k + z k d z ¯ k ) ( z ¯ l d z l + z l d z ¯ l ) ,
where d z ¯ k S d z k = d z ¯ k d z k + d z k d z ¯ k . On the other hand, the ( 1 , 1 ) -tensor of the following:
J FS = 1 i d z k z k d z ¯ k z ¯ k
defines a complex structure on CP ( H 0 ) , such that the quadruple ( CP ( H 0 ) , ω FS , g FS , J FS ) is a Kähler manifold [25,27,28,31].

2.2. Dynamics: Symplectic Hamiltonian Systems on CP ( H 0 )

Having defined the canonical projection π from H 0 to CP ( H 0 ) , one may project the dynamics of the system as follows: First, we recall that the Schrödinger dynamics may be viewed as a classical Hamiltonian system on H 0 with the symplectic structure given by the imaginary part of the Hermitian scalar product and the generating Hamiltonian function being the expectation value e H = ψ | H | ψ of the Hamiltonian operator H . The corresponding Hamiltonian vector field then reads as follows:
X H = i ψ ¯ j H k j ψ ¯ k i H j k ψ j ψ k ,
where H k j is the entry of H [25]. Now, because [ X H , Δ ] = [ X H , Γ ] = 0 , i.e., Δ and Γ are symmetries of the Schrödinger dynamics, it is possible to project X H onto the dynamics on CP ( H 0 ) [25,28,33]. It turns out that the projected dynamics X e H is again a Hamiltonian vector field with respect to the symplectic structure ω FS , with generating Hamiltonian function e H F ( CP ( H 0 ) ) , where here and in the following F ( M ) stands for the smooth functions on M for any smooth manifold M; that is, it satisfies the following.
ω FS ( X e H , · ) = d e H .
In order to express Equation (15) in complex homogeneous coordinates on CP ( H 0 ) , we start by writing the expectation value of the observable H in coordinates such as the following:
e H = [ ψ | H | ψ ] = 1 1 + | z | 2 ( z , 1 ) H 1 V V H 2 z 1 = 1 1 + | z | 2 z H 1 z + z V + V z + H 2 ,
where H 1 is an ( n 1 ) × ( n 1 ) –dimensional matrix, V is an ( n 1 ) –dimensional column vector and H 2 a real quantity. Thus, by using Equations (11) and (15), we obtain the following:
X e H = X z k z k + X z ¯ k z ¯ k ,
where the component X z ¯ k is the complex conjugate of X z k , and the following is the case.
X z k = i z k V ¯ l z l | H 1 | l k z l + H 2 z k V k .
Therefore, the integral curves of this Hamiltonian vector field are solutions to the Hamiltonian equations of motion:
z ˙ k = i ( 1 + | z | 2 ) e H z ¯ k + z k z ¯ l e H z ¯ l = i z k V ¯ l z l | H 1 | l k z l + H 2 z k V k ,
which is also called the matrix Riccati equation [34]. Hence, the matrix Riccati equation is simply the coordinate expression of the projection of the Schrödinger equation onto the complex projective space.
In addition, one may use the symplectic structure to introduce a Poisson bracket on CP ( H 0 ) [25,28]: given the expectation values e A and e B associated with the observables A and B , one defines the following.
{ e A , e B } ω FS = ω FS ( X e A , X e B ) .
In complex homogeneous coordinates this bracket reads as follows.
{ e A , e B } ω FS = i ( 1 + | z | 2 ) e A z k e B z ¯ k e A z ¯ k e B z k + z k e A z k z ¯ l e B z ¯ l z ¯ l e A z ¯ l z k e B z k .
Furthermore, after some calculation, it is possible to prove that the Poisson bracket satisfies { e A , e B } ω FS = e 1 i [ A , B ] , where [ A , B ] = A B B A . Therefore, one has a clear connection between the Poisson bracket and the quantum commutator. In particular, considering the Hamiltonian of the system H with expectation value e H , the evolution of the expectation value e A of an arbitrary observable A is given by the following.
d e A d t = { e A , e H } ω FS = e 1 i [ A , H ] .
This result implies immediately that, in the time-independent case, e H is a first integral of the flow, i.e., the expectation value of the Hamiltonian is conserved. In addition, the expectation value of any observable commuting with H is a first integral too.
Proceeding in parallel with the above construction of Hamiltonian vector fields and of the Poisson bracket by means of the symplectic structure, we are going to define now the gradient vector fields Y e H and the Jordan bracket by using the Fubini–Study metric.
Gradient vector fields are defined as follows:
g FS ( Y e H , · ) = d e H .
and it is direct to verify that
J FS ( X e H ) = Y e H ,
with J FS given in (13). Therefore, by employing Equations (17) and (18), one finds the explicit coordinate expression for the gradient vector field, which reads as follows.
Y e H = 1 z k V ¯ l z l | H 1 | l k z l + H 2 z k V k z k + 1 z ¯ l V l z ¯ k z ¯ l | H 1 | k l + H 2 z ¯ k V ¯ k z ¯ k .
On the other hand, the Jordan bracket between the expectation values e A and e B is defined as follows:
{ e A , e B } g FS = g FS ( Y e A , Y e B ) ,
and in the coordinates, we have the following.
{ e A , e B } g FS = 1 ( 1 + z ¯ l z l ) e A z k e B z ¯ k + e A z ¯ k e B z k + z k e A z k z ¯ l e B z ¯ l + z ¯ l e A z ¯ l z k e B z k .
Finally, the Jordan bracket is connected with the dispersion and the correlation of the observables. That is, for every couple of observables A and B , their uncertainties and correlations are given by the following:
σ A 2 = e A 2 e A 2 = 2 { e A , e A } g FS
and
σ A B = 1 2 e [ A , B ] + e A e B = 2 { e A , e B } g FS ,
with [ A , B ] + = A B + B A being the anti-commutator. Thus, the Riemannian metric on CP ( H 0 ) takes into account part of the probabilistic character of quantum mechanics [35].

Example: The Conservative Qubit

As an example, let us consider a qubit in which CP ( H 0 ) is two-dimensional, and each point is given by [ ψ 1 , ψ 2 ] . For this case, the complex projective space may be thought of as the unit sphere:
S 2 = { ( x 1 , x 2 , x 3 ) R 3 | ( x 1 ) 2 + ( x 2 ) 2 + ( x 3 ) 2 = 1 } .
One may use the homogeneous coordinates on CP ( H 0 ) to induce coordinates on S 2 as follows: take each coordinate patch U j CP ( H 0 ) where ψ j 0 , j = 1 , 2 and introduce the homogeneous coordinates on U 1 (resp. U 2 ) as defined above. Thus, we have the following.
ϕ 1 : [ ψ 1 , ψ 2 ] z = ψ 1 ψ 2 resp . ϕ 2 : [ ψ 1 , ψ 2 ] ζ = ψ 2 ψ 1 .
Then, by writing z = z R + i z I and using the stereographic projection from the north pole of the sphere, one obtains the corresponding point ( x 1 , x 2 , x 3 ) S 2 by the following.
x 1 = 2 z R 1 + | z | 2 , x 2 = 2 z I 1 + | z | 2 , x 3 = 1 + | z | 2 1 + | z | 2 .
For this example, let us consider the Hamiltonian operator as the matrix.
H q = H 1 V V ¯ H 2 ,
Thus, the expectation value of H q is given by the following:
e H q = [ ψ | H q | ψ ] = 1 1 + | z | 2 ( z ¯ , 1 ) H 1 V V ¯ H 2 z 1 = 1 1 + | z | 2 H 1 | z | 2 + z ¯ V + V ¯ z + H 2 ,
and Hamilton’s Equation (19) reduces to the following.
z ˙ = i V ¯ z 2 ( H 1 H 2 ) z V .
To provide a qualitative description of the behavior of this nonlinear system, one first realizes that the critical points of the flow are located at the following.
z ± = ( H 1 H 2 ) ± ( H 1 H 2 ) 2 + 4 | V | 2 2 V ¯ ,
They are both centers. The phase portrait in the complex plane is depicted in Figure 1a, while the corresponding vector field on S 2 obtained by means of (32) is displayed in Figure 1b.
Another important quantity that may be easily calculated by means of Equation (28) is the uncertainty of the energy given by the following.
σ H q 2 = 1 ( 1 + | z | 2 ) 2 | V z ¯ 2 ( H 1 H 2 ) z ¯ V ¯ | 2 .
Clearly this is a positive quantity that only vanishes at the singular points z ± .

3. Contact Geometry and Dynamics of Dissipative n-Level Quantum Systems

In this section, we consider a particular class of dissipative quantum systems, which are those that admit a contact Hamiltonian description (see also [7,11,12,21,26,36,37,38,39] for detailed discussions on the strengths and limitations of this approach both in the classical and quantum settings).

3.1. Kinematics: From CP ( H 0 ) to CP ( H 0 ) × R

In order to introduce dissipation, we will work on the contactification of CP ( H 0 ) . Therefore, let us first recall some basic facts of contact geometry [40,41].
A ( 2 n + 1 ) -dimensional manifold M is said to be an exact contact manifold if it is endowed with a global differentiable 1-form η such that η ( d η ) n 0 everywhere on M . Then, η is called the contact form, and the contact structure on M is given by the (non-integrable) distribution of hyperplanes D = ker ( η ) . In order to introduce a contact manifold for finite-dimensional quantum mechanical systems, we perform a contactification of the space of pure quantum states CP ( H 0 ) . Indeed, since CP ( H 0 ) is an exact symplectic manifold, one can consider the extended space CP ( H 0 ) × R , which carries a natural contact structure given as the kernel of the following contact form:
η = d S θ FS ,
where S is the global coordinate function on the fiber R , where S has unit of action, and we use a slight abuse of notation by indicating as θ FS the pullback of this 1-form on CP ( H 0 ) by means of the natural projection from CP ( H 0 ) × R onto the first argument. Let us observe that it is possible to define a Sasakian structure on CP ( H 0 ) × R starting from the Kähler structure on CP ( H 0 ) [40,42]. However, we do not need such structure in the following; therefore, we defer its analysis to a future contribution.

3.2. Dynamics: Contact Hamiltonian Systems on CP ( H 0 ) × R

The contact form allows one to associate a Hamiltonian vector field X H with every smooth function H F ( CP ( H 0 ) × R ) , defined by the following:
i X H d η = d H ( £ ξ H ) η and i X H η = H .
In this case, H is called the contact Hamiltonian function. Now, we proceed to compute the coordinate expression of X H , which may be written as the following in general:
X H = X z k z k + X z ¯ k z ¯ k + X S S ,
where the component X z ¯ k is the complex conjugate of X z k . Using Equations (38) and (39), one finds that the components of the contact Hamiltonian vector field in these coordinates are as follows:
X z k = i ( 1 + | z | 2 ) H z ¯ k + z k z ¯ l H z ¯ l + z k 2 ( 1 + | z | 2 ) H S ,
X S = H 1 2 ( 1 + | z | 2 ) z ¯ k H z ¯ k + z k H z k ,
which implies that the integral curves are solutions to the system of the differential equations.
z ˙ k = i ( 1 + | z | 2 ) H z ¯ k + z k z ¯ l H z ¯ l + z k 2 ( 1 + | z | 2 ) H S , z ¯ ˙ k = i ( 1 + | z | 2 ) H z k + z ¯ k z l H z l + z ¯ k 2 ( 1 + | z | 2 ) H S , S ˙ = H 1 2 ( 1 + | z | 2 ) z ¯ k H z ¯ k + z k H z k .
From Equations (41) and (42), we can compute the evolution of any arbitrary real function F F ( CP ( H 0 ) × R ) to obtain the following:
d F d t = X H [ F ] = H £ ξ F + Λ ( d F , d H ) ,
where the bivector Λ is given by the following.
Λ = ( 1 + | z | 2 ) i z k z ¯ k + z k z k z ¯ l z ¯ l + 1 2 S z k z k + z ¯ k z ¯ k .
Since we are interested in characterizing dissipative systems that dissipate the energy of some declared reference Hamiltonian system in this work, we shall henceforth assume that the contact Hamiltonian H may be written as follows:
H = e H + f ( S ) ,
where the first term e H is the expectation value of the Hamiltonian operator H of the conservative reference system given in Equation (16), while the second term f ( S ) is an instrumental perturbation which will be taken to be linear in order to result in projectable dynamics so that the projected dynamics can be observed as an effective description of some dissipative dynamics. Therefore, considering the expectation value e A F ( CP ( H 0 ) ) of an arbitrary observable A of the reference system, from Equation (44) we obtain the following:
d e A d t = { e A , e H } ω FS + 1 2 ( 1 + | z | 2 ) z k e A z k + z ¯ k e A z ¯ k f ( S ) ,
where { · , · } ω F S is the Poisson bracket given in Equation (21). Hence, it should be clear that the expectation value of the energy e H is not necessarily preserved along the trajectories of the contact Hamiltonian system.
Finally we remark that, from Equations (43), the contact dynamics X H X ( CP ( H 0 ) × R ) is projectable onto a vector field X X ( CP ( H 0 ) ) if and only if the contact Hamiltonian function H is linear in S. This is because the equations for the variables z k that characterize the state of the reference system defined on CP ( H 0 ) are decoupled from the equation for the additional variable S only in this case. It is for this reason that we will focus on this case only, leaving a detailed analysis of the general case to future investigations.

Example: The Dissipative Qubit

As an example, let us introduce dissipation in the qubit system described in the previous section. To perform this, we consider the following contact Hamiltonian:
H = e H γ S ,
where e H is the expectation value of the Hamiltonian of the conservative qubit (Equation (34)), and γ is a positive real constant that quantifies the strength of the coupling between the conservative system and the environment.
The Hamiltonian equations of motion (43) associated with the contact Hamiltonian (48) is read as follows.
z ˙ = i V ¯ z 2 ( H 1 H 2 ) z V γ 2 z ( 1 + | z | 2 )
S ˙ = H 1 2 1 | z | 2 1 + | z | 2 ( V ¯ z + V z ¯ ) + | z | 2 1 + | z | 2 ( H 2 H 1 ) .
We observe that, as expected, Equation (49) provides an effective equation of motion for the variable z, describing the state of the reference system, which is decoupled from the equation for the additional variable S.
As in the conservative case, we may now study the qualitative behavior of the dissipative qubit, focusing only on Equation (49). For the sake of simplicity, we consider here an example the particular case H 1 = H 2 . The critical points z s satisfythe following.
i V ¯ z s 2 V γ 2 z s ( 1 + | z s | 2 ) = 0 .
To find solutions to this algebraic equation, first we assume that | z s | = 1 in order to obtain the reduced second-order equation of the following:
V ¯ z s 2 + i γ z s V = 0 ,
with the following solutions:
z s ( 1 , 2 ) = i γ ± Δ 2 V ¯ ,
where Δ : = 4 | V | 2 2 γ 2 0 . A direct computation shows that | z s ( j ) | = 1 for j = 1 , 2 if and only if Δ 0 ; therefore, these two critical points exist only in such case. When they exist, in the conservative limit γ 0 , one recovers the analogous solutions in Equation (36) for the case H 1 = H 2 . Additionally, in the dissipative qubit, one finds the following critical point:
z s ( 3 ) = 2 i V γ ,
which always exists.
To continue the qualitative study of the dynamics, we observe that the eigenvalues of the linearized system at the critical points z s ( 1 ) and z s ( 2 ) are always the same and they are given by the following.
λ ± ( 1 , 2 ) = γ 2 ± i Δ 2 γ 2 4 ,
While at the critical point z s ( 3 ) , we have the following eigenvalues:
λ ± ( 3 ) = 2 γ 2 ± 4 | V | 2 2 2 γ ,
which are always real. From this linearization, it is possible to conclude that there is a bifurcation depending on the value of γ :
(i)
For γ such that Δ < 0 , we have only the critical point z s ( 3 ) , which is a stable node;
(ii)
For γ such that Δ = 0 , we have two critical points z s ( 1 ) = z s ( 2 ) and z s ( 3 ) , and both are non-hyperbolic (eigenvalues { 0 , γ } );
(iii)
For γ such that Δ > 0 , we have three different critical points with the following behavior: z s ( 3 ) is always a saddle; the behavior of z s ( 1 ) and z s ( 2 ) depends on the term Δ 2 γ 2 4 . If this is positive, then the critical points are stable foci. On the other hand, if this term is non-positive, then they are stable nodes.
Note that even in this simple example, we obtain very interesting dynamical behavior. Moreover, we remark that these types of critical points cannot be obtained in the case of a conservative (unitary) evolution.
By way of example, we have depicted the phase portrait for a choice of γ such that Δ > 0 and Δ 2 γ 2 4 > 0 in Figure 2. As observed in first panel in Figure 2, we have two stable foci and a saddle. However, recalling that this chart covers only the part of the Bloch sphere S 2 excluding the north pole, we also depict the corresponding vector field on S 2 in the second panel in Figure 2 from which it is clear that there is an additional critical point at the north pole, which is an unstable node.
Finally, taking into account the definition of the uncertainty of the energy operator σ H q 2 = H q 2 H q 2 , one has the following.
σ H q 2 = 1 ( 1 + | z | 2 ) 2 | V z ¯ 2 ( H 1 H 2 ) z ¯ V ¯ | 2 .
We stress that, although the functional forms of (37) and (57) are the same, they are different functions of time since, in the two cases, z ( t ) evolves according to different dynamical equations.
.

4. The Contact Master Equation

So far we have considered dissipative quantum evolutions by using contact Hamiltonian systems defined on the manifold CP ( H 0 ) × R . However, the most common description of dissipative quantum systems is by means of density operators; therefore, we shall adapt our approach to pure states considered as rank-one projectors in this section.
Given a unital C * -algebra A , a state ρ on A is a continuous linear function in the dual A * of A such that for all observables a A , one has the following:
ρ ( a ) R , ρ ( a a ) 0 , and ρ ( I ) = 1 ,
where I A is the identity, and an observable a A is a self-adjoint element. In particular, for finite-dimensional systems, we may restrict our attention to the C * -algebra of n × n complex matrices, i.e., we may consider A n = M n ( C ) with n 2 . Thus, the space of states S of A n is identified as the following:
S : = ρ D n * A n * | ρ ( a a ) 0 , for all a D n , ρ ( I ) = 1 ,
where D n A n denotes the space of observables. The pairing map  μ : D n × S R is the evaluation of the state ρ on the self-adjoint element a , and it is given explicitly by the following:
μ ( a , ρ ) ρ ( a ) : = Tr { ρ a } ,
which corresponds to the mean value for the outcome of the measurement of the observable a when the system is in the state ρ . Furthermore, since there is a one-to-one correspondence between elements in D n * and D n , it follows that the space S may be decomposed as the following:
S = k = 1 n S k ,
where S k = { ρ S | rk ( ρ ) = k } . It is proven in References [30,43] that every S k is a homogeneous space for the Lie group S L ( A n ) ; thus, every S k admits the structure of a differential manifold.
Here, we are only interested in pure states, i.e., rank-one projectors in S 1 . Then, as we have already mentioned in Section 2, there is a one-to-one correspondence between states | ψ ] CP ( H 0 ) and rank-one projectors ρ S 1 given by the relation (6). Thus, in homogeneous coordinates, we have the following.
| ψ ] = 1 1 + | z | 2 z 1 , then ρ ψ = | ψ ] [ ψ | = 1 1 + | z | 2 z z z z 1 .
Now, starting from the expression for the density matrix in the second equation in Equation (62) and taking into account the equations of motion (43), one may deduce the equation of motion for ρ , which takes the following general form.
ρ ˙ = i z j H z k z ¯ k H z ¯ j H z ¯ j z j z l H z l H z k + z ¯ k z ¯ l H z ¯ l z l H z l + z ¯ l H z ¯ l + H S z j z ¯ k 1 + | z | 2 z j 2 1 | z | 2 1 + | z | 2 z ¯ k 2 1 | z | 2 1 + | z | 2 | z | 2 1 + | z | 2 .
We first observe from the above equation that Tr { ρ ˙ } = 0 independently of the choice of the contact Hamiltonian. Therefore, the trace of ρ is preserved along the evolution, which is a fundamental statistical condition for any admissible quantum evolution.
Remarkably, the evolution (63) can be written using two brackets, which is similar to the structure of the metriplectic and GENERIC formalisms [44,45,46,47,48] (although an approach to dissipative quantum systems directly based on the analogy with the GENERIC equation results in a different type of master equation; see [49,50]). Indeed, we have the following (see Appendix A for the proof):
ρ ˙ = i [ ρ , H ] + H S ρ , A ( z ) + ,
where the dissipative bracket [ · , · ] + is the anti-commutator, and we have defined the dissipative potential to be the Hermitian operator A ( z ) given by the following:
A ( z ) = A w ( z ) w ( z ) A ( z ) ,
where A is an arbitrary ( n 1 ) × ( n 1 ) -dimensional Hermitian matrix, and w is the ( n 1 ) -dimensional column vector given by the following.
w ( z ) = 1 1 + | z | 2 1 + | z | 2 I z z 1 2 1 + | z | 2 I z A z I A z .
A ( z ) is the real quantity.
A ( z ) = z A z | z | 2 .
Let us remark that in Equation (64) there is a conservative part given by the standard von Neumann equation in addition to a contact perturbation introducing dissipation in conservative reference dynamics. Moreover, in general, Equation (64) is coupled to the equation of motion for the variable S given by the following.
S ˙ = H 1 2 ( 1 + | z | 2 ) z ¯ k H z ¯ k + z k H z k ,
They decouple if and only if the contact Hamiltonian H is linear in S. In such case, one obtains from Equation (64) a Markovian equation for ρ , which we call the contact master equation. Note that Markovianity in this case should be understood as satisfying the Chapman–Kolmogorov law for dynamical evolution. This is the case we will address in the following, leaving a detailed study of the general case for future investigations.
Let us now compare the contact master Equation (64) with the GKLS equation normally employed to describe dissipative phenomena in quantum systems. For the GKLS equation, we have the following:
ρ ˙ = i [ ρ , H ] 1 2 j ρ , V j V j + + j V j ρ V j ,
where V j and j V j V j are bounded operators. First of all, we note that, by comparing the GKLS evolution in Equation (69) with the contact evolution in Equation (64), it is clear that the absence in the latter of the “jump term” j V j ρ V j (also known as the Choi–Kraus term [2]). This is expected because this term is the one responsible for the change of rank of the density matrix, and the contact evolution considered here preserves the purity of the states. However, while in the GKLS equation the jump term is needed in order to enforce conservation of the trace and the complete-positivity of the map, we observe that it is no longer needed in the contact evolution, for these two properties are both automatically satisfied (at least for pure states). Hence, we observe that the contact master equation provides a (nonlinear) description of coherent dissipative phenomena.
Furthermore, one may also compare the contact master equation with other descriptions normally employed to describe dissipative systems. For instance, one may consider the description of an optical MASER (Microwave Amplification by Stimulated Emission of Radiation) [51,52]. Considering, for simplicity, the case of a 2-state system, one introduces the non-normalized density matrix:
ρ = ψ 1 ψ ¯ 1 ψ 1 ψ ¯ 2 ψ 2 ψ ¯ 1 ψ 2 ψ ¯ 2 ,
in which the equation of motion reads as follows:
ρ ˙ = i [ ρ , H ] 1 2 ρ , Γ + ,
where H is a time-dependent Hermitian operator, and Γ is a diagonal operator with entries γ 1 , γ 2 > 0 that describe phenomenologically the radiative decays of the eigenstates of H to the ground state. Let us observe that the positivity of Γ implies that Tr { ρ ˙ } < 0 ; therefore, the evolution is not trace-preserving. We conclude that the pathologies of Equation (71) can be removed with at least two methods: On the one side, one can introduce the jump term Γ ρ Γ to obtain an equation of GKLS type, but this results in introducing dissipation of the rank; on the other side, the contact master equation provides an alternative option, one that fixes the pathologies of Equation (71) while preserving the coherence of the states. However, this comes at the price of introducing nonlinearities.

5. Application: Radiative Decay

The simplest dissipative quantum phenomenon that one can describe by the contact evolution is the radiative decay of a 2-level atom with levels | 1 and | 2 . To show this, let us consider as the conservative system the following Hamiltonian:
H = H 1 0 0 H 2
where H 1 and H 2 are the energies of the states, with H 1 > H 2 . Then, the decay of the particle may be modelled by the contact Hamiltonian:
H = e H γ S
where e H represents the expectation value of the Hamiltonian (72), and the constant damping factor γ > 0 describes phenomenologically the radiative decay from the state | 1 to | 2 . The corresponding contact Hamiltonian equations of motion can be obtained from Equation (43) and read as follows.
z ˙ = i ( H 1 H 2 ) z γ 2 z ( 1 + | z | 2 ) ,
S ˙ = H | z | 2 1 + | z | 2 ( H 1 H 2 ) ,
In other words, we have a decoupled system of differential equations as expected; therefore, we can focus on the first equation only, for which its solution can be given explicitly by the following:
z ( t ) = e i φ 0 e γ t + 2 κ 0 1 e i ( H 1 H 2 ) t ,
where the phase φ 0 and the real constant κ 0 are defined by the initial condition z ( 0 ) = z 0 . Let us recall that because the homogeneous coordinates z CP ( H 0 ) have been mapped to the sphere via the stereographic projection, in order to consider the state with energy H 1 as the initial condition, one has to consider z 0 as the “point at infinity” in the plane, corresponding to the north pole of the sphere. In Equation (76), this corresponds to taking κ 0 0 .
As an illustration, we depict the behaviour of the solutions of Equation (76) in Figure 3. In Figure 3a, we show the phase portrait in homogeneous coordinates in the plane corresponding to the stereographic projection from the north pole of the sphere. Here, we have a stable focus at the origin; consequently, all nearby solutions evolve towards the origin. On the other hand, considering the coordinates obtained by stereographic projection from the south pole, one has an unstable focus at the origin, as shown in Figure 3b. Finally, both charts form an atlas for the Bloch sphere S 2 , and the behaviors of the vector field on such sphere are displayed in Figure 3c. All this is in agreement with our physical interpretation because the north pole corresponds to the excited state with energy H 1 and the south pole to the state with lower energy H 2 , i.e., we are dissipating energy until the system finally decays in the lower state.
To observe clearly the continuous dissipation of the energy, one may compute the evolution of the expectation value of the conservative Hamiltonian operator (72) by using the solution (76) to obtain the following.
e H ( t ) = 1 1 + | z ( t ) | 2 H 1 | z ( t ) | 2 + H 2 = H 2 + ( H 1 H 2 ) e ( γ t + 2 κ 0 ) .
Then, it is not difficult to prove that for the initial condition κ 0 0 we have e H ( t = 0 ) = H 1 and lim t + e H ( t ) = H 2 , thus describing the radiative decay from the excited to the ground state. In Figure 4a, we display the continuous transition between these states, for which its speed depends on the value of γ . Furthermore, we remark that one may equivalently describe excitations within this model by simply exchanging γ γ and then having the transition from the lower state H 2 to the excited state H 1 .
One may also compute the uncertainty of the energy operator σ H 2 = H 2 H 2 , namely the following:
σ H 2 = | z ( t ) | 2 ( 1 + | z ( t ) | 2 ) 2 ( H 1 H 2 ) 2 = 1 e ( γ t + 2 κ 0 ) e γ t 2 κ 0 ( H 1 H 2 ) 2 ,
Its evolution is plotted in Figure 4b. From this figure, one may observe that the uncertainty of the energy starts at zero, then increases up to the maximum value 1 4 ( H 1 H 2 ) 2 at t = 1 γ ln 2 . After this maximum, the curve decreases asymptotically to zero. Since σ H 2 is the statistical fluctuation around the expectation value, it follows that a measurement of e H with complete certainty is only possible in principle at the initial and final times.
Another relevant quantity is the probability of a transition between the states, denoted by P ( t ) and given by P ( t ) = | [ ψ 0 | ψ t ] | 2 , where | ψ 0 ] is the initial state and | ψ t ] corresponds to the evolved normalized state at time t. Using (62), we obtain that P ( t ) is given by the following:
P ( t ) = | z ( t ) | 2 1 + | z ( t ) | 2 ,
where we considered the initial condition κ 0 0 in this expression. Moreover, by means of Equation (47), we can compute the rate of dissipation of this probability, which is d P d t = γ P , meaning that one has an exponential decay P ( t ) = e γ t for the probability of a transition.
We remark again at this point that radiative decays are usually represented as sudden “jumps” and not as continuous transitions. However, recent experimental and theoretical investigations [5,6] have shown that it is possible to observe quantum jumps as continuous processes that preserve the coherence of the state. Here, we have observed in this example that the contact master equation yields another possibility to model quantum decays (or excitations) as coherent and continuous processes.

6. Conclusions and Perspectives

In this work, we have put forward a novel approach to the description of dissipative quantum systems based on the geometric approach to quantum mechanics and on the analogy with the description of classical dissipative systems based on contact Hamiltonian dynamics. Of special importance is the fact that the obtained evolution dissipates the expectation value of the energy of the reference system while preserving the purity of the states. In this manner, it yields a method to describe coherent dissipative dynamics, which inevitably escapes more standard approaches such as the GKSL equation. Among the possible applications, we have considered here, in particular, the important case of radiative decay for a 2-level system, both because of its theoretical importance in understanding quantum mechanics and because of recent experiments that point to an explanation in terms of the existence of coherent quantum trajectories for these systems [5,6].
At this point, there is a number of interesting questions open for future investigations. In particular, in this effort, we have considered the additional variable S as a “means to an end”, that is, as an effective tool that we have employed in order to produce a dissipative dynamics on the reference system in a geometric manner. However, the term γ S as we have used it has the dimensions of an energy, and similar terms appear in the thermodynamic literature as the actual interchange of heat between the reference system and the environment [53,54]. Therefore, a deeper connection with an energy conservation principle may be responsible for our construction.
Still, concerning the variable S, in this work we have dealt only with contact Hamiltonian functions that are linear in S. This is because we can decouple the dynamics of S from that of the quantum states and obtain an evolution that is purely quantum in the sense that it does not depend on other variables.
However, we have remarked that the evolution corresponding to the most general contact Hamiltonian still preserves the trace of the density operator. This is interesting because, in such case, one expects to obtain a non-Markovian evolution mapping pure states into pure states.
Another interesting problem is the extension of the present approach to infinite dimensional quantum systems. Solving this problem in all its generality is not an easy task; however, one may be restricted to the problem of looking for the immersion of a contact manifold into the Hilbert space such that the evolution of the states is parametrized by the contact evolution.
Furthermore, given the importance of the Schrödinger dynamics in the description of quantum systems, one may wonder whether there is a Schrödinger-like equation associated with the contact evolution introduced here. Indeed, for some particular cases, it is possible to construct a Schrödinger equation; for instance, one may reproduce the contact evolution (49) by means of the projection (8) starting from the nonlinear Schrödinger equations:
i ψ ˙ 1 ψ ˙ 2 = H 1 V V ¯ H 2 ψ 1 ψ 2 + γ 2 1 ψ 1 ψ 2 ψ ¯ 1 ψ ¯ 2 1 ψ 1 ψ 2 ,
or
i ψ ˙ 1 ψ ˙ 2 = H 1 V V ¯ H 2 ψ 1 ψ 2 + γ 2 1 / 2 0 ψ ¯ 1 ψ ¯ 2 1 / 2 ψ 1 ψ 2 .
It is not difficult to see that, in both equations, the normalization and the phase are not invariant, and they are involved in the evolution, which introduces nonlinearities in the equations.
Finally, a major motivation for our work and for further investigation is the applications of formalism and a deeper comparison with existing approaches. In the present work, we have deliberately only scratched the surface of some applications, but we expect that, as it happened in the classical case, further systems may be analyzed from this perspective, and new tantalizing results may possibly be found.

Author Contributions

Conceptualization and writing—original draft preparation, H.C.-P.; formal analysis, visualization, writing—review and editing H.C.-P., A.B. and A.G.-C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Proof of the Contact Master Equation Form

In this appendix, we show that the contact master equation in Equation (63) can be written as in Equation (64). To perform this, let us assume that there exists an operator A ( z ) of the following form:
A ( z ) = A w ( z ) w ( z ) A ( z ) ,
where A is an arbitrary ( n 1 ) × ( n 1 ) -dimensional Hermitian matrix, w ( z ) is an ( n 1 ) -dimensional column vector and A ( z ) is a real number. Our goal is to prove that the anticommutator ρ , A ( z ) + of this operator and the density operator ρ is given by the following:
ρ , A ( z ) + = 1 1 + | z | 2 z z z 2 ( 1 | z | 2 ) z 2 ( 1 | z | 2 ) | z | 2 ,
where the density matrix is of the following form.
ρ = 1 1 + | z | 2 z z z z 1 .
A direct computation of the anticommutator yields the following:
ρ , A ( z ) + = 1 1 + | z | 2 z z , A + + z w ( z ) + w ( z ) z ( z z + I ) w ( z ) + ( A + A ( z ) I ) z w ( z ) ( z z + I ) + z ( A + A ( z ) I ) z w ( z ) + w ( z ) z + 2 A ( z ) ,
where I represents the ( n 1 ) × ( n 1 ) -dimensional identity matrix. Equating the right hand sides of Equations (A2) and (A4) provides the following system of equations.
z z , A + + z w ( z ) + w ( z ) z = z z ,
( z z + I ) w ( z ) + ( A + A ( z ) I ) z = z 2 ( 1 | z | 2 ) ,
w ( z ) ( z z + I ) + z ( A + A ( z ) I ) = z 2 ( 1 | z | 2 ) ,
z w ( z ) + w ( z ) z + 2 A ( z ) = | z | 2 .
Let us impose the following conditions:
A = A , A * ( z ) = A ( z ) ,
which guarantee that (A7) is the conjugate transpose of (A6), i.e., Equation (A7) is no longer independent.
We now solve for w ( z ) from Equation (A6) and obtain the following:
w ( z ) = I z z 1 + | z | 2 1 2 ( 1 | z | 2 ) I A ( z ) I A z ,
where the inverse of the matrix I + z z is given by the following.
I + z z 1 = I z z 1 + | z | 2 .
We then insert (A10) into (A8) and solve for A ( z ) and obtain the following.
A ( z ) = z A z | z | 2 .
Equation (A12) allows us to express w ( z ) given in (A10) in terms of the matrix A as follows.
w ( z ) = I z z 1 + | z | 2 1 2 ( 1 + | z | 2 ) I ( z A z ) I A z .
So far, we have obtained expressions for A ( z ) and w ( z ) and Equations (A12) and (A13) respectively, for which Equations (A6)–(A8) are satisfied. The remaining Equation (A5) is trivially satisfied by using the following condition.
z w ( z ) + w ( z ) z = [ z z , A ] + + z z ,
When the above is inserted into (A5), we yield the following.
z z , A + + z w ( z ) + w ( z ) z = z z , A + [ z z , A ] + + z z = z z .
As a result, Equation (A5) is trivially satisfied, without any further requirements on the Hermitian matrix A . To summarize, we have the operator A ( z ) in (A1), with A being any Hermitian matrix, and A ( z ) and w ( z ) given in (A12) and (A13), respectively, satisfies (A2), thus concluding the proof.

References

  1. Landau, L. Das dämpfungsproblem in der wellenmechanik. Z. Phys. 1927, 45, 430–441. [Google Scholar] [CrossRef]
  2. Chruściński, D.; Pascazio, S. A brief history of the GKLS equation. arXiv 2017, arXiv:1710.05993. [Google Scholar] [CrossRef]
  3. Gorini, V.; Kossakowski, A.; Sudarshan, E.C.G. Completely positive dynamical semigroups of n-level systems. J. Math. Phys. 1976, 17, 821–825. [Google Scholar] [CrossRef]
  4. Lindblad, G. On the generators of quantum dynamical semigroups. Commun. Math. Phys. 1976, 48, 119–130. [Google Scholar] [CrossRef]
  5. Minev, Z.K.; Mundhada, S.O.; Shankar, S.; Reinhold, P.; Gutiérrez-Jáuregui, R.; Schoelkopf, R.J.; Mirrahimi, M.; Carmichael, H.J.; Devoret, M.H. To catch and reverse a quantum jump mid-flight. Nature 2019, 570, 200–204. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Snizhko, K.; Kumar, P.; Romito, A. Quantum Zeno effect appears in stages. Phys. Rev. Res. 2020, 2, 033512. [Google Scholar] [CrossRef]
  7. Bravetti, A.; Cruz, H.; Tapias, D. Contact Hamiltonian mechanics. Ann. Phys. 2017, 376, 17–39. [Google Scholar] [CrossRef] [Green Version]
  8. Bravetti, A.; Garcia-Chung, A. A geometric approach to the generalized Noether theorem. J. Phys. Math. Theor. 2021, 54, 095205. [Google Scholar] [CrossRef]
  9. Bravetti, A.; Tapias, D. Liouville’s theorem and the canonical measure for nonconservative systems from contact geometry. J. Phys. Math. Theor. 2015, 48, 245001. [Google Scholar] [CrossRef] [Green Version]
  10. Bravetti, A.; de León, M.; Marrero, J.C.; Padrón, E. Invariant measures for contact Hamiltonian systems: Symplectic sandwiches with contact bread. J. Phys. Math. Theor. 2020, 53, 455205. [Google Scholar] [CrossRef]
  11. De León, M.; Valcázar, M.L. Infinitesimal symmetries in contact Hamiltonian systems. J. Geom. Phys. 2020, 153, 103651. [Google Scholar] [CrossRef] [Green Version]
  12. Gaset, J.; Gràcia, X.; Muñoz-Lecanda, M.C.; Rivas, X.; Román-Roy, N. New contributions to the Hamiltonian and Lagrangian contact formalisms for dissipative mechanical systems and their symmetries. Int. J. Geom. Methods Mod. Phys. 2020, 17, 2050090. [Google Scholar] [CrossRef]
  13. Lazo, M.; Paiva, J.; Frederico, G. Noether theorem for action-dependent Lagrangian functions: Conservation laws for non-conservative systems. Nonlinear Dyn. 2019, 97, 1125–1136. [Google Scholar] [CrossRef] [Green Version]
  14. Cannarsa, P.; Cheng, W.; Wang, K.; Yan, J. Herglotz’ generalized variational principle and contact type Hamilton-Jacobi equations. In Trends in Control Theory and Partial Differential Equations; Springer: Berlin/Heidelberg, Germany, 2019; pp. 39–67. [Google Scholar]
  15. Georgieva, B.; Guenther, R. First Noether-type theorem for the generalized variational principle of Herglotz. Topol. Methods Nonlinear Anal. 2002, 20, 261–273. [Google Scholar] [CrossRef] [Green Version]
  16. Georgieva, B.; Guenther, R.; Bodurov, T. Generalized variational principle of Herglotz for several independent variables. First Noether-type theorem. J. Math. Phys. 2003, 44, 3911–3927. [Google Scholar] [CrossRef]
  17. Liu, Q.; Torres, P.J.; Wang, C. Contact Hamiltonian dynamics: Variational principles, invariants, completeness and periodic behavior. Ann. Phys. 2018, 395, 26–44. [Google Scholar] [CrossRef]
  18. Vermeeren, M.; Bravetti, A.; Seri, M. Contact variational integrators. J. Phys. A Math. Theor. 2019, 52, 445206. [Google Scholar] [CrossRef] [Green Version]
  19. Wang, K.; Wang, L.; Yan, J. Implicit variational principle for contact Hamiltonian systems. Nonlinearity 2016, 30, 492. [Google Scholar] [CrossRef]
  20. Wang, K.; Wang, L.; Yan, J. Aubry–Mather theory for contact Hamiltonian systems. Commun. Math. Phys. 2019, 366, 981–1023. [Google Scholar] [CrossRef]
  21. Gaset, J.; Gràcia, X.; Muñoz-Lecanda, M.C.; Rivas, X.; Román-Roy, N. A contact geometry framework for field theories with dissipation. Ann. Phys. 2020, 414, 168092. [Google Scholar] [CrossRef] [Green Version]
  22. Ashtekar, A.; Schilling, T.A. Geometrical Formulation of Quantum Mechanics. In On Einstein’s Path; Springer: Berlin/Heidelberg, Germany, 1999; pp. 23–65. [Google Scholar]
  23. Brody, D.C.; Hughston, L.P. Geometric quantum mechanics. J. Geom. Phys. 2001, 38, 19–53. [Google Scholar] [CrossRef] [Green Version]
  24. Cruz-Prado, H.; Marmo, G.; Schuch, D. Nonlinear description of quantum dynamics. N-level quantum systems. J. Phys. Conf. Ser. 2020, 1612, 012010. [Google Scholar] [CrossRef]
  25. Ercolessi, E.; Marmo, G.; Mor, I.G. From the equations of motion to the canonical commutation relations. arXiv 2010, arXiv:1005.1164. [Google Scholar]
  26. Ciaglia, F.M.; Cruz, H.; Marmo, G. Contact manifolds and dissipation, classical and quantum. Ann. Phys. 2018, 398, 159–179. [Google Scholar] [CrossRef] [Green Version]
  27. Bengtsson, I.; Życzkowski, K. Geometry of Quantum States: An Introduction to Quantum Entanglement; Cambridge University Press: Cambridge, UK, 2017. [Google Scholar]
  28. Cariñena, J.F.; Ibort, A.; Marmo, G.; Morandi, G. Geometry from Dynamics, Classical and Quantum; Springer: Berlin/Heidelberg, Germany, 2015. [Google Scholar]
  29. Cariñena, J.; Clemente-Gallardo, J.; Jover-Galtier, J.; Marmo, G. Tensorial dynamics on the space of quantum states. J. Phys. Math. Theor. 2017, 50, 365301. [Google Scholar] [CrossRef] [Green Version]
  30. Ciaglia, F.M.; Di Cosmo, F.; Ibort, A.; Laudato, M.; Marmo, G. Dynamical vector fields on the manifold of quantum states. Open Syst. Inf. Dyn. 2017, 24, 1740003. [Google Scholar] [CrossRef] [Green Version]
  31. McDuff, D.; Salamon, D. Introduction to Symplectic Topology; Oxford University Press: Oxford, UK, 2017. [Google Scholar]
  32. Wootters, W.K. Statistical distance and Hilbert space. Phys. Rev. D 1980, 23, 357. [Google Scholar] [CrossRef]
  33. Marmo, G.; Saletan, E.J.; Simoni, A.; Vitale, B. Dynamical Systems: A Differential Geometric Approach to Symmetry and Reduction; John Wiley & Sons: Chichester, UK, 1985. [Google Scholar]
  34. Chaturvedi, S.; Ercolessi, E.; Marmo, G.; Mukunda, N.; Simon, R. Ray space ’riccati’ evolution and geometric phases for N-level quantum systems. Pramana 2007, 69, 317–327. [Google Scholar] [CrossRef] [Green Version]
  35. Cirelli, R.; Mania, A.; Pizzocchero, L. Quantum mechanics as an infinite-dimensional Hamiltonian system with uncertainty structure: Part I. J. Math. Phys. 1990, 31, 2891–2897. [Google Scholar] [CrossRef]
  36. Cariñena, J.F.; Guha, P. Nonstandard Hamiltonian structures of the liénard equation and contact geometry. Int. J. Geom. Methods Mod. Phys. 2019, 16, 1940001. [Google Scholar] [CrossRef]
  37. De León, M.; Lainz, M. A review on contact Hamiltonian and Lagrangian systems. arXiv 2020, arXiv:2011.05579. [Google Scholar]
  38. Guha, P.; Block, J.; Ghose-Choudhury, A. Generalized conformal Hamiltonian dynamics and the pattern formation equations. J. Geom. Phys. 2018, 134, 195–208. [Google Scholar] [CrossRef]
  39. Sloan, D. Dynamical similarity. Phys. Rev. D 2018, 97, 123541. [Google Scholar] [CrossRef] [Green Version]
  40. Blair, D.E. Riemannian Geometry of Contact and Symplectic Manifolds; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  41. Geiges, H. An Introduction to Contact Topology; Volume 109 of Cambridge Studies in Advanced Mathematics; Cambridge University Press: Cambridge, UK, 2008; ISBN 978-0-521-86585-2. [Google Scholar]
  42. Boyer, C.; Galicki, K. Sasakian Geometry; Oxford University Press: Oxford, UK, 2008. [Google Scholar]
  43. Chruściński, D.; Ciaglia, F.M.; Ibort, A.; Marmo, G.; Ventriglia, F. Stratified manifold of quantum states, actions of the complex special linear group. Ann. Phys. 2019, 400, 221–245. [Google Scholar] [CrossRef] [Green Version]
  44. Grmela, M.; Öttinger, H.C. Dynamics and thermodynamics of complex fluids. I. Development of a general formalism. Phys. Rev. E 1997, 56, 6620. [Google Scholar] [CrossRef]
  45. Guha, P. Metriplectic structure, leibniz dynamics and dissipative systems. J. Math. Anal. Appl. 2007, 326, 121–136. [Google Scholar] [CrossRef] [Green Version]
  46. Morrison, P. Thoughts on brackets and dissipation: Old and new. J. Phys. Conf. Ser. 2009, 169, 012006. [Google Scholar] [CrossRef] [Green Version]
  47. Morrison, P.J. Bracket formulation for irreversible classical fields. Phys. Lett. A 1984, 100, 423–427. [Google Scholar] [CrossRef]
  48. Pavelka, M.; Klika, V.; Grmela, M. Multiscale Thermo-Dynamics: Introduction to GENERIC; Walter de Gruyter GmbH & Co KG: Berlin, Germany, 2018. [Google Scholar]
  49. Öttinger, H.C. Nonlinear thermodynamic quantum master equation: Properties and examples. Phys. Rev. A 2010, 82, 052119. [Google Scholar] [CrossRef] [Green Version]
  50. Öttinger, H.C. The geometry and thermodynamics of dissipative quantum systems. EPL Europhys. Lett. 2011, 94, 10006. [Google Scholar] [CrossRef] [Green Version]
  51. Carmichael, H. An Open Systems Approach to Quantum Optics: Lectures Presented at the Université Libre de Bruxelles, October 28 to November 4, 1991; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2009. [Google Scholar]
  52. Lamb, W.E., Jr. Theory of an optical maser. Phys. Rev. 1964, 134, A1429. [Google Scholar] [CrossRef]
  53. Eberard, D.; Maschke, B.; Van Der Schaft, A. An extension of Hamiltonian systems to the thermodynamic phase space: Towards a geometry of nonreversible processes. Rep. Math. Phys. 2007, 60, 175–198. [Google Scholar] [CrossRef] [Green Version]
  54. Simoes, A.A.; De León, M.; Valcázar, M.L.; DeDiego, D.M. Contact geometry for simple thermodynamical systems with friction. Proc. R. Soc. A 2020, 476, 20200244. [Google Scholar] [CrossRef]
Figure 1. (a) Phase portrait of the Hamiltonian vector field (35) on the plane C . (b) The same vector field on the sphere S 2 using the map (32). In both cases, the values of the parameters are H 1 = 4 , H 2 = 2 and V = 1 + i .
Figure 1. (a) Phase portrait of the Hamiltonian vector field (35) on the plane C . (b) The same vector field on the sphere S 2 using the map (32). In both cases, the values of the parameters are H 1 = 4 , H 2 = 2 and V = 1 + i .
Quantumrep 03 00042 g001
Figure 2. Phase portrait of the dissipative qubit with H 1 = H 2 both in the complex plane and on S 2 . The parameters chosen for this representation are the following: V = 1 + i and γ = 1 . One may observe the two stable foci and the saddle point both in the plane and in S 2 together with the unstable node at the north pole of the sphere, which cannot be observed in this chart on the plane. More details in the main text.
Figure 2. Phase portrait of the dissipative qubit with H 1 = H 2 both in the complex plane and on S 2 . The parameters chosen for this representation are the following: V = 1 + i and γ = 1 . One may observe the two stable foci and the saddle point both in the plane and in S 2 together with the unstable node at the north pole of the sphere, which cannot be observed in this chart on the plane. More details in the main text.
Quantumrep 03 00042 g002
Figure 3. Phase portrait of the vector field associated with the quantum decay from the higher energy H 1 = 4 to the lower energy H 2 = 2 with damping parameter γ = 1 . In (a) the vector field is plotted using the chart ϕ 1 with coordinates z C defined in Equation (32). In (b) this vector field is plotted using the chart ϕ 2 with coordinates ζ = 1 / z defined in Equation (32). Finally in (c) the vector field is plotted in the Bloch sphere.
Figure 3. Phase portrait of the vector field associated with the quantum decay from the higher energy H 1 = 4 to the lower energy H 2 = 2 with damping parameter γ = 1 . In (a) the vector field is plotted using the chart ϕ 1 with coordinates z C defined in Equation (32). In (b) this vector field is plotted using the chart ϕ 2 with coordinates ζ = 1 / z defined in Equation (32). Finally in (c) the vector field is plotted in the Bloch sphere.
Quantumrep 03 00042 g003
Figure 4. (a) Continuous transition from the state with energy H 1 = 4 to the state H 2 = 2 with damping parameter γ = 1 . (b) Time evolution of the uncertainty of the energy operator with the same parameters.
Figure 4. (a) Continuous transition from the state with energy H 1 = 4 to the state H 2 = 2 with damping parameter γ = 1 . (b) Time evolution of the uncertainty of the energy operator with the same parameters.
Quantumrep 03 00042 g004
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cruz-Prado, H.; Bravetti, A.; Garcia-Chung, A. From Geometry to Coherent Dissipative Dynamics in Quantum Mechanics. Quantum Rep. 2021, 3, 664-683. https://doi.org/10.3390/quantum3040042

AMA Style

Cruz-Prado H, Bravetti A, Garcia-Chung A. From Geometry to Coherent Dissipative Dynamics in Quantum Mechanics. Quantum Reports. 2021; 3(4):664-683. https://doi.org/10.3390/quantum3040042

Chicago/Turabian Style

Cruz-Prado, Hans, Alessandro Bravetti, and Angel Garcia-Chung. 2021. "From Geometry to Coherent Dissipative Dynamics in Quantum Mechanics" Quantum Reports 3, no. 4: 664-683. https://doi.org/10.3390/quantum3040042

APA Style

Cruz-Prado, H., Bravetti, A., & Garcia-Chung, A. (2021). From Geometry to Coherent Dissipative Dynamics in Quantum Mechanics. Quantum Reports, 3(4), 664-683. https://doi.org/10.3390/quantum3040042

Article Metrics

Back to TopTop