Home Nonclassical correlation dynamics of Heisenberg XYZ states with (x, y)-spin--orbit interaction, x-magnetic field, and intrinsic decoherence effects
Article Open Access

Nonclassical correlation dynamics of Heisenberg XYZ states with (x, y)-spin--orbit interaction, x-magnetic field, and intrinsic decoherence effects

  • Fahad Aljuaydi EMAIL logo , Sattam N. Almutairi , Heba Allhibi and Abdel-Baset A. Mohamed
Published/Copyright: September 16, 2025

Abstract

In quantum information, it is important to recognize the effects of additional interactions, such as spin–orbit interactions, on the quantum information resources of two-qubit Heisenberg states. Therefore, we study the nonlocal correlation dynamics affected by the intrinsic decoherence of the spin–spin-Heisenberg-XYZ interaction, which is supported by spin–orbit interactions (Dzyaloshinsky–Moriya) of the x and y directions together. The two spins are coupled to an external inhomogeneous magnetic field (EIMF) in the x -direction. We investigate and compare the nonclassical correlation dynamics of local quantum Fisher information, local quantum uncertainty, and log-negativity. The results show that spin–spin and spin–orbit interactions have a high capability to enhance non-local correlations in the presence of an external magnetic field. The enhanced non-local correlation can be further improved by strengthening the spin–spin and spin–orbit interactions, as well as by increasing the EIMF’s inhomogeneity and uniformity, which increases the amplitudes and fluctuations of the generated non-local correlation oscillations. The degradation of non-local correlations due to intrinsic decoherence can be controlled by spin–spin interactions. These degradation correlations can be enhanced by increasing the intensities of spin–orbit interactions, as well as by increasing the EIMF’s inhomogeneity and uniformity.

1 Introduction

Among the various quantum systems proposed for implementing quantum information and computation [1,2], superconducting circuits, trapped ions, and semiconductor quantum dots are crucial techniques for realizing quantum bits (qubits). Based on electron spins trapped in quantum dots, a quantum computer protocol has been initially proposed [35]. The electron, having a spin of ( 1 2 ), is the simplest natural qubit. Recently, quantum computation with electron spins (as a single-spin–qubit geometric gate) has been realized in quantum dots [6,7]. Due to electron tunneling from one dot to another, spin–spin coupling and spin–orbit coupling interactions between two qubits can be realized by considering a two-qubit system represented by two electrons in coupled quantum dots. Therefore, Heisenberg XYZ models describing spin–spin interactions are among the important proposed qubit systems. Two-qubit Heisenberg XYZ models have been realized in various systems, including bosonic atoms inside an optical lattice [8], trapped ions [9], superconductor systems [10], and linear molecules [11]. These models have been updated to include the first order of spin–orbit coupling known as Dzyaloshinsky–Moriya interactions [1214], realized through an antisymmetric superexchange interaction in La2CuO4 [15], and the second order of spin–orbit coupling known as the Kaplan–Shekhtman–Entin–Wohlman–Aharony interaction [16]. Additionally, Heisenberg XYZ models have been updated to include dipole–dipole interactions [17] and inhomogeneous external magnetic fields [18,19]. Recently, the spin–orbit interaction was experimentally realized in two magnetic cobalt layers [20] (with Co/Ag/Co system). Moreover, the spin–orbit interaction was experimentally implemented in the prototypical ferromagnet by polarized neutrons [21] and in ferroelectrics and antiferroelectrics [22] systems as well as in the epitaxial Ni/Cu ( 001 ) system [23]. The Heisenberg XYZ qubit models have shown several potential applications in teleportation, [24], quantum dense coding [25], thermodynamics [26], and quantum correlation generations [27].

Exploring two-qubit information dynamics in various proposed qubit systems, based on different types of nonlocal correlations (NLCs) (such as entanglement and quantum discord), is one of the most critical research fields in implementing quantum information and computation [28]. Quantum entanglement (QE), quantified by measures such as entropy [29], concurrence [30], negativity, and log-negativity [31], is a significant type of two-qubit NLCs [32,33]. It has a wide range of applications in quantum information fields, including quantum computation, teleportation [34,35], quantum optical memory [36], and quantum key distribution [37]. After implementing quantum discord as another type of qubits’ NLCs beyond entanglement [38], several quantifiers have been introduced to address other NLCs [39] using Wigner–Yanase (WY) skew information [40] and quantum Fisher information (QFI) [41]. WY-skew-information minimization (local quantum uncertainty, LQU) [42] and WY-skew-information maximization (uncertainty-induced nonlocality) [43] have been introduced to quantify other NLCs beyond entanglement. Additionally, the minimization of QFI (local quantum Fisher information, LQFI) has been used to implement another two-qubit NLC [44,45]. LQU has a direct connection to LQFI [46,47], establishing more two-qubit NLCs in several proposed qubit systems [48], such as hybrid–spin systems (under random noise [49] and intrinsic decoherence [50]), two-coupled double quantum dots [51], the mixed-spin Heisenberg model [52], and the Heisenberg system [53].

The information dynamics of two-spin Heisenberg XYZ states have been investigated using the Milburn intrinsic decoherence model [54]. This includes studies on entanglement teleportation based on the Heisenberg XYZ chain [24,55], the LQFI of Heisenberg XXX states beyond IEMF effects [56], and quantum correlations of concurrence and LQU [57]. Previous works have focused on exploring the time evolution of the two-spin Heisenberg XYZ states’ NLCs under limited conditions on spin–spin and spin–orbit interactions, as well as applied magnetic fields, to ensure residing quantum information resources of two-qubit X -states (having an X density matrix) [5862].

Motivated by the aforementioned experimental evidence for realizations for the spin–orbit interaction having a high ability to support the generating NLCs, and the importance of general two-qubit Heisenberg states, this study employs the Milburn intrinsic decoherence and Heisenberg XYZ models to explore the NLC dynamics of LQFI, LQU, and log-negativity for general two-qubit Heisenberg XYZ states with non- X density matrices, influenced by specific conditions on spin–spin and spin–orbit interactions, as well as applied magnetic fields.

The manuscript structure includes the Milburn intrinsic decoherence equation, the Heisenberg XYZ model, and its solution in Section 2. In Section 3, we introduce the definitions of the NLCs’ quantifiers: LQFI, LQU, and LN. Section 4 presents the outcomes of the dependence of these quantifiers on the physical parameters. Our conclusions are provided in Section 5.

2 Heisenberg spin model

Here, the Milburn intrinsic decoherence and Heisenberg XYZ models are used to examine the capabilities embedded in spin–spin interaction supported by the spin–orbit (Dzyaloshinsky–Moriya) interactions in the x and y directions (described by the first-order of spin–orbit couplings D x and D x ), to generate essential NLCs between the two spin qubits under the effects of the uniformity and the inhomogeneity of an applied external inhomogeneous magnetic field (EIMF). For two spins (each described by the upper 1 k and lower 1 k states, where ( k = A , B )), the Hamiltonian of the Heisenberg XYZ model with a spin–orbit interaction and an applied EIMF (Figure 1) is written as

(1) H ˆ = α = x , y , z J α σ ˆ A α σ ˆ B α + k = A , B B k σ k + D A B ( σ A × σ B ) ,

where σ k = ( σ ˆ k x , σ ˆ k y , σ ˆ k z ) represents the vector of Pauli matrices of the k -spin, and B k = ( B k x , B k y , B k z ) represents the vector of the external magnetic field applying on k -spin. In our work, we consider that the EIMF is applied only in the x -direction: B k = ( B k x , 0,0 ) , B A x = B m + b m , and B B x = B m b m . Here, B m and b m represent the degree of the uniformity and the inhomogeneity of the applied EIMF, respectively. For the spin–orbit interaction vector D A B = ( D x , D y , D z ) , we have D A B . ( σ A × σ B ) = D x C ˆ x + D y C ˆ y + D z C ˆ z with C ˆ α = σ ˆ A α + 1 σ ˆ B α + 2 σ ˆ A α + 2 σ ˆ B α + 1 ( α = x , y , z ) . After considering only the spin–orbit interactions of the x and y directions D A B = ( D x , D y , 0 ) with an applied EIMF in the x direction, to support the spin–spin interaction in generating two-spin–qubit correlations, the considered Hamiltonian is written as

(2) H ˆ = i = x , y , z J i σ ˆ A i σ ˆ B i + D x ( σ ˆ A y σ ˆ B z σ ˆ A z σ ˆ B y ) + D y ( σ ˆ A z σ ˆ B x σ ˆ A x σ ˆ B z ) + ( B m + b m ) σ ˆ A x + ( B m b m ) σ ˆ B x .

In the two-spin–qubits basis: { ψ 1 = 1 A 1 B , ψ 2 = 1 A 0 B , ψ 3 = 0 A 1 B , ψ 4 = 0 A 0 B } , the two-spin system’s Hamiltonian in (1) can read as a non- X matrix of

(3) H ˆ = J z β β + * J x J y β * J z J x + J y β + β + J x + J y J z β * J x J y β + * β J z

