{"id":"c03684b6-9f0a-44aa-ac0d-ae111092f0ff","arxiv_id":"2502.06310","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Exact closed-form natural orbital occupancies and collective angular-momentum occupancies are derived for the 2D Moshinsky model, with a new identity K = κ_l K_η.","lead":"This paper derives exact formulas for how many particles occupy each angular momentum state in the ground state of a solvable 2D model of trapped interacting bosons. The formulas provide exact benchmarks for approximate many-body methods and reveal a structural identity linking two measures of correlation.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The exact λ_nl and η_l formulas are correct, but the abstract's claim that natural orbitals are 'uniformly distributed across all significant l components' is unsupported: η_l decays geometrically with |l|.","rationale":"The reader identified the imported 1-RDM kernel as the weakest assumption, but I verified it independently via normal-mode decomposition for general N: S = ωI + ((1-ω)/N)J, and the resulting Gaussian after integrating out N-1 particles matches Eqs. (5)-(7) exactly, including the normalization A = ω/(πγ). The Hardy-Hille application is also correct: the matching conditions (15) follow from the exponent and Bessel argument, and the eigenvalues λ_nl = Aπ t^{n+|l|/2}(1-t)/z² satisfy the trace condition Σ_{l=-∞}^{∞}η_l = 1. The reader's claimed numerical flaws are not supported. For fixed Λ and large N, λ_00 ≈ 1 - √(Λ/(2N)) is the correct leading depletion: Aπ/z² ≈ (ω/γ)/(ω/γ + C/2) ≈ 1 - Λ/ω ≈ 1 - √(Λ/(2N)), and the (1-t) correction is subleading. For Λ=1, numerical evaluation of K_η from the exact formulas gives K_η(2)≈1.17, K_η(3)≈1.22, K_η(4)≈1.24, K_η(5)≈1.25, K_η(6)≈1.25, K_η(10)≈1.24, K_η(100)≈1.13, so a maximum near N≈5-7 clearly exists. The genuine issue is the interpretation: Eq. (30) only establishes l-independent radial participation, not uniform population across l. Since the central exact formulas stand but the advertised 'unique feature' is overstated, the CONDITIONAL verdict remains appropriate, though the required revision is in the abstract and discussion rather than in the derivations.","tokens_in":6566,"tokens_out":61947,"duration_ms":426919,"concrete_test":"Using Eqs. (23) and (25), compute η_l and K_η for a representative strong-interaction case, e.g., N=500 and Λ=10^5. Since η_{l+1}/η_l = √t < 1, the occupancies decay geometrically; check whether the first K_η angular-momentum sectors each carry approximately 1/K_η of the population as 'uniform distribution' would imply. If they do not (e.g., η_0 is several times η_{K_η-1}), the abstract's claim is contradicted.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The core mathematical result is sound: I re-derived the 1-RDM kernel from the normal-mode decomposition and confirmed Eqs. (5)-(7), and the Hardy-Hille matching in Eqs. (12)-(18) is algebraically consistent and satisfies the trace normalization. The load-bearing concern is the paper's headline interpretation. The abstract and Section 3 claim that natural orbitals are 'uniformly distributed across all significant l components' and that 'natural orbitals contributing to the correlations are derived equally from all significant l components.' The only supporting result, K = κ_l K_η (Eq. 30), shows that the effective number of radial modes within each l sector, κ_l = (1+t)/(1-t), is independent of l. It does not mean the collective occupancies η_l are spread uniformly: by Eq. (23), η_l = πA t^{|l|/2}/z² decays geometrically in |l| because 0<t<1. Thus the population is concentrated at small |l|, and the significant l components do not contribute equally. The reader's two specific objections (that λ_00 ≈ 1 - √(Λ/2N) is false and that K_η(N) has no maximum for Λ=1) are not valid: the λ_00 formula is the correct leading-order depletion for fixed Λ and large N, and a direct calculation for Λ=1 gives K_η increasing from N=2 to a maximum near N≈6 before decreasing, so the maximum does exist. The real flaw is the overstatement of the 'uniform distribution' result, which requires correction in the abstract and discussion.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript considers N bosons in a two-dimensional isotropic harmonic trap with pairwise harmonic interactions (the Moshinsky model). Starting from the polar-coordinate one-particle density matrix (5) with coefficients (6)-(7), it derives a Schmidt decomposition (19)-(20) by Fourier-Lagrange expansion and the Hardy-Hille formula. The central output is the occupancy formula λ_nl = Aπ t^{n+|l|/2}(1-t)/z^2 and the collective occupancy η_l = πA t^{|l|/2}/z^2, together with closed forms for the participations K, K_η, and κ_l. Section 3 uses these formulas to discuss fragmentation and argues that correlations are uniformly distributed over significant angular-momentum sectors.","tokens_in":6820,"tokens_out":44831,"duration_ms":356556,"significance":"If correct, the closed-form natural orbitals and occupancies would provide a useful benchmark for the harmonic-interaction model, and the Hardy-Hille diagonalization is elegant, explicit, and trace-normalized. The derivation is fully analytic and involves no fitted parameters. However, the coefficient C in Eq. (6) appears not to be the exact one-particle reduced density-matrix coefficient of the stated Hamiltonian, and the headline interpretation of Eq. (30) is not supported by Eq. (23). The significance is therefore conditional on a correction of the input kernel and a revision of the abstract and discussion.","major_comments":[{"comment":"The coefficient C in Eq. (6) is not the exact 1-RDM coefficient of the Moshinsky Hamiltonian. For N=2, the exact ground state after separating center-of-mass and relative coordinates is Ψ ∝ exp[-(1+ω)/4(r1²+r2²)+(ω-1)/2 r1·r2] with ω=√(1+4Λ); integrating out one particle gives B=(1+6ω+ω²)/(4(1+ω)) and C=(ω-1)²/(2(1+ω)) in the notation of Eq. (5). The printed Eq. (6) gives instead C=(ω-2)²/(2(1+ω)) and B=(ω+2)²/(4(1+ω)) for N=2, which agree with the exact values only at isolated points. In particular, at Λ=0 the printed C is 1/4 rather than 0, so the kernel does not reduce to the noninteracting harmonic-oscillator 1-RDM. The general corrected form appears to be C=(ω-1)²(N-1)/(N²γ). This error affects all numerical results and the physical conclusions in Section 3.","section":"2.2, Eq. (6)"},{"comment":"The statement λ00 ≈ 1 - √(Λ/(2N)) is inconsistent with the manuscript's own equations as printed. Using Eq. (18) with the printed coefficient C from Eq. (6), a fixed-Λ, large-N expansion gives λ00 of order N^{-1/2}, tending to zero, not to 1. If the coefficient C is corrected as noted in the first major comment, the formula may become the correct leading-order depletion, but as submitted the text and the equations disagree. This should be reconciled explicitly.","section":"3, condensation statement"},{"comment":"The claim that natural orbitals are 'uniformly distributed across all significant l components' and that they 'originate equally from all relevant l fragments' is not supported. Eq. (23) gives η_l = πA t^{|l|/2}/z² with 0<t<1, a geometric decay in |l|, so the population is concentrated at small |l|. Eq. (30), K = κ_l K_η, only shows that the radial participation within each l-sector is independent of l; it does not imply that different l sectors contribute equal total populations. The abstract and the discussion should either remove this claim or redefine 'uniform' precisely as uniformity of the radial participation within each l sector.","section":"Abstract and Sec. 3"}],"minor_comments":[{"comment":"The heading 'Discusions' should be corrected to 'Discussion'.","section":"3, heading"},{"comment":"The spelling 'Moshynsky' appears in the introduction and should be 'Moshinsky'.","section":"1, Introduction"},{"comment":"Reference [8] is incomplete: the article number for Phys. Rev. A 101 (2020) is missing.","section":"References"},{"comment":"The asymptotic expression (27) can exceed 1 for moderate values of Λ even when the exact η_l cannot; for example, at N=500 and Λ=100 the leading term is above 1. The range of validity of the approximation should be stated.","section":"3, Eq. (27)"},{"comment":"After Eq. (13), it would be helpful to state explicitly that the root t is chosen in (0,1) so that the geometric sums leading to Eq. (23) converge and the collective occupancies are positive.","section":"2.3, Eq. (13)"}],"recommendation":"major_revision","confidential_remarks":"The error in Eq. (6) may be a typo, but it is load-bearing because the exactness claim of the paper depends on the input kernel. If the authors correct C to C=(ω-1)²(N-1)/(N²γ), the Hardy-Hille machinery and the closed-form formulas will survive, but all figures and numerical statements must be recomputed. In addition, the 'uniform distribution' claim in the abstract and Section 3 should be substantially softened, since the collective occupancies decay geometrically with |l|."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: the exact results are correct and worth publishing; the abstract overclaims. The polar-coordinate Schmidt decomposition of the 1-RDM, the closed-form λ_nl, and the collective-occupancy formula η_l = πA t^{|l|/2}/z² are new relative to the cited Cartesian treatment, and they are derived cleanly from the Hardy–Hille formula. I re-derived the input kernel from normal modes and checked the trace and normalization identities; Section 2 is sound. The paper also does well to report the asymptotic K_η ~ Λ^{1/4} growth and the finite-N behavior, which are useful benchmarks for approximate methods.\n\nThe reader's conditional verdict is reasonable in spirit, but two of the specific objections in the report do not survive contact with the equations. The λ_00 ≈ 1 − sqrt(Λ/2N) formula is not in conflict with the paper's own derivation; for fixed Λ and N→∞ it is the correct leading depletion. And for Λ=1, K_η(N) does have a maximum near N≈6 before decreasing, so that criticism misses. The actual soft spot is the 'uniform distribution' language. The abstract and Section 3 claim that natural orbitals contributing to correlations are uniformly distributed across significant l components. The evidence is K = κ_l K_η with κ_l = (1+t)/(1-t) independent of l. That says the effective number of radial modes per l-sector is the same across sectors; it does not imply the collective occupancies η_l are uniformly spread. Since η_l decays geometrically in |l| at fixed Λ, the angular population is concentrated at small |l|. The word 'uniformly' is not supported and should be removed, with the conclusions rephrased as: within each l-sector, radial participation is l-independent, while the angular weight still falls off geometrically. This is a wording and interpretation defect, not a load-bearing mathematical flaw.\n\nRecommendation: send it to peer review. The core derivation is solid, the model is explicitly harmonic so its physical reach is modest, but exact 2D fragmentation benchmarks are genuinely scarce. A referee can require the abstract and discussion to be recalibrated, and the paper will then be a clean, citable contribution.","headline":"Core formulas are right and worth publishing, but the abstract's 'uniform distribution' claim overstates what K = κ_l K_η actually shows.","tokens_in":7445,"tokens_out":3799,"would_cite":false,"duration_ms":596267,"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":"Collective occupancies of the 2D Moshinsky model derived exactly","keywords":["Moshinsky model","harmonium","natural orbitals","collective occupancy","boson fragmentation","participation ratio","Hardy-Hille formula","two-dimensional harmonic trap"],"falsifier":"Take a specific $N$ and interaction strength $\\Lambda$, build the many-body Hamiltonian in a sufficiently large truncated single-particle basis, diagonalize it numerically, extract the one-particle reduced density matrix from the ground state, and compare its collective occupancies $\\eta_0$, $\\eta_1$, $\\eta_2$ with the closed formula (23); a systematic disagreement beyond the truncation error would falsify the claimed exact result.","tokens_in":6307,"feed_emoji":"⚛️","tokens_out":8822,"duration_ms":66531,"temperature":0.7,"pith_summary":"This paper proves that for the ground state of $N$ bosons in a two-dimensional harmonic trap with harmonic interactions, the one-particle reduced density matrix can be diagonalized exactly in polar coordinates. Each natural orbital is labelled by a radial quantum number $n$ and an angular momentum $l$, and its occupancy is exactly $\\lambda_{nl} = A\\pi t^{n+|l|/2}(1-t)/z^2$. Summing over $n$ gives the collective occupancy $\\eta_l = \\pi A t^{|l|/2}/z^2$, the fraction of particles carrying angular momentum $l$. These closed forms make the full interaction regime and any particle number analytically accessible, and they show that fragmentation into angular-momentum components grows with interaction strength.","feed_headline":"Collective occupancies of the 2D Moshinsky model derived exactly","feed_subtitle":"Closed-form formula for η_l shows fragmentation into angular-momentum sectors grows with interaction strength.","key_machinery":"The workhorse is the Hardy--Hille formula, the bilinear generating function for generalized Laguerre polynomials, which reproduces the Gaussian-times-Bessel kernel appearing in the polar-coordinate 1-RDM. Matching the kernel parameters via $B=(1+2t/(1-t))z^2$ and $C=4\\sqrt{t}(1-t)^{-1}z^2$, where $z^2=(4B^2-C^2)^{1/2}/2$, turns each angular-momentum partial wave into a diagonal sum over Laguerre-based natural orbitals, so occupancies come out as explicit powers of $t$. The same geometric structure then yields every collective quantity by summation.","core_discovery":"The central discovery is that the polar-coordinate one-particle reduced density matrix, a Gaussian kernel with an angular Bessel factor, is diagonalized by the Hardy--Hille formula for Laguerre polynomials, yielding explicit natural orbitals $u_{nl}(r,\\phi)\\propto (zr)^{|l|} e^{-z^2r^2/2} L_n^{|l|}(z^2r^2) e^{il\\phi}$ and occupancies $\\lambda_{nl}=A\\pi t^{n+|l|/2}(1-t)/z^2$. Because the occupancies are geometric in $n$, the collective occupancy $\\eta_l=\\pi A t^{|l|/2}/z^2$ and the participation measures $K_\\eta=z^4(1-t)/(\\pi^2 A^2(1+t))$, $K=z^4/(\\pi^2 A^2)$, and $\\kappa_l=(1+t)/(1-t)$ follow in closed form. The paper further establishes the identity $K=\\kappa_l K_\\eta$, showing that the natural orbitals contributing to the correlations are drawn uniformly from all significant angular-momentum sectors. In the large-interaction limit the effective number of $l$-fragments grows as $\\Lambda^{1/4}$, while for a macroscopic number of particles the zero-angular-momentum orbital becomes macroscopically occupied, signalling condensation.","pith_inferences":["Because the derivation only needs the Gaussian-Bessel form of the 1-RDM, the same Hardy-Hille diagonalization should go through for excited states or anisotropic variants of the Moshinsky model whose reduced density matrix retains that form; this is a direct extension the paper does not pursue.","The $\\Lambda^{1/4}$ growth of $K_\\eta$ suggests a measurable signature: in a trapped ultracold gas, the number of populated angular-momentum sectors could be inferred from momentum-space noise correlations, providing an experimental test of the fragmentation law.","The uniform-fragmentation identity $K=\\kappa_l K_\\eta$ implies that truncating the natural-orbital expansion to a few $l$-sectors introduces errors of the same relative size in every kept sector; numerical simulations of such systems should therefore include complete $l$-shells rather than individual orbitals."],"forward_implications":["For any particle number $N$, the collective occupancy $\\eta_l$ decreases geometrically with $|l|$, so the angular-momentum distribution of the gas is fully characterized by one number $t$.","The participation $K_\\eta$ grows like $\\Lambda^{1/4}$, so the number of significantly populated angular-momentum fragments increases with interaction strength at a quarter-power rate.","At strong interactions the state fragments into many macroscopically occupied $l$ components, while at large $N$ the $l=0$ orbital approaches unit occupancy and the system condenses.","The identity $K=\\kappa_l K_\\eta$ means that every angular-momentum sector contributes the same effective number of natural orbitals to the correlations, a property that distinguishes the Moshinsky model from generic many-body systems.","Closed-form occupancies provide exact benchmarks for approximate methods in two-dimensional trapped boson systems."],"supporting_citations":[{"why":"Supplies the exact ground-state 1-RDM in Cartesian coordinates, which the paper rewrites in polar form as Eq. (5); the entire diagonalization starts from this kernel.","marker":"[18]"},{"why":"Provides the Hardy-Hille bilinear generating function for Laguerre polynomials used to diagonalize each angular-momentum partial wave.","marker":"[20]"},{"why":"Supplies the integral identity for the modified Bessel function $I_l(z)$ used to extract the partial-wave components of the 1-RDM.","marker":"[19]"},{"why":"Defines the collective occupancy as the sum of occupancies with a given angular momentum, the central quantity the paper evaluates exactly.","marker":"[22]"},{"why":"Introduces the participation ratio $K$ as an indicator of single-particle correlations, which the paper evaluates in closed form.","marker":"[23]"}],"fun_headline_variants":["Exact collective occupancies of 2D Moshinsky model","Hardy-Hille diagonalization gives exact occupancies","Angular-momentum fragmentation exactly solved","Uniform distribution across angular momentum sectors","Closed-form collective occupancy for 2D trapped bosons"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire derivation assumes that the Gaussian kernel of Eq. (5), imported from reference [18], is the exact one-particle reduced density matrix of the ground state; if that kernel were wrong, the Hardy-Hille diagonalization would produce orbitals and occupancies for the wrong state.","fun_headline_variants_meta":{"raw":{"variants":["Exact collective occupancies of 2D Moshinsky model","Hardy-Hille diagonalization gives exact occupancies","Angular-momentum fragmentation exactly solved","Uniform distribution across angular momentum sectors","Closed-form collective occupancy for 2D trapped bosons"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00031,"raw_usage":{"total_tokens":1785,"prompt_tokens":976,"completion_tokens":809,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":592,"completion_tokens_details":{"reasoning_tokens":735}},"tokens_in":592,"tokens_out":809,"duration_ms":7042,"temperature":1.0,"reasoning_tokens":735,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T15:54:09.033023+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a specific $N$ and interaction strength $\\Lambda$, build the many-body Hamiltonian in a sufficiently large truncated single-particle basis, diagonalize it numerically, extract the one-particle reduced density matrix from the ground state, and compare its collective occupancies $\\eta_0$, $\\eta_1$, $\\eta_2$ with the closed formula (23); a systematic disagreement beyond the truncation error would falsify the claimed exact result.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the exact ground-state 1-RDM in Cartesian coordinates, which the paper rewrites in polar form as Eq. (5); the entire diagonalization starts from this kernel."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Hardy-Hille bilinear generating function for Laguerre polynomials used to diagonalize each angular-momentum partial wave."},{"cited_title":"Abramowitz, I","cited_arxiv_id":null,"evidence_quote":"Supplies the integral identity for the modified Bessel function $I_l(z)$ used to extract the partial-wave components of the 1-RDM."},{"cited_title":"Ko´ scik, Quantum Information Processing, 23(7), 260 (2024)","cited_arxiv_id":null,"evidence_quote":"Defines the collective occupancy as the sum of occupancies with a given angular momentum, the central quantity the paper evaluates exactly."},{"cited_title":"Grobe, K","cited_arxiv_id":null,"evidence_quote":"Introduces the participation ratio $K$ as an indicator of single-particle correlations, which the paper evaluates in closed form."}],"review_version":1}