{"id":"1cebf585-124c-46ce-8507-e03da0dcf48d","arxiv_id":"1908.02186","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Doping the Z2 Bose-Hubbard model creates Z_n solitons that bind 1/2 or 1/3 fractionalized bosons, with a generalized bulk-defect correspondence established via inter-soliton pumping.","lead":"This paper shows that a chain of interacting bosons coupled to Ising spins can spontaneously form topological walls called solitons when doped away from special fillings, and that these walls bind fractional pieces of bosons. It offers a concrete cold-atom model for observing fractionalization of matter in a clean system, avoiding the disorder problems that have plagued polymer experiments.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Pumping-cycle gap is unchecked: the Chern-number assignment in Eq. (14) and the generalized bulk-defect correspondence assume the interacting ground state of Eq. (13) remains gapped for all φ, but no gap spectrum is shown and δ is fixed to t despite claiming δ ≫ t.","rationale":"The reader's weakest assumption matches the concern I would raise: the interacting pumping cycle in Sec. III B 2 is used to define Chern numbers via Eq. (14), but the paper never shows that the many-body gap remains open along the cycle. This is not a stylistic omission; the generalized bulk-defect correspondence is precisely what elevates the observed fractionalized bound states from a numerical observation to a topologically robust phenomenon. Without the gap check, the Chern numbers could be artifacts of a level crossing, and the statement that the 1/3 bosons are remnants of 2D edge states would lack its quantitative foundation. I do not think this invalidates the equilibrium fractionalization evidence, which is supported by the integrated-density plateaus in Fig. 5 and the exponential repulsion in Fig. 7. The concern is therefore conditional: the central claim stands if the gap is verified, and requires revision if it is not. Since the reader already assigned CONDITIONAL for essentially this reason, I recommend no change to the verdict. A secondary issue, undocumented DMRG convergence, is real but less load-bearing because the fractionalization plateaus are fairly direct observables; the gap issue is the one that the topological interpretation cannot do without.","tokens_in":26296,"tokens_out":5333,"duration_ms":62891,"concrete_test":"Recompute the many-body energy gap E_1(φ)-E_0(φ) for the Hamiltonian (13) with δ=t, using the same L=42 chain and the same U, Δ, β, and pinning parameters as in Figs. 11-12. Sample at least 100 φ points across [0,2π) and verify the gap remains positive; repeat for L=60 to check finite-size scaling. If the gap closes at any φ, rerun the Berry-phase integration with a smaller δ or a modified cycle and demonstrate that the pumped charge is quantized only when the gap is open.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing step is the assumption that the interacting pumping cycle in Eq. (13) keeps the many-body gap open. The Chern numbers extracted from Eq. (14), and the resulting generalized bulk-defect correspondence for the Z4 solitons, require a smooth, non-degenerate evolution of the many-body ground state over the full 2π cycle. The paper states 'By choosing δ≫t, we guarantee that the spins rotate periodically...' and then immediately says 'In practice, it is enough to fix δ=t.' With δ=t the superlattice modulation is not parametrically large, so the adiabatic assumption is not automatically justified. No many-body gap spectrum, level-crossing check, or finite-size scaling along φ is presented. If the gap closes at any φ, the Berry-phase integral in Eq. (14) is not a well-defined quantized Chern number, the assignments ν_A=-1, ν_B=1, ν_C=1 in Fig. 12 are not justified, and the claim that the 1/3-bound bosons are protected by a generalized bulk-defect correspondence does not follow. The equilibrium fractionalization evidence in Figs. 5-7 is independent of this pumping argument, so the fractionalization results may survive, but the topological-origin claim pivots on the unchecked gap.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the Z2 Bose-Hubbard model (Eq. (1)) and claims that, when doped away from commensurate fillings, the ground state contains dynamical Z_n topological solitons that bind fractionalized bosonic quasi-particles (1/2 bosons for Z2 solitons around half filling, 1/3 bosons for Z4 solitons around one- and two-third fillings). It further argues that these bound quasi-particles repel, forming a fractional soliton lattice at higher doping, and that their topological origin can be understood through an adiabatic inter-soliton pumping protocol whose quantized transport establishes a generalized bulk-defect correspondence. The evidence for the equilibrium fractionalization includes DMRG density plateaus at 1/2 and 1/3, fits to tanh and sech^2 profiles, and an energy-distance curve that excludes a polaron minimum. The pumping argument is used to assign Chern numbers to the homogeneous symmetry-broken sectors and to explain why Z4 solitons, which separate regions of equal Berry phase, nevertheless host protected bound states.","tokens_in":26534,"tokens_out":5428,"duration_ms":57906,"significance":"If the central claims hold, this is a significant contribution to the physics of symmetry-protected topological defects in strongly correlated systems. The equilibrium results—spontaneous Zn solitons, fractionalized bosonic bound states, and a numerically exact fractional soliton lattice—are well supported by the presented DMRG data and are interesting in their own right, particularly in the context of cold-atom implementations. The paper also makes a conceptual claim: that interacting solitons can be understood via a generalized bulk-defect correspondence based on quantized inter-soliton pumping. The presentation of explicit numerical fits to known quantum-field-theory profiles and the quantitative density-plateau diagnostics are strengths. However, the pumping-based topological origin is the least supported part of the manuscript, and it is load-bearing for the paper's central narrative.","major_comments":[{"comment":"The interacting pumping cycle is assumed to keep the many-body gap open for all φ, but no gap spectrum, level-crossing check, or finite-size scaling of the gap along the cycle is presented. The text states 'By choosing δ≫t, we guarantee that the spins rotate periodically' and then immediately 'In practice, it is enough to fix δ=t'; with δ=t the superlattice modulation is not parametrically large, so the adiabatic assumption is not automatically justified. If the gap closes at any φ, the Berry-phase integral in Eq. (14) is not a well-defined quantized Chern number, and the assignments ν_Ā=-1, ν_B=1, ν_C=1 in Fig. 12 do not follow. Since the generalized bulk-defect correspondence for the Z4 solitons is a central claim, this missing gap analysis is a load-bearing issue. The equilibrium fractionalization results of Figs. 5-7 are independent of this pumping argument and are not affected by this comment.","section":"Sec. III.B.2, Eq. (13)"},{"comment":"The pumping argument is developed for the topological configurations Ā, ¯B, ¯C (the U=15t TBOW phases in Fig. 8), whereas the 1/3-fractionalized bound states whose origin the argument is meant to explain were obtained in Sec. II for the trivial BOW phase at U=10t (Figs. 5-7). The manuscript does not specify the interaction strength used in the pumping calculations, nor does it demonstrate that the same A-B-C-A defect structure with 1/3 bound bosons persists in the topological phase. Without an explicit connection between these parameter regimes, the Chern numbers computed for the topological sectors do not explain the fractionalization observed in the trivial BOW phase.","section":"Sec. III.B.2 and Figs. 11-12"}],"minor_comments":[{"comment":"There are numerous typos and misspellings, including 'relatisvistic' (Sec. II.A), 'inebitably' (Sec. II), 'stablished' (Sec. III.B), 'distnace' and 'miniminum' (Sec. II.C), and 'Peielrs' (Conclusions). A careful proofread is needed.","section":"Throughout"},{"comment":"The bond dimension is stated as D=100, but no convergence analysis in D is shown for the key quantities, such as the 1/2 and 1/3 density plateaus in Fig. 5 or the energy-distance curve in Fig. 7(a). Providing such a check would strengthen the numerical claims.","section":"Sec. II, DMRG parameters"},{"comment":"The notation ⟨: n_j :⟩ is used for the density deviation from ρ*, but this normal-ordering symbol is never defined; it should be explicitly stated as ⟨n_j⟩ - ρ*.","section":"Eq. (10)"},{"comment":"The x-axis label 'β0/β' is related to the pinning strength ε only through the definition β0 = β(1+ε) below Eq. (2); the text should make this relationship explicit in the figure caption.","section":"Fig. 3(c) and pinning definition"},{"comment":"Reference [1] is cited as unpublished and is used to motivate the extension beyond half filling; the final version should clarify its status (e.g., companion paper) or provide the relevant results explicitly.","section":"Reference [1]"}],"recommendation":"major_revision","confidential_remarks":"The main risk to the paper's central narrative is the unchecked gap in the interacting pumping cycle. I would ask the authors to provide a many-body gap spectrum along φ for the parameters used in Figs. 11-12, or at least a finite-size scaling of the gap, and to clarify the parameter regime connecting the pumping argument to the fractionalization results of Sec. II. The equilibrium fractionalization and soliton-lattice results are solid and could support a strong paper even if the pumping claim is weakened or reframed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a good paper with one under-supported load-bearing claim. The equilibrium results – Z2 and Z4 solitons, fractionalized boson densities, soliton lattice – are backed by clean DMRG data and convincing fits to the expected QFT profiles. The new content beyond the companion paper is the Z4 case at 1/3 and 2/3 fillings plus the generalized bulk-defect correspondence via inter-soliton pumping. The half-filling fractionalization is attributed to Ref. [1], which is awkward but not disqualifying if the paper directs to an available companion.\n\nWhat's genuinely good: the density profiles in Figs. 5 and 7 match sech^2 and the integrated densities give clean 1/2 and 1/3 jumps. The energy-versus-distance curve in Fig. 7(a) reasonably excludes a polaron minimum. The local Berry phase calculations are a nice way to visualize the topology of the different SSB sectors. This is careful, reproducible-looking numerics, though no code/data is shipped.\n\nThe soft spot is the pumping argument in Sec. III.B.2. The paper states that δ≫t guarantees periodic spin rotation, then fixes δ=t. That alone is fine if the gap stays open, but no many-body gap spectrum or level-crossing check along φ is shown for the interacting cycle. The Chern numbers from Eq. (14) and the generalized bulk-defect correspondence for the Z4 solitons depend on this. If the gap closes anywhere, the assignment νA=-1, νB=1, νC=1 in Fig. 12 is not established. The fractionalization evidence does not depend on this, so the central physics likely survives, but the topological-origin claim should come with a gap analysis or at least a finite-size scaling statement.\n\nThe other weakness, minor by comparison, is the absence of bond-dimension convergence tests. D=100 with n0=2 is plausible for these parameters, but the paper should show it for the key observables.\n\nOverall: the paper deserves refereeing. It addresses a real question, the main numerical results look credible, and the gap issue is addressable in revision. I'd send it to a serious referee and ask for the gap check plus convergence data.","headline":"The fractionalization results are solid and the Z4 soliton extension is real, but the pumping-based bulk-defect correspondence needs a gap check before the topological-origin claim is fully justified.","tokens_in":27116,"tokens_out":3626,"would_cite":false,"duration_ms":36517,"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":"Doped ground states of the $\\mathbb{Z}_2$ Bose-Hubbard model host topological solitons that split each added boson into 1/2 or 1/3 charges, with the protection explained by a generalized bulk-defect correspondence from quantized…","keywords":["Z2 Bose-Hubbard model","topological solitons","boson fractionalization","symmetry-protected topological defects","Thouless pumping","generalized bulk-defect correspondence","fractional soliton lattice","cold-atom quantum simulation"],"falsifier":"Compute, with a matrix-product-state algorithm, the lowest two many-body energies of Hamiltonian (13) for $\\delta=t$ and the parameters used for the two-third-filling solitons, scanning $\\phi\\in[0,2\\pi]$: any gap closing makes the Chern numbers from Eq.~(14) ill-defined and falsifies the generalized bulk-defect correspondence as stated. A complementary check is to measure the boson number transported between two pinned solitons after one full adiabatic cycle; it must be an integer equal to the Chern-number difference, and a non-integer result would signal gap closure or non-adiabaticity.","tokens_in":26098,"feed_emoji":"⚛️","tokens_out":13390,"duration_ms":136536,"temperature":0.7,"pith_summary":"This paper argues that the $\\mathbb{Z}_2$ Bose-Hubbard model---one-dimensional interacting bosons whose tunnelling is dressed by Ising spins on the bonds---spontaneously forms topological solitons in its ground state when doped away from special fillings, and that each added boson splits into a quantized fraction: $1/2$ above or below half filling, and $\\pm 1/3$ around two-third filling. The solitons are dynamical objects, not fixed backgrounds: they move unless pinned, repel each other, and at sufficient density order into a fractional soliton lattice. The deeper claim is that the bound fractional bosons are topologically protected, even when a soliton joins two symmetry-broken sectors with the same equilibrium Berry phase. That robustness is established through a generalized bulk-defect correspondence built on quantized inter-soliton Thouless pumping, in which the pumping phase acts as a synthetic dimension and Chern-number differences across a defect dictate the pumped charge.","feed_headline":"Ground-state solitons split extra bosons into halves and thirds","feed_subtitle":"The Z2 Bose-Hubbard model's pumped Chern numbers show why these fractional bound states are protected.","key_machinery":"The load-bearing machinery is the inter-soliton Thouless pumping construction: the staggered field couplings are modulated as $\\Delta_{\\phi,i}=2(-1)^i\\delta\\cos\\phi$ and $\\beta_{\\phi,i}=(-1)^i\\delta\\sin\\phi$ (Eq.~13), so that $\\phi$ acts as momentum along a synthetic dimension, and each symmetry-broken sector is assigned a many-body Chern number $\\nu$ from the flow of the Berry phase $\\gamma(\\phi)$ through $\\nu=\\frac{1}{2\\pi}\\int_0^{2\\pi}d\\phi\\,\\partial_\\phi\\gamma(\\phi)$ (Eq.~14). The Chern-number difference between the sectors joined by a soliton determines how many bosons are pumped between solitons over a cycle, and the presence or absence of spectral flow decides whether a bound mode is topologically protected. Complementing this, the soliton order parameter takes the universal $\\tanh$ kink profile and the bound boson density the $\\mathrm{sech}^2$ zero-mode profile of relativistic field theory, connecting the numerical ground states to established soliton physics.","core_discovery":"The paper's central discovery is that the $\\mathbb{Z}_2$ Bose-Hubbard model has doped ground states in which the $\\mathbb{Z}_2$ fields spontaneously form solitons that interpolate between symmetry-broken bond-order sectors. Doping one boson above or below half filling creates a soliton-antisoliton pair, and the integrated local density changes by $1/2$ at each defect; around two-third filling, an added particle or hole creates three defects with $1/3$ or $-1/3$ bound bosons. The paper establishes that these are not polarons: the energy of the pair decreases monotonically with separation, so the fractionalized pair is the ground state, and at higher doping the repelling solitons crystallize into a fractional soliton lattice. It then shows that a slow cyclic modulation of the staggered couplings, interpreted as a Thouless pump, assigns each symmetry-broken sector a many-body Chern number through the integrated Berry phase; the difference of Chern numbers across a defect fixes how many bosons are pumped between solitons. This generalized bulk-defect correspondence explains the bound fractional bosons even for $\\mathbb{Z}_4$ solitons that separate sectors with identical equilibrium Berry phases, identifying them as remnants of protected edge states of a synthetic two-dimensional system.","pith_inferences":["The construction explicitly predicts that at other commensurate fillings, such as one-third filling, solitons should also have Chern-number differences (0, $\\pm 2$, and so on) and hence quantized pumped charges even when all sectors share the same Berry phase; this is a testable extension the paper does not carry out.","Because the transport calculations pin the solitons and identify the ground state with the adiabatic evolution, a direct time-dependent simulation of unpinned defects during the cycle would test the paper's claim that the pumped charge remains quantized when the defects move.","The statistical-interaction parameters the paper mentions as an outlook, $g=1/2$ at half filling and $g=1/3$ at two-third filling, imply fractional exclusion statistics with measurable thermodynamic signatures in a cold-atom realization of the model.","The absence of polarons, unlike fermionic SSH models, suggests that either bosonic statistics or the dynamical nature of the $\\mathbb{Z}_2$ field is what selects fractionalization over polaron formation; a bosonic SSH chain without spins would isolate the responsible ingredient."],"forward_implications":["Above or below half filling, a single added boson does not form a polaron; the ground state contains a soliton-antisoliton pair and the integrated density jumps by exactly $1/2$ at each defect.","Around two-third filling, adding or removing one boson creates three defects with $\\pm 1/3$ bound bosons, even though the defects separate sectors that have the same equilibrium Berry phase.","Solitons repel with an exponentially small energy, so a finite density of extra bosons self-assembles into a fractional soliton lattice whose boson density is a periodic array of fractional charges.","The Chern-number difference across a defect, computed through the many-body Berry-phase flow, fixes both the number of bound modes and the direction of inter-soliton transport, generalizing the bulk-defect correspondence to cases where one-dimensional invariants alone cannot explain the bound states.","In the hard-core limit the fractional mode has support on a single sublattice due to chiral symmetry, while at finite $U$ the inversion-symmetry-quantized Berry phase and the pumping argument provide the topological explanation."],"supporting_citations":[{"why":"It supplies the prior result on dynamical solitons and boson fractionalization near half filling that this paper extends to other densities.","marker":"[1]"},{"why":"It provides the relativistic-field-theory zero-mode fractionalization mechanism whose $\\tanh$/ $\\mathrm{sech}^2$ profiles the paper reproduces for bosons.","marker":"[17]"},{"why":"It introduces the Su-Schrieffer-Heeger model that supplies the soliton and bond-order-wave paradigm being bosonized here.","marker":"[19]"},{"why":"It establishes the symmetry-protected bulk-defect correspondence that the paper generalizes through pumping.","marker":"[27]"},{"why":"It defines the $\\mathbb{Z}_2$ Bose-Hubbard model and its Peierls-type phase diagram including the bond-ordered wave.","marker":"[30]"},{"why":"It shows that the half-filled topological bond-ordered wave is a symmetry-protected phase with quantized Berry phase.","marker":"[31]"},{"why":"It identifies the trivial and topological threefold-degenerate phases at one-third and two-third fillings that produce the $\\mathbb{Z}_4$ solitons.","marker":"[32]"},{"why":"It supplies the real-space local Berry phase used to assign quantized invariants to each symmetry-broken sector.","marker":"[63]"},{"why":"It supplies Thouless pumping, the quantization mechanism used to build the synthetic-dimension Chern numbers and the generalized correspondence.","marker":"[68]"}],"fun_headline_variants":["Doping Z2 Bose-Hubbard chain splits bosons into 1/2 and 1/3 solitons","Fractional solitons appear in Z2 Bose-Hubbard model at fractional fillings","Dynamical solitons with 1/2 and 1/3 charges from topological phase intertwining","Pumping reveals why fractional solitons are topologically protected in Z2 model","Spontaneous solitons split bosons into fractional charges in strongly correlated chain"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The generalized bulk-defect correspondence rests on the assumption that the many-body gap of the interacting pumping Hamiltonian stays open for every value of $\\phi$ with $\\delta=t$, so that the Berry-phase integral yields well-defined quantized Chern numbers; the paper computes the resulting Berry phases but does not display the gap spectrum along the cycle.","fun_headline_variants_meta":{"raw":{"variants":["Doping Z2 Bose-Hubbard chain splits bosons into 1/2 and 1/3 solitons","Fractional solitons appear in Z2 Bose-Hubbard model at fractional fillings","Dynamical solitons with 1/2 and 1/3 charges from topological phase intertwining","Pumping reveals why fractional solitons are topologically protected in Z2 model","Spontaneous solitons split bosons into fractional charges in strongly correlated chain"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000795,"raw_usage":{"total_tokens":3571,"prompt_tokens":1089,"completion_tokens":2482,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":705,"completion_tokens_details":{"reasoning_tokens":2360}},"tokens_in":705,"tokens_out":2482,"duration_ms":17023,"temperature":1.0,"reasoning_tokens":2360,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:52:10.609709+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute, with a matrix-product-state algorithm, the lowest two many-body energies of Hamiltonian (13) for $\\delta=t$ and the parameters used for the two-third-filling solitons, scanning $\\phi\\in[0,2\\pi]$: any gap closing makes the Chern numbers from Eq.~(14) ill-defined and falsifies the generalized bulk-defect correspondence as stated. A complementary check is to measure the boson number transported between two pinned solitons after one full adiabatic cycle; it must be an integer equal to the Chern-number difference, and a non-integer result would signal gap closure or non-adiabaticity.","supporting_citations":[{"cited_title":"As described above, there is a bosonic Peierls’ transition where the Ising spins de- velop an antiferromagnetic N´eel-type order [see Fig","cited_arxiv_id":null,"evidence_quote":"It supplies the prior result on dynamical solitons and boson fractionalization near half filling that this paper extends to other densities."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It provides the relativistic-field-theory zero-mode fractionalization mechanism whose $\\tanh$/ $\\mathrm{sech}^2$ profiles the paper reproduces for bosons."},{"cited_title":"Hsieh, D","cited_arxiv_id":null,"evidence_quote":"It defines the $\\mathbb{Z}_2$ Bose-Hubbard model and its Peierls-type phase diagram including the bond-ordered wave."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It shows that the half-filled topological bond-ordered wave is a symmetry-protected phase with quantized Berry phase."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It identifies the trivial and topological threefold-degenerate phases at one-third and two-third fillings that produce the $\\mathbb{Z}_4$ solitons."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It supplies Thouless pumping, the quantization mechanism used to build the synthetic-dimension Chern numbers and the generalized correspondence."}],"review_version":1}