with β ± = B m ± b m + D y + i D x . Generally, with a motion equation, this non- X Hamiltonian matrix generates two-qubit non- X states due to spin–spin interaction combined with x and y spin–orbit interactions. However, if we consider only the spin–orbit interaction in the z direction, D = ( 0 , 0 , D z ) , as discussed in previous studies [5862], the two-spin system’s Hamiltonian (1) generates X -states characterized by a density X -matrix. Because deriving the analytical expressions for the eigenvalues and eigenvectors of the non- X Hamiltonian matrix (Eq. (3)) is very difficult, the eigenvalues are computed numerically.

Figure 1 
               Diagram of a Heisenberg XYZ chain model, where two arbitrary spin–qubits (
                     
                        
                        
                           A
                        
                        A
                     
                   and 
                     
                        
                        
                           B
                        
                        B
                     
                  ) are selected with spin–orbit interaction vector 
                     
                        
                        
                           
                              
                                 
                                    
                                       D
                                    
                                    
                                       →
                                    
                                 
                              
                              
                                 A
                                 B
                              
                           
                           =
                           
                              (
                              
                                 
                                    
                                       D
                                    
                                    
                                       x
                                    
                                 
                                 ,
                                 
                                    
                                       D
                                    
                                    
                                       y
                                    
                                 
                                 ,
                                 
                                    
                                       D
                                    
                                    
                                       z
                                    
                                 
                              
                              )
                           
                        
                        {\overrightarrow{D}}_{AB}=\left({D}_{x},{D}_{y},{D}_{z})
                     
                  , and an EIMFs 
                     
                        
                        
                           
                              
                                 
                                    
                                       B
                                    
                                    
                                       →
                                    
                                 
                              
                              
                                 k
                              
                           
                        
                        {\overrightarrow{B}}_{k}
                     
                   in the 
                     
                        
                        
                           x
                        
                        x
                     
                  -direction.
Figure 1

Diagram of a Heisenberg XYZ chain model, where two arbitrary spin–qubits ( A and B ) are selected with spin–orbit interaction vector D A B = ( D x , D y , D z ) , and an EIMFs B k in the x -direction.

The time evolution of the NLCs in the generated two spin–qubit states, represented by the density matrix M ˆ ( t ) , will be explored using the Milburn intrinsic decoherence model [54], which is given by

(4) d d t M ˆ ( t ) = i [ H ˆ , M ˆ ] γ 2 [ H ˆ , [ H ˆ , M ˆ ] ] ,

where γ is the intrinsic spin–spin decoherence (ISSD) coupling.

After calculating the eigenvalues V k ( k = 1, 2, 3, 4) and the eigenstates V k of the Hamiltonian in Eq. (3), the two-spin density matrix M ˆ ( t ) of Eq. (4) can be obtained numerically using the following solution, given by

(5) M ˆ ( t ) = m , n = 1 4 U m n ( t ) S m n ( t ) V m M ˆ ( 0 ) V n V m V n .

This solution depends on the unitary interaction U m n ( t ) and the ISSD coupling S m n ( t ) , given by the following terms:

(6) U m n ( t ) = e i ( V m V n ) t , S m n ( t ) = e γ 2 ( V m V n ) 2 t .

Eq. (5) is used to numerically calculate and explore the dynamics of the NLCs within the two-spin–qubit states’ Heisenberg XYZ model under the effects of spin–orbit interactions along the x and y directions and an applied external magnetic field in the x direction.

3 NLC quantifiers

Here, the two-spin NLCs will be measured by the following quantifiers: LQFI, LQU, and LN.

  • LQFI

    Here, we use LQFI to quantify another type of two-spin Heisenberg-XYZ correlation beyond entanglement. After calculating the two-spin eigenvalues π k ( k = 1, 2, 3, 4) and the eigenstates Π k of the density matrix of Eq. (5), which has the representation matrix M ( t ) = m π m Π m Π m with π m 0 and m π m = 1 , the LQFI is calculated using the closed expression given by [41,44,45]

    F ( t ) = 1 π R max ,

    which depends on the highest eigenvalue π R max of the symmetric matrix R = [ r i j ] . Based on the Pauli spin- 1 2 matrices σ i ( i = 1, 2, 3) and the elements ξ m n i = Π m I σ i Π n , the symmetric matrix elements r i j are given by

    r i j = π m + π n 0 2 π m π n π m + π n ξ m n i ( ξ n m j ) .

    For a maximally correlated two-spin–qubit state, the LQFI function converges to F ( t ) = 1 . Otherwise, the LQFI function oscillates and is bounded by the inequality 0 < F ( t ) < 1 , indicating that the states have partial LQFI NLC.

  • LQU

    Also, we use LQU of WY skew information [40] to realize another type of two-spin–qubits’ NLC [40,42,43]. For the two-spin density matrix M ( t ) of Eq. (5), the LQU can be calculated by the following closed expression [42]:

    (7) U ( t ) = 1 λ max ( Λ A B ) ,

    which depends on the largest eigenvalue λ max of the 3 × 3 -matrix Λ = [ a i j ] , which have the following elements:

    a i j = T r { M ( t ) ( σ i I ) M ( t ) ( σ j I ) } .

    The LQU function oscillates and is bounded by the inequality 0 U ( t ) 1 . It converges to U ( t ) = 1 , otherwise indicating that the states have partial LQU NLC.

  • Logarithmic negativity (LN)

    We employ LN [31] to measure the generated two-spin–qubit entanglement. The LN expression is based on the negativity’s definition μ t , which is defined as the absolute sum of the negative eigenvalues of the partial transposition matrix ( M ( t ) ) T of the two-spin–qubit density matrix M ( t ) of Eq. (5). The LN can be expressed as

    (8) N ( t ) = log 2 [ 1 + 2 μ t ] .

    The LN function vanishes, N ( t ) = 0 , for a disentangled two-spin state. It converges to its maximum value, N ( t ) = 1 , for a maximally entangled two-spin state. Otherwise, LN oscillates and is bounded by the inequality 0 N ( t ) 1 , indicating that the two-spin states have partial entanglement.

In the following, we work in a system of units where = 1 , and employ the nondimensionalized parameter method as described in previous studies [24,63,64]. We also consider the case of spin–spin interactions with antiferromagnetic couplings satisfying J α > 0 . Meanwhile, the other physical parameters, including the spin–orbit couplings and the degree of uniformity and inhomogeneity of the magnetic field, satisfy D x , D y , B m , b m 0 . Small values of these parameters indicate weak spin–orbit interaction and a weak applied magnetic field.

4 Two-spin qubit dynamics

To explore the generation of non-local correlations between two spin qubits, we consider that the two spins are initially in their uncorrelated upper states 1 A 1 B . In this state, the density matrix has no non-local correlations according to the considered quantifiers. Our focus is on the effects of J α spin–spin interactions ( D x and D y ) and inhomogeneous x -direction magnetic field parameters ( B m and b m ) in the presence of ISSD coupling.

Our first analysis, starting from Figure 2, illustrates the dynamics of non-local correlations (LQFI, LQU, and LN) between two spin qubits. These correlations are generated by the couplings ( J x , J y , J z ) = (0.8, 0.8, 0.8), supported by varying intensities of x and y spin–orbit interactions. This is done in the presence of an inhomogeneous x -direction magnetic field with small uniformity and inhomogeneity ( B m , b m ) = (0.3, 0.5), and in the absence of intrinsic spin–spin decoherence ( γ = 0 ). Figure 2(a) with ( D x , D y ) = (0.0, 0.0) shows that the LQFI, LQU, and log-negativity grow and reach their maximum values. They are subject to slow quasi-regular oscillations with the same frequencies and different amplitudes. LQFI and LQU have the same behavior, i.e., the two spin qubit correlation is called “LQFI–LQU correlation.” The amplitude of the LN is always larger than that of the LQFI and LQU. Under these circumstances of a weak coupling regime ( J α = 0.8 ) and the applied inhomogeneous x -direction magnetic field (with weak uniformity and inhomogeneity), the initial pure uncorrelated two-spin state evolves into various time-dependent partially correlated states. At specific times, it transforms into maximally correlated states. The two-spin states exhibit maximal LQFI–LQU correlation ( F ( t ) = U ( t ) = 1 ) and log-negativity ( N ( t ) = 1 ) simultaneously. At particular times, we observe that partially entangled two-spin states have neither LQFI nor LQU correlation.

