Next Article in Journal
Testing Primordial Black Hole Dark Matter with Atacama Large Millimeter Array Observations of the Gravitational Lens B1422+231
Previous Article in Journal
Arrival Time from Hamiltonian with Non-Hermitian Boundary Term
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Generalization of the Schrödinger Equation for Open Systems Based on the Quantum-Statistical Approach

by
Konstantin G. Zloshchastiev
Institute of Systems Science, Durban University of Technology, Durban 4000, South Africa
Universe 2024, 10(1), 36; https://doi.org/10.3390/universe10010036
Submission received: 29 November 2023 / Revised: 8 January 2024 / Accepted: 10 January 2024 / Published: 12 January 2024
(This article belongs to the Section Foundations of Quantum Mechanics and Quantum Gravity)

Abstract

:
Within the framework of the quantum-statistical approach, utilizing both non-Hermitian Hamiltonian and Lindblad’s jump operators, one can derive various generalizations of the von Neumann equation for reduced density operators, also known as hybrid master equations. If one considers the evolution of pure states only, i.e., disregarding the coherence between states and spontaneous transitions from pure to mixed states, then one can resort to quantum-mechanical equations of the Schrödinger type. We derive them from the hybrid master equations and study their main properties, which indicate that our equations have a larger range of applicability compared to other generalized Schrödinger equations proposed hitherto. Among other features, they can describe not only systems which remain in the stationary eigenstates of the Hamiltonian as time passes, but also those which evolve from those eigenstates. As an example, we consider a simple but important model, a quantum harmonic oscillator driven by both Hamiltonian and non-Hamiltonian terms, and derive its classical limit, which turns out to be the damped harmonic oscillator. Using this model, we demonstrate that the effects of dissipative environments of different types can cancel each other, thus resulting in an effectively dissipation-free classical system. Another discussed phenomenon is whether a non-trivial quantum system can reduce to a classical system in free motion, i.e., without experiencing any classical Newtonian forces. This uncovers a large class of quantum-mechanical non-Hamiltonian systems whose dynamics are not determined by conventional mechanics’ potentials and forces, but rather come about through quantum statistical effects caused by the system’s environment.

1. Introduction

Any realistic quantum system can be influenced not only by its natural environment, but can also be affected by the process of measurement, i.e., the interaction with an external measuring apparatus. Therefore, any observable quantum system is, or eventually becomes, open. The history of the theory of open quantum systems (OQS) goes back a few decades, counting from the early studies of dissipative systems based on quantizing classical models with damping, viscosity, and other phenomenologically introduced characteristics of dissipation [1].
It soon became clear that these initial attempts were quite limited in terms of their physical relevance because the role, complexity, and diversity of dissipative phenomena drastically increase, both qualitatively and quantitatively, in quantum systems, compared to classical ones. Therefore, one cannot take them all into account by simply starting from classical dissipative models and quantizing them in a “textbook” manner.
Moreover, quantum systems and experiments exist whose complete description cannot be achieved using a single-state vector because of the additional randomness or uncertainty regarding precisely which of their possible quantum states is prepared. These additional uncertainties could be caused by imperfectness of the devices used in experiments, which inevitably induce the classical randomness, or they could occur as a result of correlations of states due to quantum entanglement.
Thus, a consistent description of open quantum systems requires extensive usage of quantum statistical approaches [2,3], due to the occurrence of mixed states, i.e., statistical ensembles of pure states, as soon as a system “opens up”. This requires the OQS formalism to be upgraded from state vectors to density operators because the latter contain information not only about quantum-mechanical states, but also about the classical statistical probabilities thereof and correlations between them.
In conventional quantum mechanics, dealing with conservative systems, the transition from the state-vector to density-matrix descriptions is straightforward, and leads to the von Neumann equation, but, in the case of open systems, the situation is obviously more complicated, as a result of the above-mentioned diversity in dissipative processes and environments. This results in various generalizations of the von Neumann equation, usually aggregated into the OQS approach commonly referred to as the master equations for reduced density operators.
The Gorini–Kossakowski–Sudarshan–Lindblad (GKSL) master equation, or, simply, the Lindblad equation, became popular first (chronologically) [4,5], due to the relative simplicity of its foundations, which are based on the Markovian and rotating-wave approximations [3]. Despite the achieved substantial technical simplifications in open systems’ descriptions, it is obvious that these approximations significantly reduce the range of applicability of the GKSL approach.
An alternative framework and a range of models, which do not use the above-mentioned approximation, are based on non-Hermitian Hamiltonians (NH). Probably the earliest example can be found in the works of Feshbach, who showed that Hamiltonian operators can acquire anti-Hermitian terms as a result of partition between a subsystem and its environment [6,7].
Although master equations with non-Hermitian Hamiltonians have long been used—for example, in quantum optics [2]—the modern version of the NH approach was proposed relatively recently: when the post-selecting normalization procedure of a reduced density matrix was introduced [8,9] and transitions from pure to mixed states were considered [10]. This immediately resulted in a vast number of applications in nearly all branches of quantum physics; for examples from only the very recent literature, see [11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33].
In this paper, we derive various generalizations of the Schrödinger equation for open systems while using the quantum-statistical approach as a starting point and motivation. In Section 2, we give a brief introduction to the modern quantum-statistical approach based on unification of the NH- and GKSL-type master equations. In Section 3, we perform the dephasing and rank reduction procedures on density matrices, thus resorting from the density-matrix to the state-vector description. We discuss three types of generalized Schrödinger equations which occur as a result of dephasing, determine their equivalence class, and discuss special cases: NH-driven versus GKSL-driven evolution. Their basic properties are discussed in Section 4, where we focus on the evolution of stationary states, invariances, and relations between mean values. An instructive example is considered in Section 5, where we consider OQS models that reduce to the damped oscillator upon dephasing and averaging in the coherent states basis. Conclusions are drawn in Section 6.

2. Generalized von Neumann Equations

