{"id":"b3d14edc-ebe2-45b4-aa40-7defcc9bbcca","arxiv_id":"2506.10254","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"DFT+U and spin-dynamics simulations find an antiferromagnetic J1-J2 Heisenberg model for CaMn2P2 whose ground state is a spin spiral with q=(1/6,1/6,0), and predict a spin-liquid-like phase at large frustration.","lead":"Using density functional theory plus spin dynamics, the authors show that the magnetic behavior of the compound CaMn2P2 can be described by two competing antiferromagnetic interactions, and that this competition produces a spiral arrangement of magnetic moments. The paper is useful because it reproduces the experimentally observed spin spiral and maps out what happens as the competition is tuned, including a possible spin-liquid-like phase at strong frustration.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Figure 3 compares only three trial ordered states; if the SpinW/Luttinger-Tisza minimization was not global over q, the claimed q=(1/6,1/6,0) ground state is a selection among candidates rather than a prediction from the fitted J1-J2-J3 model.","rationale":"The paper has a genuine quantitative success: GGA+U with two independent codes gives a localized Mn2+ state, exchange couplings that produce a spiral at q=(1/6,1/6,0), and an ASD transition temperature (66 K) close to experiment (69.8 K). The Heisenberg truncation with J3 only 4% of J1 is plausible, and the fitted J1-J2 ratio falls in the reported spiral window. I considered the reader's flagged assumptions about omitted longer-range or anisotropic exchanges and about the classical-spin-liquid claim. The omitted longer-range exchange is a standard fitting risk but would require additional DFT data to test; the spin-liquid label is speculative but is presented in the text as a 'possible' phase in a model regime rather than as the measured ground state. The more precise, internal, and immediately testable weakness is the q-space minimization: Figure 3 and the accompanying text suggest that only three commensurate q-vectors were compared when constructing the phase diagram. If true, the match to neutron diffraction is partly built into the choice of trial states. This is not an accusation of misconduct; it is a request to see the unconstrained minimization or to run it. The proposed check is cheap and decisive: scan a dense q mesh for the same Hamiltonian. I therefore recommend no change to the reader's CONDITIONAL verdict, with the added condition that the authors report the global q-space search and, for the spin-liquid region, system-size and U-dependence checks.","tokens_in":14816,"tokens_out":8831,"duration_ms":111870,"concrete_test":"Re-run the classical ground-state search for the same Hamiltonian (J1 = -39.96 meV, J2 = -16.02 meV, J3 = -1.72 meV) and for J2/J1 in [0.2, 1.0], using an unconstrained Luttinger-Tisza or energy minimization over a dense q mesh (e.g., at least 101x101 points in the a-b plane) that also allows an arbitrary phase between the two Mn sublattices. Compare the energy at q=(1/6,1/6,0) with the global minimum. If the global minimum at J2/J1=0.4 is q=(1/6,1/6,0), the central spiral claim survives; if a neighboring incommensurate q is lower by more than numerical tolerance, the predicted ground state is not the reported one.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's strongest quantitative success is matching the neutron-diffraction spiral q=(1/6,1/6,0) using the DFT-derived Heisenberg parameters. The text states that the Luttinger-Tisza method in SpinW was used to identify wavevectors minimizing the Hamiltonian, but the phase-diagram figure (Fig. 3, Sec. III, \"Magnetic ground state using linear spin wave theory\") plots energies for exactly three commensurate q vectors: (0,0,0), (1/6,1/6,0), and (1/3,1/3,0). The text then reports q=(1/6,1/6,0) as the ground state over the entire interval 0.23 < J2/J1 < 0.52. For a classical Heisenberg model, J(q) is a sum of cosines, and its global minimum generically moves continuously with J2/J1. A plateau in q over a wide parameter range therefore suggests that only a discrete set of trial q values was compared, so earlier statements about an unconstrained minimization are not clearly supported by the data presented. If the search was restricted to these candidates, the agreement with neutron diffraction is weaker than claimed: the experimental q was one of the three trial vectors, and the calculation only shows it beats two competitors. A nearby incommensurate q or a different spiral pitch could be lower in energy. This concern is internal to the spin-wave step and does not depend on questioning the fitted exchange parameters. A separate weakness remains for the spin-liquid-like phase (J2/J1 > 0.91), which is inferred from classical sLLG dynamics with no system-size or U-sensitivity checks, but the q-space issue is more load-bearing because it affects the paper's central experimental match.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript combines DFT+U electronic structure calculations (WIEN2k and RSPt) with exchange coupling extraction via the magnetic force theorem, linear spin wave theory (SpinW), and atomistic spin dynamics (UppASD) to study the magnetism of CaMn2P2. The authors find an indirect-gap semiconductor with localized Mn2+ moments, and compute nearest-, next-nearest-, and third-nearest-neighbor exchanges J1 = -39.96 meV, J2 = -16.02 meV, and J3 = -1.72 meV. From a J1-J2-J3 Heisenberg model they identify a spin-spiral ground state with propagation vector q = (1/6, 1/6, 0), consistent with neutron diffraction, and reproduce the experimental ordering temperature (calculated TN = 66 K vs measured ~70 K). By tuning J2/J1 they construct a phase diagram with Néel order, two spiral phases, and, for J2/J1 > 0.91, a disordered slow-relaxing low-temperature phase that they interpret as a possible spin-liquid-like state. The central claims are that CaMn2P2 is a 3D realization of the J1-J2 model and that the computed exchanges quantitatively explain the magnetic properties.","tokens_in":15203,"tokens_out":5562,"duration_ms":60424,"significance":"If the claims are confirmed, the paper provides a valuable first-principles characterization of a three-dimensionally frustrated magnet, with parameter-free (up to U and Hund's coupling) extraction of exchange constants and a direct prediction of the experimental spiral wavevector and transition temperature. The use of two independent all-electron methods and the direct comparison with neutron diffraction and specific-heat data are strengths. The phase diagram for J2/J1, including the predicted Néel-to-spiral transitions and the disordered regime, is a useful guide for future chemical substitution or pressure experiments in CaAl2Si2-type compounds. However, the significance is moderated by two gaps: the global-minimization evidence for the spiral wavevector is not shown, and the spin-liquid-like phase is deduced from classical spin dynamics without finite-size or sensitivity analysis. The spin-liquid claim, as presently supported, is not yet commensurate with the strength of the wording in the title and abstract.","major_comments":[{"comment":"The claim that the Luttinger-Tisza minimization identifies q=(1/6,1/6,0) as the global ground state is not supported by the data shown. Figure 3 compares energies for only three commensurate q vectors, (0,0,0), (1/6,1/6,0), and (1/3,1/3,0). For a classical Heisenberg model with J1-J2-J3, the Fourier transform J(q) is a sum of cosines, so the global minimum should vary continuously with J2/J1; a plateau in q over the range 0.23 < J2/J1 < 0.52 is therefore suspicious and suggests that only a discrete set of trial vectors was considered. Because the experimental q was one of the trial vectors, the agreement with neutron diffraction is weaker than claimed unless the authors demonstrate that no other q in the full Brillouin zone is lower in energy. I recommend showing a dense q-space scan of the classical energy for the ab initio J2/J1, providing the analytic minimum of J(q), or otherwise presenting evidence that the SpinW relaxation was not seeded by the experimental wavevector.","section":"Section III, 'Magnetic ground state using linear spin wave theory', Fig. 3"},{"comment":"The assignment of a spin-liquid-like ground state for J2/J1 > 0.91 is not adequately supported. The evidence consists of (i) a broad, peak-less Cmag(T) at low temperature, (ii) a visually disordered spin texture in Fig. 4(h), (iii) waiting-time-dependent spin autocorrelations in Fig. 5(a), and (iv) a broad, incoherent dynamical structure factor in Fig. 5(b). All of these are expected for a classical, disordered, possibly glassy spin state on a finite lattice; they do not distinguish a spin liquid from a conventional frustrated classical paramagnet or a slowly relaxing glass. No system-size scaling, no cooling-rate dependence, no equilibration test, and no sensitivity of the phase boundary to the Hubbard U are reported. Moreover, the nominal S=5/2 moments make quantum spin-liquid behavior unlikely, so the manuscript should either supply substantially stronger evidence (e.g., scaling of the correlation length and absence of long-range order in the thermodynamic limit, comparison of classical and quantum results, or a clear statement that 'spin-liquid-like' means only 'classical disordered regime') or soften the title, abstract, and conclusion accordingly.","section":"Section III, 'Magnetic transitions using atomistic spin dynamics simulations', Figs. 4(d) and 5"},{"comment":"The assertion that the isotropic Heisenberg model truncated at third-nearest neighbors fully captures the magnetism is not quantitatively demonstrated. The text states that Dzyaloshinskii-Moriya, symmetric anisotropic exchange, and single-ion anisotropy terms are negligible, but no computed values or comparison to J1 are given, and no convergence of Jij with distance beyond J3 is shown. Since the phase diagram and the predicted q vector depend on the ratio J2/J1 and on the possible presence of longer-range couplings, this missing support is load-bearing. I ask the authors to report the values (or bounds) of the anisotropic terms and of the next few exchange couplings, or to provide a convergence test with respect to the real-space cutoff.","section":"Section III, 'Inter-site exchange couplings and Heiseberg Spin-Hamiltonian', Table I"}],"minor_comments":[{"comment":"The phrase 'a isotropic Heisenberg Hamiltonian' should be 'an isotropic Heisenberg Hamiltonian'.","section":"Abstract"},{"comment":"The notation '− →S i' is malformed; use \\vec{S}_i and define the site average over the lattice.","section":"Section II, Eq. (5)"},{"comment":"The sentence containing 'and and N is the number of lattice sites' has a duplicated 'and'.","section":"Section III, text near Fig. 5"},{"comment":"'The structure parameters used in our calculations' should read 'The structural parameters used in our calculations'.","section":"Section II"},{"comment":"The caption 'three di fferent ordering vector' should be 'three different ordering vectors'.","section":"Fig. 3 caption"},{"comment":"The manuscript does not state the system size (number of spins) and boundary conditions used in the UppASD simulations; these details are needed to judge the finite-size effects in Figs. 4 and 5.","section":"Section III, 'Magnetic transitions using atomistic spin dynamics simulations'"}],"recommendation":"major_revision","confidential_remarks":"The paper is within scope for a condensed-matter journal. The main value is the ab initio exchange parameters and the reproduction of TN. I would be comfortable with acceptance after the authors provide a full q-space minimization and either substantially strengthen or explicitly downgrade the spin-liquid claim. Please ensure the title and abstract align with the evidence level."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis is a genuine first-principles result: DFT+U with the magnetic force theorem gives J1, J2, J3 for CaMn2P2, and the Heisenberg model built from these exchanges reproduces the experimentally observed spiral q=(1/6,1/6,0) and a transition temperature of 66 K against 69.8 K experimental. That is not trivial. The minimal J1-J2 model with J2/J1=0.4 falling in the spiral window is a nice structural insight, and the phase diagram connecting the isostructural compounds is a useful contribution.\n\nThe soft spots are two.\n\nFirst, the spin-wave step. The text says Luttinger-Tisza was used to minimize over q, but Fig. 3 only shows energies for three commensurate q vectors: (0,0,0), (1/6,1/6,0), and (1/3,1/3,0). A plateau in q over 0.23<J2/J1<0.52 looks odd if the minimum is continuous. The authors need to show that the search was global over q, e.g. a plot of J(q) or energy versus q, otherwise the match to experiment is partly a selection among three candidates. This is the most load-bearing issue and it is fixable.\n\nSecond, the spin-liquid-like phase for J2/J1>0.91 is inferred from classical sLLG dynamics on a finite lattice. No system-size analysis, no U-sensitivity check, and the autocorrelation aging analysis is qualitative. The paper's own text is appropriately cautious ('suggesting', 'possible'), but the title says 'spin liquid ground state'. That is an overreach. A classical spin liquid is not a quantum spin liquid; at most this is a frustrated classical disordered phase.\n\nAlso, no code or data artifacts are provided, which makes reproduction of the sLLG results harder. That is minor but worth a referee note.\n\nThe citation pattern looks fine: they cite the experimental neutron work and prior theoretical studies of J1-J2 models. No obvious missing references.\n\nOverall: the core spiral result is credible and worth publishing after revision. The spin-liquid claim should be softened and the q-space search clarified.\n\nRecommendation: send it to peer review. It deserves referee time.\n\nBest,","headline":"Credible first-principles model of CaMn2P2's spiral order; the spin-liquid claim is overreached and the q-space search needs clarification, but the core physics holds.","tokens_in":15735,"tokens_out":2416,"would_cite":true,"duration_ms":28592,"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":"This paper argues that the magnetism of CaMn2P2 reduces to a two-coupling Heisenberg model whose competition yields the observed spin spiral and, at strong frustration, a spin-liquid-like disordered phase.","keywords":["CaMn2P2","frustrated magnetism","J1-J2 Heisenberg model","spin spiral","spin liquid","DFT+U","atomistic spin dynamics","magnetic force theorem"],"falsifier":"A DFT+U scan varying the Hubbard $U$ between roughly 2 and 6 eV would settle the ground-state prediction: if the ratio $J_2/J_1$ leaves the 0.23 to 0.52 interval under a reasonable choice of $U$, the computed $q=(1/6,1/6,0)$ spiral is not a robust consequence of the model.","tokens_in":14608,"feed_emoji":"🧲","tokens_out":16669,"duration_ms":160401,"temperature":0.7,"pith_summary":"This paper argues that the magnetism of the layered compound CaMn$_2$P$_2$ is governed by two competing antiferromagnetic exchange couplings: a strong coupling $J_1$ between adjacent Mn layers and a frustrated coupling $J_2$ within each triangular Mn layer. The authors extract both couplings from density functional theory with a Hubbard $U$ correction, then solve the resulting isotropic Heisenberg spin Hamiltonian by spin-wave minimization and by stochastic spin dynamics. The same minimal model reproduces the spin-spiral ordering vector $q=(1/6,1/6,0)$ seen in neutron diffraction and a magnetic transition temperature near 66 K, close to the measured 69.8 K. Tuning the ratio $J_2/J_1$, the model yields a phase diagram that runs from Néel order through two different spin spirals to a disordered spin-liquid-like state at large frustration. If correct, this would make CaMn$_2$P$_2$ a three-dimensional realization of the $J_1$-$J_2$ model, which until now has mostly been studied in two dimensions.","feed_headline":"A two-coupling spin model reproduces CaMn2P2's spiral magnet","feed_subtitle":"If right, CaMn2P2 is the first 3D case of the classic two-coupling spin model, including a possible spin-liquid regime.","key_machinery":"The central object is the Heisenberg spin Hamiltonian $H = -J_1 \\sum_{\\langle ij\\rangle} \\mathbf{S}_i \\cdot \\mathbf{S}_j - J_2 \\sum_{\\langle ik\\rangle} \\mathbf{S}_i \\cdot \\mathbf{S}_k - J_3 \\sum_{\\langle il\\rangle} \\mathbf{S}_i \\cdot \\mathbf{S}_l$, with exchange parameters extracted from the DFT+U electronic structure through the magnetic force theorem, a linear-response method that maps total-energy differences onto pairwise couplings. The essential competition is between the non-frustrated interlayer coupling $J_1$ and the frustrating intralayer coupling $J_2$ on the triangular Mn network. That Hamiltonian, truncated at the tiny $J_3$, is then used in two ways: a reciprocal-space spin-wave minimization selects the ordering vector of the classical ground state, and stochastic spin-dynamics simulations generate the finite-temperature phase behavior, including the ordering temperatures and the candidate spin-liquid signatures.","core_discovery":"On its own terms, the central claim is that a classical isotropic Heisenberg Hamiltonian with only two significant exchange constants describes the magnetism of CaMn$_2$P$_2$: $J_1 = -39.96$ meV couples Mn spins across the bilayer along the $c$-axis and $J_2 = -16.02$ meV couples spins within the triangular $a$-$b$ plane, while the third-neighbor coupling $J_3$ is only about 4 percent of $J_1$ and can be neglected. Because $J_1$ is a non-frustrated antiferromagnet and $J_2$ is a frustrating antiferromagnet on a triangular lattice, the energy-minimizing spin configuration is an in-plane spiral with propagation vector $q=(1/6,1/6,0)$ and with adjacent layers coupled antiferromagnetically, exactly as measured by neutron diffraction. Finite-temperature spin dynamics on the same Hamiltonian give a magnetic specific-heat peak at 66 K, matching the observed transition temperature, and the computed spin autocorrelation and dynamical structure factor at $J_2/J_1 > 0.91$ show slow relaxation, aging, and broad momentum-selective spectral weight that the authors read as signatures of a spin-liquid-like ground state. The paper also claims a general phase diagram in the $J_2/J_1$ plane: Néel order below 0.23, the $q=(1/6,1/6,0)$ spiral between 0.23 and 0.52, a $q=(1/3,1/3,0)$ spiral above 0.52, and a disordered phase above 0.91.","pith_inferences":["A decisive check neither reported nor implied by the paper is finite-size scaling: repeating the $J_2/J_1 > 0.91$ simulations on progressively larger lattices would show whether the slow relaxation and broad spectral features persist or anneal into a conventional ordered or glassy state. (This check is an editorial suggestion, not a paper claim.)","The phase diagram implies a materials-search strategy: extracting $J_1$ and $J_2$ for the related compounds CaMn$_2$As$_2$, CaMn$_2$Sb$_2$, and CaMn$_2$Bi$_2$ would place each on the same ratio axis and could identify which, if any, already sits in the spin-liquid window. The paper discusses these compounds but does not compute their couplings.","The magnetic specific-heat peak shape from the simulations at the experimental ratio could be compared quantitatively with the measured anomaly as a natural next step; the paper reports the peak temperature and shows the measured data in an inset but does not carry out that line-shape comparison."],"forward_implications":["If the model is right, tuning $J_2/J_1$ by chemical substitution, pressure, or strain should move CaMn$_2$P$_2$ and structurally related 122 compounds through the predicted Néel-to-spiral-to-disordered sequence.","The agreement between the computed 66 K transition and the measured 69.8 K supports the quantitative accuracy of the DFT+U-derived exchange couplings and makes the same extraction method applicable to other frustrated 122 pnictides.","The $q=(1/6,1/6,0)$ spiral is the stable classical ground state of the extracted Hamiltonian, giving a microscopic explanation of the neutron diffraction structure without requiring an additional Potts-nematic ordering mechanism.","A material with $J_2/J_1$ above 0.91 would be predicted to show no long-range magnetic order down to zero temperature, with slow relaxation and persistent spin fluctuations characteristic of a spin-liquid-like phase.","Because $J_1$ acts along the $c$-axis while $J_2$ acts in the $a$-$b$ plane, the model represents a genuinely three-dimensional variant of the $J_1$-$J_2$ model, so its phase diagram is a distinct prediction rather than a copy of the well-studied two-dimensional results."],"supporting_citations":[{"why":"Reports the neutron diffraction spin-spiral structure with propagation vector $q=(1/6,1/6,0)$ that the paper's spin-wave calculation reproduces.","marker":"[22]"},{"why":"Provides the experimental crystal structure, the 69.8 K transition, and the specific-heat data used for comparison.","marker":"[16]"},{"why":"Supplies the magnetic force theorem used to extract the exchange couplings $J_1$, $J_2$, and $J_3$ from the DFT+U electronic structure.","marker":"[36,37]"},{"why":"Provides the linear spin-wave theory machinery used to minimize the Heisenberg Hamiltonian and find the ground-state $q$ vector.","marker":"[42]"},{"why":"Provides the atomistic spin-dynamics simulation method used for finite-temperature dynamics and the spin-liquid-like signatures.","marker":"[43,44]"},{"why":"Gives the calculated dynamical structure factors of frustrated triangular-lattice spin models that the paper compares with its broad $S(q,\\omega)$ spectrum.","marker":"[56,57]"},{"why":"Supplies the reciprocal-space energy-minimization technique used to scan candidate ordering vectors in the spin-wave calculation.","marker":"[48]"}],"fun_headline_variants":["Two couplings explain CaMn2P2's spiral and possible spin liquid","CaMn2P2: a 3D frustrated magnet with spiral and spin-liquid states","Model with just J1 and J2 reproduces CaMn2P2's spiral magnet","Frustrated two-coupling model yields spiral and spin-liquid in CaMn2P2","CaMn2P2's magnetism captured by two exchange couplings and a spiral"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument rests on the assumption that a classical isotropic Heisenberg model with only two significant exchange couplings, and no longer-range, biquadratic, or anisotropic terms, completely describes the magnetism of CaMn$_2$P$_2$, and that spin dynamics on a finite lattice can distinguish a spin-liquid-like state from a merely disordered or glassy one.","fun_headline_variants_meta":{"raw":{"variants":["Two couplings explain CaMn2P2's spiral and possible spin liquid","CaMn2P2: a 3D frustrated magnet with spiral and spin-liquid states","Model with just J1 and J2 reproduces CaMn2P2's spiral magnet","Frustrated two-coupling model yields spiral and spin-liquid in CaMn2P2","CaMn2P2's magnetism captured by two exchange couplings and a spiral"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000704,"raw_usage":{"total_tokens":3339,"prompt_tokens":1274,"completion_tokens":2065,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":890,"completion_tokens_details":{"reasoning_tokens":1956}},"tokens_in":890,"tokens_out":2065,"duration_ms":15090,"temperature":1.0,"reasoning_tokens":1956,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T04:31:45.237357+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A DFT+U scan varying the Hubbard $U$ between roughly 2 and 6 eV would settle the ground-state prediction: if the ratio $J_2/J_1$ leaves the 0.23 to 0.52 interval under a reasonable choice of $U$, the computed $q=(1/6,1/6,0)$ spiral is not a robust consequence of the model.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports the neutron diffraction spin-spiral structure with propagation vector $q=(1/6,1/6,0)$ that the paper's spin-wave calculation reproduces."},{"cited_title":"Saparov and A","cited_arxiv_id":null,"evidence_quote":"Provides the experimental crystal structure, the 69.8 K transition, and the specific-heat data used for comparison."},{"cited_title":"Kargeti, A","cited_arxiv_id":null,"evidence_quote":"Provides the linear spin-wave theory machinery used to minimize the Heisenberg Hamiltonian and find the ground-state $q$ vector."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the reciprocal-space energy-minimization technique used to scan candidate ordering vectors in the spin-wave calculation."}],"review_version":1}