pith. sign in

arxiv: 2603.20106 · v2 · submitted 2026-03-20 · ❄️ cond-mat.mes-hall · cond-mat.stat-mech· physics.comp-ph· quant-ph

Micromagnetic Modeling of Surface Acoustic Wave Driven Dynamics: Interplay of Strain, Magnetorotation, and Magnetic Anisotropy

Florian Millo , Pauline Rovillain , Massimiliano Marangolo , Daniel Stoeffler This is my paper

Pith reviewed 2026-05-15 08:02 UTC · model grok-4.3

classification ❄️ cond-mat.mes-hall cond-mat.stat-mechphysics.comp-phquant-ph
keywords surface acoustic wavesspin wavesmagnetoelastic couplingmicromagnetic modelingCoFeBuniaxial anisotropyparallel propagationmagnetorotation
0
0 comments X

The pith

Anisotropy orientation tunes the resonant coupling between surface acoustic waves and spin waves in parallel geometry.

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

Micromagnetic simulations implement the full SAW magnetoacoustic field, including every strain and lattice-rotation contribution, to model the interaction with spin waves in a realistic CoFeB film that carries only weak in-plane uniaxial anisotropy. The analysis focuses on the parallel configuration in which the SAW travels along the external magnetic field, a geometry favored for device integration. The simulations show that the direction of the anisotropy axis directly controls the strength of the resonant interaction in this parallel case. A reader would care because this supplies a concrete handle for designing efficient magnetoacoustic devices without requiring changes to the applied field or film thickness.

Core claim

In micromagnetic simulations of a CoFeB film, the complete implementation of the SAW-induced magnetoelastic field, encompassing both strain and magnetorotation terms, reveals that the orientation of the weak in-plane uniaxial anisotropy acts as a tuning parameter for the resonant coupling between SAW and spin waves propagating parallel to the applied magnetic field.

What carries the argument

The full SAW magnetoacoustic excitation field that incorporates all strain-induced and lattice-rotation contributions, implemented inside the micromagnetic model for a CoFeB film with weak in-plane uniaxial anisotropy.

Load-bearing premise

The micromagnetic model with every strain and lattice-rotation term included accurately represents the physical SAW-SW coupling inside a realistic CoFeB film that has only weak in-plane uniaxial anisotropy.

What would settle it

An experiment that measures the amplitude or frequency shift of the parallel SAW-SW resonance while rotating the direction of the film's uniaxial anisotropy axis relative to the SAW propagation direction.

Figures

Figures reproduced from arXiv: 2603.20106 by Daniel Stoeffler, Florian Millo, Massimiliano Marangolo, Pauline Rovillain.

Figure 1
Figure 1. Figure 1: a) Geometry of the problem and definitions. The slab has dimensions of ℓx = ℓy = λSAW, ℓz = 34 nm and is meshed into {Nx, Ny, Nz} = {128, 128, 1}. An in-plane static field B⃗0 is applied with an angle ψ with respect to ⃗kSAW = kSAW. ˆ⃗x. A weak uniaxial anisotropy field B⃗ u with an angle φu is introduced in the system. EA stands for easy axis. b) Symmetry of the SAW–SW coupling strength ∆P(ψ) at fSAW = 1.… view at source ↗
Figure 2
Figure 2. Figure 2: 2D maps of the SAW–SW coupling strength ∆P(B0, ψ) highlighting the interplay of anisotropy with magnetoacoustic excitation. Left column: εxx (& εxz). Right column: full magnetoacoustic drive including strain and rotation (εµν and ωµν). (a,b) In the absence of anisotropy Bu = 0. (c,d) Weak anisotropy Bu = 1.5 mT oriented at φu = 105◦ . The horizontal dashed line marks ψ = 180◦ . clude the shear strain εxz, … view at source ↗
Figure 3
Figure 3. Figure 3: 2D maps of the SAW–SW coupling strength ∆P(ψ, φu), where φu is varied in increments of 15◦ and ψ in increments of 1 ◦ , at fixed uniaxial anisotropy field Bu = 1.5 mT, including full magnetoacoustic drive, with strain and rotation (εµν and ωµν). Polynomial nth-order interpolation is applied between simulated φu-curves. Left column: B0 = 1.5 mT. Right column: B0 = 5 mT. (a,b) fSAW = 1.72 GHz. (c,d) fSAW = 3… view at source ↗
Figure 4
Figure 4. Figure 4: 2D maps of the SAW–SW coupling strength ∆P(fSAW, ψ) for Bu = 1.5 mT and φu = 120◦ . a) field strength B0 = 1.5 mT. b) B0 = 5 mT. Parameters used in this simulation are: | ⃗kSAW| = 2πfSAW/vSAW, full magnetoacoustic drive including strain and rotation (εµν and ωµν). The horizontal dashed line marks ψ = 180◦ . The arrows show the intersection at which the parallel SAW–SW interaction occurs. tion shifts both ∆… view at source ↗
read the original abstract

