Next Article in Journal
Market Choices Driven by Reference Groups: A Comparison of Analytical and Simulation Results on Random Networks
Next Article in Special Issue
Para-Hamiltonian form for General Autonomous ODE Systems: Introductory Results
Previous Article in Journal
Noise-Robust Image Reconstruction Based on Minimizing Extended Class of Power-Divergence Measures
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

General Non-Markovian Quantum Dynamics

1
Skobeltsyn Institute of Nuclear Physics, Lomonosov Moscow State University, 119991 Moscow, Russia
2
Faculty Information Technologies and Applied Mathematics, Moscow Aviation Institute, National Research University, 125993 Moscow, Russia
Submission received: 24 June 2021 / Revised: 27 July 2021 / Accepted: 29 July 2021 / Published: 31 July 2021
(This article belongs to the Special Issue Non-Hamiltonian Dynamics, Open Systems and Entropy)

Abstract

:
A general approach to the construction of non-Markovian quantum theory is proposed. Non-Markovian equations for quantum observables and states are suggested by using general fractional calculus. In the proposed approach, the non-locality in time is represented by operator kernels of the Sonin type. A wide class of the exactly solvable models of non-Markovian quantum dynamics is suggested. These models describe open (non-Hamiltonian) quantum systems with general form of nonlocality in time. To describe these systems, the Lindblad equations for quantum observable and states are generalized by taking into account a general form of nonlocality. The non-Markovian quantum dynamics is described by using integro-differential equations with general fractional derivatives and integrals with respect to time. The exact solutions of these equations are derived by using the operational calculus that is proposed by Yu. Luchko for general fractional differential equations. Properties of bi-positivity, complete positivity, dissipativity, and generalized dissipativity in general non-Markovian quantum dynamics are discussed. Examples of a quantum oscillator and two-level quantum system with a general form of nonlocality in time are suggested.

1. Introduction

