pith. sign in

arxiv: 1907.05739 · v1 · pith:MBF5G7C2new · submitted 2019-07-12 · ❄️ cond-mat.mtrl-sci · cond-mat.mes-hall

Spin Seebeck effect and phonon energy transfer in heterostructures containing layers of a normal metal and a ferromagnetic insulator

A. I. Bezuglyj , V. A. Shklovskij , V. V. Kruglyak , R. V. Vovk This is my paper

Pith reviewed 2026-05-24 22:33 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci cond-mat.mes-hall
keywords spin Seebeck effectphonon heat transferN/F heterostructuresBoltzmann equationinverse spin Hall voltageelectron-magnon temperature differenceferromagnetic insulator
0
0 comments X

The pith

The spin Seebeck voltage arises from the difference between electron temperature in the metal and magnon temperature in the insulator, both set by phonon heat flow.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper applies the Boltzmann kinetic equation to phonons in an N/F bilayer to find the steady-state temperatures when a gradient is imposed either by current in the normal metal or by external dielectrics at different temperatures. It shows that the inverse spin Hall voltage on the normal-metal side is fixed by the resulting Te minus Tm. The calculated voltage versus bath temperature and versus N and F thicknesses reproduces the trends seen in existing experiments.

Core claim

In the kinetic approach based on the Boltzmann equation for the phonon distribution function, the electron temperature Te in the N-layer and the magnon temperature Tm in the F-layer are calculated for heterostructures with current heating or external temperature differences, with the difference Te - Tm determining the ISHE voltage VISHE whose dependence on temperature and thicknesses agrees with experiment.

What carries the argument

Boltzmann-equation solution for the phonon distribution that transfers energy at the N/F interface to electrons and to magnons, fixing the steady-state Te and Tm.

If this is right

  • VISHE is proportional to the phonon-limited temperature offset between the two layers.
  • Increasing the thickness of either layer raises the thermal resistance and therefore changes the voltage at fixed heat flow.
  • Current heating of the N layer produces a different Te-Tm profile than uniform external heating of the whole bilayer.
  • The voltage falls with rising bath temperature because phonon scattering rates increase.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same phonon-balance method could be applied to trilayers or multilayers to predict how additional interfaces alter the effective Te-Tm difference.
  • If the model were applied to materials with strong magnon-phonon coupling inside the F layer, the assumption that Tm is uniform would need re-examination.
  • The framework supplies a parameter-free route to estimate the phonon contribution in other spin-caloritronic geometries once the interface conductances are known.

Load-bearing premise

The interface coupling of phonons to electrons and magnons, together with the bulk scattering rates, fully accounts for all energy transfer channels between the layers.

What would settle it

An experiment that measures VISHE while systematically varying the thicknesses of both layers and finds a thickness dependence that deviates from the phonon-resistance prediction would falsify the model.

Figures

Figures reproduced from arXiv: 1907.05739 by A. I. Bezuglyj, R. V. Vovk, V. A. Shklovskij, V. V. Kruglyak.

