Protecting quantum resources via frequency modulation of qubits in leaky cavities

Abstract

Finding strategies to preserve quantum resources in open systems is nowadays a main requirement for reliable quantum-enhanced technologies. We address this issue by considering structured cavities embedding qubits driven by a control technique known as frequency modulation. We first study a single qubit in a lossy cavity to determine optimal modulation parameters and qubit-cavity coupling regime allowing a gain of four orders of magnitude concerning coherence lifetimes. We relate this behavior to the inhibition of the qubit effective decay rate rather than to stronger memory effects (non-Markovianity) of the system. We then exploit these findings in a system of noninteracting qubits embedded in separated cavities to gain basic information about scalability of the procedure. We show that the determined modulation parameters enable lifetimes of quantum resources, such as entanglement, discord and coherence, three orders of magnitude longer than their natural (uncontrolled) decay times. We discuss the feasibility of the system within the circuit-QED scenario, typically employed in the current quantum computer prototypes. These results provide new insights towards efficient experimental strategies against decoherence.

Introduction

Quantum coherence, stemming from the superposition principle, is the key concept that distinguishes quantum mechanics from classical mechanics. It is well-known that quantum coherence among different bodies of a multipartite systems is the basic ingredient for the formation of quantum correlations^{1,2,3,4,5,6,7,8,9,10}. Among these, entanglement represents the part due to non-separability of states, utilized in many paradigmatic quantum information processes such as teleportation^11,12,13, quantum error correction^14,15, quantum key distribution¹⁶ and quantum dense coding¹⁷. Useful quantum correlations beyond entanglement, occurring even for separable states, have been identified by the so-called discord^18,19, which is in turn employed for specific quantum information tasks^20,21,22. As a very fundamental trait, quantum coherence itself can be quantified^{23,24,25,26,27,28} and exploited as a resource in various quantum protocols^{23,28,29,30,31}. One of the cutting-edge fields of application of these quantum features is, for instance, quantum metrology, aiming at achieving high precision measurements characterized by the quantum Fisher information^{32,33,34,35,36,37,38,39,40,41,42,43}, with recent developments in open quantum systems^{39,41,44,45,46,47}.

Open quantum systems undergo decoherence due to the interaction with the surrounding environment, which usually destroys coherence and correlations thus limiting their practical use in quantum information⁴⁸ and metrology^{36,45,49,50,51}. The dynamics of open quantum systems can be classified in Markovian (memoryless) and non-Markovian (memory-keeping) regimes^52,53. Markovian regime is typically associated to the lack of memory effects, a situation where the information irreversibly flows from the system to the environment. Differently, non-Markovian regime implies that the past history of the system affects the present one: such a memory effect results in an information feedback from the environment to the system^{48,52,53,54,55,56,57,58,59,60,61,62}. In this sense, non-Markovianity has been quantified by a variety of measures^{61,63,64,65,66,67,68,69,70,71,72} and regarded as a precious resource for certain applications^61,70,73. Although memory effects permit coherence, entanglement and discord of noninteracting subsystems to partially revive after disappearing during the time evolution^{48,53,65,74,75,76,77}, these revivals eventually decay. Decoherence remains one of the main drawbacks to overcome towards the implementation of quantum-enhanced technologies, which require qubits with long-lived quantum features. To this purpose, many strategies have been proposed in order to harnessing and protecting quantum resources in systems of qubits under different environmental conditions^{78,79,80,81,82,83,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98,99,100,101,102,103,104,105,106,107,108,109,110,111,112,113,114,115,116,117,118,119}.

In this paper, we address this issue by using qubit frequency modulation as control technique. Generally, a quantum system is frequency modulated when its energy levels are shifted by an external driving. Frequency modulation in an atomic qubit can be performed by applying an external off-resonant field^120,121,122. Moreover, the most recent experimental progress in fabrication and control of quantum circuit-QED devices enables frequency modulation in superconducting Josephson qubits (artificial atoms)^{122,123,124,125,126}, which are the preferred building blocks in current quantum computer prototypes¹²⁷. It has been reported that external control of the qubit frequency can induce sidebands transitions^128,129, modify its fluorescence spectrum¹³⁰ and population dynamics^{131,132,133,134,135}, as well as amplify non-Markovianity¹³⁶. Motivated by these considerations, here we aim at providing a thorough quantitative analysis of the role frequency modulation can play in prolonging lifetimes of desired quantum resources. Differently from previous studies, we mainly focus on determining the values of parameters of the driving field and of the qubit-cavity coupling which can efficiently maintain coherence and correlations for times orders of magnitude longer than their natural lifetimes in absence of external control. We first consider a single frequency-modulated qubit inside a structured leaky cavity, which allows us to determine the desired values of parameters. We also analyze how the choice of these parameters plays a role in preserving quantum Fisher information so to enhance the precision of phase estimation in the system during the evolution. The reliability of the theoretical results is checked by taking into account experimental parameters typical of the circuit-QED scenario for a transmon qubit in a coplanar resonator. We then extend the analysis to a system of two noninteracting qubits in separated cavities for establishing at which extent entanglement, discord and coherence can be shielded from decay thanks to individual control of the qubit frequencies. This novel investigation of simple systems constitutes the essential step for acquiring useful insights which can be straightforwardly generalized to a system of many separated qubits for scalability.

The paper is organized as follows. In Sec. 1 we describe the evolution of the frequency-modulated open qubit, giving the dynamics of coherence and quantum Fisher information (QFI). In Sec. 2 we discuss the dynamics of quantum resources in the two-qubit system for some classes of initially correlated states. Finally, in Sec. 2.3 we summarize the main conclusions.

Results

Single-qubit system

The system consists of a qubit (two-level system) coupled to a zero-temperature reservoir formed by the quantized modes of a high-Q cavity. The excited and ground states of the qubit are labeled by $|e\rangle$ and $|g\rangle$, respectively. It is assumed that the transition frequency of the qubit ω₀ is modulated sinusoidally by an external driving field, as depicted in Fig. 1. The Hamiltonian of the system is

$$\hat{H}={\hat{H}}_{{\rm{q}}}+{\hat{H}}_{{\rm{r}}}+{\hat{H}}_{{\rm{in}}},$$

(1)

where ${\hat{H}}_{{\rm{q}}}$ is the qubit Hamiltonian (ħ = 1)

$${\hat{H}}_{{\rm{q}}}=\frac{1}{2}[{\omega }_{0}+\delta \,\cos ({\rm{\Omega }}t)]{\hat{\sigma }}_{z},$$

(2)

