## Ergotropy and capacity optimization in Heisenberg spin-chain quantum batteries

Asad Ali , 1, ∗ Saif Al-Kuwari , 1, † M. I. Hussain , 1 Tim Byrnes, 2, 3, 4, 5 M. T. ‡

Rahim , 1 James Q. Quach, 6 Mehrdad Ghominejad , 7 and Saeed Haddadi 7,

1 Qatar Centre for Quantum Computing, College of Science and Engineering, Hamad Bin Khalifa University, Doha, Qatar 2 New York University Shanghai, NYU-ECNU Institute of Physics at NYU Shanghai,

Shanghai Frontiers Science Center of Artificial Intelligence and Deep Learning,

567 West Yangsi Road, Shanghai, 200126, China

3 State Key Laboratory of Precision Spectroscopy, School of Physical and Material Sciences,

East China Normal University, Shanghai 200062, China

4 Center for Quantum and Topological Systems (CQTS),

NYUAD Research Institute, New York University Abu Dhabi, UAE

5 Department of Physics, New York University, New York, NY 10003, USA

6 The University of Adelaide, SA 5005, Australia

7 Faculty of Physics, Semnan University, P.O. Box 35195-363, Semnan, Iran

(Dated: November 4, 2024)

This study examines the performance of finite spin quantum batteries (QBs) using Heisenberg XX XY XXZ , , , and XYZ spin models with Dzyaloshinsky-Moriya (DM) and Kaplan-ShekhtmanEntin-Wohlman-Aharony (KSEA) interactions. The QBs are modeled as interacting quantum spins in local inhomogeneous magnetic fields, inducing variable Zeeman splitting . We derive analytical expressions for the maximal extractable work, ergotropy and the capacity of QBs, as recently examined by Yang et al. [Phys. Rev. Lett. 131, 030402 (2023)]. These quantities are analytically linked through certain quantum correlations, as posited in the aforementioned study. Different Heisenberg spin chain models exhibit distinct behaviors under varying conditions, emphasizing the importance of model selection for optimizing QB performance. In antiferromagnetic (AFM) systems, maximum ergotropy occurs with a Zeeman splitting field applied to either spin, while ferromagnetic (FM) systems benefit from a uniform Zeeman field. Temperature significantly impacts QB performance such that ergotropy in the AFM case is generally more robust against increasing temperature compared to the FM case. Incorporating DM and KSEA couplings can significantly enhance the capacity and ergotropy extraction of QBs. However, there exists a threshold beyond which additional increases in these interactions cause a sharp decline in capacity and ergotropy. This behavior is influenced by temperature and quantum coherence, which signal the occurrence of a sudden phase transition. The resource theory of quantum coherence proposed by Baumgratz et al. [Phys. Rev. Lett. 113, 140401 (2014)] plays a crucial role in enhancing ergotropy and capacity. However, ergotropy is limited by both the system's capacity and the amount of coherence. These findings support the theoretical framework of spin-based QBs and may benefit future research on quantum energy storage devices.

Keywords: Quantum battery, Ergotropy, Capacity, Spin chain, Quantum coherence

## I. INTRODUCTION

Ongoing battery research faces critical challenges, such as low energy density, slow charging speeds, limited shelf life, and environmental concerns. This introduces a pressing demand for innovative solutions to address these issues [1, 2]. The operation and architecture of chemical batteries are purely classical, completely disregarding the underlying useful quantum resource with the potential to leverage a new paradigm in the energy sector. Recent theoretical proposals on quantum battery (QB), based on the principles of quantum mechanics, present compelling heuristics [3-5]. These proposals attempt to harness quantum traits such as superposition and quantum correlations [6-9] to tackle the challenges faced by conventional classical batteries, e.g., by offering superior energy densities and fast charging rates. This could potentially offer crucial insights into quantum thermodynamics, particularly in overcoming challenges related to environmental-induced superselection (einselection) in quantum mechanics [10, 11]. Recent findings also emphasize the importance of QBs for truly reversible quantum gates on quantum computers [12]. Moreover, studying QBs can illuminate the quantum mechanical origins of thermodynamic quantities such as heat, internal energy, ergotropy, and entropy.

∗ asal68826@hbku.edu.qa

† smalkuwari@hbku.edu.qa

‡ haddadi@semnan.ac.ir

Although experimental realization of QBs in the literature is scarce, a recent study [13] reported experimental findings showing that a QB uses superconducting qutrits, which are optimized to maintain stable charging states. Furthermore, the authors highlighted a selfdischarge mechanism reminiscent of supercapacitors, suggesting efficient energy storage capabilities in superconducting circuits. Similarly, recent experimental tests confirm a quantum advantage in charging QBs using NMR star-topology spin systems [14] and verified the charging

FIG. 1. (a) Non-interacting QC-based QB with global charging field provided by the common cavity field without coupling between cells. (b) Interacting QC-based QB with local magnetic field provided by the field stored inside each cavity with couplings between cells.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00133v2-figure-698c7b88-a602-4e58-89da-d3a924410eac.jpg)

behavior in organic QBs [15]. Moreover, the authors in [16] realized Dicke QBs, showcasing superextensive power storage capabilities in their experimental setup. More recently, QB charging utilizing single photons within a linear optics setup has been demonstrated [17].

QBs are fundamentally realized through two distinct designs utilizing quantum cells (QCs) as their basic units [18], as depicted in Fig. 1. The first design, the isolated QC-based QB, conceptualizes each QC as an independent, non-interacting entity, all subjected to a uniform global field for charging. These QCs only remain entangled during the charging process via a global entangled charger field, as shown in Fig. 1(a). This approach primarily focuses on the intrinsic properties of each QC, such as quantum superposition and a non-degenerate energy spectrum. In contrast, the second design, the interacting QC-based QB, illustrated in Fig. 1(b), considers each QC as part of a system where inter-QC couplings are pivotal in determining the QB's overall performance in terms of energy storage and retrieval, acknowledging that no quantum system operates in isolation. This design evolves from isolated systems to a coupled many-body quantum system. QBs can be modeled using spin-chain models, which represent many-body quantum systems with couplings between constituent spins that function as QCs. Unlike non-interacting QC-based QBs, the interacting QC-based QB model necessitates at least two QCs or spins and emphasizes the importance of interspin couplings, as seen in one-dimensional spin chains with realistic two-body couplings [19-23].

Heisenberg spin chains (HSCs), which serve as foundational models in condensed matter physics and quantum information technologies, exhibit a variety of multipartite couplings and potential non-classical correlations such as quantum entanglement and quantum discord. These HSCs have been used extensively in quantum technologies including quantum transistors [24], cold atoms [25], interferometry [26], optical lattices [27], and quantum communications [28-42]. Since the introduction of the QB concept, researchers have been actively pursuing the- oretical and experimental implementations across diverse quantum platforms, with the spin chain-based QB emerging as particularly promising [43-65]. By leveraging spin chains as QBs, we can examine the interplay between exchange couplings and the charging process, potentially unlocking quantum advantages in terms of power and efficiency. Unlike earlier studies that focused on global entangling operations during charging, spin chain-based QBs facilitate an in-depth analysis of the role of local charging fields and the contribution of intra-system couplings to the charging process.

In this paper, we focus on QBs based on the HSC with antisymmetric Dzyaloshinsky-Moriya (DM) exchange coupling and symmetric KaplanShekhtman-EntinWohlman-Aharony (KSEA) coupling [20, 22, 66, 67] along with variable Zeeman splitting, which has not been addressed in the literature yet to the best of our knowledge. There are different sub-categories of 1D Heisenberg spin models, such as the Ising, XX XY , , XXX , XXZ , XYZ , and different research groups have addressed these specific cases of HSC models with intrinsic spin-spin couplings without DM and KSEA couplings [21, 43-45, 56, 68-70]. Therefore, our major emphasis is on exploiting DM and KSEA couplings for extracting ergotropy production since in previous papers the spinspin coupling in different subcategories of HSC such as Ising, XX XY , , XXZ XYZ , have been reported. By exploring the spin-chain QB with added DM and KSEA couplings, we initiate a new avenue for experimentalists to design experimental schemes or working substances for spin chain-based efficient QBs, paving the way for the development of efficient and practical QBs tailored for quantum technologies. The rest of this paper is organized as follows: Sec. II covers the modeling of QBs in a spin-1 / 2 chain, their initialization, and the charging process. It also introduces the performance indicators for QBs, including work done, ergotropy, power, and efficiency. Sec. III presents the results and discussion. Finally, Sec. IV concludes this work.

## II. MODELING QUANTUM BATTERIES WITH HEISENBERG SPIN CHAIN

## A. Classes of Heisenberg spin models

A spin-based QB is typically a system of quantum mechanical spins' collection with associated couplings to each other. The general 1D Heisenberg spin-1 / 2 chain with nearest neighbor coupling can be modeled via a Hamiltonian, given by

