The quantum Zeno and anti-Zeno effects with strong system-environment coupling

TL;DR

This paper explores the quantum Zeno and anti-Zeno effects in strongly coupled systems, revealing that strong interactions lead to non-linear decay rates and indirect atom interactions, contrasting with weak coupling behaviors.

Contribution

It introduces a polaron transformation approach to analyze the Zeno effects under strong system-environment coupling, highlighting qualitative and quantitative differences from weak coupling scenarios.

Findings

Strong coupling alters the decay rate dependence on spectral density.

Increasing coupling strength can decrease the effective decay rate.

Collective interactions emerge among atoms via the common environment.

Abstract

The fact that repeated projective measurements can slow down (the Zeno effect) or speed up (the anti-Zeno effect) quantum evolution is well-known. However, to date, studies of these effects focus on quantum systems that are weakly interacting with their environment. In this paper, we investigate what happens to a quantum system under the action of repeated measurements if the quantum system is strongly interacting with its environment. We consider as the quantum system a single two-level system coupled strongly to a collection of harmonic oscillators. A so-called polaron transformation is then used to make the problem in the strong system-environment coupling regime tractable. We find that the strong coupling case exhibits quantitative and qualitative differences as compared with the weak coupling case. In particular, the effective decay rate does not depend linearly on the spectral density of the environment. This then means that, in the strong coupling regime that we investigate, increasing the system-environment coupling strength can actually decrease the effective decay rate. We also consider a collection of two-level atoms coupled strongly with a common environment. In this case, we find that there are further differences between the weak and strong coupling cases since the two-level atoms can now indirectly interact with one another due to the common environment.

Equations

$$H_L=\frac{\epsilon }{2}{\sigma }_z+\frac{\Delta }{2}{\sigma }_x+\sum _k{\omega }_k{b}_k^{†}{b}_k+{\sigma }_z\sum _k(g_k^{\ast }{b}_k+g_k{b}_k^{†}),$$

$$s(\tau )={\text{Tr}}_{S,B}[|\uparrow ⟩⟨\uparrow |e^{-i{H}_L\tau }{\rho }_L(0)e^{i{H}_L\tau }].$$

$$s(\tau )={\text{Tr}}_{S,B}[|\uparrow ⟩⟨\uparrow |e^{-iH\tau }{P}_{\uparrow }\frac{e^{-\beta H}}{Z}{P}_{\uparrow }{e}^{iH\tau }],$$

$$s(\tau )={\text{Tr}}_{S,B}[|\uparrow ⟩⟨\uparrow |e^{-iH\tau }(|\uparrow ⟩⟨\uparrow |\otimes e^{-\beta H_B}/Z_B)e^{iH\tau }]={\text{Tr}}_S[|\uparrow ⟩⟨\uparrow |{\rho }_S(\tau )],$$

$${\rho }_S(\tau )=$$

$$\displaystyle\sum_{\mu\nu}\int_{0}^{\tau}dt_{1}\int_{0}^{t_{1}}dt_{2}\bigl{\{}C_{\mu\nu}(t_{1},t_{2})[\widetilde{F}_{\nu}(t_{2})\rho_{S}(0),\widetilde{F}_{\mu}(t_{1})]+\text{h.c.}\bigr{\}}\biggr{)}U_{S}^{\dagger}(\tau).$$

$$\Gamma (\tau )$$

$${\Phi }_R(t)={\int }_0^{\infty }d\omega J(\omega )\frac{1-\cos (\omega t)}{{\omega }^2}\coth \left(\frac{\beta \omega }{2}\right),{\Phi }_I(t)={\int }_0^{\infty }d\omega J(\omega )\frac{\sin (\omega t)}{{\omega }^2},$$

$$\Gamma (\tau )$$

$$s(\tau )={\text{Tr}}_{S,B}[(|\uparrow ⟩⟨\uparrow |)e^{i{H}_{S,L}\tau }{e}^{-i{H}_L\tau }{\rho }_{\text{L}}(0)e^{i{H}_L\tau }{e}^{-i{H}_{S,L}\tau }].$$

$$s(\tau )=1-{\text{Tr}}_{S,B}[(|\downarrow ⟩⟨\downarrow |)e^{i{H}_{S,P}\tau }{e}^{-iH\tau }(|\uparrow ⟩⟨\uparrow |\otimes {\rho }_B)e^{iH\tau }{e}^{-i{H}_{S,P}\tau }],$$

$${\Gamma }_n(\tau )=\Gamma (\tau )+{\Gamma }_{\text{mod}}(\tau ),$$

$${\Gamma }_{\text{mod}}(\tau )=\frac{{\Delta }^2}{\tau }\left\{\frac{1}{{\epsilon }^2}{\sin }^2\left(\frac{\epsilon \tau }{2}\right)e^{-{\Phi }_R(0)}{e}^{-i{\Phi }_I(0)}-\frac{1}{\epsilon }\sin \left(\frac{\epsilon \tau }{2}\right){\int }_0^{\tau }dt{e}^{-{\Phi }_R(t)}\cos \left[\epsilon \left(t-\frac{\tau }{2}\right)-{\Phi }_I(t)\right]\right\}.$$

$${\Gamma }_n(\tau )={\int }_0^{\infty }d\omega J(\omega )Q(\omega ,\tau ),$$

$$Q(\omega ,\tau )=\frac{2}{\tau }\left\{\coth \left(\frac{\beta \omega }{2}\right)D_1(\omega ,\tau )+D_2(\omega ,\tau )\right\},$$

$$H_L=\epsilon J_z+\Delta J_x+\sum _k{\omega }_k{b}_k^{†}{b}_k+2J_z\sum _k(g_k^{\ast }{b}_k+g_k{b}_k^{†}),$$

$$s(\tau )={\text{Tr}}_{S,B}[|j⟩⟨j|e^{-i{H}_L\tau }{\rho }_L(0)e^{i{H}_L\tau }].$$

$$H=\epsilon J_z+\sum _k{\omega }_k{b}_k^{†}{b}_k-\kappa J_z^2+\frac{\Delta }{2}(J_{+}X+J_{-}{X}^{†}),$$

