{"id":"611cab33-a82a-4736-867c-2a8e353921d5","arxiv_id":"2602.06169","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A mobile impurity in a 1D open Hubbard chain with a fermionic bath phase-separates from particles under strong repulsion and from holes under strong attraction, with Friedel oscillations shaping intermediate regimes.","lead":"This paper studies one mobile impurity inside a small chain of interacting fermions, modeled with a Hubbard Hamiltonian and solved by exact diagonalization. It finds new phase-separation and oscillatory density patterns for the impurity that could be probed with ultracold atoms in optical lattices.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 'five regimes for any U/t>0' claim is undercut by the size-dependent FR regime: absent for M=8 and M=10, with no M>12 evidence.","rationale":"The reader's weakest assumption correctly identified the finite-size limitation, specifically noting that FR disappears for M=8 and M=10. My stress-test focuses on this as the single most load-bearing issue: the paper's headline assertion of five regimes for any finite U/t>0 is contradicted by its own Appendix B, and no evidence is given that FR survives for M>12. The hPS is indeed forced by particle-hole symmetry and is not the weak point; the Friedel-induced regimes are. The discrete peak classifier and the continuous-model check (which uses M=8, where FR is absent) reinforce rather than resolve this concern. The appropriate verdict remains CONDITIONAL, as the claim is plausible but requires a systematic finite-size and particle-number study. No change to the reader's verdict is needed.","tokens_in":20368,"tokens_out":5041,"duration_ms":50394,"concrete_test":"Repeat the phase diagram computation of Fig. 4b for M=16 and M=20 at quarter filling (Nf=4 and Nf=5 per spin), and for M=12 with Nf=4 at filling 1/3, using the same site(max)_I classifier. Check whether the FR region persists and whether the transition lines converge as M increases. Also compute M=12 with Nf=2 (filling 1/6) to verify that FR is absent when Nf=2, confirming the size/particle-number dependence. If FR cannot be identified for larger M or for Nf>3, the claim of five regimes for any finite U/t>0 fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim (Sec. III C, Fig. 4b, and Sec. V) is that five distinct interaction regimes exist 'for any finite U/t > 0'. This is supported almost entirely by exact diagonalization at M=12 with Nf=3 per spin. However, Appendix B explicitly states that M=8 and M=10 do not support the Friedel-repulsive (FR) regime, because with only two fermions per spin the bath lacks the non-centered minima needed to induce the impurity pattern. Thus the existence of FR is tied to having Nf>=3, i.e., to a specific small-system particle number, and no finite-size scaling or larger-M data are provided to show FR persists for M=16, 20, or in the thermodynamic limit. The hole phase separation (hPS) is protected by particle-hole symmetry, but the Friedel-induced regimes are not; they are boundary-and-filling sensitive. The continuous-model check in Appendix C uses M=8, which is exactly the case where FR is absent, so it cannot validate the FR claim. Without evidence that the phase diagram in Fig. 4b stabilizes for larger M and other Nf, the assertion that the five regimes are generic is not established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the ground state of a single mobile impurity immersed in a balanced spin-1/2 Fermi gas on a one-dimensional open Hubbard chain, using exact diagonalization for M = 8, 10, and 12 sites. By varying the impurity-bath interaction U_fI for fixed bath interaction U, the authors identify five regimes: a weakly interacting miscible phase, particle phase separation (pPS), hole phase separation (hPS), and two intermediate regimes (F_A, FR) in which the impurity density is claimed to be shaped by Friedel oscillations of the finite bath. The hPS regime is derived from the particle-hole symmetry relations Eqs. (7)-(8), and the phase diagram is summarized in Fig. 4(b). The paper also examines two-body correlations, von Neumann entanglement, and a continuum-lattice check in Appendix C.","tokens_in":20619,"tokens_out":6229,"duration_ms":73235,"significance":"If the claims are established, the paper would add a genuinely new symmetry-enforced phase-separation channel (hPS) to the impurity-in-lattice literature and would propose the impurity as a local probe of Friedel oscillations. The manuscript has clear strengths: the particle-hole symmetry relations (7)-(8) are exact and provide an internal consistency check that the reported quarter/three-quarter data respect; the numerics are standard, transparent ED with explicit Hilbert-space construction; and the diagnostics (double occupations, Schmidt gap, von Neumann entropy) are appropriate for few-body systems. The continuous-model appendix provides a useful check that the tight-binding conclusions survive in deep lattices. However, the central qualitative claims are extracted from very small systems, and the universality statement 'the five configurations appear for any finite U/t > 0' is not supported by any finite-size scaling and is contradicted by the paper's own Appendix B for M = 8 and M = 10.","major_comments":[{"comment":"The statement in Sec. III C that 'the five configurations appear for any finite U/t > 0' is load-bearing but not established. The phase diagram in Fig. 4(b) is built almost entirely from M = 12 with N_f = 3 fermions per spin, while Appendix B explicitly states that M = 8 and M = 10 do not support the FR regime because they have only two fermions per spin. These two sizes mix a change in M with a change in N_f (and at M = 10 also a change in filling from ν_f = 1/4 to 1/5), so they cannot be used to argue either that FR is generic or that it is an N_f = 3 artifact. A clean finite-size study at fixed ν_f = 1/4 would require M = 16, 20, ... with N_f = 4, 5, ..., or a DMRG/iDMRG calculation. Without such data, the claim that the Friedel-induced regimes are generic, and not a special few-fermion effect, remains unsupported. This is central because the abstract and title advertise Friedel-induc","section":"Sec. III C / Fig. 4(b) / App. B"},{"comment":"The causal attribution to Friedel oscillations is not uniquely supported by the data. In Fig. 3(bottom), for U/t = 0.5 the non-interacting bath has its largest Friedel maxima at i = 3 and M-2, yet the impurity at intermediate attraction still localizes at i = 2 and M-1, and the fermionic peaks move there. The text itself states in Sec. III B 4 that 'the impurity changes the pattern of the Friedel oscillations.' This is more naturally read as a self-consistent impurity-induced rearrangement of the bath, rather than the impurity passively following pre-existing Friedel oscillations. The F_A regime is also called a 'phase separation induced by Friedel oscillations,' but Fig. 5(a) shows the total impurity-fermion double occupation in F_A is about 1.5, not 0 as in pPS, so the overlap is substantial and the terminology overstates the case. Please separate the existence of the density-localizat","section":"Sec. III B 3-4, Figs. 2-3"},{"comment":"The phase boundaries in Fig. 4 are defined by a discrete observable, site_max(I), whose jumps can be artifacts of discreteness, as the authors acknowledge in Sec. III E. However, the entanglement diagnostic only partly rescues the identification: the von Neumann entropy shows a change in slope at the sA/FA, FA/w, and FR/sR transitions, but the w/FR transition has no entropy signature and the correlations (Fig. 5) are smooth across it. Thus the FR regime is identified by one discrete criterion and is exactly the regime whose existence is most size-sensitive (App. B). The paper should provide an additional, continuously varying order parameter (e.g., participation ratio of n_I, impurity-hole overlap, or low-lying energy crossings) that marks the FR region, or explicitly state that this regime is defined only by the position of the maximum of n_I.","section":"Sec. III C/E, Fig. 4(a)/Fig. 6"},{"comment":"The continuum-model check in Appendix C uses M = 8, the one lattice size for which Appendix B reports no FR regime. Therefore the continuous-model data cannot validate the FR claim; they can only validate F_A, pPS, and hPS (and indeed hPS is shown to require V_0/E_R ≳ 2). This limitation should be stated where the continuous check is summarized in Sec. V. In addition, the abstract and conclusions present hPS as a general finding without the tight-binding/deep-lattice caveat that the authors themselves establish in App. C; since particle-hole symmetry is a property of the Hubbard model, the hPS claim should be qualified in the abstract by 'in the tight-binding regime' to avoid overgeneralization.","section":"App. C / Sec. V"}],"minor_comments":[{"comment":"The definition of \\bar{x}_I in Eq. (13) is not normalized explicitly; since n_I is normalized to one, the formula is fine, but a parenthetical noting \\sum_i n_I(i)=1 would improve clarity.","section":"Sec. II B, Eq. (13)"},{"comment":"The statement that the change of the largest Friedel peak from i = 3 to i = 2 occurs around U ≈ 1.8t would benefit from a reference or a brief derivation; as written it appears as an unexplained numerical observation.","section":"Sec. III A"},{"comment":"The convergence test in Appendix C is described only by the statement that M_σ = n M modes are needed; please specify how many bands n were used for the data in Figs. 13-14 and whether convergence was checked in U_fI and V_0 as well as in the densities.","section":"App. C"},{"comment":"In the color density plots, the color scale is not defined in the caption. Adding a colorbar or stating the normalized density range would make the figures interpretable.","section":"Figs. 3 and 8"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is carefully done and the hPS mechanism is a nice, symmetry-protected result, but the headline universality of the five regimes is currently overreaching. The finite-size issue is the single reason I am not recommending acceptance: the FR regime is absent in M = 8 and M = 10, and no M > 12 data are provided. If the authors can supply a proper finite-size analysis for the Friedel regimes (or substantially recast the claims as M- and N_f-dependent few-body physics), the paper would be a solid contribution to the specialist literature."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is a straightforward, honest exact-diagonalization study of a mobile impurity in a 1D Fermi-Hubbard chain with open boundaries. The genuinely new thing is the hole phase separation (hPS): for strong attraction, the impurity localizes at the borders together with the fermions, which is counterintuitive if you think of attraction as collapsing everything to the center. The explanation via particle-hole symmetry is convincing, and the symmetry relations (7)-(8) are a useful, checkable consistency constraint that the quarter- and three-quarter-filling results respect. That alone is worth noting.\n\nThe paper also does several things well. The phase classification is not just based on one observable: density profiles, average position, double occupations, and the impurity von Neumann entropy all corroborate the boundaries, and the entropy shows genuine changes at the transitions. The continuous-lattice check in Appendix C is a good faith effort to test the tight-binding assumption, and the authors are explicit that hPS disappears in shallow lattices, as it should since the particle-hole symmetry is a Hubbard-model property.\n\nThe soft spots are real but not fatal. The central claim in Sec. III C that five regimes appear for any finite U/t > 0 is too strong. Appendix B itself says that the Friedel-repulsive (FR) regime does not exist for M=8 or M=10 because those systems have only two fermions per spin. So FR is tied to having Nf >= 3, and there is no finite-size scaling or data beyond M=12 showing it persists. The continuous-model check also uses M=8, which is exactly the case where FR is absent, so it cannot validate that regime. The hPS, in contrast, is protected by the symmetry relation and should survive to larger systems, so I would separate the robustness of hPS from the more fragile Friedel-induced localizations. The FA regime also appears only in a narrow interaction window, and the discrete site(max) classifier is admittedly heuristic, though the entropy data helps.\n\nOverall, I trust the numerics and the symmetry argument. The hPS result is new and solid; the Friedel-induced localization is a plausible suggestion that needs larger-system evidence. The paper would benefit from one systematic finite-size study, even just M=14 or M=16, and from toning down the \"any finite U/t\" phrasing. It is not circular, and the cited literature is handled fairly.\n\nI would send it to a serious referee. It is a small-system numerical survey, but the hPS mechanism is a genuine contribution to the lattice-polaron literature, and the paper is honest about its limitations. For a reading group, it would spark a good discussion about finite-size effects and symmetry in polaron problems.","headline":"The hPS result is real and clean; the Friedel regimes are suggestive but the paper overclaims its 'any finite U/t' generality given the small-system evidence.","tokens_in":21135,"tokens_out":1965,"would_cite":true,"duration_ms":21507,"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":"In a small open Hubbard chain, a mobile impurity's ground state splits into five regimes, including a hole phase separation driven by particle-hole symmetry and localizations induced by Friedel oscillations.","keywords":["Hubbard model","mobile impurity","phase separation","Friedel oscillations","exact diagonalization","fermionic bath","particle-hole symmetry","optical lattices"],"falsifier":"A systematic finite-size scaling of the phase diagram for M = 6, 8, 10, 12, 14, 16 at quarter filling would settle the claim: if the hPS or FA regions shrink and vanish as M grows, the central assertion fails. Alternatively, an ultracold-atom experiment in a 1D optical lattice with V0/ER ≈ 4 and adjustable impurity-bath scattering length could directly image the impurity density; observing the impurity at the edges for strong attraction (hPS) and at sites 2 and M-1 for intermediate attraction (FA) would confirm the paper's prediction.","tokens_in":20252,"feed_emoji":"⚛","tokens_out":6522,"duration_ms":59509,"temperature":0.7,"pith_summary":"The paper establishes that a single mobile impurity in a one-dimensional open Hubbard chain with a balanced spin-1/2 fermionic bath does not simply stay mobile or phase-separate from the bath. It shows five distinct ground-state regimes as the impurity-bath interaction is varied: a miscible weak-coupling regime, a repulsive phase separation between impurity and fermions, a hole phase separation at strong attraction where the impurity binds to empty sites rather than fermions, and two intermediate regimes where the impurity's density pattern is guided by the bath's Friedel oscillations. The hole-phase-separation regime follows from the particle-hole symmetry of the balanced Hubbard model and is absent for bosonic baths. If these results hold beyond the small lattices studied, they would make the impurity a practical probe of Friedel physics in fermionic systems and enlarge the known phase diagram of lattice polarons.","feed_headline":"Impurity separates from holes, not particles, at strong attraction","feed_subtitle":"Five regimes in a 12-site chain: impurity density tracks Friedel oscillations, and strong attraction pairs impurity with holes.","key_machinery":"The central object is the three-component Fermi-Hubbard Hamiltonian on an open chain, with one mobile impurity and N_up = N_down fermions, plus its particle-hole symmetry. In the tight-binding limit the balanced Hubbard model obeys n_sigma(i; UfI, nu_f) + n_sigma(i; -UfI, 1-nu_f) = 1 for each spin and n_I(i; UfI, nu_f) = n_I(i; -UfI, 1-nu_f) for the impurity, which forces attractive and repulsive regimes at complementary fillings. The second ingredient is the bath's Friedel oscillations — boundary-induced density waves with peaks (for nu_f < 1/2) or dips (for nu_f > 1/2) whose position depends on U — which act as an effective periodic potential for the impurity. Exact diagonalization of the","core_discovery":"Using exact diagonalization of a three-component Fermi-Hubbard Hamiltonian with open boundaries, the authors find that for quarter filling (and, by particle-hole symmetry, three-quarter filling) the ground state of one mobile impurity plus a balanced spin-1/2 fermion bath crosses through five interaction regimes for any finite U/t > 0. For strong repulsive impurity-bath coupling the impurity phase-separates from the fermions (pPS), while for strong attraction it phase-separates from the fermionic holes (hPS) — the impurity concentrates at the two edge sites while fermions and holes sort themselves accordingly, a direct consequence of the particle-hole symmetry of the Hubbard model. At interm","pith_inferences":["If the hPS and FA regimes persist under finite-size scaling, similar impurity-hole separation should appear in other particle-hole-symmetric lattice fermion models beyond the single-band Hubbard chain — e.g., in the t-J model or in fermionic mixtures with longer-range interactions, where the particle-hole relation is modified.","The role of the open boundary as the source of the Friedel oscillations suggests that in a closed ring or in a bulk system the intermediate regimes would vanish, leaving only miscible and phase-separated phases; a direct comparison between open and periodic boundary conditions would cleanly test this.","Since the FA regime behaves like a local hPS with the impurity pinned next to the boundary, the paper implicitly predicts that the local double-occupation maximum at sites 2 and M-1 reaches about 1.5; an experimental measurement of the on-site pair correlation at those sites would distinguish FA from hPS.","The phase diagram is computed for balanced baths; an imbalanced bath or higher filling (e.g., nu_f = 1/3) would likely shift the Friedel peak positions and could either widen or destroy the FR/FA windows, providing a control knob for the intermediate regimes."],"forward_implications":["The same five regimes appear at both quarter and three-quarter filling, with attractively and repulsively induced phases exchanged, because of particle-hole symmetry; experiments can accordingly look for hPS at either sign of the interaction depending on filling.","The impurity's localizations at sites i=2 and i=M-1 for intermediate attraction (FA) and at the sites between Friedel peaks (FR) offer a direct way to image the bath's Friedel oscillations: the impurity density pattern is a readout of the bath density.","A continuous optical-lattice calculation confirms the hPS and FA regimes survive for lattice depths V0/ER ≳ 2, meaning the results are not purely an artifact of the tight-binding model, and suggests feasible parameters for ultracold-atom experiments.","The von Neumann entropy of the impurity peaks at the transition points between regimes, giving a measurable entanglement signature of the phase boundaries beyond density measurements."],"fun_headline_variants":["Impurity bonds with holes at strong attraction","Friedel oscillations pattern impurity in Hubbard chain","Five regimes for a mobile impurity in a fermionic bath","Attraction pairs impurity with holes, not fermions","Particle-hole symmetry flips impurity separation"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The few-site exact-diagonalization results (M = 8, 10, 12 with at most three fermions per spin) are assumed to represent the physics of larger or continuous systems, with no systematic finite-size scaling shown for the phase boundaries; in particular the FR regime disappears for M = 8 and M = 10.","fun_headline_variants_meta":{"raw":{"variants":["Impurity bonds with holes at strong attraction","Friedel oscillations pattern impurity in Hubbard chain","Five regimes for a mobile impurity in a fermionic bath","Attraction pairs impurity with holes, not fermions","Particle-hole symmetry flips impurity separation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001039,"raw_usage":{"total_tokens":4170,"prompt_tokens":668,"completion_tokens":3502,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":412,"completion_tokens_details":{"reasoning_tokens":3431}},"tokens_in":412,"tokens_out":3502,"duration_ms":24668,"temperature":1.0,"reasoning_tokens":3431,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T03:58:18.851148+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A systematic finite-size scaling of the phase diagram for M = 6, 8, 10, 12, 14, 16 at quarter filling would settle the claim: if the hPS or FA regions shrink and vanish as M grows, the central assertion fails. Alternatively, an ultracold-atom experiment in a 1D optical lattice with V0/ER ≈ 4 and adjustable impurity-bath scattering length could directly image the impurity density; observing the impurity at the edges for strong attraction (hPS) and at sites 2 and M-1 for intermediate attraction (FA) would confirm the paper's prediction.","supporting_citations":[],"review_version":1}