Geometric quantum discord and coherence in a dipolar interacting magnetic system

Page created by Leroy Patterson
 
CONTINUE READING
Geometric quantum discord and coherence in a dipolar interacting magnetic system

                                                              Clebson Cruz,1, ∗ Maron F. Anka,2, † Hamid-Reza Rastegar-Sedehi,3, ‡ and Cleidson Castro4, §
                                                                1
                                                                    Grupo de Informação Quântica e Fı́sica Estatı́stica, Centro de Ciências Exatas e das Tecnologias,
                                                                       Universidade Federal do Oeste da Bahia - Campus Reitor Edgard Santos. Rua Bertioga,
                                                                                       892, Morada Nobre I, 47810-059 Barreiras, Bahia, Brasil.
                                                                                         2
                                                                                           Instituto de Fı́sica, Universidade Federal Fluminense,
                                                                            Av. Gal. Milton Tavares de Souza s/n, 24210-346 Niterói, Rio de Janeiro, Brasil.
                                                                       3
                                                                         Department of Physics, College of Sciences, Jahrom University, Jahrom 74135-111, Iran
                                                                          4
                                                                            Centro de Formação de Professores, Universidade Federal do Recôncavo da Bahia,
                                                                                      Avenida Nestor de Mello Pita, 535 Amargosa, Bahia, Brazil.
                                                                                                          (Dated: January 10, 2023)
                                                               The study of low-dimensional metal complexes has revealed fascinating characteristics regarding the ground-
                                                            state crossover shown by spin-gaped systems. In this context, this work explores the effect of the quantum-level
arXiv:2301.02891v1 [quant-ph] 7 Jan 2023

                                                            crossing, induced by the magnetic anisotropies of dipolar interacting systems, on the quantum discord and
                                                            coherence of the system. The analytical expressions for the quantum discord, based on Schatten 1-norm, and
                                                            the l1 trace-norm quantum coherence for dinuclear spin-1/2 systems, are provided in terms of the magnetic
                                                            anisotropies. The results show that, while the quantum discord has a clear signature of the quantum level-
                                                            crossing, the basis dependence of the quantum coherence hides the crossover regarding the measured basis. In
                                                            addition, the global quantum coherence is wholly stored within the correlations of the system, regardless of its
                                                            reference basis.
                                                            Keywords: Dipolar Interaction; Qunatum discord; Quantum Coherence; Quantum-level crossing.

                                                                                                         I.   INTRODUCTION

                                              The study of the quantum properties of composite systems has led to a revolution in the development of emerging quantum
                                           technologies [1–3]. The new generation of quantum devices explores physical properties associated with quantum correlations
                                           between particles [4, 5] and superposition principle for the system states [1, 6, 7]. In this scenario, the characterization of the
                                           quantumness of the physical systems is of paramount importance since the existence of quantum correlations and coherence are
                                           a valuable resource for several quantum tasks [4, 8, 9].
                                              However, the characterization of quantum correlations is a rather complicated task from the theoretical [10] and experimental
                                           [11] point of view. This scenario is aggravated in Condensed Matter systems, where the number of interacting components in
                                           the system is usually on the order of the Avogadro number [12]. Nevertheless, there are a few exceptions, like low-dimensional
                                           metal complexes (LDMC), for which full knowledge about their quantum properties can be obtained through the corresponding
                                           analytical solutions [4, 6, 13–20]. In such solid-state systems, intra-molecular interactions are strong enough to suppress ex-
                                           trinsic and intermolecular interactions [4, 13, 14, 21]. Therefore, their quantum features exhibit high stability against external
                                           perturbations such as temperatures [4, 13, 14, 16, 21], magnetic fields [4, 6, 16, 22], and pressures [6, 15]. These characteris-
                                           tics make these systems promising platforms for the development of emerging quantum technologies [4, 23–26]. In regard to
                                           these possible applications, the study of dipolar interacting magnetic systems has received considerable attention in the quantum
                                           information literature [17, 27–31].
                                              In this work, we present a theoretical study of quantum correlations and coherence for a dipolar interacting magnetic system,
                                           exploring the effects of magnetic anisotropies on the quantumness of the system. As a result, this study provides to the literature
                                           analytical expressions, in terms of the magnetic anisotropies, for the quantum discord, based on Schatten 1-norm, and the l1
                                           trace-norm quantum coherence written in an arbitrary basis, defined by the co-latitude and longitude angles of the Bloch sphere
                                           representation. According to the findings, the behavior of the quantum discord carries a noteworthy signature of the quantum
                                           level-crossing, caused by population changes resulting from the alteration of Boltzmann weights arising from the change of the
                                           magnetic anisotropies of the system. On the other hand, the basis dependency of quantum coherence is detrimental in terms of
                                           recognizing this crossover. In this regard, the measurement of the average quantum coherence is numerically obtained in order
                                           to obtain a basis-independent perspective for this quantum resource. The results not only demonstrate that the average coherence

                                           ∗ Electronic address: clebson.cruz@ufob.edu.br
                                           † Electronic address: maronanka@id.uff.br
                                           ‡ Electronic address: h.rastegar@jahromu.ac.ir
                                           § Electronic address: ccastro@ufrb.edu.br
2

is completely stored within the correlations of the system, yet they also demonstrate that it is possible to retrieve the signature of
the energy-level crossover present on the quantum discord measurement. Furthermore, the findings show how dipolar interaction
coupling magnetic anisotropies impact quantum correlations and coherence in a dinuclear spin-1/2 system. Thus, the dipolar
interaction model is a viable foundation for quantum technologies based on quantum discord and coherence.

                            II.   DINUCLEAR METAL COMPLEX WITH DIPOLAR INTERACTION

  The class of dinuclear metal complexes undergoes several types of magnetic coupling [12]. Among these are Heisenberg
exchange, which is isotropic under rotations in spin space [4, 6, 20, 21], and Dzyaloshinskii-Moriya (DM) interaction [32–
34], which accounts for weak ferromagnetism in some antiferromagnetic materials [12]. A ubiquitous example of anisotropic
coupling in LDMCs is the dipolar interaction [17, 27–31]. This coupling arises from the influence of a magnetic field yielded
by one of the magnetic moments in the other ones [12]. In particular, for dinuclear metal complexes, the Hamiltonian which
describes this interaction is given by:
                                                               1        ↔
                                                          H = − S~ TA · D · S~ B ,                                                  (1)
                                                               3
                                                       ↔