$$\mathfrak { e } _ { \substack { \mathfrak { e } \\ \mathfrak { l } = \\ \mathfrak { e } = \\ \mathfrak { h } = \\ \mathfrak { h } = \\ \mathfrak { h } = \\ \mathfrak { h } = \\ \mathfrak { h } = \\ \mathfrak { h } = \\ \mathfrak { h } = \\ \mathfrak { h } = \\ \mathfrak { h } = \\ \mathfrak { h } = \\ \mathfrak { h } = \\ \mathfrak { h } = \\ \mathfrak { h } = \\ \mathfrak { h } = \\ \mathfrak { h } = \\ \mathfrak { h } = \\ \mathfrak { h } = \\ \mathfrak { h } = \\ \mathfrak { h } = \\ \mathfrak { h } = \\ \mathfrak { h } = \\ \mathfrak { h } = \\ \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfra \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfring { \mathfrak { h } \\ \mathfrak { h } = \\ \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \\ \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathi \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathfrak { h } \mathf ]$$

where γ is called the anisotropy parameter, and J and ∆ are the nearest neighbor coupling parameters in the xy -and z -plane respectively. When both J and ∆ are positive, this model exhibits antiferromagnetic (AFM) behavior, whereas when both J and ∆ are negative, it demonstrates ferromagnetic (FM) characteristics. The standard Pauli spin operators denoted as ˆ σ k i (for k = x, y, z ) define the spin vector on the i th site as -→ σ i = (ˆ σ , σ x i ˆ y i , σ ˆ z i ). The different ranges of γ and ∆ determine how the model falls into different categories of the onedimensional Heisenberg model, as shown in Table I.

In this paper, we focus on four out of the six potential cases. Specifically, we analyze the XX XY , , XXZ , XYZ models. This selection allows us to thoroughly investigate the key dynamics and characteristics pertinent to our research objectives.

For a more comprehensive understanding of the other cases not considered in this study, we direct the reader to the interesting work presented in [18, 43, 45, 46, 51, 55, 56, 66, 71]. We note that, to date, no research has incorporated DM and/or KSEA coupling(s) in the QBs based on various subclasses of HSCs.

## B. The QB Hamiltonian

The Hamiltonian ˆ H S for working substance of our system reads

$$\hat { \mathcal { H } } _ { S } & = \frac { 1 } { 4 } \sum _ { j = 1 } ^ { N - 1 } \left \{ J \left [ ( 1 + \gamma ) \hat { \sigma } _ { j } ^ { x } \otimes \hat { \sigma } _ { j + 1 } ^ { x } + ( 1 - \gamma ) \hat { \sigma } _ { j } ^ { y } \otimes \hat { \sigma } _ { j + 1 } ^ { y } \right ] \quad & \text{behav} \\ & + \Delta \hat { \sigma } _ { j } ^ { z } \otimes \hat { \sigma } _ { j + 1 } ^ { z } + \overrightarrow { \mathcal { D } }. ( \overrightarrow { \hat { \sigma } } _ { j } \times \overrightarrow { \hat { \sigma } } _ { j + 1 } ) + \overrightarrow { \hat { \sigma } } _ { j }. \overleftarrow { \Gamma }. \overrightarrow { \sigma } _ { j + 1 } \right \}, \\ & \quad \sim.$$

where the DM coupling is denoted by -→ D = ( D ,D ,D x y z ) and the KSEA coupling ← → Γ is a symmetric traceless tensor which can be written as

$$\leftrightarrow \left ( \begin{matrix} 0 & G _ { z } & G _ { y } \\ G _ { z } & 0 & G _ { x } \\ G _ { y } & G _ { x } & 0 \end{matrix} \right ). \quad \quad \left ( 3 \right ) \quad \text{$\text{$this$k$} } \\ \text{$\text{on$this$} }$$

The DM coupling, stemming from spin-orbit interaction, introduces an antisymmetric exchange coupling between neighboring spins, which is crucial for understanding phenomena like weak ferromagnetism [72-75]. Simultaneously, the often-overlooked KSEA coupling restores symmetry disturbed by the DM term, providing a more accurate depiction of spin dynamics [76]. Neglecting these couplings would obscure crucial aspects of quantum coherence and correlation, hindering our understanding of spin systems under external magnetic fields. Thus, their inclusion forms a fundamental basis for advancing QB technologies and deepening our understanding of correlated quantum materials-based QBs where DM and KSEA couplings are dominant [77, 78].

TABLE I. Different cases of 1D Heisenberg spin models.

̸

̸

| Anisotropy in xy -plane ( γ   | z -axis coupling (∆)   | Case   |
|-------------------------------|------------------------|--------|
| ± 1                           | 0                      | Ising  |
| 0                             | 0                      | XX     |
| 0                             | = 0                    | XXZ    |
| 0                             | J                      | XXX    |
| 0 < γ < 1                     | 0                      | XY     |
| 0 < γ < 1                     | = 0                    | XYZ    |

In order to guarantee the non-degeneracy in the energy state of the working substance to act as QB, we apply Zeeman splitting to the spins using the following Hamiltonian

$$\text{ier} \\ \underset { \text{as} } { \text{ier} } \\ \text{is}$$

where -→ B j = ( B ,B ,B x j y j z j ) represents the magnetic field applied in an arbitrary direction to induce Zeeman splitting at any j th and ( j +1)th site. We consider inhomogeneous magnetic fields along the chain, necessitating a site-dependent Hamiltonian formulation to capture this variation accurately. This approach is crucial for comprehensively understanding the system's dynamics and behaviors.

Having the Hamiltonian of the working substance (2) and the Zemaan splitting Hamiltonian (4), we can establish the Hamiltonian for a QB as

$$\hat { \mathcal { H } } _ { \mathcal { Q } } = \hat { \mathcal { H } } _ { \mathcal { S } } + \hat { \mathcal { H } } _ { \mathcal { Z } }.$$

̸

Given the intricate nature of the model, characterized by many parameters, maintaining complete generality becomes daunting. We add that not all subcategories of Heisenberg spin models are exactly solvable such that the analytical spectrum of energies and eigenstates is known. Thus, we opt for a tractable approach, focusing on a scenario involving two neighboring spins. Within this framework, we set the magnitude of the magnetic field along x and y directions to be zero and investigate the influence of a Zeeman splitting due to an inhomogeneous sinusoidal magnetic field acting along the conventional z -direction, denoted as -→ B 1 = (0 0 , , B cos θ ) and -→ B 2 = (0 0 , , B sin θ ). Let us define 0 ≤ θ ≤ π 2 such that when θ = 0 (or equivalently θ = π 2 ), the Zeeman splitting is applied to only one spin. When θ = π 4 , the magnetic field is applied equally to both spins, and the same amount of Zeeman splitting is applied to both QCs. Furthermore, we consider the DM and KSEA couplings to be aligned with the z -direction, specified by -→ D = (0 0 , , D z ), with the constraint G x = G y = 0 and G z = 0. Consequently, Eq. (5) through Eqs. (2) and (4) can be written to accommodate a pair of two spins, namely

$$\hat { \mathcal { H } } _ { \mathcal { Q } } & = \frac { J } { 4 } [ ( 1 + \gamma ) \hat { \sigma } _ { 1 } ^ { x } \otimes \hat { \sigma } _ { 2 } ^ { x } + ( 1 - \gamma ) \hat { \sigma } _ { 1 } ^ { y } \otimes \hat { \sigma } _ { 2 } ^ { y } ] \\ & + \Delta \hat { \sigma } _ { 1 } ^ { z } \otimes \hat { \sigma } _ { 2 } ^ { z } + B \cos \theta ( \hat { \sigma } _ { 1 } ^ { z } \otimes \hat { \mathbb { I } } _ { 2 } ) + B \sin \theta ( \hat { \mathbb { I } } _ { 2 } \otimes \hat { \sigma } _ { 2 } ^ { z } ) \\ & + D _ { z } ( \hat { \sigma } _ { 1 } ^ { x } \otimes \hat { \sigma } _ { 2 } ^ { y } - \hat { \sigma } _ { 1 } ^ { y } \otimes \hat { \sigma } _ { 2 } ^ { x } ) + G _ { z } ( \hat { \sigma } _ { 1 } ^ { x } \otimes \hat { \sigma } _ { 2 } ^ { y } + \hat { \sigma } _ { 1 } ^ { y } \otimes \hat { \sigma } _ { 2 } ^ { x } ). \quad \text{with}$$

## C. The initial state of QB

In the context of QBs, we may consider two initial states: the uncharged state and the charged state. For the uncharged QB's initial condition, it can either be in its ground state | 00 ⟩⟨ 00 , | representing absolute zero temperature, or in a thermal state with finite temperature [45, 79]. In the latter case, the system is initially in its ground state, and through a unitary operator, it undergoes a transition to an excited state, effectively fully charging the battery over time for subsequent power or ergotropy extraction. Since in the real world no system can reach absolute zero, as the third law of thermodynamics and the Heisenberg uncertainty principle do not allow, we choose the second option and find Gibbs thermal state at equilibrium temperature T from the Hamiltonian expressed in Eq. (6). Hence, we get

$$\hat { \varrho } _ { t h } = \frac { e ^ { - \hat { \mathcal { H } } _ { \mathfrak { Q } } / k _ { B } T } } { Z } = \frac { 1 } { Z } \sum _ { \mu = 1 } ^ { 4 } \exp \left ( - \nu _ { \mu } / k _ { B } T \right ) | \psi _ { \mu } \rangle \left \langle \psi _ { \mu } | \,. \quad \text{$\quad$\quad$(\quad$(7)$\quad$\ \text{energy}$} \\ \tau _ { \mu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \nu \numeric}$$

In the above expression, Z = Tr[exp( -H ˆ Q /k B T )] is the partition function and ν µ stands for the eigenenergy corresponding to the eigenstate | ψ µ ⟩ . Here, k B signifies the Boltzmann constant (set to one for convenience) and T is the absolute temperature. In the matrix representation, Eq. (7) can be written as below

$$\hat { \varrho } _ { t h } = \frac { 1 } { Z } \left ( \begin{array} { c c c c } \varrho _ { 1 1 } & 0 & 0 & \varrho _ { 1 4 } \\ 0 & \varrho _ { 2 2 } & \varrho _ { 2 3 } & 0 \\ 0 & \varrho _ { 3 2 } & \varrho _ { 3 3 } & 0 \\ \varrho _ { 4 1 } & 0 & 0 & \varrho _ { 4 4 } \end{array} \right ), \quad \quad ( 8 ) \quad \text{init} \\ \text{pur} \\ \text{rate} \\ \text{red} \end{array}$$

and its nonzero elements read (setting B = 1)

$$\i d i t s \, \text{nonzero elements read (setting $B=1$)} \\ \varrho _ { 1 1 } = \frac { e ^ { ( \mathcal { E } - \Delta ) / T } [ \mathcal { E } - \sqrt { 2 } \sin ( \theta + \pi / 4 ) ] } { 2 \mathcal { E } } \\ + \frac { e ^ { - ( \Delta + \mathcal { E } ) / T } [ \mathcal { E } + \sqrt { 2 } \sin ( \theta + \pi / 4 ) ] } { 2 \mathcal { E } },$$

$$\varrho \\ \varrho$$

$$& \text{and its nonzero elements read (setting $B=1$} ) & \text{satisfy} \\ & \varrho _ { 1 1 } = \frac { e ^ { ( \mathcal { E } - \Delta ) / T } [ \mathcal { E } - \sqrt { 2 } \sin ( \theta + \pi / 4 ) ] } { 2 \mathcal { E } } \\ & \quad + \frac { e ^ { - ( \Delta + \mathcal { E } ) / T } [ \mathcal { E } + \sqrt { 2 } \sin ( \theta + \pi / 4 ) ] } { 2 \mathcal { E } }, \\ & \varrho _ { 2 2 } = e ^ { \Delta / T } \left ( \frac { \sqrt { 2 } \sinh ( \mathcal { R } / T ) \sin ( \theta - \pi / 4 ) } { \mathcal { R } } + \cosh ( \mathcal { R } / T ) \right ), \quad \text{which} \ \! \ \text{For} \! } { \text{tary of} } \\ & \varrho _ { 3 3 } = e ^ { \Delta / T } \left ( \frac { \sqrt { 2 } \sinh ( \mathcal { R } / T ) \cos ( \theta - \pi / 4 ) } { \mathcal { R } } + \cosh ( \mathcal { R } / T ) \right ), \\ & \varrho _ { 4 4 } = \frac { e ^ { ( \mathcal { E } - \Delta ) / T } [ \mathcal { E } + \sqrt { 2 } \sin ( \theta + \pi / 4 ) ] } { 2 \mathcal { E } } \\ & \quad - \frac { e ^ { - ( \Delta + \mathcal { E } ) / T } [ - \mathcal { E } + \sqrt { 2 } \sin ( \theta + \pi / 4 ) ] } { 2 \mathcal { E } } \quad \text{where} \ \ - \frac { i } { 2 } \sin ( \ \ }$$

$$\varrho _ { 1 4 } = \varrho _ { 4 1 } ^ { * } = \frac { 2 i \chi e ^ { - \Delta / T } \sinh ( \mathcal { E } / T ) } { \mathcal { E } }, \\ \varrho _ { 2 3 } = \varrho _ { 3 2 } ^ { * } = - \frac { 2 i \phi ^ { * } e ^ { \Delta / T } \sinh ( \mathcal { R } / T ) } { \mathcal { R } },$$

$$\begin{matrix} \end{matrix}$$

with

$$\stackrel { \prime \cdot } { \Phi } \\ Z & = 2 \left \{ \cosh \left ( \frac { \Delta } { T } \right ) \left [ \cosh \left ( \frac { \mathcal { E } } { T } \right ) + \cosh \left ( \frac { \mathcal { R } } { T } \right ) \right ] \\ & \quad + \sinh \left ( \frac { \Delta } { T } \right ) \left [ \cosh \left ( \frac { \mathcal { R } } { T } \right ) - \cosh \left ( \frac { \mathcal { E } } { T } \right ) \right ] \right \}, \quad ( 9 )$$

where R = √ 4 D 2 z -sin(2 θ ) + 4 J 2 +1, ϕ = D z + iJ , E = √ 1 + 4 G 2 z +4 J γ 2 2 +sin(2 ), and θ χ = G z + iγJ .

## D. Charging the QB

The QB can be charged by applying a local external magnetic field either in the x -direction or the y -direction, with the charging field strength Ω = Ω( ), which can be t time-dependent or time-independent. For simplicity, we assume a time-independent charging field. The charging Hamiltonian is given by:

$$\hat { \mathcal { H } } _ { \mathcal { C } } ^ { ( x, y ) } = \Omega \left ( \hat { \sigma } _ { 1 } ^ { ( x, y ) } \otimes \hat { \mathbb { I } } _ { 2 } + \hat { \mathbb { I } } _ { 2 } \otimes \hat { \sigma } _ { 2 } ^ { ( x, y ) } \right ). \quad ( 1 0 )$$

The magnetic field is applied uniformly to all qubits in the system. As the battery reaches its maximum stored energy at time t = τ , the charging field must be disconnected to prevent the battery from reverting back to an uncharged state due to continuous cycling.

Both the PauliX and PauliY gates can equivalently serve as charging operators for the QB. The PauliX gate acts as a bit-flip, while the PauliY gate performs both bit-flip and phase-flip operations. One can use either PauliX or PauliY gates for charging because both act as unitary operations that drive the QB system from its initial state to a charged state. The PauliX gate is a pure bit-flip operator, while the PauliY gate incorporates both bit-flip and phase-flip, providing more versatility depending on the QB's configuration and field interactions.

In a closed system, the charging process is implemented as a unitary evolution:

$$\hat { U } _ { \mathcal { C } } ( t ) = \exp ( - i \hat { \mathcal { H } } _ { \mathcal { C } } t ), \text{ \quad \ \ } ( 1 1 )$$

which governs the cyclic charging of the QB over time.

For the PauliY gate charging Hamiltonian, the unitary operator takes the form:

$$\hat { \mathcal { U } } _ { c } ^ { ( y ) } ( t ) = \left ( \begin{array} { c c c c } \alpha & \lambda & \lambda & \beta \\ \lambda & \alpha & \beta & \lambda \\ \lambda & \beta & \alpha & \lambda \\ \beta & \lambda & \lambda & \alpha \end{array} \right ), \quad \quad ( 1 2 )$$

where α = cos 2 (Ω ), t β = -sin 2 (Ω ), t and λ = -i 2 sin(2Ω ). t

Similarly, for the PauliX gate charging Hamiltonian, the unitary operator is:

$$\hat { U } _ { \mathcal { C } } ^ { ( x ) } ( t ) = \left ( \begin{array} { c c c c } \alpha & \vartheta _ { - } & \vartheta _ { - } & 1 - \alpha \\ \vartheta _ { + } & \alpha & \alpha - 1 & \vartheta _ { - } \\ \vartheta _ { + } & \alpha - 1 & \alpha & \vartheta _ { - } \\ 1 - \alpha & \vartheta _ { + } & \vartheta _ { + } & \alpha \end{array} \right ), \quad ( 1 3 ) \quad \text{(12)} \quad \text{a so-} \quad \text{is un} \quad \text{sive $\varepsilon$} \quad \text{.}$$

where ϑ ± = ± 1 2 sin(2Ω ). t

During the charging process, the evolution of the system is determined solely by the charging Hamiltonian ˆ H ( x,y ) C , while the Hamiltonian ˆ H Q is relevant before and after charging.

## E. Performance indicators for QBs

Let ˆ be ϱ a quantum state of QB with Hamiltonian ˆ H Q . The question we are interested in is: how much work can be extracted from this QB through a cyclic unitary process? In this context, cyclicity implies that the Hamiltonian of the QB must be the same at the start and end of the process, i.e., ˆ H Q = ˆ H Q (0) = ˆ H Q ( τ ) [80]. Since the evolution is unitary, any change in the internal energy ⟨H ˆ Q ⟩ is attributed to the work done on or by QB.

To analyze the QB's Hamiltonian, we express it in terms of its increasing spectral decomposition, namely

$$\hat { \mathcal { H } } _ { \mathcal { Q } } \coloneqq \sum _ { \mu = 1 } ^ { 4 } \nu _ { \mu } | \psi _ { \mu } \rangle \langle \psi _ { \mu } |, \text{ where } \nu _ { \mu + 1 } \geq \nu _ { \mu } \, \forall \mu. \quad ( 1 4 ) \quad \max _ { t } \max _ { \tau }$$

In this expression, ν µ denotes the energy eigenvalues, and | ψ µ ⟩ represents the corresponding eigenstates. This formulation enables the evaluation of energy variations and potential for work extraction in a cyclic unitary process, where the QB's Hamiltonian remains unchanged at the start and end of the cycle.

The state ˆ, ϱ on the other hand, is expressed in its decreasing eigen-decomposition

$$\hat { \varrho } \coloneqq \sum _ { \kappa = 1 } ^ { 4 } r _ { \kappa } | r _ { \kappa } \rangle \langle r _ { \kappa } |, \text{ with } r _ { \kappa + 1 } \leq r _ { \kappa } \, \forall \kappa. \quad ( 1 5 ) \quad \overset { m \kappa } { \text{ergo} }$$

The goal is to transform ˆ into a state with lower inϱ ternal energy, extracting the difference in internal energy in the process.

After maximal cyclic, unitary work extraction, no further work can be extracted, and the system ends up in a so-called passive state ˆ [45, 80-82]. π A passive state is unique up to degeneracy in the Hamiltonian. Passive states are diagonal in the Hamiltonian's eigenbasis, with decreasing populations for increasing energy levels [80, 82]. That is, a state ˆ, ϱ as defined above, is passive if | r n ⟩ = | ψ n ⟩ ∀ n . Gibbs thermal states are consequently passive [80]; therefore, we set ˆ = ˆ π ϱ th .

The maximum work that can be extracted from a nonpassive state ˆ for a Hamiltonian ϱ ˆ H Q via a cyclic unitary process ˆ ϱ th → ˆ (ˆ = ϱ ϱ ˆ U ϱ C ˆ th ˆ U † C ) [45], called ergotropy, is found to be [82, 83]

$$\xi ( t ) \coloneqq \sum _ { m, n } r _ { m } \nu _ { n } \left [ | \langle \psi _ { n } | r _ { m } \rangle | ^ { 2 } - \delta _ { m n } \right ], \quad ( 1 6 )$$

where δ mn is the Kronecker delta function. Equivalently, one can write the above expression as [51]

$$\xi ( t ) \coloneqq \text{Tr} [ ( \hat { \varrho } - \hat { \varrho } _ { t h } ) \hat { \mathcal { H } } _ { \mathcal { Q } } ].$$

If the final state is another non-passive state with relatively higher energy than ˆ, say Π, then work done by π ˆ the QB can be defined as [45]

$$\mathcal { W } ( t ) = \text{Tr} [ ( \hat { \varrho } - \hat { \Pi } ) \hat { \mathcal { H } } _ { \mathcal { Q } } ].$$

Notice that the average power is defined as ⟨ p t ( ) ⟩ = W ( ) t /t , and the peak average power is ⟨ p max ( ) t ⟩ = max t W ( ) t /t [45]. Besides, the efficiency η = W ( ) t /ξ t ( ) quantifies the fraction of stored energy W ( ) that can be t extracted as ergotropy ξ t ( ). Efficiency is typically constrained to η ≤ 1, with η = 1 indicating W ( ) = t ξ t ( ) for which the initial state should be a passive state.

## III. RESULTS AND DISCUSSION

The study of temporal ergotropy evolution in various Heisenberg model subclasses, detailed in Section II, provides key insights into the charging and discharging dynamics of spin-based QB systems, facilitating their output optimization. The closed-form expression for ergotropy based on our considered model with Y -gate charging Hamiltonian is given by

$$\xi ( t ) = & \frac { 2 \sin ^ { 2 } ( \Omega t ) } { \sqrt { a b } \, ( 2 h c _ { a } + d ^ { 2 } + 1 ) } \left \{ c \left [ - 2 \sqrt { a b } h \epsilon _ { a } + \sqrt { a } \cos ( 2 \Omega t ) \left ( - 2 \sqrt { b } d g \epsilon _ { a } + \sqrt { b } \, ( d ^ { 2 } + 1 ) + 2 \gamma \, ( d ^ { 2 } - 1 ) \, J \right ) + 2 \sqrt { a b } d \epsilon _ { b } \right ] \\ & + 2 \sqrt { b } f _ { a } \left [ h ( a + 2 J ( c - 2 J ) ) + 2 c d g \, J \cos ( 2 \Omega t ) + \sin ( 2 \theta ) ( h - d g ) \right ] + 2 \sqrt { a d } f _ { b } \left [ b + 2 \gamma \, J ( c - 2 \gamma \, J ) \right ] \right \},$$

where a = R 2 , b = E 2 , c = -∆+ γJ + J , d = e √ b T , ϵ a = cosh ( √ a T ) , ϵ b = cosh ( √ b T ) , f a = sinh ( √ a T ) , f b = sinh ( √ b T ) , g = e 2∆ T , and h = e 2∆+ √ b T .

## A. Effect of variable Zeeman splitting

We investigate the temporal evolution of ergotropy (19) under an inhomogeneous magnetic field parameterized by θ in Eq. (6), focusing on QB performance over time Ω t and θ , with DM and KSEA couplings set to zero for various Heisenberg model subclasses. Setting B = 1 and 0 ≤ θ ≤ π 2 , Zeeman splitting affects one spin at θ = 0 or θ = π 2 , while both spins experience equal splitting at θ = π 4 . We constrain Ω t to 0 ≤ Ω t ≤ 2 π , allowing the QB to undergo multiple charging and discharging cycles at T = 0 1, illustrating cyclic unitary charging operation. . Acommon strategy, often adopted in experimental protocols, involves deactivating the NOT-gate-based charging field once maximum ergotropy is reached. This prevents the quantum NOT-gate-based charging Hamiltonian, as described in Eq. (10), from returning the system to its thermal state. Without deactivation, the system would oscillate between charged and passive states over time. It is important to note that this approach is well-established and not a novel execution of this work.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00133v2-figure-9f5d048a-15a5-48b1-8fb1-afd33ed7da39.jpg)

The density plots in Fig. 2 show the variation of ergotropy ξ t ( ) in the AFM scenario ( J = 1), with γ ≥ 0, ∆ ≥ 0, D z = 0, G z = 0, and T = 0 1, . as functions of time Ω t and Zeeman splitting inhomogeneity parameter θ . Maximum ξ t ( → τ ) occurs when the magnetic field affects one spin at θ = 0 or θ = π 2 , while the lowest ξ t ( ) is observed when both spins are equally affected at θ = π 4 . Although trends appear similar across HSC cases, the XY model exhibits distinct behavior around θ = π 2 , maintaining nonzero ξ t ( ) at θ = π 4 and suggesting slightly better extractable work compared to the XX XXZ , , and XYZ models. Within the interval 0 ≤ Ω t ≤ 2 π , the QB experiences two complete chargedischarge cycles. The maximal charge states are denoted by red regions along the Ω -axis. t For the given set of parameters, the optimal charging time is determined to be τ = π 2Ω . The peak value of ξ max ( ) remains invariant in t time due to the unitary nature of the charging process.

We have analyzed the effect of inhomogeneous Zeeman splitting on the QB in the AFM case over time. Next, we will investigate the FM case. The density plots in Fig. 3 show the temporal evolution of ξ t ( ) as functions of Ω and t θ in the FM scenario, where J = -1 and ∆ ≤ 0 with γ ≥ 0, D z = 0, G z = 0, and T = 0 1. . Obviously, the results in the FM case differ significantly from those in the AFM case.

Unlike the AFM case illustrated in Fig. 2, where the highest values of ξ t ( → τ ) are attained at τ = π 2Ω for all considered HSC cases and at θ = 0 or π 2 , the FM scenario presents a distinct behavior. In the XX model (Fig. 3(a)) and the XY model (Fig. 3(b)), the peak value of ergotropy, or first full-charging, occurs at τ = π 4Ω . This

FIG. 2. Density plots of ξ t ( ) in the AFM scenario versus Ω t and θ , with D z = 0, G z = 0, T = 0 1, and . J = 1: (a) XX (∆ = 0, γ = 0); (b) XY (∆ = 0, γ = 0 5); . (c) XXZ (∆ = 0 5, . γ = 0); (d) XYZ (∆ = 0 5, . γ = 0 5). .

peak is reached four times within the interval 0 ≤ Ω t ≤ 2 π . This shows that, in the FM case, the charging rate for the XX and XY models in Figs. 3(a)-(b) is two times faster than in the AFM scenario for all the HSC cases depicted in Figs. 2(a)-(d). For the AFM case, the peak value of ergotropy is achieved for the Zeeman splitting parameter at θ = 0 (or equivalently π 2 ), whereas in the FM scenario, the ergotropy value remains nonzero for all values of 0 ≤ θ ≤ π 2 , with maxima at θ = 0 and θ = π 2 . Thus, the inhomogenous Zeeman splitting in the FMcase with XX and XY models does not substantially alter the maximum ergotropy.

Considering XXZ and XYZ models from Figs. 3(c) and 3(d), one can see that unlike Figs. 3(a) and 3(b) where the impact of Zeeman splitting parameter θ is roughly independent on ergotropy, here we find that the maximum value of ξ t ( → τ ) occurs around θ = π 4 wheres minimum value but not zero occurs at θ = 0 or equivalently π 2 . Furthermore, in contrast to Figs. 3(a) and 3(b), where the first maximum charge is achieved at t = π/ 4Ω, the first peak of ξ t ( ) in Figs. 3(c) and 3(d) occurs at t = π/ 2Ω.

A key observation is that the overall peak ergotropy in all HSC categories under the FM scenario is nearly three times greater than the peak ergotropy in the AFM case. Notably, the Heisenberg XY model demonstrates greater robustness compared to the other three models for QB in both FM and AFM scenarios.

FIG. 3. Density plots of ξ t ( ) in the FM scenario as functions of Ω t and θ , with D z = 0, G z = 0, T = 0 1, and . J = -1: (a) XX (∆ = 0, γ = 0); (b) XY (∆ = 0, γ = 0 5); . (c) XXZ (∆ = -0 5, . γ = 0); (d) XYZ (∆ = -0 5, . γ = 0 5) .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00133v2-figure-6002ce67-6764-44b0-82b9-d00bdd6d7dfd.jpg)

When optimizing a QB with HSC models, the strategies differ significantly between FM and AFM scenarios. For the FM case, where spin polarizations are parallel, the XX and XY models achieve the fastest charging rates, with peak ergotropy occurring at τ = π 4Ω and repeating four times within 0 ≤ Ω t ≤ 2 π . In contrast, the XXZ and XYZ models reach maximum ergotropy at θ = π 4 and have slower charging, peaking at τ = π 2Ω . This indicates efficient energy transfer in the FM scenarios, but the choice of coupling model affects the timing of optimal charging. For AFM scenario, where spins are oppositely aligned, maximum ergotropy for all models occurs at θ = 0 or π 2 . The XY model shows minimal ergotropy around θ = π 4 , suggesting that the AFM setting benefits from specific spin chain couplings to enhance performance. Thus, according to obtained results from Figs. 2 and 3, we will use θ = π 4 for the FM case and θ = 0 or π 2 for the AFM case in the continuation of our analysis.

## B. Effect of Temperature

Let's study how temperature T affects ξ t ( ) without DM and KSEA couplings in FM and AFM cases.

In the AFM scenario depicted in Fig. 4, the behavior of ξ t ( ) as functions of temperature T and Ω t exhibits a nonmonotonic pattern across all subcategories of the HSC models. For instance, in Fig. 4(a), ξ t ( ) starts from a nonzero value at T = 0, rises as T increases, peaks around T ≈ 1, and then gradually decreases as T continues to rise, eventually approaching zero at asymptot-

FIG. 4. Density plots of ξ t ( ) in the AFM scenario versus Ω t and T , with D z = 0, G z = 0, θ = 0 (or π 2 ), and J = 1: (a) XX (∆ = 0, γ = 0); (b) XY (∆ = 0, γ = 0 5); . (c) XXZ (∆ = 0 5, . γ = 0); (d) XYZ (∆ = 0 5, . γ = 0 5). .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00133v2-figure-3bed3f34-2a98-491a-b811-afea9d2d2693.jpg)