Here, we adopt the quantum-statistical approach that deals exclusively with the degrees of freedom of a subsystem, wherein all degrees of freedom of the environment have been averaged (“traced out”). We thus regard the (reduced) density operator of a subsystem as a primary operator for quantum-statistical description: the mean value of a quantum observable defined on a Hilbert space associated with the subsystem’s degrees of freedom is computed as the trace of the product of its operator with the density operator.
If density operator is normalized, i.e., its trace’s value is constant (can thus be set to 1) and preserved during the system’s time evolution, then we denote it by ρ ^ ; otherwise, we denote it by ϱ ^ . The relation between these two density matrices is straightforward:
ρ ^ = ϱ ^ / Tr ϱ ^ ,
therefore, definitions for statistical mean values are
O ^ Tr ( ρ ^ O ^ ) , O ^ ϱ Tr ( ϱ ^ O ^ ) = O ^ Tr ϱ ^ ,
where we use the subscript ϱ to differentiate the notations.
Furthermore, we assume that the dissipative environment’s influence upon a subsystem can be encoded in two different, but not mutually exclusive, ways.
The first approach is to use the anti-Hermitian terms of the subsystem’s Hamiltonian. Any non-Hermitian Hamiltonian operator can always be written as the sum of self- and skew-adjoint components:
H ^ = H ^ + + H ^ = H ^ + i Γ ^ ,
where we denote H ^ ± = 1 2 ( H ^ ± H ^ ) = ± H ^ ± and introduce a self-adjoint operator Γ ^ i H ^ , which is commonly referred to as the decay rate operator. It should be noted that, strictly speaking, the operator H ^ is not a true Hamiltonian, but rather an analytical way of writing down a set of two physically relevant operators H ^ ( H ^ + , Γ ^ ) , as we shall demonstrate below.
A second type of dissipation can be described using the non-Hamiltonian (Liouvillian) term, dubbed the Lindblad’s dissipator:
D ^ ρ ^ = k γ k L ^ k ρ ^ L ^ k 1 2 { ρ ^ , L ^ k L ^ k } = 1 2 k γ k [ L ^ k , ρ ^ L ^ k ] + [ L ^ k ρ ^ , L ^ k ] ,
where L ^ k is called the Lindblad, collapse, or jump operator; γ s are non-negative couplings which can be derived from the correlation functions of the environment (usually playing the role of relaxation rates for different decay modes); and the summation index runs k = 1 , , N 2 1 , where N is the dimension of the system [4,5]. Thus, D ^ is a Liouvillian-type operator in the quantum dynamical semigroup form, which is traceless, linear with respect to the density matrix, and quadratic with respect to the jump operators.
Furthermore, the general equation of motion for the density matrix was proposed in [9]. For the normalized density matrix, it reads as
d d t ρ ^ = i H ^ ρ ^ ρ ^ H ^ + i H ^ H ^ ρ ^ + k γ k L ^ k ρ ^ L ^ k 1 2 { ρ ^ , L ^ k L ^ k } = i H ^ + , ρ ^ 1 Γ ^ , ρ ^ + 2 Γ ^ ρ ^ + k γ k L ^ k ρ ^ L ^ k 1 2 { ρ ^ , L ^ k L ^ k } ,
where square and curly brackets denote, respectively, the commutator and anticommutator, Γ ^ = tr ( ρ ^ Γ ^ ) , according to our notations for mean values in (2).
Furthermore, the right-hand side of Equation (5) is essentially a sum of the following three terms:
i H ^ + , ρ ^ , 1 Γ ^ , ρ ^ + 2 Γ ^ , k γ k L ^ k ρ ^ L ^ k 1 2 { ρ ^ , L ^ k L ^ k } ,
which means that it is the generalization of the quantum von Neumann equation (represented by the first term) based on the unification of the quantum-statistical equation for non-Hermitian Hamiltonians (second term) with the Lindblad master equation approach (third term). As discussed in Section 1, the Lindblad’s dissipator and non-Hermitian Hamiltonian parts have different physical origins and, thus, should describe different types of dissipation. Comparing the structure of Equation (5) to the conventional von Neumann equation, one can see that it is the operator H ^ + that plays the role of a Hamiltonian of the system, whereas the anti-Hermitian part of H ^ contributes, together with the Lindblad’s dissipator, to the non-Hamiltonian (Liouvillian) part. This reaffirms our earlier remark after Equation (3) about the physical meaning of H ^ .
Each of the three terms in (6) independently has a vanishing trace, which means that the trace of ρ ^ is conserved at all times, as it is necessary for the probabilistic interpretation. It should also be noted that the Lindblad term cannot be entirely absorbed into the non-Hermitian Hamiltonian term, in general, because it contains not only the anticommutator { ρ ^ , L ^ k L ^ k } , but also the term L ^ k ρ ^ L ^ k , which is required by the semigroup structure and tracelessness of the Lindblad’s dissipator.
Since the density operator is capable of describing not only pure but also mixed states, it is an important feature when dealing with dissipative systems. However, the case of pure states is still fairly special: using the master equation in (5), one can show that the evolution equation for the quantum purity reads as
d d t P = R = 4 Γ ^ P Tr ( ρ ^ 2 Γ ^ ) + 2 k γ k Tr ( ρ ^ L ^ k ρ ^ L ^ k ) Tr ( ρ ^ 2 L ^ k L ^ k ) ,
where P = tr ( ρ ^ 2 ) is the purity. This equation indicates that the purity evolution rate R is determined by the NH terms and Lindblad terms independently of each other and in different ways. If evolution starts from an initially pure state, for which ρ ^ 2 ( t 0 ) = ρ ^ ( t 0 ) = | Ψ 0 Ψ 0 | , then one obtains
R t = t 0 = 2 k γ k L ^ k Ψ 0 L ^ k Ψ 0 L ^ k L ^ k Ψ 0 ,
where we used the notation . Ψ = Ψ | . | Ψ for mean values in the Hilbert space of | Ψ s. This formula indicates that initially pure states continuously evolve into mixed ones due to the Lindblad term in the master equation, while its non-Hermitian Hamiltonian part does not affect this process (this, however, does not preclude non-Hermitian Hamiltonians alone from triggering spontaneous transitions from pure to mixed states in the presence of quantum fluctuations, cf. [10]).
Let us now consider the non-normalized density matrix. Its evolution equation can be derived similarly:
d d t ϱ ^ = i H ^ ϱ ^ ϱ ^ H ^ + k γ k L ^ k ϱ ^ L ^ k 1 2 { ϱ ^ , L ^ k L ^ k } = i H ^ + , ϱ ^ 1 Γ ^ , ϱ ^ + k γ k L ^ k ϱ ^ L ^ k 1 2 { ϱ ^ , L ^ k L ^ k } ,
where the trace conservation breaks due to the absence of the compensating term 2 Γ ^ ϱ ^ , cf. Equation (6). The dynamical value of the density matrix’s trace is no longer an identity but a solution to the first-order differential equation
d d t Tr ϱ ^ = 2 i Tr ϱ ^ H ^ = 2 Γ ^ ϱ ,
where, in the last step, we assume the notations in (2).
The existence of two types of density operators, normalized and non-normalized, can be useful in certain applications, for example, when a system can evolve in either a sustainable or a non-sustainable way [34]. Then, the non-normalized density operator describes the system which inevitably decays (or diverges), whereas the normalized density operator indicates that gain and loss somehow balance each other out; therefore, the system sustains itself.

3. Generalized Schrödinger Equations

In this section, we derive relations between quantum-statistical master equations from the previous section and the quantum-mechanical Schrödinger-type equations for state vectors.

3.1. Dephasing and State-Vector Reduction