where S~ j = {S xj , S yj , S zj } are the spin operators and D = diag(∆ − 3, ∆ + 3, −2∆) is a diagonal tensor, with  and ∆ being the
rhombic and axial parameters, respectively, related to the magnetic anisotropies in the dipolar model [12]. In particular, ∆ is
related to the zero-field splitting of the energy levels [35]. Considering the Hamiltonian, Eq. (1), written in the S (z) eigenbasis,
∆ > 0 becomes a signature that the spins are on z-axis, while ∆ < 0 indicates that the spins will be on the x − y plane.
   Considering a dinuclear metal complex in with d9 electronic configuration, Eq. (1) can describe two coupled spin 1/2 particles
in the corresponding S z eigenbasis {|00i, |01i, |10i, |11i}

                                                                      ∆
                                                                                     
                                                                                  3 
                                                                 1       −∆ −∆
                                                                                      
                                                            H =                        ,                                      (2)
                                                                                        
                                                                 6       −∆ −∆
                                                                                  ∆
                                                                                            
                                                                         3
The energy levels of the system from the coupling parameters are composed by the
                                      E1 = 0 ,    E2 = −2∆ ,          E3 = ∆ + 3 ,                  E4 = ∆ − 3 .                  (3)
  From the thermal equilibrium, the density matrix for the coupled system is described by the Gibbs form ρAB = Z−1 e−H/kB T ,
where
                                                                 β∆                    β 
                                                                     !
                                 Z = Tr(e−H/kB T ) = 2eβ∆/6 cosh       + 2e−β∆/6 cosh        .                          (4)
                                                                  6                     2
is the canonical the partition function, with kB representing the Boltzmann’s constant. Thus, the dinuclear density matrix at sites
labeled by A and B can be written in the S z eigenbasis as the so-called X-shaped mixed state

                                                cosh β2                                       − sinh β2 
                                                                                                       

                                                                   e 6 cosh β∆               β∆
                                        β∆                         2β∆            2β∆
                                                                                                              
                                      e− 6                                     6 sinh
                                                 
                                                                              6 e
                               ρAB =                                                                                   .
                                                                                                                    
                                                                                              6 
                                                                  2β∆     β∆    2β∆      β∆                                 (5)
                                       Z       
                                                                   e 6 sinh 6 e 6 cosh 6                                
                                                           β                                        β  
                                                     − sinh 2                                      cosh 2

The density matrix eigenvalues (population) and their corresponding eigenvectors can be written as:
                                                                      1
                                             PΨ− =          β∆            β∆         β  → |Ψ− i,                                 (6)
                                                     1+e     3   + 2e− 6 cosh           2
                                                                      1
                                             PΨ+ =                                       β  → |Ψ+ i,                             (7)
                                                            − β∆           β∆
                                                     1+e       3   + 2e− 6 cosh             2
                                                                      1
                                             PΦ+ =                                      ∆+3
                                                                                                    → |Φ+ i,                        (8)
                                                                    β( ∆+
                                                     1+    eβ   +e     2 )     + eβ(    6      )
                                                                 eβ
                                             PΦ− =                ∆+       ∆+3
                                                                                 → |Φ− i,                                           (9)
                                                     1 + eβ + eβ( 2 ) + eβ( 6 )
3

where
                                                           1                                   1
                                                   |Ψ± i = √ (|01i ± |10i) ,           |Φ± i = √ (|00i ± |11i) .                                    (10)
                                                            2                                   2
are the so-called Bell states, which represent the maximally entangled states for a bipartite system [36].
   The study of LDMC has attracted the attention of both theoretical and experimental condensed matter physics communities
due to the fascinating properties of their ground states [6, 37]. In the presence of an external magnetic field, this systems
typically show a quantum level-crossing between its ground state and the first excited one when the field reaches a critical value
since the external magnetic field splits its energy levels, changing their corresponding populations. However, since the dipolar
interaction arises from the influence of the magnetic field created by one of the magnetic moments in the other, the splitting in
energy levels is ruled by the axial (∆) and rhombic () parameters, as can be seen in Eq. (3). In this regard, Fig. 1 shows the
populations as a function of the ratio between the magnetic anisotropies and the energy scale factor kB T .
                        1.0                                                                   1.0
                                                                                                                                              PΨ-
                                                            ϵ/ kB= 5K                                                            ϵ/ kB=-5K    PΨ+
                        0.8                                                                   0.8
                                                                                                                                              PΦ+
          Populations

                                                                                Populations
                                                                                                                                              PΦ-
                        0.6                                                                   0.6

                        0.4                                                                   0.4

                        0.2                                                                   0.2

                        0.0                                                                   0.0
                              - 20   - 10      0       10        20                                 - 20   - 10      0      10        20

                                            Δ/kBT                                                                 Δ /kBT
                                                                          (a)

                        1.0                                                                   1.0

                        0.8                                                                   0.8
         Populations

                                                                                Populations

                        0.6                                                                   0.6

                        0.4                                                                   0.4

                        0.2                                                                   0.2
                                                            Δ /kB = 5K                                                           Δ /kB =-5K
                        0.0                                                                   0.0
                              - 20   - 10      0       10        20                                 - 20   - 10      0      10        20

                                            ϵ /kB T                                                               ϵ /kB T

                                                                          (b)

FIG. 1: (Color online) Populations, Eqs. (6)–(9), as a function of the ratio between the magnetic anisotropies and the energy scale factor kB T .
(a) Axial dependence considering the rhombic parameter /kB = 5 K (left) and /kB = −5 K (right). (b) Rhombic dependence considering the
axial parameter ∆/kB = 5 K (left) and ∆/kB = −5 K (right). The inset shows the magnetic anisotropy dependence on the energy levels.

   As can be seen, in agreement with Eq. (3), when the spins are in the x − y plane (∆ < 0) with positive rhombic parameter
