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.
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.
We define non-Markovian quantum dynamics as evolution in which the current state depends on all past states, and completely characterize its structure under the assumptions of complete positivity and non-signalling. The resulting continuous-time dynamics is an integro-differential equation that augments the Gorini-Kossakowski-Sudarshan-Lindblad equation with a memory integral, and is capable of describing the quantum state of systems exposed to noise with any integrable power spectral density with no further approximations. We then establish a formalism to evaluate multi-time correlations of measurement outcomes in this general setting, obviating the need for a regression theorem. As an application, we derive the emission spectrum of a driven two-level system coupled to a non-Markovian bath: the familiar Mollow triplet acquires a frequency-dependent linewidth that encodes the memory of the bath. Our work provides a rigorous yet transparent description of the quantum state of non-Markovian systems, opening the door for state estimation and state-based quantum control beyond the Markovian regime.
In this paper, we proposed the exactly solvable model of non-Markovian dynamics of open quantum systems. This model describes open quantum systems with memory and periodic sequence of kicks by environment. To describe these systems, the Lindblad equation for quantum observable is generalized by taking into account power-law fading memory. Dynamics of open quantum systems with power-law memory are considered. The proposed generalized Lindblad equations describe non-Markovian quantum dynamics. The quantum dynamics with power-law memory are described by using integrations and differentiation of non-integer orders, as well as fractional calculus. An example of a quantum oscillator with linear friction and power-law memory is considered. In this paper, discrete-time quantum maps with memory, which are derived from generalized Lindblad equations without any approximations, are suggested. These maps exactly correspond to the generalized Lindblad equations, which are fractional differential equations with the Caputo derivatives of non-integer orders and periodic sequence of kicks that are represented by the Dirac delta-functions. The solution of these equations for coordinates and momenta are derived. The solutions of the generalized Lindblad equations for coordinate and momentum operators are obtained for open quantum systems with memory and kicks. Using these solutions, linear and nonlinear quantum discrete-time maps are derived.
Existing quantum algorithms for simulating dynamical systems -- from Hamiltonian simulation to linear and nonlinear differential equations solvers -- simulate Markovian dynamics, in which the system's future evolution depends solely on its current state. We turn our attention to developing quantum algorithms for non-Markovian dynamical systems where the system's future evolution depends on its past history and thus has memory. Specifically, we develop efficient algorithms for linear Volterra integro-differential equations (VIDEs) with a convolution memory kernel that output a quantum state encoding the state description over a time interval or at a particular time. Given efficient circuits for the problem inputs, our algorithms achieve an exponential speedup in system size over existing classical algorithms.
We develop an algorithm for general kernels assuming that M<1, where M characterizes the strength of the memory term relative to the dissipation of the Markovian part of the dynamics. We complement this with lower bounds for general-kernel VIDEs when M≥1, showing that the problem becomes intractable for a family of systems. However, by specializing to structured kernels which admit concise decompositions over exponentials, we develop efficient quantum algorithms even when M≥1 by converting the VIDE into a larger set of ODEs, a procedure which we call Markovianization. As an application of the overall framework, we discuss the Mori-Zwanzig formalism used in open quantum systems and fluid dynamics. Overall, our results expand the range of dynamical systems that quantum computers can simulate efficiently.
Abstract We present a general numerical simulation method to solve non-Markovian dynamics of an open quantum system influenced by quantum Brownian motion. Based on the determined memory kernel function, this method enables the resolution of non-Markovian dynamics for a wide range of system Hamiltonians and spectral densities. The system dynamics are described by exact integro-differential operator equations without any common approximations and they are simulated in this work by definite-number equations with stochastic initial conditions. This approach ensures the applicability of mature numerical methods and maintains computational complexity that remains largely invariant, even when dealing with more complex models. The high accuracy of our simulation is evident from a comparison with the results obtained from corresponding exact master equations, underscoring the reliability and precision of our method.
Open quantum systems exhibit dynamics ranging from unitary evolution to irreversible dissipation. While the Gorini--Kossakowski--Sudarshan--Lindblad (GKSL) equation uniquely characterizes Markovian CPTP evolution, many physical platforms display non-Markovian features such as algebraic relaxation and coherence backflow. Fractional calculus provides a natural way to model such long-memory behavior through power-law temporal kernels introduced by fractional time derivatives. Here we develop a unified framework that embeds fractional master equations within the broader hierarchy of open-system formalisms. The fractional equation forms a structured subclass of memory-kernel models, reduces to the Lindblad form at unit order, and, through Bochner--Phillips subordination, admits a CPTP representation as an average over Lindblad semigroups. Its resolvent structure further connects fractional dynamics to established non-Markovian approaches, including Nakajima--Zwanzig kernels and hierarchical equations of motion, providing a compact surrogate for long-memory effects. This formulation positions fractional calculus as a rigorous and practical language for quantum dynamics with intrinsic memory, supporting both analytical insight and efficient quantum simulation.
Applications of Time-Fractional Schrodinger Equations (TFSEs) to quantum processes are instructive for understanding and describing the time behavior of real physical systems. By applying three popular TFSEs, namely Naber's TFSE I, Naber's TFSE II, and XGF's TFSE, to a basic open system model of a two-level system (qubit) coupled resonantly to a dissipative environment, we solve exactly for Time-Fractional Single Qubit Open Systems (TFSQOSs). However, the three TFSEs perform badly for the following reasons. On the other hand, in the respective frameworks of the three TFSEs, the total probability for obtaining the system in a single-qubit state is not equal to one with time at fractional order, implying that time-fractional quantum mechanics violates quantum mechanical probability conservation. On the other hand, the latter two TFSEs are not capable of describing the non-Markovian dynamics of the system at all fractional order, only at some fractional order. To address this, we introduce a well-performed TFSE by constructing a new analytic continuation of time combined with the conformable fractional derivative, in which for all fractional order, not only does the total probability for the system equal one at all times but also the non-Markovian features can be observed throughout the time evolution of the system. Furthermore, we study the performances of the four TFSEs applying to an open system model of two isolated qubits each locally interacting with its dissipative environment. By deriving the exact solutions for time-fractional two qubits open systems, we show that our TFSE still possesses the above two advantages compared with the other three TFSEs.