Off-resonant Dicke Quantum Battery: Charging by Virtual Photons

TL;DR

This paper explores an off-resonant Dicke quantum battery where virtual photons enable charging through matter-radiation coupling modulation, revealing new charging mechanisms and performance scaling in the dispersive regime.

Contribution

It introduces the concept of charging a quantum battery via virtual photons in the off-resonant regime, expanding understanding of quantum energy storage mechanisms.

Findings

Charging occurs even with an empty cavity due to virtual photon interactions.

Performance scales favorably with the number of two-level systems in strong coupling.

Efficient energy extraction is possible across various parameter ranges.

Abstract

We investigate a Dicke quantum battery in the dispersive regime, where the photons trapped into a resonant cavity are way more energetic with respect to the two-level systems embedded into it. Under such off-resonant conditions, even an empty cavity can lead to the charging of the quantum battery through a proper modulation of the matter-radiation coupling. This counterintuitive behaviour has its roots in the effective interaction between two-level systems mediated by virtual photons emerging from the fluctuations of the quantum electromagnetic field. In order to properly characterize it, we address relevant figures of merit such as the stored energy, the time required to reach the maximum charging, and the averaged charging power. Moreover, the possibility of efficiently extracting energy in various ranges of parameters is discussed. The scaling of stored energy and power as a function of the number $N$ of two-level systems and for different values of the matter-radiation coupling is also discussed showing, in the strong coupling regime, performances in line with what reported for the Dicke quantum battery in the resonant regime.

Tables

Table 1. Table 1: Scaling behaviour for some figures of merit in the weak and strong coupling regimes.

Coupling regime	$E_{max}$	$P_{max}$	$t_E,t_P$
Weak coupling	$N^2$	$N^2$	Constant
Strong coupling	$N$	$N^{\frac{3}{2}}$	$N^{-\frac{1}{2}}$

Equations

$${\hat{H}}_{\mathrm{Dicke}}$$

$$=$$

$${\hat{S}}_{\alpha }=\frac{1}{2}\sum _{i=1}^N{\hat{\sigma }}_{\alpha }^i,$$

$${\hat{H}}_{\mathrm{eff}}={\omega }_z{\hat{S}}_z+{\omega }_c{\hat{a}}^{†}\hat{a}-\frac{2{\lambda }^2{\omega }_z}{{\omega }_c^2-{\omega }_z^2}{\hat{S}}_z{\left({\hat{a}}^{†}+\hat{a}\right)}^2-\frac{4{\lambda }^2{\omega }_c}{{\omega }_c^2-{\omega }_z^2}{\hat{S}}_x^2,$$

$${\hat{H}}_{\mathrm{eff}}^{\prime }={\omega }_z{\hat{S}}_z+{\omega }_c{\hat{a}}^{†}\hat{a}-\frac{4{\lambda }^2}{{\omega }_c}{\hat{S}}_x^2.$$

$$\hat{H}$$

$$=$$

$$g=\frac{4{\lambda }^2}{{\omega }_c}.$$

$$\hat{H}(t)={\hat{H}}_B+f(t){\hat{H}}_{\mathrm{int}},$$

$${\hat{S}}^2={\hat{S}}_x^2+{\hat{S}}_y^2+{\hat{S}}_z^2$$

$${\hat{S}}^2|s,m⟩$$

$${\hat{S}}_z|s,m⟩$$

$$\hat{\Pi }=e^{-i\pi {\hat{S}}_z}$$

$$|\psi (0)⟩=|\frac{N}{2},-\frac{N}{2}⟩.$$

$$|{\psi }_{max}⟩=|\frac{N}{2},+\frac{N}{2}⟩$$

$$E_B(t)=⟨\psi (t)|{\hat{H}}_B|\psi (t)⟩-⟨\psi (0)|{\hat{H}}_B|\psi (0)⟩,$$

$$P(t)=\frac{E_B(t)}{t}.$$

$$\mathcal{E}\left(\rho (t),\hat{\mathcal{H}}\right)=\mathrm{tr}[\hat{\mathcal{H}}\hat{\rho }(t)]-\mathrm{tr}[\hat{\mathcal{H}}{\hat{\pi }}_{\rho }(t)],$$

$$\hat{\rho }(t)=\sum _n{r}_n(t)|r_n(t)⟩⟨r_n(t)|,$$

