{"id":"3961b7ae-7a56-4833-9eab-1f782faf7f6e","arxiv_id":"2412.09670","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A quaternion reformulation of composite fermion theory allows accurate calculations for hundreds of electrons, and finds no magneto-roton instability in the Jain sequence up to nu = 15/31.","lead":"Physicists studying the fractional quantum Hall effect use a powerful method called composite fermions, but its calculations become slow and unstable for large systems. This paper rewrites the method using quaternions, making it fast enough to simulate hundreds of electrons, and finds that the fractional states remain stable as they approach the half-filled Landau level.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The no-instability claim for n up to 15 rests on the untested accuracy of the JK projected CF wave functions at large n; a small-system exact-diagonalization benchmark would settle it.","rationale":"The reader identifies the same core concern: the large-n physics conclusion depends on CF trial wave function accuracy that has not been independently benchmarked for n > 7. This is the right weakest link because the method's mathematical derivation appears self-consistent, the complexity gain is real, and the error plots in the SM validate the new implementation against the old one but not the underlying JK ansatz against exact states. The no-instability result is a variational upper-bound calculation; without a benchmark, the conclusion of 'extraordinary robustness' could be an artifact of the trial states. The proposed ED test is practical because the Hilbert space dimension for large n and moderate N is small (near half filling, the number of empty orbitals is N/n), so N=30 at nu=15/31 is feasible and directly probes the untested regime. The reader's conditional verdict is therefore appropriate, and no change is needed.","tokens_in":29754,"tokens_out":10879,"duration_ms":103296,"concrete_test":"Perform exact diagonalization for N=30 electrons at nu=15/31 (2Q=32, Hilbert space dimension C(33,3)=5456) and compute the squared overlap and energy difference between the quaternion JK ground state and the exact Coulomb ground state, plus the overlap between the quaternion JK state and the exact LLL projection of the unprojected CF wave function. Compare with the known n<=7 benchmarks. If the overlap drops substantially (e.g., below ~0.9) or the energy difference grows significantly relative to n=7, the large-n robustness conclusion is not supported. The same test can be repeated for n=10,20 at the largest N with dimension <=10^5.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The single most load-bearing assumption is that the Jain-Kamilla (JK) projected composite-fermion wave functions remain quantitatively accurate representations of the exact Coulomb ground and neutral excited states for the newly accessed regime n > 7, up to n = 15 and N = 390. The paper's own error analysis (SM Sec. II, Fig. 3) demonstrates only that the quaternion implementation matches the JK ansatz to high precision; it does not test whether the JK ansatz itself approximates the exact LLL-projected state or the exact ground state at large n. Because the magneto-roton dispersion is computed as the expectation value of the Coulomb interaction in a single CF-exciton trial state, a positive Delta E is a variational upper bound on the true neutral gap and does not exclude a different soft mode or a negative true gap. The paper flags this accuracy caveat for the second Landau level but not for the lowest Landau level at large n. The finite-size extrapolation in Fig. 3a (Delta E at k -> 0 for N = n^2, n = 8..18) is also a delicate sign test, and magic-number shell effects at N = n^2 could bias the slope, but the primary unresolved issue remains the unbenchmarked accuracy of the trial wave functions.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript introduces a quaternion-based reformulation of the Jain-Kamilla (JK) lowest-Landau-level projection for composite-fermion wave functions on the sphere, claims access to Jain states with n in the range -20 <= n <= 20 and with more than 400 electrons, and uses the method to compute magneto-roton dispersions for filling factors up to nu = 15/31 as well as the k -> 0 CF-exciton energy for the CF Fermi sea at nu = 1/2. The authors report no magneto-roton instability along the sequence n/(2n+1), interpret this as evidence for extraordinary robustness of the FQHE approaching half filling, and contrast it with an instability at 2/5 and 3/7 in the second Landau level, where they explicitly note that the JK ansatz is less accurate.","tokens_in":29875,"tokens_out":4076,"duration_ms":41061,"significance":"The numerical method is a potentially important technical advance: the Supplementary Material contains a detailed derivation of the quaternion JK projection, shows that the cost scales linearly in l_max - Q*, and validates the implementation against higher-precision arithmetic to about 10^-14 accuracy. If the physics conclusion is correct, the paper provides a strong statement about the stability of the Jain sequence and the HLR Fermi sea in the lowest Landau level. However, the central physics claim rests on the untested accuracy of the JK trial wave functions in the newly accessed large-n regime, and the k -> 0 extrapolation is a delicate sign test. These issues are load-bearing and need to be addressed before the conclusion can be accepted as stated.","major_comments":[{"comment":"The accuracy validation in SM Fig. 3 benchmarks the new quaternion implementation against the same JK ansatz evaluated in higher precision; it does not establish that the JK ansatz itself approximates the exact Coulomb ground and neutral excited states for n > 7. For n up to about 7 this is supported by previous exact-diagonalization comparisons, but for nu = 8/17 through 15/31 with N up to 390 no independent benchmark is provided. Because Delta E is computed as the expectation value of the Coulomb interaction in a single CF-exciton trial state, a positive Delta E is a variational upper bound on the true neutral gap and cannot exclude a different soft mode or a negative true gap. The paper flags this accuracy caveat for the second Landau level but not for the lowest Landau level at large n; a small-system exact-diagonalization comparison of the lowest magneto-roton branch at, for example, nu = 8/17 or 9/19 would directly test the load-bearing assumption.","section":"SM Sec. II; main text Fig. 2"},{"comment":"The conclusion that the CF Fermi sea at nu = 1/2 has no nematic instability is based on the N -> infinity extrapolation of Delta E(k -> 0) for N = n^2 with n = 8..18. The plotted energies decrease slowly, and the published figure does not show a fitted functional form, a confidence interval for the thermodynamic limit, or alternative extrapolation forms. The shell-closed choice N = n^2 may also introduce systematic shell effects. The statement that the energy 'approaches zero (rather than a negative value)' needs quantitative support, such as a linear or quadratic extrapolation with an uncertainty estimate; otherwise the claim should be softened to 'no evidence for instability in the studied systems.'","section":"Main text Fig. 3a; SM Sec. III"}],"minor_comments":[{"comment":"The displayed expression for N_l^{m',Q*,Q1} contains typesetting artifacts ('radicaltp/radicalvertex'), which should be corrected to a properly typeset square root.","section":"Eq. (11), main text"},{"comment":"The abstract contains the typo 'Suplementary', and SM Sec. II contains the repeated phrase 'called Λ levels called Λ levels'; both should be fixed.","section":"Abstract and SM Sec. II"},{"comment":"The claim that the method handles 'at least -20 <= n <= 20' is supported by the SM timing and accuracy figures, but the main text states this without referencing those figures; adding a cross-reference would help readers locate the evidence.","section":"SM Sec. II, after Eq. (12)"},{"comment":"The Monte Carlo sample size and error estimation procedure are described in the SM but not in the main text; a brief sentence stating the number of samples and the jackknife procedure would improve the reproducibility of Fig. 2.","section":"SM Sec. III and main text Fig. 2"}],"recommendation":"major_revision","confidential_remarks":"The manuscript describes a genuine technical advance that will interest the FQH community. The main obstacle to acceptance is not the implementation, which appears carefully validated internally, but the absence of an independent physical benchmark for the newly accessed large-n regime. If the authors add a small-system exact-diagonalization comparison for one or two previously inaccessible fillings, or clearly reframe the robustness conclusion as conditional on the JK ansatz, I would support publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the quaternion reformulation of the JK projection is a genuine methodological advance, and the paper is honest about what it does and does not show. The physics conclusion (no instability up to nu=15/31) is a prediction of the CF ansatz, not an established fact, because there is no independent benchmark at large n. That said, the method itself is solid and deserves a serious referee.\n\nThe genuinely new piece is the reformulation of LLL projection in terms of elementary symmetric polynomials, with linear scaling in the number of filled Lambda levels and numerical stability in double precision. The supplement's derivation is detailed and self-consistent; the error plots show the new method matches higher-precision arithmetic to machine precision where the traditional derivative-based method fails for n>=8. The speed-up ratios are dramatic. The second-LL control, where the method finds an instability at 2/5 and 3/7, is a good sanity check: the method can detect a magneto-roton instability when one is present. No fitting, no circularity - the exciton energy is a variational number and is allowed to go negative.\n\nThe soft spot is exactly what your reader flagged: the accuracy of the JK wave functions at large n is not benchmarked against anything independent. The paper cites earlier ED tests for small n, but for the new regime (n up to 15, N up to 390) the only validation is agreement with the same ansatz at higher precision. That validates the implementation, not the physics. Combined with the variational nature of Delta E, the \"extraordinary robustness\" claim is really a statement about the CF exciton ansatz, not a theorem about the Coulomb ground state. The paper does flag this caveat for the second LL but not for the LLL at large n, which is an asymmetry the authors should address.\n\nTwo smaller points. The advertised reach (-20<=n<=20, more than 400 electrons) is not demonstrated by the shown data - the largest state shown is 15/31 with N=390. And the k->0 extrapolation for N=n^2 (Fig. 3a) is a delicate sign test, so the \"approaches zero from positive side\" conclusion needs more care, though the data do look persuasive. These are minor; the primary caveat is the unbenchmarked ansatz.\n\nWho should read this: anyone working on composite fermion numerics or the stability of Jain states. The method will likely become the standard tool for this subfield. I would send it to review; I would ask referees to press on the LLL accuracy question - either a small ED benchmark at n=8 or 9, or a clear statement that the robustness conclusion is a consequence of the CF ansatz and not established for the true Coulomb problem. I would not cite the robustness claim in my own work without an independent check, but I would cite the method.","headline":"A real computational advance wrapped in a physics claim that overreaches its evidence; the method deserves review, the robustness conclusion needs qualification.","tokens_in":30499,"tokens_out":2678,"would_cite":true,"duration_ms":26291,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["73.43.-f","71.10.Pm"],"model":"deepseek-v4-flash","headline":"A quaternion reformulation of composite-fermion wave functions extends quantitative fractional quantum Hall calculations to hundreds of electrons and reveals no instability along the Jain sequence up to ν=15/31.","keywords":["fractional quantum Hall effect","composite fermions","quaternions","monopole harmonics","lowest Landau level projection","Jain sequence","magneto-roton","nematic instability"],"falsifier":"Compute the low-energy neutral excitations of a large-$n$ Jain state such as $\\nu = 8/17$ or $\\nu = 10/21$ with an unbiased method (e.g., exact diagonalization on small systems or DMRG on a cylinder) and compare the magneto-roton dispersion with the quaternion CF result; a negative-energy roton, or a rapid decay of the CF wave-function overlap with system size, would overturn the robustness conclusion.","tokens_in":29444,"feed_emoji":"🧲","tokens_out":5575,"duration_ms":45461,"temperature":0.7,"pith_summary":"This paper introduces a quaternion-based formulation of composite-fermion (CF) wave functions for the fractional quantum Hall effect (FQHE). The new formulation makes the Jain-Kamilla lowest-Landau-level projection dramatically faster and numerically stable, extending reliable calculations from Jain states with index $n$ between $-5$ and $7$ to at least $n$ between $-20$ and $20$ for systems with more than 400 electrons. As a first application, the authors compute magneto-roton (neutral CF exciton) dispersions along the Jain sequence $\\nu = n/(2n+1)$ up to $\\nu = 15/31$ and find no sign of an instability, concluding that the FQHE along this sequence is highly robust and that the composite-fermion Fermi sea at $\\nu = 1/2$ shows no nematic or charge-density-wave instability in these calculations. If correct, this method opens many previously inaccessible questions in FQHE to quantitative study.","feed_headline":"Quaternions unlock fractional quantum Hall states at 400 electrons","feed_subtitle":"New composite-fermion calculations reach ν=15/31 and find the Jain sequence stable toward ν=1/2.","key_machinery":"The central object is the quaternion extension of the monopole harmonics, $Y_{Q,l,m}(r)$, built from the symmetric and antisymmetric projections of a unit quaternion $r$; these transform ordinarily under quaternion multiplication with no extra phase. The load-bearing identity is that the Jastrow factor $d(r_i,r_j) = (r_i^{-1} \\cdot r_j)_A$ is invariant under left multiplication, so the LLL projection can be done after rotating electron $i$ to $r_i = 1$. There the projection of $Y_{Q^*,l,m}$ times the Jastrow product becomes a sum of Wigner-$D$ matrices times (regularized) elementary symmetric polynomials in the variables $(r_i^{-1} \\cdot r_j)_S / (r_i^{-1} \\cdot r_j)_A$, and only $l_{\\max} - Q^* + 1$ such polynomials are needed per column. This replaces the quadratic number of numerically unstable mixed derivatives in the standard Jain-Kamilla projection with a linear, stable computation.","core_discovery":"The central discovery is that rewriting the Jain CF wave functions in terms of quaternions (unit quaternions parametrized as $r = e^{\\phi k/2} e^{\\theta j/2} e^{\\psi k/2}$) removes the phase-factor obstruction that makes monopole harmonics awkward under rotations, and reveals that the Jastrow factor $(u_i v_j - u_j v_i)$ is a quaternion displacement invariant under left quaternion multiplication. This invariance lets one rotate each electron to the identity quaternion before performing the lowest-Landau-level projection, where the projection reduces to computing elementary symmetric polynomials. The required number of such polynomials scales linearly with the CF Landau-level index $l_{\\max} - Q^*$ rather than quadratically, and they can be evaluated accurately in double precision. The authors thereby reach Jain states $\\nu = n/(2n+1)$ with $n$ up to 15 (and state the method works for about $-20$ to $20$) and hundreds of electrons, e.g., $N = 390$ at $\\nu = 15/31$ whose Fock dimension is $\\sim 10^{236}$. Using this tool, they compute magneto-roton dispersions and find that the excitation energy never goes negative, indicating no instability of the FQHE along $n/(2n+1)$ and no instability of the HLR Fermi sea at $\\nu = 1/2$, while the same calculation in the second Landau level does show instabilities at $\\nu = 2/5$ and $3/7$, consistent with the Read-Rezayi physics there.","pith_inferences":["If the quaternion formulation is as stable as claimed, it may be extended to other projected wave functions beyond Jain states, such as parton states or paired CF states, where LLL projection is also a bottleneck.","The linear scaling in $l_{\\max} - Q^*$ suggests the method could be pushed to even higher $n$ or to systems with larger $Q^*$ (higher Landau levels), possibly testing the stability of FQHE arbitrarily close to $\\nu = 1/2$.","The paper's identification of many approximately equally spaced magneto-roton minima ($n$ minima for $\\nu = n/(2n+1)$) could be connected to the higher-spin description of magneto-rotons; the method could be used to extract the chirality or spin content of these modes.","A testable prediction implied by the robustness result: in ultra-clean samples, Jain states with $n$ as high as $15/31$ (and beyond) should be experimentally observable in the lowest Landau level, and their neutral mode dispersions should show the predicted multi-minimum structure."],"forward_implications":["Reliable CF calculations for Jain states with hundreds of electrons, enabling quantitative study of thermodynamic limits of energy gaps, dispersions, and phase boundaries.","Ability to probe the approach to $\\nu = 1/2$: the magneto-roton gap at $k \\to 0$ approaches zero but does not become negative, consistent with a stable CF Fermi sea.","New access to reverse-flux Jain states ($\\nu = n/(2n-1)$, negative $n$) and larger system sizes for previously accessible fractions, improving benchmarks against experiments.","The second Landau level results show the method can detect genuine instabilities (at $\\nu = 2/5$ and $3/7$), demonstrating that the absence of instability in the LLL is a nontrivial finding."],"supporting_citations":[{"why":"Supplies the quaternion representation of spin-weighted spherical harmonics that inspired the reformulation of the monopole harmonics.","marker":"[72]"},{"why":"Introduces the Jain-Kamilla projection ansatz that this paper reformulates in quaternion language.","marker":"[7]"},{"why":"Establishes the JK method's quantitative accuracy and demonstrates the original limits on system size and filling.","marker":"[8]"},{"why":"Provides the spherical geometry and monopole harmonic framework on which the whole calculation is built.","marker":"[76]"},{"why":"Defines the monopole harmonics used in the wave functions and their transformation properties.","marker":"[73]"},{"why":"Introduces the composite-fermion theory of the FQHE that the paper extends.","marker":"[2]"},{"why":"Establishes the CF exciton model of magneto-roton modes used for the instability search.","marker":"[16]"},{"why":"Identifies the $k \\to 0$ magneto-roton minimum relevant to the nematic instability question.","marker":"[91]"}],"fun_headline_variants":["Quaternions propel fractional quantum Hall theory to new heights","Quaternion math cracks fractional quantum Hall states at 390 electrons","Rewriting CF theory with quaternions reaches ν=15/31","Quaternions unlock massive FQH simulations, test nematic instabilities","Quaternion twist stabilizes composite-fermion Fermi sea at half filling"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire calculation treats the Jain-Kamilla projected composite-fermion wave functions as quantitatively accurate representations of the true Coulomb ground and excited states, a benchmarked assumption only for fillings up to about $n = 7$; for the new large-$n$ systems there is no independent check.","fun_headline_variants_meta":{"raw":{"variants":["Quaternions propel fractional quantum Hall theory to new heights","Quaternion math cracks fractional quantum Hall states at 390 electrons","Rewriting CF theory with quaternions reaches ν=15/31","Quaternions unlock massive FQH simulations, test nematic instabilities","Quaternion twist stabilizes composite-fermion Fermi sea at half filling"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000676,"raw_usage":{"total_tokens":3110,"prompt_tokens":1012,"completion_tokens":2098,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":628,"completion_tokens_details":{"reasoning_tokens":2006}},"tokens_in":628,"tokens_out":2098,"duration_ms":13862,"temperature":1.0,"reasoning_tokens":2006,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T16:54:36.677209+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the low-energy neutral excitations of a large-$n$ Jain state such as $\\nu = 8/17$ or $\\nu = 10/21$ with an unbiased method (e.g., exact diagonalization on small systems or DMRG on a cylinder) and compare the magneto-roton dispersion with the quaternion CF result; a negative-energy roton, or a rapid decay of the CF wave-function overlap with system size, would overturn the robustness conclusion.","supporting_citations":[{"cited_title":"How should spin-weighted spheri- cal functions be defined? Journal of Mathematical Physics, 57(9), September 2016","cited_arxiv_id":null,"evidence_quote":"Supplies the quaternion representation of spin-weighted spherical harmonics that inspired the reformulation of the monopole harmonics."},{"cited_title":"Scarola, Kwon Park, and Jainendra K","cited_arxiv_id":null,"evidence_quote":"Establishes the CF exciton model of magneto-roton modes used for the instability search."},{"cited_title":"Balram, G","cited_arxiv_id":null,"evidence_quote":"Identifies the $k \\to 0$ magneto-roton minimum relevant to the nematic instability question."}],"review_version":1}