Spin Seebeck effect of correlated magnetic molecules

In this paper we investigate the spin-resolved thermoelectric properties of strongly correlated molecular junctions in the linear response regime. The magnetic molecule is modeled by a single orbital level to which the molecular core spin is attached by an exchange interaction. Using the numerical renormalization group method we analyze the behavior of the (spin) Seebeck effect, heat conductance and figure of merit for different model parameters of the molecule. We show that the thermopower strongly depends on the strength and type of the exchange interaction as well as the molecule’s magnetic anisotropy. When the molecule is coupled to ferromagnetic leads, the thermoelectric properties reveal an interplay between the spin-resolved tunneling processes and intrinsic magnetic properties of the molecule. Moreover, in the case of finite spin accumulation in the leads, the system exhibits the spin Seebeck effect. We demonstrate that a considerable spin Seebeck effect can develop when the molecule exhibits an easy-plane magnetic anisotropy, while the sign of the spin thermopower depends on the type and magnitude of the molecule’s exchange interaction.


Results
The schematic of the studied system is presented in Fig. 1. A high-spin magnetic molecule, of spin S = S c + s , is coupled to external magnetic leads, whose magnetizations point in the same direction. The molecule's easy axis is assumed to coincide with the direction of leads' magnetizations. It is further assumed that a single molecular energy level is active in transport and this orbital level is coupled through an exchange interaction J to the core spin S c of the molecule. Thus, the Hamiltonian of the molecule reads as 52-54 , where ε and U denote the energy of the molecule's orbital level and Coulomb correlation energy between two electrons of opposite spin occupying that level. n σ ≡d † σd σ is the occupation number operator for an electron of spin σ and d † σ ( d σ ) is the corresponding creation (annihilation) operator. The spin operator for an electron occupying the orbital level is denoted by ŝ ≡ (1/2) σ σ ′d † σ σ σ σ ′d σ ′ , where σ ≡ (σ x , σ y , σ z ) denotes the vector of the Pauli matrices, and Ŝ c is the operator for the core spin of the molecule. The two spins are coupled by the exchange interaction J, which can be either of ferromagnetic ( J > 0 ) or antiferromagnetic ( J < 0 ) type, depending on the sign of J. The molecule can be subject to magnetic anisotropy denoted by D and Ŝ z is the zth component of the molecule spin operator Ŝ =Ŝ c +ŝ . In calculations, we assume S c = 1 , while the spin of electrons on the orbital level is given by s = 1/2 or s = 0 , depending on its occupancy. Consequently, the total molecule's spin is S = 3/2 for singly occupied orbital level or S = 1 in the case when the occupation is even.
The tunneling processes between the molecule and the leads are described by the following Hamiltonian where q = L for the left and q = R for the right electrode, the operator ĉ † qkσ ( ĉ qkσ ) creates (annihilates) an electron with energy ε qkσ , momentum k and spin σ in the q-th lead, and v qkσ denotes the corresponding tunnel matrix elements. The leads are described within the non-interacting quasi-particle approximation by Having defined the three parts of the Hamiltonian, the total Hamiltonian is given by, Ĥ =Ĥ mol +Ĥ tun +Ĥ leads .
The tunnel coupling between the molecule and the leads gives rise to the broadening of the orbital level, which can be described by, Ŵ qσ = πρ qσ v 2 qσ , where ρ qσ is the spin-dependent density of states at the Fermi level in the lead q and we assumed momentum-independent tunnel matrix elements v qkσ ≡ v qσ . We then define the full broadening function, which for spin σ can be written as, Ŵ σ = (1 + ηp)Ŵ , where Ŵ = Ŵ L + Ŵ R   Schematic of the considered molecular junction. It consists of a magnetic molecule of spin S = S c + s tunnel-coupled to external leads, with spin-dependent coupling strengths, Ŵ Lσ and Ŵ Rσ , for the left and right lead. S c is the spin of the molecule's magnetic core, while s denotes the spin of electrons occupying the orbital level. The molecule is assumed to effectively possess one orbital level, through which transport takes place, that is exchange-coupled (with coupling strength J) to the molecule's core spin S c . The orbital level is characterized by on-site energy ε and the Coulomb correlations U. There is a voltage ( �µ σ ) and temperature ( T ) gradient applied to the system. In the case of magnetic contacts, the voltage gradient may be spin-dependent. In considerations we assume S c = 1 , while s = 1/2 ( s = 0 ) when the orbital level occupancy is odd (even). www.nature.com/scientificreports/ and Ŵ q = Ŵ q↑ + Ŵ q↓ , whereas p is the effective spin polarization of the left and right ferromagnetic lead, p = (p L + p R )/2 , and η = 1 ( η = −1 ) for spin-up (spin-down) electrons.
Thermopower in the case of nonmagnetic leads. We focus on the thermoelectric transport properties of the considered molecular junction in the linear response regime. The interesting thermoelectric phenomena happening across the magnetic molecule can be quantified using the transport coefficients, such as the electrical conductance G, the thermopower (Seebeck coefficient) S, as well as the thermal conductance κ and the thermoelectric figure of merit ZT. In the linear response regime, these quantities can be expressed in terms of Onsager integrals L nσ 55 where T σ (ω) is the energy-dependent transmission coefficient for the spin channel σ , f is the Fermi-Dirac distribution function and µ denotes the electrochemical potential. The electrical conductance G and the electronic contribution to the thermal conductance κ can be then found from where e is the electron charge, T denotes the temperature and L n = σ L nσ . On the other hand, the thermopower S is defined as S = −[�V /�T] J=0 , on the condition of vanishing of the charge current J. Hence, S can be expressed in the form Having defined G, S and κ , one can obtain the thermoelectric figure of merit ZT ≡ GS 2 T/κ. Let us first consider the case of nonmagnetic molecular junction. The linear conductance and the Seebeck coefficient as a function of the molecule's orbital level energy ε and temperature T are shown in Fig. 2. The first column shows the results for J = 0 , while the second (third) column corresponds to the case of the ferromagnetic (antiferromagnetic) exchange interaction J. The case of J = 0 is shown just for reference and allows us to clearly reveal the effects stemming from the presence of large-spin molecule. To begin with, we consider the behavior of the linear conductance. The largest changes with lowering the temperature are visible when the orbital level is singly occupied, i.e. for −U ε 0 . In this case, the Kondo effect can develop at sufficiently low temperatures,  www.nature.com/scientificreports/ such that T T K , where T K is the Kondo temperature 17 . Once T ≪ T K , the conductance reveals a plateau as a function of ε of height G = 2e 2 /h in the case of J = 0 18,20 . However, for magnetic molecules described by the Hamiltonian (1), the low-temperature behavior strongly depends on the type of exchange interaction J 56 . For ferromagnetic exchange, the Kondo effect always develops with lowering T, however, T K becomes reduced compared to the case of J = 0 . On the other hand, in the antiferromagnetic-J case, once |J| T K , the spin on the orbital level strongly binds with the magnetic core spin, which results in the suppression of the conductance through the system. In this case, in the singly occupied orbital regime, one only observes a small enhancement (a local maximum) of G with decreasing T followed by its strong suppression 56 . The different scenarios discussed above give rise to a unique behavior of the Seebeck coefficient, which is presented in the bottom row of Fig. 2. The first observation is that, as expected 9,12 , the thermopower changes sign with respect to the particle-hole symmetry point, ε = −U/2 . Moreover, one can see that the behavior of S for T Ŵ is hardly affected by the type of exchange interaction J. This is because in our considerations |J| < Ŵ and the effects of finite J can be visible only when T |J| . Let us anyway summarize the main features of the high-temperature behavior. One observes two pronounced maxima for T ≈ Ŵ in the case when the orbital level is either empty or doubly occupied. On the other hand, when moving to the Coulomb blockade regime where the orbital level is singly occupied, the thermopower changes sign and two local extrema develop antisymmetrically around ε = −U/2 for T ≈ Ŵ 21,36 . With lowering the temperature, distinct features appear, resulting from the interplay of the correlations driving the Kondo effect and the molecule's exchange interaction. First of all, one can note that in the case of J = 0 , the regions of large |S| extend from the empty and doubly occupied regimes for T ≈ Ŵ downwards to low temperatures in the single occupancy regime. For singly occupied orbital level, the thermopower exhibits then a sign change as a function of T, see Fig. 2d, which signals the relevance of the Kondo correlations 21 . A qualitatively similar behavior can be observed in the case of ferromagnetic exchange interaction J, see Fig. 2e, with the main difference associated with smaller temperatures at which the corresponding sign change occurs. This is associated with the fact that the sign change occurs at the onset of the Kondo correlations and, because the Kondo effect develops at much lower temperatures in the case of ferromagnetic J compared to the case of J = 0 (see Fig. 2a,b), one observes that the sign change in S is also shifted to lower temperatures. However, this shift is not proportional to the corresponding shift visible in the behavior of G. It is because while much smaller temperatures are needed for the full development of the Kondo effect in the case of ferromagnetic J, the temperature associated with the onset of the Kondo correlations only weakly decreases with J. This is why the crossover for J > 0 is only slightly shifted to lower temperatures compared to the case of J = 0 , cf. Fig. 2d,e. Interestingly, qualitatively new features compared to the case of J ≥ 0 can be observed in the case of antiferromagnetic exchange interaction, where an additional sign change at low temperatures is present, see Fig. 2f.
Further insight into the behavior of the thermoelectric properties can be obtained from the inspection of Fig. 3, which presents the temperature dependence of G, S, κ and ZT for different values of the exchange interaction J, as indicated. This figure is generated for the case when the orbital level is detuned from the particlehole symmetry point, such that a considerable thermopower can be observed. In addition, the orbital level is assumed to be singly occupied ( ε = −U/3 ), such that the system is in the local moment regime and the Kondo correlations are relevant at sufficiently low temperatures. Let us start with the analysis of the linear conductance. In the case of ferromagnetic J (see the solid lines in the figure), the enhancement of J results in a decrease of the Kondo temperature. In this case the ground state is always two-fold degenerate with total spin given by S = S c + 1/2 and the Kondo effect develops irrespective of J, though T K becomes very cryogenic with increasing J 56 . Consequently, we observe mainly quantitative changes in Fig. 3a, while qualitative behavior is the same. On the other hand, the situation is completely different in the case of antiferromagnetic exchange interaction J < 0 , see the dashed lines in Fig. 3. Now, with increasing |J|, the ground state becomes S = S c − 1/2 , since the spin on the orbital level binds anti-ferromagnetically with the molecule's core spin. Because of that, the Kondo effect is quenched once the temperature becomes lower than the energy scale responsible for this antiferromagnetic state. As a result, the temperature dependence of conductance exhibits a nonmonotonic behavior 56 , see Fig. 3a. We note that such behavior is similar to the two-stage Kondo effect observed in side-attached double quantum dots where the hopping induces an antiferromagnetic interaction between the dots 57 .
The behavior of the conductance has a strong influence on the Seebeck coefficient. This is because the temperature dependence of conductance reflects the energy dependence of the transmission coefficient T (ω) and, from the Sommerfeld expansion, the thermopower at low temperatures can be estimated from (note that temperature T is in units of energy) where T (ω) = σ T σ (ω) . One can thus see that the thermopower is related to the monotonicity of the variation of the spectral function with energy. This is why in the case of ferromagnetic exchange interaction we only observe qualitative changes in S, see the solid lines in Fig. 3b. The enhancement of J suppresses then the Kondo temperature, which is seen in the behavior of S as a shift of the local minimum towards smaller temperatures. However, a completely different scenario develops for antiferromagnetic J, where an additional sign change occurs and S exhibits an extra maximum for temperatures corresponding to the energy scale at which G(T) starts decreasing with lowering T, see the dashed curves in Fig. 3b.
Despite a spectacular impact of exchange interaction on the thermopower, its effect on the thermal conductance is less pronounced, see Fig. 3c. This is because κ is considerable only at energy scales corresponding to the coupling strength and Coulomb correlations and thus, as long as |J| ≪ Ŵ , κ hardly depends on the type and magnitude of exchange interaction. On the other hand, the figure of merit ZT displays new peaks in the case of antiferromagnetic J, see Fig. 3d, which are associated with the above-discussed maxima emerging in the www.nature.com/scientificreports/ temperature dependence of S. We also note that the influence of J on κ is more visible when one plots κ/T , see the inset in Fig. 3c. It is nicely visible that the qualitative behavior of κ/T resembles that of the linear conductance, which is a direct consequence of the Wiedemann-Franz law. A similar behavior has been recently observed for T-shaped double quantum dots 23 .
Effect of magnetic anisotropy. We now focus on elucidating the role of magnetic anisotropy on the thermoelectric properties of the considered molecular junction. We consider both an easy-axis ( D > 0 ) and easy-plane ( D < 0 ) types of magnetic anisotropy. First we note that in the case of ferromagnetic exchange interaction between the orbital level and molecule's core spin the magnetic anisotropy has a very moderate influence on the thermoelectric properties. It does not lead to new qualitative behavior as long as |D| is smaller than the corresponding Kondo temperature, therefore in the following we just analyze the case of antiferromagnetic exchange coupling J. The temperature dependence of the conductance and the Seebeck coefficient for this situation is shown in Fig. 4, where the left (right) column corresponds to the easy-plane (easy-axis) magnetic anisotropy case. This figure was generated for the same orbital level position as in Fig. 3, such that the orbital level is detuned from the particle-hole symmetry point, while the system stays in the local moment regime. In the absence of anisotropy the conductance displays a nonmonotonic dependence, characteristic of the antiferromagnetic exchange coupling. When an easy plane anisotropy arises in the system and the molecule possesses a half-integer spin, it results in a two-fold degenerate ground state of S = 1/2 , such that the Kondo effect can be restored. This is clearly seen in Fig. 4a, where one observes an upturn of the conductance with lowering T. This in turn has a considerable impact on the Seebeck coefficient, which exhibits an additional sign change, see Fig. 4b. On the other hand, in the case of easy plane anisotropy, the low-temperature conductance again becomes increased with D, as can be seen in Fig. 4c. This is however just associated with a decreased exchange interaction between the molecule's core spin and the spin of the orbital level, and not with the reinstatement of the Kondo effect. Consequently, while the Seebeck coefficient strongly depends on D, no additional sign changes are present. In fact, the maximum in S for D = 0 becomes suppressed with increasing D and smears out completely once D ≈ T K , see Fig. 4d. www.nature.com/scientificreports/ Thermopower in the case of ferromagnetic leads. We now turn to the discussion of thermoelectric properties in the case of ferromagnetic electrodes. For the ferromagnetic contacts, we assume a moderate spin polarization 58-60 , p = 20% . In this section we assume that the spin relaxation in the contacts is relatively fast, such that no spin accumulation develops and the induced potential gradient does not depend on spin, �µ ↑ = �µ ↓ . In this regime, although the spin Seebeck effect does not develop, the spin-dependence of tunneling processes greatly modifies the thermoelectric transport properties of the system as compared to the nonmagnetic case. We also note that in the absence of spin accumulation the formulas for spin-dependent thermoelectric coefficients are the same as in case of nonmagnetic leads 35,36 .
The linear conductance and the Seebeck coefficient as a function of ε and T calculated for different values of J are displayed in Fig. 5. In the behavior of the conductance one can clearly observe the signatures of an effective exchange field that develops in the molecule coupled to ferromagnetic electrodes 54,61 . Such an exchange field, which within the perturbation theory at zero temperature and for J = 0 can be described as 61 , results in a spin-splitting of the molecule's orbital level, when it is detuned from the particle-hole symmetry point of the model ε = −U/2 . In the case of considered molecule, this field depends in a nontrivial way on the properties of the molecule, such as J, D and S c , however, it still vanishes whenever ε = −U/2 54,62 . If the exchange field splitting becomes larger than the Kondo temperature, it suppresses the Kondo resonance. As a consequence, the low-temperature conductance in the Coulomb blockade regime is generally decreased except for ε = −U/2 , where a local maximum as a function of ε is present, see Fig. 5a. A similar conductance suppression can be seen in the case of ferromagnetic exchange interaction J presented in Fig. 5b. However, now the suppression of G is larger as compared to the case of J = 0 , which is associated with a smaller Kondo temperature when J > 0 . On the other hand, when J is antiferromagnetic, mainly quantitative changes can be observed in the conductance behavior, cf. Figs. 2c and 5c. Namely, the region of suppression of low-temperature conductance is now smaller, which is an www.nature.com/scientificreports/ indication that the exchange field hinders the formation of an antiparallel spin state between the orbital level and the core spin. Because of that, the decrease of conductance due to that formation is correspondingly weakened. The signatures of the interplay of exchange field, Kondo correlations, and the molecule's exchange interaction are also visible in the behavior of the thermopower, which is presented in the bottom row of Fig. 5. One can see that the main changes are visible in the low-temperature behavior, which reveals the energy scale associated with the exchange field. More specifically, in the case of ferromagnetic exchange interaction, an additional sign change in S occurs as compared to the case of nonmagnetic leads, whereas for antiferromagnetic J the Seebeck coefficient becomes generally reduced, cf. Figs. 2 and 5.
To shed more light onto the thermoelectric behavior of the considered magnetic molecular junction, in Fig. 6 we present the temperature dependence of the conductance, Seebeck coefficient, heat conductance and figure of merit calculated for different values of exchange interaction, as indicated. This figure was determined for a relatively low value of detuning from the particle-hole symmetry point, ε = −0.48U , such that the exchange field effects rather compete than surpass other energy scales, and the system's behavior is most interesting. First of all, one can note that below a certain temperature the conductance stops changing any more and retains its lowtemperature value. This characteristic energy scale is set by the exchange field �ε exch -when this field is larger than thermal energy, it determines the transport behavior and no further dependence on lowering T is visible, see Fig. 6. The Seebeck effect becomes then suppressed and so does the figure of merit. This can be understood by referring to the Sommerfeld expansion: once the exchange field is the dominant energy scale the conductance becomes constant and so does the low-energy transmission coefficient. Consequently, [∂T (ω)/∂ω] ω=0 ≈ 0 and thus S ≈ 0 . This explains why S ≈ 0 and ZT ≈ 0 for T/ Ŵ �ε exch / Ŵ ≈ 10 −4 , see Fig. 6b,d.
Let us now focus on the most interesting behavior, which is present for temperatures of the order and larger than the exchange field, and let us start with the case of ferromagnetic exchange interaction J > 0 . One can see that increasing J results in suppression of the low-temperature conductance, see Fig. 6a. This results from the fact that increasing J leads to lowering of T K , and once �ε exch T K , the Kondo peak becomes suppressed. This is also visible in the thermopower and the figure of merit where a new maximum emerges at energy scale corresponding to the exchange field T ≈ �ε exch , see Fig. 6. While in the case of ferromagnetic exchange interaction mainly qualitative effects are visible, the case of antiferromagnetic exchange is completely different, see the dashed lines in Fig. 6. In this case, increasing |J| gives rise to the suppression of the conductance due to the formation of an antiferromagnetic spin state between the molecule's orbital level and its magnetic core. When |J| is relatively low, the exchange field wins over the antiferromagnetic interaction and only a small suppression of G as a function of T is present, see e.g. the curve for J = −0.3T K in Fig. 6a. However, further enhancement of |J|, stabilizes the antiferromagnetic state of the molecule and a full suppression of G is obtained once J −2T K . This behavior gives rise to a new maximum visible both in S and ZT. Moreover, while the maximum associated with the conductance drop as T is lowered moves to higher energies with decreasing J ( J < 0 ), there is an extra maximum www.nature.com/scientificreports/ visible just at the energy scale corresponding to T ≈ �ε exch . As a consequence, the temperature dependence of the Seebeck coefficient displays an interesting triple-peak structure, see Fig. 6b.
To complete the picture, in Fig. 7 we present the temperature dependence of G and S calculated for different values of magnetic anisotropy. Similarly to the case of nonmagnetic leads, we display the data for antiferromagnetic exchange interaction, which shows the most interesting behavior. Let us first analyze the case of uniaxial anisotropy, which is presented in the right column of Fig. 7. One can see that finite anisotropy gives rise to an enhancement of the low-temperature conductance. This is a consequence of the fact that anisotropy breaks the symmetry of the antiferromagnetic state of the molecule responsible for the conductance suppression. This effect gives rise to a local maximum in the Seebeck coefficient that develops at the energy scale of the order of magnetic anisotropy, see Fig. 7d. The case when the molecule exhibits easy-plane type of anisotropy is shown in the left column of Fig. 7. Now, one can observe a very strong dependence of both G and S on the magnitude of magnetic anisotropy. First of all, the low-T conductance exhibits a nonmonotonic dependence on D < 0 . Once the easy-plane anisotropy is present in the system and |D| 0.1T K , G becomes suppressed. However, this tendency becomes reversed when D −0.1T K , such that one observes an enhancement of G at low temperatures, see Fig. 7a. This is associated with the formation of a doublet ground state in the molecule in the case of considerable easy-plane anisotropy. Now, however, one witnesses a subtle interplay between the antiferromagnetic exchange interaction, the easy-plane magnetic anisotropy, the exchange field that splits the doublet state of the molecule and the Kondo correlations. The antiferromagnetic J gives rise to the suppression of G at low temperatures, which is however slightly hindered by the exchange field. On the other hand, turning on D ( D < 0 ), results initially in a larger suppression of the conductance, nevertheless, increased values of |D| eventually make the doublet state the ground state of the molecule, enhancing thus G due to the Kondo effect. Consequently, for sufficiently large |D|, the conductance shows a pronounced Kondo resonance, see Fig. 7a. The behavior of the conductance is clearly revealed in the temperature dependence of the Seebeck coefficient, which is shown in Fig. 7b. One can see that for values of D such that the conductance starts increasing, the thermopower exhibits a considerable maximum, which actually develops for T ≈ �ε exch . With further increase of |D|, this maximum becomes however decreased and its position moves towards higher temperatures of the order of the Kondo temperature, see Fig. 7b. This large enhancement of the Seebeck coefficient is a result of interplay between the intrinsic properties of the molecule, such as its magnetic anisotropy, and the ferromagnetism of the leads. Unfortunately, this effect On the other hand, the Seebeck coefficient in the case of finite spin bias is given by 35,36 , S = − 1 2|e|T L 1↑ /L 0↑ + L 1↓ /L 0↓ , whereas the heat conductance can be expressed as, κ = 1 T σ L 2σ − L 2 1σ /L 0σ . We would like to emphasize that in the case of ferromagnetic leads in the absence of spin accumulation, the Seebeck effect is due to a generated voltage difference that is not spin dependent. However, the tunneling processes themselves do depend on spin, therefore one then observes the spin-dependent Seebeck effect. On the other hand, when the spin accumulation is relevant in the leads, a spin bias becomes generated, which gives rise to the spin Seebeck effect. Such spin thermopower results from the spin splitting of the chemical potentials in the contacts and is associated with the corresponding spin current. Moreover, we would like to note that spin Seebeck effect due to electron transport is also referred to as the spin-dependent Seebeck effect, whereas the typical Seebeck effect in the case of magnetic contacts is referred to as spin-resolved Seebeck effect 29,32 . However, in our paper we adopt the notion introduced by Świrkowicz et al. 35 .  Fig. 8. First of all, we note that the Seebeck effect at higher temperatures T ≈ Ŵ behaves generally very similarly as in the case of no spin accumulation in the leads, however, its low-temperature behavior is changed. While in the absence of spin accumulation, S exhibits an additional sign change for −U ε 0 and T 0.01Ŵ , in the case of long spin relaxation these features are not seen anymore. Instead, the Seebeck coefficient changes sign only once around T ≈ Ŵ and extends to low temperatures till it is quenched by the exchange field, see the first row of Fig. 8. On the other hand, the spin Seebeck effect displays completely different behavior. As can be seen, the only sign change occurs across the particle-hole symmetry point when tuning the position of the molecule's orbital level. Moreover, the overall behavior is reversed as compared to the Seebeck coefficient. In the regions where S is generally negative S S is positive and vice versa. This is directly associated with the definition of S and S S -while S captures the spin-resolved contributions due to hole and electron processes, S S presents mainly the difference between the spin-dependent components. As can be seen in the bottom row of Fig. 8, the spin-up contribution dominates the thermopower for ε > −U/2 , while for ε < −U/2 the spin-down thermopower is dominant. Such behavior is visible in all considered cases, namely, for J = 0 as well as for ferromagnetic and antiferromagnetic exchange interaction. However, there are some subtle differences. First of all, the spin Seebeck effect is finite for lower temperatures in the case of finite J, as compared to the case with J = 0 . Moreover, while for ferromagnetic J, S S exhibits a local maximum for T ≈ 0.01Ŵ , for antiferromagnetic J, there is a small local minimum in the temperature dependence of spin Seebeck effect, see Fig. 8f.
The above-mentioned features are better visible in the temperature dependence of the thermopower and spin thermopower presented in Fig. 9. Moreover, it turns out that the above discussion is not fully complete, since one can now clearly observe that a sign change of both the Seebeck as well as spin Seebeck effect can develop, provided that the exchange interaction is sufficiently large. In the case of ferromagnetic exchange, a large minimum in S develops with increasing J at energy scale corresponding to the exchange field, see Fig. 9a. On the other hand, for the spin Seebeck effect a maximum forms at the same temperature for which S exhibits a minimum. An interesting behavior occurs for antiferromagnetic exchange interaction, when with increasing |J|, both the thermopower and spin thermopower change sign once J −T K . More specifically, S exhibits then a pronounced maximum, whereas a considerable minimum develops in S S for T ≈ �ε exch , see Fig. 9. The sign change of the thermopower is associated with the fact that the exchange field effects, which determine the sign of S and S S for low values of J, become overwhelmed by antiferromagnetic J, once the exchange interaction becomes sufficiently large, i.e. |J| �ε exch .
To further understand the spin thermopower of magnetic molecules, we present the impact of magnetic anisotropy on the Seebeck and spin Seebeck effects for the case of antiferromagnetic exchange interaction in Fig. 10. First of all, for the case of uniaxial type of magnetic anisotropy, as shown in the right column of Fig. 10, there is a suppression in the minima (maxima) present in the Seebeck (spin Seebeck) effect with increasing anisotropy. This behavior is in fact similar to that observed in the case of no spin accumulation discussed in previous section. However, the case of an easy plane type of anisotropy presented in the left column of Fig. 10 shows far more interesting behavior. The thermopower and spin thermopower exhibit a change in sign around T ≈ �ε exch ≈ 10 −4 for values of D −0.1T K . For temperatures T �ε exch , a considerable peak, positive (negative) for the spin   www.nature.com/scientificreports/ develops for S ( S S ), see the left column of Fig. 10. Further increase in the magnitude of anisotropy D suppresses the local extrema in the thermopower, as seen in the case of no spin accumulation. Consequently, one observes a strong nonmonotonic dependence of the (spin) Seebeck effect on the value of the easy-plane anisotropy. Large values of S and S S develop for such intrinsic parameters of the molecule that the temperature dependence of G is most spectacular, i.e. when the Kondo effect starts being restored with increasing |D|. However, once a plateau in G is formed at low temperatures due to the Kondo effect, which happens for considerable values of D < 0 , the thermopower becomes suppressed, see Fig. 10.

