Ancilla-Assisted Generation of Photons from Vacuum via Time-Modulation of Extracavity Qubit

We propose a scheme for the generation of photons from a vacuum via time-modulation of a quantum system indirectly coupled to the cavity field through some ancilla quantum subsystem. We consider the simplest case when the modulation is applied to an artificial two-level atom (we call ‘t-qubit’, that can be located even outside the cavity), while the ancilla is a stationary qubit coupled via the dipole interaction both to the cavity and t-qubit. We find that tripartite entangled states with a small number of photons can be generated from the system ground state under resonant modulations, even when the t-qubit is far detuned from both the ancilla and the cavity, provided its bare and modulation frequencies are properly adjusted. We attest our approximate analytic results by numeric simulations and show that photon generation from vacuum persists in the presence of common dissipation mechanisms.


Introduction
The dynamical Casimir effect (DCE) designates a plethora of phenomena characterized by the generation of photons (or quanta of some other field) from vacuum due to time-dependent variations of the geometry (dimensions) or material properties (e.g., the dielectric constant or conductivity) of some macroscopic system (see, e.g., the reviews [1][2][3][4][5]). It was initially investigated for Electromagnetic (EM) field in the presence of non-uniformly accelerating mirrors and cavities with oscillating boundaries or time-dependent material properties [6][7][8][9][10], but the concept was later extended to optomechanical systems [11,12], Bose-Einstein condensates and ultracold gases [13][14][15], polariton condensates [16], and spinor condensates [17,18]. Recently, DCE was implemented experimentally via periodical fast changes of the boundary conditions in circuit Quantum Electrodynamics architecture (circuit QED) [19][20][21][22] and Bose-Einstein condensates [23]. In addition to serving as a direct proof of the vacuum fluctuations [5], from the practical point of view DCE can be employed to generate non-classical states of light or of an ensemble of atoms [20,24,25].
The circuit QED architecture [26][27][28][29][30] is a handy platform for the implementation of DCE and its generalizations, since both the cavity's and artificial atoms' properties can be rapidly modulated by external bias, e.g., magnetic flux [31,32]. In particular, when the atom is directly coupled to the field via the dipole interaction (described by the Quantum Rabi Model [33]), a resonant time-modulation of the atomic transition frequency or the atomfield coupling strength can be used to generate photons and light-matter entangled states from the initial vacuum state [34][35][36][37]. In this case, one can view the atom as a microscopic constituent of the intracavity medium that shifts its effective frequency; moreover, such scheme benefits from leaving the Fock states of the cavity field time-independent (as opposed to the standard case of time-varying cavity frequency, when the annihilation operator and the Fock states depend explicitly on time [9]). These non-stationary circuit QED setups exhibit several important phenomena beside photon generation from vacuum, e.g., generation of atom-field entangled states and novel non-classical states of light [38,39], quantum simulation [40][41][42], implementation of quantum gates [43], engineering of effective interactions [44], implementation of quantum thermal engines [45,46], photon generation and atom-field effective coupling via multi-photon transitions [47,48], anti-dynamical Casimir effect (coherent annihilation of excitations due to external modulation) [49][50][51][52], photon generation by both temporal and spatial modulation in metamaterials [53], vacuum Casimir-Rabi oscillations in optomechanical systems [54], etc. [5].
In this paper, we investigate whether photons can also be generated from the vacuum by modulating an artificial atom that does not interact directly with the cavity, but instead is indirectly coupled to the field through some auxiliary subsystem-the ancilla. Such a coupling scheme may have several reasons and applications. For instance, the artificial atom can be designed specifically to withstand fast external modulation of arbitrary format, at the expense of null coupling to the cavity field but large coupling to other subsystems (possibly to different kinds of artificial atoms); or the atom can be placed outside or at the end of the cavity (at the node of the electric field) to minimize the influence of external driving on the cavity field and increase the cavity quality factor. In addition, the modulated artificial atom can be designed to couple selectively to multiple cavities by means of different stationary ancillas, which are constructed with reduced dissipative losses and enhanced atom-field coupling strengths (ultrastrong coupling [33,55], for instance). Independently of the concrete scenario, it seems timely to investigate whether such indirectly coupled time-modulated atom can be harnessed to generate photons from vacuum or engineer some useful effective interactions, and under which conditions these processes are optimized.
We address analytically and numerically this issue by considering the simplest scenario in which the time-modulated artificial atom is a qubit ("t-qubit", for shortness) and the ancilla is a stationary qubit dipole-coupled to both the cavity field and the t-qubit. We find that photon generation with sufficiently large transition rates is possible provided there is a fine tuning of both the modulation frequency and the bare frequency of the t-qubit, which depend on all other system parameters.
This paper is organized as follows. The mathematical formulation of our proposal is given in Section 2, and in Section 3 we present a closed analytic description of the dynamics in terms of the system dressed-states. In Section 4, we confirm our analytic predictions by exact numeric simulations and illustrate typical system behavior in different regimes of operation. In particular, we show that the initial vacuum state can be deterministically driven either to states with only two excitations or states with multiple excitations, even in the presence of weak dissipative effects. Section 5 contains the conclusions.

