{"id":"aa26f653-6044-4109-998c-da5492c13011","arxiv_id":"2607.29165","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A universal single-parameter law describes the eigenstate entanglement transition when two chaotic subsystems with local conservation laws are coupled, tested against the kicked Bose-Hubbard model.","lead":"This paper shows that two chaotic quantum systems with separate particle-conservation laws, once coupled, develop eigenstate entanglement along one universal curve set by a single parameter combining system size and coupling strength. The result is derived in a random-matrix model, checked against numerics, and applied to a kicked Bose-Hubbard system.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Full-transition formula Eq. (73) is derived only for equal sector dimensions; no test against CQTE with Bose-Hubbard dimensions is provided, leaving single-parameter universality unverified.","rationale":"The reader's conditional verdict already identifies Appendix A approximations (equal sector dimensions, 1/2 probabilities, independence assumptions) as a load-bearing weakness, and my concern is a direct consequence of those approximations being violated for the Bose-Hubbard dimension structure. However, I sharpen it: the paper provides no numerical test of Eq. (73) against the CQTE with non-uniform sector dimensions, which is the minimal setting needed to establish that the full-transition formula is universal in the single parameter √Λ_q. The physical Bose-Hubbard comparison is confounded by the small-κ̃ caveat, so it cannot serve as that test. This does not change the conditional verdict but provides a concrete route to either confirm or falsify the central universality claim. I do not find a more fundamental flaw: the perturbative derivation of λ₁ and the numerical collapse in Fig. 6 are solid, and the parameter-free derivation of Λ_N for the Bose-Hubbard model is a genuine contribution. The missing check is therefore the decisive gap.","tokens_in":22365,"tokens_out":8769,"duration_ms":73654,"concrete_test":"Run the CQTE with the four-site Bose-Hubbard dimension structure (e.g., N=15 or N=24, as in Fig. 6) and compute the numerically averaged S2/S∞2 as a function of √Λ_q using the definition in Eq. (60) with per-sector Λ_{qA} from Eq. (51). Overlay the analytic prediction of Eq. (73) computed with the same per-sector Λ_{qA} and k_{qA}. If the analytic curve deviates from the numerics by more than the statistical uncertainty (using, say, 100+ realizations), then Eq. (73) does not generalize to non-equal sector dimensions and the single-parameter universality of the full transition is not established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Appendix A derives the purity μ_{qA}(Λ_{qA}) and hence Eq. (73) under three explicit assumptions: all sector dimensions and variances are equal, the probabilities that a coupling matrix element connects q_A to q_A+1 versus q_A−1 are each 1/2, and the Schmidt values are independent of the matrix elements. These are exact for the equal-dimension CQTE but are stated to be 'not exact' for Bose-Hubbard-type dimension structures. The paper nevertheless applies Eq. (73) to the kicked Bose-Hubbard model, where k_{N_A} = (N_A+1)(N−N_A+1) varies strongly across sectors. The only validation of Eq. (73) in a non-equal-dimension setting is against the physical Bose-Hubbard model (Fig. 11), which involves additional uncontrolled assumptions—notably the assertion that small-κ̃ sectors contribute negligibly, which the paper admits lacks a sector-resolved analysis. There is no check of Eq. (73) against the CQTE with Bose-Hubbard dimensions. If Eq. (73) fails there, the full-transition curve is not universal in the single parameter √Λ_q; the collapse in Fig. 6 only demonstrates universality of the transition parameter, not of the functional form of the full transition, and the Bose-Hubbard agreement could be accidental. This is the load-bearing gap in the central claim.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies average eigenstate entanglement in bipartite Floquet systems where two local conserved quantities are broken by a coupling into one global conserved quantity. It introduces the 'conserved quantity transition ensemble' (CQTE), built from random matrices in each symmetry sector plus a nearest-sector Gaussian coupling, and derives a single transition parameter Λ. Using regularized perturbation theory, the authors obtain a leading-order law for the largest Schmidt value, λ₁ = 1 − C√Λ, and a non-perturbative interpolation, Eq. (73), for the average linear entanglement entropy S₂ as a weighted sum over sectors. The paper then applies this to a kicked Bose-Hubbard model with L=4 sites, deriving Λ_N analytically from the Hamiltonian. Numerical comparisons against the equal-dimension CQTE and the kicked Bose-Hubbard model are presented.","tokens_in":22582,"tokens_out":7480,"duration_ms":70031,"significance":"If the central claim holds, this is a significant step: it provides a universal, single-parameter description of eigenstate entanglement across the breaking of two U(1) conservation laws into one, and it cleanly separates the number entropy from the configurational entropy. The paper's strengths include an explicit random-matrix ensemble, an analytic derivation of Λ_N (Eq. (86)) directly from the kicked Bose-Hubbard Hamiltonian rather than by fitting, and explicit numerical data in Figs. 8–11 supporting the perturbative law and the universal collapse. The main caveats are that the full-transition formula Eq. (73) is derived under equal-sector assumptions in Appendix A and that its application to the Bose-Hubbard model relies on an unverified sector-resolved dominance argument. These gaps affect the strength of the claimed universality of the full transition, not the perturbative small-Λ result.","major_comments":[{"comment":"The full-transition formula (73) is derived under the explicit assumptions stated in Appendix A: all sector dimensions and variances are equal, the probabilities for coupling q_A→q_A+1 and q_A→q_A−1 are each 1/2, and the Schmidt values are independent of the matrix elements. The paper itself notes these are 'not exact' for Bose-Hubbard-type dimension structures. Nevertheless Eq. (73) is applied to the kicked Bose-Hubbard model via Eq. (86), where k_{N_A}=(N_A+1)(N−N_A+1) varies strongly across sectors. The only validation of Eq. (73) in a non-equal-dimension setting is against the physical Bose-Hubbard model (Fig. 11), which introduces additional uncontrolled approximations. There is no numerical test of Eq. (73) against the CQTE with Bose-Hubbard dimensions. Such a test is necessary to establish that the full transition curve, not just the transition parameter, is universal in √Λ_q alon","section":"Appendix A / Eq. (73); Sec. IV"},{"comment":"The justification for applying the COE-based CQTE to the kicked Bose-Hubbard model is the assertion that 'the dominant contributions to the entanglement entropy come from subspaces with sufficiently large effective κ̃=UN_A/N'. No sector-resolved analysis is given. The authors also state that the fraction of small-κ̃ operators does not diminish as N increases, so this is not a standard large-N suppression. A quantitative check, e.g., comparing the contribution of each N_A sector to S₂ against Λ_{N_A}, or replacing the small-κ̃ U_{N_A} operators by integrable (Poisson) evolution while keeping the large-κ̃ sectors COE, is needed to support the agreement in Fig. 11.","section":"Sec. IV, final paragraph"},{"comment":"The interpolation step leading to Eq. (73) uses the exponent 1/S∞₂ without a derivation from the recursive equation (67)–(69). The exponent is chosen so that the formula approaches the RMT asymptotic value, and the numerical agreement for equal-dimension CQTE in Figs. 9–10 is good. However, because the same interpolation is then exported to the Bose-Hubbard model, the paper should state more clearly that the functional form of the full transition is partly a heuristic interpolation, not a consequence of the perturbation theory. This does not invalidate the small-Λ result, but it bears on the strength of the claim that Eq. (73) is a universal analytic description of the full transition.","section":"Sec. IIIB7 / Eq. (73)"}],"minor_comments":[{"comment":"The typography for the COE coefficients is ambiguous: '4√(2π)' appears where the perturbative coefficient should be 4/√(2π). Please check all COE coefficients in Eqs. (58), (59), (61), and (63) and ensure the division by √(2π) is printed unambiguously.","section":"Eqs. (58)–(59), (61)–(63)"},{"comment":"The notation for the large-M limit is confusing: Eq. (44) writes k̃_{q_A}=k_{q_A}/√M, while Eq. (46) writes D_{q_A±1}=2π/(M k̃_{q_A±1}). Clarify whether M is the total Hilbert-space dimension or the number of blocks, and how k relates to M. This will help readers verify the M-independence of Λ.","section":"Sec. IIB / Eqs. (44)–(46)"},{"comment":"The expression '(2π)^{1/4} 3√{2+2}' is poorly typeset; it should presumably be the cube root of (2+2), i.e., (2π)^{1/4}∛(4). Please reformat all such coefficients in Eq. (71) and Eq. (A21)–(A22) so that the derivation can be checked.","section":"Eq. (71) / App. A"},{"comment":"Figure 11 shows only two system sizes (N=20,24) and 10 disorder realizations, with no error bars. Given that the central claim is quantitative agreement, an estimate of the statistical uncertainty (e.g., with error bars or shaded band) would be helpful.","section":"Sec. IV / Fig. 11"},{"comment":"There is a typo in the sentence 'Them, the normalized distribution of the matrix elements...'; it should read 'Then'. Also, the assumptions of statistical independence used in Eq. (A11) should be stated explicitly at the point where they are first used.","section":"App. A"}],"recommendation":"major_revision","confidential_remarks":"The core perturbative result is well supported and likely correct. The main issue is that the full-transition universality claim, Eq. (73), is pushed beyond the regime in which it is actually validated: it is derived under equal-sector assumptions and applied to Bose-Hubbard dimensions without a direct CQTE test or sector-resolved Bose-Hubbard analysis. I recommend major revision. A straightforward numerical check of Eq. (73) for the CQTE with Bose-Hubbard dimensions, plus the sector-resolved analysis suggested above, would substantially strengthen the paper."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Honest take: this is a real extension of the Tomsovic–Lakshminarayan–Bäcker program, not a repackaging. The CQTE ensemble, the sector-resolved transition parameter Λ_q, the perturbative λ₁ law, and the ab initio Λ_N for the kicked Bose-Hubbard model are genuinely new, and the Bose-Hubbard comparison in Fig. 11 is the strongest part—no fitted constants against the physical system. The number/configurational entropy split is well motivated, and the use of Bianchi–Donà RMT predictions is appropriate.\n\nThe soft spots are real but not fatal to the core claim. The stress-test is the important one: Eq. (73), the “analytic description of the full transition,” is derived in Appendix A under equal sector dimensions, equal variances, 1/2 mixing probabilities, and independence of Schmidt values and matrix elements—assumptions the authors explicitly say are not exact for Bose-Hubbard structures. They then apply Eq. (73) to the kicked Bose-Hubbard model without ever testing it against the CQTE with Bose-Hubbard dimensions. So the collapse in Fig. 6 demonstrates universality of the transition parameter, but not of the functional form of the full transition. The agreement in Fig. 11 could be accidental. This is a genuine gap, and it should be fixable: run the CQTE with the Bose-Hubbard sector dimensions and see if Eq. (73) still holds.\n\nSecond, the abstract overstates the result: Eq. (73) is a heuristic interpolation, and the authors themselves note s(Λ) decreases at large Λ, which is an artifact. That should be softened.\n\nMinor issues: error bars are absent, only 10 disorder realizations, and S_conf’s failure to collapse is explained post hoc. The “manufactured universality” point about the variance scaling is not a real flaw—choosing the scaling to define the right variable is standard RMT practice, and the numerics justify it.\n\nOverall: the central perturbative result and the parameter derivation hold up. The full-transition claim needs a further test. This deserves a serious referee; an editor should send it out and ask for that test and a sector-resolved check of the κ̃-dominance assumption.","headline":"A genuinely new RMT ensemble and a clean Λ_N derivation for kicked Bose-Hubbard, but the full-transition interpolation is untested in unequal sector dimensions—worth a serious referee.","tokens_in":23243,"tokens_out":2781,"would_cite":true,"duration_ms":26803,"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":"For a bipartite system in which two chaotic subsystems each conserve a local quantity and the coupling leaves only the sum conserved, the average eigenstate entanglement is claimed to follow a universal transition governed by one parameter,","keywords":["entanglement transition","symmetry-resolved entanglement","number entropy","conserved quantity","random matrix model","kicked Bose-Hubbard model","eigenstate entanglement","quantum chaos"],"falsifier":"Compute the mean largest Schmidt value and the mean linear entropy as functions of the rescaled coupling for two independent dimension choices of the CQTE and for the kicked Bose-Hubbard model: the paper predicts the specific slope λ₁ ≈ 1 − (4/√(2π))√Λ_q for COE in the small-coupling limit and the universal collapse of S₂/S∞₂ onto Eq. (73); a different slope or a failure of the curves to collapse would refute the central claim.","tokens_in":22042,"feed_emoji":"⚛️","tokens_out":6424,"duration_ms":56657,"temperature":0.7,"pith_summary":"The paper sets out to show that when two quantum-chaotic subsystems, each carrying its own conserved quantity, are coupled so that only the total quantity survives, the average entanglement of the resulting eigenstates is controlled by a single dimensionless parameter rather than by the microscopic details. To make this precise, the authors construct a structured random-matrix ensemble—the conserved-quantity transition ensemble—that encodes the sector structure of the broken conservation law, and they derive both a leading-order perturbative law and an extrapolated analytic curve for the full transition. The symmetry-resolved decomposition splits entanglement into a number-entropy part, which tracks delocalization along the conserved sectors, and a configurational part, which is the genuine subsystem entanglement; in the weak-coupling limit the total is essentially number entropy. Applied to a kicked Bose-Hubbard model, the theory yields an explicit formula for the transition parameter and matches numerically computed eigenstate entropies, supporting the claim of universality.","feed_headline":"Universal entanglement law appears when two conservation laws merge","feed_subtitle":"One number, derived from sector sizes and coupling, predicts the whole transition in the kicked Bose-Hubbard model.","key_machinery":"The key object is the conserved-quantity transition ensemble (CQTE): a statistical model in which the time-evolution operator has the form U = E e^{iV}, with E a diagonal random unitary carrying the uncoupled COE/CUE spectra and V a Gaussian coupling that only connects adjacent sectors of the local conserved charge. The ensemble is designed so that, after a large-dimension scaling, all statistical inputs (level spacings, coupling fluctuations) combine into a sector-dependent parameter Λ_{qA}, whose dimension-weighted average √Λ_q controls the transition. The analytical work combines two tools: standard perturbation theory for the Schmidt values, regularized by assuming immediate avoided cros","core_discovery":"The central claim is that the average linear entanglement entropy S2 of eigenstates in such a symmetry-breaking bipartite system is a universal function of a single transition parameter √Λ_q, defined as the dimension-weighted average of sector transition parameters. In leading order for the time-reversal-invariant (COE) case, the largest Schmidt value obeys λ₁ = 1 − (4/√(2π))√Λ_q, and the extrapolated full transition is given by S₂ = S∞₂ (1 − (1/K) Σ_{qA} k_{qA} ν_{qA}(Λ_{qA})). The same parameter also organizes the number entropy, which dominates the total entropy for weak coupling and exhibits a slow, power-law approach to the random-matrix limit. For the kicked Bose-Hubbard model the tran","pith_inferences":["The same single-parameter description likely extends to the growth of entanglement after a sudden quench, since the eigenstate results and the local Fock-space coupling structure are closely tied; the paper does not treat dynamics, but the machinery appears transferable.","Because the number entropy is measurable through number fluctuations, the predicted Λ_N could act as a probe of the onset of quantum chaos in interacting bosonic systems without needing full state tomography.","The banded-structure localization resembles Fock-space Anderson localization; testing whether the slow power-law approach to the random-matrix limit persists at larger particle numbers would sharpen that analogy and could connect to many-body localization studies.","A natural extension is to non-Abelian conserved charges (for example spin), where the sector decomposition and the nearest-sector coupling assumption would need revision; the paper lists this as a future direction."],"forward_implications":["Any weakly chaotic bipartite system with two local conservation laws merging into one should show the same rescaled entanglement transition, independent of its microscopic realization.","For weak coupling, total entanglement is essentially the number entropy, so measuring particle-number fluctuations in one subsystem approximates the entanglement entropy.","The explicit Λ_N for the kicked Bose-Hubbard model turns the universal curve into a quantitative prediction that can be checked in current cold-atom or optical-lattice setups.","The recursive extrapolation yields a closed-form description of the entire transition, not just the small-coupling limit, allowing predictions in the region where the largest Schmidt value is no longer near unity."],"fun_headline_variants":["One parameter predicts entanglement in kicked Bose-Hubbard model","Universal entanglement transition from merged conservation laws","Symmetry breaking yields universal entanglement law","Single number captures full entanglement transition","Kicked Bose-Hubbard: entanglement set by one parameter"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The derivation of the full transition curve hinges on the assumption that resonant level crossings immediately turn into avoided crossings once the coupling is switched on, together with the sector-averaging approximations (equal sector dimensions and variances, equal probability for ±1 sector coupling, and statistical independence of matrix elements and Schmidt values) used in the extrapolation.","fun_headline_variants_meta":{"raw":{"variants":["One parameter predicts entanglement in kicked Bose-Hubbard model","Universal entanglement transition from merged conservation laws","Symmetry breaking yields universal entanglement law","Single number captures full entanglement transition","Kicked Bose-Hubbard: entanglement set by one parameter"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000146,"raw_usage":{"total_tokens":1005,"prompt_tokens":713,"completion_tokens":292,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":457,"completion_tokens_details":{"reasoning_tokens":223}},"tokens_in":457,"tokens_out":292,"duration_ms":3667,"temperature":1.0,"reasoning_tokens":223,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T12:26:00.946923+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the mean largest Schmidt value and the mean linear entropy as functions of the rescaled coupling for two independent dimension choices of the CQTE and for the kicked Bose-Hubbard model: the paper predicts the specific slope λ₁ ≈ 1 − (4/√(2π))√Λ_q for COE in the small-coupling limit and the universal collapse of S₂/S∞₂ onto Eq. (73); a different slope or a failure of the curves to collapse would refute the central claim.","supporting_citations":[],"review_version":1}