Pith. sign in

REVIEW 2 major objections 6 minor 7 cited by

Sunny.jl: A Julia Package for Spin Dynamics

T0 review · 2 major / 6 minor · reviewed 2026-08-10 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Sunny.jl puts linear spin wave theory and classical spin dynamics on one SU(N) coherent-state footing, letting a single package predict scattering data for a broad class of quantum magnets.

desk verdict Sunny.jl is a genuinely useful and well-built software package, but the abstract oversells local-entanglement capabilities that the paper itself lists as future work. read the letter →

arxiv 2501.13095 v2 pith:4NM3HU3M submitted 2025-01-22 quant-ph cond-mat.str-elphysics.comp-ph

David Dahlbom , Hao Zhang , Cole Miles , Sam Quinn , Alin Niraula , Bhushan Thipe , Matthew Wilson , Sakib Matin
show 11 more authors
Het Mankad Steven Hahn Daniel Pajerowski Steve Johnston Zhentao Wang Harry Lane Ying Wai Li Xiaojian Bai Martin Mourigal Cristian D. Batista Kipton Barros
This is my paper · ORCID
classification quant-phcond-mat.str-elphysics.comp-ph
keywords JuliaquantummagnetismspindynamicsSU(N)coherentstateslinearwavetheoryclassicaldynamicalstructurefactorscatteringexperiments
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

Sunny is a Julia package whose central claim is that the two standard workhorses of quantum magnetism—linear spin wave theory near zero temperature and classical spin dynamics at higher temperatures—can be placed on one footing and implemented in a single, user-friendly code. It does so by replacing the classical dipole on each lattice site with an SU(N) coherent state, so that the same system specification supports both semiclassical and classical solvers. The paper argues this matters because scattering experiments (neutron and X-ray) measure the dynamical spin structure factor, and Sunny can compute that quantity for a broad class of spin models, including systems with strong spin-orbit coupling or local entanglement where ordinary dipole descriptions break down. If the claim holds, a researcher can go from a crystal structure and a hypothesized spin Hamiltonian to predicted scattering intensities—and back to the model—without writing bespoke numerical code.

What carries the argument

The SU(N) coherent state formalism. Instead of assigning a classical dipole to each site, each site carries a coherent state in an N-dimensional Hilbert space, generalizing the familiar spin-coherent-state picture. With N = 2s+1 it represents a spin-s ion including crystal-field levels; grouping sites can represent inter-site entanglement. It supplies the mathematical object on which both the generalized linear spin wave theory and the generalized classical Landau-Lifshitz dynamics are built, so a single model specification can be run in either regime.

What would settle it

Take a small spin-1 or spin-3/2 cluster (or a one-dimensional chain) with single-ion anisotropy and/or spin-orbit coupling, compute its dynamical spin structure factor with Sunny's :SUN linear spin wave and classical-dynamics modes, and compare against exact diagonalization at the same parameters; if the predicted mode positions and spectral weights deviate beyond controlled finite-size effects as N grows, the faithfulness of the SU(N) representation is refuted.

Watch

Extended reading notes

Core claim

On its own terms, the paper establishes Sunny as the first package that combines linear spin wave theory and classical spin dynamics and extends both through SU(N) coherent states. Setting a system's mode to :SUN promotes each site's state from a dipole to a full N-level quantum state, with N = 2s+1 for a spin-s ion; this captures crystal-field levels and, by grouping sites, local entanglement. The same formalism powers the LSWT solver and the classical equations of motion, and the package couples the latter to a thermal bath through Langevin dynamics. Alongside this, Sunny provides symmetry-guided model construction, optimizers and Monte Carlo samplers, and tools that turn the computed spin correlations into experimental S(q,ω).

Load-bearing premise

The software's distinctive reach depends on the assumption that SU(N) coherent states faithfully represent the local quantum state of the material—including crystal-field levels and spin-orbit or cluster entanglement—well enough that classical and semiclassical dynamics on them reproduce the real quantum dynamics.

