{"id":"78572268-f5e4-49fd-a1b4-ae4c0c1068c9","arxiv_id":"2608.06233","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A mobile impurity binds to a pinned quasihole in a lattice fractional Chern insulator, and the binding-energy ratio gives a direct readout of the quasihole's fractional charge.","lead":"An impurity atom in a lattice fractional Chern insulator binds to a pinned quasihole, and the binding energy scales with the quasihole's fractional charge in the weak-coupling limit. The paper shows numerically how to read out that charge and drag the composite object, a step toward anyon braiding in cold-atom simulators.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Weak-coupling binding-energy ratio appears to be a finite-size uniform-density shift (E_B ≈ U Q / L²); without L-scaling at small U, the 'direct fractional-charge measure' claim is not established in the thermodynamic limit.","rationale":"The paper is a serious numerical and analytical study of impurity–quasihole binding in lattice FCIs, and the transport protocol is a useful independent contribution. The reader's highlighted assumption is the mean-field replacement H_U → U∑⟨n⟩a†a, with back-action validated indirectly through δE_B ≈ 0 for U/J < 0.5. I agree that concern is real, but I see a more fundamental and less appreciated issue: even if the mean-field replacement is exact, the weak-coupling binding energy in a finite system is dominated by the first-order overlap of the extended impurity with the host density depletion, giving E_B ≈ U Q / L². This produces the clean charge ratio that the paper reports, but it is a finite-size volume effect rather than a true two-body bound-state energy. In the thermodynamic limit, a two-dimensional finite-range potential binds for arbitrarily weak attraction, but with an exponentially small energy in 1/U, so the ratio of binding energies for different charges has no reason to equal the charge ratio. The paper's own finite-size study of E_B is performed at U = 2J_imp (a strong-coupling, genuinely bound regime), not in the weak-coupling window on which the charge-extraction claim rests. Therefore the central assertion requires an explicit demonstration that the weak-coupling ratio survives as L → ∞ (or a restatement that the protocol is inherently finite-size). The proposed single-particle L-scaling test at fixed small U settles this directly; if E_B ∝ 1/L², the claim must be reworded and the practical utility for larger simulators re-evaluated. The paper's strengths—the MPS/DMRG numerics, the window-MPS analysis, and the dragging demonstration—are not undermined by this concern, which is why a conditional verdict, rather than rejection, is appropriate.","tokens_in":17239,"tokens_out":21931,"duration_ms":234000,"concrete_test":"Compute E_B for the single-particle Hamiltonian H''_imp with the integer-charge (Q = −1) and the renormalized-charge (Q′ = −1/2) density profiles, at a fixed weak coupling U/J = 0.1, for system sizes L = 10, 14, 20, 28, 40 with open boundaries. Plot E_B versus 1/L² and E_B(−1/2)/E_B(−1) versus L. If E_B scales as 1/L² and the ratio remains 0.5 for all L, the weak-coupling result is the finite-size uniform-density shift; the charge-extraction protocol should be reframed as a finite-size energy-shift measurement, and its extension to larger simulators requires careful restatement. If E_B saturates to a finite value as L grows (i.e., E_B·L² increases) and the ratio still approaches 0.5, the genuine bound-state mechanism is confirmed. A complementary many-body check: repeat the FCI DMRG binding-energy ratio at U/J = 0.1 for 10×10 and 12×12 (or 14×14 if feasible) lattices at ν = 1/2.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The single-particle analysis (Fig. 1(d)) that underlies the charge-extraction claim is performed on finite lattices (20×20 in Fig. 1, 10×10 in Fig. 2). In the weak-coupling limit U→0, the binding energy defined as E_gs(V0) − E_gs(0) is dominated by first-order perturbation theory: for an unperturbed impurity ground state with uniform density 1/L² per site, E_B^(1) = U ∑_x n_{Vh}(x)/L² = U Q / L². This is linear in U and proportional to the integrated host charge Q, so the ratio E_B(Q′)/E_B(Q) trivially equals Q′/Q for any finite L. It is not the binding energy of a localized bound state; the impurity wavefunction remains extended over the whole system. In the thermodynamic limit, a finite-range attractive potential in 2D always binds for arbitrarily small U, but the binding energy is exponentially small in 1/U (e.g., ~ z0 exp(−2/z0²) with z0 ∝ √U for a shallow well), so the ratio of two such energies is not Q′/Q. The paper's finite-size scaling of E_B (inset of Fig. 1(c)) is shown only at the stronger coupling U = 2J_imp, where the bound state is deep and E_B saturates; no L-scaling is reported for the weak-coupling regime where the charge-extraction claim operates. Thus the central claim that the binding-energy ratio 'provides a direct measure of the fractional charge' is currently supported only as a finite-size volume effect, not as a bound-state property. If this is the intended mechanism, the claim should be restated; if not, the missing L-dependence at small U must be supplied.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This Letter studies a single mobile impurity coupled to a pinned quasihole in the interacting Harper–Hofstadter model in the fractional Chern insulator (FCI) regime. The authors first analyze a single particle in a 2D potential well and then map the impurity–host interaction onto an effective single-particle potential U<n> at the mean-field level, arguing that the binding energy of the impurity to the quasihole is controlled by the integrated host density defect. Using DMRG on a 10×10 lattice for hard-core bosons at ν=1/2 and a window-MPS approach on infinite cylinders for ν=1/2 and ν=1/3, they compute the many-body binding energy E_B and find that in the weak-coupling limit the ratio E_B(-1/2)/E_B(-1) approaches 1/2 and E_B(-1/3)/E_B(-1/2) approaches 2/3. They also simulate dragging of the pinning potential with TDVP and show that an impurity at U=3J follows the moving quasihole. The appendix contains the 2D binding-energy derivation, host characterization via the Streda marker, charge extraction from density profiles, window-MPS details, and the drag-dynamics analysis.","tokens_in":17660,"tokens_out":30767,"duration_ms":287046,"significance":"The paper is a serious numerical study with several concrete strengths: the many-body DMRG and window-MPS results go beyond the mean-field illustration and break its circularity by computing E_B from the full Hamiltonian while extracting the charge independently from the host density; the δE_B diagnostic (Eq. (9) and Fig. 2(c)) offers a quantitative handle on back-action; the window-MPS comparison of ν=1/2 and ν=1/3 is a nontrivial two-system test; and the dragging protocol connects directly to ongoing cold-atom quantum-gas-microscope experiments, with a self-contained appendix. If the charge-extraction claim holds in the stated regime, the scheme would be a useful in-situ probe of fractional charge in engineered FCI simulators. However, as detailed below, the weak-coupling 'binding energy' that carries the claim is currently a finite-size linear-response energy shift rather than a demonstrated bound-state property, the calibration and error budget are not quantified, and the drag simulation operates in a different (strong-coupling) regime. The significance of the result is therefore real but more modest than the abstract suggests.","major_comments":[{"comment":"The charge-extraction claim is currently supported only as a finite-size first-order energy shift, not as a bound-state property. For U→0, the ground-state energy of H''_imp = H_imp + U Σ n_Vh a†a is shifted by the first-order term E_B ≈ −U Σ_x δn(x)|ψ_0(x)|², with δn = n_0 − n_Vh and ψ_0 the unperturbed impurity ground state. For the delocalized ψ_0 of a finite lattice this shift is of order U|Q|/L², linear in U and proportional to the integrated charge Q; the ratio E_B(Q')/E_B(Q) therefore approaches Q'/Q for any finite L (Fig. 1(d)) regardless of whether the impurity wavefunction is localized. In the thermodynamic limit at fixed U, the 2D bound-state binding energy is exponentially small in 1/(U|Q|) and its ratio is not Q'/Q. The only L-dependence shown (inset of Fig. 1(c)) is at U = 2J_imp, where the well is deep, not in the weak-coupling regime where the claim operates. I ask for either (a) an L-scaling study of E_B and its ratio at small U (U ≲ 0.5J) that demonstrates the behavior beyond the 10–20-site systems used here, or (b) an explicit restatement of the result as a finite-size linear-response probe with the regime of validity and the operational meaning of 'binding energy' spelled out.","section":"Mean field analysis; Many-body calculations (Figs. 1(d), 2(d))"},{"comment":"The ratio E_B(ν=1/2)/E_B(ν=1) in Fig. 2(d) compares binding energies obtained in two different host systems (10 hard-core bosons in a ν=1/2 FCI and 21 free fermions in a CI), whose quasihole density profiles differ (Fig. 2(a)). At first order the ratio equals (Q_FCI/Q_CI) times the ratio of the impurity-density overlaps with the two defects, so the extracted 'fractional charge' carries a shape-dependent systematic error. The manuscript does not report the residual deviation of the MPS ratio from 1/2, the magnitude of the overlap correction, or error bars, so the precision implied by 'direct measure' is not quantified. Please state the numerical values of the ratios at the smallest U, the corresponding uncertainties, and a justification that the shape correction is negligible or controlled.","section":"Many-body calculations, Eq. (7), Fig. 2(d)"},{"comment":"The dragging demonstration is performed at U/J = 3.0, outside the weak-coupling window (U/J ≲ 0.5) in which δE_B ≈ 0 and the mean-field mapping was validated (Fig. 2(c)); at U = 3J the back-action dressing is significant, so the transported density minimum is not guaranteed to be the clean Q = −1/2 quasihole whose binding physics was analyzed. Appendix Fig. 9(e) restricts U ≤ 3J to 'isolate a single fractional charge of −1/2', but the text does not state the charge value at U = 3J during the drag or how much of the impurity's displacement (≈1 lattice constant after τ = 25ℏ/J, about one third of the host displacement) reflects coherent binding rather than the direct effect of the moving Gaussian potential on the impurity. Please quantify these points or soften the transport claim accordingly.","section":"Dragging the anyon-impurity composite (Fig. 4)"},{"comment":"None of the numerical sections reports bond dimensions, truncation errors, convergence of E_B with bond dimension, or TDVP time-step/truncation checks, and no error bars appear in Figs. 2–3. Since the quantitative claims (E_B ratios approaching 1/2 and 2/3) rest on these MPS simulations, at least a per-point estimate of the truncated weight and a statement of the bond dimension used are needed to establish that the reported ratios are converged.","section":"Numerical sections (Figs. 2–3)"}],"minor_comments":[{"comment":"The sentence 'In both the weak and strong confinement limits, the binding energy exhibits an approximately linear dependence on V0' is not correct for the continuum result Eq. (5): in the weak limit the 2D well binding energy is exponentially small in 1/z0², and the near-linear behavior seen in the lattice numerics reflects the finite level spacing; please rephrase.","section":"Bound states in 2D"},{"comment":"Fig. 1(d) and the surrounding text put the fractional charge in by hand as a renormalized trap depth Q'; the paper should state explicitly that this panel is an illustrative consistency check, and that the physical claim rests on the many-body calculations (Figs. 2(d), 3(e)) where the charge is obtained independently from the density.","section":"Mean field analysis (Fig. 1(d))"},{"comment":"Please clarify the impurity hopping J_imp used in Figs. 2–4; the inset of Fig. 1(a) uses J_imp/J = 0.6 and 1, and the binding energy scales with J_imp, so the value used in the many-body simulations should be stated explicitly.","section":"Many-body calculations (Fig. 2(d))"},{"comment":"For the ν=1/2 window the impurity is confined by an infinite potential on the last rung, while for ν=1/3 it is not; the resulting different impurity ground-state envelopes affect the first-order shifts and hence the L_w-dependent ratio in Fig. 3(e). Please state this explicitly and check that it does not bias the extrapolation to large L_w.","section":"Domain wall excitation and infinite MPS (Fig. 3)"},{"comment":"The conclusion states that the fractional charge is to be 'inferred directly from in-situ density measurements', while the proposed protocol in the paper infers it from the binding energy; please correct this sentence to refer to the binding-energy measurement.","section":"Conclusion"},{"comment":"Since E_B is negative for all charges considered, the ratio in Fig. 1(d) is a ratio of magnitudes; please state the sign convention explicitly in the text and captions.","section":"Fig. 1(d) caption"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid numerical study with a useful finite-size probe, but the headline claims are ahead of the evidence: the weak-coupling 'binding energy' is essentially a first-order linear-response shift, and the title-level claim ('bound state', 'direct measure of the fractional charge') needs either new L-scaling data at small U or a more careful statement of the finite-size regime. The transport section is a proof-of-principle at U=3J, a regime where back-action is non-negligible and the charge-extraction claim was not validated. The absence of convergence data and error bars is a fit-to-journal issue for a numerical Letter. I would also note that the braiding-motivation language in the abstract goes beyond what a dragging simulation can support. None of this appears fatal: with the requested additions (L-scaling at small U, error bars, charge during the drag, honest restatement of the regime), the paper could be acceptable for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a solid numerical study of impurity-quasihole binding in the interacting Harper-Hofstadter model, with a nice window-MPS treatment on infinite cylinders and a dragging simulation. The catch is that the paper's headline claim—that the binding-energy ratio directly measures fractional charge—is currently supported only by a finite-size volume effect, not by a localized bound-state property.\n\nWhat's genuinely useful: the lattice implementation, the UQ matching argument, and the transport protocol are practical additions to the cold-atom FCI toolbox. The window-MPS results for nu=1/2 and 1/3 on the cylinder are a good technical step, and the authors are upfront that a similar continuum calculation exists (Ref. [61]). The many-body binding energies are computed from the full Hamiltonian, so the circularity concern is mostly addressed: Q is extracted independently from the density.\n\nThe soft spot is the weak-coupling limit. For U->0, the binding energy defined as E_gs(V0)-E_gs(0) is first-order in U: E_B ~ U Q / L^2, because the impurity ground state is essentially uniform on the finite lattice. The ratio E_B(Q')/E_B(Q) then equals Q'/Q for any finite L, independent of any bound-state structure. The impurity wavefunction remains extended; it has not localized onto the quasihole. The paper's Fig. 1(d) shows exactly this trivial limit. The inset of Fig. 1(c) shows L-scaling only at strong coupling (U=2J_imp), where the bound state is deep. No L-scaling is reported in the weak-coupling regime where the charge-extraction claim operates. So the central claim is currently a restatement of the integrated density deficit, not a new probe. In the thermodynamic limit the shallow-well binding energy is exponentially small in 1/U, and the ratio would not equal the charge ratio. This is a load-bearing gap.\n\nThe dragging simulation is interesting but also relies on the same binding; for U/J=3 the impurity displacement is only about one lattice constant, and there are no error bars or convergence data reported for the MPS calculations anywhere.\n\nBottom line: this is a credible, well-executed numerics paper with an overstated interpretation. It deserves peer review—a good referee can push for the missing L-scaling or a careful restatement of what is being measured. I would not cite the charge-extraction claim until that is fixed, but I would cite the window-MPS and dragging methods.","headline":"A careful numerics paper whose headline claim—binding-energy ratios measure fractional charge—is currently a finite-size volume effect, not a bound-state property.","tokens_in":18166,"tokens_out":4611,"would_cite":false,"duration_ms":45899,"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":"The paper claims that in a fractional Chern insulator, the binding energy of a mobile impurity to a pinned quasihole is proportional to the quasihole's fractional charge in the weak-coupling limit, making binding-energy ratios a direct…","keywords":["fractional Chern insulator","quasihole","fractional charge","mobile impurity","anyon-impurity bound state","Harper-Hofstadter model","binding energy","matrix product states"],"falsifier":"In a $10\\times 10$ hard-core-boson Harper-Hofstadter system at $\\nu=1/2$, compute the full many-body binding-energy ratio $E_B(-1/2)/E_B(-1)$ at $U/J<0.5$; if the ratio is not $1/2$ within finite-size error, or if $\\delta E_B$ does not vanish in that window, the linear charge-extraction claim fails.","tokens_in":17077,"feed_emoji":"⚛️","tokens_out":7809,"duration_ms":68485,"temperature":0.7,"pith_summary":"Mobile impurities are established probes of correlated matter, and this paper asks whether they can reveal the fractional charge of a quasihole in an engineered fractional Chern insulator. The authors argue that they can: in the interacting Harper-Hofstadter model, a repulsive impurity binds to a pinned quasihole, and in the weak-coupling limit the binding energy is proportional to the quasihole's charge. They support this with mean-field analysis and large-scale matrix-product-state simulations for bosonic ν=1/2 and fermionic ν=1/3 Laughlin-type states, finding that the ratios E_B(−1/2)/E_B(−1) and E_B(−1/3)/E_B(−1/2) approach 1/2 and 2/3. If correct, this gives a practical, in-situ way to measure fractional charge and to manipulate anyon-impurity composites by dragging the pinning potential.","feed_headline":"Impurity binding energy exposes anyon charge in Chern insulators","feed_subtitle":"In fractional Chern insulators, binding-energy ratios approach 1/2 and 2/3, a practical charge probe.","key_machinery":"The load-bearing mechanism is the mean-field replacement of the impurity-host repulsion by a single-particle potential $U\\sum_{x,y}\\langle \\hat n_{x,y}\\rangle \\hat a^\\dagger_{x,y}\\hat a_{x,y}$, built from the host ground-state density. Since a two-dimensional attractive potential always binds a particle, this maps the anyon-impurity problem onto a shallow trap whose integrated strength is $UQ$; the binding energy depends on trap depth and grows linearly with $Q$ in the weak-coupling limit. The paper validates this mapping by comparing mean-field curves with full many-body matrix-product-state results and by using $\\delta E_B$, the difference between bare and effective binding energy, to estimate the onset of collective-mode dressing.","core_discovery":"The paper's central claim is that the binding energy of a charge-neutral mobile impurity to a single pinned quasihole is set by the quasihole's charge in the weak-interaction regime, so the ratio of binding energies for two quasihole charges equals the ratio of the charges themselves. This is demonstrated in the interacting Harper-Hofstadter model: hard-core bosons at ν=1/2 host a quasihole with charge −1/2, a Chern insulator supplies an integer charge −1, and an infinite-cylinder fermionic ν=1/3 state supplies a quasihole of charge −1/3. The numerical binding-energy ratios track the expected fractional values once the impurity-host coupling is weak enough that the impurity's back-action on the host, estimated by the difference between bare and effective binding energies, nearly vanishes. The same bound state can be moved by translating the pinning potential, giving a physical route to steer anyonic excitations.","pith_inferences":["Because only ratios of binding energies enter, the protocol could work as a calibration-free diagnostic in systems where absolute energy scales are hard to fix; the paper computes absolute energies but does not make this inference.","The same mean-field trap picture suggests that the ratio test could be applied to other localized excitations, such as domain walls or quasiparticles in non-Laughlin states, but the paper only tests Laughlin-type ν=1/2 and ν=1/3 quasiholes.","A natural experimental follow-up is to use the dragging protocol to move one composite around another and look for braiding phases; the paper demonstrates transport but not braiding."],"forward_implications":["If the linear ratio holds, binding-energy spectroscopy gives a direct measurement of quasihole charge without relying on transport or density-integration protocols.","The bound composite can be transported by moving the pinning potential, so anyon-impurity composites become movable objects suitable for braiding-type manipulations in cold-atom simulators.","Monitoring the breakdown of adiabatic following during a drag provides a dynamical estimate of the binding energy.","The same protocol works for both bosonic ν=1/2 and fermionic ν=1/3 Laughlin-type states, suggesting generality across fractional Chern insulator platforms."],"supporting_citations":[{"why":"Defines the anyon-trion binding-energy framework and the binding-energy definition that this work transfers to lattice fractional Chern insulators.","marker":"[15]"},{"why":"Earlier proposal to extract fractional charge and anyonic statistics from impurities binding to fractional quantum Hall excitations.","marker":"[9]"},{"why":"Provides spectral fingerprints of exciton-anyon binding in fractional Chern insulators, supplying the finite-size and binding-energy methodology.","marker":"[17]"},{"why":"Reports the observation of anyon-trions in twisted MoTe2, the experimental motivation for engineered anyon-impurity composites.","marker":"[23]"},{"why":"Shows how lattice quasiholes are created and their charge measured, grounding the pinned-quasihole setup used here.","marker":"[35]"},{"why":"Provides a patchwork construction for growing extended Laughlin states in quantum gas microscopes, supporting experimental feasibility.","marker":"[37]"},{"why":"Interferometric braiding protocol for anyons in Chern insulators that motivates steering anyon-impurity composites.","marker":"[20]"},{"why":"The DMRG algorithm used for the finite-system many-body ground-state calculations.","marker":"[51]"},{"why":"The time-dependent variational principle algorithm used for the real-time dragging simulations.","marker":"[65]"}],"fun_headline_variants":["Anyon-impurity binding exposes quasihole charge","Impurity binding energy measures anyon charge","Steered anyon-impurity states in Chern insulators"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument assumes the impurity acts on the host only through the host's ground-state density—the mean-field replacement $\\hat H_U \\to U\\sum_{x,y}\\langle \\hat n_{x,y}\\rangle \\hat a^\\dagger_{x,y}\\hat a_{x,y}$—so impurity back-action on the fractional Chern insulator is negligible where the charge is read out.","fun_headline_variants_meta":{"raw":{"variants":["Anyon-impurity binding exposes quasihole charge","Impurity binding energy measures anyon charge","Steered anyon-impurity states in Chern insulators"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000213,"raw_usage":{"total_tokens":1410,"prompt_tokens":921,"completion_tokens":489,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":537,"completion_tokens_details":{"reasoning_tokens":438}},"tokens_in":537,"tokens_out":489,"duration_ms":4976,"temperature":1.0,"reasoning_tokens":438,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T11:52:10.376511+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"In a $10\\times 10$ hard-core-boson Harper-Hofstadter system at $\\nu=1/2$, compute the full many-body binding-energy ratio $E_B(-1/2)/E_B(-1)$ at $U/J<0.5$; if the ratio is not $1/2$ within finite-size error, or if $\\delta E_B$ does not vanish in that window, the linear charge-extraction claim fails.","supporting_citations":[{"cited_title":"Mostaan, N","cited_arxiv_id":null,"evidence_quote":"Defines the anyon-trion binding-energy framework and the binding-energy definition that this work transfers to lattice fractional Chern insulators."},{"cited_title":"Mu˜ noz de las Heras, E","cited_arxiv_id":null,"evidence_quote":"Earlier proposal to extract fractional charge and anyonic statistics from impurities binding to fractional quantum Hall excitations."},{"cited_title":"Lu and L","cited_arxiv_id":null,"evidence_quote":"Provides spectral fingerprints of exciton-anyon binding in fractional Chern insulators, supplying the finite-size and binding-energy methodology."},{"cited_title":"Macaluso, T","cited_arxiv_id":null,"evidence_quote":"Shows how lattice quasiholes are created and their charge measured, grounding the pinned-quasihole setup used here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides a patchwork construction for growing extended Laughlin states in quantum gas microscopes, supporting experimental feasibility."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Interferometric braiding protocol for anyons in Chern insulators that motivates steering anyon-impurity composites."}],"review_version":1}