ically high temperatures. A similar trend is observed in Figs. 4(c) and 4(d), where the nonmonotonic dependence of ξ t ( ) on T is also evident. However, in Fig. 4(b), the highly charged regions split from a single nonzero ergotropy region, depicted by the yellowish-red area, into two distinct red regions as T increases. Despite this bifurcation, the overall nonmonotonic effect of T on ξ t ( ) persists. This bifurcation behavior also indicates a faster charging rate in the case illustrated in Fig. 4(b). As previously discussed, in the XY model, ξ t ( ) does not peak at T = 0 but instead starts from a lower value and increases with temperature, reaching its maximum around T ≈ 1. This signifies a higher charging rate compared to the XX XXZ , , and XYZ models, where the influence of temperature is less pronounced. In particular, the charging rate in the XY model is nearly double that of the other three models for given fixed values of parameters. Moreover, the peak value of ergotropy is largest in the XY case ( ξ t ( ) ≈ 1 25), followed by . XX ξ t ( ( ) ≈ 1), with XXZ and XYZ having comparatively lower peak values ( ξ t ( ) ≈ 0 9). .

In contrast to the AFM scenario depicted in Fig. 4, the FM case illustrated in Fig. 5 reveals a significant decrease in ξ t ( ) with increasing T , displaying a monotonic behavior. Unlike the AFM case, where the variation of ξ t ( ) with T is nonmonotonic, the FM scenario shows a steady decline as T rises. This monotonic decrease in ξ t ( ) is attributed to the more fragile nature of quantum coherence and correlations in FM configurations of the HSC models, which are prone to degradation with increasing T compared to the stronger AFM scenario. This fragility is consistently reflected in the ξ t ( ) plots across all T .

