{"id":"6705e816-63fe-4398-acad-a3b641db638b","arxiv_id":"1908.02821","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A two-species p-wave bosonic pairing Hamiltonian is exactly solved via Richardson-Gaudin integrability, exhibiting a third-order quantum phase transition and a Bose Moore-Read critical condensate.","lead":"This paper builds an exactly solvable quantum model of two bosonic species pairing through p-wave interactions, and solves its phase diagram exactly. The authors find a third-order quantum phase transition between a gapless fragmented condensate and a gapped pair superfluid, with an exact critical state analogous to the fermionic Moore-Read state.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The thermodynamic-limit reduction of the Richardson equations to Eqs. (13) is asserted without derivation or displayed verification, and the third-order transition and phase diagram rest on this unproven step.","rationale":"The reader's weakest assumption identifies the same load-bearing concern: the thermodynamic-limit reduction of the Richardson equations to the gap/number equations (13), and the numerically asserted large-g branch. My independent reading confirms that this step is the least secure part of the central claim. The exact finite-size integrability construction is a standard hyperbolic su(1,1) Richardson-Gaudin model, and the Richardson equations and eigenstates are coherent. However, the thermodynamic limit in Section IV is presented as a transformation without derivation, and the large-g solution is only described as numerically verified. The additional observation that Eq. (13) has no smooth finite-density solution for Delta = 0, mu = 0 makes the reduction's status in the gapless phase unclear, and no prescription is given for the singular condensate contribution that would carry the density. This does not mean the claimed phase diagram is wrong, but it is currently under-supported. A finite-size scaling check of the exact Richardson energy against the thermodynamic-limit formula would directly test whether the reduction is correct. Since the reader already assigned a CONDITIONAL verdict based on this same concern, my assessment does not move the verdict.","tokens_in":10495,"tokens_out":35040,"duration_ms":380592,"concrete_test":"Compute the exact Richardson ground-state energy density for fixed density rho = 0.2 at L = 100, 200, and 400 over a grid of g spanning the superfluid interval (gc, g_infinity), and compare with the energy density from Eq. (15) evaluated with the numerical solution of Eq. (13). If the difference does not decrease as 1/L toward zero, the claimed thermodynamic-limit reduction is incorrect; if it converges, the reduction is supported and the main phase-diagram predictions are confirmed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central phase diagram (Fig. 1) and the third-order nonanalyticity of the energy density (Fig. 3) are computed from the boson gap and number equations (13)-(15). The claimed reduction of the Richardson equations (3) to these equations as N,L tend to infinity is not derived, and the reduction is not uniform: for g < gc the natural solution with Delta = 0 and mu = 0 makes the right-hand side of the number equation vanish, so the finite density of the gapless fragmented phase must be carried by a singular k = 0 (or kmin) contribution that is not present in Eq. (13). The text does not specify how this singular term emerges from the pairon distribution. On the superfluid side, the large-g branch mu approximately -gamma1 g, Delta approximately gamma2 g with 4 gamma2^2 < gamma1^2 is stated to be 'numerically verified' (paragraph after Eq. (14)), but no numerical data or analytic argument are provided. Since the Volovik line, the quasi-boson dispersion, and the third-order transition are all derived from these equations, a failure of the reduction or a different large-g branch would invalidate the headline conclusions. The finite-size exact solution appears internally consistent, but the thermodynamic-limit step is the load-bearing assumption that is least secure.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper constructs an exactly solvable two-species bosonic p-wave pairing Hamiltonian by taking a linear combination of known hyperbolic su(1,1) Richardson-Gaudin integrals of motion. The finite-size exact solution is obtained from the Richardson equations, and the ground state is tracked as a function of coupling. The authors claim that for 0≤G<Gc the system is a gapless fragmented singlet-pair condensate, that at the critical coupling Gc=2/(2L+M−1) the exact ground state is a Bose analog of the Moore-Read state, and that for G>Gc the system is a gapped pair Bose superfluid. In the thermodynamic limit, the Richardson equations are asserted to reduce to boson gap and number equations (13), from which the phase diagram, the Volovik line, and a third-order quantum phase transition are derived. The paper also claims that the construction extends to any spatial dimension.","tokens_in":10786,"tokens_out":4009,"duration_ms":46605,"significance":"If the claims hold, this is a valuable addition to the small family of exactly solvable bosonic pairing models. The finite-size solution is an honest Richardson-Gaudin solution: the Hamiltonian is built from commuting su(1,1) invariants, the eigenstates are explicitly constructed, and the energy formula (8) follows directly from the Richardson equations. The identification of the ground state at Gc as a bosonic Moore-Read-type pair condensate, and the contrast with the fermionic p+ip model, are conceptually interesting and falsifiable. The main novelty—the phase diagram with a gapless fragmented phase and a gapped pair superfluid separated by a third-order transition—rests on the thermodynamic-limit reduction in Section IV, which is asserted rather than derived. The paper also gives credit where due by comparing with known fermionic results and by providing a finite-size check in Fig. 3. However, the absence of a derivation or verification of the thermodynamic-limit reduction and of the large-coupling asymptotic branch leaves the central quantitative claims under-supported.","major_comments":[{"comment":"The reduction of the Richardson equations (3) to the boson gap and number equations (13) is stated without derivation or supporting data. The reduction is not uniform in the gapless phase: for g<gc with Δ=0 and μ=0, Eq. (13) gives v_k^2=0 for all k and hence ρ=0, whereas the phase diagram in Fig. 1 assigns a finite density ρ to this region. The macroscopic occupation of the kmin state in the fragmented condensate must contribute a singular term to the number equation that is not present in Eq. (13). The paper needs to specify how this singular contribution emerges from the Richardson pairon distribution and to justify the claimed limit.","section":"Section IV, Eq. (13)"},{"comment":"The large-coupling branch μ≈−γ1g, Δ≈γ2g with 4γ2^2<γ1^2 is described as 'numerically verified,' but no numerical data, convergence criterion, or analytic argument is shown. This branch is load-bearing because it guarantees real quasi-boson energies for all k, determines the Volovik line μ+2Δ^2=0, the g∞ divergence, and the overall shape of the phase diagram. The authors should provide either a derivation of the asymptotic behavior from Eqs. (13) or a reproducible numerical verification with residuals.","section":"Section IV, paragraph after Eq. (14)"},{"comment":"The third-order nonanalyticity of E(g) is derived from the thermodynamic-limit equations (A1)-(A9), so it inherits the gap identified in the previous two comments. In addition, the expansion (A7)-(A8) assumes a specific ordering of μ, Δ, and a=μ+2Δ^2 as δ→0+; the paper should show that this ordering is the unique physical branch selected by the original finite-size equations as L→∞. Without this, the claimed universal third derivative −2π^2/gc^6 is not established.","section":"Appendix A and Eq. (16)"},{"comment":"The claim that the model is solvable in any spatial dimension is asserted in the abstract and in Section II ('It is straightforward to extend our model to higher dimensions') but no higher-dimensional Hamiltonian, single-particle dispersion, or su(1,1) construction is provided. Since the entire analysis uses η_k=sin(k/2) on a one-dimensional chain, the dimensional generalization should either be demonstrated explicitly or the claim should be softened to a conjecture. This is part of the paper's advertised scope and needs support.","section":"Abstract and Section II"}],"minor_comments":[{"comment":"In the sentence 'This latter condition guarantees that the quasi-boson energies ... are always real, even in the limit g→∞,' the text should state whether the condition 4γ2^2<γ1^2 is verified for all ρ or only for the numerically checked cases.","section":"Section IV, after Eq. (13)"},{"comment":"The pair condensate at G∞ is written as (∑_k η_k K_k^+)^M|0⟩, but the earlier sums are restricted to k>0; the notation in Eq. (12) should be made consistent with the earlier k>0 convention.","section":"Section III, Eq. (12)"},{"comment":"The exponential factors in Eq. (16) and Eq. (A9) are typeset ambiguously; a reader cannot tell whether the last factor is e^{2(gc−1)}/g̃ or e^{2(gc−1)/g̃}. The formula should be rewritten with clear parentheses.","section":"Section IV, Eq. (16) and Appendix A"},{"comment":"The Supplemental Material reference contains a publisher URL placeholder; the authors should provide an accessible link or archive identifier.","section":"Supplemental material, Ref. [31]"},{"comment":"The caption mentions 'g = 1.1' in the inset but the list of five couplings in the main text is g = 0.5, 1.2, 1.8, 5.0, 6.8; the caption should be aligned with the displayed data.","section":"Section V, Fig. 4"}],"recommendation":"major_revision","confidential_remarks":"The paper is a compact Letter whose central quantitative claims depend on the thermodynamic-limit reduction of the Richardson equations. The finite-size exact solution is sound, and the third-order transition is plausibly correct, but the current manuscript does not provide enough support for the reduction or for the large-coupling asymptotics. This is fixable in revision, so I recommend major revision rather than rejection. I would also encourage the editor to ask for at least one explicit finite-size scaling test of Eq. (13), e.g., a plot of the numerical solution of the finite system approaching the thermodynamic-limit curves for ρ=0.2."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"What you should know: this is a genuinely new integrable model, not a repackaging. The authors take the known hyperbolic su(1,1) Richardson-Gaudin integrals and construct a two-species p-wave boson Hamiltonian (Eq. 6) that is exactly solvable by design. That construction is not circular—the phase diagram, critical coupling, and critical state are derived consequences of the Richardson equations, not inputs. The finite-size solution is standard machinery and appears internally consistent. The Bose Moore-Read state at G_c = 2/(2L+M-1) is a new object, and a third-order quantum phase transition in the thermodynamic limit is a solid result if the reduction to Eqs. (13)-(15) holds.\n\nThe soft spot is exactly there. The text says that in the thermodynamic limit the Richardson equations \"transform into\" the boson gap and number equations, with no derivation. The stress-test concern is on point: for g < g_c, with Delta = 0 and mu = 0, the right-hand side of the number equation vanishes identically, so Eq. (13) cannot produce the finite density of the fragmented BEC phase. That density must be carried by a singular k = 0 (or k_min -> 0) condensate contribution that is never written down. The paper is aware that the gapless phase is a fragmented condensate, but it never reconciles that with the integral equations. On the superfluid side, the large-g asymptotics mu ~ -gamma1 g, Delta ~ gamma2 g with 4 gamma2^2 < gamma1^2 is stated as \"numerically verified,\" but no numerical data or analytic argument is shown. These are load-bearing for the phase diagram, the Volovik line, and the third-order discontinuity.\n\nA minor issue: the claim of solvability in any spatial dimension is stated but only the one-dimensional chain is worked out. In the fermionic case that extension is known, so this is likely fine, but it is asserted rather than shown.\n\nThis paper deserves a serious referee. The finite-size exact solution and the critical state are worth publishing even if the thermodynamic limit needs tightening. The referee should ask for a derivation of the reduction to Eqs. (13)-(15) or an explicit treatment of the singular condensate contribution in the gapless phase, plus numerical support for the large-g branch.\n\nI would bring it to our reading group and cite it if I worked on Richardson-Gaudin models or ultracold bosonic mixtures. Send it to peer review with the expectation of a revision.","headline":"New exactly solvable p-wave bosonic pairing model with a third-order transition and a Bose Moore-Read critical state; the finite-size solution is solid, but the thermodynamic-limit reduction needs real work.","tokens_in":11283,"tokens_out":3603,"would_cite":true,"duration_ms":39839,"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":"An exactly solvable two-species p-wave bosonic pairing model has a third-order quantum phase transition from a fragmented singlet-pair condensate to a gapped pair Bose superfluid, with a Bose Moore-Read pair condensate at the critical…","keywords":["p-wave bosonic pairing","Richardson-Gaudin integrable models","su(1,1) algebra","Bose Moore-Read state","quantum phase transition","pair Bose superfluid","fragmented condensate"],"falsifier":"Solve the full Richardson equations (3) for increasing $M$ and $L$ at fixed density $\\rho$ and large $gL$, and compare the resulting ground-state energy density with the thermodynamic formula (15) using the $\\mu \\approx -\\gamma_1 g$, $\\Delta \\approx \\gamma_2 g$ branch. If no solution with $4\\gamma_2^2 < \\gamma_1^2$ exists for some $g$, or if the exact energy density departs from the third-order behavior of Eq. (16), the predicted pair superfluid phase and transition order would be ruled out.","tokens_in":10292,"feed_emoji":"🧊","tokens_out":9930,"duration_ms":96067,"temperature":0.7,"pith_summary":"The paper introduces a two-species bosonic model with p-wave pairing and shows it is exactly solvable in any spatial dimension through an integrable su(1,1) pairing algebra. The central result is a quantum phase diagram with two phases: for weak attraction the ground state is a gapless fragmented condensate of singlet pairs, and for strong attraction it is a gapped pair Bose superfluid. Exactly at the critical coupling the ground state is a Bose Moore-Read pair condensate, the bosonic analogue of the paired quantum Hall state. If the construction holds, the model provides one of the few exactly solvable descriptions of bosonic pairing and a concrete starting point for understanding p-wave bosonic quantum liquids.","feed_headline":"Bosonic p-wave superfluid solved exactly in any dimension","feed_subtitle":"A third-order transition separates a fragmented pair condensate from a pair Bose superfluid.","key_machinery":"The engine is the hyperbolic su(1,1) Richardson-Gaudin integrability structure built from singlet-pair operators $K^+_k = b^\\dagger_k a^\\dagger_{-k} - a^\\dagger_k b^\\dagger_{-k}$, which create pairs with opposite momenta and opposite pseudo-spin. The commuting integrals of motion in Eq. (2) combine linearly into the Hamiltonian of Eq. (7), so every eigenstate is a Richardson state (4) labelled by spectral parameters $e_\\alpha$ (the pairons) that solve Richardson equations (3). In the thermodynamic limit these equations reduce to the boson gap and number equations (13)-(14), and it is this reduction that yields the critical coupling, the Volovik line, and the third-order discontinuity of the energy density.","core_discovery":"The central claim is that the Hamiltonian of Eq. (6), with single-particle dispersion $\\sin^2(k/2)$ and two-body p-wave attraction $G$, is exactly integrable and its ground state changes character at a critical coupling. In the balanced, zero-momentum sector the weak-coupling phase, $0 \\le g < g_c$, is a gapless fragmented singlet-pair condensate whose occupation sits at the lowest finite momentum $\\pm k_{\\min}$; at $g_c = 2/(2+\\rho)$ in the thermodynamic limit ($G_c = 2/(2L+M-1)$ in finite size) all pairons collapse to zero energy and the exact ground state is the Bose Moore-Read pair condensate of Eq. (9). Above $g_c$ the ground state is a gapped pair Bose superfluid, and the ground-state energy density is non-analytic at $g_c$ with vanishing first and second derivatives and a discontinuous third derivative. The superfluid quasiboson dispersion has three regimes, with the minimum at $k=0$, at $0<k<\\pi$, or at $k=\\pi$, separated by the Volovik line and the line $\\mu + 2\\Delta^2 = 1$.","pith_inferences":["The same hyperbolic su(1,1) machinery with $Q \\neq 0$ likely yields exactly solvable finite-momentum-pair (Larkin-Ovchinnikov-type) bosonic phases; the paper notes the algebra is available but does not map that regime.","The Bose Moore-Read critical state could serve as a fixed-point ansatz for perturbative studies of p-wave bosonic droplets; the paper says the pair superfluid phase is a candidate quantum liquid but does not establish self-binding.","In the fermionic $p+ip$ case the analogous $k=0$ occupation jump is tied to a topological transition, so an open question this paper leaves implicit is whether the bosonic jump at $g_c$ also carries topological or geometric meaning.","A direct experimental test would be Bragg spectroscopy of the quasiboson minimum: crossing the Volovik line should shift $k_{\\min}$ discontinuously, a signature that does not depend on the fine details of the trap."],"forward_implications":["Below the critical coupling $g_c$, the ground state is gapless and macroscopically occupied at the lowest nonzero momentum pair states, not at zero momentum.","At $g_c$ the full set of pairons collapses to zero energy and the exact ground state is the Bose Moore-Read pair condensate, an algebraic bosonic counterpart of the fermionic Moore-Read state.","Above $g_c$ the system is a gapped pair Bose superfluid, and the energy density is non-analytic at $g_c$: the first two derivatives are continuous and the third derivative jumps.","The quasiboson dispersion in the superfluid phase has three shapes, with its minimum at $k=0$ for weak densities, at intermediate $k$ in a broad region, and at $k=\\pi$ for strong coupling; the boundary of the first region is the Volovik line.","Because the construction is integrable in any dimension and includes imbalanced mixtures and finite center-of-mass momentum pairs, the phase diagram is not restricted to the one-dimensional balanced case studied in detail."],"supporting_citations":[{"why":"Introduces the hyperbolic su(1,1) Richardson-Gaudin integrals of motion that the Hamiltonian is built from.","marker":"[2]"},{"why":"Establishes the integrability structure used to construct the commuting operators $R_k$.","marker":"[6]"},{"why":"Supplies the thermodynamic-limit boson gap and number equations that become Eqs. (13).","marker":"[8]"},{"why":"Provides the fermionic p+ip model whose Moore-Read point the Bose Moore-Read state mirrors.","marker":"[13]"},{"why":"Provides the thermodynamic-limit analysis and third-order-transition framework the bosonic derivation follows.","marker":"[14]"},{"why":"Supports the fragmented singlet-pair condensate as the weak-coupling ground state.","marker":"[32]"},{"why":"Supports the mesoscopic degeneracy of spin states and selection of the singlet pair condensate.","marker":"[33]"},{"why":"Context for the third-order quantum phase transition in hyperbolic pairing models.","marker":"[34]"}],"fun_headline_variants":["Exactly solvable p-wave bosonic superfluid","Bosonic p-wave superfluid: exact in any dimension","Third-order transition in exact bosonic p-wave model","Fragmented to pair superfluid: exact p-wave bosonic model"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the exact Richardson equations reduce to the boson gap and number equations (13) as $N$ and $L$ go to infinity at fixed density, and that for large attractive coupling those equations have the real branch $\\mu \\approx -\\gamma_1 g$, $\\Delta \\approx \\gamma_2 g$ with $4\\gamma_2^2 < \\gamma_1^2$; the paper reports numerical verification but no proof of this branch.","fun_headline_variants_meta":{"raw":{"variants":["Exactly solvable p-wave bosonic superfluid","Bosonic p-wave superfluid: exact in any dimension","Third-order transition in exact bosonic p-wave model","Fragmented to pair superfluid: exact p-wave bosonic model"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000217,"raw_usage":{"total_tokens":1402,"prompt_tokens":881,"completion_tokens":521,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":497,"completion_tokens_details":{"reasoning_tokens":454}},"tokens_in":497,"tokens_out":521,"duration_ms":5621,"temperature":1.0,"reasoning_tokens":454,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:32:52.631589+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Solve the full Richardson equations (3) for increasing $M$ and $L$ at fixed density $\\rho$ and large $gL$, and compare the resulting ground-state energy density with the thermodynamic formula (15) using the $\\mu \\approx -\\gamma_1 g$, $\\Delta \\approx \\gamma_2 g$ branch. If no solution with $4\\gamma_2^2 < \\gamma_1^2$ exists for some $g$, or if the exact energy density departs from the third-order behavior of Eq. (16), the predicted pair superfluid phase and transition order would be ruled out.","supporting_citations":[{"cited_title":"Sierra, J","cited_arxiv_id":null,"evidence_quote":"Establishes the integrability structure used to construct the commuting operators $R_k$."},{"cited_title":"Ortiz and J","cited_arxiv_id":null,"evidence_quote":"Supplies the thermodynamic-limit boson gap and number equations that become Eqs. (13)."},{"cited_title":"Bortz, S","cited_arxiv_id":null,"evidence_quote":"Provides the fermionic p+ip model whose Moore-Read point the Bose Moore-Read state mirrors."},{"cited_title":"Dukelsky, G","cited_arxiv_id":null,"evidence_quote":"Provides the thermodynamic-limit analysis and third-order-transition framework the bosonic derivation follows."},{"cited_title":"Radzihovsky and S","cited_arxiv_id":null,"evidence_quote":"Supports the fragmented singlet-pair condensate as the weak-coupling ground state."},{"cited_title":"Spontaneous formation of polar superfluid droplets in a p-wave interacting Bose gas","cited_arxiv_id":"1905.08463","evidence_quote":"Supports the mesoscopic degeneracy of spin states and selection of the singlet pair condensate."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Context for the third-order quantum phase transition in hyperbolic pairing models."}],"review_version":1}