Numerical Interchain Mean-Field Theory for the Specific Heat of the Bimetallic Ferromagnetically Coupled Chain Compound MnNi(NO2)4(en)2 (en = Ethylenediamine).

We present a detailed study of the field-dependent specific heat of the bimetallic ferromagnetically coupled chain compound MnNi(NO2)4(en)2, en = ethylenediamine. For this material, which in zero field orders antiferromagnetically below TN=2.45 K, small fields suppress magnetic order. Instead, in such fields, a double-peak-like structure in the temperature dependence of the specific heat is observed. We attribute this behavior to the existence of an acoustic and an optical mode in the spin-wave dispersion as a result of the existence of two different spins per unit cell. We compare our experimental data to numerical results for the specific heat obtained by exact diagonalization and Quantum Monte Carlo simulations for the alternating spin-chain model, using parameters that have been derived from the high-temperature behavior of the magnetic susceptibility. The interchain coupling is included in the numerical treatment at the mean-field level. We observe remarkable agreement between experiment and theory, including the ordering transition, using previously determined parameters. Furthermore, the observed strong effect of an applied magnetic field on the ordered state of MnNi(NO2)4(en)2 promises interesting magnetocaloric properties.

1. Introduction

Alternation in spin systems, be it of the magnetic coupling, the local symmetry, or the spin value, induces new and exotic types of magnetic ground states and excitations [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18]. In particular, this is exemplified in novel bimetallic chain systems, viz. molecule-based chain systems with alternately arranged magnetic units carrying quantum spins and of different sizes. The ability to synthesize mixed-spin-chain materials [9,19,20,21,22,23,24,25,26,27] has stimulated theoretical investigations [10,11,12,13,14,15,16,17,18,28,29,30,31,32,33,34,35]. The magnon dispersion relation of such chains splits into an optical and an acoustic mode because of the two differently sized quantum spins and per unit cell, both for antiferromagnetic coupling along the chain [10,11,12,13,14,15,16,17,35] as well as ferromagnetic coupling [18,28,29,32,33,34]. Although ground state and fundamental excitations of a Heisenberg ferromagnet are simple, thermodynamic properties are very sensitive to interactions of the magnon excitations, as is evidenced by the ferromagnetic uniform spin-1/2 Heisenberg chain (compare chapter 11.3 of [36] and references therein). Computations for bimetallic Heisenberg chains show that the two energy scales associated to the acoustic and the optical spin excitation modes are reflected by “double-peak” kind of features in the specific heat for both antiferro- [10,11,12,13,14,15,16,17] and ferromagnetic [18,28,29] coupling.

Similar predictions have been made for chains of mixed classical and quantum spins as far back as 1975 [37]. In spite of this long history, experimental verification of the features expected in the specific heat is lacking. In fact, experimental tests of mixed-spin-chain models are scarce [38,39,40], since most materials available contain elements with larger spins [20,21,22,41,42,43,44], which are difficult to be treated adequately in theoretical calculations [10,11,12,13,14,15,16,17].

Here, we present a verification of the two energy scale prediction via a detailed study of the specific heat of MnNi(NO)(en), en = ethylenediamine = CNH, in zero and applied fields. After the field-induced suppression of long-range antiferromagnetic order, we observe a double-peak-like structure in the temperature dependence of for MnNi(NO)(en). We compare our findings with the results of numerical calculations for an , mixed spin chain in zero and external fields. We demonstrate that the in-field calculations, for which finite-size effects are negligible, fully reproduce the double-peak structure of the experimentally observed in-field specific heat. This shows that the optical and acoustic spin excitation mode are reflected by the thermodynamics of this bimetallic chain system. Quantum Monte Carlo (QMC) simulations of the individual chains augmented by a self-consistent mean-field treatment of interchain coupling even yields a remarkably accurate description of the ordering transition in a vanishing magnetic field.

The remainder of this manuscript is organized as follows: Section 2 presents more details on MnNi(NO)(en) and, in particular, a measurement of its specific heat. We then proceed in Section 3 with a detailed theoretical analysis based on exact diagonalization and QMC simulations combined with a mean-field treatment of interchain coupling; some complementary details are provided in Appendix B. In Section 4, we briefly comment on the magnetocaloric properties of MnNi(NO)(en) before we summarize our findings in Section 5. Appendix C contains a summary of a complementary single-site mean-field treatment.