Discussion
We have determined the thermoelectric properties of large-spin magnetic molecules attached to both nonmagnetic and ferromagnetic electrodes. Our analysis was focused on the strong correlation regime where the Kondo effect can emerge. To accurately address the system's behavior in this nonperturbative regime, we have employed the numerical renormalization group method, which allowed us to study the electrical and heat conductances, the Seebeck effect and the figure of merit in the full parameter space of considered molecular junction. In particular, we have considered the cases of both ferromagnetic and antiferromagnetic exchange interaction J between the orbital level of the molecule and its magnetic core. Moreover, we have also analyzed the effect of finite magnetic anisotropy of the molecule.
In the case of nonmagnetic contacts, we have shown that the behavior of the Seebeck effect strongly depends on the type of exchange interaction J. In the case of ferromagnetic exchange, we have shown that the Seebeck coefficient displays features qualitatively similar to the case of J = 0 (quantum dot case). However, due to a reduced Kondo temperature, the sign change in S occurs at a lower temperature and the thermopower was found to retain finite values extending to much lower T compared to the quantum dot case. Interestingly, in the case of antiferromagnetic exchange interaction, we have found a new sign change of the thermopower, at energy scale corresponding to the exchange interaction between the orbital level and molecule's core spin. Moreover, we have also determined the influence of finite magnetic anisotropy on the thermoelectric properties of magnetic molecule. Finite magnetic anisotropy gives rise to new qualitative features especially in the case of antiferromagnetic J, where a new sign change of thermopower occurs in the case of easy-plane magnetic anisotropy.
On the other hand, when the leads are ferromagnetic, depending on the spin relaxation time, spin accumulation may be generated in the contacts giving rise to the spin Seebeck effect. First, we have focused on the impact of the spin-resolved tunneling on the thermopower in the absence of spin accumulation. We have shown that in the case of ferromagnetic exchange interaction there is an additional sign change in the temperature dependence of the thermopower. On the other hand, the low-temperature Seebeck effect has been found to strongly depend on the position of the molecule's orbital level, which conditioned the strength of the exchange field in the molecule. Generally, the Seebeck effect becomes quenched once the temperature gets smaller than the corresponding exchange field.
Finally, we have assumed that the spin relaxation time in the leads is long, so that a spin bias can be generated in the system, resulting in the spin Seebeck effect. We have shown that for relatively low values of molecule's exchange interaction the spin Seebeck effect changes sign only when tuning the orbital level across the particlehole symmetry point. This is associated with the fact, that the sign of the spin thermopower is conditioned by the sign of the exchange field, which changes only when crossing ε = −U/2 . However, if the molecule's antiferromagnetic exchange interaction becomes larger than the Kondo temperature, the temperature dependence of the spin Seebeck effect can exhibit a sign change. This feature is associated with the interplay of exchange field, antiferromagnetic interaction between the orbital level and molecule's core spin and the Kondo correlations. A similarly nontrivial behavior have been observed in the case of finite easy-plane magnetic anisotropy. We have found that the spin Seebeck effect exhibits then a nonmonotonic dependence on the magnitude of anisotropy. This effect is related to a revival of the Kondo effect when the anisotropy becomes large enough to bring about the two-fold degenerate ground state of the molecule, i.e. when the easy-plane anisotropy wins over both the exchange field as well as the molecule's antiferromagnetic exchange interaction.
As far as the experimental progress is concerned, the field of molecular spin caloritronics is rather at its initial stage of development. The model studied here can however also describe quantum dots or impurities coupled to a large spin. It is therefore worth mentioning that, recently, there have been successful measurements of thermopower of Kondo-correlated quantum dots 10,25,26 . Moreover, spin-resolved electronic transport properties of molecular junctions have already been extensively studied 60,63-65 . Therefore, it seems that all the necessary ingredients are at hand and, with the state-of-the-art apparatus, it should be possible to explore the effects presented in this paper. We do hope that our work will foster further experimental efforts in this direction.
Methods. The thermoelectric properties of the system in the linear response regime can be characterized by the Onsager integrals, L nσ , which depend on the spin-resolved transmission coefficient T σ (ω) . The transmission coefficient T σ (ω) can be related to the spin-dependent spectral function A σ (ω) using the relation where A σ (ω) = −Im[G R σ (ω)]/π and G R σ (ω) is the Fourier transform of the retarded Green's function of the molecule's orbital level, G R σ (t) = −iθ(t)�{d σ (t),d † σ (0)}� . The main task is thus to accurately determine the spectral function of the system. One of the most powerful methods in this regard is the Wilson's numerical renormalization group method 50,51,66,67 , which allows for nonperturbative treatment of all correlations in the system. In this method one performs a discretization of the conduction band with discretization parameter and, consecutively, www.nature.com/scientificreports/ a tridiagonalization of the Hamiltonian describing such discretized system is performed. Eventually, one obtains the following NRG Hamiltonian 50 in which the molecule is coupled to the first site of the tight-binding chain with exponentially decaying hoppings ξ n . Here, f † nσ creates a spin-σ electron at the Wilson site n and W denotes the band halfwidth, which is used as energy unit W ≡ 1 . This Hamiltonian is then solved in an iterative fashion by retaining a fixed number of states N K . While the kept states are used to construct the statespace for the next iteration, the states that are discarded during calculation play a vital role in the problem as these states are used to construct the complete many-body basis of the full NRG Hamiltonian 68 . The discarded states are then used for the calculation of quantities of interest with the aid of the full density matrix 69 . In our calculations we have exploited the Abelian symmetries for the system's spin zth component and charge. We have used = 2 and kept at least N K = 2000 states in the iteration. Moreover, we have determined the Onsager transport coefficients from the raw NRG data without the need of broadening the Dirac delta peaks 36  www.nature.com/scientificreports/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.