If one considers the evolution of only the state vectors, then one narrows the physical picture by neglecting various quantum-statistical phenomena, such as coherence, because they are predominantly described by the off-diagonal components of a density matrix. In other words, a forced dephasing can be performed, which turns a system into a plain statistical mixture, thus making its description closer to the classical. If this is acceptable, if an initial state can be prepared to be pure, and also if “mixing” quantum fluctuations can be neglected, then one can indeed resort to a state-vector description, i.e., consider, for each state, the time evolution of only one diagonal element of a density matrix. The rank of the latter thus effectively reduces to one because its other elements are assumed to be zeros.
We thus consider the special case of a statistical operator:
ρ ^ | Φ Φ | , ϱ ^ | Ω Ω | ,
for which Equation (5) reduces to a Hermitian conjugate pair of equations of a Schrödinger type:
i t | Φ = H ^ + i 2 k γ k L ^ k Φ Φ | Φ L ^ k L ^ k | Φ 1 2 H ^ H ^ Φ | Φ = H ^ + + i 4 k γ k Φ | Φ L ^ k Φ L ^ k L ^ k L ^ k Φ | Φ i Γ ^ Γ ^ Φ 1 4 k γ k L ^ k Φ L ^ k + L ^ k L ^ k Φ Φ | Φ 2 L ^ k L ^ k | Φ ,
and, similarly, for the non-normalized density matrix evolution (9):
i t | Ω = H ^ + i 2 k γ k L ^ k Ω Ω | Ω L ^ k L ^ k | Ω = H ^ + + i 4 k γ k Ω | Ω L ^ k Ω L ^ k L ^ k L ^ k Ω | Ω i Γ ^ + 1 4 k γ k 2 L ^ k L ^ k L ^ k Ω L ^ k + L ^ k L ^ k Ω Ω | Ω | Ω ,
where, in the second and subsequent lines of either Equation (12) or (13), we explicitly separate Hermitian and anti-Hermitian components of the Lindblad part:
H ^ L ( Φ ) = i 2 k γ k L ^ k Φ Φ | Φ L ^ k L ^ k , H ^ L ( Ω ) = i 2 k γ k L ^ k Ω Ω | Ω L ^ k L ^ k ,
H ^ L ± = 1 2 ( H ^ L ± H ^ L ) = ± H ^ L ± , Γ ^ L = i H ^ L ,
assuming our usual notations.
One can see that the Lindblad term can be now viewed as a part of the non-Hermitian Hamiltonian in Equation (12) or Equation (13), which was not possible in the underlying master equations (due to the presence of the term containing a density matrix “sandwiched” between two Lindblad operators). This is a side effect of the imposed dephasing procedure in (11), which turns some of the Lindblad operators into their averages.
Another side effect of the dephasing and rank-truncating of the original density matrices is that the norms of the corresponding state vectors | Φ and | Ω are not conserved in general (which is not surprising for | Ω but is for | Φ ). From Equations (12) and (13), one obtains the evolution equation for the norms:
t Φ | Φ = 2 Γ ^ Φ 1 Φ | Φ + k γ k L ^ k Φ L ^ k Φ Φ | Φ L ^ k L ^ k Φ ,
t Ω | Ω = 2 Γ ^ Ω + k γ k L ^ k Ω L ^ k Ω Ω | Ω L ^ k L ^ k Ω ,
which can be viewed as first-order differential equations for, respectively, Φ | Φ and Ω | Ω as functions of time. Each of these equations is in fact a generalization of the normalization condition of a state vector in terms of an inner product. They also indicate that the GKSL-type dissipation dynamically breaks the normalization of both wavefunctions Φ and Ω .
The norm (non-)conservation can be used in certain applications, for example, in establishing the classical limit of a dissipative system. According to Equation (16), the classical regime is necessary to stop the ultimate decay of the system because the Lindblad operators tend to become uncorrelated in the classical limit; therefore, the right-hand side of Equation (16) tends to zero, and Φ | Φ = 1 becomes a solution. This feature can be useful in decoherence models and in theories of mesoscopic systems and quantum chaos [35,36].
For other applications, it is more robust to use models with a state vector whose norm is conserved at all times, cf. the generalized Schrödinger equations proposed in [37,38,39]. These models also cover situations when the measuring process is restarted or calibrated to compensate for the “leakage” of probability.
Therefore, let us introduce the normalized state vector:
| Ψ Ω | Ω 1 / 2 | Ω ,
by analogy with the normalization procedure for density matrices (1). Its dynamics are governed by the following equation:
i t | Ψ = H ^ + i 2 k γ k L ^ k Ψ L ^ k L ^ k L ^ k Ψ L ^ k Ψ + L ^ k L ^ k Ψ | Ψ 1 2 H ^ H ^ Ψ | Ψ = H ^ + + i 4 k γ k L ^ k Ψ L ^ k L ^ k L ^ k Ψ | Ψ i Γ ^ Γ ^ Ψ 1 4 k γ k L ^ k Ψ L ^ k + L ^ k L ^ k Ψ 2 L ^ k L ^ k 2 L ^ k Ψ L ^ k Ψ + 2 L ^ k L ^ k Ψ | Ψ ,
where we remember the notation in (3). One can directly check that the norm of a state vector is conserved at all times: t Ψ | Ψ = 0 . Although the corresponding density matrix | Ψ Ψ | would no longer be a solution to either of the original master equations—(5) and (9)—this equation belongs to the same equivalence class as Equations (12) and (13).
Indeed, three state vectors are related through the transformations
| Ψ = Ω | Ω 1 / 2 | Ω = Φ | Φ 1 / 2 | Φ , | Ω = exp 1 Γ ^ Ψ + 1 2 k γ k L ^ k L ^ k Ψ L ^ k Ψ L ^ k Ψ d t | Ψ = exp 1 Γ ^ Φ d t | Φ ,
defined up to a standard unitary transformation. These relations mean that the three generalized Schrödinger equations we have derived, (12), (13), and (19), belong to the same equivalence class as defined above.
To summarize, by performing dephasing and rank reduction of the quantum-statistical equations, we derived three forms of the generalized Schrödinger equation: (12), (13), and (19). The former two were derived from, respectively, Equations (5) and (9), whereas the latter is the one with a conserved norm. These equations belong to the same equivalence class as (20); therefore, the choice between them for a specific model is a matter of physical relevance, cf. the remarks at the end of Section 2.

3.2. Special Cases of a Normalized State-Vector Equation

For reasons given in the previous section, we choose Equation (19) as the main object of our studies. For brevity, we omit the Lindblad operators’ summation and related indices from now on, i.e.,
k γ k L ^ k ρ ^ L ^ k 1 2 { ρ ^ , L ^ k L ^ k } γ L ^ ρ ^ L ^ 1 2 { ρ ^ , L ^ L ^ } ,
while implying the summation whenever the Lindblad operators occur in a quadratic combination.
Since the NH and GKSL parts of Equation (19) have quite different physical origins and applications, it is convenient to consider them separately:
  • Generic NH-driven evolution. In this case, we neglect the Lindblad part; therefore, Equation (19), as well as Equations (12) and (13), reduce to
    i t | Ψ = H ^ | Ψ 1 2 H ^ H ^ Ψ | Ψ = H ^ + | Ψ i Γ ^ Γ ^ Ψ | Ψ ,
    where the anti-Hermitian part H ^ = i Γ ^ is assumed to account for all dissipative effects except those coming from the Lindblad’s jump operators;
  • GKSL-driven evolution. In this case, we neglect the non-Lindblad dissipation part; therefore, Equation (19) reduces to
    i t | Ψ = H ^ + + i 2 γ L ^ Ψ L ^ L ^ + L ^ L ^ Ψ L ^ Ψ L ^ Ψ | Ψ = H ^ + + H ^ L + | Ψ i ( Γ ^ L Γ ^ L ) | Ψ , H ^ L + i 4 γ L ^ Ψ L ^ L ^ Ψ L ^ , Γ ^ L 1 4 γ L ^ Ψ L ^ + L ^ Ψ L ^ 2 L ^ L ^ ,
    where H ^ + is a Hermitian part of the original subsystem’s Hamiltonian, and we assume the notations in (21).
    This equation confirms that the GKSL-driven evolution of a state vector belongs to a class of non-Hermitian Hamiltonian evolutions with a conserved norm; see also the remark after Equation (15).