We study the coupling mechanism of surface acoustic waves (SAW) with spin waves (SW) using micromagnetic analysis. The SAW magnetoacoustic excitation field is fully implemented, i.e., all strain and lattice-rotation terms are included. A realistic CoFeB film with a weak in-plane uniaxial anisotropy is considered. We investigate the conditions for efficient SAW--SW coupling, with particular emphasis on the case where the SAW propagates parallel to the external magnetic field, a configuration of special interest for magnonic applications. Remarkably, we find that the anisotropy orientation serves as a knob to tune the parallel resonant interaction. Overall, this work provides a unified and practical picture of SAW--SW coupling in thin magnetized films.

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

3 major / 2 minor

Summary. The manuscript presents a micromagnetic modeling study of SAW-SW coupling in thin films. It fully implements the magnetoacoustic driving field with all strain and lattice-rotation contributions, models a realistic CoFeB film possessing weak in-plane uniaxial anisotropy, and examines coupling efficiency across field orientations. The central result is that the anisotropy axis orientation acts as a tunable parameter that enables resonant interaction even when the SAW propagates parallel to the external field.

Significance. If the numerical implementation is shown to be free of discretization artifacts, the work supplies a practical design rule for magnonic devices that exploit acoustic driving, particularly in geometries where parallel coupling was previously symmetry-forbidden. The emphasis on a complete strain-plus-rotation torque model and the identification of anisotropy as a control knob are potentially useful for device engineering.

major comments (3)
  1. [§3] §3 (Micromagnetic model): the manuscript states that all strain and magnetorotation terms are included, yet provides no limiting-case analytic benchmark (e.g., zero-anisotropy dispersion or torque balance) against which the simulated resonance shift for parallel geometry can be compared. Without this cross-check it is unclear whether the reported tuning originates from the intended magnetoelastic physics or from numerical choices.
  2. [§4.2] §4.2 (Parallel-coupling results): the claim that anisotropy orientation 'serves as a knob' rests on a single set of CoFeB parameters; no sensitivity analysis or variation of the uniaxial anisotropy strength (the only free parameter listed) is shown to confirm that the effect survives changes in damping, mesh size, or film thickness.
  3. [Fig. 5] Fig. 5 (resonance spectra): the parallel-geometry spectra lack error bars or ensemble averaging over initial conditions, making it impossible to judge whether the observed frequency shift exceeds numerical noise.
minor comments (2)
  1. [§2] Notation for the magnetorotation vector is introduced without an explicit definition of its relation to the lattice displacement gradient; a short appendix equation would improve clarity.
  2. [Abstract] The abstract states 'all strain and lattice-rotation terms are included' but the main text does not list the explicit torque expressions; adding them would aid reproducibility.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive and detailed report. We have revised the manuscript to address the concerns and provide additional validation and robustness checks as outlined below.

