1. Introduction
The contemporary development of novel quantum technologies requires expertise from various scientific disciplines spanning quantum physics [
1], materials engineering [
2] to quantum information science [
3,
4]. One of the most challenging tasks in this rapidly evolving and highly interdisciplinary research area is to identify a suitable platform for realizing quantum computers [
5]. Although several alternative paradigms are currently being explored to develop novel quantum technologies using diverse physical platforms such as photons, trapped atoms, quantum dots, nuclear and electron spins, or superconducting circuits [
6,
7], none has yet satisfactorily addressed all fundamental issues related to a practical realization of quantum computation and quantum information processing [
5]. The main obstacles in overcoming this intricate problem are associated with the inherent nature of quantum materials, which typically retain their genuine quantum features only at relatively low temperatures or for relatively short coherence times.
Molecular magnetic materials are compelling candidates for encoding qubits due to their collective spin degrees of freedom, which arise from exchange-coupled electron spins inherent in the individual spin carriers they are composed of [
8,
9,
10,
11]. The main advantage of molecular qubits lies in the capability of targeted design of magnetic molecules with well-controlled number of magnetic centers and interactions between them through a rational bottom-up chemical synthesis [
12,
13]. Another important benefit of qubits encoded through molecular magnets is the relatively effective decoupling of their electronic degrees of freedom from the environment. This feature allows for tuning their coherence time to be sufficiently long for implementation in quantum computing and quantum information processing [
14,
15,
16,
17]. It has been suggested that molecular materials with a carefully selected magnetic core may provide a suitable platform for operating as a storage unit of a dynamic random access memory device [
18], single- and multiqubit quantum gates [
19,
20], or a quantum battery capable of storing a finite amount of extractable work [
21].
A smart choice of molecular magnetic materials with sufficiently strong exchange coupling constants additionally provides the opportunity to suppress the loss of their genuine quantum features due to thermal effects. The bipartite thermal entanglement may indeed persist up relatively high temperatures as for instance 40 K reported for the spin-1 Ni
dinuclear complex [Ni
(Medpt)
(
-ox)(H
O)
](ClO
)
[
22], 50 K for the spin-1 Ni
tetranuclear complex [Ni
(
-CO
)
(aetpy)
](ClO
)
[
23], 200 K for the spin-
Cu
coordination polymer Na
Cu
Si
O
[
24], 680 K for the spin-
Cu
dinuclear complex [Cu
(
-HCOO)
(HCOO)
]C
H
N
[
25] or even 730 K for the spin-
Fe
dinuclear complex [Fe
(
-oxo)(C
H
N
)
(C
O
)
] [
26]. However, none of the aforementioned molecular magnetic materials can be regarded as an experimental realization of the molecular qubit, because they suffer from a limitation of their singlet ground state.
One possible route towards a physical realization of the molecular qubit is offered by heterobimetallic dinuclear complexes, which represents paradigmatic realizations of the mixed spin-(
) Heisenberg dimer with an antiferromagnetic exchange interaction. For this reason, the mixed spin-(
) Heisenberg dimer has recently attracted considerable attention as the simplest member of this class of quantum spin models [
27,
28,
29,
30,
31,
32,
33,
34,
35]. Among other matters, it has been evidenced that the mixed spin-(
) Heisenberg dimer has a higher threshold temperature for thermal entanglement compared to the pure spin-
Heisenberg dimer [
27], the strength of thermal entanglement is significantly influenced by a nonuniform magnetic field resulting from different g-factors of individual metal centers [
28,
29], the Zeeman’s splitting of energy levels may substantially enhance the strength of thermal entanglement at low magnetic fields [
30,
31,
32], the electric field can serve as another driving force controlling the strength of thermal entanglement [
33], and the thermal entanglement can be quite robust against decoherence [
34,
35]. From a materials science perspective, the thermal entanglement of experimental realizations of the antiferromagnetic mixed spin-(
) Heisenberg dimer has been studied only for the particular case of the heterodinuclear complex [Ni(dpt)(H
O)Cu(pba)] · 2H
O [
36], which exhibits thermal entanglement up to a threshold temperature of about 140 K comparable to the magnitude of the exchange coupling constant between the spin-
Cu
and spin-1 Ni
magnetic ions [
30,
32].
In the present paper we will comprehensively examine the strength of bipartite entanglement in the mononuclear molecular complex (Et
NH)[Ni(hfac)
L] (HL denotes 2-(2-hydroxy-3-methoxy-5-nitrophenyl)-4,4,5,5-tetramethyl-4,5-dihydro-1H-imidazol-3-oxide-1-oxyl and hfacH stands for hexafluoroacetylacetone) [
37], which provides completely different platform for an experimental realization of the mixed spin-(
) Heisenberg dimer. Unlike the previous case involving the molecular magnet composed from two distinct metal centers, the magnetic structure of the mononuclear molecular compound (Et
NH)[Ni(hfac)
L] consists of the spin-
nitronyl-nitroxide radical substituted nitrophenol exchange coupled to the spin-1 Ni
magnetic ion (see
Figure 1). It will be demonstrated hereafter that this paradigmatic realization of the mixed spin-(
) Heisenberg dimer may showcase significant thermal entanglement beyond room temperature owing to the sizable exchange interaction
J/
= 505 K (
is Boltzmann’s constant) between the spin-
nitronyl-nitroxide radical and the spin-1 Ni
magnetic ion.
The organization of this paper is as follows. In the
Section 2 we will adapt the mixed spin-(
) Heisenberg dimer for theoretical modeling of the molecular complex (Et
NH)[Ni(hfac)
L and describe basic steps of the procedure used for calculating the negativity serving as a measure of the bipartite entanglement. The most interesting results for the quantum and thermal entanglement of the magnetic molecule (Et
NH)[Ni(hfac)
L are then comprehensively discussed in
Section 3. Finally, a brief summary of the most important findings and future outlooks are provided in
Section 4.
2. Model and Method
In this section, let us first introduce the respective theoretical model of the molecular complex (Et
NH)[Ni(hfac)
L and subsequently describe the basic steps of the procedure allowing calculation of the strength of bipartite thermal entanglement. Taking into consideration an isotropic nature of the exchange coupling between the spin-1 Ni
magnetic ion and the spin-
nitronyl-nitroxide radical, the molecular complex (Et
NH)[Ni(hfac)
L can be theoretically modeled by the mixed spin-(
,1) Heisenberg dimer defined by the Hamiltonian:
where
and
denote the spin-
and spin-1 operators ascribed to the nitronyl-nitroxide radical and Ni
magnetic ion, respectively. The coupling constant
J determines the isotropic exchange interaction between the spin-
nitronyl-nitroxide radical and the spin-1 Ni
magnetic ion, the parameter
B denotes a static external magnetic field,
is Bohr magneton, while
and
are Landé
g-factors of the nitronyl-nitroxide radical and Ni
magnetic ion, respectively.
The Hamiltonian (
1) acquires the following matrix form in the standard basis spanned over the complete set of eigenvectors
of
z-components of the spin-
and spin-1 operators:
For brevity, we have adapted the following notation
and
for ’local’ Zeeman’s terms acting on the spin-
and spin-1 magnetic entities, which may be different due to difference of the gyromagnetic g-factors the nitronyl-nitroxide radical and Ni
magnetic ion
. An exact analytical diagonalization of the Hamiltonian matrix (
2) affords the following set of eigenvalues:
and the respective eigenvectors:
The eigenvectors () and () are expressed in terms of the probability amplitudes:
The strength of quantum and thermal entanglement in pure and mixed states of the mixed spin-(
,1) Heisenberg dimer can be quantified through the quantity
negativity, which was introduced by Vidal and Werner [
38] based on the famous Peres-Horodecki separability criterion [
39,
40]:
The negativity is defined through the eigenvalues
of a partially transposed density matrix
derived from the density matrix
upon partial transposition
with respect to the spin-
subsystem. The negativity is zero for separable states but becomes non-zero for entangled states. It is quite evident from equation (
10) that the bipartite entanglement exists if at least one eigenvalue
of the partially transposed density matrix
becomes negative [
39,
40].
A spectral decomposition, which implements the eigenvalues (
3)-() and eigenvectors (
6)-() of the mixed spin-(1/2,1) Heisenberg dimer, provides a useful tool for calculation of the density operator
:
whereby the partition function
acquires the following explicit form:
The matrix representation of the density operator (
11) gives the following density matrix:
which is defined through the elements
explicitly quoted in
Appendix A. The partially transposed density matrix
obtained after a partial transposition
with respect to the spin-
subsystem reads:
After a straightforward diagonalization of the partially transposed density matrix (
14) one obtains the following spectrum of eigenvalues:
The eigenvalues (
15) and () with a negative sign before the square root may become negative, which implies according to the formula (
10) the possible existence of bipartite entanglement (nonzero negativity). In the next section we will adapt the theoretical results derived for the mixed spin-(
,1) Heisenberg dimer in order to provide a theoretical modeling of bipartite entanglement emerging within the nickel-radical molecular complex (Et
NH)[Ni(hfac)
L].
3. Results and discussions
Before proceeding to a discussion of the most interesting results, let us first briefly review experimental findings reported previously for the nickel-radical molecular complex (Et
NH)[Ni(hfac)
L] in reference [
37]. Electron paramagnetic resonance spectra revealed two different signals characterized through g-factors
and
, which can be attributed to the spin-
nitronyl-nitroxide radical and the spin-1 Ni
magnetic ion, respectively. It turned out, moreover, that both Landé g-factors are almost completely independent of temperature with only negligible changes of the relevant g-factors on the third decimal place. Besides, theoretical modeling of the temperature dependence of magnetic susceptibility and the low-temperature magnetization curve unveiled a substantial antiferromagnetic exchange coupling
K between the spin-
nitronyl-nitroxide radical and the spin-1 Ni
magnetic ion. The paramount finding of significant antiferromagnetic exchange coupling thus establishes the mononuclear complex (Et
NH)[Ni(hfac)
L] as a compelling candidate in magnetic materials research capable of showcasing bipartite entanglement persisting well above room temperature. In the following part, we will adjust the parameters of the exactly solved mixed spin-(
,1) Heisenberg dimer given by the Hamiltonian (
1) to the previously mentioned values
K,
, and
, which allows us to conduct theoretical modeling of the bipartite entanglement in the molecular compound (Et
NH)[Ni(hfac)
L].
3.1. Quantum entanglement in the nickel-radical molecular complex (EtNH)[Ni(hfac)L]
First, let us examine the strength of bipartite quantum entanglement in the molecular complex (Et
NH)[Ni(hfac)
L] at absolute zero temperature. To achieve this goal, we need to determine the ground state of the mixed spin-(
,1) Heisenberg dimer and construct the relevant density operator. In the absence of an external magnetic field, the ground state of the mixed spin-(
,1) Heisenberg dimer with an isotropic exchange interaction consists of a two-fold degenerate quantum ferrimagnetic state represented by the eigenvectors:
The density operator
ascribed to the ferrimagnetic doublet (
18) implies a moderately strong zero-temperature value of the negativity
, which acquires according to equation (
10) two-thirds of its maximal value
corresponding to perfectly entangled Bell states. It should be pointed out that the sizable energy gap
K between a quantum ferrimagnetic ground-state doublet and a spin-quartet excited state implies that the molecular complex (Et
NH)[Ni(hfac)
L] can be regarded, up to extremely high temperatures and magnetic fields, as a compelling experimental realization of the molecular qubit.
An even more intriguing situation emerges when the external magnetic field lifts the degeneracy of the quantum ferrimagnetic doublet (
18) through Zeeman’s splitting of energy levels. The Zeeman effect consequently leads to a nondegenerate quantum ferrimagnetic ground state given by the eigenvector:
which is characterized through the probability amplitudes depending on a mutual interplay of the exchange-coupling constant
J and the magnetic-field term
involving the difference of the relevant g-factors:
The assessment of bipartite quantum entanglement within the nondegenerate ferrimagnetic ground state
given by equations (
19)-(
20) can be readily accomplished with the help of the respective density operator
being consistent with the following zero-temperature value of the negativity:
A low-field asymptotic value of the negativity
calculated according to equation (
21) for the nondegenerate ferrimagnetic ground state
in the limit
is indicative of an unexpected enhancement of the bipartite quantum entanglement when compared to the zero-field value
corresponding to the two-fold degenerate ferrimagnetic ground states
emergent strictly at zero magnetic field
. Note that this phenomenon contradicts the conventional expectation of magnetic-field-induced suppression of the bipartite quantum entanglement, because it suggests instead its significant rise (enhancement approximately by 41 %) at nonzero magnetic fields due to the Zeeman effect. Although a subsequent rise of the magnetic field causes a gradual reduction in the strength of bipartite quantum entanglement, this effect is quite weak for achievable magnetic-field strengths due to the negligible size of the magnetic-field-dependent term
relative to the sizable exchange interaction
K. It will be demonstrated hereafter that the sudden increase in bipartite quantum entanglement due to the magnetic field persists even at low temperatures what underscores its potential significance for experimental verification.
3.2. Thermal entanglement in nickel-radical molecular complex (EtNH)[Ni(hfac)L]
Next, let us investigate the strength of bipartite thermal entanglement in the molecular complex (Et
NH)[Ni(hfac)
L] at nonzero temperatures. Typical temperature dependencies of the negativity of the mixed spin-(
,1) Heisenberg dimer with the parameter set
K,
, and
adjusted to a theoretical modeling of the molecular compound (Et
NH)[Ni(hfac)
L] are depicted in
Figure 2 for a few selected values of the external magnetic field. In agreement with the previous ground-state analysis, the negativity starts at zero magnetic field from the zero-temperature asymptotic value
and then it shows a gradual temperature-induced decline until the negativity completely vanishes at the threshold temperature
K. Contrary to this, it is evident from
Figure 2 and its inset that the negativity starts at relatively small magnetic fields (
T) from the higher zero-temperature asymptotic value
in accordance with the previous arguments concerning with the magnetic-field-driven enhancement of bipartite quantum entanglement. The negativity at relatively small magnetic fields subsequently shows a sudden drop upon increasing of temperature until it tends towards the respective zero-field value before the negativity finally displays a gradual temperature-induced decline terminating at the same threshold temperature
K as the zero-field case. Note that the initial thermally-driven reduction of the negativity is the steeper, the smaller the external magnetic field is, whereas this phenomenon gradually vanishes at higher magnetic fields. A detectable reduction of zero-temperature value of the negativity
due to the external magnetic field occurs merely at extremely high magnetic fields
T.
Furthermore, a few typical isothermal dependencies of the negativity on a magnetic field are plotted in
Figure 3 for the mixed spin-(
,1) Heisenberg dimer with the parameter set
K,
, and
adjusted to a theoretical modeling of the molecular compound (Et
NH)[Ni(hfac)
L]. At low enough temperatures
K the negativity shows an abrupt magnetic-field-induced rise from the zero-field value
towards the much higher value
, which is subsequently followed by a gradual quasi-linear decrease achieved upon further strengthening of the magnetic field until the negativity suddenly drops down to zero at the extremely high magnetic field
T. While a small magnetic field significantly strengthens the thermal entanglement roughly by 41% due to the Zeeman splitting of the ferrimagnetic ground-state doublet, the subsequent reduction of thermal entanglement observable upon further increase of the magnetic field is quite negligible. Namely, the magnetic-field change of 500 T decreases the negativity only by approximately 2%. It is noteworthy that the unconventional magnetic-field-driven strengthening of thermal entanglement is gradually suppressed upon increasing of temperature as exemplified by the isothermal dependencies of the negativity calculated at temperatures
K and 200 K. At and above room temperature, the negativity displays a conventional monotonous decrease upon increasing the magnetic field, whereby the strength of thermal entanglement is kept almost constant in the magnetic-field range
T effective for technological applications. From this perspective, it is worthwhile to remark that the strength of bipartite entanglement of the molecular compound (Et
NH)[Ni(hfac)
L] at room temperature acquires roughly 40% of its maximum value corresponding to the fully entangled Bell states independently of the magnetic-field strengths accessible for contemporary real-world technologies (
T).
Finally, let us conclude our discussion with a detailed analysis of a density plot of the negativity in the temperature versus magnetic field plane, which was obtained for the mixed spin-(
,1) Heisenberg dimer with the parameter set
K,
,
adjusted for a theoretical modeling of the molecular complex (Et
NH)[Ni(hfac)
L].
Figure 4(a) shows a detail of the relevant density plot up to the magnetic field
T, which represents a typical upper limit of the magnetic-field strength for commercially accessible devices generating static magnetic fields. It is obvious from a low-temperature part of the density plot that this magnetic-field range is capable of stabilizing almost perfect bipartite entanglement
, which nearly approaches fully entangled Bell states, only up to temperature
K. In fact, the rise in temperature necessitates a stronger magnetic field to restore almost perfect bipartite entanglement as reflected by a relatively steep gradient of the dark region in the given density plot. A color shading with character of vertical stripes observed in a high-temperature part of the density plot (
K) contrarily verifies the independence of the strength of bipartite thermal entanglement on the magnetic field in the range of sufficiently small magnetic fields
T. Although the strengthening of bipartite thermal entanglement due to the Zeeman effect cannot be utilized in the high-temperature region, the increasing magnetic field does not diminish the strength of bipartite thermal entanglement either. For completeness, the density plot of the negativity in a full range of temperatures and magnetic fields is plotted in
Figure 4(b) in order to provide the phase diagram, which delimits a parameter region with nonzero thermal entanglement. The displayed phase diagram convincingly evidences that the molecular compound (Et
NH)[Ni(hfac)
L] exhibits the thermal entanglement, which is confined to temperatures below 540 K and magnetic fields below 540 T. The respective limit set for temperature evidently implies a possible implementation of the nickel-radical molecular complex (Et
NH)[Ni(hfac)
L] for room-temperature technological applications exploiting the intense thermal entanglement for quantum computation and quantum information processing.
4. Conclusion
In this paper, we have adapted exact results derived for the negativity of the mixed spin-(
,1) Heisenberg dimer to gain insight into the strength of bipartite quantum and thermal entanglement of the mononuclear molecular complex (Et
NH)[Ni(hfac)
L], which is characterized by a sizable exchange coupling constant between the spin-
nitronyl-nitroxide radical substituted nitrophenol and the spin-1 Ni
magnetic ion [
37]. The significant energy gap
K between the quantum ferrimagnetic ground-state doublet and a quartet excited state suggests that the molecular complex (Et
NH)[Ni(hfac)
L] provides a suitable platform of a molecular qubit in a wide range of temperature and magnetic fields. It has been demonstrated that the Zeeman splitting of the two-fold degenerate quantum ferrimagnetic ground state significantly reinforces the strength of bipartite entanglement at small magnitudes of the magnetic field, whereas a further enhancement of the magnetic field gradually suppresses the strength of bipartite entanglement. The latter reduction of the bipartite entanglement is however quite weak effect for achievable magnetic-field strengths, because it is inhibited by a competition of the extremely strong exchange coupling constant
K relative to the magnetic-field-dependent term
involving a rather small difference between g-factors g
=2.005 and g
=2.275 of the spin-
nitronyl-nitroxide radical and the spin-1 Ni
magnetic ion, respectively. It could be thus concluded that the nickel-radical molecular complex (Et
NH)[Ni(hfac)
L] affords a prominent example of the molecular qubit whose entanglement strength can be substantially enhanced (cca. 41%) by a small driving magnetic field. This significant strengthening of bipartite entanglement in response to the application of small magnetic fields
T indeed persists in the molecular complex (Et
NH)[Ni(hfac)
L] up to temperatures
K highlighting the potential significance of this magnetic compound for experimental verification of this remarkable phenomenon predicted for experimental realizations of the mixed spin-
Heisenberg dimers [
30,
31].
Finally, it is worthwhile to remark that the molecular compound (EtNH)[Ni(hfac)L] maintains even at room temperature a sufficiently strong bipartite entanglement retaining approximately 40% of the maximum value corresponding to perfectly entangled Bell states. This intriguing feature persists even at extremely high magnetic fields considerably exceeding commercially available static magnetic fields. From this perspective, the molecular complex (EtNH)[Ni(hfac)L] can be regarded as a compelling example of a molecular qubit, which may be implemented for prospective technological applications exploiting quantum computing and quantum information processing at room temperature.
Author Contributions
Conceptualization, J.S.; methodology, J.S.; software, J.S. and E.S.S.; validation, J.S., and E.S.S.; formal analysis, J.S., and E.S.S.; investigation, J.S., and E.S.S.; resources, J.S.; data curation, J.S. and E.S.S.; writing—original draft preparation, J.S. and E.S.S.; writing—review and editing, J.S. and E.S.S.; visualization, J.S. and E.S.S.; supervision, J.S.; project administration, J.S.; funding acquisition, J.S. All authors have read and agreed to the published version of the manuscript.