Quantum Magnetism in Fe2Cu2 Polymeric Branched Chains: Insights from Exactly Solved Ising-Heisenberg Model

The spin-1/2 Ising-Heisenberg branched chain, inspired by the magnetic structure of three isostructural polymeric coordination compounds [(Tp)2Fe2(CN)6X(bdmap)Cu2(H2O)].H2O to be further denoted as Fe2Cu2 (Tp = tris(pyrazolyl)hydroborate, bdmapH = 1,3-bis(dimethylamino)-2-propanol, HX = acetic acid, propionic acid or trifluoroacetic acid), is rigorously studied using the transfer-matrix method. The overall ground-state phase diagram reveals three distinct phases: a quantum antiferromagnetic phase, a quantum ferrimagnetic phase and a classical ferromagnetic phase. In the zero-temperature magnetization curve, two quantum ground states are manifested as intermediate plateaus at zero and half of the saturation magnetization, while the magnetization reaches its saturated value within the classical ferromagnetic phase. The bipartite entanglement between nearest-neighbor Heisenberg spins is more pronounced in the quantum ferrimagnetic phase compared to the quantum antiferromagnetic phase due to a fully polarized nature of the Ising spins. A reasonable agreement between theoretical predictions for the spin-1/2 Ising-Heisenberg branched chain and experimental data measured for a temperature dependence of the magnetic susceptibility and a low-temperature magnetization curve suggests strong antiferromagnetic coupling between nearest-neighbor CuII-CuII magnetic ions and moderately strong ferromagnetic coupling between nearest-neighbor CuII-FeIII magnetic ions in the polymeric compounds Fe2Cu2. A thermal entanglement between nearest-neighbor CuII-CuII magnetic ions persists up to a relatively high threshold temperature T = 224 K and undergoes a transient magnetic-field-driven strengthening.


INTRODUCTION
Although one-dimensional (1D) quantum spin models may seem relatively simple at first glance, achieving an exact solution for these lattice-statistical models is often accompanied with formidable mathematical difficulties stemming from noncommuting algebra [1].As a matter of fact, there exists only a few paradigmatic examples of exactly solved 1D quantum spin models such as the quantum spin-1/2 Heisenberg chain solved by the Bethe-ansatz method [2,3,4,5] or the quantum spin-1/2 XY model solved using the Jordan-Wigner fermionization approach [6,7].It should be stressed that the exploration of 1D quantum spin models is of particular importance, because several low-dimensional magnetic materials may have negligible exchange interactions in the remaining two spatial dimensions [8].It has been demonstrated that several 1D spin chains may display intriguing quantum features such quantized magnetization plateaus [9] and quantum entanglement [10].
Figure 2. A schematic illustration of the spin-1/2 Ising-Heisenberg branched chain with the isotropic exchange interaction J between nearest-neighbor Heisenberg spins and the anisotropic exchange interaction J 1 between the nearest-neighbor Heisenberg and Ising spins, respectively.The Heisenberg spins capture lattice positions of Cu 2+ magnetic ions (yellow circles), while the Ising spins capture lattice positions of Fe 3+ magnetic ions (blue circles).
ions (S = 1/2) described by the notion of quantum Heisenberg spins, which is in contrast to highly magnetically anisotropic spin-1/2 Fe 3+ magnetic ions (σ = 1/2) described by the notion of Ising spins.Consequently, the exchange coupling between nearest-neighbor Cu 2+ -Cu 2+ magnetic ions will be assumed to be of the Heisenberg type, whereas the exchange coupling between nearest-neighbor Fe 3+ -Cu 2+ magnetic ions will be treated as the Ising-type exchange interaction (see Fig. 2).
The structure of this paper is as follows.In Section 2, we introduce the spin-1/2 Ising-Heisenberg branched chain along with a brief outline of the method employed for its calculation.Section 3 is dedicated to discussing the theoretical findings regarding the ground-state phase diagram, magnetization curves, entropy, and concurrence.In Section 4, we compare the exact theoretical results for the spin-1/2 Ising-Heisenberg branched chain with available experimental data for the magnetization and susceptibility of Siv ý et al.
Quantum Magnetism in Fe 2 Cu 2 Polymeric Branched Chains the polymeric compound Fe 2 Cu 2 together with the respective theoretical prediction for the concurrence.Finally, Section 5 provides a summary of the most important scientific findings obtained in this study.

