{"id":"7b9f5340-506b-4eb0-9b73-39f866d2942e","arxiv_id":"2602.22315","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A graph-based generalization of the Jastrow ansatz yields parent Hamiltonians with two-body edge interactions and three-body 2-path interactions, specifying exact ground states and energies for many new many-body models.","lead":"The paper shows that for any graph-shaped network of interacting particles, one can build a Hamiltonian whose exact ground state is a product of pair functions over the graph's edges. This gives a unified way to generate many exactly solvable many-body models, recovering known examples and adding new ones.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Domain/self-adjointness gap is load-bearing: the 'any pair function' claim fails for singular f (e.g., |x|) without an explicit boundary condition; the formal V2+V3 miss coincidence-manifold terms.","rationale":"The reader's weakest assumption—domain/self-adjointness—is precisely the load-bearing issue. The algebraic factorization H0=ΣQ_i†Q_i is a valid formal construction for smooth, positive pair functions, but the paper states it for arbitrary f. Singular f such as |x|^g or exp(g|x|) introduce distributional/boundary terms at particle coincidences that the explicit V2+V3 formulas miss unless the Hamiltonian is defined through a quadratic-form closure with an implicit boundary condition. This is not merely a technicality: for f=|x| in two-particle one-dimensional free space, the formal Hamiltonian is the free Laplacian, which does not annihilate |x|; the correct form-closure gives a Dirichlet Laplacian. Thus the central universal claim is false as stated, though it can be repaired by restricting to C² f or by specifying the self-adjoint extension/form domain for each singular case. This supports the reader's CONDITIONAL verdict rather than requiring rejection: the construction is valuable and largely correct for its main examples, but the missing domain analysis is necessary for the advertised generality.","tokens_in":17479,"tokens_out":19317,"duration_ms":178422,"concrete_test":"Fix N=2, D=1, f(x)=|x| on a finite interval [-L,L]. Compute the operator associated with the closed quadratic form h(ψ)=||ψ'−ψ/x||² on C_c∞((-L,0)∪(0,L)); identify it as −d²/dx² with Dirichlet boundary at x=0. Then verify that |x| lies in the form domain and h(|x|,|x|)=0, while the formal Hamiltonian −d²/dx² does not annihilate |x| in the distributional sense. If the two operators differ, the paper must specify the self-adjoint extension for each singular f before the ground-state claim can be accepted.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Sec. III's proof that H0Φ0=0 and H0≥0 (Eqs. 8–11) treats Φ0 and the Qi as smooth objects. For pair functions that are not C² at r=0, the Laplacian of Φ0 acquires singular/boundary contributions on the coincidence manifolds x_i=x_j that are not contained in the formal V2+V3 of Eqs. (6)–(7). Example: N=2, D=1, f=|x| (g=1). The formulas give V2=V3=0, so H0=−ℏ²/(2m)d²/dx² on R. But Φ0=|x| is not in the domain of the free Laplacian and H0Φ0 = −(ℏ²/m)δ(x) (up to constants), not zero. If one instead closes the quadratic form Σ||Q_iψ||², Q=∂−1/x, the associated self-adjoint operator is the Dirichlet Laplacian on R\\{0}, not the free Laplacian; the factorization implicitly imposes boundary conditions at coincidences that are never stated. The same mechanism affects |x|^g for 0<g≤1 and exp(g|x|), and its omission is load-bearing because the claim 'for any pair function f' and the ground-state status of Φ0 fail for the explicit Hamiltonian unless such boundary conditions are specified.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper constructs parent Hamiltonians for many-body systems whose ground state is a product of pair-correlation functions over the edge set of a graph (a ``Graph-Jastrow'' wavefunction). The parent Hamiltonian is shown to contain two-body interactions supported on graph edges and three-body interactions supported on 2-paths. The authors derive the general formula in arbitrary dimension, specialize to one dimension, and catalogue examples for complete, path, cycle, regular, star, and product graphs, recovering Calogero-Sutherland, Lieb-Liniger, Jain-Khare, and truncated models. They also show how graph joins and products generate composite solvable models, with explicit tables of Hamiltonians, wavefunctions, and energies.","tokens_in":17838,"tokens_out":16495,"duration_ms":153676,"significance":"If the constructions are made rigorous, the paper provides a useful unifying framework and taxonomy for exactly solvable many-body Hamiltonians with Jastrow-type ground states. The factorization H0 = ∑ Qi† Qi is elegant and gives a clear sufficient condition for non-negativity; the graph-theoretic organization is transparent and yields many explicit new examples. A particular strength is the large set of closed-form results in Tables I–VI, which makes the formal claims easy to test. The recovery of known models (Calogero-Sutherland, Lieb-Liniger, Jain-Khare, truncated Calogero-Sutherland) is a valuable sanity check. However, the central derivation and the stated ground-state claims have gaps for singular pair functions, and one step in the D-dimensional proof is incorrect as written. These issues are fixable but require careful revision.","major_comments":[{"comment":"The formal factorization H0 = ∑ Qi† Qi does not by itself establish that Φ0 is a ground state of the operator defined by Eqs. (5)-(7) when f is singular. Example: D=1, N=2, f=|x_ij|. Equations (6)-(7) give V2=V3=0, so the formal Hamiltonian is the free Laplacian; but Φ0=|x1−x2| has ΔΦ0 = 2δ(x1−x2) and is not in its domain. The factorization actually defines the Dirichlet Laplacian on R\\{0} for the relative coordinate, i.e., it imposes an unstated boundary condition at the coincidence manifold. The same issue affects |x_ij|^g (0<g≤1) and exp(g|x_ij|) in Tables I–V. The paper must either restrict f so that the formal potentials and the wavefunction share a domain, or define H0 through the closure of the form (11) and state the resulting boundary conditions explicitly.","section":"Sec. III, Eqs. (5)-(10)"},{"comment":"The statement ``Noting that Δ_{S^{D-1}}^i Φ0 = 0'' is false for the Graph-Jastrow wavefunction (2). Φ0 depends on r_i through the relative vectors r_i−r_j and therefore has nontrivial angular dependence in hyperspherical coordinates centered at the origin. The final D-dimensional formulas (5)-(7) are nevertheless correct; they follow from direct differentiation of the edge product, and the proof should be rewritten without the angular-Laplacian assertion.","section":"Sec. III, Eq. (4)"},{"comment":"For f = exp(g|x_ij|), the constants in V2 and V3 are positive; together with the kinetic term they cancel, so the parent Hamiltonian H0 has zero ground-state energy by construction. The quoted E0 = −ℏ²g²/(6m) N(N²−1) is the energy of the attractive Lieb-Liniger Hamiltonian, which differs from H0 by a constant shift. Please clarify which Hamiltonian this energy belongs to, or correct the sign/energy shift. This affects the exact-energy claim in the complete-graph row.","section":"Sec. VI, after Eq. (24) and Table II"}],"minor_comments":[{"comment":"The header for the cycle graph writes Φ0 = ∏_{i=1}^{N−1} f(x_i,i+1), but a cycle with periodic boundary conditions should include the edge (N,1); the product should run to N with indices modulo N.","section":"Table IV"},{"comment":"The statement ``for r=(N−1)/2 one recovers K_N'' requires N odd. For even N the complete graph is (N−1)-regular and is not of the form 2r with integer r.","section":"Sec. VIII"},{"comment":"The title and introduction use ``integrable'' but the paper only establishes ground-state solvability; no complete set of conserved quantities is constructed for the general graph-Jastrow Hamiltonians. Suggest softening the terminology to ``ground-state solvable''.","section":"Title/Introduction"},{"comment":"Typos and small presentation issues: ``Arbitary'' in Table III, ``Anstze'' in Sec. V, inconsistent sub/superscript notation for the Laplace-Beltrami operator in Eq. (4), and a few awkward hyphenations.","section":"Various"}],"recommendation":"major_revision","confidential_remarks":"The paper is likely to be of interest to the quantum many-body and mathematical physics communities. The central framework is valuable and the algebraic core is sound for smooth pair functions, but the missing domain/self-adjointness discussion for singular f, the incorrect angular-Laplacian step in the D-dimensional proof, and the energy-shift inconsistency must be addressed before publication. These are correctable within the manuscript's scope, so I do not recommend rejection. The authors should also be asked to clarify the distinction between ground-state solvability and integrability in the title and abstract."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper gives a genuinely useful way to generate exact ground states for many-body systems from graph adjacency matrices, and the two-path three-body structure is a nice observation. But the central claim is overbroad: for singular pair functions like |x|, the stated Hamiltonian is not the operator that has Φ0 as ground state. The examples in the tables include this case, so the error is load-bearing. The title also overreaches: 'integrable' is used loosely, since no conserved quantities or Lax pair are provided for the new graph models.\n\nThe derivation of the parent Hamiltonian is clean: acting with the Laplacian on the edge product yields the two-body and three-body terms, and the Q_i^†Q_i factorization establishes positivity. This is correct when f is smooth and the operator domain is not an issue. The graph-product constructions (star, wheel, ladder, Creutz) are a nice addition, and the tables organize a large family of models. The citation pattern is fair: the complete-graph results from the authors' own prior work are appropriately referenced, and the external literature on Jain-Khare and truncated CSM is covered.\n\nThe soft spots, in order of importance. First, the domain problem is not a technicality. For N=2, D=1, f=|x|, Eqs. (6)-(7) give V2=V3=0, so the Hamiltonian is the free Laplacian. But |x| is not in its domain, and the action is a delta function. The factorization Q^†Q with Q=∂−1/x corresponds to the Dirichlet Laplacian on R\\{0}, not the free Laplacian. The same issue affects |x|^g for g≤1 and |sinh|^g at coincidences. The paper does include a delta in the exp(g|x|) case, so the authors are aware of distributional terms; they just miss them for the power-law cases. The claim 'any pair function' needs to be replaced with a precise statement about domains and boundary conditions. Second, 'integrable' in the title is not earned. The paper establishes ground-state solvability, not integrability; either prove integrability or change the title. Third, the confined-case discussion repeats earlier results and is fine, but nothing more.\n\nWho is this for? People who like exact ground states and Calogero-type models. The framework is a useful generator of toy models, and with a careful treatment of singular interactions it could be solid. I would not cite the current version without the caveat.\n\nRecommendation: send to peer review, but choose a referee who knows the self-adjoint extension literature for Calogero models. The authors should be asked to fix the boundary-condition issue and to calibrate the integrability language. The general construction is worth preserving.","headline":"A useful framework for generating graph-based Jastrow models, but the domain issue for singular pair functions is real and the 'integrable' label overreaches; referee it, but expect revisions.","tokens_in":18283,"tokens_out":7378,"would_cite":true,"duration_ms":69435,"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":"For every graph and pair function, an edge-product wavefunction is the exact zero-energy ground state of an explicit Hamiltonian.","keywords":["Graph-Jastrow wavefunction","parent Hamiltonian","exactly solvable many-body models","adjacency matrix","three-body interactions","Calogero-Sutherland model","ground-state solvability","integrable systems"],"falsifier":"Take N=2 particles on a one-edge graph with f(r)=e^{g r}. The formal Hamiltonian predicts a repulsive delta interaction and zero ground-state energy. If a numerical diagonalization of a regularized version on a fine grid finds a negative energy eigenvalue, or if the delta contribution must be renormalized differently to produce a normalizable zero-energy state, the universal ground-state claim fails for this pair function.","tokens_in":17403,"feed_emoji":"🔗","tokens_out":7578,"duration_ms":72997,"temperature":0.7,"pith_summary":"This paper establishes a general recipe: any undirected graph G and any pair function f(r) define a many-body Hamiltonian for which the Jastrow-like product of f(r_ij) over the graph's edges is the exact zero-energy ground state. The parent Hamiltonian contains two-body interactions on the graph's edges and three-body interactions over its length-two paths, and the proof rests on factoring the Hamiltonian as a sum of squares of first-order operators that annihilate the wavefunction. This turns graph structure into a taxonomy of exactly solvable models, recovering known limits such as complete-graph quantum fluids and nearest-neighbor Jain-Khare-type systems, while generating new models on stars, wheels, ladders, and other graph products. The value is a systematic way to produce closed-form ground states, Hamiltonians, and energies for a broad class of distinguishable-particle systems.","feed_headline":"Every graph yields an exact many-body ground state","feed_subtitle":"A Jastrow product over graph edges is the zero-energy ground state of a Hamiltonian with two- and three-body interactions.","key_machinery":"The load-bearing object is the factorization H0 = ∑_i Q_i^† Q_i into first-order differential operators Q_i defined through the logarithmic derivative of the pair function and the graph's adjacency matrix. Each Q_i annihilates Φ0 because the gradient of an edge product picks up exactly the adjacency matrix entries, so H0 is positive semidefinite and Φ0 is automatically the ground state. The adjacency matrix A_ij controls which pairs enter the wavefunction and the potentials: edges generate two-body terms, and length-two paths (consecutive incident edges) generate three-body terms. Thus the graph's combinatorics determines the operator structure, and the pair function dictates the detailed fu","core_discovery":"The central claim is that the graph-Jastrow ansatz Φ0 = ∏_{(i,j)∈E} f(r_ij) is the exact ground state (with energy zero up to constant shifts) of the parent Hamiltonian H0 = -ħ²/2m ∑∇_i² + V2 + V3, where V2 couples pairs connected by an edge and V3 couples triples forming a two-edge path. The argument uses the identity H0 = ∑_i Q_i^† Q_i with Q_i = (ħ/√(2m))(∇_i − ∑_j A_ij (f'_ij/f_ij) r̂_ij), and Q_i Φ0 = 0, which makes H0 positive semidefinite. In one dimension the construction simplifies: the complete graph reproduces known exactly solvable gases (inverse-square, contact, and oscillator types), paths and cycles generalize the Jain-Khare model to arbitrary pair functions, 2r-regular graphs","pith_inferences":["The factorization suggests that ground-state properties such as correlation decay or entanglement may be controlled by graph invariants (degree distribution, number of 2-paths, spectral gap); comparing the same pair function on a path versus a complete graph would be a concrete test of this connectivity-controlled universality.","For singular pair functions like |x|^g or e^{g|x|}, the formal Hamiltonian contains delta functions and inverse-square terms; verifying the two-particle spectrum for a single-edge graph, where the self-adjoint extension issue is tractable, would determine whether the ground-state claim survives beyond formal manipulations.","The Q_i operators could seed a variational family: for graphs that are not exactly solvable, optimizing the pair function f would yield approximate ground states, extending the construction from exact to variational.","Since the three-body terms on 2-paths are required for exactness, a realistic simulator that implements only the two-body graph interactions will see a finite energy penalty; computing this penalty as a function of graph connectivity gives a testable prediction for analogue quantum simulators."],"forward_implications":["Every undirected graph and pair function yields a closed-form ground state and parent Hamiltonian, producing a broad class of exactly solvable many-body models without permutation symmetry.","The complete-graph limit recovers known exactly solvable gases (inverse-square, contact, harmonic) as special cases, showing a single construction unifies them.","Path, cycle, and 2r-regular graph models implement interactions truncated by adjacency rather than spatial range, interpolating between nearest-neighbor and all-to-all limits.","Graph joins and products generate composite system-plus-environment models (star, wheel, ladder, prism, Creutz ladder) whose ground states are exactly known, relevant for impurity and decoherence studies.","Weighted and random graph ensembles extend the framework to multi-species and disordered systems with site-dependent interaction strengths."],"fun_headline_variants":["Jastrow factor on graphs gives exact ground states","Every graph defines a solvable quantum many-body model","Graph Jastrow ansatz yields zero-energy ground states","Mapping graphs to exact many-body ground states","Graph edges become exact ground-state wave functions"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The proof that Φ0 is the ground state assumes H0 = ∑ Q_i^† Q_i defines a positive self-adjoint operator on a domain containing Φ0; for singular pair functions like |x|^g or e^{g|x|}, the derivatives and delta functions are only formal, so boundary contributions at particle coincidences could spoil positivity.","fun_headline_variants_meta":{"raw":{"variants":["Jastrow factor on graphs gives exact ground states","Every graph defines a solvable quantum many-body model","Graph Jastrow ansatz yields zero-energy ground states","Mapping graphs to exact many-body ground states","Graph edges become exact ground-state wave functions"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000458,"raw_usage":{"total_tokens":2146,"prompt_tokens":767,"completion_tokens":1379,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":511,"completion_tokens_details":{"reasoning_tokens":1306}},"tokens_in":511,"tokens_out":1379,"duration_ms":9687,"temperature":1.0,"reasoning_tokens":1306,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T20:43:19.874739+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take N=2 particles on a one-edge graph with f(r)=e^{g r}. The formal Hamiltonian predicts a repulsive delta interaction and zero ground-state energy. If a numerical diagonalization of a regularized version on a fine grid finds a negative energy eigenvalue, or if the delta contribution must be renormalized differently to produce a normalizable zero-energy state, the universal ground-state claim fails for this pair function.","supporting_citations":[],"review_version":1}