Unconventional Thermal and Magnetic-Field-Driven Changes of a Bipartite Entanglement of a Mixed Spin-(1/2,S) Heisenberg Dimer with an Uniaxial Single-Ion Anisotropy

The concept of negativity is adapted in order to explore the quantum and thermal entanglement of the mixed spin-(1/2,S) Heisenberg dimers in presence of an external magnetic field. The mutual interplay between the spin size S, XXZ exchange and uniaxial single-ion anisotropy is thoroughly examined with a goal to tune the degree and thermal stability of the pairwise entanglement. It turns out that the antiferromagnetic spin-(1/2,S) Heisenberg dimers exhibit higher degree of entanglement and higher threshold temperature in comparison with their ferromagnetic counterparts when assuming the same set of model parameters. The increasing spin magnitude S accompanied with an easy-plane uniaxial single-ion anisotropy can enhance not only the thermal stability but simultaneously the degree of entanglement. It is additionally shown that the further enhancement of a bipartite entanglement can be achieved in the mixed spin-(1/2,S) Heisenberg dimers, involving half-odd-integer spins S. Under this condition the thermal negativity saturates at low-enough temperatures in its maximal value regardless of the magnitude of half-odd-integer spin S. The magnetic field induces consecutive discontinuous phase transitions in the mixed spin-(1/2,S) Heisenberg dimers with S>1, which are manifested in a surprising oscillating magnetic-field dependence of the negativity observed at low enough temperature.


Introduction
Extraordinary correlations between subsystems of a quantum-mechanical ensemble, known as entanglement, belongs to the most fascinating phenomena attracting a lot of attention during the last few decades. A huge concerment in this field of study closely relates to perspective applications of this phenomenon in a quantum computing [1][2][3][4], quantum information [5,6], and quantum memory circuits decoding [7][8][9]. However, a long period before, it was believed that the entanglement could exist exclusively on the atomic scale and completely vanishes at macroscopic scales as a consequence of the decoherence arising from the interactions between a large number of matter constituents and their environment. It was presumed, in addition, that arbitrarily small thermal fluctuations rapidly smear out the quantum correlations and thus, the entanglement cannot exist at nonzero temperatures. Some theoretical predictions [10][11][12][13][14][15] have presented relevant arguments that the entanglement could surprisingly exist even at finite temperatures, but rapidly falls down as temperature increases.
From a theoretical perspective low-dimensional Heisenberg spin models seem to be a reasonable theoretical ground, which allows an exact study of the quantum and thermal entanglement depending on external stimuli such as the magnetic field and/or temperature [16][17][18][19][20][21][22][23][24][25][26][27][28]. Besides the most intensively analysed spin-1/2 case, a few theoretical works were focused on the entanglement of the mixed-spin Heisenberg systems . The particular interest in this field of study has been motivated by the pioneering work by Wang et al. [26], which demonstrated that a higher difference between dissimilar spin constituents can slightly enhance the thermal entanglement at higher temperatures due to the respective shift of its threshold temperature. Unfortunately, the enlargement of a threshold temperature is simultaneously accompanied by the reduction of the strength of mutual quantum correlations. In order to minimize the reduction of the degree of entanglement, many of subsequent studies were concentrated on an extended mixedspin Heisenberg chain involving the Dzyaloshinskii-Moriya interaction (DMI) [34][35][36][37], the effect of nonuniform magnetic field [37][38][39][40][41][42][43], long-range interaction [44] and uniaxial single-ion anisotropy [45,46], respectively. It was verified for the mixed spin-(1/2,1), spin-(1/2,3/2) and spin-(1/2,5/2) Heisenberg chains that the inhomogeneity of the external magnetic field can be suitable tuning parameter for enhancing the thermal entanglement in a high-temperature region. On the other hand, the DMI can enhance the low-temperature entanglement of the antiferromagnetic mixed spin-(1/2,3/2) Heisenberg dimer [36], nevertheless the enhancement of entanglement in the mixed spin-(1/2,1) Heisenberg dimer is possible only for the ferromagnetic exchange coupling [34]. Based on the results obtained for the mixed spin-(1/2,1) Heisenberg dimer [45,46], the uniaxial single-ion anisotropy seems to be another relevant driving force for an enhancement of entanglement in presence of thermal fluctuations. However, the comprehensive analysis of other mixed-spin Heisenberg dimers with higher spins S > 1 is still absent. This fact motivated us to study the mixed spin-(1/2,S) Heisenberg dimers involving the uniaxial single-ion anisotropy and magnetic field simultaneously with the special goal to verify whether the interplay between the uniaxial single-ion anisotropy, magnetic field and the spin magnitude S can enhance thermal entanglement at high enough temperatures. It should be emphasized that the uniaxial single-ion anisotropy may be relevant for several heterodinuclear complexes as for instance in magnetic compounds [49] such as the [MnCu(pbaOH)(H 2 O) 3 ]·nH 2 O (pbaOH = 2-hydroxy-1,3-propylenebis (oxamato)) [50], [Ni(dpt)(H 2 O) 3 Cu(pba)]·H 2 O (pba = 1,3-propylenebis(oxamato) and dpt = bis-(3aminopropyl)amine) [51] or [NiCu(pba)(D 2 O) 3 ]·2D 2 O [52].
The paper is organized as follows. The investigated mixed spin-(1/2,S) Heisenberg dimer will be defined in Section 2 together with a few details of the calculation procedure used in our rigorous study. The most interesting results concerned with the quantum and thermal entanglement under the influence of increasing spin magnitude S will be discussed in Section 3. Besides the effect of an applied external magnetic field on the thermal entanglement and the threshold temperature of the mixed spin-(1/2,S) Heisenberg dimers will be also discussed for an arbitrary spin-S. Finally, some concluding remarks are given in Section 4 and a few details of analytical derivations are presented in Appendices A-C.