FIG. 5. Density plots of ξ t ( ) in the FM scenario versus Ω t and T , with D z = 0, G z = 0, θ = π 4 , and J = -1: (a) XX (∆ = 0, γ = 0); (b) XY (∆ = 0, γ = 0 5); . (c) XXZ (∆ = -0 5, . γ = 0); (d) XYZ (∆ = -0 5, . γ = 0 5). .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00133v2-figure-7acc56ac-3b39-430e-b0bf-abf1119440ea.jpg)

Furthermore, the charging rates differ significantly between the FM and AFM cases. In the FM scenario, the peak charging rates for the XX and XY models are higher than those for the XXZ and XYZ models. The XY model exhibits the highest ξ t ( ), followed by the XX model, with the XXZ model showing lower peak of ξ t ( ), and the XYZ model having the lowest. This is in contrast to the AFM case, where the nonmonotonic behavior allows for temperature-dependent variations that affect the peak values differently. In both scenarios, however, the XY model consistently demonstrates the highest peak of ξ t ( ), highlighting its superior performance in both FM and AFM cases.

Overall, temperature has a more detrimental impact on ergotropy in the FM case across all models compared to the AFM scenario, indicating reduced performance of QBs as temperature increases in the FM configuration. This suggests that QBs operating under FM conditions are more sensitive to thermal effects. For optimal performance, the strategy should focus on maintaining low temperatures and leveraging the specific strengths of each model, whether it be maximizing charging rates or achieving higher peak ergotropy. In particular, models like XY consistently show better ergotropy production, both in FM and AFM cases, making them more suitable for efficient QB performance for given fixed values of other parameters.

