{"id":"b5f4b9f4-e0ec-4074-8094-89ea81e4e25f","arxiv_id":"2501.13095","paper_version":2,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Sunny.jl unifies linear spin wave theory and classical spin dynamics with SU(N) coherent states to simulate spin dynamics and scattering spectra.","lead":"Sunny.jl is a new Julia software package that helps scientists simulate the magnetic behavior of materials, including quantum effects. It combines several calculation methods into one tool so researchers can compare simulations directly with neutron scattering experiments.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Abstract claims SU(N) local-entanglement modeling, but future-work section defers it; uniqueness over-scoped.","rationale":"The reader's identified weakest assumption concerns the physical fidelity of the SU(N) coherent-state formalism, which is a reasonable domain assumption. My concern is more concrete: the paper's own future-work section indicates that a key advertised capability—modeling local entanglement within strongly coupled clusters—is not yet supported by the released software. This is a missing-support issue located in the manuscript text, not a disagreement with the overarching value of the package. The software is otherwise well supported by tests, documentation, and a substantial list of applications, so a rejection is not warranted. However, the abstract and the central uniqueness claim overstate the current scope. Acceptance should be conditional on revising these claims to distinguish shipped features (single-site SU(N), crystal-field levels, LSWT, classical dynamics) from planned extensions (cluster-based local entanglement, non-perturbative corrections). The concrete check of whether a composite local Hilbert space is supported would settle whether the condition is already met or whether a text revision is required.","tokens_in":7150,"tokens_out":9813,"duration_ms":107339,"concrete_test":"Search the Sunny.jl documentation and test suite for any 'system mode' or 'site' with a composite local Hilbert space (e.g., a spin-1/2 dimer defined as a single site with local dimension 4). Run a minimal dimer example in :SUN mode and verify that the API accepts the 4-dimensional local space, and that the resulting SU(N) dynamics or LSWT spectrum differ from two independent SU(2) sites. If no such mode exists or the API raises an error, the local-entanglement cluster capability is not implemented, and the abstract and uniqueness claim must be revised.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The paper's central claim is that Sunny uniquely generalizes LSWT and classical spin dynamics through SU(N) coherent states, 'useful for studying systems exhibiting strong spin-orbit coupling or local entanglement effects' (Abstract). This is supported by the 'SU(N) Formalism and System Modes' section, which 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' 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.' This is a self-referential inconsistency: a capability advertised as an available, unique feature is simultaneously listed as future work. The paper provides no documentation, test, or tutorial demonstrating a cluster/supersite mode with local dimension larger than 2s+1. The released software appears to support SU(N) for single-ion crystal-field states (N=2s+1), but the cluster-based extension for inter-site local entanglement seems unimplemented. If so, the 'local entanglement effects' portion of the strongest claim is unsupported, and the uniqueness claim should be narrowed to single-site SU(N) coherent states until the cluster functionality ships.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":7340,"tokens_out":3745,"duration_ms":40242,"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":[{"comment":"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.","section":"Abstract, 'SU(N) Formalism and System Modes', and 'Sunny as a Platform for Future Developments'"},{"comment":"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.","section":"Statement of need"}],"minor_comments":[{"comment":"The affiliation for Sam Quinn contains a typo: 'Univeriy of California' should be 'University of California'.","section":"Affiliation list"},{"comment":"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.","section":"References"},{"comment":"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.","section":"Figure 2"},{"comment":"The sentence 'The data shown in Figure 2 was collected at the the Spallation Neutron Source' contains a duplicated article 'the'.","section":"Acknowledgements"},{"comment":"The word 'inhomogenous' should be spelled 'inhomogeneous' in the sentence describing models that may be made 'inhomogenous'.","section":"Feature Overview, 'Symmetry analysis'"},{"comment":"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.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The paper is appropriate for JOSS and the software appears to be real, used, and maintained. The central issue is not the existence of the package but the precision of its advertised capabilities: the cluster/spin-orbital local-entanglement functionality is claimed as available in the abstract and formalism sections while being listed as future work later in the paper. If the authors clarify the release scope and adjust the uniqueness claim accordingly, I would expect the paper to be acceptable for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know two things about this paper. First, the software is real and useful: Sunny.jl integrates LSWT, classical spin dynamics, SU(N) coherent states, Monte Carlo, and symmetry analysis into one Julia package, with a public repo, tests, tutorials, and a long list of experimental papers already using it. That combination genuinely is new. Second, the stress-test note is right about the overclaim: the abstract says the SU(N) formalism is \"useful for studying systems exhibiting strong spin-orbit coupling or local entanglement effects,\" and the formalism section says it \"can also be adapted\" to local entanglement within clusters, but the \"Future Developments\" section explicitly lists \"modeling of local entanglement effects generated by spin-orbit coupling or strongly coupled clusters of spins\" as current work. No documentation, test, or tutorial demonstrates a cluster/supersite mode with local dimension larger than 2s+1. So the uniqueness claim should be narrowed to single-site SU(N) coherent states—which already covers crystal-field levels and single-ion spin-orbit physics—and the cluster extension should be labeled as planned, not shipped.\n\nWhat the paper does well: it gives a clear overview of the package's capabilities, situates it against existing codes like SpinW and Spirit, and backs the central claims with a public repository and correctness tests. The qualitative comparison to FeI2 data in Figure 2 is suggestive rather than quantitative, but for a JOSS software paper that is minor. The citation list is long and includes many self-citations, but those are mostly prior theoretical papers (Muniz 2014, Zhang & Batista 2021) and applications that genuinely used the code; that is evidence of adoption, not circularity.\n\nThe soft spots are the overclaim I already described and a somewhat vague treatment of what \"local entanglement\" means operationally. The formalism section is honest that the cluster adaptation is theoretical, but the abstract and the uniqueness sentence go further than the code delivers. That is a fixable wording issue, not a fundamental flaw.\n\nWho should read this: anyone in quantum magnetism who wants to compute S(q,ω) without rebuilding the pipeline from scratch, and anyone comparing semiclassical and SU(N) methods on real materials. It deserves a serious referee, and I would send it to review but ask for a revision that aligns the abstract and the uniqueness claim with the actually implemented features. I would cite the package if I were doing scattering-modeling work.\n\nRecommendation: accept after minor revision, with the local-entanglement language tightened.","headline":"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.","tokens_in":7929,"tokens_out":1735,"would_cite":true,"duration_ms":18676,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"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.","keywords":["Julia","quantum magnetism","spin dynamics","SU(N) coherent states","linear spin wave theory","classical spin dynamics","dynamical spin structure factor","scattering experiments"],"falsifier":"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.","tokens_in":6955,"feed_emoji":"🧲","tokens_out":6243,"duration_ms":57675,"temperature":0.7,"pith_summary":"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.","feed_headline":"One Julia package now unifies two ways to simulate quantum magnets","feed_subtitle":"Its SU(N) coherent-state formalism captures spin-orbit coupling and local entanglement, with direct scattering predictions.","key_machinery":"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.","core_discovery":"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,ω).","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the generalized spin-wave theory that extends LSWT to SU(N) coherent states, the basis of Sunny's :SUN LSWT mode.","marker":"Muniz et al. (2014)"},{"why":"Derives classical spin dynamics from SU(N) coherent states, the foundation of Sunny's generalized classical mode.","marker":"H. Zhang & Batista (2021)"},{"why":"Generalizes Langevin dynamics to SU(N) coherent states, giving Sunny's thermal-bath coupling for finite-temperature and out-of-equilibrium simulations.","marker":"Dahlbom, Miles, et al. (2022)"},{"why":"Provides the symplectic geometric integrator that Sunny uses to evolve dissipationless classical spin trajectories.","marker":"Dahlbom, Zhang, et al. (2022)"},{"why":"Supplies quantum renormalizations of biquadratic and single-ion anisotropy terms used in Sunny's dipole mode.","marker":"Dahlbom et al. (2023)"},{"why":"Introduces the kernel polynomial method that Sunny uses for efficient linear spin wave theory on large magnetic cells.","marker":"Lane et al. (2024)"},{"why":"Provides the reference method for linear spin wave theory on incommensurate spiral orderings, which Sunny extends.","marker":"Toth & Lake (2015)"},{"why":"Provides the FeI2 scattering data shown as a direct experimental comparison in the paper, demonstrating the package's predictive workflow.","marker":"Bai et al. (2021)"}],"fun_headline_variants":["Sunny.jl unifies two spin simulation approaches in Julia","One Julia package for spin dynamics and scattering predictions","Sunny.jl: SU(N) coherent states for quantum magnets","Sunny.jl merges linear spin wave and classical dynamics"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Sunny.jl unifies two spin simulation approaches in Julia","One Julia package for spin dynamics and scattering predictions","Sunny.jl: SU(N) coherent states for quantum magnets","Sunny.jl merges linear spin wave and classical dynamics"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000496,"raw_usage":{"total_tokens":2357,"prompt_tokens":793,"completion_tokens":1564,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":409,"completion_tokens_details":{"reasoning_tokens":1495}},"tokens_in":409,"tokens_out":1564,"duration_ms":12667,"temperature":1.0,"reasoning_tokens":1495,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T16:25:22.046072+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[],"review_version":1}