Entanglement of superconducting qubits via acceleration radiation

We show that simulated relativistic motion can generate entanglement between artificial atoms and protect them from spontaneous emission. We consider a pair of superconducting qubits coupled to a resonator mode, where the modulation of the coupling strength can mimic the harmonic motion of the qubits at relativistic speeds, generating acceleration radiation. We find the optimal feasible conditions for generating a stationary entangled state between the qubits when they are initially prepared in their ground state. Furthermore, we analyse the effects of motion on the probability of spontaneous emission in the standard scenarios of single-atom and two-atom superradiance, where one or two excitations are initially present. Finally, we show that relativistic motion induces sub-radiance and can generate a Zeno-like effect, preserving the excitations from radiative decay.

We show that simulated relativistic motion can generate entanglement between artificial atoms and protect them from spontaneous emission.We consider a pair of superconducting qubits coupled to a resonator mode, where the modulation of the coupling strength can mimic the harmonic motion of the qubits at relativistic speeds, generating acceleration radiation.We find the optimal feasible conditions for generating a stationary entangled state between the qubits when they are initially prepared in their ground state.Furthermore, we analyze the effects of motion on the probability of spontaneous emission in the standard scenarios of single-atom and twoatom superradiance, where one or two excitations are initially present.Finally, we show that relativistic motion induces sub-radiance and can generate a Zeno-like effect, preserving the excitations from radiative decay.
Circuit Quantum Electrodynamics (cQED) [1][2][3][4] offers both a promising architecture for quantum technologies, such as quantum computers [5,6] and simulators [7][8][9], and a natural arena for the study of quantum field theory and relativistic effects, either in a direct or simulated fashion [10][11][12][13][14].For instance, the dynamical Casimir effect (DCE), produced by the modulation of the boundary conditions of the electromagnetic field at relativistic speeds, has been observed in superconducting devices [15][16][17].Along these lines, it has been shown that DCE radiation possesses several forms of quantum correlations [18][19][20][21][22] that can be transferred to superconducting qubits [23,24].A related phenomenon is the Unruh effect, where an accelerated detector in vacuum should detect thermal radiation [25,26].Recently, some of us have shown how to mimic the generation of acceleration radiation by means of the modulation of the coupling strength of a superconducting qubit [27], a phenomenon resembling the cavity-enhanced Unruh effect [28][29][30].The simulation in a superconducting architecture of both phenomena, DCE and acceleration radiation, relies on the possibility of performing an ultrafast variation of the magnetic flux threading a superconducting quantum interferometric device (SQUID) [31][32][33].
In this paper, we consider a superconducting circuit setup in which two superconducting qubits interact with the same resonator mode and effectively move at relativistic speeds, see Fig. (1).The simulation of the relativistic motion of the qubits comes from the modulation of the coupling strength between the qubits and the resonator, which can be interpreted as the qubits movement and activates the counterrotating terms of the quantum Rabi Hamiltonian.We analyse the role of the generated acceleration radiation in several collective properties of the qubits.First, we consider an initial state with no excitations and find the conditions for an efficient generation of stationary entangled states.We find several optimal scenarios for entanglement production: either both qubits move resonantly with the (or half the) natural frequency of the cavity, or one qubit remains static while the other moves at twice the cavity frequency.Second, we analyse the effect of relativis- tic motion on the spontaneous emission rate when one or two qubit excitations are initially present.Namely, we add the ingredient of relativistic motion to the celebrated Dicke scenario of single-atom and two-atom superradiance [34], which has been recently implemented in a circuit QED architecture [35].We will show that the counterrotating dynamics generated by the motion of the qubits tends to suppress the superradiance.Moreover, we find experimental conditions under which the qubit decay is completely frozen, giving rise to a Zeno-like effect induced by the continuous modulation of the coupling strength [36][37][38].In this second case, the optimal scenario for the appearance of Zeno-like effect correspond to a synchronised motion of the qubits at twice the frequency of the cavity, which generates an anti-Jaynes-Cummings dynamics in both qubits that prevents them from spontaneous emission.
Entanglement and acceleration radiation.-TheHamiltonian of the system describes two superconducting qubits of frequency gaps ω q coupled to a single resonator mode of fre-arXiv:1606.06170v2[quant-ph] 7 Apr 2017 Here, σ z are Pauli matrices for the qubits, and a (a † ) is an annihilation (creation) operator for the resonator mode.The interaction Hamiltonian depends on the qubits position as with g the coupling strength and x q the qubit position [39].
In order to simulate the motion of the qubits, which are located in fixed positions, we modulate the coupling with external drivings, such that the interaction Hamiltonian for a qubit reads H I (x q ) = g cos (f 0 + ∆f cos(ω d t))σ x (a † + a), and kx q = f 0 + ∆f cos(ω d t).The velocity of the qubits vary harmonically in time, with the maximum value of ≈ λω d .For λ = 2L c = 1cm and ω d = 10 GHz, we reach values of ≈ 10 8 m/s = c/3.
Initially, we consider the system in the ground state for the qubits and the resonator mode |g 1 g 2 0 .In order to determine the degree of entanglement between the qubits after a certain interaction time T , we compute the concurrence, which up to second order in perturbation theory with respect to g /ω reads C = 2 Max{|X| − P e , 0} [40].Here, X is the amplitude for photon exchange between the qubits, X = 0|T (S + 1 S + 2 )|0 , with T the time-ordering operator.P e is the probability of emitting a photon, P e = 0|S − 1 S + 1 |0 , with S + = − i g T 0 e iω q t dt cos(kx q )(e iωt a † + e −iωt a)dt = −(S − ) † .In this configuration, if both qubits are at fixed positions, the emission and photon exchange are counterrotating processes, that is, related to the breakdown of the rotatingwave approximation (RWA).These processes will be significant only for ultrastrong couplings or short interaction times.However, the motion of the qubits can excite these counterrotating terms of the Hamiltonian, giving rise to a sizeable emission of photons by a sort of cavity-enhanced Unruh effect [28].
We will analyse now under which conditions this phenomenon can be exploited to efficiently generate entanglement between the qubits.We can gain first insights by using analytical techniques.For the sake of simplicity, we assume that ω q 1 = ω q 2 = ω q and g 1 = g 2 = g.Moreover, for the sake of simplicity, we consider that the harmonic oscillations of the qubits with frequency ω d preserve its relative distance, , and λ is the wavelength associated with the cavity frequency FIG.2: Concurrence C of two qubits initially located in the centre of the resonator and oscillating from mirror to mirror with frequencies ω d1 and ω d2 , respectively.For coupling constants g1 = g2 = g = 0.02, initial state |g1 g2 0 , qubit decay parameter Γ = 0.002, and T2/T1 = 0.67, we consider two cavity decay rates: a) κ = 0.002, and b) κ = 0.2 (bad-cavity limit), in units of ω.These numerical results correspond to a broader parameter range, not limited by the perturbative approximation gT < 1, which for this case of g = 0.02 breaks for ≈ 8ωt/2π.ω.Then, we have where ∆ is the detuning between the qubits and the cavity, ∆ = ω q − ω, and J n (x) are Bessel functions of the first kind.By inspection of Eq. ( 3), we observe that both X and P e are, in general, oscillating functions.However, under certain resonant conditions, the oscillations are suppressed and these magnitudes grow monotonically in time.For instance, for negligible detuning ∆ = 0, and frequencies ω d = ω = ω q , we find that |X| , by keeping only the non-oscillating terms.Since 1, the entanglement grows quadratically in time as Therefore, we predict an entanglement resonance around ω d = ω.These analytical results are limited by the perturbative approximation employed, which assumes that gT < 1.Even in the weak coupling regime, the perturbative approximation would eventually break down.In Fig. (2), we plot the results of numerical simulations which generalise our analytical insights.The dynamics is governed by a master equation where we introduce a cavity decay rate κ, a decay parameter Γ accounting for dissipative processes, as well as a decay Γ ϕ for the dephasing of the qubits.The energy relaxation time and phase coherence time are denoted with T 1 = 1/Γ and T 2 = 1/Γ φ , respectively.We consider realistic parameters, achievable with present technology [41].This allows us to analyse the long-term dynamics of the system and to consider more general types of motion with ω d1 = ω d2 , in which the relative distance of the qubits is no longer preserved.Numerical simulations confirm the generation of a high degree of entanglement in the case of ω d1 = ω d2 = ω d = ω, as expected for short-time dynamics.In the long-term dynamics, we observe non-trivial entanglement oscillations, where maximum values are achieved at particular times shown in Fig. (2).Another optimal scenario for entanglement generation appears when one qubit is effectively moving with frequency ω d1 = 2ω, and the other remains static, ω d2 = 0.Under these circumstances, the first qubit is ruled by an anti-Jaynes Cummings (anti-JC) dynamics which maximises the counterrotating emission of photons [27].In this case, the concurrence reaches its maximal value, and the amplitude of initially perfect collapse-revival cycles eventually diminish.Asymptotically, entanglement exhibits small fluctuations around a mean value close to one.Moreover, if we also consider the bad-cavity limit, where κ g Γ, entanglement oscillations are smoothed out, and highly entangled stationary states are reached, see Fig. (2b).We extend our analysis of the generation of entanglement between both qubits for the case in which the cavity is out of resonance from both drivings and qubits frequencies.Although we have analytically seen that the concurrence increases when the resonant condition is fulfilled, we expect that considering a detuned cavity will enhance the quantum correlation between the qubits (see Supplementary Information).
Single-atom superradiance and Zeno-like effect.-Inhis seminal work, Dicke showed that the decay of atomic emitters is enhanced by the presence of other atomic emitters [34].The simplest case, called single-atom superradiance, consists of a single emitter in an excited state influenced by the proximity of another emitter, even if the latter is in the ground state.This gedanken experiment has been recently realised in a circuit QED architecture in the bad-cavity limit [35].Here, we analyse the effects of relativistic motion in this scenario.To this end, we consider the initial state |g 1 e 2 0 , and discuss, firstly, the effects of the relativistic motion of the second qubit in its decay, and, secondly, the effects of the presence of the first qubit in the decay of the second one.We observe that the relativistic motion of the second qubit, encoded in ω d2 , tends to inhibit its decay leading to a decreased emission FIG.3: Motion effects in single-atom superradiance, observed in the probability P of excitation of the second qubit Q2.We consider the first qubit Q1 decoupled, g1 = 0, and the second moving with frequency ω d2 , for an initial state |g1 e2 0 .We show the behaviour for different velocities of the second qubit, ranging from the static case, ω d2 = 0, to a velocity of ω d2 = 2ω.For a coupling constant g2 = 0.02, a qubit decay parameter Γ = 0.002 and T2/T1 = 0.67, we consider two cavity decay rates: a) κ = 0.002, and b) κ = 0.2 (bad-cavity limit), in units of ω. rate, known as sub-radiance, see Fig. 3. Again, we can get some insight on the system dynamics from first-order analytical computations.Since in this case the qubit is initially excited, the probability of emitting a photon is not given by Eq. (3).In particular, in the resonant case (∆ = 0), we have . Then, considering the trajectories for both qubits such that . This means that for a static qubit -ω d2 close to 0-, the probability of emission grows quadratically in time P e g 2 T 2 , while for frequencies of motion significantly different from 0 the probability oscillates with an amplitude which decreases with ω d2 .Therefore, for large enough ω d2 the probability of emission is suppressed, as can be seen in Fig. (3).Note that the maximum acceleration of the qubit motion is proportional to ω 2 d2 .Thus, the larger the acceleration is, the larger the suppression of the probability of emission.Then, the subradiance can be seen as another relativistic effect, hitherto unexplored.At first glance, this subradiant dynamics might look surprising, since relativistic accelerated motion is typically associated to the emission of photons.However, note that both FIG.4: Motion effects in single-atom superradiance, observed in the probability P of excitation of the second qubit Q2.We consider the influence of the movement of the first qubit Q1 analysing different velocities ω d1 and the decoupled case, g1 = 0.The second qubit is moving with frequency ω d2 , for an initial state |g1 e2 0 .We compute for a coupling constant g1 = g2 = 0.02 in the cases with the first qubit coupled, a qubit decay parameter Γ = 0.002 and T2/T1 = 0.67, and a cavity decay rate κ = 0.2 (bad-cavity limit), in units of ω.We show the modification in the behaviour in the case of the second qubit a) static ω d2 = 0, and b) moving with ω d2 = ω.phenomena, Unruh effect and subrradiance, are activated by the counterrotating terms of the Hamiltonian, which become dominant for large enough ω d2 associated with relativistic motion.While non-RWA dynamics gives rise to emission of photons when the qubit and cavity start in the ground state, in the present case the initial state |e 2 0 would be stationary in the presence of only non-RWA terms.The decay dynamics of the second qubit is effectively frozen, that is, we observe a Zeno-like effect generated by the continuous modulation of the qubit-cavity coupling strength, which has an effect similar to a continuous monitoring of the system [36][37][38].
Not only the relativistic motion of the second qubit, encoded in ω d2 , tends to inhibit its decay, but also the motion of the first qubit, ω d1 , has a significant effect on the emission rate for certain values of ω d1 .In order to analyse this influence, we compare in Fig. 4 the probability of excitation of the second qubit for the case in which the first qubit is decoupled, and the case in which it is coupled moving at different relativistic speeds, ω d1 .Firstly, we consider the extreme case with the second qubit static, ω d2 = 0, and we observe that a first qubit relativistic speed ω d1 = 2ω leads to a decreased emission rate, FIG.5: Zeno-like effect in the probability P of excitation of two qubits moving with the same frequency ω d1 = ω d2 = 2ω for different amplitudes of oscillation for the first qubit Q1.We consider the initial state |g1 e2 0 , a cavity decay rate κ = 0.1, coupling constants g1 = g2 = g = 0.01, and qubit decay parameter Γ = 0.001 and T2/T1 = 0.67, in units of ω, and that the first qubit is initially placed at Lc/4, with Lc the cavity length.We show the excitation probabilities of both qubits, Q1 and Q2, with frequencies ω q 1 = ω q 2 = ω, and adding a detuning ∆ = 0.1ω between the first qubit and the cavity, for an amplitude of motion of the first qubit a) Lc/4, and b) Lc/16.known as sub-radiance, whereas for other combinations of frequencies, the behaviour of the probability of excitation of the second qubit remains unaltered.Secondly, we analyse the case with the second qubit moving with ω d2 = ω, and observe a dramatic change in the decay rate of the second qubit for the same frequency ω d1 = 2ω, as in the previous case.We notice that the emission rate of the second qubit is also slightly modified for a relation of frequencies ω d1 = ω d2 = ω that generates entanglement.However, a further analysis in the relation of the generated entanglement and single-atom superradiance allows us to discard drastic influences of the former in the decay rate (see Supplementary Information).We also interpret the frozen dynamics of the second qubit for ω d2 = ω, and for ω d2 = ω and ω d1 = 2ω, as a Zeno-like effect [36][37][38].
We extend our analysis to other scenarios by considering different frequencies and initial conditions.We confirm the Zeno-like effect when we reduce the oscillation amplitude and the excitation probability of the first qubit, which enhances even more the effect on the second qubit (see Fig. 5).Twoatom superradiance is a collective effect consisting in an enhancement of the decay rate of two emitters with respect to their individual ones [34].We consider the effects of relativistic motion on the decay of two artificial atoms, |e 1 e 2 0 and observe a sub-radiance phenomenon, as in the case of singleatom superradiance (see Supplementary Information).
Implementation in superconducting circuits.-Themodel described in Eq. ( 1) can be implemented in a circuit QED architecture [42], using a single-mode transmission line resonator (TLR) interacting with two tunable-coupling superconducting qubits.In order to observe all the phenomenology so far described in a single device, it is required independent tuning of the qubit transition frequencies and of the qubitcavity coupling strengths.Tunable coupling superconducting qubits [41,43,44] coupled to a TLR and tuning of effective couplings over nanosecond time-scale [45] have been proven in circuit QED architectures.
Conclusions.-We have proposed a possible realisation in which simulated relativistic motion generates true entanglement between artificial atoms, while protecting them from spontaneous emission in a Zeno-like effect.A natural extension would be to consider the effects of multi-atom relativistic motion, with a study of superradiant phase transition.Both the ability of generating entanglement and state protection may pave the way for new applications in superconducting quantum technologies.