FIG. 6. Density plots of ξ t ( ) in the AFM scenario versus Ω t and G z , with D z = 0, T = 0 1, . θ = 0 (or π 2 ), and J = 1: (a) XX (∆ = 0, γ = 0); (b) XY (∆ = 0, γ = 0 5); . (c) XXZ (∆ = 0 5, . γ = 0); (d) XYZ (∆ = 0 5, . γ = 0 5). .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00133v2-figure-2ed107ab-cc9e-47ee-810f-834c8e2c2651.jpg)

## C. Symmetric-antisymmetric exchange coupling

## 1. Effect of KSEA Coupling

In addition to the fact that a high ergotropy in both AFM and FM cases requires low temperatures, it is crucial to apply Zeeman splitting fields effectively. As can be seen from the results, applying the magnetic field to one spin is advantageous in AFM scenarios, while in FM context, applying it to both spins is beneficial. Hence, we investigate how DM and KSEA couplings can be utilized to enhance ergotropy in HSC-based QBs using predefined parameters from previous analyses.

Figure 6 shows the density plots of ξ t ( ) in the AFMcase as functions of Ω t and the z -component of KSEA coupling at T = 0 1, . θ = π 2 , and D z = 0. It shows that G z positively affects all HSC cases. The full-charging occurs around τ = π 2Ω during the first charge. Notably, G z significantly boosts ergotropy in the AFM case with peak values around 75 near G z ≈ 20, as seen in Fig. 6. Beyond G z ≈ 20, ergotropy sharply drops to zero for G z &gt; 20, indicating that G z has maximized the QB's charge, and further increases in G z do not enhance ergotropy. A similar trend is observed in the FM case, which is omitted here for brevity.

## 2. Effect of DM Coupling

The relationship between ξ t ( ), Ω , and t D z for the AFM case with θ = π 2 , G z = 0, and T = 0 1 is shown in Fig. . 7. The density plots show that variations in DM coupling do not affect the charging rate across the cases depicted in

FIG. 7. Density plots of ξ t ( ) in the AFM scenario versus Ω t and D z , with G z = 0, T = 0 1, . θ = 0 (or π 2 ), and J = 1: (a) XX (∆ = 0, γ = 0); (b) XY (∆ = 0, γ = 0 5); . (c) XXZ (∆ = 0 5, . γ = 0); (d) XYZ (∆ = 0 5, . γ = 0 5). .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00133v2-figure-cc8df6ec-5bff-483a-8c41-0e825da6c5eb.jpg)

FIG. 8. Density plots of ξ t ( ) in the FM scenario versus Ω t and D z , with G z = 0, T = 0 1, . θ = π 4 , and J = -1: (a) XX (∆ = 0, γ = 0); (b) XY (∆ = 0, γ = 0 5); . (c) XXZ (∆ = -0 5, . γ = 0); (d) XYZ (∆ = -0 5, . γ = 0 5). .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00133v2-figure-fd884f4c-c62f-49ef-9f61-1577552f2e1c.jpg)

Figs. 7(a) to 7(d). Remarkably, increasing the DM coupling enhances ergotropy values up to a threshold near D z ≈ 35, consistent across all configurations from XX to XYZ . Beyond this threshold, ξ t ( ) suddenly drops down to zero, similar to the effect of G z on ξ t ( ) in Fig. 6. This pattern, where ξ t ( ) increases with D z up to a point and then abruptly drops to zero, persists across all AFM cases, indicating a consistent trend.

In the FM case, as shown in Fig. 8, the variation of ξ t ( ) with D z and Ω t reveals a similar behavior. Like the AFMscenarios where ξ t ( ) drops suddenly after exceeding a certain threshold of D z , the FM case also shows these sudden drops.

To further extend our analysis of why this happens, we focus on the XYZ case as an example of complete anisotropy. Both the FM and AFM scenarios exhibit similar behaviors in the XYZ configuration, including the abrupt transition of ξ t ( ) with changes in D z . Our goal is to determine the maximum ergotropy ξ max achieved by the QB in the AFM scenario and investigate its variation with D z at different temperatures. We maintain the other parameters as in Figs. 7 and 8, and perform the analysis using numerical optimization. Various temperature values have been chosen to show the dependence of ξ max on D z in Figs. 9(a)-(d).

The plots are generated under the AFM scenario with parameters G z = 0, J = 1, ∆ = 0 5, and . γ = 0 5. . These figures illustrate how the maximum ergotropy ξ max varies with D z at different temperatures.

At a low temperature such as T = 0 01 (Fig. . 9(a)), ξ max initially increases with D z up to a threshold value of D z ≈ 2 25, then sharply drops to zero. . The inset in Fig. 9(a) shows a narrow window of oscillations before the value reaches zero, and ξ max does not recover for any further increase even until D z → ∞ . Similar behavior is observed at higher temperatures, with ξ max showing larger values and a higher threshold of D z as shown in Figs. 9(b) and 9(c). The maximum value of ξ max corresponding to the threshold D z increases with temperature.

Recently, a new measure, known as the capacity of QB is proposed [84], defined as

$$K = \text{Tr} [ \hat { \mathcal { H } } _ { \mathcal { Q } } \hat { \varrho } _ { \uparrow } ] - \text{Tr} [ \hat { \mathcal { H } } _ { \mathcal { Q } } \hat { \varrho } _ { \downarrow } ], \text{ \quad \ \ } ( 2 0 )$$

where ˆ ϱ ↓ denotes the lowest energy state of the QB, which in our case is the Gibbs thermal state defined in Eq. (8), and ˆ ϱ ↑ is the highest energy state, ˆ ϱ ↑ = 11 | ⟩⟨ 11 . | It is worth noting that the ergotropy, representing the maximum extractable work from the QB, typically does not equal K , except in the specific cases where the temperature is zero ( T = 0) or when the vacuum state is used instead of a thermal state. The evaluation of the capacity of the QB is independent of time, relying solely on the lowest and highest energy states. This makes it applicable for evaluating both unitary and non-unitary QBs. This measure, which requires no optimization, can be linked to specific entanglement and quantum coherence measures, as demonstrated in [84]. However, its values are not directly comparable to the maximum ergotropy.