$${\hat{\pi }}_{\rho }(t)=\sum _n{r}_n(t)|{ϵ}_n⟩⟨{ϵ}_n|,$$

$$\mathcal{E}\left(\hat{\rho }(t),\hat{\mathcal{H}}\right)=\sum _{k,n}{r}_k(t){ϵ}_n\left(|⟨r_k(t)|{ϵ}_n⟩{|}^2-{\delta }_{k,n}\right),$$

$$\hat{\rho }(t)=|\psi (t)⟩⟨\psi (t)|$$

$${\hat{h}}_B^1=\frac{{\omega }_z}{2}{\hat{\sigma }}_z^1.$$

$${\mathcal{E}}_1^{(N)}(t)\equiv \mathcal{E}({\hat{\rho }}_{B,1}(t),{\hat{h}}_B^1).$$

$${\hat{\rho }}_{B,1}(t)=r_0(t)|r_0(t)⟩⟨r_0(t)|+r_1(t)|r_1(t)⟩⟨r_1(t)|,$$

$${\mathcal{E}}_1^{(N)}(t)=\frac{E^{(N)}(t)}{N}-r_1(t){\omega }_z,$$

$$\mathcal{E}(t)\ne N{\mathcal{E}}_1^{(N)}(t).$$

$$E_B(t)\approx \frac{g^2}{8}\frac{1-\cos 2{\omega }_{z}t}{{\omega }_z}\left(N^2-N\right),$$

$$P(t)\approx \frac{g^2}{8}\frac{1-\cos 2{\omega }_{z}t}{{\omega }_{z}t}\left(N^2-N\right).$$

$$E_{max}$$

$${\omega }_z{t}_E$$

$$P_{max}$$

Full text

\pubvolume

1 \issuenum1 \articlenumber0

\externaleditorAcademic Editor: \datereceived \dateaccepted \datepublished \hreflinkhttps://doi.org/

\TitleOff-resonant Dicke Quantum Battery: Charging by Virtual Photons \TitleCitationOff-resonant Dicke Quantum Battery: Charging by Virtual Photons

\Author Giulia Gemme1, Gian Marcello Andolina2, Francesco Maria Dimitri Pellegrino3,4,5, Maura Sassetti1,6,∗ and Dario Ferraro${}^{1,6,\orcidA{}}$

\AuthorNamesGiulia Gemme, Gian Marcello Andolina, Francesco Maria Dimitri Pellegrino, Maura Sassetti and Dario Ferraro \AuthorCitationGemme, G.; Andolina, G.M.; Pellegrino, F.M.D.; Sassetti, M.; Ferraro, D.

\corresCorrespondence: sassetti@fisica.unige.it

1 Introduction

Quantum batteries (QBs) Alicki (2013) —genuinely non-classical devices capable of storing energy and performing useful work— have been a very active topic of research in recent years Campaioli (2018); Bhattacharjee (2021). In this framework, theoretical studies have demonstrated that entangling operations can consistently accelerate the charging Binder (2015), with the Dicke QB standing out as a particularly promising candidate Ferraro18 ; Crescente20 ; Dou (2022); Dou22 ; Erdman22 . It is given by a system where the energy of photons trapped into a cavity (acting as the charger) is coherently transferred to a QB composed of $N$ quantum units, namely a collection of identical independent two-level systems (TLSs). This model describes, for example, an ensemble of real or artificial atoms embedded in a Fabry-Perot cavity, namely a resonator consisting of two parallel, highly reflective mirrors placed a certain distance apart Dicke54 .

In this system, the energy is initially stored into the cavity and flows towards the TLSs due to a proper modulation of the matter-radiation coupling Andolina (2018). In order to favour this energy flow the cavity and the TLSs are usually assumed as resonant, namely characterized by the same typical energy. Once the charging process is complete, the coupling between the two systems is turned off, leaving the TLSs in a charged state whose stability crucially depends on their relaxation and dephasing times Farina19 ; Santos19 ; Barra19 ; Carrega20 ; Morrone22 ; Gemme22 . This kind of device displays a collective speed-up of the charging process. This fact has attracted a great deal of interest and, very recently, a first step towards the realization of a Dicke QB has been experimentally implemented in an excitonic system made of fluorescent organic molecules dispersed in an inert matrix, where the Dicke superabsorption leads to a collective boost of the charging process Quanch (2022). Despite this remarkable achievement, the investigation of this kind of device is still in its infancy and new theoretical and experimental studies could lead to promising routes towards the exploitation of purely quantum effect in miniaturized devices devoted to energy storage.