with δ and Ω denoting the modulation amplitude and frequency, respectively, and ${\hat{\sigma }}_{z}$ the Pauli operator along the z-direction. ${\hat{H}}_{{\rm{r}}}={\sum }_{k}\,{\omega }_{k}{\hat{a}}_{k}^{\dagger }{\hat{a}}_{k}$ is the reservoir Hamiltonian with ${\hat{a}}_{k}^{\dagger }$ (${\hat{a}}_{k}$) being the creation (annihilation) operator of the cavity field mode k. ${\hat{H}}_{{\rm{in}}}$ describes the interaction between the qubit and the cavity modes which, in the dipole and rotating-wave approximation, can be written as

$${\hat{H}}_{{\rm{in}}}=\sum _{k}\,({g}_{k}{\hat{\sigma }}_{+}{\hat{a}}_{k}+{g}_{k}^{\ast }{\hat{a}}_{k}^{\dagger }{\hat{\sigma }}_{-}),$$

(3)

where g_k is the coupling constant between the qubit and mode k while ${\hat{\sigma }}_{+}$ (${\hat{\sigma }}_{-}$) represents the raising (lowering) operator for the qubit. We point out that, in general, when a nonstationary mode (in this case, the qubit mode) is coupled to the environment, the system-reservoir couplings g_k also become time dependent^137,138. However, here we focus on an experimentally feasible system where $\delta /{\omega }_{0}\ll 1$ and ${\rm{\Omega }}/{\omega }_{0}\ll 1$ (see subsection 1.5), which guarantee the so-called adiabatic regime: under such condition, the qubit-reservoir interaction is still modelled by the Jaynes-Cummings interaction, with the couplings g_k independent of time, despite the nonstationary qubit mode^138,139. Moreover, the qubit dimension is taken to be much smaller than the wavelength of the emitted photon so that the dipole approximation remains valid and, since ω₀ ± δ ≈ ω₀, also the rotating wave approximation holds along our analysis^132,136.

Figure 1

Sketch of the single-qubit system. A qubit with transition frequency ω₀ is embedded in a high-Q leaky cavity, whose photon losses are characterized by the spectral width λ ∝ 1/Q of the coupling, assuming a Lorentzian spectrum for the cavity modes. The transition frequency of the qubit is modulated sinusoidally by an external driving field with a modulation amplitude δ and a modulation frequency Ω.

Moving to a non-uniformly rotating frame (interaction picture) by the unitary transformation

$$\hat{U}=\exp \,[\,-\,i\{\sum _{k}\,{\omega }_{k}{\hat{a}}_{k}^{\dagger }{\hat{a}}_{k}+[{\omega }_{0}+(\delta /{\rm{\Omega }})\,\sin \,{\rm{\Omega }}t]{\hat{\sigma }}_{z}\}],$$

(4)

the effective Hamiltonian ${\hat{H}}_{{\rm{eff}}}={\hat{U}}^{\dagger }\hat{H}\hat{U}+(\partial {\hat{U}}^{\dagger }/\partial t)\hat{U}$ is

$$\begin{array}{rcl}{\hat{H}}_{{\rm{eff}}} & = & \sum _{k}\,{g}_{k}{\hat{\sigma }}_{+}{\hat{a}}_{k}{e}^{-i({\omega }_{k}-{\omega }_{0})t}{e}^{i(\delta /{\rm{\Omega }})\sin {\rm{\Omega }}t}\\ & & +\,\sum _{k}\,{g}_{k}^{\ast }{\hat{a}}_{k}^{\dagger }{\hat{\sigma }}_{-}{e}^{i({\omega }_{k}-{\omega }_{0})t}{e}^{-i(\delta /{\rm{\Omega }})\sin {\rm{\Omega }}t}.\end{array}$$

(5)

We notice that ${\hat{H}}_{{\rm{eff}}}$ shows how the system behaves as if the same frequency modulation was applied to each qubit-cavity coupling constant g_k. This is relevant to our purpose since, under this condition, it is known that decoherence can be softened^135,140. Making use of Jacobi-Anger expansion, the exponential factors in Eq. (5) can be written as

$${e}^{\pm i(\delta /{\rm{\Omega }})\sin ({\rm{\Omega }}t)}={J}_{0}(\frac{\delta }{{\rm{\Omega }}})+2\,\sum _{n=1}^{\infty }\,{(\pm i)}^{n}{J}_{n}(\frac{\delta }{{\rm{\Omega }}})\,\cos (n{\rm{\Omega }}t),$$

(6)

where ${J}_{n}(\frac{\delta }{{\rm{\Omega }}})$ is the n-th Bessel function of the first kind.

We assume the qubit is initially in a coherent superposition $\alpha |e\rangle +\beta |g\rangle$ and the reservoir modes in the vacuum state $|0\rangle$, so that the overall initial state is

$$|{\rm{\Psi }}\mathrm{(0)}\rangle =(\alpha |e\rangle +\beta |g\rangle )|0\rangle \mathrm{.}$$

(7)

Hence, at any later time t the quantum state of the whole system can be written as

$$|{\rm{\Psi }}(t)\rangle =\alpha {C}_{e}(t)|e\rangle |0\rangle +\beta |g\rangle |0\rangle +\sum _{k}\,{C}_{g,k}(t)|g\rangle |{1}_{k}\rangle ,$$

(8)

where $|{1}_{k}\rangle$ is the cavity state with a single photon in mode k and C_g,k(t) is its probability amplitude. Using the time-dependent Schrödinger equation, the differential equations for the probability amplitudes C_e(t) and C_g,k(t) are, respectively,

$${\dot{C}}_{e}(t)=-\,i{e}^{i(\delta /{\rm{\Omega }})\sin ({\rm{\Omega }}t)}\,\sum _{k}\,{g}_{k}{e}^{-i({\omega }_{k}-{\omega }_{0})t}{C}_{g,k}(t),$$

(9)

and

$${\dot{C}}_{g,k}(t)=-\,i{e}^{-i(\delta /{\rm{\Omega }})\sin ({\rm{\Omega }}t)}{g}_{k}^{\ast }{e}^{i({\omega }_{k}-{\omega }_{0})t}{C}_{e}(t\mathrm{).}$$

(10)

Solving Eq. (10) formally and substituting the solution into Eq. (9), one obtains

$${\dot{C}}_{e}(t)+{\int }_{0}^{t}\,dt^{\prime} G(t,t^{\prime} ){C}_{e}(t^{\prime} )=0,$$

(11)

where the kernel G(t, t′), that is the correlation function including the memory effects, is

