Next Article in Journal
Fabrication and Investigation of the Microwave Absorption of Nonwovens Modified by Carbon Nanotubes and Graphene Flakes
Previous Article in Journal
Accessing Properties of Molecular Compounds Involved in Cellular Metabolic Processes with Electron Paramagnetic Resonance, Raman Spectroscopy, and Differential Scanning Calorimetry
Previous Article in Special Issue
Distribution of Bipartite and Tripartite Entanglement within a Spin-1/2 Heisenberg Star in a Magnetic Field
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Phase Diagram and Quantum Entanglement Properties of a Pentamer S = 1/2 Heisenberg Spin Cluster

Department of Solid State Physics, Faculty of Physics and Applied Informatics, University of Łódź, ul. Pomorska 149/153, PL90-236 Łódź, Poland
Molecules 2023, 28(17), 6418; https://doi.org/10.3390/molecules28176418
Submission received: 24 July 2023 / Revised: 21 August 2023 / Accepted: 29 August 2023 / Published: 3 September 2023
(This article belongs to the Special Issue Molecular Magnetism: Modern Trends and Future Perspectives)

Abstract

:
Cluster molecular magnets prove their potential for applications in quantum technologies, encouraging studies of quantum entanglement in spin systems. In the paper we discuss quantum entanglement properties of pentamer cluster composed of spins S = 1 / 2 forming a tetrahedron with additional spin in its center, with geometry reproducing the smallest nonplanar graph. We model the system with isotropic Heisenberg Hamiltonian including external magnetic field and use exact diagonalization approach to explore the ground-state phase diagram and thermodynamic properties within canonical ensemble formalism. We focus the interest on two-spin entanglement quantified by Wootters concurrence. For ground state, we find two states with total cluster spin equal to 3/2 exhibiting entanglement, occurring preferably for antiferromagnetic interactions. For finite temperatures, we predict the presence of magnetic-field-induced entanglement as well as temperature-induced entanglement.

1. Introduction

Molecular magnetism constitutes a flourishing field of study within the theoretical and experimental condensed matter physics and chemistry [1,2], combining the rich fundamental physics of quantum low-dimensional magnets [3,4] with prospective applications in crucial areas such as magnetocaloric cooling [5,6]. However, another factor stimulating research in this scope is the high usefulness of molecular magnets for information storage and processing [7,8,9], both on classical [10] and quantum level [11,12,13,14,15]. The latter revokes the phenomenon of quantum entanglement [16,17,18,19,20].
The interest in quantum entanglement brings the localized spin systems of cluster geometry to the attention of the theoreticians [21,22,23,24]. In such context, a variety of geometries can be mentioned, like dimers [25,26,27,28,29,30,31,32,33,34,35,36], trimers [37,38,39,40,41,42,43] or four-site structures [44,45,46,47,48,49,50,51,52]. Attention is paid to design of cluster structures offering advantageous properties due to interplay between geometry and spin-spin interactions [53,54]. Moreover, the interest can be generalized to entanglement in nanochains of various lengths housing itinerant electrons, with emphasis on both spin and charge degree of freedom [55,56,57,58]. In addition, several experimental works on localized spin systems aimed at capturing the entanglement in zero-dimensional structures have been reported [59,60,61,62,63], supplemented with works focused on systems with higher dimensionality [64,65,66].
Among various shapes of clusters, a pentamer geometry can be highlighted. Such clusters attracted so far some attention of theorists, but only in the context unrelated to quantum entanglement, to mention the neutron scattering simulation [67,68], magnetocaloric properties [69] or general thermodynamics [70]. The experimental studies of usual magnetic properties of these structures were conducted as well [71,72,73]. Within the group, a pentamer formed by inserting a central ion into a tetrahedron can be singled out as geometrically interesting due to high symmetry of its structure and the fact that it constitutes a smallest non-planar graph. Its triangle-based geometry presages considerable magnetic frustration linked to the possible presence of degenerate ground states in the case of competing interactions. Such structures can be found in metal-organic compounds in which Co ions contribute localized spins [74], with a single Co ion possessing octahedral coordination environment and four other ions with tetrahedral environment. Let us also notice that pentanuclear structural units of similar type, tetrahedral in shape and containing Co or Ni ions, were found to build diamond-like networks [75]. Similar pentamers, called diamondoids, were also discussed in Fe-based compounds [76]. What is even more important, an alike structure was found in Cu ion-based compound, where Cu ions contributed localized spins S = 1 / 2 [77]. This particular sort of pentameric structures was subject of some model theoretical studies limited to determining the eigenenergies [78]. Therefore, the quantum entanglement properties of these structures remain unexplored.
Motivated by the synthesability of pentamer spin clusters described above in the field of molecular magnetism and by their non-trivial, highly symmetric geometry, we present a theoretical study focused on characterization of quantum entanglement in such structure. In order to capture the maximized quantum effects, we select cluster composed of spins S = 1 / 2 . We explore the full phase diagram of the system in question and identify its possible ground states, considering then ground-state and finite-temperature properties of two-spin entanglement. In the following part of the paper we describe the theoretical model and its thermodynamic solution. Next, we discuss the analytic and numerical results concerning the entanglement in the context of system phase diagram.

2. Theoretical Model and Its Thermodynamics

The system of interest is a spin pentamer cluster, composed of spins S = 1 / 2 forming a regular tetrahedron with an additional spin placed inside it. The schematic view of the considered pentamer cluster is shown in Figure 1a, with two exchange integrals, J 1 and J 2 , denoted with solid and dashed lines, respectively. The interactions between the spins forming an external tetrahedron (labelled with 2, 3, 4 and 5) are quantified by J 2 , whereas the coupling of all the tetrahedron ions with a central spin (labelled with 1) is equal to J 1 . An alternative structure with equivalent interactions is shown in Figure 1b, composed of a planar tetramer with additional spin. Let us mention, however, that in the structure shown in Figure 1a the distances between the tetramer spins coupled with J 2 are equal, which is not the case for the structure shown in Figure 1b.
An interesting feature of the studied cluster is that it forms a smallest nonplanar Kuratowski graph, being a complete graph of 5 vertices K 5 [79], each vertex containing a spin and each spin interacting with each one. As a consequence, there are only two kinds of spin pairs - those connected by J 1 or J 2 exchange integrals - and all of them are of nearest-neighbour type. Such non-planarity of the graph representing the spin cluster was also mentioned in other studies on metal-organic pentamers [80,81].
We model the system with isotropic Heisenberg Hamiltonian:
H ^ = J 1 S ^ 1 · S ^ 2 + S ^ 3 + S ^ 4 + S ^ 5 J 2 S ^ 2 · S ^ 3 + S ^ 2 · S ^ 4 + S ^ 2 · S ^ 5 + S ^ 3 · S ^ 4 + S ^ 3 · S ^ 5 + S ^ 4 · S ^ 5 H S ^ 1 z + S ^ 2 z + S ^ 3 z + S ^ 4 z + S ^ 5 z .
In the above Hamiltonian, the symbol S ^ i S ^ i x , S ^ i y , S ^ i z is a vector of quantum operators of spin S = 1 / 2 , localized at site labelled by i, moreover, S ^ i · S ^ j = S ^ i x S ^ j x + S ^ i y S ^ j y + S ^ i z S ^ j z . Operators for individual spin projections α = x , y , z are given by S ^ i α = σ ^ i α / 2 , where σ ^ α denotes appropriate Pauli matrix. The external magnetic field is introduced by a Zeeman term parametrized with energy H.
It should be noticed that the Hamiltonian (1) commutes with an operator of the square of total spin S ^ T 2 = S ^ 1 + S ^ 2 + S ^ 3 + S ^ 4 + S ^ 5 2 as well as the projection of the total spin on z axis S ^ T , z = S ^ 1 z + S ^ 2 z + S ^ 3 z + S ^ 4 z + S ^ 5 z . Moreover, it commutes also with an operator of the square of the total spin of the tetrahedron (excluding the central spin labelled with 1 in Figure 1) which we denote with s ^ 2 = S ^ 2 + S ^ 3 + S ^ 4 + S ^ 5 2 . Therefore, the Hamiltonian eigenstates can be labelled with the appropriate quantum numbers S T = 1 / 2 , 3 / 2 , 5 / 2 , S T , z = S T , . . . , S T and s = 0 , 1 , 2 , respectively. As a consequence, the eigenenergies of the Hamiltonian (1) can be written in the following form (utilizing the vector coupling method [82]):
E S T , S T , z , s = 1 2 J 1 S T ( S T + 1 ) 1 2 J 2 J 1 s ( s + 1 ) + 3 8 J 1 + 3 2 J 2 H S T , z .
However, it must be stated that some of the states labelled with three quantum numbers S T , S T , z , s are further degenerate.
The symbolic or numerically exact diagonalization of the Hamiltonian matrix provides a set of eigenvectors ψ n and eigenenergies E n (where n = 1 , , N s and N s = 2 5 = 32 ). The complete set of Hamiltonian eigenvalues and eigenstates is given in the Appendix A in Table A1 and Table A2. We label the states as ψ S T , S T , z , s and use additional superscript if the state is degenerate, which is quite often the case for the studied system.
In the description of system states in our work, the single-spin states with spin up (down) are marked with ( ), respectively. For multispin states, the order of the arrows corresponds to increasing order of site labels from 1 to 5. This convention is used in Table A1. In order to facilitate interpretation of multispin states, we decompose them also into linear combinations of products of state of a single spin (labelled with 1) and a states of two spin dimers (one of them including spins 2 and 3 and another one composed of spins 4 and 5—see Figure 1). These dimer states are, if convenient, expressed by so called Bell states [83] possessing the following form:
ϕ ± = ± 2 ψ ± = ± 2 .
Out of Bell states, the state ψ is a singlet state with dimer spin equal to 0, whereas the remaining three states constitute members of a triplet and correspond to dimer spin of 1. Let us notice that the tetrahedron composed of spins 2, 3, 4 and 5 can be decomposed into two dimers in numerous manners, so that our choice is non-unique and arbitrary. However, all the dimer selections are fully equivalent due to high symmetry of interactions within the tetrahedron. The described convention is used to present the eigenstates in Table A2.
The output of Hamitonian diagonalization can be directly used for analysis of the ground state of the system (at temperature T = 0 ) in the full Hamiltonian parameter space spanned by J 1 , J 2 and H. The thermodynamic description of the system is constructed following the rules of canonical ensemble of statistical physics. Let us put emphasis on the fact that the external magnetic field is directly included in the Hamiltonian (1), so that we deal with a field ensemble [84,85] and the thermodynamic average of the Hamiltonian yields enthalpy, not internal energy.
For T > 0 , the statistical sum can be calculated as:
Z = Tr e H ^ / k B T = n = 1 N s e E n / k B T ,
where k B denotes Boltzmann constant. The quantum thermal state of the whole cluster is then given by the density matrix:
ρ ^ = e H ^ / k B T / Z .
From the viewpoint of the purpose of the study, the quantum states of spin pairs are of particular interest. The density matrix describing the state of a pair of spins labelled by i j can be obtained by taking a partial trace from ρ ^ over remaining spins, i.e.,
ρ ^ i j = Tr k l m ρ ^ ,
with k , l , m i , j and i , j , k , l , m 1 , 2 , 3 , 4 , 5 .
In order to quantify quantum entanglement for pairs of spins i j , Wootters concurrence C is recalled [86,87]. For separable states C = 0 , whereas 0 < C 1 corresponds to entangled state. This quantity is calculated for a general mixed state on the basis of the following formula:
C = max 0 , λ 1 λ 2 λ 3 λ 4 ,
where λ 1 λ 2 λ 3 λ 4 are the eigenvalues of the matrix ρ ^ i j ρ ^ i j . The auxiliary matrix ρ ^ i j corresponds to spin flip of the state and is given by:
ρ ^ i j = σ ^ y σ ^ y ρ ^ i j * σ ^ y σ ^ y ,
where σ ^ y is appropriate Pauli matrix.
For the studied system, due to its high symmetry, only two inequivalent values of concurrence can be found, corresponding either to type 1 spin pairs coupled with an exchange integral J 1 and containing the central spin (termed C ( 1 ) ) or to type 2 spin pairs coupled with an exchange integral J 2 containing only the spins from the tetrahedron (termed C ( 2 ) ).