2. Experiment

2.1. MnNi(NO)(en) (en = Ethylenediamine)

MnNi(NO)(en) is one of the best characterized mixed spin chain compounds [45,46,47,48], crystallizing in an orthorhombic structure, space group Pccn (lattice parameters Å, Å, Å). It contains chains of alternately arranged Ni and Mn ions linked by NO ligands, which carry magnetic moments with spin and , respectively (Figure 1). The magnetic coupling along the chain, J, is ferromagnetic [47], with K [28] (we will use units such that throughout.). A finite ionic zero-field splitting D of 0.36 K is derived from the anisotropy of the susceptibility. Due to an effective antiferromagnetic interchain coupling of K, the system undergoes a transition into an antiferromagnetically (AFM) ordered state below K in zero magnetic field and at ambient pressure [28,47,48]. The long-range magnetically ordered state is suppressed by rather small magnetic fields [47].

Figure 1

Two chains in the a-b plane of MnNi(NO)(en), based on the crystal structure of [45]. The thick arrows on the Mn and Ni atoms show the zero-field ordered state: ferromagnetic along the chains and antiferromagnetic between chains.

2.2. Specific Heat

For our study, we used single crystals of MnNi(NO)(en) investigated previously [47] that were grown by slow evaporation, as described in detail in Ref. [46]. Here, we present the easy-axis data , for which AFM ordering is suppressed in T. The heat capacity was measured using commercial calorimeters in magnetic fields up to 1.6 T at temperatures T down to 0.4 K. As will be discussed below, these c axis data allow a comparison to more accurate numerical calculations than the data .

In Figure 2a, we depict the zero-field specific heat of MnNi(NO)(en) as a function of T. The AFM anomaly at K is clearly discernible. To derive the magnetic specific heat, we determine the lattice contribution . Since a single -term does not reproduce the experimental data above , we use two Debye contributions, each calculated via the full Debye integral, to parametrize . MnNi(NO)(en) is built up by chain segments -Mn-NO-Ni-NO-, with two ethylenediamine molecules and two NO groups attached to the Mn and Ni ions, respectively (Figure 1). Intramolecular oscillations of ethylenediamine or NO, because of the light atoms involved, yield Einstein contributions, which are irrelevant for the temperatures considered here. The chain segment units Mn, Ni, and NO are similar in atomic weight. Therefore, to parametrize the lattice contribution of these units we choose one Debye temperature with modes. Analogously, the four attached molecules ethylenediamine and NO per chain segment are parametrized by a second Debye temperature contributing with 12 modes. This way, we reproduce the lattice specific heat of MnNi(NO)(en) with Debye temperatures K and K (solid line in Figure 2a).

Figure 2

(a) Zero-field specific heat of MnNi(NO)(en) as function of temperature T. (b) Zero-field magnetic specific heat and associated entropy S per mole of MnNi(NO)(en) as function of temperature.

We obtain the magnetic specific heat contribution by subtracting from the total (Figure 2b). Further, by numerically integrating , we obtain the magnetic entropy S included in Figure 2b. Both quantities indicate that above there are magnetic fluctuations present over a wide temperature range. In , there is a broad anomaly ranging up to . The associated entropy reaches only at , which is less than half of the value expected for the sum of the magnetic entropies of Ni () and Mn (), (dotted line in Figure 2b). This value is reached only at 10 . Note that the saturation of S at demonstrates the consistency and adequacy of our derivation of the lattice specific heat.

AFM order in MnNi(NO)(en) is suppressed by small magnetic fields [47]. This enables us to study magnetic fluctuations in MnNi(NO)(en), as they appear in . In Figure 3, we plot as function of field. We observe a rapid suppression of the AFM state, in agreement with Ref. [47]. Moreover, after suppression of the AFM state, the broad specific heat anomaly above becomes much more pronounced in magnetic fields and is clearly visible already in the non-phonon corrected data.

Figure 3

(a) Field dependence of of MnNi(NO)(en) for fields . (b) The same data plotted as . Inset: The magnetic phase diagram of MnNi(NO)(en) for : from Ref. [47] (🟉), from the maximum in (⊕), from the maximum in (⋄); lines are guides to the eye.