MODEL AND METHOD
The spin-1/2 Ising-Heisenberg branched chain schematically illustrated in Fig. 2 can be defined through the following Hamiltonian: where σz j,i and Ŝα j,i (α = x, y, z) label the Ising and Heisenberg spins ascribed to Fe 3+ and Cu 2+ magnetic ions, respectively.The coupling constant J > 0 denotes the antiferromagnetic exchange interaction between the nearest-neighbor Heisenberg spin pairs capturing lattice positions of the dimeric Cu 2+ −Cu 2+ units in the polymeric compounds Fe 2 Cu 2 , while the coupling constant J 1 > 0 (J 1 < 0) stands for the antiferromagnetic (ferromagnetic) exchange interaction between the nearest-neighbor Ising and Heisenberg spins capturing lattice positions of Fe 3+ and Cu 2+ magnetic ions, respectively.Furthermore, the Zeeman's terms h H = g Cu µ B B and h I = g Fe µ B B account for the magnetostatic energy of magnetic moments ascribed to the Heisenberg and Ising spins in an external magnetic field B, g Cu and g Fe are the respective Landé g-factors and µ B is Bohr magneton.The parameter ∆ measures a degree of the exchange anisotropy in the XXZ Heisenberg interaction and N denotes the total number of unit cells.For simplicity, the periodic boundary conditions σ 1,N +1 = σ 1,1 are assumed.
It is convenient to rewrite the total Hamiltonian (1) as a sum of the cell Hamiltonians Ĥ = N i=1 Ĥi with the following definition of the cell Hamiltonian Ĥi : Owing to the fact that the cell Hamiltonians (2) commute with each other [ Ĥi , Ĥj ] = 0, the partition function of the spin-1/2 Ising-Heisenberg branched chain can be factorized into the following product: where β = 1/(k B T ), k B is the Boltzmann's constant, T is the absolute temperature, Tr i stands for a trace over degrees of freedom of the Heisenberg dimer S 1,i − S 2,i from the i-th unit cell, σ 2,i denotes a summation over two states of the Ising spin from the branching of the i-th unit cell, and the summation {σ 1,i } runs over all possible configurations of the Ising spins from the main backbone of the branched chain.The partition function (3) can be rewritten into a more suitable form for further calculations: where the reduced cell Hamiltonian Ĥ′ i is defined as: The reduced cell Hamiltonian Ĥ′ i given by Eq. ( 5) can be alternatively viewed as the Hamiltonian of a spin-1/2 Heisenberg dimer in nonuniform effective fields h j,i involving the interaction with the adjacent Ising spins: To perform a trace Tr i over degrees of freedom of the spin-1/2 Heisenberg dimer from the i-th unit cell in the relevant expression for the partition function ( 4), it is convenient to proceed to a matrix representation of the reduced cell Hamiltonian (5) in the standard basis: The state vectors |↑ j,i and |↓ j,i are eigenvectors of the z-component of the Heisenberg spin S z j,i = 1/2 and −1/2, respectively.After solving the relevant eigenvalue problem ϕ k,i Ĥ′ i ϕ l,i − E j,i I = 0 one gets four eigenvalues of the reduced cell Hamiltonian (5): with the corresponding eigenvectors: where the probability amplitudes A ∓ i are given by: With the use of the eigenvalues (8) one can further rewrite the partition function (4) into the following simplified form: where the effective Boltzmann's factor T(σ 1,i , σ 1,i+1 ) is defined as follows: The effective Boltzmann's weight (12) depends only on degrees of freedom of two Ising spins σ 1,i and σ 1,i+1 and it can be alternatively viewed as the transfer matrix: where T i (i = 1, 2, 3) denote three different elements of the transfer matrix obtained from Eq. ( 12) by considering all four available states of two Ising spins σ 1,i and σ 1,i+1 : Frontiers Siv ý et al.