The dynamics of open quantum systems (OQS) is usually described by Markovian equation for quantum observables and quantum states [1,2,3,4,5]. The most general explicit form of such equations were first proposed by Gorini, Kossakowski, Sudarshan, and Lindblad in 1976 [6,7,8]. Quantum mechanics of OQS has been actively developing (see [1,2,3,4,5,9, and reviews [10,11,12]). The Markovian dynamics of OQS is characterized by non-standard properties that lead to a change in the usual relations and approaches, for example, such as the uncertainty relation for open quantum systems [12,13], path integrals [5] (pp. 475–485) and [14], pure stationary states [5] (pp. 453–462) and [15,16], and quantum computer with mixed states [5] (pp. 487–520) and [17]. We can say that the theory of open quantum systems is the most general type of modern quantum mechanics such as fundamental theory [5]. Note that OQS are non-Hamiltonian quantum systems, that is, it is not enough to specify the Hamilton operator to describe the dynamics of such systems [5]. This theory has great practical importance for the creation of quantum computers and quantum informatics due to the fact that the influence of the environment on the process of quantum computation, which is realized by the dynamics of the quantum systems of qubits.
Currently, modern quantum mechanics is faced with the question of the most general form of the equations that describe the non-Markovian dynamics of quantum systems. A general form of non-Markovian character of quantum processes can be interpreted as nonlocality in time, or memory as its special form. Nonlocality in time means a dependence of the dynamics of observables and states at the current moment of time in the history of the system’s behavior on a finite interval in the past.
Attempts to construct a non-Markovian theory of open quantum systems with non-locality in time have been actively undertaken in recent years (see, for example, reviews [18,19,20], and articles [21,22,23,24,25] and references therein). Most approaches to the formulation of non-Markovian quantum dynamics were not related to fractional calculus and mathematical theory of equations with derivatives and integrals of non-integer orders (see [26,27,28,29,30,31,32]). Equations with fractional derivatives and integrals of non-integer orders are a powerful tool to take into account nonlocality in space and time in physics [33,34]. It should be emphasized that the characteristic property of fractional differential operators of non-integer orders is nonlocality in space and time [35,36].
For the first time, the use of fractional derivatives and integrals of non-integer order with respect to time to describe non-Markovian dynamics of OQS (in the form of power-law fading memory) was proposed in [37] (see Chapter 20 in book [37] (pp. 477–482) and [38,39,40]). Exact solutions of generalized Lindblad equations, which describe non-Markovian quantum dynamics, were derived in these works.
Fractional powers of Lindblad superoperators were first defined and used to describe OQS in [5] (pp. 433–444), (see also [37] (pp. 458–464, 468–477), and [41,42,43]). Solutions of equations, which describe non-Markovian quantum dynamics, were derived in works [5,37,41,42]. The generalization, which used the Grunvald–Letnikov fractional derivatives, is suggested in [44].
We also note some other formulations of the non-Markovian quantum dynamics by using fractional calculus: the non-Markovian dynamics of OQS with time-dependent parameters [39,45], relativistic open classical systems [46], and quantum systems with memory [47].
Note that the non-Markovian quantum dynamics is also considered in the framework of generalizations of the Schrodinger and Heisenberg equations (see [5,37,43,48,49,50] (pp. 457–466), respectively).
In fractional calculus, nonlocality is described by the kernel of the operators, which are fractional integrals (FI) and fractional derivatives (FD) of non-integer orders. To describe dynamical systems with various types of nonlocality in space and time, we can use operators with various types of kernels. Therefore it is important to have a general fractional calculus that allows us to describe nonlocality in the most general form.
The concept of general fractional calculus (GFC) has been suggested by Anatoly N. Kochubei in his work [51] in 2011 (see also [52,53,54]) out of connection with quantum theory. In works in [51,52] the general fractional derivatives (GFD) and general fractional integral (GFI) are defined. For these operators, the general fundamental theorems were proven in [51,52]. This approach to GFC is based on the concept of kernel pairs, which was proposed by N.Ya. Sonin (1849–1915) in their 1884 article [55] (see also [56]. Note that the name “Sonin” [57] is mistakenly used in French transliteration as “N. Sonine” from French journals [55]). Then, very important results in constructing the GFC were derived by Yuri Luchko in 2021 [58,59,60]. In works [58,59], GFD and GFI of arbitrary order are suggested, and the general fundamental theorems for the GFI and GFDs were proven. Operational calculus for equations with general fractional derivatives was proposed in [60]. The general fractional calculus was also developed in works [61,62,63,64,65,66,67,68,69,70,71,72] devoted to mathematical aspects and some applications in classical physics.
The first application of the general fractional study is proposed in this article. The general fractional calculus is suggested as a mathematical basis to formulate the non-Markovian quantum dynamics. In the framework of the general approach, it was assumed and implied to obtain, first of all, general results that do not depend on specific types (particular implementations) of nonlocality and kernels [72]. In the general approach to non-Markovian quantum dynamics, all research and results should be related to the general form of nonlocality, operator kernels of almost all types (a wide set of operator kernels). In this paper, we consider general fractional operators with kernels that belong to the Sonin set. The general form of non-Markovian dynamics of open quantum systems is described by the equations with general fractional derivatives and integrals [58,59]. The exact solutions of these equations were derived by using the general operational calculus proposed in [60].

2. Markovian Dynamics of Quantum States and Observables

In this section, we will briefly describe the basics of Markovian quantum dynamics for fixing concepts and notations.
Quantum states can be described by density operators ρ that are normalized ( T r ρ   = 1 ), self-adjoint ( ρ = ρ ), positive ( ρ 0 ) operators. All these properties must be conserved in the time evolution S t :   ρ     ρ t . Therefore to describe Markovian quantum dynamics in general form, we should find such a time evolution for the density operators ρ t = S t ρ , which satisfy the following conditions for all t   0 , .
(1)
The self-adjoint condition
( ρ t ) = ρ t , ( forallt > 0 ) .
(2)
The positivity condition
ρ t > 0 , ( forallt > 0 ) .
(3)
The normalization condition
T r ρ t   = 1 , ( forallt > 0 ) .
If we consider the Markovian dynamics, then the semigroup condition is used
S t   S τ = S t + τ
for all t , τ 0 , and the condition S 0 ρ = ρ is represented as S t ρ   ρ at t 0 in the trace norm.
The Markovian quantum dynamics of quantum observables A t   =   Φ t A is described by the dual dynamical maps Φ t , where the duality is represented by the condition
T r A ρ t   = T r A t ρ , ( forallt > 0 ) .
Let A be an algebra of bounded operators (for example, C -algebra, or algebra of B H of bounded operators on the Hilbert space H ).
Normalization condition (3) means that Φ t I   = I , where I is the identity operator.
Positivity condition (2) means that Φ t A A     0 . The positivity condition is usually replaced by the complete positivity condition for the dynamical maps Φ t :
k = 1 n l = 1 n A k Φ t B k   B l A l 0
for all A j , B j A , j = 1 , n and all n .
To describe the Markovian dynamics of quantum observables, the semigroup condition is used
Φ t   Φ τ =   Φ t + τ
for all t , τ 0 , and the condition Φ 0 A   = A is represented as Φ t A     A at t 0 in ultraweak operator topology. The ultraweak operator topology can be defined as the topology induced on A by the set of all seminorms of the form
A σ , ρ   =   T r ρ   A .
For the complete positive semigroup Φ t :   t   0 , , for which Φ t I   = I and Φ t A     A ultraweakly at t 0 , there exists a superoperator L for which
d d t Φ t A   =   L Φ t A
holds for all A D L , where D L is a ultraweakly dense domain of A . The superoperator L is called the generator of the semigroup Φ t . The dual generator Λ of the semigroup S t is connected with L through the equation
T r L A   ρ   =   T r A   Λ ρ .
For superoperator Λ , we have
d d t S t ρ   =   Λ S t ρ .
The bounded superoperators are defined on the normed operator space A : The superoperator L is called bounded, if
L A A     c A A
for some constant c and all A A . Inequality (12) means that L transforms norm bounded sets of A into the norm bounded sets. The least value of c equal to
L   =   sup A 0 L A A A A ,
is called the norm of the superoperator L . If A is a normed space, and L is a bounded superoperator, then
Λ   =   L .
Let us define the class of real superoperators: A real superoperator is a superoperator L on A , such that
[ L A ] = L A
for all A D L   A and A D L , where A be adjoint of A A . If L is a real superoperator, then Λ is real.
The Lindblad proved theorem defines the structure of the superoperators L and Λ . This theorem is proved for the completely positive dynamical semigroup Φ t , which is norm continuous
lim t 0 Φ t A     A A   = 0 .
Equation (16) is a more restrictive condition than ultraweak continuity. For such semigroup. the superoperator L is bounded.
For completely positive maps Φ t , the bi-positivity condition
Φ t A   A     Φ t A   Φ t A
holds for all t > 0 with equality at t = 0 . The differentiation of inequality (17) gives at t = 0 the property of the superoperator
L A   A     L A   A + A   L A
for all A A . The bounded superoperators L , which satisfy conditions L I   = 0 , L A   = ( L A ) and inequality (18), are called dissipative. The real superoperator L is called completely dissipative, if
L A k   A l     L A k   A l + A k   L A l ,
where k , l = 1 , , n , for all n and all A k ,   A l ,   A k A l   D L   A , where A k ,   A k A l   D L .
The Lindblad theorem gives the most general form of completely dissipative superoperators L .
The Lindblad theorem states [7] that the superoperator L is completely dissipative and ultraweakly continuous if and only if it has the form
L A   =   1 i H , A   +   1 2 k = 1 V k   A , V k   +   V k , A V k ,
where H = H A , V k A , k V   V k     A The dual superoperator Λ has the form
Λ ρ   =   1 i H , ρ   +   1 2 k = 1 V k   ρ , V k   +   V k ,   ρ   V k .
Equations (20) and (21) give explicit forms of most general Markovian dynamics of quantum observables and states.
The quantum Markovian equation, which are also called the Lindblad equations, are written in the form
d A t d t = 1 i H , A t   +   1 2 k = 1 V k A t , V k   +   V k , A t V k ,
d ρ t d t   =   1 i H , ρ t   +   1 2 k = 1 V k ρ t , V k   +   V k , ρ t V k ,
where A t is a quantum observable; ρ t is a quantum state; H is the Hamiltonian operator; and V k are the Lindblad operators [5]. Equations (22) and (23) are the standard time-local (memoryless) Markovian equations for quantum observables and states [6,7,8].
If V k = 0 for all k , then Equations (22) and (23) give the standard Heisenberg equation and von Neumann equation, respectively. In this case, Equations (22) and (23) describe Markovian dynamics of Hamiltonian quantum systems without nonlocality in time.

3. General Non-Markovian Dynamics of Quantum Observables and States

For the description of non-Markovian quantum processes, we can take into account nonlocality in time, which means that the behavior of the quantum observable A t (or state ρ t ) and its derivatives may depend on the history of previous changes of this operator. To describe this type of behavior, we cannot use differential equations of integer orders.
To take into account nonlocality in time (non-Markovality) in open quantum systems in work [37] (Chapter 20) and [38,39], it has been proposed to use the derivatives of non-integer orders instead of integer-order derivatives with respect to time.
In this section, we proposed using general fractional calculus and general fractional derivatives (GFD) as mathematical tools to allow us to take into account the general nonlocality in time for non-Markovian quantum processes.

3.1. Generalization of Lindblad Equation for Quantum Observables

The Lindblad equation for quantum observables is described by the operator differential equation of the first order
d A t d t = L   A t ,
where the superoperator L is defined by expression (20). Equation (24) can be written in the integral form
A t     A 0   =   0 t   d τ   L   A τ .
The nonlocality in time can be taken into account by using an integral kernel in generalization of Equation (25) in the form
A t     A 0   =   0 t d τ   M t τ   L   A τ ,
where the function M t τ describes nonlocality in time. If M t   = 1 for all t     0 , , then Equation (26) gives standard Equation (24), which describes the Markovian quantum dynamics.
Obviously, not all kernels M t τ can describe nonlocality in time. The nonlocality requirement can be formulated as follows. If integral Equation (26) can be written as a differential equation of an integer order (or a finite system of such equations), then the process, which is described by Equation (26), is local in time. The first obvious example of such “local” kernel is
M t   = h α t   = t α 1 Γ α
with positive integer values of α = n . Examples of “local” kernels are also the probability density functions of the exponential distribution and the gamma distribution with integer shape parameters (the Erlang distribution) [73,74].
Note that kernel (27) with non-integer values of α > 0 define “non-local” kernel of the Riemann–Liouville and Caputo fractional operators [29].
Let us assume that the functions M t belongs to the space C 1 , 0 0 , , and suppose that there exists a function K t   C 1 , 0 0 , , such that the Laplace convolution of these functions is equal to one for all t   0 , . The function X t belongs to the space C 1 , 0 0 , , if this function can be represented in the form X t   = t p   Y t , where 1 < p < 0 , and Y t     c   0 , .
Definition 1.
The functions M t , K t are a Sonin pair of kernels, if the following conditions are satisfied
(1) 
The Sonin condition for the kernels M t and K t requires that the relations
0 t M t τ   K τ   d τ = 1
holds for all t   0 , .
(2) 
The functions M t , K t belong to the space
C 1 , 0 0 ,   = { X t :   X t   = t p   Y t ,     t > 0 ,   1 < p < 0 ,     Y t     C 0 , } .  
The set of such Sonin kernels is denoted by S 1 .
If the kernel pair M t , K t belongs to Sonin set S 1 [59], the kernel K t is called the associated kernel to M t . Note that if K t is the associated kernel to M t , then M t is the associated kernel to K t . Therefore, if M t , K t belongs to set S 1 , then I M t , D K t , and I K t , D M t can be used as the general fractional integrals (GFI) and general fractional derivatives (GFD).
To define GFI and GFD, we used Luchko’s approach to general fractional calculus, which is proposed in [58,59].
Definition 2.
Let M t   S 1 and X t   C 1 0 ,   = C 1 , 0 , . The general fractional integral (GFI) with the kernel M t   C 1 , 0 0 , is the operator on the space C 1 0 , , that is
I M t :     C 1 0 ,   C 1 0 , ,  
that is defined by the equation
I M t τ X τ   = M     X t   =   0 t d τ   M t τ   X τ .
If the functions M t and K t belong to the Sonin set, then we can define general fractional derivatives D K t and D K t , that are associated with GFI I M t .
Definition 3.
Let M t , K t   S 1 and X t   C 1 1 0 , , i.e., X 1 C 1 0 , . The general fractional derivatives (GFD) with kernel K t   C 1 , 0 0 , , which is associated with GFI (31), is defined as
D K t , τ X τ   =   K     X 1 t   =   0 t d τ   K t τ   X 1 τ
for t   0 , . The GFD D K t is defined by the equation
D K t τ X τ   = d dt K     X t   =   0 t d dt d τ   K t τ   X τ
for t   0 , .
As proven in [58,59], operators (32) and (33) are connected (see Equation (47) in Definition 4 of [58] (p. 8)) by the equation
D K t , τ X τ   = D K t τ X τ     K t   X 0 .
The proposed GFI and GFD can be used to formulate non-Markovian dynamics in the general form, where the nonlocality in time is described by the kernel pairs that belong to the Sonin set S 1 .
If M t , K t   S 1 , then Equation (26) can be written through the GFI with kernel M t as
A t     A 0   = I M t τ   L   A τ .
The action of the GFD D K s t with kernel K t , which is associated to M t , in Equation (35), gives
D K s t   A t     A 0 D K s t   1 = D K s t   I M t τ   L   A τ .
For the right-hand side of Equation (36), we use the first fundamental theorem of GFC (see Theorem 3 of [58] (p. 9)). This theorem states that the equation
D K s t   I M t τ   f τ   = f s ,
holds for f t   C 1 0 , .
Using Equation (37) and the equality
D K s t   1   =   d d s 0 s d t   K s t   1   =   d d s 0 s d ξ   K ξ   =   K s ,
Equation (36) is written in the form
D K s t   A t     A 0   K t   = L   A s .
The left-hand side of Equation (39) can be expressed through GFD D K s , by using Equation (34) in the form
D K s t   A t     A 0 K s   = D K s , t   A t
for A t   C 1 1 0 , (see also Equations (47) and (49) in [58] (p. 8)).
As a result, Equation (36) can be written as
D K t , τ   A τ   = L   A t ,
where L is the Lindblad superoperator (the Lindbladian, quantum Liouvillian). Equation (41) describes non-Markovian dynamics of quantum observables in the general form, where the kernel belongs to the Sonin set S 1 .
As a result, we proved the following theorem.
Theorem 1.
The integral equation
A t     A 0   =   0 t d τ   M t τ   L   A τ ,
where M t   C 1 , 0 0 , and A t   C 1 1 0 , , L is the Lindblad superoperator (20), which can be represented in the form
D K t , τ   A τ   = L   A t ,  
if there exists the kernel K t   C 1 , 0 0 , , which is associated with M t such that the pair M t , K t belongs to the Sonin set S 1 .

3.2. Generalization of Lindblad Equation for Quantum States

Equation (23) is the standard memoryless Markovian quantum master equation [6,7,8]. Equation (23) can be written in the integral form
ρ t     ρ 0   =   0 t Λ   ρ τ   d τ ,
The nonlocality in time can be taken into account by using an integral kernel in Equation (44) in the form
ρ t     ρ 0   =   0 t M t     τ   Λ   ρ τ   d τ .
If M t   = 1 for t   0 , , then Equation (45) gives standard Equation (44). In general, the kernel M t τ can be used to describe non-Markovian quantum dynamics.
As a result, the non-Markovian master equation for quantum states takes the form
D K t , τ   ρ τ   =   Λ   ρ t ,
if M t , K t   S 1 . Equation (46) describes non-Markovian dynamics of quantum states in the general form, where the kernel K t belongs to the Sonin set S 1 .

3.3. Luchko Functions

Let us consider the triple R 1 =   C 1 0 , , , + , where the multiplication is the Laplace convolution and + the standard addition of functions. The triple R 1 is a commutative ring without divisors of zero [58,75].
The solutions of Equations (41) and (46) can be derived by using the Luchko operational calculus [60]. To describe the solution, we give the following theorem and define the Luchko function.
Theorem 2.
Let M t be a kernel from the Sonin set S 1 and the power series
f x , λ   =   j = 0   λ j 1   x j
has non-zero convergence radius r . Then the convolution series
F M , λ , t   =   j = 0   M , j t   λ j 1
is convergent for all t   0 , and the function F x , λ , t belongs to the ring R 1 .
Theorem 2 is proven in [60] (see Theorem 4.4 in [60] (p. 359) and comments on [60] (p. 360)).
Definition 4.
Let M t be a kernel from the Sonin set S 1 , and M , j t is the convolution j -power:
M , j t   : =   M 1         M j t ,  
where M k t   =   M t for all k = 1 ,   ,   j , and t   0 , .
Then, the function
F M , λ , t   =   j = 0   M , j t   λ j 1
will be called the first Luchko function.
Let us give examples of the first Luchko function [60] (p. 361).
(1)
For the Sonin kernel
M t   = h α t   = t α 1 Γ α ,   M , j t   = h j α t ,
the first Luchko function has the form
F M , λ , t   =   j = 1 λ j 1 h j α t   = t α 1 j = 1 λ j   t j α Γ j α + α = t α 1 E α , α λ   t α .
Here, E α , β z is the two-parameters Mittag–Leffler function [76] that is defined as
E α , β z   : =   k = 0 z k Γ α k + β ,
where α > 0 , β (or β ).
(2)
For the Sonin kernel
M t   = h α , β t   = h α t e β t , M , j t   = h j α , β t .
the first Luchko function takes the form
F M , λ , t   =   j = 1 λ j 1 h j α , β t   = e β t t α 1 j = 1 λ j   t j α Γ j α + α = e β t t α 1 E α , α λ   t α .
(3)
For the Sonin kernel
M t   = h 1 β + α t   +   h 1 β t ,
where 0 < α < β < 1 , the first Luchko function is given as
F M , λ , t   =   j = 1 λ j 1 i = 0 j j i h i α + j 1 β t   = 1 λ t j = 1 l 1 + l 2 = j j ! l 1 ! l 2 ! ( λ t 1 β + α ) l 1 ( λ t 1 β ) l 2 Γ l 1 1 β + α   +   l 2 1 β = 1 λ t E 1 β , 1 β + α , 0 λ t 1 β , λ t 1 β + α ,
where E 1 β , 1 β + α , 0 is the binomial (multinomial) Mittag–Leffler function [29] (p. 49) and [77], which is defined as
E α 1 , , α m , β z 1 , , z m   : = j = 0 l 1 + + l m = j j ! i = 1 m l i ! i = 1 m z i l i Γ i = 1 m α i l i + β
Using the first Luchko function F M , λ , t and the kernel K t that is associated with M t , we can define the second Luchko function.
Definition 5.
Let the functions M t , K t belong to the Sonin set S 1 , and the first Luchko function is
F M , λ , t   = j = 1   M , j t   λ j 1 .
Then, the function
L M , λ , t   : =   0 t d τ   K t τ F M , λ , τ   =   I K t τ F M , λ , τ ,
will be called the second Luchko function.
Note that Equation (60) contains the GFI with kernel K t rather than kernel M t .
The second Luchko function (60) is used [60] in solution of equations with GFD that is defined by the kernel K t associated with the kernel M t of the GFI.
To derive solutions of equations with GFD, Yu. Luchko proposed the general operational calculus [60]. Theorem 5.1 of [60] (p. 366) proves that the solution of the equation
D K t , τ   A τ   = λ   A t ,
where A t   C 1 1 0 , is expressed through the second Luchko function L M , λ , τ .
If M t , K t   S 1 , then F M , λ , t   C 1 0 , and L M , λ , t   C 1 0 , . Therefore, these Luchko functions belong to the ring R 1 . Thess statements are based on the fact that GFI I K t τ is the operator on C 1 0 , (see Equation (30) and reference [60]).
Note that the second Luchko function (60) can be considered as independent of the kernel K t since the Sonin condition K     M t   = 1 is satisfied for all t   0 , .
Using the superoperator form of the first Luchko function F M , λ , t and the second Luchko function L M , λ , t , we can propose solutions of equations for non-Markovian open quantum systems with nonlocality in time.

3.4. General Form of Solutions for Non-Markovian Equations

Let L be a bounded superoperator on the normed operator algebra A , in other words,
λ   : =   L   < ,
and L 0 = L I is the unit superoperator ( L I   A = A for all A A ). The superoperator power series
F x , L   : =   j = 1 x j   L j 1
converges in norm and the radius of convergence is equal to
r = λ 1 =   L 1 ,
if L is the bounded superoperator on the normed operator algebra A .
Using Theorem 2, we can state that the series
F M , L , t   : =   j = 1 M , j t   L j 1
converges in norm for all t   0 , , if M t   C 1 0 , and L is the bounded superoperator.
The solution of Equation (41) can be expressed through the superoperator (65).
To obtain the solution of the general non-Markovian equation for quantum observable A A , we will use Theorem 5.1 of [60] (p. 366).
Theorem 3.
Let A t   C 1 1 0 , , the pair M t , K t is Sonin pair from S 1 , and L is a bounded superoperator on A . Then the initial value problem
D K t , τ   A τ   = L   A t ,         A 0   = A 0  
has the unique solution
A t   =   Φ M t A 0 ,
where
Φ M t   =   L M , L , t   =   0 t   d τ   K t τ   F M , L , τ   =   I K t τ   F M , L , τ .
The proof of this theorem is based on Theorem 5.1, which was proven in [60] (pp. 366).
Remark 1.
Note that superoperator Φ M t is independent of the kernel K t due to the Sonin condition
M     K t   =   1 .  
Using condition (69), we get the convolution of kernel K = K t and L = F M , L , t in the form
Φ M t   =   K     L t   =   j = 1 K     M * , j t   L j 1   =   j = 1   K     M     M * , j 1 t   L j 1   =   j = 1 1     M , j 1 t   L j 1   =   L I   +   1     j = 2 M , j 1   L j 1 t   =   L I   +   1     j = 1 M * , j   L j t ,
where L 0 = L I is the unit superoperator ( L I   A = A for all A A ).
Therefore the map Φ M t is written as
Φ M t   = 0 t   K t τ   F M , L , τ   =   L I   +   0 t j = 1   M * , j τ   L j d τ .
In the next section, we give some examples of the general non-Markovian Equation (66) and solutions (67).

3.5. Example of General Non-Markovian Dynamics

Let us consider some special cases of the proposed general non-Markovian equations and its solutions that are derived by the Luchko operational calculus [60].
(1)
In the first example, we consider the Sonin pairs of the kernels
M t   = h α t   = t α 1 Γ α , K t   = h 1 α t ,
where t > 0 , and 0 < α < 1 . In this case, the GFD D K t , is the Caputo derivative D C , 0 + α with α   0 , 1 . Then
Φ M t   =   L I   +   1     j = 1   h j α t   L j   =   L I +   j = 1 h j α + 1 t   L j   =   E α t α   L ,
where the Mittag–Leffler function E α z   = E α , 1 z [76]. The fractional differential equation
D K t , τ A τ   = L A t , A 0   = A 0 ,
has the solution
A t   = E α t α   L   A 0 .
This type of non-Markovian quantum dynamics was first described in [37] (Chapter 20) and [38,39].
(2)
In the second example, we consider the Sonin pairs of the kernels
M t   = h α , β t   = h α t   e β t , K t   = h 1 α t   e β t + β α Γ 1 α   γ 1 α , β t ,
where t > 0 , α   0 , 1 , β > 0 , and γ β , t is the incomplete gamma function
γ β , t   =   0 t τ β 1 e τ   d τ .
The GFD has the form
D K t , τ   A τ   =   0 t h 1 α τ   e β τ   +   β α Γ 1 α   γ 1 α , β τ   A 1 t τ   d τ .
The superoperator Φ M t is written as
Φ M t   =   K   F M , L , τ t   =   h 1 α , β + β 1     h 1 α , β , j = 1 h j α , β   L j 1 t   =   j = 1 h j 1 α + 1 , β t L j 1   +   β 1   j = 1 h j 1 α + 1 , β   L j 1 t   =   e β t   E α t α   L   +   β   0 t d τ   e β τ E α τ α   L ,
where we use h α , β     h γ , β   = h α + γ , β .
For nonlocality (76), the equation of non-Markovian dynamics
D K t , τ A τ   = L A t , A 0   = A 0 ,
has the solution
A t   =   e β t   E α t α   L   +   β 0 t   d τ   e β τ E α τ α   L   A 0 .
(3)
In the third example, we consider the Sonin pairs of the kernels
M t   = h 1 β + α t   +   h 1 β t , K t   = t β 1   E α , β t α ,
where 0 < α < β < 1 . In this case, we use the GFD that has the Mittag–Leffler function in the kernel
D K t , τ   A τ   =   0 t   d τ   τ β 1   E α , β τ α   A 1 t τ .
Then, the non-Markovian dynamics is described by the superoperator
Φ M t   =   L I   +   1   j = 1 M , j t L j = L I   +   1   j = 1 L j i = 0 j j i h i α + j 1 β t   =   E 1 β , 1 β + α , 1 t 1 β   L , t 1 β + α   L ,
where the function E 1 β , 1 β + α , 1 is a special case (binomial) of the multinomial Mittag–Leffler function [77]. The equation with GFD (83) and initial condition
D K t , τ A τ   = L A t , A 0   = A 0 ,
has the solution
A t   = E 1 β , 1 β + α , 1 t 1 β L , t 1 β + α   L   A 0 .
These examples of equations and their solutions describe the non-Markovian dynamics of quantum observables.

4. Properties of Non-Markovian Quantum Dynamical Maps

4.1. Violation of Semigroup Property for Non-Markovian Maps

The non-Markovian maps Φ t M α , t > 0 describe dynamics of open quantum systems with power-law memory. The superoperator L can be considered as a generator of the one-parameter groupoid Φ t M α on an operator algebra of quantum observables:
D K t , τ   Φ τ M = L   Φ τ M .
The set Φ t M   | t   >   0 , is called a quantum dynamical groupoid [5,38]. Note that the following properties are realized
Φ t M   I = I ,
( Φ t M   A ) = Φ t M   A
for self-adjoint operators A ( A = A ), and
lim t 0 + Φ t M = L I ,
where L I is an identity superoperator ( L I   A = A ). As a result, the non-Markovian map Φ t M , t > 0 , are real and unit preserving maps on the operator algebra of quantum observables.
For M t   = h α t and K t   = h 1 α t with 0 < α 1 , the GFD is the Caputo fractional derivative of the order α , and the non-Markovian map has the form
Φ t M = E α t α   L
that is described in [37,38], where E α z is the Mittag–Leffler function [76]. For α = 1 , we have
Φ t M = E 1 t   L   =   exp t   L   =   Φ t .
The superoperators Φ t form a semigroup such that
Φ t   Φ s =   Φ t + s , ( t , s > 0 ) ,   Φ 0 = L I .
This property holds since
exp t   L   exp s   L   =   exp t + s   L V .
For α , we have the violation of the semigroup property [78,79,80]:
E α [ t α   L   E α s α   L     E α t + s ) α   L .
Therefore, the semigroup property is not satisfied for non-Markovian dynamics
Φ t M   Φ s M     Φ t M , ( t , s > 0 ) .
As a result, the non-Markovian maps Φ t M cannot form a semigroup. This property is a characteristic property, which means that we have a quantum process with non-locality in time. The maps Φ t M describe the quantum dynamics of open systems with non-locality in time. This non-locality in time means that their present value of quantum observable (or quantum state) of quantum system A t   = Φ t M A 0 depends on all past values of A τ for τ   0 , t .

4.2. Some Properties of Markovian Maps

In this section, we will briefly describe the properties of Markovian maps Φ t and the superoperator L , for the convenience of generalizations to non-Markovian dynamics.

4.2.1. Bi-Positivity and Dissipativity in Markovian Theory

The Markovian quantum dynamics is described by the map
Φ t   : =   n = 0 h n + 1 t   L n ,
where
h n + 1 t   = t n Γ n + 1 = t n n ! .
We will assume that D L   = A to simplify the description. To this purpose, we will also consider the bi-positivity condition instead of complete positivity condition.
The bi-positivity condition for the Markovian map Φ t can be considered in the form
Φ t A   A     Φ t A   Φ t A ,
which should be satisfied for all t > 0 and A A .
The importance of condition (88) is due to the fact that it leads to the positive condition
Φ t A   A     0 ,
for all A A , if Φ t A   = Φ t A . Let us prove this statement. If L A   = ( L A ) , then
Φ t A   =   ( Φ t A ) ,
and
Φ t A   Φ t A   =   Φ t A   ( Φ t A ) = Φ t A 2     0 .
As a result, we obtain
Φ t A   A     Φ t A 2     0 .
If the real superoperator L is completely dissipative, for which the inequality
L A k   A l     L A k   A l + A k   L A l
is satisfied for all A k , A l A , then the quantum Markovian map map Φ t is completely positive, if Φ t A   = Φ t A for all A A . This statement can be proved similarly by using the following transformations
k = 1 n l = 1 n B k Φ t A k   A l B l k = 1 n B k Φ t A k l = 1 n Φ t A l B l   = l = 1 n Φ t A l B l l = 1 n Φ t A l B l   = l = 1 n Φ t A l B l 2     0 .
Let us consider two approaches to find a condition that the real superoperator L must satisfy in order for the bi-positivity condition to be satisfied for all t > 0 in the Markovian quantum dynamics, and a possibility to generalize these approaches to the general non-Markovian maps.

4.2.2. Markovian Case: First Approach

The condition on L to have the bi-positivity of maps Φ t is usually obtained by differentiating inequality (88) with respect to time
d d t Φ t A   A     d d t Φ t A   Φ t A
and using the standard Leibniz rule
d d t Φ t A   Φ t A   = d d t Φ t A   Φ t A   +   Φ t A   d d t Φ t A .
The Markovian equations for quantum observables
d d t Φ t A   = L   Φ t A   =   Φ t L A
and Equation (91) allow us to get inequality (90) in the form
Φ t L A   A     Φ t L A   Φ t A   +   Φ t A   Φ t L A .
In the limit t     0 , we get the condition
L A   A     L A   A + A   L A .
Unfortunately, this approach cannot be used for equations with fractional derivatives and GFD because the standard Leibniz rule (the product rule) is violated
D K t , τ Φ τ A   Φ τ A     D K t , τ Φ τ A   Φ t A   +   Φ t A   D K t , τ Φ t A .
For example, in the non-Markovian quantum theory, which was proposed in [37] (pp.477–482) and [38,39], the Sonin pair of kernels is used in the form
M t   = h α t   = t α 1 Γ α , K t   = h 1 α t   = t α Γ 1 α ,
with 0 < α < 1 . In this case, the GFI and GFD are the Riemann–Liouville fractional integral and the Caputo fractional derivative. The generalized Leibniz rule (see Theorem 3.17 in [30] (p. 59)) has the form
D h 1 α t , τ   f τ   g τ   =   h 1 α t   f t     f 0   g 0   + f t   D h 1 α t , τ   g τ   +   k = 1 α k   f k t   I h k α t τ   g τ ,
where α k   : = α α 1 α k + 1 k ! .
The violation of the standard Leibniz rule is a characteristic property of fractional derivatives of non-integer order [81].
For GFD, there is no rule for differentiating the product in the general case. Therefore, we should use another approach to derive the conditions on the superoperator L . For Markovian dynamics, another method of obtaining the condition on L can be used, and this method can be generalized to the case of non-Markovian quantum theory.

4.2.3. Markovian Case: Second Approach

Let us consider the bi-positivity condition in the form
Φ t A   A     Φ t A   Φ t A ,
which should be satisfied for all t > 0 and all A A , where the average value of the quantum observable A t   = Φ t A is defined as
Φ t A = T r ρ   Φ t A .
Substitution of expression (87) into inequality (95) gives
n = 0 h n + 1 t   L n A   A     m = 0 k = 0 h m + 1 t   h k + 1 t   L m A   L k A ,
where we use the linearity of the average value. Using n = m + k , and
h m + 1 t   h k + 1 t   = t m Γ m + 1   t k Γ k + 1   = Γ m + k + 1 Γ m + 1   Γ k + 1   t m + k Γ m + k + 1   =   m + k k   h m + k + 1 t   =   n k   h n + 1 t
inequality (96) is represented as
n = o h n + 1 t   L n A   A     n = o h n + 1 t   k = o n n k   L n k A   L k A .
Since the bi-positivity condition must be satisfied for all t > 0 , we obtain
L n A   A     k = 0 n n k   L n k A   L k A ,
which should be satisfied for all n . For n = 1 , inequality (98) has the form
L A   A     L A   A + A   L A ,
which should be satisfied for all A A .
The bounded superoperators L , which satisfy conditions L I   = 0 , L A   = ( L A ) and inequality (99), are called dissipative.
Inequality (99) for real superoperator L is a necessary and sufficient condition in order for the Markovian quantum map Φ t to have a bi-positivity property.
Theorem 4.
Let L be real superoperator, which satisfies the dissipativity condition
L A   A     L A   A + A   L A  
for all A A .
Then, the Markovian map Φ t satisfies the bi-positive condition
Φ t A   A     Φ t A   Φ t A
for all t > 0 and all A A .
Proof. 
Using the series representation of Φ t A   A and inequality
L A   A     L A   A + A   L A
we get
Φ t A A   =   n = 0 h n + 1 t L n A A   n = 0 k = 0 n n k h n + 1 t L n k A L k A   = m = 0 k = 0 m + k k h m + k + 1 t L m A L k A .
Using
n k   h m + k + 1 t   = Γ m + k + 1 Γ m + 1 Γ k + 1   h m + k + 1 t   = Γ m + k + 1 Γ m + 1 Γ k + 1   t m + k Γ m + k + 1   = t m   t k Γ m + 1 Γ k + 1   = h m + 1 t   h k + 1 t
we get
Φ t A   A     m = 0 k = 0 h m + 1 t   h k + 1 t   L m A   L k A   =   m = 0   h m + 1 t   L m A k = 0   h k + 1 t   L k A   =  Φ t A  Φ t A .
 □

4.3. General Non-Markovian Maps: Bi-Positivity and Complete Positivity

In this section, we will consider the bi-positivity condition instead of the complete positivity condition to simplify the description and proofs. Condition of complete positivity proved similar to the proofs for bi-positivity, and conditions for complete dissipativity of the superoperator L are written analogously to conditions for general dissipativity. We will also assume that D L   = A to simplify the description.
Let us give the definition for the non-Markovian quantum maps that are described in the paper.
Definition 6.
Let the pair of kernels M t ,   K t belongs to the Sonin set S 1 , where the Sonin condition has the form
K     M t   =   1 ,  
and L is real Lindblad superoperator.
Then, the one-parameter superoperator
Φ t M A   =   j = 1 K     M , j t   L j 1   =   n = 0 M n + 1 t   L n ,
where
M n + 1 t   : =   K     M , n + 1 t   =   1     M , n t   = 0 t   M , n τ   d τ ,
will be called the general non-Markovian quantum map.
For the general non-Markovian quantum dynamics, the bi-positivity property of the map Φ t M can be described in the following form.
Definition 7.
The bi-positivity property of the general non-Markovian quantum dynamical map Φ t M can be defined in the form of the inequality
Φ t M A   A     Φ t M A   Φ t M A ,  
which holds for all t > 0 and all A A . The average value of the quantum observable A t   = Φ t M A is defined as
Φ t M A = T r ρ   Φ t M A .  
Lemma 1.
Let the general non-Markovian map Φ t M be bi-positive and
Φ t M A = Φ t M A  
for all t > 0 and all A A .
Then, the positivity condition
Φ t M A   A     0
holds for all t > 0 and all A A .
Proof. 
The following equalities hold
Φ t M A   Φ t M A = A t   A t = A t 2     0
for all t > 0 and for all A A . Therefore, bi-positivity condition (104) leads to the positivity
Φ t M A   A     Φ t M A   Φ t M A = Φ t M A   Φ t M A = | Φ t M A | 2     0
in the form
Φ t M A   A     0 .
 □
Definition 8.
The complete positivity condition of the general non-Markovian quantum dynamical map Φ t M can be defined in the form
k = 1 n l = 1 n B k Φ t M A k   A l B l     0 ,
which holds for all t > 0 , and for all A j , B j A , j = 1 , n and all n .
Lemma 2.
Let the general non-Markovian map Φ t M satisfy the complete positivity condition, and
Φ t M A = Φ t M A  
for all t > 0 and all A A .
Then the positivity condition
Φ t M A   A     0 .
holds for all t > 0 .
Proof. 
Using A k = A and B k = I , inequality (105) takes form (106). □

4.3.1. From Bi-Positivity to General Dissipativity

Let us find a condition that the real superoperator L must satisfy in order to bi-positivity condition (104) to be satisfied for all t > 0 for the general non-Markovian quantum maps.
This required condition is given by the following theorem.
Theorem 5.
Let the bi-positivity condition
Φ t M A   A     Φ t M A   Φ t M A  
be satisfied for all t > 0 for the general non-Markovian map
Φ t M A   =   n = 0 M n + 1 t   L n ,
where
M n t   >   0  
for all n and all t > 0 .
Then real superoperator L satisfies the inequality
L n A , A k = 1 n E k n K , M L n k A L k A ,
for all n and all t > 0 , where
E k n K , M   : = M n k + 1 t   M k + 1 t M n + 1 t .  
Proof. 
Substitution of expression (102) in the right side of inequality (104) gives
Φ t M A   Φ t M A   =   m = 0 k = 0 M m + 1 t   M k + 1 t   L m A   L k A ,
where we use the linearity of the average value. Using m = n k , Equation (112) takes the form
Φ t M A   Φ t M A   =   m = 0 k = 0 n M n k + 1 t   M k + 1 t   L n k A   L k A .
Using (102) and (113), bi-positivity condition (104) is represented by the inequality
n = 0 M n + 1 t   L n A   A     n = 0 k = 0 n M n k + 1 t   M k + 1 t   L n k A   L k A .
Since inequality (114) must hold for all t > 0 , we will look for conditions on L by using the inequalities
M n + 1 t   L n A   A     k = 0 n M n k + 1 t   M k + 1 t   L n k A   L k A .
Using the assumption that the inequality
M n t   >   0
holds for all n and all t > 0 , we obtain the condition
L n A   A     k = 0 n E k n K , M   L n k A   L k A ,
where
E k n K , M   : = M n k + 1 t   M k + 1 t M n + 1 t .
and k , n 0 and 0 k n . □
Definition 9.
Let a pair of kernels M t ,   K t belong to the Sonin set S 1 , and
M n t   >   0    
for all n and all t > 0 . Then the function
E k n K , M   : = M n k + 1 t   M k + 1 t M n + 1 t ,  
where k ,   n 0 and 0 k n , will be called the general binomial coefficients.
Remark 2.
If the kernel M t is positive
M t   >   0  
for all t > 0 , then condition (119) holds.
Definition 10.
Let the real operator L satisfy the inequalities
L n A   A   k = 0 n   E k n K , M   L n k A   L k A
for all n 0 and all A A . Then L will be called the general dissipative superoperator.
The general complete dissipativity condition is defined in the form
L n A i     , A j     k = 0 n E k n K , M   L n k A i   L k A j
which holds for all t > 0 , and for all A i , A j A , j = 1 , m and all n , m .
We will consider the bi-positivity condition instead of the complete positivity condition to simplify the description and proof in the next sections. Conditions for complete positivity were proved similarly to the proofs for bi-positivity. The condition for general complete dissipativity of the real superoperator L will be used analogously as the condition for general dissipativity.
In connection with general dissipativity condition (121), which must be satisfied for any value of all n , two questions arise: (1) What condition must the superoperator L 1 = L satisfy in order for the general dissipativity condition (121) to hold for all n ? (2) What is the connection between the general dissipativity condition, and the condition of the dissipativity? Answers to these questions will be offered in the following sections.

4.3.2. General Dissipativity for n = 1

Let us consider the general dissipativity condition for the superoperator L n with n = 1 .
Theorem 6.
Let a pair of kernels M t ,   K t belongs to the Sonin set S 1 , and
M n t   >   0  
for all n and all t > 0 .
Then the general binomial coefficients with k = 0 and k = n are equal to one
E 0 n K , M   = E n n K , M   = 1
for all n .
Proof. 
Using the definition of the general binomial coefficients, we have
E 0 n K , M   : = M n + 1 t   M 1 t M n + 1 t ,
E n n K , M   : = M 1 t   M n + 1 t M n + 1 t .
Using
M 1 t   : =   K     M t   =   1 ,
we get
E 0 n K , M   = E n n K , M   = 1   M n + 1 t   M n + 1 t   = 1 .
 □
Theorem 7.
Let a pair of kernels M t ,   K t belongs to the Sonin set S 1 , and
M n t   >   0  
for all n and all t > 0 .
Then the general dissipative superoperator L n with n = 1 satisfies the condition
L A   A     L A   A + A   L A
for all A A .
Proof. 
For n = 1 , inequality (110) takes the form
L A , A )     E 0 1 K , M   L A   A + E 1 1 K , M   A   L A ,
where
E 0 1 K , M   : = M 2 t   M 1 t M 2 t ,
E 1 1 K , M   : = M 1 t   M 2 t M 2 t .
Using Theorem 7, we have
E 0 1 K , M   = E 1 1 K , M   = 1   M 2 t M 2 t   =   1 .
Therefore, inequality (127) has the form
L A , A     L A   A + A   L A .
Theorem 8.
Let Φ t M be general non-Markovian map with the positive kernel M t   >   0 for all t > 0 such that the bi-positivity condition
Φ t M A   A     Φ t M A   Φ t M A  
is satisfied for all t > 0 , and all A A .
Then, the real superoperator L satisfies the dissipativity condition
L A   A     L A   A + A   L A q
holds for all A A .
Proof. 
Using Theorem 5, we get that the bi-positivity condition (132) for general non-Markovian map Φ t M leads to the general dissipativity condition for the superoperator L n . □
Then, using Theorem 6, the general dissipativity condition leads to the dissipativity condition (133).