read point-by-point responses

  1. Referee: [§3] §3 (Micromagnetic model): the manuscript states that all strain and magnetorotation terms are included, yet provides no limiting-case analytic benchmark (e.g., zero-anisotropy dispersion or torque balance) against which the simulated resonance shift for parallel geometry can be compared. Without this cross-check it is unclear whether the reported tuning originates from the intended magnetoelastic physics or from numerical choices.

    Authors: We agree that an explicit analytic benchmark strengthens the validation. In the revised manuscript we have added a new paragraph in §3 deriving the zero-anisotropy resonance condition from the magnetoelastic torque balance (including both strain and rotation terms) and directly comparing the analytic frequencies to micromagnetic results for the parallel geometry. The close agreement confirms that the observed tuning originates from the intended physics. revision: yes

  2. Referee: [§4.2] §4.2 (Parallel-coupling results): the claim that anisotropy orientation 'serves as a knob' rests on a single set of CoFeB parameters; no sensitivity analysis or variation of the uniaxial anisotropy strength (the only free parameter listed) is shown to confirm that the effect survives changes in damping, mesh size, or film thickness.

    Authors: We have performed additional simulations varying the uniaxial anisotropy constant (over ±50% of the nominal value), Gilbert damping, mesh cell size, and film thickness. The anisotropy-orientation tuning of the parallel resonance persists in all cases, with only quantitative shifts in coupling strength. These results are summarized in a new supplementary figure and briefly discussed in the revised §4.2. revision: yes

  3. Referee: [Fig. 5] Fig. 5 (resonance spectra): the parallel-geometry spectra lack error bars or ensemble averaging over initial conditions, making it impossible to judge whether the observed frequency shift exceeds numerical noise.

    Authors: We have recomputed the spectra in Fig. 5 using ensemble averaging over ten independent runs, each started from a slightly perturbed initial magnetization (random fluctuations of order 1°). Standard-deviation error bars are now shown; the frequency shifts induced by anisotropy orientation remain well above the numerical variability. revision: yes

Circularity Check

0 steps flagged

No circularity detected in micromagnetic SAW-SW coupling simulations

full rationale

The paper implements standard micromagnetic equations with all strain and lattice-rotation terms for a realistic CoFeB film and reports simulation results on anisotropy-tuned parallel coupling. No step reduces by construction to a fitted input renamed as prediction, a self-citation chain, or an ansatz smuggled via prior work; the tuning effect is a direct numerical output from the implemented physical model using independent parameters. The derivation chain is self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The model rests on standard micromagnetic assumptions plus the specific choice of a weak uniaxial anisotropy for CoFeB; no new entities are introduced.

free parameters (1)
  • weak in-plane uniaxial anisotropy strength
    Chosen to represent a realistic CoFeB film; its value directly affects the reported tuning behavior.
axioms (1)
  • domain assumption Standard micromagnetic equations remain valid when all strain and lattice-rotation terms from the SAW are included
    Invoked throughout the modeling of magnetoacoustic excitation.

pith-pipeline@v0.9.0 · 5445 in / 1253 out tokens · 46844 ms · 2026-05-15T08:02:50.140252+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