$$G(t,\,t^{\prime} )={e}^{i(\delta /{\rm{\Omega }})[\sin ({\rm{\Omega }}t)-\sin ({\rm{\Omega }}t^{\prime} )]}\,\sum _{k}\,{|{g}_{k}|}^{2}{e}^{-i({\omega }_{k}-{\omega }_{0})(t-t^{\prime} )}\mathrm{.}$$

(12)

In the continuous limit, the kernel above becomes

$$G(t,\,t^{\prime} )={e}^{i(\delta /{\rm{\Omega }})[\sin ({\rm{\Omega }}t)-\sin ({\rm{\Omega }}t^{\prime} )]}\,{\int }_{0}^{\infty }\,J(\omega ){e}^{-i(\omega -{\omega }_{0})(t-t^{\prime} )}d\omega ,$$

(13)

where $J(\omega )={\sum }_{k}\,|{g}_{k}{|}^{2}\delta (\omega -{\omega }_{k})$ is, as usual, the spectral density of the reservoir (cavity) modes^52,136. We choose a Lorentzian spectral density, which is typical of a structured cavity⁵², whose form is

$$J(\omega )=\frac{1}{2\pi }\frac{\gamma {\lambda }^{2}}{{({\omega }_{0}-\omega )}^{2}+{\lambda }^{2}},$$

(14)

where λ indicates the spectral width of the coupling and is related to the reservoir correlation time τ_r via τ_r = λ^−1. On the other hand, γ represents the decay rate of the excited state of the qubit in the Markovian limit of flat spectrum (i.e., the spontaneous emission decay rate) and it is linked to the qubit relaxation time τ_q by τ_q = γ^−1 ⁵². Qubit-cavity weak coupling occurs for λ > γ (τ_r < τ_q); the opposite condition λ < γ (τ_r > τ_q) thus identifies strong coupling. The larger the cavity quality factor, the smaller the spectral width λ.

With such a spectral density, the kernel of Eq. (13) becomes

$$G(t,\,t^{\prime} )=\frac{\gamma \lambda }{2}{e}^{-\lambda (t-t^{\prime} )}{e}^{i(\delta /{\rm{\Omega }})[\sin ({\rm{\Omega }}t)-\sin ({\rm{\Omega }}t^{\prime} )]}\mathrm{.}$$

(15)

Substituting it into Eq. (11), one gets

$${\dot{C}}_{e}(t)+\frac{\gamma \lambda }{2}{e}^{i(\delta /{\rm{\Omega }})\sin ({\rm{\Omega }}t)}\,{\int }_{0}^{t}\,dt^{\prime} {e}^{-i(\delta /{\rm{\Omega }})\sin ({\rm{\Omega }}t^{\prime} )}{e}^{-\lambda (t-t^{\prime} )}{C}_{e}(t^{\prime} )=0.$$

(16)

Calculating C_e(t) from this equation, the reduced density matrix of the qubit ρ_q(t) in the basis $\{|e\rangle ,|g\rangle \}$ is given by

$${\rho }_{{\rm{q}}}(t)=(\begin{array}{cc}{|\alpha |}^{2}{|{C}_{e}(t)|}^{2} & \alpha {\beta }^{\ast }{C}_{e}(t)\\ {\alpha }^{\ast }\beta {C}_{e}^{\ast }(t) & 1-{|\alpha |}^{2}{|{C}_{e}(t)|}^{2}\end{array}).$$

(17)

Taking the derivative of Eq. (16) with respect to time, we have

$$\begin{array}{rcl}{\dot{\rho }}_{q}(t) & = & -\,i\frac{{\rm{\Omega }}(t)}{2}[{\sigma }_{+}{\sigma }_{-},{\rho }_{q}(t)]+\frac{{\rm{\Gamma }}(t)}{2}(2{\sigma }_{-}{\rho }_{q}(t){\sigma }_{+}\\ & & -\,{\sigma }_{+}{\sigma }_{-}{\rho }_{q}(t)-{\rho }_{q}(t){\sigma }_{+}{\sigma }_{-}),\end{array}$$

(18)

where ${\rm{\Omega }}(t)=-\,2{\rm{Im}}\,[\frac{{\dot{C}}_{e}(t)}{{C}_{e}(t)}]$ plays the role of a time-dependent Lamb shift and ${\rm{\Gamma }}(t)=-\,2{\rm{Re}}\,[\frac{{\dot{C}}_{e}(t)}{{C}_{e}(t)}]$ can be interpreted as a time-dependent decay rate⁵².

The evolved density matrix of the qubit ρ_q(t) shall be utilized for obtaining the dynamics of the quantum properties of interest, such as coherence, quantum Fisher information and non-Markovianity. Before displaying the dynamics of these quantities, we recall their definitions and give their time-dependent expressions in the following.

Coherence

Coherence of a quantum state characterizes the property of superposition among the basis states of the system and can be regarded as a basis-dependent resource by itself ²⁸. Among the various bona-fide measures of coherence, we choose here the so-called l₁-norm defined as²⁹

$$\zeta =\mathop{{\rm{\min }}}\limits_{{\rho }_{{\rm{inc}}}\in {\mathcal I} }\,\parallel \rho -{\rho }_{{\rm{inc}}}{\parallel }_{{l}_{1}}=\sum _{i\ne j}\,|{\rho }_{ij}|,$$

(19)

where ρ is the density matrix of an arbitrary quantum state and the minimum is taken over the set ${\mathcal I}$ of incoherent states ρ_inc. Such a measure, depending only on the off-diagonal elements ρ_ij (i ≠ j) of the quantum state, is clearly related to the fundamental property of quantum interference.

For ρ_q(t) of Eq. (18) and assuming $\alpha =\beta =1/\sqrt{2}$ in the initial state of Eq. (7), we have ζ(0) = 1 (maximum initial coherence) and a time-dependent qubit coherence ζ(t) = |C_e(t)|.

Quantum Fisher information

Quantum metrology exploits quantum-mechanical effects to reach high precision measurements. In a typical metrological procedure, one first encodes the parameter of interest ϕ on a probe state ρ_in by means of a unitary process U_ϕ. Thus, the output state is ${\rho }_{\varphi }={U}_{\varphi }{\rho }_{in}{U}_{\varphi }^{\dagger }$. The output state ρ_ϕ is then measured by a set of positive operator valued measurements and the value of ϕ finally estimated from the outcomes. It is known that the precision in estimating ϕ is limited by the quantum Cramer-Rao bound inequality^141,142

$$\delta \varphi \ge 1/\sqrt{{F}_{\varphi }},$$

(20)