Figure 2 
               Time evolution of the LQFI, LQU, and LN are shown with the two-spin couplings 
                     
                        
                        
                           
                              (
                              
                                 
                                    
                                       J
                                    
                                    
                                       x
                                    
                                 
                                 ,
                                 
                                    
                                       J
                                    
                                    
                                       y
                                    
                                 
                                 ,
                                 
                                    
                                       J
                                    
                                    
                                       z
                                    
                                 
                              
                              )
                           
                        
                        \left({J}_{x},{J}_{y},{J}_{z})
                     
                   = (0.8, 0.8, 0.8) and the applied magnetic field parameters 
                     
                        
                        
                           
                              (
                              
                                 
                                    
                                       B
                                    
                                    
                                       m
                                    
                                 
                                 ,
                                 
                                    
                                       b
                                    
                                    
                                       m
                                    
                                 
                              
                              )
                           
                        
                        \left({B}_{m},{b}_{m})
                     
                   = (0.3, 0.5) for different 
                     
                        
                        
                           x
                           ,
                           y
                        
                        x,y
                     
                   spin–orbit interactions: 
                     
                        
                        
                           
                              (
                              
                                 
                                    
                                       D
                                    
                                    
                                       x
                                    
                                 
                                 ,
                                 
                                    
                                       D
                                    
                                    
                                       y
                                    
                                 
                              
                              )
                           
                        
                        \left({D}_{x},{D}_{y})
                     
                   = (0.0, 0.0) in (a), 
                     
                        
                        
                           
                              (
                              
                                 
                                    
                                       D
                                    
                                    
                                       x
                                    
                                 
                                 ,
                                 
                                    
                                       D
                                    
                                    
                                       y
                                    
                                 
                              
                              )
                           
                        
                        \left({D}_{x},{D}_{y})
                     
                   = (0.5, 0.0) in (b), and 
                     
                        
                        
                           
                              (
                              
                                 
                                    
                                       D
                                    
                                    
                                       x
                                    
                                 
                                 ,
                                 
                                    
                                       D
                                    
                                    
                                       y
                                    
                                 
                              
                              )
                           
                        
                        \left({D}_{x},{D}_{y})
                     
                   = (0.5, 0.5) in (c).
Figure 2

Time evolution of the LQFI, LQU, and LN are shown with the two-spin couplings ( J x , J y , J z ) = (0.8, 0.8, 0.8) and the applied magnetic field parameters ( B m , b m ) = (0.3, 0.5) for different x , y spin–orbit interactions: ( D x , D y ) = (0.0, 0.0) in (a), ( D x , D y ) = (0.5, 0.0) in (b), and ( D x , D y ) = (0.5, 0.5) in (c).

The D x -spin–orbit interaction ( D x , D y ) = (0.5, 0) dramatically improves the appearance of the intervals of the maximal LQFI–LQU correlation and log-negativity entanglement, as well as the intervals in which two-spin entangled states have no LQFI or LQU correlation. The effects of weak spin–orbit interactions in the x direction only are shown in Figure 2(b). As illustrated, the regularity and fluctuations of the generated LQFI–LQU correlation and log-negativity entanglement are significantly greater than previously observed in the absence of x and y spin–orbit interactions. The weak D x spin–orbit interaction dramatically enhances the intervals of maximal LQFI-LQU correlation and log-negativity entanglement, as well as the intervals where two-spin entangled states have neither LQFI nor LQU correlation. In Figure 2(c), we combined the D x and D y spin–orbit interactions ( D x , D y ) = (0.5, 0.5). As shown, the fluctuations of the two-spin NLCs between their partial and maximal values are significantly fewer than in Figure 2(a) and (b). Additionally, the NLC frequencies have been reduced, and their lower bounds have shifted upward. This indicates that the combined D x and D y spin–orbit interactions enhance the generated partial two-spin–qubit LQFI–LQU correlation and log-negativity entanglement.

Figure 3(a) and (b) illustrates that higher spin–spin interaction couplings ( ( J x , J y , J z ) = (1, 0.5, 1.5) in (a) and ( J x , J y , J z ) = (5, 1, 1.5) in (b)) significantly enhance the two-spin LQFI–LQU correlation and log-negativity entanglement. By comparing the generated spin–spin NLCs shown in Figur 2(c) and 3(a), we find that relatively strong couplings of J α -spin–spin interactions ( ( J x , J y , J z ) = (1, 0.5, 1.5)), supported by weak D x , y -spin–orbit interactions ( ( D x , D y ) = (0.5, 0.5)), increase the amplitudes and frequencies of the LQFI–LQU correlation and log-negativity entanglement oscillations. Figure 3(a) and (b) shows that higher J α -couplings lead to that the spin–spin NLCs’ oscillations have more fluctuations. The time positions of the maximal LQFI–LQU correlation and log-negativity entanglement are enhanced. Figure 3(c) is plotted to demonstrate the capability of spin–spin interactions ( J α = 0.8 ), supported by x , y -spin–orbit interactions ( D x = D y = 2 ), to enhance the generated spin–spin NLCs when an external magnetic field with weak determinants ( ( B m , b m ) = (0.3, 0.5)) is applied. By comparing the qualitative dynamics of the generated LQFI-LQU correlation and log-negativity entanglement shown in Figure 2(c) ( D x = D y = 0.5 ) with those in Figure 3(c) ( D x = D y = 2 ), we can deduce that D x , y -spin–orbit interactions play a significant role in enhancing the generated LQFI-LQU correlation and log-negativity entanglement. Their amplitudes are increased, and their oscillations exhibit more fluctuations between extreme values. Additionally, strong x , y -spin–orbit interactions potentially strengthen and accelerate the generation of LQFI–LQU correlation and log-negativity entanglement due to J α -spin–spin interactions.

Figure 3 
               Time evolution of the LQFI, LQU, and LN of Figure 2(c) are plotted for different two-spin couplings: 
                     
                        
                        
                           
                              (
                              
                                 
                                    
                                       J
                                    
                                    
                                       x
                                    
                                 
                                 ,
                                 
                                    
                                       J
                                    
                                    
                                       y
                                    
                                 
                                 ,
                                 
                                    
                                       J
                                    
                                    
                                       z
                                    
                                 
                              
                              )
                           
                        
                        \left({J}_{x},{J}_{y},{J}_{z})
                     
                   = (1, 0.5, 1.5) in (a) and 
                     
                        
                        
                           
                              (
                              
                                 
                                    
                                       J
                                    
                                    
                                       x
                                    
                                 
                                 ,
                                 
                                    
                                       J
                                    
                                    
                                       y
                                    
                                 
                                 ,
                                 
                                    
                                       J
                                    
                                    
                                       z
                                    
                                 
                              
                              )
                           
                        
                        \left({J}_{x},{J}_{y},{J}_{z})
                     
                   = (5, 1, 1.5) in (b). In (c), they are plotted for 
                     
                        
                        
                           x
                           ,
                           y
                        
                        x,y
                     
                   spin–orbit couplings 
                     
                        
                        
                           
                              
                                 D
                              
                              
                                 x
                              
                           
                           =
                           
                              
                                 D
                              
                              
                                 y
                              
                           
                           =
                           2
                        
                        {D}_{x}={D}_{y}=2
                     
                   in (c).
Figure 3

Time evolution of the LQFI, LQU, and LN of Figure 2(c) are plotted for different two-spin couplings: ( J x , J y , J z ) = (1, 0.5, 1.5) in (a) and ( J x , J y , J z ) = (5, 1, 1.5) in (b). In (c), they are plotted for x , y spin–orbit couplings D x = D y = 2 in (c).

Figure 4 illustrates the LQFI–LQU correlation and log-negativity entanglement dynamics of Figure 3(a) (where ( J x , J y , J z ) = (1, 0.5, 1.5), b m = 0.5 , and D x = D y = 0.5 ) for different uniformities of the applied EIMF. Figure 4(a) illustrates that with a large uniformity ( B f = 2 ), increasing the EIMF uniformity delays the growth of LQFI, LQU, and log-negativity. It also increases the fluctuations of the two-spin state between different partially and maximally correlated states. The generations of the LQFI–LQU correlation and log-negativity entanglement are shown in Figure 4(b) (with B m = 10 ), confirming that increasing the EIMF uniformity enhances the ability of strong J α -spin–spin interactions, supported by weak x , y -spin–orbit interactions, to create partially and maximally correlated states with greater stability. However, the generated spin–spin NLCs become more sensitive to the EIMF uniformity.

Figure 4 
               Time evolutions of LQFI, LQU, and LN of Figure 3(a) (for 
                     
                        
                        
                           
                              (
                              
                                 
                                    
                                       J
                                    
                                    
                                       x
                                    
                                 
                                 ,
                                 
                                    
                                       J
                                    
                                    
                                       y
                                    
                                 
                                 ,
                                 
                                    
                                       J
                                    
                                    
                                       z
                                    
                                 
                              
                              )
                           
                        
                        \left({J}_{x},{J}_{y},{J}_{z})
                     
                   = (1, 0.5, 1.5), 
                     
                        
                        
                           
                              
                                 b
                              
                              
                                 m
                              
                           
                           =
                           0.5
                        
                        {b}_{m}=0.5
                     
                  , and 
                     
                        
                        
                           
                              
                                 D
                              
                              
                                 x
                              
                           
                           =
                           
                              
                                 D
                              
                              
                                 y
                              
                           
                           =
                           0.5
                        
                        {D}_{x}={D}_{y}=0.5
                     
                  ) are plotted for different EIMF uniformities: 
                     
                        
                        
                           
                              
                                 B
                              
                              
                                 m
                              
                           
                           =
                           2
                        
                        {B}_{m}=2
                     
                   in (a) and 
                     
                        
                        
                           
                              
                                 B
                              
                              
                                 m
                              
                           
                           =
                           10
                        
                        {B}_{m}=10
                     
                   in (b).