Mathematical Formulation
Our tripartite system is represented schematically in Figure 1, and is described by the Hamiltonian (h = 1) (1) The terms in the first brackets describe the free cavity field and the t-qubit coupled to the ancilla, while the terms in the second brackets describe the ancilla qubit and its coupling to the cavity field. ν is the cavity frequency,n =â †â is the photon number operator andâ (â † ) is the annihilation (creation) operator. The t-qubit operators areσ e = |e e|, σ z = |e e| − |g g|,σ − = |g e|,σ + =σ † − andσ x =σ + +σ − , where |g (|e ) denotes the ground (excited) state. For the ancilla, the operators are similar and are indicated by the upper index (a), while its ground and excited states are denoted as |g a and |e a , respectively. We assume that the ancilla, with constant transition frequency Ω a , interacts directly with the cavity field via the dipole interaction with the time-independent coupling strength g. The t-qubit is not directly coupled to the cavity field; instead, it interacts with the ancilla via the dipole interaction with the coupling constant h. We shall derive closed (approximate) analytic description of the dynamics for moderate coupling strengths, g, h 0.05ν, but our numeric simulations will also explore some interesting phenomena in the ultrastrong coupling regime with g, h ∼ 0.1-0.3ν. Notice that the direct qubit-qubit coupling occurs naturally in many circuit QED setups [56][57][58][59], and the coupling constant h can be calculated in terms of the parameters of the superconducting quantum interference devices (SQUIDs) that form the qubits [60,61]. A related case in which the t-qubit is coupled simultaneously to two cavity modes was recently studied in [62], while the possibility of coupling distant qubits using a chiral ring resonator was analyzed in [63]. We assume that the transition frequency of the t-qubit is modulated externally as where Ω 0 is the bare (average) frequency, ε Ω 0 is the modulation amplitude and η is the frequency of modulation. Notice that a periodic external modulation of the system Hamiltonian is a current topic of research in many areas of physics, e.g., two-dimensional electron systems with Rashba spin-orbit coupling irradiated by an off-resonant highfrequency electromagnetic field [64]. Moreover, the generalization of our scheme to nonharmonic modulations is straightforward using the Fourier decomposition [35], while a clever choice of the time-dependence of the modulation frequency η(t) could enhance the photon generation process and originate novel dynamical behaviors [34,65].
In circuit QED, the Hamiltonian alone does not describe accurately the dynamics due to the system coupling to the environment, so instead of the Schrödinger Equation (SE) one has to use the master equation for the density operatorρ where the LiouvillianL accounts for the influence of the environment [66,67]. The form of the Liouvillian depends on the spectral density of the reservoir and the type of coupling to the system [68]. For moderate coupling strengths one can use the standard Markovian master Equation (SMME) of quantum optics [69] L =L d +L ph +L where the superoperatorL d (L ph ) describes the energy damping (pure dephasing) of the t-qubit by a thermal reservoir.L ph have similar meaning for the ancilla, and L κ describes the cavity damping. Indeed, it was shown in [50,65] that in similar setups the discrepancy between this equation and a more rigorous one (a microscopic derivation taking into account the qubit-resonator coupling [70]) is small for the coupling strengths g 0.1ν, while the qualitative agreement is excellent.
In this paper, we consider the zero-temperature limit of the SMMÊ where γ (γ ph ) is the atom relaxation (pure dephasing) rate, κ is the cavity damping rate and is the so called Lindblad superoperator that preserves the hermiticity, normalization and positivity ofρ [68]. For numeric simulations in Section 4 we shall adopt the following parameters: This represents the presumable situation in which the t-qubit is exposed to moderate dissipation due to external modulation, while the ancilla is less susceptible to dissipation by proper design. These dissipative rates were assumed sufficiently small yet readily achievable experimentally [71][72][73]. Such approximate treatment is sufficient to assess the feasibility of our scheme in realistic situations, and the parameters in Equation (7) can be viewed as a benchmark for sufficient dissipation rates. We note that an accurate description of the dissipative dynamics would require solving the microscopic master Equation [70], for which the spectral densities of the baths must be known. Moreover, the driving field could cause the system heating. Since in our scheme the t-qubit is placed outside or at the end of the cavity, such undesirable effects are minimized and could be taken into account by the additional dissipative term γ heat D[σ + ] in Equation (4).