where δϕ is the standard deviation associated to the variable ϕ and F_ϕ is the quantum Fisher information (QFI) defined as F_ϕ = Tr(ρ_ϕL²), with L being the so-called symmetric logarithmic derivative determined by ∂_ϕρ_ϕ = (Lρ(ϕ) + ρ(ϕ)L)/2 with ∂_ϕ = ∂/∂ϕ^141,142. From the spectral decomposition of the ϕ-dependent density matrix ${\rho }_{\varphi }={\sum }_{m}\,{p}_{m}|{\psi }_{m}\rangle \langle {\psi }_{m}|$, where p_m and |ψ_m〉 are, respectively, eigenvalues and eigenstates, the QFI is known to have the analytical expression³⁷

$${F}_{\varphi }=\sum _{m,n}\,\frac{2}{{p}_{m}+{p}_{n}}|\langle {\psi }_{m}|{\partial }_{\varphi }{\rho }_{\varphi }|{\psi }_{n}\rangle {|}^{2}.$$

(21)

For our dynamical system, we consider a phase-estimation problem where ${U}_{\varphi }\equiv |g\rangle \langle g|+{e}^{i\varphi }|e\rangle \langle e|$ acts on the initial maximally coherent state of the qubit $|{\psi }_{+}\rangle =(|e\rangle +|g\rangle )/\sqrt{2}$³² and successively let the system evolve under the dissipative noise and frequency modulation. The initial overall qubit-cavity state is therefore (U_ϕ|ψ₊〉)|0〉 and the evolved reduced density matrix of the qubit ρ_q,ϕ(t) has the form of Eq. (17) with a ϕ-dependence in the off-diagonal elements. Using Eq. (21) we get F_ϕ(t) = |C_e(t)|², that is the coherence squared. The associated minimum estimation error is therefore $\delta {\varphi }_{{\rm{\min }}}(t)=\mathrm{1/}\sqrt{{F}_{\varphi }(t)}=\mathrm{1/|}{C}_{e}(t)|$.

Non-Markovianity

To discuss the non-Markovian character of our system we employ, among the various quantifiers, the well-known measure based on the dynamics of the trace distance between two initially different states ρ₁(0) and ρ₂(0) of the qubit, identifying information backflows. It is defined as⁶⁴

$${\mathscr{N}}=\mathop{{\rm{\max }}}\limits_{{\rho }_{1}(0),{\rho }_{2}(0)}\,{\int }_{\sigma > 0}\,\sigma [t,{\rho }_{1}(0),{\rho }_{2}(0)]\,{\rm{d}}t,$$

(22)

where $\sigma [t,\,{\rho }_{1}\mathrm{(0)},\,{\rho }_{2}\mathrm{(0)]}={\rm{d}}D[{\rho }_{1}(t),\,{\rho }_{2}(t)]/{\rm{d}}t$, with $D[{\rho }_{1}(t),\,{\rho }_{2}(t)]=\mathrm{(1/2)}\mathrm{Tr}|{\rho }_{1}(t)-{\rho }_{2}(t)|$ being the trace distance $(|X|=\sqrt{{X}^{\dagger }X})$. In Eq. (22) the integration is taken over time intervals when σ > 0 and the maximization is made over all the possible pairs of initial states ρ₁(0) and ρ₂(0). It is noteworthy that the trace distance is related to the distinguishability between quantum states, whereas its time derivative (σ) means a flow of information between the system and its environment. While Markovian processes satisfy σ ≤ 0 for all pairs of initial states ρ_1,2(0) at any time t, non-Markovian ones admit at least a pair of initial states such that σ > 0 for some time intervals, so that the information flows from the environment back to the system⁶⁴.

For our system, which experiences a dissipative dynamics describable as an amplitude damping channel, the evolved qubit density matrix has the form of Eq. (17). Under this condition, it is known that the non-Markovianity measure ${\mathscr{N}}$ is maximized for the choice of the initial orthogonal states $(|e\rangle \pm |g\rangle )/\sqrt{2}$ of the qubit⁷¹ and assumes the expression

$${\mathscr{N}}=-\,(1/2)\,{\int }_{{\rm{\Gamma }} < 0}\,{\rm{\Gamma }}(t)|{C}_{e}(t)|\,{\rm{d}}t,$$

(23)

where Γ(t) is the effective time-dependent decay rate defined after Eq. (18).

Time evolution of single-qubit resources

We start the quantitative analysis by studying the time evolution of coherence ζ(t) under a weak coupling regime (λ = 3γ), plotted in Fig. 2 for different values of the modulation frequency Ω and fixed modulation amplitude δ = 10γ. This evolution is compared to the case when the external driving is off (δ = 0, Ω = 0), for which coherence disappears at time t* ~ 10/γ. As shown, under the given condition for δ, coherence survives for times longer than t* when Ω is smaller. Presence of oscillations in the dynamics (see Fig. 2(b)) indicates the manifestation of non-Markovian effects under weak coupling regime¹³⁶. However, frequency modulation can even produce negative effects for large values of the modulation frequency, Ω ≥ 100γ, which are such that the coherence disappears at times t < t*.

Figure 2

Qubit coherence ζ(t) as a function of scaled time γt for different values of the modulation frequency Ω. The values of other parameters are: δ = 10γ, λ = 3γ (weak coupling), $\alpha =\beta =1/\sqrt{2}$. Solid-blue line in panel (a) corresponds to the situation in which frequency modulation is off.

This fact suggests to see what happens under a strong coupling regime. Figure 3 shows the time behavior of coherence under such a condition (λ = 0.01γ) for different values of modulation frequency at a fixed modulation amplitude δ = 5γ. Interesting results are obtained for Ω = 0.001γ, Ω = 0.9γ and Ω = 2.1γ. The coherence dynamic for Ω = 0.9γ and Ω = 2.1γ is accompanied with rapid and small oscillations, whereas the latter are observed only in the range γt < 400 for Ω = 0.001γ. This implies that non-Markovianity, associated to backflows of information, should be greater for the first two values of Ω. In general, a non-monotonic behavior of non-Markovianity is expected as a function of Ω, since the intensity of oscillations appears to be very sensitive to different values of the modulation frequency (this aspect shall be treated below). We then point out that, although the plot for Ω = 0.001γ seems to perform better than those for Ω = 0.9γ and Ω = 2.1γ within the time range shown in the figures, our calculations instead show that the coherence lasts much longer for these latter values of the modulation frequency. In general, however, reducing Ω to values smaller and smaller, with a fixed nonzero δ, helps to prolonging the time when coherence vanishes. This happens because the final effect is to simply detune the central frequency of the high-Q cavity Lorentzian spectral density from the qubit transition frequency^136,143, as can be deduced from H_eff of Eq. (5) for Ω → 0. Nevertheless, this extreme condition is not relevant to our purpose, since we are interested in the modulation of the qubit frequency during the evolution¹²² and not to a detuning effect.