To support the trend observed in Fig. 9 for the peak value of ergotropy versus D z , we evaluate K for the same parameter values used in Fig. 10.

The closed-form expression for the capacity of QB turns out to be (with S = sin( θ ) + cos( θ ))

$$\stackrel { \text{sum} } { \text{point} } & \text{$\text{$\text{sum} \colon \infty \colon \infty \colon \infty \colon \infty \colon \infty \colon \infty \colon \infty \colon \infty \colon \infty \colon \infty \colon \infty \colon \infty \colon \infty \colon \infty \colon \infty \colon \infty \colon \infty \colon \infty \colon \infty \colon \infty \colon \infty \colon \infty \colon \infty \colon \infty \colon \infty \colon \infty \colon \infty \colon \infty \colon \infty \colon \infty \colon \infty \colon \infty \colon \infty \colon \infty \text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text	$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text
$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text#$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text}$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text $\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text`$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text'$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text.$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\texts$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text1$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text<$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text^{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$
$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text$$

Figure 10 shows K as a function of D z at various temperatures: T = 0 01 in Fig. . 10(a), T = 0 1 in Fig. . 10(b),

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00133v2-figure-d0a2c3f2-2820-4a84-b2ee-d2835a01152d.jpg)

FIG. 9. Maximum ergotropy ξ max versus D z for the AFM case in the Heisenberg XYZ model with parameters G z = 0, J = 1, ∆ = 0 5, . γ = 0 5, Ω = 1, and . θ = π 2 . The subfigures represent different temperatures: (a) T = 0 01, (b) . T = 0 1, and (c) . T = 1.

FIG. 10. Capacity of QB K versus D z for the AFM case in the Heisenberg XYZ model with parameters G z = 0, J = 1, ∆ = 0 5, . γ = 0 5, Ω = 1, and . θ = π 2 . The subfigures represent different temperatures: (a) T = 0 01, (b) . T = 0 1, and (c) . T = 1.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00133v2-figure-1aecd2d5-d8eb-4f6b-87cc-f28065050eca.jpg)

FIG. 11. Density plots of the maximized l 1 -norm of quantum coherence, Q max (ˆ), ϱ as functions of D z and G z for the AFM Heisenberg XYZ model with J = 1, ∆ = 0 5, . γ = 0 5, Ω = 1, and . θ = π 2 . Subfigures display different temperatures: (a) T = 0 01, (b) . T = 0 1, and (c) . T = 1.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00133v2-figure-e73c0509-669d-45ce-a9f2-8ee2baf5598a.jpg)

and T = 1 0 in Fig. . 10(c). The plots are generated under the AFM scenario with fixed parameters G z = 0, J = 1, ∆ = 0 5, and . γ = 0 5. .

K showing larger values and a higher threshold of D z as shown in Figs. 10(b) and 10(c). The maximum value of K corresponding to the threshold D z increases with temperature.

These plots reveal how K varies with D z at different temperatures. At T = 0 01, . K initially increases with D z up to a threshold value of D z ≈ 2 25, and then sharply . drops to zero. The inset plot in Fig. 10(a) highlights a narrow and sharp oscillation window before K reaches zero and does not recover for any higher values of D z . Similar behavior is observed at higher temperatures, with

We observe a similarity between the capacity of the QB and the peak value of ergotropy with respect to the threshold value of D z , where both measures reach their maximum and then suddenly drop to zero. Although the peak values of these quantifiers differ, they exhibit the same trend concerning the threshold value of D z at var-

ious fixed temperatures. This suggests that the QB has a finite capacity for supplying maximum work, which explains why infinite ergotropy cannot be achieved by continuously increasing D z . The QB's capacity also reflects this constraint, showing that while increasing D z up to a certain threshold can enhance the QB's capacity, it cannot be extended beyond that limit.

To investigate whether the observed signature in ergotropy and QB capacity is linked to an abrupt transition in quantum coherence, we evaluate the l 1 -norm of coherence Q (ˆ) for the AFM case of ϱ XYZ scenario. Specifically, let us consider the following density matrix

$$\hat { \varrho } = \hat { U } _ { \mathcal { C } } \hat { \varrho } _ { t h } \hat { U } _ { \mathcal { C } } ^ { \dagger }, \quad \quad ( 2 2 ) \quad \text{in}$$

where ˆ ϱ th and ˆ U C are defined in Eqs. (8) and (11), respectively.

The Q (ˆ) is defined by [85] ϱ

̸