The temperature of the maximum in represents a measure for an energy scale characteristic for the magnetic fluctuation spectrum (indicated for the 1.6 T data in Figure 3). In the inset of Figure 3, we record its field dependence up to 1.6 T, with a modest increase in of about 1 K/T. Further, after suppression of AFM order in the T dependence of there is additional structure. This is most clearly seen for , where one now observes a double-peak-like structure (see Figure 3b). We take as measure for a second characteristic energy scale the maximum in and include its field dependence in Figure 3. Again, we find a modest increase of by about 1 K between 0.4 and 1.6 T.

and are clearly distinct temperatures and increase at a similar rate. Therefore, they do not stem from ionic states Zeeman split in an external field. Further, extrapolating to zero field yields a finite value of about 0.7 K, implying that does not arise from Zeeman splitting of ionic degenerate states. Therefore, we associate both characteristic energy scales and with collective excitation modes of the magnetic fluctuation spectrum of MnNi(NO)(en) as the result of the existence of an acoustic and an optical magnon mode.

3. Theory

We now proceed to provide a theoretical description of the experimental findings.

3.1. Model

We start from the basic chain model where the () correspond to the spins of the Ni ions (Mn ions) and have (). Following Refs. [28,47], we take a single-ion anisotropy into account only for the Mn sites. The main role of this anisotropy is to select a preferred axis, it should not matter too much if this is due to the Mn or the Ni sites, and it is the form Equation (1) for which parameters were extracted in Ref. [28] by analyzing the high-temperature behavior of the magnetic susceptibility. Nevertheless, we refer to Appendix A for a discussion of the one-magnon dispersion for the case where both anisotropies are present. In the following discussion, we will use the parameters that have been determined in Ref. [28]—namely, K and K—or in units with : . In the latter units, and assuming magnetic g factors , the magnetic fields of T and T shown in Figure 3 are modeled by and , respectively.

3.2. Numerical Treatment of Decoupled Chains

Previously, some of the present authors have performed exact (full) diagonalization and Quantum Monte Carlo (QMC) simulations of chains with , [28,29]. The previous exact diagonalization (ED) investigations went to spins with and . When we replace a spin by , the local Hilbert space dimension increases from 2 to 6, i.e., by a factor 3. Thus, here, we have to contend ourselves with ED for chains with . Adding one unit cell would increase the total Hilbert space dimension by a factor of 18 for the case , such that the next system size remains out of reach. We use conservation of as well as spatial symmetries. The magnetic susceptibility and specific heat c can then be calculated from the eigenvalues and the associated quantum numbers.

To access longer chains, we use QMC. The present QMC simulations were carried out with the ALPS [49,50] directed loop applications [51,52] in the stochastic series expansion framework [53]. To be precise, these computations were started a while ago. Therefore, we used version 1.3 of the ALPS applications [50] rather than the more recent release 2.0 [54]. The specific heat in a magnetic field can be sensitive to the pseudorandom number generator; so, this needs to be carefully chosen. Here, we used the “Mersenne Twister 19937” pseudorandom number generator [55]. To verify reliability of our results, we performed QMC simulations for (data not shown here) and double-checked them against our ED computations for the same system size.

Figure 4, Figure 5, Figure 6 and Figure 7 show ED () and QMC () results for the specific heat. The QMC simulations become challenging at the lowest temperatures, in particular for finite D and h. This leads to visible statistical error bars at low T, in particular, in Figure 6 and Figure 7; otherwise, statistical errors are negligible. For , finite-size effects are relevant, as demonstrated by visible deviations between the and 100 data in Figure 4 and Figure 5. On the other hand, no further change is visible for larger N, i.e., can be considered as representative of the thermodynamic limit for . Finally, a field of lifts the ground-state degeneracy and opens a sufficiently large gap in the spectrum such that and become indistinguishable (see Figure 6 and Figure 7) and ED suffices to describe the thermodynamic limit.

Figure 4

Specific heat per spin c calculated for , , by exact diagonalization (ED) and Quantum Monte Carlo (QMC). The left panel shows the specific heat itself while the right panel shows the specific heat divided by temperature .

Figure 5

Specific heat per spin c calculated for , , . As in Figure 4, the left panel shows the specific heat itself while the right panel shows the specific heat divided by temperature .

Figure 6

Specific heat per spin divided by temperature calculated for , in a magnetic field .