In this direction we want to investigate the role played by exchange of virtual photons between the TLSs composing the QB. A paradigmatic example of the dynamical effect of virtual photons is the Casimir force which arises from the quantum fluctuations of the electromagnetic field. It was first theoretically investigated by Hendrik Casimir in 1948 Casimir48 and experimentally proved almost fifty years later Lamoureaux97 . This force is due to the presence of virtual photons in the vacuum and, in the case of a Fabry-Perot cavity, depends on the boundary conditions imposed by the mirrors. In particular, in the case of two parallel, uncharged, and perfect metallic plates, the virtual cavity photons are restricted in their fluctuation modes by the plates themselves. As a result, the plates experience a net attractive force.

Motivated by this, we introduce and characterize a new type of device for energy storage where the charging is mediated by virtual photons. In contrast to the conventional Dicke QB, the cavity in our proposed off-resonant Dicke QB can be also empty and the two systems (charger and QB) are far detuned in energy, preventing any direct energy flow between them. This makes the systems easier to realize with respect to the conventional Dicke QB and allows us to describe it in the so-called dispersive regime Schleich (2021), typically used in quantum computing for qubit readout Krantz (2019) and characterized by an effective infinite range interaction among the TLSs of the Lipkin-Meshkov-Glick (LMG) kind Lipkin65 ; Dou22b . Our charging protocol works as follows: initially, both the cavity and the QB are in their ground states and decoupled. The cavity mirrors are then brought close together, leading to a finite matter-radiation coupling. Thus, the virtual detuned photons excite the TLSs and charge the QB. The mirrors are then moved away, switching-off again the matter-radiation coupling leaving the atoms in an excited state. The use of virtual photons in this process gives the QB a unique quantum behaviour. Remarkably enough, in the strong coupling regime, the performances of these off-resonant Dicke quantum batteries are comparable with the ones of the resonant case.

The article is organized as follows. In Section 2 we consider the dispersive regime of the conventional Dicke QB with dipolar matter-radiation coupling. Under the assumption of highly energetic cavity photons, we derive an effective description of the off-resonant Dicke QB based on an LMG model characterized by an infinite range interaction among the TLSs and showing relevant simplifying conservation laws. The charging protocol based on switching on and off the matter-radiation coupling is also discussed. The main figures of merit, namely the energy stored into the QB, the averaged charging power, and the times required to reach their maxima are discussed in Section 3. The ergotropy, namely the maximum energy which can be extracted by means of either collective or individual unitary operation is also introduced. Section 4 reports the results concerning the time evolution of the above quantities for different values of the matter-radiation coupling and at various numbers $N$ of TLSs. The scaling of their maxima as a function of $N$ is also considered. A universal behaviour of the maximum stored energy is discussed in connection with the quantum phase transition predicted for the underlying effective model in Section 5. Section 6 is devoted to the conclusions, while technical details of the calculations are reported in three Appendices.

2 Model

Let’s start by considering a conventional Dicke-QB with a dipole matter-radiation coupling between $N$ identical TLSs and a resonant cavity Ferraro18 . It is described by the Hamiltonian (from now on we will set $\hbar=1$)

[TABLE]

with $\hat{H}_{B}$ the Hamiltonian of the QB, $\hat{H}_{C}$ the one of the cavity and $\hat{H}_{B-C}$ their coupling. In the above equation, we have introduced the collective pseudo-spin operator

[TABLE]

with $\hat{\sigma}^{i}_{\alpha}$ ($\alpha=x,y,z$) the conventional Pauli matrices associated with the $i$-th TLS and where $\hat{a}$ ($\hat{a}^{\dagger}$) are the annihilation (creation) operator for the photons in the cavity. Here, $\omega_{z}$ is the level spacing of each TLS, $\omega_{c}$ is the energy of the photons, and $\lambda$ is the intensity of the matter-radiation coupling. Notice that this configuration can be effectively engineered in conventional cavity Walther06 and circuit-quantum electrodynamics Blais04 set-ups. In the latter case, the TLSs can be realized by means of superconducting qubits coupled to an LC circuit playing the role of a resonant cavity Krantz (2019).

According to the analysis carried out in Appendix A, considering the Schrieffer–Wolff transformation Schleich (2021) up to the second order in $\lambda$ one obtains the effective Hamiltonian

