Non-Markovianity of qubit evolution under the action of spin environment

TL;DR

This paper investigates the non-Markovian behavior of a qubit interacting with a spin environment, analyzing how different coupling types and strengths influence the information flow and the transition between non-Markovian and Markovian dynamics.

Contribution

It provides a detailed analysis of non-Markovianity in qubit evolution under various coupling conditions with a spin bath environment, including explicit Kraus operators.

Findings

Time-independent and time-polynomial couplings induce non-Markovianity.

Certain parameter values of time-exponential coupling also lead to non-Markovian behavior.

Transition from non-Markovian to Markovian dynamics occurs as coupling strength varies.

Abstract

The question, whether an open system dynamics is Markovian or non-Markovian can be answered by studying the direction of the information flow in the dynamics. In Markovian dynamics, information must always flow from the system to the environment. If the environment is interacting with only one of the subsystems of a bipartite system, the dynamics of the entanglement in the bipartite system can be used to identify the direction of information flow. Here we study the dynamics of a two-level system interacting with an environment, which is also a heat bath, and consists of a large number of two-level quantum systems. Our model can be seen as a close approximation to the `spin bath' model at low temperatures. We analyze the Markovian nature of the dynamics, as we change the coupling between the system and the environment. We find the Kraus operators of the dynamics for certain classes of couplings. We show that any form of time-independent or time-polynomial coupling gives rise to non-Markovianity. Also, we witness non-Markovianity for certain parameter values of time-exponential coupling. Moreover, we study the transition from non-Markovian to Markovian dynamics as we change the value of coupling strength.

Equations

$${\Lambda }_t=V_{t,s}\circ {\Lambda }_s,$$

$$E({\rho }_{as})=\frac{||({\rho }_{as}{)}^{T_s}|{|}_1-1}{2}$$

$$\ket{\Phi^{+}}=\frac{1}{\sqrt{2}}\big{(}\ket{00}+\ket{11}\big{)}.$$

$${\tilde{H}}_{se}(t)$$

$$H_s$$

$$H_e$$

$$H_{\text{se}}(t)$$

$$\displaystyle+\ket{0}_{s}\bra{1}\otimes\sum_{n=1}^{N}\tilde{g}_{n}(t)e^{-i\delta\omega_{n}\,t}\ket{0..1_{n}..0}_{e}\bra{0..0_{n}..0}\Big{]}$$

$$\displaystyle=\hbar\alpha\Big{[}\ket{1}_{s}\bra{0}\otimes\sum_{n=1}^{N}g^{*}_{n}(t)\ket{0..0_{n}..0}_{e}\bra{0..1_{n}..0}$$

$$\displaystyle+\ket{0}_{s}\bra{1}\otimes\sum_{n=1}^{N}g_{n}(t)\ket{0..1_{n}..0}_{e}\bra{0..0_{n}..0}\Big{]}$$

$$\rho_{e}(0)=\Big{[}p\ket{0}\bra{0}+(1-p)\ket{1}\bra{1}\Big{]}^{\otimes N},$$

$$H_{\text{se}}(t)$$

$$\displaystyle+g_{n}(t)\sigma_{-}^{(s)}\otimes\big{[}\ket{0}\bra{0}\otimes..\sigma_{+}^{(n)}..\otimes\ket{0}\bra{0}\big{]}_{e}\Big{\}}$$

$$H_{spin-bath}=\hbar\alpha\sum_{n=1}^{N}\big{(}\sigma_{x}^{(s)}\sigma_{x}^{(n)}+\sigma_{y}^{(s)}\sigma_{y}^{(n)}+\sigma_{z}^{(s)}\sigma_{z}^{(n)}\big{)}$$

$${\mathcal{E}}_{\pm }(t)=\pm \hslash \alpha \sqrt{\sum _{n=1}^N|g_n(t){|}^2}=\pm \mathcal{E}(t),$$

$$\displaystyle\ket{\chi_{\pm}(t)}_{se}=\frac{1}{\sqrt{2}}\big{[}\ket{1}_{s}\otimes\ket{0}^{\otimes N}_{e}\pm\ket{\xi(t)}_{se}\big{]},$$

$${|{\beta }_0⟩}_e=\frac{\hslash \alpha }{\mathcal{E}(t)}\sum _{n=1}^N{g}_n(t){|\mathrm{0..}{1}_n..0⟩}_e.$$

$$U(t,0)=\mathcal{T}\exp \left[-\frac{i}{\hslash }{\int }_0^t{H}_{se}({\tau }^{\prime })d{\tau }^{\prime }\right],$$

$${\rho }_{\text{ase}}(0)=|{\Phi }^{+}⟩⟨{\Phi }^{+}|\otimes {\rho }_e(0),$$

$$\rho_{ase}(t)=\big{[}\mathbb{I}_{a}\otimes U(t,0)\big{]}\rho_{\text{ase}}(0)\big{[}\mathbb{I}_{a}\otimes U^{\dagger}(t,0)\big{]}.$$

$${\rho }_{as}(t)={\text{Tr}}_e{\rho }_{\text{ase}}(t).$$

$$U(t,0)=$$

$$+\left(e^{i{\omega }_{0}t}-1\right)|{\chi }_{-}⟩⟨{\chi }_{-}|+\text{openone}.$$

$${\rho }_s(0)\to {\rho }_s(t)=\sum _{m,n=1}^N{K}_{mn}(t){\rho }_s(0)K_{mn}^{†}(t),$$

$$K_{mn}(t)=\sqrt{p^{N-s_{n}}(1-p)^{s_{n}}}\left[\begin{array}[]{c c c}\frac{\cos\omega_{0}t-1}{\omega_{0}^{2}}~{}g_{N-\log m}(t)~{}g^{*}_{N-\log n}(t)+\delta_{mn}&&-\frac{i}{\omega_{0}}\sin\omega_{0}t~{}g_{N-\log m}(t)~{}\delta_{0n}\\ \\ -\frac{i}{\omega_{0}}\sin\omega_{0}t~{}g^{*}_{N-\log n}(t)~{}\delta_{0m}&&(\cos\omega_{0}t-1)~{}\delta_{0m}\delta_{0n}+\delta_{mn}\\ \end{array}\right]$$

$${\rho }_{\text{as}}(t)=$$

$$\displaystyle+(1-p)\kappa_{0}\big{(}\ket{00}\bra{00}+\ket{01}\bra{01}\big{)}$$

$$\displaystyle+\delta_{0}(\ket{00}\bra{11}+\ket{11}\bra{00})\Big{]},$$

$$\lambda (t)=\frac{1}{2}\left({\kappa }_0-\sqrt{(1-2p{\kappa }_0{)}^2+4{\delta }_0^2}\right).$$

$$\Omega (t)$$

$$=\sqrt{N}{\int }_0^t|g(\tau )|d\tau .$$

$$U_{se}(t,0)$$

Full text

I Introduction

We rarely come across systems that are completely isolated from the surrounding world. Had it been the case, dealing with quantum mechanical systems would have been lot more easier. So, although arduous to deal with, real quantum systems are mostly open quantum systems $-$ a system interacting with an environment. In these situations, information exchange between system and environment becomes an essential feature. Information that has been previously transferred to the environment may come back and affect the system, and this may appear as a memory-effect on the system. When this information backflow from the environment is negligible we have a situation analogous to the discrete Markov process, where the instantaneous state of the system depends solely on the immediately previous step, the system dynamics is called memory-less or Markovian Breuer et al. (2016); Rivas et al. (2014). On the other hand, when this information backflow affects the system significantly i.e. when some long past history of the system influences its present state, the system dynamics becomes retentive, and is called non-Markovian.

In recent years, non-Markovianity has been used as a resource in a number of information theoretic protocols, namely, channel discrimination Bylicka et al. (2014), preserving coherence and correlation Man et al. (2015a); Franco (2016) and retrieving quantum correlations in both quantum and classical environments Bellomo et al. (2012); González-Gutiérrez et al. (2016a); Xu et al. (2013); Franco et al. (2013); Bellomo et al. (2011). Non-Markovian effects also play important roles in areas ranging from fundamental physics of strong fields Schmidt et al. (1999); Bloch et al. (2000) to energy transfer process of photosynthetic complexes Rebentrost and Aspuru-Guzik (2011).

Owing to its diverse applications, various aspects of non-Markovianity are now being studied. Lately, researchers have been focusing on transition from non-Markovian to Markovian dynamics Liu et al. (2011); Bernardes et al. (2014); Brito and Werlang (2015); Garrido et al. (2016); Man et al. (2015b); Franco (2015). Some of them have dealt with bosonic bath of infinite or finite degrees of freedom, while some have considered a qudit system as the environment. But in all of these studies, system-environment interaction has been considered to be homogeneous in space, and the issue of non-Markovian to Markovian transition in terms of system-environment coupling strength has not been addressed. Note that, non-Markovian to Markovian transition is, in general, not a trivial issue, as in most cases finite dimensional environments give rise to non-Markovianity.

In our study, we attempt to analyze the problem of whether a transition from non-Markovianity to Markovianity can be engineered for the spin bath model Prokof’ev and Stamp (2000). We particularly choose the spin bath model since it has wide ranging applications in simulating real physical scenarios Hutton and Bose (2004); Breuer et al. (2004); Prokof’ev and Stamp (2000). In our attempt, we face a serious difficulty in diagonalizing the spin bath Hamiltonian, either analytically or numerically, for larger number of spins in the environment. Although, analytic solutions do exist for constant coupling Bhattacharya et al. (2017) and some special forms of time dependent coupling Jing and Wu (2018), general solution for arbitrary forms of system-environment coupling of the spin bath Hamiltonian are hard to find. We therefore, try to circumvent the problem by choosing a simple model, which we argue, is a close approximation to the spin bath model for low temperatures. We choose an exchange type of interaction between a system qubit and individual environment qubits, where for each environment qubit the coupling can be chosen to be of different time dependent forms. But unlike the spin bath case, in our model, when the exchange interaction takes place between the system and a particular environment qubit, the rest of the environment qubits remain in a ground state; which also closely resembles the state of environment for low temperatures. As we will see in the paper, this approximation helps us to calculate and analyze non-Markovian to Markovian transition for different types of system-environment coupling.

We present four scenarios here, for different forms of system-environment coupling: (i) the coupling is time-independent and homogeneous over environment qubits, (ii) the coupling is time-independent but inhomogeneous over environment qubits, (iii) the coupling is homogeneous over the environment but is time-dependent, and (iv) the coupling is both time-dependent and inhomogeneous. We find that cases (i) and (ii) always give rise to non-Markovian system dynamics. For cases (iii) and (iv), we find that some functional forms of coupling for certain ranges of coupling strengths gives rise to non-Markovianity. For example in case (iii), polynomial forms of coupling always give rise to non-Markovian system dynamics, while exponential coupling give rise to non-Markovian system dynamics only for certain ranges of parameter values. In case (iv) we find that a cross-over from non-Markovianity to Markovianity can be achieved by varying the strength of coupling. We also calculate, the extremal values of coupling parameter beyond which non-Markovianity can no longer be detected. Thus we see, these extremal values act as critical values for transition from non-Markovian to Markovian regime. It is worth mentioning here that, for the purpose of detecting non-Markovianity we use Rivas-Huelga-Plenio (RHP) measure of non-Markovianity as proposed in Rivas et al. (2010). Although there are different approaches of defining Markovianity and each approach represent different aspects of Markovianity, for the purpose of the present paper we choose, detection by the RHP measure as the definition of Markovianity.

Similar works on this line were done in Apollaro et al. (2011); Lorenzo et al. (2013); Wang et al. (2013). But in the first approach Apollaro et al. (2011), the system qubit directly interacts with a single environment qubit and the rest of the environment qubits, only have an indirect effect on the system via the environment qubit directly attached. Also, the coupling parameters involved do not have any time dependence. In the second approach Lorenzo et al. (2013), the transition from Markovianity to non-Markovianity was shown with a two tier environment; the first one being a multiple-spin system, while the second one was a bosonic bath. Also in Wang et al. (2013), the coupling between the system and individual environment qubits were constant in space and time. We take into account all these factors and present a detailed study of a spin environment and cover all the relevant cases.

In section II, we discuss the relevant background required for following the techniques used in the paper. In section III, we present our model and in section IV we introduce different types of couplings and theredy analyze them. Finally, in section V, we present the results of our analysis, before concluding in section VI.

II Background

In this section we present the relevant background of Markovian dynamics and the definitions used in the paper. We also describe the measure of entanglement for two-qubit systems, which will also be used to quantify non-Markovianity of our dynamics.

II.1 Quantum Markovian dynamics

A discrete time stochastic process is called Markovian (Markov chain) if the state of the system at time $t_{n}$ depends solely on the state of the system at time $t_{n-1}$. This concept of Markov chain can be extended to the continuous time stochastic processes as well Rivas et al. (2014). However, generalizing it to quantum dynamics is a difficult task. Numerous prescriptions have been proposed to capture different aspects of quantum Markovianity. Broadly these prescriptions can be classified into two classes: information backflow Breuer et al. (2009); Wißmann et al. (2015); Luo et al. (2012); Bylicka et al. (2017); Chruściński and Rivas (2017); Chakraborty (2018) and completely positive divisibility (CP-divisibility) Breuer et al. (2016); Chruściński and Kossakowski (2012).

Information backflow: The information backflow approach is inspired from the fact that a Markovian dynamics is characterized by unidirectional flow of information from the system to the environment. As for example, in the Lindblad master equation Breuer and Petruccione (2002), the non-negativity of the entropy production rate signifies unidirectional information flow from system to the environment, and thereby, is a signature of Markovianity. A dynamics is called Markovian from the information backflow approach, if some information quantifier decays over time in a monotonic way. Any departure from monotonicity of such quantifier is seen as a backflow of information from the environment, back to the system. Different quantifiers of information like distinguishability of states Breuer et al. (2009), measure of entanglement Rivas et al. (2010), quantum mutual information Luo et al. (2012), etc has also been suggested for this purpose. Each quantifier provides a different definition of Markovianity; all of which, are not in general equivalent. Only recently, there has been attempts to unify all these different definitions Bylicka et al. (2017); Chruściński and Rivas (2017); Chakraborty (2018) to provide a unified approach to information backflow.

CP-divisibility: Any dynamical process, given by completely-positive (CP) trace preserving (TP) map $\Lambda_{t}$ representing evolution up to time $t$ is called CP-divisible if

[TABLE]

where $V_{t,s}$ is CP for any $t\geq s\geq 0$, and $\circ$ denotes composition.

Although the most general description of Markovian dynamics is given by CP-Divisibility Breuer et al. (2016); Chruściński and Kossakowski (2012), for our purpose we consider information backflow, in terms of measure of entanglement i e. the RHP measure, as the description of Markovianity.

II.2 Detecting non-Markovianity through Entanglement

Let us first discuss entanglement measure of two qubit states. The entanglement between two two-level systems (two qubits) can be characterized by the Peres-Horodecki criterion Peres (1996); Horodecki et al. (1996) which states that a two-qubit state $\rho_{\text{as}}$, shared between a system qubit $s$ and an ancilla qubit $a$, is entangled if and only if the partial transpose of this state, i.e. $(\rho_{\text{as}})^{T_{s}}$, is not a positive-semidefinite operator i.e. $(\rho_{as})^{T_{s}}\ngeq 0$. Notably, for a two-qubit entangled state, the operator $(\rho_{\text{as}})^{T_{s}}$ has exactly one negative eigenvalue $\lambda$ Sanpera et al. (1998); Rana and Parashar (2012). Thus $|\lambda|$ may be used as a measure of entanglement for the state $\rho_{\text{as}}$. Formally, the entanglement measure can be defined as follows

[TABLE]

where, $||A||_{1}=\operatorname{Tr}\sqrt{A^{\dagger}A}$ is the trace norm of a matrix $A$. Note that, $E(\rho_{as})$ is nothing but the negativity of the bipartite state $\rho_{as}$ Vidal and Werner (2002). We will use this measure of entanglement as the quantifier for ascertaining Markovianity of the dynamics from the information backflow approach. Using entanglement to detect non-Markovianity was first used by Rivas, Huelga and Plenio in Rivas et al. (2010), and this measure has been so called the RHP measure of non-Markovianity. Following their technique we attach an ancilla to the system, on which a dynamical map $\Lambda_{t}$ is acting. Following the information backflow approach, the dynamical map $\Lambda_{t}$ is called Markovian if E\big{(}(\openone\otimes\Lambda_{t})\big{[}\ket{\Phi^{+}}\bra{\Phi^{+}}\big{]}\big{)} is a non-increasing function of time $t$, where $\ket{\Phi^{+}}\bra{\Phi^{+}}$ is the maximally entangled state, given by

[TABLE]

III The model

In this section, we present our model and discuss the motivation behind choosing it. We also describe the technique in detail, in which non-Markovianity in the system dynamics is detected. We consider two qubits, one of which is called the system ($s$) and the other, the ancilla ($a$). The system qubit is placed in an environment consisting of $N$ non-interacting qubits (see Fig. 1). We take the interaction between the system qubit and the environment in the following form,

[TABLE]

where $\ket{0}$ and $\ket{1}$, respectively represent the ground and excited states of each qubit. The coupling strength $\tilde{g}_{n}(t)$ is in general, a complex number and can also be time-dependent as well as site-dependent, and $\alpha$ is a real parameter with the dimension of frequency. The extra factor $\alpha$ is introduced to make the coupling strengths $\tilde{g}_{n}(t)$ dimensionless. For all practical purposes $\alpha$ can be assumed to be $1$. The free Hamiltonians of the system and the environment are respectively given by,

[TABLE]

It is convenient to work in the interaction picture where we replace the total Hamiltonian $H=H_{s}+H_{e}+\tilde{H}_{\text{se}}\equiv H_{0}+\tilde{H}_{\text{se}}$ by the interaction picture Hamiltonian $H_{\text{se}}(t)=\exp(iH_{0}t/\hbar)\tilde{H}_{\text{se}}\exp(-iH_{0}t/\hbar)$ which reads,

[TABLE]

where $\delta\omega_{n}=\omega_{s}-\omega_{n}$ and $g_{n}(t)=\tilde{g}_{n}(t)e^{-i\delta\omega_{n}}$. Henceforth, our discussion will be based on the Hamiltonian $H_{se}(t)$. We also consider the initial state of the environment to be in the thermal state,

[TABLE]

where p=\big{(}1+e^{-\beta}\big{)}^{-1} and $\beta$ is a positive real parameter which can be identified as the inverse of the temperature $T$ of the environment.

III.1 Motivation behind the model

Here we argue that, our model is a close approximation to the ‘spin bath’ model Hutton and Bose (2004); Breuer et al. (2004) at low temperatures. Note that the Hamiltonian in Eq. (III) can also be written as,

[TABLE]

where $\sigma_{+}=\ket{0}\bra{1}$ and $\sigma_{-}=\ket{1}\bra{0}$. When we compare Eq. (III.1) with the usual Hamiltonian of a spin bath model Hutton and Bose (2004); Breuer et al. (2004); Bhattacharya et al. (2017) in the interaction picture, given by,

[TABLE]

we find that the only difference comes from the $\ket{0}\bra{0}$ factors arising in Eq. (III.1), which are replaced by $\openone$ for the spin bath Hamiltonian. As a result of this difference, the dynamics of the spin bath model is not entirely the same as our model. In the former, an exchange of one quanta of energy takes place between the system and individual environment qubit, when the rest of the environment qubits are allowed to be in any state, whereas in the later, the exchange will only take place when the rest of the environment qubits are in their ground state. This difference, although significant in general, will not play a major role when the state of the environment is close to the ground states, or in other words, temperature of the environment is low. Note that low temperature of environment correspond to values of $p$ in Eq. (8), which are very close to $1$, and this also confirms the fact that for low temperatures $\rho_{e}$ is close to the ground state. Thus we see for low temperatures our model serves as a close approximation to the spin bath model. The main advantage of our model is the fact that our Hamiltonian is easily diagonalizable, and for certain types of couplings, as we discuss later in detail, allows for exact determination of the system dynamics in terms of Kraus operators, for any number of environment qubits.

We also stress that, although our model shows similarity to the spin bath model for low temperatures, we find solutions and analyze the dynamics of our model for any temperature whatsoever. The reason behind this is that our model being analytically solvable for certain types of couplings, allows for an opportunity to exactly solve the dynamics for any number of environment qubits, which is not often the case for systems with large number of spins. Note that, even for the spin-bath Hamiltonian, it is not easy to find the exact solution for non-zero temperature.

III.2 Diagonalizing the Hamiltonian of our model

There are only two non-zero eigenvalues of the Hamiltonian $H_{\text{s}e}(t)$ and they are,

[TABLE]

corresponding to the eigenvectors,

[TABLE]

where $\ket{\xi(t)}_{se}=\ket{0}_{s}\otimes\ket{\beta_{0}}_{e}$ and,

[TABLE]

Thus, the time evolution operator $U(t,0)$ corresponding to the Hamiltonian $H_{\text{s}e}$ is,

[TABLE]

where $\mathcal{T}$ represents time ordering.

The ancilla qubit is used as a probe to characterize the non-Markovianity of the dynamics of the system in the presence of the environment. In order to do so we prepare the system and ancilla qubits in a maximally entangled state $\Ket{\Phi^{+}}$, as given in Eq. (3). Due to the interaction of the system qubit with the environment, the entanglement between the system and the ancilla qubit will evolve with time. The deviation of this time evolution of the entanglement, from monotonic decay is used to establish the non-Markovian character of the dynamics. Note here, that this idea was used by Rivas et al Rivas et al. (2010) to devise a measure of non-Markovianity. In the present paper, we follow this technique to consider the system dynamics to be non-Markovian whenever the entanglement between system and ancilla, as described above, shows non-monotonic behaviour, otherwise we consider the dynamics to be Markovian.

The joint initial state of the system plus ancilla plus environment is of the form,

[TABLE]

which evolves to,

[TABLE]

Therefore, reduced time-evolved system-ancilla state can be calculated by tracing out the environment part,

[TABLE]

IV System - Environment couplings

In this section, we introduce various classes of system-environment coupling, and in each case, we study their effect on the evolution of the system-ancilla joint state. We classify all the couplings into four major classes : (A) when the coupling parameter $g_{n}(t)$ is independent of the site index $n$ (homogeneous) and time-independent; (B) when $g_{n}(t)$ is inhomogeneous but time-independent; (C) when $g_{n}(t)$ is homogeneous but time-dependent, and (D) when $g_{n}(t)$ is inhomogeneous and time-dependent. For each class, we calculate the entanglement of the time evolved state of system-ancilla, and thereby try to characterize the non-Markovian behaviour of the system dynamics. Henceforth, we assume $\alpha$ to be $1$.

IV.1 Homogeneous and time-independent coupling

We have here the simplest situation, where the coupling of the system with all the environment qubits are uniform and time-independent i.e. $g_{n}(t)=g$, a constant. As a result, the non-zero eigenvalues of the Hamiltonian, as given in Eq. (11), takes the form $\mathcal{E}_{\pm}(t)=\pm\mathcal{E}=\pm\hbar\sqrt{N}~{}|g|\equiv\hbar\omega_{0}$, where $\omega_{0}=\sqrt{N}~{}|g|$ is a constant with the dimension of frequency. The time-evolution operator $U(t,0)$ is of the form,

[TABLE]

Using the above form and the form of $\rho_{e}$ given in Eq. (8), we find the Kraus operators $K_{mn}(t)$ of system dynamics, which are defined in the following way,

[TABLE]

where the $N^{2}$ Kraus operators are given by,

[TABLE]

where $m,n=1,\dots,N$, $\log x$ refers to $\log_{2}x$, and $s_{n}$ is the number of $1$’s in the binary equivalent of $n$. For example, if $n=6$, then the binary equivalent of $n$ is $110$. Therefore $s_{n}=2$.

We then find time evolved state of the system-ancilla, using Eqs. (15), (16) and (17),

[TABLE]

where $\kappa_{0}=p^{N-1}\sin^{2}(\omega_{0}t)$ and $\delta_{0}=2p^{N-1}\sin^{2}\left(\frac{\omega_{0}t}{2}\right)$. The only possible negative eigenvalue of $[\rho_{\text{as}}(t)]^{T_{s}}$, if any, is of the form,

[TABLE]

We present the plot of $E(\rho_{a}s(t))=\lambda(t)$ versus time, later in the Result section.

IV.2 Inhomogeneous and time-independent coupling

Consider a system, where a single two-level system (perhaps an ion as an impurity) is placed in a spin lattice. The lattice sites, closest to the impurity interacts very strongly with the system, while, as we go away from the impurity site, the strength of interaction becomes weaker and weaker. In such cases the interaction parameter $g_{n}(t)$ is in general inhomogeneous, but there is no explicit time dependence. Therefore, $g_{n}(t)=g_{n}$. Hence, $\mathcal{E}_{\pm}(t)=\pm\mathcal{E}=\hbar\omega$ in Eq. (11) are also time-independent. Note, in this case also $\omega=\sqrt{\sum_{n=1}^{N}|g_{n}|^{2}}$ is a constant with dimensions of frequency. Following this, the analysis is same as in the last subsection. As a result, the evolution operator $U(t,0)$, the time evolved state $\rho_{\text{as}}(t)$ and the only possible negative eigenvalue, if any, of $[\rho_{\text{as}}(t)]^{T_{s}}$ for this case are of the same forms as in Eqs (IV.1), (20), (IV.1) and (22) respectively, except for $\omega_{0}$, in appropriate places, replaced by $\omega$.

IV.3 Homogeneous and time-dependent coupling

So far we have considered only couplings which are independent of time. In this section, we consider time-dependent and homogeneous couplings. We take an arbitrary real function of time, which is independent of site index $n$ i.e. $g_{n}(t)=g(t)$. Note that our coupling operator between system and individual environment qubit, as given in Eq. (III), is of the form $\sigma_{+}\otimes\sigma_{-}+\sigma_{-}\otimes\sigma_{+}$ which can also be expressed as $\sigma_{x}\otimes\sigma_{x}+\sigma_{y}\otimes\sigma_{y}$. Thus, our system-environment coupling is a special case of the $XY$ coupled Hamiltonian. Such coupling with time-dependent coefficients have been used to show non-trivial entanglement dynamics Sadiek et al. (2010); Bortz and Stolze (2007).

Fortunately, the Hamiltonians $H_{se}(t)$ in this case commutes at different times, which makes the analysis similar to the one in section IV.2. The only difference being the non-zero eigenvalues, in Eq. (11), to be of the form $\mathcal{E}_{\pm}(t)=\pm\hbar\sqrt{N}|g(t)|=\pm\mathcal{E}(t)$, which is no longer constant in time. The whole treatment of the dynamics of the system and the ancilla remains the same if we replace $\omega_{0}$ and $\omega_{0}t$, in Eqs (IV.1), (20), (IV.1) and (22), by $\sqrt{N}~{}|g(t)|$ and $\Omega(t)$, respectively, where,

[TABLE]

IV.4 Inhomogeneous and time-dependent coupling

The most general class of coupling $g_{n}(t)$ is when it depends on both the site $n$ and time $t$. The interaction Hamiltonian in such a situation does not commute at different times and this makes the calculation for solving the dynamics difficult. However, we can use numerical methods to simulate the time-evolution and get the solution for $\rho_{\text{as}}(t)$.

One can obtain the following results analytically before starting the simulation part.

Analytical Part: Two of the eigenvalues of the Hamiltonian in Eq. (III) are non-zero, as given in Eq. (11). The remaining $(2^{N+1}-2)$ of the eigenvalues are zero. A possible choice for these null space eigenvectors are found in the following way:

Step I: We first feed the eigenvectors (corresponding to non-zero eigenvalues) given in Eq. (12) as rows of a matrix $A$. Note, $A$ is a $2\times 2^{(N+1)}$ matrix

Step II: By row reduction method Hoffman and Kunze , we find out a basis $\mathcal{B}$ for the Null space of $A$. Note that $\mathcal{B}$ is not necessarily ortho-normal.

Simulation Part: Obtaining an orthonormal basis $\mathcal{B}^{\prime}$ from $\mathcal{B}$ analytically is a challenging job. We, therefore resort to numerical techniques for this case.

Step I: From $\mathcal{B}$, using Gram-Schmidt Orthonormalization procedure Hoffman and Kunze , we find an orthonormal basis $\mathcal{B}^{\prime}$. Note, $\mathcal{B}^{\prime}$ forms the set of eigenvectors of the Hamiltonian, corresponding to zero eigenvalues.

Step II: As, the eigenvectors are time-dependent, the Hamiltonian is not different-time commuting. Hence, the evolution operator may be found numerically from the following expression,

[TABLE]

where $\mathcal{T}$ represents time ordering.

Step III: We evolve the initial ancilla -system-environment state $\rho_{ase}(0)$ by the unitary operator $U_{ase}(t,0)=\mathbb{I}_{a}\otimes U_{se}(t,0)$ and get the time evolved state $\rho_{ase}(t)$.

Step IV: We trace out the environment from $\rho_{ase}(t)$ and get $\rho_{as}(t)=Tr_{e}[\rho_{ase}(t)]$. We then evaluate our entanglement measure $E(t)$ given in Eq. (2), on $\rho_{as}(t)$ and plot it as a function of time.

V Results

In this section, we show that some of the classes of the couplings that we have considered, always results in non-Markovian dynamics. However, there are also some classes for which we can tune the parameters and find a transition from non-Markovian to Markovian dynamics. In order to do so, we plot the entanglement dynamics between the system and ancilla, for each class, as a function of time and observed if there is any departure from monotonicity with time. As mentioned in Sec III, this technique helps in characterizing non-Markovianity present in the system dynamics. In Fig. 2, we present entanglement dynamics for different classes of the system-environment coupling, considered in the previous section. We now present our findings for each class of system-environment coupling.

Section IV.1 and IV.2 : In Figs. 2(a) and 2(b), we plot the entanglement as a function of time for the homogeneous time-independent, and the inhomogeneous time-independent couplings, respectively, i.e, $g_{n}(t)=g$ and $g_{n}(t)=g_{n}$, respectively. Note here, that $g$ and $g_{n}$, for all values of $n$, are arbitrary complex functions. Analytic calculations for the entanglement measure, as given in Eq. ((22)) suggests a periodic behaviour for both the classes, which can also be seen in Figs. 2(a) and 2(b). As a result, we conclude in both of these classes of couplings, the dynamics is always non-Markovian.

Section IV.3 : We consider homogeneous and time-dependent couplings and find that if $g(t)$ is some polynomial function of $t$, we will get $\Omega(t)$ as a polynomial function of t. This gives rise to a periodic function $\lambda(t)$. As a result, the dynamics is non-Markovian, in general. Interestingly, it was recently pointed out Dive et al. (2015), that the Hamiltonian dialation obtained from the Choi state of the system dynamics, diverges whenever the system dynamics is time-independent Markovian. Although, in our case, the system-environment Hamiltonian (equivalent to the dialated Hamiltonian)is not related to the Choi state of the system dynamics in the same way as in Dive et al. (2015). As a result the conclusion of Dive et al. (2015) differs from ours. If $g(t)=\exp(-\gamma t)$ then non-Markovianity can be witnessed if the real part $\gamma_{r}$ of $\gamma$ fails to be positive or violates the inequality $\alpha\sqrt{N}\geq\gamma_{r}\pi$. Fig. 2(c) shows the entanglement vs time plot for two values of $\gamma_{r}$; one of which violates the above mentioned inequality. We also consider the case $g(t)=t^{-\delta}$, and show that the dynamics is non-Markovian for $\delta=0.01$.

Section IV.4 : For inhomogeneous time-dependent couplings, the dynamics can be made both non-Markovian and Markovian by choosing the strength of the coupling appropriately. We consider two special cases of inhomogeneous time-dependent coupling: $(i)$ $g_{n}(t)=e^{-\gamma_{1}nt}$, and $(ii)$ $g_{n}(t)=t^{-n\gamma}$. For simplicity, we have assumed $\alpha=1$ and set the temperature of the environment to be zero i.e. $\beta\to\infty$. For analyzing coupling $(i)$, we plot the system-ancilla entanglement measure as a function of time in Figures 2(d), and 2(e), for different values of the coupling parameter $\gamma_{1}$, at a fixed value of $N$ i.e. number of qubits in the environment. In Figures 2(d), and 2(e), monotonically decreasing entanglement values show signs of Markovianity and non-monotonic decay are evidence of non-Markovianity. As expected, increasing the coupling parameter $\gamma_{1}$ i.e., decreasing coupling strength, leads to the transition from non-Markovian to Markovian dynamics. The figures also show an interesting feature that, after sufficient time, the entanglement in the system-ancilla state saturates to fixed values irrespective of their Markovian or non-Markovian nature. This feature can be signs of possible equilibration of the system ancilla state. For coupling $(ii)$, the dynamics showed non-Markovianity for $\gamma=0.3$, as shown in Fig 2(f).

Next, for coupling $(i)$ i.e. $g_{n}(t)=e^{-\gamma_{1}nt}$, we find the extremal values of $\gamma_{1}$ for which non-Markovianity is witnessed. These extremal values serve as transition parameters from non-Markovianity to Markovianity. On plotting these transition values as a function of $N$ (Fig. 3) it appears that a saturation is reached as $N$ is increased.

VI Conclusion

In this paper, we have addressed the question of how non-Markovianity of a dynamics changes with the interaction between the system and the environment and also with size of the environment. We have taken a simple model constituting of a few qubits, which can also be seen as a close approximation to the spin bath model, for low temperatures. Even in this minimalistic scenario, we were able to find a transition from non-Markovian to Markovian dynamics by tuning the system-environment interaction. This is somewhat counterintuitive as it is generally conceived that for having Markovian dynamics the bath/environment should have infinite degrees of freedom, although there are exceptions Pang et al. (2017). We also found, in our model, that if the interaction Hamiltonian is time-independent, the dynamics is always non-Markovian, irrespective of the size of environment. Note that in this scenario for a general interaction Hamiltonian in the weak coupling limit, we generally get to see Markovian dynamics only for a very large size of bath; at least in the case of harmonic oscillator bath. The present scenario is different as we have considered spin environment and the system-environment interaction is of very specific type. In the case of site-independent interaction, polynomial forms and certain cases of exponential forms of interaction show non-Markovianity. Lastly, we study time-dependent and site-dependent interaction for certain forms of system-environment coupling. In this last case, we also saw a transition from non-Markovian to Markovian regime. Interestingly, the transition values appear to saturate to a certain value depending on the initial temperature of the environment, as the number of environment qubits increases. Examining this type of spin environment is recently drawing some amount of interest Prokof’ev and Stamp (2000). Studies on similar lines was also done recently in Bhattacharya et al. (2017), where an analysis of a qubit system interacting with a sea of spins was given. A number of questions arise from the present study: whether such a transition can be found by considering more general forms of interaction, what happens if along with system-environment interaction there is some interaction present among the environment particles themselves, etc. One may also find it useful to check, whether the aforesaid saturation of transition parameters is a general feature of interaction that exhibits non-Markovian to Markovian transition.

Acknowledgements.

The authors would like to thank R. Lo Franco for helping in revising the earlier version of the manuscript as well as providing the references Franco (2016); Bellomo et al. (2012); González-Gutiérrez et al. (2016a); Xu et al. (2013); Franco et al. (2013); Bellomo et al. (2011); Man et al. (2015b); Franco (2015). The authors would also like to thank C. González-Gutiérrez for providing reference González-Gutiérrez et al. (2016b), F. Plastina and T. J. G. Apollaro for providing references Apollaro et al. (2011); Lorenzo et al. (2013). SC, AM, and SG gratefully acknowledge discussions with Samyadev Bhattacharya and Chiranjib Mukhopadhyay.