( > 0), the ground state is the given by the state |Φ− i (with population PΦ− ), and there is no quantum level crossing. Thus, the
system remains in the ground state. However, changing the signal of the rhombic parameter induces a quantum level crossing
between states |Φ− i and |Φ+ i (with population PΦ+ ). Moreover, for the spins oriented in the z axis (∆ > 0), it is possible to
observe a quantum level crossing between the state |Φ− i (if  > 0) or |Φ+ i (if  < 0) and the state |Ψ− i (with population PΨ+ ) by
increasing the ratio ∆/kB T to the critical point ∆ = ||.
   In reference [29], the authors study the effect of the magnetic anisotropies, described by the axial and rhombic parameters,
on the nonlocal correlations of a dipolar interacting system of two spins-1/2, identified by the Peres-Horodecki separability
criterion [38, 39]. In addition, they explore the change in the ground state on the thermal entanglement for the teleportation
process. However, although quantum entanglement provides one path toward the characterization of nonlocal correlations, it
does not encompass all quantum correlations in the system [21, 40–51]. Therefore, in order to expand this result, the following
section presents a study of the quantum correlations and coherence described by the Schatten 1-norm geometric quantum discord
and the l1 trace-norm quantum coherence.
4

                                                     III.   QUANTUM DISCORD

   Quantum discord has been defined as a measurement of the quantumness of correlations in a quantum system. It has been first
introduced as an entropic measurement of genuinely quantum correlations in a quantum state, defined as the difference between
the total and the classical correlation [40] Q(ρAB ) = I(ρA : ρB ) − C(ρAB ), where I(ρA : ρB ) = S (ρA ) + S (ρB ) − S (ρAB ) represents
the mutual information between  the subsystems      A and B, and C(ρAB ) is the classical correlation of the composite system ρAB
defined as C(ρAB ) = max{Bk } S (ρA ) − k pk S (ρk ) , with the maximization taking over positive operator-valued measurements
                                         P
(POVM’s) {Bk } performed locally only on subsystem B. However, this analytical maximization over POVMs is an arduous task
even for a two-qubit system [10, 11, 40, 41, 47, 51]. In this scenario, the class of entropic measurements of correlations, such as
the entropic quantum discord, is defined as nondeterministic polynomial time (NP-complete) problems [10]. Consequently, only
a few results for the analytical expression of entropic quantum discord, and only for certain classes of states are exact solutions
known [41, 43, 45, 46, 49, 50, 52]. Due to this fact, alternative measurements of quantum correlations have been proposed
[41, 43, 45–50, 52–60], especially quantifiers based on geometric arguments [47–49, 53, 57–60].
   Geometric approaches are widely used to characterize and quantify quantum resources in a wide variety of quantum systems
[61]. In particular, the Schatten 1-norm quantum discord [4, 21, 47, 48, 62], is a reliable geometric-based quantifier of the
amount of quantum correlations in metal complexes [4, 21, 62, 63]. The so-called geometric quantum discord can be defined in
terms of the minimal distance between a set ω of closest classical-quantum states ρc [21, 47, 48], given by:
                                                              X
                                                         ρc =    pk Π{A}
                                                                     k ⊗ ρk ,
                                                                            {B}
                                                                                                                                     (11)
                                                                k

where 0 ≤ pk ≤ 1 and k pk = 1; {Πk{A} } define a set of orthogonal projectors for a given subsystem A and ρ{B}
                         P
                                                                                                           k the reduced
density matrix for the subsystem B [47, 48]. Therefore, the geometric quantum discord can be expressed as

                                                    QG (ρAB ) = min kρAB − ρc k ,                                                   (12)
                                                                     ω
                h√      i
where kMk = Tr M † M is the so-called 1-norm, and ρAB is the given quantum state at thermal equilibrium, Eq. (5).
  Therefore, considering the given dinuclear magnetic system of spins-1/2 in a quantum spin-lattice, ruled by a dipolar Hamil-
tonian H, Eq.(1), the invariance under π rotation around a given spin axis (Z2 symmetry) [12, 46] allow us to compute the
geometric quantum discord, based on Schatten 1-norm, for the two-qubit X state, Eq. (5), as [53, 64]
                                                 s
                                               1     φ21 max{φ22 , φ23 } − φ22 min{φ21 , φ23 }
                                   QG (ρAB ) =                                                                            (13)
                                               2 max{φ22 , φ23 } − min{φ21 , φ23 } + φ21 − φ22

where
                                                     β∆
                                                   e 6 −1 + eβ∆/3 + 2 sinh β2
                                                                              
                                              φ1 =                                 ,                                                (14)
                                                      2 cosh β2 + eβ∆/6 + eβ∆/2
                                                             

                                                   e 6 −1 + eβ∆/3 − 2 sinh β2
                                                     ∆
                                                                             
                                              φ2 =                               ,                                                  (15)
                                                     2 cosh β2 + eβ∆/6 + eβ∆/2
                                                             

                                                                         2
                                              φ3 =                    β∆            β         −1.                              (16)
                                                      eβ∆/3   cosh       6    sech     2      +1

   Considering the dipolar magnetic system in thermal equilibrium described by Eq. (5), it is possible to examine how the
magnetic anisotropies, represented by the axial (∆) and rhombic () coupling parameters, affects the thermal quantum discord
in the system. Fig. 2 shows the geometric quantum discord, based on Schatten 1-norm, Eq. (13), as a function of the ratio
∆/kB T and /kB T . As expected, the quantum discord reaches its maximum (saturated) value of 1/2 as T approaches zero. As
the temperature rises, the value of quantum discord decreases inexorably and goes to zero when T  |∆| and T  ||. On the
other hand, given the spins in the x − y plane (∆ < 0), it is sufficient that only T  || to the discord reaches its minimum value.
However, if the spins are in the z-axis (∆ > 0), one can increase the quantum discord by increasing the axial parameter ∆ even
when T  ||.
   Furthermore, regarding the magnetic anisotropies, the quantum discord presents a signature of the quantum level crossing in
the dipolar interacting system, highlighted on the solid white line in Fig. 2. Considering the spins oriented in the z-axis (∆ > 0),
the zero-field splitting leads the system to a quantum level crossing in the critical boundary ∆ = ||, where it is possible to
detect a crossover between the states |Ψ+ i and |Φ− i, if  > 0, or |Φ+ i, if  < 0. Moreover, for the spins oriented in the x − y
5