Model and Method
Let us consider the mixed spin-(1/2,S) Heisenberg dimers with an uniaxial singleion anisotropy under the influence of the external magnetic field defined through the following Hamiltonian: In above, the symbolsμ α andŜ α (α = x, y, z) correspond to spatial components of spin-1/2 and spin-S (S ≥ 1) operators, J is the XXZ exchange interaction with an exchange anisotropy ∆, D is the uniaxial single-ion anisotropy acting on the spin-S magnetic ion only. Finally, the model under the investigation accounts for the effect of external magnetic field B applied along the z-direction, g denotes the gyromagnetic Landé g-factor and µ B is the Bohr magneton.
In order to study the quantum and thermal entanglement of the mixed spin-(1/2,S) Heisenberg dimers we will employ the concept of Peres-Horodecki [53,54], according to which negative eigenvalue of a partially transposed density matrix is a necessary condition for the onset of entanglement. To quantify the strength of quantum and thermal entanglement one may therefore utilize the quantity known as negativity [55]. The negativity of the mixed state given by the density matrix ρ is by definition the sum of all negative eigenvalues λ i of partially transposed density matrixρ T 1/2 It is worthwhile to remark that the negativity of the maximally entangled state is equal to one-half (N = 1/2) for the mixed spin-(1/2,S) Heisenberg dimers, whereas the negativity completely vanishes (N = 0) in the separable states without the bipartite entanglement.
Before calculating the negativity it is necessary to derive the eigenvalues and eigenvectors of the Hamiltonian (1), which can be easily calculated in the standard orthonormal basis |µ z , S z constructed from all available eigenvectors of z-components of two constituent spins with eigenvalues µ z = ±1/2 and S z = −S, −S+1, . . . , S−1, S. For this purpose, let us introduce first the notation for raising and lowering ladder operatorsŜ ± =Ŝ x ±iŜ y and µ ± =μ x ±iμ y , which allows us to rewrite the Hamiltonian (1) into the more convenient form As a result, one immediately realizes that the z-component of the spin operatorsŜ z (μ z ) trivially act on the original basis statesŜ z |S z = S z |S z andμ z |µ z = µ z |µ z , whereas the raising and lowering ladder operatorsŜ + (μ + ) andŜ − (μ − ) shift the relevant quantum spin number by unitŷ Subsequently, the application of the Hamiltonian (3) on the basis state | ± 1/2, S z leads to the identityĤ where S z = −S, −S+1, . . . , S−1, S and h = Bgµ B . The non-zero matrix elements define the block diagonal structure of the Hamiltonian, which consists of two one-by-one blocks and 2S two-by-two blocks characterized by a specific value of the z-component of the total spin S z t = S z +µ z running from −S−1/2 to S+1/2. Consequently, one can easily derive the respective eigenvalues and eigenvectors. The extremal values of S z t = ±(S+1/2) define two one-by-one blocks, whose element is identical with its eigenvalue and the respective standard basis state designates the corresponding eigenvector Based on the knowledge of a complete energy spectrum of eigenvalues ε S,±(S+1/2) , ε ∓ S,S z t and corresponding eigenvectors | ± (S+1/2) , |(S z t ) ∓ (S z t = −S+1/2, . . . , S−1/2), one is able to construct the relevant density operatorρ according to the formulâ where β = 1/(k B T), T is an absolute temperature, k B is a Boltzmann's constant and Z is the partition function The density matrixρ representing a matrix representation of the density operator (14) again has in the standard basis the same block diagonal form (classified according to the S z t value) involving two one-by-one blocks with the extremal values of the total spin S z t = ±(S+1/2) and 2S two-by-two blocks with S z t = −S+1/2, . . . , S −1/2. All non-zero elements of the density matrix can be commonly expressed through the following general formulas For a completeness, the readers can find the explicit form of the density-matrix elements for a few selected mixed spin-(1/2,S) Heisenberg dimers (S = 1, 3/2, 2, 5/2) in the Appendix A.
In order to calculate the density matrixρ T 1/2 partially transposed with respect to the spin-1/2 subsystem, it is sufficient to replace the bra and ket state vectors referred to the spin-1/2 subsystem. It is clear, that the diagonal elements derived from Equations (16) and (17) remain unchanged, whereas off-diagonal ones are moved to other positions. The partial transposition does not conserve the total spin momentum S z t , but it conserves the staggered spin momentum S z tm = S z −µ z running from −S−1/2 to S+1/2. Subsequently, the non-zero elements of the partially transposed density matrixρ T 1/2 expressed in term of S z tm have the form Note that the density matrixρ T 1/2 is a block diagonal with a maximal block's size of 2 × 2 achieving a specific S z tm value. Two one-by-one blocks with extremal S z tm = ±(S+1/2) involve a single element which directly determines two positive eigenvalues λ ∓(S+1/2) of a partially transposed density matrix Other 2S two-by-two blocks determined by basis state vectors with S z tm = −S + 1/2, . . . , S−1/2 immediately result to the remaining couple of eigenvalues Due to length of explicit form of the Equation (25), the readers can find them in Appendix B. At the same time, the complete list of partially transposed density matricesρ T 1/2 for the mixed spin-(1/2,S) Heisenberg dimers with specific spin values S = 1, 3/2, 2, 5/2 is given in Appendix C. Analysing Equation (25) in detail one identifies that only eigenvalues λ − S z tm can be negative, and hence, the respective bipartite entanglement is in accordance to the definition (2) determined by the formula

