This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Influence of initial states on memory effects: A study of early-time superradiance

S. C. Hou housc@dlmu.edu.cn School of Science, Dalian Maritime University, Dalian 116026, China    G. Q. Shuai School of Science, Dalian Maritime University, Dalian 116026, China    X. Y. Zhang School of Science, Dalian Maritime University, Dalian 116026, China    J. Shen Department of Criminal Technology, Liaoning Police College, Dalian 116036, China    X. X. Yi yixx@nenu.edu.cn Center for Quantum Sciences and School of Physics, Northeast Normal University, Changchun 130024, China
Abstract

The initial states of a quantum system can significantly influence its future dynamics, especially in non-Markovian quantum processes due to the environmental memory effects. Based on a previous work of ours, we propose a method to quantify the memory effects of a non-Markovian quantum process conditioned on a particular system initial state. With this method, we analytically study the early-time memory effects of a superradiance model consisting of NN atoms (the system) interacting with a single-mode vacuum cavity (the environment) with two types of initial states: the Dicke states and the factorized identical states. We find that the radiation intensity are closely related to the memory effects, and correspondingly, the enhancement of memory effects (from independent radiation to collective radiation) is important for the degree of superradiance, especially for the Dicke states. Furthermore, numerical results in other regimes show that the characteristics of memory effects and superradiance in longer-time dynamics can be reflected through its early-time dynamics.

I Introduction

The initial state of a quantum system can significantly influence its future dynamics [1, 2], especially in non-Markovian quantum processes due to the environmental memory effects. One trivial example is that if the system is initially in a steady state in a non-Markovian quantum process, it can hardly exhibit any non-Markovian features afterward, such as the nonmonotonic behaviors of energy and information flows [3, 4, 5, 6]. More intriguing phenomena may emerge when a system consists of a collection of subsystems, such that the properties of its initial state, such as entanglement and coherence, may significantly influence its future dynamics. A well-known example is the concept of superradiance [7, 9, 8, 10, 11, 12, 13, 14] introduced by Dicke in 1954, where the emission intensity from an ensemble of atoms interacting with a common electromagnetic field can be enhanced compared with that from independent atoms. It is also well-known that the superradiance behaviors are highly relevant to the initial states of the system. In recent years, superradiance has received a large amount of attention due to its theoretical significance and potential applications [16, 15, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31]. Under certain approximations, such as a coarse-grained timescale, the superradiance process could be regarded as Markovian [10, 11, 12, 13, 14], whereas, it is intrinsically non-Markovian. With the advances in theories and technologies, understanding the non-Markovian dynamics of superradiance becomes more demanding [28, 29, 30, 31]. In a recent work [31], the author shows that non-Markovian memory effects play an important role in superradiance beyond retardation, featuring the quadratic dynamics in the early-time (Zeno) regime, whereas, the memory effects are not quantitatively evaluated. In view of the significant influences of the system initial states on non-Markovian quantum processes, especially, a superradiance process, some interesting questions arise. For example, how to quantitatively evaluate the memory effects of a quantum process conditioned on a particular system initial state? Are there quantitative relations between the memory effects and the superradiance characteristics in the early-time regime? What is the role of the (very weak) memory effects in a superradiance process that could be well approximated as Markovian?

In recent years, a number of measures or manifestations of non-Markovianity were proposed to characterize the memory effects which are often connected to nonmonotonic behaviors [3, 4, 5, 6]. For example, the well-known BLP measure [32] and the RHP measure [33] use the increases of distinguishability and entanglement, respectively, to measure the non-Markovianity. The nonmonotonic behaviors originate from the noncomplete positivity of the intermediate dynamical maps in non-Markovian processes. These measures generally apply to quantum processes where the evolutions start at a fixed time t0t_{0} (usually t0=0t_{0}=0 for simplicity). However, a key distinction of the non-Markovian process is the existence of t0t_{0} in the generator of the system’s time-local equation, i.e., ρ˙S=(tt0)ρS\dot{\rho}_{S}=\mathcal{L}(t-t_{0})\rho_{S} [34]. Thus some features of non-Markovianity might not be characterized with only evolutions starting from a fixed time t0t_{0}. For example, in the early-time regime or a near-Markovian regime, the spontaneous radiation of an atom could be described by a time-local equation with a positive (time-dependent) damping rate, which looks like a time-dependent Markovian master equation (except t0t_{0} in its generator) [35]. Thus the nonmonotonic behaviors do not occur in these regimes, implying weaker memory effects. Meanwhile, the non-Markovianity measures based on the noncomplete positivity may give zero memory effects. On the other hand, the previously proposed non-Markovianity measures generally do not depends on a particular system initial state, but deal with a quantum process where all the potential initial states are considered. Therefore, an initial-state-dependent measure which is able to characterize weak memory effects is demanding.

In this paper, we propose such a method developed from a previous work [35] of ours. The non-Markovianity measure in Ref. [35] is based on the inequality of completely positive dynamical maps T(t2,t0)T(t2,t1)T(t1,t0)T(t_{2},t_{0})\neq T(t_{2},t_{1})T(t_{1},t_{0}) and connected with the validity of the Markov approximation. By acting the inequality on a particular initial state ρS(t0)\rho_{S}(t_{0}), the evolutions starting at t0t_{0} and t1t_{1} have the same system states at t1t_{1} but different histories before t1t_{1}. Thus the memory effects (past-future dependence) could be understood by the difference between the final states ρS(t2)\rho_{S}(t_{2}) and ρS(t2)\rho_{S}^{\prime}(t_{2}). Based on the above interpretation, we suggest measuring the memory effects conditioned on the initial state ρS(t0)\rho_{S}(t_{0}) by the maximal difference of ρS(t2)\rho_{S}(t_{2}) and ρS(t2)\rho_{S}^{\prime}(t_{2}) in a time interval of interest. Using this method, we calculate the memory effects (as well as the radiation characteristics) of a superradiance model with two types of initial states: the Dicke states and factorized identical states in different regimes. The model describes NN two-level atoms interacting with a single-mode cavity initially in vacuum. We give analytic expressions of the characteristics of memory effects (and superradiance) in the early-time regime. Moreover, we extend our work to longer time intervals in the near-Markovian regime, and to the whole radiation lifetime in a strongly non-Markovian regime with numerical results.

The main findings are as follows. 1. The radiation intensity is closely related to the memory effects, and correspondingly, the degree of superradiance is closely related to the degree of memory-effect-enhancement (from independent radiation to collective radiation), especially for the Dicke states. 2. The entanglement in the Dicke states and the single-atom coherence in the factorized identical states are necessary for the superradiance and important for the memory-effect-enhancement. 3. The (change of) environmental photon number is a main physical source of memory effects for our model. 4. The characteristics of memory effects and superradiance in long-time dynamics (especially for the near-Markovian regime) can be reflected through its early-time dynamics.

The paper is organized as follows: In Sec. II, we propose our method to evaluate the initial-state-dependent memory effects in a quantum process. In Sec. III, we obtain the early-time solution of NN two-level atoms interacting with a vacuum cavity, with which we give analytic expressions for the value of memory effects, the cavity photon number, the degree of superradiance, and so on. In Sec. IV, the influences of two types of initial states on the memory effects as well as the superradiance are calculated and analyzed using the early-time results. Our work is extended to the near-Markovian regime and a strongly non-Markovian regime in Sec. V. At last, conclusions and discussions are presented in Sec. VI.

II Initial-state-dependent memory effects

In this section, we first review the measure of non-Markovianity in Ref. [35], which quantifies the memory effects (the past-future dependence) in a quantum process. Using its physical interpretations, we then propose a method to quantify the memory effects in a quantum process conditioned on a particular initial state of the system. The object of study in Ref. [35] is a quantum process of an open quantum system described by the total Hamiltonian

H=HS+HE+HSE,\displaystyle H=H_{S}+H_{E}+H_{SE}, (1)

an arbitrary system initial state ρS(tI)\rho_{S}(t_{I}) and a fixed initial state of the environment ρE(tI)\rho_{E}(t_{I}) (independent of the system). Here tIt_{I} is the initial time of an evolution which is also arbitrary. The initial condition of an evolution is

ρSE(tI)=ρS(tI)ρE(tI)\displaystyle\rho_{SE}(t_{I})=\rho_{S}(t_{I})\otimes\rho_{E}(t_{I}) (2)

which means that the system is initially uncorrelated with the environment before the evolution. In general, ρE(tI)\rho_{E}(t_{I}) is governed by [6]

ρE(tI)=𝒯ei0tIHE(τ)𝑑τρE(0).\displaystyle\rho_{E}(t_{I})=\mathcal{T}e^{-i\int_{0}^{t_{I}}H_{E}(\tau)d\tau}\rho_{E}(0). (3)

Typically, one deals with a time-independent environment Hamiltonian HEH_{E} and assumes a steady state of HEH_{E} as the environmental initial state (e.g., a thermal state). Then, the initial condition of an evolution is

ρSE(tI)=ρS(tI)ρE\displaystyle\rho_{SE}(t_{I})=\rho_{S}(t_{I})\otimes\rho_{E} (4)

for any tIt_{I}. Since ρS(tI)\rho_{S}(t_{I}) and tIt_{I} is arbitrary, the measure of non-Markovianity (memory effects) in Ref. [35] is determined by HH and ρE(tI)\rho_{E}(t_{I}) and the time interval of interest.

In a non-Markovian quantum process, the meaning of a dynamical map given by ρS(t2)=ε(t2,t1)ρS(t1)\rho_{S}(t_{2})=\varepsilon(t_{2},t_{1})\rho_{S}(t_{1}) is ambiguous unless ρSE(t1)\rho_{SE}(t_{1}), or from another perspective, the initial time tIt_{I} is specified (tIt1t_{I}\leqslant t_{1}). For clarity, we define T(t2,t1)T(t_{2},t_{1}) as a memoryless dynamical map that transfers ρS(t1)\rho_{S}(t_{1}) to ρS(t2)\rho_{S}(t_{2}), where t1t_{1} is the initial time of the evolution (tI=t1t_{I}=t_{1}) [35], i.e.,

ρS(t2)\displaystyle\rho_{S}(t_{2}) =\displaystyle= T(t2,t1)ρS(t1)\displaystyle T(t_{2},t_{1})\rho_{S}(t_{1}) (5)
=\displaystyle= TrE[U(t2,t1)ρS(t1)ρE(t1)U(t2,t1)],\displaystyle\text{Tr}_{E}[U(t_{2},t_{1})\rho_{S}(t_{1})\otimes\rho_{E}(t_{1})U(t_{2},t_{1})^{{\dagger}}],

where ρE(t1)\rho_{E}(t_{1}) is fixed, ρS(t1)\rho_{S}(t_{1}) is arbitrary and U(t2,t1)=𝒯eit1t2H(τ)𝑑τU(t_{2},t_{1})=\mathcal{T}e^{-i\int_{t_{1}}^{t_{2}}H(\tau)d\tau} in general. Particularly, when HH is time-independent and ρE(tI)=ρE\rho_{E}(t_{I})=\rho_{E} is a steady state of HEH_{E}, Eq. (5) is simplified to

ρS(t2)\displaystyle\rho_{S}(t_{2}) =\displaystyle= T(t2,t1)ρS(t1)\displaystyle T(t_{2},t_{1})\rho_{S}(t_{1}) (6)
=\displaystyle= TrE[eiH(t2t1)ρS(t1)ρEeiH(t2t1)]\displaystyle\text{Tr}_{E}[e^{-iH(t_{2}-t_{1})}\rho_{S}(t_{1})\otimes\rho_{E}e^{iH(t_{2}-t_{1})}]

such that T(t2,t1)=T(t2t1,0)T(t_{2},t_{1})=T(t_{2}-t_{1},0) [34]. The dynamical map TT is trace-preserving and completely positive and called a universal dynamical map (UDM) that is independent of the state it acts upon [36].

Let t0t1t2t_{0}\leqslant t_{1}\leqslant t_{2}, a (time-dependent) quantum process is Markovian if the dynamical maps (denoted by ΛM\Lambda_{M}) satisfy the divisibility condition

ΛM(t2,t0)=ΛM(t2,t1)ΛM(t1,t0),\displaystyle\Lambda_{M}(t_{2},t_{0})=\Lambda_{M}(t_{2},t_{1})\Lambda_{M}(t_{1},t_{0}), (7)

where each dynamical map is uniquely defined and a UDM. Remark that the dynamical map ΛM(t2,t1)\Lambda_{M}(t_{2},t_{1}) is a UDM if and only if it is induced by

ρS(t2)\displaystyle\rho_{S}(t_{2}) =\displaystyle= ΛM(t2,t1)ρS(t1)\displaystyle\Lambda_{M}(t_{2},t_{1})\rho_{S}(t_{1}) (8)
=\displaystyle= TrE[U(t2,t1)ρS(t1)ρE(t1)U(t2,t1)],\displaystyle\text{Tr}_{E}[U(t_{2},t_{1})\rho_{S}(t_{1})\otimes\rho_{E}(t_{1})U(t_{2},t_{1})^{{\dagger}}],

where ρE(t1)\rho_{E}(t_{1}) is fixed and ρS(t1)\rho_{S}(t_{1}) is arbitrary [36]. In an open quantum system, the condition ρSE(t)=ρS(t)ρE(t)\rho_{SE}(t)=\rho_{S}(t)\otimes\rho_{E}(t) [ρE(t)\rho_{E}(t) does not dependent on the system state] may not be satisfied exactly for t>tIt>t_{I}. Thus the dynamics of an exact open quantum system is typically not Markovian [36]. However, ρSE(t)ρS(t)ρE(t)\rho_{SE}(t)\approx\rho_{S}(t)\otimes\rho_{E}(t) can be a good approximation where the correlation between the system and the environment does not affect the system’s dynamics so much [36]. It is observed that the Markovian dynamical map ΛM(t2,t1)\Lambda_{M}(t_{2},t_{1}) is unique and does not depend on the initial time of an evolution, e.g., tI=t1t_{I}=t_{1} or tI<t1t_{I}<t_{1}. Therefore, ΛM(t2,t1)=T(t2,t1)\Lambda_{M}(t_{2},t_{1})=T(t_{2},t_{1}) according to the definition of the memoryless dynamical map TT. Then, the Markovian divisibility condition (7) can also be expressed in terms of the memoryless dynamical map TT by

T(t2,t0)=T(t2,t1)T(t1,t0),\displaystyle T(t_{2},t_{0})=T(t_{2},t_{1})T(t_{1},t_{0}), (9)

whose violation is a sign of non-Markovianity. Unlike the dynamical maps used in some non-Markovianity measures [33, 37], all the dynamical maps TT are completely positive. The violation of Eq. (9) is manifested by the inequality

T(t2,t0)T(t2,t1)T(t1,t0).\displaystyle T(t_{2},t_{0})\neq T(t_{2},t_{1})T(t_{1},t_{0}). (10)