Figure 4

Time evolutions of LQFI, LQU, and LN of Figure 3(a) (for ( J x , J y , J z ) = (1, 0.5, 1.5), b m = 0.5 , and D x = D y = 0.5 ) are plotted for different EIMF uniformities: B m = 2 in (a) and B m = 10 in (b).

In the upcoming analysis of Figure 5, we maintain the same parameter values as in Figure 3a (with ( J x , J y , J z ) = (1, 0.5, 1.5), B m = 0.3 , and D x = D y = 0.5 ) and examine different magnetic field inhomogeneities: b m = 2 in (a) and b m = 10 in (b). In the case of Figure 5(a), we observe that greater EIMF uniformities improve the efficiency of generating LQFI–LQU correlations and log-negativity entanglement. The uniformity of the EIMF increases the fluctuations of the two-spin state between different partially and maximally correlated states. The timing of the maxima ( F ( t ) = U ( t ) = N ( t ) 1 ) and minima (zero-value) ( F ( t ) = U ( t ) = N ( t ) 0 ) of the generated LQFI-LQU correlations and log-negativity entanglement is enhanced. Figure 5(b) demonstrates that an increase in EIMF inhomogeneity significantly enhances the generated two-spin–qubits’ NLCs, leading to greater amplitudes and fluctuations in the NLC oscillations.

Figure 5 
               Time evolutions of LQFI, LQU, and LN of Figure 3(a) (for 
                     
                        
                        
                           
                              (
                              
                                 
                                    
                                       J
                                    
                                    
                                       x
                                    
                                 
                                 ,
                                 
                                    
                                       J
                                    
                                    
                                       y
                                    
                                 
                                 ,
                                 
                                    
                                       J
                                    
                                    
                                       z
                                    
                                 
                              
                              )
                           
                        
                        \left({J}_{x},{J}_{y},{J}_{z})
                     
                   = (1, 0.5, 1.5), 
                     
                        
                        
                           
                              
                                 B
                              
                              
                                 m
                              
                           
                           =
                           0.3
                        
                        {B}_{m}=0.3
                     
                  , and 
                     
                        
                        
                           
                              
                                 D
                              
                              
                                 x
                              
                           
                           =
                           
                              
                                 D
                              
                              
                                 y
                              
                           
                           =
                           0.5
                        
                        {D}_{x}={D}_{y}=0.5
                     
                  ) are plotted for different EIMF inhomogeneities: 
                     
                        
                        
                           
                              
                                 b
                              
                              
                                 m
                              
                           
                           =
                           2
                        
                        {b}_{m}=2
                     
                   in (a) and 
                     
                        
                        
                           
                              
                                 b
                              
                              
                                 m
                              
                           
                           =
                           10
                        
                        {b}_{m}=10
                     
                   in (b).
Figure 5

Time evolutions of LQFI, LQU, and LN of Figure 3(a) (for ( J x , J y , J z ) = (1, 0.5, 1.5), B m = 0.3 , and D x = D y = 0.5 ) are plotted for different EIMF inhomogeneities: b m = 2 in (a) and b m = 10 in (b).

The next illustrations in Figures 6, 7, 8 depict the time evolutions of NLCs of LQFI, LQU, and log-negativity in the presence of non-zero ISSD coupling. By comparing the results of Figure 2(a) ( γ = 0.0 ) with those of Figure 6(a) ( γ = 0.05 ), we observe that LQFI, LQU, and log-negativity exhibit different decaying oscillatory dynamical evolutions. The generations of Heisenberg-XYZ states’ NLCs, due to J α = 0.8 spin–spin couplings and the applied magnetic field ( B m , b m ) = (0.3, 0.5) without spin–orbit interaction, are weakened and exhibit different amplitudes, which decrease with increasing ISSD coupling. After a certain time interval, with non-zero ISSD coupling, LQFI and LQU show different NLCs with varying amplitudes but similar behaviors. Moreover, the robustness of LQFI and log-negativity against the ISSD effect is greater for LQU.

Figure 6 
               Time evolutions of LQFI, LQU, and LN of Figure 2(a) are shown in the presence of the ISSD effect (
                     
                        
                        
                           γ
                           =
                           0.05
                        
                        \gamma =0.05
                     
                  ) with EIMF uniformity and inhomogeneity 
                     
                        
                        
                           
                              (
                              
                                 
                                    
                                       B
                                    
                                    
                                       m
                                    
                                 
                                 ,
                                 
                                    
                                       b
                                    
                                    
                                       m
                                    
                                 
                              
                              )
                           
                        
                        \left({B}_{m},{b}_{m})
                     
                   = (0.3, 0.5), and two-spin couplings 
                     
                        
                        
                           
                              
                                 J
                              
                              
                                 α
                              
                           
                           =
                           0.8
                        
                        {J}_{\alpha }=0.8
                     
                   for different couplings: 
                     
                        
                        
                           
                              
                                 D
                              
                              
                                 k
                              
                           
                           =
                           0
                        
                        {D}_{k}=0
                     
                   (
                     
                        
                        
                           k
                           =
                           x
                           ,
                           y
                        
                        k=x,y
                     
                  ) in (a), 
                     
                        
                        
                           
                              
                                 D
                              
                              
                                 k
                              
                           
                           =
                           0.5
                        
                        {D}_{k}=0.5
                     
                   in (b), and 
                     
                        
                        
                           
                              
                                 D
                              
                              
                                 k
                              
                           
                           =
                           2
                        
                        {D}_{k}=2
                     
                   in (c).
Figure 6

Time evolutions of LQFI, LQU, and LN of Figure 2(a) are shown in the presence of the ISSD effect ( γ = 0.05 ) with EIMF uniformity and inhomogeneity ( B m , b m ) = (0.3, 0.5), and two-spin couplings J α = 0.8 for different couplings: D k = 0 ( k = x , y ) in (a), D k = 0.5 in (b), and D k = 2 in (c).

Figure 7 
               Time evolutions of LQFI, LQU, and LN of Figure 6(b) and (c) are shown but for strong spin–spin couplings with 
                     
                        
                        
                           
                              (
                              
                                 
                                    
                                       J
                                    
                                    
                                       x
                                    
                                 
                                 ,
                                 
                                    
                                       J
                                    
                                    
                                       y
                                    
                                 
                                 ,
                                 
                                    
                                       J
                                    
                                    
                                       z
                                    
                                 
                              
                              )
                           
                        
                        \left({J}_{x},{J}_{y},{J}_{z})
                     
                   = (1, 0.5, 1.5).
Figure 7

Time evolutions of LQFI, LQU, and LN of Figure 6(b) and (c) are shown but for strong spin–spin couplings with ( J x , J y , J z ) = (1, 0.5, 1.5).

Figure 8 
               Time evolutions of LQFI, LQU, and LN in Figure 4(a), for 
                     
                        
                        
                           
                              (
                              
                                 
                                    
                                       J
                                    
                                    
                                       x
                                    
                                 
                                 ,
                                 
                                    
                                       J
                                    
                                    
                                       y
                                    
                                 
                                 ,
                                 
                                    
                                       J
                                    
                                    
                                       z
                                    
                                 
                              
                              )
                           
                        
                        \left({J}_{x},{J}_{y},{J}_{z})
                     
                   = (1, 0.5, 1.5) and 
                     
                        
                        
                           
                              
                                 D
                              
                              
                                 x
                              
                           
                           =
                           
                              
                                 D
                              
                              
                                 y
                              
                           
                           =
                           0.5
                        
                        {D}_{x}={D}_{y}=0.5
                     
                  , are presented considering the ISSD effect (
                     
                        
                        
                           γ
                           =
                           0.05
                        
                        \gamma =0.05
                     
                  ) for EIMF uniformity with 
                     
                        
                        
                           
                              (
                              
                                 
                                    
                                       B
                                    
                                    
                                       m
                                    
                                 
                                 ,
                                 
                                    
                                       b
                                    
                                    
                                       m
                                    
                                 
                              
                              )
                           
                        
                        \left({B}_{m},{b}_{m})
                     
                   = (2, 0.5) in (a) and EIMF inhomogeneity with 
                     
                        
                        
                           
                              (
                              
                                 
                                    
                                       B
                                    
                                    
                                       m
                                    
                                 
                                 ,
                                 
                                    
                                       b
                                    
                                    
                                       m
                                    
                                 
                              
                              )
                           
                        
                        \left({B}_{m},{b}_{m})
                     
                   = (0.3, 2) in (b).
Figure 8

Time evolutions of LQFI, LQU, and LN in Figure 4(a), for ( J x , J y , J z ) = (1, 0.5, 1.5) and D x = D y = 0.5 , are presented considering the ISSD effect ( γ = 0.05 ) for EIMF uniformity with ( B m , b m ) = (2, 0.5) in (a) and EIMF inhomogeneity with ( B m , b m ) = (0.3, 2) in (b).