[TABLE]

which is a good approximation of the Dicke Hamiltonian in Eq. (2) in the dispersive regime $\lambda\ll|\omega_{c}-\omega_{z}|$ Schleich (2021). Notice that, in this regime, the possibility to increase the energy stored into a QB through virtual excitations mediated by a higher excited state has been recently discussed Santos23 .

Further considering the condition of highly energetic photons $\omega_{c}\gg\omega_{z}$ one finally obtains Roche22

[TABLE]

Getting rid of the cavity contribution it is possible to map Eq. (5) into

[TABLE]

with

[TABLE]

Notice that, according to this mapping, only positive values of the parameter $g$ are physically meaningful. Moreover, as shown in Appendix A, even a more exotic two-photon matter-radiation coupling Felicetti15 ; Felicetti18 leads to an analogous effective Hamiltonian.

The Hamiltonian in Eq. (7) is characterized by an effective infinite range interaction between the TLSs of the LMG kind Lipkin65 ; Dou22b ; Abah22 mediated by virtual photons playing the role of a quantum bus Larson10 . Such kind of infinite range interaction has been indicated as an ideal condition to achieve high power QBs Campaioli (2017); Le (2018); Gyhm (2022).

In the following, we will assume a charging protocol addressing a time-dependent version of Eq. (7)

[TABLE]

where $f(t)$ is a classical external control parameter introduced to turn on and off the effective interaction among the TLSs in Eq. (7). This can be ultimately seen as a direct consequence of the modulation of the matter-radiation coupling in the original Dicke model in Eq. (2) obtained by bringing the cavity mirrors close together. For sake of simplicity, we will assume this control function to be step-wise: equal to $1$ for $t\in[0,\tau_{c}]$, with $\tau_{c}$ a given time which will be specified in the following, and zero elsewhere. It is worth noting that smoother charging protocols are expected to lead only to marginal changes in the behaviour discussed in the following Crescente22 . Generally speaking, the possibility to optimize the time evolution of the charging protocol has been recently object of intense research for both the Dicke QB Erdman22 and other QB technologies Hu (2021); Rodriguez23 . Moreover, we further underline the fact that differently from what is usually considered Campaioli (2018); Bhattacharjee (2021), and in the same spirit of Refs. Rosa (2020); Rossini (2020), here the charger is encoded in the effective interaction among the TLSs composing the QB and not into some external classical or quantum device.

2.1 Conserved quantities

In order to investigate the behaviour of QB, for a fixed number $N$ of TLSs, a full numerical exact diagonalization is needed. In order to consistently simplify this task some symmetries can be exploited.

First of all, the operator

[TABLE]

is a conserved quantity, namely $\left[\hat{H}(t),\hat{S}^{2}\right]=0$. This allows us to work in the $(N+1)\times(N+1)$ Hilbert space with maximal value $s=N/2$ and to characterize the states of the system in terms of a pseudo-spin basis $|s,m\rangle$ such that:

[TABLE]

with $m=-s,-s+1,\dots,s-1,s$.

Moreover, also the spin number parity operator

[TABLE]

satisfies $\left[\hat{H},\hat{\Pi}\right]=0$, leading to a further reduction of the Hilbert space dimension for proper initial conditions of the QB. According to this, in the following, we will investigate the charging of the QB starting from the ground state of $\hat{H}_{B}$

[TABLE]

It is worth noticing that, according to this parity constraint, the perfect charging of the QB, namely the complete transition to the maximally charged state

[TABLE]

can only occur for an even number $N$ of TLSs. However, as will be clear in the following, the efficiency of this kind of QB is lower than $100\%$, therefore we can investigate equivalently both even and odd numbers of TLSs.

3 Figures of merit

We introduce now relevant quantities able to characterize the performances of the considered off-resonant Dicke QB (with $\omega_{c}\gg\omega_{z}$). In the following, we will focus on the energy stored into the QB itself, the time required to reach the optimal charging, the energy trapped into the interaction term, the averaged charging power, and the ergotropy, namely the maximum amount of work extractable from a charged QB using either collective or local unitary operations Alicki (2013); Allahverdyan04 .

3.1 Stored energy and charging time

First of all, we consider the energy stored into the QB. At a given time $t$ it is given by Ferraro18 ; Andolina (2018)

[TABLE]