3. Results and Discussion

The results of analytic or numerically exact calculations based on the formalism sketched above will be discussed in this section. In particular, the ground-state phase diagram will be a starting point for analysis of ground-state entanglement and further investigation of the influence of finite temperature on the system properties.

3.1. Ground-State Phase Diagram

Let us commence the study of the cluster from the discussion of the ground-state magnetic phase diagram, for T 0 . Such diagrams are shown in Figure 2, separately for the case of ferromagnetic interaction J 1 > 0 (Figure 2a) and antiferromagnetic interaction J 1 < 0 (Figure 2b). In all further discussion, the dimensionless quantities will be used, with | J 1 | as an unit of energy. Therefore, the cases of J 1 > 0 and J 1 < 0 will be discussed separately and J 1 0 is assumed (i.e., the limit of spin tetrahedron without coupling to the central spin is not discussed, as corresponding to two completely uncoupled systems).
It should be emphasized that we deal with a finite spin cluster, therefore, no actual phase transitions are expected and no symmetry breaking phenomenon occurs, contrary to infinite systems in thermodynamic limit. Instead, the diagrams termed as phase diagrams illustrate various ground states of the system and these ground states are not phases in rigorous thermodynamic meaning. However, the term ’phase diagram’ is commonly accepted in this generalized meaning.
Let us notice that due to possible degeneracy of the Hamiltonian eigenstates for the cluster, the ground states for some parameter ranges can be expected to be mixed quantum states (i.e., probabilistic mixtures of degenerate eigenstates with equal probabilities) rather than pure states. Therefore, if we deal with d states of the same energy, ψ ( q ) , q = 1 , . . . , d , the ground state when T 0 is
ρ ^ = 1 d q = 1 d ψ ( q ) ψ ( q ) .
If the case of H > 0 away from phase boundary is considered, all the states involved must have the same values of all quantum numbers S T , S T , z and s. The remark about mixed ground state also applies to the case of the phase boundary (i.e., line along which two ground states have the same energy) or a triple/quadruple point (where three/four ground states indicate the same energy) when the degeneracy appears necessarily, even if the neighbouring ground states are pure ones. In both such cases the state of the system is a mixed quantum state being a probabilistic mixture of all the involved Hamiltonian eigenstates sharing the same energy along specific line or at point (this time with different values of the quantum numbers). In the absence of the external magnetic field, for H = 0 , the additional degeneracy with respect to the value of S T , z emerges [see Equation (2)] and all the states with the same values of S T and s share the same energy.
If the ground state for a given range of parameters is non-degenerate, it is marked in Figure 2 for brevity with a ket ψ S T , s (labelled with the quantum numbers S T for total spin of the cluster and s for spin of tetrahedron). In the presence of the external magnetic field H > 0 (even arbitrarily weak), the ground state always contains the Hamiltonian eigenstates for which S T , z = S T . Therefore, it is sufficient to label the ground states in Figure 2 with just S T and s. On the contrary, in the case of a mixed state, it is described by a density matrix ρ ^ S T , s (labelled with the same numbers, as the states degenerate in energy share the same S T and s).
In both panels of Figure 2, the boundaries at which the energies corresponding to two different ground states are equal are marked with lines. The equations describing these lines (phase boundaries) are given in the inset. Moreover, the triple or quadruple points, at which the energies of three or four different ground states are equal are marked with bold circles. Let us remark that we discuss a finite spin cluster, so that the existence of quadruple point does not violate Gibbs phase rule, as we do not deal with phase transitions. Such points can be found in ground-state phase diagrams of spin clusters, see for example Refs. [70,88].

3.1.1. Ferromagnetic J 1 > 0

The ground-state phase diagram for ferromagnetic J 1 > 0 is shown in Figure 2a. In the whole range of magnetic fields for J 2 / | J 1 | > 1 / 4 and just for H > 1 2 J 1 2 J 2 otherwise, the system takes the saturated ferromagnetic pure state with S T = S T , z = 5 / 2 and s = 2 , i.e., ψ 5 / 2 , + 5 / 2 , 2 (see Table A1 and Table A2), having the energy of
E = J 1 3 2 J 2 5 2 H .
For 1 2 J 1 J 2 < H < 1 2 J 1 2 J 2 the ground state is a mixed state ρ ^ 3 / 2 , 1 , being a statistical mixture given by Equation (9) involving d = 3 states ψ 3 / 2 , + 3 / 2 , 1 with S T = S T , z = 3 / 2 and s = 1 , listed in Table A1 and Table A2, all sharing the same energy,
E = 1 2 J 1 + 1 2 J 2 3 2 H .
In this state, the central spin inside the tetrahedron (labelled with 1) takes the pure state .
For H < 1 2 J 1 J 2 (possible only for J 2 / | J 1 | < 1 / 2 ), the ground state is a mixed state ρ ^ 1 / 2 , 0 of the form given by Equation (9), composed of d = 2 states ψ 1 / 2 , + 1 / 2 , 0 with S T = S T , z = 1 / 2 and s = 0 , as listed in Table A1 and Table A2, having the energy equal to:
E = 3 2 J 2 1 2 H .
Also in this state, the central spin inside the tetrahedron (labelled with 1) takes the pure state .
Let us separately mention the form of the ground states exactly at H = 0 , as under this condition the system symmetry is higher and more degeneracies appear because the state energy is S T , z -independent.
If J 2 / | J 1 | > 1 / 4 , the quantum states ψ 5 / 2 , S T , z , 2 with S T = 5 / 2 and all possible S T , z = ± 5 / 2 , ± 3 / 2 , ± 1 / 2 share the same energy of E = J 1 3 2 J 2 and correspond to s = 2 . The ground state is a mixed state expressed by Equation (9) with d = 6 and the appropriate states are listed in Table A1 and Table A2.
If the coupling fulfils the condition 1 / 2 < J 2 / | J 1 | < 1 / 4 , the ground state is composed of d = 12 eigenstates ψ 3 / 2 , S T , z , 1 sharing the energy equal to E = 1 2 J 1 + 1 2 J 2 , having S T = 3 / 2 and S T , z = ± 1 / 2 , ± 3 / 2 as well as s = 1 , possessing the form given by Equation (9) with the relevant states from Table A1 and Table A2.
Some of the above states, corresponding to S T , z = + 3 / 2 , contribute to the ground state ρ ^ 3 / 2 , 1 for H > 0 .
Finally, if J 2 / | J 1 | < 1 / 2 , the mixed ground state has energy equal to E = 3 2 J 2 and it takes the form predicted by Equation (9) with d = 4 states ψ 1 / 2 , S T , z , 0 having S T = 1 / 2 , S T , z = ± 1 / 2 and s = 0 , given in Table A1 and Table A2 (note that the states for S T , z = + 1 / 2 contribute to the state ρ ^ 1 / 2 , 0 for H > 0 .).

3.1.2. Antiferromagnetic J 1 < 0

The ground-state phase diagram for antiferromagnetic J 1 < 0 is depicted in Figure 2b. If J 2 / | J 1 | < 1 , the same states as those found for J 1 > 0 are present in the diagram. If J 2 / | J 1 | > 1 , two more quantum states can emerge as ground states. In this range, if the magnetic field exceeds H / | J 1 | = 5 / 2 , the saturated ferromagnetic state occurs. However, below this field, but for H > 1 2 J 1 2 J 2 , a pure state with S T = S T , z = 3 / 2 and s = 2 , ψ 3 / 2 , + 3 / 2 , 2 , is a ground state of the system (see Table A1 and Table A2); it has energy equal to
E = 3 2 J 1 3 2 J 2 3 2 H .
In this state, the central spin inside a tetrahedron, labelled with 1, is in mixed state of the form 1 5 + 4 5 .
If H < 1 2 J 1 2 J 2 (only for J 2 / | J 1 | < 1 / 4 ), the ground state is mixed state ρ ^ 1 / 2 , 1 [see Equation (9)] composed of d = 3 states with S T = S T , z = 1 / 2 and s = 1 , listed in Table A1 and Table A2 as ψ 1 / 2 , + 1 / 2 , 1 , sharing the energy of
E = 1 2 J 1 + 1 2 J 2 1 2 H .
In this state, the spin labelled with 1 is in mixed state of the form 1 3 + 2 3 .
Again, for H = 0 the degeneracy of the Hamiltonian eigenvalues increases significantly and the ground states are discussed separately.
If J 2 / | J 1 | > 1 / 4 , the ground state is a mixed state [Equation (9)] containing d = 4 eigenstates ψ 3 / 2 , S T , z , 2 of energy E = 3 2 J 1 1 2 J 2 , with S T = 3 / 2 , S T , z = ± 3 / 2 , ± 1 / 2 and s = 2 , taking the form indicated in Table A1 and Table A2.
In the case of 1 < J 2 / | J 1 | < 1 / 4 , the ground state [Equation (9)] is formed out of d = 6 eigenstates ψ 1 / 2 , S T , z , 1 with energy E = J 1 + 1 2 J 2 and S T = 1 / 2 , S T , z = ± 1 / 2 and s = 1 , which are expressed in Table A1 and Table A2.
Finally, if J 2 / | J 1 | < 1 , the quantum state is a mixed state discussed already for ferromagnetic J 1 > 0 , possessing the energy of E = 3 2 J 2 (it is notable that the energy is independent on J 1 ) and taking the form predicted by Equation (9) with d = 4 states ψ 1 / 2 , S T , z , 0 having S T = 1 / 2 , S T , z = ± 1 / 2 and s = 0 , of the form expressed in Table A1 and Table A2.
After analysis of the ground-state phase diagram it can be concluded that relatively large area is filled with states showing degeneracy, due to high symmetry of the studied cluster.

3.2. Ground-State Entanglement