The physical meaning of Eq. (10) can be explained with Fig. 1. Let the left-hand and the right-hand sides of Eq. (10) act on a system initial state ρS(t0)\rho_{S}(t_{0}). On the left-hand side of Eq. (10), ρS(t0)\rho_{S}(t_{0}) is mapped to ρS(t2)\rho_{S}(t_{2}) by T(t2,t0)T(t_{2},t_{0}) in evolution A. On the right-hand side, ρS(t0)\rho_{S}(t_{0}) is first mapped to ρS(t1)\rho_{S}(t_{1}) by T(t1,t0)T(t_{1},t_{0}) in evolution B. Then, as the initial state of evolution C, ρS(t1)\rho_{S}(t_{1}) is mapped to ρS(t2)\rho^{\prime}_{S}(t_{2}) by T(t2,t1)T(t_{2},t_{1}), which means that the initial condition of evolution C is ρSE(t1)=ρS(t1)ρE(t1)\rho_{SE}(t_{1})=\rho_{S}(t_{1})\otimes\rho_{E}(t_{1}) with fixed ρE(t1)\rho_{E}(t_{1}) defined by Eq. (3) or a steady one ρE\rho_{E}. At moment t1t_{1}, evolution A and C have the same system state ρS(t1)\rho_{S}(t_{1}) but different histories: evolution A has a history in time interval [t0,t1][t_{0},t_{1}], which is encoded in ρSE(t1)=U(t1,t0)ρS(t0)ρE(t0)U(t1,t0)\rho_{SE}(t_{1})=U(t_{1},t_{0})\rho_{S}(t_{0})\otimes\rho_{E}(t_{0})U(t_{1},t_{0})^{{\dagger}}, while evolution C (starting at t1t_{1}) has no history before t1t_{1}. Therefore, ρS(t2)ρS(t2)\rho_{S}(t_{2})\neq\rho^{\prime}_{S}(t_{2}) is an evidence that the future state (after t1t_{1}) of the system depends on its history (in [t0,t1][t_{0},t_{1}]) in this quantum process. This is the physical meaning of the memory effects in this paper and [35]. Otherwise, if the divisibility (9) holds, the process is Markovian and ρS(t2)=ρS(t2)\rho_{S}(t_{2})=\rho^{\prime}_{S}(t_{2}) for any ρS(t0)\rho_{S}(t_{0}).

Refer to caption
Figure 1: Physical interpretation of Eq. (10) as memory effects. At the moment t1t_{1}, the system states in evolution A and C are the same. Recall that evolution A has a history from t0t_{0} to t1t_{1} while evolution C does not have any history before t1t_{1}. Therefore, ρS(t2)ρS(t2)\rho_{S}(t_{2})\neq\rho^{\prime}_{S}(t_{2}) is an evidence that the future (after t1t_{1}) state of the system depends on its history (from t0t_{0} to t1t_{1}) in a quantum process.

On the other hand, the inequality could be understood by focusing on the change of the environment state between evolution B and C. At the end of evolution B, ρSE(t1)=U(t1,t0)ρS(t0)ρE(t0)U(t1,t0)\rho_{SE}(t_{1})=U(t_{1},t_{0})\rho_{S}(t_{0})\otimes\rho_{E}(t_{0})U(t_{1},t_{0})^{{\dagger}}. Then, at the beginning of evolution C, the environment is initialized by T(t2,t1)T(t_{2},t_{1}) such that ρSE(t1)ρS(t1)ρE(t1)\rho_{SE}(t_{1})\rightarrow\rho_{S}(t_{1})\otimes\rho_{E}(t_{1}), where ρE(t1)\rho_{E}(t_{1}) is independent of the system. The information of the system’s history in [t0,t1][t_{0},t_{1}] is erased by T(t2,t1)T(t_{2},t_{1}). In contrast, such an initialization never happens in evolution A. Therefore, ρS(t2)ρS(t2)\rho_{S}(t_{2})\neq\rho^{\prime}_{S}(t_{2}) signifies that the environment (as well as the correlations between the system and environment [36]) remembers the system’s history and this memory can influence the future of the system.

In Ref. [35], we define the maximal difference of T(t2,t0)T(t_{2},t_{0}) and T(t2,t1)T(t1,t0)T(t_{2},t_{1})T(t_{1},t_{0}) as the non-Markovianity where t1t_{1} and t2t_{2} are optimized. Based on the dynamical maps TT, all the initial states are potentially considered, thus the measure of non-Markovianity in Ref. [35] does not depend on the initial state of the system. In this paper, our goal is to evaluate the influence of the system initial state on the memory effects in a quantum process. Using the physical interpretations discussed above, we define the value of memory effects conditioned on ρS(t0)\rho_{S}(t_{0}) as the maximal trace distance between ρS(t2)\rho_{S}(t_{2}) and ρS(t2)\rho^{\prime}_{S}(t_{2}):

NM[ρS(t0)]=maxt1,t212ρS(t2)ρS(t2).\displaystyle N_{M}[\rho_{S}(t_{0})]=\max_{t_{1},t_{2}}\frac{1}{2}\|\rho_{S}(t_{2})-\rho^{\prime}_{S}(t_{2})\|. (11)

Here A=Tr(AA)||A||=\text{Tr}(\sqrt{A^{\dagger}A}) is the trace norm of an operator AA. Assuming a fixed t0t_{0} for convenience, NM[ρS(t0)]N_{M}[\rho_{S}(t_{0})] could be calculated by optimizing t1t_{1} and t2t_{2} in a time interval [t0,t][t_{0},t] of interest where t0t1t2tt_{0}\leqslant t_{1}\leqslant t_{2}\leqslant t. For example, the time interval might be [t0,][t_{0},\infty] or [t0,t0+τ][t_{0},t_{0}+\tau], where τ\tau is finite (as done later in this paper). Similar to the measure in Ref. [35], 0NM[ρS(t0)]10\leqslant N_{M}[\rho_{S}(t_{0})]\leqslant 1 is satisfied due to the properties of the trace distance.

Using Eq. (11), the influence of initial states on the memory effects of a quantum process could be quantitatively evaluated. Note that NM[ρS(t0)]>0N_{M}[\rho_{S}(t_{0})]>0 is a sufficient condition for the inequality Eq. (10), but not a necessary one. The theoretical calculation and experimental observation of Eq. (11) might be easier than those in Ref. [35] since the determination of the dynamical maps TT is not compulsory. Moreover, as discussed in Ref. [35], Tr[Pρ(t2)]Tr[Pρ(t2)]\text{Tr}[P\rho(t_{2})]\neq\text{Tr}[P\rho^{\prime}(t_{2})] is a sufficient condition for ρS(t2)ρS(t2)\rho_{S}(t_{2})\neq\rho_{S}^{\prime}(t_{2}), where PP is an operator of a physical quantity. Therefore, when the full information of the system state is not easily obtainable, ΔP[ρS(t0)]=|Tr[Pρ(t2)]Tr[Pρ(t2)]|\Delta P[\rho_{S}(t_{0})]=|\text{Tr}[P\rho(t_{2})]-\text{Tr}[P\rho^{\prime}(t_{2})]| (maxt1,t2ΔP[ρS(t0)]\max_{t_{1},t_{2}}\Delta P[\rho_{S}(t_{0})]) could be used as an evidence (manifestation) of the memory effects for simplicity and intuitiveness both theoretically and experimentally.

III Early-time superradiance

III.1 Theoretical model

We consider a fundamental model that describes NN two-level atoms (the system, denoted by “SS”) interacting with a cavity (the environment, denoted by “EE” ) initially in a vacuum state. The Hamiltonian is given by

H\displaystyle H\! =\displaystyle= HS+HE+HI\displaystyle\!H_{S}+H_{E}+H_{I} (12)
=\displaystyle= ωAn=1Nσn+σn+ωBbb+n=1Ng(σn+b+σnb).\displaystyle\!\omega_{A}\!\sum_{n=1}^{N}\!\sigma_{n}^{+}\sigma_{n}^{-}\!+\!\omega_{B}b^{{\dagger}}b\!+\!\sum_{n=1}^{N}g(\sigma_{n}^{+}b+\sigma_{n}^{-}b^{{\dagger}}).

Here HSH_{S} and HEH_{E} represent NN noninteracting two-level atoms and a single-mode electromagnetic field in the cavity, respectively. HI=n=1Ng(σn+b+σnb)H_{I}=\sum_{n=1}^{N}g(\sigma_{n}^{+}b+\sigma_{n}^{-}b^{{\dagger}}) describes the interactions between the atoms and the cavity with the rotating wave approximation, and the coupling strength gg is a real constant. The lowering and raising operators for the nnth atom is defined as σn=|gne|n\sigma_{n}^{-}=|g\rangle_{n}\langle e|_{n} and σn+=|eng|n\sigma_{n}^{+}=|e\rangle_{n}\langle g|_{n}.

The initial condition of the model is assumed to be ρSE(t0)=ρS(t0)|00|\rho_{S\!E}(t_{0})=\rho_{S}(t_{0})\otimes|0\rangle\langle 0| where ρS(t0)\rho_{S}(t_{0}) is arbitrary and |0|0\rangle is the vacuum state of the cavity. For simplicity, we use t0=0t_{0}=0 in the remainder of this paper without loss of generality. The density matrix of the composite system ρSE\rho_{S\!E} is described by the master equation,

ρ˙SE\displaystyle\dot{\rho}_{S\!E} =\displaystyle= i[HS+HE+HI,ρSE]\displaystyle-i[H_{S}+H_{E}+H_{I},\rho_{S\!E}] (13)
+γ(bρSEb12bbρSE12ρSEbb),\displaystyle+\gamma(b\rho_{S\!E}b^{{\dagger}}-\frac{1}{2}b^{{\dagger}}b\rho_{S\!E}-\frac{1}{2}\rho_{S\!E}b^{{\dagger}}b),

where γ\gamma is the dissipation strength of the cavity. Assume the cavity is in resonance with the atoms (ωA=ωB\omega_{A}=\omega_{B}), then the density matrix ρSE\rho_{S\!E} in the interaction picture is described by

ρ˙SE\displaystyle\dot{\rho}_{S\!E} =\displaystyle= i[n=1Ng(σn+b+σnb),ρSE]\displaystyle-i[\sum_{n=1}^{N}g(\sigma_{n}^{+}b+\sigma_{n}^{-}b^{{\dagger}}),\rho_{S\!E}] (14)
+γ(bρSEb12bbρSE12ρSEbb).\displaystyle+\gamma(b\rho_{S\!E}b^{{\dagger}}-\frac{1}{2}b^{{\dagger}}b\rho_{S\!E}-\frac{1}{2}\rho_{S\!E}b^{{\dagger}}b).

The superradiance phenomenon could be understood as the enhancement of collective spontaneous radiation caused by a common environment (compared with independent environments). To understand the role of memory effects in the superradiance process, we are also interested in the dynamics where the atoms radiate independently, e.g., each atom radiates in its own cavity. This dynamics could be described by the Hamiltonians HI=n=1Ng(σn+bn+σnbn)H_{I}=\sum_{n=1}^{N}g(\sigma_{n}^{+}b_{n}+\sigma_{n}^{-}b_{n}^{{\dagger}}), HE=n=1NωBbnbnH_{E}=\sum_{n=1}^{N}\omega_{B}b_{n}^{{\dagger}}b_{n} and the environment initial state ρE=Πn=1N|0n0|n\rho_{E}=\Pi_{n=1}^{N}|0\rangle_{n}\langle 0|_{n}. Under the resonance condition, the master equation for the independently radiating atoms and the cavities in the interaction picture is given by

ρ˙SE\displaystyle\dot{\rho}_{S\!E} =\displaystyle= i[n=1Ng(σn+bn+σnbn),ρSE]\displaystyle-i[\sum_{n=1}^{N}g(\sigma_{n}^{+}b_{n}+\sigma_{n}^{-}b_{n}^{{\dagger}}),\rho_{S\!E}] (15)
+γn=1N(bnρSEbn12bnbnρSE12ρSEbnbn).\displaystyle+\gamma\sum_{n=1}^{N}(b_{n}\rho_{S\!E}b_{n}^{{\dagger}}\!-\!\frac{1}{2}b_{n}^{{\dagger}}b_{n}\rho_{S\!E}\!-\!\frac{1}{2}\rho_{S\!E}b_{n}^{{\dagger}}b_{n}).

III.2 Early-time solution

In this section, we focus on the early-time dynamics where gt1gt\ll 1 due to the following reasons: First, non-Markovian characters are non-negligible on such a short timescale. Second, it stresses the influence of a particular initial state since the state hardly changes in this time interval. Furthermore, it is helpful for understanding the creation mechanism of the superradiance. In the early-time limit gt0gt\rightarrow 0, one might ignore the influence of cavity dissipation γ\gamma, providing that γ\gamma is not infinitely large. In this case, ρSE\rho_{S\!E} evolves unitarily via

ρSE(t)=U(t)ρS(0)ρEU(t),\displaystyle\rho_{S\!E}(t)=U(t)\rho_{S}(0)\otimes\rho_{E}U^{{\dagger}}(t), (16)

where U(t)=eiHItU(t)=e^{-iH_{I}t} and ρE=|00|\rho_{E}=|0\rangle\langle 0|. The influence of ignoring γ\gamma on the early-time dynamics will be discussed at the end of this section.

The reduced dynamics of the atoms could be determined by tracing out the degrees of the environment:

ρS(t)\displaystyle\rho_{S}(t) =\displaystyle= TrE[U(t)ρS(0)ρEU(t)]\displaystyle\text{Tr}_{E}[U(t)\rho_{S}(0)\otimes\rho_{E}U^{{\dagger}}(t)] (17)
=\displaystyle= kk|EU(t)ρS(0)ρEU(t)|kE,\displaystyle\sum_{k}\langle k|_{E}U(t)\rho_{S}(0)\otimes\rho_{E}U^{{\dagger}}(t)|k\rangle_{E},

where |kE|k\rangle_{E} are a set of basis in E\mathcal{H}_{E}. With the help of the Baker-Campbell-Hausdorff formula

eαABeαA=B+α[A,B]+α22![A,[A,B]]+,\displaystyle e^{\alpha A}Be^{-\alpha A}=B+\alpha[A,B]+\frac{\alpha^{2}}{2!}[A,[A,B]]+\cdots, (18)

after the replacement A=HIA=H_{I}, B=ρS(0)ρEB=\rho_{S}(0)\otimes\rho_{E}, and α=it\alpha=-it, we have