RELATION OF CONCURRENCE AND SINGLE-ATOM SUPERRADIANCE
We have explored numerically the influence of the movement of the first qubit in the emission rate of the second qubit in Fig. 4 of the main text.In order to better interpret the possible relation between the generated entanglement and the influence of the first qubit, we show in Fig. 7 the concurrence for the relation of frequencies ω d1 and ω d2 analyzed in Fig. 4. We conclude that a combination of frequencies generating entanglement is not relevant for observing significant changes in the probability of the decay of the second qubit, although for ω d1 = ω d2 = ω it results in a slight modification of the decay rate.FIG.7: Concurrence C of two qubits moving with different velocities ω d1 and ω d2 , previously analyzed in the light of single-atom superradiance.We consider a coupling constant g1 = g2 = 0.02, a qubit decay parameter Γ = 0.002 and T2/T1 = 0.67, and a cavity decay rate κ = 0.2 (bad-cavity limit), in units of ω.We show results for different velocities ω d1 of the first qubit, when the second qubit is a) static ω d2 = 0, and b) moving with ω d2 = ω.FIG.8: Probability P of excitation of two qubits, Q1 and Q2, moving with different frequencies ω d1 , ω d2 , and a qubit decay parameter Γ = 0.001 and T2/T1 = 0.67, in units of ω. a), b) Both qubits are initially excited.c), d) The initial state of the qubits is | +1 +2 .a) , c) κ = 0.001, and b), d) κ = 0.1 (bad-cavity limit), in units of ω.All figures show the case of both qubits coupled to the cavity with coupling strengths g1 = g2 = g = 0.01, ω d1 = ω d2 = 2ω, and the case of an uncoupled first qubit g1 = 0, with a static second qubit ω d2 = 0 with g2 = 0.01.