Quantum Magnetism in Fe 2 Cu 2 Polymeric Branched Chains
In the spirit of the transfer-matrix method, the final expression for the partition function can be finally expressed in terms of the transfer-matrix eigenvalues: The transfer matrix (13) can be readily diagonalized in order to obtain the eigenvalues λ ± defined through the transfer-matrix elements (14): In the thermodynamic limit N → ∞, the Gibbs free energy of the spin-1/2 Ising-Heisenberg branched chain per unit cell is given by the largest transfer-matrix eigenvalue: From the Gibbs free energy (17) one can subsequently derive the local magnetization of the Heisenberg and Ising spins as well as the total magnetization (per spin) of the spin-1/2 Ising-Heisenberg branched chain: From the exact result (17) for the Gibbs free energy one can eventually derive other quantities such as magnetic susceptibility, entropy and specific heat using the standard relations of thermodynamics.Besides, one may also investigate a strength of the bipartite entanglement inside of the Heisenberg dimers using the quantity concurrence [38], which can be calculated from the pair correlation functions and the local magnetization [39,40]: Here, Ŝz 2,i denote two spatial components of the pair correlation function between the nearest-neighbor Heisenberg spins and m H is the respective local magnetization given by Eq. (18).The pair correlation functions can be straightforwardly derived from the Gibbs free energy (17) according to the relations: but their explicit form is too cumbersome to write it down here explicitly.

THEORETICAL RESULTS
Let us proceed to a discussion of the most interesting results obtained for the spin-1/2 Ising-Heisenberg branched chain by assuming the isotropic antiferromagnetic Heisenberg coupling constant J > 0, while keeping the anisotropy parameter fixed at the constant value of ∆ = 1 for illustrative purposes.Throughout this section, the coupling constant J > 0 will serve as an energy unit when defining a relative strength of the magnetic field h/J, the dimensionless temperature k B T /J and the interaction ratio J 1 /J.
To reduce the number of free parameters, the Landé g-factors of the Heisenberg and Ising spins assigned to Cu 2+ and Fe 3+ magnetic ions were set equal to each other, thereby ensuring that the 'local' magnetic fields exert an equal strength h ≡ h H = h I .
By comparing energies of all available eigenstates of the spin-1/2 Ising-Heisenberg branched chain one finds just three different ground states schematically shown in Fig. 3 and mathematically given by the eigenvectors: While the factorized form of the third ground state ( 23) is consistent with the classical nature of the fully polarized ferromagnetic state |III , the quantum antiferromagnetic phase |I, I ′ and the quantum ferrimagnetic phase |II, II ′ described by Eqs. ( 21) and ( 22) exhibit bipartite quantum entanglement among the Heisenberg spin pairs S 1,i -S 2,i unambiguously characterized through the probability amplitudes a ± and b ± : From this perspective, it would be insightful to compare the strength of the bipartite quantum entanglement emerging within the quantum antiferromagnetic (21) and ferrimagnetic (22) phases using the concept of the concurrence (19), which acquires within these two ground states the following values: It is obvious from Eq. ( 25) that the bipartite quantum entanglement of the Heisenberg spin pairs is consistently stronger within the quantum ferrimagnetic phase |II, II ′ compared to the quantum For the sake of completeness, let us also quote eigenenergies (per unit cell) corresponding to the quantum antiferromagnetic phase |I, the quantum ferrimagnetic phase |II, the classical ferromagnetic phase |III It should be stressed that the ground states |I and |I ′ have the same energy and they just differ by a relative orientation of the singlet-like state of the Heisenberg spin pairs with respect to their neighboring Ising spins [see Fig. 3(a)-(b)] depending on whether the Ising coupling constant is antiferromagnetic J 1 > 0 or ferromagnetic J 1 < 0. The same statement holds true also for the quantum ferrimagnetic phases |II and |II ′ schematically drawn in Fig. 3(c)-(d), respectively.A direct comparison of all three ground-state energies ( 26)-( 28) affords exact expressions for the relevant ground-state phase boundaries: 1. phase boundary |I /|II and |I ′ /|II ′ : 2. phase boundary |I ′ /|III : 3. phase boundary |II /|III and |II ′ /|III : 4. phase boundary |II /|II ′ : Siv ý et al.