4.3.3. From General Dissipativity to Bi-Positivity

Let us prove the theorem that is converse to Theorem 5.
Theorem 9.
Let L be real superoperator, for which the inequalities
L n A   A   k = 0 n   E k n K , M   L n k A   L k A  
are satisfied for all n , and all A A , where the functions M n + 1 t are positive
M n + 1 t   : =   K     M , n + 1 t   =   1     M , n t   >   0
for all n 0 .
Then the bi-positivity condition
Φ t M A   A     Φ t M A   Φ t M A  
is satisfied for all t > 0 and all A A for the general non-Markovian map
Φ t M A   =   n = 0 M n + 1 t   L n .
Proof. 
Let us use condition (134) in the form
L n A   A     k = 0 n M n k + 1 t   M k + 1 t M n + 1 t   L n k A   L k A .
Using assumption that M n t   >   0 holds for all n and all t   >   0 , inequality (138) can be written as
M n + 1 t   L n A   A     k = 0 n M n k + 1 t   M k + 1 t   L n k A   L k A .
Since, by condition of the theorem, inequality (138) holds for any n 0 . Then, summing from 0 to inequality (139), we obtain
n = 0 M n + 1 t   L n A   A     n = 0 k = 0 n M n k + 1 t   M k + 1 t   L n k A   L k A .
Using m = n k , condition (140) takes the form
n = 0 M n + 1 t   L n A   A     m = 0 k = 0 M m + 1 t   M k + 1 t   L m A   L k A .
Then
n = 0 M n + 1 t   L n A   A     m = 0 k = 0   M m + 1 t   M k + 1 t   L m A   L k A .
n = 0 M n + 1 t   L n A   A     m = 0 M m + 1 t   L m A     k = 0 M k + 1 t   L k A .
Φ t M A   A     Φ t M A   Φ t M A .
 □