Figure 3

Qubit coherence ζ(t) as a function of scaled time γt for different values of the modulation frequency, for a fixed modulation amplitude δ = 5γ. Other parameters are: λ = 0.01γ (strong coupling), $\alpha =\beta =1/\sqrt{2}$. Solid-blue line in panel (a) corresponds to the situation in which frequency modulation is off.

Therefore, to gain a more comprehensive understanding of the system evolution, we now study the interplay of the two driving parameters δ and Ω, which is expected to play a significant role in affecting the decoherence process¹⁴⁰. Figure 4 displays the dynamics of coherence, under strong coupling (λ = 0.01γ), when the ratio δ/Ω is tuned such as to assume values which make the n-th Bessel function J_n of Eq. (6) vanish. As seen in Fig. 4(a), the initial coherence can be strongly protected against the noise by increasing the modulation frequency when δ/Ω is such that J₀ vanishes. However, this behavior is not general for larger values of the ratio δ/Ω, as shown in Fig. 4(b–d). The interplay between δ and Ω is thus non-trivial, since it is not sufficient to fix δ/Ω as a zero of the Bessel functions to have long-lasting coherence. In particular, we find that no other settings of δ/Ω and Ω supply a result superseding that occurring for δ = 2.40483Ω and Ω = 5γ (solid blue line of Fig. 4(a)).

Figure 4

Qubit coherence ζ(t) as a function of scaled time γt for different values of the modulation frequency: Ω = 0.05γ (dashed green line), Ω = 0.5γ (dotted red line), Ω = 5γ (solid blue line). The panels correspond to various values of the modulation amplitudes: (a) δ = 2.40483Ω [J₀(2.40483Ω) = 0], (b) δ = 3.83170Ω [J₁(3.83170Ω) = 0], (c) δ = 5.13562Ω [J₂(5.13562Ω) = 0], (d) δ = 6.38016Ω [J₃(6.38016Ω) = 0]. Other parameters are: λ = 0.01γ (strong coupling), $\alpha =\beta =1/\sqrt{2}$.

We then examine the effect of frequency modulation on the dynamics of both QFI F_ϕ(t) and optimal phase estimation δϕ_min(t) for the ratio δ/Ω = 2.40483 under strong coupling (λ = 0.01γ). From Fig. 5 it is immediately seen, as expected, that maintenance of QFI and a significant improvement in phase estimation can be reached during the evolution by increasing Ω up to the value Ω = 5γ.

Figure 5

(a) Quantum Fisher information F_ϕ(t) and (b) optimal phase estimation δϕ_min(t) as a function of scaled time γt for various values of modulation frequency: Ω = 0.05γ (dashed light green line), Ω = 0.2γ (dash-dotted dark green line), Ω = 0.5γ (dotted red line), Ω = 5γ (solid blue line). Modulation amplitude is fixed at δ = 2.40483Ω. Other parameters are: λ = 0.01γ (strong coupling), $\alpha =\beta =1/\sqrt{2}$.

At this stage, one may ask whether the high coherence protection found above is ultimately linked to strong memory effects. The answer requires the knowledge of the degree of non-Markovianity ${\mathscr{N}}$ as a function of Ω under different parameter conditions. This is first done under the weak coupling regime (λ = 3γ) for some modulation amplitudes δ in Fig. 6. These plots disclose that the frequency modulation process can induce non-Markovian features even in the weak-coupling regime¹³⁶, which justifies the existence of oscillations in the time evolution of coherence for Ω = 0.5γ and Ω = 1γ (see Fig. 2(b)). Moreover, this induced non-Markovianity can be reinforced by increasing the modulation amplitude, while no memory effect is observed for small values of Ω tending to turn off the frequency modulation. In Fig. 7 we then display the behavior of non-Markovianity ${\mathscr{N}}$ as a function of Ω under the strong coupling regime (λ = 0.01γ). Panels (a) and (b) of this figure consider the cases when the ratio δ/Ω corresponds to a zero of the Bessel functions, as in Fig. 4: one sees that non-Markovian characteristics can be amplified only for Ω < 1γ, with maximum amplification occurring when δ/Ω is the zero of the Bessel function J₀. Panels (c) and (d) of Fig. 7 instead consider fixed modulation amplitudes δ as Ω changes: it is evident that an increase of δ enlarges the range of Ω where non-Markovianity can be enriched, albeit this is achieved only when Ω < δ. As a general behavior, all plots of Fig. 7 show that ${\mathscr{N}}$ tends to reach approximately the same value as Ω/γ increases, independently of the values of the modulation amplitude δ. Thus, on the basis of our above analysis of qubit dynamics, a more efficient preservation of coherence cannot be related to a higher degree of non-Markovianity (stronger memory effects).

Figure 6

Non-Markovianity ${\mathscr{N}}$ as a function of Ω/γ. Values of other parameters are: δ = 10γ, λ = 3γ (weak coupling), $\alpha =\beta =1/\sqrt{2}$.

Figure 7

Non-Markovianity ${\mathscr{N}}$ as a function of Ω/γ for: (a) δ = 2.40483Ω (J₀(2.40483Ω) = 0) (solid-blue line), δ = 3.83170Ω (J₁(3.83170Ω) = 0) (dotted-red line); (b) δ = 5.13562Ω (J₂(5.13562Ω) = 0) (solid-blue line); (c) δ = 0.1Ω (solid-blue line), δ = 1Ω (dotted-red line); (d) δ = 5Ω (solid-blue line), δ = 10Ω (dotted-red line). The values of other parameters are: λ = 0.01γ (strong coupling), $\alpha =\beta =1/\sqrt{2}$.

To have a deeper physical interpretation behind the effective coherence protection, especially exhibited in Fig. 4(a), we plot in Fig. 8 the time-dependent decay rate Γ(t) (in units of γ), appearing in the qubit master equation of Eq. (18), for some relevant values of Ω. From these plots (see, in particular, panel (d)), it is clear that the main cause of long-lasting coherence is the inhibition of the effective decay rate Γ(t) of the qubit.

Figure 8

Decay rate Γ(t) as a function of scaled time γt for different values of the modulation frequency: (a) Ω = 0.001γ, (b) Ω = 0.05γ, (c) Ω = 0.5γ, and (d) Ω = 5γ. Other parameters are those of Fig. 3(a), that is: δ = 2.40483Ω, λ = 0.01γ (strong coupling), $\alpha =\beta =1/\sqrt{2}$.