37 extracted references · 37 canonical work pages

  1. [1]

    Kittel, Phys

    C. Kittel, Phys. Rev. 110, 1295 (1958)

    work page 1958

  2. [2]

    A. I. Akhiezer, V . G. B. Iakhtar, and S. V . Peletminskii, Soviet Physics JETP (1958)

    work page 1958

  3. [3]

    H. F. Tiersten, Journal of Mathematical Physics 5, 1298 (1964)

    work page 1964

  4. [4]

    Bömmel and K

    H. Bömmel and K. Dransfeld, Phys. Rev. Lett. 3, 83 (1959)

    work page 1959

  5. [5]

    Kobayashi, R

    T. Kobayashi, R. C. Barker, and A. Y elon,Phys. Rev. B 7, 3286 (1973)

    work page 1973

  6. [6]

    Flebus, D

    B. Flebus, D. Grundler, B. Rana, Y . Otani, I. Barsukov, A. Bar- man, G. Gubbiotti, P . Landeros, J. Akerman, U. Ebels, P . Pirro, V . E. Demidov, K. Schultheiss, G. Csaba, Q. Wang, F. Ciubo- taru, D. E. Nikonov, P . Che, R. Hertel, T. Ono, D. Afanasiev, J. Mentink, T. Rasing, B. Hillebrands, S. V . Kusminskiy, W. Zhang, C. R. Du, A. Finco, T. van der Sar...

    work page 2024

  7. [7]

    H. J. Krenner, P . V . Santos, C. Westerhausen, G. Ander- sson, A. N. Cleland, H. Sellier, S. Takada, C. Bäuerle, D. Wigger, T. Kuhn, P . Machnikowski, M. WeiSS, G. Moody, A. Hernández-Mínguez, S. Lazi, A. S. Kuznetsov, M. KüSS, M. Albrecht, M. Weiler, J. Puebla, Y . Hwang, Y . Otani, K. C. Balram, I.-T. Chen, K. Lai, M. Li, G. R. Nash, E. D. S. Nys- ten,...

    work page 2026

  8. [8]

    Rayleigh, Proceedings of the London Mathematical Society s1-17, 4 (1885)

    L. Rayleigh, Proceedings of the London Mathematical Society s1-17, 4 (1885)

    work page

  9. [9]

    Weiler, L

    M. Weiler, L. Dreher, C. Heeg, H. Huebl, R. Gross, M. S. Brandt, and S. T. B. Goennenwein, Phys. Rev. Lett. 106, 117601 (2011)

    work page 2011

  10. [10]

    B. A. Kalinikos and A. N. Slavin, J. Phys. C: Solid State Phys. 19, 7013 (1986)

    work page 1986

  11. [11]

    Dreher, M

    L. Dreher, M. Weiler, M. Pernpeintner, H. Huebl, R. Gross, M. S. Brandt, and S. T. B. Goennenwein, Phys. Rev. B 86, 134415 (2012)

    work page 2012

  12. [12]

    Thevenard, C

    L. Thevenard, C. Gourdon, J. Y . Prieur, H. J. V on Bardeleben, S. Vincent, L. Becerra, L. Largeau, and J.-Y . Duquesne, Phys. Rev. B 90, 094401 (2014)

    work page 2014

  13. [13]

    Labanowski, A

    D. Labanowski, A. Jung, and S. Salahuddin, Appl. Phys. Lett. 108, 022905 (2016)

    work page 2016

  14. [14]

    Sasaki, Y

    R. Sasaki, Y . Nii, Y . Iguchi, and Y . Onose, Phys. Rev. B 95, 020407 (2017)

    work page 2017

  15. [15]

    M. Xu, K. Y amamoto, J. Puebla, K. Baumgaertl, B. Rana, K. Miura, H. Takahashi, D. Grundler, S. Maekawa, and Y . Otani,Sci. Adv. 6, eabb1724 (2020)

    work page 2020

  16. [16]

    Hernández-Mínguez, F

    A. Hernández-Mínguez, F. Macià, J. M. Hernàndez, J. Herfort, and P . V . Santos,Phys. Rev. Applied 13, 044018 (2020)

    work page 2020

  17. [17]

    Puebla, M

    J. Puebla, M. Xu, B. Rana, K. Y amamoto, S. Maekawa, and Y . Otani, Journal of Physics D: Applied Physics 53, 264002 (2020)

    work page 2020

  18. [18]

    N. K. P . Babu, A. Trzaskowska, P . Graczyk, G. Centaa, S. Mieszczak, H. Gowiski, M. Zdunek, S. Mielcarek, and J. W. Kos, Nano Lett. 21, 946 (2021)

    work page 2021

  19. [19]

    R. Gao, Y . Y e, H. Wu, X. Li, X. Liu, J. Wang, and Q. Liu, Applied Physics Letters 121, 042401 (2022)

    work page 2022

  20. [20]

    F. Chen, J. Lu, S. Liang, Y . Otani, X. Y ang, Y . Zhang, and W. Luo, Journal of Alloys and Compounds 980, 173635 (2024). 7

    work page 2024

  21. [21]

    Lopes Seeger, L

    R. Lopes Seeger, L. La Spina, V . Laude, F. Millo, A. Barta- syte, S. Margueron, A. Solignac, G. De Loubens, L. Thevenard, C. Gourdon, C. Chappert, and T. Devolder, Phys. Rev. B 109, 104416 (2024)

    work page 2024

  22. [22]

    Millo, R

    F. Millo, R. L. Seeger, C. Chappert, A. Solignac, and T. De- volder, AIP Advances 15, 045230 (2025)

    work page 2025

  23. [23]

    P . G. Gowtham, T. Moriyama, D. C. Ralph, and R. A. Buhrman, Journal of Applied Physics 118, 233910 (2015)

    work page 2015

  24. [24]

    Rovillain, J.-Y

    P . Rovillain, J.-Y . Duquesne, L. Christienne, M. Eddrief, M. G. Pini, A. Rettori, S. Tacchi, and M. Marangolo, Phys. Rev. Ap- plied 18, 064043 (2022)

    work page 2022

  25. [25]

    Y amamoto, M

    K. Y amamoto, M. Xu, J. Puebla, Y . Otani, and S. Maekawa, Journal of Magnetism and Magnetic Materials 545, 168672 (2022)

    work page 2022

  26. [26]

    W. Hu, M. Huang, H. Xie, H. Zhang, and F. Bai, Phys. Rev. Appl. 19, 014010 (2023)

    work page 2023

  27. [27]

    Y amamoto and S

    K. Y amamoto and S. Maekawa, Annalen der Physik , 2300395 (2023)

    work page 2023

  28. [28]

    Huang, Y

    M. Huang, Y . Liu, W. Hu, Y . Wu, W. Wang, W. He, H. Zhang, and F. Bai, Phys. Rev. Applied 21, 014035 (2024)

    work page 2024

  29. [29]

    A. K. Vythelingum, T. Tremblais, D. B. Nguyen, D. Ourdani, E. Dandeu, V . Laude, C. Gourdon, and L. Thevenard, ACS Applied Electronic Materials 7, 4055 (2025)

    work page 2025

  30. [30]

    Sharma, L

    A. Sharma, L. Christienne, F. Millo, M. Eddrief, E. Dandeu, J.-Y . Duquesne, C. Gourdon, L. Thevenard, P . Rovillain, and M. Marangolo, Applied Physics Letters 128, 062404 (2026)

    work page 2026

  31. [31]

    V anderveken, J

    F. V anderveken, J. Mulkers, J. Leliaert, B. V an Waeyenberge, B. Sorée, O. Zografos, F. Ciubotaru, and C. Adelmann, Phys. Rev. B 103, 054439 (2021)

    work page 2021

  32. [32]

    Jander, P

    A. Jander, P . Dhagat, and C. Rivard,Journal of Applied Physics 137 (2025), 10.1063/5.0251252

    work page doi:10.1063/5.0251252 2025

  33. [33]

    V ansteenkiste, J

    A. V ansteenkiste, J. Leliaert, M. Dvornik, M. Helsen, F. Garcia- Sanchez, and B. V an Waeyenberge, AIP Advances 4, 107133 (2014)

    work page 2014

  34. [34]

    mumax+: ex- tensible gpu-accelerated micromagnetics and beyond,

    L. Moreels, I. Lateur, D. D. Gusem, J. Mulkers, J. Maes, M. V . Milosevic, J. Leliaert, and B. V . Waeyenberge, “mumax+: ex- tensible gpu-accelerated micromagnetics and beyond,” (2025), arXiv:2411.18194 [cond-mat.mes-hall]

    work page arXiv 2025

  35. [35]

    Ngouagnia Y emeli, L

    I. Ngouagnia Y emeli, L. Christienne, P . Rovillain, J.-Y . Duquesne, A. Anane, M. Marangolo, and D. Stoeffler, Jour- nal of Applied Physics 137, 153908 (2025)

    work page 2025

  36. [36]

    We verified that discretizing the thickness with Nz = 8 cells does not improve the result as compared to Nz = 1

    work page

  37. [37]

    The choice of SAW frequency ( fSA W = 1 .72 GHz) follows from experimental measurements in ref. [ 22]. Appendix A: Magnetic power absorption used in MuMax3 simulations The SAW attenuation into the magnetic film is quantified by the magnetic power absorption per unit area [W/m 2]. In linear response, for a harmonic excitation at frequency f = fSA W, the magn...

    work page