Quantum Magnetism in Fe 2 Cu 2 Polymeric Branched Chains
Using these exact expressions, we have plotted in Fig. 4(a) the overall ground-state phase diagram of the spin-1/2 Ising-Heisenberg branched chain in the J 1 /J − h/J plane.The parameter regions denoted in Fig. 4(a) as I(I'), II(II') and III correspond to the quantum antiferromagnetic (QAF), the quantum ferrimagnetic (QFI) and the classical ferromagnetic (CF) phase, respectively.Note furthermore that the coexistence lines (29)- (31) correspond to discontinuous magnetic-field-driven phase transitions.All three phases QAF, QFI and CF coexist together in Fig. 4(a) at a triple coexistence point emerging for the specific value of the interaction ratio J 1 /J ≈ −0.97.Below this value, the QFI ground state |II ′ is absent and only the two ground states |I ′ and |III appear.
To verify this observation, we present magnetic-field dependencies of the local and total magnetization in Fig. 4(b) for the interaction ratio J 1 /J = −1.5 and three distinct temperature values.Under this condition, the local magnetization m H and m I defined as a single-site magnetization of the Heisenberg and Ising spins are zero until the saturation field h s /J ≈ 0.3 is reached.This consequently leads to a zero magnetization plateau of the total magnetization m, which is consistent with the antiferromagnetic alignment of the Ising spins and the singlet-like ordering of the Heisenberg spin pairs realized within the QAF ground state |I ′ [see Fig. 3(b)].The observed single magnetization jump can thus be attributed to a discontinuous magnetic-field-driven phase transition from the QAF ground state |I ′ to the CF ground state |III .It should be also noted that the true magnetization jump exists at zero temperature only, which however changes to a steep but continuous rise upon increasing of temperature.
Next, the local and total magnetization of the spin-1/2 Ising-Heisenberg branched chain are plotted against the magnetic field in Fig. 4(c) for another value of the interaction ratio J 1 /J = −0.5 and three distinct temperatures.The spin-1/2 Ising-Heisenberg branched chain undergoes in this particular case two magnetic-field-driven phase transitions: initially from the QAF ground state |I ′ towards the QFI ground state |II ′ subsequently followed by a transition from the QFI ground state |II ′ towards the CF ground state |III .In agreement with this statement, the local magnetization of the Ising spins m I does not contribute to the total magnetization m only at sufficiently low magnetic fields h/J 0.11 stabilizing the QAF ground state |I ′ and then it undergoes sudden rise towards the saturated value acquired within the QFI ground state |II ′ and the CF ground state |III .On the contrary, the local magnetization of the Heisenberg spins m H does not contribute to the total magnetization m unless the magnetic field exceeds its saturation value h s /J ≈ 0.64.
Finally, we illustrate in Fig. 4(d) the local and total magnetization of the spin-1/2 Ising-Heisenberg branched chain for a fixed value of the interaction ratio J 1 /J = 1.5.The behavior observed here closely resembles that of the previous cases, with two intermediate magnetization plateaus and two magnetic-fieldinduced transitions: first from the QAF phase |I towards the QFI phase |II , followed by a transition from the QFI phase |II towards the CF phase |III .The primary quantitative difference lies in a somewhat wider intermediate one-half magnetization plateau corresponding to the QFI phase |II , which generally broadens with an increase in the interaction ratio J 1 /J.
An independent validation of the aforementioned magnetization scenario is provided through a comprehensive analysis of the magnetic entropy, as depicted in Fig. 5 in the form of a density plot in the magnetic field versus temperature plane (h/J -k B T /J).The isoentropy lines shown in Fig. 5 bring further insight into the adiabatic demagnetization process of the spin-1/2 Ising-Heisenberg branched chain.It is evident from Fig. 5 that a sharp decrease (or increase) in temperature is observed upon decreasing (or increasing) the magnetic field above (below) each magnetic-field-driven phase transition.Consistent with this observation, the isoentropy lines tend towards absolute zero temperature at a single magnetic field for the specific case J 1 /J = −1.5, owing to the absence of the quantum ferrimagnetic phase |II ′ at Siv ý et al.  moderate magnetic fields [refer to Fig. 5(a)].Conversely, for the other two specific cases, J 1 /J = −0.5 and 1.5, the isoentropy lines approach absolute zero temperature at two distinct magnetic fields due to the presence of the quantum ferrimagnetic phase |II ′ at moderate magnetic fields.The temperature dependence of the initial magnetic susceptibility times temperature product is illustrated in Fig. 6 for various values of the interaction ratio J 1 /J.It is noteworthy that the spin-1/2 Ising-Heisenberg branched chain exhibits the quantum antiferromagnetic (QAF) ground state in the zero-field limit for any value of the interaction ratio J 1 /J [refer to Fig. 4(a)].This aligns with the vanishing behavior of the susceptibility times temperature product as temperature approaches zero.It is apparent from Fig. 6 that the product of susceptibility and temperature gradually tends to zero as temperature decreases for the interaction ratio J 1 /J = 1.5.For the interaction ratio J 1 /J = −0.25, the product χT contrarily initially declines, followed by maintaining a nearly constant value below the temperature k B T /J 0.25, before rapidly decreasing to zero as absolute zero temperature is approached.The notable temperature dependence with a nearly constant value of the susceptibility times temperature product observed at moderate temperatures is attributed to the ferromagnetic nature of the Ising coupling constant J 1 < 0, which has opposite tendency of increasing of the product χT upon decreasing of temperature.Indeed, the reinforcement of the ferromagnetic coupling constant J 1 leads to an increase in the susceptibility times temperature product as temperature decreases, as evidenced in Fig. 6 for specific values of the interaction ratio J 1 /J = −1.5 and −2.3, where the product χT initially reaches a rounded maximum before ultimately approaching zero as temperature decreases.