In the next section, we study these special cases separately, while keeping OQS applications in mind.

4. Basic Properties

Let us study the main properties of the generalized Schrödinger equation in (19) while focusing on its limit cases, as discussed in Section 3.2 above.

4.1. Generic NH-Driven Evolution

In this special case, evolution is governed by Equation (22). It conserves the norm of the state vector because
t Ψ | Ψ = 0
at all times. The evolution equation is invariant under a complex-valued shift of the zero of the operator H ^ :
H ^ H ^ + ( ε 0 i γ 0 ) I ^ , | Ψ e i ε 0 t / | Ψ ,
with I ^ being the identity operator, as per usual. Remember that the conventional Schrödinger equation for conservative systems is invariant only under a real-valued shift in the zero of the Hamiltonian operator; therefore, the invariance in (25) is something which is more general. This seems plausible from the phenomenological point of view because observables should be independent not only from the real-valued shifts in the energy’s zero but also from the imaginary-valued shifts that could possibly occur due to quantum fluctuations and virtual particles.
Norm conservation and shift invariance are the features that simultaneously occur in the master equation in (5).
Furthermore, if | Ψ is an eigenvector of the decay rate operator Γ ^ , then our equation reduces to the ordinary Schrödinger equation with the Hermitian Hamiltonian H ^ + .
If H ^ is time-independent, then mean values of its components evolve according to the following formulae:
d d t H ^ + Ψ = 2 H ^ + Ψ Γ ^ Ψ 1 2 { H ^ + , Γ ^ } Ψ ,
d d t Γ ^ Ψ = 2 Γ ^ Ψ 2 Γ ^ 2 Ψ + i 2 [ H ^ + , Γ ^ ] Ψ = i [ H ^ + , Γ ^ ] Ψ 2 Var ( Γ ) ,
where Var ( X ) X 2 Ψ X Ψ 2 is a variance of a physical variable described by the operator X ^ . Therefore, if | Ψ is an eigenvector of H ^ + , then the average energy is conserved, and, if it is an eigenvector of the decay operator Γ ^ , then the mean values of both energy and decay rate are conserved.
Furthermore, the evolution Equation (22) allows the usual stationary states, i.e, the eigenvectors of the Hamiltonian H ^ + , only if the latter commutes with the decay operator (which corresponds to H ^ being a normal operator).
This means that our model, even in its truncated form (without the Lindblad part), can describe two types of dissipative systems: those with balanced gain and loss, for which H ^ + , Γ ^ = ( i / 2 ) H ^ , H ^ = 0 (hence, the usual stationary states can exist [13]), and those with strong dissipation, for which H ^ + , Γ ^ 0 (hence, no usual stationary states can exist in general). This feature significantly expands the range of applicability of our generalized Schrödinger equation approach compared to some other generalizations, cf. [37], and places the proverbial pyramid on its base.

4.2. GKSL-Driven Evolution of L-Eigenvector

Let us consider the special case of purely GKSL-driven evolution in (23), with two additional assumptions: we restrict ourselves to only one species of Lindblad jump operators (hence, the index k takes only one value), and we assume that our state vector is an eigenvector of the jump operator:
| Ψ = | L : L ^ | L = L | L , L | L = 1 ,
where an eigenvalue L = L | L ^ | L is a complex-valued function of time in general.
Under these assumptions, Equation (23) yields the following evolution equation for the state vector | L :
i t | L = H ^ + i 2 γ 1 L ^ L * I ^ L ^ | L = H ^ + + H ^ L + | Ψ i ( Γ ^ L Γ ^ L L ) | Ψ , H ^ L + i 4 γ 1 L * L ^ L L ^ i 4 γ 1 | L | 2 L L ^ , Γ ^ L 1 4 γ 1 L * L ^ + L L ^ 2 L ^ L ^ 1 4 γ 1 | L | 2 L L ^ ,
where the notation ‘↦’ indicates that the projection of the jump operator on the eigen(ket)vector is implied. Notice also that, for this model,
H ^ L + L = 0 , Γ ^ L L = 0 ,
which can be directly verified.
Otherwise, this model belongs to a class of NH models discussed in Section 4.1, as one can clearly see from the form of its evolution equation. This means that the properties in (24) and (25), as well as other generic properties of NH-driven evolution equations, would be valid for this case too.
If operators H ^ + and L ^ are themselves time-independent, then their mean values evolve according to the following equations:
d d t H ^ + L = γ 1 2 | L | 2 2 H ^ + L 1 L L ^ H ^ + L 1 L * H ^ + L ^ L = γ 1 2 | L | 2 C ^ L + + C ^ L + L = γ 1 2 | L | 2 C ^ Δ + L ,
1 L d d t L = i 1 L L ^ H ^ + L H ^ + L + γ 1 2 | L | 2 L ^ L ^ L = i C ^ L + L γ 1 2 C ^ L L ,
1 L * d d t L * = i 1 L * H ^ + L ^ L H ^ + L + γ 1 2 | L | 2 L ^ L ^ L = i C ^ L + L γ 1 2 C ^ L L ,
where we denoted C ^ L L ^ , L ^ = C ^ L , C ^ L + H ^ + , L ^ / L , C ^ Δ + H ^ + , Δ ^ L and Δ ^ L L ^ / L L ^ / L * = Δ ^ L .
From the first of these equations, one can deduce that the original subsystem’s average energy is conserved if operators H ^ + and i Δ ^ L commute, i.e., they have common eigenvectors.
Furthermore, by rearranging those equations one can obtain the relation
d d t H ^ + L = γ 1 | L | 2 H ^ + L i γ 1 2 1 L d L d t 1 L * d L * d t ,
which can be regarded as a first-order differential equation for H ^ + L as a function of time. Solving it, we obtain the relation between mean values
H ^ + L = H ^ + 0 e γ 1 t 0 t | L | 2 d t i γ 1 2 ln L / L * i γ 1 2 2 e γ 1 t 0 t | L | 2 d t t 0 t e γ 1 t 0 t | L | 2 d t ln L / L * | L | 2 d t ,
where H ^ + 0 = H ^ + L | t = t 0 is an initial value of average energy as a function of time.
Furthermore, let us study the eigenvectors of the Hamiltonian H ^ + , which are solutions to the evolution equation (Equation (29)). For simplicity, let us consider a special case when H ^ + commutes with the Lindblad operator:
H ^ + , L ^ = 0 ,
thus, having common eigenvectors, H ^ + | L ε = ε | L ε . In this case, mean value of H ^ + is conserved in time, which corresponds to an open system with balanced gain and loss.
Let us also perform the decomposition into the projected and orthogonal parts with respect to the vector | L ε :
L ^ | L ε = L * | L ε   +   L | ,
where L ε | = 0 and L = | L ^ | L ε / | . Then, Equation (29) reduces to the differential equation
t | L ε = i ε | L ε     γ 1 2 L L | ,
whose formal solution can be written in the form
| L ε = e i ε ( t t 0 ) / | L 0 γ 1 2 t 0 t d t e i ε t / L L | ,
where | L 0 = | L ε | t = t 0 is an initial value of | L ε as a vector function of time.
Equation (38) and the solution in (39) demonstrate that the eigenstates of the Hamiltonian H ^ + , which are solutions of Equation (29), are not stationary in general, even in the balanced case with conserved average energy. This reaffirms that the evolution in (29) is indeed of a dissipative type, with the time-dependent source/sink vector L L | .

