{"id":"3cef2cff-8278-418d-bc13-800c2edbc16d","arxiv_id":"1909.01097","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Rydberg Composites, a Rydberg atom coupled to many lattice atoms, show derived scaling laws and a two-dimensional band structure, connecting few-body trilobite molecules to dense matter.","lead":"One giant Rydberg atom is placed inside a dense lattice of ordinary atoms, and the surrounding atoms sculpt the electron's wave function into new shapes with energy bands. The paper derives how those states scale with atom size and shows that a flat two-dimensional sheet gives the richest, most controllable structure.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The analytic 2D band structure and the headline ν^{-11/2} scaling rest on the unproven continuum replacement of the discrete lattice (Eq. 10); the heuristic dc criterion does not by itself guarantee the discrete spectrum converges to the continuum result at the densities where bands are predicted.","rationale":"The reader's weakest_assumption identified exactly the same load-bearing concern that I would flag: the coarse-graining replacement of the discrete lattice by a homogeneous density, introduced in Sec. IV and formalized in Eq. 10. This is the bridge that converts the discrete many-center Hamiltonian (Eq. 4) into the solvable integral model that produces the analytic band dispersion (Eq. 16), the band-edge energies (Eq. 17), the ν^{-11/2} scaling, and the crossover to ν^{-6}. If this continuum replacement fails at the finite lattice spacings where the bands are claimed to appear, the central quantitative content of the paper is modified. The paper explicitly labels the replacement an assumption and supplies only the heuristic dc criterion of Eq. 31, which is a wave-function-resolution estimate rather than a controlled expansion of the difference between the lattice sum and the integral. The numerical support for the coarse-graining is present but limited: Fig. 2 shows convergence of the density of states for ν=30 and Fig. 5 shows the qualitative Pm transition, but neither verifies the analytic band-edge energies or the scaling exponents directly against the discrete model for multiple ν. I therefore agree with the reader that this is the key soft spot. The single-ν manifold truncation is also a model assumption, but it is explicitly stated and argued to improve with increasing ν; the central band-structure claim is about the model itself, whereas the continuum replacement is the mathematical step that produces the claimed analytic results. For this reason the continuum coarse-graining is the most load-bearing concern. The reader's conditional verdict already accounts for this issue, so I recommend no change to the verdict.","tokens_in":22890,"tokens_out":26490,"duration_ms":258264,"concrete_test":"Diagonalize the full discrete-lattice matrix (Eq. 4 with Eq. 3b) for a square lattice at d/ν = 0.25 and 0.5 for ν = 30, 50, 80, and 100, including all scatterers within r = 2ν^2. Compare the lowest few band-edge eigenvalues of each β band with the continuum eigenvalues of Eq. 14. If the root-mean-square deviation δ = sqrt((1/N) Σ (E_disc - E_cont)^2) does not decrease faster than (d/ν)^2, or if the fitted power-law exponent in ν for the lowest band edge deviates from 11/2 by more than 0.1, the continuum-derived scaling law is not a quantitatively faithful description of the discrete composite in the band regime.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The 2D band structure and the headline scaling laws (ν^{-11/2} for the lower bands, ν^{-6} in the overlap region) are derived entirely in the continuum limit defined by Eq. 10, where the discrete lattice sum of Eq. 4 is replaced by an integral over a uniform plane. The paper states this replacement is 'taken as fact' (Sec. IV) and justifies it with the heuristic dc criterion of Eq. 31, which is based on the angular node spacing of the unperturbed wave function at the inner turning point rather than derived from the full discrete Hamiltonian. The m-block diagonal structure of Eq. 14, and therefore the non-degenerate band structure of Fig. 6, depends on this continuum isotropy; a finite lattice spacing breaks continuous rotational symmetry via lattice harmonics. If the discrete lattice remains partially resolvable at the densities where the bands are predicted (d/ν between roughly 0.5 and 3.5 in Fig. 5), the band energies and the extracted power-law exponents would be modified. The numerical evidence in Figs. 2 and 5 is for a single value ν=30 and validates the DoS, not the analytic band-edge eigenvalues or scaling exponents against the discrete model across ν. Hence the central quantitative claims rest on an unquantified coarse-graining step.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript introduces a model of a Rydberg atom coupled through a Fermi pseudopotential to many ground-state scatterers arranged on one-, two-, and three-dimensional lattices. Working in a single-ν-manifold basis, it identifies the number of shifted states and characteristic lattice spacings for each dimensionality, and it numerically computes densities of states and wave functions. In the homogeneous-density limit the paper derives an analytic 2D band structure with energy scaling ν^{-11/2} for the lower bands and ν^{-6} in the overlap region, together with degenerate ν^{-5} and ν^{-6} shifts in 1D and 3D. It also constructs a supposedly ν-independent 'universal' density of states for the 2D case and uses the adjacent gap ratio to argue that partially filled lattices display quantum chaos with GOE-like statistics once symmetry blocks are accounted for.","tokens_in":23080,"tokens_out":8526,"duration_ms":92833,"significance":"If the central claims hold, the paper provides a valuable systematic bridge between few-scatterer trilobite physics and dense atomic environments, with explicit, analytically transparent scaling predictions that are in principle testable in ultracold-atom experiments. The algebraic derivation of the band-edge exponents, the m-block reduction, and the symmetric C4v decomposition of the spectral statistics are concrete strengths. The RMT analysis is careful in using an unfolding-free diagnostic and in modeling the symmetry-blocked GOE reference, which strengthens the evidence for chaos in the disordered regime.","major_comments":[{"comment":"The analytic band structure and the headline scaling exponents ν^{-11/2} and ν^{-6} are obtained in the continuum limit in which the discrete lattice sum of Eq. (4) is replaced by the uniform-plane integral of Eq. (10). The text explicitly takes this coarse-graining 'as fact,' and the rotational symmetry that produces the m-block-diagonal structure of Eq. (14), and hence the non-degenerate bands of Fig. 6, is a property of the continuum limit rather than of the discrete lattice. The critical-spacing criterion in Eq. (31) is heuristic: it is based on the angular node spacing of the unperturbed wave function at the inner turning point and is not derived from the full discrete Hamiltonian. The numerical evidence in Figs. 2 and 5 is for ν=30 and validates the DoS and Pm structure, but it does not test the band-edge eigenvalues or the scaling exponents against diagonalization of the discrete Eq. (4) as a function of ν. I request a direct convergence study comparing the discrete-lattice spectrum with the continuum eigenvalues of Eq. (14) for several ν values over the d/ν range where the bands are predicted, and a report of whether the exponents ν^{-11/2} and ν^{-6} are recovered before rotational symmetry is broken by lattice harmonics. Without this, the central quantitative claims rest on an unquantified assumption.","section":"Sec. IV B 4 and Appendix D, Eq. (D5)"},{"comment":"The 'universal' density of states in Fig. 8(c) is constructed with an interpolating tanh function whose parameters x0=-0.011 and w=0.0028 are fit to the numerical spectra, and with a compression factor b=0.1 that the text itself describes as 'somewhat arbitrary.' The integrated normalizations in Eqs. (D3) and (D4) also differ by factors of b^2. As presented, this is an empirical collapse and not a parameter-free universal prediction. The main text should state this distinction explicitly, and the stability of the collapse with respect to b, x0, and w (for instance over the range of ν shown in Fig. 8) should be quantified.","section":"Sec. IV B 4 and Appendix D, Eq. (D5)"}],"minor_comments":[{"comment":"In the sentence beginning 'We will devote much of the reminder of this paper,' 'reminder' should be 'remainder.'","section":"Sec. III C"},{"comment":"Equation (31) uses the magnetic quantum number m in sin(π/m), which gives a negative spacing for negative m; the formula should use |m|, and the approximate relation l=2(β-1)-m in the text should be written with |m| so that it remains valid for the full spectrum in Fig. 6.","section":"Sec. V, Eq. (31)"},{"comment":"The symbol b is used both for the band index in Sec. IV A and for the compression factor in Appendix D; this overloading makes the formulas in Appendix D confusing and should be resolved by renaming one of them.","section":"Appendix D, Eq. (D1)"},{"comment":"The parameter g is introduced without an explicit definition in the main text; the reader is left to infer from context that g=1 in the band region and g=1/2 in the overlap region, and this should be stated.","section":"Appendix D, Eq. (D1)"},{"comment":"The caption of Fig. 8 reports a FWHM of approximately 0.0235, but the corresponding Gaussian width σ used in Eq. (D1) is not stated for panel (c); reporting σ would improve reproducibility.","section":"Fig. 8"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is well within the journal's scope and the citation practice appears appropriate. My main concern is the continuum replacement behind the analytic scaling laws; I am recommending major revision rather than rejection because the required additional evidence is a computational convergence study that is natural to carry out within the manuscript's existing framework."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. First, this paper actually delivers a new synthesis: it takes a single Rydberg atom in a dense scatterer environment, treats 1D, 2D, and 3D lattices, and derives analytic scaling laws (ν^-5, ν^-11/2, ν^-6) plus a non-degenerate band structure in the 2D homogeneous limit. Second, the 2D band structure and those scaling exponents rest entirely on a continuum replacement of the discrete lattice (Eq. 10), which the authors say they 'take as fact.' The heuristic dc criterion doesn't bridge that gap, and the numerical diagonalization is only shown for ν=30. That makes the central quantitative claims conditional, not wrong.\n\nWhat's genuinely new is the systematic interpolation between few-scatterer trilobite molecules and dense homogeneous environments, and the geometry-dependent classification: degenerate shifts in 1D and 3D, bands in 2D. The analytic diagonal approximation leading to the scaling laws is clean, and the appendices are useful—the trilobite-basis integral equation, parabolic coordinates, and symmetry-adapted orbitals are real technical contributions. The paper is honest about its soft spots: it flags the coarse-graining assumption, admits the compression factor in the universal DoS is 'somewhat arbitrary,' and uses an auxiliary block-diagonal GOE model for the filled-lattice chaos rather than claiming raw GOE.\n\nThe soft spots are in proportion. First, the coarse-graining: if the Rydberg wave function can still resolve lattice structure at the densities where bands are predicted, the isotropy that produces the m-block structure is broken, and band energies and exponents would shift. A multi-ν check of discrete band edges against the analytic continuum result would settle this. It is a missing validation, not a disproven claim. Second, the 'universal' DoS depends on fitted tanh parameters and an arbitrary compression factor; the universality is partly a fit. Third, the quantum-chaos classification for the filled lattice is symmetry-consistent rather than direct: the AGR matches a block-diagonal GOE, so it's evidence of chaotic dynamics within symmetry sectors, not a raw GOE statement.\n\nWho gets value: ultracold Rydberg experimentalists, optical lattice people, and the Rydberg molecule theory community. The experimental section is realistic about lattice spacing and ν requirements.\n\nRecommendation: send it to a serious referee. The framework is important, the analytic core is sound within its stated model, and the main weakness is a validation gap rather than a fatal flaw. A good referee can press for multi-ν tests and a sharper treatment of the continuum limit.","headline":"Genuinely new framework with clean analytic scaling laws, but the 2D band structure rests on an unproven continuum limit and single-ν numerics; still deserves a serious referee.","tokens_in":23716,"tokens_out":2760,"would_cite":true,"duration_ms":27056,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A flat atomic sheet gives one Rydberg atom a band spectrum","keywords":["Rydberg atoms","Rydberg molecules","optical lattices","ultracold gases","Fermi pseudopotential","band structure","quantum chaos","level statistics"],"falsifier":"Diagonalize the full discrete-lattice matrix (Eq. 4) without any coarse-graining for several $\\nu$ at lattice spacings at or below the critical spacing $d_c$, and compare the eigenvalues with the continuum-integral result of Eq. (14); if the spectra fail to converge to the predicted bands and to the $\\nu^{-11/2}$/ $\\nu^{-6}$ scalings, the homogeneous-density assumption is where the argument breaks down.","tokens_in":22588,"feed_emoji":"⚡️","tokens_out":8899,"duration_ms":87658,"temperature":0.7,"pith_summary":"This paper introduces the Rydberg Composite: a single highly excited Rydberg atom enveloped by many ground-state atoms arranged on a one-, two-, or three-dimensional lattice. The authors aim to show that this object has a systematic spectral theory that connects the familiar few-scatterer “trilobite” molecules to the dense, homogeneous limit, with properties set by the principal quantum number $\\nu$, the lattice spacing $d$, and the fill factor $F$. Their central result is that a two-dimensional monolayer keeps a rich, non-degenerate band structure in the dense limit, with band-edge energies scaling as $\\nu^{-11/2}$ and crossing over to $\\nu^{-6}$, while one- and three-dimensional composites shift uniformly and become degenerate again, with energies $\\nu^{-5}$ and $\\nu^{-6}$. In partially filled, disordered lattices the same Hamiltonian produces chaotic spectra whose level statistics match the Gaussian orthogonal ensemble. The appeal is that the planar geometry turns the lattice into a sculpting tool for the electronic wave function, promising controllable dense Rydberg systems.","feed_headline":"A flat atomic sheet gives one Rydberg atom a band spectrum","feed_subtitle":"The 2D composite keeps a non-degenerate tunable spectrum; wires and crystals only shift the whole level set.","key_machinery":"The central object is the single-manifold pseudopotential matrix $V_{lm,l'm'} = 2\\pi a_s \\sum_i \\langle\\nu lm|R_i\\rangle\\langle R_i|\\nu l'm'\\rangle$, evaluated inside one Rydberg manifold. In the homogeneous 2D limit this becomes an integral over the plane, and the factorization of the integrand into an azimuthal integral, a radial overlap integral, and a projection of spherical harmonics into the plane is what produces tractable structure: the azimuthal part gives $\\delta_{mm'}$, the radial overlap gives near-diagonal dominance, and the spherical-harmonic projection gives the circular-state selectivity that organizes states into $b$-bands. The band index $b$, defined by $l=2b-2+m$, is the identity that carries the argument: states along a fixed $b$ diagonal have nearly identical overlap with the plane, so the problem becomes quasi-one-dimensional per band and yields the analytic dispersion relation $\\tilde E^{2D}_{bk}\\approx (8\\nu+4(b-k)-1)\\Gamma(b-1/2)/(8\\pi^2\\nu^{13/2}\\Gamma(b))$. The same machinery, applied to 1D with the quantization axis along the chain and to 3D with the full volume integral, produces the degenerate shifts of Eqs. (27) and (29).","core_discovery":"At the paper’s core is the claim that the symmetry of the scatterer geometry, not the number of scatterers, determines the dense-limit spectrum. In the homogeneous limit, where the discrete lattice sum is replaced by an integral over a continuous density, the 2D monolayer Hamiltonian becomes block-diagonal in the angular momentum projection $m$ because the azimuthal integral forces $\\delta_{mm'}$. Because only states with $l+m$ even have amplitude in the plane, each $m$-block contains states labeled by band indices $b$ and $k$, with $l=2b-2+m$ and $k=\\nu-l$, and the dominant diagonal elements produce a dispersion-like spectrum. The paper derives analytically that the band lower edges scale as $\\nu^{-11/2}$, the level spacing within a band as $b^{-1/2}\\nu^{-13/2}$, and that the scaling crosses to $\\nu^{-6}$ when bands overlap for large $b$; a universal density of states is obtained by interpolating between the two regimes with a hyperbolic tangent. For 1D and 3D the same homogeneous replacement leaves all states degenerate: in 1D a radial selection rule makes every state shift identically ($\\nu^{-5}$), and in 3D the normalization integral gives a uniform shift ($\\nu^{-6}$). For disordered partial filling, the adjacent gap ratio approaches the Gaussian orthogonal ensemble value $0.530$, indicating quantum chaos in the spectral fluctuations.","pith_inferences":["The paper does not explore it, but the dispersion-like dependence of energy on $m$ within each band suggests that patterned lattice densities could encode a desired spectral function, effectively using the lattice as a template for the composite’s level structure.","A direct testable extension is to measure the adjacent-gap ratio across the predicted critical spacing $d_c$ for several $\\nu$ values: Eq. (31) predicts a crossover from chaotic to symmetry-dominated statistics whose location should scale as $\\nu^2$, a scaling the paper derives heuristically but does not verify numerically.","The same coarse-graining technique could be applied to other central potentials whose wave functions have a planar node structure, such as excitonic systems, as long as the scatterer distribution has a symmetry plane; this would generalize the band-formation mechanism beyond atomic Rydberg physics."],"forward_implications":["A monolayer Rydberg Composite at lattice spacings below the critical value has a predictable, non-degenerate sequence of energy levels separated by a band gap, in contrast to dense 1D and 3D cases where all levels collapse onto one shifted energy.","Dense one- and three-dimensional composites give a single uniform shift of the whole Rydberg manifold, scaling as $\\nu^{-5}$ or $\\nu^{-6}$, so the spectrum itself stays simple while the shift encodes the scatterer density.","At partial filling, spectral fluctuations follow Gaussian orthogonal ensemble statistics, so a disordered 2D lattice offers a tunable laboratory for quantum chaos with $\\nu$, $d$, and $F$ as control parameters.","At sufficiently high $\\nu$ the composite manifold is energetically isolated from neighboring Rydberg manifolds, because its level shifts fall faster than the $\\nu^{-3}$ inter-manifold spacing, making the single-manifold description self-consistent.","The known few-scatterer trilobite states and the dense homogeneous limit are connected by one Hamiltonian, so the same formalism interpolates between molecular and condensed-matter treatments of Rydberg matter."],"supporting_citations":[{"why":"Introduces the pseudopotential contact interaction between the Rydberg electron and a ground-state atom, which becomes the coupling term in the Hamiltonian (Eq. 2).","marker":"[1]"},{"why":"Supplies the trilobite-basis representation and molecular point-group projection method used to classify states under the $C_{4v}$ symmetry and separate real level crossings.","marker":"[19]"},{"why":"Provides the hydrogen radial matrix element whose selection rule makes the 1D homogeneous shift identical for all states, underpinning the $\\nu^{-5}$ degeneracy claim.","marker":"[33]"},{"why":"Gives the exact adjacent-gap-ratio distributions for Gaussian orthogonal and Poisson ensembles, the reference values used to identify quantum chaos in the disordered regime.","marker":"[35]"},{"why":"Introduces the adjacent-gap-ratio statistic adopted in the paper as an unfolding-free measure of level repulsion.","marker":"[37]"}],"fun_headline_variants":["2D Rydberg composite turns a single atom into a band structure","Symmetry, not scatterer count, sets Rydberg composite's spectrum","Flat array gives one Rydberg atom a tunable band spectrum","Monolayer Rydberg composite: bands from a single excited atom"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that below a critical lattice spacing the discrete lattice can be exchanged for a continuous, uniform density; if the Rydberg wave function still resolves individual scatterers at those densities, the band structure and the $\\nu$ scaling exponents derived from that continuum replacement would change.","fun_headline_variants_meta":{"raw":{"variants":["2D Rydberg composite turns a single atom into a band structure","Symmetry, not scatterer count, sets Rydberg composite's spectrum","Flat array gives one Rydberg atom a tunable band spectrum","Monolayer Rydberg composite: bands from a single excited atom"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000337,"raw_usage":{"total_tokens":1900,"prompt_tokens":1014,"completion_tokens":886,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":630,"completion_tokens_details":{"reasoning_tokens":807}},"tokens_in":630,"tokens_out":886,"duration_ms":8546,"temperature":1.0,"reasoning_tokens":807,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:27:23.174621+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Diagonalize the full discrete-lattice matrix (Eq. 4) without any coarse-graining for several $\\nu$ at lattice spacings at or below the critical spacing $d_c$, and compare the eigenvalues with the continuum-integral result of Eq. (14); if the spectra fail to converge to the predicted bands and to the $\\nu^{-11/2}$/ $\\nu^{-6}$ scalings, the homogeneous-density assumption is where the argument breaks down.","supporting_citations":[{"cited_title":"trilobite","cited_arxiv_id":null,"evidence_quote":"Introduces the pseudopotential contact interaction between the Rydberg electron and a ground-state atom, which becomes the coupling term in the Hamiltonian (Eq. 2)."},{"cited_title":"Electric field control in ultralong-range tri- atomic polar Rydberg molecules,","cited_arxiv_id":null,"evidence_quote":"Supplies the trilobite-basis representation and molecular point-group projection method used to classify states under the $C_{4v}$ symmetry and separate real level crossings."},{"cited_title":"Wave functions with localizations on classical periodic orbits in weakly perturbed quantum billiards,","cited_arxiv_id":null,"evidence_quote":"Provides the hydrogen radial matrix element whose selection rule makes the 1D homogeneous shift identical for all states, underpinning the $\\nu^{-5}$ degeneracy claim."},{"cited_title":"The application of general zero-range po- tentials to multi-center problems,","cited_arxiv_id":null,"evidence_quote":"Gives the exact adjacent-gap-ratio distributions for Gaussian orthogonal and Poisson ensembles, the reference values used to identify quantum chaos in the disordered regime."},{"cited_title":"Bose-Einstein condensate in a harmonic trap decorated with Diracδ functions,","cited_arxiv_id":null,"evidence_quote":"Introduces the adjacent-gap-ratio statistic adopted in the paper as an unfolding-free measure of level repulsion."}],"review_version":1}