$$\Gamma (\tau )=\frac{{\Delta }^{2}j}{\tau }{\int }_0^{\tau }dt{\int }_0^{t}d{t}^{\prime }{e}^{-{\Phi }_R(t^{\prime })}\cos [\epsilon t^{\prime }+\kappa (1-2j)t^{\prime }-{\Phi }_I(t^{\prime })].$$

$${\Gamma }_n(\tau )=\Gamma (\tau )+{\Gamma }_{\text{mod}}(\tau ),$$

$${\Gamma }_{\text{mod}}(\tau )=\frac{{\Delta }^2}{\tau }(2j)\left\{\frac{1}{{\epsilon }^2}{\sin }^2\left(\frac{\epsilon \tau }{2}\right)e^{-{\Phi }_R(0)-i{\Phi }_I(0)}-\frac{1}{\epsilon }\sin \left(\frac{\epsilon \tau }{2}\right){\int }_0^{\tau }dt{e}^{-{\Phi }_R(t)}\cos \left[\kappa (2j-1)(2\tau -t)+\epsilon \left(t-\frac{\tau }{2}\right)-{\Phi }_I(t)\right]\right\}.$$

$$e^{\theta A}B{e}^{-\theta A}=B+\theta [A,B]+\frac{{\theta }^2}{2!}[A,[A,B]]+\dots$$

$$e^{\chi {\sigma }_z/2}\left[\frac{\epsilon }{2}{\sigma }_z+\sum _k{\omega }_k{b}_k^{†}{b}_k+{\sigma }_z\sum _k(g_k^{\ast }{b}_k+g_k{b}_k^{†})\right]e^{-\chi {\sigma }_z/2}=\frac{\epsilon }{2}{\sigma }_z+\sum _k{\omega }_k{b}_k^{†}{b}_k.$$

$$[\chi J_z,\epsilon J_z+\sum _k{\omega }_k{b}_k^{†}{b}_k+2J_z\sum _k(g_k^{\ast }{b}_k+g_k{b}_k^{†})]=-2J_z\sum _k(g_k{b}_k^{†}+g_k^{\ast }{b}_k)-8J_z^2\sum _k\frac{|g_k{|}^2}{{\omega }_k}.$$

$$[\chi J_z,-2J_z\sum _k(g_k{b}_k^{†}+g_k^{\ast }{b}_k)-8J_z^2\sum _k\frac{|g_k{|}^2}{{\omega }_k}]=8J_z^2\sum _k\frac{|g_k{|}^2}{{\omega }_k}.$$

$${\rho }_S(\tau )\approx$$

$${\text{Tr}}_B\{U_0(\tau )A_2{\rho }_{\text{tot}}(0)U_0^{†}(\tau )\}=-\sum _{\mu \nu }{\int }_0^{\tau }d{t}_1{\int }_0^{t_1}d{t}_2{U}_S(\tau ){\tilde{F}}_{\mu }(t_1){\tilde{F}}_{\nu }(t_2){\rho }_S(0)U_S^{†}(\tau )C_{\mu \nu }(t_1,t_2),$$

$${\text{Tr}}_B\{U_0(\tau ){\rho }_{\text{tot}}(0)A_2^{†}{U}_0^{†}(\tau )\}=-\sum _{\mu \nu }{\int }_0^{\tau }d{t}_1{\int }_0^{t_1}d{t}_2{U}_S(\tau ){\rho }_S(0){\tilde{F}}_{\nu }(t_2){\tilde{F}}_{\mu }(t_1)U_S^{†}(\tau )C_{\nu \mu }(t_2,t_1).$$

$${\text{Tr}}_B\{U_0(\tau )A_1{\rho }_{\text{tot}}{A}_1^{†}{U}_0^{†}(\tau )\}=\sum _{\mu \nu }{\int }_0^{\tau }d{t}_1{\int }_0^{\tau }d{t}_2{U}_S(\tau ){\tilde{F}}_{\mu }(t_1){\rho }_S(0){\tilde{F}}_{\nu }(t_2)U_S^{†}(\tau )C_{\nu \mu }(t_2,t_1).$$

$${\text{Tr}}_B\{U_0(\tau )A_1{\rho }_{\text{tot}}{A}_1^{†}{U}_0^{†}(\tau )\}=\sum _{\mu \nu }{\int }_0^{\tau }d{t}_1{\int }_0^{t_1}d{t}_2{U}_S(\tau ){\tilde{F}}_{\mu }(t_1){\rho }_S(0){\tilde{F}}_{\nu }(t_2)U_S^{†}(\tau )C_{\nu \mu }(t_2,t_1)+\text{h.c.},$$

$$s(\tau )$$

Full text

By repeatedly measuring a quantum system very frequently, the evolution of the quantum system can be slowed down, an effect that has been dubbed as the Quantum Zeno effect (QZE)[1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20]. On the other hand, if the quantum system is measured repeatedly not very rapidly, the measurements can actually speed up the temporal evolution. This effect, the opposite of the QZE, is known as the Quantum anti-Zeno effect (QAZE)[21, 22, 23, 24, 25, 26]. Both the QZE and the QAZE have attracted tremendous theoretical and experimental interest due to their great importance for emerging quantum technologies as well as their fundamental theoretical interest. However, it is worth noting that the emphasis in studies performed on the QZE and the QAZE to date has been on the population decay of quantum systems. In these studies, the quantum system is prepared in an excited state, and then the system is repeatedly checked to see if the system is still in the excited state or not[21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32]. It is well-known then that the decay rate of the quantum system depends on the overlap of the spectral density of the environment and a measurement-induced level width[21]. Depending on this overlap, decreasing the measurement interval can lead to a decrease (the QZE) or an increase (the QAZE) of the decay rate.

While studies of the QZE and the QAZE performed to date by and large focus on the population decay model where only decay takes place, we also know from the study of open quantum systems that, in general, quantum systems interacting with their environment also undergo dephasing. To this end, the QZE and the QAZE were studied for the exactly solvable pure dephasing model in Ref. [33], where it was shown that the QZE and the QAZE are significantly different for the pure dephasing case as compared with the population decay case. This study was then extended to arbitrary system-environment models in Ref. [34] where a general framework for calculating the effective decay rate of the system for an arbitrary system-environment model was presented. It was found that the effective decay rate can be written as an overlap integral of the spectral density of the environment and an effective ‘filter function’ that depends on the system-environment model at hand, the measurement interval, and the measurement being repeatedly performed. This general formalism was then used to study the QZE and the QAZE when both dephasing and population decay are present. For example, repeated measurements for the paradigmatic spin-boson model [35] were considered and it was shown that the presence of both population decay and dephasing make the results differ considerably both quantitatively and qualitatively as compared to the pure population decay case.

