Quick Charging of A Quantum Battery With Superposed Trajecotries
Quick Charging of A Quantum Battery With Superposed Trajecotries
Quick Charging of A Quantum Battery With Superposed Trajecotries
We propose novel charging protocols for quantum batteries based on quantum superpositions of
trajectories. Specifically, we consider that a qubit (the battery) interacts with multiple cavities or a
single cavity at various positions, where the cavities act as chargers. Further, we introduce a quan-
tum control prepared in a quantum superposition state, allowing the battery to be simultaneously
charged by multiple cavities or a single cavity with different entry positions. To assess the battery’s
performance, we evaluate the maximum extractable work, referred to as ergotropy. Our main result
is that the proposed protocols can utilize quantum interference effects to speed up the charging
process. For the protocol involving multiple cavities, we observe a substantial increase in ergotropy
as the number of superposed trajectories increases. In the case of the single-cavity protocol, we
arXiv:2307.09010v1 [quant-ph] 18 Jul 2023
show that two superposed trajectories (entry positions) are sufficient to achieve the upper limit of
the ergotropy throughout the entire charging process. Furthermore, we propose circuit models for
these charging protocols and conduct proof-of-principle demonstrations on IBMQ and IonQ quan-
tum processors. The results validate our theoretical predictions, demonstrating a clear enhancement
in ergotropy.
aˆ † .
Here, aˆ j aˆ † annihilates (creates) a photon in Cj with
j 1 Σ
frequency ωc, ¯h ω a represents the energy splitting between |m⟩ ⟨n|D . (9)
the ground state |g⟩ and the excited state |e⟩ of Q. The |ξk=1⟩ ⟨ξk=1|D ≡
N
m ,
n=1
| ⟩e⟨ |−|
Pauli operators are therefore given by σˆ z = e g⟩ ⟨ g
| According to the assumption that all chargers are iden-
and σˆ x = |e⟩ ⟨g |+ g| ⟩e⟨ . Moreover, the dimensionless
tical as well as the orthonormality of the projectors, in
constant λ |represents the coupling strength between Q
and all the chargers. Throughout this work, we focus on Appendix A, we show that the explicit form of the rest
the resonant regime ωa = ωc. of the projectors is irrelevant, enabling us to further sim-
We first send the battery Q into the multiport beam plify the analysis.
splitter (MPBS1 in Fig. 1). In general, the beam splitter Here, we evaluate the performance of a QB by con-
can prepare the trajectories in a quantum superposition sidering the ergotropy, which quantifies the maximum
state so that Q can be charged by these N chargers extractable work. Given a charged state of the QB ρ(τ),
simul- taneously. For simplicity, we assume the the ergotropy is defined as
superposition state of the trajectories is W (ρ(τ)) ≡ Tr (ρ(τ)HQ) − min Tr Uρ(τ)U † H Q
N U (10)
Σ = Tr (ρ(τ)HQ) − Tr (φ(τ)HQ) ,
|ψ⟩D = √1 |j⟩D . (3)
N i=1 where U represents the unitary operation for work ex-
We consider that the battery and chargers are initialized traction. Also, φ is known as the passive state [21] (as-
in the ground state and the single-photon Fock state, sociated with ρ), which cannot provide useful work for
respectively. Therefore, the total initial state all possible work extraction operations U. According to
O reads Ref. [21], the passive state of the battery can be written
1 Σ as
N N
|ψ (0)⟩DQC = √ |j⟩D ⊗ |g⟩Q |1⟩Cj . (4)
N j=1 φ = s0 |e⟩ ⟨e| + s1 |g⟩ ⟨g| . (11)
j=1
where |ϕj (τ)⟩QC is defined as E (τ) ≡ Σ pk (τ) Tr [HQρk (τ)] − Tr (HQ |g⟩
⟨g|) (12)
k=1
Finally, we make these trajectories interfere with pro- jectors {Pk}k acting on D, namely
one another by using another beam splitter (MPBS2
in Fig. 1). In principle, MPBS2 has N different outputs, Pk = |ξk⟩ ⟨ξk|D ,
which can be described using a set of orthonormal
Σ
4
σk(τ). We can obtain W ≤ E, be-
cause Tr(HQ g g ) Tr(HQφ) in general. In addition, | ⟩ ⟨outputs
are N different | ≤ (labeled as k ). In principle,
the inequality saturates if and only if φ = g g , which | ⟩ ⟨ for each
one can find the optimal work extraction strategies
implies that ρ is a pure state. output and obtain the average ergotropy, i.e.,|
As aforementioned, in our charging protocol, there Σ
N W= pkW (ρk(τ)) { }
Pk = 11, (7) Σk
(13)
k=1
′ = pk [Tr (ρk(τ)HQ) − Tr (φk(τ)HQ)] ,
⟨ξk|ξk′ ⟩ = δk,k′ ∀ k, k . k
5
where φk denotes the passive state associated with ρk. Here, for the case k = 1, the post-measurement state
Following similar reasoning as mentioned earlier, the av-
is passive during the time√period τ ∈ [0, TN ] with the inver-
erage ergotropy is also upper bonded by the stored en-
ergy, and the optimal extractable work can be obtained, sion time TN = tan−1( N)T/2π. Remarkably, for the
i.e., W (τ) = E(τ), if and only if pk, ρk forms a pure cases of k ̸ = 1, the post-measurement states are
{ } exactly the excited state, implying that the maximal
state decomposition of ρ.
extractable work h¯ ω c can be obtained. Therefore,
To gain some analytical insight, we now consider the the average er- gotropy can be expressed as
rotating-wave approximation, which is usually valid for
coupling strength λ ≤0.1 [52], so that the interaction W (τ)
Hamiltonian of the QB and the chargers in Eq. (2) can = W (σk=1(τ)) + (N − 1)W (σk̸=1(τ))
T
¯hω 2 sin2 (ωcλτ) − 1 if TN ≤ τ ≤ .
c
( 6
be reduced to the Jaynes-Cummings model, namely ¯h ω N−1 2
sin (ω λτ) if 0 ≤ τ ≤ T (20)
′ †
c N c N
H˜ =
j
σˆ + aˆ j + , (14) =
h¯ωcλ
σˆ− aˆj 4
QC
where σˆ + =| e⟩ ⟨g|and σˆ − =| g⟩ ⟨e represent the creation In Fig. 2, we present the time-dependent stored energy
|
and annihilation operators of Q, respectively. We can and the average ergotropy for different values of N. In
then evaluate Eq. (6) in this case: contrast to the case of N = 1, we observe non-zero av-
O erage ergotropy for the entire interval of interest
N
|ϕj (τ)⟩QC = − i sin (ωcλτ) |e⟩Q ⊗ aˆ j |1⟩C ′ because the states with k ̸= 1 are non-passive right after
the QB-
j′=1
j
chargers interaction is turned on (τ > 0). Therefore, the
(15) protocol can be used for “quick charging”, enabling im-
N
O mediate storage of useful work after the charging
+ cos (ωcλτ) |g⟩Q |1⟩Cj . process begins. Furthermore, the result indicates that
j=1
increasing N delays the inversion time TN and
enhances the av-
Let us start from the simplest case with only one charger erage ergotropy before TN . According to Eq. (20), this
(i.e., N = 1), where the reduced state of Q is expressed enhancement originates from the increase of the purity
as for the output k = 1 as its excited state population de-
ρ(τ) = sin2 (ωcλτ) |e⟩ ⟨e| + cos2 (ωcλτ) |g⟩ ⟨g| . (16) creases with N. This result aligns with the Zeno-like
state freezing effect described in Ref. [45]. In the asymp-
In this case, the stored energy is totic limit (N → ), we can further obtain a pure state
E(τ) = h¯ ω c sin2 (ωcλτ) , (17) decomposition,∞i.e., σk=1(τ) ∝ g | ⟩ ⟨g | σk̸=1(τ) ∝ e| ⟩ e⟨ ,
implying that the stored energy| can be fully converted
which oscillates with a period T = 2π/(ωcλ). We now into extractable work for the whole time interval, i.e,
focus on the time interval τ [0, T/4] ∈ (such that ωcλτ [0, E(τ) = W (τ).
π/2]), ∈where the stored energy monotonically increases
from 0 to its maximum value ¯hω c . Note that ρ is
diagonalized under basis e{|, g ⟩ |. ⟩Thus,
} according to Eq. 1 E(τ )
(10), the criterion for obtaining non-zero ergotropy is the W [N = 4]
moment that population inversion occurs, where the W [N = 2]
excited state population becomes larger than the ground W [N = 1]
state population (⟨ e| ρ(τ) | e⟩ > ⟨g| ρ(τ) |g⟩ ). The time-
dependence ( of the ergotropy can then be derived as
[h¯ωc ]
E, W
0 2 if 0 ≤ τ < T
8T
,
W (τ) = ¯hωc (2 T (18)
(ωcλτ) − 1) if 8 ≤τ ≤ 4 .
sin
One can observe that in the duration τ [0, T/8], al- 0
∈
though the stored energy E monotonically increases,
0 π/4 π/2
there is no extractable work, W = 0, for the battery
because ρ remains a passive state during this period. τ [ωc−1λ−1]
We now consider the scenario involving N chargers.
When the selective measurements satisfy Eq. (7) and FIG. 2. The stored energy E and average ergotropy W
Eq. (9), the unnormalized post-measurement states (both in units of ¯hω c ) on time τ (in units of 1/ωc) for λ =
can be written as (see Appendix A for detailed 0.05. The black solid curve plots the stored energy while the
derivations): dashed curves plot the average ergotropy. From bottom to
1 top, the blue, red and green dashed curves show the results
σk=1(τ) = c |e⟩ ⟨e| + cos (ω cλτ) |g⟩ ⟨g| ,
sin2 (ω λτ) 2
for N = 1,
N
1 N = 2 and N = 4, respectively.
σ (τ) = sin2 (ω λτ) |e⟩ ⟨e| .
k̸=1 c We now extend our scope of discussion into the
N (19) ultrastrong–coupling regime [52], where the rotating–
7
E, W k /= 1
[h¯ωc]
E, W
[h¯ωc]
0.0 0.0
0 0.63 1.00 π/2 0 π/4 π/2
—1 —1
τ [ωc λ ] τ [ωc—1λ—1]
0.8 (b) (d)
1.0
P
E, Wk=1
[h¯ωc]
0.0 0.5
0 0.63 1.00 π/2 0 π/4 π/2
—1 —1
τ [ωc λ ] τ [ωc—1λ—1]
8
FIG. 3. (a) The stored energy E and average ergotropy W (both in units of ¯hω c ) as functions of time τ (in units of 1/ωc)
for λ = 0.5. (b) The stored energy E and the average ergotropy contributed by states of k = 1 Wk=1. Two brown vertical
lines at τ = 0.63, 1.00 indicate the inversion points TN for N = 1(blue) and N = 2(red). Here, Wk=1 is 0 when N =
7(green) and N = ∞(magenta). (c) The stored energy E and the average ergotropy contributed by states of k ̸= 1 Wk̸=1 on
time τ . Except for N = 1 (blue), these states contribute average ergotropy when τ > 0. (d) The dashed curves plot the
change in average purity P against time τ . The cutoff photon number is set to 9 in the above results.
wave approximation is no longer valid. In Fig. 3 (a), we contribution of k̸ = 1, since the corresponding average
present the dynamics of the average ergotropy and the er- gotropy becomes non-zero as soon as the charging
stored energy. We can still observe the quick charging process begins when N > 2. However, in contrast to the
effect, a delay of the inversion point, and an enhance- previous results with the rotating–wave approximation,
ment in average ergotropy as N increases. As shown the aver- age ergotropy cannot reach the upper bound
in Figs. 3 (b) and (c), we further present the individual even in the asymptotic limit, implying that the post-
contributions of the average ergotropy from the outputs measurement states are not pure.
k = 1 and k ̸= 1, which are respectively defined by In Fig. 3(d), we present the average purityP associated
W k=1 = pk=1W (ρk=1(τ)) with the post-measurement states, which is defined by
and " #
Σ σk 2
Σ (21) Tr[σk]Tr . (22)
W k̸=1 P= Tr[σ ]
= p W (ρ
k
(τ)).
k k k
k̸=1
We can observe that the inversion points come from the We can observe the overall average purity increasing as
contribution of k = 1. Further, we can observe that its N increases, thereby leading to the enhancement of
contribution decreases as N increases. When N ≥ 7, aver- age ergotropy. However, the average purity cannot
the contribution vanishes for the entire charging period. reach unity even in the asymptotic limit. Thus, the stored
In Appendix A, we provide analytical analysis, showing en- ergy cannot be fully converted into extractable work.
that the decrease in average ergotropy can also attribute In Figs. 4 (a) and (b), we compare the maximum average
to the Zeno-like state freezing effect. More specifically, ergotropy with the average purity with respect to the
we demonstrate that the excited state population cou- pling strength λ, where the maximum average
decreases when N increases. In the asymptotic limie ergotropy is defined by
(N → ∞),
the state can be frozen in the ground state. We further
W max = max W (τ). (23)
observe that the quick charging effect originates from the τ
∈[0,T/4]
9
1.0 (a)
where HC = h¯ ω c aˆ † aˆ . Here, H′ represents the
interaction
Hamiltonian when QB is located at the position rj. Its
explicit form is given by
[N = 1]
j c rj x
[N = 10] H′ = h¯ ω λ cos πL σˆ aˆ + aˆ † , (25)
max
W
1 Σ
N
[N = 1] |j⟩D ⊗ |g⟩Q ⊗ |1⟩C , (27)
P
|ψ(0)⟩DQC = √
[N = 10] N j=1
[N = 100]
[N = 1]
10
and allow it to evolve according to the total Hamiltonian
Htot, namely
0.6
1 Σ
N
0.05 0.25 0.50 |ψ(τ)⟩DQC = √ |j⟩D ⊗ |ϕj(τ)⟩QC , (28)
∏ N j=1
c
j=1
Therefore, one can observe that the total evolution (in- dt1 · · · dtnT HI (t1) · · · HI (tn) |g⟩Q ⊗ |1⟩C .
cluding the qudit, the QB, and the charger) is described 0 0
(40)
by a linear combination of the position-dependent prop-
In this case, all the even terms in the Dyson series vanish,
agators, which characterizes the collective quantum in-
implying the QB is in the excited state |e⟩ for all τ > 0
terference effect among different positions [45].
(with zero probability of obtaining the output k = 2 at
We now show that two superposed trajectories (posi-
tions) can lead to the saturation of the ergotropy to its τ = 0). Because the post-measurement QB states for
upper bound with an appropriate adjustment of the col- these two outputs are pure states (i.e, the average pu-
lective interference effect. More specifically, we consider rity is one), we can conclude that the stored energy can
that the two positions satisfy r1 + r2 = L, such that be fully converted to the extractable work, i.e., E = W ,
throughout the whole charging process. Note that the
presented analysis does not rely on the rotating wave
r2 r1 r1 approximation, thus indicating the saturation of the er-
cos π L = cos π − π L = − cos π L . (36)
gotropy to its upper limit generally holds for all regimes
Therefore, the coupling strengths share the same magni- of the QB-charger coupling strength.
tude but are completely out of phase. Furthermore, we
consider
1
|ξk=1⟩ = √ (|1⟩D + | IV. IMPLEMENTATION ON
2⟩D), 2 QUANTUM DEVICES
(37)
1
|ξk=2⟩ = √ (|1⟩D − | In this section, we provide circuit models for the pro-
2⟩D), 2
and, according to Eq. (31), the coefficients in this case posed charging protocols and perform proof-of-concept
are experiments on the quantum processors provided by
IBMQ and IonQ, which involves two superposed trajec-
1 tories (N = 2).
c1,1 = = = 2,2 = . (38)
−c 2 The quantum circuit for the multiple-charging setup is
c1,2 c2,1
Therefore, the time evolution of the post-measurement described by Fig. 5 (a). The circuit consists of 4 qubits,
state of the QB and the charger for k = 1 is given by representing the control qubit D, the quantum battery
2 Q, and the two charging cavities C1 and C2, respec-
Σ
c tively. The circuit can be divided into three parts: state
1,j Uj,I (τ, 0) |g⟩Q ⊗ |1⟩C preparation, charging process, and measurements on the
j=1
control qubit and quantum battery. In the state prepa-
1 ration part, the qubits are prepared in the initial state
= [U1,I (τ, 0) + U2,I (τ, 0)] |g⟩Q ⊗ | specified in Eq. (4) using single-qubit gates. The charg-
1⟩2C
∞
Σ 1 ing process involves the utilization of two controlled-
cos r1 n
= − π unitary gates to simulate the simultaneous charging of
i
12
n! ¯h L
n even the qubit by the two cavities through Jaynes-Cummings
∫ τ ∫ τ
interactions. In Fig. 5 (b), we present the decomposi-
× dt1 · · · dtnTˆH′I(t1) · · · H′ (t
I n) |g⟩ Q⊗ |1⟩ C .
0 0
tion of the controlled–unitaries into bit-flip (X) gates,
(39) controlled-z gates (CZ), and Ising coupling gates [XX(θ)
13
D Q C1
C1 : |1ih1 T
exp[—iHQC1 ⌧/~] r Z
Q : |0ih0 C2
exp[—iHQC2 ⌧/~] (d)
C2 : |1ih1 T XX [✓(⌧)]
(c) r
H H
exp[—iHQC1 ⌧/~] XX [—✓(⌧)] YY [—✓(⌧)] XX [✓(⌧)] YY [✓(⌧)] H Rz [✓(⌧)] H
FIG. 5. (a) Quantum circuit for multiple-chargers protocol. Here, D, Q, C1, C2 represents the control qubit, battery qubit, first
charger and second charger, respectively. (b) Decomposition of a controlled unitary in (a). (c) Qubit configuration used on
ibmq algiers. (d) Decomposition of a XX gate into CNOT gates.
Here, we map the charging time τ into the angle θ FIG. 6. Charging process of the single-charger setup, where
using the following relation the coupling strengths of the two superposed trajectories have
the same magnitude but are out of phase.
θ(τ) = ωcλτ/2. (42)
Finally, in the measurement part, we measure the con- average of 1000 experiment repetitions. The experimen-
trol qubit D in the x-direction, aligned with the projec- tal results demonstrate a notable increase in average er-
tors described by Eq. (37). Furthermore, as indicated gotropy, aligning with the theoretical predictions. Fur-
in Eq. (19), the QB’s post-measurement states are di- thermore, we can observe that the deviation between
agonalized under the energy eigenstates. Consequently, the experimental and theoretical results is correlated by
we can only measure Q in the z-direction to determine the circuit size, primarily determined by the number of
the stored energy as well as the ergotropy. We utilize two- qubit gates involved. Therefore, comparing the
the ibmq algiers and IonQ-Aria 1 devices. Note that the results from the IonQ and IBMQ devices, we find that the
qubit configuration for ibmq algiers is illustrated in Fig. 5 ex- perimental data from the IonQ device exhibit a closer
(c), while the qubits in IonQ-Aria 1 are fully connected. match to the theoretical curves compared to those from
Also, since the Ising coupling gates are not native gates the IBMQ device. In addition, the errors associated with
for the IBMQ device, we need to further decompose the single-charger protocol are smaller in magnitude
them into CNOT gates, which is shown in Fig. 5 (d). than those of the multiple-chargers protocol.
Conse- quently, the circuits for the IBMQ and IonQ
devices con- sist of 20 and 12 two-qubit gates,
respectively. Here, we use the Hadamard (H) gate and V. SUMMARY AND OUTLOOK
the rotation-z [Rz(θ)] gate, which are defined by
In this work, we utilize superposition of trajectories
H= to propose two charging protocols for quantum batteries
,
2 1 11 −11 (QBs), leading to improved storage of extractable work.
The first protocol, called the multiple-charger protocol,
(43)
Rz(θ) = exp(−iσˆ z θ/2). allows for simultaneous interaction between the QB and
multiple chargers. Leveraging the interference effect, we
For the single-charger setup, the circuit model for the observe an increase in extractable work, attributed to
charging process is presented in Fig. 6, which consists of the Zeno-like state freezing effect. The second protocol,
6 and 4 two-qubit gates in the circuits for the IBMQ and re- ferred to as the single-charger protocol, enables the
the IonQ devices, respectively. QB to interact with a single charger from multiple
Figures 7 and 8 illustrate the results obtained from positions concurrently. We demonstrate that the
the two protocols, with each data point representing the collective inter- ference effect generated by this setup
allows the QB to
14
0 ACKNOWLEDGEMENT
0 π/4 π/2
τ [ωc−1λ−1] We acknowledge the NTU-IBM Q Hub and the IBM
quantum experience for providing us a platform to im-
FIG. 7. The stored energy E and average ergotropy W
plement the experiment. This work is supported by the
(both in units of ¯hω c ) on time τ (in units of 1/ωc) National Center for Theoretical Sciences and National
for λ = 0.05, N = 2. The black solid and dashed curves rep- Science and Technology Council, Taiwan, Grant Nos.
resent the stored energy and average ergotorpy predicted
MOST 111-2123-M-006-001and NSTC 111-2627-M-006-
by numerical simulations. The green circles and blue “x”s 008.
repre- sent experimental results performed on ibmq algiers.
The red triangles and magenta diamonds represent
experimental re- sults performed on IonQ Aria 1. Each data Appendix A: Effects of the choice of projectors on
point is obtained after averaging 1000 experimental the post-measurement quantum battery states
repetitions.
Here, we prove that the choice of projectors for
ξ| k̸=1 ξ⟩ k̸⟨=1 D is irrelevant to the post-measurement
1 E, W [Sim.] states| σk̸=1(τ) as long as they are orthonormal to
E [IBMQ Exp.] |ξk=1 ⟩ ⟨ξk=1 D.
E [IonQ Exp.] Let| us start from Eq. (8), which can be expanded into
W [IBMQ Exp.]
W [IonQ Exp.] σk(τ) =
[h¯ωc]
1 N h i
E, W
Σ (A1)
N ⟨ξk|j⟩ ⟨f|ξk⟩ TrC
j,f =1 |ϕj(τ)⟩ ⟨ϕf (τ)|QC
h¯
Inserting this into Eq. (A1), we can obtain Recall that for the case of k = 1, we require the
Σ
Σ 1 projectors to be orthonormal to̸ that of k = 1, i.e.,
2
σ (τ) = ∞ | (τ)| σn |e⟩ ⟨e| σn |⟨ξ |j⟩ ⟨ξk̸=1|ξk=1⟩ = 0, implying that the ket |ξk̸=1⟩ can be
2 N
α |
k n x x k expressed by
n=0 j=f
N (A6) N
Σ Σ Σ
1 N N
+ |α (τ)|2|g⟩ ⟨g| ⟨ξ |j⟩ ⟨f|ξ ⟩ .
1 k |ξk̸=1⟩ = βk,j |j⟩D ,
N k j̸=f f =1 j=1
N (A9)
For the case of k = 1, according to Eq. (7), we can Σ
obtain with β k,j = 0 ∀ k ̸= 1.
j=1
1
⟨ξk=1|j⟩ = √ ∀ j = 1, · · · , N. (A7)
N
Therefore, the corresponding post-measurement state Therefore, the post-measurement states can be written
can be written as as
∞
∞ Σ
Σ 1 |αn(τ)| σ
2 x n|e⟩ ⟨e| σxn
σk=1(τ) = 1 σk̸=1(τ) =
2
|αn(τ)| σ xn|e⟩ ⟨e| σ xn
N N n=0 (A10)
n=0
N (A8) 1
−1 2
2
[1] R. Uzdin, A. Levy, and R. Kosloff, Equivalence of quan- localized quantum batteries, Phys. Rev. B 100, 115142
tum heat machines, and quantum-thermodynamic signa- (2019).
tures, Phys. Rev. X 5, 031044 (2015). [12] A. Crescente, M. Carrega, M. Sasetti, and D. Ferraro,
[2] W. Niedenzu, V. Mukherjee, A. Ghosh, A. G. Kofman, Charging and energy fluctuations of a driven quantum
and G. Kurizki, Quantum engine efficiency bound battery, New J. Phys. 22, 063057 (2020).
beyond the second law of thermodynamics, Nat. [13] A. Crescente, M. Carrega, M. Sassetti, and D. Ferraro,
Commun. 9, 165 (2018). Ultrafast charging in a two-photon dicke quantum bat-
[3] R. Horodecki, P. Horodecki, M. Horodecki, and tery, Phys. Rev. B 102, 245407 (2020).
K. Horodecki, Quantum entanglement, Rev. Mod. Phys. [14] S. Ghosh, T. Chanda, S. Mal, and A. Sen(De), Fast
81, 865 (2009). charging of a quantum battery assisted by noise,
[4] D. Bruß, Characterizing entanglement, Journal of Math- Phys. Rev. A 104, 032207 (2021).
ematical Physics 43, 4237 (2002). [15] A. C. Santos, B. i. e. i. f. m. c. C¸ akmak, S. Campbell,
[5] A. Streltsov, G. Adesso, and M. B. Plenio, Colloquium: and
Quantum coherence as a resource, Rev. Mod. Phys. 89, N. T. Zinner, Stable adiabatic quantum batteries,
041003 (2017). Phys. Rev. E 100, 032107 (2019).
[6] O. Nairz, M. Arndt, and A. Zeilinger, Quantum interfer- [16] F. Pirmoradian and K. Mølmer, Aging of a quantum bat-
ence experiments with large molecules, American Journal tery, Phys. Rev. A 100, 043833 (2019).
of Physics 71, 319 (2003). [17] J. Liu, D. Segal, and G. Hanna, Loss-free excitonic
[7] F. C. Binder, S. Vinjanampathy, K. Modi, and J. Goold, quan- tum battery, J. Phys. Chem. C 123, 18303
Quantacell: powerful charging of quantum batteries, New (2019).
J. Phys. 17, 075015 (2015). [18] J. Q. Quach and W. J. Munro, Using dark states to
[8] F. Campaioli, F. A. Pollock, F. C. Binder, L. C´eleri, charge and stabilize open quantum batteries, Phys. Rev.
J. Goold, S. Vinjanampathy, and K. Modi, Enhancing the Applied 14, 024092 (2020).
charging power of quantum batteries, Phys. Rev. Lett. [19] S. Gherardini, F. Campaioli, F. Caruso, and F. C. Binder,
118, 150601 (2017). Stabilizing open quantum batteries by sequential mea-
[9] D. Ferraro, M. Campisi, G. M. Andolina, V. Pellegrini, surements, Phys. Rev. Res. 2, 013095 (2020).
and M. Polini, High-power collective charging of a solid- [20] D. Rosa, D. Rossini, G. M. Andolina, M. Polini, and
state quantum battery, Phys. Rev. Lett. 120, 117702 M. Carrega, Ultra-stable charging of fast-scrambling syk
(2018). quantum batteries, J. High Energy Phys. 2020 (11), 67.
[10] T. P. Le, J. Levinsen, K. Modi, M. M. Parish, and F. A. [21] R. Alicki and M. Fannes, Entanglement boost for ex-
Pollock, Spin-chain model of a many-body quantum bat- tractable work from ensembles of quantum batteries,
tery, Phys. Rev. A 97, 022106 (2018). Phys. Rev. E 87, 042123 (2013).
[11] D. Rossini, G. M. Andolina, and M. Polini, Many-body
17
[22] G. Francica, J. Goold, F. Plastina, and M. Paternos-
tro, Daemonic ergotropy: enhanced work extraction from
quantum correlations, npj Quantum Inf. 3, 12 (2017).
18