The system ground states can exhibit entanglement for the case of ferromagnetic J 1 (Figure 3a) and of antiferromagnetic J 1 (Figure 3b); both cases are analysed in the further part of the paper.
In the context of the ground state being a mixed state given by Equation (9) for the case of degeneracy, it should be mentioned that the creation of probabilistic mixture of entangled states may lead to entangled or to separable state. Therefore, even if the member states ψ ( q ) in Equation (9) are entangled, the state ρ ^ is not necessarily entangled. This is an important observation for the case of cluster systems with degeneracy of ground states. Some results regarding this phenomenon were, for example, discussed in Ref. [89], for the case of probabilistic mixture of two pure states. The results discussed in Ref. [89] show the overall complicacy of the relation of the entanglement of pure states and their probabilistic mixture.
The values of the Wootters concurrence for the ground state, at T = 0 , are shown in Figure 3 in a form of pairs C ( 1 ) , C ( 2 ) , indicating the concurrence of inequivalent spin pairs of both types, as discussed previously. Apart from the values valid for given phase, the values of concurrence are also assigned to the boundaries between the phases and to triple point, as shown by the arrows. Let us notice here that at the phase boundary or triple point all the Hamiltonian eigenstates corresponding to two or three phases take the same energy, so that the quantum state is a probabilistic mixture of all that states. Moreover, also the concurrence value at H = 0 for some range of interactions is indicated with an arrow, as the ground state degeneracy in the absence of the magnetic field is increased and the concurrence value may be different than that achieved for H > 0 .
It is visible that three ground states are separable: the ferromagnetic saturated spin-5/2 one (given by ψ 5 / 2 , 2 (see Table A1 or Table A2) and two states with spin-1/2, given by ρ ^ 1 / 2 , 0 and ρ ^ 1 / 2 , 1 . As a consequence, only two ground states with spin-3/2 exhibit quantum entanglement.
For the cluster state described by ρ ^ 3 / 2 , 1 , only the spin pairs belonging to the outer tetrahedron are entangled and the concurrence value is 1/6. The form of pair state is:
ρ ^ = 1 2 + 1 4 + 1 12 +
or
ρ ^ = 1 4 ϕ + ϕ + + ϕ ϕ + ϕ + ϕ + ϕ ϕ + + 1 6 ψ + ψ + + 1 3 ψ ψ .
It can be mentioned that the state (15) is an example of probabilistic mixture of all four Bell (maximally entangled) states with additional admixture of non-diagonal states. It is visible that the concurrence value is significantly decreased with respect to the component states (as shown, for example, in Ref. [89] for the specific case of two states in probabilistic mixture).
Moreover, exactly at the boundary of this state with saturated ferromagnetic state, weaker entanglement demonstrated by concurrence equal to 1/8 is present. Due to the topology of the phase diagram for the ferromagnetic J 1 > 0 , the entanglement is present only for a finite range of magnetic fields for J 2 strong enough and of antiferromagnetic sign, as this is the only spin-3/2 state found for J 1 > 0 (see Figure 3a). The same spin-3/2 state is also present in the phase diagram for J 1 < 0 (Figure 3b), but its presence is limited to J 2 / | J 1 | < 1 . For J 2 / | J 1 | > 1 it is replaced with a pure spin-3/2 state ψ 3 / 2 , 2 , for which both kinds of spin pairs are entangled. The pairs involving the central spin indicate quite noticeable concurrence of 2/5 and their state takes the form of:
ρ ^ = 3 20 + 1 20 + 4 5 1 5 +
or
ρ ^ = 3 40 ϕ + ϕ + + ϕ ϕ + ϕ + ϕ + ϕ ϕ + + 9 40 ψ + ψ + + 5 8 ψ ψ 3 8 ψ + ψ + ψ ψ + .
The pairs built of tetrahedron (outer) spins are less entangled, with concurrence value of 1/10 and quantum state:
ρ ^ = 9 10 + 1 20 + + 1 20 +
or
ρ ^ = 9 20 ϕ + ϕ + + ϕ ϕ + ϕ + ϕ + ϕ ϕ + + 1 10 ψ + ψ + .
Let us note that also the boundary of both spin-3/2 states shows the presence of entanglement, as well as the triple point at which two spin-3/2 states and spin-5/2 state have the same energy. Also the boundary between spin-3/2 and saturated ferromagnetic state corresponds to non-vanishing entanglement of both kinds of spin pairs. On the other hand, at the boundary between the pure spin-3/2 state and spin-1/2 state ρ ^ 1 / 2 , 1 only the first kind pairs are weakly entangled, with concurrence value of 14 130 / 40 0.065 . For completeness, we also mention that for J 2 / | J 1 | > 1 / 4 exactly at H = 0 , the spin pairs from the tetrahedron exhibit concurrence value of 1/4 and their state is:
ρ ^ = 1 8 + + 3 8 + 1 4 +
or
ρ ^ = 1 8 ϕ + ϕ + + ϕ ϕ + 1 8 ψ + ψ + + 5 8 ψ ψ .

3.3. Finite Temperature Entanglement

The effect of the finite temperature T > 0 on the quantum entanglement generally consists in creation of a thermal mixture of all available system states [Equation (5)]. In the regions of phase diagram where the ground state is entangled, this thermal admixture of separable states may tend to reduce the concurrence (and in the limit of T the thermal state is always separable). Let us notice that the robustness of entanglement in thermal states focused some general attention [90,91,92] and even for the case of two qubits (spins) the thermal behaviour of entanglement can be rich [93]; noticeably, the separability of the ground state does not imply absence of entanglement at finite temperatures (so that a sort of temperature-induced entanglement can emerge, as observed for example in Ref. [25]). Again, it recalls the problem of entanglement of a probabilistic mixture of pure Hamiltonian states according to Equation (5) in canonical ensemble.
The discussion of the finite-temperature properties will first involve the case of J 1 > 0 (when only the spin pairs of second kind can be entangled) and then the case of J 1 < 0 (when both kinds of spin pairs can exhibit entanglement).

3.3.1. Ferromagnetic J 1 > 0

The effect of finite temperature on the concurrence values for J 1 > 0 can be first tracked in Figure 4a, for spin pairs of second kind, at the normalized temperature of k B T / | J 1 | = 0.2 (compare the ground-state phase diagram in Figure 3a). It is visible that close to the boundaries between the states the entanglement is most effectively reduced (as visible for weakly antiferromagnetic J 2 ). This is due to the fact that close to the boundary between the entangled and separable ground state their energy difference is small, so that the coefficients describing the amount of both states in a thermal mixture become comparable. On the other hand, a range of weak entanglement appears for J 2 / | J 1 | < 1 / 2 for the span of stronger magnetic fields without upper critical magnetic field (above which the ground-state would be separable).
The effect of variable temperature and magnetic field on the concurrence in the case of J 1 > 0 can be followed in Figure 4b, in a contour plot prepared for constant J 2 / | J 1 | = 1 . The influence of the temperature comprises reduction of concurrence and also shifting the critical magnetic field between separable and entangled states to higher values. At finite temperatures, the entanglement is present above the critical magnetic field, exceeding always the value of H / | J 1 | = 1 / 2 . This sort of boundary dividing the quantum thermal states to separable and entangled ones can be plotted separately in normalized temperature-normalized magnetic field plane for various values of J 2 / | J 1 , as it is done in Figure 5. For the values of J 2 / | J 1 | > 1 / 2 , the critical magnetic field increases very fast with the temperature, with the value of critical field equal to 0 at T = 0 . For couplings J 2 / | J 1 | < 1 / 2 the nonzero critical magnetic field for the onset of entanglement occurs; it increases with rising temperature and also with antiferromagnetic coupling J 2 magnitude.
The cross-sections of Figure 4b are shown in Figure 6a for various selections of constant magnetic field and variable temperature and Figure 6b for various selections of constant temperature and variable magnetic field. In Figure 6a, for lower magnetic fields, the concurrence decreases monotonically with temperature from the initial, ground-state value of 1/6 characteristic of the state ρ ^ 3 / 2 , 1 . Exactly at H / | J 1 | = 3 / 2 the initial value switches to 1/8 when the state ρ ^ 3 / 2 , 1 gains statistical admixture of separable saturated ferromagnetic state ψ 5 / 2 , 2 , but the overall behaviour as a function of temperature is unchanged. On the contrary, for the higher fields, the temperature is a factor causing the emergence of entanglement, which is absent in the ground state ρ ^ 1 / 2 , 0 . For this regime, a maximum of concurrence is formed when T increases (due to the thermal admixture of spin-3/2 states), but its height is reduced when the magnetic field rises (which is the factor favouring energetically the saturated spin-5/2 separable state). The effect of the increasing magnetic field can be tracked in Figure 6b. In all studied temperatures the state is separable below certain critical magnetic field, which switches on the entanglement (due to transition to the state with significant participation of the entangled state ρ ^ 3 / 2 , 1 at finite T) and a maximum is formed (with values up to 1/6 for the lowest temperatures, as in ρ ^ 3 / 2 , 1 state). At higher fields the concurrence tends to vanish gradually as the separable state with spin 5/2 dominates. This is an example of competing influence of temperature and magnetic field on the entanglement. For lower magnetic fields, the magnetic field increases the amount of entangled ρ ^ 3 / 2 , 1 state in the thermal state, boosting the concurrence, whereas the thermal fluctuations at elevated temperature restore the higher amount of separable ground state ρ ^ 1 / 2 , 0 lowering the concurrence value. At higher fields the situation is opposite: further increase of the field enriches the thermal state with saturated ferromagnetic separable state reducing the concurrence, whereas the thermal fluctuations at elevated temperature restore the higher amount of entangled state ρ ^ 3 / 2 , 1 boosting the concurrence value.

3.3.2. Antiferromagnetic J 1 < 0

For antiferromagnetic J 1 , both types of spin pairs can exhibit entanglement. A finite temperature ( k B T / | J 1 | = 0.2 ) contour plot of concurrence is shown in the J 2 / | J 1 | H / | J 1 | plane in Figure 7, for type 1 (Figure 7a) and type 2 (Figure 7b) spin pairs. Again, it resembles a thermally diffused version of the relevant ground-state phase diagram shown in Figure 3b. For spin pairs connected with J 1 , the most pronounced entanglement is present for the parameters range corresponding to ψ 3 / 2 , 2 ground state in Figure 2b, with values reduced close to the boundaries of this range due to thermal admixture of separable states. On the other hand, the thermal state for J 2 > 0 and considerably high magnetic field acquires some residual entanglement due to thermal effects (as opposed to separable saturated ferromagnetic ground state). A similar qualitative picture can be tracked in Figure 7b for concurrence in type 2 spin pairs. Both entangled ranges of the phase diagram (Figure 3b) corresponding to ground states with S T = 3 / 2 reduce their entanglement and become separated when the temperature is elevated, whereas the saturated ferromagnetic state acquires residual concurrence by admixture of other states.
The influence of the temperature and magnetic field on the entanglement for both types of spin pairs can be followed in contour plots in Figure 8, prepared for J 2 / | J 1 | = 1. For type 1 pairs (Figure 8a), the maximum concurrence is achieved for a finite magnetic field which does not shift significantly when the temperature is elevated (and there is no minimum field necessary to trigger the entanglement). Moreover, the boundary dividing the phase diagram into separable and entangled states corresponds to relatively high temperature and is relatively insensitive to the field. For type 2 pairs (Figure 8b), the magnetic field at which the concurrence peaks shifts to higher values when the temperature rises and the entanglement is absent below certain critical field. Contrary to the case of type 1 pairs, the boundary separating the entangled and separable states is linear at lower temperatures and its slope increases for higher values of temperature.
The mentioned boundary lines dividing entangled and separable states can be plotted complementarily in temperature-magnetic field plane in Figure 9 for various selected values of J 2 / | J 1 | . For type-1 pairs (Figure 9a), the boundary for J 2 > 0 is rather weakly dependent on the applied field and the non-zero critical temperature above which the concurrence vanishes decreases when J 2 coupling tends to 0. Exactly for J 2 = 0 (i.e., spins belonging to the tetrahedron uncoupled to each other) the boundary commences at T = 0 and the critical temperature rises when stronger field is applied. For antiferromagnetic J 2 < 0 this tendency remains, but the entanglement is present only above certain critical field which is increased when J 2 becomes stronger. For type-2 pairs (Figure 9b), for ferromagnetic or weakly antiferromagnetic J 2 , the entanglement is present only at finite magnetic fields and this minimum field increases when the temperature is elevated; in addition also reducing the coupling J 2 results in increase of the field. For J 2 / | J 1 | > 1 / 4 all the boundaries commence at T = 0 and H = 0 . If the coupling J 2 is more antiferromagnetic, the critical magnetic field necessary to trigger the entanglement is non-zero even for T = 0 and it is increased when the temperature rises (the overall temperature range for which concurrence is positive gets reduced when the antiferromagnetic coupling J 2 becomes stronger). For even stronger antiferromagnetic J 2 , the temperature dependence of the critical field tends to flatten.
The detailed temperature dependence of concurrence for both considered spin pairs for selected magnetic fields can be tracked in Figure 10. The case of ferromagnetic J 2 / | J 1 | = 1 is depicted in Figure 10a for type 1 pairs and Figure 10b for type 2 pairs. For type 1 pairs, in the absence of the field C ( 1 ) takes the ground state value of 1 / 4 at T = 0 (characteristic of probabilistic mixture of states ψ 3 / 2 , S T , z , 2 with all allowed values of S T , z ) and monotonically decreases to 0 with increasing slope when the temperature rises (with a long plateau at the lowest temperatures). For H > 0 , the initial ground-state value of concurrence is equal to 2 / 5 (as the ground state is ψ 3 / 2 , 2 ), but C ( 1 ) first drops down fast to approximately 1 / 4 (because for weak field the states ψ 3 / 2 , S T , z , 2 with all allowed values of S T , z lie close in energy and form probabilistic mixture as thermal fluctuations arise) and then further decreases when T is elevated. For stronger fields the temperature dependence of C ( 1 ) becomes more regular. Exactly for H / | J 1 | = 5/2 the initial value of concurrence is 1/5 (as the mixture of states with spin-3/2 and spin-5/2 is created) with further monotonic decrease. At stronger fields the behaviour changes qualitatively, as the entanglement vanishes at T = 0 (as the ground state is saturated ferromagnetic one) and a local maximum builds up when the temperature increases (so that we deal with thermally induced entanglement). For the strongest field the maximum height becomes reduced and its position is shifted to higher temperatures. The behaviour of C ( 2 ) is shown in Figure 10b. For this kind of spin pairs the initial value of concurrence for low temperatures is 1/10 for zero field or low fields (characteristic of ground state ψ 3 / 2 , 2 ), with monotonic behaviour as a function of the temperature. The same sort of behaviour is seen for H / | J 1 | = 5/2 but with the initial value of 1/20 (as the ground state involves entangled state with spin 3/2 and separable ferromagnetic one with spin 5/2). For stronger fields a temperature-induced local maximum of concurrence is formed, like in the case of type 1 spin pair.
When the coupling J 2 is antiferromagnetic, the situation is qualitatively different. For type 1 pairs (as shown in Figure 10c) at weak magnetic field the entanglement is only induced by finite temperature, with a local maximum (as the ground state is separable ρ ^ 1 / 2 , 1 and only the thermal admixture of spin-3/2 state induces the entanglement). When the field increases, the ground-state entanglement emerges (either due to ground state being ψ 3 / 2 , 2 or owing to its admixture to other states at phase boundaries) and the temperature dependence takes the form of monotonic decrease (with varying initial values in concert with the diagram in Figure 3b). For even stronger fields the picture switches back to the scenario with temperature-induced local maximum of concurrence (as the ground state becomes saturated separable ferromagnetic state). In the case of type 2 pairs (Figure 10d), the monotonic decrease from the initial value of 1/10 (characteristic of entangled ground state ψ 3 / 2 , 2 ) is noticeable for weaker fields, whereas for stronger ones a local maximum of temperature-induced concurrence is present (and the entanglement is thermally induced by admixing spin-3/2 states to saturated ferromagnetic ground state).
Tracking of the detailed magnetic field dependence of concurrence for various temperatures is possible in Figure 11. For type 1 spin pairs and ferromagnetic J 2 / | J 1 | = 1 (Figure 11a), at lower temperatures the concurrence takes in the absence of the field the value of 1/4 (characteristic of mixed state involving ψ 3 / 2 , S T , z , 2 with all allowed values of S T , z , see Figure 3b) and fast jumps to 2/5 (stemming from ψ 3 / 2 , 2 ground state), reaching a plateau persisting up to certain critical magnetic field. The increasing temperature first tends to reduce the value at the plateau (as the thermal state acquires contribution from saturated ferromagnetic state) and then also diminishes the initial concurrence value at H = 0. At the same time the range of entangled states is extended to higher magnetic fields, as the thermal admixture of entangled spin-3/2 state builds up entanglement in the range where the ground state has spin 5/2. For type 2 pairs (Figure 11b), the entanglement is induced by the finite field H. At the lowest temperatures its field dependence reaches the highest plateau with the value of 1/10 (emerging for ψ 3 / 2 , 2 ground state), and the increasing temperature smears this maximum, increases a threshold magnetic field for the onset of entanglement and makes the entanglement more robust in the range of higher fields. This behaviour exemplifies again the competing effect of temperature and magnetic field on the entanglement, as discussed in the context of Figure 6b. The magnetic field influence on the concurrence of type 1 spin pairs for antiferromagnetic J 2 / | J 1 | = −0.5 (Figure 11c) is somehow similar as in the previous case (the entanglement is field-induced); however, the maximum corresponds to the values reaching 2/5. The situation illustrated in Figure 11d for type 2 pairs resembles even more the scenario for ferromagnetic J 2 .

4. Conclusions

In the paper we have reported computational study of a pentamer spin cluster (with single spin located in the center of spin tetrahedron) composed of spins S = 1 / 2 and described with isotropic Heisenberg model with external magnetic field and two exchange integrals: one for spin pairs belonging to tetrahedron and another one quantifying the coupling between tetrahedron spins and a central spin. The cluster has non-trivial geometry characterized by high symmetry (being the smallest non-planar Kuratowski graph K 5 ). The selected geometry offers the possibility of influencing the ratio of exchange integrals by chemical composition and of creating networks of weakly interacting clusters [75]. The interest of our study was focused on the quantum two-spin entanglement properties. Both ground-state properties and the effect of finite temperature on the system behaviour were discussed, preceded by the analysis of the ground-state phase diagram as a function of exchange integrals and external magnetic field. A specific feature of the selected cluster is relatively frequent occurrence of ground state degeneracy due to high system symmetry. We found the presence of two-spin entanglement for two ground states (a pure one and a mixed one), both with total cluster spin number equal to S T = 3 / 2 . The studied cluster exhibits the entanglement preferably for all-antiferromagnetic couplings, which is an usual pattern in the case of magnetic materials [94]. For the case of finite temperature, the phenomenon of temperature-induced and magnetic field-induced entanglement were predicted. The behaviour of the Wootters concurrence was extensively discussed, developing the interest in thermal entanglement and magnetic entanglement, as called in pioneering work, Ref. [25] exposing the interplay of the temperature and magnetic field when influencing the entanglement of naturally occurring thermal quantum states in magnetic nanoclusters.
The system in question was inspired by synthesis of Cu-based cluster molecular magnets containing spins 1/2 [77] and Co-based structures [74] built of higher spins. Our selection of cluster containing localized spins S = 1 / 2 (equivalent to qubits) was additionally justified by expectation of most pronounced quantum behaviour for the lowest possible spin. However, the extension of the study to higher spins would further augment the parameter space, leading to enriched phase diagram and possibility to explore entanglement for more general qudit case. Another direction of extending the study would involve incorporation of intercluster couplings to correlate the model better to cluster molecular magnets in which such couplings occur (see a general model study [95] or particular examples of coupled dimers [96,97]).

Funding

This research received no external funding.

Data Availability Statement

The data presented in this study are available on reasonable request from the corresponding author.

Conflicts of Interest

The author declares no conflict of interest.

Sample Availability

No samples were used in this study.

Appendix A. Complete Set of the Eigenvalues and Eigenstates of the Cluster Hamiltonian

The appendix contains the complete set of eigenstates of cluster Hamiltonian given by Equation (1) with relevant eigenenergies and sets of quantum numbers S T , S T , z and s [related to the eigenenergies by Equation (2)]. The states are listed in two forms: Table A1 contains the eigenstates written in the basis using only spin-up and spin-down states for single spins, whereas Table A2 manifestly emphasizes division of the system to central spin and two dimers using Bell states for spin pairs if possible.
Table A1. Eigenenergies and eigenstates of the Hamiltonian [Equation (1)] with quantum numbers.
Table A1. Eigenenergies and eigenstates of the Hamiltonian [Equation (1)] with quantum numbers.
S T , S T , z , s EigenenergyEigenstate
5 2 , + 5 2 , 2 J 1 3 2 J 2 5 2 H ψ 5 2 , + 5 2 , 2 =
5 2 , + 3 2 , 2 J 1 3 2 J 2 3 2 H ψ 5 2 , + 3 2 , 2 = 1 5 + + + +
5 2 , 1 2 , 2 J 1 3 2 J 2 1 2 H ψ 5 2 , + 1 2 , 2 = 1 10 + + + + +
+ + + +
5 2 , 1 2 , 2 J 1 3 2 J 2 + 1 2 H ψ 5 2 , 1 2 , 2 = 1 10 + + + + +
+ + + +
5 2 , 3 2 , 2 J 1 3 2 J 2 + 3 2 H ψ 5 2 , 3 2 , 2 = 1 5 + + + +
5 2 , 5 2 , 2 J 1 3 2 J 2 + 5 2 H ψ 5 2 , 5 2 , 2 =
3 2 , + 3 2 , 2 3 2 J 1 3 2 J 2 3 2 H ψ 3 2 , + 3 2 , 2 = 1 2 5 + + + 4
3 2 , + 1 2 , 2 3 2 J 1 3 2 J 2 1 2 H ψ 3 2 , + 1 2 , 2 = 1 2 15 2 + 2 + 2 + 2 + 2 + 2
3 3 3 3
3 2 , 1 2 , 2 3 2 J 1 3 2 J 2 + 1 2 H ψ 3 2 , 1 2 , 2 = 1 2 15 2 + 2 + 2 + 2 + 2 + 2
3 3 3 3
3 2 , 3 2 , 2 3 2 J 1 3 2 J 2 + 3 2 H ψ 3 2 , 3 2 , 2 = 1 2 5 + + + 4
3 2 , + 3 2 , 1 1 2 J 1 + 1 2 J 2 3 2 H ψ 3 2 , + 3 2 , 1 ( 1 ) = 1 2 +
ψ 3 2 , + 3 2 , 1 ( 2 ) = 1 2 +
ψ 3 2 , + 3 2 , 1 ( 3 ) = 1 2 +
3 2 , + 1 2 , 1 1 2 J 1 + 1 2 J 2 1 2 H ψ 3 2 , + 1 2 , 1 ( 1 ) = 1 2 3 + + 2 2
ψ 3 2 , + 1 2 , 1 ( 2 ) = 1 2 3 + + 2 2
ψ 3 2 , + 1 2 , 1 ( 3 ) = 1 2 3 + 2 + 2
3 2 , 1 2 , 1 1 2 J 1 + 1 2 J 2 + 1 2 H ψ 3 2 , 1 2 , 1 ( 3 ) = 1 2 3 + 2 + 2
ψ 3 2 , 1 2 , 1 ( 2 ) = 1 2 3 + + 2 2
ψ 3 2 , 1 2 , 1 ( 1 ) = 1 2 3 + + 2 2
3 2 , 3 2 , 1 1 2 J 1 + 1 2 J 2 + 3 2 H ψ 3 2 , 3 2 , 1 ( 3 ) = 1 2 +
ψ 3 2 , 3 2 , 1 ( 2 ) = 1 2 +
ψ 3 2 , 3 2 , 1 ( 1 ) = 1 2 +
1 2 , + 1 2 , 1 J 1 + 1 2 J 2 1 2 H ψ 1 2 , + 1 2 , 1 ( 1 ) = 1 6 + +
ψ 1 2 , + 1 2 , 1 ( 2 ) = 1 6 + + +
ψ 1 2 , + 1 2 , 1 ( 3 ) = 1 6 + +
1 2 , 1 2 , 1 J 1 + 1 2 J 2 + 1 2 H ψ 1 2 , 1 2 , 1 ( 3 ) = 1 6 + +
ψ 1 2 , 1 2 , 1 ( 2 ) = 1 6 + + +
ψ 1 2 , 1 2 , 1 ( 1 ) = 1 6 + +
1 2 , + 1 2 , 0 3 2 J 2 1 2 H ψ 1 2 , + 1 2 , 0 ( 1 ) = 1 2 +
ψ 1 2 , + 1 2 , 0 ( 2 ) = 1 2 3 2 + 2
1 2 , 1 2 , 0 3 2 J 2 + 1 2 H ψ 1 2 , 1 2 , 0 ( 2 ) = 1 2 3 2 + 2
ψ 1 2 , 1 2 , 0 ( 1 ) = 1 2 +
Table A2. Eigenenergies and eigenstates of the Hamiltonian [Equation (1)] with quantum numbers, expressed in a basis emphasizing the system division into central spin and two dimers.
Table A2. Eigenenergies and eigenstates of the Hamiltonian [Equation (1)] with quantum numbers, expressed in a basis emphasizing the system division into central spin and two dimers.
S T , S T , z , s EigenenergyEigenstate
5 2 , + 5 2 , 2 J 1 3 2 J 2 5 2 H ψ 5 2 , + 5 2 , 2 =
5 2 , + 3 2 , 2 J 1 3 2 J 2 3 2 H ψ 5 2 , + 3 2 , 2 = 2 5 ψ + + ψ + + 1 5
5 2 , + 1 2 , 2 J 1 3 2 J 2 1 2 H ψ 5 2 , + 1 2 , 2 = 1 10 ϕ + ϕ + ϕ ϕ + 2 ψ + ψ + + 1 5 ψ + + ψ +
5 2 , 1 2 , 2 J 1 3 2 J 2 + 1 2 H ψ 5 2 , 1 2 , 2 = 1 10 ϕ + ϕ + ϕ ϕ + 2 ψ + ψ + + 1 5 ψ + + ψ +
5 2 , 3 2 , 2 J 1 3 2 J 2 + 3 2 H ψ 5 2 , 3 2 , 2 = 2 5 ψ + + ψ + + 1 5
5 2 , 5 2 , 2 J 1 3 2 J 2 + 5 2 H ψ 5 2 , 5 2 , 2 =
3 2 , + 3 2 , 2 3 2 J 1 3 2 J 2 3 2 H ψ 3 2 , + 3 2 , 2 = 1 10 ψ + + ψ + 2 5
3 2 , + 1 2 , 2 3 2 J 1 3 2 J 2 1 2 H ψ 3 2 , + 1 2 , 2 = 1 15 ϕ + ϕ + ϕ ϕ + 2 ψ + ψ + 3 30 ψ + + ψ +
3 2 , 1 2 , 2 3 2 J 1 3 2 J 2 + 1 2 H ψ 3 2 , 1 2 , 2 = 1 15 ϕ + ϕ + ϕ ϕ + 2 ψ + ψ + 3 30 ψ + + ψ +
3 2 , 3 2 , 2 3 2 J 1 3 2 J 2 + 3 2 H ψ 3 2 , 3 2 , 2 = 1 10 ψ + + ψ + 2 5
3 2 , + 3 2 , 1 1 2 J 1 + 1 2 J 2 3 2 H ψ 3 2 , + 3 2 , 1 ( 1 ) = 1 2 ψ + ψ
ψ 3 2 , + 3 2 , 1 ( 2 ) = 1 2 ψ ψ
ψ 3 2 , + 3 2 , 1 ( 3 ) = 1 2 ψ + ψ +
3 2 , + 1 2 , 1 1 2 J 1 + 1 2 J 2 1 2 H ψ 3 2 , + 1 2 , 1 ( 1 ) = 1 6 ψ + ψ + 1 3 ϕ ϕ + ϕ + ϕ
ψ 3 2 , + 1 2 , 1 ( 2 ) = 1 6 ψ ψ + 1 3 ψ ψ + ψ + ψ
ψ 3 2 , + 1 2 , 1 ( 3 ) = 1 6 ψ + ψ + 1 3 ψ ψ + + ψ + ψ
3 2 , 1 2 , 1 1 2 J 1 + 1 2 J 2 + 1 2 H ψ 3 2 , 1 2 , 1 ( 3 ) = 1 6 ψ + ψ + 1 3 ψ ψ + + ψ + ψ
ψ 3 2 , 1 2 , 1 ( 2 ) = 1 6 ψ ψ + 1 3 ψ ψ + ψ + ψ
ψ 3 2 , 1 2 , 1 ( 1 ) = 1 6 ψ + ψ + + 1 3 ϕ ϕ + ϕ + ϕ
3 2 , 3 2 , 1 1 2 J 1 + 1 2 J 2 + 3 2 H ψ 3 2 , 3 2 , 1 ( 3 ) = 1 2 ψ + ψ +
ψ 3 2 , 3 2 , 1 ( 2 ) = 1 2 ψ ψ
ψ 3 2 , 3 2 , 1 ( 1 ) = 1 2 ψ + ψ
1 2 , + 1 2 , 1 J 1 + 1 2 J 2 1 2 H ψ 1 2 , + 1 2 , 1 ( 1 ) = 1 6 ψ + ψ + ψ ψ + 1 3 ψ + ψ
ψ 1 2 , + 1 2 , 1 ( 2 ) = 1 6 ψ + ψ + ψ ψ + + 1 3 ψ ψ
ψ 1 2 , + 1 2 , 1 ( 3 ) = 1 6 ψ + ψ + ψ ψ + 1 3 ψ + ψ +
1 2 , 1 2 , 1 J 1 + 1 2 J 2 + 1 2 H ψ 1 2 , 1 2 , 1 ( 3 ) = 1 6 ψ + ψ + ψ ψ + + 1 3 ψ + ψ +
ψ 1 2 , 1 2 , 1 ( 2 ) = 1 6 ψ + ψ + ψ ψ + + 1 3 ψ ψ
ψ 1 2 , 1 2 , 1 ( 1 ) = 1 6 ψ + ψ + ψ ψ + 1 3 ψ + ψ
1 2 , + 1 2 , 0 3 2 J 2 1 2 H ψ 1 2 , + 1 2 , 0 ( 1 ) = ψ ψ
ψ 1 2 , + 1 2 , 0 ( 2 ) = 1 3 ϕ + ϕ + ϕ ϕ ψ + ψ +
1 2 , 1 2 , 0 3 2 J 2 + 1 2 H ψ 1 2 , 1 2 , 0 ( 2 ) = 1 3 ϕ + ϕ + ϕ ϕ ψ + ψ +
ψ 1 2 , 1 2 , 0 ( 1 ) = ψ ψ

References

  1. Gatteschi, D. Molecular Magnetism: A Basis for New Materials. Adv. Mater. 1994, 6, 635–645. [Google Scholar] [CrossRef]
  2. Gatteschi, D.; Bogani, L.; Cornia, A.; Mannini, M.; Sorace, L.; Sessoli, R. Molecular Magnetism, Status and Perspectives. Solid State Sci. 2008, 10, 1701–1709. [Google Scholar] [CrossRef]
  3. Dei, A.; Gatteschi, D. Molecular (Nano) Magnets as Test Grounds of Quantum Mechanics. Angew. Chem. Int. Ed. 2011, 50, 11852–11858. [Google Scholar] [CrossRef] [PubMed]
  4. Moreno-Pineda, E.; Wernsdorfer, W. Measuring Molecular Magnets for Quantum Technologies. Nat. Rev. Phys. 2021, 3, 645–659. [Google Scholar] [CrossRef]
  5. Evangelisti, M.; Brechin, E.K. Recipes for Enhanced Molecular Cooling. Dalton Trans. 2010, 39, 4672–4676. [Google Scholar] [CrossRef] [PubMed]
  6. Konieczny, P.; Sas, W.; Czernia, D.; Pacanowska, A.; Fitta, M.; Pełka, R. Magnetic Cooling: A Molecular Perspective. Dalton Trans. 2022, 51, 12762–12780. [Google Scholar] [CrossRef]
  7. Coronado, E. Molecular Magnetism: From Chemical Design to Spin Control in Molecules, Materials and Devices. Nat. Rev. Mater. 2020, 5, 87–104. [Google Scholar] [CrossRef]
  8. Chilton, N.F. Molecular Magnetism. Annu. Rev. Mater. Res. 2022, 52, 79–101. [Google Scholar] [CrossRef]
  9. Cruz, C. Quantum Entanglement in Low-Dimensional Metal Complexes: An Outlook. Phys. B Condens. Matter 2023, 653, 414485. [Google Scholar] [CrossRef]
  10. Moneo-Corcuera, A.; Nieto-Castro, D.; Cirera, J.; Gómez, V.; Sanjosé-Orduna, J.; Casadevall, C.; Molnár, G.; Bousseksou, A.; Parella, T.; Martínez-Agudo, J.M.; et al. Molecular Memory near Room Temperature in an Iron Polyanionic Complex. Chem 2023, 9, 377–393. [Google Scholar] [CrossRef]
  11. Leuenberger, M.N.; Loss, D. Quantum Computing in Molecular Magnets. Nature 2001, 410, 789–793. [Google Scholar] [CrossRef]
  12. Stepanenko, D.; Trif, M.; Loss, D. Quantum Computing with Molecular Magnets. Inorganica Chim. Acta 2008, 361, 3740–3745. [Google Scholar] [CrossRef]
  13. Ghirri, A.; Candini, A.; Affronte, M. Molecular Spins in the Context of Quantum Technologies. Magnetochemistry 2017, 3, 12. [Google Scholar] [CrossRef]
  14. Gaita-Ariño, A.; Luis, F.; Hill, S.; Coronado, E. Molecular Spins for Quantum Computation. Nat. Chem. 2019, 11, 301–309. [Google Scholar] [CrossRef] [PubMed]
  15. Wasielewski, M.R.; Forbes, M.D.E.; Frank, N.L.; Kowalski, K.; Scholes, G.D.; Yuen-Zhou, J.; Baldo, M.A.; Freedman, D.E.; Goldsmith, R.H.; Goodson, T.; et al. Exploiting Chemistry and Molecular Systems for Quantum Information Science. Nat. Rev. Chem. 2020, 4, 490–504. [Google Scholar] [CrossRef]
  16. Galindo, A.; Martín-Delgado, M.A. Information and Computation: Classical and Quantum Aspects. Rev. Mod. Phys. 2002, 74, 347–423. [Google Scholar] [CrossRef]
  17. Amico, L.; Fazio, R.; Osterloh, A.; Vedral, V. Entanglement in Many-Body Systems. Rev. Mod. Phys. 2008, 80, 517–576. [Google Scholar] [CrossRef]
  18. Horodecki, R.; Horodecki, P.; Horodecki, M.; Horodecki, K. Quantum Entanglement. Rev. Mod. Phys. 2009, 81, 865–942. [Google Scholar] [CrossRef]
  19. Latorre, J.I.; Riera, A. A Short Review on Entanglement in Quantum Spin Systems. J. Phys. A Math. Theor. 2009, 42, 504002. [Google Scholar] [CrossRef]
  20. Laflorencie, N. Quantum Entanglement in Condensed Matter Systems. Phys. Rep. 2016, 646, 1–59. [Google Scholar] [CrossRef]
  21. Bose, I.; Tribedi, A. Thermal Entanglement Properties of Small Spin Clusters. Phys. Rev. A 2005, 72, 022314. [Google Scholar] [CrossRef]
  22. Tribedi, A.; Bose, I. Entangled Spin Clusters: Some Special Features. Phys. Rev. A 2006, 74, 012314. [Google Scholar] [CrossRef]
  23. Souza, A.M.; Reis, M.S.; Soares-Pinto, D.O.; Oliveira, I.S.; Sarthour, R.S. Experimental Determination of Thermal Entanglement in Spin Clusters Using Magnetic Susceptibility Measurements. Phys. Rev. B 2008, 77, 104402. [Google Scholar] [CrossRef]
  24. Pal, A.K.; Bose, I. Quantum Discord in the Ground and Thermal States of Spin Clusters. J. Phys. B At. Mol. Opt. Phys. 2011, 44, 045101. [Google Scholar] [CrossRef]
  25. Arnesen, M.C.; Bose, S.; Vedral, V. Natural Thermal and Magnetic Entanglement in the 1D Heisenberg Model. Phys. Rev. Lett. 2001, 87, 017901. [Google Scholar] [CrossRef] [PubMed]
  26. Lagmago Kamta, G.; Starace, A.F. Anisotropy and Magnetic Field Effects on the Entanglement of a Two Qubit Heisenberg XY Chain. Phys. Rev. Lett. 2002, 88, 107901. [Google Scholar] [CrossRef]
  27. Guo-Feng, Z.; Jiu-Qing, L.; Qi-Wei, Y. Thermal Entanglement in Spin-Dimer V4+ with a Strong Magnetic Field. Chin. Phys. Lett. 2003, 20, 452. [Google Scholar] [CrossRef]
  28. Terzis, A.F.; Paspalakis, E. Entanglement in a Two-Qubit Ising Model under a Site-Dependent External Magnetic Field. Phys. Lett. A 2004, 333, 438–445. [Google Scholar] [CrossRef]
  29. Soares-Pinto, D.O.; Souza, A.M.; Sarthour, R.S.; Oliveira, I.S.; Reis, M.S.; Brandão, P.; Rocha, J.; dos Santos, A.M. Entanglement Temperature in Molecular Magnets Composed of S-spin Dimers. Europhys. Lett. 2009, 87, 40008. [Google Scholar] [CrossRef]
  30. Del Cima, O.M.; Franco, D.H.T.; Silva, M.M. Magnetic Shielding of Quantum Entanglement States. Quantum Stud. Math. Found. 2019, 6, 141–150. [Google Scholar] [CrossRef]
  31. Čenčariková, H.; Strečka, J. Unconventional Strengthening of the Bipartite Entanglement of a Mixed Spin-(1/2,1) Heisenberg Dimer Achieved through Zeeman Splitting. Phys. Rev. B 2020, 102, 184419. [Google Scholar] [CrossRef]
  32. Ghannadan, A.; Strečka, J. Magnetic-Field-Orientation Dependent Thermal Entanglement of a Spin-1 Heisenberg Dimer: The Case Study of Dinuclear Nickel Complex with an Uniaxial Single-Ion Anisotropy. Molecules 2021, 26, 3420. [Google Scholar] [CrossRef] [PubMed]
  33. Vargová, H.; Strečka, J. 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. Nanomaterials 2021, 11, 3096. [Google Scholar] [CrossRef]
  34. Vargová, H.; Strečka, J. Conventional and Inverse Magnetocaloric and Electrocaloric Effects of a Mixed Spin-(1/2, 1) Heisenberg Dimer. Eur. Phys. J. Plus 2022, 137, 490. [Google Scholar] [CrossRef]
  35. Naveena, P.; Bhuvaneswari, S.; Radha, R.; Muthuganesan, R. Quantum Correlations in a Mixed Spin-(1/2,1) Heisenberg Dimer. J. Magn. Magn. Mater. 2022, 563, 169863. [Google Scholar] [CrossRef]
  36. Oumennana, M.; Mansour, M. Quantum Coherence versus Quantum-Memory-Assisted Entropic Uncertainty Relation in a Mixed Spin-(1/2,1) Heisenberg Dimer. Opt. Quantum Electron. 2023, 55, 594. [Google Scholar] [CrossRef]
  37. Wang, X.; Fu, H.; Solomon, A.I. Thermal Entanglement in Three-Qubit Heisenberg Models. J. Phys. A Math. Gen. 2001, 34, 11307. [Google Scholar] [CrossRef]
  38. Li, C.X.; Wang, C.Z.; Lin, X.M.; Guo, G.C. Thermal Entanglement versus Mixture in a Spin Chain: Generation of Maximally Entangled Mixed States. Phys. A Stat. Mech. Its Appl. 2006, 370, 430–440. [Google Scholar] [CrossRef]
  39. Liu, R.; Liang, M.L.; Yuan, B. Entanglement in the Three-Qubit Heisenberg Model with next Nearest Neighbor Interaction and a Nonuniform Magnetic Field. Eur. Phys. J. D 2007, 41, 571–578. [Google Scholar] [CrossRef]
  40. Pal, A.K.; Bose, I. Entanglement in a Molecular Three-Qubit System. J. Phys. Condens. Matter 2009, 22, 016004. [Google Scholar] [CrossRef]
  41. Cima, O.M.D.; Franco, D.H.T.; da Silva, S.L.L. Quantum Entanglement in Trimer Spin-1/2 Heisenberg Chains with Antiferromagnetic Coupling. Quantum Stud. Math. Found. 2016, 3, 57–63. [Google Scholar] [CrossRef]
  42. Milivojević, M. Maximal Thermal Entanglement Using Three-Spin Interactions. Quantum Inf. Process. 2019, 18, 48. [Google Scholar] [CrossRef]
  43. Benabdallah, F.; Anouz, K.E.; Strecka, J.; Daoud, M. Thermal Non-Classical Correlation via Skew Information, Quantum Fisher Information, and Quantum Teleportation of a Spin-1/2 Heisenberg Trimer System. Eur. Phys. J. Plus 2022, 137, 1096. [Google Scholar] [CrossRef]
  44. Wu, K.D.; Zhou, B.; Cao, W.Q. Thermal Entanglement in a Four-Qubit Heisenberg Spin Model with External Magnetic Fields. Phys. Lett. A 2007, 362, 381–389. [Google Scholar] [CrossRef]
  45. Irons, H.R.; Quintanilla, J.; Perring, T.G.; Amico, L.; Aeppli, G. Control of Entanglement Transitions in Quantum Spin Clusters. Phys. Rev. B 2017, 96, 224408. [Google Scholar] [CrossRef]
  46. Karľová, K.; Strečka, J. Interplay of Bipartite Entanglement between Two Geometrically Inequivalent Spin Pairs of a Spin-1/2 Heisenberg Distorted Tetrahedron. Acta Phys. Pol. A 2020, 137, 595–597. [Google Scholar] [CrossRef]
  47. Karlóvá, K.; Strečka, J. Distribution of Bipartite and Tripartite Entanglement within a Spin-1/2 Heisenberg Star in a Magnetic Field. Molecules 2023, 28, 4037. [Google Scholar] [CrossRef] [PubMed]
  48. Szałowski, K. Two-Spin and Multi-Spin Quantum Entanglement in V12 Polyoxovanadate Molecular Nanomagnet. J. Magn. Magn. Mater. 2022, 546, 168782. [Google Scholar] [CrossRef]
  49. Zad, H.A.; Zoshki, A.; Ananikian, N.; Jaščur, M. Robust Quantum Entanglement and Teleportation in the Tetrapartite Spin-1/2 Square Clusters: Theoretical Study on the Effect of a Cyclic Four-Spin Exchange. J. Magn. Magn. Mater. 2022, 559, 169533. [Google Scholar] [CrossRef]
  50. Ghannadan, A.; Karl’ová, K.; Strečka, J. On the Concurrent Bipartite Entanglement of a Spin-1 Heisenberg Diamond Cluster Developed for Tetranuclear Nickel Complexes. Magnetochemistry 2022, 8, 156. [Google Scholar] [CrossRef]
  51. Kuzmak, A.R. Entanglement of the Ising–Heisenberg Diamond Spin- Cluster in Evolution. J. Phys. A Math. Theor. 2023, 56, 165302. [Google Scholar] [CrossRef]
  52. Vargová, H.; Strečka, J. Distribution of a Bipartite Entanglement in a Mixed Spin-(1/2,1) Heisenberg Tetramer. Phys. A Stat. Mech. Its Appl. 2023, 625, 129046. [Google Scholar] [CrossRef]
  53. Chiesa, A.; Petiziol, F.; Chizzini, M.; Santini, P.; Carretta, S. Theoretical Design of Optimal Molecular Qudits for Quantum Error Correction. J. Phys. Chem. Lett. 2022, 13, 6468–6474. [Google Scholar] [CrossRef] [PubMed]
  54. Chizzini, M.; Crippa, L.; Chiesa, A.; Tacchino, F.; Petiziol, F.; Tavernelli, I.; Santini, P.; Carretta, S. Molecular Nanomagnets with Competing Interactions as Optimal Units for Qudit-Based Quantum Computation. Phys. Rev. Res. 2022, 4, 043135. [Google Scholar] [CrossRef]
  55. Rycerz, A. Pairwise Entanglement and the Mott Transition for Correlated Electrons in Nanochains. New J. Phys. 2017, 19, 053025. [Google Scholar] [CrossRef]
  56. Sousa, H.S.; Pereira, M.S.S.; de Oliveira, I.N.; Strečka, J.; Lyra, M.L. Phase Diagram and Re-Entrant Fermionic Entanglement in a Hybrid Ising-Hubbard Ladder. Phys. Rev. E 2018, 97, 052115. [Google Scholar] [CrossRef]
  57. Souza, F.; Almeida, G.M.A.; Lyra, M.L.; Pereira, M.S.S. Interplay between Charge and Spin Thermal Entanglement in Hubbard Dimers. Phys. Rev. A 2020, 102, 032421. [Google Scholar] [CrossRef]
  58. Lo Schiavo, M.; Cipriani, F.; De Riso, G.; Romano, A.; Noce, C. Quantum entanglement in an extended Hubbard model as evaluated from a spin concurrence measure. J. Magn. Magn. Mater. 2023, 584, 171066. [Google Scholar] [CrossRef]
  59. Ansbro, S.; Moreno-Pineda, E.; Yu, W.; Ollivier, J.; Mutka, H.; Ruben, M.; Chiesa, A. Magnetic Properties of Transition Metal Dimers Probed by Inelastic Neutron Scattering. Dalton Trans. 2018, 47, 11953–11959. [Google Scholar] [CrossRef]
  60. Garlatti, E.; Guidi, T.; Ansbro, S.; Santini, P.; Amoretti, G.; Ollivier, J.; Mutka, H.; Timco, G.; Vitorica-Yrezabal, I.J.; Whitehead, G.F.S.; et al. Portraying Entanglement between Molecular Qubits with Four-Dimensional Inelastic Neutron Scattering. Nat. Commun. 2017, 8, 14543. [Google Scholar] [CrossRef]
  61. Chakraborty, T.; Singh, H.; Mitra, C. Experimental Evidences of Singlet to Triplet Transition in a Spin Cluster Compound. J. Magn. Magn. Mater. 2015, 396, 247–253. [Google Scholar] [CrossRef]
  62. Das, D.; Singh, H.; Chakraborty, T.; Gopal, R.K.; Mitra, C. Experimental Detection of Quantum Information Sharing and Its Quantification in Quantum Spin Systems. New J. Phys. 2013, 15, 013047. [Google Scholar] [CrossRef]
  63. Reis, M.S.; Soriano, S.; dos Santos, A.M.; Sales, B.C.; Soares-Pinto, D.O.; Brandão, P. Evidence for Entanglement at High Temperatures in an Engineered Molecular Magnet. Europhys. Lett. 2012, 100, 50001. [Google Scholar] [CrossRef]
  64. Sahling, S.; Remenyi, G.; Paulsen, C.; Monceau, P.; Saligrama, V.; Marin, C.; Revcolevschi, A.; Regnault, L.P.; Raymond, S.; Lorenzo, J.E. Experimental Realization of Long-Distance Entanglement between Spins in Antiferromagnetic Quantum Spin Chains. Nat. Phys. 2015, 11, 255–260. [Google Scholar] [CrossRef]
  65. Laurell, P.; Scheie, A.; Mukherjee, C.J.; Koza, M.M.; Enderle, M.; Tylczynski, Z.; Okamoto, S.; Coldea, R.; Tennant, D.A.; Alvarez, G. Quantifying and Controlling Entanglement in the Quantum Magnet Cs2CoCl4. Phys. Rev. Lett. 2021, 127, 037201. [Google Scholar] [CrossRef]
  66. Scheie, A.; Laurell, P.; Samarakoon, A.M.; Lake, B.; Nagler, S.E.; Granroth, G.E.; Okamoto, S.; Alvarez, G.; Tennant, D.A. Witnessing Entanglement in Quantum Magnets Using Neutron Scattering. Phys. Rev. B 2021, 103, 224434. [Google Scholar] [CrossRef]
  67. Haraldsen, J.T. Heisenberg Pentamer: Insights into Inelastic Neutron Scattering on Magnetic Clusters. Phys. Rev. Lett. 2011, 107, 037205. [Google Scholar] [CrossRef] [PubMed]
  68. Haraldsen, J.T. Evolution of Thermodynamic Properties and Inelastic Neutron Scattering Intensities for Spin-1/2 Antiferromagnetic Quantum Rings. Phys. Rev. B 2016, 94, 054436. [Google Scholar] [CrossRef]
  69. Szałowski, K.; Kowalewska, P. Magnetocaloric Effect in Cu5-NIPA Molecular Magnet: A Theoretical Study. Materials 2020, 13, 485. [Google Scholar] [CrossRef] [PubMed]
  70. Torrico, J.; Plascak, J.A. Study of the Ground State and Thermodynamic Properties of Cu5-NIPA-like Molecular Nanomagnets. J. Magn. Magn. Mater. 2022, 552, 169151. [Google Scholar] [CrossRef]
  71. Rajca, A.; Wongsriratanakul, J.; Rajca, S. Organic Spin Clusters: Macrocyclic-Macrocyclic Polyarylmethyl Polyradicals with Very High Spin S = 5–13. J. Am. Chem. Soc. 2004, 126, 6608–6626. [Google Scholar] [CrossRef]
  72. Shatruk, M.; Chambers, K.E.; Prosvirin, A.V.; Dunbar, K.R. Systematic Investigation of Trigonal-Bipyramidal Cyanide-Bridged Clusters of the First-Row Transition Metals. Inorg. Chem. 2007, 46, 5155–5165. [Google Scholar] [CrossRef] [PubMed]
  73. Nath, R.; Tsirlin, A.A.; Khuntia, P.; Janson, O.; Förster, T.; Padmanabhan, M.; Li, J.; Skourski, Y.; Baenitz, M.; Rosner, H.; et al. Magnetization and Spin Dynamics of the Spin S = 1/2 Hourglass Nanomagnet Cu5(OH)2(NIPA)4·10H2O. Phys. Rev. B 2013, 87, 214417. [Google Scholar] [CrossRef]
  74. Werner, T.W.; Reschke, S.; Bunzen, H.; von Nidda, H.A.K.; Deisenhofer, J.; Loidl, A.; Volkmer, D. [Co5Tp*4(Me2bta)6]: A Highly Symmetrical Pentanuclear Kuratowski Complex Featuring Tris(Pyrazolyl)Borate and Benzotriazolate Ligands. Inorg. Chem. 2016, 55, 1053–1060. [Google Scholar] [CrossRef]
  75. Bai, Y.L.; Tao, J.; Huang, R.B.; Zheng, L.S. The Designed Assembly of Augmented Diamond Networks From Predetermined Pentanuclear Tetrahedral Units. Angew. Chem. Int. Ed. 2008, 47, 5344–5347. [Google Scholar] [CrossRef] [PubMed]
  76. Fedulin, A.; Gupta, S.K.; Rüter, I.; Meyer, F.; Jacobi von Wangelin, A. Polynuclear Iron(II) Pyridonates: Synthesis and Reactivity of Fe4 and Fe5 Clusters. Inorg. Chem. 2022, 61, 6149–6159. [Google Scholar] [CrossRef]
  77. Cui, L.; Geng, Y.F.; Leong, C.F.; Ma, Q.; D’Alessandro, D.M.; Deng, K.; Zeng, Q.D.; Zuo, J.L. Synthesis, Properties and Surface Self-Assembly of a Pentanuclear Cluster Based on the New π-Conjugated TTF-triazole Ligand. Sci. Rep. 2016, 6, 25544. [Google Scholar] [CrossRef]
  78. Wang, F.; Chen, Z. Geometric Spin Frustration for Isolated Plaquettes of the Lattices: An Extended Irreducible Tensor Operator Method. Chem. Phys. 2006, 327, 427–433. [Google Scholar] [CrossRef]
  79. Kuratowski, C. Sur le problème des courbes gauches en Topologie. Fundam. Math. 1930, 15, 271–283. [Google Scholar] [CrossRef]
  80. Biswas, S.; Tonigold, M.; Speldrich, M.; Kögerler, P.; Weil, M.; Volkmer, D. Syntheses and Magnetostructural Investigations on Kuratowski-Type Homo- and Heteropentanuclear Coordination Compounds [MZn4Cl4(L)6] (MII = Zn, Fe, Co, Ni, or Cu; L = 5,6-Dimethyl-1,2,3-Benzotriazolate) Represented by the Nonplanar K3,3 Graph. Inorg. Chem. 2010, 49, 7424–7434. [Google Scholar] [CrossRef]
  81. Buvaylo, E.A.; Nesterova, O.V.; Kokozay, V.N.; Vassilyeva, O.Y.; Skelton, B.W.; Boča, R.; Nesterov, D.S. Discussion of Planarity of Molecular Structures Using Novel Pentanuclear Cu/Ni Complexes as an Example. Cryst. Growth Des. 2012, 12, 3200–3208. [Google Scholar] [CrossRef]
  82. Kambe, K. On the Paramagnetic Susceptibilities of Some Polynuclear Complex Salts. J. Phys. Soc. Jpn. 1950, 5, 48–51. [Google Scholar] [CrossRef]
  83. Nielsen, M.A.; Chuang, I.L. Quantum Computation and Quantum Information; Cambridge University Press: Cambridge, UK, 2000. [Google Scholar]
  84. Castellano, G. Thermodynamic Potentials for Simple Magnetic Systems. J. Magn. Magn. Mater. 2003, 260, 146–150. [Google Scholar] [CrossRef]
  85. Plascak, J.A. Ensemble Thermodynamic Potentials of Magnetic Systems. J. Magn. Magn. Mater. 2018, 468, 224–229. [Google Scholar] [CrossRef]
  86. Hill, S.A.; Wootters, W.K. Entanglement of a Pair of Quantum Bits. Phys. Rev. Lett. 1997, 78, 5022–5025. [Google Scholar] [CrossRef]
  87. Wootters, W.K. Entanglement of Formation of an Arbitrary State of Two Qubits. Phys. Rev. Lett. 1998, 80, 2245–2248. [Google Scholar] [CrossRef]
  88. Torrico, J.; Plascak, J.A. Ground State and Thermodynamic Properties of Spin-1/2 Isosceles Heisenberg Triangles for V6-like Magnetic Molecules. Phys. Rev. E 2020, 102, 062116. [Google Scholar] [CrossRef]
  89. Wang, A.-M. A Simplified and Obvious Expression of Concurrence in Wootters’ Measure of Entanglement of a Pair of Qubits. Chin. Phys. Lett. 2003, 20, 1907–1909. [Google Scholar] [CrossRef]
  90. Osenda, O.; Raggio, G.A. Entanglement in Thermal Equilibrium States. Phys. Rev. A 2005, 72, 064102. [Google Scholar] [CrossRef]
  91. Markham, D.; Anders, J.; Vedral, V.; Murao, M.; Miyake, A. Survival of Entanglement in Thermal States. Europhys. Lett. 2008, 81, 40006. [Google Scholar] [CrossRef]
  92. Sadiek, G.; Kais, S. Persistence of Entanglement in Thermal States of Spin Systems. J. Phys. B At. Mol. Opt. Phys. 2013, 46, 245501. [Google Scholar] [CrossRef]
  93. Berry, D.W.; Dowling, M.R. Two Qubits Can Be Entangled in Two Distinct Temperature Regions. Phys. Rev. A 2006, 74, 062301. [Google Scholar] [CrossRef]
  94. Sharma, S.; Kothari, R. Antiferromagnetic Materials Pave the Path for Quantum Entanglement. Mater. Today Proc. 2021, 43, 2428–2431. [Google Scholar] [CrossRef]
  95. Schnack, J. Influence of Intermolecular Interactions on Magnetic Observables. Phys. Rev. B 2016, 93, 054421. [Google Scholar] [CrossRef]
  96. Gálisová, L.; Strečka, J.; Verkholyak, T.; Havadej, S. Magnetization Plateaus and Bipartite Entanglement of an Exactly Solved Spin-1/2 Ising-Heisenberg Orthogonal-Dimer Chain. Phys. E Low-Dimens. Syst. Nanostruct. 2021, 125, 114089. [Google Scholar] [CrossRef]
  97. Gálisová, L. Insight into Ground-State Spin Arrangement and Bipartite Entanglement of the Polymeric Coordination Compound [Dy2Cu2]n through the Symmetric Spin-1/2 Ising–Heisenberg Orthogonal-Dimer Chain. J. Magn. Magn. Mater. 2022, 561, 169721. [Google Scholar] [CrossRef]
Figure 1. Schematic views of the considered spin cluster with exchange integrals of two kinds marked with different lines: J 1 with solid lines and J 2 with dashed lines: (a) a tetrahedral structure for spin tetramer with additional spin in the center; (b) a planar structure for the tetramer with additional spin interacting with all its members. The numbers from 1 to 5 label the spins.
Figure 1. Schematic views of the considered spin cluster with exchange integrals of two kinds marked with different lines: J 1 with solid lines and J 2 with dashed lines: (a) a tetrahedral structure for spin tetramer with additional spin in the center; (b) a planar structure for the tetramer with additional spin interacting with all its members. The numbers from 1 to 5 label the spins.
Molecules 28 06418 g001
Figure 2. Ground-state magnetic phase diagram of considered spin cluster in the plane J 2 / | J 1 | H / | J 1 | , for (a) ferromagnetic J 1 and (b) antiferromagnetic J 1 exchange integral. The (linear) phase boundaries are marked with lines and their equations are listed in the inset. The filled circles denote the triple/quadruple points. The pure ground state is denoted by vector ψ S T , s , whereas the mixed (degenerate) ground state is denoted by density matrix ρ ^ S T , s .
Figure 2. Ground-state magnetic phase diagram of considered spin cluster in the plane J 2 / | J 1 | H / | J 1 | , for (a) ferromagnetic J 1 and (b) antiferromagnetic J 1 exchange integral. The (linear) phase boundaries are marked with lines and their equations are listed in the inset. The filled circles denote the triple/quadruple points. The pure ground state is denoted by vector ψ S T , s , whereas the mixed (degenerate) ground state is denoted by density matrix ρ ^ S T , s .
Molecules 28 06418 g002
Figure 3. Ground-state phase diagram showing the values of Wootters concurrences C ( 1 ) , C ( 2 ) for spin pairs interacting with J 1 and J 2 exchange integral, respectively (see Figure 1), for (a) ferromagnetic J 1 and (b) antiferromagnetic J 1 exchange integral. The color lines mark the phase boundaries as in Figure 2. The values of concurrences exactly along the phase boundaries or at the triple point are shown with arrows. If no arrow points to the phase boundary or a point, both concurrences vanish there.
Figure 3. Ground-state phase diagram showing the values of Wootters concurrences C ( 1 ) , C ( 2 ) for spin pairs interacting with J 1 and J 2 exchange integral, respectively (see Figure 1), for (a) ferromagnetic J 1 and (b) antiferromagnetic J 1 exchange integral. The color lines mark the phase boundaries as in Figure 2. The values of concurrences exactly along the phase boundaries or at the triple point are shown with arrows. If no arrow points to the phase boundary or a point, both concurrences vanish there.
Molecules 28 06418 g003
Figure 4. Contour plot of C ( 2 ) concurrence value for J 1 > 0 : (a) at the finite normalized temperature of k B T / | J 1 | = 0.2, in the plane J 2 / | J 1 | H / | J 1 | ; (b) for J 2 / | J 1 | = −1, in the plane k B T / | J 1 | H / | J 1 | . Color bar right to each panel shows the concurrence values for each contour. The white range corresponds to separable states.
Figure 4. Contour plot of C ( 2 ) concurrence value for J 1 > 0 : (a) at the finite normalized temperature of k B T / | J 1 | = 0.2, in the plane J 2 / | J 1 | H / | J 1 | ; (b) for J 2 / | J 1 | = −1, in the plane k B T / | J 1 | H / | J 1 | . Color bar right to each panel shows the concurrence values for each contour. The white range corresponds to separable states.
Molecules 28 06418 g004
Figure 5. The boundary between the entangled and the separable state of type 2 spin pairs in the plane k B T / | J 1 | H / | J 1 | , for J 1 > 0 and for various values of J 2 / | J 1 | .
Figure 5. The boundary between the entangled and the separable state of type 2 spin pairs in the plane k B T / | J 1 | H / | J 1 | , for J 1 > 0 and for various values of J 2 / | J 1 | .
Molecules 28 06418 g005
Figure 6. Dependence of the concurrence values C ( 2 ) (a) on the normalized temperature, for various values of normalized magnetic field and (b) on the normalized magnetic field, for various values of normalized temperature, for J 1 > 0 and J 2 / | J 1 | = −1.
Figure 6. Dependence of the concurrence values C ( 2 ) (a) on the normalized temperature, for various values of normalized magnetic field and (b) on the normalized magnetic field, for various values of normalized temperature, for J 1 > 0 and J 2 / | J 1 | = −1.
Molecules 28 06418 g006
Figure 7. Contour plot of (a) C ( 1 ) and (b) C ( 2 ) concurrence value at the finite normalized temperature of k B T / | J 1 | = 0.2, in the plane J 2 / | J 1 | H / | J 1 | , for J 1 < 0 . Color bar right to each panel shows the concurrence values for each contour. The white range corresponds to separable states.
Figure 7. Contour plot of (a) C ( 1 ) and (b) C ( 2 ) concurrence value at the finite normalized temperature of k B T / | J 1 | = 0.2, in the plane J 2 / | J 1 | H / | J 1 | , for J 1 < 0 . Color bar right to each panel shows the concurrence values for each contour. The white range corresponds to separable states.
Molecules 28 06418 g007
Figure 8. Contour plot of (a) C ( 1 ) and (b) C ( 2 ) concurrence value for J 2 / | J 1 | = 1 and J 1 < 0 , in the plane k B T / | J 1 | H / | J 1 | . Color bar right to each panel shows the concurrence values for each contour. The white range corresponds to separable states.
Figure 8. Contour plot of (a) C ( 1 ) and (b) C ( 2 ) concurrence value for J 2 / | J 1 | = 1 and J 1 < 0 , in the plane k B T / | J 1 | H / | J 1 | . Color bar right to each panel shows the concurrence values for each contour. The white range corresponds to separable states.
Molecules 28 06418 g008
Figure 9. The boundary between the entangled and the separable state of type 1 spin pairs (a) and type 2 spin pairs (b) in the plane k B T / | J 1 | H / | J 1 | , for J 1 < 0 and for various values of J 2 / | J 1 | .
Figure 9. The boundary between the entangled and the separable state of type 1 spin pairs (a) and type 2 spin pairs (b) in the plane k B T / | J 1 | H / | J 1 | , for J 1 < 0 and for various values of J 2 / | J 1 | .
Molecules 28 06418 g009
Figure 10. Dependence of the concurrence values C ( 1 ) [for (a,c)] and C ( 2 ) [for (b,d)] on the normalized temperature, for various values of normalized magnetic field, for J 1 < 0 and for J 2 / | J 1 | = 1 [(a,b)] or for J 2 / | J 1 | = −0.5 [(c,d)].
Figure 10. Dependence of the concurrence values C ( 1 ) [for (a,c)] and C ( 2 ) [for (b,d)] on the normalized temperature, for various values of normalized magnetic field, for J 1 < 0 and for J 2 / | J 1 | = 1 [(a,b)] or for J 2 / | J 1 | = −0.5 [(c,d)].
Molecules 28 06418 g010
Figure 11. Dependence of the concurrence values C ( 1 ) [for (a,c)] and C ( 2 ) [for (b,d)] on the normalized magnetic field, for various values of normalized temperature, for J 1 < 0 and for J 2 / | J 1 | = 1 [(a,b)] or for J 2 / | J 1 | = −0.5 [(c,d)].
Figure 11. Dependence of the concurrence values C ( 1 ) [for (a,c)] and C ( 2 ) [for (b,d)] on the normalized magnetic field, for various values of normalized temperature, for J 1 < 0 and for J 2 / | J 1 | = 1 [(a,b)] or for J 2 / | J 1 | = −0.5 [(c,d)].
Molecules 28 06418 g011
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Szałowski, K. Phase Diagram and Quantum Entanglement Properties of a Pentamer S = 1/2 Heisenberg Spin Cluster. Molecules 2023, 28, 6418. https://doi.org/10.3390/molecules28176418

AMA Style

Szałowski K. Phase Diagram and Quantum Entanglement Properties of a Pentamer S = 1/2 Heisenberg Spin Cluster. Molecules. 2023; 28(17):6418. https://doi.org/10.3390/molecules28176418

Chicago/Turabian Style

Szałowski, Karol. 2023. "Phase Diagram and Quantum Entanglement Properties of a Pentamer S = 1/2 Heisenberg Spin Cluster" Molecules 28, no. 17: 6418. https://doi.org/10.3390/molecules28176418

APA Style

Szałowski, K. (2023). Phase Diagram and Quantum Entanglement Properties of a Pentamer S = 1/2 Heisenberg Spin Cluster. Molecules, 28(17), 6418. https://doi.org/10.3390/molecules28176418

Article Metrics

Back to TopTop