It should be pointed out, however, that the results presented in Ref. [34] were derived under the assumption that the system-environment coupling is weak. This is consistent with studies performed for the population decay models, where the effective decay rate can be derived using time-dependent perturbation theory [22]. On the other hand, the behavior of a quantum system, subjected to repeated measurements, that is interacting strongly with its environment is not well understood. For instance, one could ask whether or not the effective decay rate is still an overlap integral of the spectral density function and a ‘filter’ function. This paper intends to answer precisely such questions by looking at what happens to the spin-boson model under the action of repeated measurements if the central two-level system is interacting strongly with a surrounding environment of harmonic oscillators. Since the system-environment coupling is strong, the system-environment interaction cannot be treated perturbatively, and thus the treatment given in Ref. [34] is no longer applicable. Our strategy then is to perform a unitary transformation, known as the polaron transformation, on the system-environment Hamiltonian [36, 37, 38, 39, 40, 41]. One then finds that the system and the environment can end up interacting weakly in this new ‘polaron’ frame. Perturbation theory can then be applied and the effect of repeated measurements is analyzed. We find that the analysis of the QZE and QAZE are in general very different compared to the population decay case. For example, it is clear that for the usual population decay case, increasing the system-environment strength increases the effective decay rate. However, for the strong system-environment regime that we investigate, we find that increasing the system-environment coupling regime can actually decrease the effective decay rate. We also study the QZE and the QAZE for more than one two-level system interacting with a common environment. For the weak coupling regime, the effective decay rate is directly proportional to the number of two-level systems coupled to the common environment [34]. On the other hand, for the strong system-environment coupling regime, we find that the effective decay rate for more than one two-level system is very different compared to the single two-level system case. The indirect interaction between the two-level systems due to their interaction with a common environment now plays a very important role, and the effective decay rate is no longer simply proportional to the number of two-level systems coupled to the common environment.

Results

Spin-boson model with strong system-environment coupling

We start with the paradigmatic spin-boson model Hamiltonian [35, 42, 43] which we write as (we set $\hbar=1$ throughout)

[TABLE]

where the system Hamiltonian is $H_{S,L}=\frac{\varepsilon}{2}\sigma_{z}+\frac{\Delta}{2}\sigma_{x}$, the environment Hamiltonian is $H_{B}=\sum_{k}\omega_{k}b_{k}^{\dagger}b_{k}$, and the system-environment coupling is $V_{L}=\sigma_{z}\sum_{k}(g_{k}^{*}b_{k}+g_{k}b_{k}^{\dagger})$. $\varepsilon$ is the energy level difference of the two-level system, $\Delta$ is the tunneling amplitude, $\omega_{k}$ are the frequencies of the harmonic oscillators, $b_{k}$ and $b_{k}^{\dagger}$ are the annihilation and creation operators for the harmonic oscillators, and $\sigma_{x}$ and $\sigma_{z}$ are the standard Pauli operators. The ‘$L$’ denotes the ‘lab’ frame. If the system-environment coupling is strong, we cannot treat the system-environment coupling perturbatively. Furthermore, the system-environment correlation effects are significant as well in general. To motivate our basic approach in this strong coupling regime, we note that if the system tunneling amplitude is negligible and the initial system state is an eigenstate of $\sigma_{z}$, then, even though the system and the environment are strongly interacting, the evolution of the system state is negligible. This then means that we should look to unitarily transform $H_{L}$ such that the effective system-environment coupling contains the tunneling amplitude $\Delta$. This unitary transformation is provided by the ‘polaron’ transformation, whereby the system-environment Hamiltonian in this new ‘polaron’ frame becomes $H=e^{\chi\sigma_{z}/2}H_{L}e^{-\chi\sigma_{z}/2}$, where $\chi=\sum_{k}\left[\frac{2g_{k}}{\omega_{k}}b_{k}^{\dagger}-\frac{2g_{k}^{*}}{\omega_{k}}b_{k}\right]$ [36, 37, 38, 39, 40, 41]. The system-environment Hamiltonian in the polaron frame is then $H=H_{S}+H_{B}+V$, where $H_{S}=\frac{\varepsilon}{2}\sigma_{z}$, $H_{B}=\sum_{k}\omega_{k}b_{k}^{\dagger}b_{k}$, and $V=\frac{\Delta}{2}\left[\sigma_{+}X+\sigma_{-}X^{\dagger}\right]$, with $X=e^{\chi}$ (see Methods for details). Now, if the tunneling amplitude is small, we can use time-dependent perturbation theory, treating $V$ as the perturbation. This is the key idea to deal with the strong system-environment coupling regime. Although the system and the environment are strongly interacting, in the polaron frame, they are effectively interacting weakly. Let us now use this fact in order to calculate the survival probability, and thereby the effective decay rate. For concreteness, we assume that the initial state prepared is $\left|\uparrow\right\rangle$, where $\sigma_{z}\left|\uparrow\right\rangle=\left|\uparrow\right\rangle$. In other words, we consider the same initial state as that considered in the analysis of the usual population decay model [21, 22]. At time $t=0$, we prepare the system state $\left|\uparrow\right\rangle$, and we subsequently perform measurements with time interval $\tau$ to check if the system state is still $\left|\uparrow\right\rangle$ or not. The survival probability after time interval $\tau$ is then $s(\tau)=\text{Tr}_{S,B}[\left|\uparrow\right\rangle\left\langle\uparrow\right|\rho_{L}(\tau)]$, where $\rho_{L}(\tau)$ is the combined density matrix of the system and the environment at time $\tau$ just before the projective measurement. Then,

[TABLE]