with $|\psi(t)\rangle$ the quantum state of the systems for $t\leq\tau_{c}$, evolved according to $\hat{H}(t)$.

The charging time, namely the time $t_{E}$ at which the quantity in Eq. (16) reaches is first maximum $E_{max}$, will be also investigated in the following. Of particular relevance will be the scaling of both $E_{max}$ and $t_{E}$ as a function of the number $N$ of TLSs composing the QB for different values of the effective coupling $g$.

It is worth underlining the fact that the system under investigation is considered as closed. Under this approximation, which holds as long as the time evolution is shorter than the typical time scale associated with relaxation and dephasing of the TLSs Dou (2022); Carrega20 , one expects a net energy is injected into the QB via the modulation of the external control $f(t)$.

3.2 Averaged charging power

Another relevant figure of merit is the averaged charging power, namely the energy stored in a given time, defined as Binder (2015); Campaioli (2017)

[TABLE]

In analogy with what was done above, we are interested in the scaling of the maximum charging power $P_{max}$ and of the time $t_{P}$ requited to reach it as a function of $N$ and for various values of the effective coupling $g$.

3.3 Ergotropy

Here we introduce the ergotropy, namely the maximum work extractable from a QB Allahverdyan04 . We will focus both on the total ergotropy and the one associated with only one of the TLS composing the QB. On a very general ground at a given time $t$, for an Hamiltonian $\hat{\mathcal{H}}$ and a time dependent density matrix $\hat{\rho}(t)$ describing a quantum state at a given time, it is possible to define the quantity

[TABLE]

with $\hat{\pi}_{\rho}(t)$ the passive state associated to $\hat{\rho}(t)$ at a given time. Given the diagonal representation of the density matrix

[TABLE]

with $r_{0}(t)\geq r_{1}(t)\geq\cdots$, the associated passive state is diagonal on the eigenbasis of $\hat{\mathcal{H}}$ ($|\epsilon_{n}\rangle$) and reads

[TABLE]

with energy eigestates ordered in such a way that $\epsilon_{0}\leq\epsilon_{1}\leq\cdots$. This finally leads to the general definition

[TABLE]

which can be directly applied to our particular case considering

[TABLE]

and $\hat{\mathcal{H}}=\hat{H}_{B}$.

Due to the fact that we are considering only unitary operations acting on a closed system, the ergotropy defined above coincides with the energy stored into the QB Alicki (2013).

For what it concerns the ergotropy of a single TLS of the QB in Eq.(2) one can proceed as follows. The $2\times 2$ reduced density matrix describing the first TLS at a given instant is indicated by $\hat{\rho}_{{B},1}$, while the energy is measured with respect to the Hamiltonian

[TABLE]

Notice that, due to the invariance under TLSs permutations of the Hamiltonian in Eq. (7), we can focus on the first TLS without any loss of generality. Thus, the maximum energy that can be extracted from a single battery unit in an $N$ TLSs device is given by

[TABLE]

Consistently with what is done above, this expression can be further simplified by writing $\hat{\rho}_{{B},1}(t)$ in the diagonal basis,

[TABLE]

where the eigenvalues are again ordered in such a way that $r_{0}(t)\geq r_{1}(t)$. In this case, the ergotropy $\mathcal{E}_{1}^{(N)}(t)$ reduces to

[TABLE]

where we used the fact that ${\rm tr}[\hat{\rho}_{{B},1}(t)\hat{h}_{1}^{B}]=({E^{(N)}(t)}/{N})$, with $E^{(N)}={\rm tr}[\hat{H}_{B}\hat{\rho}(t)]$, due to permutation symmetry. We have evaluated the above expression numerically, further details about this analysis can be found in Ref. Erdman22 .

We conclude this Section by pointing out that for a fixed number $N$ of TLSs, in general,

[TABLE]

This is a consequence of the fact that local unitary operations acting on a unique TLS are not enough in order to extract all the energy stored into the QB, and non-local unitary operations collectively acting on the whole QB are needed in order to extract also the energy trapped in the correlations among the elementary building blocks composing the device Alicki (2013).

4 Results and scaling laws

In the following, we will investigate in detail the behaviour of the previously discussed figures of merit for different values of the coupling $g$ and as a function of the number $N$ of TLSs composing the off-resonant Dicke QB.

4.1 Weak coupling

