Characterization of a Two-Photon Quantum Battery: Initial Conditions, Stability and Work Extraction

We consider a quantum battery that is based on a two-level system coupled with a cavity radiation by means of a two-photon interaction. Various figures of merit, such as stored energy, average charging power, energy fluctuations, and extractable work are investigated, considering, as possible initial conditions for the cavity, a Fock state, a coherent state, and a squeezed state. We show that the first state leads to better performances for the battery. However, a coherent state with the same average number of photons, even if it is affected by stronger fluctuations in the stored energy, results in quite interesting performance, in particular since it allows for almost completely extracting the stored energy as usable work at short enough times.


Introduction
Quantum thermodynamics is a growing field of research that aims at extending concepts, such as heat and work to the realm of quantum physics [1][2][3][4][5][6][7][8][9]. In this framework, several conventional assumptions of classical thermodynamics should be reconsidered. In particular, one needs to depart from the concept of thermodynamic equilibrium to properly include quantum coherences or possible external drives acting on individual quantum systems. In recent years, a very active branch of this new field of research focused on the study of energy transfer and storage in quantum devices leading to the concept of Quantum Battery (QB) [7,10,11].
Various strategies have been developed in order to exploit quantum features to outperform classical devices, in particular enhancing the charging power (the energy that is stored in a given time interval) and the extractable work [10,[12][13][14][15][16][17][18][19][20][21][22]. Theoretical proposals for possible actual implementations of QBs are often based on engineered two-level systems (TLSs) such as trapped ions [23][24][25], superconducting qubits [26] or quantum dot in semiconductors [27]. These platforms, which are usually exploited for qubit implementations, can be used as QB and can be coupled to chargers of a different nature able to coherently transfer energy into them. The proposed charging mechanisms range from classical external drives [28][29][30][31] to a proper control of the interaction in an array of TLSs [32][33][34][35][36].
In this respect, promising platforms are devices where a collection of TLSs is coupled to a monochromatic cavity radiation in the same spirit of what has been realized in cavity and circuit quantum electrodynamics [37,38]. Indeed, it has been shown that, for such kind of QBs, an enhancement of the averaged charging power, which scales as √ M, with M being the total number of TLSs, can be achieved, due to the collective interaction among the TLSs [39][40][41][42]. Interestingly, a similar scaling, compatible with theoretical predictions, has been reported in a very recent experiment [43], demonstrating, for the first time, the collective quantum advantage first discussed in Ref. [13].
In addition, recent theoretical proposals have shown that it is possible to engineer, in both trapped ions [44] and flux qubit systems [45], the suppression of the conventional dipole coupling, linear in the quantum field radiation, in order to access the peculiar phenomenologies, such as super-radiance and spectral collapse, which are associated to the two-photon coupling (quadratic in the radiation field) [46][47][48][49][50][51]. These ideas have also recently been considered in the context of QBs by some of us, showing that two-photon coupling, once realized, can lead to an even greater (proportional to M instead of √ M) enhancement of the averaged charging power [52]. This analysis has been carried out considering as initial state for the cavity radiation a Fock state where the number of photons is exactly the double of the number of TLSs in such a way that, at resonance, all the radiation energy is transferred to the QB (charging) and back (discharging). However, this form for the initial state is quite ad hoc and not so easy to be implemented experimentally. According to this, it is relevant to investigate also other possible initial states for the cavity radiation, such as coherent states and squeezed states, which are easier to be realized from an experimental point of view [37,53].
In the present paper, we address this aspect to achieve a full characterization of the two-photon based charging mechanism. For the sake of clarity, we will consider the case of a single TLS as the QB to be charged, avoiding other competing mechanisms that are based on collective behaviour [28,33]. We will also adopt the so-called rotating-wave approximation (RWA), in the same spirit of Ref. [54] for the single-photon case, which allows for a simple analytical treatment. We will describe the charging of a QB investigating different figures of merit. In particular, we will study the energy storage and associated charging time. Moreover, we will consider the average charging power and the fluctuations of the stored energy [29,52,55,56]. Finally, a discussion regarding the maximum extractable work, namely the maximum fraction of the energy stored into the battery, which can be actually extracted and used to further purposes, will be addressed. This quantity, which is known in the literature as ergotropy [57], in general does not coincide with the total energy that is stored in the QB due to quantum correlations [15,19,40]. This analysis has already been carried out for a conventional single-photon coupling [40] identifying the Fock and coherent state as the best choices for extracting almost all of the stored energy, but is still lacking for a two-photon coupling.
All of these different aspects will be analyzed while taking various initial conditions for the cavity states into account. In particular, we will consider a Fock state, a coherent state, and a squeezed state with the same average number of photons as representative of the more conventional quantum state for the cavity radiation. Our study confirms that a properly designed Fock state is optimal for the functioning of the QB leading to maximal stored energy, short charging time, and a good amount of extractable work. However, interesting performances can also be found for a coherent state with the same averaged number of photons. In particular, this state shows an important ratio between extractable work and stored energy at short enough times. An opposite behaviour is observed for a squeezed state whose performances are all around very poor. The present analysis could give relevant hints for actual experimental implementations of these kind of QBs.
The paper is organized, as follows. In Section 2, we introduce the model of QB, where a single TLS is coupled to a unique cavity mode via a two-photon interaction. Here, we consider different initial states for the cavity such as a Fock state, a coherent state and a squeezed state. Section 3 is devoted to the characterization of the QB for these three different cavity states. In particular, we study the charging of the QB when considering the stored energy and corresponding charging times, as well as the related charging power. Moreover, we investigate the energy fluctuations and maximum work that can be extracted from the QB. All of these figures of merit are relevant in view of future development of this kind of QBs. In Section 4, we draw the conclusions of our work.

Model
We consider a QB described as a single TLS coupled to the radiation of a cavity through a two-photon coupling (quadratic in the electron field) with matter-radiation coupling λ. This simple single cell can be replicated M times to form a QB that works in a parallel charging configuration [39,54]. The theoretical possibility of engineering such a kind of coupling in realistic devices has been recently discussed for trapped ions subject to bichromatic driving [44] and for flux qubit coupled with a symmetric dc superconducting quantum interference device (SQUID) [45]. The corresponding Hamiltonian can be written aŝ whereâ (â † ) are the annihilation (creation) operators of the cavity radiation,σ x ,σ z are the Pauli matrices, ω a the level spacing of the TLS, and ω c the frequency of the radiation in the cavity. We underline that the TLS and cavity do not interact at t < 0 and we assume that the two-photon interaction is switched on at time t ≥ 0, as indicated by the θ(t) step function. According to this, it is possible to transfer energy from the cavity to the TLS and observe the charging of the QB.
In the limit where the matter-radiation coupling λ satisfies λ ω a , ω c , one can consider the RWA of Equation (1) and describe the system in terms of the two-photon Jaynes-Cummings model [52,58,59] with Hamiltonian being given bŷ where the first term is associated with QB, the second with the cavity micro-wave radiation used as charger, and the last one is the two-photon interaction term in the RWA [44,45] and withσ ± =σ whereσ y is one of the Pauli matrices.
In the following, we will focus on the resonant case ω a = 2ω c (see Figure 1), where all of the energies associated to the radiation can be transferred to the TLS and viceversa [60]. This condition guarantees an optimal functioning of the QB, which is not achieved in other configurations of the device [52]. The initial condition at time t = 0 for the total system is chosen as a tensor product state where the first factor represents the ground state of the TLS, while the second characterizes the charger with |n , a state with n photons in the cavity. We will consider, as paradigmatic examples, a Fock (fixed number of photons), a coherent, and a squeezed state for the cavity. Thus, we introduce the probability amplitudes α n , such that where the superscripts F, C, S denote the Fock, coherent, and squeezed state, respectively. We underline that N represents exactly the number of photons for the Fock state, while it is the average number N of photons for the other two states where fluctuations in the photon numbers are present. Moreover, the probability amplitudes α (S) n are valid only for even values of n and are zero otherwise.
It is worth mentioning that the QB+charger is a closed system. Interactions with possible external environments can lead to energy relaxation of the TLS and the loss of photons in the cavity. These processes can be characterized in terms of typical time scales t r and t γ respectively [37,[61][62][63]. According to this we need to restrict our analysis to times, such that t t r , t γ , in such a way to safely neglect the dissipation effects. However, we note that these constraints are comparable with the ones that are required by quantum information protocols [26]. Moreover, a protocol that is able to overcome the detrimental effects of energy relaxation in TLS have been theoretically discussed very recently in Ref. [16].

Figures of Merit of the Quantum Battery
We now characterize the performances of a two-photon QB comparing the three cavity states that were introduced in the previous Section.
In order to investigate the performance of the QB, it is necessary to study the time evolution of the initial state in Equation (4). We observe that, over the basis given by the states |g ⊗ |n and |e ⊗ |n − 2 (n ≥ 2), with |e excited state of the TLS, the Hamiltonian in Equation (2) where the notation keeps track of the number of photons n and the dynamics remains confined in a two-dimensional space. This 2 × 2 matrix can be diagonalized in terms of the eigenstates with eigenvalues In terms of these states, one can rewrite Equation (4) as that evolves in time as It is instructive to also consider the time dependent behaviour of the reduced density matrix that was obtained by tracing out the photon degrees of freedom when we consider the initial state in Equation (11) ρ TLS (t) = ∑ n n|ρ(t)|n = ∑ n p n sin 2 (λ n(n − 1)t)|e e| + p n cos 2 (λ n(n − 1)t)|g g| + e iω a t |g e| i √ p n+2 p n sin(λ (n + 1)(n + 2)t) cos(λ n(n − 1)t) − e −iω a t |e g| i √ p n+2 p n sin(λ (n + 1)(n + 2)t) cos(λ n(n − 1)t) where p n = |α n | 2 andρ is the total density matrix at time t with |ψ(t) the state in Equation (12). Equation (13) can be written in the conventional form using the Bloch vector. Being written in terms of the identity matrix and the Pauli matrices, it reads [64,65] ρ TLS (t) = Here, starting from Equation (13) and recalling that we can obtain the expressions for the components u(t), v(t), and w(t) u(t) = − sin (ω a t) ∑ n √ p n p n+2 sin(λ (n + 1)(n + 2)t) cos (λ n(n − 1)t), (20) v(t) = cos (ω a t) ∑ n √ p n p n+2 sin(λ (n + 1)(n + 2)t) cos (λ n(n − 1)t), Notice that, for the simple case of a Fock state, where the probability amplitude α n is given by Equation (5), the above expression assumes the simple form Instead, the explicit form for u(t), v(t), and w(t) for the coherent and squeezed states can be obtained by replacing the α n given in Equations (6) and (7) into Equations (20)- (22), and cannot be further simplified.

Stored Energy and Average Charging Power
The energy stored in the QB at time t is given by [7,29,39,52,66] the contribution to the total Hamiltonian in Equation (2) that is associated to the QB. Inserting the time evolved state of the system [see Equation (12)] into Equation (24) one has An important task in the context of QB is to store the maximum amount of energy in the fastest time. In this perspective, we define the maximum of the stored energy as [52,54] where t E is the time at which the maximum occurs. In the case of a Fock state (α (F) n = δ n,N ), it is possible to analytically find the value of t (F) E imposing the condition sin 2 (λ n(n − 1)t) = 1.
Consequently the value of the maxima of the energy occurs when where k ∈ Z indicates which maximum is considered. Instead, in the following, for the coherent and squeezed state, the charging time t E is obtained numerically.
In Figure 2, we report the behavior of the energy E(t) (see Equation (26)  Here, one can see that a complete charging is only reached if the cavity is in a Fock state (E (F) max = ω a ). However, a coherent state with the same average number of photons (N = 8 in the considered case) can also reach a quite large fraction of the maximal stored energy (the first maximum is E (C) max ≈ 0.76ω a ). This aspect can be particularly relevant for the actual experimental realization of this kind of QB due to the fact that a coherent state is usually easier to realize in a cavity with respect to a well defined Fock state [53]. Conversely, the squeezed state only reaches a very poor maximum charging E (S) max ≈ 0.52ω a . Moreover, the minimal charging times for the Fock and coherent state are almost comparable and both are faster with respect to the squeezed state. Indeed, from Equation (29), we have that λt 19. Instead, the maximum of the energy of the squeezed state is achieved for the longer time λt (S) E ≈ 1.53. It is also worth pointing out that, for a given quantum state of the charger, the charging times in the two-photon interaction model are usually shorter than the corresponding ones for a conventional dipolar singlephoton coupling with the same average value of photons [39,52,54]. In fact, as reported in Ref. [54], the energy for a single photon process is given by E 1ph (t) = ∑ n p n sin 2 (λ √ nt), meaning that the time at which the maximum of the energy is reached scales as 1/ √ n. Instead, in the two-photon case, it scales as 1/ n(n − 1), leading to a charging time, which is √ n − 1 times faster fixing all other parameters. From the previous discussion, one can infer that the Fock state is the best choice for the initial state of the charger, in order to store energy into a QB, in agreement with Ref. [54]. Conversely, the squeezed one is the less efficient, due to both the longer charging time and smaller maximum stored energy when compared to the other ones. Notice that the Fock states still show the best charging performances in terms of stored energy and charging times also in comparison with other states of the general form in Equation (4).
The charging behaviour can be better understood by looking at the evolution of the TLS state on the Bloch sphere. In Figure 3, we show the path that is followed by the TLS state on the Bloch sphere for the three different cavity states, up to the first maximum of the stored energy. Notice that, here, the |g and the |e states are respectively represented by the south and north pole. As discussed above, we initialize the TLS in the ground state. Here, we can observe that, only for the Fock state, the TLS state reaches the excited state, corresponding to the complete charging of the QB, while, for the coherent and the squeezed states, this never happens and the path of the state vector on the Bloch sphere is much more complicated. It is worth pointing out the fact that more involved initial states for the TLS, such as coherent superpositions of ground and excited state, are also characterized by complicated evolution in the Bloch sphere.
Another relevant figure of merit is the average charging power, which is defined as [7,13,39,66] Regarding the energy, also in this case we are interested in achieving the maximum value of the charging power in the fastest possible time. Therefore, one needs to consider the maximum charging power [39,52,54] where t P represents the time at which the maximum occurs.
In Figure 4, we report the behaviour of P(t) in Equation (30) as a function of time for the three considered initial states of the cavity. Additionally, in this case, one finds that this quantity is maximal for the Fock state (P

Energy Fluctuations
In order to have a complete characterization of the QB, we now evaluate the quantum fluctuations that are associated to the stored energy, since they can induce detrimental effects on the charging performances [29,52,55]. We discuss the stability of the charging process in terms of the fluctuations of the stored energy at equal times, as represented by the correlator [29,52,55] with |ψ(0) the initial state in Equation (11) andĤ QB (t) the TLS Hamiltonian evolved in time in the Heisenberg representation according toĤ.
In Figure 5, we show the evolution of Ξ(t) for the three initial states of the cavity. In the following, we are interested in studying the value of the correlator at time t E , where the stored energy has its maximum, to understand how it affects the functionality of the QB. We then define Here, it clearly emerges that the Fock state has no energy fluctuations in correspondence of the maximum of the charging energy (Ξ (F) = 0), as also reported in Ref. [52]. This can be understood from Equation (32) and it is a consequence of the fact that, at time t E , the quantum state of the system (QB+charger) is separable (|ψ(t E ) = |e ⊗ |n − 2 ). This is not true for both the coherent and the squeezed states that show relevant fluctuations in the stored energy. Because neither of these states achieve E max = ω a , we will never observeΞ = 0. In fact, for the coherent state, where the maximum energy was higher, we obtainΞ (C) ≈ 0.43ω a , while, for the squeezed state, which is even worse in terms of energy storage, we haveΞ (S) ≈ 0.49ω a . Again, this confirms that the Fock state appears as the most convenient cavity state to build a good QB. Indeed, this particular initial state is not subject to stored energy fluctuations.

Ergotropy
Another relevant quantity to look at to determine the efficiency of a QB is the so-called ergotropy [10,57]. It consists in the maximal stored energy that can be converted into usable work. It can be extracted from the QB at a given time t of its evolution (charging). In general, this quantity is different from the stored energy due to the fact that part of the energy may be locked into correlations and, therefore, cannot be extracted for further purposes [67]. A general derivation of this quantity starts from the Hamiltonian of the QB, written as [10]Ĥ QB = ∑ n n | n n | (34) with the energy eigenvalues ordered, such that n < n+1 and with | n associated eigenvectors, as well as the density matrix at a given timê ρ TLS (t) = ∑ n r n (t)|r n (t) r n (t)| (35) ordered, such that r n > r n+1 and with |r n eigenvectors. According to this, the work done on the system after a time t is given by The maximum work that can be extracted from the battery after a time t, called ergotropy, is defined as E = max(−W). Exploiting Equations (34) and (35), the ergotropy can be written as [57] It is worth pointing out that E = 0 when the initial state of the system is passive [10]. Moreover, in general, the maximum extractable work E th can be obtained when the final state of the system is thermal, i. e. ρ TLS (t) = e −βĤ QB /Z (with β the inverse temperature of the system and Z = Tr{e −βĤ QB }). Subsequently, the bound on the maximum and minimum extractable work is given by 0 ≤ E ≤ E th [57]. We now want to find the explicit form of the ergotropy for our system. To do so, we need to diagonalize the density matrix of the TLS in Equation (13). Its eigenvalues are where s = 0, 1. Moreover, the eigenvalues ofĤ QB in Equation (25) are ε s = (−1) s ω a /2. Consequently, starting from Equation (37), the ergotropy at time t can be written as Figure 6 shows its time-evolution. Here, one has a qualitatively similar behaviour with respect to the one that is observed for the stored energy. We observe that only the Fock state reaches the maximum of the ergotropy E max = ω a . This occurs at a time t where the energy has its maxima and the state of the TLS is pure, confirming the relevant performances of this cavity state.
To better quantify the actual fraction of extractable energy, we now evaluate the ratio between the ergotropy and stored energy, as given by In Figure 7, we compare the behavior of the ratio η for the three considered initial states as a function of time. As one can observe, at quite short times the coherent state allows an almost complete extraction of the energy as usable work. However, in this region, the energy that is stored in the system can be quite small [see Figure 2]. For the Fock state it is possible to completely extract the stored energy only for narrow time windows in correspondence of t (F) E , with the advantage that the stored energy approaches or is equal to the maximum value E This result is a consequence of the nature of the three considered states. In fact, as also stated in Ref. [40], the more there is a mixing between the TLS and the photonic part, the less energy we can extract, since the entanglement between the QB and charger has a negative impact on the possibility to extract work.

Conclusions
We have characterized various figures of merit for a quantum battery given by a qubit, namely a two-level system, coupled with a cavity radiation through a two-photon coupling (quadratic in the quantum electromagnetic field). As possible initial conditions for the quantum radiation in the cavity, we have inspected a Fock state, a coherent state, and squeezed state. According to our analysis, the Fock state emerges as the ideal candidate for optimizing the performances of the quantum battery, being the only state able to reach a complete charging of the quantum battery, in short times and without showing fluctuations. Quite remarkably, this kind of state also allows for completely extracting its maximal stored energy. A coherent state with the same average number of photons also shows quite interesting performances, in particular for what concerns the fraction of extractable work at short enough times, even if it is affected by fluctuations in the stored energy. This is not true for a squeezed state whose performances are very poor, which makes it the worst state for implementing a quantum battery. Indeed, it can only store a fraction of the optimal energy and it is affected by very strong fluctuations. The present analysis gives important hints towards the possible implementations of quantum batteries coupled to cavity radiation.

Data Availability Statement:
The data and analysis used in this work are available from the corresponding author upon reasonable request.