It is important to note that the initial state that we have prepared cannot simply be taken as the usual product state $\left|\uparrow\right\rangle\left\langle\uparrow\right|\otimes e^{-\beta H_{B}}/Z_{B}$, with $Z_{B}=\text{Tr}_{B}[e^{-\beta H_{B}}]$ since the system and the environment are strongly interacting and consquently there will be significant initial system-environment correlations [44, 45] . Rather, the initial state that we should consider is $\rho_{L}(0)=P_{\uparrow}e^{-\beta H_{L}}P_{\uparrow}/Z$, where $P_{\uparrow}=\left|\uparrow\right\rangle\left\langle\uparrow\right|$, and $Z=\text{Tr}_{S,B}[P_{\uparrow}e^{-\beta H_{L}}]$. Keeping this in mind, we use the polaron transformation to cast the expression for the survival probability $s(\tau)$ after the measurement at time $\tau$ in terms of quantities in the polaron frame. Doing so leads us to

[TABLE]

with the Hamiltonian $H$ now in the polaron frame. Now, for small $\Delta$, to a first approximation, the initial state in the polaron frame can be written as $P_{\uparrow}\otimes e^{-\beta H_{B}}/Z_{B}$. This is a similar approximation as the usual assumption that the initial system-environment state is $\rho_{S}(0)\otimes\rho_{B}$ since, in the polaron frame, the system and the environment are weakly interacting. We thus get

[TABLE]

where $\rho_{S}(\tau)=\text{Tr}_{B}[e^{-iH\tau}(\left|\uparrow\right\rangle\left\langle\uparrow\right|\otimes e^{-\beta H_{B}}/Z_{B})e^{iH\tau}]$. Our objective then is to find $\rho_{S}(\tau)$, given the initial system-environment state $\rho(0)=\left|\uparrow\right\rangle\left\langle\uparrow\right|\otimes e^{-\beta H_{B}}/Z_{B}$. We find that (see the Methods section)

[TABLE]

Here $U_{S}(\tau)=e^{-iH_{S}\tau}$, $F_{1}=\frac{\Delta}{2}\sigma_{+}$, $B_{1}=X$, $F_{2}=\frac{\Delta}{2}\sigma_{-}$, $B_{2}=X^{\dagger}$, $\widetilde{F}_{\mu}(t)=U_{S}^{\dagger}(t)F_{\mu}U_{S}(t)$, $\widetilde{B}_{\mu}(t)=U_{B}^{\dagger}(t)B_{\mu}U_{B}(t)$ with $U_{B}(t)=e^{-iH_{B}t}$, $\langle\ldots\rangle_{B}=\text{Tr}_{B}[\rho_{B}(\ldots)]$ where $\text{Tr}_{B}$ denotes taking trace over the environment, the environment correlation functions are defined as $C_{\mu\nu}(t_{1},t_{2})=\langle\widetilde{B}_{\mu}(t_{1})\widetilde{B}_{\nu}(t_{2})\rangle_{B}$, and h.c. denotes the hermitian conjugate. Now, since the system-environment coupling in the polaron frame is weak, we can neglect the build up of correlations between the system and the environment. Thus, we can write the survival probability after time $t=N\tau$, where $N$ is the number of measurements performed after time $t=0$, as $S(t=N\tau)=[s(\tau)]^{N}\equiv e^{-\Gamma(\tau)N\tau}$, thereby defining the effective decay rate $\Gamma(\tau)$. It then follows that $\Gamma(\tau)=-\frac{1}{\tau}\ln s(\tau)$. Since we have the system density matrix in the polaron frame, we can work out the survival probability $s(\tau)$ and hence the effective decay rate $\Gamma(\tau)$. The result is that (see the Methods section for details)

[TABLE]

where

[TABLE]

and the spectral density of the environment has been introduced as $\sum_{k}|g_{k}|^{2}(\ldots)\rightarrow\int_{0}^{\infty}d\omega J(\omega)(\ldots)$. At this point, it is useful to compare this expression for the effective decay rate for the case of strong system-environment coupling with the case of the usual population decay model where the effective decay rate is $\Gamma(\tau)=\tau\int_{0}^{\infty}\,d\omega\,J(\omega)\text{sinc}^{2}\left[\frac{(\varepsilon-\omega)\tau}{2}\right]$[21, 22]. It should be clear that for the strong system-environment coupling case, the effective decay rate given by Eq. (3) has a very different qualitative behavior. In particular, the effective decay rate can no longer be regarded as simply an overlap integral of the spectral density of the environment with a sinc-squared function. Rather, the effective decay rate now has a very prominent non-linear dependence on the spectral density, leading to very different behavior as compared with the population decay case. For example, as the system-environment coupling strength increases, $\Phi_{R}(t)$ increases, and thus we expect $\Gamma(\tau)$ to decrease. To make this claim concrete, let us model the spectral density as $J(\omega)=G\omega^{s}\omega_{c}^{1-s}e^{-\omega/\omega_{c}}$, where $G$ is a dimensionless parameter characterizing the system-environment coupling strength, $\omega_{c}$ is the cutoff frequency, and $s$ is the Ohmicity parameter [43]. For concreteness, we look at the Ohmic case $(s=1)$. In this case, $\Phi_{R}(t)=\frac{G}{2}\ln(1+\omega_{c}^{2}t^{2})$, while $\Phi_{I}(t)=G\tan^{-1}(\omega_{c}t)$, leading to

[TABLE]

The double integral can be worked out numerically. Results are shown in Fig. 1(a) for different system-environment coupling strengths $G$. For the strong system-environment regime that we are dealing with, it is clear that increasing the system-environment coupling strength $G$ actually decreases the effective decay rate. This is in contrast with what happens in the weak system-environment regime for the paradigmatic population decay model [see Fig. 1(b)]. Here it is clear that increasing the system-environment coupling strength increases the effective decay rate as expected. It should also be noted that the behaviour of $\Gamma(\tau)$ as a function of $\tau$ allows us to identify the Zeno and anti-Zeno regimes. One approach is to simply say that if $\Gamma(\tau)$ decreases when $\tau$ decreases, we are in the Zeno regime, while if $\Gamma(\tau)$ increases if $\tau$ decreases, then we are in the anti-Zeno regime [21, 28, 31, 33]. From Fig. 1(b), it should also be noted that increasing the coupling strength does not change the qualitative behavior of the Zeno to anti-Zeno transition, but for the strong coupling regime [see Fig. 1(a)], while we only observe the Zeno effect for $G=1$, both the Zeno and anti-Zeno effects are observed for $G=2.5$. Similarly, as shown in Fig. 2(a), increasing the cutoff frequency for the strong coupling case decreases the effective decay rate, but the opposite behaviour is observed for the weak coupling case [see Fig. 2(b)].