plane (∆ < 0), it is possible to observe a quantum level crossing between the state |Φ− i and |Φ− i in the critical boundary  = 0.
On the other hand, the degree of quantum discord in the system can be increased by gradient ascent of the function QG (ρAB ),
perpendicularly to the crossing boundary, which occurs for values in which ||  kB T (for ∆ < 0), ∆  kB T (for ||  ∆), and
||  ∆), corresponding to the lightest region in Fig. 2.Therefore, by controlling the axial (∆) and rhombic () anisotropies is
possible to manage the degree of quantum discord in the dipolar interacting system.
   In addition, in order to compare quantum discord to the level of entanglement in the system under investigation, we use the
concurrence measure. Typically, concurrence is used to assess entanglement in bipartite systems, and it can be easily computed
for any two-qubit system. The thermal concurrence examines the resemblance between the considered quantum state in thermal
equilibrium and its bit-flipped density matrix, ρ̄ = ρAB (σy ⊗ σy )ρ∗AB (σy ⊗ σy ). In particular, for the X- shaped density matrix, Eq.
(5), the concurrence is analytically defined as

                                                         C(ρAB ) := max{0, A, B},                                                    (17)

where
                                                         β∆
                                                  e− 6               β∆            β #
                                                              "         !
                                                                  2β∆
                                              A =        e sinh    6       − cosh         ,                                          (18)
                                                   Z                  6             2
                                                    β∆ "       β 
                                                  e− 6                             β∆
                                                                                       !#
                                                                         2β∆
                                              B =        sinh        − e 6 cosh           .                                          (19)
                                                   Z            2                   6

   Dashed green line in Fig. 2 denotes the boundary given by C(ρAB ) = 0. Inside this region, the concurrence is zero, and the
state of the system is separable. However, within the region where entanglement is absent, the quantum discord of the system
is still considerably more than zero, ensuring the presence of quantum-correlated states even when the system is in a separable
syaye. On the other hand, for low temperatures, the entanglement is zero in the quantum level crossing boundary alongside
the quantum discord at the quantum level crossing boundary. In this scenario, the existence or absence of entanglement and,
therefore, quantum correlations, is dependent on its ground state, which might vary in response to magnetic anisotropies. Thus,
the variation of Boltzmann’s weights, Eqs. (6)-(9), associated with the occupancy of the energy levels, is the physical mechanism
responsible for the abrupt change in the quantum correlations near the energy-level crossover.

                                                              PΦ-

                                          ϵ
                                                                                    PΨ+
                                       kB T

                                                              PΦ+

                                                                         Δ
                                                                        kB T

FIG. 2: (Color online) Quantum Discord, based on Schatten 1-norm, for a dipolar interacting magnetic system, Eq. (13), as a function of
the ratios ∆/kB T and /kB T . The solid white line denotes the boundary between the quantum level crossings. The dashed green line is the
boundary given by the concurrence, Eq. (17), C(ρAB ) = 0, inside which the entanglement of the system is absent.

                                                   IV.        QUANTUM COHERENCE

   Similar to the approach proposed for the entanglement theory, where the quantum entanglement can be characterized by the
distance between a state of interest (ρ) and a set of states closed under local operations, and classical communication (separable
states) [38, 61, 65], Baumgratz et al. [65] provided the mathematical tools for quantifying the amount of quantum coherence
in a quantum system. Considering a d-dimensional Hilbert space, quantum coherence can be obtained from the minimal value
6

of a distance measurement D(ρ, σ), between the considered quantum state ρ and a set {σ = dk |kihk| ∈ I} of incoherent states,
                                                                                                P
where the reference basis {|ki}{k=1,...,d} can be adequately defined considering the physics of the problem under investigation or
the task that requires this quantum resource [6, 65, 66]. In this scenario, since the non-vanishing off-diagonal terms of the
density operator ρ, which characterizes the quantum state of the system of interest, constitute the superposition from the chosen
reference basis [61, 65], the authors established a reliable measurement of quantum coherence through the l1 trace norm as [65]
                                                                         X
                                                  Cl1 = min kρ − σkl1 =      |hi|ρ| ji| .                                     (20)
                                                                σ∈I
                                                                                   i, j

   Since coherence is a quantity that is reliant on the basis on which it is measured, it is essential to choose a reference basis for
the system within a metrology setting [61, 65]. In this scenario, the basis of an arbitrary quantum state can be altered by means
of unitary operations [36, 61]. In particular, for two-level systems such as spin-1/2, any reference basis can be obtained from the
unitary transformation

                                                            cos θ    −eiφ sin 2θ 
                                                                                
                                               U(θ, φ) =  −iφ   2             ,                                          (21)
                                                              e sin 2θ     cos 2θ
                                                                     

where the θ and φ angles are the spherical equivalents of the co-latitude with respect to the z-axis, and the longitude concerning
the x-axis in a Bloch sphere representation, respectively [36, 67]. In this regard, the unitary transformation for the bipartite state
given by Eq. (5) is given by ρ{θ,φ}
                              AB = ÛAB (θ, φ)ρAB ÛAB (θ, φ) [67], where

                                                       cos2 2θ         −eiφ sin 2θ cos 2θ −eiφ sin 2θ cos 2θ      e2iφ sin2 2θ
                                                                                                                           
                                                                                                                                     
                                                                                                                                   
                                                            θ        θ               θ                  θ                θ        θ
                                             −iφ                                                             
                                                                                 2                   2           iφ
                                               e sin         cos 2         cos 2             − sin 2         −e sin 2 cos 2 
    ÛAB (θ, φ) = U(θ, φ) ⊗ U(θ, φ) =  −iφ           2                                                                       (22)
                                                e sin 2θ cos 2θ          − sin2 2θ           cos2 2θ        −eiφ sin 2θ cos 2θ 
                                                                                                                  

                                                     e−2iφ sin2 2θ     e−iφ sin 2θ cos 2θ e−iφ sin 2θ cos 2θ         cos2 2θ
                                                                                                                          

  By varying the co-latitude and longitude angles {θ, φ}, one can obtain the bipartite state ρAB , Eq. (5), in any reference basis.