Quantum Magnetism in
To better understand the quantum nature of the spin-1/2 Ising-Heisenberg branched chain, we conduct a detailed analysis of the concurrence (19), which quantifies the strength of bipartite entanglement within the Heisenberg spin pairs.In Fig. 7(a), we illustrate the magnetic-field dependence of the concurrence for a fixed interaction ratio J 1 /J = −1.5 and three distinct temperatures.At sufficiently low magnetic fields, the concurrence is non-zero due to the presence of the quantum antiferromagnetic (QAF) ground state |I ′ , characterized by singlet-like states of the Heisenberg dimers.This state breaks down at the magnetic-field-induced phase transition h c /J ≈ 0.3, leading to the emergence of the fully polarized classical ferromagnetic (CF) ground state |III .Fig. 7(a) demonstrates that increasing temperature generally smears out the stepwise magnetic-field dependence of the concurrence.The corresponding thermal variations of the concurrence are depicted in Fig. 7(d) for three selected magnetic field values.At zero magnetic field, the concurrence exhibits a standard monotonous decline with increasing temperature until it completely vanishes at a threshold temperature.However, a more intricate temperature dependence is observed for the magnetic field h/J = 0.3, around which the magnetic-fielddriven phase transition occurs at zero temperature.Under this condition, the concurrence undergoes an abrupt decrease to a local minimum at the lowest temperatures, followed by a gradual rise to a local maximum at moderate temperatures, before eventually tending to zero with further temperature increase.Moreover, the concurrence initially starts from zero due to the CF ground state |III at the magnetic field h/J = 1, which exceeds the saturation value.However, the concurrence exhibits a striking doublereentrant behavior, developing first at a lower threshold temperature k B T /J ≈ 0.2, then reaching a rounded maximum at a moderate temperature, before returning to zero at an upper threshold temperature k B T /J ≈ 0.85.This double-reentrant behavior is attributed to the thermal population of the QAF phase |I ′ .
Next, we depict in Fig. 7(b) the concurrence as a function of the magnetic field for a fixed interaction ratio J 1 /J = −0.5 and three distinct temperatures.In contrast to the previous case, the spin-1/2 Ising-Heisenberg branched chain encompasses at zero temperature all three possible ground states: QAF, QFI, and CF.Consequently, the concurrence initially attains a relatively large value (C ≈ 0.8) due to the QAF ground state |I ′ realized at sufficiently low magnetic fields (h/J < 0.11).This is followed by an almost maximal value (C ≈ 0.97) attributed to the QFI ground state |II ′ , before the concurrence abruptly drops to zero near the saturation field (h s /J ≈ 0.64) at very low temperature (k B T /J = 0.01), as observed in Fig. 7(b).Furthermore, Fig. 7(b) reveals that the concurrence may become non-zero above the saturation field (h s /J ≈ 0.64) at slightly higher temperatures (k B T /J = 0.1 and 0.3) due to thermal excitations from the fully polarized CF ground state |III to the quantum states |II ′ and |I ′ .To gain a comprehensive understanding, Fig. 7(e) illustrates the concurrence against temperature for the same value of the interaction ratio J 1 /J = −0.5 and three different magnetic field values.At zero magnetic field, the concurrence starts from a finite value (C ≈ 0.8) corresponding to the QAF ground state, then steadily rises to reach a rounded maximum before gradually tending to zero.The prominent rounded peak of the concurrence can be attributed to thermal excitations from the less entangled QAF phase to the QFI phase with stronger entanglement.If the magnetic field h/J = 0.6 favors the QFI ground state, the concurrence sharply drops from its almost maximal value within a narrow temperature range (k B T /J ∈ (0; 0.08)) and then decreases steadily before showing a more abrupt thermally-induced decline.At a higher magnetic field (h/J = 1) exceeding the saturation value, the concurrence exhibits a double-reentrant behavior quite analogous to the previously described case for J 1 /J = −1.5.
Last but not least, we have depicted in Fig. 7(c) and (f) the magnetic-field and temperature dependencies of the concurrence for the fixed value of the interaction ratio J 1 /J = 1.5.It is evident that the results shown in Fig. 7(c) and (f) are qualitatively very similar to those displayed in Fig. 7(b) and (e) for the other specific case J 1 /J = −0.5.The main quantitative difference lies in the strength of bipartite entanglement, i.e., the magnitude of concurrence, which generally weakens with an increase in the absolute value of the interaction ratio |J 1 |/J.Notably, the concurrence attains the same zero-temperature asymptotic value within the QAF ground states |I and |I ′ for two specific values of the interaction ratio J 1 /J = 1.5 and −1.5, respectively. .Blue squares represent experimental data reported in Ref. [37], while green curves are the respective theoretical fits obtained with the help of the spin-1/2 Ising-Heisenberg branched chain model.