In our treatment until now, we have considered the change in the system state due to the tunneling term. This tunneling term, due to its presence in $H_{S,L}=\frac{\varepsilon}{2}\sigma_{z}+\frac{\Delta}{2}\sigma_{x}$, leads to the system state changing even if the system and the environment are not coupled to each other. Thus, an alternative way to quantify the effective decay rate would be to remove the evolution due to the system Hamiltonian (in the ‘lab’ frame) $H_{S,L}$ before performing each measurement since what we are really interested in is the change in the system state due to the system-environment interaction. A similar approach has been followed in Refs. [46, 33, 34]. Therefore, we now derive an expression for the effective decay rate of the system state when, just before each measurement, we remove the system evolution due to $H_{S,L}$. The survival probability, after one measurement, is now (starting from the state $\left|\uparrow\right\rangle$)

[TABLE]

Notice now the presence of $e^{iH_{S,L}\tau}$ and $e^{-iH_{S,L}\tau}$ which remove the evolution of the system due to the system Hamiltonian before performing the measurement. Once again transforming to the polaron frame, we obtain

[TABLE]

where $H_{S,P}=\frac{\varepsilon}{2}\sigma_{z}+\frac{\Delta}{2}(\sigma_{+}X+\sigma_{-}X^{\dagger})$ and $\rho_{B}=e^{-\beta H_{B}}/Z_{B}$. Since we are assuming that the tunneling amplitude is small, the unitary operator $e^{-iH_{S,P}\tau}$ can be expanded as a perturbation series. At the same time, $e^{-iH\tau}$ can also expanded as a perturbation series. Keeping terms to second order in the tunneling amplitude (see the Methods section), we find that now the modified decay rate $\Gamma_{n}(\tau)$ is

[TABLE]

where the modification to the previous decay rate is

[TABLE]

Using these expressions, we have plotted the behavior of $\Gamma_{n}(\tau)$ for the strong system-environment coupling regime in Fig. 3(a). It should be clear that once again increasing the system-environment coupling strength generally decreases the effective decay rate $\Gamma_{n}(\tau)$. This is in sharp contrast with what happens in the weak coupling regime. For the weak coupling case, it is known that [34]

[TABLE]

where the filter function $Q(\omega,\tau)$ is

[TABLE]

with $D_{1}(\omega,\tau)=\int_{0}^{\tau}dt\int_{0}^{t}dt^{\prime}\cos(\omega t^{\prime})[a_{x}(t-t^{\prime})a_{x}(t)+a_{y}(t-t^{\prime})a_{y}(t)]$, and $D_{2}(\omega,\tau)=\int_{0}^{\tau}dt\int_{0}^{t}dt^{\prime}\sin(\omega t^{\prime})[-a_{x}(t-t^{\prime})a_{y}(t)+a_{x}(t)a_{y}(t-t^{\prime})]$. Here $a_{x}(t)=\frac{2\varepsilon\Delta}{\Omega^{2}}\sin^{2}\left(\frac{\Omega t}{2}\right)$ and $a_{y}(t)=\frac{\Delta}{\Omega}\sin(\Omega t)$ with $\Omega^{2}=\varepsilon^{2}+\Delta^{2}$. Using these expressions, we can investigate how the decay rate varies as the measurement interval changes for different system-environment coupling strengths in the weak coupling regime. Typical results are illustrated in Fig. 3(b) from which it should be clear that increasing the coupling strength in the weak coupling regime increases the effective decay rate. Furthermore, changing the coupling strength has no effect on the measurement time interval at which the Zeno to anti-Zeno transition takes place for the weak coupling regime as the three curves in Fig. 3(b) achieve their maximum value for the same value of $\tau$. This is not the case for the strong coupling regime [see Fig. 3(a)].

At this point, it is worth pausing to consider where the qualitative difference in the behavior of the effective decay rate in the weak and the strong coupling regime comes from. The effective decay rate is derived from the survival probability after one measurement $s(\tau)$. For both the weak and the strong coupling regimes, the survival probability after one measurement is given by Eq. (6). For both cases, the Hamiltonian $H_{S,L}$ and $H_{L}$ are the same. The only difference is in the choice of the system-environment state $\rho_{L}(0)$. For the weak coupling case, this state is simply the product state $\left|\uparrow\right\rangle\left\langle\uparrow\right|\otimes e^{-\beta H_{B}}/Z_{B}$. This is not the case for the strong coupling due to the significant system-environment correlations. Thus, we can say that the qualitative difference in the behavior of the effective decay rate is because of the presence of the system-environment correlations. It seems that these correlations can protect the quantum state of the system - as the coupling strength increases, these correlations become more and more significant, and at the same time, the effective decay rate goes down.

Large spin-boson model with strong system-environment coupling

Let us now generalize the usual spin-boson model to deal with $N_{S}$ two-level systems interacting with a common environment. In this case, the system-environment Hamiltonian (in the ‘lab’ frame) is given by [37, 45, 34]

[TABLE]

where $J_{x,y,z}$ are the usual angular momentum operators obeying the commutation relations $[J_{k},J_{l}]=i\varepsilon_{klm}J_{m}$. We now start from the spin coherent state $\left|j\right\rangle$ such that $J_{z}\left|j\right\rangle=j\left|j\right\rangle$ with $j=N_{S}/2$. Other eigenstates of $J_{z}$ can be considered as the initial state in a similar manner. Our objective is to again perform repeated projective measurements, described by the projector $\left|j\right\rangle\left\langle j\right|$, with time interval $\tau$ and thereby investigate what happens to the effective decay rate. As before, the survival probability after one measurement is

[TABLE]

Since we consider the system and the environment to be strongly interacting, we once again perform the polaron tranformation given by $H=e^{\chi J_{z}}H_{L}e^{-\chi J_{z}}$, with $\chi$ the same as before. Then, we find that