Using the unitary transformation for the bipartite states, Eq. (22), in Eq. (5), one can obtain the representation of the density
operator for the dipolar interacting magnetic system of two spins-1/2 written in an arbitrary basis as

                                                             %11 %12 %12 %14 
                                                                                   
                                                       β∆ 
                                                    e− 6  %∗12 %22 %23 −%12 
                                            ρ{θ,φ}
                                              AB =
                                                                                     ,                                      (23)
                                                     4Z  %12 %23 %22 −%12 
                                                              ∗
                                                               %∗14 −%∗12 −%∗12 %11

where
                                     β∆            β∆               β                         β 
                                        !               !                           !
           %11 = 2 sin2 (θ) sinh          + cosh          − sinh          cos(2φ) + cosh               (cos(2θ) + 3) ,                        (24)
                                      2             2                2                            2
                                                    β                            β                           β               
           %12 =    e−3iφ sin(θ) 2e2iφ cos(θ) cosh           − eβ∆/2 + e4iφ sinh             (cos(θ) + 1) + sinh             (cos(θ) − 1) ,   (25)
                                                       2                                2                               2
                                          β                              β 
                                                                                      4 θ
                                                                                                           β 
                                                                                                                    4 θ
                                                                                                                  
                                                  β∆/2
           %14 =    2e −2iφ    2
                            sin (θ) cosh        −e        − 4e −4iφ
                                                                     sinh         sin         − 4 sinh           cos       ,                  (26)
                                           2                                 2           2                  2          2
                                                               β                       β 
           %22 =    2 eβ∆/2 cos2 (θ) + eβ∆/6 + sin2 (θ) sinh            cos(2φ) + cosh                ,                                       (27)
                                                                  2                           2
                                                                  β                       β 
           %23 =    2eβ∆/2 cos2 (θ) − 2eβ∆/6 + 2 sin2 (θ) sinh           cos(2φ) + cosh                 .                                     (28)
                                                                     2                          2
   The diagonal entries of Eq. (23) are real, and the trace is 1. In addition, to ensure real eigenvalues, hermiticity restricts
off-diagonal elements to two complex numbers, i.e., %i j is the complex conjugate of % ji .
   Thus, from Eqs. (20) and (23), it is possible to write an analytical expression for the normalized quantum coherence in an
arbitrary basis, defined by the co-latitude and longitude angles {θ, φ}, as:
                                                                   β∆
                                                               e− 6 
                                                    C{θ,φ}   =        4 |%12 | + |%14 | + |%23 | .
                                                                                                
                                                     l1                                                                                       (29)
                                                               6 |Z|
   In order to examine the relationship between quantum coherence and quantum correlations, a new metric known as corre-
lated coherence was established recently [67–69]. Quantum correlated coherence is a measure of coherence in which all local
7

components have been eliminated, i.e., all coherence in the system is totally recorded in the quantum correlations. For any
given quantum state ρ, the correlated contribution to quantum coherence may be calculated by subtracting the local coherence
of subsystems ρA = TrB (ρ) and ρB = TrA (ρ) from the overall coherence [68, 69]. Thus, the definition of correlated coherence
according to the l1 -norm of coherence is:

                                              Ccc (ρ{θ,φ}
                                                    AB ) := Cl1 (ρ) − Cl1 (ρA ) − Cl1 (ρ B ).                                     (30)

   Considering the density matrix of the dipolar interacting magnetic system written in an arbitrary basis, Eq. (23), the reduced
density matrices of local subsystems are ρA = ρB = I/2, the maximally mixed state. Thus, regardless of the basis, the local
subsystems will remain in the maximally mixed state, since it is basis invariant [36]. Consequently, the local contribution for
the quantum coherence in this dipolar interacting system is always null, and the global coherence of the system, Eq. (29),
is totally recorded in the quantum correlations of the system, regardless of its reference basis. Therefore, for a number of
different combinations of values for the co-latitude and longitude angles {θ, φ}, the unitary transformation, Eq. (22), gives a
direct connection between the overall and the correlated degrees of coherence.

                                                          A.    Axial Coherence

   In particular, due to the rotation symmetry of the dipolar interaction, the density matrix will be invariant when rotated both
spins by an angle π along any given spin axis. Thus, choosing the co-latitude angle as θ = nπ (n = {0, 1, 2, ...}), regardless of
the longitude angle φ, one can obtain the density matrix in the X-shaped form as described in Eq. (5). On the other hand, by
applying the unitary transformation for the bipartite states, Eq. (22), for {θ = π/2; φ = nπ}, and {θ = π/2; φ = nπ/2}, in Eq. (5)
one can obtain the density matrix S (x) and S (y) eigenbasis, respectively.

                                         eβ∆/2 + e∓β/2                                       ∓ eβ∆/2 − e∓β/2 
                                                                                                               
                                                                        0              0
                                   β∆ 
                                e− 6                         eβ∆/6 + e±β/2 e±β/2 − eβ∆/6
                                                                                                                    
                                                      0                                                  0
                       ρAB =
                         {X,Y}
                                                                                                                         . (31)
                                                                                                                      
                                                                           β∆/6   β∆/6
                                                                  ±β/2
                                                                                       +   ±β/2
                                          
                                 2Z               0          e      − e     e        e               0
                                                 β∆/2                                               β∆/2
                                                               
                                                                                                         +e
                                                         ∓β/2                                             ∓β/2
                                                                                                                           
                                              ∓ e     −e
                                            
                                                                        0              0          e

   As can be seen, due to the symmetry of the X-shaped density matrices [70], the X-structure of the operator is preserved.
Therefore, from Eqs. (5), (20) and (31), one can obtain the analytical expressions for the normalized axial quantum coherences
as
                                                                  β∆                  β  #
                                                     "               !
                                                 2      β∆/6
                                      Cl1 =
                                        {Z}
                                                      e      sinh      +e −β∆/6
                                                                                sinh                                      (32)
                                               3 |Z|               6                   2
                                                    β∆
                                                  e− 6 h β∆/2
                                                              − e∓β/2 + eβ∆/6 − e±β/2
                                                                                        i
                                      Cl{X,Y} =         e                                                                         (33)
                                         1        3 |Z|
    Fig. 3 shows the axial quantum coherence in S (i) spin eigenbasis, where i = {x, y, z}. Different from the behavior observed