THEORY VS. EXPERIMENT
In this section, we will take advantage of the spin-1/2 Ising-Heisenberg branched-chain model described by the Hamiltonian (1) to theoretically interpret available magnetic data experimentally measured for the heterobimetallic coordination polymer [(Tp) 2 Fe 2 (CN) 6 (CH 3 COO)(bdmap)Cu 2 (H 2 O)] to be further abbreviated as Fe 2 Cu 2 , which serves as its experimental realization [37].The experimental data for the temperature dependence of the magnetic susceptibility times temperature product χ M T recorded at a magnetic field of B = 0.2 T [37] are confronted in Fig. 8(a) with the best theoretical fit obtained using the spin-1/2 Ising-Heisenberg branched-chain model.The least-squares fitting procedure yielded a rather strong antiferromagnetic coupling constant J/k B = 247.8K between nearest-neighbor Cu 2+ -Cu 2+ magnetic ions, a moderately strong ferromagnetic coupling constant J 1 /k B = −64.2K between nearest-neighbor Cu 2+ -Fe 3+ magnetic ions, and Landé g-factors g Cu = 2.16 and g Fe = 2.23 for Cu 2+ and Fe 3+ magnetic ions, respectively.A moderate value of the interaction ratio J 1 /J ≈ −0.26 contributes to a relatively rapid decline in the susceptibility times temperature product as temperature decreases in a high-temperature regime T 50 K.This is successively followed by a less steep decline observed in a range of moderate temperatures 20 K T 50 K before it finally rapidly falls to zero with further temperature reduction.Note furthermore that the steep reduction of the susceptibility times temperature product observed in the low-temperature range is consistent with the quantum antiferromagnetic (QAF) ground state |I ′ .
Utilizing the determined set of fitting interaction parameters, we have also conducted theoretical modeling of the low-temperature magnetization curve of the polymeric compound Fe 2 Cu 2 , as reported in Ref. [37], for sufficiently low temperatures T = 1.8K up to a magnetic field of B = 5 T. The magnetization data recorded for the coordination polymer Fe 2 Cu 2 are confronted in Fig. 8(b) with the corresponding theoretical fit based on the spin-1/2 Ising-Heisenberg branched chain model, calculated up to an ultrahigh magnetic field of 150 T surpassing the saturation field.The spin-1/2 Ising-Heisenberg branched chain model provides a relatively satisfactory fit to the experimental magnetization data.Moreover, an extrapolation to higher magnetic fields predicts the presence of a sizable intermediate onehalf magnetization plateau in the magnetic field range between B ≈ 10 T and 140 T.This intermediate magnetization plateau reflects the emergence of the quantum ferrimagnetic phase |II ′ accompanied by a stronger bipartite entanglement of the Heisenberg spin pairs.It is worthwhile to mention that a small deviation between the theoretical and experimental data can be attributed to the Ising approximation of a highly (but not infinitely) anisotropic exchange coupling between Cu 2+ and Fe 3+ magnetic ions.
To shed light on the strength of bipartite entanglement within the Cu 2+ -Cu 2+ dimeric units of the polymeric chain Fe 2 Cu 2 , we utilized the determined fitting parameters of the corresponding spin-1/2 Ising-Heisenberg branched chain to calculate the concurrence (19).Temperature dependencies of the concurrence (19) calculated for the Heisenberg spin pairs of the spin-1/2 Ising-Heisenberg branched chain are illustrated in Fig. 9(a) for a few selected values of the magnetic field B. At zero magnetic field, the concurrence exhibits a striking temperature-induced rise almost to its maximum possible value, after which it gradually declines with increasing temperature until it completely vanishes at the threshold temperature T ≈ 224 K.A small magnetic field can be considered as another driving force that may reinforce the strength of bipartite entanglement.Indeed, it turns out that the concurrence approaches almost its maximum value in the zero-temperature asymptotic limit when a small magnetic field is applied.If the magnetic field exceeds the saturation field, a double reentrant behavior of the concurrence is observed.
Last but not least, a density plot of the concurrence (19) calculated for the Heisenberg spin pairs of the spin-1/2 Ising-Heisenberg branched chain is presented in Fig. 9(b) in the magnetic field versus temperature plane.This density plot can be interpreted as a phase diagram, delineating the entangled parameter region from the disentangled one.Remarkably, the density plot in Fig. 9(b) reveals that a relatively high value of the concurrence (C = 0.1) can still be detected even up to relatively high temperatures of around T ≈ 188 K irrespective of the accessible magnetic-field strength for experimental testing (B < 10 T).This suggests the persistence of significant bipartite entanglement in the system over a wide temperature range, which could be of interest for experimental investigations.

