{"id":"80a91cd1-4bea-4080-a79e-6329294f6c26","arxiv_id":"2502.09448","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"In a trapped 1D lattice, a single repulsive impurity in a Mott-insulator bath forms a correlated counterflow state with a combined unity-filling profile and slowly decaying anti-pair correlations.","lead":"A single impurity atom in a 1D optical lattice with a repulsively interacting boson bath can, at intermediate interaction strengths, form a 'counterflow' state where the impurity and a bath hole move in opposite directions through a Mott-insulator region. The paper maps this behavior with numerical simulations and a simple particle-hole model, suggesting counterflow order can appear in highly imbalanced atomic mixtures, which is relevant to ultracold atom experiments.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The counterflow phase claim rests on a single-particle cosine fit to CAP in finite DMRG samples; the analytic model inserts this profile by construction, so thermodynamic long-range order is not established.","rationale":"The reader's weakest_assumption identifies the core issue, and my independent read of Eqs. (3)-(7) confirms it. The paper's own criterion for counterflow order is algebraic decay of CAP, but the presented CAP is a finite-box cosine with a zero at a finite distance. The analytic impurity-hole model recovers nb + nI = 1 and CAP = alpha0*alpha_i by construction from the single-particle ansatz, so it cannot serve as independent evidence of a phase. The DMRG and ED cross-checks are valuable and the qualitative profile observations are likely robust, but they are finite-size observations and do not establish a thermodynamic transition. The proposed scaling test would directly settle whether the correlation tail becomes algebraic. The existing CONDITIONAL verdict therefore remains appropriate, with the condition being the finite-size scaling of CAP.","tokens_in":12895,"tokens_out":10334,"duration_ms":103565,"concrete_test":"At fixed rho_tilde = 3.58 and t/Ubb = 0.1, UbI/Ubb = 0.6, run DMRG for Nb = 25, 40, and 60 (with Vho = 0.0205t, 0.008t, and 0.0036t respectively, and M = Nb + 20), then plot CAP versus xi on log-log and semi-log scales. If the counterflow is thermodynamic, CAP at a fixed rescaled distance should not be controlled by a zero at xi = l_CF/2 and the tail should approach an algebraic decay as Nb grows. If instead the curves collapse onto n0 cos(pi (xi - <x_I>)/l_CF) with the same l_CF and no power-law tail emerges, region v is a finite-size trap coherence and the phase claim should be downgraded.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires long-range counterflow order, which the paper itself defines (Eq. 3) as slow algebraic decay of CAP. The numerical evidence in Fig. 3(c) is instead fitted by Eq. (7), CAP(i) = n0 cos(i pi / M_CF), a single-particle finite-box wavefunction. A cosine is not algebraic and vanishes at the box edge, so it cannot by itself certify thermodynamic order. The analytical model (Eqs. 5-6) assumes |Psi_CF> = sum_i alpha_i |i_IH>, from which nb + nI = 1 and CAP = alpha0*alpha_i follow exactly; the cos^2 profile of Eq. (4) is inserted as a fit to the numerical profile, not derived. Thus the model does not provide independent confirmation of order. Moreover, because the calculations are performed at fixed characteristic density rho_tilde = 3.58, increasing Nb at fixed rho_tilde increases M_CF in lattice units but keeps the rescaled box width l_CF approximately constant. CAP then remains the same finite-box cosine in rescaled distance rather than developing an algebraic tail. The DMRG benchmarks with Nb = 25 and the ED runs collapse onto the same finite-size curves and do not supply the missing scaling. The data are therefore consistent with a trap-induced single-particle coherence, not with a thermodynamic counterflow phase.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a single mobile impurity immersed in a one-dimensional harmonically confined bosonic lattice bath described by a two-component Bose-Hubbard model. Using DMRG for Nb=40 bath atoms (with ED and smaller-Nb benchmarks), it maps the phase diagram in t/Ubb and UbI/Ubb. For baths with a central Mott-insulator domain (t/Ubb < 0.155) and intermediate impurity-bath repulsion (roughly 0.4 < UbI/Ubb < 1), it identifies a new 'counterflow phase' in which the sum of bath and impurity densities forms a unity-filling domain, the impurity profile takes a cos^2 shape, the anti-pair correlator CAP decays slowly, and the polaron residue drops sharply. The authors propose an impurity-hole pair model to explain these features and discuss experimental probes.","tokens_in":13286,"tokens_out":3820,"duration_ms":41412,"significance":"If the counterflow phase is genuine, the result would be significant: it would extend counterflow/Mott-insulator physics to highly imbalanced mixtures and demonstrate that a single impurity can form a delocalized correlated structure with the bath over the trap scale, with a sudden orthogonality. The numerical work is careful: the DMRG parameters and convergence checks are described, the constant-characteristic-density scaling is tested with Nb=25 and Nb=40, and ED benchmarks support the raw density profiles. The bath-only characterization and the phase diagram based on average positions and cloud size are useful. However, the central claim of long-range counterflow order currently rests on a finite-size cosine fit and a circular analytical model, so the significance cannot be fully assessed until the thermodynamic-limit evidence is supplied.","major_comments":[{"comment":"The paper defines long-range counterflow order as a slow algebraic decay of CAP, yet the numerical evidence is fitted by Eq. (7), CAP(i) = n0^I cos(i pi/M_CF), which is a single-particle finite-box function, not algebraic, and it vanishes at the edge of the box. The DMRG results are presented for a single system size in the main text (Nb=40, M=60); the supplemental benchmarks with Nb=25 and ED with Nb=8 collapse onto essentially the same rescaled curves, but that is precisely what a finite-size single-particle coherence would do when the characteristic density is held fixed. No finite-size scaling is provided, so the data are consistent with trap-induced coherence rather than thermodynamic counterflow order. Please provide a scaling analysis at fixed rho_tilde (e.g., several Nb with corresponding Vho and M), and show whether the decay of CAP becomes algebraic in the thermodynamic limit, or at least whether the correlation length in lattice units grows linearly with the system size while the vanishing at the box edge moves out.","section":"Counterflow phase, Eq. (3) and Fig. 3(c)"},{"comment":"The analytical model is not an independent confirmation of the counterflow phase. The ansatz |Psi_CF> = sum_i alpha_i |i_IH> enforces by construction nb(i) = 1 - |alpha_i|^2 and CAP(i) = alpha_0^* alpha_i, as derived in the supplemental material. Then Eq. (4), the cos^2 impurity profile, is taken directly from the numerical fit, and Eq. (7) follows immediately. Thus the model merely restates the numerical observations; it does not explain why the cosine profile is energetically selected or why the combined system forms an insulator. To support the central claim, the model should derive alpha_i from the Hamiltonian, for example by a variational minimization over the coefficients, and show that the cosine-like delocalized state is the ground-state solution rather than an input.","section":"Counterflow phase, Eqs. (5)-(6) and supplemental material 'Impurity-hole model'"},{"comment":"The claim of an orthogonality catastrophe in the thermodynamic limit is an extrapolation from finite-size data. The paper states that the residue 'does not show a fully discontinuous' drop and that the discontinuity 'becomes more abrupt with a larger number of particles,' but the evidence shown comprises only Nb=25 and Nb=40 DMRG and Nb=8 ED, all at the same rho_tilde. No scaling of the residue with Nb (e.g., Z at the transition vs. 1/Nb) is presented, so the thermodynamic-limit inference is not supported. Either provide such scaling or phrase the conclusion more cautiously as a sharp finite-size drop.","section":"Polaron properties, Eq. (8) and Fig. 4(c)"}],"minor_comments":[{"comment":"The sentence 'In contrast, in baths with MI domains, the residue remains constant and equal to one for small UbI [region ii]' refers to the wrong phase label: the miscible low-UbI region in MI baths is region iii, while region ii is the phase-separated region for SF baths. Please correct the label.","section":"Polaron properties, paragraph after Eq. (9)"},{"comment":"The legends contain the typo 'DRMG' in several places; it should read 'DMRG'.","section":"Supplemental material, Figs. 5 and 6"},{"comment":"The correlation length symbol used in 'decays as exp(-|xi|/xi_e)' is not defined; please introduce xi_e or use a different notation to avoid confusion with the rescaled coordinate xi.","section":"Counterflow phase, text below Eq. (3)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a letter that reports an interesting and potentially important phase, but the evidence for thermodynamic counterflow order is currently incomplete. The finite-size scaling required to distinguish a genuine long-range-ordered phase from a trap-induced single-particle coherence is well within the authors' numerical capabilities, as they already have a DMRG setup with multiple system sizes. The analytical model is elegant but circular in its present form; it should be reframed as a variational ansatz whose coefficients are optimized, or the authors should clearly state that it is a fitting model rather than a derivation. The paper's scope (a PRL-style letter) is appropriate for the topic, but the central claim needs stronger support before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis paper is worth reading. It maps a single impurity in a harmonically trapped 1D Bose-Hubbard bath with full back-action and finds a new region at intermediate UbI/Ubb where the combined bath+impurity profile forms a unity-filling plateau even though neither species does individually. That observation is new and likely robust, and the numerics are careful: DMRG with convergence checks, ED cross-checks, and consistent results for N=25 and N=40. The polaron residue and energy data across the phase diagram are also useful.\n\nThe soft spot is the central claim. The paper defines counterflow order as algebraic decay of CAP (Eq. 3), but the evidence in Fig. 3(c) is a fit to a cosine, CAP(i) = n0 cos(i pi / M_CF), which is a finite-box single-particle wavefunction. A cosine is not algebraic, and no finite-size scaling is provided. Because the system is run at fixed characteristic density, increasing Nb at fixed rho_tilde leaves the rescaled box width l_CF essentially unchanged, so the same cosine persists; the N=25 benchmarks collapse onto the same curves and do not supply the missing scaling. The analytic model does not rescue this: the ansatz |Psi_CF> = sum_i alpha_i |i_IH> builds in nb+nI=1 and CAP = alpha0*alpha_i, and the cos^2 profile is inserted as a fit. The supplemental candidly notes the model does not conserve particle number, which further limits its use for a thermodynamic phase assignment. So the \"counterflow phase\" is currently a finite-size/trap-induced coherence, not the long-range ordered state claimed.\n\nThe observed profiles and phase boundaries are still valuable; the issue is the interpretation. This is fixable with extra analysis—finite-size scaling of CAP, algebraic-decay diagnostics, or a reframed claim about a correlated finite-size regime.\n\nWho it's for: the cold-atom polaron community and people working on trapped lattice mixtures. A serious referee should engage; the paper should not be desk-rejected. But it needs major revision before acceptance, and the counterflow terminology should be conditional on the scaling.","headline":"A carefully computed phase diagram with a genuinely new combined-unity-filling domain, but the counterflow phase is not backed by its own algebraic-decay criterion—the CAP evidence is a finite-size cosine and the analytic model assumes the conclusion.","tokens_in":13745,"tokens_out":3291,"would_cite":false,"duration_ms":33467,"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":"A single impurity in a trapped lattice can pair with a bath hole to form a delocalized counterflow phase of unity filling.","keywords":["Bose polaron","Bose-Hubbard model","counterflow order","Mott insulator","optical lattice","harmonic confinement","anti-pair correlator","impurity-hole pair"],"falsifier":"For fixed characteristic density $\\tilde{\\rho}_b=3.58$, compute $C_{\\mathrm{AP}}$ at $t/U_{bb}=0.1$ and $U_{bI}/U_{bb}=0.6$ for increasing system sizes, for example $N_b=25,40,80,160$ with $M=N_b+20$. If the fitted cosine's zero crossing or decay length fails to converge, or if $C_{\\mathrm{AP}}$ crosses over to exponential decay beyond a finite $x_i$, then the slow decay is a finite-size effect and the counterflow phase as defined collapses.","tokens_in":12725,"feed_emoji":"🔄","tokens_out":12349,"duration_ms":150236,"temperature":0.7,"pith_summary":"This paper studies a single mobile impurity repulsively coupled to a bosonic bath in a one-dimensional optical lattice with a harmonic trap, and it identifies a new phase in the bath-impurity phase diagram. When the bath develops a central Mott-insulator domain ($t/U_{bb}$ below about 0.155) and the impurity-bath repulsion is intermediate (roughly $0.4 < U_{bI}/U_{bb} < 1$), the impurity no longer behaves as a localized polaron. Instead, the impurity and a bath hole form a delocalized bound structure: the combined density $n_b+n_I=1$ across the central region, the impurity profile becomes the ground-state density of a particle in an infinite square well, and the anti-pair correlator decays slowly, signaling long-range counterflow order. The paper argues this is an unconventional counterflow phase, distinct from the miscible and phase-separated regimes, and explains it with an impurity-hole wavefunction ansatz. If correct, this would show that counterflow order is not limited to balanced mixtures and offers a sharp probe of the trapped superfluid-to-Mott-insulator transition.","feed_headline":"Impurity and bath hole form a delocalized counterflow phase","feed_subtitle":"An impurity and a bath hole lock into a unity-filling insulator, and the polaron suddenly loses its overlap.","key_machinery":"The load-bearing object is the impurity-hole state $|i_{\\mathrm{IH}}\\rangle = \\hat{a}^\\dagger_{i,I}\\hat{a}_{i,b}|\\mathrm{MI}\\rangle$, where $|\\mathrm{MI}\\rangle$ is a unit-filled Mott-insulator state of the bath, and the variational wavefunction $|\\Psi_{\\mathrm{CF}}\\rangle = \\sum_i \\alpha_i |i_{\\mathrm{IH}}\\rangle$. With this ansatz the bath density automatically becomes $n_b(i)=1-|\\alpha_i|^2$, so $n_b+n_I=1$, and the anti-pair correlator becomes $C_{\\mathrm{AP}}(i)=\\alpha_0^*\\alpha_i$. Imposing the observed $\\cos^2$ impurity profile fixes $\\alpha_i=\\sqrt{n_0}\\,\\cos(i\\pi/M_{\\mathrm{CF}})$ and turns the correlator into a cosine, which is the signature of delocalized counterflow. The model also explains why the harmonic trap is essential: in a homogeneous lattice the Mott filling condition leaves no room for an impurity-hole pair, whereas the trap's compressible edges supply the extra site.","core_discovery":"The central claim is that, at intermediate impurity-bath repulsion in a harmonically confined lattice with a central Mott-insulator domain, the system enters a counterflow phase in which the bath and impurity together form an extended insulating state of unity filling, even though neither species alone shows a density plateau. In this phase the impurity's density takes the profile $n_I^{\\mathrm{CF}}(x_i)=n_0\\cos^2(\\pi(x_i-\\langle x_I\\rangle)/\\ell_{\\mathrm{CF}})$, the ground-state density of a free particle in an infinite square well of width $\\ell_{\\mathrm{CF}}$, and the anti-pair correlator decays as $C_{\\mathrm{AP}}(i)=n_0\\cos(i\\pi/M_{\\mathrm{CF}})$ rather than exponentially. The transition into this phase is marked by a sudden vanishing of the polaron residue, indicating an orthogonality catastrophe in the thermodynamic limit. The paper supports the phase with a minimal model in which the bath starts in a Mott state and each configuration contains exactly one impurity-hole pair; that model yields $n_b=1-|\\alpha_i|^2$, $n_I=|\\alpha_i|^2$, and $C_{\\mathrm{AP}}(i)=\\alpha_0^*\\alpha_i$, reproducing the combined-insulator condition and the slow correlator decay.","pith_inferences":["A direct check the paper leaves implicit: the effective well width $\\ell_{\\mathrm{CF}}$ should scale with the width of the central Mott plateau as $t/U_{bb}$ is varied at fixed $U_{bI}/U_{bb}$; measuring this scaling would separate a trap-geometry explanation from an interaction-driven one.","The predicted orthogonality catastrophe has a dynamical signature the paper does not compute: after a quench of $U_{bI}$ across the transition, the polaron overlap should decay with a diverging timescale, observable in Ramsey or radio-frequency spectroscopy.","The one-sided impurity location seen in the numerical profiles is a symmetry-breaking artifact the paper notes; a symmetric two-sided superposition should be restored in the thermodynamic limit, and larger exact-diagonalization checks could confirm the phase diagram is not an artifact of that symmetry breaking."],"forward_implications":["For baths with a central Mott-insulator domain and intermediate $U_{bI}/U_{bb}$, the impurity and bath lock into a combined unity-filling insulator in which neither species alone shows a constant-density plateau.","The impurity's cloud shape changes from a trap-controlled Gaussian to a $\\cos^2$ profile with a new length $\\ell_{\\mathrm{CF}}$, so strong correlations dominate over the harmonic confinement.","The anti-pair correlator $C_{\\mathrm{AP}}$ changes from exponential to slow trigonometric decay, meaning one impurity forms particle-hole correlations with the entire Mott domain rather than locally.","The polaron residue drops sharply at the counterflow transition, so the interacting ground state becomes orthogonal to the non-interacting one and the standard polaron quasiparticle picture fails in the thermodynamic limit.","Counterflow of this type requires a harmonic trap; a homogeneous lattice with integer filling cannot host the impurity-hole pair."],"supporting_citations":[{"why":"Supplies the characteristic-density scaling and locates the bath's superfluid-to-Mott transition at $t/U_{bb}\\approx 0.155$, defining the bath regime in which the counterflow phase appears.","marker":"[59]"},{"why":"Introduces the anti-pair correlator and the slow algebraic decay criterion that the paper uses to identify counterflow order.","marker":"[66]"},{"why":"Reports the recent experimental realization of supercounterflows in binary Mott insulators, providing the observable context for the anti-pair correlation signature.","marker":"[61]"},{"why":"Gives the non-lattice harmonically confined Bose polaron results used as the baseline for the miscible and immiscible regimes.","marker":"[51]"},{"why":"Provides the numerical description of harmonically confined Bose-Hubbard density profiles and the trap parameters used to create the central Mott domain.","marker":"[58]"},{"why":"Establishes that lattice polaron properties are controlled by $t/U_{bb}$ and $U_{bI}/U_{bb}$, the parameter set on which the phase diagram is built.","marker":"[33]"}],"fun_headline_variants":["Counterflow phase: impurity and bath form a single unity-filling insulator","Orthogonality catastrophe signals new counterflow phase in lattices","Impurity and bath hole pair into a delocalized unity-filling insulator","New phase: polaron and bath hole move together as a single insulator"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The counterflow phase rests on reading the slow decay of the anti-pair correlator in finite numerical samples as a thermodynamic algebraic order; without a finite-size scaling that rules out a finite-size cosine coherence, the phase could be an artifact of the trapped finite chain.","fun_headline_variants_meta":{"raw":{"variants":["Counterflow phase: impurity and bath form a single unity-filling insulator","Orthogonality catastrophe signals new counterflow phase in lattices","Impurity and bath hole pair into a delocalized unity-filling insulator","New phase: polaron and bath hole move together as a single insulator"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000676,"raw_usage":{"total_tokens":3077,"prompt_tokens":948,"completion_tokens":2129,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":564,"completion_tokens_details":{"reasoning_tokens":2051}},"tokens_in":564,"tokens_out":2129,"duration_ms":14361,"temperature":1.0,"reasoning_tokens":2051,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T21:26:31.706573+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"For fixed characteristic density $\\tilde{\\rho}_b=3.58$, compute $C_{\\mathrm{AP}}$ at $t/U_{bb}=0.1$ and $U_{bI}/U_{bb}=0.6$ for increasing system sizes, for example $N_b=25,40,80,160$ with $M=N_b+20$. If the fitted cosine's zero crossing or decay length fails to converge, or if $C_{\\mathrm{AP}}$ crosses over to exponential decay beyond a finite $x_i$, then the slow decay is a finite-size effect and the counterflow phase as defined collapses.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the characteristic-density scaling and locates the bath's superfluid-to-Mott transition at $t/U_{bb}\\approx 0.155$, defining the bath regime in which the counterflow phase appears."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the anti-pair correlator and the slow algebraic decay criterion that the paper uses to identify counterflow order."},{"cited_title":"Zheng, A","cited_arxiv_id":null,"evidence_quote":"Reports the recent experimental realization of supercounterflows in binary Mott insulators, providing the observable context for the anti-pair correlation signature."},{"cited_title":"Mistakidis, G","cited_arxiv_id":null,"evidence_quote":"Gives the non-lattice harmonically confined Bose polaron results used as the baseline for the miscible and immiscible regimes."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the numerical description of harmonically confined Bose-Hubbard density profiles and the trap parameters used to create the central Mott domain."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes that lattice polaron properties are controlled by $t/U_{bb}$ and $U_{bI}/U_{bb}$, the parameter set on which the phase diagram is built."}],"review_version":1}