Editorial extensions

If this is right

  • A single model definition in Sunny can be studied with LSWT near zero temperature and with classical dynamics at elevated or out-of-equilibrium conditions, so results from both regimes can be compared without reimplementing the Hamiltonian.
  • Materials with strong spin-orbit coupling, multipolar interactions, or local entanglement—where a simple dipole is inadequate—become accessible to semiclassical and classical simulation through the :SUN mode.
  • The dynamical spin structure factor computed by Sunny can be compared directly with neutron and RIXS scattering data, including instrument-specific region-of-interest integrations, making model validation a closed loop.
  • Symmetry-guided model construction and automatic propagation of interactions over symmetry-equivalent bonds lower the barrier to building realistic models from a crystallographic CIF file.
  • Because LSWT on large cells is accelerated by iterative matrix-vector products, chemically disordered or complex ordered magnets can be simulated at scales previously impractical.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the SU(N) formalism is as faithful as claimed, a natural testable extension is to use Sunny's classical dynamics to extract finite-temperature linewidths and compare them with the same material measured across a temperature sweep; the paper reports FeI2 as a showcase but does not claim this as a general proxy for quantum effects.
  • The same model-specification layer could serve as a common input for other solvers—such as exact diagonalization or tensor-network methods—so that classical, semiclassical, and fully quantum results are compared on identical Hamiltonians; Sunny's data-retrieval layer already pushes in this direction.
  • The planned RIXS observables would extend the scattering comparison beyond neutron data, which is a natural next testbed since RIXS measures S(q,ω)-like responses but with different matrix elements.
  • A broader implication is that the 'mode' abstraction (dipole versus SU(N)) may become a template for software that lets users choose the level of quantum description per site, potentially including mixed cluster sites, without changing the Hamiltonian specification.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 6 minor

Summary. The paper describes Sunny.jl, an open-source Julia package for simulating spin dynamics in magnetic materials. The authors report that Sunny supports symmetry-guided model construction from crystallographic data, classical ground-state optimization and Monte Carlo sampling, linear spin wave theory (LSWT), classical Landau-Lifshitz dynamics and its generalization to SU(N) coherent states, and calculations of the dynamical spin structure factor for comparison with neutron and X-ray scattering experiments. The central claimed novel contribution is that Sunny uniquely unifies LSWT and classical spin dynamics within an SU(N) coherent-state formalism, with claimed applicability to systems with strong spin-orbit coupling or local entanglement effects. The paper also documents a substantial list of experimental and theoretical studies that have used Sunny.

Significance. If the feature claims are accurate, Sunny is a valuable community resource for quantum magnetism: it integrates a broad workflow from crystal symmetry analysis through spin-model construction, simulation, and scattering-data comparison in a single, documented, Julia-based package. The manuscript's strongest evidence is the extensive list of independent applications and the claim of a public repository with correctness tests and tutorials. These are concrete signs of adoption and practical utility. However, the paper's most distinctive scientific claim—uniqueness and generality of the SU(N) coherent-state approach for local entanglement effects—is weakened by an internal inconsistency about what is already implemented versus planned. Because the manuscript is a software-description paper rather than a methods paper, its soundness largely rests on the accuracy and scope of its feature claims, which need to be clarified before the paper can be accepted.