Let’s start by considering the regime $g\ll\omega_{z}$. The time evolution of the energy stored in the QB and the averaged charging power are reported in Fig. 1.

This behavior can be understood considering the analytical expressions for both the above quantities using a time dependent perturbative approach (see Appendix B for more details about the derivation). This leads to

[TABLE]

and

[TABLE]

These approximated expressions correspond to the black dashed curves reported in Fig. 1. For the reported coupling $g=10^{-3}\omega_{z}$ they show an excellent agreement with the numerical curves up to $N\approx 30$, further validating the correctness of the considered exact diagonalization and supporting the observation that there is only one relevant frequency of the state dynamics in this regime, leading to very smooth curves with no beats. It is worth noting that the amount of energy stored in this regime is extremely small (see the scale in Fig. 1). As a consequence of this, the ergotropy for a single TLS is zero (dashed orange curve in Fig. 1), namely no energy can be extracted for a single TLS using local unitary operations.

Starting from Eq. (28) one can obtain the values of the maximal stored energy and of the charging time associated with the first maximum, namely

[TABLE]

respectively. Analogously for the averaged charging power in Eq. (29) one has

[TABLE]

where the values in this case need to be extracted numerically. The validity of the scalings of stored energy and charging power as a function of the number of TLSs has been checked in Fig. 2. In this framework, we underline the fact that throughout all the paper we will consider the scaling at a finite number of TLSs, without addressing the thermodynamical limit, differently from what is discussed for example in Ref. JuliaFarre (2020). This is justified by the fact that, for example, circuit quantum-electrodynamics devices are usually far from the thermodynamic limit Fink09 .

4.2 Strong coupling

We investigate here the strong coupling regime, namely a situation in which $g$ is of the same order of $\omega_{z}$ Giannelli22 . Even if difficult to be realized in the dissipative Dicke model with dipolar coupling, this situation could be reached by exploiting more exotic matter-radiation couplings such as the two-photon interaction and properly changing the initial photon distribution in the cavity (see Appendix A for more details). The time evolution of the figures of merits, reported in Fig. 3, appears again quite regular at small $N$ but shows a complicated beating structure at greater values of $N$. Notice that, in this regime, the first maximum of the energy stored into the QB is only a local maximum, with higher values of the energy possibly appearing at greater times. However, we focus on it due to the fact that it is characterized by a greater value of the power (see Fig. 3 (d)). In this regime, the charging of the QB can exceed the $50\%$ near the maxima. In correspondence to these points, one observes the emergence of a non-zero single TLS ergotropy (see the dashed orange curve in Fig. 3 (a) and (b)), the signature of the possibility to extract energy from the QB also by using local unitary operations.

According to the full numerical analysis reported in Fig. 4 it is possible to deduce the scaling behaviour of the various figures of merit in the strong coupling regime, namely

[TABLE]

Notice that these are reminiscent of the ones reported for the Dicke model with conventional dipole matter-radiation coupling Ferraro18 . This can be justified by considering the fact that in this case the expectation value of the operator $\hat{S}_{x}$ over both the ground and the excited states of the system becomes extensively large, similarly to what happens for both $\hat{S}_{x}$ and the operator $(\hat{a}^{\dagger}+\hat{a})$ in the original Dicke model in Eq. (2) after the superradiant quantum phase transition Emary03 .

5 Considerations about universality

The previous analysis showed two drastically different scaling behaviours of the considered quantities for different regimes of the coupling, with the weak coupling showing a better scaling with $N$ and the strong coupling regime more promising in terms of absolute values of the energy storage, charging power, and ergotropy. These are summarized in Table 1. Notice that for intermediate values of the coupling, a crossover between the two behaviour emerges (not shown).

This change of behaviours can be seen as a consequence of the Excited State Quantum Phase Transition (ESQPT) shown by both the Dicke and the effective LMG model (see Ref. Cejnar21 for a review on the subject). This represents the footprint of a quantum phase transition at the level of the excited states of the system. Its role in the present QB is justified by the fact that it is initialized in the ground state and charged by means of a protocol which naturally involves the excited states. In order to better clarify this point, we can consider the evolution of $E_{max}$ as a function of the renormalized coupling $G=gN$ (see Fig. 5). The choice of this renormalized parameter to characterize the evolution of the system can be justified by considering the classical limit of the effective LMG model (see Appendix C for more details) and reminiscent of the renormalization $\Lambda=\lambda\sqrt{N}$ observed for the Dicke model Ferraro18 . Notice that the emergence of a crossover behaviour for the maximum stored energy is motivated by the fact that it is closely related to the order parameter of the system, namely the total magnetization, which at a given time $t$ can be written as