As shown in Figure 6(b) and (c), increasing the intensities of x , y -spin–orbit interactions reduces the robustness of NLCs against the ISSD effect. The amplitudes of NLCs significantly decrease as the x , y -spin–orbit interactions increase. Moreover, LQFI and LQU exhibit sudden changes at different times. The phenomenon of sudden changes has been studied both theoretically [65] and experimentally [66]. For very strong x , y -spin–orbit interactions ( D k = 2 ) Figure 6(c), we observe that the log-negativity of the two-spin qubit drops instantly to zero at a specific time and remains zero for an extended period (i.e., the sudden-death LN-entanglement phenomenon occurs). After this, the disentangled two-spin states exhibit only different stable NLCs of LQFI and LQU. We can conclude that the decay of NLCs due to ISSD can be intensified by increasing the intensities of x , y -spin–orbit interactions.

Figure 7 illustrates the time evolutions of LQFI, LQU, and LN of Figure 6(b) and (c), but for strong spin–spin couplings with ( J x , J y , J z ) = (1, 0.5, 1.5). By comparing Figures 6(b), (c) and 7(a), (b), we find that strong spin–spin couplings ( J x , J y , J z ) = (1, 0.5, 1.5) reduce the ISSD effect and improve the robustness of the NLCs (against the ISSD effect) of LQFI, LQU, and LN. For very strong spin–orbit interactions D k = 2 (Figure 7(b)), the sudden-death LN-entanglement phenomenon does not occur, except instantaneously at t 0.5 π . The generated two-spin states have different stable partial NLCs of LQFI, LQU, and LN. In this case, the decay of NLCs due to ISSD can be mitigated by enhancing the spin–spin interactions.

Figure 8(a) shows the time evolutions of LQFI, LQU, and LN from Figure 4(a) for ( J x , J y , J z ) = (1, 0.5, 1.5) and D x = D y = 0.5 , considering the ISSD effect ( γ = 0.05 ) after strengthening the EIMF uniformity with ( B m , b m ) = (2, 0.5). In this case, we observe that the large EIMF uniformity B m = 2 increases the NLCs’ decay resulting from ISSD. Time intervals appear in which the disentangled two-spin states have only different stable NLCs of LQFI and LQU. Moreover, the robustness of LQFI and LN NLCs, against the ISSD effect, is reduced by increasing EIMF uniformity. The results shown in Figure 8(b) demonstrate that increasing the inhomogeneity of the EIMF to b m = 2 enhances the degradation of the NLC functions. Under the parameters b m = 2 , ( J x , J y , J z ) = (1, 0.5, 1.5), and D k = 0.5 , we observe that the generated NLCs (LQFI, LQU, and entanglement) in Figure 4(a) degrade (due to the ISSD effect) and quickly reach their partially stable oscillatory behaviors, compared to the case with a small value of b m = 0.5 in Figure 7(a). We find that the EIMF’s inhomogeneity has a lesser ability to enhance the ISSD effect compared to its uniformity.

5 Conclusion

In this study, we use the Milburn intrinsic decoherence and the Heisenberg XYZ models to explore the capabilities of spin–spin and spin–orbit interactions (in the x and y directions) to generate NLCs (measured by LQFI, LQU, and LN) under the influence of the uniformity and inhomogeneity of the EIMF in the x direction. The generated NLCs are examined in the absence of the ISSD, as the parameters of spin–spin and spin–orbit interactions, as well as the EIMF’s uniformity and inhomogeneity, are increased. It is found that the spin–spin Heisenberg XYZ and x , y -spin–orbit interactions have a strong ability to enhance non-local correlations with small external magnetic field parameters. The spin–orbit interactions significantly contribute to the enhancement of the generated two-spin–qubits LQFI–LQU correlation and log-negativity entanglement, increasing their oscillation amplitudes and fluctuations. In the presence the ISSD, the generation of NLCs is weakened and exhibits varying amplitudes, decreasing as the ISSD coupling increases. The robustness of LQFI and log-negativity against the ISSD effect is greater than that of LQU. Sudden changes occur during the dynamics of LQU and LQFI, while sudden death occurs during the dynamics of log-negativity entanglement. The decay of NLCs due to ISSD can be enhanced by increasing the intensities of x , y -spin–orbit interactions. Strengthening the spin–spin interactions, however, weakens the NLCs’ decay resulting from ISSD. The generated NLCs degrade quickly due to the ISSD effect with a large IMF’s inhomogeneity, reaching their partially stable oscillatory behaviors. The ability of the IMF’s inhomogeneity to increase the ISSD effect is small compared to that of the EIMF’s uniformity. Our findings on the generated two-spin qubit NLCs (measured by LQFI, LQU, and LN) resulting from spin–spin interactions, supported by spin–orbit interactions in the x and y directions and an applied EIMF in the x direction, pave the way for exploring other quantum effects with the considered interaction directions or with different interaction directions. Additionally, these quantum effects can be explored using other differential motion equations with the considered interaction directions. Furthermore, fractional quantum calculus and Riemann–Liouville integration [67,68] can be employed to investigate fractional quantum phenomena in the considered Heisenberg XYZ model using fractional differential motion equations [69,70].

Acknowledgments

The author is very grateful to the referees and the associate editor for their important remarks which have helped him to improve the manuscript. This study is supported via funding from Prince sattam bin Abdulaziz University Project Number (PSAU/2025/R/1446).

  1. Funding information: This study was funded by Princess Nourah bint Abdulrahman University Researchers Supporting Project Number (PNURSP2025R906), Princess Nourah bint Abdulrahman University, Riyadh, Saudi, Arabia.

  2. Author contributions: Sattam N. Almutairi, Heba Allhibi, and Fahad Aljuaydi prepared all the figures and performed the mathematical calculations. Heba Allhibi, Sattam N. Almutairi, and Fahad Aljuaydi wrote the original draft. Abdel-Baset A. Mohamed, Sattam N. Almutairi, and Fahad Aljuaydi reviewed and edited the draft. All authors have accepted responsibility for the entire content of this manuscript and approved its submission.

  3. Conflict of interest: The authors state no conflict of interest.

  4. Data availability statement: The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.

References

[1] Horodecki R, Horodecki P, Horodecki M, Horodecki K. Quantum entanglement. Rev Mod Phys. 2009;81:865. 10.1103/RevModPhys.81.865Search in Google Scholar

[2] Nielsen MA, Chuang IL. Quantum computation and quantum information. Cambridge: Cambridge University Press; 2000. Search in Google Scholar

[3] Loss D, DiVincenzo DP. Quantum computation with quantum dots. Phys Rev A 1998;57:120. 10.1103/PhysRevA.57.120Search in Google Scholar

[4] Burkard G, Loss D, DiVincenzo DP. Coupled quantum dots as quantum gates. Phys Rev B. 1999;59:2070. 10.1103/PhysRevB.59.2070Search in Google Scholar

[5] Vandersypen LMK, Hanson R, van Willems Beveren LH, Elzerman JM, Greidanus JS, De Franceschi S, et al. Quantum computing and quantum bits in mesoscopic systems. Boston, M A.: Springer; 2004. Search in Google Scholar

[6] Ma R-L, Li A-R, Wang C, Kong ZZ, Liao W-Z, Ni M, et al. Single-spin–qubit geometric gate in a silicon quantum dot. Phys Rev Appl. 2024;21:014044. 10.1103/PhysRevApplied.21.014044Search in Google Scholar

[7] Zwolak JP, Taylor JM. Colloquium: Advances in automation of quantum dot devices control. Rev Mod Phys. 2023;95:011006. 10.1103/RevModPhys.95.011006Search in Google Scholar PubMed PubMed Central

[8] Pinheiro F, Bruun GM, Martikainen J-P, Larson J. XYZ Quantum Heisenberg Models with p-Orbital Bosons. Phys Rev Lett. 2013;111:205302. 10.1103/PhysRevLett.111.205302Search in Google Scholar PubMed

[9] Bermudez A, Tagliacozzo L, Sierra G, Richerme P. Long-range Heisenberg models in quasiperiodically driven crystals of trapped ions. Phys Rev B Condens Matter Mater Phys. 2017;95:024431. 10.1103/PhysRevB.95.024431Search in Google Scholar

[10] Nishiyama M, Inada Y, Zheng GQ. Spin triplet superconducting state due to broken inversion symmetry in Li2Pt3B. Phys Rev Lett. 2007;98:047002. 10.1103/PhysRevLett.98.047002Search in Google Scholar PubMed

[11] Yue W, Wei Q, Kais S, Friedrichc B, Herschbach D. Realization of Heisenberg models of spin systems with polar molecules in pendular states. Phys Chem Chem Phys. 2022;24:25270. 10.1039/D2CP00380ESearch in Google Scholar PubMed

[12] Dzyaloshinski I. A thermodynamic theory of “weak” ferromagnetism of antiferromagnetics. J Phys Chem Solids. 1958;4:241. 10.1016/0022-3697(58)90076-3Search in Google Scholar