major comments (2)
  1. [Abstract, 'SU(N) Formalism and System Modes', and 'Sunny as a Platform for Future Developments'] The abstract and summary state that the SU(N) coherent-state approach is 'useful for studying systems exhibiting strong spin-orbit coupling or local entanglement effects,' and the section 'SU(N) Formalism and System Modes' says the formalism 'can also be adapted to model local entanglement effects... within a cluster of spins on different sites.' However, the later section 'Sunny as a Platform for Future Developments' lists 'the modeling of local entanglement effects generated by spin-orbit coupling or strongly coupled clusters of spins' as a current development effort. This is an internal inconsistency: a capability presented as an available, unique feature is simultaneously described as future work. The manuscript should state clearly whether the released version of Sunny supports cluster/supersite SU(N) states with local dimension larger than 2s+1. If it does not, the claims should be narrowed to single-site SU(N) coherent states for crystal-field levels (N=2s+1), and the cluster-based spin-orbital entanglement extension should be explicitly labeled as planned work.
  2. [Statement of need] The paper claims 'Sunny is unique in offering both approaches and generalizing them through a formalism based on SU(N) coherent states.' This uniqueness claim is not substantiated by a comparative feature table or by citations demonstrating the absence of other SU(N)-based LSWT and classical-dynamics codes. Given that the previous major comment identifies an unimplemented component (local entanglement in clusters) within the claimed scope, the uniqueness claim is currently over-scoped. Please either provide concrete comparative evidence (for example, a feature matrix covering SpinW, SpinWaveGenie, pyLiSW, UppASD, Spirit, and any SU(N)-capable codes) or restrict the claim to the specific combination of features that are demonstrably present in the released version.
minor comments (6)
  1. [Affiliation list] The affiliation for Sam Quinn contains a typo: 'Univeriy of California' should be 'University of California'.
  2. [References] Several references are placeholders rather than formatted citations: '(2024a)', '(2024b)', '(2024c)', '(2024d)', '(2024e)', '(2024f)', and 'weber:2016'. These need to be replaced with proper author-year entries.
  3. [Figure 2] The qualitative comparison in Figure 2 would be more informative with a shared color scale and an explicit statement of the intensity normalization used for the experimental and calculated panels.
  4. [Acknowledgements] The sentence 'The data shown in Figure 2 was collected at the the Spallation Neutron Source' contains a duplicated article 'the'.
  5. [Feature Overview, 'Symmetry analysis'] The word 'inhomogenous' should be spelled 'inhomogeneous' in the sentence describing models that may be made 'inhomogenous'.
  6. [General] The manuscript does not specify the reviewed version of Sunny or its archive DOI; please include the exact version number and a stable software archive reference in the metadata.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity; self-citations document prior theory and adoption, and no claim reduces to its own input.

full rationale

This is a software description paper, not a derivation: the central claim is that Sunny implements a set of numerical methods and that the implementation works. The SU(N) coherent-state formalism is imported from prior published work, including Muniz et al. (2014) and H. Zhang & Batista (2021), some of whose authors overlap with this paper. That reliance is legitimate foundation rather than circularity: the formalism is not being derived here, and the package's correctness is supported by an independent test suite, documented tutorials, and comparisons to experimental scattering data (e.g., FeI2 in Figure 2). The extended list of publications using Sunny includes self-citations, but these are evidence of adoption, not load-bearing arguments that make a prediction equivalent to a fit. The uniqueness claim ('Sunny is unique in offering both approaches and generalizing them through a formalism based on SU(N) coherent states') is a scoping statement about competing software, not a theorem imported from the authors' prior work. One internal inconsistency is worth noting but is not circularity: the Abstract and Feature Overview say the formalism 'can also be adapted to model local entanglement effects ... within a cluster of spins on different sites,' while the Future Developments section states that 'current efforts are directed at supporting ... the modeling of local entanglement effects generated by spin-orbit coupling or strongly coupled clusters of spins.' That discrepancy affects how broadly the abstract's 'local entanglement effects' claim is supported by the released software, but it does not make the paper's central claim reduce to its own input. No fitted parameter is renamed as a prediction, no result is defined in terms of the quantity it claims to predict, and no load-bearing uniqueness theorem is imported solely from the authors' own work. Accordingly, the circularity score is 0.

Assumptions & free parameters 0 free parameters · 3 assumptions · 0 invented entities

The central claim of the paper is that Sunny is a useful, unique software package. This rests on the validity of the underlying theoretical approximations (SU(N) coherent states, LSWT, classical dynamics), all taken from prior literature, and on the correctness of the implementation, which is supported by tests but not formally verified here.