for quantum discord (see Fig. 2), the axial quantum coherence is not sensible to the quantum level crossing. The quantum
coherence in each axis {x, y, z} is minimized in only one energy-level crossover. As can be seen, considering the spins oriented
in the z-axis (∆ > 0), the axial coherence in the S (x) eigenbasis is minimized on the critical boundary ∆ = − (with  < 0), where
it is possible to detect a crossover between the states |Ψ+ i and |Φ+ i, while the coherence in S (y) eigenbasis is minimized on the
critical boundary ∆ =  (with  > 0), where it is possible to detect a crossover between the states |Ψ+ i and |Φ− i (see Fig. 2). As
expected from Eq. 33, if the rhombic parameter is null ( = 0), C{X}  l1 = Cl1 . On the other hand, for the spins oriented in the x − y
                                                                            {Y}

plane, (∆ < 0), it is possible to observe that C{Z}                                                                     −
                                                 l1 is minimized in the quantum level crossing between the state |Φ i and |Φ i in
                                                                                                                                  −

the critical boundary  = 0.
    As shown in Fig. 3, the basis dependence of the quantum coherence hides the energy-level crossover in this dipolar interacting
system regarding the measured basis. Therefore, the basis dependence on the quantum coherence defined by Baumgratz et al.
[65], can be unfavorable to recognizing the quantum level crossing caused by population changes resulting from the alteration
of Boltzman weights, Eqs. (6)-(9), arising from the change of the magnetic anisotropies of the dipolar interacting system.

                                                         B.    Average Coherence

  Since the coherence formulated in the quantum resource theory is a basis-dependent measurement [8, 61, 66], it is natural to
define a basis-independent measurement [71–75]. Recent research has shown, via the use of relative entropies, as distance mea-
surements of quantum correlations, that basis-independent measurements of entropic quantum coherence are precisely identical
8

        10                                             10                                                  10

                 SX eigenbasis                                  SY eigenbasis                                       SZ eigenbasis
                                                                                                                                                 0.30
         5                                              5                                                   5
                                                                                                                                                 0.25

  ϵ                                              ϵ                                                   ϵ                                           0.20
         0                                              0                                                   0
kB T                                           kB T                                                kB T                                          0.15

                                                                                                                                                 0.10
        -5                                             -5                                                  -5
                                                                                                                                                 0.05

       - 10                                           - 10                                                - 10
          - 10      -5     0      5      10              - 10        -5         0      5      10             - 10      -5     0     5   10
                           Δ                                                   Δ                                             Δ
                          kB T                                                kB T                                          kB T
                          (a)                                                  (b)                                           (c)

FIG. 3: (Color online) Axial quantum coherence based on l1 trace norm, for a dipolar interacting magnetic system, as a function of the ratios
∆/kB T and /kB T . The dashed white line represents the minimum value for the axial quantum coherence.

to entropic discord [74]. On the other hand, a possible basis-free measurement of quantum coherence for a quantum system
can be obtained from a geometrical standpoint by averaging the coherence of a state across all reference bases [71–73, 75].
From a theoretical point of view, this measurement corresponds to averaging the coherence on a standard basis across all equiv-
alent states ρ{θ,φ}
              AB = ÛAB (θ, φ)ρAB ÛAB (θ, φ). Therefore, as any two-qubit reference base can be created by applying the unitary
operation described in Eq. (22), the average quantum coherence can be obtained from Eq. (29) as

                                                                     Z2π Zπ
                                                             1
                                                   hCl1 i =                   sin (θ)C{θ,φ}
                                                                                      l1 dθdφ .                                         (34)
                                                            4π
                                                                     0    0

   It is worth mentioning that these integrals are not trivial to solve, and an analytical expression for the average coherence is not
presented. However, it can be numerically integrated by any quadrature method [76]. In this scenario, Eq. (34) is estimated by
using the Clenshaw-Curtis rule on adaptively refined subintervals of the integration area [76, 77] since the numerical integration
algorithms are often equally efficient and effective as conventional algorithms for well-behaved integrands such as Eqs. (29) and
(30) [76].
   Fig. 4 shows the average quantum coherence for the dipolar magnetic interacting system. The solid white line represents the
threshold at which the quantum-level crossing, described in previous sections, actually occurs. As expected, based on Fig. 3,
when the temperature rises reaching the threshold T  |∆| and T  ||, the value of coherence reaches its lowest point and
will be equal to zero. However, the behavior of the average coherence is completely different from that observed in the axial
(basis-dependent) coherence shown in Fig. 3.
   Moreover, besides unified frameworks from relative entropic measurements has shown that basis-independent entropic quan-
tum coherence is equivalent to entropic discord [74], this is not true for this geometrical approach. However, although the
contour lines of the average coherence are quite different from that shown in the discord presented in Fig. 2, it is still able to
identify the signature of the energy-level crossing that was seen during the measurement of the quantum discord. This result is
due to the fact that the global coherence is totally stored within the correlations of the system, and its average behavior is affected
by the presence of genuine quantum correlations measured by the quantum discord.
   In addition, the entanglement of the system is absent within the area shown by the dashed green line that denotes the boundary
supplied by the concurrence, which is denoted by Eq. (17), C(ρAB ) = 0. Thus, as one would anticipate based on the observation
of the quantum discord in Fig. 2, even in the absence of entanglement, the average coherence that is completely stored on the
correlations of the system is noticeably distinct from zero.

                                                                V.   CONCLUSIONS

   In summary, this paper explored the influence of magnetic anisotropies on the quantumness of a dipolar interacting magnetic
system via a theoretical examination of the geometric quantum discord, measured by Schatten 1-norm, and the l1 trace-norm
quantum coherence. The analytical formulations for these quantum information quantifiers were obtained in terms of magnetic
anisotropies. In this scenario, the effects of dipolar coupling constants on these quantifiers are highlighted. It is demonstrated
that the presence of dipolar anisotropies increases the degree to which the system possesses quantum correlation and coherence.
9

                                                 10

                                                             PΦ-                                   0.7
                                                  5                                                0.6

                                                                                                   0.5
                                          ϵ
                                                  0                                    PΨ+         0.4
                                        kB T
                                                                                                   0.3

                                                -5                                                 0.2

                                                             PΦ+                                   0.1

                                               - 10
                                                  - 10       -5         0          5         10
                                                                       Δ
                                                                      kB T