[TABLE]

where $\kappa=4\sum_{k}\frac{|g_{k}|^{2}}{\omega_{k}}$, and $J_{\pm}=J_{x}\pm iJ_{y}$ are the standard raising and lowering operators. Interestingly, the transformed Hamiltonian now contains a term proportional to $J_{z}^{2}$. This term arises because the collection of two-level systems interacting with the collective environment are indirectly interacting with each other. This term is obviously proportional to the identity operator for a single two-level system, and thus has no influence for a single two-level system. If the tunneling amplitude is small, then we again use perturbation theory and assume that, in the polaron frame, the system-environment correlations can be neglected. We find that now the effective decay rate is (see the Methods section)

[TABLE]

Here $\Phi_{R}(t)$ and $\Phi_{I}(t)$ are the same as defined before. This result obviously agrees with the result that we obtained for a single two-level system. Moreover, it is clear from Eq. (10) that increasing the system-environment coupling strength $G$ should reduce the effective decay rate due to the $e^{-\Phi_{R}(t^{\prime})}$ factor in the integrand. This is precisely what we observe in Fig. 4(a). Furthermore, it may be thought that increasing $j$ (or, equivalently, $N_{S}$) increases the effective decay rate. On the other hand, the dependence on $j$ is not so clear because of the presence of the indirect interaction. Namely, increasing $j$ increases the oscillatory behavior of the integrand due to the dependence of the integrand on $\cos[\varepsilon t^{\prime}+\kappa(1-2j)t^{\prime}-\Phi_{I}(t^{\prime})]$. Thus, once the integral over this rapidly oscillating integrand is taken, we can again get a small number. Such a prediction is borne out by Fig. 4(b) where the effective decay rate has been plotted for different values of $j$. It is obvious that there is a big difference between the single two-level system case and the more than one two-level system case. Furthermore, it seems that increasing $j$ can largely reduce the value of the effective decay rate, meaning that in the strong coupling regime, the indirect interaction helps in keeping the quantum state alive.

Let us now consider the situation where the evolution to the system Hamiltonian $H_{S,L}=\varepsilon J_{z}+\Delta J_{x}$ is removed before each measurement. In the polaron frame $H_{S,L}$ becomes $H_{S,P}=\varepsilon J_{z}+\frac{\Delta}{2}(J_{+}X+J_{-}X^{\dagger})$. The major difference now compared to the previous single two-level system case is that the total system-environment Hamiltonian in the polaron frame $H=H_{S,P}+H_{B}-\kappa J_{z}^{2}$ contains a term (namely, $-\kappa J_{z}^{2}$) that is not part of the system Hamiltonian in the polaron frame. As a result, when the system evolution is removed just before performing each measurement, the evolution induced by this extra term survives. Keeping this fact in mind, the effective decay rate $\Gamma_{n}(\tau)$ is now

[TABLE]

where $\Gamma(\tau)$ is given by Eq. (10) and

[TABLE]

In Fig. 5(a), we have shown the behavior of $\Gamma_{n}(\tau)$ when the system-environment coupling strength is increased for $N_{S}=2$. It should be obvious that we observe multiple Zeno-anti Zeno regimes. Also, increasing the coupling strength does not generally increase the effective decay rate $\Gamma_{n}(\tau)$. This behavior should be contrasted with the weak coupling scenario. For weak coupling, it has been found that the effective decay rate is still given by Eq. (8), but now the filter function is $N_{S}$ times the filter function given by Eq. (9)[34]. Thus, increasing the coupling strength should now increase the effective decay rate. This is precisely what is observed in Fig.5(b). Consequently, the weak coupling and the strong coupling regimes are very different for the strong and the weak coupling regimes. The difference is again due to the system-environment correlations.

Discussion

We have investigated the quantum Zeno and anti-Zeno effects for a single two-level system interacting strongly with an environment of harmonic oscillators. Although it seems that perturbation theory cannot be applied, we have applied a polaron transformation that can make the coupling strength effectively small in the transformed frame and thereby validate the use of perturbation theory. We have obtained general expressions for the effective decay rate, independent of any particular form of the spectral density of the environment. Thereafter, we have shown that the strong coupling regime shows both qualitative and quantitative differences in the behavior of the effective decay rate as a function of the measurement interval and the QZE to QAZE transitions as compared with the weak system-environment coupling scenario. The effective decay rate is no longer an overlap integral of the spectral density of the environment and some other function. Rather, there is a very pronounced non-linear dependence on the spectral density of the environment. Most importantly, increasing the coupling strength in the strong coupling regime can actually reduce the effective decay rate. These differences can be understood in terms of the significant role played by the system-environment correlations. Moreover, we have extended our results to many two-level systems interacting with a common environment. Once again, we obtained expressions for the effective decay rate that are independent of the spectral density of the environment. We illustrated that in this case as well the behavior of the effective decay rate is very different from the commonly considered weak coupling regime. Our results should be important for understanding the role of repeated measurements in quantum systems that are interacting strongly with their environment.

Methods

The polaron transformation

For completeness, let us sketch how to transform the spin-boson Hamiltonian to the polaron frame [36, 37, 38, 39, 40, 41]. We need to find $H=e^{\chi\sigma_{z}/2}H_{L}e^{-\chi\sigma_{z}/2}$. We use the identity

[TABLE]

Now, it is clear that $[\chi\sigma_{z}/2,\sigma_{z}/2]=0$. Also, $\left[\sum_{k}\left(\frac{g_{k}}{\omega_{k}}b_{k}^{\dagger}-\frac{g_{k}^{*}}{\omega_{k}}b_{k}\right),\sum_{k}\omega_{k}b_{k}^{\dagger}b_{k}\right]=-\sigma_{z}\sum_{k}(g_{k}b_{k}^{\dagger}+g_{k}^{*}b_{k}).$ Carrying on, we find that $\left[\sigma_{z}\chi/2,\sigma_{z}\sum_{k}(g_{k}^{*}b_{k}+g_{k}b_{k}^{\dagger})\right]=-2\sum_{k}\frac{|g_{k}|^{2}}{\omega_{k}}.$ This is simply a c-number, so the higher-order commutators are zero. Furthermore, this c-number leads to a constant shift in the transformed Hamiltonian, and can thus be dropped. Putting all the commutators together, we find that