5. Example: Dissipative Harmonic Oscillator

In the conventional quantum mechanics of conservative systems, the harmonic oscillator is a benchmark model with a clearly derivable classical limit and plenitude of applications. Its Hamiltonian can be written in the standard form:
H ^ HO = p ^ 2 / 2 + m ω 2 q ^ 2 / 2 = ω ( a ^ a ^ + 1 / 2 ) ,
where a ^ = ( m ω / q ^ + i p ^ / m ω ) / 2 and a ^ are ladder operators whose commutator [ a ^ , a ^ ] = 1 , as per usual.
By taking quantum averages in the eigenbasis of the operator a ^ , usually referred to as coherent states, one can obtain a set of ordinary differential equations:
d d t α = i ω α ,
d d t p ^ α q ^ α = 0 m ω 2 1 / m 0 p ^ α q ^ α ,
where a ^ | α = α | α . From these equations, a second-order differential equation immediately follows:
d 2 q d t 2 + ω 2 q = 0 ,
where q = q ^ α , which describes the classical harmonic oscillator.
On the other hand, let us recall the classical damped oscillator
d 2 q d t 2 + 2 δ d q d t + ω 0 2 q = 0 ,
where δ = μ / 2 m ˜ , ω 0 = κ / m ˜ , and m ˜ , κ , and μ are, respectively, mass, spring constant, and friction. This model was extensively used in the early quantization attempts of dissipative systems.
The straightforward quantization of the system in (44) via the Lagrangian approach has been shown to be somewhat inconsistent [1,40]. The currently prevailing point of view is that the damped oscillator can be a classical limit of a number of quantum models; one example is the model based on Gisin’s generalized Schrödinger equation.
In our case, the question arises: what model, as defined by means of both Hamiltonian and non-Hamiltonian (Liouvillian) terms, recovers Equation (44) in the classical limit?
Similarly to the previous section, let us consider cases of generic NH-driven and GKSL-driven evolutions separately first; then, we consider the unified NH–GKSL evolution.

5.1. Generic NH-Driven Evolution

Let us generalize the Hamiltonian (40) by moving its overall coupling constant into the complex domain. From the quantum-mechanical point of view, it is a natural generalization because coupling constants are not observables per se and, therefore, do not have to be real numbers, especially in the presence of quantum fluctuations and virtual particles. We thus assume that
H ^ = ( 1 i η / ω ) H ^ HO = ( ω i η ) ( a ^ a ^ + 1 / 2 ) ,
while the notations in Equation (40) remain intact. Correspondingly, the evolution of this system is described by Equation (22), with the following associations:
H ^ + = ω ( a ^ a ^ + 1 / 2 ) , Γ ^ = η a ^ a ^ ,
while mean values are assumed to be computed in the eigenbasis | α of the operator a ^ . Note that solutions to other generalized Schrödinger Equations (12) and (13), in this basis, are related to | Ψ via the transformations
exp η | α | 2 d t | α ( Ω ) = | α ( Ψ ) = | α ( Φ ) ,
which can be easily derived from Equation (20); it is, thus, sufficient to study a case of a normalized state vector | Ψ , as mentioned before.
Upon averaging in the basis | α , we obtain the analogues of Equations (41) and (42) for the model in (45):
d d t α = ( η + i ω ) α ,
d d t p ^ α q ^ α = η m ω 2 1 / m η p ^ α q ^ α ,
from which a second-order differential equation immediately follows:
d 2 q d t 2 + 2 η d q d t + ( ω 2 + η 2 ) q = 0 ,
where q = q ^ α . Thus, the classical limit of the model (45) is the damped oscillator (44), with
δ η , ω 0 2 ω 2 + η 2 ,
which means that the NH coupling induces friction in the classical limit. Solving this, we reveal the following cases:
  • If ω 0 and η > 0 , then the damping is subcritical, which is an exponential decay with oscillations with the frequency ω :
    q ( t ) = q 0 exp ( η t ) cos ( ω t ) + η + q ˙ 0 / q 0 ω sin ( ω t ) ,
    where q 0 and q ˙ 0 are, respectively, oscillator’s position and velocity at t = 0 ;
  • If ω = 0 and η > 0 , then the damping is critical:
    q ( t ) = lim ω 0 q 0 exp ( η t ) cos ( ω t ) + η + q ˙ 0 / q 0 ω sin ( ω t ) = q 0 exp ( η t ) 1 + ( η + q ˙ 0 / q 0 ) t ,
    which is rather surprising because the “Hamiltonian” (45) is completely skew-adjoint in this case, yet it still yields a realistic classical system upon averaging.
Thus, in the classical limit, the model described by Equations (45) and (22) reduces to either an underdamped or a critically damped harmonic oscillator, depending on the value of frequency.

5.2. GKSL-Driven Evolution