[TABLE]

It is possible to observe now two different regimes for $G<\omega_{z}$ and $G>\omega_{z}$ respectively, with a crossover in correspondence of the ESQPT of the model, occurring at $G=\omega_{z}$. Incidentally, by increasing the value of $G$ (Fig. 5 (b)) one has that the maximum value of energy per TLS reachable in the considered QB can exceed, the $60\%$. Remarkably enough, this value is compatible with what is observed for the Dicke model in the strong coupling regime and under resonant conditions Ferraro18 . This makes the present off-resonant Dicke QB competitive for practical experimental realizations, also taking into account the fact that it can be realized also using a vacuum cavity which is easier to be engineered and more stable with respect to what happens in the resonant Dicke QB. Moreover, the dependence on the renormalized coupling $G$ allows to extend the validity of the previous analysis carried out for the strong coupling regime by increasing the number of TLS composing the QB.

6 Conclusions

We have investigated an off-resonant Dicke QB, namely a device composed by $N$ TLSs embedded into a highly detuned resonant cavity. Under this condition, the TLSs feel an effective infinite range interaction among them. By properly switching on and off this coupling it is possible to promote at least a part of the TLSs from the ground to the excited states, thus charging the QB, even in the case of an empty cavity. This phenomenology can be ultimately related to the exchange of virtual photons whose coupling with the TLSs is varied bringing the mirrors of the cavity close together.

By looking at relevant figures of merits such as the stored energy, the averaged charging power, and the time required to reach their respective maxima, we have characterized the performance of this kind of system. We have identified a small coupling regime, showing a very poor charging and no possibility of energy extraction (single TLS ergotropy) despite the relevant scaling in both the energy and the averaged charging power, and a strong coupling limit showing the same collective advantage in averaged charging power of the conventional Dicke QB Ferraro18 with dipole matter-radiation interaction and in resonant conditions. The crossover between these two regimes is related to the quantum phase transition occurring in the effective model describing the QB Cejnar21 .

The present device shows, above the critical coupling, performances comparable to the ones of the Dicke QB in the resonant regime. In addition, the observed physics can be obtained also working with an empty cavity, thus avoiding the need of controlling the initial quantum state of radiation, the physically relevant parameters of the system, or exotic matter-radiation interactions Crescente20 . These remarkable properties make it a very convincing candidate for future practical implementations of a functioning QB in various platforms already considered in the framework of the Dicke QB, such as superconducting qubits or array of semiconducting quantum dots coupled to an LC circuit through a tunable capacitance Xiang13 ; Stockklauser17 .

\authorcontributions

Conceptualization, G. M. A and D. F.; methodology, G. M. A. and D. F.; software, G. G. and G. M. A.; validation, F. M. D. P. and G. M. A.; formal analysis, D. F.; investigation, G. M. A. and G. G.; data curation, G. G.; writing—original draft preparation, D. F. and G. G.; writing—review and editing, G. M. A., F. M. D. P. and M. S.; supervision, M. S. and All authors have read and agreed to the published version of the manuscript.

\funding

D. F. would like to thank the funding of the European Union-NextGenerationEU through the "Quantum Busses for Coherent Energy Transfer" (QUBERT) project, in the framework of Curiosity Driven 2021 initiative of the University of Genova. F.M.D.P. was supported by the Università degli Studi di Catania, Piano di Incentivi per la Ricerca di Ateneo 2020/2022 (progetto QUAPHENE), and PNRR MUR project PE0000023-NQSTI.

\institutionalreview

Not applicable.

\informedconsent

Not applicable.

\dataavailability

Data are available from the authors upon request.

Acknowledgements.

The authors thank A. Crescente, G. Falci, L. Giannelli, and E. Paladino for illuminating discussions and fruitful comments. \conflictsofinterestThe authors declare no conflict of interest. \abbreviationsAbbreviations The following abbreviations are used in this manuscript:

[TABLE]

\abbreviationsNomenclature List of the Relevant Parameters in the Model:

\appendixtitlesyes \appendixstart

Appendix A Derivation of the effective LMG Hamiltonian

A.1 Basics of the Schrieffer–Wolff transformation