Figure 7

Specific heat per spin divided by temperature calculated for , in a magnetic field .

For and , the ground state is an multiplet with components. This leads to a difference between the zero-temperature entropies per site for and of Accordingly, the entropy integral , i.e., the corresponding area under the curve of the right panel of Figure 4, is expected to be bigger than that of the corresponding curve by this amount . The QMC data for the specific heat c not only exhibit a maximum at but also a shoulder at (see left panel of Figure 4), corresponding to the two expected features [28,29].

Figure 5 shows the result with the single-ion anisotropy included, still at . The presence of the single-ion anisotropy reduces the ground-state degeneracy to two and opens a gap in the one-magnon spectrum, see Appendix A for details. For , the resulting ground-state entropy is still almost 5% of the total entropy. This leads to a difference between the zero-temperature entropies per site for and of While this is smaller than in the case of , the difference is still visible in the ED data compared with those of QMC, shown in the right panel of Figure 5. From the point of view of physics, the specific heat c in the left panel of Figure 5 may be more instructive. The shoulder-like feature for developed into a sharp peak around for the value while, in turn, the previous global maximum of c became a shoulder around . In any case, these two features can be traced from to finite D.

Finally, we add a magnetic field , corresponding to the experimental case where we actually observed two features in the specific heat (see Figure 3b). Application of a finite field not only lifts the remaining ground-state degeneracy, but opens a sufficiently large gap in the spectrum such that finite-size effects are negligible already for , as mentioned before and shown in Figure 6 and Figure 7. As in the experiment, we observe the emergence of a double-peak structure where both the feature at and in particular the one at shift to higher temperatures with increasing magnetic field (compare Figure 6 and Figure 7).

3.3. Mean-Field Treatment of the Interchain Coupling

In zero external magnetic field, an antiferromagnetic phase transition with a Néel temperature is observed experimentally, as discussed in Section 2. This demonstrates that interchain coupling should be included in a quantitative description, at least for and , even if the numerical results of Section 3.2 already qualitatively reproduce the experiment in a finite magnetic field.

Since the chains are ferromagnetic, we assume that only the total magnetization of one chain acts via an effective field on the neighboring chains. The assumption of only average magnetizations of one chain affecting the neighboring ones is motivated by the exact exchange paths between chains in MnNi(NO)(en) being unknown (compare the crystal structure of Figure 1) and was also made in Ref. [28]. To be precise, one starts from a coupling between chains i and j of the form which one replaces by

We drop the term for the time being, but one should remember to add this term for total energy computations and in particular if one wants to write expectation values as derivatives of the free energy, see also Ref. [56]. This leads to a family of interchain mean-field Hamiltonians where the magnetization of the ith chain should satisfy the self-consistency condition with (recall that we chose units such that ).

We now consider two cases. Firstly, for , we expect antiferromagnetic order that should be described by two types of chains . Furthermore, by symmetry, one expects that . This sign difference can be absorbed by a spin inversion on every other chain, which also flips the sign of the interchain coupling. Therefore, we introduce an effective interchain coupling , where the minus sign will allow us to treat all chains as having the same magnetization . Secondly, for , one stays in a paramagnetic phase where we expect all chain magnetizations to be equal . Now, we straightforwardly set the effective interchain coupling .

Under either of these assumptions, the family of mean-field Hamiltonians (4) reduces to a single interchain mean-field Hamiltonian with

The magnetization should now satisfy the modified self-consistency condition

Recall that in order to cast both the antiferromagnetic case at and the paramagnetic case at in the same single-chain form, it was necessary to introduce different signs for the effective interchain coupling in the two cases. Still, the absolute value of is the same in both cases.

Since the magnetization is easily evaluated even within QMC, it is possible to run a self-consistency loop using a numerical evaluation of the chain magnetization , i.e., one starts with an initial guess such as , recomputes from Equation (8), and iterates until a desired level of accuracy is reached (compare Appendix B for further details). Some ED and QMC results for the self-consistent for are shown in Figure 8. The vertical line in Figure 8 shows an estimate of the Néel temperature that will be discussed in the following Section 3.3.1. One observes in Figure 8 that the estimated Néel temperature varies by almost a factor two as one goes from to 100 spins in a chain. Even with , one still deviates by about 25% from the estimate obtained with . On the other hand, analysis of the data shown in Figure 5 and Figure 9 below indicates that should indeed be sufficient to represent the thermodynamic limit along the chains.