Results and Discussion
In order to minimize the huge parametric space, all further discussions will be limited to the physically most interesting case with an isotropic exchange interaction defined through the parameter ∆ = 1. For simplicity, the gyromagnetic factor g of both types of magnetic ions is set equal to two (g = 2). It is worthwhile to note that other choice of the unequal g factor has only the quantitative, but not qualitative, impact on all obtained observations.

Quantum Negativity
The behaviour of the quantum negativity of the mixed spin-(1/2,S) Heisenberg dimers with S = 1, 3/2, 2, 5/2, 3, 7/2 is illustrated in Figure 1 in the D/J − µ B B/J plane by considering the antiferromagnetic exchange coupling J > 0. The density plots of quantum negativity simultaneously illustrate stability regions of all relevant ground states |(S z t ) − and |S+1/2 .
It is worthwhile to remark that the negativity at zero magnetic field was comprehensively analysed in our preceding paper [28] and thus, the case µ B B/J = 0 will be just marginally explored in our subsequent discussions. It has been found that the negativity at µ B B/J = 0 exhibits qualitatively different behaviour for integer and half-odd-integer spin-S constituents, if the uniaxial single-ion anisotropy D/J > 0 of easy-plane type is taken into account. It was surprisingly detected that the enhancement of a degree of entan-glement for the mixed-spin Heisenberg dimers involving an integer spin S emerges as a consequence of interplay between the increasing spin magnitude S and uniaxial single-ion anisotropy D/J > 1/2. Nevertheless, the highest negativity N = ( √ 5−1)/4 reached for the mixed-spin Heisenberg dimers with integer spin S at the specific value D/J = 1/2 is significantly smaller than the maximal negativity N = 1/2 detected for the mixed-spin Heisenberg dimers with an arbitrary half-odd-integer spin S. The external magnetic field reduces the degeneracy of the mixed spin-(1/2,S) Heisenberg dimers due to the Zeeman's splitting of energy levels and the quantum negativity of an arbitrary mixed spin-(1/2,S) Heisenberg dimer in a state with the total spin S z t can be expressed through the general formula It is evident from Equation (27) and Figure 1 that the quantum negativity of each |(S z t ) − ground state is fully independent of the external magnetic field and its magnitude decreases towards to the completely separable (N = 0) ferromagnetic state |1/2+S with S z t = S+1/2. However, the increasing magnetic field is responsible for existence of discontinuous changes of the quantum negativity at all field-driven magnetic phase transitions. It follows from Figure 1 that maximal bipartite entanglement is reached for the specific value of the uniaxial single-ion anisotropy D/J = 1/2 (similar as in the µ B B/J = 0 case), at which an arbitrary mixed spin-(1/2,S) Heisenberg dimer shows, in agreement with Equation (27), the highest quantum negativity N = 1/2 until the sufficiently high magnetic field reorients both spins into its direction.
The behaviour of the quantum negativity in the regime of easy-axis uniaxial single-ion anisotropy D/J ≤ 0 confirms previously reported findings [30,33,34] that the increasing spin magnitude S enlarges the stability of entangled state with respect to the magnetic field. Nevertheless, the degree of respective bipartite entanglement is gradually reduced. Contrary to this, the increasing spin magnitudes S induces the enhancement of a quantum negativity for an arbitrary |(S z t ) − (S z t ≤ S−1/2) ground state of the mixed spin-(1/2,S) Heisenberg dimers assuming easy-plane single-ion anisotropy D/J > 0. In contradiction to the zero-field case, the enhancement of a negativity is observed even for 0 < D/J < 1/2 as a consequence of reduction of ground-state degeneracy in respective parametric space. This is a very important observation from the application perspective, because variation of a magnetic ion in the mixed spin-(1/2,S) Heisenberg dimers, offers a relative simple alternative how to enhance the bipartite entanglement. The origin of qualitatively different behaviour of the negativity below and above D/J = 0 can be explained through the respective variation of the total spin value S z t . In the |(S−1/2) − ground state emergent in easy-axis regime D/J < 0, the total spin S z t = S−1/2 is gradually enhanced with increasing spin size S, but the difference 4S(S+1)−(2S z t −1)(2S z t +1) entering into the Equation (27) Considering the fixed value of the total spin S z t and the easy- (27) enlarges with an increasing spin size S, which means that the denominator in rewritten form of Equation (27) increases and the respective quantum negativity is thus naturally enhanced. In the special case of D/J = 0 the negativity is an inverse function of spin magnitude S and thus the quantum entanglement reduces upon strengthening of the spin size S. For the ferromagnetic exchange coupling J < 0 (Figure 2), the quantum entanglement can be achieved only for the mixed-spin Heisenberg dimers with an easy-plane single-ion anisotropy D/|J| > 0 due to possible preference of various ferrimagnetic or antiferromagnetic ground states |(S z t ) + (S z t ≤ S−1/2). Since the previously derived relation remains in force, the quantum negativity of each |(S z t ) + ground state with the total spin S z t ≤ S−1/2 always increases as the spin magnitudes increase. It should be pointed out that the degree of bipartite entanglement for an arbitrary |(S z t ) + ground state is significantly smaller in comparison to its antiferromagnetic counterpart, which makes the ferromagnetic quantum mixed-spin Heisenberg dimers less attractive for practical utilizations. For a completeness, it should be emphasized, that the invariant point at D/|J| = 1/2 is completely absent in the ferromagnetic case J < 0 and the highest negativity (excluding the ground state |(0) + ) can be found in a proximity of the isotropic point D/|J| = 0. Strictly at D/|J| = 0 and µ B B/J = 0 all |(S z t ) + (|S z t | ≤ S−1/2) ground states and | ± (S+1/2) ones are degenerate and the respective negativity follows the simple relation

