Fractional Langevin Equation in Quantum Systems with Memory Effect

In this paper, we introduce the fractional generalized Langevin equation (FGLE) in quantum systems with memory effect. For a particular form of memory kernel that characterizes the quantum system, we obtain the analytical solution of the FGLE in terms of the two-parameter Mittag-Leffler function. Based on this solution, we study the time evolution of this system including the qubit excited-state energy, polarization and von Neumann entropy. Memory effect of this system is observed directly through the trapping states of these dynamics.


Introduction
Application of the fractional calculus on physics has attracted an increasing attention this decade including physical kinetics [1], anomalous transport theory in solid-state physics [2], and nonlinear dynamics [3].These applications aim to explore the nonlocal quantum phenomena found for either long-range interactions or time-dependent processes with long-time memory effect.Fractional space and time derivatives are used to describe these systems with nonlocal dynamics, e.g.anomalous diffusion or fractional Brownian motion [4] [5].Recently, the fractional time derivative is given a physical interpretation by Iomin which describes an effective interaction of a quantum system with its environment [6].Here we will introduce the practical quantum system with long-time memory effect [7].
Quantum information systems are currently attracting enormous interest for their fundamental nature and potential applications to computation and secure communication.Many novel methods have been proposed to generate controllable qubit states through the environment of the qubits [8]- [10].When a qubit made of a twolevel atom is embedded inside a structured reservoir with memory effect, the correlation between the qubit and environment will affect the dynamics of the qubit [7] [11] [12].In this case, the effective interaction of the qubit with the photon modes of the environment is expressed as a memory kernel.In this paper, we take a particular memory kernel of an anisotropic photonic crystal (PhC) as an example to illustrate the method and suitability of applying fractional calculus to this qubit system.We find that the kinetic equation can be expressed as a fractional generalized Langevin equation (FGLE) when the fractional time derivative is used to express memory kernel term.The expression of FGLE opens a new route for the application of fractional calculus to quantum information systems.
This paper is organized as follows.In Section 2, we present a FGLE for the quantum system with memory effect.The general solution of this FGLE is expressed in terms of two-parameter Mittag-Leffler function through the methodology of Laplace transform.In Section 3, a particular memory kernel characterizing the quantum system of a qubit in an anisotropic PhC is used to illustrate the solving procedures of the FGLE.Based on the analytical solution of this FGLE, we study the dynamics of the qubit energy, polarization and von Neumann entropy.Finally, we summarize our results in Section 4.

Fractional Langevin Equation
When the quantum system with memory effect is considered, the general form of the kinetic equation can be derived from the time dependent Schrődinger equation as where ( ) A t denotes the time evolution of the quantum system and ( ) the memory kernel of the reservoir.This kinetic equation reads that the future of the system is determined by the memory of the reservoir in its previous state.
To solve Equation (1.1), we introduce the Laplace transform.When the Laplace image of the memory kernel ( ) K t is considered, we need to construct a Cantor set with infinite number of divisions of time interval because no characteristic time scale exists in this memory kernel.As the Cantor set is constructed, we could choose the entire time interval as T with unity height.When the central part of the time interval is removed, the interval leaves two intervals with length T ξ ( ) . In order to keep the integral memory, the heights of the two intervals must be increased from unity to the value ( ) 2 T ξ .In the next stage, each remaining interval with length T ξ is subjected to the same division process.As the division process is performed n times, the memory kernel in Laplace image will be represented by a set of 2 n intervals with length n T ξ and height ( ) T ξ .The integral of the memory kernel in Laplace image is approached as the division time n is taken to be infinity.This approach leads to the Laplace transform of the memory kernel has the form of ( ) ( ) ( ) with ν being the fractal dimension of time blocks participating in the construct of the Cantor set.
As the inverse Laplace transform is performed on the kernel ( ) , we obtain the memory kernel approximately as with Gamma function ( ) . Substituting this memory kernel into Equation (1.1), we have the kinetic equation as Comparing the right-hand-side term of this equation with Riemann-Liouville fractional derivative we could express the kinetic equation as a differential equation with fractional order, i.e., ( ) ( ) In order not to lose the initial conditions of the quantum systems, we apply the integral operator ( This differential equation with negative fractional order could be further processed through manipulating the fractional differential operator.This manipulation yields ( ) ( ) ( ) ( ) Here we have expressed the kinetic equation of this quantum system with memory effect as a FGLE through using the fractional time derivative for the memory kernel.
When we proceed to solve the FGLE through Laplace transform, we need the formula of Laplace transform in fractional order with Laplace variable s .Here the operator of fractional calculus a t D ν is defined as with

[ ]
Re ν being the real part of the order ν and the fractional derivative d dt ν ν being defined through the Riemann-Liouvile form in Equation (1.4).And the Laplace transform of exponential order 1 t ν − − could be obtained from elementary calculus as With these two formula, the Laplace transform of the FGLE has the algebraic form as ( ) ( ) ( ) This algebraic form could be expressed as a sum of partial fractions as with expansion coefficients n a and n X being the roots of the indicial Equation (1.11).
When the inverse Laplace transform is performed on these fractional expansions of ( )

A s 
, we need the formula ( ) ( )  ( ) which is a linear combination of the two-parameter Mittag-Leffler functions.These two-parameter Mittag-Leffler functions are defined as with the derivative formula They are related to the one-parameter Mittag-Leffler function

A Particular Memory Kernel
In this section, we consider a quantum system of a qubit in an anisotropic photonic crystal (PhC) with a particular memory kernel as shown in Figure 1.The Hamiltonian of this quantum system is ( ) .
Here the atomic operators  As we assume only one photon is created or annihilated for one atomic transition (single-photon sector), the wave function of the system has the form with the initial condition ( ) 0 0 0 e cos 2 (1.20) and ( ) with detuning frequency . The two Equations of (1.19) and (1.20) can be combined as with the memory kernel . This evolution equation relates the excited amplitude ( ) ω from the band edge frequency c ω [7].With this special form of memory kernel, the kinetic equation of the excited amplitude ( ) if the transformation ( ) ( ) is performed.Here we express this memory effect as a fractional time derivative which describes the effective interaction of the qubit with the environmental PhC reservoir.That is, the right-hand-side term of the kinetic Equation (1.24) is written as a Riemann-Louville fractional time derivative with order 1 2 This expression leads to the fractional form of the kinetic equation as We define this fractional differential equation as a fractional generalized Langevin equation (FGLE) of this quantum system.We solve this FGLE by applying Laplace transform and obtain ( ) Here ( )

C s 
is the Laplace transform of ( ) C t with the initial condition ( ) ( ) In order to find the solution of this excited amplitude, standard procedures of expressing this algebraic equation as a sum of partial fractions and performing inverse Laplace transform on these partial fractions are taken.In the first step, we need to find the roots of the indicial equation 2 For the degenerate case, we have β = ∆ which leads to the indicial equation as ( ) The partial fractions of ( ) As the inverse Laplace transform is applied to these partial-fractional forms of ( ) C s  with fractional powers of the variables s , we use the formulas of ( ) ( ) with the fractional exponential function , whose definition and properties are listed in Equations (1.15) and (1.16).This procedure leads to the solution of the fractional kinetic equation being expressed as the linear combination of the fractional exponential functions such that for the different-root case β ≠ ∆ ; and for the degenerate case β = ∆ .Based on this analytic solution, we can obtain the reduced density matrix of the qubit from the wave function in Equation (1.18) through tracing over the reservoir degrees of freedom as The elements in this matrix are associated with the information of the qubit energy and coherence.In the following, we will study the dynamics of the excited-state energy, polarization and von Neumann entropy of the qubit.

Excited-State Probability
We show the time evolution of the qubit excited-state energy through the probability ( ) ( ) ( ) in Figure 2. The energy dynamics exhibits oscillatory behavior and does not decay with time as the qubit frequency lies inside the photonic bang gap (PBG) region ( 0 β ∆ < ).The memory effect of this quantum system is observed directly through these trapping states of the qubit energy.

Polarization Dynamics
As a qubit interacts with the environment, it will randomize the polarization of the qubit.The quantum phase information of a qubit carried by the qubit polarization will thus escape from the qubit into the environment through this randomization of polarization and lead to the quantum decoherence.Here we show the polarization dynamics in Figure 3 through the expression of qubit polarization Re ).The polarization dynamics of the qubit with frequency lying inside the PBG region ( 0 β ∆ < ) exhibits non-decaying oscillation.The qubit loses partial of its polarization in the very beginning period of time and then preserves the remaining polarization through the steady oscillation.This trapping state of the qubit polarization reveals the memory effect of the system which leads to the preservation of the qubit phase information.