TWO-ATOM SUPERRADIANCE AND COLLECTIVE ZENO-LIKE EFFECT
We consider the phenomenon of two-atom superradiance in our setup.We compare in Fig. 8 the decay rate of both qubits moving with relativistic velocities, encoded in ω d1 and ω d2 , with respect to their individual decay rates, i.e., with respect to the decay of a single static qubit.We observe that the relativistic motion alters the emission rate and reduces drastically for both qubits the probability of decay, leading to a sub-radiance phenomenon.Both in the good cavity and bad cavity limit, we observe a collective Zeno-like effect, namely the effective freezing of the decay of both qubits.This ties in with the prediction of the expected anti-JC dynamics at this particular driving frequency [27].Finally, if we consider a different initial state where both qubits are in a superposition of their ground and excited states, | + 1 + 2 0 , with |+ = 1/ √ 2(|e + |g ), we observe that the qubits are driven fast to their excited states, and the ensuing dynamics follows then the one already observed for an initial |e 1 e 2 .

IMPLEMENTATION DETAILS
With the aim of providing a specific implementation proposal, we focus on a three island superconducting qubit [43,44].The qubit scheme consists in two shunted SQUID loops described by a model composed of two interacting anharmonic oscillators, whose dynamics can be effectively restricted to their two lowest eigenstates.Independent control over the corresponding transition frequencies results in a completely tunable qubit.We denote |0 a the ground and |1 a the first excited state of the oscillator a, and we call ω a the corresponding transition frequency.Same notation will apply for the oscillator b.The qubit logical levels are given by the ground |E 0 and first-excited |E 1 states of the collective system, whose structure depends on the ratio between the frequencies ω i .When the two anharmonic oscillators are detuned, the lowest collective eigenstates are given by |E 0 = |0 a |0 b and |E 1 = |1 a |0 b , where we set ω a < ω b .In this configuration, the qubit is strongly coupled to the TLR.On the other hand, when the two anharmonic oscillators are nearly degenerate ω a ≈ ω b , the first collective excited state is given by Such state corresponds to an antiparallel configuration of the dipoles of the SQUID loops, hence creating a quadrupolar moment that does not couple with the resonator.Swapping between the two collective energyconfigurations, the effective coupling can be continuously tuned in real time, without exciting higher energy levels [46], which are detuned by at least 1 GHz from the primary qubit transition Using this scheme, tuning of the qubit-cavity coupling strength has been experimentally achieved in the range from 40 MHz to less than 200 KHz [44].
System initialization is trivial in all the cases considered in this manuscript, as only ground-state cooling and single-qubit gates are required.Fast read-out of the qubit state can be implemented [41] driving a transition to the second-excited collective state |E 1 → |E 2 .In the resonant case, ω a = ω b , such state is given by |E 2 = (|1 a |0 b + |0 a |1 b ) / √ 2 and it is strongly coupled to the resonator.The qubit excitation states can be then individually obtained measuring the state-dependent shift on the TLR resonant frequency.

FIG. 1 :
FIG.1: Two superconducting qubits strongly coupled to a single resonator mode and driven with frequencies ω d1 and ω d2 simulating harmonic relativistic speeds.The resonator of length Lc is initially in the vacuum while the qubits, which are located in the middle of the cavity x = 0, are initially a) both in the ground state, b) one in the excited state and the other in the ground state, c) both in the excited state.The red wavy arrows indicate emission or absence of emission of photons from the qubits, showing a subradiant behaviour.