{"id":"6c88bfed-aaea-457c-9dea-b2e4e58bbf37","arxiv_id":"2411.14792","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"The authors derive closed-form resonance frequency formulas for two-sublattice ferrimagnets valid at all temperatures, including near magnetic compensation, and fit them to garnet data.","lead":"This paper derives analytic formulas for the magnetic resonance frequencies of two-sublattice ferrimagnets near the compensation temperature, for both in-plane and out-of-plane magnetization. The formulas reproduce experimental resonance data for two bismuth-doped iron garnets when the exchange stiffness and molecular field coefficient are fitted.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central claim that Eqs. (5)–(6) reproduce experiment depends on an unverified, interpolated K(T) assigned equally to both sublattices; an independent measurement of K(T) is needed to separate model prediction from parameter reconstruction.","rationale":"The theoretical derivation is standard: linearizing coupled LLG equations around a collinear equilibrium gives secular equations whose solutions reduce to Kittel and exchange-resonance forms, so I do not object to the algebra. The load-bearing step is the experimental validation: the paper's central claim is that Eqs. (5) and (6) reproduce measured frequencies, but the reproduction uses K(T) reconstructed from coercivity via a relation that diverges at T_M, with interpolation bridging the compensation region, and assumes equal sublattice anisotropy. λ and A are fit parameters, so the comparison has effectively two free parameters plus a reconstructed input. The authors themselves note the equal-K assumption is a limitation. An independent K(T) measurement would break the circularity and show whether the model has predictive content. This is the same weak spot identified by the reader, so the conditional verdict is appropriate; no change to the verdict is needed.","tokens_in":12309,"tokens_out":14831,"duration_ms":144927,"concrete_test":"Independently determine K(T) for the same GdYb-BIG and Gd-BIG crystals using torque magnetometry or angle-resolved FMR in the regime where M||H, obtaining K_Fe(T) and K_RE(T) separately. Then recompute the curves in Fig. 4 with these measured anisotropy values, allowing only λ and A to float. If the agreement with the pump-probe and BLS data degrades substantially (e.g., frequency shifts beyond the scatter of the data), the original agreement depends on the post-hoc K interpolation rather than on Eqs. (5)–(6). Reporting the raw data and error bars for the measured frequencies would make this comparison quantitative.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The weakest point is not the algebra of Eqs. (5)–(6) but the input K(T) used to compare with experiment. Section 3.1 sets K1=K2=K and obtains K from H_FMR^u = 4K/(M1-M2), where H_FMR^u is inferred from coercivity and M1(T), M2(T) from molecular-field theory. Because H_FMR^u diverges at T_M while M1-M2 vanishes, K(T) near T_M is reconstructed by interpolation, not measured. The resulting K(T) enters both the resonance formulas and the equilibrium-angle calculation that defines where the in-plane formula is applicable. Since λ and A are then fitted to the same datasets, the agreement in Fig. 4 can be partly an artifact of choosing K(T) and the fit parameters to match the target curves. The assumption K1=K2 is also acknowledged by the authors as a limitation, yet no estimate of the error it introduces is given. This does not invalidate the derivation, but it does mean the claim 'successfully reproduce experimental data' is not independently tested.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper derives closed-form expressions, Eqs. (5) and (6), for the two magnetic-resonance frequencies of a two-sublattice ferrimagnet with uniaxial anisotropy, for in-plane and out-of-plane equilibrium magnetizations, starting from linearized undamped Landau-Lifshitz-Gilbert equations. The formulas include Zeeman, exchange, demagnetizing, anisotropy, and exchange-stiffness terms and are presented as valid at all temperatures including near compensation. The authors test the formulas against time-resolved Faraday-rotation (in-plane) and Brillouin light scattering (out-of-plane) data on two Bi-substituted garnets, GdYb-BIG and Gd-BIG, using sublattice magnetizations from molecular-field theory, a reconstructed temperature-dependent anisotropy K(T), and two adjustable parameters lambda and A. They further show that far from T_M the expressions reduce to conventional FMR and exchange-resonance formulas and note the extension to gamma1 != gamma2.","tokens_in":12558,"tokens_out":8519,"duration_ms":82694,"significance":"The result, if fully supported, would be useful: it consolidates the two-mode resonance description of compensated ferrimagnets into practical analytic formulas that reduce to the known Kittel and Geschwind-Walker limits, and it demonstrates that one parameter set can describe both in-plane and out-of-plane modes when exchange stiffness is included. The paper's analytic consistency checks (convergence to conventional formulas and to the antiferromagnetic-resonance limit) are valuable, and the simultaneous reproduction of two independent measurement geometries with common lambda and A is a non-trivial success. The main reservation is that the comparison is not an independent test: K(T) is reconstructed from a divergent quantity with interpolation, and lambda and A are fit to the same experimental frequencies. Thus the significance is real but conditional on independent validation or a sensitivity analysis.","major_comments":[{"comment":"The validation is not independent. K(T) is obtained by converting the coercivity-derived anisotropy field via H_FMR^u = 4K/(M1-M2) and interpolating through the divergence at T_M; the molecular-field parameters are chosen to match the same crystals' magnetization; and lambda and A in Table I are fit using Eqs. (5) and (6) to the very frequency-vs-temperature data shown in Fig. 4. Therefore the agreement of the solid curves with the points is an in-sample fit, not a prediction, and the abstract's 'successfully reproduce' overstates the test. The discrepancy below T_M for the LF mode of Gd-BIG in Fig. 4(c) is acknowledged but not quantitatively accounted for. Please provide independent K(T) data (torque or FMR on the same crystals), or at minimum a sensitivity analysis over K and over the uncertainties in lambda and A, with residuals and error bars on the fitted parameters.","section":"Section 3.1 and Figs. 2(b), 2(d), 4"},{"comment":"The key in-plane result is not verifiable from the manuscript. The text says a quartic secular equation is obtained and then gives Eq. (5), with all details in 'Supplemental Material 1', which is not included. The approximations needed to reduce the quartic to Eq. (5), and the handling of the equilibrium angles, are therefore not stated. Please include the secular equation and the reduction in the main text or in a self-contained appendix; the same applies to the equilibrium-angle calculation in Section 4.1 (Supplemental Material 2).","section":"Section 2, Eq. (5)"},{"comment":"The assumption K1 = K2 = K is acknowledged as a limitation, but its impact on the central claim is not quantified. The anisotropy terms appear with opposite signs for the two sublattices in Eqs. (5), (6), (10), and (11), so unequal sublattice anisotropies are not a priori negligible. The authors should estimate the error, for example by repeating the fit with published Fe- and RE-sublattice anisotropy values for these garnets or by showing that the predicted frequencies near T_M are insensitive to the difference.","section":"Section 3.1 and Section 5"},{"comment":"The claim that Eqs. (5) and (6) are valid 'over all temperature ranges' should be delimited. The derivation assumes the collinear equilibrium of Eq. (4) with M1 and M2 along H0. The paper itself states that the magnetization tilts out of plane near T_M (no precession observed at 70-145 K for GdYb-BIG and 140-390 K for Gd-BIG) and that Eq. (4) is inapplicable there. The formulas may still be valid wherever the collinear state is the energy minimum, but this restriction should be stated explicitly in the abstract and conclusion.","section":"Section 4.1 and Conclusion"}],"minor_comments":[{"comment":"'bithmuth-doped' should be 'bismuth-doped'.","section":"Section 3.1, first paragraph"},{"comment":"The interpolation procedure used to obtain K(T) is not specified; please state the interpolating function (e.g., spline or polynomial) and show the raw H_FMR^u data points.","section":"Section 3.1"},{"comment":"H_FMR^u is used both for the measured anisotropy field and for the effective field 2(K1+K2)/(M1-M2); distinct symbols would prevent confusion.","section":"Sections 3.1 and 4.3"},{"comment":"'the contribution of the BSM mode to the LF mode is large enough to reach 40 GHz' is ambiguous; the LF mode frequency reaches tens of GHz because of the large k, so the sentence should be rephrased as a frequency shift or mode frequency, not a 'contribution'.","section":"Section 4.2"},{"comment":"'in the simplified case at TM' is imprecise; at the compensation temperature M1 = M2 by definition, so this is a limit taken at the compensation point, not a temperature interval. Rephrase accordingly.","section":"Section 4.3, Eq. (12)"},{"comment":"The data points are shown without error bars or a statement of frequency uncertainty; at least a representative uncertainty is needed to judge the quality of the fit.","section":"Figure 4"}],"recommendation":"major_revision","confidential_remarks":"The central algebra appears plausible and the convergence checks are a useful asset. The main issue is the in-sample validation; this is fixable with independent K(T) data or a sensitivity analysis, so I recommend major revision rather than rejection. The authors should also be required to supply the missing Supplemental Material, since the derivation of Eq. (5) is currently inaccessible to the reader."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the new thing here is the in-plane analytic resonance formula, Eq. (5), and the approximate forms Eqs. (8)-(11). The out-of-plane solution matches Kamra et al., which is good consistency. The derivation is standard two-sublattice LLG linearization, and the resulting expressions are compact enough to be practically useful across the compensation temperature. That alone is a solid incremental contribution to the ferrimagnetic spintronics toolkit.\n\nThe experimental section is the weaker part, and the weakness is exactly where the stress-test note points. K(T) is not measured independently; it is extracted from H_FMR^u = 4K/(M1-M2), which diverges at T_M, so K(T) near compensation is obtained by interpolation. The authors state the assumption K1=K2 and say the parameters might differ, but they do not quantify the error this introduces. On top of that, lambda and A are fitted using the same resonance data that the formulas are then compared with. So Fig. 4 is an in-sample parameter reconstruction more than an out-of-sample prediction. The shape of the temperature dependence still carries some weight - the formulas reproduce the dips and the HF mode rise - but a skeptic can reasonably say the agreement is not a strong test of the model.\n\nThat said, the core analytic result does not stand or fall with this experiment. Equations (5) and (6) are derived from the model, and the convergence to known FMR/exchange limits far from T_M is verified. The missing supplement prevents me from checking the algebra in detail, but the structure is standard and the consistency checks are reassuring. The paper would benefit from showing the full derivation in the main text or a publicly available supplement, and from providing error bars on the extracted frequencies and some estimate of the K(T) uncertainty.\n\nWho is this for? People modeling resonance in compensated ferrimagnets, especially experimentalists who want a closed-form expression instead of solving coupled LLG numerically. It is a genuine convenience result. I would send it to review; a referee should push on the K(T) reconstruction and the fitting protocol, but the analytic contribution deserves attention.","headline":"A useful analytic addition to ferrimagnetic resonance phenomenology, with an experimental comparison that is suggestive rather than decisive because key inputs are reconstructed or fitted.","tokens_in":13104,"tokens_out":1542,"would_cite":true,"duration_ms":15952,"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":"The paper derives closed-form resonance-frequency formulas for a two-sublattice ferrimagnet, valid for in-plane and out-of-plane magnetizations over all temperatures, and reproduces pump-probe and Brillouin scattering data on two garnets…","keywords":["ferrimagnetic resonance","compensation temperature","two-sublattice ferrimagnet","Landau-Lifshitz-Gilbert equation","exchange resonance mode","Brillouin light scattering","time-resolved magneto-optical Faraday rotation","rare-earth iron garnet"],"falsifier":"Measure $K_1(T)$ and $K_2(T)$ independently on the same garnet crystals using a technique such as torque magnetometry, then recompute Eqs. (5) and (6); if the independently measured anisotropies do not reproduce the observed low- and high-frequency modes near $T_M$, the equal-anisotropy interpolation is falsified.","tokens_in":12121,"feed_emoji":"🧲","tokens_out":17345,"duration_ms":137220,"temperature":0.7,"pith_summary":"This paper derives two analytic formulas, Eqs. (5) and (6), for the magnetic resonance frequencies of a two-sublattice ferrimagnet—a magnet with two oppositely aligned sublattice moments—and argues they remain valid at every temperature, including at the magnetization-compensation point $T_M$ where the net moment vanishes. The formulas are obtained by linearizing the coupled Landau–Lifshitz–Gilbert equations around the equilibrium configuration for both in-plane and out-of-plane magnetizations. The paper shows they reproduce temperature-dependent resonance frequencies measured by time-resolved magneto-optical Faraday rotation (in-plane) and Brillouin light scattering (out-of-plane) on two bismuth-doped rare-earth iron garnets, GdYb-BIG and Gd-BIG, using one parameter set for the molecular field and exchange stiffness. Far from $T_M$, the formulas reduce to the conventional ferromagnetic-resonance and exchange-resonance approximations, so the new expressions extend rather than replace existing models.","feed_headline":"Two equations describe ferrimagnet resonance at any temperature","feed_subtitle":"Analytic in-plane and out-of-plane formulas match data on two garnets, including near compensation.","key_machinery":"The central machinery is the two-sublattice Landau–Lifshitz–Gilbert equation driven by a free-energy density $\\Phi$ (Eq. (2)) that includes the Zeeman energy, the exchange coupling $\\lambda\\mathbf{M}_1\\cdot\\mathbf{M}_2$, uniaxial anisotropy terms, demagnetizing energy, and exchange-stiffness terms. Linearizing these coupled equations for small oscillations around equilibrium and solving the resulting secular determinant yields a quartic equation whose roots are the resonance frequencies; Eqs. (5) and (6) are the closed-form solutions. The exchange-stiffness terms are what let the model describe the backscattering magnon mode observed in Brillouin scattering, and the temperature-dependent anisotropy, reconstructed from $H_u^{\\mathrm{FMR}}=4K/(M_1-M_2)$ by interpolation, is what lets the same formulas track the modes through the compensation region.","core_discovery":"The central claim is that Eqs. (5) and (6) are closed-form analytical solutions for the two resonance branches (low-frequency and high-frequency) of a two-sublattice ferrimagnet with uniaxial anisotropy, for magnetization in the plane and perpendicular to the plane respectively. The solutions contain the molecular-field coefficient $\\lambda$, sublattice magnetizations $M_1$ and $M_2$, gyromagnetic ratios $\\gamma_1$ and $\\gamma_2$, uniaxial anisotropies $K_1$ and $K_2$, demagnetizing fields, and exchange-stiffness terms $A_i k^2$. The paper demonstrates that these formulas reproduce the measured temperature dependence of both the low- and high-frequency modes in GdYb-BIG and Gd-BIG, that they reduce to the conventional ferromagnetic-resonance and exchange-resonance equations far from $T_M$, and that the out-of-plane expression is consistent with a previously derived analytical result for a perpendicular applied field.","pith_inferences":["A natural extension would be to invert Eqs. (5) and (6) as a fitting protocol for extracting the molecular-field coefficient and exchange stiffness from resonance data on other compensated ferrimagnets; the paper only demonstrates parameter extraction for the two garnets it studies.","Because the formulas depend on wavevector $k$, they are ready-made for computing spin-wave dispersion across the compensation temperature; the paper itself reports only the uniform and backscattering modes, not a full dispersion analysis.","The equal-anisotropy assumption ($K_1=K_2=K$) could be tested by allowing distinct sublattice anisotropies; a systematic shift in the predicted low- or high-frequency mode near $T_M$ would show how much of the agreement depends on that simplification.","In the tilted equilibrium region near $T_M$ the paper switches to numerical solution; a fully analytic treatment of the noncollinear phase would be a direct follow-up."],"forward_implications":["Far from $T_M$, Eqs. (5) and (6) reduce to the conventional ferromagnetic-resonance and exchange-resonance formulas, so existing analyses remain valid and the new expressions show precisely where those approximations break down.","The formulas retain wavevector dependence through the exchange-stiffness terms $A_i k^2$, so they can describe the backscattering magnon mode in Brillouin light scattering as well as the uniform precession seen in pump-probe experiments.","The same molecular-field and exchange-stiffness parameters reproduce both the in-plane and out-of-plane data for each garnet, giving a single consistent parameter set for both geometries.","At $T_M$, with equal sublattice moments and equal gyromagnetic ratios, the out-of-plane formula reduces to the antiferromagnetic resonance expression $\\omega=\\gamma[\\pm H_0+\\sqrt{H_u(H_u+2\\lambda M)}]$, connecting compensated ferrimagnets to antiferromagnetic resonance physics.","Because the formulas do not require identical gyromagnetic ratios, they extend to ferrimagnets whose magnetization and angular-momentum compensation temperatures do not coincide, the systems used for fast domain-wall motion and skyrmion racetracks."],"supporting_citations":[{"why":"provides the general two-sublattice resonance expression that Eqs. (5) and (6) simplify.","marker":"[34]"},{"why":"derives an analytical solution for a perpendicular applied field, with which the out-of-plane formula Eq. (6) is consistent.","marker":"[35]"},{"why":"derives the exchange-resonance mode in ferrites, which the high-frequency branch reproduces far from the compensation temperature.","marker":"[6]"},{"why":"establishes exchange resonances in gadolinium iron garnet near compensation, giving the dominant term of the exchange-resonance approximation.","marker":"[7]"},{"why":"performs numerical spin-wave calculations for in-plane fields that the paper's in-plane analytic result is compared with.","marker":"[15]"},{"why":"supplies experimental resonance frequencies across the compensation point in GdFeCo that motivate the effective gyromagnetic-ratio description.","marker":"[16]"},{"why":"gives an earlier magnetization-dynamics analysis of CoGd near compensation, which the paper's all-temperature solutions extend.","marker":"[33]"},{"why":"determines the exchange stiffness constant in an iron garnet and provides the backscattering-magnon wavenumber used to fit the Brillouin data.","marker":"[13]"},{"why":"demonstrates nonthermal optical excitation of spin precession in a magneto-optical garnet, the pump-probe technique used for the in-plane measurements.","marker":"[18]"}],"fun_headline_variants":["Closed-form formulas nail ferrimagnet resonance at all temperatures","Exact resonance formulas for ferrimagnets, no temperature gaps","Two closed-form equations solve ferrimagnet resonance everywhere","Analytic fits for ferrimagnet resonance across all temperatures"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the two sublattices share one temperature-dependent uniaxial anisotropy, recovered by interpolating a measured anisotropy field that diverges at the compensation temperature; if that reconstructed $K(T)$ is wrong, the apparent agreement with the resonance data does not test the model.","fun_headline_variants_meta":{"raw":{"variants":["Closed-form formulas nail ferrimagnet resonance at all temperatures","Exact resonance formulas for ferrimagnets, no temperature gaps","Two closed-form equations solve ferrimagnet resonance everywhere","Analytic fits for ferrimagnet resonance across all temperatures"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000654,"raw_usage":{"total_tokens":2975,"prompt_tokens":905,"completion_tokens":2070,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":521,"completion_tokens_details":{"reasoning_tokens":2002}},"tokens_in":521,"tokens_out":2070,"duration_ms":53927,"temperature":1.0,"reasoning_tokens":2002,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T14:52:54.660631+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure $K_1(T)$ and $K_2(T)$ independently on the same garnet crystals using a technique such as torque magnetometry, then recompute Eqs. (5) and (6); if the independently measured anisotropies do not reproduce the observed low- and high-frequency modes near $T_M$, the equal-anisotropy interpolation is falsified.","supporting_citations":[{"cited_title":"THz-Scale Field-Induced Spin Dynamics in Ferrimagnetic Iron Garnets","cited_arxiv_id":null,"evidence_quote":"provides the general two-sublattice resonance expression that Eqs. (5) and (6) simplify."},{"cited_title":"Unconventional spin dy- namics in the noncollinear phase of a ferrimagnet","cited_arxiv_id":null,"evidence_quote":"derives an analytical solution for a perpendicular applied field, with which the out-of-plane formula Eq. (6) is consistent."},{"cited_title":"Ferrimagnetism","cited_arxiv_id":null,"evidence_quote":"derives the exchange-resonance mode in ferrites, which the high-frequency branch reproduces far from the compensation temperature."},{"cited_title":"Ultrafast spin dynamics and spintronics for ferrimagnets close to the spin compensation point (Review)","cited_arxiv_id":null,"evidence_quote":"establishes exchange resonances in gadolinium iron garnet near compensation, giving the dominant term of the exchange-resonance approximation."},{"cited_title":"Far-infrared spectra of the magnetic exchange resonances and optical phonons and their connection to magnetic and dielectric properties of Dy3Fe5O12 garnet","cited_arxiv_id":null,"evidence_quote":"performs numerical spin-wave calculations for in-plane fields that the paper's in-plane analytic result is compared with."},{"cited_title":"Evidence of relativistic field-derivative torque in nonlinear THz response of magnetization dynamics","cited_arxiv_id":"2408.05510","evidence_quote":"supplies experimental resonance frequencies across the compensation point in GdFeCo that motivate the effective gyromagnetic-ratio description."},{"cited_title":"H-T phase diagram of rare-earth–transition-metal alloys in the vicinity of the compensation point","cited_arxiv_id":null,"evidence_quote":"gives an earlier magnetization-dynamics analysis of CoGd near compensation, which the paper's all-temperature solutions extend."},{"cited_title":"Far Infrared Spectra of Rare-Earth Iron Garnets","cited_arxiv_id":null,"evidence_quote":"determines the exchange stiffness constant in an iron garnet and provides the backscattering-magnon wavenumber used to fit the Brillouin data."},{"cited_title":"Optical determination of the exchange stiffness constant in an iron garnet","cited_arxiv_id":null,"evidence_quote":"demonstrates nonthermal optical excitation of spin precession in a magneto-optical garnet, the pump-probe technique used for the in-plane measurements."}],"review_version":1}