Figure 8

Self-consistent mean-field results for the magnetization with , , , . The vertical line indicates the estimate for the Néel temperature obtained from chains with sites (for details, see text).

Figure 9

Inverse magnetic susceptibility calculated by QMC for a single chain with , , . Furthermore, also shown are the estimated value of the interchain coupling K [28] as well as the experimental Néel temperature K [28,47] divided by K.

3.3.1. Magnetic Susceptibility and Ordering Temperature

The numerical treatment of a single chain yields direct access to where may be included in the self-consistent effective field but is considered to be fixed, i.e., contributions from the self-consistent field are not included in Equation (9).

The magnetic susceptibility should be defined by within the interchain mean-field approximation. Insertion of the definition Equation (8) for the magnetization and some straightforward algebra leads to

The result (11) can be solved for and one finds (As a consequence of the spin inversion that we have applied to half of the chains at , the following is actually not the uniform, but a staggered susceptibility in the case of a vanishing external field.)

This approximation is widely used in the literature (see for example [57,58,59]) and also known under the name “random phase approximation”. Since there are some similarities with the Stoner model of ferromagnetism (see, e.g., chapter 7.4 of [60]), one can also call a “Stoner factor”. Note that the above derivation is essentially the same as the computation on page 66 of [61], but the linearizing assumption has been dropped. Accordingly, we see that Equation (12) also applies for a finite magnetization of a single chain.

A zero of the denominator in Equation (12) signals a second-order phase transition. This yields the standard condition for the Néel temperature

Let us use this condition to take a look at the ordering transition in zero external field where for such that Equation (13) can be evaluated without running a self-consistency cycle. Our QMC results for at are shown in Figure 9 for , 200, and 400. The fact that these three system sizes are essentially indistinguishable on the scale of the figure shows that suffices to represent the thermodynamic limit .

If one assumes the value K (horizontal line in Figure 9) that has been deduced in [28] by fitting the magnetic susceptibility for K, one reads off an ordering temperature K. This deviates by about 12% from the experimental value K, which is remarkably good for a mean-field theory. Conversely, if one insists on the experimental value , one infers an interchain coupling K, which is about 50% larger than the estimate of [28]. In fact, varies quite strongly in this temperature range. Therefore, is not very sensitive to the interchain coupling .

In any case, an interchain coupling of a few percent suffices to yield an antiferromagnetic ordering temperature at that is of the same order as the coupling in an individual chain, reflecting strong ferromagnetic ordering tendencies of the decoupled chains.

3.3.2. Specific Heat

Let us now take a closer look at the specific heat in interchain mean-field theory. As in the case of the magnetic susceptibility, the numerical treatment of the individual chains provides convenient access to where may again be included in the self-consistent effective field but is considered to be fixed.

The self-consistent magnetization is also temperature-dependent such that the specific heat should be written as a first derivative of the internal energy

The temperature derivative can in principle be calculated numerically. For reasons of numerical stability, particularly in a Monte-Carlo setting, it is nevertheless preferable to carry the derivatives out analytically. Since we are not aware of such an analysis having been presented before, we present it here in some detail. With the help of , we find from Equation (15) that

This expression contains another derivative for which we can find an expression that is very similar to Equation (12) (including a “Stoner factor” ):

Noting the relation the combination of Equations (16) and (17) can also be written in the following form:

In this form, the sign of the second term is evident. This form is also useful for the purpose of evaluation since Equation (19) contains only quantities that can be related to static expectation values for a single chain with a fixed value of via Equations (9), (14), and (18). The only object that is non-standard is the crosscorrelator in Equation (18); however, it represents exactly the same observable as was used in Ref. [62] to compute the adiabatic cooling rate by QMC.

3.3.3. Comparison with Experimental Specific Heat

We are now in a position to perform a comparison with the experimental results for the specific heat of Figure 3. Figure 10, Figure 11 and Figure 12 show the results for , , and (corresponding to the experimental magnetic fields , T, and T, respectively). For , we used QMC with , while for and , we used ED with . These systems sizes should be sufficiently large to render finite-size effects negligible according to the discussions in Section 3.2. From a technical point of view, we note that at and in the paramagnetic phase, such that and the correction term in Equation (19) vanishes, i.e., , and the blue circles are identical to the green diamonds in Figure 10 for .