$$Q ( \hat { \varrho } ) = \sum _ { i \neq j } | \langle i | \hat { \varrho } | j \rangle | = \sum _ { i, j } | \hat { \varrho } _ { i j } | - \sum _ { i } | \hat { \varrho } _ { i i } |, \quad ( 2 3 ) \quad \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{                                                                                                                                                                                                        
\text{ }$$

which is the sum of the absolute off-diagonal elements of a quantum state ˆ in the reference basis. ϱ

Figure 11 illustrates the variations in the maximum l 1 -norm of quantum coherence for the AFM Heisenberg XYZ model versus D z and G z , using the same fixed parameter values as in Figs. 9 and 10 for comparison. Similar to the plots of ξ max and the QB's capacity shown in Figs. 9 and 10, Q max (ˆ) displays a comparable trend. ϱ A sudden drop in Q max (ˆ) ϱ around the same threshold values of DM and KSEA interactions is observed, analogous to the behaviors of ergotropy and capacity. This clearly indicates that the amount of ergotropy is explicitly bounded by the capacity of the QB, and capacity is constrained by the availability of quantum coherence, implying that ergotropy, or the maximum work extraction, is ultimately (implicitly) limited by the QB's coherence.

Furthermore, an increase in temperature requires even larger values of D z to achieve the same maximum quantum coherence ( ≈ 3 as shown in Fig. 11(a)-(c)), as also observed in the QB's capacity plot in Fig. 10. This suggests that indefinitely extracting work from a spin-based QB by merely increasing the values of DM or KSEA interactions is not feasible. The extractable work is constrained by the QB's capacity and ultimately limited by the availability of quantum coherence.

The occurrence of unexpected oscillation followed by non-analytical behavior in ξ max and K in Fig. 9 and Fig. 10, especially within the insets, indicates a kind of sudden phase transition. Initially, the ergotropy ξ t ( ) increases with D z , reflecting enhanced work capacity due to quantum coherence and correlations. However, beyond a critical D z , ξ max abruptly drops to zero, signifying a transition where the essential quantum coherence for work extraction is lost. This behavior is mirrored in the QB's capacity K and the l 1 -norm of coherence Q (ˆ), ϱ which also sharply decline at the same threshold. This consistent pattern suggests a fundamental reorganization of the system's quantum state, highlighting the critical role of quantum coherence in work extraction and underscoring the influence of quantum correlations on the system's energetic properties.

## IV. CONCLUSION AND OUTLOOK

In this paper, we provided a thorough investigation of finite spins quantum batteries based on four different subcategories of Heisenberg spin chains ( XX XY , , XXZ XYZ , ), integrating the variable Zeeman splitting and DM and KSEA interactions. By exploring the maximum extractable work from diverse spin chain models in both ferromagnetic and antiferromagnetic settings, we uncovered significant insights into the mechanisms of energy storage and extraction in these quantum systems.

Our results indicate that the performance of these quantum batteries is notably affected by several factors:

- · Inhomogenous Zeeman splitting and temperature : In antiferromagnetic systems, optimal energy extraction is achieved when the external magnetic field is applied to one of the two spins. In contrast, ferromagnetic systems benefit from applying the magnetic field uniformly to both spins for certain fixed parameter values. Generally, lower temperatures enhance performance across most models. However, the XY model in antiferromagnetic configurations exhibits considerably improved robustness at elevated temperatures. Conversely, in ferromagnetic systems, energy extraction diminishes with increasing temperatures due to heightened thermal fluctuations.
- · Symmetric-antisymmetric exchange interactions : The incorporation of specific spin interactions, such as Dzyaloshinsky-Moriya and Kaplan-Shekhtman-Entin-Wohlman-Aharony couplings, can substantially boost the energy storage and extraction capabilities of these quantum batteries. Nevertheless, there exists a threshold beyond which further increases in interaction strength result in a sharp decline in energy storage capacity.
- · Model-specific behavior : Different Heisenberg spin chain models ( XX , XY , XXZ , XYZ ) exhibit distinct behaviors under varying conditions, highlighting the importance of model selection for optimizing battery performance in specific applications.
- · Quantum phenomena : The interplay between various quantum phenomena, including coherence and entanglement, significantly influences the performance of the battery. The maximum amount of extractable energy is limited by both the battery's inherent capacity and the level of quantum coherence present within the system.

These findings highlight the interplay between energy extraction, quantum battery capacity, and quantum coherence under varying conditions. Our findings reinforce the theoretical framework for finite spin chain-based quantum batteries and provide practical insights into material selection and experimental setups. By recognizing the limitations of quantum resources, this study underscores the need for strategies to optimize their use in advanced quantum systems. Moreover, one can improve the efficiency of energy storage and extraction in quantum batteries by controlling magnetic fields, symmetric and antisymmetric exchange interactions, temperature, and quantum effects.

In summary, this research offers valuable insights into

- [1] R. Dell and D. A. J. Rand, Understanding batteries , Vol. 28 (Royal society of chemistry, 2001).
- [2] X. Luo, J. Wang, M. Dooner, and J. Clarke, Overview of current development in electrical energy storage technologies and the application potential in power system operation, Applied energy 137 , 511 (2015).
- [3] D. Ferraro, M. Campisi, G. M. Andolina, V. Pellegrini, and M. Polini, High-power collective charging of a solidstate quantum battery, Physical Review Letters 120 , 117702 (2018).
- [4] F. Campaioli, F. A. Pollock, and S. Vinjanampathy, Quantum batteries, Thermodynamics in the Quantum Regime: Fundamental Aspects and New Directions , 207 (2018).
- [5] F. C. Binder, S. Vinjanampathy, K. Modi, and J. Goold, Quantacell: powerful charging of quantum batteries, New Journal of Physics 17 , 075015 (2015).
- [6] H.-L. Shi, S. Ding, Q.-K. Wan, X.-H. Wang, and W.-L. Yang, Entanglement, coherence, and extractable work in quantum batteries, Physical Review Letters 129 , 130602 (2022).
- [7] F. Kamin, F. Tabesh, S. Salimi, and A. C. Santos, Entanglement, coherence, and charging process of quantum batteries, Physical Review E 102 , 052109 (2020).
- [8] R. Alicki and M. Fannes, Entanglement boost for extractable work from ensembles of quantum batteries, Physical Review E 87 , 042123 (2013).
- [9] I. Henao and R. M. Serra, Role of quantum coherence in the thermodynamics of energy transfer, Physical Review E 97 , 062105 (2018).
- [10] A. Albrecht, R. Baunach, and A. Arrasmith, Einselection, equilibrium, and cosmology, Physical Review D 106 , 123507 (2022).
- [11] W. H. Zurek, Decoherence, einselection, and the quantum origins of the classical, Reviews of Modern Physics 75 , 715 (2003).
- [12] J. Q. Quach, G. Cerullo, and T. Virgili, Quantum batteries: The future of energy storage?, Joule 7 , 2195 (2023).
- [13] C.-K. Hu, J. Qiu, P. J. P. Souza, J. Yuan, Y. Zhou, L. Zhang, J. Chu, X. Pan, L. Hu, J. Li, Y. Xu, Y. Zhong, S. Liu, F. Yan, D. Tan, R. Bachelard, C. J. Villas-Boas, A. C. Santos, and D. Yu, Optimal charging of a superconducting quantum battery, Quantum Science and Technol-

optimizing quantum batteries derived from spin systems, contributing to the development of efficient energy storage devices in quantum technologies.

## ACKNOWLEDGEMENTS

M.G. and S.H. were supported by Semnan University under Contract No. 21270.

## DISCLOSURES

The authors declare that they have no known competing financial interests.

ogy 7 , 045018 (2022).

- [14] J. Joshi and T. Mahesh, Experimental investigation of a quantum battery using star-topology NMR spin systems, Physical Review A 106 , 042601 (2022).
- [15] J. Q. Quach, T. Virgili, G. Cerullo, K. McGhee, L. Ganzer, and D. Lidzey, Organic quantum batteries, in International Conference on Ultrafast Phenomena (Optica Publishing Group, 2020) pp. Tu4A-1.
- [16] J. Q. Quach, K. E. McGhee, L. Ganzer, D. M. Rouse, B. W. Lovett, E. M. Gauger, J. Keeling, G. Cerullo, D. G. Lidzey, and T. Virgili, Superabsorption in an organic microcavity: Toward a quantum battery, Science Advances 8 , eabk3160 (2022).
- [17] X. Huang, K. Wang, L. Xiao, L. Gao, H. Lin, and P. Xue, Demonstration of the charging progress of quantum batteries, Physical Review A 107 , L030201 (2023).
- [18] T. P. Le, J. Levinsen, K. Modi, M. M. Parish, and F. A. Pollock, Spin-chain model of a many-body quantum battery, Physical Review A 97 , 022106 (2018).
- [19] R. Grimaudo, Y. Belousov, H. Nakazato, and A. Messina, Time evolution of a pair of distinguishable interacting spins subjected to controllable and noisy magnetic fields, Annals of Physics 392 , 242 (2018).
- [20] A. Ali, S. Al-Kuwari, and S. Haddadi, Trade-off relations of quantum resource theory in Heisenberg models, Physica Scripta 99 , 055111 (2024).
- [21] R. Grimaudo, A. Messina, and H. Nakazato, Exactly solvable time-dependent models of two interacting twolevel systems, Physical Review A 94 , 022108 (2016).
- [22] A. Ali, S. Al-Kuwari, M. Rahim, M. Ghominejad, H. Ali, and S. Haddadi, A study on thermal quantum resources and probabilistic teleportation in spin-1/2 Heisenberg XYZ+DM+KSEA model under variable Zeeman splitting, Applied Physics B 130 , 177 (2024).
- [23] R. Grimaudo, N. Vitanov, and A. Messina, Couplingassisted Landau-Majorana-St¨ckelberg-Zener u transition in a system of two interacting spin qubits, Physical Review B 99 , 174416 (2019).
- [24] O. V. Marchukov, A. G. Volosniev, M. Valiente, D. Petrosyan, and N. Zinner, Quantum spin transistor with a Heisenberg spin chain, Nature communications 7 , 13070 (2016).

- [25] S. Murmann, F. Deuretzbacher, G. Z¨rn, J. Bjerlin, S. M. u Reimann, L. Santos, T. Lompe, and S. Jochim, Antiferromagnetic heisenberg spin chain of a few cold atoms in a one-dimensional trap, Physical Review Letters 115 , 215301 (2015).
- [26] G. De Chiara, T. Calarco, S. Fishman, and G. Morigi, Ramsey interferometry with a spin embedded in a coulomb chain, Physical Review A 78 , 043414 (2008).
- [27] J. Simon, W. S. Bakr, R. Ma, M. E. Tai, P. M. Preiss, and M. Greiner, Quantum simulation of antiferromagnetic spin chains in an optical lattice, Nature 472 , 307 (2011).
- [28] S. Bose, Quantum communication through spin chain dynamics: an introductory overview, Contemporary Physics 48 , 13 (2007).
- [29] T. Yildirim, A. B. Harris, A. Aharony, and O. EntinWohlman, Anisotropic spin Hamiltonians due to spinorbit and Coulomb exchange interactions, Physical Review B 52 , 10239 (1995).
- [30] J. Lee and M. Kim, Entanglement teleportation via werner states, Physical Review Letters 84 , 4236 (2000).
- [31] T. E. Lee, S. Gopalakrishnan, and M. D. Lukin, Unconventional magnetism via optical pumping of interacting spin systems, Physical Review Letters 110 , 257204 (2013).
- [32] J. Gaudet, E. Smith, J. Dudemaine, J. Beare, C. Buhariwalla, N. Butch, M. Stone, A. Kolesnikov, G. Xu, D. Yahne, et al. , Quantum spin ice dynamics in the dipole-octupole pyrochlore magnet Ce Zr O , Physical 2 2 7 Review Letters 122 , 187201 (2019).
- [33] A. S. Patri, M. Hosoi, and Y. B. Kim, Distinguishing dipolar and octupolar quantum spin ices using contrasting magnetostriction signatures, Physical Review Research 2 , 023253 (2020).
- [34] A. Signoles, T. Franz, R. F. Alves, M. G¨rttner, S. Whita lock, G. Z¨rn, and M. Weidem¨ller, Glassy dynamics in u u a disordered Heisenberg quantum spin system, Physical Review X 11 , 011011 (2021).
- [35] Y. Yeo, Teleportation with a mixed state of four qubits and the generalized singlet fraction, Physical Review A 74 , 052305 (2006).
- [36] G.-F. Zhang and S.-S. Li, Thermal entanglement in a two-qubit Heisenberg XXZ spin chain under an inhomogeneous magnetic field, Physical Review A 72 , 034302 (2005).
- [37] F. Ozaydin and A. A. Altintas, Quantum metrology: Surpassing the shot-noise limit with Dzyaloshinskii-Moriya interaction, Scientific Reports 5 , 16360 (2015).
- [38] C. J. van Diepen, T.-K. Hsiao, U. Mukhopadhyay, C. Reichl, W. Wegscheider, and L. M. Vandersypen, Quantum simulation of antiferromagnetic Heisenberg chain with gate-defined quantum dots, Physical Review X 11 , 041025 (2021).
- [39] P. Cappellaro, C. Ramanathan, and D. G. Cory, Simulations of information transport in spin chains, Physical Review Letters 99 , 250506 (2007).
- [40] F. Meier, J. Levy, and D. Loss, Quantum computing with spin cluster qubits, Physical Review Letters 90 , 047901 (2003).
- [41] H. Yu, Y. Zhao, T.-C. Wei, et al. , Simulating largesize quantum spin chains on cloud-based superconducting quantum computers, Physical Review Research 5 , 013183 (2023).
- [42] R. T. Chepuri and I. A. Kov´cs, Complex quantum neta work models from spin clusters, Communications Physics 6 , 271 (2023).
- [43] F.-Q. Dou, H. Zhou, and J.-A. Sun, Cavity Heisenbergspin-chain quantum battery, Physical Review A 106 , 032212 (2022).
- [44] Y. Yao and X. Shao, Optimal charging of open spin-chain quantum batteries via homodyne-based feedback control, Physical Review E 106 , 014138 (2022).
- [45] S. Ghosh, T. Chanda, A. Sen, et al. , Enhancement in the performance of a quantum battery by ordered and disordered interactions, Physical Review A 101 , 032115 (2020).
- [46] D. Rossini, G. M. Andolina, and M. Polini, Many-body localized quantum batteries, Physical Review B 100 , 115142 (2019).
- [47] E. Mascarenhas, G. Giudice, and V. Savona, A nonequilibrium quantum phase transition in strongly coupled spin chains, Quantum 1 , 40 (2017).
- [48] Y. Huangfu and J. Jing, High-capacity and high-power collective charging with spin chargers, Physical Review E 104 , 024129 (2021).
- [49] M. B. Arjmandi, H. Mohammadi, A. Saguia, M. S. Sarandy, and A. C. Santos, Localization effects in disordered quantum batteries, Physical Review E 108 , 064106 (2023).
- [50] S.-f. Qi and J. Jing, Magnon-mediated quantum battery under systematic errors, Physical Review A 104 , 032606 (2021).
- [51] S. Zakavati, F. T. Tabesh, and S. Salimi, Bounds on charging power of open quantum batteries, Physical Review E 104 , 054117 (2021).
- [52] M. B. Arjmandi, H. Mohammadi, and A. C. Santos, Enhancing self-discharging process with disordered quantum batteries, Physical Review E 105 , 054115 (2022).
- [53] I. Medina, A. Drinko, G. I. Correr, P. C. Azado, and D. O. Soares-Pinto, Variational-quantum-eigensolverinspired optimization for spin-chain work extraction, Phys. Rev. A 110 , 012443 (2024).
- [54] M. Hadipour, S. Haseli, H. Dolatkhah, and M. Rashidi, Study the charging process of moving quantum batteries inside cavity, Scientific Reports 13 , 10672 (2023).
- [55] X.-L. Zhang, X.-K. Song, and D. Wang, Quantum battery in the Heisenberg spin chain models with Dzyaloshinskii-Moriya interaction, Advanced Quantum Technologies 7 , 2400114 (2024).
- [56] B. Mojaveri, R. Jafarzadeh Bahrbeig, and M. Fasihi, Extracting ergotropy from nonequilibrium steady states of an XXZ spin-chain quantum battery, Physical Review A 109 , 042619 (2024).
- [57] M.-L. Song, L.-J. Li, X.-K. Song, L. Ye, and D. Wang, Environment-mediated entropic uncertainty in charging quantum batteries, Phys. Rev. E 106 , 054107 (2022).
- [58] M.-L. Song, X.-K. Song, L. Ye, and D. Wang, Evaluating extractable work of quantum batteries via entropic uncertainty relations, Phys. Rev. E 109 , 064103 (2024).
- [59] Z.-G. Lu, G. Tian, X.-Y. L¨, and C. Shang, Topological u quantum batteries (2024), arXiv:2405.03675 [quant-ph].
- [60] A. Catalano, S. Giampaolo, O. Morsch, V. Giovannetti, and F. Franchini, Frustrating quantum batteries, PRX Quantum 5 , 030319 (2024).
- [61] B. Ahmadi, P. Mazurek, P. Horodecki, and S. Barzanjeh, Nonreciprocal quantum batteries, Phys. Rev. Lett. 132 , 210402 (2024).

- [62] R. R. Rodr´ ıguez, B. Ahmadi, G. Su´ arez, P. Mazurek, S. Barzanjeh, and P. Horodecki, Optimal quantum control of charging quantum batteries, New Journal of Physics 26 , 043004 (2024).
- [63] J.-Y. Gyhm and U. R. Fischer, Beneficial and detrimental entanglement for quantum battery charging, AVS Quantum Science 6 , 012001 (2024).
- [64] M. Hadipour and S. Haseli, Enhancing the efficiency of open quantum batteries via adjusting the classical driving field, Results in Physics 64 , 107928 (2024).
- [65] M. Hadipour, S. Haseli, D. Wang, and S. Haddadi, Proposed scheme for a cavity-based quantum battery, Advanced Quantum Technologies 7 , 2400115 (2024).
- [66] C. Radhakrishnan, M. Parthasarathy, S. Jambulingam, and T. Byrnes, Quantum coherence of the heisenberg spin models with Dzyaloshinsky-Moriya interactions, Scientific Reports 7 , 13865 (2017).
- [67] C. Radhakrishnan, I. Ermakov, and T. Byrnes, Quantum coherence of planar spin models with DzyaloshinskyMoriya interaction, Physical Review A 96 , 012341 (2017).
- [68] F. Zhao, F.-Q. Dou, and Q. Zhao, Quantum battery of interacting spins with environmental noise, Physical Review A 103 , 033715 (2021).
- [69] S.-Q. Liu, L. Wang, H. Fan, F.-L. Wu, and S.-Y. Liu, Better performance of quantum batteries in different environments compared to closed batteries, Physical Review A 109 , 042411 (2024).
- [70] T. K. Konar, L. G. C. Lakkaraju, and A. Sen, Quantum battery with non-hermitian charging, Physical Review A 109 , 042207 (2024).
- [71] A. Ullah, M. T. Naseem, and ¨ . O E. M¨ ustecaplıo˘lu, g Low-temperature quantum thermometry boosted by coherence generation, Physical Review Research 5 , 043184 (2023).
- [72] L. Shekhtman, O. Entin-Wohlman, and A. Aharony, Moriya's anisotropic superexchange interaction, frustration, and Dzyaloshinsky's weak ferromagnetism, Physical Review Letters 69 , 836 (1992).
- [73] A. Szilva, Y. Kvashnin, E. A. Stepanov, L. Nordstr¨m, o O. Eriksson, A. I. Lichtenstein, and M. I. Katsnelson, Quantitative theory of magnetic interactions in solids,
- Reviews of Modern Physics 95 , 035004 (2023).
- [74] M. Kuepferling, A. Casiraghi, G. Soares, G. Durin, F. Garcia-Sanchez, L. Chen, C. H. Back, C. H. Marrows, S. Tacchi, and G. Carlotti, Measuring interfacial dzyaloshinskii-moriya interaction in ultrathin magnetic films, Reviews of Modern Physics 95 , 015003 (2023).
- [75] Q. Zhang, J. Liang, K. Bi, L. Zhao, H. Bai, Q. Cui, H.-A. Zhou, H. Bai, H. Feng, W. Song, et al. , Quantifying the Dzyaloshinskii-Moriya interaction induced by the bulk magnetic asymmetry, Physical Review Letters 128 , 167202 (2022).
- [76] T. Kaplan, Single-band hubbard model with spin-orbit coupling, Zeitschrift f¨r Physik B Condensed Matter u 49 , 313 (1983).
- [77] M. A. Yurischev and S. Haddadi, Local quantum Fisher information and local quantum uncertainty for general X states, Phys. Lett. A 476 , 128868 (2023).
- [78] M. A. Yurischev, On the quantum correlations in two-qubit XYZ spin chains with Dzyaloshinsky-Moriya and Kaplan-Shekhtman-Entin-Wohlman-Aharony interactions, Quantum Inf. Process. 19 , 336 (2020).
- [79] A. Ullah and S. J. Anwar, Characterization of squeezed thermal baths via qubit probe, The European Physical Journal B 94 , 134 (2021).
- [80] F. Binder, S. Vinjanampathy, K. Modi, and J. Goold, Quantum thermodynamics of general quantum processes, Physical Review E 91 , 032119 (2015).
- [81] W. Pusz and S. L. Woronowicz, Passive states and kms states for general quantum systems, Communications in Mathematical Physics 58 , 273 (1978).
- [82] A. E. Allahverdyan, R. Balian, and T. M. Nieuwenhuizen, Maximal work extraction from finite quantum systems, Europhysics Letters 67 , 565 (2004).
- [83] J. Carrasco, J. R. Maze, C. Hermann-Avigliano, and F. Barra, Collective enhancement in dissipative quantum batteries, Phys. Rev. E 105 , 064119 (2022).
- [84] X. Yang, Y.-H. Yang, M. Alimuddin, R. Salvia, S.-M. Fei, L.-M. Zhao, S. Nimmrichter, and M.-X. Luo, Battery capacity of energy-storing quantum systems, Physical Review Letters 131 , 030402 (2023).
- [85] T. Baumgratz, M. Cramer, and M. B. Plenio, Quantifying coherence, Physical Review Letters 113 , 140401 (2014).
Asad Ali Saif Al-Kuwari M. I. Hussain Tim Byrnes M. T. Rahim James Q. Quach Mehrdad Ghominejad Saeed Haddadi
Ergotropy and capacity optimization in Heisenberg spin-chain quantum batteries
This study examines the performance of finite spin quantum batteries (QBs)
using Heisenberg spin models with Dzyaloshinsky-Moriya (DM) and
Kaplan--Shekhtman--Entin-Wohlman--Aharony (KSEA) interactions. The QBs are
modeled as interacting quantum spins in local inhomogeneous magnetic fields,
inducing variable Zeeman splitting. We derive analytical expressions for the
maximal extractable work, ergotropy and the capacity of QBs, as recently
examined by Yang et al. [Phys. Rev. Lett. 131, 030402 (2023)]. These quantities
are analytically linked through certain quantum correlations, as posited in the
aforementioned study. Different Heisenberg spin chain models exhibit distinct
behaviors under varying conditions, emphasizing the importance of model
selection for optimizing QB performance. In antiferromagnetic (AFM) systems,
maximum ergotropy occurs with a Zeeman splitting field applied to either spin,
while ferromagnetic (FM) systems benefit from a uniform Zeeman field.
Temperature significantly impacts QB performance, with ergotropy in the AFM
case being generally more robust against temperature increases compared to the
FM case. Incorporating DM and KSEA couplings can significantly enhance the
capacity and ergotropy extraction of QBs. However, there exists a threshold
beyond which additional increases in these interactions cause a sharp decline
in capacity and ergotropy. This behavior is influenced by temperature and
quantum coherence, which signal the occurrence of a sudden phase transition.
The resource theory of quantum coherence proposed by Baumgratz et al. [Phys.
Rev. Lett. 113, 140401 (2014)] plays a crucial role in enhancing ergotropy and
capacity. However, ergotropy is limited by both the system's capacity and the
amount of coherence. These findings support the theoretical framework of
spin-based QBs and may benefit future research on quantum energy storage
devices.