assumptions (3)
  • domain assumption SU(N) coherent states faithfully represent the local quantum state of spins, crystal field levels, and local inter-site entanglement (Muniz et al. 2014; H. Zhang & Batista 2021).
    This is the theoretical foundation for Sunny's unique features, adopted from prior work without re-derivation in this paper. Its validity determines whether the software's central selling point is meaningful.
  • domain assumption Linear spin wave theory (LSWT) is a valid approximation for calculating S(q,ω) at low temperatures.
    Sunny implements LSWT following established literature; the paper does not justify the approximation, relying on standard practice.
  • domain assumption Classical Landau-Lifshitz dynamics and its SU(N) generalization capture the relevant spin dynamics at finite temperature and out of equilibrium.
    The software's classical dynamics solvers rest on this assumption, which is conventional in the field but not proven in the paper.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Sunny.jl: A Julia Package for Spin Dynamics." pith.science (2026). https://pith.science/paper/4NM3HU3M

@misc{pith2026250113095,
  author       = {Pith},
  title        = {Pith review of: Sunny.jl: A Julia Package for Spin Dynamics},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/4NM3HU3M}},
  note         = {Machine review of arXiv:2501.13095}
}
read the original abstract

Sunny is a Julia package designed to serve the needs of the quantum magnetism community. It supports the specification of a very broad class of spin models and a diverse suite of numerical solvers. These include powerful methods for simulating spin dynamics both in and out of equilibrium. Uniquely, it features a broad generalization of classical and semiclassical approaches to SU(N) coherent states, which is useful for studying systems exhibiting strong spin-orbit coupling or local entanglement effects. Sunny also offers a well-developed framework for calculating the dynamical spin structure factor, enabling direct comparison with scattering experiments. Ease of use is a priority, with tools for symmetry-guided modeling and interactive visualization.

Figures

Figures reproduced from arXiv: 2501.13095 by the authors.

Figure 1
Figure 1. a) Ground state of FeI2, found using Sunny’s minimize_energy! function and visualized with plot_spins. b) The crystal of FeI2 visualized with the view_crystal function. Hovering the cursor over a bond reveals the exchange interaction, if already assigned, or a general expression for all symmetry-allowed interactions. Visualization Both the symmetry analysis and data retrieval features of Sunny include 3D visualizati… view at source ↗
Figure 2
Figure 2. Left: Scattering intensities of FeI2 as measured on the SEQUOIA instrument at the Spallation Neutron Source, Oak Ridge National Laboratory (Bai et al., 2021). Right: Predicted scattering intensities calculated with Sunny’s SU(N) linear spin wave solver. The figure was generated with Sunny’s data retrieval and plotting functions. Optimization and Monte Carlo Tools Identifying a classical ground state is often the fir… view at source ↗

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 7 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Skyrmions of Frustrated Quantum Dimer Systems

    cond-mat.str-el 2025-06 conditional novelty 7.0 of 10

    A classical-limit study of a frustrated spin-dimer bilayer with only isotropic Heisenberg exchanges finds two field-induced CP^3 skyrmion crystal phases at zero temperature.

  2. Magnetic Order in bilayer Ruddlesden-Popper Nickelates

    cond-mat.supr-con 2026-07 conditional novelty 6.0 of 10

    Combining superexchange with RKKY interactions between orbital-selective local moments reproduces the (π/2,π/2) magnetic order and ~80 meV spin excitations of bilayer nickelate La₃Ni₂O₇.

  3. Anisotropic magnons in a layered honeycomb ferromagnet

    cond-mat.str-el 2026-02 conditional novelty 6.0 of 10

    Detailed neutron spectroscopy shows the reported Dirac magnon gap in CrSiTe3 shrinks with integration volume and is likely a resolution artifact, while revealing frustrated out-of-plane couplings and unexplained aniso...

  4. Altermagnetic and dipolar splitting of magnons in FeF$_2$

    cond-mat.str-el 2026-01 conditional novelty 6.0 of 10

    In FeF2, dipolar interactions dominate the observed magnon splitting; a small altermagnetic splitting of about 35 μeV is inferred from line broadening at half-integer L.

  5. Dynamics of Quantum Chiral Solitons

    cond-mat.str-el 2025-12 conditional novelty 6.0 of 10

    Chiral solitons in a one-dimensional spin chain can be quantized on the lattice, producing fermionic bands whose nearest-neighbor hopping changes sign between integer and half-integer spin.

  6. Microscopic Origin of Reduced Magnetic Order in a Frustrated Metal

    cond-mat.str-el 2025-02 conditional novelty 6.0 of 10

    In the frustrated fcc metal HoInCu4, measured exchanges J1 and J2 with J2/J1 = 0.45(5) match between two independent neutron experiments, and the zero-field state shows overdamped excitations plus a ~30% quantum-drive...

  7. Wannier states and spin supersolid physics in the triangular antiferromagnet K$_2$Co(SeO$_3$)$_2$

    cond-mat.str-el 2024-12 conditional novelty 6.0 of 10

    K2Co(SeO3)2 shows Wannier entropy and continuum spin excitations consistent with a spin supersolid described by a triangular-lattice XXZ model, as supported by neutron scattering and quantum Monte Carlo.

