{"id":"5ca70bc8-7a30-4955-8d42-9b92828cb643","arxiv_id":"2412.11686","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"An analytic continuation via contour integrals gives the emptiness instanton and emptiness formation probability for all polytropic indices gamma > 1 in a 1D quantum gas.","lead":"The authors find the exact shape of the emptiness instanton, the rare density ripple that empties a long interval in a one-dimensional quantum polytropic gas, valid for every polytropic index gamma > 1. Generalists might read this because it reduces a previously special-case large-deviation problem to a closed integral formula and confirms a conjecture about interacting quantum fluids.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Branch mismatch admitted in Appendix B leaves the global spatiotemporal instanton and the x=-1 boundary condition unproven; the full-profile claim needs a consistency check.","rationale":"The reader's weakest_assumption already identifies the branch-patching gap, and I agree that it is the most load-bearing issue. I weighed whether the EFP exponent itself could be invalidated: the amplitude α in Eq. (39) is computed from Eq. (38) at λ=-\\bar λ=iμ, μ→1+, which is the boundary of the domain Re λ=v>0 where Eq. (36) was derived; the integral in Eq. (37) is regular for μ>1 and the hypergeometric asymptotic is standard, so the EFP exponent f(n) is credible. But the paper's central claim explicitly includes the analytic spatiotemporal profile (Fig. 1), and that claim depends on patching two representations whose global equivalence the authors explicitly disclaim in Appendix B. Without a consistency check, a reader cannot tell whether the plotted interior is the physical instanton or an artifact of branch choices. The concern is concrete and testable, hence conditional rather than rejection. Credit is due for the independent boundary/asymptotic verifications and for the exact axis profiles at x=0 and τ=0, which support the EFP part of the claim.","tokens_in":14473,"tokens_out":23587,"duration_ms":213188,"concrete_test":"Numerically evaluate ∂λV_n from Eq. (36) and from Eq. (43) for n=1/2 and n=√2 at a set of conjugate points λ=\\bar λ* with v=Re λ<0 and Im λ from 0.5 to 2, covering the ρ=1 crossing in the left half-plane. Insert each value into Eqs. (41)-(42) and check whether the two representations give the same (x,τ). If they differ anywhere on this grid, the Fig. 1 profile is not a single-valued solution and the branch-patching claim fails; if they agree up to numerical precision, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Eqs. (34)-(35) are claimed to give the emptiness instanton for all n>-1/2, including the full spatiotemporal profile plotted in Fig. 1. Appendix B ends by stating that the two representations used for the profile, Eq. (36) (said to be suitable for ρ>1) and Eq. (43) (for ρ<1), can lie on different branches and that their equivalence 'can be violated' outside Re λ = Re \\bar λ = v > 0. The physical profile, however, includes the left half of the astroid where v<0, i.e., Re λ<0; the paper does not prove that the patched branches agree where the density crosses ρ=1, nor does it prove the parity/branch choice that would give boundary condition (22) at x=-1 (Re λ<0). Thus the off-axis density shapes, the n-dependent exponents at x=±1, and the claim that Eqs. (34)-(35) define a single global solution are not established. The EFP exponent f(n), Eq. (3), is less exposed because it follows from the μ→1+ amplitude on the axis λ=-\\bar λ=iμ, but the global instanton claim is load-bearing on the unproven branch consistency.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper constructs a Pochhammer contour integral, Eq. (34), for the potential V_n governing the imaginary-time hydrodynamic description of a one-dimensional polytropic gas with equation of state P ~ rho^gamma, gamma = 1 + 2/(2n+1), and claims that this potential solves the emptiness instanton problem for all n > -1/2. The construction verifies the Euler-Poisson equation by linearity, checks the endpoint and quadrupole boundary conditions using Eq. (36) and hypergeometric asymptotics, and rederives the emptiness formation amplitude alpha in Eq. (39), and therefore the EFP exponent f(n) in Eq. (3), for non-integer n. The paper then uses Eqs. (41)-(42) to compute the spatiotemporal profile in Fig. 1, including the astroid-shaped empty region, the critical time Eq. (4), and n-dependent singular exponents near x = ±1 and tau = ±tau_c.","tokens_in":14729,"tokens_out":10653,"duration_ms":102442,"significance":"If the construction is fully established, this is a substantial result: it gives the first analytic emptiness instanton for non-integer polytropic indices, proves the conjectured EFP exponent of Ref. [3] for all gamma > 1, and connects the hydrodynamic limit-shape problem to Dotsenko-Fateev integrals. The paper has genuine strengths: there are no fitted parameters; the Euler-Poisson equation is satisfied by linearity; the boundary conditions and the quadrupole amplitude are checked through explicit hypergeometric manipulations; and the amplitude reproduces the integer-n result exactly. The result would be of interest to the statistical mechanics and cold-atom communities. However, the global spatiotemporal profile claim is currently weakened by an admitted branch ambiguity in the patching of the two integral representations, so the paper needs revision before the full claim can be accepted.","major_comments":[{"comment":"The full spatiotemporal profile is obtained by patching two integral representations whose equivalence is explicitly not guaranteed. The derivation of Eq. (36) and the identity leading to Eq. (43) assume Re lambda = Re bar lambda = v > 0, as stated in Appendix B. The physical profile includes the left half of the astroid, where v < 0, and the density crosses rho = 1 where the text switches between Eq. (36) and Eq. (43). Appendix B closes by stating that outside v > 0 the equivalence between the two representations 'can be violated.' Consequently the boundary condition (22) at x = -1, the asymptotic density (47), and the boundary scaling (56) are not established for the continued solution. The EFP exponent obtained from Eq. (38) on the axis is less exposed, but the global instanton claim requires either a branch-continuation proof or a direct numerical check that both representations agree on the overlap and on the v < 0 part of the physical domain.","section":"Section 5 and Appendix B, Eqs. (36), (43), and Fig. 1"},{"comment":"The manuscript highlights gamma = 2, n = 1/2, as a case studied numerically in Ref. [2], but it contains no quantitative comparison with that numerical solution. Since the new step is analytic continuation in n, a comparison of at least rho(x,0), rho(0,tau), and the boundary shape with the numerical data of Ref. [2] is the natural falsifiable check of the branch choice. Without it, the claim that the construction correctly covers the weakly interacting Bose gas rests on internal consistency alone.","section":"Sections 1 and 5, gamma = 2 case"},{"comment":"The prefactor 2/(1 + e^{2 pi i n})^2 in Eq. (34) is singular at n = 1/2, which is the value corresponding to gamma = 2 and is displayed in Fig. 1. The text explains the factor only for integer n, where it equals 1/2. If the Pochhammer integral vanishes at the same values so that the limit is finite, that cancellation should be shown explicitly; otherwise Eq. (34) is not a valid representation for all n > -1/2 as claimed.","section":"Section 4, Eq. (34)"}],"minor_comments":[{"comment":"The stated expansion near x = 0, |tau| -> tau_c^+ is not the expansion of Eq. (45); the correct leading behavior is (2(|tau| - tau_c)/tau_c)^(n + 1/2), not ((2|tau| - tau_c)/tau_c)^(n + 1/2).","section":"Section 5, Eq. (46)"},{"comment":"The numerical evaluation of Eq. (43) uses a principal-value prescription at q = 0; the implementation should be described briefly so that the profiles in Fig. 1 are reproducible and it is clear how the branch switch between Eq. (36) and Eq. (43) is implemented in practice.","section":"Section 5, Eq. (43)"},{"comment":"The caption should indicate which region of the profile is computed with Eq. (36) and which with Eq. (43), since the switching between representations is part of the construction and is relevant to interpreting the plotted density.","section":"Figure 1 caption"},{"comment":"There are minor typographical errors, such as 'indpendent' in Appendix A and 'densties' in Section 5, which should be corrected.","section":"Appendix A and Section 5"}],"recommendation":"major_revision","confidential_remarks":"The main risk is the branch mismatch that the authors themselves acknowledge at the end of Appendix B. I do not regard this as a fatal defect: the axis amplitude and the EFP exponent appear to be on much firmer ground, and a consistency check of the two representations on the overlap, plus a comparison with Ref. [2] for gamma = 2, could resolve it. The self-citation of Ref. [3] is appropriate because the present derivation is independent and goes beyond that work. I would not recommend rejection if the branch issue is addressed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: the paper does something real—it extends the emptiness instanton from integer n to all n > -1/2 and proves the conjectured EFP exponent f(n) from Ref [3]. The Dotsenko-Fateev/Pochhammer representation is an elegant tool, and the verification by linearity plus explicit hypergeometric asymptotics is convincing for the central result. I would send this to a referee.\n\nWhat is genuinely new: the analytic continuation via the Pochhammer contour, Eq. (34), which satisfies the Euler-Poisson equation by linearity; the explicit check of the boundary condition at |lambda|→infty and the quadrupole amplitude alpha, Eq. (39), matching Ref [3]; and the resulting proof of f(n) for arbitrary polytropic index. The critical time tau_c and the density on the axes x=0 and tau=0 are also clean and explicit. The self-citation to Ref [3] is legitimate; the paper is explicitly building on that construction.\n\nThe soft spot is the branch issue, exactly where the stress-test note lands. Appendix B ends by stating that the two representations used for the profile, Eqs. (36) and (43), can lie on different branches outside Re lambda = Re lambdabar = v > 0, and that their equivalence 'can be violated'. The paper patches them to draw the full spatiotemporal profile, including the rho<1 region and the vicinity of x=-1, but never proves the patched branches agree across rho=1, nor that the branch choice gives the boundary condition at x=-1. So the off-axis density shapes and the n-dependent exponents at x=±1, Eqs. (46)-(47), are not rigorously established. They are probably right—the construction is natural—but the proof as written does not close the argument. The EFP exponent itself is safe: it follows from the mu→1+ amplitude on the axis, which is verified, so the main quantitative claim stands.\n\nOne more gap: no comparison with the numerical solution of Ref [2] for gamma=2 (n=1/2), the most physically accessible interacting case. A quick overlay of the analytic profile against those numerics would be a strong check and might also illuminate the branch question.\n\nWho is this for: people working on large deviations in 1D fluids, emptiness probabilities, and limit shapes. The f(n) proof is worth citing on its own. The paper deserves a serious referee; the revision should either prove the branch consistency or explicitly qualify the full-profile claims as a formal patch with a conjecture for the remaining regions.","headline":"Proves the conjectured emptiness exponent for arbitrary polytropic index via a neat analytic continuation, but the admitted branch mismatch leaves the full spacetime profile and the x=-1 endpoint unproven.","tokens_in":15241,"tokens_out":6650,"would_cite":true,"duration_ms":60217,"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":"This paper proves that the potential $V_n$ defined by a Pochhammer contour integral solves the emptiness instanton problem for every polytropic index $\\gamma > 1$, confirming the conjectured emptiness formation probability exponent $f(n)$…","keywords":["emptiness formation probability","polytropic gas","instanton","hydrodynamics in imaginary time","Pochhammer contour integral","Dotsenko-Fateev integrals","Euler-Poisson equation","large deviations"],"falsifier":"Compute the density profile along a curve in the $(x,\\tau)$ plane that crosses the boundary $\\rho = 1$ away from the axes using both Eq. (36) and Eq. (43); if the two branches disagree there, the off-axis part of the instanton and the claimed exponent at $x = \\pm 1$ are unsupported.","tokens_in":14255,"feed_emoji":"🕳️","tokens_out":9053,"duration_ms":71040,"temperature":0.7,"pith_summary":"The paper addresses the emptiness formation problem: how likely it is that a macroscopically empty interval spontaneously appears in the ground state of a one-dimensional quantum gas with a polytropic equation of state $P \\sim \\rho^{\\gamma}$. For a long interval, the probability is exponentially small, $P_{\\mathrm{EFP}} \\sim e^{-(R^2 \\rho_0 m c_0 / \\hbar) f(n)}$, and is dominated by a classical instanton solution of the hydrodynamic equations in imaginary time. Previous work obtained the instanton only for integer $n$ and conjectured the closed form $f(n) = \\frac{2}{n+1}\\left[\\frac{\\Gamma(n+3/2)}{\\Gamma(n+1)}\\right]^2$ for all $n$. This paper constructs the potential $V_n(\\lambda, \\bar{\\lambda})$ through a Pochhammer contour integral, verifies the Euler-Poisson equation, the endpoint boundary conditions, and the quadrupole asymptotic amplitude, and thereby proves the instanton and the exponent $f(n)$ for all $n > -1/2$, that is, for all $\\gamma > 1$. The result makes the emptiness probability analytically known across a continuous range of polytropic indices, including the weakly interacting Bose gas at $\\gamma = 2$ ($n = 1/2$), which was previously accessible only numerically.","feed_headline":"One contour integral solves emptiness for any polytropic gas","feed_subtitle":"The probability exponent f(n), proven only for integer n before, now holds for all γ > 1 — including weakly interacting Bose gas at γ = 2.","key_machinery":"The load-bearing object is the potential $V_n(\\lambda,\\bar{\\lambda})$ appearing in the hodograph equation $x - w\\tau = \\partial_\\lambda V_n$ and its complex conjugate. For integer $n$ this potential was a finite sum; here it is written as a contour integral, Eq. (34), whose integrand contains $z(z^2+1)^{n-1/2} / [(z-\\lambda)^n(z-\\bar{\\lambda})^n]$ integrated around a Pochhammer contour that winds around the branch points $i$ and $-i$ in opposite senses. That contour is what makes the analytic continuation to non-integer $n$ single-valued. The same integrand satisfies the Euler-Poisson equation (21) for every $z$, so $V_n$ inherits it by linearity; then the boundary condition at the interval endpoints and the quadrupole asymptotics are extracted from the equivalent one-dimensional integrals (36) and (43). From this machinery the paper derives the empty-region boundary via the function $g(v)$ in Eq. (50) and the explicit singularity exponents of the density at the axes.","core_discovery":"The central claim is that Eq. (34), the potential $V_n(\\lambda, \\bar{\\lambda})$ defined by a Pochhammer contour integral together with the prefactor $2/(1+e^{2\\pi i n})^2$, is the correct analytic continuation of the emptiness instanton of Ref. [3] to non-integer $n$, and is valid for every real $n > -1/2$ (equivalently $\\gamma > 1$). The argument proceeds in three steps: the integrand obeys the Euler-Poisson equation (21) at every point $z$, so $V_n$ satisfies it by linearity; the collapsed-contour representation (36) gives $\\partial_\\lambda V_n \\to \\pm 1$ as $|\\lambda| \\to \\infty$, matching the required singularities at the interval endpoints; and the evaluation (38) on the symmetry line $\\lambda = -\\bar{\\lambda} = i\\mu$ reproduces the inverse-square-root quadrupole asymptotics (24) with amplitude $\\alpha = \\frac{1}{2}(2n+1)\\left[\\frac{\\Gamma(n+3/2)}{\\Gamma(3/2)\\Gamma(n+1)}\\right]^2$, which feeds into $f(n)$ through Eq. (25). Thus the emptiness formation probability exponent (3), the critical time (4), and the spatiotemporal instanton profile (41)-(42) hold for arbitrary polytropic index $\\gamma > 1$, not just the integer values treated before.","pith_inferences":["If the branch-patching between Eqs. (36) and (43) is benign on the overlap, the same Pochhammer-contour construction should produce exact instanton solutions for other large-deviation observables in polytropic gases, such as full counting statistics, by changing only the boundary conditions.","The Dotsenko-Fateev form of the potential hints that emptiness formation in an interacting gas may be governed by a conformal field theory with a $\\gamma$-dependent central charge; if so, exact microscopic universality would extend beyond free fermions.","The coexistence of a universal $3/2$ exponent at $\\tau = \\pm\\tau_c$ with $n$-dependent exponents at $x = \\pm 1$ suggests that fluctuations around the emptiness boundary may exhibit an $n$-dependent dynamical exponent, a question the authors explicitly leave open and that could be tested by instanton fluctuation calculations.","Formally continuing to $n = -1$ (Chaplygin gas) makes $\\tau_c = f(n) = 0$; the comment in the paper that this continuation solves a different problem can be tested by checking whether the analytic potential (34) at $n = -1$ yields a real, positive hydrodynamic density profile at all."],"forward_implications":["The emptiness formation probability exponent is now proven in closed form for all $\\gamma > 1$: $f(n) = \\frac{2}{n+1}\\left[\\frac{\\Gamma(n+3/2)}{\\Gamma(n+1)}\\right]^2$, so the exponential suppression of emptiness is known for every polytropic index.","The analytic continuation includes the weakly interacting Bose gas case $\\gamma = 2$, $n = 1/2$, for which the instanton profile was previously obtained only by numerical solution of the hydrodynamic equations.","The spacetime shape of the empty region is astroid-like for all $n$; near $x=0$, $\\tau = \\pm\\tau_c$ the boundary always scales as $|x| \\sim (|\\tau|-\\tau_c)^{3/2}$, while near $x = \\pm 1$, $\\tau = 0$ the scaling exponent is $n$-dependent, $|\\tau| \\sim (1-|x|)^{(2n+3)/(2n+2)}$.","Explicit density profiles are available on the symmetry axes: $\\rho(0,\\tau) = [1 - (\\tau_c/\\tau)^2]^{n+1/2}$ and $\\rho(x,0) \\sim (A_n/(|x|-1))^{(2n+1)/(2n+2)}$ near the endpoint.","The closed-form hydrodynamic profile provides concrete predictions that can be checked by direct numerical simulation of the imaginary-time hydrodynamic equations for any $\\gamma > 1$."],"supporting_citations":[{"why":"Supplies the integer-n emptiness instanton solution and the conjectured formula for f(n) that this paper extends to all n > -1/2.","marker":"[3]"},{"why":"Introduces the imaginary-time hydrodynamic method and gives the free-fermion (n=0) astroid solution on which the integral construction is built.","marker":"[24]"},{"why":"Provides the numerical hydrodynamic solution for gamma=2 (n=1/2), the case the analytic continuation now covers exactly.","marker":"[2]"},{"why":"Classical treatment of the Pochhammer contour and beta-function continuation that justifies the single-valued contour for non-integer n.","marker":"[28]"},{"why":"Supplies the hypergeometric transformation used to extract the inverse-square-root asymptotics and the quadrupole amplitude alpha.","marker":"[29]"},{"why":"Identifies the integrals as Dotsenko-Fateev integrals, connecting the instanton potential to correlation functions in conformal field theory.","marker":"[30]"}],"fun_headline_variants":["Emptiness instanton now solved for all polytropic gases","Non-integer polytropic index tamed by contour integral","Any gamma > 1: emptiness probability from one integral","Analytic continuation unlocks emptiness for any polytropic gas"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The two integral representations used for the plotted profile must give the same branch of the multi-valued potential in the overlap region where the density crosses $\\rho = 1$; the paper asserts but does not prove their equivalence there.","fun_headline_variants_meta":{"raw":{"variants":["Emptiness instanton now solved for all polytropic gases","Non-integer polytropic index tamed by contour integral","Any gamma > 1: emptiness probability from one integral","Analytic continuation unlocks emptiness for any polytropic gas"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000318,"raw_usage":{"total_tokens":1821,"prompt_tokens":997,"completion_tokens":824,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":613,"completion_tokens_details":{"reasoning_tokens":754}},"tokens_in":613,"tokens_out":824,"duration_ms":7110,"temperature":1.0,"reasoning_tokens":754,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T14:40:39.963428+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the density profile along a curve in the $(x,\\tau)$ plane that crosses the boundary $\\rho = 1$ away from the axes using both Eq. (36) and Eq. (43); if the two branches disagree there, the off-axis part of the instanton and the claimed exponent at $x = \\pm 1$ are unsupported.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the integer-n emptiness instanton solution and the conjectured formula for f(n) that this paper extends to all n > -1/2."},{"cited_title":"Hydrodynamics of correlated systems. Emptiness Formation Probability and Random Matrices","cited_arxiv_id":"cond-mat/0504307","evidence_quote":"Introduces the imaginary-time hydrodynamic method and gives the free-fermion (n=0) astroid solution on which the integral construction is built."},{"cited_title":"Yeh and A","cited_arxiv_id":null,"evidence_quote":"Provides the numerical hydrodynamic solution for gamma=2 (n=1/2), the case the analytic continuation now covers exactly."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Classical treatment of the Pochhammer contour and beta-function continuation that justifies the single-valued contour for non-integer n."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the hypergeometric transformation used to extract the inverse-square-root asymptotics and the quadrupole amplitude alpha."}],"review_version":1}