Let us consider the conventional harmonic oscillator in (40) and assume that its evolution is affected by the dissipative effects described by two jump operators:
L ^ 1 = a ^ , L ^ 2 = a ^ ,
with the non-negative relaxation rates γ 1 and γ 2 , respectively. Then, Equation (23) takes the form
i t | Ψ = H ^ HO | Ψ + i 2 γ 1 ( a ^ Ψ a ^ ) a ^ + a ^ a ^ Ψ a ^ Ψ a ^ Ψ + γ 2 a ^ a ^ | Ψ = H ^ HO + H ^ L + | Ψ i ( Γ ^ L Γ ^ L ) | Ψ , H ^ L + i 4 γ 1 a ^ Ψ a ^ a ^ Ψ a ^ + γ 2 a ^ Ψ a ^ a ^ Ψ a ^ , Γ ^ L 1 4 γ 1 a ^ Ψ a ^ + a ^ Ψ a ^ 2 a ^ a ^ + γ 2 a ^ Ψ a ^ + a ^ Ψ a ^ 2 a ^ a ^ ,
where H ^ HO is borrowed from Equation (40), and a shorthand notation ↔ means an interchange of symbols.
In the eigenbasis of the operator a ^ , solutions to other generalized Schrödinger Equations (12) and (13) are related via the transformations
| α ( Ω ) = exp 1 2 γ 2 t | α ( Ψ ) = | α ( Φ ) ,
which can be easily derived from Equation (20); it is, thus, sufficient to study a case of a normalized state vector | Ψ only.
Upon averaging in the basis | α , Equation (55) yields the following analogues of Equations (41) and (42):
d d t α = 1 2 γ 1 + i ω α ,
d d t p ^ α q ^ α = γ 1 / 2 m ω 2 1 / m γ 1 / 2 p ^ α q ^ α ,
from which a second-order differential equation immediately follows:
d 2 q d t 2 + γ 1 d q d t + ω 2 + 1 4 γ 1 2 q = 0 ,
where q = q ^ α , as per usual.
Thus, the classical limit of the model in (55) is the damped oscillator (44) with
δ 1 2 γ 1 , ω 0 2 ω 2 + 1 4 γ 1 2 ,
which means that the Lindblad evolution also induces friction in the classical limit.
It is worth noticing that the γ 2 -coupled term, γ 2 ( a ^ Ψ a ^ ) a ^ + a ^ a ^ Ψ a ^ Ψ a ^ Ψ , does not affect the evolution of coherent-state mean values based on Equation (19). This indicates that this term represents higher-order corrections that can be taken into account when returning to the more general (density operator) formalism.
Furthermore, by solving Equation (59), we reveal the following special cases:
  • If ω 0 , then the damping is subcritical, which is an exponential decay with oscillations with the frequency ω :
    q ( t ) = q 0 exp ( γ 1 t / 2 ) cos ( ω t ) + γ 1 / 2 + q ˙ 0 / q 0 ω sin ( ω t ) ,
    where q 0 and q ˙ 0 are, respectively, oscillator’s position and velocity at t = 0 ;
  • If ω 0 , then the damping turns critical:
    q ( t ) = lim ω 0 q 0 exp ( γ 1 t / 2 ) cos ( ω t ) + γ 1 / 2 + q ˙ 0 / q 0 ω sin ( ω t ) = q 0 exp ( γ 1 t / 2 ) 1 + ( γ 1 / 2 + q ˙ 0 / q 0 ) t ,
    which is also surprising because the Hamiltonian can be neglected. In this case, system’s evolution is driven solely by the non-Hamiltonian (Liouvillian) term, yet we obtain a physically relevant classical system upon averaging.

5.3. NH–GKSL Evolution and Cancellation of Dissipations

Now let us unify the models described in Section 5.1 and Section 5.2: we assume that the dissipative effects are encoded both in the non-Hermitian Hamiltonian (45) and in the Lindblad term with two jump operators (54).
In this case, Equation (19) takes the form
i t | Ψ = 1 i η ω H ^ HO | Ψ + i η ω H ^ HO Ψ | Ψ + i 2 γ 1 ( a ^ Ψ a ^ ) a ^ + a ^ a ^ Ψ a ^ Ψ a ^ Ψ + γ 2 a ^ a ^ | Ψ = H ^ HO + H ^ L + | Ψ i Γ ^ L Γ ^ L Ψ + Γ ^ η Γ ^ η Ψ | Ψ , H ^ L + i 4 γ 1 a ^ Ψ a ^ a ^ Ψ a ^ + γ 2 a ^ Ψ a ^ a ^ Ψ a ^ , Γ ^ L 1 4 γ 1 a ^ Ψ a ^ + a ^ Ψ a ^ 2 a ^ a ^ + γ 2 a ^ Ψ a ^ + a ^ Ψ a ^ 2 a ^ a ^ , Γ ^ η η a ^ a ^ ,
where H ^ HO has been defined in Equation (40), and shorthand notation ‘↔’ means an interchange of symbols.
In the eigenbasis of the operator a ^ , solutions to the other generalized Schrödinger Equations (12) and (13) are related via the transformations in (20):
exp η | α | 2 d t | α ( Ω ) = exp 1 2 γ 2 t | α ( Ψ ) = | α ( Φ ) ,
which supports our decision to study the normalized state-vector case (see the end of Section 3.1).
Furthermore, upon averaging in the basis | α , Equation (63) yields the following analogues of Equations (41) and (42):
d d t α = 1 2 γ 1 + η + i ω α ,
d d t p ^ α q ^ α = η γ 1 / 2 m ω 2 1 / m η γ 1 / 2 p ^ α q ^ α ,
from which a second-order differential equation immediately follows:
d 2 q d t 2 + 2 ( η + γ 1 / 2 ) d q d t + ω 2 + ( η + γ 1 / 2 ) 2 q = 0 ,
where q = q ^ α , as per usual.
Thus, the classical limit of the model in (55) is the damped oscillator in (44) with
δ η + γ 1 / 2 , ω 0 2 ω 2 + ( η + γ 1 / 2 ) 2 ,
which means that the NH–GKSL evolution also induces friction in the classical limit. Its solutions reveal the following special cases:
  • If ω 0 and η = γ 1 / 2 (the latter is possible because η can be both positive and negative, unlike the strictly non-negative Lindblad relaxation constants), then the damping term disappears in Equation (67). In this case, evolution becomes a plain oscillation with the frequency ω :
    q ( t ) = q 0 cos ( ω t ) + q ˙ 0 q 0 ω sin ( ω t ) ,
    where q 0 and q ˙ 0 are, respectively, the oscillator’s position and velocity at t = 0 . This indicates that the dissipative effects associated with NH and GKSL terms can act against each other, which leads to their mutual cancellation;
  • If ω 0 and η γ 1 / 2 , then the damping is subcritical, which is an exponential decay with oscillations at the frequency ω :
    q ( t ) = q 0 exp ( ( η + γ 1 / 2 ) t ) cos ( ω t ) + η + γ 1 / 2 + q ˙ 0 / q 0 ω sin ( ω t ) ,
    where q 0 and q ˙ 0 are, respectively, oscillator’s position and velocity at t = 0 ;
  • If ω 0 and η + γ 1 / 2 > 0 , then the damping becomes critical:
    q ( t ) = lim ω 0 q 0 exp ( ( η + γ 1 / 2 ) t ) cos ( ω t ) + η + γ 1 / 2 + q ˙ 0 / q 0 ω sin ( ω t ) = q 0 exp ( ( η + γ 1 / 2 ) t ) 1 + ( η + γ 1 / 2 + q ˙ 0 / q 0 ) t ,
    which is somewhat surprising because the Hamiltonian is neglected in this case; therefore, the evolution becomes dominated by the Liouvillian term, yet we obtain a physically realistic classical system upon averaging;
  • If ω 0 and η + γ 1 / 2 = 0 , then the damping disappears, while the average value of q ^ α grows linearly in time:
    q ( t ) = lim ω 0 q 0 cos ( ω t ) + q ˙ 0 q 0 ω sin ( ω t ) = q 0 + q ˙ 0 t ,
    which corresponds to a system in free motion. Notice that this classically free system descends from a quantum system that does not have a Hamiltonian, hence it experiences no forces in a sense of Newton’s first law of motion, yet its underlying dynamics are non-trivial (due to the presence of non-Hamiltonian terms).