Reference graph

Works this paper leans on

13 extracted references · 12 canonical work pages · cited by 7 Pith papers

  1. [1]

    In GitHub repository

    (2024d). In GitHub repository. GitHub. https://github.com/SunnySuite/Sunny.jl (2024a). In GitHub repository. GitHub. https://github.com/SpinWaveGenie/SpinWaveGenie (2024c). In GitHub. GitHub. https://github.com/UppASD/UppASD (2024b). In GitHub repository. GitHub. https://github.com/bingli621/pyLiSW (2024f). In GitHub repository. GitHub. https://github.com...

    work page 2024

  2. [5]

    Renormalized Classical Theory of Quantum Magnets

    Journal of Open Source Software , ¿VOL?(¿ISSUE?), ¿PAGE?. https://doi.org/DOIunavailable 5 Dahlbom, D. A., Thomas, J., Johnston, S., Barros, K., & Batista, C. D. (2024). Classical dynamics of the antiferromagnetic Heisenberg S=1/2 spin ladder.Physical Review B , 110(10), 104403. Dahlbom, D. A., Zhang, H., Laraib, Z., Pajerowski, D. M., Barros, K., & Batis...

    work page Pith review arXiv 2024

  3. [6]

    Evans, R. F. L., Fan, W. J., Chureemart, P., Ostler, T. A., Ellis, M. O. A., & Chantrell, R. W. (2014). Atomistic spin model simulations of magnetic nanomaterials.Journal of Physics: Condensed Matter, 26(10), 103202. Hall, S. R., Allen, F. H., & Brown, I. D. (1991). The crystallographic information file (CIF): A new standard archive file for crystallograp...

    work page 2014

  4. [7]

    https://doi.org/DOIunavailable 6 Nagl, J., Flavián, D., Hayashida, S., Povarov, K

    Journal of Open Source Software , ¿VOL?(¿ISSUE?), ¿PAGE?. https://doi.org/DOIunavailable 6 Nagl, J., Flavián, D., Hayashida, S., Povarov, K. Y., Yan, M., Murai, N., Ohira-Kawamura, S., Simutis, G., Hicken, T. J., Luetkens, H., & others. (2024). Excitation spectrum and spin Hamiltonian of the frustrated quantum Ising magnet Pr3BWO9. Physical Review Researc...

    work page 2024

  5. [48]

    D., & Park, J.-G

    Park, P., Cho, W., Kim, C., An, Y., Iida, K., Kajimoto, R., Matin, S., Zhang, S.-S., Batista, C. D., & Park, J.-G. (2024). Contrasting dynamical properties of single-Q and triple-Q magnetic orderings in a triangular lattice antiferromagnet.arXiv Preprint arXiv:2410.02180. Park, P., Cho, W., Kim, C., An, Y., Kang, Y.-G., Avdeev, M., Sibille, R., Iida, K., ...

    arXiv 2024

  6. [145]

    D., Cheong, S.-W., Park, J.-H., Brosha, E., & others

    Lee, M., Schönemann, R., Zhang, H., Dahlbom, D., Jang, T.-H., Do, S.-H., Christianson, A. D., Cheong, S.-W., Park, J.-H., Brosha, E., & others. (2023). Field-induced spin level crossings within a quasi-XY antiferromagnetic state in Ba2FeSi2O7. Physical Review B , 107(14), 144427. Legros, A., Zhang, S.-S., Bai, X., Zhang, H., Dun, Z., Phelan, W. A., Batist...

    work page 2023

  7. [2050]

    R., Cage, G., Fortney, J., Granroth, G

    Watson, G. R., Cage, G., Fortney, J., Granroth, G. E., Hughes, H., Maier, T., McDonnell, M., Ramirez-Cuesta, A., Smith, R., Yakubov, S., & others. (2022). Calvera: A platform for the interpretation and analysis of neutron scattering data.Smoky Mountains Computational Sciences and Engineering Conference , 137–154. Zhang, H., & Batista, C. D. (2021). Classi...

    work page 2022

  8. [2607]

    Togo, A., Shinohara, K., & Tanaka, I. (2024). Spglib: A software library for crystal symmetry search. Science and Technology of Advanced Materials: Methods , 4(1), 2384822. Toth, S., & Lake, B. (2015). Linear spin wave theory for single-Q incommensurate magnetic structures. Journal of Physics: Condensed Matter , 27(16), 166002. Wang, F., & Landau, D. P. (...

    work page 2024

Show all 13 references
  1. [3349]

    A., Williams, T

    Do, S.-H., Zhang, H., Dahlbom, D. A., Williams, T. J., Garlea, V. O., Hong, T., Jang, T.-H., Cheong, S.-W., Park, J.-H., Barros, K., & others. (2023). Understanding temperature- dependent SU(3) spin dynamics in the S=1 antiferromagnet Ba2FeSi2O7. Npj Quantum Materials, 8(1),

    2023

  2. [3626]

    Sunny.jl: A Julia Package for Spin Dynamics

    ¿citation_author?, (2025). Sunny.jl: A Julia Package for Spin Dynamics

    2025

  3. [4199]

    A., Brooks, F

    Dahlbom, D. A., Brooks, F. T., Wilson, M. S., Chi, S., Kolesnikov, A. I., Stone, M. B., Cao, H., Li, Y. W., Barros, K., Mourigal, M., & others. (2024). Quantum-to-classical crossover in generalized spin systems: Temperature-dependent spin dynamics ofFeI2. Physical Review B, 10...

    2024

  4. [7264]

    M., May, A

    Park, P., Sala, G., Pajerowski, D. M., May, A. F., Kolopus, J. A., Dahlbom, D. A., Stone, M. B., Halász, G. B., & Christianson, A. D. (2024). Quantum and classical spin dynamics across temperature scales in theS=1/2 Heisenberg antiferromagnet. Physical Review Research, 6(3), 0...

    2024

  5. [8346]

    A., May, A

    Park, P., Ghioldi, E. A., May, A. F., Kolopus, J. A., Podlesnyak, A. A., Calder, S., Paddison, J. A. M., Trumper, A. E., Manuel, L. O., Batista, C. D., & others. (2024). Anomalous continuum scattering and higher-order van Hove singularity in the strongly anisotropic S=1/2 tria...

    2024

Pith tools

Reviewed August 10, 2026 · model on record in the stance chip above.