CONCLUSION
In this paper, we have investigated the ground-state phase diagram, magnetization process, magnetocaloric effect, and bipartite entanglement of the spin-1/2 Ising-Heisenberg branched chain with two different coupling constants using the transfer-matrix method.Our analysis reveals that the spin-1/2 Ising-Heisenberg branched chain exhibits five distinct ground states: two quantum antiferromagnetic phases |I and |I ′ , two quantum ferrimagnetic phases |II and |II ′ , and the classical ferromagnetic phase |III .These ground states are separated by discontinuous magnetic-field-driven phase transitions, whereby the ground states |I, I ′ and |II, II ′ are responsible for zero and one-half magnetization plateaus, respectively.Furthermore, our exact calculation of the concurrence has confirmed the quantum character of the ground states |I, I ′ and |II, II ′ .We find that the strength of bipartite entanglement within the quantum ferrimagnetic phases |II, II ′ is stronger compared to that within the quantum antiferromagnetic phases |I, I ′ .This insight contributes to a deeper understanding of the entanglement properties of the system and its magnetic behavior.
It has been demonstrated that the spin-1/2 Ising-Heisenberg branched chain with two different coupling constants provides a satisfactory description of the magnetic behavior observed in the heterobimetallic polymeric chain Fe 2 Cu 2 .By fitting available experimental data for the temperature dependence of the susceptibility times temperature product and the low-temperature magnetization curve, we have determined a rather strong coupling constant J/k B = 247.8K between the Cu 2+ -Cu 2+ magnetic ions and a moderately strong coupling constant J 1 /k B = −64.2K between the Cu 2+ -Fe 3+ magnetic ions with subtle differences in the Landé g-factors g Cu = 2.16 and g F e = 2.23.Our analysis reveals that the quantum antiferromagnetic phase is realized at sufficiently low magnetic fields (B 10 T), while the quantum ferrimagnetic phase, characterized by a sizable intermediate one-half magnetization plateau, can be detected in the magnetic-field range between B ≈ 10 T and 140 T.The substantial coupling constant between the Cu 2+ -Cu 2+ dimeric units further implies a robust bipartite entanglement between the nearest-neighbor Cu 2+ -Cu 2+ magnetic ions.This entanglement persists up to a relatively high threshold temperature (T ≈ 224 K) and undergoes a transient magnetic-field-driven strengthening.

