{"id":"d8bf7f21-0e4f-4550-8109-b4b07a383701","arxiv_id":"1908.08714","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"For strongly interacting trapped 1D bosons, the canonical Tan's contact at any temperature is conjectured to factor into a two-boson scaling function and a particle-number dependent Tonks-Girardeau contact.","lead":"This paper studies the Tan's contact, a measure of short-range pair correlations, for small numbers of repulsive bosons in a one-dimensional trap at finite temperature. It finds that in the strongly interacting regime the contact divided by its infinite-interaction value depends only on rescaled interaction and temperature, independent of particle number.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (23)'s conjectured N-dependence of the TG contact is unverified beyond N=5 and, taken literally, contradicts the paper's own exact high-T pair limit (Eq. 19); the 'any N, any T' extrapolation needs a direct check.","rationale":"The reader identified Eq. (23) as the weakest load-bearing premise, and I agree that the conjectured N-scaling of the TG contact is the soft spot for the paper's 'arbitrary N and T' formula. My stress-test adds a sharper, internal-consistency red flag: the formula as printed has a high-temperature limit incompatible with the paper's own exact pair-counting result Eq. (19). This makes the need for a direct check concrete rather than merely a request for more data. I do not regard this as fatal to the central scaling relation Eq. (25), because f_N(z,τ) is tested using the exact TG contact from Eq. (7), not the analytic s(N); the collapse in Figs. 4(d) and 5(d) is evidence for universality independent of Eq. (23). But the abstract and conclusion claim an explicit expression for the canonical contact at any N and T, and that claim is only as good as Eq. (23). Since the reader already rendered CONDITIONAL, I recommend no change in verdict, with the condition being a direct verification of Eq. (23) for larger N and over the full τ range, including the large-τ asymptotics. If Eq. (23) fails, the universal ratio picture may survive, but the explicit formula and the 'all N-dependence embedded in the TG contact' quantitative statement would need revision.","tokens_in":12663,"tokens_out":21439,"duration_ms":203182,"concrete_test":"Compute the exact canonical Tonks-Girardeau contact from Eqs. (7)-(9) for N=6 and N=7 over τ∈[0.05,20], and compare Cc_N(∞,τ)/h2(∞,τ) with the proposed s(N), checking in particular that the large-τ limit reproduces Eq. (19). If the printed s(N) fails at large τ, redo Figs. 3-5 with the corrected interpolation; if it passes for N=6 and N=7, the extrapolation is substantially supported. Optionally, run QMC for N=6 at z=1 and z=2.5 over τ∈[0.1,5] and test whether f_N(z,τ)/f_2(z,τ)=1 within error bars, which directly probes Eq. (25) without relying on the analytic s(N).","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's practical claim that the canonical contact is known for arbitrary N and T rests on Eq. (23), the conjectured factorization Cc_N(∞,τ)=h2(∞,τ)s(N). This conjecture is verified only for N=2..5 and for the QMC window τ≲5. As printed, s(N)=N^{5/2}-N^{3/4}(1+exp(-2/τ)) has the wrong high-temperature limit: at τ→∞ it tends to N^{5/2}-2N^{3/4}, whereas the exact pair-counting result Eq. (19) requires N^{5/2}-N^{3/2}. If the intended interpolation is instead s(N)=N^{5/2}-N^{3/4}-N^{3/2}exp(-2/τ), the high-T contradiction disappears, but the choice of the crossover function and the N^{3/4} exponent are still hand-fitted and have been tested only for N≤5. In either reading, the quantitative embedding statement is only as secure as this s(N): Fig. 3 uses s(N) to rescale the TG data, and the final formula uses s(N) to extrapolate to arbitrary N and T. A failure of s(N) for N>5, or an unverified crossover form, would not necessarily destroy the f_N universality claim of Eq. (25), but it would invalidate the paper's advertised explicit expression for the contact at arbitrary N and T.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the canonical Tan's contact of N repulsive Lieb-Liniger bosons in a one-dimensional harmonic trap at finite temperature, with a focus on small N. The authors compute the contact exactly for two bosons at arbitrary interaction and for N bosons in the Tonks-Girardeau (TG) limit, and use quantum Monte Carlo simulations for intermediate interactions. Their central claim is that in the strongly interacting regime (z > 1) the ratio f_N(z, τ) = C_N(g(z), T(τ))/C_N(∞, T(τ)) is approximately independent of N, so that all nontrivial particle-number dependence is carried by the TG contact. To make this explicit they conjecture an analytic N-dependence for the TG contact, Eq. (23), with an interpolating factor s(N) that bridges the zero-temperature N^{3/4} correction and the high-temperature pair-counting behavior. The paper also compares canonical and grand-canonical contacts and shows that the proposed scaling fails in the grand-canonical ensemble at intermediate interactions.","tokens_in":13098,"tokens_out":9128,"duration_ms":92793,"significance":"If the central universality claim were established, it would be a practically useful result: the canonical contact for a trapped Lieb-Liniger gas in the strong-coupling regime could be obtained from the exactly computable two-boson curve and the N-body TG contact. The paper has genuine strengths: the two-boson calculation and the TG contact evaluation are clean, the QMC data appear to have small errors, and the comparison between canonical and grand-canonical ensembles is informative. The proposed scaling relation is interesting and plausible. However, the main quantitative formula for the TG contact, Eq. (23), is internally inconsistent with the paper's own high-temperature limit, and the central universality claim is tested only in a narrow parameter window. These issues currently prevent the paper from supporting the advertised \"any N, any T\" conclusion.","major_comments":[{"comment":"The central scaling hypothesis (25) is tested only for N ≤ 5, z = 1 and 2.5, and rescaled temperatures roughly in the range 0.1 ≲ τ ≲ 1 in Figs. 4(d) and 5(d). The abstract and introduction state \"we show\" and \"for any number of particles and temperature,\" which overstates the evidence. Since Eq. (25) is presented as a conjecture, the language should be qualified throughout, and the conclusion should clearly delimit the tested regime. A direct test at larger N or wider τ, even for the TG contact where exact results are available, would materially strengthen the claim.","section":"Sec. III.C, Eq. (23)"}],"minor_comments":[{"comment":"The dashed line in Fig. 3 is labeled as the high-temperature limit h_2(∞, τ ≫ 1), but the plotted object appears to be h_N(∞, τ ≫ 1) of Eq. (20), which is not the τ → ∞ limit of h_2(∞, τ) as defined in Eq. (24). This labeling should be corrected to avoid confusion, especially given the issue raised in the first major comment.","section":"Sec. III.C, Eq. (24) and Fig. 3"},{"comment":"There are several typos: \"g is ininite\" should be \"g is infinite,\" and \"can de derived\" should be \"can be derived.\" These should be corrected.","section":"Sec. II.B and Sec. III.C"},{"comment":"The horizontal axes of Figs. 4 and 5 are labeled only \"0.1 1\"; the tick labels and the intended τ range should be made explicit so the reader can see exactly which temperatures are covered.","section":"Figs. 4 and 5"},{"comment":"Reference [3] (Yang and Yang) appears to have an incorrect journal and volume: the standard citation is J. Math. Phys. 10, 1115 (1969). Please verify all bibliographic entries.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper contains solid exact two-body and TG calculations, and the scaling idea is worth publishing if the presentation is made accurate. The main problem is that Eq. (23) has a large-τ limit that contradicts the paper's own exact pair-counting result; this looks like a fixable interpolation error, but as written it undermines the advertised explicit formula. The abstract and conclusion should also be brought in line with the status of Eq. (25) as a conjecture supported by limited numerical evidence."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The useful core is the finite-temperature canonical scaling collapse for small N, and it is worth a careful look. The exact two-boson contact, the exact Tonks-Girardeau contact from the fermionic two-body density matrix, and the QMC data with honest error bars all look clean. The collapse in Figs. 4 and 5 for z = 1 and 2.5, N = 2..5, over τ ≈ 0.1–1, is a genuine numerical check of the embedding idea. What is new relative to [37] and the grand-canonical work is precisely this finite-T canonical statement and the explicit conjectured form for the TG contact.\n\nThe soft spot is Eq. (23), and the stress-test note is right about it. As printed, s(N) = N^{5/2} − N^{3/4}(1 + exp(−2/τ)) tends to N^{5/2} − 2N^{3/4} at large τ, whereas the paper's own exact pair-counting result, Eq. (19), requires N^{5/2} − N^{3/2}. The likely intended form is s(N) = N^{5/2} − N^{3/4} − N^{3/2} exp(−2/τ), but that is not what is written, and either way the crossover is hand-fitted and tested only for N ≤ 5. This is load-bearing for the advertised explicit 'any N, any T' expression, but not for the f_N collapse itself, because the denominators in Figs. 4(d) and 5(d) are exact C_N(∞,τ) values, not Eq. (23). The mild circularity—the N = 2 branch is built into h_2 = C_2/s(2), and s(N) was chosen to reproduce the TG data before testing—does not sink the main collapse, but it does weaken the extrapolation.\n\nThe paper itself calls Eq. (23) a conjecture, but the abstract says 'we show', which overstates what is established. A derivation or a direct check for N = 6–8, plus a corrected high-temperature limit, would materially strengthen the claim. Proportional to its soft spots, this is a solid paper with one real inconsistency in an auxiliary formula.\n\nWho is it for? People working on Tan's contact in 1D traps, few-body ultracold gases, and canonical versus grand-canonical ensemble effects. I would send it to peer review; a serious referee should ask for the Eq. (23) fix, a larger-N check if feasible, and a more careful abstract. I would not cite the explicit formula in my own work until it is corrected.","headline":"Useful, mostly solid paper with a real high-temperature inconsistency in the advertised explicit formula; the canonical scaling collapse for N≤5 is credible.","tokens_in":13579,"tokens_out":3704,"would_cite":false,"duration_ms":38257,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81V70","82B10","82B80"],"pacs":["67.85.-d","05.30.Jp"],"model":"deepseek-v4-flash","headline":"In the strongly interacting regime, every N-particle contact ratio collapses onto the two-boson curve.","keywords":["Tan's contact","Lieb-Liniger gas","Tonks-Girardeau limit","canonical ensemble","finite temperature","one-dimensional bosons","scaling function","quantum Monte Carlo"],"falsifier":"Evaluate the exact fermionic density-matrix expression for the canonical Tonks-Girardeau contact, Eq. (7), for $N=6$ and $N=7$ over a range of $\\tau$; if $C_N^c(\\infty,\\tau)/h_2(\\infty,\\tau)$ departs from $s(N)$ or gains a visible $\\tau$ dependence, Eq. (23) is wrong. A second check is to compute $f_N(z=2.5,\\tau)$ for $N=6$ by quantum Monte Carlo and ask whether it still lies on the $N=2$ curve.","tokens_in":12483,"feed_emoji":"⚛️","tokens_out":8792,"duration_ms":77564,"temperature":0.7,"pith_summary":"Strongly interacting one-dimensional trapped Bose gases are usually studied one system at a time, with a fixed particle number, interaction, and temperature. This paper claims that in the strongly interacting regime the canonical Tan's contact for any number $N$ of bosons is fixed by two simpler quantities: the exactly computable two-boson contact and the $N$-body contact at infinite repulsion, the Tonks-Girardeau limit. Concretely, the ratio $C_N(g,T)/C_N(\\infty,T)$ is claimed to be a universal function of the rescaled interaction strength $z>1$ and rescaled temperature $\\tau$, with no remaining $N$ dependence. The paper also conjectures an explicit formula for the $N$ dependence of the Tonks-Girardeau contact at any temperature, verified numerically for $N=2$ to $5$. If the claim holds, the full contact of a few-atom experiment is obtained from the two-boson curve and the known Tonks limit, with correlations and pair effects separated cleanly.","feed_headline":"One curve fixes the contact of any strongly interacting 1D Bose gas","feed_subtitle":"For fixed particle number, the Tan contact is set by the two-boson curve and the Tonks-Girardeau limit.","key_machinery":"The mechanism is a chain of exact reductions. The canonical contact is obtained from Tan's sweep relation, $C_N^c = -(m^2/\\pi\\hbar^4)\\partial F/\\partial g^{-1}$, equivalently $C_N^c = (g m^2/\\pi\\hbar^4)\\langle H_{\\mathrm{int}}\\rangle$. For two bosons the relative spectrum comes from the implicit equation $\\Gamma(-\\nu/2)/\\Gamma(-\\nu/2+1/2) = -\\sqrt{2}\\,|a_{1\\mathrm D}|/a_{\\mathrm{ho}}$, giving the exact curve $C_2^c(z,\\tau)$ and its Tonks-Girardeau limit. For $N$ bosons at infinite repulsion the contact is written through the fermionic two-body density matrix, $C_N^c(\\infty,T) = (2/\\pi)\\int dx\\,F(x)$, and for finite interactions the authors use quantum Monte Carlo on a discretized Hubbard model. The identity that carries the claim is $f_N(z>1,\\tau)\\simeq f_2(z>1,\\tau)$, and the conjectured $N$-dependent factor $s(N)=N^{5/2}-N^{3/4}(1+e^{-2/\\tau})$ is what makes the Tonks-Girardeau data for different $N$ collapse onto a single curve.","core_discovery":"The central claim is the scaling law $f_N(z>1,\\tau) \\simeq f_2(z>1,\\tau)$, where $f_N = C_N^c(z,\\tau)/C_N^c(\\infty,\\tau)$ is the canonical contact rescaled by its value at infinite repulsion, $z = a_{\\mathrm{ho}}/(|a_{1\\mathrm D}|\\sqrt{N})$ is the rescaled interaction, and $\\tau=T/T_F$ the rescaled temperature. Equivalently, all the non-trivial particle-number dependence of the contact is embedded in the Tonks-Girardeau contact $C_N^c(\\infty,\\tau)$. For that quantity the paper conjectures $C_N^c(\\infty,\\tau) = h_2(\\infty,\\tau)\\,[N^{5/2}-N^{3/4}(1+e^{-2/\\tau})]$, with $h_2$ fixed by the exact two-boson formula; Monte Carlo data for $N=2$ to $5$ collapse on this curve within about $5\\%$ at $z=1$ and $1\\%$ at $z=2.5$. The same collapse is shown to fail in the grand-canonical ensemble at intermediate temperatures, so the universal form is specific to fixed particle number.","pith_inferences":["I would expect the $3/4$ exponent to be a fingerprint of the harmonic confinement, so a box-trapped version of this construction should show a different finite-$N$ factor; checking the homogeneous analogue cited in the paper would settle that.","The factor $1+e^{-2/\\tau}$ is an interpolation between the low- and high-temperature limits; deriving it from the fermionic density-matrix sums might extend the collapse to weaker interactions ($z\\le 1$) or to multi-component mixtures.","A practical experimental test is to measure the momentum-tail coefficient in a few-atom 1D tube for two or three particle numbers at fixed $z$ and $\\tau$, and compare the ratio $C_N/C_N(\\infty)$ with the two-boson ratio."],"forward_implications":["For every $N\\ge 2$ and every temperature in the regime $z>1$, the canonical contact is $C_N^c(z,\\tau)=f_2(z,\\tau)\\,C_N^c(\\infty,\\tau)$, so no separate $N$-body calculation is required once the two-boson curve and the Tonks-Girardeau contact are known.","The conjectured formula $C_N^c(\\infty,\\tau)=h_2(\\infty,\\tau)(N^{5/2}-N^{3/4}(1+e^{-2/\\tau}))$ supplies an explicit bridge from few-atom samples to the thermodynamic limit, where $N^{5/2}$ scaling is recovered.","At high temperature the formula reduces to the pair-counting result $N(N-1)/2$ times the two-boson contact, while at low temperature it carries the correlation-induced $N^{5/2}-N^{3/4}$ dependence.","The universal collapse does not survive in the grand-canonical ensemble at intermediate temperatures, so fixed-$N$ (canonical) experiments are required to observe it.","The authors note that an analogous universal ratio follows from recent results for the homogeneous gas, suggesting the scaling extends beyond harmonic traps."],"supporting_citations":[{"why":"Tan's sweep relation defines the contact as $-(m^2/\\pi\\hbar^4)\\partial F/\\partial g^{-1}$, the starting point of the calculation.","marker":"[20]"},{"why":"Supplies the exact two-boson relative spectrum through the implicit equation for $\\nu$, from which $C_2^c(z,\\tau)$ is computed.","marker":"[12]"},{"why":"Establishes the zero-temperature scaling $f_N(z,0)\\simeq f_2(z,0)$ and the $N^{5/2}-N^{3/4}$ behaviour that the conjecture extends to finite temperature.","marker":"[37]"},{"why":"Gives the Tonks-Girardeau contact as an integral of the fermionic two-body density matrix, used in Eqs. (7)-(9).","marker":"[45]"},{"why":"Supplies the grand-canonical Tonks-Girardeau contact formulas used to compare statistical ensembles.","marker":"[46]"},{"why":"Provides the large-temperature grand-canonical contact and virial scaling used to identify the pair-counting limit.","marker":"[16]"},{"why":"Presents the stochastic Green function algorithm used for the canonical quantum Monte Carlo simulations.","marker":"[47]"},{"why":"Extends the stochastic Green function algorithm to the parameter regimes needed for the finite-temperature contact calculations.","marker":"[48]"}],"fun_headline_variants":["Tan contact scaling: all N collapse onto two-boson curve","One curve fixes Tan contact for any N in 1D Bose gas","Universal Tan contact: N-dependence lives in Tonks limit","Fixed-N Tan contact reduces to two-particle scaling","Strongly interacting 1D gas: contact scales with two-body curve"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is the conjectured formula for the Tonks-Girardeau contact, $C_N^c(\\infty,\\tau)=h_2(\\infty,\\tau)(N^{5/2}-N^{3/4}(1+e^{-2/\\tau}))$, whose exponents come from a zero-temperature fit and large-temperature pair counting, whose exponential crossover is put in by hand, and which is verified only for $N=2$ to $5$.","fun_headline_variants_meta":{"raw":{"variants":["Tan contact scaling: all N collapse onto two-boson curve","One curve fixes Tan contact for any N in 1D Bose gas","Universal Tan contact: N-dependence lives in Tonks limit","Fixed-N Tan contact reduces to two-particle scaling","Strongly interacting 1D gas: contact scales with two-body curve"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000237,"raw_usage":{"total_tokens":1494,"prompt_tokens":922,"completion_tokens":572,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":538,"completion_tokens_details":{"reasoning_tokens":484}},"tokens_in":538,"tokens_out":572,"duration_ms":6611,"temperature":1.0,"reasoning_tokens":484,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:30:52.522526+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate the exact fermionic density-matrix expression for the canonical Tonks-Girardeau contact, Eq. (7), for $N=6$ and $N=7$ over a range of $\\tau$; if $C_N^c(\\infty,\\tau)/h_2(\\infty,\\tau)$ departs from $s(N)$ or gains a visible $\\tau$ dependence, Eq. (23) is wrong. A second check is to compute $f_N(z=2.5,\\tau)$ for $N=6$ by quantum Monte Carlo and ask whether it still lies on the $N=2$ curve.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the exact two-boson relative spectrum through the implicit equation for $\\nu$, from which $C_2^c(z,\\tau)$ is computed."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the zero-temperature scaling $f_N(z,0)\\simeq f_2(z,0)$ and the $N^{5/2}-N^{3/4}$ behaviour that the conjecture extends to finite temperature."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the Tonks-Girardeau contact as an integral of the fermionic two-body density matrix, used in Eqs. (7)-(9)."},{"cited_title":"Yan et al","cited_arxiv_id":null,"evidence_quote":"Supplies the grand-canonical Tonks-Girardeau contact formulas used to compare statistical ensembles."},{"cited_title":"5, τ )a3 ho τ FIG","cited_arxiv_id":null,"evidence_quote":"Provides the large-temperature grand-canonical contact and virial scaling used to identify the pair-counting limit."},{"cited_title":"Hoinka et al","cited_arxiv_id":null,"evidence_quote":"Presents the stochastic Green function algorithm used for the canonical quantum Monte Carlo simulations."},{"cited_title":"Laurent et al","cited_arxiv_id":null,"evidence_quote":"Extends the stochastic Green function algorithm to the parameter regimes needed for the finite-temperature contact calculations."}],"review_version":1}