As a result, Theorems 5 and 9 allow us to formulate the following statement.
Corollary 1.
A necessary and sufficient condition for the bi-positivity property of the general non-Markovian quantum dynamical map Φ t M with positive kernel M t in the form
Φ t M A   A     Φ t M A   Φ t M A  
that holds for all t > 0 , is the general dissipativity of the real superoperator L in the form
L n A   A   k = 0 n   E k n K , M   L n k A   L k A  
that satisfies for all n .
Remark 3.
Let the inequality
L n A i   A j   k = 0 n   E k n K , M   L n k A i   L k A j .
be satisfied for all n and all A i , A j A , where i , j = 1 , , m .
Then, the completely positive condition holds in the form
i = 1 m i = 1 m B i Φ t M A i   A j B j 0 .
for all n and all A i , B i A , where i , j = 1 , , m .
The proof of this statement is realized similarly to Theorem 9.
Corollary 2.
The condition that a superoperator is dissipative is a necessary condition for non-Markovian maps to be bi-positive, but it is not a sufficient condition.

4.3.4. From General Dissipativity to Dissipativity

Let us consider the relationship between the concepts of dissipativity and general dissipativity.
Theorem 10.
Let real superoperator L satisfy the general dissipativity condition
L n A   A   k = 0 n   E k n K , M   L n k A   L k A
for all n and all A A .
Then, the superoperator L satisfies the dissipativity condition
L n A   A   k = 0 n   n k   L n k A   L k A
for all n .
Proof. 
Using Theorem 7, we get that general dissipativity condition (148) with n = 1 has the form
L A   A     L A   A + A   L A .
Using inequality (150) two times for the expression L 2 A   A and the linearity property of the superoperator L , we obtain
L 2 A   A = L L A   A     L L A   A + A   L A = L L A   A   +   L A   L A     L 2 A   A + 2   L A   L A   +   A   L 2 A .
Using inequality (150) n -times for the expression L n A   A and the linearity property of the superoperator L in a similar way, we obtain inequality (149). □
Corollary 3.
Let the general non-Markovian map Φ t M with positive kernel M t satisfy the bi-positivity condition in the form of the inequality
Φ t M A   A     Φ t M A   Φ t M A ,  
which holds for all t > 0 and all A A .
Then, the superoperator L satisfies the dissipativity condition
L n A   A     k = 0 n n k   L n k A   L k A
for all n , and all A A .
Remark 4.
The dissipativity condition
L n A   A     k = 0 n n k   L n k A   L k A
which holds for all n does not lead to the fact that the general non-Markov map Φ t M is bi-positive.
If we assume in non-Markovian quantum theory that L is dissipative superoperators only, then the standard bi-positive condition
Φ τ A   A   Φ τ A   Φ τ A     0 .
is replaced by the inequality
0 t   d τ   Φ τ M A   A     Φ t τ M A   Φ τ M A     0 .
that should be satisfied for all t 0 . Let us prove this statement.
Theorem 11.
Let real superoperator L satisfy the inequality
L A   A     L A   A + A   L A .  
Then, the inequality
0 t Φ τ M A   A   d τ     0 t Φ t τ M A   Φ τ M A   d τ
holds for all t > 0 and all A A , where Φ τ M is general non-Markovian map.
Proof. 
Using the proof of Theorem 10, we get that the repeated action of the inequality (154) gives the condition
L n A   A     k = 0 n n k   L n k A   L k A .
Using the series representation of Φ t A   A and inequality (18), we get
Φ t M A   A   =   k = 0 M n + 1 t   L n A   A     n = 0 k = 0 n   M n + 1 t   n k   L n k A   L k A   =   m = 0 k = 0   M m + k + 1 t   m + k k   L m A   L k A .
Using
m + k k     1
for k , n 0 , 0 k n , and M n + 1 t   >   0 for all n 0 , we can write the inequality
Φ t M A   A     m = 0 k = 0 M m + k + 1 t   L m A   L k A .
The convolution with 1 means the integration in the form
0 t f τ   d τ   =   1     f t ,
Therefore, multiplying the left and right sides of the inequality by 1 and using Sonin’s condition
K     M t   =   1 ,
We obtain
1     Φ t M A   A t     1     m = 0 k = 0 M m + k + 1 t   L m A   L k A t   =   m = 0 k = 0   1     M , m     1     M , k t   L m A   L k A   =   m = 0   1     M , m t   L m A     k = 0   1     M , k t   L k A   =   Φ t M A     Φ t M A   =   0 t Φ t τ M A Φ τ M A   d τ .
 □