Dynamics of von Neumann Entropy
Entropy, a measurement of information amount stored in a qubit, will be changed as the qubit is correlated to the environment.The correlation between the environment and the state will transform the initially pure state of the qubit into a finally mixed state where the amount of information of the qubit is changed.For a qubit state with density matrix ( )  4 for the initially excited qubit ( 0 0 θ = ).The entropy has its minimal value zero at 0 t = and reaches its maximal value log 2 0.693 = at the very beginning of time.After a period of time on the order of the decay timescale, the entropy becomes steady with nonzero value for the qubit frequency in PBG region ( 0 β ∆ < ).This result shows that the initially pure system becomes maximally mixed in the very beginning period of time.As the qubit equilibrating with the PhC reservoir, the system becomes steady with less mixed state.The trapping state of the von Neumann entropy reveals that the amount of information stored in the qubit   is preserved through the steady mixed state.The memory effect of the system is observed directly through the preservation of the qubit information.

Conclusion
We have used the fractional time derivative to express the kinetic equation of the quantum system with memory effect as a FGLE.For a particular memory kernel, we obtain the solution of the FGLE in terms of the twoparameter Mittag-Leffler function.In the study of the qubit dynamics in the particular memory kernel of an anisotropic PhC, we observe the memory effect directly through the trapping states of the qubit energy, polarization and von Neumann entropy.

.
13) with positive integer 1 q ν = and the two-parameter Mittag-Leffler function ( ) The analytical solution of the FGLE is obtained as ( )

10 ω
stands for the atomic transition frequency from excited state 1 to ground state 0 and k ω for the photon mode frequency of the reservoir with wave- vector k  .The coupling strength between the atom and the photon (electromagnetic field) is characterized by V , dielectric constant 0  and polarization unit vector ˆk e of the photon mode with frequency k ω .

Figure 1 .
Figure 1.(Color online) (a) A qubit with excited state 1 and ground state 0 .The transition frequency 10 ω is nearly resonant with the frequency range of the PhC reservoir; (b) Directional dependent dispersion relations near band edge expressed by the effective-mass approximation with the edge frequency c ω by solid and dashed curves; (c) Photon DOS the single-photon sector of the Hilbert space as A t of the qubit to the reservoir memory through the memory kernel ( ) K t τ − .For the anisotropic PhC reservoir, this memory kernel has a particular form of

=
, where the variable 1 2 s has been converted into Y .Two kinds of roots exist in this indicial equation: one with different roots .For the case of different roots, ( )

Figure 4 .
Figure 4. (Color online) Dynamics of the von Neumann entropy for different qubit frequencies.