[13] Moriya T. Theory of Magnetism of NiF2. Phys Rev. 1960;117:635. 10.1103/PhysRev.117.635Search in Google Scholar

[14] Moriya T. New mechanism of anisotropic superexchange interaction. Phys Rev Lett. 1960;4:228. 10.1103/PhysRevLett.4.228Search in Google Scholar

[15] Shekhtman L, Entin-Wohlman O, Aharony A. Bond-dependent symmetric and antisymmetric superexchange interactions in La2CuO4. Phys Rev B. 1993;47:174. 10.1103/PhysRevB.47.174Search in Google Scholar

[16] Moriya T. Anisotropic superexchange interaction and weak Ferromagnetism. Phys Rev. 1960;120:91. 10.1103/PhysRev.120.91Search in Google Scholar

[17] Kuznetsova EI, Yurischev MA. Quantum discord in spin systems with dipole-dipole interaction. Quantum Inf Process. 2013;12:3587–605. 10.1007/s11128-013-0617-6Search in Google Scholar

[18] Arnesen MC, Bose S, Vedral V. Natural thermal and magnetic entanglement in the 1D Heisenberg model. Phys Rev Lett. 2001;87:017901. 10.1103/PhysRevLett.87.017901Search in Google Scholar PubMed

[19] Li D-C, Cao Z-L. Thermal entanglement in the anisotropic Heisenberg (XYZ) model with different inhomogeneous magnetic fields. Optics Commun. 2009;282:1226–30. 10.1016/j.optcom.2008.11.087Search in Google Scholar

[20] Matthies T, Rózsa L, Szunyogh L, Wiesendanger R, Vedmedenko EY. Interlayer and interfacial Dzyaloshinskii-Moriya interaction in magnetic trilayers: First-principles calculations. Phys Rev Res. 2024;6:043158. 10.1103/PhysRevResearch.6.043158Search in Google Scholar

[21] Thoma H, Hutanu V, Deng H, Dmitrienko VE, Brown PJ, Gukasov A, et al. Revealing the absolute direction of the Dzyaloshinskii-Moriya interaction in prototypical weak ferromagnets by polarized neutrons. Phys Rev X. 2021;11:011060. 10.1103/PhysRevX.11.011060Search in Google Scholar

[22] Zhao HJ, Chen P, Prosandeev S, Artyukhin S, Bellaiche L. Dzyaloshinskii-Moriya-like interaction in ferroelectrics and antiferroelectrics. Nature Materials. 2021;20:341. 10.1038/s41563-020-00821-3Search in Google Scholar PubMed

[23] Allenspach R, Bischof A, Boehm B, Drechsler U, Reich O, Sousa M, et al. Dzyaloshinskii-Moriya interaction in Ni/Cu(001). Phys Rev B. 2024;110:014402. 10.1103/PhysRevB.110.014402Search in Google Scholar

[24] Hosseiny SM, Seyed-Yazdi J, Norouzi M, Livreri P. Quantum teleportation in Heisenberg chain with magnetic-field gradient under intrinsic decoherence. Sci Rep. 2024;14:9607. 10.1038/s41598-024-60321-1Search in Google Scholar PubMed PubMed Central

[25] Khedif Y, Muthuganesan R. Intrinsic decoherence dynamics and dense coding in dipolar spin system. Appl Phys B. 2023;129:19. 10.1007/s00340-022-07956-ySearch in Google Scholar

[26] Abd-Rabbou MY, ur Rahman A, Yurischev MA, Haddadi S. Comparative study of quantum Otto and Carnot engines powered by a spin working substance. Phys Rev E. 2023;108:034106. 10.1103/PhysRevE.108.034106Search in Google Scholar PubMed

[27] Yurischev MA. On the quantum correlations in two-qubit XYZ spin chains with Dzyaloshinsky–Moriya and Kaplan–Shekhtman–Entin-Wohlman–Aharony interactions. Quantum Inf Process. 2020;19:336. 10.1007/s11128-020-02835-xSearch in Google Scholar

[28] Le Bellac M. A short introduction to quantum information and quantum computation. Cambridge: Cambridge University Press; 2006. 10.1017/CBO9780511755361Search in Google Scholar

[29] Eisert J, Cramer M, Plenio MB. Colloquium: Area laws for the entanglement entropy. Rev Mod Phys. 2010;82:277. 10.1103/RevModPhys.82.277Search in Google Scholar

[30] Wootters WK. Entanglement of formation of an arbitrary state of two qubits. Phys Rev Lett. 1998;80:2245. 10.1103/PhysRevLett.80.2245Search in Google Scholar

[31] Vidal G, Werner RF. Computable measure of entanglement. Phys Rev A. 2002;65:032314. 10.1103/PhysRevA.65.032314Search in Google Scholar

[32] Eftekhari F, Tavassoly MK, Behjat A, Faghihi MJ. Entanglement and atomic inversion in a dissipative two-atom-optomechanical system. Optics Laser Tech. 2024;168:109934. 10.1016/j.optlastec.2023.109934Search in Google Scholar

[33] Abdel-Wahab NH, Zangi SM, Seoudy TA, Haddadi S. Dynamical evolution of a five-level atom interacting with an intensity-dependent coupling regime influenced by a nonlinear Kerr-like medium. Sci Rep. 2024;14:25211. 10.1038/s41598-024-76629-xSearch in Google Scholar PubMed PubMed Central

[34] Braunstein SL, Kimble HJ. Teleportation of continuous quantum variables. Phys Rev Lett. 1998;80:869. 10.1103/PhysRevLett.80.869Search in Google Scholar

[35] Bouwmeester D, Pan JW, Mattle K, Eibl M, Weinfurter H, Zeilinger A. Experimental quantum teleportation. Nature. 1997;390:575. 10.1038/37539Search in Google Scholar

[36] Lei Y, Asadi FK, Zhong T, Kuzmich A, Simon C, Hosseini M. Quantum optical memory for entanglement distribution. Optica. 2023;10:1511–28. 10.1364/OPTICA.493732Search in Google Scholar

[37] Ekert AK. Quantum cryptography based on Bell’s theorem. Phys Rev Lett. 1991;67:661. 10.1103/PhysRevLett.67.661Search in Google Scholar PubMed

[38] Ollivier H, Zurek WH. Quantum discord: A measure of the quantumness of correlations. Phys Rev Lett. 2001;88:017901. 10.1103/PhysRevLett.88.017901Search in Google Scholar PubMed

[39] Hu M-L, Hu X, Wang J, Peng Y, Zhang Y-R, Fan H. Quantum coherence and geometric quantum discord. Phys Rep. 2018;762:1. 10.1016/j.physrep.2018.07.004Search in Google Scholar

[40] Wigner EP, Yanase MM. Information contents of distributions. Proc Natl Acad Sci. 1963;49:910. 10.1073/pnas.49.6.910Search in Google Scholar PubMed PubMed Central

[41] Girolami D, Souza AM, Giovannetti V, Tufarelli T, Filgueiras JG, Sarthour RS, et al. Quantum discord determines the interferometric power of quantum states. Phys Rev Lett. 2014;112:210401. 10.1103/PhysRevLett.112.210401Search in Google Scholar

[42] Girolami D, Tufarelli T, Adesso G. Characterizing nonclassical correlations via local quantum uncertainty. Phys Rev Lett. 2013:110:240402. 10.1103/PhysRevLett.110.240402Search in Google Scholar PubMed

[43] Wu S-X, Zhang J, Yuuu C-S, Song H-S. Uncertainty-induced quantum nonlocality. Phys Lett A. 2014;378:344. 10.1016/j.physleta.2013.11.047Search in Google Scholar

[44] Dhar HS, Bera MN, Adesso G. Characterizing non-Markovianity via quantum interferometric power. Phys Rev A. 2015;991:032115. 10.1103/PhysRevA.91.032115Search in Google Scholar

[45] Kim S, Li L, Kumar A, Wu J. Characterizing nonclassical correlations via local quantum Fisher information. Phys Rev A. 2018;97:032326. 10.1103/PhysRevA.97.032326Search in Google Scholar

[46] Helstrom CW. Quantum detection and estimation theory. New York: Academic; 1976. Search in Google Scholar

[47] Paris MGA. Quantum estimation for quantum technology. Int J Quant Inf. 2009;7:125. 10.1142/S0219749909004839Search in Google Scholar

[48] Mohamed A-B A, Farouk A, Yassen MF, Eleuch H. Fisher and skew information correlations of two coupled trapped ions: Intrinsic decoherence and Lamb-Dicke nonlinearity. Symmetry. 2021;13:2243. 10.3390/sym13122243Search in Google Scholar

[49] Benabdallah F, Ur Rahman A, Haddadi S, Daoud M. Long-time protection of thermal correlations in a hybrid-spin system under random telegraph noise. Phys Rev E. 2022;106:034122. 10.1103/PhysRevE.106.034122Search in Google Scholar PubMed

[50] Benabdallah F, El Anouz K, Ur Rahman A, Daoud M, El Allati A, Haddadi S. Witnessing quantum correlations in a hybrid Qubit-Qutrit system under intrinsic decoherence. Fortschr Phys. 2023;71:2300032. 10.1002/prop.202300032Search in Google Scholar