The Thermal Negativity
In order to analyse the thermal behaviour of the negativity we have chosen the specific sets of model parameters under the influence of magnetic field being consistent with |(S−1) − ground state (Figure 3), |(0) ± or the |(1/2) ± ground state ( Figure 4) and finally with the |(1) ± or |(3/2) ± ones ( Figure 5).
Focusing on Figure 3 we can generalize our previous zero-temperature conclusion [28], which states that the increasing spin magnitude S reduces the quantum as well as lowtemperature thermal negativity. On the other hand, the increasing spin magnitude S enlarges the threshold temperature, which subsequently allows us to detect a subtle enhancement of the thermal entanglement at larger temperatures if the spin size S increases. The obtained results are in a perfect quantitative agreement with previous observations for the Heisenberg dimers without the uniaxial single-ion anisotropy D/J [26,29,33,34].  The most significant finding follows from Figures 4 and 5. It is evident from these figures that the increasing spin magnitude S can enhance not only the threshold temperature, but it can also enhance the degree of thermal entanglement in contrast to previous knowledge. It should be emphasized that the above statement holds provided that the ensemble of integer or half-odd-integer spins S is taken into account separately. At low magnetic fields, the thermal negativity of Heisenberg dimers consisting of both half-oddinteger-spins always saturate in the maximal value N = 1/2 in the asymptotic limit of absolute zero temperature, whereas the maximum of negativity of Heisenberg dimers composed of integer and one-half spin converges to the value N = 1 2 4S(S+1) (1−2D/J) 2 +4S(S+1) . Hence, one can immediately conclude that the negativity of the Heisenberg dimers with integer spins S can reach the maximum value N = 1/2 just for the special case D/J = 1/2.  Furthermore, in Figure 6, we present, the behaviour of the thermal negativity under the changes of magnetic field. The same model parameters have been used as in the above analysis. It should be emphasized that absence of Zeeman's term leads to the ground-state degeneracy in the zero-temperature limit and thus, the zero-field negativity is always smaller than that in an arbitrary small but non-zero magnetic field. This fact is visualized through the symbols on the y-axis determining the respective asymptotic values of the negativity in zero-field limit.
In agreement with general expectation the increasing temperature reduces the bipartite entanglement with a significant drop of the negativity emergent in a proximity of all magnetic-field-induced phase transitions associated with crossing of energy levels. Around each level-crossing field the interplay between thermal and quantum fluctuations is the most pronounced and the negativity shows a marked local maximum located between two neighbouring level-crossing fields. In a consequence of that, the mixed spin-(1/2,S) Heisenberg dimers exhibit very specific oscillating changes of the negativity at low and moderate temperatures (Figure 6b). It is worthwhile to remark that such oscillating behaviour is possible only for the spin-(1/2,S) Heisenberg dimers with higher spin magnitude S > 1, because existence of at least two level-crossing fields has to be guaranteed. In addition it turns out that the distance between two local maxima can be tuned through the uniaxial single-ion anisotropy D/J, which unambiguously determines the stability region (field range) of a given magnetic ground state (see Figures 1 and 2). Finally, let us turn our attention to the dependence of the threshold temperature on the spin magnitude S as well as the external magnetic field. The results presented in Figure 7 in the form of the threshold temperature versus magnetic field plot for a few selected spin sizes S confirm previous conclusions [26,29,33,34] that the threshold temperature gradually enlarges with an enhancement of the spin magnitude S. It is also quite evident from Figure 7b that the antiferromagnetic mixed spin-(1/2,S) Heisenberg dimers are more persistent against rising temperature and magnetic field than their ferromagnetic counterparts (the negative values of k B T c /J in Figure 7b correspond to the mixed-spin Heisenberg dimers with the ferromagnetic exchange coupling J < 0). Another remarkable finding is that all investigated mixed spin-(1/2, S) Heisenberg dimers exhibit a striking reentrant behaviour of the threshold temperature regardless of the character and size of the exchange coupling and the uniaxial single-ion anisotropy. The origin of the unconventional reentrant phenomenon could be explained by incapability of the magnetic field to suppress the thermally induced population of the entangled excited states, which is reflected in a thermally stimulated rise of the negativity. It appears worthwhile to remark that existence of magnetic-field-driven phase transitions gives rise to a stepwise dependence of the respective threshold temperature. The threshold temperature slightly decreases at each level-crossing field for the antiferromagnetic Heisenberg dimers (J > 0), while an opposite effect is observed for the ferromagnetic Heisenberg dimers (J < 0), see Figure 7b.