Summarizing, according to our analysis, the best performance of coherence preservation controlled by qubit frequency modulation occurs under strong coupling and for values of modulation amplitude δ and modulation frequency Ω given, respectively, by

$$\delta =2.40483{\rm{\Omega }},\,{\rm{\Omega }}=5\gamma ,$$

(24)

which are expressed in units of the spontaneous decay rate of the qubit $\gamma ={\tau }_{{\rm{q}}}^{-1}$.

Experimental context

We discuss here how the previous results may impact in experimental contexts where quantum computing platforms allow for accurate preparation of initial states and control of the qubit. This is particularly the case of superconducting qubits, which have been imposing as the building blocks of quantum computer prototypes based on architectures of circuit quantum electrodynamics (circuit-QED)^127,144.

In order to give a realistic account of our theoretical results, we choose a superconducting transmon qubit embedded in a coplanar resonator. A typical transition frequency for such a qubit is ω₀ ≈ 50 GHz, with a spontaneous (free space) decay rate γ ≈ 10 MHz and size l ≈ 200 μm¹⁴⁵. Within the strong coupling regime, which is the condition of major interest in our study, the considered range of amplitude and frequency modulation parameters is 0 ≤ δ ≤ 10γ and 0 ≤ Ω ≤ 50γ, that guarantee the adiabatic (stationary) regime for the interaction Hamiltonian occurring when $\delta /{\omega }_{0},\,{\rm{\Omega }}/{\omega }_{0}\ll 1$^138,139. Notice that this adiabatic regime appears to be already experimentally achieved for a transmon qubit system in a coplanar resonator¹²⁶. These values also imply ω₀ ± δ ≈ ω₀ and ${\bf{k}}\cdot {\bf{r}}=\mathrm{(2}\pi /{\lambda }_{{\rm{e}}})l\approx ({\omega }_{0}/c)l\ll 1$, where λ_e is the wavelength of the photon emitted by the qubit and c the speed of light. The latter conditions justify the dipole and rotating wave approximations for the system Hamiltonian along the analysis. Moreover, it is known¹³⁶ that the effects due to a condition of non-rotating wave approximation are negligible when ${\omega }_{0}\gg \lambda$, λ being the width of the spectral density J(ω), which is largely satisfied in our study.

From our dynamical results reported above, we know that the coherence of a qubit in a high-Q cavity would naturally (i.e., without frequency modulation) disappear completely at a time ${t}_{{\rm{wc}}}^{\ast }\approx 10{\gamma }^{-1}=10{\tau }_{{\rm{q}}}$ under weak coupling (see Fig. 2(a)) and at ${t}_{{\rm{sc}}}^{\ast }\approx {10}^{3}{\tau }_{{\rm{q}}}$ (after oscillations) under strong coupling (see Fig. 3(a)). These times are meant as lifetimes of coherent superpositions. For the typical superconducting transmon qubits here considered, which have (free space) relaxation times τ_q ~ 100 ns¹⁴⁵, one would obtain ${t}_{{\rm{sc}}}^{\ast }\approx 100$ μs, that is comparable to the current performance achieved in the IBM-Q quantum processor, where a single qubit has average coherence lifetimes of 90 μs¹⁴⁶.

Let us now see to which extent our results with frequency-modulated (FM) qubit are able to prolong coherence lifetimes. Looking at the time behavior obtained under strong coupling with the optimal driving parameters of Eq. (24), as displayed by the solid (top) curve of Fig. 4(a), our calculations estimate the dimensionless time when coherence vanishes at $\gamma {t}_{{\rm{FM}}}^{\ast }\sim {10}^{7}$. This implies a coherence lifetime ${t}_{{\rm{FM}}}^{\ast }\sim {10}^{7}{\tau }_{{\rm{q}}}$ and thus ${t}_{{\rm{FM}}}^{\ast }\sim {10}^{4}{\tau }_{{\rm{sc}}}^{\ast }\approx 1$ s, producing an extension of four orders of magnitude compared to the case of uncontrolled qubit. This achievement seems feasible with the current technology. In fact, quality factors of cavities in circuit quantum electrodynamics can be adjusted to high values such that the cavity spectral bandwidth (or photon decay rate) λ is smaller than the spontaneous emission decay rate γ, thus entering the strong coupling regime λ < γ¹⁴⁵. Moreover, frequency modulation of individual superconducting qubits is already realized^122,126, with the possibility to suitably fix modulation amplitude and frequency to the desired optimal values.

Two-qubit system

After individuating the optimal parameters of the driving external field such as to efficiently shield single-qubit coherence from the detrimental effects of the environment, we now extend our study to a composite system of two separated noninteracting qubit-cavity subsystems, namely A and B. Each qubit-cavity subsystem is structured like that considered in the above Sec. 1, with the transition frequency of each qubit individually modulated.

Initial states and evolved density matrix

The qubits are initially prepared in the extended Werner-like (EWL) states⁷⁵

$${\rho }_{{\rm{\Psi }}}(0)=r|{\rm{\Psi }}\rangle \langle {\rm{\Psi }}|+\frac{1-r}{4}I,\,{\rho }_{{\rm{\Phi }}}(0)=r|{\rm{\Phi }}\rangle \langle {\rm{\Phi }}|+\frac{1-r}{4}I,$$

(25)

where r is a measure of the state purity P, being for both states P = Tr(ρ²) = (1 + 3r²)/4, I is the 4 × 4 identity matrix and

$$|{\rm{\Psi }}\rangle =\mu |{e}_{{\rm{A}}}{e}_{{\rm{B}}}\rangle +\nu |{g}_{{\rm{A}}}{g}_{{\rm{B}}}\rangle ,\,|{\rm{\Phi }}\rangle =\mu |{e}_{{\rm{A}}}{g}_{{\rm{B}}}\rangle +\nu |{g}_{{\rm{A}}}{e}_{{\rm{B}}}\rangle ,$$

(26)

represent the pure part of the state in the form of Bell-like states with |μ|²â€‰+ |ν|²â€‰= 1. This class of states is quite general, including both Bell states and Werner states.