[TABLE]

Next, we observe that $e^{\chi\sigma_{z}/2}\frac{\Delta}{2}\sigma_{x}e^{-\chi\sigma_{z}/2}=e^{\chi\sigma_{z}/2}\frac{\Delta}{2}(\sigma_{+}+\sigma_{-})e^{-\chi\sigma_{z}/2}$, where $\sigma_{+}$ and $\sigma_{-}$ are the standard spin half raising and lowering operators. Furthermore, $[\chi\sigma_{z}/2,\sigma_{+}]=\sigma_{+}\chi$, leading to $e^{\chi\sigma_{z}/2}\sigma_{+}e^{-\chi\sigma_{z}/2}=\sigma_{+}e^{\chi}$. Similarly, $e^{\chi\sigma_{z}/2}\sigma_{-}e^{-\chi\sigma_{z}/2}=\sigma_{-}e^{-\chi}$. Thus, we finally have the required Hamiltonian in the polaron frame.

For the large spin case, the calculation is very similar [37]. The major difference is that now the c-number term that we dropped before cannot be dropped any longer since this term is proportional to $J_{z}^{2}$ (for the spin half case, this is proportional to the identity operator, so this is just a constant shift for the spin half case). Namely, we now find that

[TABLE]

Also,

[TABLE]

The rest of the calculation is very similar to the spin half case, and leads to the Hamiltonian in the polaron frame.

Finding the system density matrix in the polaron frame

Here we describe how to obtain the system density matrix in the polaron frame $\rho_{S}(\tau)$ just before performing the measurement at time $\tau$. We define $U_{\text{tot}}(\tau)=e^{-iH\tau}=U_{0}(\tau)U_{I}(\tau)$, where $U_{0}(\tau)$ is the unitary time-evolution operator corresponding to $H_{S}$ and $H_{B}$, while $U_{I}(\tau)$ is the ‘left over’ part that we can find using time-dependent perturbation theory. Writing the system-environment coupling in the polaron frame as $\sum_{\mu}F_{\mu}\otimes B_{\mu}$, with $F_{1}=\frac{\Delta}{2}\sigma_{+}$, $B_{1}=X$, $F_{2}=\frac{\Delta}{2}\sigma_{-}$ and $B_{2}=X^{\dagger}$, $U_{I}(\tau)$ can be found to be $U_{I}(\tau)=\mathds{1}+A_{1}+A_{2}$, with $A_{1}=-i\sum_{\mu}\int_{0}^{\tau}\widetilde{F}_{\mu}(t_{1})\widetilde{B}_{\mu}(t_{1})dt_{1}$ and $A_{2}=-\sum_{\mu\nu}\int_{0}^{\tau}dt_{1}\int_{0}^{t_{1}}dt_{2}\widetilde{F}_{\mu}(t_{1})\widetilde{F}_{\nu}(t_{2})\widetilde{B}_{\mu}(t_{1})\widetilde{B}_{\nu}(t_{2})$. Correct to second order in the tunneling amplitude $\Delta$, we can then write

[TABLE]

where $\rho_{\text{tot}}(0)=\rho_{S}(0)\otimes\rho_{B}$. Eq. (12) can now be simplified term by term. First, we find that $\text{Tr}_{B}\{U_{0}(\tau)\rho_{\text{tot}}(0)U_{0}^{\dagger}(\tau)\}=\widetilde{\rho}_{S}(\tau)$, where $\widetilde{\rho}_{S}(\tau)=U_{S}(\tau)\rho_{S}(0)U_{S}^{\dagger}(\tau)$ is the system density matrix if the tunneling amplitude is zero. Next, we find that $\text{Tr}_{B}\{U_{0}(\tau)\rho_{\text{tot}}(0)A_{1}^{\dagger}U_{0}^{\dagger}(\tau)\}=i\sum_{\mu}\int_{0}^{\tau}dt_{1}\,U_{S}(\tau)\rho_{S}(0)\widetilde{F}_{\mu}(t_{1})U_{S}^{\dagger}(\tau)\langle\widetilde{B}_{\mu}(t_{1})\rangle_{B}$, where $\langle\widetilde{B}_{\mu}(t_{1})\rangle_{B}=\text{Tr}_{B}\{U_{B}(\tau)\rho_{B}\widetilde{B}_{\mu}(t_{1})U_{B}^{\dagger}(\tau)\}$. Similarly, $\text{Tr}_{B}\{U_{0}(\tau)A_{1}\rho_{\text{tot}}(0)U_{0}^{\dagger}(\tau)\}=-i\sum_{\mu}\int_{0}^{\tau}dt_{1}\,U_{S}(\tau)\widetilde{F}_{\mu}(t_{1})\rho_{S}(0)U_{S}^{\dagger}(\tau)\langle\widetilde{B}_{\mu}(t_{1})\rangle_{B}$. Carrying on,

[TABLE]

with the environment correlation function $C_{\mu\nu}(t_{1},t_{2})$ defined as $C_{\mu\nu}(t_{1},t_{2})=\langle\widetilde{B}_{\mu}(t_{1})\widetilde{B}_{\nu}(t_{2})\rangle_{B}=\text{Tr}_{B}\{\widetilde{B}_{\mu}(t_{1})\widetilde{B}_{\nu}(t_{2})\rho_{B}\}$. Similarly,

[TABLE]

Finally,

[TABLE]

Using the fact that $\int_{0}^{\tau}dt_{1}\int_{0}^{\tau}dt_{2}=\int_{0}^{\tau}dt_{1}\int_{0}^{t_{1}}dt_{2}+\int_{0}^{\tau}dt_{2}\int_{0}^{t_{2}}dt_{1}$,

[TABLE]

where h.c. denotes hermitian conjugate. Putting all the terms back together, the system density matrix can be written as Eq. (2).

Finding the effective decay rate

We now explain how to find the effective decay rate given by Eq. (3). With the system density matrix at time $\tau$ available, we first calculate the survival probability $s(\tau)$. This can be done via $s(\tau)=1-\langle\downarrow|\rho_{S}(\tau)|\downarrow\rangle$. Since the state $\left|\downarrow\right\rangle$ is an eigenstate of $H_{S}$, and $\rho_{S}(0)=\left|\uparrow\right\rangle\left\langle\uparrow\right|$, it is straightforward to see that