Figure 1
Figure 1. Figure 1: FIG. 1. Refraction and reflection of phonon modes at layer [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Refraction and reflection of phonon modes at in [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. The dependence of the difference between the electron [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. The dimensionless difference between the electron [PITH_FULL_IMAGE:figures/full_fig_p010_4.png] view at source ↗
read the original abstract

In the framework of the kinetic approach based on the Boltzmann equation for the phonon distribution function, we analyze phonon heat transfer in a heterostructure containing a layer of a normal metal ($ N $) and a layer of a ferromagnetic insulator ($ F $). Two realistic methods for creating a temperature gradient in such a heterostructure are considered: by heating of the $N$-layer by an electric current and by placing the $N/F$-bilayer between massive dielectrics with different temperatures. The electron temperature $ T_e $ in the $ N $-layer and the magnon temperature $ T_m $ in the $ F $-layer are calculated. The difference in these temperatures determines the voltage $ V_{ISHE} $ on the $ N $-layer in the Seebeck spin effect regime. The dependence of $ V_{ISHE} $ on the bath temperature and on the thickness of the $ N $ and $ F $ layers is compared with the available experimental data.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript develops a kinetic model based on the Boltzmann equation for the phonon distribution in N/F heterostructures. It treats two heating methods (Joule heating of the N layer by current, and placement between dielectrics at different temperatures), solves for the steady-state electron temperature Te in N and magnon temperature Tm in F, and asserts that the difference Te − Tm sets the inverse spin Hall voltage VISHE. The resulting dependence of VISHE on bath temperature and on N/F layer thicknesses is compared with existing experimental data.

Significance. If the interface phonon-electron and phonon-magnon couplings can be fixed by material parameters rather than adjusted to data, the calculation would supply a concrete microscopic route from phonon kinetics to the observed VISHE(T, dN, dF) trends. Such a derivation would be useful for interpreting spin-Seebeck experiments and for estimating the relative importance of phonon-mediated versus direct channels.

major comments (2)
  1. [model setup / collision integrals] The Boltzmann-equation treatment requires explicit interface scattering rates (or coupling constants) between phonons and electrons in N and between phonons and magnons in F. The abstract and model description give no indication whether these rates are computed from first-principles matrix elements or chosen to reproduce the cited experiments; if the latter, the claim that the kinetic treatment determines Te and Tm without additional fitting parameters is not supported.
  2. [results / comparison with experiment] The two heating scenarios are stated to be treated, yet the manuscript does not demonstrate that the same set of interface parameters reproduces both the current-heating and the dielectric-boundary-heating data sets. A single-parameter-set comparison would be required to substantiate that the phonon channel alone accounts for the observed thickness and temperature trends.
minor comments (2)
  1. [kinetic equation] Notation for the phonon distribution function and the form of the interface collision terms should be written explicitly (e.g., as an equation) rather than described only in words.
  2. [figures] The manuscript should state the numerical values or functional forms adopted for the layer thicknesses and the bath-temperature range when the VISHE curves are plotted.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of the manuscript and the constructive comments. We address each major comment below and indicate the revisions we will make.

read point-by-point responses

  1. Referee: [model setup / collision integrals] The Boltzmann-equation treatment requires explicit interface scattering rates (or coupling constants) between phonons and electrons in N and between phonons and magnons in F. The abstract and model description give no indication whether these rates are computed from first-principles matrix elements or chosen to reproduce the cited experiments; if the latter, the claim that the kinetic treatment determines Te and Tm without additional fitting parameters is not supported.

    Authors: The interface coupling constants are phenomenological parameters chosen to reproduce the cited experimental trends; they are not derived from first-principles matrix elements. The manuscript does not claim that the kinetic treatment is free of fitting parameters. The Boltzmann approach determines the steady-state phonon distribution and the resulting Te and Tm once the couplings are specified, but the values of the couplings are adjusted to data. We will revise the model section and abstract to state this explicitly and remove any implication that the model is parameter-free. revision: yes

  2. Referee: [results / comparison with experiment] The two heating scenarios are stated to be treated, yet the manuscript does not demonstrate that the same set of interface parameters reproduces both the current-heating and the dielectric-boundary-heating data sets. A single-parameter-set comparison would be required to substantiate that the phonon channel alone accounts for the observed thickness and temperature trends.

    Authors: The manuscript formulates the kinetic model for both heating methods but performs the explicit comparison with experiment only for the current-heating (Joule) case, as that matches the available experimental datasets referenced. The dielectric-boundary heating is solved theoretically but not matched to a separate experimental series with the identical parameter set. We agree that a unified demonstration would be stronger; we will add a clarifying paragraph noting the parameter choice and the limitation that the two experimental configurations are not compared with exactly the same couplings in the present work. revision: partial

Circularity Check

0 steps flagged

Boltzmann phonon model yields independent Te/Tm difference for VISHE without self-referential reduction

full rationale

The provided abstract and description present a kinetic Boltzmann treatment of phonons with interface couplings to electrons and magnons, solving for steady-state Te in N and Tm in F layers under two heating methods. VISHE is stated to be determined by the resulting Te-Tm difference, with model outputs then compared to external experimental data on bath temperature and layer thicknesses. No equations are exhibited that define VISHE directly from fitted conductances or that rename a fit as a prediction. No self-citations are quoted as load-bearing for uniqueness or ansatz. The derivation chain is therefore self-contained against the cited experiments and does not reduce to its inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review supplies no explicit free parameters, axioms, or invented entities; the central claim rests on the unstated assumption that the Boltzmann phonon treatment plus interface coupling fully determines Te and Tm.

pith-pipeline@v0.9.0 · 5733 in / 1064 out tokens · 17736 ms · 2026-05-24T22:33:09.812437+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

39 extracted references · 39 canonical work pages

  1. [1]

    150, 459 (2010)

    G.E.W.Bauer, A.H.MacDonald, and S.Maekawa, Solid State Commun. 150, 459 (2010)

    work page 2010

  2. [2]

    Bauer, Spin caloritronics

    G.E.W. Bauer, Spin caloritronics. In: Spin Current [Edited by Sadamichi Maekawa, Sergio O. Valenzuela, Eiji Saitoh, Takashi Kimura] Oxford University Press, 2012

    work page 2012

  3. [3]

    Bauer, E.Saitoh, B.J

    G.E.W. Bauer, E.Saitoh, B.J. Van Wees, Nature materi- als, 11, 391 (2012)

    work page 2012

  4. [4]

    Sci., 7, 885 (2014)

    S.R.Boona, R.C.Myers, and J.P.Heremans, Energy Env- iron. Sci., 7, 885 (2014)

    work page 2014

  5. [5]

    Uchida, S

    K. Uchida, S. Takahashi, K. Hari, J. Ieda, W. Kishibae, K. Ando, S. Maekawa and E. Saitoh, Nature, 455, 778 (2008)

    work page 2008

  6. [6]

    Bauer, S.Maekawa, and E

    K.Uchida, J.Xiao, H.Adachi, J.Ohe, S.Takahashi1, J.Ieda, T.Ota, Y.Kajiwara, H.Umezawa, H.Kawai, G.E.W. Bauer, S.Maekawa, and E. Saitoh, Nature mate- rials, 9, 894 (2010)

    work page 2010

  7. [7]

    K.Uchida, T.Nonaka, T.Ota, and E.Saitoh, Appl. Phys. Lett. 97, 262504 (2010)

    work page 2010

  8. [8]

    Serga, V.Lauer, E.Th.Papaioannou, B.Hillebrands, and V.I.Vasyuchka, Appl

    M.Agrawal, A.A. Serga, V.Lauer, E.Th.Papaioannou, B.Hillebrands, and V.I.Vasyuchka, Appl. Phys. Lett. 105, 092404 (2014)

    work page 2014

  9. [9]

    M.Agrawal, V.I.Vasyuchka, A.A.Serga, A.Kirihara, P.Pirro, T.Langner, M.B.Jungfleisch, A.V.Chumak, E.Th.Papaioannou, and B.Hillebrands, Phys. Rev. B 89 , 224414 (2014)

    work page 2014

  10. [10]

    M.Schreier, F.Kramer, H.Huebl, S.Geprags, R.Gross, S.T.B.Goennenwein, T.Noack, T.Langner, A.A.Serga, B.Hillebrands, and V.I.Vasyuchka, Phys. Rev. B 93 , 224430 (2016)

    work page 2016

  11. [11]

    Azevedo, L

    A. Azevedo, L. H. Vilela Leao, R. L. Rodriguez-Suarez, A. B. Oliveira, and S. M. Rezende, J. Appl. Phys. 97, 10C715 (2005)

    work page 2005

  12. [12]

    E.-J.Guo, J.Cramer, A.Kehlberger, C.A.Ferguson, D.A.MacLaren, G.Jakob, and M.Klaui, Phys. Rev. X 6 , 031012 (2016)

    work page 2016

  13. [13]

    Xiao, G.E.W.Bauer, R.Gross, and S.T.B

    M.Schreier, A.Kamra, M.Weiler, J. Xiao, G.E.W.Bauer, R.Gross, and S.T.B. Goennenwein, Phys. Rev. B 88 , 094410 (2013)

    work page 2013

  14. [14]

    S.Y.Huang, W.G.Wang, S.F.Lee, J.Kwo, and C.L.Chien, Phys. Rev. Lett. 107, 216604 (2011)

    work page 2011

  15. [15]

    M.Schreier, N.Roschewsky, E.Dobler, S.Meyer, H.Hueb l, R.Gross, and S.T.B.Goennenwein, Appl. Phys. Lett. 103, 242404 (2013)

    work page 2013

  16. [16]

    W. X. Wang, S. H. Wang, L. K. Zou, J. W. Cai, Z. G. Sun, and J. R. Sun, Appl. Phys. Lett. 105, 182403 (2014)

    work page 2014

  17. [17]

    A.W.Little, Can. J. Phys. 37, 334 (1959)

    work page 1959

  18. [18]

    Shklovskij, V.V

    V.A. Shklovskij, V.V. Kruglyak, R.V. Vovk, and O.V.Dobrovolskiy, to be published

    work page

  19. [19]

    A. I. Bezuglyj and V. A. Shklovskij, Fiz. Nizk. Temp. 39, 459 (2013) [Low Temp. Phys. 39, 357 (2013)]

    work page 2013

  20. [20]

    A.I.Bezuglyj and V.A.Shklovskij, Phys. Rev. B 89 , 214303 (2014)

    work page 2014

  21. [21]

    Low Temp

    S.B.Kaplan, J. Low Temp. Phys. 37, 343 (1979)

    work page 1979

  22. [22]

    E.T.Swartz, R.O.Pohl, Rev. Mod. Phys. 61, 605 (1989)

    work page 1989

  23. [23]

    A.I.Akhieser and L.A.Shishkin, Zh. Eksp. Teor. Fiz. 34, 875 (1958) [Sov. Phys. JETP 34, 1267 (1958)]

    work page 1958

  24. [24]

    A.Shklovskii, Zh

    V. A.Shklovskii, Zh. Eksp. Teor. Fiz. 78, 1281 (1980) [Sov. Phys. JETP 51, 646 (1980)]

    work page 1980

  25. [25]

    A. I. Bezuglyj and V. A. Shklovskij, Zh. Eksp. Teor. Fiz. 12 111, 2106 (1997) [JETP 84, 1149 (1997)]

    work page 1997

  26. [26]

    J.Xiao, G.E.W.Bauer, K.Uchida, E.Saitoh, and S.Maekawa, Phys. Rev. B 81 , 214418 (2010)

    work page 2010

  27. [27]

    A. I. Bezuglyj and V. A. Shklovskij, J.Phys.:Condens.Matter, 29, 295001 (2018)

    work page 2018

  28. [28]

    J.Kimling, G.-M.Choi, J.T.Brangham, T.Matalla- Wagner, T.Huebner, T.Kuschel, F.Yang, and D.G.Cahill, Phys. Rev. Lett. 118, 057201 (2017)

    work page 2017

  29. [29]

    A.I.Bezuglyj, V.A.Shklovskij, V.V.Kruglyak, R.V.Vo vk to be published

    work page

  30. [30]

    Kikkawa, K.Uchida, Y.Shiomi, Z.Qiu, D.Hou, D.Tian, H.Nakayama, X.-F.Jin, and E.Saitoh, Phys

    T. Kikkawa, K.Uchida, Y.Shiomi, Z.Qiu, D.Hou, D.Tian, H.Nakayama, X.-F.Jin, and E.Saitoh, Phys. Rev. Lett. 110, 067207 (2013)

    work page 2013

  31. [31]

    S.M.Rezende, R.L.Rodriguez-Suarez, R.O.Cunha, A.R.Rodrigues, F.L.A.Machado, G.A.Fonseca Guerra, J.C.Lopez Ortiz, and A.Azevedo, Phys. Rev. B 89 , 014416 (2014)

    work page 2014

  32. [32]

    Y. H. Shen, X. S. Wang, and X. R. Wang, Phys. Rev. B 94, 014403 (2016)

    work page 2016

  33. [33]

    H.Adachi, K.Uchida, E.Saitoh, and S.Maekawa, Rep. Prog. Phys. 76, 036501 (2013)

    work page 2013

  34. [34]

    S.A.Bender and Y.Tserkovnyak, Phys. Rev. B 91, 140402 (2015)

    work page 2015

  35. [35]

    E.Saitoh, M.Ueda, H.Miyajima, and G.Tatara, Appl. Phys. Lett. 88, 182509 (2006)

    work page 2006

  36. [36]

    A.Prakash, B.Flebus, J.Brangham, F.Yang, Y.Tserkovnyak, and J.P.Heremans Phys. Rev. B 97, 020408(R) (2018)

    work page 2018

  37. [37]

    A. I. Bezuglyj and V. A. Shklovskij, Fiz. Nizk. Temp. 42, 809 (2016) [Low Temp. Phys. 42, 636 (2016)]

    work page 2016

  38. [38]

    E. M. Gershenzon, M. E. Gershenzon, G. N. Gol’tsman, A. D. Semenov, and A. V. Sergeev, Pis’ma Zh. Eksp. Teor. Fiz. 36, 241 (1982) [Sov. Phys. JETP Letters 36, 296 (1982)]

    work page 1982

  39. [39]

    Bar’yakhtar, and S.V

    Akhiezer, A.I., V.G. Bar’yakhtar, and S.V. Peletminsk ii, Spin Waves (North Holland, Amsterdam, 1968)

    work page 1968