Figure 10

Magnetic specific heat per spin c for K, K in zero magnetic field in comparison with the experimental results for . Theoretical results are obtained by QMC with both for decoupled chains () and with a self-consistent mean-field treatment for K.

Figure 11

Magnetic specific heat per spin divided by temperature for K, K in a magnetic field in comparison with the experimental results for T. Theoretical results are obtained by ED with both for decoupled chains () and with a self-consistent mean-field treatment for K.

Figure 12

Magnetic specific heat per spin divided by temperature for K, K in a magnetic field in comparison with the experimental results for T. As in Figure 11, theoretical results are obtained by ED with both for decoupled chains () and with a self-consistent mean-field treatment for K.

Figure 11 and Figure 12 show that the interchain coupling leads only to small corrections for a magnetic field ; the trend is towards the experimental data, but the shift by interchain coupling does not change the situation significantly. Nevertheless, the two theory curves and the experimental one in Figure 11 and Figure 12 exhibit double-peak structures where the two peaks are located at very similar temperatures between theory and experiment.

Figure 10 demonstrates that in zero field (), interchain coupling is not only essential for reproducing the ordering transition to good accuracy, as we have seen before, but that thanks to the “Stoner factor”, the correction term in (19) dominates the specific heat just below the ordering transition and thus gives rise to the characteristic ordering peak. We note, however, that the singularity in the denominator of Equation (19) is deceptive since the numerator (18) also vanishes such that c has a finite limit for . Consequently, our interchain mean-field theory remains in the universality class of Landau theory [63] with a specific heat exponent .

4. Magnetocaloric Properties

The strong dependence of the specific heat of MnNi(NO)(en) on an applied magnetic field promises a strong magnetocaloric effect and potential relevance to low-temperature magnetic refrigeration by adiabatic demagnetization, e.g., see Refs. [64,65]. Therefore, let us have a closer look at its magnetocaloric properties.

Figure 13 shows the experimental magnetic entropy that is obtained by integrating the experimental results for the specific heat of Figure 3b with respect to temperature T. The curve corresponds to the one shown already in Figure 2b. Figure 13 shows that the magnetic entropy is significantly reduced by applying a magnetic field of T, or even T, corresponding to polarization of the spin system by the applied magnetic field. Consequently, we expect cooling of the spin system during adiabatic demagnetization. Let us consider, for example, an ideal adiabatic process that starts with K for T. We read off from Figure 13 that the same entropy is found at for K, i.e., adiabatic demagnetization from T to would cool from an initial temperature K to a final temperature of K. Likewise, an ideal adiabatic process starting with K at T would cool to K during a single ideal adiabatic demagnetization process. These are relatively large effects in the liquid Helium range, which is also remarkable since one is cooling through a phase transition into a magnetically ordered state. The main caveat is that the processes of the two examples exploit only 8% or 19% of the total magnetic entropy in the first and second case, respectively.

Figure 13

Magnetic entropy per mole of MnNi(NO)(en) in magnetic fields of , , and T, respectively.

Next, let us comment on a numerical description. The entropy is not directly accessible in QMC simulations such that we resort to ED even if this leads to stronger finite-size effects. Furthermore, for the full h and T dependence of the magnetic entropy S, we would have to model the ordered state in an external magnetic field (gray shaded region in the inset of Figure 3b). However, this is expected to correspond to a canted spin configuration and is thus beyond the present investigation. Therefore, we also neglect interchain coupling, i.e., we focus on a situation corresponding to the one discussed in Section 3.2 (see, however, Appendix C for a discussion of simple single-site mean-field theory). Figure 14 shows the corresponding result for the entropy (now normalized per spin) of an chain. This density plot of permits to immediately read off the magnetocaloric effect. In particular, the isentropes, corresponding to the white lines in Figure 14, directly show the behavior under an adiabatic process. Finite-size effects are expected to be small for T (corresponding to , compare Figure 6 and Figure 7), but they are known to be relevant throughout the temperature range of Figure 14 for (, compare Figure 5).

Figure 14

Entropy S per spin calculated by ED for an isolated chain with spins and K, K.

