Most quantum entanglement investigations are focused on two qubits or some finite (small) chain structure, since the infinite chain structure is a considerably cumbersome task. Therefore, the quantum entanglement properties involving an infinite chain structure is quite important, not only because the mathematical calculation is cumbersome but also because real materials are well represented by an infinite chain. Thus, in this paper we consider an entangled diamond chain with Ising and anisotropic Heisenberg (Ising-XXZ) coupling. Two interstitial particles are coupled through Heisenberg coupling or simply two-qubit Heisenberg, which could be responsible for the emergence of entanglement. These two-qubit Heisenberg operators are interacted with two nodal Ising spins. An infinite diamond chain is organized by interstitial-interstitial and nodal-interstitial (dimer-monomer) site couplings. We are able to get the thermal average of the two-qubit operator, called the reduced two-qubit density operator. Since these density operators are spatially separated, we could obtain the concurrence (entanglement) directly in the thermodynamic limit. The thermal entanglement (concurrence) is constructed for different values of the anisotropic Heisenberg parameter, magnetic field and temperature. We also observed the threshold temperature via the parameter of anisotropy, Heisenberg and Ising interaction, external magnetic field, and temperature.

Quantum entanglement is one of the most fascinating types of correlation that can be shared only among quantum systems. The Heisenberg chain is one of the simplest quantum chains which exhibits a reach entanglement feature, due to the Heisenberg interaction is quantum coupling in the spin system. The two particles were coupled trough XYZ coupling or simply called as two-qubit XYZ spin, which are the responsible for the emergence of thermal entanglement. These two-qubit operators are bonded to two nodal Ising spins, and this process is repeated infinitely resulting in a diamond chain structure. We will discuss two-qubit thermal entanglement effect on Ising-XYZ diamond chain structure. The concurrence could be obtained straightforwardly in terms of two-qubit density operator elements, using this result, we study the thermal entanglement, as well as the threshold temperature where entangled state vanishes. The present model displays a quite unusual concurrence behavior, such as, the boundary of two entangled regions becomes a disentangled region, this is intrinsically related to the XY-anisotropy in the Heisenberg coupling. Despite a similar property had been found for only two-qubit, here we show in the case of a diamond chain structure, which reasonably represents real materials. S i S z J J 0 J Figure 1: (Color Online) Schematic representation of Ising-XYZ chain on diamond structure, σa,i and σ i,b are Heisenberg spins, while Si corresponds to Ising spins.

The entanglement quantum properties of a spin-1/2 Ising-Heisenberg model on a symmetrical diamond chain were analyzed. Due to the separable nature of the Ising-type exchange interactions between neighboring Heisenberg dimers, calculation of the entanglement can be performed exactly for each individual dimer. Pairwise thermal entanglement was studied in terms of the isotropic Ising-Heisenberg model and analytical expressions for the concurrence (as a measure of bipartite entanglement) were obtained. The effects of external magnetic field H and next-nearest neighbor interaction J(m) between nodal Ising sites were considered. The ground state structure and entanglement properties of the system were studied in a wide range of coupling constant values. Various regimes with different values of ground state entanglement were revealed, depending on the relation between competing interaction strengths. Finally, some novel effects, such as the two-peak behavior of concurrence versus temperature and coexistence of phases with different values of magnetic entanglement, were observed.

The Ising chain consisting of the antiferromagneticaly coupled ferromagnetic trimer is considered in the external magnetic field. In the framework of the transfer-matrix formalism the thermodynamics of the system is described. The magnetization per site (m) is obtained as the explicit function of the external magnetic field (H). The corresponding plots of m(H) are drawn. Two qualitatively different regions of the values of coupling constants are established: weak antiferromagnetic coupling (J A < 3J F ) and the strong antiferromagnetic coupling (J A ≥ 3J F ). For the latter case the magnetization curve with plateau at m/m sat = 1/3 is obtained. It is proven that the plateau is caused by the stability of spatially modulated spin structure 3111 . The values of magnetic field determining the width of the plateau are obtained in the limit of zero temperature. * E-mail address: vohanyan@www.physdep.r.am Fig. 4 The magnetization curves for κ = 6 at different temperatures: (a) T = 2.5J F ; (b) T = 0.8J F ; (c) T = 0.1J F ; (d) T = 0.001J F .