Conclusions
In the present paper, we have exactly examined the effect of the spin magnitude S, magnetic field and uniaxial single-ion anisotropy on the quantum and thermal entanglement of the mixed spin-(1/2,S) Heisenberg dimers. In particular, it has been verified that the concurrent interplay of the uniaxial single-ion anisotropy and the external magnetic field basically influences bipartite entanglement of the mixed-spin Heisenberg dimers. To quantify the degree of bipartite entanglement we have derived the exact analytical expression for the negativity in terms of Peres-Horodecki criterion [53,54] followed by the mathematical formulation due to Vidal and Werner [55]. In the present study, we have provided first an exhaustive analysis of all possible ground states of the mixed spin-(1/2,S) Heisenberg dimers as a necessarily prerequisite for further entanglement analysis. Two different scenarios were observed for both the antiferromagnetic (J > 0) and ferromagnetic (J < 0) coupling constants depending on the character of an uniaxial single-ion anisotropy D/J. For the easy-axis single-ion anisotropy D/J < 0 the mixed-spin Heisenberg dimers exhibit either one or none magnetic-field-driven phase transition, whereas the increasing magnetic field generates S+(2S mod 2)/2 consecutive field-driven phase transitions between the ground states with the total spin |S z t | ≤ S−1/2 for the easy-plane single-ion anisotropy D/J > 0.
As a direct consequence of different effect of easy-axis and easy-plane uniaxial singleion anisotropy one detects very different influence of increasing spin magnitude S on the bipartite entanglement. In the case of easy-axis uniaxial single-ion anisotropy, an increasing spin S always suppresses the degree of quantum entanglement as dictated by the formula Equation (28). In contrast to this, the increasing spin magnitude S for the easy-plane single-ion anisotropy D/J > 0 leads to the coincidence of regions with a fixed number of the total spin S z t ≤ S−1/2, and in accordance with the Formula (29), the enhancement of a quantum entanglement can be observed. Interestingly, two specific conditions with maximal entanglement invariant on external stimuli have been identified for: (i) antiferromagnetic mixed spin-(1/2,S) Heisenberg dimers with an arbitrary spin-S magnitude for the particular value of the uniaxial single-ion anisotropy D/J = 1/2, or (ii) the antiferromagnetic ground state |(0) − exclusively existing only in the mixed spin-(1/2,S) Heisenberg dimers with both half-odd-integer spin constituents.
The comprehensive analysis of a thermal entanglement above a unique ground state |(S−1/2) − confirms previously reported findings that the increasing spin magnitude S enlarges the threshold temperature, but unfortunately reduces the degree of thermal entanglement. However, a completely novel and unexpected behaviour has been observed for the easy-plane single-ion anisotropy D/J > 0, where the increasing spin S may simultaneously enlarge the threshold temperature as well as the degree of thermal entanglement if one compares solely integer or half-odd-integer spin-S cases. It has been evidenced that except the particular case with D/J = 1/2 the mixed spin-(1/2,S) Heisenberg dimers with both half-odd-integer spins generally achieve higher degree of entanglement, which makes them more attractive for a practical utilization. Last but not least, a fascinating oscillating changes of the negativity were observed upon the variation of the external magnetic field at low and moderate temperatures. The unconventional oscillating behaviour of the negativity originates from existence of consecutive field-driven phase transitions emergent at level-crossing fields, at which the negativity rapidly falls down before it produces significant local maxima localized in between two level-crossing fields. It is noteworthy that the local maxima of the negativity of the mixed spin-(1/2,S) Heisenberg dimers with integer and half-odd-integer spin S are not equal due to existence of different ground states in the parameter region of the easy-plane single-ion anisotropy D/J > 0. Of course, the envelope of such oscillations is gradually suppressed upon strengthening of the magnetic field until the fully polarized state is reached. The density matrix of the mixed spin-(1/2,1) Heisenberg dimer can be recast into the following block-diagonal form The individual elements of the density matrix are defined as The partition function Z of the mixed spin-(1/2,1) Heisenberg dimer reads as follows (A10) Finally, the probability coefficients c ± S,S z t and the respective energies ε ± S,S z t (S z t = −3/2, −1/2, 1/2, 3/2) have the explicit forms

. Density Matrix of the Mixed Spin-(1/2,3/2) Heisenberg Dimer
The density matrix of the mixed spin-(1/2,3/2) Heisenberg dimer can be recast into the following block-diagonal form The individual elements of the density matrix are defined as The partition function Z of the mixed spin-(1/2,3/2) Heisenberg dimer reads as follows Finally, the probability coefficients c ± S,S z t and the respective energies ε ± S,S z t (S z t = −2, −1, 0, 1, 2) have the explicit forms The density matrix of the mixed spin-(1/2,2) Heisenberg dimer can be recast into the following block-diagonal form The individual elements of the density matrix are defined as (A46) Finally, the probability coefficients c ± S,S z t and the respective energies ε ± S,S z t (S z t = −5/2, −3/2, −1/2, 1/2, 3/2, 5/2) have the explicit forms  The individual elements of the density matrix are defined as     The non-zero elements of ρ T 1/2 i,j are explicitly defined in Equations (A32)-(A45).