6. Conclusions

In the present study, we started out with the quantum-statistical approach, using both the non-Hermitian Hamiltonian and Lindblad’s jump operators: we presented various generalizations of the von Neumann equation for density operator, which did not assume the Markovian approximation in a general case. We demonstrated that, if one performs the dephasing procedure by neglecting off-diagonal components of the density matrix that are responsible for coherence and spontaneous transitions from pure to mixed states, followed by neglecting all but one diagonal component thereof, then one reduces quantum-statistical equations to quantum-mechanical Schrödinger-type equations for state vectors. We thus restricted ourselves to dealing with pure states only, but it can be a robust approximation in a number of cases where one has the full classical statistical information about a system (which does not imply the full knowledge of the possible outcomes of quantum measurement).
We first obtained generalized Schrödinger Equations (12) and (13), which directly originated from the master equations, i.e., Equations (5) and (9), for normalized and non-normalized density operators, respectively. We showed that both evolutions, (12) and (13), do not preserve the norm of a state vector, due to the fact that the dephasing and rank reduction procedures replace some jump operators with their mean values. Therefore, in addition to those equations, we derived the generalized Schrödinger Equation (19), which directly descends from neither of the original master equations, but describes the time evolution of the state vector with a conserved norm. In a sense, Equation (19) can be viewed as a state-vector analogue of Equation (5) because it is also derived using the post-selecting procedure of normalization.
Furthermore, we demonstrated that all three Equations, (12), (13), and (19), belong to the same equivalence class. For this reason, we focused our studies on the evolution in (19), while keeping in mind applications to open quantum systems. We demonstrated that equations from the equivalence class in (20) can describe not only systems which can stay in the stationary eigenstates of the Hamiltonian, but also those which can evolve from those eigenstates. This indicates that our generalizations have a broader range of applicability than some of the generalized Schrödinger equations proposed hitherto. Another advantage is that they do not contain any stochastic terms or other random functions.
From a mathematical viewpoint, each of these equations in general is both nonlinear (with respect to wavefunction) and non-local because various mean values inside their structure contain wavefunction and integration over the whole configuration space. These two features indicate that our equations form a modeling framework that not only has a rich mathematical structure, but is also a priori more natural for describing diverse quantum phenomena and systems, including those involving measuring apparati.
Finally, as an instructive example, we considered the model of the quantum harmonic oscillator in the presence of dissipative effects described by non-Hermitian Hamiltonian and Lindblad’s terms. Using coherent states’ basis, we derived its classical limit, which turned out to be the damped harmonic oscillator, though with different regimes. In particular, we considered the case when non-Hermitian Hamiltonian and Lindblad’s dissipative effects cancel each other out, thus resulting in the decay-free evolution of an oscillating system, at least in the leading-order approximation after dephasing and averaging.
Yet another interesting example shown is the classically free system which descends from a non-trivial (not free from interactions) quantum system. A seeming contradiction to Newton’s first law of motion is resolved due to the fact that the quantum system in question does not have a Hamiltonian, therefore no Newtonian forces are present therein; instead, its quantum dynamics are governed by NH and GKSL terms, which are both non-Hamiltonian by nature. This uncovers a way to formulate a large class of quantum-mechanical non-Hamiltonian models whose dynamics are not determined by conventional potentials and forces, but which come about through the quantum statistical effects of the system’s environment.

Funding

This work is based on research supported by the Department of Higher Education and Training of South Africa and in part by the National Research Foundation of South Africa (Grants Nos. 95965 and 132202).

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Acknowledgments

The author gratefully acknowledges an invitation and a full waiver on article publication charges by Multidisciplinary Digital Publishing Institute. Proofreading of the manuscript by P. Stannard is greatly appreciated.

Conflicts of Interest

The author declares no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:
GKSLGorini–Kossakowski–Sudarshan–Lindblad
HOHarmonic oscillator (quantum)
NHNon-Hermitian Hamiltonian
OQSOpen quantum system

