{"id":"0c3ee72a-9e5f-426a-a86b-22d9096f0949","arxiv_id":"2411.18198","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A U(1)-symmetric pseudo-Majorana functional renormalization scheme for finite-temperature XXZ magnets with Z-fields is derived and validated against exact, quantum Monte Carlo, and experimental results.","lead":"This paper extends a computational method, pseudo-Majorana functional renormalization, to spin-1/2 XXZ magnets with a magnetic field along the spin-Z axis at finite temperature. The method reproduces measured magnetization for CeMgAl11O19 and the ordering temperature for Na2BaCo(PO4)2, two materials quantum Monte Carlo cannot easily handle.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The Na2BaCo Tc claim relies on an uncontrolled one-loop truncation in a regime where the only error benchmark already shows strong deviations; no convergence test at the application parameters is given.","rationale":"Read in good faith: the paper provides a clear derivation, a public code, and several honest benchmarks, including a consistency check between susceptibility and magnetization derivative. The methodological advance, adapting pm-fRG to a U(1)-symmetric mixed representation, appears sound. The central application-level claim, however, is quantitative: reproducing measured magnetizations and transition temperatures. For CeMgAl11O19 the temperatures are high enough (T/J ~ 3) that the benchmark regime supports the calculation. The Na2BaCo claim is the load-bearing one, and it sits at the edge of the demonstrated validity range. The dimer benchmark is the only place where the controlled error is quantified, and it shows the onset of strong deviations at T ~ J/3; the application has T_c/J_z ~ 0.3. Frustration on the triangular lattice and the spontaneous UUD order parameter introduce additional uncontrolled ingredients: the frequency-box projection, the finite seed field needed to stabilize the order, and the finite correlation radius. The paper gives no direct convergence test in this regime, so the agreement with XTRG/experiment is a single-point comparison without an error estimate. The proposed check, varying nmax and (Delta_h, L) at the application parameters, would directly settle whether the one-loop truncation is the bottleneck; a stable peak position would convert the conditional acceptance into a firm one, while a shift would show that the agreement is coincidental. This agrees with the reader's weakest_assumption; no stronger internal inconsistency or non-faithfulness issue was found.","tokens_in":21512,"tokens_out":5111,"duration_ms":48209,"concrete_test":"Recompute the Na2BaCo(PO4)2 susceptibility chi~zz(T) from Eq. (37) at parameters J_z = 1.48 K, J_perp = 0.8 K, h = 2.465 K with larger frequency boxes (e.g., nmax = 16 and 24 bosonic with correspondingly larger fermionic cutoffs) and with at least one additional pair (Delta_h, L) such as (0.02 K, 12). If the peak position shifts by more than ~10% of the XTRG value (~0.45 K) or fails to trend toward the XTRG line as nmax increases and Delta_h decreases, the one-loop truncation is not converged in this regime and the agreement is not quantitatively established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The decisive application is Sec. VII B: the claim that U(1)-pm-fRG determines the transition temperature into the up-up-down spin solid of Na2BaCo(PO4)2. At the quoted parameters (J_z = 1.48 K, J_perp = 0.8 K, h = 2.465 K) the extracted peak sits near T_c/J_z ~ 0.3, which is precisely the scale where the dimer benchmark in Sec. VI A (Fig. 4) begins to deviate strongly from exact results. The square-lattice FM benchmark (Fig. 6) also shows an unphysical low-field dip, and the cubic-lattice benchmark (Fig. 7) admits a small overshoot artifact from frequency-sum truncation. No benchmark is provided for a frustrated triangular lattice in the relevant low-temperature regime. The numerical parameters are fixed at nmax = 10 bosonic Matsubara frequencies with a projection-to-boundary rule (Appendix C), and no convergence test in nmax is reported for the application. In addition, the transition is located from a susceptibility peak computed at finite seed fields Delta_h = 0.05 K and 0.0375 K and finite correlation radii L = 6, 8, 10, with no extrapolation to Delta_h -> 0 or L -> infinity; the authors state that smaller Delta_h causes numerical instabilities. Consequently, the single-point agreement with the XTRG/experimental Tc could be coincidental.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper extends the pseudo-Majorana functional renormalization group (pm-fRG) to spin-1/2 XXZ models with a magnetic field or magnetization along the spin-Z axis at finite temperature. The authors use a mixed Majorana-complex fermion ('drone-fermion') representation to diagonalize the non-interacting part, derive the U(1)-symmetry-adapted flow equations for self-energies and four-point vertices with the Katanin truncation, and express observables such as magnetization and static spin susceptibilities in terms of the resulting Green's functions and vertices. The method is benchmarked against exact dimer solutions and QMC results for the square-lattice ferromagnetic and cubic-lattice XXZ models. It is then applied to two triangular-lattice materials: CeMgAl11O19, where magnetization curves at T = 2 K and 5 K are reproduced with parameters from neutron-scattering work, and Na2BaCo(PO4)2, where a susceptibility peak is used to identify the critical temperature of the three-sublattice up-up-down transition.","tokens_in":21831,"tokens_out":5095,"duration_ms":45947,"significance":"The principal contribution is a systematic derivation of the U(1)-pm-fRG flow equations and their numerical implementation, together with reproducible code and a set of nontrivial benchmarks. If the application claims hold, the method provides a sign-problem-free tool for frustrated XXZ magnets in a field at finite temperature, which is otherwise difficult to access. The benchmarks against exact dimer results and QMC in unfrustrated systems are meaningful, and the internal consistency check in Eq. (36) and Fig. 5 is a useful addition. However, the decisive application to Na2BaCo(PO4)2 is not yet supported by a convergence test in the regime where the method is not otherwise benchmarked, and the finite-seed-field/correlation-radius extrapolation is absent; these gaps make the magnetic-field applications exploratory rather than conclusive.","major_comments":[{"comment":"The central claim that U(1)-pm-fRG determines the transition temperature of Na2BaCo(PO4)2 is not supported by a convergence test in the relevant regime. The dimer benchmark in Sec. VI A (Fig. 4) shows that the method deviates strongly from exact results for T ≲ J/3, while the susceptibility peak in Fig. 9 lies near T/J_z ≈ 0.3 for the quoted parameters J_z = 1.48 K, J_perp = 0.8 K, h = 2.465 K. The authors do not provide a benchmark for a frustrated triangular-lattice XXZ model in this low-temperature regime, nor a check of the nmax = 10 frequency box and projection-to-boundary rule (Appendix C) at these parameters. The single-point agreement with XTRG and experiment may therefore be coincidental, and the paper should demonstrate convergence or provide an independent error estimate before claiming that the method determines this transition temperature.","section":"VII B, Fig. 9"},{"comment":"The transition temperature is extracted from a susceptibility peak computed at finite symmetry-breaking fields Δh = 0.05 K and 0.0375 K and correlation radii L = 6, 8, 10. No extrapolation to Δh → 0 and L → ∞ is reported; the text states that smaller Δh causes numerical instabilities. Since the peak position is the sole basis for the transition-temperature estimate, its sensitivity to these numerical parameters must be quantified before the agreement with XTRG/experiment can be taken as a quantitative validation. A plot showing the peak position as a function of Δh and L, or a systematic extrapolation, would be necessary to support the stated accuracy.","section":"VII B, Fig. 9"},{"comment":"The unphysical low-field dip in the magnetization for the square-lattice ferromagnet at h = 0.2 is attributed to a numerical instability in vertex components. This is an admitted failure of the truncation in a regime of low field and low temperature, and the paper does not provide a diagnostic to distinguish physical from unstable flows in the material applications. Since the CeMgAl11O19 application (Sec. VII A) involves low-field and finite-temperature data, the agreement in Fig. 8 would be strengthened by showing that the flows in that parameter range are free of the instability seen in Fig. 6, or by providing a separate controlled benchmark in the relevant low-field regime.","section":"VI B, Fig. 6"}],"minor_comments":[{"comment":"The caption contains a stray symbol '□0.5' before the y-axis label; this appears to be a rendering error and should be corrected.","section":"Fig. 4 caption"},{"comment":"The derivation of the initial self-energy condition in Eq. (B7) would benefit from a more explicit statement that the constant J^z sum is the Hartree term and that it cancels the bare chemical-potential-like term from Eq. (16b); currently the cancellation is only sketched.","section":"Appendix B 2"},{"comment":"The notation 'Φj ≡ (cj, c†j, ηj)T' is introduced without a clear statement of the transpose convention for Grassmann vectors; a brief comment on the superfield ordering would reduce the chance of confusion in the subsequent derivation.","section":"Sec. II A"},{"comment":"The paper states that 'better accuracy of the fRG can be expected in higher spatial dimensions by full incorporation of the mean-field equations' (citing Ref. 16), but this expectation is not demonstrated for the triangular-lattice frustrated case that is the main application. A supporting reference or a brief numerical check would make this statement more precise.","section":"Sec. VIII"},{"comment":"The susceptibility formulas depend on vertex functions at specific bosonic transfer frequencies, but the paper does not test the effect of the finite frequency box (nmax = 10) and the projection-to-boundary rule on these particular frequency combinations. A sensitivity check for at least one of the computed susceptibilities would be useful.","section":"Eqs. (31)-(32)"}],"recommendation":"major_revision","confidential_remarks":"The stress-test concern from the reader is justified: the Na2BaCo(PO4)2 claim is a single-point agreement in a regime where the only available error benchmark already shows strong deviations, and no convergence test at the application parameters is provided. The remaining issues (finite seed fields, correlation-radius extrapolation, and the acknowledged low-field instability) are fixable within the manuscript's scope, but they are load-bearing for the headline application claims. I therefore recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is a real extension of pm-fRG to XXZ models with a field or magnetization along Z. The mixed Majorana/complex-fermion representation is not new, but the derived flow equations, the self-energy and vertex structures, and the magnetization and susceptibility formulas in this representation are. The benchmarks against exact dimer results and QMC for the square-lattice FM and cubic-lattice XXZ models are honest and mostly reassuring. The code is on GitLab, which is a plus.\n\nThe main soft spot is the Na2BaCo(PO4)2 application. The transition temperature is extracted from a susceptibility peak at T_c/J_z around 0.3, which is precisely where the dimer benchmark (Fig. 4) starts to deviate strongly from exact results. The paper gives no convergence test at those parameters: nmax=10 bosonic frequencies with projection to the boundary, and no sweep of nmax or the seed field/correlation radius. The agreement with XTRG/experiment could be coincidental. I would not call this a fatal flaw, because the method is new and the derivation is sound, but the claim in the abstract that the results \"confirm\" the model parameters is too strong. Reproducing one magnetization curve at fixed parameters is a consistency check, not a confirmation.\n\nTwo smaller points: the abstract's phrasing overstates what a single observable can establish, and the coupling-ratio statement for CeMgAl11O19 (J_perp/J_z = -0.5) looks like a typo given the quoted J_perp=0.6469 K and J_z=-0.2784 K. The experimental agreements also lack error bars or convergence estimates, so the reader cannot judge how meaningful the deviations are.\n\nThat said, the central derivation appears correct and the method will be useful to people working on frustrated XXZ magnets. I would take the paper seriously as a method paper, but the applications should be framed as demonstrations rather than confirmations. A referee should ask for a convergence test in the frustrated low-temperature regime, or at least a clear statement that the one-loop truncation is uncontrolled there. This is not a desk reject; it needs a careful referee and probably minor revisions.","headline":"A genuine and useful extension of pm-fRG to U(1)-symmetric field problems, but the headline Na2BaCo transition temperature rests on an uncontrolled truncation in exactly the regime where the method's own benchmark shows trouble.","tokens_in":22357,"tokens_out":1765,"would_cite":true,"duration_ms":17757,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A U(1)-symmetry-adapted pseudo-Majorana functional renormalization group is derived for spin-1/2 XXZ models with a field along the Z axis, and it reproduces magnetization and transition-temperature data for two frustrated…","keywords":["pseudo-Majorana fRG","XXZ model","triangular lattice","frustrated magnetism","finite temperature","magnetization","spin solid","drone-fermion representation"],"falsifier":"Run the Na2BaCo(PO4)2 calculation with a doubled Matsubara frequency box ($n_{\\mathrm{max}} = 20$) and correlation radius $L = 12$; if the peak of the sublattice susceptibility $\\tilde{\\chi}^{zz}$ shifts by more than a few percent or fails to keep sharpening with $L$, the claimed agreement with the experimental $T_c$ is not numerically converged. A second check would compare against XTRG on a cylinder of width 8 or 10 at the same field.","tokens_in":21250,"feed_emoji":"🧲","tokens_out":7592,"duration_ms":60478,"temperature":0.7,"pith_summary":"This paper extends the pseudo-Majorana functional renormalization group to spin-1/2 XXZ models at finite temperature when a magnetic field or a net magnetization points along the spin-Z direction. By switching to a mixed drone-fermion representation with one complex fermion and one Majorana per site, the Zeeman term becomes diagonal and the U(1) spin-rotation symmetry reduces to fermion number conservation, making the fRG flow tractable. The paper derives the one-loop Katanin-truncated flow equations and the expressions for magnetization and static susceptibilities, and benchmarks them against exact dimer results and quantum Monte Carlo. As demonstrations, the method reproduces the experimental magnetization curves of CeMgAl11O19 and locates the transition temperature into the three-sublattice up-up-down spin solid of Na2BaCo(PO4)2 in agreement with XTRG and experiment. A sympathetic reader would care because frustrated triangular-lattice magnets are out of reach of quantum Monte Carlo due to the sign problem, and this scheme offers a diagrammatic alternative that can handle the field-driven ordered phase.","feed_headline":"New fRG matches magnetization and Tc for two frustrated magnets","feed_subtitle":"Pseudo-Majorana fRG now handles Zeeman fields, matching data for CeMgAl11O19 and Na2BaCo(PO4)2.","key_machinery":"The drone-fermion (mixed) representation, which pairs one complex fermion $c_j$ with one Majorana fermion $\\eta_j$ per site, carries the argument: it maps the XXZ–Z Hamiltonian into a form where the Zeeman term is diagonal while the U(1) spin-rotation symmetry becomes conservation of $c$-particle number. This restricts the four-point vertices to five independent bi-local types and yields a closed set of one-loop flow equations with a Katanin-modified single-scale propagator, solved with a Lorentzian cutoff on Matsubara frequencies and a finite frequency box.","core_discovery":"The central claim is that a U(1)-symmetry-adapted pm-fRG, based on the drone-fermion representation $S^+_j = -i\\sqrt{2}\\, c_j \\eta_j$, $S^-_j = -i\\sqrt{2}\\, c^\\dagger_j \\eta_j$, $S^z_j = \\tfrac12 - c^\\dagger_j c_j$, correctly describes finite-temperature XXZ models with a Z-field, including the spontaneously magnetized side of a transition. Because the field term $h_j S^z_j$ becomes a simple on-site potential for the complex fermion, the non-interacting Green's function is diagonal in flavor without the mixing that previously blocked pm-fRG in a field. The resulting flow equations, truncated at one loop with a Katanin substitution and a Lorentzian frequency cutoff, compute magnetization from the complex-fermion self-energy and static susceptibilities from four-point vertices. Applied to CeMgAl11O19, the method reproduces measured magnetization curves at 2 K and 5 K, supporting the parameter set put forward from inelastic neutron scattering. Applied to Na2BaCo(PO4)2, the susceptibility peak associated with the up-up-down order parameter gives a critical temperature consistent with experiment and XTRG, and the result sharpens with correlation distance $L$ up to 10 on an infinite lattice.","pith_inferences":["Because the Zeeman term is diagonal for arbitrary site-dependent $h_j$, the method should handle inhomogeneous fields or local impurity fields at no extra structural cost, a natural test being a single flipped site or a staggered seed field in the ordered phase.","The Na2BaCo(PO4)2 agreement occurs at temperatures where the dimer benchmark already shows strong truncation error, suggesting that higher-dimensional frustration improves the mean-field-like cancellation; comparing against XTRG on cylinders of width larger than 6 would test this directly.","The unphysical low-temperature dip seen in the FM square-lattice benchmark at small $h$ indicates that spontaneous symmetry breaking with an infinitesimal seed field remains delicate; a temperature-flow version of U(1)-pm-fRG could mitigate this.","The same diagrammatic structure should carry over to $S = 3/2$ via faithful pseudo-Majorana representations, potentially opening finite-temperature field-dependent studies of higher-spin frustrated magnets."],"forward_implications":["Finite-temperature fRG can now be applied to frustrated XXZ magnets in a magnetic field along Z, including regions on both sides of a spontaneous U(1)-preserving ordering transition.","The magnetization data of CeMgAl11O19 are reproduced from the neutron-scattering parameter set, strengthening the case that this material sits close to the exactly solvable $J_\\perp/J_z = -0.5$ spin-liquid point.","For Na2BaCo(PO4)2, the up-up-down spin solid transition temperature at $h = 2.465\\,\\mathrm{K}$ is obtained by a system-size converged ($L \\le 10$) calculation, where the XTRG reference is limited to cylinder width 6.","The same machinery, with essentially no modification, extends to three-dimensional frustrated XXZ magnets and to other U(1)-symmetric spin models such as retarded spin-spin interactions."],"supporting_citations":[{"why":"Defines the pseudo-Majorana fRG at finite temperature and the SO(3) Majorana representation that the mixed representation is built from.","marker":"[18]"},{"why":"Prior adaptation of pm-fRG to XXZ-type anisotropy, providing the convention for XXZ couplings and benchmarks that the new scheme extends.","marker":"[24]"},{"why":"Introduced the drone-fermion representation used here to diagonalize the Zeeman term.","marker":"[26]"},{"why":"Supplies the CeMgAl11O19 model parameters, g-factor, and the experimental magnetization data being reproduced.","marker":"[29]"},{"why":"Supplies the Na2BaCo(PO4)2 model parameters, experimental transition temperature, and XTRG reference values.","marker":"[37]"},{"why":"Provides the Katanin truncation that partially includes six-point vertex effects in the flow equations.","marker":"[47]"},{"why":"Error-controlled QMC results for the ferromagnetic square-lattice Heisenberg model used as a benchmark.","marker":"[48]"},{"why":"Worm QMC code used to benchmark the cubic-lattice XXZ magnetization across the transition.","marker":"[49]"}],"fun_headline_variants":["Pseudo-Majorana fRG now handles Zeeman fields, matching two experiments","pm-fRG extended to Z-direction fields, reproduces two frustrated magnets","Zeroing in on spins: fRG matches magnetization and Tc for two materials","U(1)-adapted fRG now matches two frustrated magnet experiments","Pseudo-Majorana fRG: from magnetization to phase transition for two magnets"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The flow equations, truncated to one loop with the Katanin substitution and a finite Matsubara frequency box, remain quantitatively accurate at the low temperatures of the two material applications, although the same truncation already deviates strongly from the exact dimer solution below roughly $T \\approx J/3$.","fun_headline_variants_meta":{"raw":{"variants":["Pseudo-Majorana fRG now handles Zeeman fields, matching two experiments","pm-fRG extended to Z-direction fields, reproduces two frustrated magnets","Zeroing in on spins: fRG matches magnetization and Tc for two materials","U(1)-adapted fRG now matches two frustrated magnet experiments","Pseudo-Majorana fRG: from magnetization to phase transition for two magnets"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001123,"raw_usage":{"total_tokens":4704,"prompt_tokens":1011,"completion_tokens":3693,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":627,"completion_tokens_details":{"reasoning_tokens":3593}},"tokens_in":627,"tokens_out":3693,"duration_ms":25134,"temperature":1.0,"reasoning_tokens":3593,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T11:25:45.174052+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the Na2BaCo(PO4)2 calculation with a doubled Matsubara frequency box ($n_{\\mathrm{max}} = 20$) and correlation radius $L = 12$; if the peak of the sublattice susceptibility $\\tilde{\\chi}^{zz}$ shifts by more than a few percent or fails to keep sharpening with $L$, the claimed agreement with the experimental $T_c$ is not numerically converged. A second check would compare against XTRG on a cylinder of width 8 or 10 at the same field.","supporting_citations":[{"cited_title":"Sbierski , author M","cited_arxiv_id":null,"evidence_quote":"Prior adaptation of pm-fRG to XXZ-type anisotropy, providing the convention for XXZ couplings and benchmarks that the new scheme extends."},{"cited_title":"\\ Cottam \\ and\\ author R.B","cited_arxiv_id":null,"evidence_quote":"Introduced the drone-fermion representation used here to diagonalize the Zeeman term."},{"cited_title":"Xiang , author C","cited_arxiv_id":null,"evidence_quote":"Supplies the Na2BaCo(PO4)2 model parameters, experimental transition temperature, and XTRG reference values."},{"cited_title":"Henelius , author A.W","cited_arxiv_id":null,"evidence_quote":"Error-controlled QMC results for the ferromagnetic square-lattice Heisenberg model used as a benchmark."}],"review_version":1}