Analytic Description
The analysis is simplified by introducing a new "conjoint" atomic basis {|A i , i = 0, . . . , 3}, in which |A i are the eigenstates of the two-atom time-independent HamiltonianĤ a = Ω 0σe + Ω aσ (a) e + hσ xσ (a) x containing the bare t-qubit frequency Ω 0 : Here W ± = (Ω a ± Ω 0 )/2, D ± = W 2 ± + h 2 and N i are the normalization constants (with the dimension 1/frequency); some useful approximate formulae for the eigenstates |A i are given in the Appendix A. The corresponding eigenenergies are In the basis {|A i } the time-independent system Hamiltonian readŝ and the total Hamiltonian isĤ =Ĥ 0 + ε sin(ηt)σ e .

Spectrum Far from Degeneracies
Far from any degeneracy between the eigenenergies ofĤ 0 (g = 0) (system Hamiltonian without the matter-field coupling) the spectrum ofĤ 0 can be found from the nondegenerate perturbation theory. We use the orthonormal complete basis |A n i ≡ |A i ⊗ |n (where |n is the Fock state of the field with n ≥ 0 and i = 0, . . . , 3) and denote the eigenstate (eigenvalue) ofĤ 0 with the dominant contribution of the state |A n i as |A i , n (λ i,n ). For example, to the second order in g one obtains the (non-normalized) state x |A j and the coefficients Ξ (k) n are given in the Appendix A. The corresponding eigenenergy reads Expressions for other states |A i , n and eigenenergies λ i,n can be obtained similarly, but they are not required for the present work.

