{"id":"c86b028a-d094-4abd-b512-a4744c17fff7","arxiv_id":"1908.01824","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"In the trapped one-dimensional Bose-Hubbard model, bipartite entanglement entropy acts as an order parameter for the formation of a Mott-insulating domain and for the emergence of a central superfluid domain.","lead":"This paper uses numerical simulations to show that entanglement entropy can detect two local quantum phase transitions in a trapped cloud of bosons on a one-dimensional lattice. The result suggests a practical way to locate these transitions in ultracold-atom experiments that can already measure entanglement.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Lower-transition extrapolation relies on a possibly underdetermined quartic fit in 1/N*, so the claimed LDA agreement is not yet quantitatively established.","rationale":"The reader's verdict is CONDITIONAL, and my read supports keeping it: the paper contains genuine numerical evidence (DMRG data, LDA comparisons, scaling collapses) that entanglement entropies track both local transitions, and the energy-derivative kink/discontinuity is convincingly demonstrated against LDA. The concern I identify is a sharper version of the reader's weakest assumption: the lower-transition thermodynamic-limit density is extracted from very few small-system crossing points via a quartic fit that, as described, has more parameters than data points. This makes the 'agreement with LDA' for the lower transition fragile. It does not warrant rejection, because the discrepancy is fixable by reanalysis and by adding error bars, and because the upper-transition extrapolation and the LDA reference values are less severely underdetermined. The verdict remains CONDITIONAL: the central claim about quantitative extraction of rho_l^c needs either a corrected, well-constrained fit or an independent large-system extraction of rho_l^c from the scaling collapse.","tokens_in":15170,"tokens_out":7016,"duration_ms":83049,"concrete_test":"Re-analyze the Fig. 6 data by extracting all available crossing points rho*(N1,N2) for pairs with 0<|N1-N2|<=4 and fit them with quadratic, cubic, and quartic polynomials in 1/N*, reporting the extrapolated rho_c with bootstrap confidence intervals. If the number of crossing points is <=3, a quartic fit is underdetermined; show instead e.g. a quadratic fit and a leave-one-out stability check. The lower-transition claim is settled only if the extrapolated rho_c is stable (within, say, 0.005 of the LDA value) across fit orders and across removal of the smallest-N point.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central quantitative claim is that entanglement entropies locate both local transitions at the LDA critical densities in the thermodynamic limit. For the lower transition, the only systematic extraction from entanglement data is the crossing-point analysis of S2 in Fig. 6. The text states that crossing points are taken between curves with N and N+4 particles, and for the values listed (U=5: N=16,20,24,28; U=6.25,8,10: N=8,12,16,20) that yields at most three crossing points per panel. The extrapolating function is a quartic polynomial in 1/N*, i.e., five parameters. Three data points cannot determine five parameters; the reported extrapolated rho_c is consequently not a well-defined result of the data, and the apparent agreement with LDA in the insets cannot be evaluated. The upper-transition linear fits in 1/N and the LDA reference extrapolations (cubic in 1/L0) are also empirical, but at least they are not formally underdetermined; the lower-transition quartic is the weakest link in the chain connecting finite-system entanglement data to thermodynamic-limit critical densities.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies the ground-state bipartite entanglement entropies (von Neumann and second Rényi) in the one-dimensional Bose-Hubbard model in a harmonic trap. It identifies two local quantum phase transitions upon increasing the characteristic density ρ=N/R: the formation of an n=1 Mott-insulating domain at the trap center (lower transition) and the emergence of an n>1 superfluid domain at the center of that Mott domain (upper transition). Using DMRG, the authors present evidence that the second derivative of the energy density with respect to ρ has a kink at the lower transition and a discontinuity at the upper transition, and they propose that entanglement entropies serve as order parameters. They extract critical densities from crossing points and midpoints of entropy curves for finite systems, extrapolate them to the thermodynamic limit, and compare with local-density-approximation (LDA) predictions.","tokens_in":15460,"tokens_out":4973,"duration_ms":50161,"significance":"If the claims hold, the paper provides a practical, experimentally accessible route to locate both local transitions using the second-order Rényi entropy, which has been measured in cold-atom experiments. The DMRG calculations appear carefully converged (bond dimension 3200, truncation 1e-12), and the homogeneous-limit entanglement scalings in Fig. 4 are checked against perturbative results. The data collapse in Fig. 5 and the linear extrapolations for the upper transition in Fig. 7 are suggestive and constitute a useful scaling framework for trapped systems. However, the quantitative lower-transition extrapolation is not well determined, so the central agreement claim is not yet established.","major_comments":[{"comment":"The extrapolation of the lower-transition crossing points ρ× via a quartic polynomial in 1/N* is underdetermined. For U=5, the curves for N=16 and 20 do not cross (as stated in the text), leaving only two crossing points (20/24 and 24/28). For U=6.25, 8, and 10, the four particle numbers N=8,12,16,20 give only three crossing points. A quartic polynomial has five free parameters, so these fits are not uniquely defined, and the extrapolated values shown in the insets are not well-determined results of the data. The apparent agreement with LDA therefore cannot be evaluated. The authors should either use more system sizes, reduce the polynomial order to a determined fit, or provide a constrained extrapolation with an error estimate.","section":"Sec. III A, Fig. 6 and inset"},{"comment":"The scaling collapse in the insets of Fig. 5 uses the LDA value ρ_l^c to define the rescaled variable \\tildeρ=(ρ−ρ_l^c)N. Thus the collapse cannot independently confirm the LDA prediction; it only demonstrates consistency with it. To support the claim that entanglement measures determine the transition point, the collapse should be used with ρ_l^c treated as an adjustable parameter or otherwise validated by a method that does not use the LDA value as input.","section":"Sec. III A, Eq. (9) and Fig. 5 insets"},{"comment":"The abstract states that the second derivative of the total energy is continuous with a kink at the lower transition and discontinuous at the upper transition. The evidence in Fig. 3 consists of finite-size DMRG curves for two values of R plus LDA curves. No finite-size scaling or extrapolation of \\bar E'' to the thermodynamic limit is presented, and LDA itself is an approximation whose accuracy is not quantified at these parameters. As written, this asserts a stronger result than the numerics establish. Please soften the claim or add a quantitative finite-size analysis of the derivative.","section":"Sec. III, Fig. 3 and abstract"}],"minor_comments":[{"comment":"Equation (9) writes SvN = F((ρ−ρ_l^c)N), but the insets of Fig. 5 use \\tildeρ=(ρ−ρ_l^c)N; please make the notation consistent and clarify the collapse criterion.","section":"Eq. (9) and Fig. 5"},{"comment":"There are several typos: 'superfulid' in Section I, 'in presented' in Section IV should be 'is presented', and 'Renyi' should be 'Rényi'.","section":"Throughout"},{"comment":"'Quartic order polynomial' should be 'quartic polynomial'; also, please state how many crossing points enter each fit and why a quartic form was chosen.","section":"Fig. 6 caption"},{"comment":"The LDA reference extrapolations are shown without error bars; please state the number of L0 values used in each cubic fit and report the resulting uncertainty on ρ_l^c and ρ_u^c.","section":"Fig. 2 and Sec. II"}],"recommendation":"major_revision","confidential_remarks":"The underdetermined quartic fit in Fig. 6 is a genuine technical flaw that affects the paper's central quantitative claim that entanglement entropies locate the lower transition at the LDA density in the thermodynamic limit. I would ask the authors to redo the lower-transition extraction with a determined fitting scheme, add more system sizes if possible, and report uncertainties. The rest of the paper is largely sound and the upper-transition analysis is more convincing."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a genuinely useful extension of the authors' fermionic work to bosons, and the qualitative picture is almost certainly right. The paper's real contributions are (i) showing that the second derivative of the energy has a kink at the lower local transition and a discontinuity at the upper one, and (ii) demonstrating that the von Neumann and Renyi entropies behave like order parameters, with a data collapse for the lower transition and a clean finite-size trend for the upper one. The DMRG is carefully converged, and the comparison with LDA is the right check.\n\nNow the soft spots. The stress-test note is correct. The crossing-point extrapolation for the lower transition (Fig. 6) uses three crossing points (from N and N+4 curves, so at most 3 points) fitted to a quartic polynomial in 1/N* — that is five parameters. The fit is formally underdetermined, so the extrapolated rho_c is not a well-defined result of the data, and the apparent agreement with LDA in the insets cannot be evaluated. This is the weakest link connecting finite-system entropies to thermodynamic-limit critical densities. The fix is straightforward: use a lower-order fit or more crossing points, and report error bars.\n\nThe other extrapolations are also empirical but at least not underdetermined: cubic in 1/L0 for the LDA reference values and linear in 1/N for the upper-transition midpoint. The scaling collapse in Fig. 5 uses the LDA value rho_lc as input, so it cannot independently confirm rho_lc; the authors don't overclaim this, but the circularity should be acknowledged more explicitly. The claim about the thermodynamic-limit nonanalytic behavior of E'' is inferred from finite-size DMRG and LDA curves; it is plausible but not proven. Minor: no quantitative comparison with local-property methods is given, though the qualitative statement that entanglement is more accurate is fine.\n\nWho is this for? People working on trapped ultracold gases and entanglement measures. It deserves a serious referee; the central idea is solid and the paper is mostly careful. But the underdetermined fit should be fixed or downplayed before I would trust the lower-transition quantitative claim.","headline":"Entanglement order parameters for local transitions in trapped Bose-Hubbard is a genuinely useful idea with careful DMRG, but the lower-transition extrapolation rests on an underdetermined quartic fit and needs a fix.","tokens_in":15960,"tokens_out":1945,"would_cite":true,"duration_ms":20371,"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":"Bipartite entanglement entropies are order parameters for both local quantum phase transitions in trapped one-dimensional lattice bosons, and finite-size scaling recovers the thermodynamic-limit critical densities predicted by the local…","keywords":["Bose-Hubbard model","local quantum phase transition","entanglement entropy","Renyi entanglement entropy","harmonic trap","local density approximation","density-matrix renormalization group","Mott insulator"],"falsifier":"Compute the transition densities for the same trap and interaction parameters directly from a trapped-system observable that does not rely on the local density approximation—for example, the kink or jump in the second energy derivative or the peak of the local compressibility—and extrapolate those to the thermodynamic limit; if the resulting $\\rho_l^c$ and $\\rho_u^c$ disagree with the entanglement-extrapolated values beyond the numerical uncertainty of the extrapolations, the claim that entanglement entropies order these transitions would be falsified.","tokens_in":14952,"feed_emoji":"⚛️","tokens_out":11566,"duration_ms":106194,"temperature":0.7,"pith_summary":"The paper asks whether entanglement, rather than only local densities, can detect the two local quantum phase transitions that occur when bosons in a one-dimensional harmonic trap are compressed by increasing the characteristic density $\\rho=N/R$. It argues that the second derivative of the total energy is continuous with a kink at the lower transition (formation of an $n=1$ Mott-insulating core) and discontinuous at the upper transition (emergence of an $n>1$ superfluid core). It then shows that bipartite von Neumann and second Renyi entanglement entropies act as order parameters: the entropy decreases smoothly toward the homogeneous Mott value at the lower transition and rises sharply at the upper transition. Using density-matrix renormalization group data, the paper shows that crossing-point and midpoint extrapolations from small systems recover the local-density-approximation transition densities in the thermodynamic limit. This matters because the second Renyi entropy is measurable in current cold-atom experiments, so the transitions could be located directly from entanglement measurements.","feed_headline":"Entanglement entropy pins down both trapped-boson transitions","feed_subtitle":"Entropy signals the Mott core and the superfluid core, matching LDA densities in the thermodynamic limit.","key_machinery":"The load-bearing object is the bipartite entanglement entropy—the von Neumann entropy $S_{\\mathrm{vN}}=-\\mathrm{Tr}\\,\\rho_A\\ln\\rho_A$ and the Renyi entropies $S_\\alpha=(1-\\alpha)^{-1}\\ln\\mathrm{Tr}\\,\\rho_A^\\alpha$, with $S_2$ the experimentally measurable case—computed from the reduced density matrix of half the trapped chain by density-matrix renormalization group. Around the lower transition, a universal scaling function $S_{\\mathrm{vN}}(U)=F([\\rho-\\rho_l^c(U)]N)$ supplies data collapse that locates the critical density; around the upper transition, the midpoint of the sharp entropy rise, defined by Eq. (10) as $S_{\\mathrm{vN}}(\\bar\\rho)=(S_{\\mathrm{vN}}^{(0)}+S_{\\mathrm{vN}}^{\\max})/2$, serves as the finite-size estimator. The local density approximation, Eq. (3), provides the reference thermodynamic-limit transition densities that these extrapolations are checked against.","core_discovery":"The central claim is that nonlocal information measures carry the same ordering information as the energy nonanalyticities in an inhomogeneous Bose-Hubbard system. In the trapped one-dimensional Bose-Hubbard model, the paper finds a universal scaling form $S_{\\mathrm{vN}}(U)=F([\\rho-\\rho_l^c(U)]N)$ for the lower transition, with data collapse for different particle numbers, and a sharp rise in $S_{\\mathrm{vN}}$ at the upper transition whose midpoint approaches the local-density-approximation value $\\rho_u^c$ as $1/N$ decreases. Extrapolating those features with polynomial fits in inverse system size reproduces the LDA critical densities for the interaction strengths studied, and the same analysis works for the experimentally accessible second Renyi entropy in systems with a few tens of particles. The implication is that entanglement entropy is a quantitative order parameter for local quantum phase transitions in trapped bosons, not merely a qualitative indicator.","pith_inferences":["If the low-temperature expectation stated in the paper holds, the same $S_2$ crossing and midpoint analysis could be applied directly to finite-temperature experimental snapshots, turning the Renyi entropy into a practical transition locator; the paper only states this as an expectation, not a demonstrated result.","The pattern of nonanalyticities—a kink at the lower transition and a jump at the upper one—suggests the two local transitions may belong to different classes of critical behavior, but the paper does not make that classification.","Because the local-density-approximation reference values are themselves extrapolated from homogeneous finite-size data, an independent determination of $\\rho_l^c$ and $\\rho_u^c$ from trapped-system observables alone would test whether the agreement with LDA is robust or partly an artifact of shared finite-size assumptions."],"forward_implications":["The second derivative of the total energy distinguishes the two local transitions: a kink at the lower transition and a jump at the upper one, so energy measurements alone can already separate them.","Both the von Neumann entropy and the second Renyi entropy $S_2$ work as order parameters, and $S_2$ is the quantity measured in current ultracold-atom experiments.","For sufficiently large $U$, the critical characteristic densities for both transitions can be extracted from systems with only tens of particles, within reach of optical-lattice experiments.","The universal scaling collapse of $S_{\\mathrm{vN}}$ versus $(\\rho-\\rho_l^c)N$ implies a system-size-independent order-parameter shape at fixed $U$, paralleling the fermionic case.","As $U$ approaches the homogeneous critical value $U_c^{n=1}=3.28$ from above, the lower and upper transitions are no longer well separated, and larger systems are needed both theoretically and experimentally."],"supporting_citations":[{"why":"Companion study of trapped spinless fermions; supplies the universal scaling and entanglement-entropy order-parameter approach that this paper extends to bosons.","marker":"[43]"},{"why":"State diagrams for trapped bosons; identifies the lower and upper transitions from local observables and gives the baseline critical densities.","marker":"[33]"},{"why":"Optical-lattice experiment that measured the second-order Renyi entanglement entropy; motivates the small-system extrapolation and experimental relevance.","marker":"[17]"},{"why":"Perturbative entanglement-spectrum analysis of gapped systems; provides the analytical homogeneous Mott entropies used as reference values.","marker":"[57]"},{"why":"Homogeneous Bose-Hubbard data establishing the critical interaction $U_c^{n=1}=3.28$ and the ground-state energies used in the local-density-approximation state diagram.","marker":"[72, 73]"}],"fun_headline_variants":["Entanglement entropy orders trapped-boson phase transitions","Entropy measures pinpoint trapped-boson transitions","Renyi entropy detects both trapped-boson critical points","Quantum information unveils local transitions in trapped bosons"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument assumes the finite-size corrections to the transition densities follow the specific polynomial forms chosen for the extrapolations, and that the local density approximation accurately describes the trapped system at the interaction strengths studied.","fun_headline_variants_meta":{"raw":{"variants":["Entanglement entropy orders trapped-boson phase transitions","Entropy measures pinpoint trapped-boson transitions","Renyi entropy detects both trapped-boson critical points","Quantum information unveils local transitions in trapped bosons"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000593,"raw_usage":{"total_tokens":2780,"prompt_tokens":947,"completion_tokens":1833,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":563,"completion_tokens_details":{"reasoning_tokens":1769}},"tokens_in":563,"tokens_out":1833,"duration_ms":13835,"temperature":1.0,"reasoning_tokens":1769,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:01:33.289270+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the transition densities for the same trap and interaction parameters directly from a trapped-system observable that does not rely on the local density approximation—for example, the kink or jump in the second energy derivative or the peak of the local compressibility—and extrapolate those to the thermodynamic limit; if the resulting $\\rho_l^c$ and $\\rho_u^c$ disagree with the entanglement-extrapolated values beyond the numerical uncertainty of the extrapolations, the claim that entanglement entropies order these transitions would be falsified.","supporting_citations":[{"cited_title":"Zhang, L","cited_arxiv_id":null,"evidence_quote":"Companion study of trapped spinless fermions; supplies the universal scaling and entanglement-entropy order-parameter approach that this paper extends to bosons."},{"cited_title":"Rigol, G","cited_arxiv_id":null,"evidence_quote":"State diagrams for trapped bosons; identifies the lower and upper transitions from local observables and gives the baseline critical densities."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Perturbative entanglement-spectrum analysis of gapped systems; provides the analytical homogeneous Mott entropies used as reference values."}],"review_version":1}