FIG. 4: (Color online) Average quantum coherence based on l1 trace norm, for a dipolar interacting magnetic system, as a function of the ratios
∆/kB T and /kB T . The solid white line represents the boundary between the quantum level crossings. The dashed green line is the boundary
given by the concurrence, Eq. (17), C(ρAB ) = 0, inside which the entanglement of the system is absent.

   As another remarkable result, it is proved that the global coherence, expressed in an arbitrary reference basis, determined by
the co-latitude and longitude angles of the Bloch sphere representation, is totally stored within the correlations of the system.
Moreover, according to the results, the behavior of quantum discord contains a notable hallmark of quantum level-crossing in
the system, in contrast to the basis-dependent axial quantum coherence, which hides the energy-level crossover regarding the
measured basis.
   Therefore, the dependency of the base on the quantum coherence specified by Baumgratz might be deleterious in identifying
the crossing of levels owing to population changes originating from the changing of Boltzman weights due to the modification
of the magnetic anisotropies of the studied system. In this regard, the average quantum coherence was measured numerically
obtained in order to gain a viewpoint independent of the reference basis, unraveling that the average coherence is able to extract
the signature of the energy-level crossover present in the measurement of quantum discord.
   Finally, the findings that were given provide light on the ways in which magnetic anisotropies caused by the dipolar interaction
coupling of a dinuclear spin-1/2 system influence quantum correlations and coherence. Therefore, the dipolar interaction model
is an excellent option for usage as a platform for quantum technologies that are based on quantum resources such as quantum
coherence and quantum discord.

                                                         ACKNOWLEDGEMENTS

  C. Cruz gratefully acknowledges Mario Reis for the valuable discussions. M. F. Anka thanks FAPERJ for financial support.

 [1]   M. Mohseni, P. Read, H. Neven, S. Boixo, V. Denchev, R. Babbush, A. Fowler, V. Smelyanskiy, and J. Martinis, Nature 543, 171 (2017).
 [2]   M. Atzori and R. Sessoli, Journal of the American Chemical Society 141, 11339 (2019).
 [3]   I. H. Deutsch, PRX Quantum 1, 020101 (2020).
 [4]   C. Cruz, M. F. Anka, M. S. Reis, R. Bachelard, and A. C. Santos, Quantum Science and Technology 7, 025020 (2022).
 [5]   G. L. Giorgi and S. Campbell, J. Phys. B: At. Mol. Opt. Phys. 48, 035501 (2015).
 [6]   C. Cruz and M. Anka, EPL (Europhysics Letters) 130, 30006 (2020).
 [7]   F. Caravelli, G. Coulter-De Wit, L. P. Garcı́a-Pintos, and A. Hamma, Physical Review Research 2, 023095 (2020).
 [8]   A. Streltsov, G. Adesso, and M. B. Plenio, Reviews of Modern Physics 89, 041003 (2017).
 [9]   F. Sapienza, F. Cerisola, and A. J. Roncaglia, Nature communications 10, 1 (2019).
[10]   Y. Huang, New journal of physics 16, 033027 (2014).
[11]   M. Cramer, M. Plenio, and H. Wunderlich, Physical review letters 106, 020401 (2011).
[12]   M. Reis (Academic Press, Boston, 2013), ISBN 978-0-12-405545-2.
[13]   A. M. Souza, D. O. Soares-Pinto, R. S. Sarthour, I. S. Oliveira, M. S. Reis, P. Brandao, and A. M. dos Santos, Physical Review B 79,
       054408 (2009).
10

[14]   M. S. Reis, S. Soriano, A. M. dos Santos, B. C. Sales, D. Soares-Pinto, and P. Brandao, EPL (Europhysics Letters) 100, 50001 (2012).
[15]   C. Cruz, Á. Alves, R. dos Santos, D. Soares-Pinto, J. de Jesus, J. de Almeida, and M. Reis, EPL (Europhysics Letters) 117, 20004 (2017).
[16]   H. Čenčariková and J. Strečka, Physical Review B 102, 184419 (2020).
[17]   E. I. Kuznetsova and M. A. Yurischev, Quantum Information Processing 12, 3587 (2013).
[18]   M. A. Yurishchev, Physical Review B 84, 024418 (2011).
[19]   S. Aldoshin, E. Fel’dman, and M. Yurishchev, Low Temperature Physics 40, 3 (2014).
[20]   C. Cruz, H.-R. Rastegar-Sedehi, M. F. Anka, T. R. de Oliveira, and M. Reis, arXiv preprint arXiv:2208.14548 (2022).
[21]   C. Cruz, D. O. Soares-Pinto, P. Brand?o, A. M. dos Santos, and M. S. Reis, EPL (Europhysics Letters) 113, 40004 (2016).
[22]   A. M. Souza, M. S. Reis, D. O. Soares-Pinto, I. S. Oliveira, and R. S. Sarthour, Physical Review B 77, 104402 (2008).
[23]   M. R. Wasielewski, M. D. Forbes, N. L. Frank, K. Kowalski, G. D. Scholes, J. Yuen-Zhou, M. A. Baldo, D. E. Freedman, R. H. Goldsmith,
       T. Goodson, et al., Nature Reviews Chemistry pp. 1–15 (2020).