The two subsystems, being separated, evolve independently and are individually governed by the Hamiltonain of Eq. (1), so that the total Hamiltonian is simply H_AB = H_A + H_B. For the EWL initial states above, the density matrix of the two qubits at time t, in the standard (tensor product) computational basis $\{|1\rangle \equiv |{e}_{A}{e}_{B}\rangle ,\,|2\rangle \equiv$$|{e}_{A}{g}_{B}\rangle ,\,|3\rangle \equiv |{g}_{A}{e}_{B}\rangle ,\,|4\rangle \equiv |{g}_{A}{g}_{B}\rangle \}$, takes the X-type structure^75,108

$$\rho (t)=(\begin{array}{cccc}{\rho }_{11}(t) & 0 & 0 & {\rho }_{14}(t)\\ 0 & {\rho }_{22}(t) & {\rho }_{23}(t) & 0\\ 0 & {\rho }_{32}(t) & {\rho }_{33}(t) & 0\\ {\rho }_{41}(t) & 0 & 0 & {\rho }_{44}(t)\end{array}),$$

(27)

where ${\rho }_{ji}^{\ast }(t)={\rho }_{ij}(t)$. The elements of this evolved density matrix can be straightforwardly calculated by the well-known method introduced in ref.⁷⁴ for separated qubits, based on the knowledge of the single-qubit density matrix. We are interested in the dynamics of quantum resources such as entanglement, discord and coherence associated to the evolved two-qubit state ρ(t). In the following we recall the definitions of their quantifiers.

Quantification of quantum resources

Entanglement between subsystems of any bipartite quantum system can be measured by concurrence. For the density matrix given by Eq. (27), the concurrence is⁷⁸

$${\mathscr{C}}(t)=2\,{\rm{\max }}\,\mathrm{\{0},\,{{\rm{\Lambda }}}_{1}(t),\,{{\rm{\Lambda }}}_{2}(t)\},$$

(28)

where ${{\rm{\Lambda }}}_{1}(t)=|{\rho }_{14}(t)|-\sqrt{{\rho }_{22}(t){\rho }_{33}(t)}$ and ${{\rm{\Lambda }}}_{2}(t)=|{\rho }_{23}(t)|-\sqrt{{\rho }_{11}(t){\rho }_{44}(t)}$. The two initial EWL states of Eq. (25) have the same concurrence ${\mathscr{C}}\mathrm{(0)}=2\,{\rm{\max }}\,\mathrm{\{0},(|\mu \nu |+\mathrm{1/4)}r-\mathrm{1/4}\}$ and are thus entangled for r > 1/(1 + 4|μν|).

Discord captures all nonclassical correlations between two qubits beyond entanglement^18,19,22. For a X-structure density matrix, like that of Eq. (27), the analytic expression of discord can be written as¹⁴⁷

$${\mathscr{D}}(t)=\,{\rm{\min }}\,\{{{\mathscr{D}}}_{1}(t),\,{{\mathscr{D}}}_{2}(t)\},$$

(29)

where ${{\mathscr{D}}}_{j}(t)=h({\rho }_{11}+{\rho }_{33})+{\sum }_{i=1}^{4}\,{\lambda }_{i}\,{\mathrm{log}}_{2}\,{\lambda }_{i}+{d}_{j}$ (j = 1, 2), with λ_i being the eigenvalues of the density matrix ρ(t), $h(x)=-\,x\,{\mathrm{log}}_{2}\,x-\mathrm{(1}-x)\,{\mathrm{log}}_{2}\,\mathrm{(1}-x)$, d₁ = h(τ) with $\tau =(1+$$\sqrt{{({\rho }_{11}+{\rho }_{22}-{\rho }_{33}-{\rho }_{44})}^{2}+4{(|{\rho }_{14}|+|{\rho }_{23}|)}^{2}})/2$, and ${d}_{2}=-\,h({\rho }_{11}+{\rho }_{33})-{\sum }_{i=1}^{4}\,{\rho }_{ii}\,{\mathrm{log}}_{2}\,{\rho }_{ii}$ (the explicit time-dependence of the density matrix elements has been omitted for simplicity). The two initial EWL states have the same discord ${\mathscr{D}}(0)$ as a function of r and μ whose explicit expression, obtainable by the previous formulas, is not reported here.

Concerning coherence, the definition given in Eq. (19) is immediately applicable to the two-qubit evolved density matrix of Eq. (27), so to give the two-qubit coherence with respect to the computational basis

$${\zeta }_{2}(t)=|{\rho }_{14}(t)|+|{\rho }_{41}(t)|+|{\rho }_{23}(t)|+|{\rho }_{32}(t)|.$$

(30)

Both the initial EWL states of Eq. (25) have the same coherence ζ₂(0) = 2r|μν| and are thus coherent for any r > 0 provided that μ ≠ 0. Notice that for $\mu =\nu =1/\sqrt{2}$, that is the pure part of the EWL states is a (maximally entangled) Bell state, one has ζ₂(0) = r.

We are now ready to display the dynamics of the relevant quantities defined above.

Time evolution of two-qubit resources

We assume that the subsystems are identical, that is characterized by the same values of qubit-environment parameters, with the optimal frequency modulation parameters fixed as in Eq. (24), under the strong coupling regime with λ = 0.01γ. Such a choice has the advantage to make us directly focus on the best quantitative protection of two-qubit resources attainable by individual qubit frequency modulation under a given strong coupling condition. To acquire information about the protection efficiency of individual frequency modulation, we report the dynamics of ${\mathscr{C}}(t)$, ${\mathscr{D}}(t)$ and ζ₂(t) starting from the initial EWL states of Eq. (25) with r = 1 and $\mu =\nu =1/\sqrt{2}$, which reduces them to the Bell states, respectively, $|{\rm{\Psi }}\rangle =(|{e}_{{\rm{A}}}{e}_{{\rm{B}}}\rangle +\nu |{g}_{{\rm{A}}}{g}_{{\rm{B}}}\rangle )/\sqrt{2}$ and $|{\rm{\Phi }}\rangle =(\mu |{e}_{{\rm{A}}}{g}_{{\rm{B}}}\rangle +\nu |{g}_{{\rm{A}}}{e}_{{\rm{B}}}\rangle )/\sqrt{2}$. Figure 9 immediately shows that individual qubit control by frequency modulation enables lifetimes of quantum resources $\bar{t}$ (that is, the time when they completely disappear) which are three orders of magnitude longer than the lifetimes when modulation is off. In particular, we find ${\bar{t}}_{{\rm{FM}}}\sim {10}^{6}{\tau }_{{\rm{q}}}$ against $\bar{t}\approx 500{\tau }_{{\rm{q}}}$ without frequency modulation. Concerning dynamical details, one sees that while the evolution of coherence is the same for both initial states, discord and entanglement vanish earlier for the initial state ρ_Ψ(0) = |Ψ〉 〈Ψ| compared to ρ_Φ(0) = |Φ〉 〈Φ|. Another difference in the time behavior of quantum resources for the two initial states is that entanglement and coherence coincide at any time for ρ_Φ(0) = |Φ〉 〈Φ| (panels (b) and (d) of Fig. 9), whereas this is not the case for ρ_Ψ(0) = |Ψ〉 〈Ψ|.