Spectrum near Degeneracies
Most accessible applications of our scheme explore the regime of parameters in which two states of the subspace A n = {|A n 0 ,|A n−1 In this case, the non-degenerate perturbation theory fails, but one can obtain excellent analytic results by expandingĤ 0 in the subspace A n as where X = νn + λ 0 ,Î is the 4 × 4 identity matrix and (the parameters of M n are given in the Appendix A). The eigenvalues Λ n,i of M n (with i = 1, . . . , 4 for a given n, where n labels the subspace A n ) can be found exactly using the Ferrari's method. The normalized eigenstate corresponding to the eigenvalue Λ n,i is where φ (k) n,i are the probability amplitudes of the conjoint atomic state |A k (see the Appendix A for the derivation). The eigenenergy of the state |φ n,i is denoted as λ φ n,i and reads Entropy 2023, 25, 901 6 of 20 λ φ n,i = X + Λ n,i . For example, according to our notation, near the degeneracy the eigenstate |A 0 , n of Υ n is the state |φ n,i=J for which φ (0) n,i=J is the largest among the four states {|φ n,i }; the eigenenergy λ 0,n is the corresponding function λ φ n,i=J (so λ 0,n can be a discontinuous function of parameters near the degeneracy).
The above diagonalization procedure has one drawback, it neglects the coupling of the subspace A n to the neighboring subspaces A n±2 , carried by the counter-rotating terms g(âσ in the HamiltonianĤ 0 (notice that the counter-rotating terms were fully taken into account within the subspace A n ). For small values of g the main effect of the neglected contributions is to introduce small frequency shifts ("Bloch-Siegert" shifts [35]) to the eigenfrequencies of the bare HamiltonianĤ 0 (g = 0). To include these cor- has the same structure as M n (with ν D ≡ ν + D + ), so its eigenvalues can be found as previously. Denoting the eigenvalues ofM n in the increasing order asΛ n,i (i = 1, . . . , 4), the frequency shifts of the states in the subspace B n are found as δ n−2 . Now one can insert these shifts into the matrix (13) by replacing 0 , obtaining analytically the eigenvalues and eigenstates of the HamiltonianĤ 0 near the degeneracy points. As will be shown in Section 4, for moderate coupling strengths g 0.05ν, this procedure is in excellent agreement with exact numeric results.
For higher values of g or direct transitions involving n > 2 photons, the above diagonalization procedure becomes improper (since A n only contains states differing by at most two photons) and must be generalized by adding more states to the subspace A n (forming a larger subspace A gen n ). However, since there is no algebraic solution to general polynomial equations of degree five or higher with arbitrary coefficients (Abel-Ruffini theorem), it is easier to perform the diagonalization numerically, as will be completed in Section 4.1.

Dynamics in the Dressed-Basis
The unitary system dynamics is straightforward in terms of the eigenstates (dressedstates) ofĤ 0 , which can be found either numerically or analytically as in the previous subsections. Expanding the wavefunction as where |ϕ l and E l are the eigenstates and eigenvalues ofĤ 0 in the increasing order (E l+1 ≥ E l ), we obtain for the probability amplitudes where E lm = E l − E m is the energy difference. For the realistic case of weak modulations, ε η, one can neglect the rapidly oscillating terms to obtaiṅ where Therefore, for the resonant modulation frequency η = |E lm | the external perturbation induces transition between the dressed-states |ϕ m and |ϕ l with the transition rate |R m;l |. From the practical standpoint, the numeric results are obtained by finding the eigenvalues E m and eigenstates |ϕ m ofĤ 0 in the basis {|A n k }, where k = 0, . . . , 3 and n = 0, . . . , n tr , and then evaluating R m;l and E lm . n tr stands for the maximum number of photons fixed by the truncation procedure, which does not affect the results for eigenstates with photon numbers n n tr (here the value n tr = 30 was enough). For analytic calculations, one simply uses the closed form expressions for dressed-states and eigenenergies found in Section 3.1 (for the states |A i , n far from degeneracy points, with eigenenergies λ i,n ) and Section 3.2 (for dressed-states |φ n,i near degeneracy points, with eigenenergies λ φ n,i ) to evaluate the transition rate and energy differences. Since all the states were expanded in the common basis {|A n k }, such calculations are long but straightforward. In this work, we are primarily interested in photon generation from the initial (nondegenerate) ground state of the system |A 0 , 0 . Near the degeneracy points one can obtain approximate analytic expression for the transition rate using the formulae (11) and (15). To the first order in g, the transition rate to the two-excitations state |φ 2,i of subspace A 2 is where Analogously, the transition rate between the dressed-states |φ n,i and |φ n+2,j is When the modulation frequency matches only one possible value |E lm | (with non-zero transition rate R m;l between the respective eigenstates), the dressed-states |ϕ m and |ϕ l become resonantly coupled and exhibit sinusoidal behaviors (see Figures 2 and 3 below), |A m (t)| 2 = cos 2 (R m;l t) and |A l (t)| 2 = sin 2 (R m;l t) (assuming that only A m was initially nonzero). On the other hand, when several values |E lm | are close to the modulation frequency (with the mismatches η − |E lm | smaller or of the order of |R m;l |), then several dressedstates can become simultaneously coupled and the dynamics becomes more intricate (see . We also note that the neglected rapidly oscillating terms introduce small corrections to the resonant frequencies |E lm | [35,49], which are found numerically in this paper. The ground state can also be directly coupled to dressed-states with more than two excitations. To obtain reliable analytic formulae for the resonant modulation frequencies and transition rates, one would need to generalize the results of Section 3.2 for larger subspaces (more than four states in A gen n ). However, we can assure that these transitions do take place by substituting the formulae (11) and (15) into Equation (20) to obtain the (overly underestimated) 4-excitations transition rate 6,i for the transition |A 0 , 0 → |φ 6,i of the subspace A 6 . In Section 4.1, we shall calculate the transition rates for the transition |A 0 , 0 → |A 0 , n with n = 4 and 6 by exact numeric diagonalization of the Hamiltonian H 0 , showing that these transition rates are strongly enhanced in the vicinity of degeneracy between |A 0 , n and the states |A 2 , 1 or |A 3 , 0 . Since such multi-photon transitions are weaker than two-photons ones, we shall study their implementation in the ultrastrong coupling regime [33,55] with g ∼ 0.2-0.3ν.
Our scheme can be readily generalized to simultaneous modulation of other system parameters, one merely needs to add the corresponding matrix element on the RHS of Equation (18), which would produce an additive contribution to the transition rate (20). The inclusion of terms proportional to ε 2 in Equation (19) is also possible [74], but is not considered here because the resulting transition rates are roughly η/ε times smaller than Equation (20) (although in this case one benefits from halved resonant modulation frequencies).

Numeric Results and Discussion
To assess the feasibility of our scheme for photon generation from vacuum in circuit QED, we first assume conservative experimental parameters g = h = 0.05ν. In the following, we set the modulation amplitude as ε = 0.1Ω 0 .
In Figure 2a, we plot the dimensionless transition rate r ≡ |R 0,0;1,1 |/ν for transition from the system ground state |A 0 , 0 ≈ |g, g a , 0 to the state |A 1 , 1 ≈ −|g, e a , 1 as function of the t-qubit's bare frequency Ω 0 for the ancilla frequency Ω a = 0.6ν. Figure 2b shows the resonant modulation frequency η r = E 1,1 − E 0,0 (r and η r were obtained by diagonalizing numericallyĤ 0 ). This transition can be roughly interpreted as the Anti-Jaynes-Cummings ancilla-field regime in which one photon and one ancilla excitation are created, while the t-qubit remains approximately in the ground state. We verified that the relative errors between the exact numeric results and the analytic expressions of Section 3 is below 0.6% for the modulation frequency η r and below 5% for the transition rate r (data not shown). In Figure 2c, we show the fidelities F k,n (t) = | A k , n|ψ(t) | 2 of the dressed-states as function of time for (k = 0, n = 0) and (k = 1, n = 1), obtained by solving numerically the SE for the original Hamiltonian (1) with parameters Ω 0 = 0.95ν, η = 1.586ν, and the initial state |g, g a , 0 . For these parameters, the weights in the state |A 1 , 1 [Equation (15) with 2,i = 0.021 and φ (3) 2,i = 0.132 (so the contribution of the bare state |A 1 1 is the largest) and the transition rate is r = 1.93 × 10 −4 . As predicted in Section 3, there is a periodic population exchange between the dressed-states |A 0 , 0 and |A 1 , 1 with period π/(νr), while other states remain practically unpopulated. In Figure 2d, we present the numeric solution of the master Equation (3). We plot the average photon number n , the qubits excitations σ e and σ is also shown. This figure attests that photon generation from vacuum persists in the presence of dissipation. The order of magnitude of maximal allowed dissipation rates (denoted as γ max for the atoms and κ max for the cavity dampings) can be estimated from the condition γ max /ν, κ max n max /ν ∼ r (where |n max is highest Fock state significantly populated during the evolution, equal to n max = 1 in this example), which for r ∼ 10 −4 ν are ubiquitous in several circuit QED setups [71][72][73].

t-Qubit Frequency
This hint is confirmed in Figure 4c, where we solved numerically the SE and the master equation for parameters Ω a = 0.6ν, Ω 0 = 1.405ν and η = 2.0086ν, when the transition rates are approximately r 0 = 1.8 × 10 −4 , r 2 = 1.1 × 10 −4 and r 4 = 9.6 × 10 −5 . This plot displays the average photon number, the average total excitation number and the Mandel's Q-factor Q = (∆n) 2 / n − 1; bold (thin) lines depict the unitary (dissipative) evolution. Recall that Q quantifies the spread of the photon number distribution; some common values are Q = −1 (Fock state), Q = 0 (coherent state), and Q = n (thermal state). We see that a small number of photons can be created from vacuum even in the presence of dissipation, and the qubits remain approximately in the ground states because N tot is always close to n . The behavior of n is better understood by looking at Figure 4d, in which we plot the photon number probabilities of the Fock states, P n = Tr[ρ|n n|], with occupation probabilities above 10% under the unitary evolution. We see that up to six photons can be generated with significant probabilities, and the populations of the Fock states exhibit irregular oscillations due to the simultaneous coupling between the states |A 0 , 2n with n = 0, . . . , 3. The created field state is very different from the squeezed vacuum state generated in standard single-mode cavity DCE with vibrating walls or timedependent permittivity [1], for which Q = 1 + 2 n . As expected, the dissipation destroys the oscillating behavior of N tot , n , and Q after some time (r 0 νt π in Figure 4c) and these quantities tend to non-zero stationary values, loosing any resemblance with the unitary evolution.

Multi-Photon Transitions
For larger values of the cavity-ancilla coupling constant, g/ν ∼ 0.2-0.3 (ultrastrong coupling regime [33,55]), generation of photons from vacuum becomes feasible via direct driving of the ground state of the HamiltonianĤ 0 to excited eigenstates with n > 2 excitations. Here, we illustrate this rich variety of phenomena by considering the transitions from |A 0 , 0 to the states |A 0 , 4 and |A 0 , 6 . For the reasons already mentioned at the end of Section 3.2, it is more advantageous to obtain the transition rates and resonant modulation frequencies via numeric diagonalization ofĤ 0 . To assess the experimental feasibility of these phenomena in open quantum systems we solve numerically the SMME, but warn that the obtained results serve only to attest that the corresponding transition rates are sufficiently high to overcome the dissipation for initial times. Although the quantitative (and even qualitative) accuracy of SMME cannot be guaranteed in the ultrastrong coupling regime, the obtained results are useful for indicating that during the time interval 0 ≤ rνt π/2 the dissipative dynamics is qualitatively similar and not too different from the unitary one, so the generation of several photons should be possible.
In Figure 5a, we set the parameters Ω a = 0.6ν, g = 0.2ν, and h = 0.1ν, and evaluate numerically the dimensionless transition rate r = |R 0,0;0,4 |/ν for the transition |A 0 , 0 → |A 0 , 4 . As in previous figures, the transition rate presents sharp peaks near the degeneracies between |A 0 , 4 and the states |A 2 , 1 or |A 3 , 0 . The corresponding resonant modulation frequency η r = E 0,4 − E 0,0 is shown in Figure 5b. To confirm the feasibility of such multi-photon transition, we solved numerically the SE for Ω 0 = 3.12ν and η = 4.1873ν, when the transition rate is r ≈ 3.2 × 10 −4 (hence larger than the dissipation rates). For these parameters, the ground state is |A 0 , 0 ≈ 0.99|A 0 0 + 0.13|A 1 1 + 0.02|A 2 0 and the near degenerate dressed-states are |0, 4 ≈ 0.72|A 4 0 − 0.57|A 3 1 − 0.27|A 1 2 + 0.25|A 5 1 + 0.12|A 0 3 (with approximate energy 4.1868ν above the ground state energy) and |2, 1 ≈ 0.86|A 1 2 − 0.39|A 0 3 + 0.22|A 4 0 − 0.18|A 3 1 + 0.16|A 2 3 (with the corresponding energy 4.1899ν). Figure 5c shows the behavior of n , N tot and Q under unitary evolution, as well as the total excitation number N (dis) tot assuming the standard dissipation kernel (6). We see that for rνt 2, N tot , and N (dis) tot have similar behavior, confirming that the transition rate is sufficiently high to overcome the dissipation for initial times. Figure 5d shows the excitation probability of t-qubit and ancilla for the initial zeroexcitation state |g, g a , 0 , together with the fidelities F k,l = | A k , l|ψ(t) | 2 during the unitary evolution. This intricate behavior occurs because the modulation couples the ground state |A 0 , 0 to both the states |A 0 , 4 and |A 2 , 1 , and we verified that F 0,0 + F 0,4 + F 2,1 > 0.98, proving that only the states |A 0 , 0 , |A 0 , 4 , and |A 2 , 1 become significantly populated throughout the evolution. We see that nearly four excitations are created during the expected time interval rνt = π/2, and both qubits can be found in excited states with significant probabilities. The excitation of the ancilla comes mainly from the contribution of the state |A 3 1 , and the excitation of the t-qubit comes mainly from the contribution of |A 1 2 in the dressed-states (since |A 1 ≈ −|g, e a and |A 2 ≈ |e, g a for the chosen parameters). Up to five photons can be detected with probabilities above 5%, and the contribution of the vacuum state becomes nearly zero for rνt ≈ π/2, as illustrated in Figure 5e. Figure 6 is analogous to Figure 5 but for the direct transition |A 0 , 0 → |A 0 , 6 . We performed numerical simulations for larger coupling strengths, g = 0.3ν and h = 0.2ν, but the same ancilla frequency Ω a = 0.6ν as in Figures 4 and 5. In panel Figure 6a, we verify a strong enhancement of the transition rate r = |R 0,0;0,6 |/ν in the vicinity of degeneracy between |A 0 , 6 and the states |A 2 , 1 or |A 3 , 0 . The resonant modulation frequency η r = E 0,6 − E 0,0 is shown in Figure 6b. In panels Figure 6c-e, we solved numerically the SE for the initial state |g, g a , 0 and parameters Ω 0 = 4.057ν and η = 5.201ν, when the transition rate is r = 4.6 × 10 −4 . The near degenerate states are |A 0 , 6 ≈ 0.68|A 6 0 + 0.52|A 5 1 − 0.39|A 4 0 + 0.25|A 7 1 + 0.16|A 3 1 − 0.14|A 1 2 (with approximate energy 5.2025ν above the ground state energy) and |A 2 , 1 ≈ 0.8|A 1 2 − 0.52|A 0 3 + 0.25|A 2 3 + 0.11|A 6 0 (with the corresponding energy 5.1978ν), while the ground state is |A 0 , 0 ≈ 0.98|A 0 0 + 0.19|A 1 1 + 0.04|A 2 0 . Figure 6c confirms the generation of photons with n ≤ 4.5 and N tot ≤ 5 during the unitary evolution. Under dissipation, the plot of N (dis) tot shows that for rνt 1 the influence of dissipation is small, and generation of a few photons is feasible. The populations of Fock states with occupation probabilities above 10% are shown in Figure 6e (for unitary evolution). As expected, up to six photons are generated with significant probabilities (while the Fock states |2 and |3 are practically unpopulated), and the substantial population of the one-photon state comes from the partial excitation of the dressed-state |A 2 , 1 . Figure 6d shows the fidelities of the dressed states |A 0 , 6 , |A 2 , 1 and |A 0 , 0 , whose sum is always above 96%, attesting that only these states become significantly populated during the evolution.
To conclude, we emphasize that in the regimes studied above the ancilla-cavity counter-rotating term g(âσ (a) − + h.c.) in the HamiltonianĤ plays an essential role. In the two-excitations transitions studied in Figures 2-4, its major effect is to alter the resonant modulation frequencies and the position of peaks in transition rates, so in the first approximation it could be neglected via the Rotating Wave Approximation (RWA) (the transition rate would be wrong by roughly 30% in this case, provided one adjusted by hand the positions of the peaks). However, for 4-and 6-excitations transitions studied in Figures 5 and 6 that term is indispensable, and the transition rates would be orders of magnitude smaller if it were neglected. This is illustrated in Figure 7, where we plot the transition rates of Figures 4a, 5a and 6a (corresponding to the direct transition |A 0 , 0 → |A 0 , n , n = 2, 4, 6) with and without the counter-rotating term and indicate the origin of each peak (degeneracy of |A 0 , n with |A 2 , 1 or |A 3 , 0 ).

Conclusions
We showed that photons can be generated from a vacuum inside a stationary cavity due to the resonant time-modulation of a quantum system that is indirectly coupled to the field via some ancilla quantum subsystem. Such coupling may be advantageous when the modulated system is located outside or at the end of the resonator (when its coupling to the field is zero) to reduce the influence of the external driving on the cavity field, while it is strongly coupled to one or several ancilla subsystems. Moreover, this setup permits fabrication of quantum systems designed specifically for rapid external modulation at the cost of increased dissipation, while the stationary ancillas and cavities can be optimized for minimal dissipative losses.
We considered the simplest scenario of a single-mode cavity and a time-modulated qubit (we named "t-qubit") playing the role of the indirectly coupled system, while the ancilla consisted of a stationary qubit dipole-coupled to both the cavity and the t-qubit. Our scheme and mathematical analysis can be readily generalized to more complex systems with alternative coupling mechanisms. We deduced analytically the system dynamics under unitary evolution and showed that the rate of photon generation can be drastically enhanced near certain resonances of the tripartite system (related to anti-crossings in the energy spectrum). We exemplified our scheme by studying photon generation from the initial vacuum state of the system, considering regimes in which either a single dressed-state or a small number of dressed-states with a few photons are populated under sinusoidal modulation. Moreover, we demonstrated that these phenomena can withstand dissipation with damping rates of the order of 10 −5 -10 −4 ν, since the typical photon generation rates are of the order of 10 −4 ν (where ν is the cavity frequency). Although the present study focused on prospects of photon creation from vacuum, our scheme could find applications in the engineering of effective interactions and generation of useful multipartite entangled states and non-classical states of light.