Figure 3 .
Figure 3.A schematic illustration of typical configurations realized within all possible ground states of the spin-1/2 Ising-Heisenberg branched chain: (a)-(b) the quantum antiferromagnetic phases |I and |I ′ ; (c)-(d) the quantum ferrimagnetic phases |II and |II ′ ; (e) the classical ferromagnetic phase |III .The pair of arrows situated at the Heisenberg spins symbolizes a singlet-like quantum superposition, whereby the size of arrows reflects different magnitudes of the respective probability amplitudes (24).

Figure 4 .
Figure 4. (a) The ground-state phase diagram for the spin-1/2 Ising-Heisenberg branched chain in the J 1 /J − h/J plane by assuming the isotropic value of the parameter ∆ = 1.The regions I(I'), II(II') and III correspond to a quantum antiferromagnetic, a quantum ferrimagnetic and a classical ferromagnetic phase, respectively; (b)-(d) Magnetic-field dependencies of the local and total magnetization given by Eq. (18) for the spin-1/2 Ising-Heisenberg branched chain at three different temperatures and three selected values of the interaction ratio: (b) J 1 /J = −1.5;(c) J 1 /J = −0.5;(d) J 1 /J = 1.5.

Fe 2 2 5 Figure 5 .
Figure 5.A density plot of the reduced entropy S/Nk B of the spin-1/2 Ising-Heisenberg branched chain in the magnetic field versus temperature plane h/J − k B T /J for three selected values of the interaction ratio: (a) J 1 /J = −1.5,(b) J 1 /J = −0.5, and (c) J 1 /J = 1.5.

Figure 6 .
Figure 6.Temperature dependence of the initial magnetic susceptibility times the temperature χk B T /J of the spin-1/2 Ising-Heisenberg branched chain for a few different values of the interaction ratio J 1 /J.

Siv ý et al. Quantum Magnetism in Fe 2
Cu 2 Polymeric Branched Chains

Figure 8 .
Figure 8. Temperature dependence of the molar magnetic susceptibility χ M times temperature T data at the magnetic field B = 0.2 T (a) and the isothermal magnetization curve of the polymeric coordination compound [(Tp) 2 Fe 2 (CN) 6 (CH 3 COO)(bdmap)Cu 2 (H 2 O)] recorded at temperature T = 1.8K (b).Blue squares represent experimental data reported in Ref.[37], while green curves are the respective theoretical fits obtained with the help of the spin-1/2 Ising-Heisenberg branched chain model.

Figure 9 .
Figure 9. (a) Temperatures dependencies of the concurrence (19) calculated for the Cu 2+ -Cu 2+ dimeric units using the spin-1/2 Ising-Heisenberg branched chain for several values of the magnetic field B (the inset shows a detail from the low-field region); (b) A density plot of the concurrence (19) calculated for the Cu 2+ -Cu 2+ dimeric units of the spin-1/2 Ising-Heisenberg branched chain in the magnetic field vs. temperature plane.Displayed lines correspond to isocontour lines of the concurrence.All calculations are presented for the set of fitting parameters J/k B = 247.8K, J 1 /k B = −64.2K, g Cu = 2.16 and g F e = 2.23 obtained for the polymeric coordination compound [(Tp) 2 Fe 2 (CN) 6 (CH 3 COO)(bdmap)Cu 2 (H 2 O)].

Siv ý et al. Quantum Magnetism in Fe 2
Cu 2 Polymeric Branched Chains

Siv ý et al. Quantum Magnetism in Fe 2 Cu 2 Polymeric Branched Chains antiferromagnetic
phase |I, I ′ for a given set of the interaction ratio J 1 /J and the exchange anisotropy ∆, because a fully polarized nature of the enclosing Ising spins tends to bring a singlet-like state closer to a perfect singlet-dimer state as one of the Bell states.