[51] Elghaayda S, Dahbi Z, Mansour M. Local quantum uncertainty and local quantum Fisher information in two-coupled double quantum dots. Opt Quant Electron 2022;54:419. 10.1007/s11082-022-03829-ySearch in Google Scholar

[52] Wei P-F, Luo Q, Wang H-Q-C, Xiong S-J, Liu B, Sun Z. Local quantum Fisher information and quantum correlation in the mixed-spin Heisenberg XXZ chain. Front Phys. 2024;19:21201. 10.1007/s11467-023-1336-9Search in Google Scholar

[53] Fedorova AV, Yurischev MA. Behavior of quantum discord, local quantum uncertainty, and local quantum Fisher information in two-spin-1/2 Heisenberg chain with DM and KSEA interactions. Quantum Inf Process. 2022;21:92. 10.1007/s11128-022-03427-7Search in Google Scholar

[54] Milburn GJ. Intrinsic decoherence in quantum mechanics. Phys Rev A. 1991;44:5401. 10.1103/PhysRevA.44.5401Search in Google Scholar PubMed

[55] Qin M, Ren Z-Z. Influence of intrinsic decoherence on entanglement teleportation via a Heisenberg XYZ model with different Dzyaloshinskii-Moriya interactions. Quantum Inf Process. 2015;14:2055–66. 10.1007/s11128-015-0978-0Search in Google Scholar

[56] Mohamed A-B A, Aldosari FM, Eleuch H. Two-qubit-Heisenberg local quantum Fisher information dynamics induced by intrinsic decoherence model. Results Phys. 2023;49:106470. 10.1016/j.rinp.2023.106470Search in Google Scholar

[57] Ait Chlih A, Habiballah N, Nassik M. Dynamics of quantum correlations under intrinsic decoherence in a Heisenberg spin chain model with Dzyaloshinskii–Moriya interaction. Quantum Inf Process. 2021;20:92. 10.1007/s11128-021-03030-2Search in Google Scholar

[58] Jafari R, Akbari A. Dynamics of quantum coherence and quantum Fisher information after a sudden quench. Phys Rev A. 2020;101:062105. 10.1103/PhysRevA.101.062105Search in Google Scholar

[59] Mo C, Zhang G-F. The effect of acceleration parameter on thermal entanglement and teleportation of a two-qubit Heisenberg XXX model with Dzyaloshinski-Moriya interaction in non-inertial frames. Results Phys. 2021;21:103759. 10.1016/j.rinp.2020.103759Search in Google Scholar

[60] Indrajith VS, Sankaranarayanan R. Effect of environment on quantum correlation in Heisenberg XYZ spin model. Phys A. 2021;582:126250. 10.1016/j.physa.2021.126250Search in Google Scholar

[61] El Aroui A, Khedif Y, Habiballah N, Nassik M, Aroui AE, Khedif Y, et al. Characterizing the thermal quantum correlations in a two-qubit Heisenberg XXZ spin-chain 1/2 under Dzyaloshinskii-Moriya and Kaplan-Shekhtman-Entin-Wohlman-Aharony interactions. Opt Quantum Elec. 2022;54:694. 10.1007/s11082-022-04042-7Search in Google Scholar

[62] Yurischev MA, Haddadi S. Local quantum Fisher information and local quantum uncertainty for general X states. Phys Lett A. 2023;476:128868. 10.1016/j.physleta.2023.128868Search in Google Scholar

[63] Hosseiny SM. Quantum teleportation and phase quantum estimation in a two-qubit state influenced by dipole and symmetric cross interactions. Phys Scr. 2023;98(11):115101. 10.1088/1402-4896/acfc7aSearch in Google Scholar

[64] Mohr PJ, Phillips WD. Dimensionless units in the SI. Metrologia. 2014;52:40. 10.1088/0026-1394/52/1/40Search in Google Scholar

[65] Maziero J, Celeri LC, Serra RM, Vedral V. Classical and quantum correlations under decoherence. Phys Rev A. 2009;80:044102. 10.1103/PhysRevA.80.044102Search in Google Scholar

[66] Xu J-S, Xu X-Y, Li C-F, Zhang C-J, Zou X-B, Guo G-C. Experimental investigation of classical and quantum correlations under decoherence. Nature Commun. 2010;1:7. 10.1038/ncomms1005Search in Google Scholar PubMed

[67] Tariboon J, Ntouyas SK, Agarwal P. New concepts of fractional quantum calculus and applications to impulsive fractional q-difference equations. Adv Differ Equ. 2015;2015:18. 10.1186/s13662-014-0348-8Search in Google Scholar

[68] Tuncay Gençoǧlu M, Agarwal P. Use of quantum differential equations in sonic processes. App Math Nonl Sci. 2021;6(1):21–8. 10.2478/amns.2020.2.00003Search in Google Scholar

[69] El Allati A, Bukbech S, El Anouz K, El Allali Z. Entanglement versus Bell non-locality via solving the fractional Schrödinger equation using the twisting model. Chaos Solitons Fractals. 2024;179:114446. 10.1016/j.chaos.2023.114446Search in Google Scholar

[70] El Anouz K, El Allati A, Metwally N, Obada AS. The efficiency of fractional channels in the Heisenberg XYZ model. Chaos Solitons Fractals. 2023;172:113581. 10.1016/j.chaos.2023.113581Search in Google Scholar

Received: 2024-11-15
Revised: 2025-01-10
Accepted: 2025-01-28
Published Online: 2025-09-16