References

  1. Hasse, R.W. On the Quantum Mechanical Treatment of Dissipative Systems. J. Math. Phys. 1975, 16, 2005–2011. [Google Scholar] [CrossRef]
  2. Faisal, F.H.M. Theory of Multiphoton Processes; Plenum Press: New York, NY, USA, 1987; p. 408. [Google Scholar]
  3. Breuer, H.-P.; Petruccione, F. The Theory of Open Quantum Systems; Oxford University Press: Oxford, UK, 2002; p. 636. [Google Scholar]
  4. 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]
  5. Lindblad, G. On the Generator of Quantum Dynamical Semigroups. Commun. Math. Phys. 1976, 48, 119–130. [Google Scholar] [CrossRef]
  6. Feshbach, H. Unified Theory of Nuclear Reactions. Ann. Phys. 1958, 5, 357–390. [Google Scholar] [CrossRef]
  7. Feshbach, H. A Unified Theory of Nuclear Reactions. II. Ann. Phys. 1962, 19, 287–313. [Google Scholar] [CrossRef]
  8. Sergi, A.; Zloshchastiev, K.G. Non-Hermitian Quantum Dynamics of a Two-level System and Models of Dissipative Environments. Int. J. Mod. Phys. B 2013, 27, 1350163. [Google Scholar] [CrossRef]
  9. Zloshchastiev, K.G.; Sergi, A. Comparison and Unification of Non-Hermitian and Lindblad Approaches with Applications to Open Quantum Optical Systems. J. Mod. Optics 2014, 61, 1298. [Google Scholar] [CrossRef]
  10. Zloshchastiev, K.G. Non-Hermitian Hamiltonians and Stability of Pure States in Quantum Mechanics. Eur. Phys. J. D 2015, 69, 253. [Google Scholar] [CrossRef]
  11. Liao, Q.-H.; Deng, W.-C.; Wen, J.; Zhou, N.-R.; Liu, N.-H. Phonon Blockade Induced by a Non-Hermitian Hamiltonian in a Nanomechanical Resonator Coupled with a Qubit. Acta Phys. Sin. 2019, 68, 114203. [Google Scholar] [CrossRef]
  12. Zloshchastiev, K.G. Density Operator Approach to Turbulent Flows in Plasma and Atmospheric Fluids. Universe 2020, 6, 216. [Google Scholar] [CrossRef]
  13. Zloshchastiev, K.G. Model Hamiltonians of Open Quantum Optical Systems: Evolvement from Hermiticity to Adjoint Commutativity. J. Phys. Conf. Ser. 2022, 2407, 012011. [Google Scholar] [CrossRef]
  14. Parkavi, J.R.; Muthuganesan, R.; Chandrasekar, V.K. Temporal Evolution of Quantum Correlations under Non-Hermitian Operation. Opt. Quant. Electron. 2022, 54, 729. [Google Scholar] [CrossRef]
  15. Mosquera, M.A.; Marmolejo-Tejada, J.M.; Borys, N.J. Theoretical Quantum Model of Two-dimensional Propagating Plexcitons. J. Chem. Phys. 2022, 157, 124103. [Google Scholar] [CrossRef]
  16. Li, D.; Zheng, C. Non-Hermitian Generalization of Rényi Entropy. Entropy 2022, 24, 1563. [Google Scholar] [CrossRef] [PubMed]
  17. Rosso, L.; Biella, A.; De Nardis, J.; Mazza, L. Dynamical Theory for One-dimensional Fermions with Strong Two-body Losses: Universal Non-Hermitian Zeno Physics and Spin-charge Separation. Phys. Rev. A 2023, 107, 013303. [Google Scholar] [CrossRef]
  18. Grimaudo, R.; Valenti, D.; Sergi, A.; Messina, A. Superradiant Quantum Phase Transition for an Exactly Solvable Two-qubit Spin-boson Model. Entropy 2023, 25, 187. [Google Scholar] [CrossRef]
  19. Itable, G.M.M.; Paraan, F.N.C. Entanglement Entropy Distinguishes PT-symmetry and Topological Phases in a Class of Non-unitary Quantum Walks. Quantum Inf. Process. 2023, 22, 106. [Google Scholar] [CrossRef]
  20. Pires, D.P.; Macrì, T. Mixedness Timescale in Non-Hermitian Quantum Systems. Phys. Rev. A 2023, 107, 022202. [Google Scholar] [CrossRef]
  21. Chakrabarty, A.; Datta, S. Skin Effect and Dynamical Delocalization in Non-Hermitian Quasicrystals with Spin-orbit Interaction. Phys. Rev. B 2023, 107, 064305. [Google Scholar] [CrossRef]
  22. Geller, M.R. Fast Quantum State Discrimination with Nonlinear Positive Trace-preserving Channels. Adv. Quantum Technol. 2023, 2023, 2200156. [Google Scholar] [CrossRef]
  23. Sun, J.Y.; Shen, H.Z. Photon Blockade in Non-Hermitian Optomechanical Systems with Nonreciprocal Couplings. Phys. Rev. A 2023, 107, 043715. [Google Scholar] [CrossRef]
  24. Luan, T.Z.; Sun, J.Y.; Shen, H.Z. Dynamical Approach to Shortcuts to Adiabaticity for General Two-level Non-Hermitian Systems. Europhys. Lett. 2023, 142, 58001. [Google Scholar] [CrossRef]
  25. Zhao, W.-L.; Wang, R.-R.; Ke, H.; Liu, J. Scaling Laws of the Out-of-time-order Correlators at the Transition to the Spontaneous PT-symmetry Breaking in a Floquet System. Phys. Rev. A 2023, 107, 062201. [Google Scholar] [CrossRef]
  26. Cheng, D.; Lustig, E.; Wang, K.; Fan, S. Multi-dimensional Band Structure Spectroscopy in the Synthetic Frequency Dimension. Light Sci. Appl. 2023, 12, 158. [Google Scholar] [CrossRef] [PubMed]
  27. Lima, L.S. Spin Transport in Non-Hermitian Quantum Systems. Sci. Rep. 2023, 13, 11112. [Google Scholar] [CrossRef]
  28. Saleem, Z.H.; Shaji, A.; Gray, S.K. Optimal Time for Sensing in Open Quantum Systems. Phys. Rev. A 2023, 108, 022413. [Google Scholar] [CrossRef]
  29. Niu, X.; Li, J.; Wu, S.L.; Yi, X.X. Effect of Quantum Jumps on Non-Hermitian Systems. Phys. Rev. A 2023, 108, 032214. [Google Scholar] [CrossRef]
  30. Yuan, Z.-H.; Chen, Y.-J.; Han, J.-X.; Wu, J.-L.; Li, W.-Q.; Xia, Y.; Jiang, Y.-Y.; Song, J. Periodic Photon-magnon Blockade in an Optomagnonic System with Chiral Exceptional Points. Phys. Rev. B 2023, 108, 134409. [Google Scholar] [CrossRef]
  31. Leng, Y.; Chen, F. Improving the Capacity of Quantum Dense Coding via Non-Hermitian Operation. Laser Phys. 2023, 33, 125203. [Google Scholar] [CrossRef]
  32. Meden, V.; Grunwald, L.; Kennes, D.M. P T -symmetric, Non-Hermitian Quantum Many-body Physics—A Methodological Perspective. Rep. Prog. Phys. 2023, 86, 124501. [Google Scholar] [CrossRef]
  33. Singh, G.; Singh, R.K.; Banerjee, S. Embedding of a Non-Hermitian Hamiltonian to Emulate the Von Neumann Measurement Scheme. J. Phys. A Math. Theor. 2024, 57, 035301. [Google Scholar] [CrossRef]
  34. Zloshchastiev, K.G. Sustainability of Environment-assisted Energy Transfer in Quantum Photobiological Complexes. Ann. Phys. 2017, 529, 1600185. [Google Scholar] [CrossRef]
  35. Imry, Y. Introduction to Mesoscopic Physics; Oxford University Press: New York, NY, USA, 2002; p. 236. [Google Scholar]
  36. Nakamura, K. Quantum versus Chaos: Questions Emerging from Mesoscopic Cosmos; Kluwer: Dordrecht, The Netherlands, 2002; p. 214. [Google Scholar]
  37. Gisin, N. A Simple Nonlinear Dissipative Quantum Evolution Equation. J. Phys. A 1981, 14, 2259–2267. [Google Scholar] [CrossRef]
  38. Dalibard, J.; Castin, Y.; Mølmer, K. Wave-Function Approach to Dissipative Processes in Quantum Optics. Phys. Rev. Lett. 1992, 68, 580–583. [Google Scholar] [CrossRef]
  39. Gisin, N.; Percival, I.C. The Quantum-state Diffusion Model Applied to Open Systems. J. Phys. A 1992, 25, 5677–5691. [Google Scholar] [CrossRef]
  40. Ray, J.R. Lagrangians and Systems They Describe—How not to Treat Dissipation in Quantum Mechanics. Am. J. Phys. A 1979, 47, 626–629. [Google Scholar] [CrossRef]
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.

Share and Cite

MDPI and ACS Style

Zloshchastiev, K.G. Generalization of the Schrödinger Equation for Open Systems Based on the Quantum-Statistical Approach. Universe 2024, 10, 36. https://doi.org/10.3390/universe10010036

AMA Style

Zloshchastiev KG. Generalization of the Schrödinger Equation for Open Systems Based on the Quantum-Statistical Approach. Universe. 2024; 10(1):36. https://doi.org/10.3390/universe10010036

Chicago/Turabian Style

Zloshchastiev, Konstantin G. 2024. "Generalization of the Schrödinger Equation for Open Systems Based on the Quantum-Statistical Approach" Universe 10, no. 1: 36. https://doi.org/10.3390/universe10010036

APA Style

Zloshchastiev, K. G. (2024). Generalization of the Schrödinger Equation for Open Systems Based on the Quantum-Statistical Approach. Universe, 10(1), 36. https://doi.org/10.3390/universe10010036

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop