Quantum Non-Markovian Environment-to-System Backflows of Information: Nonoperational vs. Operational Approaches.

Quantum memory effects can be qualitatively understood as a consequence of an environment-to-system backflow of information. Here, we analyze and compare how this concept is interpreted and implemented in different approaches to quantum non-Markovianity. We study a nonoperational approach, defined by the distinguishability between two system states characterized by different initial conditions, and an operational approach, which is defined by the correlation between different outcomes associated to successive measurement processes performed over the system of interest. The differences, limitations, and vantages of each approach are characterized in detail by considering diverse system-environment models and dynamics. As a specific example, we study a non-Markovian depolarizing map induced by the interaction of the system of interest with an environment characterized by incoherent and coherent self-dynamics.

Emergence of stationary multimodality under two-timescaled dichotomic noise.

We study a linear Langevin dynamics driven by an additive non-Markovian symmetrical dichotomic noise. It is shown that when the statistics of the time intervals between noise transitions is characterized by two well differentiated timescales, the stationary distribution may develop multimodality (bi- and trimodality). The underlying effects that lead to a probability concentration in different points include intermittence and also a dynamical locking of realizations. Our results are supported by numerical simulations as well as by an exact treatment obtained from a Markovian embedding of the full dynamics, which leads to a third-order differential equation for the stationary distribution.

First-passage time for superstatistical Fokker-Planck models.

The first-passage-time (FPT) problem is studied for superstatistical models assuming that the mesoscopic system dynamics is described by a Fokker-Planck equation. We show that all moments of the random intensive parameter associated to the superstatistical approach can be put in one-to-one correspondence with the moments of the FPT. For systems subjected to an additional uncorrelated external force, the same statistical information is obtained from the dependence of the FPT moments on the external force. These results provide an alternative technique for checking the validity of superstatistical models. As an example, we characterize the mean FPT for a forced Brownian particle.

Quantum Non-Markovian Processes Break Conditional Past-Future Independence.

For classical Markovian stochastic systems, past and future events become statistically independent when conditioned to a given state at the present time. Memory non-Markovian effects break this condition, inducing a nonvanishing conditional past-future correlation. Here, this classical memory indicator is extended to a quantum regime, which provides an operational definition of quantum non-Markovianity based on a minimal set of three time-ordered quantum system measurements and postselection. The detection of memory effects through the measurement scheme is univocally related to departures from Born-Markov and white noise approximations in quantum and classical environments, respectively.

Memory-induced diffusive-superdiffusive transition: Ensemble and time-averaged observables.

The ensemble properties and time-averaged observables of a memory-induced diffusive-superdiffusive transition are studied. The model consists in a random walker whose transitions in a given direction depend on a weighted linear combination of the number of both right and left previous transitions. The diffusion process is nonstationary, and its probability develops the phenomenon of aging. Depending on the characteristic memory parameters, the ensemble behavior may be normal, superdiffusive, or ballistic. In contrast, the time-averaged mean squared displacement is equal to that of a normal undriven random walk, which renders the process nonergodic. In addition, and similarly to Lévy walks [Godec and Metzler, Phys. Rev. Lett. 110, 020603 (2013)PRLTAO0031-900710.1103/PhysRevLett.110.020603], for trajectories of finite duration the time-averaged displacement apparently become random with properties that depend on the measurement time and also on the memory properties. These features are related to the nonstationary power-law decay of the transition probabilities to their stationary values. Time-averaged response to a bias is also calculated. In contrast with Lévy walks [Froemberg and Barkai, Phys. Rev. E 87, 030104(R) (2013)PLEEE81539-375510.1103/PhysRevE.87.030104], the response always vanishes asymptotically.

Inhomogeneous diffusion and ergodicity breaking induced by global memory effects.

We introduce a class of discrete random-walk model driven by global memory effects. At any time, the right-left transitions depend on the whole previous history of the walker, being defined by an urnlike memory mechanism. The characteristic function is calculated in an exact way, which allows us to demonstrate that the ensemble of realizations is ballistic. Asymptotically, each realization is equivalent to that of a biased Markovian diffusion process with transition rates that strongly differs from one trajectory to another. Using this "inhomogeneous diffusion" feature, the ergodic properties of the dynamics are analytically studied through the time-averaged moments. Even in the long-time regime, they remain random objects. While their average over realizations recovers the corresponding ensemble averages, departure between time and ensemble averages is explicitly shown through their probability densities. For the density of the second time-averaged moment, an ergodic limit and the limit of infinite lag times do not commutate. All these effects are induced by the memory effects. A generalized Einstein fluctuation-dissipation relation is also obtained for the time-averaged moments.

Weak ergodicity breaking induced by global memory effects.

We study the phenomenon of weak ergodicity breaking for a class of globally correlated random walk dynamics defined over a finite set of states. The persistence in a given state or the transition to another one depends on the whole previous temporal history of the system. A set of waiting time distributions, associated to each state, sets the random times between consecutive steps. Their mean value is finite for all states. The probability density of time-averaged observables is obtained for different memory mechanisms. This statistical object explicitly shows departures between time and ensemble averages. While the residence time in each state may have a divergent mean value, we demonstrate that this condition is in general not necessary for breaking ergodicity. Hence, we conclude that global memory effects are an alternative mechanism able to induce ergodicity breaking without involving power-law statistics. Analytical and numerical calculations support these results.