Given a generic Hamiltonian $\hat{\mathcal{H}}$ and an operator $\hat{\mathcal{S}}$ it is possible to define

[TABLE]

According to the Baker-Campbell-Hausdorff relation, one can develop the exponential operators in such a way that

[TABLE]

On a general ground, one can now assume that it is possible to decompose $\hat{\mathcal{H}}$ in a free and an interacting term, namely

[TABLE]

By properly choosing the operator $\hat{\mathcal{S}}$ in such a way to satisfy

[TABLE]

one obtains

[TABLE]

Assuming now that both $\hat{\mathcal{S}}$ and $\hat{\mathcal{V}}$ depend linearly on a given coupling constant one can consider a second order expansion leading to an effective Hamiltonian of the form

[TABLE]

In the following, we will apply this general scheme to the Dicke model assuming both single- and two-photon interaction.

A.2 Dicke model with single-photon coupling

In the case of the Dicke model with conventional dipole matter-radiation coupling Dicke54 ; Schleich (2021) in Eq. (2) one can identify (using the notation of the main text)

[TABLE]

and consequently

[TABLE]

with

[TABLE]

Notice that both Eq. (46) and (47) are linear in the matter-radiation coupling $\lambda$ as required by the previous discussion. Using the Schrieffer–Wolff transformation discussed above one obtains the effective Hamiltonian

[TABLE]

which is a good approximation of the Dicke Hamiltonian in the so called dispersive regime $\lambda\ll|\omega_{c}-\omega_{z}|$ Schleich (2021). Further considering the limit $\omega_{c}\gg\omega_{z}$ one finally has

[TABLE]

Getting now rid of the cavity contribution it is possible to map Eq. (50) into Eq. (7) by defining

[TABLE]

A.3 Dicke model with two-photon coupling

This more exotic interaction can be obtained by properly engineering both atomic and solid state systems Felicetti15 ; Felicetti18 . Here, one can write the Hamiltonian

[TABLE]

identifying

[TABLE]

and consequently

[TABLE]

According to the previous analysis, one obtains the effective Hamiltonian

[TABLE]

which again represents a good approximation of the two-photon Dicke Hamiltonian in the dispersive limit $\lambda\ll|\omega_{c}-\omega_{z}|$.

Further considering the limit $\omega_{c}\gg\omega_{z}$ one finally has

[TABLE]

For proper initial states of the photon in the cavity, it is possible to map Eq. (57) into Eq. (7) (up to a constant) by defining

[TABLE]

with $n$ the averaged number of photons into the cavity. Notice that in this case the renormalization of the effective coupling due to the presence of the cavity photons can be exploited to increase the value of $g$ with respect to the previous case.

Appendix B Perturbative approach to the weak coupling case

It is possible to study the dynamics induced by Eq. (7) and the figures of merits discussed in the main text by applying the time-dependent perturbation theory in terms of the Dyson series.

Working in the Interaction picture (indicated in the following by the apex $I$), the time evolution operator $\hat{U}_{I}(t,0)$ from the initial time $t=0$ to a generic time $t$ satisfies the differential equation

[TABLE]

with

[TABLE]

Imposing the initial condition $\hat{U}_{I}(0,0)=1$ and formally integrating the above equation we have

[TABLE]

We can solve iteratively this equation in terms of Dyson’s series Sakurai_Book

[TABLE]

In this picture, the energy stored in the QB can be written as

[TABLE]

After some long and tedious algebra, one finally obtains

[TABLE]

and consequently

[TABLE]

Appendix C Comments on the classical limit

In the main text, we have considered the scaling at a finite number of TLSs, without addressing the thermodynamical limit needed for a proper study of the ESQPT. Some hints in this direction can be derived starting from the classical version of the Hamiltonian in Eq. (7). This can be obtained replacing the operators $\hat{S}_{x,y,z}$ with the classical angular momentum components

[TABLE]

This leads to the classical Hamiltonian

[TABLE]

One can identify now two conjugated variables Andolina19

[TABLE]

[TABLE]

in such a way that

[TABLE]

Recalling the equations of motion

[TABLE]

Further manipulating the above Equations introducing rescaled and shifted variables

[TABLE]

[TABLE]

one has

[TABLE]

where clearly emerge the dependence on the renormalized interaction parameter

[TABLE]

which plays a major role in the discussion about the universality reported in the main text.

\reftitle

References