ρS(t)\displaystyle\rho_{S}(t) \displaystyle\approx kk|(ρS(0)ρEit[HI,ρS(0)ρE]\displaystyle\sum_{k}\langle k|(\rho_{S}(0)\otimes\rho_{E}-it[H_{I},\rho_{S}(0)\otimes\rho_{E}] (19)
12t2[HI,[HI,ρS(0)ρE]]|k\displaystyle-\frac{1}{2}t^{2}[H_{I},[H_{I},\rho_{S}(0)\otimes\rho_{E}]]|k\rangle

where the terms with higher orders of gtgt have been omitted due to the early-time limit (gt0gt\rightarrow 0). Using the number state basis of the cavity |k=|0,|1,|k\rangle=|0\rangle,|1\rangle,\cdots and the assumption ρE=|00|\rho_{E}=|0\rangle\langle 0|, it is straightforward to derive the solution of the system evolution in the early-time limit gt0gt\rightarrow 0,

ρS(t)=ρS(0)+(gt)2σ[ρS(0)].\displaystyle\rho_{S}(t)=\rho_{S}(0)+(gt)^{2}\mathcal{L}_{\sigma^{-}}[\rho_{S}(0)]. (20)

Here σ=nσn\sigma^{-}=\sum_{n}\sigma_{n}^{-} is the collective lowering operator and K(ρ)=KρK+12KKρ12ρKK\mathcal{L}_{K}(\rho)=K\rho K^{+}-\frac{1}{2}K^{\dagger}K\rho-\frac{1}{2}\rho K^{\dagger}K is the Lindblad superoperator. b|0=|1b^{{\dagger}}|0\rangle=|1\rangle and 0|b=1|\langle 0|b=\langle 1| have been used to obtain Eq. (20). Notice that the partial trace over the second term of Eq. (19) is zero, thus the change of ρS(0)\rho_{S}(0) is quadratic in tt in the early-time limit. Using similar procedures done above, we obtain the early-time solution for independently radiating atoms:

ρS(t)=ρS(0)+(gt)2n=1Nσn[ρS(0)].\displaystyle\rho_{S}(t)=\rho_{S}(0)+(gt)^{2}\sum_{n=1}^{N}\mathcal{L}_{\sigma_{n}^{-}}[\rho_{S}(0)]. (21)

III.3 Memory effects

With the above results, we now calculate the initial-state-dependent memory effects defined in Sec.II. We consider the dynamics in a short time interval [0,t][0,t] where t0=0t1t2tt_{0}=0\leqslant t_{1}\leqslant t_{2}\leqslant t. In “evolution A” (with a history from t0t_{0} to t1t_{1}) as mentioned before, the initial state is ρS(0)\rho_{S}(0) and the final state ρS(t2)\rho_{S}(t_{2}) is given by

ρS(t2)\displaystyle\rho_{S}(t_{2}) =\displaystyle= ρS(0)+[g(τ10+τ21)]2σ[ρS(0)].\displaystyle\rho_{S}(0)+[g(\tau_{10}+\tau_{21})]^{2}\mathcal{L}_{\sigma^{-}}[\rho_{S}(0)]. (22)

Here τ10=t1t0\tau_{10}=t_{1}-t_{0} and τ21=t2t1\tau_{21}=t_{2}-t_{1} are used for convenience. The initial state of “evolution C” at t1t_{1} is

ρS(t1)\displaystyle\rho_{S}(t_{1}) =\displaystyle= ρS(0)+(gτ10)2σ[ρS(0)].\displaystyle\rho_{S}(0)+(g\tau_{10})^{2}\mathcal{L}_{\sigma^{-}}[\rho_{S}(0)]. (23)

In “evolution C” (without any history before t1t_{1}), the final state ρS(t2)\rho^{\prime}_{S}(t_{2}) is given by

ρS(t2)\displaystyle\rho^{\prime}_{S}(t_{2}) =\displaystyle= ρS(t1)+(gτ21)2σ[ρS(t1)].\displaystyle\rho_{S}(t_{1})+(g\tau_{21})^{2}\mathcal{L}_{\sigma^{-}}[\rho_{S}(t_{1})]. (24)

Substituting Eq. (23) into Eq. (24), we have

ρS(t2)\displaystyle\rho^{\prime}_{S}(t_{2}) =\displaystyle= ρS(0)+(gτ10)2σ[ρS(0)]\displaystyle\rho_{S}(0)+(g\tau_{10})^{2}\mathcal{L}_{\sigma^{-}}[\rho_{S}(0)] (25)
+(gτ21)2σ{ρS(0)+(gτ10)2σ[ρS(0)]}\displaystyle+(g\tau_{21})^{2}\mathcal{L}_{\sigma^{-}}\{\rho_{S}(0)+(g\tau_{10})^{2}\mathcal{L}_{\sigma^{-}}[\rho_{S}(0)]\}
=\displaystyle= ρS(0)+(gτ10)2σ[ρS(0)]\displaystyle\rho_{S}(0)+(g\tau_{10})^{2}\mathcal{L}_{\sigma^{-}}[\rho_{S}(0)]
+(gτ21)2σ[ρS(0)]\displaystyle+(g\tau_{21})^{2}\mathcal{L}_{\sigma^{-}}[\rho_{S}(0)]

where the high-order term g4τ102τ212σ{σ[ρS(0)]}g^{4}\tau_{10}^{2}\tau_{21}^{2}\mathcal{L}_{\sigma^{-}}\{\mathcal{L}_{\sigma^{-}}[\rho_{S}(0)]\} has been omitted due to the early-time limit. The memory effects are manifested by the difference between Eq. (22) and Eq. (25):

ρS(t2)ρS(t2)\displaystyle\rho_{S}(t_{2})\!-\!\rho^{\prime}_{S}(t_{2})\! =\displaystyle= g2[(τ10+τ21)2τ102τ212]σ[ρS(0)]\displaystyle\!g^{2}[(\tau_{10}+\tau_{21})^{2}\!-\!\tau_{10}^{2}\!-\!\tau_{21}^{2}]\mathcal{L}_{\sigma^{-}}[\rho_{S}(0)] (26)
=\displaystyle= 2g2τ10τ21σ[ρS(0)].\displaystyle\!2g^{2}\tau_{10}\tau_{21}\mathcal{L}_{\sigma^{-}}[\rho_{S}(0)].

According to Eq. (11), the value of the memory effects for the initial state ρS(0)\rho_{S}(0) is

NM[ρS(0)]\displaystyle N_{M}[\rho_{S}(0)] =\displaystyle= maxτ10,τ2112ρS(t2)ρS(t2)\displaystyle\max_{\tau_{10},\tau_{21}}\frac{1}{2}\|\rho_{S}(t_{2})-\rho^{\prime}_{S}(t_{2})\| (27)
=\displaystyle= g2σ[ρS(0)]maxτ10,τ21(τ10τ21).\displaystyle g^{2}\|\mathcal{L}_{\sigma^{-}}[\rho_{S}(0)]\|\max_{\tau_{10},\tau_{21}}(\tau_{10}\tau_{21}).

Since g2σ[ρS(0)]g^{2}\|\mathcal{L}_{\sigma^{-}}[\rho_{S}(0)]\| is a constant for a given ρS(0)\rho_{S}(0), Eq. (27) can be calculated simply by maximizing τ10τ21\tau_{10}\tau_{21} with the constraint 0τ10+τ21=τ20t0\leqslant\tau_{10}+\tau_{21}=\tau_{20}\leqslant t. Let τ10+τ21=t\tau_{10}+\tau_{21}=t^{\prime}, then τ10τ21=(τ10t2)2+t24t24\tau_{10}\tau_{21}=-(\tau_{10}-\frac{t^{\prime}}{2})^{2}+\frac{t^{\prime 2}}{4}\leqslant\frac{t^{\prime 2}}{4}. Therefore, the maximum of τ10τ21\tau_{10}\tau_{21} is t24\frac{t^{\prime 2}}{4} when τ10=τ21=t2\tau_{10}=\tau_{21}=\frac{t^{\prime}}{2}. It is easy to see that the maximum of τ10τ21\tau_{10}\tau_{21} in the time interval [0,t][0,t] happens when t=tt^{\prime}=t. Eventually, NM[ρ(t0)]N_{M}[\rho(t_{0})] in the early-time limit is

NM[ρS(0)]\displaystyle N_{M}[\rho_{S}(0)] =\displaystyle= g2σ[ρS(0)]maxτ10,τ21(τ10τ21)\displaystyle g^{2}\|\mathcal{L}_{\sigma^{-}}[\rho_{S}(0)]\|\max_{\tau_{10},\tau_{21}}(\tau_{10}\tau_{21}) (28)
=\displaystyle= g2σ[ρS(0)](t2)2\displaystyle g^{2}\|\mathcal{L}_{\sigma^{-}}[\rho_{S}(0)]\|(\frac{t}{2})^{2}
=\displaystyle= 14(gt)2σ[ρS(0)].\displaystyle\frac{1}{4}(gt)^{2}\|\mathcal{L}_{\sigma^{-}}[\rho_{S}(0)]\|.

Equation (28) demonstrates that in a short time interval [0,t][0,t], NM[ρ(0)]N_{M}[\rho(0)] grows quadratically with tt. To focus on the influence of initial states (rather than the time tt or the atom-cavity coupling gg), it is convenient to discuss the normalized value of memory effects

NM[ρS(0)](gt)2=14σ[ρS(0)].\displaystyle\frac{N_{M}[\rho_{S}(0)]}{(gt)^{2}}=\frac{1}{4}\|\mathcal{L}_{\sigma^{-}}[\rho_{S}(0)]\|. (29)

It represents the strength of memory effects in the early-time limit as a function of only ρS(0)\rho_{S}(0). Notice that NM[ρS(0)](gt)2\frac{N_{M}[\rho_{S}(0)]}{(gt)^{2}} is not bounded as NM[ρS(0)]N_{M}[\rho_{S}(0)]. In principle, there is 0NM[ρS(0)](gt)2<0\leqslant\frac{N_{M}[\rho_{S}(0)]}{(gt)^{2}}<\infty. Similarly, the normalized memory effects for independent cavities (denoted by NMindN_{M}^{ind}) can be derived from Eq. (21). The result can be expressed simply by replacing σ[ρS(0)]\mathcal{L}_{\sigma^{-}}[\rho_{S}(0)] in Eq. (29) by n=1Nσn[ρS(0)]\sum_{n=1}^{N}\mathcal{L}_{\sigma_{n}^{-}}[\rho_{S}(0)], i.e.,

NMind[ρS(0)](gt)2=14n=1Nσn[ρS(0)].\displaystyle\frac{N_{M}^{\textrm{ind}}[\rho_{S}(0)]}{(gt)^{2}}=\frac{1}{4}\|\sum_{n=1}^{N}\mathcal{L}_{\sigma_{n}^{-}}[\rho_{S}(0)]\|. (30)

III.4 Cavity photon number

In a non-Markovian proceess, the relation ρSE(t)ρS(t)ρE\rho_{SE}(t)\approx\rho_{S}(t)\otimes\rho_{E} does not hold in general. Particularly, for the superradiance problem, the change of photon number in the environment (the cavity) might be an important source of the memory effects so that ρSE(t)ρS(t)ρE\rho_{SE}(t)\neq\rho_{S}(t)\otimes\rho_{E}. Thus it is desirable to know the cavity photon number (denoted by NPN_{P} in this paper) in the early-time limit in order to understand the physics of the memory effects. The total excitation number of our model represented by N^ex=bb+nσn+σn\hat{N}_{ex}=b^{\dagger}b+\sum_{n}\sigma_{n}^{+}\sigma_{n}^{-} is conserved since [N^ex,H]=0[\hat{N}_{ex},H]=0. Besides, there are no photon in the cavity initially. Therefore, the cavity photon number is equal to the loss of excitations of the atoms that represents the emission intensity of superradiance (in the early-time limit).

According to Eq. (20), the cavity photon number at tt for ρS(0)\rho_{S}(0) is given by,

NP[ρS(0)]\displaystyle N_{P}[\rho_{S}(0)] =\displaystyle= Tr[nσn+σnρS(0)]Tr[nσn+σnρS(t)]\displaystyle\text{Tr}[\sum_{n}\sigma_{n}^{+}\sigma_{n}^{-}\rho_{S}(0)]-\text{Tr}[\sum_{n}\sigma_{n}^{+}\sigma_{n}^{-}\rho_{S}(t)] (31)
=\displaystyle= Tr[nσn+σnρS(0)]Tr[nσn+σnρS(0)]\displaystyle\text{Tr}[\sum_{n}\sigma_{n}^{+}\sigma_{n}^{-}\rho_{S}(0)]-\text{Tr}[\sum_{n}\sigma_{n}^{+}\sigma_{n}^{-}\rho_{S}(0)]
Tr{nσn+σn(gt)2σ[ρS(0)]}\displaystyle-\text{Tr}\{\sum_{n}\sigma_{n}^{+}\sigma_{n}^{-}(gt)^{2}\mathcal{L}_{\sigma^{-}}[\rho_{S}(0)]\}
=\displaystyle= (gt)2Tr{nσn+σnσ[ρS(0)]}.\displaystyle(gt)^{2}\text{Tr}\{-\sum_{n}\sigma_{n}^{+}\sigma_{n}^{-}\mathcal{L}_{\sigma^{-}}[\rho_{S}(0)]\}.

Using [σ+σ,nσn+σn]=0[\sigma^{+}\sigma^{-},\sum_{n}\sigma_{n}^{+}\sigma_{n}^{-}]=0 and [nσn+σn,σ+]=σ+[\sum_{n}\sigma_{n}^{+}\sigma_{n}^{-},\sigma^{+}]=\sigma^{+}, Eq. (31) can be simplified to a more concise form

NP[ρS(0)]=(gt)2Tr[σ+σρS(0)].\displaystyle N_{P}[\rho_{S}(0)]=(gt)^{2}\text{Tr}[\sigma^{+}\sigma^{-}\rho_{S}(0)]. (32)

It is observed that in the early-time limit, the cavity number increases quadratically with time. As mentioned in the last section, it is convenient to discuss the normalized cavity photon number

NP[ρS(0)](gt)2=Tr[σ+σρS(0)]\displaystyle\frac{N_{P}[\rho_{S}(0)]}{(gt)^{2}}=\text{Tr}[\sigma^{+}\sigma^{-}\rho_{S}(0)] (33)

such that it does not depend on the evolution time tt and atom-cavity coupling gg. Equation (33) represents the emission intensity in the early-time limit as a function of only the initial state. When each atom radiates independently, e.g., each atom radiates in its own cavity, the normalized cavity photon number for the nnth atom is reduced to

NP(n)[ρSn(0)](gt)2=Tr[σn+σnρSn(0)],\displaystyle\frac{N_{P}^{(n)}[\rho_{S_{n}}(0)]}{(gt)^{2}}=\text{Tr}[\sigma_{n}^{+}\sigma_{n}^{-}\rho_{S_{n}}(0)], (34)

where ρSn(0)\rho_{S_{n}}(0) is the reduced density matrix of the nnth atom. Then the degree of superradiant in the early-time regime might be measured by the ratio [9, 15]

S[ρS(0)]\displaystyle S[\rho_{S}(0)] =\displaystyle= NP[ρS(0)]nNP(n)[ρS(0)]\displaystyle\frac{N_{P}[\rho_{S}(0)]}{\sum_{n}N_{P}^{(n)}[\rho_{S}(0)]} (35)
=\displaystyle= Tr[σ+σρS(0)]nTr[σn+σnρSn(0)]\displaystyle\frac{\text{Tr}[\sigma^{+}\sigma^{-}\rho_{S}(0)]}{\sum_{n}\text{Tr}[\sigma_{n}^{+}\sigma_{n}^{-}\rho_{S_{n}}(0)]}

for nonzero denominator. If S[ρS(0)]S[\rho_{S}(0)] is greater (less) than one, the state ρS(0)\rho_{S}(0) is superradiant (subradiant).

III.5 Other manifestations of memory effects

As discussed in Sec. II, the difference between Tr[Pρ(t2)]\text{Tr}[P\rho(t_{2})] and Tr[Pρ(t2)]\text{Tr}[P\rho^{\prime}(t_{2})] can be used as an evidence or a manifestation of the memory effects. For the superradiance problem, the atom excitation number operator P=nσn+σnP=\sum_{n}\sigma_{n}^{+}\sigma_{n}^{-} serves as a nature choice. Following similar treatments in Eqs. (26)-(28), it is seen that the maximum of the atom excitation difference ΔNexatom[ρS(0)]=|Tr[nσn+σnρS(t2)]Tr[nσn+σnρS(t2)]|\Delta N_{ex}^{\textrm{atom}}[\rho_{S}(0)]=|\text{Tr}[\sum_{n}\sigma_{n}^{+}\sigma_{n}^{-}\rho_{S}(t_{2})]-\text{Tr}[\sum_{n}\sigma_{n}^{+}\sigma_{n}^{-}\rho^{\prime}_{S}(t_{2})]| in the time interval [0,t][0,t] also happens at τ10=τ21=t2\tau_{10}=\tau_{21}=\frac{t}{2} and

maxτ10,τ21ΔNexatom[ρS(0)]\displaystyle\max_{\tau_{10},\tau_{21}}\Delta N_{ex}^{\textrm{atom}}[\rho_{S}(0)] (36)
=\displaystyle= |2g2Tr{nσn+σn(t2)2σ[ρS(0)]}|.\displaystyle|2g^{2}\text{Tr}\{\sum_{n}\sigma_{n}^{+}\sigma_{n}^{-}(\frac{t}{2})^{2}\mathcal{L}_{\sigma^{-}}[\rho_{S}(0)]\}|.

Using the results given by Eq. (31) and (32), the manifestation can be further simplified by

maxτ10,τ21ΔNexatom[ρS(0)]\displaystyle\max_{\tau_{10},\tau_{21}}\Delta N_{ex}^{\textrm{atom}}[\rho_{S}(0)] =\displaystyle= (gt)22Tr[σ+σρS(0)]\displaystyle\frac{(gt)^{2}}{2}\text{Tr}[\sigma^{+}\sigma^{-}\rho_{S}(0)] (37)
=\displaystyle= 12NP[ρS(0)].\displaystyle\frac{1}{2}N_{P}[\rho_{S}(0)].

It is seen that the maximal difference of the atom excitation numbers in [0,t][0,t] is exactly half the cavity photon number at tt in the early-time limit. Therefore, the nonzero cavity photon number itself is also a manifestation of memory effects.

III.6 Early-time regime

Mathematically, Eq. (20) holds true as gt0gt\rightarrow 0. One may wonder the validity of the early-time solution Eq. (20) in a longer time interval and how the cavity dissipation deteriorates the validity. In this section, we discuss this problem by comparing the systems dynamics by Eq. (20) with that by numerically solving the full dynamics in SE\mathcal{H}_{S}\otimes\mathcal{H}_{E} with Eq. (14) and tracing out the environment. A dynamics could be represented by the dynamical maps T(t,0)T(t,0) that turns all the possible initial states at 0 to their final states at tt. Furthermore, the difference of two dynamical maps can be evaluated through their Choi-Jamiółkowski matrices [38, 39]. Here we evaluate the error of Eq. (20) at instant tt by the trace distance of ρT(t,0)quad\rho_{T(t,0)}^{\textrm{quad}} and ρT(t,0)exact\rho_{T(t,0)}^{\textrm{exact}}, i.e.,

Error=12ρT(t,0)quadρT(t,0)exact\displaystyle Error=\frac{1}{2}\|\rho_{T(t,0)}^{\textrm{quad}}-\rho_{T(t,0)}^{\textrm{exact}}\| (38)

where ρT(t,0)quad\rho_{T(t,0)}^{\textrm{quad}} and ρT(t,0)exact\rho_{T(t,0)}^{\textrm{exact}} are the Choi-Jamiółkowski matrices calculated by Eq. (20) and (14), respectively. The Choi-Jamiółkowski matrix is defined as ρT(t,0)=T(t,0)𝕀(|ΨΨ|)\rho_{T}(t,0)=T(t,0)\otimes\mathbb{I}(|\Psi\rangle\langle\Psi|) with |Ψ=12Ni=12N|i|i|\Psi\rangle=\frac{1}{\sqrt{2^{N}}}\sum_{i=1}^{2^{N}}|i\rangle|i\rangle a maximally entangled state of the NN-atom system and an ancillary system of the same dimension.

Refer to caption
Figure 2: The error of Eq. (20) as a function of gtgt for (a) 2 atoms and (b) 6 atoms with different strengths of cavity dissipation. The error represents the difference of the dynamical maps T(t,0)T(t,0) corresponding to Eq. (20) (early-time assumption) and Eq. (14) (exact). The error vanishes as gt0gt\rightarrow 0 and increases with the cavity dissipation γ\gamma or the atom number NN.
Refer to caption
Figure 3: The error of Eq. (20) in terms of the atom excitation number difference for (a) 2 and (b) 6 fully excited atoms. Other parameters are the same with those in Fig. 2.

We numerically calculated the error of Eq. (20) for γ/g=0,1,10\gamma/g=0,1,10 and N=2,6N=2,6. The results are illustrated in Fig. 2 in logarithmic scale. It is seen that the order of magnitude of error decreases linearly with that of tt. Although the cavity dissipation and the increasing of atom number can increase the error, there exists a time interval where Eq. (20) is a good approximation as long as γ\gamma and NN are finite. For example, in a time interval [0,0.01/g][0,0.01/g], the error is less than 10510^{-5} even when γ=10g\gamma=10g for both N=2N=2 and N=6N=6, implying that Eq. (20) is a good approximation in this scenario. The above discussion implies that the early-time limit gt0gt\rightarrow 0 could be relaxed to an early-time regime, i.e., a time interval [0,τ][0,\tau] where the decay is quadratic. Here τ\tau is linked to the Zeno time [40, 41, 42, 43]. We also illustrate the difference of the atom excitation numbers from Eqs. (20) and (14), i.e., Errorex=|nσn+σnquadnσn+σnexact|Error_{ex}=|\langle\sum_{n}\sigma_{n}^{+}\sigma_{n}^{-}\rangle_{quad}-\langle\sum_{n}\sigma_{n}^{+}\sigma_{n}^{-}\rangle_{exact}|, in Fig. 3 as another manifestation of the error. The initial states in Figs. 3(a) and 3(b) are 2 and 6 fully excited atoms, respectively. The result is similar to that in Fig. 2 but easier to be verified experimentally.

IV Influence of initial states

In this section, the results Eqs. (29), (30), (33) and (35) are applied to several types of NN-atom initial states. Our first aim is to study the influence of initial states on the memory effects. Another aim is to reveal the role of memory effect in creating superradiance. Conversely, it also helps to understand the physical source of memory effects in our model.

IV.1 Dicke states

The Dicke states, written as |JM|JM\rangle, are extensively studied in the field of superradiance. It is defined as the common eigenstate of the pseudospin operators D2=12(D+D+DD+)+Dz2D^{2}=\frac{1}{2}(D^{+}D^{-}+D^{-}D^{+})+D_{z}^{2} and DzD_{z} with J=N/2J=N/2. The eigenvalues are given by

D2|JM\displaystyle D^{2}|JM\rangle =\displaystyle= J(J+1)|JM\displaystyle J(J+1)|JM\rangle (39)
Dz|JM\displaystyle D_{z}|JM\rangle =\displaystyle= M|JM(M=J,,J)\displaystyle M|JM\rangle\qquad(M=-J,\cdots,J) (40)

where the operators are defined by D+=nσn+D^{+}=\sum_{n}\sigma_{n}^{+}, D=nσnD^{-}=\sum_{n}\sigma_{n}^{-}, and Dz=12n(|ene|n|gng|n)D_{z}=\frac{1}{2}\sum_{n}(|e\rangle_{n}\langle e|_{n}-|g\rangle_{n}\langle g|_{n}). A Dicke state is symmetrical and invariant by atom permutation with Ne=J+MN_{e}=J+M excited atoms and Ng=JMN_{g}=J-M ground-state atoms. It can be constructed by the following formula [10]:

|JM=(J+M)!N!(JM)!(nσn)JM|e,e,,e.\displaystyle|JM\rangle=\sqrt{\frac{(J+M)!}{N!(J-M)!}}(\sum_{n}\sigma_{n}^{-})^{J-M}|e,e,\cdots,e\rangle. (41)

For example, for 3 two-level atoms, there is |3/2,1/2=13(|egg+|geg+|gge)|3/2,-1/2\rangle=\frac{1}{\sqrt{3}}(|egg\rangle+|geg\rangle+|gge\rangle).

A well-known result of the free-space spontaneous emission rate for two-level atoms is WN=Γσ+σW_{N}=\Gamma\langle\sigma^{+}\sigma^{-}\rangle, where Γ\Gamma is the natural linewidth of a single atom. For the Dicke states, there is

WN\displaystyle W_{N} =\displaystyle= ΓJM|σ+σ|JM\displaystyle\Gamma\langle JM|\sigma^{+}\sigma^{-}|JM\rangle (42)
=\displaystyle= Γ(J+M)(JM+1).\displaystyle\Gamma(J+M)(J-M+1).

When M=JM=J (fully excited), WN=NΓW_{N}=N\Gamma. Then the emission rate is proportional to NN, or say, the NN-atom emission rate is equal to the summation of those from independent atoms. Thus there is no superradiance. When M=0M=0 (half excited), WN=ΓN2(N2+1)W_{N}=\Gamma\frac{N}{2}(\frac{N}{2}+1), the emission rate increase with N2N^{2}, or say, the NN-atom emission rate is greater than the summation of those from independent atoms (N2Γ\frac{N}{2}\Gamma). The superradiance with the most intensive emission happens in this case. Note that the results holds under the Born-Markov approximation and on a coarse-grained timescale [10, 11, 12, 13, 14].

Refer to caption
Figure 4: (a) Normalized memory effects NMN_{M} and (b) normalized cavity photon number NPN_{P} in the early-time regime for different Dicke states |JM|JM\rangle. The number of atoms NN varies from 1 to 15, correspondingly, J=N/2=1/2,3/2,,15/2J=N/2=1/2,3/2,\cdots,15/2. The results demonstrate that if the atoms are initially in |JM|JM\rangle, the normalized value of memory effects is half the normalized cavity photon number in the early-time regime. The strongest memory effects as well as photon emission happens when MM is 0 or next to 0.

In the non-Markovian early-time regime, we calculated the normalized value of memory effects and cavity photon number with Eqs. (29) and (33) for the Dicke states of different NN. The results are shown in Fig. 4 where N=1,2,,15N=1,2,\cdots,15 (J=N/2J=N/2) and M=J,J+1,,JM=-J,-J+1,\cdots,J. It is observed that when the number of atoms NN is fixed, the strongest memory effects as well as photon emission happens when MM is 0 or next to 0, which is in agreement with Eq. (42). Interestingly, we find that the value of memory effects for the Dicke states can be fully determined by the cavity photon number in the early-time regime:

NM(|JM)\displaystyle N_{M}(|JM\rangle) =\displaystyle= 12NP(|JM).\displaystyle\frac{1}{2}N_{P}(|JM\rangle). (43)

Considering the result in Eq. (37), there is also

NM(|JM)=maxτ10,τ21ΔNexatom[|JM].\displaystyle N_{M}(|JM\rangle)=\max_{\tau_{10},\tau_{21}}\Delta N_{ex}^{\textrm{atom}}[|JM\rangle]. (44)

Thus the value of memory effects can be directly measured through the atom excitation number difference or the cavity photon number. The relation Eq. (43) is proved in the following. Using the property σ±|JM=J(J+1)M(M±1)|J(M±1)\sigma^{\pm}|JM\rangle=\sqrt{J(J+1)-M(M\pm 1)}|J(M\pm 1)\rangle, it is straightforward to see that

NM(|JM)(gt)2\displaystyle\frac{N_{M}(|JM\rangle)}{(gt)^{2}}\! =\displaystyle= 14σ(|JM)\displaystyle\!\frac{1}{4}\|\mathcal{L}_{\sigma^{-}}(|JM\rangle)\| (45)
=\displaystyle= 14σ|JMJM|σ+12σ+σ|JMJM|\displaystyle\!\frac{1}{4}\|\sigma^{-}|JM\rangle\langle JM|\sigma^{+}\!-\!\frac{1}{2}\sigma^{+}\sigma^{-}|JM\rangle\langle JM|
12|JMJM|σ+σ\displaystyle\!-\frac{1}{2}|JM\rangle\langle JM|\sigma^{+}\sigma^{-}\|
=\displaystyle= 14(J+M)(JM+1)|JMJM|\displaystyle\!\frac{1}{4}(J\!+\!M)(J\!-\!M\!+\!1)\big{\|}|JM^{\prime}\rangle\langle JM^{\prime}|
|JMJM|\displaystyle-|JM\rangle\langle JM|\big{\|}

where M=M1M^{\prime}=M-1. Remind that the trace norm could be calculated by A=TrAA=m|λm|\|A\|=\text{Tr}\sqrt{A^{\dagger}A}=\sum_{m}|\lambda_{m}|, where λm\lambda_{m} is the eigenvalue of operator AA. Meanwhile, the eigenvalues of |JMJM||JMJM||JM^{\prime}\rangle\langle JM^{\prime}|-|JM\rangle\langle JM| is 11 and 1-1 whose absolute values sum up to 2. Therefore,

NM(|JM)(gt)2=12(J+M)(JM+1).\displaystyle\frac{N_{M}(|JM\rangle)}{(gt)^{2}}=\frac{1}{2}(J+M)(J-M+1). (46)

Considering that the normalized cavity photon number for the Dicke states is

NP(|JM)(gt)2\displaystyle\frac{N_{P}(|JM\rangle)}{(gt)^{2}} =\displaystyle= Tr(σ+σ|JMJM|)\displaystyle\text{Tr}(\sigma^{+}\sigma^{-}|JM\rangle\langle JM|) (47)
=\displaystyle= Tr[(J+M)(JM+1)|JMJM|]\displaystyle\text{Tr}[(J+M)(J-M+1)|JM\rangle\langle JM|]
=\displaystyle= (J+M)(JM+1),\displaystyle(J+M)(J-M+1),

Eq. (43) is proved. The relation Eq. (43) demonstrates that the (change of) the cavity photon number is a fundamental reason of the memory effects in the early-time regime for the Dicke states.

Refer to caption
Figure 5: (a) The degree of superradiance SS and (b) the degree of memory-effect-enhancement NM/NMindN_{M}/N_{M}^{\textrm{ind}} (from independent radiation to collective radiation) in the early-time regime for different Dicke states. The states are the same with those in Fig. 4 except that SS and NM/NMindN_{M}/N_{M}^{\textrm{ind}} for M=JM=-J are undefined and not shown. For a given JJ, the highest degree of superradiance happens when M=J+1M=-J+1 (single excitation).

We now discuss the degree of superadiance for the Dicke states in the early-time regime. The normalized cavity photon number for one independent atom in its excited state is given by NP(|e)/(gt)2=NP(|12,12)/(gt)2=1N_{P}(|e\rangle)/(gt)^{2}=N_{P}(|\frac{1}{2},\frac{1}{2}\rangle)/(gt)^{2}=1 by Eq. (34). Using Eq. (47), the degree of superradiance for the Dicke states could be calculated by

S(|JM)\displaystyle S(|JM\rangle) =\displaystyle= NP(|JM)n=1J+MNP(|e)\displaystyle\frac{N_{P}(|JM\rangle)}{\sum_{n=1}^{J+M}N_{P}(|e\rangle)} (48)
=\displaystyle= (J+M)(JM+1)J+M\displaystyle\frac{(J+M)(J-M+1)}{J+M}
=\displaystyle= JM+1\displaystyle J-M+1

for J+M>0J+M>0 as shown in Fig. 5(a). Here the denominator represents the contribution of the emission from the J+MJ+M independently excited atoms in |JM|JM\rangle. Note that the physical meaning of the denominator in Eq. (48) is different from that in Eq. (35), but their values are the same. In Fig. 5(a), SS for |J,J|J,-J\rangle (no atom excited) are undefined and not shown. Another group of Dicke states without superradiance are |J,J|J,J\rangle (fully excited) corresponding to the right row in Fig. 5(a). Except for the two groups, other Dicke states are superradiant which could be examined by S(|JM)>1S(|JM\rangle)>1 with Eq. (48). For example, for |J,J+1|J,-J+1\rangle (one atom excited) represented by the left row in Fig. 5(a), there is S(|J,J+1)=NS(|J,-J+1\rangle)=N. The single-photon superradiance is called “the greatest radiation anomaly” by Dicke. For the states in the right row in Fig. 5(a), there is S(|J,J1)=2S(|J,J-1\rangle)=2 for N>1N>1, showing a weaker degree of superradiance.

The degree of superradiance represents the magnification of the radiation intensity caused by a common environment (compared with NN independent environments). One interesting question is whether the memory effects are simultaneously enhanced by a common environment. To uncover the role of memory effects in the superradiance process, we next examine the degree of memory-effect-enhancement, i.e., NM(|JM)/NMind(|JM)N_{M}(|JM\rangle)/N_{M}^{\textrm{ind}}(|JM\rangle), where NMind(|JM)N_{M}^{\textrm{ind}}(|JM\rangle) is the value of memory effects with NN independent cavities calculated by Eq. (30). We find by mathematical induction that the value of memory effects for independently radiating atoms in a Dicke state is NMind(|JM)=12(gt)2(J+M)N_{M}^{\textrm{ind}}(|JM\rangle)=\frac{1}{2}(gt)^{2}(J+M). Therefore, the degree of memory-effect-enhancement in the early-time regime is

NM(|JM)NMind(|JM)=JM+1=S(|JM)\displaystyle\frac{N_{M}(|JM\rangle)}{N_{M}^{\textrm{ind}}(|JM\rangle)}=J-M+1=S(|JM\rangle) (49)

as shown in Fig. 5(b). The enhancement of radiation intensity (caused by a common environment) is equal to the enhancement of memory effects for the Dicke states, implying that the memory effects of a common environment contribute to the establishment of superradiance.

IV.2 Factorized identical states

In this section, we investigate NN-atom initial states in a factorized from

ρSfact(0)\displaystyle\rho_{S}^{\textrm{fact}}(0) =\displaystyle= ρS1ρS2ρSN\displaystyle\rho_{S_{1}}\otimes\rho_{S_{2}}\otimes\cdots\otimes\rho_{S_{N}} (50)

with identical single-atom states

ρS1=ρS2==ρSN=(ρeeρegρegρgg).\displaystyle\rho_{S_{1}}=\rho_{S_{2}}=\cdots=\rho_{S_{N}}=\left(\begin{array}[]{cc}\rho_{ee}&\rho_{eg}\\ \rho_{eg}^{*}&\rho_{gg}\\ \end{array}\right). (53)

Such states are fully characterized by ρee\rho_{ee} and ρeg\rho_{eg} and may also be superradiant. For example, states of this form are experimentally studied in Ref. [15] to realize the single-atom superradiance where a serious of atoms enter a cavity one by one. In Ref. [23], we show that such a model can be used to achieve the Heisenberg Limit in a quantum metrology scheme for measuring the atom-cavity coupling. The superradiance of two artificial atoms with initial state 12(|e+|g)(|e+|g)\frac{1}{2}(|e\rangle+|g\rangle)\otimes(|e\rangle+|g\rangle) are experimentally observed in Ref. [17]. In the above works, the single-atom coherence (|ρeg||\rho_{eg}|) are crucial for creating superradiance. Thus it is desirable to know the influence of such states on the memory effects in superradiance, especially the role of the single-atom coherence.

Refer to caption
Figure 6: Normalized memory effects NMN_{M} (left column) and normalized cavity photon number NPN_{P} (right column) for ρSfact(0)\rho_{S}^{\textrm{fact}}(0) as functions of ρee\rho_{ee} and |ρeg||\rho_{eg}| in the early-time regime. The atom numbers are N=1N=1, 22, and 1111 from top to bottom. For any ρee\rho_{ee} and NN (N>1N>1), NMN_{M} (NPN_{P}) can always be enhanced nonlinearly by the single-atom coherence |ρeg||\rho_{eg}|.

In the following, the normalized value of memory effects and the normalized cavity photon number for states described by Eq. (50) are calculated using Eqs. (29) and (33). When Eq. (53) is a pure state, i.e., ρSn=|ϕϕ|\rho_{Sn}=|\phi\rangle\langle\phi|, ρSfact(0)\rho_{S}^{\textrm{fact}}(0) can be decomposed into the superpositions of different Dicke states [15]. Using the analytical result of ϕ|σ+σ|ϕ\langle\phi|\sigma^{+}\sigma^{-}|\phi\rangle in Ref. [15], we obtain the normalized cavity photon number for ρSfact(0)\rho_{S}^{\textrm{fact}}(0):

NP[ρSfact(0)](gt)2=N(N1)|ρeg|2+Nρee.\displaystyle\frac{N_{P}[\rho_{S}^{\textrm{fact}}(0)]}{(gt)^{2}}=N(N-1)|\rho_{eg}|^{2}+N\rho_{ee}. (54)

In contrast with Eq. (46), the analytical solution of Eq. (29) for the factorized initial states is more complicated in general. Therefore, Eq. (29) is solved numerically here. The results of the normalized NM[ρSfact(0)]N_{M}[\rho_{S}^{\textrm{fact}}(0)] and NP[ρSfact(0)]N_{P}[\rho_{S}^{\textrm{fact}}(0)] are depicted in Fig. 6 as functions of ρee\rho_{ee} and |ρeg||\rho_{eg}| with N=1N=1, 2 and 1111. The condition |ρeg|2(1ρee)ρee|\rho_{eg}|^{2}\leqslant(1-\rho_{ee})\rho_{ee} are satisfied in Fig. 6 to guarantee the semipositivity of ρSn\rho_{S_{n}}. The single-atom coherence ρeg\rho_{eg} are chosen to be real in the simulations for simplicity. Numerical results shows that the phase of ρeg\rho_{eg} does not affect the value of memory effects. Unlike the case for the Dicke states, NM[ρSfact(0)]/(gt)2N_{M}[\rho_{S}^{\textrm{fact}}(0)]/(gt)^{2} is not proportional to NP[ρSfact(0)]/(gt)2N_{P}[\rho_{S}^{\textrm{fact}}(0)]/(gt)^{2} in general. It is observed that the normalized value of memory effects can always be enhanced by |ρeg||\rho_{eg}| for any NN. Similar situations happen for NP[ρSfact(0)]/(gt)2N_{P}[\rho_{S}^{\textrm{fact}}(0)]/(gt)^{2} when N>1N>1 as seen in Eq. (54) and Fig. 6 (right column). Besides, as NN increases, NM[ρSfact(0)]/(gt)2N_{M}[\rho_{S}^{\textrm{fact}}(0)]/(gt)^{2} grows more and more fast with |ρeg||\rho_{eg}|, so does NP[ρSfact(0)]/(gt)2N_{P}[\rho_{S}^{\textrm{fact}}(0)]/(gt)^{2}. Thus we infer that for states with large NN, the cavity photon number NPN_{P} are dominated by |ρeg||\rho_{eg}|, and the memory effects are dominated by NPN_{P}. In Fig. 6(c) and (e), it is found that the states around ρeg=0\rho_{eg}=0 and ρee=0.5\rho_{ee}=0.5 tend to take lower values of memory effects compared with their neighboring states, especially when NN is large. This might be intuitively interpreted by the fact that ρeg=0\rho_{eg}=0 and ρee=0.5\rho_{ee}=0.5 leads to a maximally mixed state of the NN-atom system. Besides, Fig. 6 shows that the (change of) the cavity photon number is not the only factor to determine the value of the memory effects. Other factors might include the (change of) coherence of the cavity state between the basis |n|n\rangle.

Refer to caption
Figure 7: The degree of superradiance SS (left column) and the degree of memory-effect-enhancement NM/NMindN_{M}/N_{M}^{\textrm{ind}} (from independent radiation to collective radiation, right column) for ρSfact(0)\rho_{S}^{\textrm{fact}}(0) as functions of ρee\rho_{ee} and |ρeg||\rho_{eg}| in the early-time regime. The values of SS and NM/NMindN_{M}/N_{M}^{\textrm{ind}} for ρee=0\rho_{ee}=0 are not defined and not shown. The atom numbers are N=2N=2, 33, and 1212 from top to bottom. It is observed that large values of SS and NM/NMindN_{M}/N_{M}^{\textrm{ind}} happen for low-excitation states with high coherence.

According to Eqs. (35) and (54), the degree of superradiance for ρSfact(0)\rho_{S}^{\textrm{fact}}(0) is given by

S[ρSfact(0)]\displaystyle S[\rho_{S}^{\textrm{fact}}(0)] =\displaystyle= NP[ρSfact(0)]n=1NNP(n)(ρSn)\displaystyle\frac{N_{P}[\rho_{S}^{\textrm{fact}}(0)]}{\sum_{n=1}^{N}N_{P}^{(n)}(\rho_{S_{n}})} (55)
=\displaystyle= N(N1)|ρeg|2+NρeeNρee\displaystyle\frac{N(N-1)|\rho_{eg}|^{2}+N\rho_{ee}}{N\rho_{ee}}
=\displaystyle= 1+(N1)|ρeg|2ρee\displaystyle 1+(N-1)\frac{|\rho_{eg}|^{2}}{\rho_{ee}}

for state with ρee>0\rho_{ee}>0. The state ρSfact(0)\rho_{S}^{\textrm{fact}}(0) is superradiant when |ρeg||\rho_{eg}| is nonzero and N>1N>1. Besides, the degree of superradiance grow quadratically with the single-atom coherence |ρeg||\rho_{eg}| for N>1N>1. The maximum of S[ρSfact(0)]S[\rho_{S}^{\textrm{fact}}(0)] is NN when ρee0\rho_{ee}\rightarrow 0 and |ρeg|=(1ρee)ρee|\rho_{eg}|=\sqrt{(1-\rho_{ee})\rho_{ee}} (its maximum under the semi-positivity constraint). We illustrate the degrees of superradiance for different atom numbers in Fig. 7 (left column). As done in last section, we also calculate the degree of memory-effect-enhancement NM[ρSfact(0)]/NMind[ρSfact(0)]N_{M}[\rho_{S}^{\textrm{fact}}(0)]/N_{M}^{\textrm{ind}}[\rho_{S}^{\textrm{fact}}(0)] in the early-time regime with Eqs. (29) and (30) for different atom numbers. These results are shown in Fig. 7 (right column). Unlike the case for the Dicke states, the dependence of NM/NMindN_{M}/N_{M}^{\textrm{ind}} on ρee\rho_{ee} and ρeg\rho_{eg} is not the same as that of SS. But there are many similarities between them. When there are only one atom (N=1N=1), NM[ρSfact(0)]/NMind[ρSfact(0)]=S[ρSfact(0)]=1N_{M}[\rho_{S}^{\textrm{fact}}(0)]/N_{M}^{\textrm{ind}}[\rho_{S}^{\textrm{fact}}(0)]=S[\rho_{S}^{\textrm{fact}}(0)]=1 according to our definitions. For N=2N=2, SS and NM[ρSfact(0)]/NMind[ρSfact(0)]N_{M}[\rho_{S}^{\textrm{fact}}(0)]/N_{M}^{\textrm{ind}}[\rho_{S}^{\textrm{fact}}(0)] are slightly different, implying that the memory effects are not merely determined by the cavity photon numbers for the factorized identical states. As the atom number increases (for N3N\geqslant 3), the surface shapes in the right-side figures become smoother and more similar to those in the left column, especially when SS is large. As shown in Fig. 7, strong memory-effect-enhancement happens for low-excitation states with high coherence (as SS does). Similar to Eq. (55), numerical simulations demonstrate that the maximum of NM[ρSfact(0)]/NMind[ρSfact(0)]N_{M}[\rho_{S}^{\textrm{fact}}(0)]/N_{M}^{\textrm{ind}}[\rho_{S}^{\textrm{fact}}(0)] is also NN when ρee0\rho_{ee}\rightarrow 0 and |ρeg|=(1ρee)ρee|\rho_{eg}|=\sqrt{(1-\rho_{ee})\rho_{ee}}. These results imply that for states with high degree of superradiance, the (enhancement of) memory effects are dominated by the cavity photon number. Besides, numerical simulations show that for initial states whose spontaneous radiation can be enhanced by a common environment (S>1S>1), their memory effects can always be enhanced simultaneously (NM[ρSfact(0)]/NMind[ρSfact(0)]>1N_{M}[\rho_{S}^{\textrm{fact}}(0)]/N_{M}^{\textrm{ind}}[\rho_{S}^{\textrm{fact}}(0)]>1). Therefore, we infer that the memory effects also contribute to the establishment of the superradiance for the factorized identical initial states.

IV.3 Other initial states

Entanglement and coherence of the initial states may play important roles in superradiance [44, 24, 25, 26, 27], especially for the Dicke states. Next, we examine the influence of the coherence (or the entanglement among the atoms) in the dephased Dicke states ρSdeph(0)\rho^{\textrm{deph}}_{S}(0) on the memory effects as well as the superradiance. The initial state we consider is given by

ρSdeph(0)=λ|JMJM|+(1λ)𝒟(|JMJM|)\displaystyle\rho^{\textrm{deph}}_{S}(0)=\lambda|JM\rangle\langle JM|+(1-\lambda)\mathcal{D}(|JM\rangle\langle JM|) (56)

where 0λ10\leqslant\lambda\leqslant 1. Here 𝒟\mathcal{D} turns |JMJM||JM\rangle\langle JM| into a fully dephased state with vanishing off-diagonal elements in the basis |e(g),e(g),,e(g)|e(g),e(g),\cdots,e(g)\rangle. Meanwhile, 𝒟(|JMJM|)\mathcal{D}(|JM\rangle\langle JM|) is a mixed state of separable states whose entanglement is zero. Thus λ\lambda reflects the strength of entanglement or the value of coherence in ρSdeph(0)\rho^{\textrm{deph}}_{S}(0). With the help of numerical solutions, we find by mathematical induction that NM[𝒟(|JMJM|)](gt)2=12(J+M)\frac{N_{M}[\mathcal{D}(|JM\rangle\langle JM|)]}{(gt)^{2}}=\frac{1}{2}(J+M) and

NM[ρSdeph(0)](gt)2\displaystyle\frac{N_{M}[\rho^{\textrm{deph}}_{S}(0)]}{(gt)^{2}} =\displaystyle= λNM(|JMJM|)(gt)2\displaystyle\lambda\frac{N_{M}(|JM\rangle\langle JM|)}{(gt)^{2}} (57)
+(1λ)NM[𝒟(|JMJM|)](gt)2\displaystyle+(1-\lambda)\frac{N_{M}[\mathcal{D}(|JM\rangle\langle JM|)]}{(gt)^{2}}
=\displaystyle= 12(J+M)[(JM)λ+1].\displaystyle\frac{1}{2}(J+M)[(J-M)\lambda+1].

It is seen that the normalized value of memory effects increases linearly with λ\lambda when MJ,JM\neq J,-J. The normalized cavity photon number can be calculated by Eq. (33). We obtain

NP[ρSdeph(0)](gt)2=\displaystyle\frac{N_{P}[\rho^{\textrm{deph}}_{S}(0)]}{(gt)^{2}}= =\displaystyle= λNP(|JM)(gt)2\displaystyle\lambda\frac{N_{P}(|JM\rangle)}{(gt)^{2}} (58)
+(1λ)NP[𝒟(|JMJM|)](gt)2\displaystyle+(1-\lambda)\frac{N_{P}[\mathcal{D}(|JM\rangle\langle JM|)]}{(gt)^{2}}
=\displaystyle= (J+M)[(JM)λ+1].\displaystyle(J+M)[(J-M)\lambda+1].

Similarly, NM[ρSdeph(0)]=12NP[ρSdeph(0)]=maxτ10,τ21ΔNexatom[ρSdeph(0)]N_{M}[\rho^{\textrm{deph}}_{S}(0)]=\frac{1}{2}N_{P}[\rho^{\textrm{deph}}_{S}(0)]=\max_{\tau_{10},\tau_{21}}\Delta N_{ex}^{\textrm{atom}}[\rho^{\textrm{deph}}_{S}(0)] is satisfied for the dephased Dicke states in the early-time regime.

The degree of superradiance for the dephased Dicke states can be calculated from Eq. (58) by

S[ρSdeph(0)]\displaystyle S[\rho^{\textrm{deph}}_{S}(0)] =\displaystyle= NP[ρSdeph(0)]n=1J+MNP(|e)\displaystyle\frac{N_{P}[\rho^{\textrm{deph}}_{S}(0)]}{\sum_{n=1}^{J+M}N_{P}(|e\rangle)} (59)
=\displaystyle= (J+M)[(JM)λ+1]J+M\displaystyle\frac{(J+M)[(J-M)\lambda+1]}{J+M}
=\displaystyle= (JM)λ+1\displaystyle(J-M)\lambda+1

for J+M>0J+M>0. It is observed that the entanglement or coherence (represented by λ\lambda) in the depahsed Dicke state is necessary for superradiance. The degree of superradiance increases linearly with λ\lambda for superradiant Dicke states. Next we discuss the degree of memory-effect-enhancement caused by a common environment. When the initial state is 𝒟(|JMJM|)\mathcal{D}(|JM\rangle\langle JM|) (fully dephased), we find that NMind[𝒟(|JMJM|]/(gt)2=NMind[|JMJM|]/(gt)2=12(J+M)N_{M}^{\textrm{ind}}[\mathcal{D}(|JM\rangle\langle JM|]/(gt)^{2}=N_{M}^{\textrm{ind}}[|JM\rangle\langle JM|]/(gt)^{2}=\frac{1}{2}(J+M). Moreover, for a partially dephased Dicke state described by Eq. (56), there is NMind[ρSdeph(0)]/(gt)2=12(J+M)N_{M}^{\textrm{ind}}[\rho^{\textrm{deph}}_{S}(0)]/(gt)^{2}=\frac{1}{2}(J+M) regardless of λ\lambda. Therefore, the degree of memory-effect-enhancement for ρSdeph(0)\rho^{\textrm{deph}}_{S}(0) is given by

NM[ρSdeph(0)]NMind[ρSdeph(0)]=(JM)λ+1=S[ρSdeph(0)]\displaystyle\frac{N_{M}[\rho^{\textrm{deph}}_{S}(0)]}{N_{M}^{\textrm{ind}}[\rho^{\textrm{deph}}_{S}(0)]}=(J-M)\lambda+1=S[\rho^{\textrm{deph}}_{S}(0)] (60)

as the case in Sec.IV A.

At the end of this section, we calculate the value of the normalized value of memory effects of a probabilistic mixture of Dicke states ρSmix(0)=MpM|JMJM|\rho_{S}^{\textrm{mix}}(0)=\sum_{M}p_{M}|JM\rangle\langle JM| where J=N/2J=N/2, M=J,J+1,,JM=-J,-J+1,\cdots,J, and MpM=1\sum_{M}p_{M}=1. Using similar treatments as done in Eq. (45), we find that

NM[ρSmix(0)](gt)2\displaystyle\frac{N_{M}[\rho_{S}^{\textrm{mix}}(0)]}{(gt)^{2}}\! =\displaystyle= 14[M=JJ1|pM+1f(M+1)pMf(M)|\displaystyle\!\frac{1}{4}[\!\sum_{M=-J}^{J-1}\!|p_{M+1}f(M+1)-p_{M}f(M)| (61)
+pJf(J)]\displaystyle+p_{J}f(J)]

where f(X)=(J+X)(JX+1)f(X)=(J+X)(J-X+1). The result reduces to Eq. (46) when ρSmix(0)\rho_{S}^{\textrm{mix}}(0) is a pure state |JMJM||JM\rangle\langle JM|.

V Results in other regimes

V.1 Near-Markovian regime

In the previous sections, we investigate the superradiance and memory effect characteristics of our model in the early-time regime with quadratic expressions. The early-time dynamics can be understood by setting γ=0\gamma=0 in Eq. (14) and consider the dynamics in a short time interval (gt1gt\ll 1). The definition of memory effects in Ref. [35] reflects the quality of Markovian approximation in a quantum process. Correspondingly, the initial-state-dependent memory effect in this paper reflects the quality of Markovian approximation conditioned on a particular initial state, which could be understood by 1NM[ρS(0)]1-N_{M}[\rho_{S}(0)]. By our definition, NM[ρS(0)]N_{M}[\rho_{S}(0)] is generally nonzero (although small) for nontrivial exact dynamics. Therefore, it could be used to investigate another interesting regime, i.e., the near-Markovian regime, by setting γ/g1\gamma/g\gg 1 and considering the exact dynamics in 0<t<0<t<\infty or as least a longer time interval far beyond the environment memory time τE=1/γ\tau_{E}=1/\gamma. In this section, for similar reasons, we consider the early-stage (far beyond τE\tau_{E}) dynamics in the near-Markovian regime where the system’s state does not change so much. In this way, we stress the influence of the initial states on the superradiance and the memory effects, and compare the results with those in the early-time regime. The specific scale of the early-stage time interval will be discussed in the following.

Refer to caption
Figure 8: An example of the system and environment dynamics as well as the memory effect manifestations in a near-Markovian regime with γ/g=103\gamma/g=10^{3} and a common cavity. The initial state is |JM=|2,0|JM\rangle=|2,0\rangle. (a) and (b) show the evolutions of the excitation number of the atoms NexatomN_{ex}^{\textrm{atom}} and the cavity photon number NPN_{P}, respectively, where the insets display their early-stage dynamics. (c) and (d) show the difference of the excitation numbers ΔNexatom\Delta N_{ex}^{\textrm{atom}} and the trace distance DD, respectively, as functions of τ10\tau_{10} and τ21\tau_{21} (both in units of 1/γ1/\gamma).

The near-Markovian dynamics of the atoms can be calculated by numerically solving Eqs. (14) or (15) and tracing out the environment under the condition γ/g1\gamma/g\gg 1. We first present a near-Markovian (γ/g=103\gamma/g=10^{3}) example showing typical system and environment dynamics and the memory effect manifestations in Fig. 8, where 4 atoms initially in the Dicke state |2,0|2,0\rangle radiate in a common cavity. The system dynamics and environment dynamics are represented by the evolutions of the atom excitation number Nexatom=nσn+σnN_{ex}^{\textrm{atom}}=\langle\sum_{n}\sigma_{n}^{+}\sigma_{n}^{-}\rangle and the cavity photon number NP=bbN_{P}=\langle b^{\dagger}b\rangle, respectively, as shown in Figs. 8(a) and 8((b). The memory effects are based on the trace distance of the final states in evolution A and C, i.e., D=12ρS(t2)ρS(t2)D=\frac{1}{2}\|\rho_{S}(t_{2})-\rho^{\prime}_{S}(t_{2})\| which is shown in Fig. 8(d) as a function of τ10=t1t0\tau_{10}=t_{1}-t_{0} and τ21=t2t1\tau_{21}=t_{2}-t_{1}. Similarly, we also illustrate the difference of the atom excitation numbers denoted by ΔNexatom=|Tr[nσn+σnρS(t2)]Tr[nσn+σnρS(t2)]|\Delta N_{ex}^{\textrm{atom}}=|\text{Tr}[\sum_{n}\sigma_{n}^{+}\sigma_{n}^{-}\rho_{S}(t_{2})]-\text{Tr}[\sum_{n}\sigma_{n}^{+}\sigma_{n}^{-}\rho^{\prime}_{S}(t_{2})]| as a manifestation of memory effects in Fig. 8(c).

Table 1: Characteristics of superradiance and memory effects for the dephased Dicke states (λ=1\lambda=1 for the Dicke states) in different regimes (τE=1/γ\tau_{E}=1/\gamma).
Characteristics  Early-time result  Characteristics  Early-stage Markovian result
NP/(gt)2N_{P}/(gt)^{2} (J+M)[(JM)λ+1](J+M)[(J-M)\lambda+1] NPsteady/(gτE)2N_{P}^{\textrm{steady}}/(g\tau_{E})^{2} 4(J+M)[(JM)λ+1]4(J+M)[(J-M)\lambda+1]
NM/(gt)2N_{M}/(gt)^{2} 12(J+M)[(JM)λ+1]\frac{1}{2}(J+M)[(J-M)\lambda+1] NM/(gτE)2N_{M}/(g\tau_{E})^{2} 8(J+M)[(JM)λ+1]8(J+M)[(J-M)\lambda+1]
NPind/(gt)2N_{P}^{\textrm{ind}}/(gt)^{2} J+MJ+M NPsteady,ind/(gτE)2N_{P}^{\textrm{steady},\textrm{ind}}/(g\tau_{E})^{2} 4(J+M)4(J+M)
NMind/(gt)2N_{M}^{\textrm{ind}}/(gt)^{2} 12(J+M)\frac{1}{2}(J+M) NMind/(gτE)2N_{M}^{\textrm{ind}}/(g\tau_{E})^{2} 8(J+M)8(J+M)
S=NP/NPindS=N_{P}/N_{P}^{\textrm{ind}} (JM)λ+1(J-M)\lambda+1 S=NPsteady/NPstead,indS=N_{P}^{\textrm{steady}}/N_{P}^{\textrm{stead},\textrm{ind}} (JM)λ+1(J-M)\lambda+1
NM/NMindN_{M}/N_{M}^{\textrm{ind}} (JM)λ+1(J-M)\lambda+1 NM/NMindN_{M}/N_{M}^{\textrm{ind}} (JM)λ+1(J-M)\lambda+1
Table 2: Characteristics of superradiance and memory effects for the factorized identical states in different regimes (τE=1/γ\tau_{E}=1/\gamma).
Characteristics  Early-time result  Characteristics  Early-stage Markovian result
NP/(gt)2N_{P}/(gt)^{2} N(N1)|ρeg|2+NρeeN(N-1)|\rho_{eg}|^{2}+N\rho_{ee} NPsteady/(gτE)2N_{P}^{\textrm{steady}}/(g\tau_{E})^{2} 4[N(N1)|ρeg|2+Nρee]4[N(N-1)|\rho_{eg}|^{2}+N\rho_{ee}]
NM/(gt)2N_{M}/(gt)^{2} 14σ[ρS(0)]\frac{1}{4}\|\mathcal{L}_{\sigma^{-}}[\rho_{S}(0)]\| NM/(gτE)2N_{M}/(g\tau_{E})^{2} 4σ[ρS(0)]4\|\mathcal{L}_{\sigma^{-}}[\rho_{S}(0)]\|
NPind/(gt)2N_{P}^{\textrm{ind}}/(gt)^{2} NρeeN\rho_{ee} NPsteady,ind/(gτE)2N_{P}^{\textrm{steady},\textrm{ind}}/(g\tau_{E})^{2} 4Nρee4N\rho_{ee}
NMind/(gt)2N_{M}^{\textrm{ind}}/(gt)^{2} 14n=1Nσn[ρS(0)]\frac{1}{4}\|\sum_{n=1}^{N}\mathcal{L}_{\sigma_{n}^{-}}[\rho_{S}(0)]\| NMind/(gτE)2N_{M}^{\textrm{ind}}/(g\tau_{E})^{2} 4n=1Nσn[ρS(0)]4\|\sum_{n=1}^{N}\mathcal{L}_{\sigma_{n}^{-}}[\rho_{S}(0)]\|
S=NP/NPindS=N_{P}/N_{P}^{\textrm{ind}} 1+(N1)|ρeg|2/ρee1+(N-1)|\rho_{eg}|^{2}/\rho_{ee} S=NPsteady/NPstead,indS=N_{P}^{\textrm{steady}}/N_{P}^{\textrm{stead},\textrm{ind}} 1+(N1)|ρeg|2/ρee1+(N-1)|\rho_{eg}|^{2}/\rho_{ee}
NM/NMindN_{M}/N_{M}^{\textrm{ind}} σ[ρS(0)]/n=1Nσn[ρS(0)]\|\mathcal{L}_{\sigma^{-}}[\rho_{S}(0)]\|/\|\sum_{n=1}^{N}\mathcal{L}_{\sigma_{n}^{-}}[\rho_{S}(0)]\| NM/NMindN_{M}/N_{M}^{\textrm{ind}} σ[ρS(0)]/n=1Nσn[ρS(0)]\|\mathcal{L}_{\sigma^{-}}[\rho_{S}(0)]\|/\|\sum_{n=1}^{N}\mathcal{L}_{\sigma_{n}^{-}}[\rho_{S}(0)]\|

As in Fig. 8, further numerical results shows that in a near-Markovian regime (γ/g1\gamma/g\gg 1) with initially excited Dicke states or factorized identical states, the cavity photon number NPN_{P} first experiences a rapid increase before t10/γ=10τEt\sim 10/\gamma=10\tau_{E} and then reaches a plateau while the atom emission rate is almost constant as shown in the insets of Figs. 8(b) and 8(a). After the plateau, NPN_{P} evolves on a larger timescale given by gtγ/ggt\propto\gamma/g where NPN_{P} may monotonically decay for a highly superradiant initial state, or first increase and then decay, for a less superradiant one. Thus the early stage of a near-Markovian dynamics in our work could be understood as a time interval [0,t~][0,\tilde{t}] where τEt~τL\tau_{E}\ll\tilde{t}\ll\tau_{L}. Here τL\tau_{L} is the lifetime of the radiation that is found to be proportional to γ/g\gamma/g and initial-state-dependent. In other words, in the early stage of a near-Markovian dynamics, the environment photon number and the atom emission rate becomes almost steady while the system does not change so much. Therefore, the almost-steady cavity photon number NPsteadyN_{P}^{\textrm{steady}}, as shown by the plateau in the inset of Fig. 8(b), represents the typical cavity photon number during the early stage of a near-Markovian process. In the near-Markovian regime, the memory effects NM[ρS(0)]N_{M}[\rho_{S}(0)] can be calculated by Eq. (11) where 0t1t2t~0\leqslant t_{1}\leqslant t_{2}\leqslant\tilde{t}. Numerical evidences imply that there exists only one local maximum of D=12ρS(t2)ρS(t2)D=\frac{1}{2}\|\rho_{S}(t_{2})-\rho^{\prime}_{S}(t_{2})\| when τ10\tau_{10} and τ21\tau_{21} are optimized (at least) in the early stage of a near-Markovian dynamics. For example, in Fig. 8, the optimal τ10\tau_{10} and τ21\tau_{21} are about τ10=τ2120/γ\tau_{10}=\tau_{21}\sim 20/\gamma in 0t1t20\leqslant t_{1}\leqslant t_{2}\leqslant\infty. In contrast, the optimal τ10\tau_{10} and τ21\tau_{21} that maximizing ΔNexatom\Delta N_{ex}^{\textrm{atom}} might be beyond the early-stage time interval for a less superadiant initial state. However, in practice, the values of DD as well as ΔNexatom\Delta N_{ex}^{\textrm{atom}} in the early stage of a near-Markovian dynamics are nearly steady and robust against τ10\tau_{10} and τ21\tau_{21} after τ10=τ21>10τE\tau_{10}=\tau_{21}>10\tau_{E}, especially when γ/g\gamma/g is large. This makes it easy to define and numerically calculate the near-steady value of the excitation number difference ΔNex,steadyatom\Delta N_{ex,\textrm{steady}}^{\textrm{atom}}.

We now discuss quantitatively the characteristics of the superadiance and the memory effects of the early-stage dynamics in near-Markovian regime. In Fig. 8, there are NM(|2,0)=4.7954×1054.8×105N_{M}(|2,0\rangle)=4.7954\times 10^{-5}\approx 4.8\times 10^{-5} and NPsteady(|2,0)=Rsteady(|2,0)=2.3996×1052.4×105N_{P}^{\textrm{steady}}(|2,0\rangle)=R^{\textrm{steady}}(|2,0\rangle)=2.3996\times 10^{-5}\approx 2.4\times 10^{-5}, where R=dNexatom/dtR=dN_{ex}^{\textrm{atom}}/dt is the photon emission rate of the atoms. Meanwhile, the difference of cavity photon number ΔNexatom(τ10,τ21)\Delta N_{ex}^{\textrm{atom}}(\tau_{10},\tau_{21}) in Fig. 8(c) and the trace distance D(τ10,τ21)D(\tau_{10},\tau_{21}) in Fig. 8(d) are approximately the same. This implies that the difference of excitation numbers can be used to measure the memory effects for the Dicke states instead of the trace distance for simplicity. Further simulations imply that NM(|2,0)8(g/γ)2(J+M)[(JM)+1]N_{M}(|2,0\rangle)\rightarrow 8(g/\gamma)^{2}(J+M)[(J-M)+1] and NPsteady(|2,0)=Rsteady4(g/γ)2(J+M)[(JM)+1]N_{P}^{\textrm{steady}}(|2,0\rangle)=R^{\textrm{steady}}\rightarrow 4(g/\gamma)^{2}(J+M)[(J-M)+1] as γ/g\gamma/g goes larger and larger (the Markovian limit). For example, for γ/g=104\gamma/g=10^{4}, there are NM(|2,0)=4.7999×107N_{M}(|2,0\rangle)=4.7999\times 10^{-7} and NPsteady(|2,0)=Rsteady(|2,0)=2.39999×105N_{P}^{\textrm{steady}}(|2,0\rangle)=R^{\textrm{steady}}(|2,0\rangle)=2.39999\times 10^{-5} for the same time interval as in Fig. 8.

Inspired by the above results, we conducted a number of simulations in the near-Markovian regime (with different γ/g\gamma/g) where the initial states included the (dephased) Dicke states with different JJ, MM, and λ\lambda, and the factorized identical states with different NN, ρee\rho_{ee}, and ρeg\rho_{eg}. As in former sections, we find that the characteristics of superradiance and memory effects in the early stage of a near-Markovian dynamics can be well approximated by analytic expressions. These characteristics include the steady cavity photon number of a common cavity (NPsteadyN_{P}^{\textrm{steady}}) and independent cavities (NPsteady,ind=nbnbnsteadyN_{P}^{\textrm{steady},\textrm{ind}}=\langle\sum_{n}b_{n}^{\dagger}b_{n}\rangle_{steady}), the value of memory effects for a common cavity (NMN_{M}) and independent cavities (NMindN_{M}^{\textrm{ind}}), the degree of superradiance S=σ+σ/nσn+σn=NPsteady/NPstead,ind=Rsteady/Rsteady,indS=\langle\sigma^{+}\sigma^{-}\rangle/\sum_{n}\langle\sigma_{n}^{+}\sigma_{n}^{-}\rangle=N_{P}^{\textrm{steady}}/N_{P}^{\textrm{stead},\textrm{ind}}=R^{\textrm{steady}}/R^{\textrm{steady},\textrm{ind}}, and the degree of memory-effect-enhancement NM/NMindN_{M}/N_{M}^{\textrm{ind}}. Using τE\tau_{E} instead of 1/γ1/\gamma for formal correspondence, these early-stage characteristics of the superradiance and the memory effects in the near-Markovian regime are listed in Table 1 (dephased Dicke states) and 2 (factorized identical states) and compared with the early-time results. As shown in Table 1 and 2, the relations between the early-stage Markovian results and the early-time results can be summarized by

NPsteady[ρS(0)](gτE)2\displaystyle\frac{N_{P}^{\textrm{steady}}[\rho_{S}(0)]}{(g\tau_{E})^{2}} =\displaystyle= 4NP[ρS(0)](gt)2,\displaystyle 4\frac{N_{P}[\rho_{S}(0)]}{(gt)^{2}}, (62)
NM[ρS(0)](gτE)2\displaystyle\frac{N_{M}[\rho_{S}(0)]}{(g\tau_{E})^{2}} =\displaystyle= 16NM[ρS(0)](gt)2\displaystyle 16\frac{N_{M}[\rho_{S}(0)]}{(gt)^{2}} (63)

for both collective and independent radiations and the two types of initial states. Therefore, the degree of superradiance SS in the early stage of a Markovian regime is consistent with that in the early-time regime, so does the degree of memory-effect-enhancement NM/NMindN_{M}/N_{M}^{\textrm{ind}}. Besides, for the Dicke states, there are ΔNexatom,steady2NPsteadyNM\Delta N_{ex}^{\textrm{atom},\textrm{steady}}\approx 2N_{P}^{\textrm{steady}}\approx N_{M} and ΔNexatom(τ10,τ21)D(τ10,τ21)\Delta N_{ex}^{\textrm{atom}}(\tau_{10},\tau_{21})\approx D(\tau_{10},\tau_{21}) in the early stage of the near-Markovian dynamics. In comparison, for the factorized identical states, the following relations holds approximately in the early stage of the near-Markovian dynamics: ΔNexatom,steady2NPsteady\Delta N_{ex}^{\textrm{atom},\textrm{steady}}\approx 2N_{P}^{\textrm{steady}} and ΔNexatom(τ10,τ21)D(τ10,τ21)\Delta N_{ex}^{\textrm{atom}}(\tau_{10},\tau_{21})\propto D(\tau_{10},\tau_{21}). The results in this section demonstrate that the characteristics of memory effects and superradiance in the early stage of a near-Markovian dynamics can be well reflected by its early-time dynamics.

V.2 Strongly non-Markovian regime

In this section, we continue to discuss the influence of the initial states on the superradiance and memory effects in a strongly non-Markovian regime where gg and γ\gamma is comparable and 0t<0\leqslant t<\infty are considered. In this regime, typically there are backflows of photons from the cavity to the atoms and the values of memory effects are much higher than those in the early-time regime or the near-Markovian regime. It is worth noting that, in this regime, the superradiance characteristics, such as the atom emission rate R=dNexatom/dtR=dN_{ex}^{\textrm{atom}}/dt and the cavity photon number NPN_{P}, are neither quadratic nor steady in time. Therefore, we may discuss their maxima (denoted by RmaxR^{\textrm{max}} and NPmaxN_{P}^{\textrm{max}}) and use S=Rmax/Rmax,indS=R^{\textrm{max}}/R^{\textrm{max},\textrm{ind}} as the degree of superradiance. Moreover (unlike what happens in other regimes), the initial states usually have evolved dramatically when the atom emission rate RR or other quantities reach their maxima.

Refer to caption
Figure 9: An example of the system and environment dynamics as well as the memory effect manifestations in a strongly non-Markovian regime with γ/g=0.5\gamma/g=0.5 and a common cavity. The initial state is |JM=|2,0|JM\rangle=|2,0\rangle. (a) shows the evolutions of the atom excitation number NexatomN_{ex}^{\textrm{atom}} (solid line) and the emission rate of atoms R=dNexatom/dtR=dN_{ex}^{\textrm{atom}}/dt (dashed line). (b) shows the evolution of the cavity photon number NPN_{P}. (c) and (d) show the difference of the excitation numbers ΔNexatom\Delta N_{ex}^{\textrm{atom}} and the trace distance DD, respectively, as functions of τ10\tau_{10} and τ21\tau_{21} (both in units of 1/g1/g).

As in last section, we provide a typical example in Fig. 9 showing the dynamics of the atoms and one common cavity as well as the memory effect manifestations in a strongly non-Markovian regime with γ/g=0.5\gamma/g=0.5. The results are obtained by numerically solving Eqs. (14) and (15). The initial state is the same as that in Fig. 8. As shown in Fig. 9(a), there are backflows of photons from the cavity to the atoms and the maximal atoms emission rate happens during the first emission pulse. Besides, the maximal cavity photon number NPmaxN_{P}^{\textrm{max}} happens at the end of the first emission pulse as shown in Fig. 9(b). The trace distance D=12ρS(t2)ρS(t2)D=\frac{1}{2}\|\rho_{S}(t_{2})-\rho^{\prime}_{S}(t_{2})\| as a function of τ10\tau_{10} and τ21\tau_{21} are shown in Fig. 9(d) where NM(|2,0)=maxτ10,τ21D(τ10,τ21)0.89N_{M}(|2,0\rangle)=\max_{\tau_{10},\tau_{21}}D(\tau_{10},\tau_{21})\approx 0.89, implying that the dynamics is far from Markovian. Figure 9(d) also demonstrate that the optimal τ10\tau_{10} and τ21\tau_{21} for the global maximum of D(τ10,τ21)D(\tau_{10},\tau_{21}) can be found in a limited time interval. The manifestation of memory effects by ΔNexatom\Delta N_{ex}^{\textrm{atom}} as a function of τ10\tau_{10} and τ21\tau_{21} is shown in Fig. 9(c). Although the surface shapes in Figs. 9(c) and 9(d) are similar, the value of ΔNexatom\Delta N_{ex}^{\textrm{atom}} are not bounded as DD is. Further simulations with other initially excited Dicke states and factorized identical states (with γ/g=0.5\gamma/g=0.5) show similar results as shown in Fig. 9, e.g., there are backflows of photons and RmaxR^{\textrm{max}} (NPmaxN_{P}^{\textrm{max}}) happens during (after) the first emission pulse. Thus the degree of superradiance S=Rmax/Rmax,indS=R^{\textrm{max}}/R^{\textrm{max},\textrm{ind}} can be clearly defined in our simulations.

Refer to caption
Figure 10: Some memory effect and superradiance characteristics in the strongly non-Markovian regime with γ/g=0.5\gamma/g=0.5 for the Dicke states. The number of atoms NN varies from 1 to 10, correspondingly, J=1/2,3/2,,5J=1/2,3/2,\cdots,5. (a) and (b) show the value of memory effects NMN_{M} and the maximal atoms emission rate RmaxR^{\textrm{max}}, respectively, as functions of JJ and MM. Similarly, (c) and (d) show the degree of superradiance S=Rmax/Rmax,indS=R^{\textrm{max}}/R^{\textrm{max},\textrm{ind}} and the degree of the memory-effect-enhancement NM/NMindN_{M}/N_{M}^{\textrm{ind}}, respectively, where M=J+1,J+2,,JM=-J+1,-J+2,\cdots,J

We numerically calculated the value of memory effects NMN_{M}, the maximal atom emission rate RmaxR^{\textrm{max}}, the degree of superradiance SS and the degree of memory-effect-enhancement NM/NMindN_{M}/N_{M}^{\textrm{ind}} for the Dicke states and the factorized identical states with γ/g=0.5\gamma/g=0.5. The results for the Dicke states are shown in Fig. 10 with no more than 1010 atoms. In the strongly non-Markovian regime, the values of memory effect are close to 1 for most of the Dicke states as shown in Fig. 10(a), implying that the Markovian approximation are invalid for their dynamics. Unlike the results in the early-time regime or in the near-Markovian regime, NM(|JM)N_{M}(|JM\rangle) is not proportional to (J+M)(JM+1)(J+M)(J-M+1) here. The first reason is that it tends to saturate (NM[ρS(0)]1N_{M}[\rho_{S}(0)]\leqslant 1) for highly non-Markovian dynamics. The second reason is as follows: on the timescale of the optimal τ10\tau_{10} (τ21\tau_{21}) that maximizing D(τ10,τ21)D(\tau_{10},\tau_{21}), the initial states have evolved dramatically. Thus the connection between NM[ρS(0)]N_{M}[\rho_{S}(0)] and ρS(0)\rho_{S}(0) is weaker compared with those in the other two regimes we discussed. However, the dependence of NM(|JM)N_{M}(|JM\rangle) on JJ and MM is still roughly similar to that in Fig. 4. For example, for a fixed JJ, the strongest memory effects happens for M=0M=0 or near 0. In the strongly non-Markovian regime, the radiation intensity is represented by the maximal emission rate RmaxR^{\textrm{max}} as shown in Fig. 9(a) (dashed line). For J5/2J\leqslant 5/2, the state |J,J|J,J\rangle have the strongest emission intensity with a given JJ. As JJ increases, the states with the strongest emission intensity become |J,J1|J,J-1\rangle as shown in Fig. 9(b). We infer that a large JJ, the Dicke states with M0M\sim 0 might be the most radiative ones. Compared with the results in Fig. 4, the Dicke states with small JJ and large MM are more radiative than those with M0M\sim 0. The reason is as mentioned before, i.e., the states with large MM may have evolved to more superradiant ones when RmaxR^{\textrm{max}} happens. Figures 10(c) and 10(d) show the degree of superradiance and the degree of memory-effect-enhancement. They are not proportional to JM+1J-M+1 as in Fig. 5 in this regime. However, the highest degree of superradiance as well as memory-effect-enhancement in Figs. 10(c) and 10(d) still happens for the state with the largest JJ and M=J+1M=-J+1. Other simulations show that, for the Dicke state |JM=|J,J+1|JM\rangle=|J,-J+1\rangle (one excitation) in a common cavity, the memory effects for the dephased Dicke states satisfies NM[ρSdeph(0)]=λNM(|JMJM|)+(1λ)NM[𝒟(|J,MJM|)]N_{M}[\rho^{\textrm{deph}}_{S}(0)]=\lambda N_{M}(|JM\rangle\langle JM|)+(1-\lambda)N_{M}[\mathcal{D}(|J,M\rangle\langle JM|)], i.e., NM[ρSdeph(0)]N_{M}[\rho^{\textrm{deph}}_{S}(0)] increases linearly with λ\lambda. For other Dicke states, this relation holds approximately at least for a few atoms.

Refer to caption
Figure 11: Memory effects NMN_{M} (left column) and the maximal atom emission rate RmaxR^{\textrm{max}} (right column) for ρSfact(0)\rho_{S}^{\textrm{fact}}(0) as functions of ρee\rho_{ee} and |ρeg||\rho_{eg}| in the strongly non-Markovian regime. The atom numbers are N=2N=2 and 77.

The value of memory effects and the maximal atom emission rate for the factorized identical states are shown in Fig. 11 with two and seven atoms as functions of ρee\rho_{ee} and |ρeg||\rho_{eg}|. In the strongly non-Markovian regime, the population ρee\rho_{ee}, the coherence |ρeg||\rho_{eg}| and the atom number NN can increase the memory effects as well as the maximal atom emission rate. The results show some resemblance to those in Fig. 6 but with stronger memory effects and emission intensities in general. When ρee\rho_{ee} is small, the memory effects are still weak as the system almost stays in its ground state, which can be approximated by a Markovian dynamics in a trivial sense. The degree of superradiance S=Rmax/Rmax,indS=R^{\textrm{max}}/R^{\textrm{max},\textrm{ind}} and the degree of memory-effect-enhancement NM/NMindN_{M}/N_{M}^{\textrm{ind}} are shown in Fig. 12. In contrast with the results in Fig. 7, SS may be greater (less) than 1 for states with ρeg=0\rho_{eg}=0. Because when RmaxR^{\textrm{max}} or Rmax,indR^{\textrm{max},\textrm{ind}} happen, these initial states may evolve to more (less) superradiant ones. Similarly, NM/NMindN_{M}/N_{M}^{\textrm{ind}} may also be less than 1 in the strongly non-Markovian regime. However, as in Fig. 7, the coherence |ρeg||\rho_{eg}| can still increase the degree of superradiance, and initial states with low-excitations and strong coherence have large values of SS and NM/NMindN_{M}/N_{M}^{\textrm{ind}}. Besides, the surface shapes in Fig. 12 (right column) are roughly similar to those in Fig. 12 (left column). Other calculations demonstrate that, for both |JM|JM\rangle and ρSfact(0)\rho_{S}^{\textrm{fact}}(0) as initial states, the following relation holds roughly: NPmaxΔNexatom,maxRmaxN_{P}^{\textrm{max}}\propto\Delta N_{ex}^{\textrm{atom},\textrm{max}}\propto R^{\textrm{max}} in our simulations. Considering the weaker connection between the discussed characteristics and the initial states, we may say that the results in the strongly non-Markovian regime are consistent with those in the early-time regime.

Refer to caption
Figure 12: The degree of superradiance SS (left column) and the degree of memory-effect-enhancement NM/NMindN_{M}/N_{M}^{\textrm{ind}} (right column) for ρSfact(0)\rho_{S}^{\textrm{fact}}(0) as functions of ρee\rho_{ee} and |ρeg||\rho_{eg}| in the strongly non-Markovian regime. The atom numbers are N=2N=2 and 77. The values of SS and NM/NMindN_{M}/N_{M}^{\textrm{ind}} for ρee=0\rho_{ee}=0 are not defined and not shown.

VI Discussions and Conclusions

In this paper, we propose a method to evaluate the memory effects in a quantum process conditioned on a particular system initial state. The method is based on the physical interpretations of the inequality T(t2,t0)T(t2,t1)T(t1,t0)T(t_{2},t_{0})\neq T(t_{2},t_{1})T(t_{1},t_{0}) [35]. Some features of the non-Markovianity measure in Ref. [35] are inherited. For example, nonzero memory effects can be characterized even in regimes where nonmonotonic behaviors do not occur, or have not occurred yet. Besides, 0NM[ρS(t0)]10\leqslant N_{M}[\rho_{S}(t_{0})]\leqslant 1 is satisfied without additional normalizations. This allows us to compare the memory effects of quantum systems with different dimensions (e.g., different numbers of atoms in our model). Using our method, we calculate the influence of two types of initial states on the memory effects as well as the superradiance in different regimes. The characteristics of the memory effects and the superradiance are compared to reveal the role of memory effects in superradiance, and conversely, to understand the physical sources of the memory effects.

In the early-time (Zeno) regime and the (early-stage) near-Markovian regime where the initial states do not change so much, the main observations and conclusions are as follows. The memory effects for the (dephased) Dicke states are fully determined by the cavity photon number (connected to the atom radiation intensity) and they are proportional to (gt)2(J+M)[(JM)λ+1](gt)^{2}(J+M)[(J-M)\lambda+1] for quadratic dynamics, or (g/γ)2(J+M)[(JM)λ+1](g/\gamma)^{2}(J+M)[(J-M)\lambda+1] in the Markovian limit. Meanwhile, the degree of superradiance is equal to the degree memory-effect-enhancement (from independent radiation to collective radiation) given by (JM)λ+1(J-M)\lambda+1. For the factorized identical states, the memory effects, the radiation intensity and the degree of superradiance can all be enhanced by the single-atom coherence. The memory effects are closely related to the cavity photon number, meanwhile, the superradiance (S>1S>1) of an initial state is accompanied by the enhancement of memory effects (by a common cavity). These results demonstrate that the memory effects are closely related to the radiation intensity. Correspondingly, the enhancement of memory effects by a common environment plays an important role on the establishment of superradiance, especially for the Dicke states. On the other hand, these results imply that the (change of) cavity photon number is one fundamental source of memory effects, especially for the Dicke states. Besides, the entanglement in the Dicke states and the single-atom coherence in the factorized identical state are necessary for the existence of superradiance and important for the memory-effect-enhancement.

For the whole dynamics in a strongly non-Markovian regime, an initial state may evolves dramatically during the radiation lifetime, making its influence weaker. However, simulations demonstrate that the relationships between the radiation intensity (degree of superradiance) and the memory effects (memory-effect-enhancement) are basically consistent with those in the other two regimes. The results in different regimes demonstrate that the characteristics of memory effects and superradiance in long-time dynamics can be reflected by its early-time behaviours. In addition, our work also provides other ways to measure the memory effects through the atom excitation number or the cavity photon number.

VII ACKNOWLEDGMENTS

The authors thank X. L. Huang, X. L. Zhao, J. Cheng and B. Cui for helpful discussions. The work is supported by the National Natural Science Foundation of China under Grant No. 11705026, the Fundamental Research Funds for the Central Universities in China under Grant No. 3132020178 and the Natural Science Foundation of Liaoning Province under Grant No. 2023-MS-333.

References

  • [1] S. Cordero and G. García–Calderón, Analytical study of quadratic and nonquadratic short-time behavior of quantum decay, Phys. Rev. A 86, 062116 (2012).
  • [2] S. Wenderoth, H.-P. Breuer, and M. Thoss, Quantifying the influence of the initial state on the dynamics of an open quantum system, Phys. Rev. A 107, 022211 (2023).
  • [3] Á. Rivas, S. F. Huelga, and M. B. Plenio, Quantum non-Markovianity: Characterization, quantification and detection, Rep. Prog. Phys. 77, 094001 (2014).
  • [4] H.-P. Breuer, E.-M. Laine, J. Piilo, and B. Vacchini, Colloquium: Non-Markovian dynamics in open quantum systems, Rev. Mod. Phys. 88, 021002 (2016).
  • [5] I. de Vega and D. Alonso, Dynamics of non-Markovian open quantum systems, Rev. Mod. Phys. 89, 015001 (2017).
  • [6] L. Li, M. J. W. Hall, and H. M. Wiseman, Concepts of quantum non-Markovianity: A hierarchy, Phys. Rep. 759, 1 (2018).
  • [7] R. H. Dicke, Coherence in spontaneous radiation processes, Phys. Rev. 93, 99 (1954).
  • [8] A. V. Andreev, V. I. Emel’yanov, and Y. A. Il’inskiĭ, Collective spontaneous emission (Dicke superradiance), Sov. Phys. Usp. 23, 493 (1980).
  • [9] N. E. Rehler and J. H. Eberly, Superradiance, Phys. Rev. A 3, 1735 (1971).
  • [10] M. Gross and S. Haroche, Superradiance: An essay on the theory of collective spontaneous emission, Phys. Rep. 93, 301 (1982).
  • [11] G. S. Agarwal, Master-Equation Approach to Spontaneous Emission, Phys. Rev. A 2, 2038 (1970).
  • [12] B. Bonifacio, P. Schwendimann, and F. Haake, Quantum statistical theory of superradiance. I, Phys. Rev. A 4, 302 (1971).
  • [13] B. Bonifacio, P. Schwendimann, and F. Haake, Quantum statistical theory of superradiance. II, Phys. Rev. A 4, 854 (1971).
  • [14] H. J. Carmichael, Statistical methods in quantum optics 1: master equations and Fokker-Planck equations (Springer, Berlin, 1999).
  • [15] J. Kim, D. Yang, S. Oh, and K. An, Coherent single-atom superradiance, Science 359, 662 (2018).
  • [16] M. O. Scully and A. A. Svidzinsky, The super of superradiance, Science 325, 1510 (2009).
  • [17] J. A. Mlynek, A. A. Abdumalikov, C. Eichler, and A. Wallraff, Observation of Dicke superradiance for two artificial atoms in a cavity with high decay rate, Nat. Commun. 5, 5186 (2014).
  • [18] A. Rastogi, E. Saglamyurek, T. Hrushevskyi, and L. J. LeBlanc, Superradiance-Mediated Photon Storage for Broadband Quantum Memory, Phys. Rev. Lett. 129, 120502 (2022).
  • [19] F. Robicheaux, Theoretical study of early-time superradiance for atom clouds and arrays, Phys. Rev. A 104, 063706 (2021).
  • [20] D. Mogilevtsev, E. Garusov, M. V. Korolkov, V. N. Shatokhin, and S. B. Cavalcanti, Restoring the Heisenberg limit via collective non-Markovian dephasing, Phys. Rev. A 98, 042116 (2018).
  • [21] M. Bojer and J. von Zanthier, Dicke-like superradiance of distant noninteracting atoms, Phys. Rev. A 106, 053712 (2022).
  • [22] V. Paulisch, M. Perarnau-Llobet, A. González-Tudela, and J. I. Cirac, Quantum metrology with one-dimensional superradiant photonic states, Phys. Rev. A 99, 043807 (2019).
  • [23] W. Cheng, S. C. Hou, Z. Wang, and X. X. Yi, Phys. Rev. A 100, 053825 (2019).
  • [24] K. C. Tan, S. Choi, H. Kwon, and H. Jeong, Coherence, quantum Fisher information, superradiance, and entanglement as interconvertible resources, Phys. Rev. A 97, 052304 (2018).
  • [25] E. Wolfe and S. F. Yelin, Certifying Separability in Symmetric Mixed States of N Qubits, Phys. Rev. Lett. 112, 140402 (2014).
  • [26] M. E. Tasgin, Many-Particle Entanglement Criterion for Superradiantlike States, Phys. Rev. Lett. 119, 033601 (2017).
  • [27] F. Lohof, D. Schumayer, D. A. W. Hutchinson, and C. Gies, Signatures of Superradiance as a Witness to Multipartite Entanglement, Phys. Rev. Lett. 131, 063601 (2023).
  • [28] F. Dinc and A. M. Brańczyk, Non-Markovian super-superradiance in a linear chain of up to 100 qubits, Phys. Rev. Research 1, 032042(R) (2019).
  • [29] K. Sinha, P. Meystre, E. A. Goldschmidt, F. K. Fatemi, S. L. Rolston, and P. Solano, Non-Markovian collective emission from macroscopically separated emitters, Phys. Rev. Lett. 124, 043603 (2020).
  • [30] Q-.Y. Qiu, Y. Wu, and X-.Y. Lü, Collective radiance of giant atoms in non-Markovian regime, Sci. China Phys. Mech. Astron. 66, 224212 (2023).
  • [31] Y. X. Zhang, Zeno regime of collective emission: Non-Markovianity beyond retardation, Phys. Rev. Lett. 131, 193603 (2023).
  • [32] H.-P. Breuer, E.-M. Laine, and J. Piilo, Measure for the Degree of Non-Markovian Behavior of Quantum Processes in Open Systems, Phys. Rev. Lett. 103, 210401 (2009).
  • [33] Á. Rivas, S. F. Huelga, and M. B. Plenio, Entanglement and Non-Markovianity of Quantum Evolutions, Phys. Rev. Lett. 105, 050403 (2010).
  • [34] D. Chruściński and A. Kossakowski, Non-Markovian Quantum Dynamics: Local versus Nonlocal, Phys. Rev. lett. 104, 070406 (2010).
  • [35] S. C. Hou, S. L. Liang, and X. X. Yi, Non-Markovianity and memory effects in quantum open systems, Phys. Rev. A. 91, 012109 (2015).
  • [36] Á. Rivas and S. F. Huelga, Open Quantum Systems An Introduction (Springer, Heidelberg, 2012).
  • [37] S. C. Hou, X. X. Yi, S.X.Yu, and C.H.Oh, Alternative non-Markovianity measure by divisibility of dynamical maps, Phys. Rev. A 83, 062115 (2011).
  • [38] M.-D. Choi, Completely positive linear maps on complex matrices, Lin. Alg. and Appl. 10, 285 (1975).
  • [39] A. Jamiołkowski, Linear transformations which preserve trace and positive semidefiniteness of operators, Rep. Math. Phys. 3, 275 (1972).
  • [40] P. Facchi, H. Nakazato, and S. Pascazio, From the quantum Zeno to the inverse quantum Zeno effect, Phys. Rev. Lett. 86, 2699 (2001).
  • [41] L. S. Schulman, Observational line broadening and the duration of a quantum jump, J. Phys. A 30, L293 (1997);
  • [42] A. Crespi, F.V. Pepe, P. Facchi, F. Sciarrino, P. Mataloni, H. Nakazato, S. Pascazio, and R. Osellame, Experimental investigation of quantum decay at short, intermediate and long times via integrated photonics, Phys. Rev. Lett. 122, 130401 (2019).
  • [43] W. Wu and H.-Q. Lin, Quantum Zeno and anti-Zeno effects in quantum dissipative systems, Phys. Rev. A 95, 042132 (2017).
  • [44] G. S. Agarwal, Quantum Optics (Cambridge University Press, Cambridge, 2013).