Central limit theorem for a class of globally correlated random variables.

The standard central limit theorem with a Gaussian attractor for the sum of independent random variables may lose its validity in the presence of strong correlations between the added random contributions. Here, we study this problem for similar interchangeable globally correlated random variables. Under these conditions, a hierarchical set of equations is derived for the conditional transition probabilities. This result allows us to define different classes of memory mechanisms that depend on a symmetric way on all involved variables. Depending on the correlation mechanisms and statistics of the single variables, the corresponding sums are characterized by distinct probability densities. For a class of urn models it is also possible to characterize their domain of attraction, which, as in the standard case, is parametrized by the probability density of each random variable. Symmetric and asymmetric q-Gaussian attractors (q<1) arise in a particular two-state case of these urn models.

Extended q-Gaussian and q-exponential distributions from gamma random variables.

The family of q-Gaussian and q-exponential probability densities fit the statistical behavior of diverse complex self-similar nonequilibrium systems. These distributions, independently of the underlying dynamics, can rigorously be obtained by maximizing Tsallis "nonextensive" entropy under appropriate constraints, as well as from superstatistical models. In this paper we provide an alternative and complementary scheme for deriving these objects. We show that q-Gaussian and q-exponential random variables can always be expressed as a function of two statistically independent gamma random variables with the same scale parameter. Their shape index determines the complexity q parameter. This result also allows us to define an extended family of asymmetric q-Gaussian and modified q-exponential densities, which reduce to the standard ones when the shape parameters are the same. Furthermore, we demonstrate that a simple change of variables always allows relating any of these distributions with a beta stochastic variable. The extended distributions are applied in the statistical description of different complex dynamics such as log-return signals in financial markets and motion of point defects in a fluid flow.

Post-Markovian quantum master equations from classical environment fluctuations.

In this paper we demonstrate that two commonly used phenomenological post-Markovian quantum master equations can be derived without using any perturbative approximation. A system coupled to an environment characterized by self-classical configurational fluctuations, the latter obeying a Markovian dynamics, defines the underlying physical model. Both Shabani-Lidar equation [A. Shabani and D. A. Lidar, Phys. Rev. A 71, 020101(R) (2005)] and its associated approximated integrodifferential kernel master equation are obtained by tracing out two different bipartite Markovian Lindblad dynamics where the environment fluctuations are taken into account by an ancilla system. Furthermore, conditions under which the non-Markovian system dynamics can be unraveled in terms of an ensemble of measurement trajectories are found. In addition, a non-Markovian quantum jump approach is formulated. Contrary to recent analysis [L. Mazzola, E. M. Laine, H. P. Breuer, S. Maniscalco, and J. Piilo, Phys. Rev. A 81, 062120 (2010)], we also demonstrate that these master equations, even with exponential memory functions, may lead to non-Markovian effects such as an environment-to-system backflow of information if the Hamiltonian system does not commutate with the dissipative dynamics.

Generalized fluctuation relation for power-law distributions.

Strong violations of existing fluctuation theorems may arise in nonequilibrium steady states characterized by distributions with power-law tails. The ratio of the probabilities of positive and negative fluctuations of equal magnitude behaves in an anomalous nonmonotonic way [H. Touchette and E. G. D. Cohen, Phys. Rev. E 76, 020101(R) (2007)]. Here, we propose an alternative definition of fluctuation relation (FR) symmetry that, in the power-law regime, is characterized by a monotonic linear behavior. The proposal is consistent with a large deviationlike principle. As an example, we study the fluctuations of the work done on a dragged particle immersed in a complex environment able to induce power-law tails. When the environment is characterized by spatiotemporal temperature fluctuations, distributions arising in nonextensive statistical mechanics define the work statistics. In that situation, we find that the FR symmetry is solely defined by the average bath temperature. The case of a dragged particle subjected to a Lévy noise is also analyzed in detail.

Large deviations of ergodic counting processes: a statistical mechanics approach.

The large-deviation method allows to characterize an ergodic counting process in terms of a thermodynamic frame where a free energy function determines the asymptotic nonstationary statistical properties of its fluctuations. Here we study this formalism through a statistical mechanics approach, that is, with an auxiliary counting process that maximizes an entropy function associated with the thermodynamic potential. We show that the realizations of this auxiliary process can be obtained after applying a conditional measurement scheme to the original ones, providing is this way an alternative measurement interpretation of the thermodynamic approach. General results are obtained for renewal counting processes, that is, those where the time intervals between consecutive events are independent and defined by a unique waiting time distribution. The underlying statistical mechanics is controlled by the same waiting time distribution, rescaled by an exponential decay measured by the free energy function. A scale invariance, shift closure, and intermittence phenomena are obtained and interpreted in this context. Similar conclusions apply for nonrenewal processes when the memory between successive events is induced by a stochastic waiting time distribution.

Fluctuation relations with intermittent non-Gaussian variables.

Nonequilibrium stationary fluctuations may exhibit a special symmetry called fluctuation relations (FRs). Here, we show that this property is always satisfied by the subtraction of two random and independent variables related by a thermodynamiclike change of measure. Taking one of them as a modulated Poisson process, it is demonstrated that intermittence and FRs are compatible properties that may coexist naturally. Strong non-Gaussian features characterize the probability distribution and its generating function. Their associated large deviation functions develop a "kink" at the origin and a plateau regime respectively. Application of this model in different stationary nonequilibrium situations is discussed.

Thermodynamics of quantum jump trajectories in systems driven by classical fluctuations.

The large-deviation method can be used to study the measurement trajectories of open quantum systems. For optical arrangements this formalism allows to describe the long time properties of the (nonequilibrium) photon counting statistics in the context of a (equilibrium) thermodynamic approach defined in terms of dynamical phases and transitions between them in the trajectory space [J. P. Garrahan and I. Lesanovsky, Phys. Rev. Lett. 104, 160601 (2010)]. In this paper, we study the thermodynamic approach for fluorescent systems coupled to complex reservoirs that induce stochastic fluctuations in their dynamical parameters. In a fast modulation limit the thermodynamics corresponds to that of a Markovian two-level system. In a slow modulation limit, the thermodynamic properties are equivalent to those of a finite system that in an infinite-size limit is characterized by a first-order transition. The dynamical phases correspond to different intensity regimes, while the size of the system is measured by the transition rate of the bath fluctuations. As a function of a dimensionless intensive variable, the first and second derivatives of the thermodynamic potential develop an abrupt change and a narrow peak, respectively. Their scaling properties are consistent with a double-Gaussian probability distribution of the associated extensive variable.

From power law intermittence to macroscopic coherent regime.

We address the problem of establishing which is the proper form of quantum master equation generating a survival probability identical to that corresponding to the nonergodic sequence of "light on" and "light off" fluorescence fluctuations in blinking quantum dots. We adopt a theoretical perspective based on the assumption that the abrupt transitions from the light on to light off state are the results of many collisions between system and environment, properly described by the Lindblad equation, and that between two consecutive collisions the system dynamics are frozen. This generates a quantum master equation belonging to the recently proposed class of generalized Lindblad equations, with a time convoluted structure, involving in the specific case of this paper both the unitary and the nonunitary contribution of the Lindblad equation. This is the property that under the low-frequency condition makes the new class of generalized Lindblad equation generates the required survival probability. We make the conjecture that this equation corresponds to the cooperative dynamics of many units that, in isolation, are described by the ordinary Lindblad equation. When the time scale of the unitary term of the Lindblad equation is shorter than the dephasing time, the cooperation generates a surprisingly extended macroscopic coherence.

Generalized non-Markovian optical Bloch equations.

By considering single chromophore systems whose radiative decay can be written in terms of a nonlocal Lindblad-type evolution, the authors extend the formalism of generalized optical Bloch equations [Y. Zheng and F. L. H. Brown, Phys. Rev. Lett. 90, 238305 (2003)] to non-Markovian dynamics. They demonstrate that photon statistical properties such as bunching and antibunching, as well as sub- and super-Poissonian photon statistics can be fitted in the context of non-Markovian dynamics. The nonlocal effects may arise due to the interaction with a complex structured environment. In this case, the photon statistics can be related with the parameters that define the microscopic system-environment interaction. Alternatively, the authors demonstrate that effective dynamics such as triplet blinking, where the system is coupled via incoherent transitions to an extra dark state, can also be worked out in terms of generalized non-Markovian optical Bloch equations. The corresponding memory contributions are mapped with those that arise from the microscopic approach.

Random Lindblad equations from complex environments.

In this paper we demonstrate that Lindblad equations characterized by a random rate variable arise after tracing out a complex structured reservoir. Our results follows from a generalization of the Born-Markov approximation, which relies on the possibility of splitting the complex environment into a direct sum of subreservoirs, each one being able to induce by itself a Markovian system evolution. Strong non-Markovian effects, which microscopically originate from the entanglement with the different subreservoirs, characterize the average system decay dynamics. As an example, we study the anomalous irreversible behavior of a quantum tunneling system described in an effective two-level approximation. Stretched exponential and power law decay behaviors arise from the interplay between the dissipative and unitary hopping dynamics.

Functional characterization of linear delay Langevin equations.

We present an exact functional characterization of linear delay Langevin equations driven by any noise structure defined through its characteristic functional. This method relies on the possibility of finding an explicitly analytical expression for each realization of the delayed stochastic process in terms of those of the driving noise. General properties of the transient dissipative dynamics are analyzed. The corresponding interplay with a color Gaussian noise is presented. As a full application of our functional method we study a model for population growth with non-Gaussian fluctuations: the Gompertz model driven by multiplicative white shot noise.