We also read off cooling by adiabatic magnetization from Figure 14 with a size of the effect corresponding to the experimental data of Figure 13. Note that the entropy is normalized to mole in Figure 13 and per spin in Figure 14, amounting to a factor two difference in addition to the factor R. Since we have ferromagnetic chains, the strongest cooling occurs in Figure 14 upon approaching a zero external field. If one adds antiferromagnetic interchain coupling, we expect to recover the magnetically ordered phase that is observed experimentally for T (compare inset of Figure 3b) and then cooling might actually occur when entering this phase (compare Figure A4 and the related discussion in Appendix C for the behavior in single-site mean-field theory). Indeed, Figure 10 demonstrates that interchain coupling reshuffles entropy from low temperatures to the ordering transition such that the most significant cooling probably occurs around it.

Figure A4

Single-site mean-field results for the entropy per spin at , , and (a) and at (b). Red lines are the phase boundaries from Figure A3.

5. Conclusions and Perspectives

We carried out specific heat measurement in zero and applied fields on the bimetallic chain compound MnNi(NO)(en). By determining the lattice contribution of the specific heat, we have extracted the magnetic specific heat . For the first time, in its temperature dependence, we verify a long-predicted double-peak-like structure. Comparison with numerical calculations for the bimetallic , ferromagnetic spin chain yields a very close resemblance on a semi-quantitative level.

Alternating spins are not the only mechanism that may give rise to a double-peak structure in the specific heat. For example, also a ferromagnetic chain alone can give rise to such structures when subjected to a strong single-ion anisotropy D [66]. However, the numerical data of Section 3.2 (and further results that we do not show here) demonstrate that these two features are already present at in the , chain and can be traced to finite D even if the presence of a single-ion anisotropy does affect the behavior of the specific heat at a quantitative level. Hence, we conclude that our experimental observation of a double-peak-like structure in the specific heat directly reflects the alternating spins and along the chains. The application of an external magnetic field to MnNi(NO)(en) is essential to suppress magnetic order and, thus, reveal this double-peak feature experimentally.

The ordered phase that is observed in MnNi(NO)(en) for low temperatures and small applied magnetic fields is due to an antiferromagnetic interchain coupling. Although its absolute value is much smaller than the ferromagnetic coupling along the chains, it has a strong effect at low temperatures and in the absence of a magnetic field. In order to describe this ordered phase, we developed a mean-field treatment of interchain coupling. The combination of QMC simulations for isolated chains and such an interchain mean-field theory not only yields a remarkably accurate value for the ordering transition temperature using previously determined parameters [28] but also yields excellent agreement for the full temperature dependence of the magnetic specific heat. For fields , the mean-field corrections are small, reflecting the smallness of the interchain coupling constant .

Beyond the very close resemblance on a qualitative level, there are some quantitative differences between experiment and theory. For instance, while in the calculations the maximum of is found close to , in the experiments it is observed at . These small differences may be due to the single-ion anisotropy being located on both the Ni and Mn sites, and not just the Mn ones, or effects of interchain coupling beyond mean-field theory. However, a further refinement of the model would require additional information about the excitation spectrum such as inelastic neutron scattering.

Another theoretical challenge concerns the theoretical description of the ordered state in a magnetic field. For and strong fields along the anisotropy axis (), one may restrict the discussion to magnetization along the z-axis only. However, for a magnetic field applied at an angle to the anisotropy axis, and also for ordered phases where the ordered moment cants away from the field/anisotropy axis, it will in general be necessary to replace the last term in (4) by vectors, i.e., by . This generalization can be implemented in single-site mean-field theory; however, such a strong approximation fails to be quantitatively accurate for the present situation (compare Appendix C). By contrast, such generic field directions in the interchain mean-field theory of Section 3.3 break conservation of total , render the computations even more challenging, and thus go beyond the present investigation.

Finally, we have shown that the strong sensitivity of MnNi(NO)(en) to even small applied magnetic fields gives rise to a strong magnetocaloric effect, i.e., large cooling by adiabatic demagnetization from initial fields B on the order of 1 T. Even if the magnetic entropy of MnNi(NO)(en) may be a bit small for practical applications in the temperature range of interest, this observation suggests materials with competing strong ferromagnetic and weaker antiferromagnetic interactions as promising candidates for efficient low-temperature refrigeration.