[TABLE]

We now note that, since $F_{1}=\frac{\Delta}{2}\sigma_{+}$, $F_{2}=\frac{\Delta}{2}\sigma_{-}$, and $H_{S}=\frac{\varepsilon}{2}\sigma_{z}$, $\widetilde{F}_{1}(t)=\frac{\Delta}{2}\sigma_{+}e^{i\varepsilon t}$ and $\widetilde{F}_{2}(t)=\frac{\Delta}{2}\sigma_{-}e^{-i\varepsilon t}$. Therefore,

[TABLE]

What remains to be worked out is the environment correlation function $C_{12}(t_{1},t_{2})=\text{Tr}_{B}[\rho_{B}\widetilde{X}(t_{1})\widetilde{X}^{\dagger}(t_{2})]$. Using the cyclic invariance of the trace, it is clear that this correlation function is actually only a function of $t_{1}-t_{2}$ only, since $\text{Tr}_{B}[\rho_{B}\widetilde{X}(t_{1})\widetilde{X}^{\dagger}(t_{2})]=\text{Tr}_{B}[\rho_{B}e^{iH_{B}(t_{1}-t_{2})}Xe^{-iH_{B}(t_{1}-t_{2})}X^{\dagger}]=C_{12}(t_{1}-t_{2})$. Thus,

[TABLE]

where we have introduced $t^{\prime}=t_{1}-t_{2}$. The calculation of $C_{12}(t^{\prime})$ can be performed as follows. First, we use the useful fact that for $\text{Tr}_{B}[\rho_{B}e^{Z}]=e^{\langle Z^{2}\rangle/2}$ where $Z$ is a linear function of the creation and annihilation operators. Second, to obtain a single exponential so that the previous identity can be used, we use the identity that for any two operators $X$ and $Y$, $e^{X}e^{Y}=e^{X+Y+\frac{1}{2}[X,Y]\ldots}$. Fortunately for us, the series terminates for our case, so the higher order terms are zero. Using these two identities, we find that $C_{12}(t^{\prime})=e^{-\Phi_{R}(t^{\prime})}e^{-i\Phi_{I}(t^{\prime})}$ where $\Phi_{R}(t)$ and $\Phi_{I}(t)$ have been defined in Eq. (4), and the spectral density of the environment has been introduced as $\sum_{k}|g_{k}|^{2}(\ldots)\rightarrow\int_{0}^{\infty}d\omega J(\omega)(\ldots)$. This finally leads to $s(\tau)=1-\frac{\Delta^{2}}{2}\int_{0}^{\tau}dt\int_{0}^{t}dt^{\prime}e^{-\Phi_{R}(t^{\prime})}\cos[\varepsilon t^{\prime}-\Phi_{I}(t^{\prime})]$. We can define an effective decay rate $\Gamma(\tau)=-\frac{1}{\tau}\ln s(\tau)$. For small $\Delta$, we expect the deviation of the survival probability from one to be small. Thus, we end up with Eq. (3).

Calculating the modified decay rate

Let us now briefly sketch how to obtain Eq. (7). We note that the system Hamiltonian $H_{S,L}$ becomes in the polaron frame $H_{S,P}=\frac{\varepsilon}{2}\sigma_{z}+\frac{\Delta}{2}(\sigma_{+}X+\sigma_{-}X^{\dagger})$. Then, to second order in $\Delta$, $e^{-iH_{S,P}\tau}\approx\mathds{1}+A_{SP}^{(1)}+A_{SP}^{(2)}$, with $A_{SP}^{(1)}=-i\int_{0}^{\tau}dt_{1}\sum_{\mu}\widetilde{F}_{\mu}(t_{1})B_{\mu}$, and $A_{SP}^{(2)}=-\int_{0}^{\tau}dt_{1}\int_{0}^{\tau}dt_{2}\sum_{\mu\nu}\widetilde{F}_{\mu}(t_{1})B_{\mu}\widetilde{F}_{\nu}(t_{2})B_{\nu}$. Here, $\widetilde{F}_{\mu}(t)=e^{i\varepsilon\sigma_{z}t/2}F_{\mu}e^{-i\varepsilon\sigma_{z}t/2}$. Substituting these expressions in the expression for the survival probability as well as the perturbation expansions for $e^{iH\tau}$ and $e^{-iH\tau}$, and keeping terms up to second order in $\Delta$, we find that the new survival probability consists of the previous survival probability plus some additional terms. It can be be easily seen that most of these additional terms, once the trace with the projector $\left|\downarrow\right\rangle\left\langle\downarrow\right|$ is taken, give zero. The additional terms that need to be worked out are $\text{Tr}_{B}[\langle\downarrow|U_{B}(\tau)A_{1}\rho_{S}(0)\rho_{B}U_{B}^{\dagger}(\tau)A_{SP}^{(1)}|\downarrow\rangle]$, $\text{Tr}_{B}[\langle\downarrow|A_{SP}^{(1)\dagger}U_{B}(\tau)\rho_{S}(0)\rho_{B}A_{1}^{\dagger}U_{B}^{\dagger}(\tau)|\downarrow\rangle]$, and $\text{Tr}_{B}[\langle\downarrow|A_{SP}^{(1)\dagger}U_{B}(\tau)\rho_{S}(0)\rho_{B}U_{B}^{\dagger}A_{SP}^{(1)}|\downarrow\rangle]$. The first of these terms is equal to

[TABLE]

while the second is simply the hermitian conjugate of the first. On the other hand, the last term is equal to

[TABLE]

Next, we use the fact that $F_{1}=\frac{\Delta}{2}\sigma_{+}$ and $F_{2}=\frac{\Delta}{2}\sigma_{-}$ to simply the inner products. Putting all the pieces together, we arrive at Eq. (7). The calculation of $\Gamma_{n}(\tau)$ for the large spin case is quite similar. One only needs to be careful about the fact that the system-environment Hamiltonian, in the polaron frame, contains a term, namely $-\kappa J_{z}^{2}$, that is not a part of the transformed system Hamiltonian $H_{S,P}$.