Remark 5.
Alternative forms of writing this inequality
0 t Φ τ M A   A     Φ t τ M A Φ τ M A   d τ     0 , Φ t M A   A     d d t 0 t Φ t τ M A   Φ τ M A   d τ .
Remark 6.
Similarly, we can prove the inequality
0 t d τ k = 1 n l = 1 n B k Φ t M A k   A l B l   0 t d τ k = 1 n Φ t τ M A k B k )   ( l = 1 n Φ τ M A l B l ,
where we use Φ t M A   = ( Φ t M A ) .

4.3.5. From Dissipativity to Bi-Positivity

Let us consider the relationship between the dissipativity and bi-positivity in general non-Markovian dynamics.
Theorem 12.
Let real superoperator L satisfy the dissipativity condition
L n A   A     k = 0 n n k   L n k A   L k A
for all n and all A A , and the general binomial coefficients satisfy the condition
n k     E k n K , M  
for all t > 0 and all k , n 0 , 0 k n .
Then, the superoperator L satisfies the general dissipativity condition
L n A   A   k = 0 n   E k n   L n k A   L k A
for all n and all A A .
Then, the general non-Markovian map Φ t M satisfies the bi-positivity condition
Φ t M A   A     Φ t M A   Φ t M A ,
which holds for all t > 0 .
Proof. 
Using inequality (156), we get
k = 0 n n k   L n k A   L k A     k = 0 n E k n K , M   L n k A   L k A .
Then, using inequalities (155) and (159), we obtain the general dissipativity condition
L n A   A     k = 0 n E k n K , M   L n k A   L k A
for all n .
Using Theorem 9, condition (160) gives the bi-positivity condition (158). □
Theorem 13.
Let the general non-Markovian map Φ t M with positive kernel M t satisfy the bi-positivity condition in the form of the inequality
Φ t M A   A     Φ t M A   Φ t M A ,  
which holds for all t > 0 and all A A , and the general binomial coefficients satisfy the condition
E k n K , M     n k  
for all t > 0 and all k , n 0 , 0 k n .
Then, the superoperator L satisfies the dissipativity condition
L n A   A     k = 0 n n k   L n k A   L k A
for all n and all A A .
Proof. 
Using Theorem 5, bi-positivity condition (161) gives the general dissipativity condition
L n A   A     k = 0 n E k n K , M   L n k A   L k A
for all n and all A A . Then using the condition (162) for the general binomial coefficients, we obtain
k = 0 n E k n K , M   L n k A   L k A     k = 0 n n k   L n k A   L k A
for all n and all A A . Then, using conditions (164) and (165), we derive the dissipativity condition (163).

4.3.6. Examples of General Binomial Coefficients

Let us consider the general binomial coefficients for the non-Markovian quantum maps with the pair of kernels
M t   = h α t   = t α 1 Γ α , K t   = h 1 α t   = t α Γ 1 α
that belongs to the Sonin set for α   0 , 1 . We see that the condition M t   > 0 holds for all t > 0 .
Let us give a definition of the generalized binomial coefficients (see [26] (p. 14–15) and [29] (p. 26–27)). These well-known “generalized” coefficients should not be confused with the “general” coefficients suggested in this article.
Definition 11.
The generalized binomial coefficients are defined by the equation
α β   : = Γ α + 1 Γ β + 1 Γ α β + 1 ,  
where α , β , α 1 , 2 ,
Let us prove that the general binomial coefficients for kernel pair (166) are expressed through the generalized binomial coefficients.
Theorem 14.
Let the pair of kernels M t ,   K t be defined by Equation (166).
Then, we have the equation
M n + 1 t   =   h α n + 1 t ,  
which satisfies the condition M n t   > 0 for all n and all t > 0 .
The general binomial coefficients are defined by the equation
E k n K , M   =   E k n h 1 α , h α   =   α n α k ,  
where α   0 , 1 , k , n 0 , 0 k n .
Proof. 
Using the property
h α     h β t   = h α + β t
that holds for α , β > 0 , we obtain
M , n t   = ( h α ) , n   = h α n t .
Then,
M n + 1 t   =   K     M , n + 1 t   =   h 1 α     h α n + 1 t   = t α n Γ α n + 1   = h α n + 1 + 1 α t   = h α n + 1 t .
We see that the inequality M n + 1 t   >   0 holds for all n and for all t > 0 .
Then, using (169), we get
M n k + 1 t   M k + 1 t   = h α n k + 1 t   h α k + 1 t   = t α n k Γ α n k   +   1   t α k Γ α k + 1   = Γ α n + 1 Γ α k + 1   Γ α n k   +   1   t α n Γ α n + 1   =   α n α k   h α n + 1 t ,
where the generalized binomial coefficients α n α k are defined by (167).
As a result, using (170), we derive the equation
E k n K , M   = E k n h 1 α , h α   = M n k + 1 t   M k + 1 t M n + 1 t   =   α n α k ,
where α   0 , 1 , k , n 0 , 0 k n . □
Let us prove the following statement about the properties of the general binomial coefficients (171) for the kernels (166).
Theorem 15.
Let us consider the function
f n , k α   : =   α n α k   = Γ α n + 1 Γ α k + 1 Γ α n k   +   1 ,  
where α > 0 , 0 k n , k , n 0 . This function is increasing with respect to α > 0 with fixed parameters k , n 0 , and the following inequality holds
f n , k β     f n , k 1     f n , k α ,  
if β 1 α > 0 .
Proof. 
Let us prove that
d d α f n , k α     0 .
Using the equation
d d α Γ α   =   Γ α   ψ α ,
where ψ α is the digamma function of real argument α > 0 , we obtain
d d α f n , k α   = f n , k α   n   ψ n α + 1     k   ψ k α + 1     n k   ψ n k α + 1 .
Using n = m + k , Equation (172) can be written in the form
d d α f n , k α   = f m + k , k α   ( m   ψ m + k α + 1     ψ m α + 1   +   k   ψ m + k α + 1     ψ k α + 1 ) .
Using the digamma function ψ α of real argument α > 0 , increases function with respect to α > 0 :
ψ m + k α + 1     ψ m α + 1     0
for k 0 , and that the generalized binomial coefficients are positive functions
f n , k α   : =   α n α k   >   0
for α > 0 , 0 k n , k , n 0 , we derive that the generalized binomial coefficients f n , k α as a function of the variable α > 0 has a non-negative derivative d f n , k α / d α     0 with respect to α . □
Remark 7.
For general binomial coefficients
E k n K , M   =   α n α k ,  
we have the inequalities
n k     α n α k ,         ( 0 < α < 1 ) , α n α k     n k ,         ( α > 1 ) .  

4.3.7. Examples of Inequalities for General Binomial Coefficients

Let us give some examples of inequalities for general binomial coefficients.
Example 1.
Let us consider the pair of the kernels
M t   = h α t   = t α 1 Γ α ,         K t   = h 1 α t   = t α Γ 1 α .  
This pair belongs to the Sonin set S 1 , if α   0 , 1 . For these kernels, the general binomial coefficients are given by the equations
E k n K , M   = E k n h 1 α , h α   =   α n α k   = Γ α n + 1 Γ α k + 1   Γ α n k   +   1 .  
If the parameter α satisfies the condition α   0 , 1 , then the inequalities
n k     E k n K , M  
hold for all k , n 0 , when 0 k n and all t > 0 .
Example 2.
If the general binomial coefficients are given by equations
E k n K , M   =   α n α k  
with α   1 , 2 , then the inequalities
E k n K , M     n k  
hold for all k , n 0 , when 0 k n .
Note that the expressions of the general non-Markovian map Φ t M and the general binomial coefficients E k n K , M are derived for the kernels M t , K t that belongs to the Sonin set. Kernels (173) with α   1 , 2 do not belong to the Sonin set S 1 .
We can assume that within the framework of the GFC of arbitrary order [59,71], the expression of the general non-Markovian map Φ t M and the general binomial coefficients E k n K , M can also be derived.
For the kernels
M t   = h α t   = t α 1 Γ α , K t   = h 2 α t   = t 1 α Γ 2 α .
with α   1 , 2 , the non-Markovian equation for quantum observables has the form
D h 2 α t , τ A τ   = L   A t ,
where D h 2 α t τ is the Caputo fractional derivative of the order α   1 , 2 that is defined as
D h 2 α t τ A τ   : = I h α t τ   A 2 τ .
To solve Equation (175) and derive general binomial coefficients, we can use Theorem 4.3 and Example 4.10 of [29] (p. 231–231). The solution is described by the equation
A t   = A 0   E α , 1 t α L   +   A 1 0   t   E α , 2 t α L ,
where E α , β z is the two-parameter Mittag–Leffler function
E α , β t α L   =   n = 0 t α n Γ α n + β   L n
For A 0   = A and A 1 0   = 0 , we have the non-Markovian map
Φ t M A   =   E α , 1 t α L   A   =   n = 0 h α n + 1   L n A .
Corollary 4.
Let the general dissipativity condition
L n A   A     k = 0 n E k n h 2 α , h α   L n k A   L k ( A )
with the kernels (174), be satisfied for all n and A A .
Then, the general binomial coefficients have the form
E k n h 2 α , h α   =   α n α k       n k ,  
where α   1 , 2 .
For the non-Markovian map (176), the bi-positivity condition
Φ t M A   A     Φ t M A   Φ t M A
is satisfied for all t > 0 and all A A .

5. Non-Markovian Quantum Oscillator with Nonlocality in Time

Let us consider a non-Markovian quantum oscillator with nonlocality in time.
As is usually assumed for oscillators that are open quantum system [10,11,12,37], the general form of a bounded completely dissipative superoperator holds for an unbounded superoperator L . Then, the general non-Markovian dynamics of coordinate Q and momentum P is described by the equations
D K t , τ Q τ   =   L Q t ,
D K t , τ Q τ   =   L P t ,
where L is defined by Equation (20). For the linear quantum oscillator, the operators V k = V k Q , P and H = H Q , P are the functions of the coordinate and momentum operators in the form
H = 1 2 m P 2 + m ω 2 2 Q 2 + μ 2 Q P + P Q ,
V k = a k   P + b k   Q ,
where a k , and b k , k = 1 , 2 , are complex numbers. The term with parameter μ can be interpreted as friction, for which force is proportional to the velocity m 1 P .
Remark 8.
In general, the operators coordinate Q and momentum P are unbounded operators. Due to this, instead of the Hilbert space, one can use the so-called rigged Hilbert space (the Gelfand triplet). A rigged Hilbert space is the ordered triplet
B H = H B ,  
where H is a Hilbert space, B is a Banach space, and B is dual of B . The term “rigged Hilbert space” is also used to describe the dual pairs ( B , B ) generated from a Hilbert space H . The term “Gelfand triplet” is sometimes used instead of the term “rigged Hilbert space”. Example of a rigged Hilbert space is the triple of spaces that consists of the Banach space J n of test functions, the Hilbert space L 2 n of square integrable functions, and the Banach space J n of the linear functionals on J n . For details, see Chapter 2 in [5].
Remark 9.
The Lindblad result has been extended by E.B. Davies [9] to a class of quantum dynamical semi-group with unbounded generating superoperators.
Remark 10.
In this linear model, the parameters a k , b k cannot be arbitrary. Let us consider the real parameters
d a a   : =   2 | a 1 | 2   +   | a 2 | 2 ,         d b b   : =   2 | b 1 | 2   +   | b 2 | 2 ,          
d a b   : =   2 R e a 1 b 1 + a 2 b 2 ,       c a b =   2   I m a 1 b 1 + a 2 b 2 .
There is a fundamental constraint [12] on the parameters in the form
d a a   d b b d a b 2     c a b 2 ,  
which follows from the Schwartz inequality
( R e k a k   b k ) 2   +   ( I m k a k   b k ) 2 k | a k | 2   k | b l | 2 .
Using the canonical commutation relations for operators Q and P , we obtain Equations (180) and (181) for operators Q t and P t in the form
D K t , τ Q τ   = 1 m P t   +   μ λ Q t ,
D K t , τ P τ   =   m ω 2 Q t     μ + λ P t ,
where
λ = I m a 1 b 1 + a 2 b 2 ,
where D K t , is the general fractional derivative with the kernel K t     C 1 0 , , for which M t is the associated kernel so that the pair M t , K t belongs to the Sonin set S 1 .
For the kernel K t   = h 1 α t   = t α / Γ 1 α with 0 < α 1 , Equations (184) and (185) describe the exactly solvable model of non-Markovian dynamics that was first proposed in [37,38,39]. The exact solutions of equation for this case are derived in these works. For α = 1 , this model gives the standard Markovian quantum model [8,10,11,12].
Let us represent Equations (184) and (185) in the matrix form. Using the matrices
A t   =   Q t P t , N   =   μ λ m 1 m ω 2 μ λ ,
Equations (184) and (185) are written as
D K t , τ A τ   = N   A t ,
where we used L A t   = N   A t .
Theorem 16.
Let the function K t , M t   C 1 0 , belong to the Sonin set S 1 . Then, the initial value problem for the equation
D K t , τ A τ   = N   A t ,  
and the condition A 0   = A 0 , where A t and N are defined by (187), has the solution in the form
A t   = Φ t M A 0 ,
with the quantum dynamical map
Φ t M   =   L M , N , t   =   j = 1 K     M , j t   N j 1
where L M , N , t is the second Luchko function with the matrix argument N .
The statement of Theorem 16 follows directly from Theorem 3.
Example 3.
For the kernel K t   = h 1 α t   = t α / Γ 1 α with 0 < α 1 , the GFD is the Caputo fractional derivative of the order α   0 , 1 , and solution (191) has the form
Φ t h α   =   E α t α N   =   k = 0 t k α Γ α k + 1   N n .
For α = 1 , solution (192) gives
Φ t   =   e tN   =   n = 0 t n n !   N n
that describes the Markovian quantum dynamics of open system without nonlocality in time.
Let us prove the following theorem.
Theorem 17.
Equations (184) and (185) that describe the general non-Markovian quantum dynamics for coordinate and momentum have the solutions
Q t   =   C h M λ , ν , t   +   μ ν   S h M λ , ν , t Q 0 + 1 m ν   S h M λ , ν , t P 0 ,  
P t   = m ω 2 ν   S h M λ , ν , t Q 0   +   C h M λ , ν , t     μ ν   S h M λ , ν , t P 0 ,
where the functions S h M λ , ν , t and C h M λ , ν , t are defined by the expressions
S h M λ , ν , t   = 1 2 L M , λ + ν , t     L M , λ ν , t ,  
C h M λ , ν , t   = 1 2 L M , λ + ν , t   +   L M , λ ν , t ,
and ν is the complex parameter such that ν 2 = μ 2 ω 2 .
Proof. 
To get exact expression of the solution for coordinate and momentum operators, we represent the matrix N in the form
N = R   D   R 1 ,
where
D   =   λ + ν 0 0 λ ν ,
R   =   μ     ν c μ   +   ν c + m ω 2 c m ω 2 c + ,
and
c ± = μ ± ν 2 + ( m ω 2 ) 2 , ν = μ 2 ω 2 .
Using (198), the non-Markovian quantum dynamical map Φ t M is given as
Φ t M   =   L M , N , t   =   j = 1 K     M , j t   ( R   D   R 1 ) j 1   =   j = 1   K M , j t   R   D j 1   R 1   =   R   L M , D , t   R 1 .
As a result, we have
Φ t M = R   L M , D , t   R 1 .
Substituting expression (199) and (200) into Equation (203), we get the dynamical map
Φ t M α   =   Ch M λ , ν , t   +   μ / ν   Sh M λ , ν , t 1 / m ν   Sh M λ , ν , t m ω 2 / ν   Sh M λ , ν , t Ch M λ , ν , t     μ / ν   Sh M λ , ν , t ,
where we use the functions (196) and (197).  □
Remark 11.
For M t   = h α t , we get
Φ t h α = R   E α t α D   R 1 .  
For α = 1 , map (205) is given by the standard expression
Φ t h α = Φ t = R   e t D   R 1 .  
Theorem 17 describes the non-Markovian dynamics of coordinate and momentum of open quantum system (linear oscillator with friction) with general form on the nonlocality in time.
For λ = 0 , expressions (194) and (195) describe solutions of non-Markovian generalization of the Heisenberg equation, which is the equation for linear oscillator with Hamiltonian (182) and general form of nonlocality in time.
The solutions of Equations (184) and (185) are described by expressions (194) and (195). Depending on the specific form of the operator kernels M t , K t , the general trigonometric functions Ch M and Sh M , which are given by (196) and (197), will differ. Let us give examples.
Example 4.
For M t   = h α t , the solution of Equations (184) and (185), which describe non-Markovian dynamics of quantum system with power-law memory, were first derived in [37,38,39], where
Sh M λ , ν , t   = S α λ , ν , t   = 1 2 E α , 1 λ + ν t α     E α , 1 λ     ν t α ,  
Ch M λ , ν , t   = C α λ , ν , t   = 1 2 E α , 1 λ + ν t α   +   E α , 1 λ ν t α .
For the case with α = 1 , equations describe the Markovian dynamics of open quantum systems without nonlocality in time ( α = 1 ), since E 1 z   = e x p z , and
Sh M λ , ν , t   = e λ t s i n h ν t ,      Ch M λ , ν , t   = e λ t c o s h ν t ,
where s i n h and c o s h are hyperbolic sine and cosine.
Example 5.
For the kernel
M t   = h α t   e β t = t α 1 Γ α   e β t ,  
where t > 0 , α   0 , 1 , β > 0 , we have
Sh M λ , ν , t   = e β t S α λ , ν , t   +   β   0 t d τ   e β τ S α λ , ν , τ ,
Ch M λ , ν , t   = e β t C α λ , ν , t   +   β   0 t d τ   e β τ C α λ , ν , τ .
Example 6.
For the kernel
M t   = h 1 β + α t   +   h 1 β t ,  
where 0 < α < β < 1 , we get
Sh M λ , ν , t   = 1 2 ( E 1 β , 1 β + α , 1 λ + ν t 1 β ,   λ + ν t 1 β + α   E 1 β , 1 β + α , 1 λ ν t 1 β ,   λ ν t 1 β + α ) ,  
Ch M λ , ν , t   = 1 2 ( E 1 β , 1 β + α , 1 λ + ν t 1 β , λ + ν t 1 β + α   + E 1 β , 1 β + α , 1 λ ν t 1 β , λ ν t 1 β + α ) ,  
where E 1 β , 1 β + α , 0 is the multinomial Mittag–Leffler function (58).

6. Non-Markovian Quantum Dynamics of Two-Level System

In the case of an N -level open quantum system, the problem was investigated by V. Gorini, A. Kossakowski, and E.C.G. Sudarshan [6]. The general form of the generating superoperator of a completely positive dynamical semi-group of this system has been established [6]. In the case of the N -level quantum system, the Hilbert space H has the dimension d i m H = N . Each N -dimensional separable Hilbert space over is isomorphic to n .
Let us consider a general non-Markovian dynamics of quantum states. The non-Markovian dynamics of the density operator ρ t can be described by the equations with GFD in the form
D K t , τ ρ τ   =   1 i H , ρ t   +   1 k = 1 V k ρ t V k     1 2 V k   V k , ρ t ,
where A , B   = A   B + B   A .
Let us consider the non-Markovian two-level quantum systems with general nonlocality in time, which is described by the Sonin kernel K t . The Hamiltonian will be considered in the form
H = 1 2     ω 0   σ 3 = 1 2     ω 0   1 1     0 0 ,
where ω 0 > 0 is the transition frequency. Then, the Hamiltonian of the two-level quantum system is diagonal in the basis | 0 , | 1 .
The operators V k 0 with k = 1 , 2 will be considered in the form
V 1 = η 1   σ , V 2 = η 2   σ + ,
V 1 = η 1   σ + , V 2 = η 2   σ ,
where η 1 , η 2 , and
σ +   =   1 0   = 1 2 σ 1 + i   σ 2 , σ   =   0 1   = 1 2 σ 1 i   σ 2 .
For operators (217), (218), (219), the general non-Markovian master Equation (216) has the form
D K t , τ ρ τ   = i H , ρ   +   γ 1 σ ρ σ + 1 2 σ + σ , ρ   +   γ 2 σ + ρ σ 1 2 σ σ + , ρ ,
where
γ 1 = 1 η 1 2 , γ 2 = 1 η 2 2 .
We can consider the parameters γ 1 and γ 2 in the form
γ 1 = γ 0   N 0 + 1 , γ 2 = γ 0   N 0 , γ 0 = 4 ω 0 3 | d | 2 3 c 3 ,
where γ 0 is the spontaneous emission rate. Then, Equation (221) with K t   = δ t takes the form of Equation (3.219) in [3] (p. 148). In non-Markovian master Equation (221), the terms with γ 1 and γ 2 describe spontaneous emission with rate γ 0 , the thermally induced emission and absorption with the rate γ 0   N 0 . The total transition rate is equal to
γ = γ 1 + γ 2 = γ 0   2 N 0 + 1 ,
where N 0 = N 0 ω 0 is the Planck distribution at the transition frequency.
Using the density operator of two-level system in the form
ρ = k , l = 0 1 ρ k l   k l ,
the general non-Markovian equation of this two-level quantum system is represented by the equations with GFD:
D K t , τ ρ 00 τ   = γ 2 ρ 00 + γ 1 ρ 11
D K t , τ ρ 01 τ   = i ω 0 ρ 01 1 2 γ 1 + γ 2 ρ 01
D K t , τ ρ 10 τ , = i ω 0 ρ 10 1 2 γ 1 + γ 2 ρ 10
D K t , τ ρ 11 τ   = γ 2 ρ 00 γ 1 ρ 11
The interaction representation of quantum theory cannot be used to solve equations with GFD, which describe non-Markovian dynamics. This fact follows from the violation of the standard Leibniz rule (product rule) for GFD. Note that the time-ordered product (chronological product) for non-Markovian quantum dynamics with nonlocality in type (in the form of memory) was proposed in [45].
To obtain an explicit form of the solution for the components of the density operator of two-level system, Equations (226)–(229) are considered by two pairs.
First pair is the equations
D K t , τ ρ 01 τ   = λ 1   ρ 01 ,
D K t , τ ρ 10 τ   = λ 2   ρ 10 ,
where
λ 1 = + i ω 0 1 2 γ 1 + γ 2 , λ 2 = i ω 0 1 2 γ 1 + γ 2 .
For 0 < α 1 , the solutions of Equations (230) and (231) have the form
ρ 01 t   = L M , λ 1 , t   ρ 01 0 ,
ρ 10 t   = L M , λ 2 , t   ρ 10 0 .
The second pair is the equations
D K t , τ ρ 00 τ   = γ 2 ρ 00 + γ 1 ρ 11 ,
D K t , τ ρ 11 τ   = γ 2 ρ 00 γ 1 ρ 11 .
To get solutions of Equations (235) and (236), we considered these equations in the matrix form
D K t , τ ρ 00 τ ρ 11 τ   =   G   ρ 00 ρ 11 ,
where
G   =   γ 2 γ 1 γ 2 γ 1 ,
and γ 1 > 0 , and γ 2 > 0 .
The solution of Equation (237) has the form
ρ 00 t ρ 11 t   = S t α   ρ 00 0 ρ 11 0 ,
where
S t α   = L M , G , t   +   S γ 1 , γ 2 ,
S γ 1 , γ 2   =   γ 1 / γ γ 1 / γ γ 2 / γ γ 1 / γ
with γ = γ 1 + γ 2 .
The matrix G can be diagonalized as
G   =     K 1   D 1   K 1 1 ,
where
D 1   =   0 0 0 γ 1 γ 2 ,
K 1   =   γ 1 γ 1 2   +   γ 2 2 1 2 γ 2 γ 1 2   +   γ 2 2 1 2 ,
K 1 1   =   γ 1 2   +   γ 2 2 γ 1   +   γ 2 γ 1 2   +   γ 2 2 γ 1   +   γ 2 2 γ 2 γ 1   +   γ 2 2 γ 1 γ 1   +   γ 2 .
Exact expression for solution (239) of Equation (237) is derived by the transformations
L M , K 1   D 1   K 1 1 , t   =   K 1   L M , D 1 , t   K 1 1 =
γ 1 γ 1 2   +   γ 2 2 1 2 γ 2 γ 1 2   +   γ 2 2 1 2   0 0 0 L M , γ , t   γ 1 2   +   γ 2 2 γ 1   +   γ 2 γ 1 2   +   γ 2 2 γ 1   +   γ 2 2 γ 2 γ 1   +   γ 2 2 γ 1 γ 1   +   γ 2   =
γ 2 γ   L M , γ , t γ 1 γ   L M , γ , t γ 2 γ   L M , γ , t γ 1 γ   L M , γ , t .
Therefore, we get
S t α   = L M , G , t   +   S γ 1 , γ 2   =   γ 1 γ + γ 2 γ   L M , γ , t γ 1 γ γ 1 γ   L M , γ , t γ 2 γ γ 2 γ   L M , γ , t γ 2 γ + γ 1 γ   L M , γ , t ,
where γ = γ 1 + γ 2 , and the solution for ρ 00 t and ρ 11 t has in the form
ρ 00 t ρ 11 t   =   γ 1 γ + γ 2 γ   L M , γ , t γ 1 γ γ 1 γ   L M , γ , t γ 2 γ γ 2 γ   L M , γ , t γ 2 γ + γ 1 γ   L M , γ , t ·   ρ 00 0 ρ 11 0 .
As a result, we obtain the solution for components of the density matrix ρ k l t in the form
ρ 01 t   = L M , λ 2 , t   ρ 01 0 ,
ρ 10 t   = L M , λ 1 , t   ρ 10 0 ,
and
ρ 00 t   =   γ 1 γ + γ 2 γ L M , γ , t   ρ 00 0   +   γ 1 γ γ 1 γ L M , γ , t   ρ 11 0 ,
ρ 11 t   =   γ 1 γ γ 2 γ L M , γ , t   ρ 00 0   +   γ 1 γ + γ 1 γ L M , γ , t   ρ 11 0 ,
where γ = γ 1 + γ 2 .

7. Entropy for General Non-Markovian Quantum Dynamics

In quantum mechanics and quantum statistics, the concept of entropy S is defined through the density operator ρ , which is a positive, normalized self-adjoint linear operator. John von Neumann defines [82] the entropy as an extension of the Gibbs entropy concepts from classical mechanics to the quantum mechanics.
For a quantum-mechanical system described by a density operator, the von Neumann entropy is defined (see Sections V.2 and V.3 in [82]) by the equation
S t   = T r ρ t   ln ρ t ,
where T r is the trace and ln denotes the (natural) matrix logarithm.
For Markovian dynamics of Hamiltonian quantum systems, the von Neumann entropy does not change
S t   = T r ρ t   ln ρ t   = T r ρ 0   ln ρ 0   = S 0 .
For non-Markovian and non-Hamiltonian systems, the von Neumann entropy changes in the general case.
Let us consider the von Neumann entropy for general non-Markovian dynamics of two-level quantum systems.
For two-level quantum systems, the density operator is given as
ρ t   =   ρ 00 t ρ 01 t ρ 10 t ρ 11 t .
Matrix (257) can be diagonalized
ρ t   = K 2   D 2   K 2 1 .
The diagonal matrix D 2 has the form
D 2   =   D 0 0 D + ,
where
D ± = 1 2 ρ 00 t   +   ρ 11 t   ±   ( ρ 00 t     ρ 11 t ) 2 + 4 ρ 01 t ρ 10 t .
Using the normalization condition T r ρ t   = 1 for all t   0 , , we get
ρ 00 t   +   ρ 11 t   =   1 .
Using representation (258) of density operator (257), we can derive an explicit form of the von Neumann entropy for non-Markovian two-level quantum system.
The von Neumann entropy has the form
S t   = T r D 2   ln D 2 ,
where we use
T r ρ t ln ρ t   = T r K 2   D 2   K 2 1 ln K 2   D 2   K 2 1   = T r K 2   D 2   K 2 1 K 2   ln D 2   K 2 1   = T r K 2   D 2   ln D 2   K 2 1   = T r K 2 1   K 2   D 2   ln D 2   = T r D 2   ln D 2 .
The von Neumann entropy is given by the equation
S t   = 1 2 1 B t   ln 1 2 1 B t     1 2 1 + B t   ln 1 2 1 + B t ,
where B t is the Bloch vector
B t   = ( ρ 00 t     ρ 11 t ) 2 + 4   ρ 01 t   ρ 10 t .
Using solutions (251), (252), and (253), (254), we derive B t in the form
B 2 t   = 2 γ 2 γ L M , γ , t   ρ 00 0     2 γ 1 γ L M , γ , t   ρ 11 0   +   γ 1 γ 2 γ 2   + 4   L M , λ 2 , t   L M , λ 1 , t   ρ 01 0   ρ 10 0 ,
where γ = γ 1 + γ 2 and ρ 00 t   +   ρ 11 t   = 1 for all t 0 since T r ρ t   = 1 .
Let us give some examples of the function L M , γ , t .
(1)
For M t   = h α t   = t α 1 / Γ α , the second Luchko function has the form
L M , γ , t   = E α γ t α .
(2)
For the kernel
M t   = t α 1 Γ α   e β t ,
we have
L M , γ , t   =   e β t   E α γ t α   +   β   0 t d τ   e β τ E α γ τ α ,
where α   0 , 1 and β > 0 .
(3)
For the kernel
M t   =   t β + α Γ 1 β + α   +   t β Γ 1 β ,
the function is
L M , γ , t   =   E 1 β , 1 β + α , 1 γ t 1 β , γ t 1 β + α ,
where 0 < α < β < 1 .

8. Conclusions

In this paper, we proposed the formulation of non-Markovian quantum theory in the general form. The non-locality in time is represented by kernels of integral and integro-differential operators. These kernels are described by functions that belong to the Sonin set of kernel pairs. The results can be derived in the general form without using special realization of these kernels. Therefore, these results are valid for any operator kernels from the Sonin set. This approach to non-Markovian quantum theory is directly connected with the concept of general fractional dynamics suggested in [72].
Non-Markovian equations for quantum observables and states are suggested by using general fractional calculus. The exact solutions of these equations are derived by using the operational calculus, which is proposed by Luchko in [60] for equations with general fractional derivatives. A wide class of the exactly solvable models of non-Markovian quantum dynamics is suggested. These models describe non-Markovian open quantum systems with the general form of nonlocality in time. The non-Markovian models of quantum oscillator and two-level quantum system with general form of nonlocality in time are described. The exact solutions of equations for these models are proposed.
This paper proposes a general approach to describing non-Markov quantum dynamics. Many important issues are not covered in this article. This work does not offer solutions to all the problems of constructing general non-Markov dynamics of open quantum systems. Let us note some unresolved questions that await their solution in future research.
(1)
A quantum system can be embedded in some environments and therefore the system is not isolated. The environment of a quantum system is in principle unobservable or is unknown. This would render the non-Markovian theory of open quantum systems a fundamental generalization of quantum mechanics. However, for practical applications, it is useful to have models of open quantum systems that can be derived from some closed systems including the system under study and some environments. In this regard, the problem arises of constructing models of such closed systems and obtaining general non-Markov dynamics, for example, within the framework of the Caldeira–Leggett approach [83]. At the moment, this problem has not been solved, and the question remains open. We think that the construction of such models is possible. This opinion is based on the following: in the framework of the simplest models, the kernels of fractional derivatives, which describe nonlocality in time, were obtained in [84].
(2)
For open quantum systems, its “reduced” dynamics not to violate thermodynamics must not decrease entropy of the evolving state [85]. In this regard, the problem arises of a detailed study of the behavior of entropy for general non-Markov dynamics. At the moment, this problem has not been solved. This question is interesting for further research and computer simulation of the behavior of entropy.
(3)
The form of the superoperator L was determined by the Lindblad theorem, which describes the relationship between a completely positive semigroup and a completely dissipative superoperator. The condition for the dissipativity of the superoperator is in fact the standard Leibniz rule, in which equality is replaced by inequality. In non-Markovian dynamics, the semigroup property is violated, and the fractional derivative violates the standard Leibniz rule. In this regard, the question arises about the existence of a generalization of the Lindblad superoperator, in the framework of the proposed general non-Markovian dynamics. In our opinion, such a possibility exists and may be associated with the fractional powers of Lindblad superoperators and the models proposed in the works [5] (pp. 433–444) and [37] (pp. 458–464, 468–477), and [41,42].

Funding

This research received no external funding.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Davies, E.B. Quantum Theory of Open Systems; Academic Press: Waltham, MA, USA, 1976; p. 179. ISBN 0-12-206150-0. [Google Scholar]
  2. Ingarden, R.S.; Kossakowski, A.; Ohya, M. Information Dynamics and Open Systems: Classical and Quantum Approach; Kluwer: New York, NY, USA, 1997; p. 310. [Google Scholar] [CrossRef]
  3. Breuer, H.-P.; Petruccione, F. Theory of Open Quantum Systems; Oxford University Press: Oxford, UK, 2002; p. 625. ISBN 0-19-852063-8. [Google Scholar]
  4. Attal, A.; Joye, A.; Pillet, C.A. (Eds.) Open Quantum Systems: The Markovian Approach; Springer: Berlin/Heidelberg, Germany, 2006; p. 244. [Google Scholar] [CrossRef]
  5. Tarasov, V.E. Quantum Mechanics of Non-Hamiltonian and Dissipative Systems; Elsevier: Amsterdam, The Netherlands, 2008; p. 540. ISBN 9780444530912. [Google Scholar]
  6. 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]
  7. Lindblad, G. On the generators of quantum dynamical semigroups. Commun. Math. Phys. 1976, 48, 119–130. [Google Scholar] [CrossRef]
  8. Lindblad, G. Brownian motion of a quantum harmonic oscillator. Rep. Math. Phys. 1976, 10, 393–406. [Google Scholar] [CrossRef]
  9. Davies, E.B. Quantum dynamical semigroups and neutron diffusion equation. Rep. Math. Phys. 1977, 11, 169–188. [Google Scholar] [CrossRef]
  10. Isar, A.; Sandulescu, A.; Scutaru, H.; Stefanescu, E.; Scheid, W. Open quantum systems. Int. J. Mod. Phys. E 1994, 3, 635–714. [Google Scholar] [CrossRef] [Green Version]
  11. Isar, A.; Sandulescu, A.; Scheid, W. Phase space representation for open quantum systems with the Lindblad theory. Int. J. Mod. Phys. B 1996, 10, 2767–2779. [Google Scholar] [CrossRef] [Green Version]
  12. Sandulescu, A.; Scutaru, H. Open quantum systems and the damping of collective models in deep inelastic collisions. Ann. Phys. 1987, 173, 277–317. [Google Scholar] [CrossRef]
  13. Tarasov, V.E. Uncertainty relation for non-Hamiltonian quantum systems. J. Math. Phys. 2013, 54, 13. [Google Scholar] [CrossRef] [Green Version]
  14. Tarasov, V.E. Path integral for quantum operations. J. Phys. A 2004, 37, 3241–3257. [Google Scholar] [CrossRef] [Green Version]
  15. Tarasov, V.E. Pure stationary states of open quantum systems. Phys. Rev. E 2002, 66, 056116. [Google Scholar] [CrossRef] [Green Version]
  16. Tarasov, V.E. Stationary states of dissipative quantum systems. Phys. Lett. A 2002, 299, 173–178. [Google Scholar] [CrossRef] [Green Version]
  17. Tarasov, V.E. Quantum computer with mixed states and four-valued logic. J. Phys. A 2002, 35, 5207–5235. [Google Scholar] [CrossRef] [Green Version]
  18. Chruscinski, D.; Kossakowski, A. From Markovian semigroup to non-Markovian quantum evolution. Europhys. Lett. 2012, 97, 20005. [Google Scholar] [CrossRef] [Green Version]
  19. Vacchini, B.; Smirne, A.; Laine, E.-M.; Piilo, J.; Breuer, H.-P. Markovianity and non-Markovianity in quantum and classical systems. New J. Phys. 2011, 13, 26. [Google Scholar] [CrossRef]
  20. Breuer, H.-P. Foundations and measures of quantum non-Markovianity. J. Phys. B At. Mol. Opt. Phys. 2012, 45, 154001. [Google Scholar] [CrossRef] [Green Version]
  21. Rivas, A.; Huelga, S.F.; Plenio, M.B. Quantum non-Markovianity: Characterization, quantification and detection. Rep. Prog. Phys. 2014, 77, 094001. [Google Scholar] [CrossRef] [PubMed]
  22. Breuer, H.-P.; Laine, E.-M.; Piilo, J.; Vacchini, B. Non-Markovian dynamics in open quantum systems. Rev. Mod. Phys. 2016, 88, 021002. [Google Scholar] [CrossRef] [Green Version]
  23. De Vega, I.; Alonso, D. Dynamics of non-Markovian open quantum systems. Rev. Mod. Phys. 2017, 89, 015001. [Google Scholar] [CrossRef] [Green Version]
  24. Li, L.; Hall, M.J.W.; Wiseman, H.M. Concepts of quantum non-Markovianity: A hierarchy. Phys. Rep. 2018, 759, 1–51. [Google Scholar] [CrossRef] [Green Version]
  25. Megier, N.; Smirne, A.; Vacchini, B. The interplay between local and non-local master equations: Exact and approximated dynamics. New J. Phys. 2020, 22, 083011. [Google Scholar] [CrossRef]
  26. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives Theory and Applications; Gordon and Breach: New York, NY, USA, 1993; p. 1006. ISBN 9782881248641. [Google Scholar]
  27. Kiryakova, V. Generalized Fractional Calculus and Applications; Longman and John Wiley: New York, NY, USA, 1994; p. 360. ISBN 9780582219779. [Google Scholar]
  28. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA, 1998; p. 340. [Google Scholar]
  29. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006; p. 540. [Google Scholar]
  30. Diethelm, K. The Analysis of Fractional Differential Equations: An Application-Oriented Exposition Using Differential Operators of Caputo Type; Springer: Berlin/Heidelberg, Germany, 2010; p. 247. [Google Scholar] [CrossRef] [Green Version]
  31. Kochubei, A.; Luchko, Y. Handbook of Fractional Calculus with Applications. Volume 1 Basic Theory; Walter de Gruyter: Berlin, Germany, 2019; p. 481. [Google Scholar] [CrossRef]
  32. Kochubei, A.; Luchko, Y. Handbook of Fractional Calculus with Applications. Volume 2 Fractional Differential Equations; Walter de Gruyter: Berlin, Germany, 2019; p. 519. [Google Scholar] [CrossRef] [Green Version]
  33. Tarasov, V.E. Handbook of Fractional Calculus with Applications. Volume 4 Application in Physics. Part A; Walter de Gruyter: Berlin, Germany, 2019; p. 306. [Google Scholar] [CrossRef]
  34. Tarasov, V.E. Handbook of Fractional Calculus with Applications. Volume 5. Application in Physics. Part B; Walter de Gruyter: Berlin, Germany, 2019; p. 319. [Google Scholar] [CrossRef] [Green Version]
  35. Tarasov, V.E. No nonlocality. No fractional derivative. Commun. Nonlinear Sci. Numer. Simul. 2018, 62, 157–163. [Google Scholar] [CrossRef] [Green Version]
  36. Tarasov, V.E. Rules for fractional-dynamic generalizations: Difficulties of constructing fractional dynamic models. Mathematics 2019, 7, 554. [Google Scholar] [CrossRef] [Green Version]
  37. Tarasov, V.E. Fractional Dynamics: Applications of Fractional Calculus to Dynamics of Particles, Fields and Media; Springer: Berlin/Heidelberg, Germany, 2010; p. 516. [Google Scholar] [CrossRef] [Green Version]
  38. Tarasov, V.E. Quantum dissipation from power-law memory. Ann. Phys. 2012, 327, 1719–1729. [Google Scholar] [CrossRef]
  39. Tarasov, V.E. Fractional quantum mechanics of open quantum systems. In Handbook of Fractional Calculus with Applications. Volume 5: Applications in Physics, Part B; Walter de Gruyter: Berlin, Germany, 2019; Chapter 11; pp. 257–277. [Google Scholar] [CrossRef]
  40. Tarasov, V.E. Quantum maps with memory from generalized Lindblad equation. Entropy 2021, 23, 544. [Google Scholar] [CrossRef] [PubMed]
  41. Tarasov, V.E. Fractional generalization of the quantum Markovian master equation. Theor. Math. Physi. 2009, 158, 179–195. [Google Scholar] [CrossRef] [Green Version]
  42. Tarasov, V.E. Fractional Dynamics of Open Quantum Systems. In Fractional Dynamics: Recent Advances; Klafter, J., Lim, S.C., Metzler, R., Eds.; World Scientific: Singapore, 2012; Chapter 19; pp. 449–482. [Google Scholar] [CrossRef]
  43. Tarasov, V.E. Fractional Heisenberg equation. Phys. Lett. A 2008, 372, 2984–2988. [Google Scholar] [CrossRef] [Green Version]
  44. Tarasov, V.E. Fractional diffusion equations for open quantum systems. Nonlinear Dyn. 2013, 71, 663–670. [Google Scholar] [CrossRef]
  45. Tarasov, V.E.; Tarasova, V.V. Time-dependent fractional dynamics with memory in quantum and economic physics. Ann. Phys. 2017, 383, 579–599. [Google Scholar] [CrossRef]
  46. Tarasov, V.E. Fractional dynamics of relativistic particle. Int. J. Theor. Phys. 2010, 49, 293–303. [Google Scholar] [CrossRef] [Green Version]
  47. Tarasov, V.E. Dirac particle with memory: Proper time non-locality. Phys. Lett. A. 2020, 384, 8. [Google Scholar] [CrossRef] [Green Version]
  48. Iomin, A. Fractional-time quantum dynamics. Phys. Rev. E 2009, 80, 022103. [Google Scholar] [CrossRef] [Green Version]
  49. Iomin, A. Fractional time quantum mechanics. In Handbook of Fractional Calculus with Applications. Volume 5. Applications in Physics, Part B; Tarasov, V.E., Ed.; De Gruyter: Berlin, Germany, 2019; pp. 299–315. [Google Scholar] [CrossRef]
  50. Iomin, A.; Mendez, V.; Horsthemke, W. Comb model: Non-Markovian versus Markovian. Fractal Fract. 2019, 3, 54. [Google Scholar] [CrossRef] [Green Version]
  51. Kochubei, A.N. General fractional calculus, evolution equations, and renewal processes. Integral Equ. Oper. Theory 2011, 71, 583–600. [Google Scholar] [CrossRef] [Green Version]
  52. Kochubei, A.N. General fractional calculus. In Handbook of Fractional Calculus with Applications. Volume 1. Basic Theory; Kochubei, A., Luchko, Y., Eds.; De Gruyter: Berlin, Germany, 2019; Chapter 5; pp. 111–126. [Google Scholar] [CrossRef]
  53. Kochubei, A.N. Equations with general fractional time derivatives. Cauchy problem. In Handbook of Fractional Calculus with Applications. Volume 2. Fractional Differential Equations; Tenreiro Machado, J.A., Ed.; De Gruyter: Berlin, Germany, 2019; Chapter 11; pp. 223–234. [Google Scholar] [CrossRef]
  54. Kochubei, A.N.; Kondratiev, Y. Growth equation of the general fractional calculus. Mathematics 2019, 7, 615. [Google Scholar] [CrossRef] [Green Version]
  55. Sonine, N. On the generalization of an Abel formula. (Sur la generalisation d’une formule d’Abel). Acta Math. 1884, 4, 171–176. (In French) [Google Scholar] [CrossRef]
  56. Sonin, N.Y. On the generalization of an Abel formula. In Investigations of Cylinder Functions and Special Polynomials; Sonin, N.Y., Ed.; GTTI: Moscow, Russia, 1954; pp. 148–154. [Google Scholar]
  57. Sonin, N.Y. Math-Net.Ru. 2021. Available online: www.mathnet.ru/eng/person40324 (accessed on 12 June 2021).
  58. Luchko, Y. General fractional integrals and derivatives with the Sonine kernels. Mathematics 2021, 9, 594. [Google Scholar] [CrossRef]
  59. Luchko, Y. General fractional integrals and derivatives of arbitrary order. Symmetry 2021, 13, 755. [Google Scholar] [CrossRef]
  60. Luchko, Y. Operational calculus for the general fractional derivatives with the Sonine kernels. Fract. Calc. Appl. Anal. 2021, 24, 338–375. [Google Scholar] [CrossRef]
  61. Kochubei, A.N.; Kondratiev, Y. Fractional kinetic hierarchies and intermittency. Kinetic and related models. Am. Inst. Math. Sci. 2017, 10, 725–740. [Google Scholar] [CrossRef]
  62. Sin, C.-S. Well-posedness of general Caputo-type fractional differential equations. Fract. Calc. Appl. Anal. 2018, 21, 819–832. [Google Scholar] [CrossRef]
  63. Luchko, Y.; Yamamoto, M. General time-fractional diffusion equation: Some uniqueness and existence results for the initial-boundary-value problems. Fract. Calc. Appl. Anal. 2016, 19, 675–695. [Google Scholar] [CrossRef]
  64. Luchko, Y.; Yamamoto, M. The general fractional derivative and related fractional differential equations. Mathematics 2020, 8, 2115. [Google Scholar] [CrossRef]
  65. Samko, S.G.; Cardoso, R.P. Integral equations of the first kind of Sonine type. Int. J. Math. Math. Sci. 2003, 203, 3609–3632. [Google Scholar] [CrossRef] [Green Version]
  66. Samko, S.G.; Cardoso, R.P. Sonine integral equations of the first kind in L~p (0, b). Fract. Calc. Appl. Anal. 2003, 6, 235–258. [Google Scholar]
  67. Kinash, N.; Janno, J. Inverse problems for a generalized subdiffusion equation with final overdetermination. Math. Model. Anal. 2019, 24, 236–262. [Google Scholar] [CrossRef]
  68. Kinash, N.; Janno, J. An inverse problem for a generalized fractional derivative with an application in reconstruction of time- and space-dependent sources in fractional diffusion and wave equations. Mathematics 2019, 7, 1138. [Google Scholar] [CrossRef] [Green Version]
  69. Hanyga, A. A comment on a controversial issue: A generalized fractional derivative cannot have a regular kernel. Fract. Calc. Appl. Anal. 2020, 23, 211–223. [Google Scholar] [CrossRef] [Green Version]
  70. Giusti, A. General fractional calculus and Prabhakar’s theory. Commun. Nonlinear Sci. Numer. Simul. 2020, 83, 105114. [Google Scholar] [CrossRef] [Green Version]
  71. Tarasov, V.E. General fractional calculus: Multi-kernel approach. Mathematics 2021, 9, 1501. [Google Scholar] [CrossRef]
  72. Tarasov, V.E. General fractional dynamics. Mathematics 2021, 9, 1464. [Google Scholar] [CrossRef]
  73. Tarasov, V.E. Caputo-Fabrizio operator in terms of integer derivatives: Memory or distributed lag? Comput. Appl. Math. 2019, 38, 15. [Google Scholar] [CrossRef]
  74. Tarasov, V.E.; Tarasova, V.V. Logistic equation with continuously distributed lag and application in economics. Nonlinear Dyn. 2019, 97, 1313–1328. [Google Scholar] [CrossRef]
  75. Luchko, Y.; Gorenflo, R. An operational method for solving fractional differential equations. Acta Math. Vietnam 1999, 24, 207–234. [Google Scholar]
  76. Gorenflo, R.; Kilbas, A.A.; Mainardi, F.; Rogosin, S.V. Mittag-Leffler Functions, Related Topics and Applications; Springer: Berlin/Heidelberg, Germany, 2020; p. 540. [Google Scholar] [CrossRef]
  77. Hadid, S.B.; Luchko, Y. An operational method for solving fractional differential equations of an arbitrary real order. Panam. Math. J. 1996, 6, 57–73. [Google Scholar]
  78. Peng, J.; Li, K. A note on property of the Mittag-Leffler function. J. Math. Anal. Appl. 2010, 370, 635–638. [Google Scholar] [CrossRef] [Green Version]
  79. Elagan, S.K. On the invalidity of semigroup property for the Mittag-Leffler function with two parameters. J. Egypt. Math. Soc. 2016, 24, 200–203. [Google Scholar] [CrossRef]
  80. Sadeghi, A.; Cardoso, J.R. Some notes on properties of the matrix Mittag-Leffler function. Appl. Math. Comput. 2018, 338, 733–738. [Google Scholar] [CrossRef]
  81. Tarasov, V.E. No violation of the Leibniz rule. No fractional derivative. Commun. Nonlinear Sci. Numer. Simul. 2013, 18, 2945–2948. [Google Scholar] [CrossRef] [Green Version]
  82. Von Neumann, J. Mathematical Foundations of Quantum Mechanics; Princeton University Press: Princeton, NJ, USA, 2018; p. 328. ISBN 978-0-691-7856-1. [Google Scholar]
  83. Caldeira, A.O.; Leggett, A.J. Path integral approach to quantum Brownian motion. Phys. A 1983, 121, 587–616. [Google Scholar] [CrossRef]
  84. Tarasov, V.E. The fractional oscillator as an open system. Cent. Eur. J. Phys. 2012, 10, 382–389. [Google Scholar] [CrossRef] [Green Version]
  85. Lindblad, G. Non-Equilibrium Entropy and Irreversibility; Springer: Dordrecht, The Netherlands, 1983; p. 166. ISBN 978-90-277-1640-8. [Google Scholar]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tarasov, V.E. General Non-Markovian Quantum Dynamics. Entropy 2021, 23, 1006. https://0-doi-org.brum.beds.ac.uk/10.3390/e23081006

AMA Style

Tarasov VE. General Non-Markovian Quantum Dynamics. Entropy. 2021; 23(8):1006. https://0-doi-org.brum.beds.ac.uk/10.3390/e23081006

Chicago/Turabian Style

Tarasov, Vasily E. 2021. "General Non-Markovian Quantum Dynamics" Entropy 23, no. 8: 1006. https://0-doi-org.brum.beds.ac.uk/10.3390/e23081006

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