© 2025 the author(s), published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

  1. Research Articles
  2. Single-step fabrication of Ag2S/poly-2-mercaptoaniline nanoribbon photocathodes for green hydrogen generation from artificial and natural red-sea water
  3. Abundant new interaction solutions and nonlinear dynamics for the (3+1)-dimensional Hirota–Satsuma–Ito-like equation
  4. A novel gold and SiO2 material based planar 5-element high HPBW end-fire antenna array for 300 GHz applications
  5. Explicit exact solutions and bifurcation analysis for the mZK equation with truncated M-fractional derivatives utilizing two reliable methods
  6. Optical and laser damage resistance: Role of periodic cylindrical surfaces
  7. Numerical study of flow and heat transfer in the air-side metal foam partially filled channels of panel-type radiator under forced convection
  8. Water-based hybrid nanofluid flow containing CNT nanoparticles over an extending surface with velocity slips, thermal convective, and zero-mass flux conditions
  9. Dynamical wave structures for some diffusion--reaction equations with quadratic and quartic nonlinearities
  10. Solving an isotropic grey matter tumour model via a heat transfer equation
  11. Study on the penetration protection of a fiber-reinforced composite structure with CNTs/GFP clip STF/3DKevlar
  12. Influence of Hall current and acoustic pressure on nanostructured DPL thermoelastic plates under ramp heating in a double-temperature model
  13. Applications of the Belousov–Zhabotinsky reaction–diffusion system: Analytical and numerical approaches
  14. AC electroosmotic flow of Maxwell fluid in a pH-regulated parallel-plate silica nanochannel
  15. Interpreting optical effects with relativistic transformations adopting one-way synchronization to conserve simultaneity and space–time continuity
  16. Modeling and analysis of quantum communication channel in airborne platforms with boundary layer effects
  17. Theoretical and numerical investigation of a memristor system with a piecewise memductance under fractal–fractional derivatives
  18. Tuning the structure and electro-optical properties of α-Cr2O3 films by heat treatment/La doping for optoelectronic applications
  19. High-speed multi-spectral explosion temperature measurement using golden-section accelerated Pearson correlation algorithm
  20. Dynamic behavior and modulation instability of the generalized coupled fractional nonlinear Helmholtz equation with cubic–quintic term
  21. Study on the duration of laser-induced air plasma flash near thin film surface
  22. Exploring the dynamics of fractional-order nonlinear dispersive wave system through homotopy technique
  23. The mechanism of carbon monoxide fluorescence inside a femtosecond laser-induced plasma
  24. Numerical solution of a nonconstant coefficient advection diffusion equation in an irregular domain and analyses of numerical dispersion and dissipation
  25. Numerical examination of the chemically reactive MHD flow of hybrid nanofluids over a two-dimensional stretching surface with the Cattaneo–Christov model and slip conditions
  26. Impacts of sinusoidal heat flux and embraced heated rectangular cavity on natural convection within a square enclosure partially filled with porous medium and Casson-hybrid nanofluid
  27. Stability analysis of unsteady ternary nanofluid flow past a stretching/shrinking wedge
  28. Solitonic wave solutions of a Hamiltonian nonlinear atom chain model through the Hirota bilinear transformation method
  29. Bilinear form and soltion solutions for (3+1)-dimensional negative-order KdV-CBS equation
  30. Solitary chirp pulses and soliton control for variable coefficients cubic–quintic nonlinear Schrödinger equation in nonuniform management system
  31. Influence of decaying heat source and temperature-dependent thermal conductivity on photo-hydro-elasto semiconductor media
  32. Dissipative disorder optimization in the radiative thin film flow of partially ionized non-Newtonian hybrid nanofluid with second-order slip condition
  33. Bifurcation, chaotic behavior, and traveling wave solutions for the fractional (4+1)-dimensional Davey–Stewartson–Kadomtsev–Petviashvili model
  34. New investigation on soliton solutions of two nonlinear PDEs in mathematical physics with a dynamical property: Bifurcation analysis
  35. Mathematical analysis of nanoparticle type and volume fraction on heat transfer efficiency of nanofluids
  36. Creation of single-wing Lorenz-like attractors via a ten-ninths-degree term
  37. Optical soliton solutions, bifurcation analysis, chaotic behaviors of nonlinear Schrödinger equation and modulation instability in optical fiber
  38. Chaotic dynamics and some solutions for the (n + 1)-dimensional modified Zakharov–Kuznetsov equation in plasma physics
  39. Fractal formation and chaotic soliton phenomena in nonlinear conformable Heisenberg ferromagnetic spin chain equation
  40. Single-step fabrication of Mn(iv) oxide-Mn(ii) sulfide/poly-2-mercaptoaniline porous network nanocomposite for pseudo-supercapacitors and charge storage
  41. Novel constructed dynamical analytical solutions and conserved quantities of the new (2+1)-dimensional KdV model describing acoustic wave propagation
  42. Tavis–Cummings model in the presence of a deformed field and time-dependent coupling
  43. Spinning dynamics of stress-dependent viscosity of generalized Cross-nonlinear materials affected by gravitationally swirling disk
  44. Design and prediction of high optical density photovoltaic polymers using machine learning-DFT studies
  45. Robust control and preservation of quantum steering, nonlocality, and coherence in open atomic systems
  46. Coating thickness and process efficiency of reverse roll coating using a magnetized hybrid nanomaterial flow
  47. Dynamic analysis, circuit realization, and its synchronization of a new chaotic hyperjerk system
  48. Decoherence of steerability and coherence dynamics induced by nonlinear qubit–cavity interactions
  49. Finite element analysis of turbulent thermal enhancement in grooved channels with flat- and plus-shaped fins
  50. Modulational instability and associated ion-acoustic modulated envelope solitons in a quantum plasma having ion beams
  51. Statistical inference of constant-stress partially accelerated life tests under type II generalized hybrid censored data from Burr III distribution
  52. On solutions of the Dirac equation for 1D hydrogenic atoms or ions
  53. Entropy optimization for chemically reactive magnetized unsteady thin film hybrid nanofluid flow on inclined surface subject to nonlinear mixed convection and variable temperature
  54. Stability analysis, circuit simulation, and color image encryption of a novel four-dimensional hyperchaotic model with hidden and self-excited attractors
  55. A high-accuracy exponential time integration scheme for the Darcy–Forchheimer Williamson fluid flow with temperature-dependent conductivity
  56. Novel analysis of fractional regularized long-wave equation in plasma dynamics
  57. Development of a photoelectrode based on a bismuth(iii) oxyiodide/intercalated iodide-poly(1H-pyrrole) rough spherical nanocomposite for green hydrogen generation
  58. Investigation of solar radiation effects on the energy performance of the (Al2O3–CuO–Cu)/H2O ternary nanofluidic system through a convectively heated cylinder
  59. Quantum resources for a system of two atoms interacting with a deformed field in the presence of intensity-dependent coupling
  60. Studying bifurcations and chaotic dynamics in the generalized hyperelastic-rod wave equation through Hamiltonian mechanics
  61. A new numerical technique for the solution of time-fractional nonlinear Klein–Gordon equation involving Atangana–Baleanu derivative using cubic B-spline functions
  62. Interaction solutions of high-order breathers and lumps for a (3+1)-dimensional conformable fractional potential-YTSF-like model
  63. Hydraulic fracturing radioactive source tracing technology based on hydraulic fracturing tracing mechanics model
  64. Numerical solution and stability analysis of non-Newtonian hybrid nanofluid flow subject to exponential heat source/sink over a Riga sheet
  65. Numerical investigation of mixed convection and viscous dissipation in couple stress nanofluid flow: A merged Adomian decomposition method and Mohand transform
  66. Effectual quintic B-spline functions for solving the time fractional coupled Boussinesq–Burgers equation arising in shallow water waves
  67. Analysis of MHD hybrid nanofluid flow over cone and wedge with exponential and thermal heat source and activation energy
  68. Solitons and travelling waves structure for M-fractional Kairat-II equation using three explicit methods
  69. Impact of nanoparticle shapes on the heat transfer properties of Cu and CuO nanofluids flowing over a stretching surface with slip effects: A computational study
  70. Computational simulation of heat transfer and nanofluid flow for two-sided lid-driven square cavity under the influence of magnetic field
  71. Irreversibility analysis of a bioconvective two-phase nanofluid in a Maxwell (non-Newtonian) flow induced by a rotating disk with thermal radiation
  72. Hydrodynamic and sensitivity analysis of a polymeric calendering process for non-Newtonian fluids with temperature-dependent viscosity
  73. Exploring the peakon solitons molecules and solitary wave structure to the nonlinear damped Kortewege–de Vries equation through efficient technique
  74. Modeling and heat transfer analysis of magnetized hybrid micropolar blood-based nanofluid flow in Darcy–Forchheimer porous stenosis narrow arteries
  75. Activation energy and cross-diffusion effects on 3D rotating nanofluid flow in a Darcy–Forchheimer porous medium with radiation and convective heating
  76. Insights into chemical reactions occurring in generalized nanomaterials due to spinning surface with melting constraints
  77. Influence of a magnetic field on double-porosity photo-thermoelastic materials under Lord–Shulman theory
  78. Soliton-like solutions for a nonlinear doubly dispersive equation in an elastic Murnaghan's rod via Hirota's bilinear method
  79. Analytical and numerical investigation of exact wave patterns and chaotic dynamics in the extended improved Boussinesq equation
  80. Nonclassical correlation dynamics of Heisenberg XYZ states with (x, y)-spin--orbit interaction, x-magnetic field, and intrinsic decoherence effects
  81. Exact traveling wave and soliton solutions for chemotaxis model and (3+1)-dimensional Boiti–Leon–Manna–Pempinelli equation
  82. Unveiling the transformative role of samarium in ZnO: Exploring structural and optical modifications for advanced functional applications
  83. On the derivation of solitary wave solutions for the time-fractional Rosenau equation through two analytical techniques
  84. Analyzing the role of length and radius of MWCNTs in a nanofluid flow influenced by variable thermal conductivity and viscosity considering Marangoni convection
  85. Advanced mathematical analysis of heat and mass transfer in oscillatory micropolar bio-nanofluid flows via peristaltic waves and electroosmotic effects
  86. Exact bound state solutions of the radial Schrödinger equation for the Coulomb potential by conformable Nikiforov–Uvarov approach
  87. Some anisotropic and perfect fluid plane symmetric solutions of Einstein's field equations using killing symmetries
  88. Nonlinear dynamics of the dissipative ion-acoustic solitary waves in anisotropic rotating magnetoplasmas
  89. Curves in multiplicative equiaffine plane
  90. Exact solution of the three-dimensional (3D) Z2 lattice gauge theory
  91. Propagation properties of Airyprime pulses in relaxing nonlinear media
  92. Symbolic computation: Analytical solutions and dynamics of a shallow water wave equation in coastal engineering
  93. Wave propagation in nonlocal piezo-photo-hygrothermoelastic semiconductors subjected to heat and moisture flux
  94. Comparative reaction dynamics in rotating nanofluid systems: Quartic and cubic kinetics under MHD influence
  95. Laplace transform technique and probabilistic analysis-based hypothesis testing in medical and engineering applications
  96. Physical properties of ternary chloro-perovskites KTCl3 (T = Ge, Al) for optoelectronic applications
  97. Gravitational length stretching: Curvature-induced modulation of quantum probability densities
  98. Review Article
  99. Examination of the gamma radiation shielding properties of different clay and sand materials in the Adrar region
  100. Special Issue on Fundamental Physics from Atoms to Cosmos - Part II
  101. Possible explanation for the neutron lifetime puzzle
  102. Special Issue on Nanomaterial utilization and structural optimization - Part III
  103. Numerical investigation on fluid-thermal-electric performance of a thermoelectric-integrated helically coiled tube heat exchanger for coal mine air cooling
  104. Special Issue on Nonlinear Dynamics and Chaos in Physical Systems
  105. Analysis of the fractional relativistic isothermal gas sphere with application to neutron stars
  106. Abundant wave symmetries in the (3+1)-dimensional Chafee–Infante equation through the Hirota bilinear transformation technique
  107. Successive midpoint method for fractional differential equations with nonlocal kernels: Error analysis, stability, and applications
  108. Novel exact solitons to the fractional modified mixed-Korteweg--de Vries model with a stability analysis
Downloaded on 2.11.2025 from https://www.degruyterbrill.com/document/doi/10.1515/phys-2025-0143/html
Scroll to top button