Figure 9

Time evolution of two-qubit coherence (TQC) ζ₂(t) (solid blue line), quantum discord (QD) ${\mathscr{D}}(t)$ (dotted red line) and concurrence (CON) ${\mathscr{C}}(t)$ (dashed green line) for the initial states ρ_Ψ(0) (column I) and ρ_Φ(0) (column II), with r = 1 and $\mu =\nu =1/\sqrt{2}$ (Bell states), under λ = 0.01γ (strong coupling). Panels (a and b) uncontrolled dynamics (no frequency modulation). Panels (c and d) controlled dynamics (individual frequency modulation) with optimal parameters δ = 2.40483Ω, Ω = 5γ.

Finally, we analyze an attenuated strong coupling condition λ = 0.1γ for decreasing values of the purity parameter r in the EWL states of Eq. (25), in order to take into account both lower quality factors of cavities and possible imperfections in the initial state preparation. From Figs 10 and 11, as expected, one sees that the amount of all the quantum resources diminishes for smaller r and tends to zero faster. For r = 0.3 entanglement is always zero, since ${\mathscr{C}}(0)=0$. We also observe that two-qubit coherence is the more robust resource of the three ones, entanglement being instead the more fragile. Once again, the general emerging result is that all the quantum resources are efficiently shielded from the noise by qubit frequency modulation, their lifetime $\bar{t}$ being extended of about three orders of magnitude with respect to the case without modulation. This is seen, for instance, in the case r = 1 of Figs 10(a) and 11(a) where ${\bar{t}}_{{\rm{FM}}}\approx {10}^{4}{\tau }_{{\rm{q}}}$, against $\bar{t}\approx 50{\tau }_{{\rm{q}}}$ found in absence of frequency modulation.

Figure 10

Two-qubit coherence (TQC) ζ₂(t) (solid blue line), quantum discord (QD) ${\mathscr{D}}(t)$ (dotted red line) and concurrence (CON) ${\mathscr{C}}(t)$ (dashed green line) for the two-qubit initial state ρ_Ψ(0) as functions of scaled time γt for different degrees of purity: (a) r = 1, (b) r = 0.8, (c) r = 0.5, and (d) r = 0.3. Other parameters are: $\mu =\nu =1/\sqrt{2}$, δ = 2.40483Ω, Ω = 5γ and λ = 0.1γ (strong coupling).

Figure 11

Two-qubit coherence (TQC) ζ₂(t) (solid blue line), quantum discord (QD) ${\mathscr{D}}(t)$ (dotted red line) and concurrence (CON) ${\mathscr{C}}(t)$ (dashed green line) for the two-qubit initial state ρ_Φ(0) as functions of scaled time γt for different degrees of purity: (a) r = 1, (b) r = 0.8, (c) r = 0.5, and (d) r = 0.3. Other parameters are taken as in Fig. 10.

These results assume experimental significance in the context of circuit quantum electrodynamics, where the required individual control of the subsystems has been already discussed above (Sec. 1.5) and a reliable preparation of entangled states of superconducting qubits has been implemented^148,149. For example, entangled states with purity ≈0.87 and fidelity to ideal Bell states ≈0.90 have been generated in the laboratory by using a two-qubit interaction mediated by a cavity bus in a circuit quantum electrodynamics architecture¹⁴⁸. These states may be approximately described as the EWL states of Eq. (25) with r = r_exp ≈ 0.91 and $\mu =\nu =1/\sqrt{2}$¹⁰², which is just the configuration of the initial states we have assumed in our study.

Discussion

In this work we have investigated in detail the effect of individual modulation of qubit transition frequency in protecting quantum resources from the detrimental effect of the environment, constituted by leaky high-Q cavities. Compared to previous studies about qubit frequency modulation, our work has the merit (i) to provide an extensive quantitative analysis of the dynamics of quantum resources under a known control technique of qubit systems and (ii) to supply its potential impact on a cutting-edge technology employed for quantum computer prototypes, such as circuit quantum electrodynamics with superconducting qubits.

We have first put our attention on a single frequency-modulated qubit embedded in a cavity, determining optimal modulation parameters δ (modulation amplitude) and Ω (modulation frequency) and the qubit-cavity coupling regime allowing a long-time preservation of coherence. We have found that, under a strong coupling regime, qubit coherence lifetimes can be extended of orders of magnitude with respect to the case when modulation is off. In particular, when the ratio between the cavity spectral bandwidth λ and the spontaneous emission rate of the qubit γ is λ/γ = 10^−2 we have shown that this coherence lifetime extension can be of four orders of magnitude. We have also seen that the same conditions guarantee a maintenance of quantum Fisher information with a consequent relevant advantage in phase estimation processes of the single-qubit state during the evolution. We have individuated the inhibition of the effective decay rate of the qubit as the mechanism underlying the efficient dynamical preservation of quantum coherence, which is therefore not due to stronger non-Markovianity (memory effects) of the system.

We have then exploited the findings of the single-qubit case to control resources like entanglement, discord and coherence in a composite system of two separated qubits. Such a study, so far unexplored, has revealed itself very useful as a first step towards understanding the scalability of the control procedure to individually addressable subsystems. We have found that individual qubit frequency modulation is still capable to increase the lifetimes of the desired quantum resources of orders of magnitude compared to their natural (uncontrolled) disappearance times. In fact, albeit the gain in the two-qubit system is weaker than in the single-qubit case, with the optimal parameters one may reach lifetime prolongations of three orders of magnitude.

We have discussed the experimental feasibility of the systems, particularly considering the state-of-art achievements in the context of circuit quantum electrodynamics. Setups of superconducting qubits can indeed modulate the qubit transition frequency by external control, reach strong qubit-cavity coupling conditions and permit high-fidelity initial state preparation. The explicit values of the control parameters provided here to achieve the high preservation of quantum resources appear to be feasible and constitute a useful practical information for experimental applications within the so-called noisy intermediate-scale quantum (NISQ) devices¹⁵⁰. Moreover, the technique is straightforwardly applicable to an array of many separated qubit-cavity subsystems to further investigate scalability. Our results provide novel insights for efficient experimental strategies against decoherence towards reliable quantum-enhanced technologies.