[24]   A. Gaita-Ariño, F. Luis, S. Hill, and E. Coronado, Nature chemistry 11, 301 (2019).
[25]   Y. A. Mezenov, A. A. Krasilin, V. P. Dzyuba, A. Nominé, and V. A. Milichko, Advanced Science 6, 1900506 (2019).
[26]   E. Moreno-Pineda, C. Godfrin, F. Balestro, W. Wernsdorfer, and M. Ruben, Chemical Society Reviews 47, 501 (2018).
[27]   D. F. Pinto and J. Maziero, Quantum Information Processing 17, 1 (2018).
[28]   D. F. Pinto and J. Maziero, Quantum Information Processing 20, 1 (2021).
[29]   C. Castro, O. Duarte, D. Pires, D. Soares-Pinto, and M. Reis, Physics Letters A 380, 1571 (2016).
[30]   A. A. Mohamed, H. Hessian, and H. Eleuch, Physica Scripta 95, 075104 (2020).
[31]   R. Muthuganesan and V. Chandrasekar, Physica Scripta 96, 125113 (2021).
[32]   R. Hoshikawa, K. Yoshida, R. Mitsuhashi, M. Mikuriya, T. Okuno, and H. Sakiyama, Molecules 26, 897 (2021).
[33]   M.-A. Bouammali, N. Suaud, R. Maurice, and N. Guihéry, The Journal of Chemical Physics 155, 164305 (2021).
[34]   M.-A. Bouammali, N. Suaud, C. Martins, R. Maurice, and N. Guihéry, The Journal of Chemical Physics 154, 134301 (2021).
[35]   E. Moreno-Pineda and W. Wernsdorfer, Nature Reviews Physics 3, 645 (2021).
[36]   M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information: 10th Anniversary Edition (Cambridge University
       Press, New York, NY, USA, 2011), 10th ed., ISBN 1107002176, 9781107002173.
[37]   T. Chakraborty and C. Mitra, Journal of Physics: Condensed Matter 31, 475802 (2019).
[38]   R. Horodecki, P. Horodecki, M. Horodecki, and K. Horodecki, Reviews of modern physics 81, 865 (2009).
[39]   A. Peres, Phys. Rev. Lett. 77, 1413 (1996).
[40]   H. Ollivier and W. H. Zurek, Physical Review Letters 88, 017901 (2001).
[41]   V. Vedral, Physical Review Letters 90, 050401 (2003).
[42]   B.-Q. Liu, L.-A. Wu, G.-M. Zeng, J.-M. Song, W. Luo, Y. Lei, G.-A. Sun, B. Chen, and S.-M. Peng, Physics Letters A 378, 3441 (2014).
[43]   Z. Ma, Z. Chen, F. F. Fanchini, and S.-M. Fei, Scientific Reports 5 (2015).
[44]   D. Girolami and G. Adesso, Physical Review A 83, 052108 (2011).
[45]   T. Nakano, M. Piani, and G. Adesso, Physical Review A 88, 012117 (2013).
[46]   M. Sarandy, Physical Review A 80, 022108 (2009).
[47]   F. Paula, T. R. de Oliveira, and M. Sarandy, Physical Review A 87, 064101 (2013).
[48]   J. Montealegre, F. Paula, A. Saguia, and M. Sarandy, Physical Review A 87, 042115 (2013).
[49]   S. Luo, Physical Review A 77, 042303 (2008).
[50]   A. Datta, A. Shaji, and C. M. Caves, Physical Review Letters 100, 050502 (2008).
[51]   L. Henderson and V. Vedral, Journal of Physics A: Mathematical and General 34, 6899 (2001).
[52]   A. Brodutch and D. R. Terno, Physical Review A 81, 062103 (2010).
[53]   P. C. Obando, F. M. Paula, and M. S. Sarandy, Physical Review A 92, 032307 (2015).
[54]   D. Girolami, T. Tufarelli, and G. Adesso, Physical review letters 110, 240402 (2013).
[55]   D. Girolami, A. M. Souza, V. Giovannetti, T. Tufarelli, J. G. Filgueiras, R. S. Sarthour, D. O. Soares-Pinto, I. S. Oliveira, and G. Adesso,
       Physical Review Letters 112, 210401 (2014).
[56]   D. Girolami, A. M. Souza, V. Giovannetti, T. Tufarelli, J. G. Filgueiras, R. S. Sarthour, D. O. Soares-Pinto, I. S. Oliveira, and G. Adesso,
       Physical Review Letters 112, 210401 (2014).
[57]   B. Dakić, V. Vedral, and Č. Brukner, Physical Review Letters 105, 190502 (2010).
[58]   M. Piani, Physical Review A 86, 034101 (2012).
[59]   F. Paula, A. Saguia, T. R. de Oliveira, and M. Sarandy, EPL (Europhysics Letters) 108, 10003 (2014).
[60]   D. Spehner, F. Illuminati, M. Orszag, and W. Roga, arXiv preprint arXiv:1611.03449 (2016).
[61]   M.-L. Hu, X. Hu, J. Wang, Y. Peng, Y.-R. Zhang, and H. Fan, Physics Reports (2018).
[62]   Y. Khedif, S. Haddadi, M. Daoud, H. Dolatkhah, and M. R. Pourkarimi, Quantum Information Processing 21, 1 (2022).
[63]   C. Cruz, arXiv preprint arXiv:1610.05255 (2016).
[64]   F. Ciccarello, T. Tufarelli, and V. Giovannetti, New Journal of Physics 16, 013038 (2014).
[65]   T. Baumgratz, M. Cramer, and M. Plenio, Physical review letters 113, 140401 (2014).
[66]   A. Streltsov, G. Adesso, and M. B. Plenio, Reviews of Modern Physics 89, 041003 (2017).
[67]   C. Filgueiras, O. Rojas, and M. Rojas, Annalen der Physik 532, 2000207 (2020).
[68]   T. Kraft and M. Piani, Journal of Physics A: Mathematical and Theoretical 51, 414013 (2018).
[69]   K. C. Tan, H. Kwon, C.-Y. Park, and H. Jeong, Physical Review A 94, 022329 (2016).
[70]   A. Rau, Journal of Physics A: Mathematical and Theoretical 42, 412002 (2009).
[71]   X.-Y. Liu and M.-L. Hu, Physica A: Statistical Mechanics and its Applications 609, 128308 (2023).
[72]   S. Luo and Y. Sun, Physics Letters A 383, 2869 (2019).
[73]   S. Cheng and M. J. Hall, Physical Review A 92, 042101 (2015).
11

[74] Y. Yao, X. Xiao, L. Ge, and C. Sun, Physical Review A 92, 022112 (2015).
[75] S. Designolle, R. Uola, K. Luoma, and N. Brunner, Physical Review Letters 126, 220404 (2021).
[76] W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, Numerical recipes 3rd edition: The art of scientific computing
     (Cambridge university press, 2007).
[77] G. Liu and S. Xiang, Applied Mathematics and Computation 340, 251 (2019).
You can also read