{"id":"0c6653e2-e737-4244-bd57-c6e9ed68fad8","arxiv_id":"1908.05260","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A continuum dielectric model reproduces the full coupled-dipole photon scattering lineshape of a cold Gaussian atom cloud for large atom numbers, until dense packing induces dipole-dipole correlations.","lead":"This paper compares two ways of calculating how a cold cloud of atoms scatters light: treating each atom as a tiny antenna, or treating the cloud as a smooth glass-like material. It finds that the smooth model works surprisingly well for large clouds, and fails when atoms are packed densely enough to interact strongly with each other.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"For the largest clouds (N=2^17), the point-dipole average uses only 4 configurations and no error bars are shown, so the claimed 'indistinguishable' agreement with the continuum model may not be robust.","rationale":"The paper's core claim is a numerical comparison between a full point-dipole simulation and a continuum dielectric model. The full point-dipole results are presented as ensemble-averaged curves, but the ensemble size for the largest N is tiny (4 samples), and no statistical uncertainty is reported. Since the paper highlights that the smallest-N case (N=2^11, with 256 samples) already shows deviations in forward scattering, the credibility of the 'indistinguishable' agreement at N=2^17 hinges entirely on the reliability of that average. This is exactly the weakness the reader identified. An independent check with more configurations and explicit error bars would settle whether the agreement is quantitative or an artifact of insufficient sampling. No issue with the internal derivations or the iterative method was found that would supersede this concern, so the appropriate recommendation remains CONDITIONAL.","tokens_in":13988,"tokens_out":5865,"duration_ms":57182,"concrete_test":"Repeat the N=2^17, b0=40, ξ=1 total-scattering calculation for at least 20 independent random configurations, compute the mean and standard error of γ(Δ) at a set of detunings (e.g., Δ = -9Γ, -6Γ, ..., 9Γ), and overlay error bars on Fig. 3 together with the paraxial continuum curve. If the standard error is larger than a few percent of γ or the mean shifts by an amount visible at the line-width scale of the figure, the 'indistinguishable' agreement claim is not established; if the error bars are negligible, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Sec. III A the paper states, 'The calculations were averaged over many runs until a total of 2^19 atoms were included in the calculation.' For N=2^17, this means only 4 independent random configurations contribute to the average, while the continuum paraxial curves in Figs. 3, 4, 7, and 8 are plotted without any error bars or convergence metric. The central claim rests on the assertion that the continuum model quantitatively reproduces the full point-dipole result for large N, up to b0=40. Since the paper itself reports a discernible difference for N=2^11 (which has 256 configurations) in the forward scattering (Fig. 8), the large-N comparison is the critical supporting evidence. If the 4-configuration average has significant configuration-to-configuration fluctuations, the displayed curves could shift, potentially changing the lineshape features (e.g., the plateau and the |Δ| dependence) and the inferred regime of validity. Without a variance estimate, the agreement is not fully quantified.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper compares two theoretical descriptions of weakly driven photon scattering from a frozen, Gaussian-distributed cloud of two-level atoms: a microscopic point-dipole model solved numerically for up to 2^17 atoms, and a continuum dielectric model solved in the eikonal or paraxial approximation. The authors identify three polarization regimes as the on-resonance optical depth b0 increases: nearly uniform polarization for dilute clouds, front-localized absorption for intermediate densities, and back-of-cloud refocusing for red detuning at large b0. They find that the continuum paraxial model quantitatively reproduces the total and forward scattering rates of the point-dipole model without adjustable parameters for b0 up to 40, that the eikonal approximation works in a lower-density regime, and that the Clausius-Mossotti local-field correction does not improve and often worsens agreement. They also present an iterative numerical method for solving the point-dipole equations with many atoms.","tokens_in":14198,"tokens_out":3068,"duration_ms":34820,"significance":"If the claimed agreement is robust, the paper is a valuable benchmark: it shows that a simple, parameter-free continuum dielectric calculation can capture non-Lorentzian lineshapes and forward-scattering plateaus in dense cold-atom clouds, and it clearly delineates the regime where mean-field local-field corrections fail because of atom-atom correlations. The paper also provides a concrete iterative algorithm for point-dipole simulations with more than 10^5 atoms and gives a systematic derivation of higher-order paraxial corrections. The false-prediction potential is high, because the continuum model is compared directly to the microscopic equations with nothing fit to make agreement happen.","major_comments":[{"comment":"The central quantitative claim—that the continuum paraxial model is indistinguishable from the point-dipole result for N=2^17—rests on disorder averages of only 4 configurations, since the text states the total number of atoms included across runs is 2^19. No error bars, standard deviations, or convergence metric are shown for these averages, so the visual claim of agreement is not quantitatively supported. Please add configuration-to-configuration variance estimates, or demonstrate convergence by increasing the number of runs, and state the resulting uncertainty in the plotted lineshapes.","section":"Sec. III A, Figs. 3, 4, 7, 8"},{"comment":"In Fig. 8 the N=2^11 point-dipole result, which is averaged over 256 configurations, shows a noticeable difference from the continuum model at small |Δ| for ξ=2, while the N=2^17 result, averaged over only 4 configurations, is claimed to be in excellent agreement. Because the finite-N comparison is the only one with a reasonably large number of configurations, the text should quantify whether the N=2^11 discrepancy is a genuine finite-density continuum breakdown or partly a statistical fluctuation, and error bars on the N=2^17 curve are needed to support the assertion that the large-N curve is the more reliable comparison.","section":"Sec. III B, Fig. 8"},{"comment":"The forward scattering rate is defined by integrating Eq. (4) up to an angle θmax chosen a posteriori so that the differential rate has decreased by two orders of magnitude from its maximum. This is a reasonable convention, but the robustness of the plateau and of the model comparison to the exact cutoff value is not demonstrated. Please show that the conclusions are insensitive to the specific choice of θmax, or state the sensitivity explicitly.","section":"Sec. III A, forward scattering definition"}],"minor_comments":[{"comment":"Several captions refer to the total scattering rate as 'Eq. (4)', but Eq. (4) is the differential rate dγ/dΩ; the total rate is Eq. (6). Please correct these cross-references.","section":"Figure captions 2, 3, 6, 7, 9"},{"comment":"The figure caption contains the LaTeX artifact 'resizebox80mm!'; this should be removed.","section":"Sec. III B, Fig. 6 caption"},{"comment":"The word 'poralization' in the concluding paragraph should be 'polarization'.","section":"Sec. IV, Conclusions"},{"comment":"The spelling 'Claussius-Mossotti' appears in the figure caption but the standard spelling 'Clausius-Mossotti' is used elsewhere in the text; please make the spelling consistent.","section":"Sec. III C, Fig. 9"}],"recommendation":"major_revision","confidential_remarks":"The paper is well within scope for a journal on atomic physics and light-matter interaction. The main scientific claim is interesting, but the quantitative comparison needs statistical support before publication, especially for the largest N values where the number of disorder realizations is very small. I would encourage the authors to add error bars or convergence data, and if this is not possible, to temper the 'indistinguishable' language accordingly."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: for large cold-atom clouds, you can often replace a full coupled-dipole simulation with a continuous dielectric using the low-density susceptibility, and this paper maps out exactly where that shortcut works and where it breaks down. That is a genuinely useful result for the community, and the iterative solver that reaches N=2^17 (with a test at 2^19) is a real technical advance.\n\nThe three-polarization-regime classification is new, as is the systematic comparison of eikonal vs. paraxial vs. full dipole across N and ξ. The negative result on Clausius-Mossotti isn't surprising—it matches earlier work—but it's cleanly demonstrated here. The derivations in Secs. II and V are careful, the paraxial higher-order corrections are checked, and there are no hidden fit parameters; the forward-scattering cutoff is a defined observable rather than a tuneable.\n\nThe soft spot is statistical. The largest runs (N=2^17) are averaged over only 4 configurations, and the figures show no error bars or convergence metric. The paper says the continuum and full-dipole curves are 'indistinguishable,' which is a strong claim to rest on four random draws. I don't think this is fatal—the agreement holds across multiple b0 values and for both total and forward scattering, and Fig. 2 shows the total scattering is nearly N-independent, which suggests the observables are self-averaging. But a reader cannot verify the central quantitative claim without some measure of run-to-run variation. A single sentence on the spread of the four configurations would have fixed this. The paper also ships no code or data, though the method section is detailed enough to reimplement.\n\nThis paper is for cold-atom experimentalists who want to model large clouds without brute-force dipole codes, and for theorists working on cooperative scattering. It deserves a serious referee. I'd recommend accepting it for review with a request for error bars or a convergence statement on the large-N averages. The physics is sound; the presentation just needs to let the noise show.","headline":"Solid computational paper with a practical message for cold-atom experimenters; the only real weakness is missing error bars on the 4-configuration large-N averages.","tokens_in":14714,"tokens_out":4347,"would_cite":true,"duration_ms":40394,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper establishes that for sufficiently large Gaussian atomic clouds, the total and forward photon scattering rates from the full coupled point-dipole calculation are quantitatively reproduced by a continuous dielectric with the…","keywords":["light scattering","cold atomic clouds","point dipole model","continuum dielectric","paraxial approximation","eikonal approximation","Clausius-Mossotti susceptibility","superradiance"],"falsifier":"Recompute the total and forward scattering curves at $b_0 = 40$, $N = 2^{17}$ with tens of independent random configurations and plot the run-to-run spread; if the spread around the mean is comparable to or larger than the gap between the point-dipole and continuum curves, or if the mean shifts away from the continuum result, the central claim fails.","tokens_in":13804,"feed_emoji":"⚛️","tokens_out":8429,"duration_ms":82829,"temperature":0.7,"pith_summary":"This paper asks when a cold atomic cloud can be treated as a smooth chunk of dielectric rather than as a collection of point dipoles. For a weakly driven, frozen Gaussian cloud of two-level atoms, the authors argue that once the cloud is large enough, the total and forward scattering rates from the full coupled-dipole calculation are reproduced by Maxwell's equations with the low-density susceptibility $\\chi_e^{(\\mathrm{ld})}$, solved in the eikonal or paraxial approximation, with no adjustable parameters, up to $b_0 = 40$. This matters because it ties a many-body scattering problem to classical optics, and it sharply delimits where that tie holds: the Clausius-Mossotti local-field correction does not help and usually hurts. The paper also charts three polarization regimes, uniform, front-loaded, and refocusing, each with its own valid mean-field description, and shows where mean-field models stop working altogether.","feed_headline":"A smooth dielectric reproduces cold-atom scattering up to b0 = 40","feed_subtitle":"A plain dielectric reproduces atom-cloud scattering; adding the Clausius-Mossotti correction makes it worse.","key_machinery":"The load-bearing object is the continuum electric susceptibility $\\chi_e^{(\\mathrm{ld})}(\\Delta) = \\frac{i\\rho\\sigma/k}{1-2i\\Delta/\\Gamma}$, with $\\sigma = 6\\pi/k^2$, placed into Maxwell's equation for the electric field. The paper solves that equation in the paraxial approximation, writing the field as $\\vec{E} \\simeq \\hat{e}_x e^{ikz} E_0 \\psi_x$, with $i\\partial_z \\psi_x = -\\frac{1}{2k}\\nabla_T^2\\psi_x - \\frac{k}{2}\\chi_e\\psi_x$; dropping the transverse Laplacian gives the eikonal solution $\\psi_x = \\exp\\!\\left(i\\frac{k}{2}\\int_{-\\infty}^{z}\\chi_e\\,dz'\\right)$, which is why the optical depth OD alone controls the lineshape whenever the eikonal approximation holds. The reference calculation is the coupled-dipole system with the point-dipole Green's function, solved by an iterative update that reaches $N = 2^{17}$ atoms, and the comparison is made on total scattering, integrated forward scattering, and small-angle angular scattering.","core_discovery":"The central claim is that the photon scattering lineshape of a cold Gaussian cloud is controlled by the spatial distribution of induced polarization, and that a classical continuum dielectric correctly describes that distribution in three regimes. For dilute clouds the polarization is nearly uniform, and the lineshape is close to the broadened timed-Dicke Lorentzian; for denser clouds near-resonant light is absorbed at the front, producing a non-Lorentzian, $|\\Delta|$-like lineshape that the eikonal approximation captures; and for wavelength-scale clouds red-detuned light refocuses toward the back, requiring the higher-order paraxial approximation. In all three regimes the continuum model with $\\chi_e^{(\\mathrm{ld})}$ reproduces the full point-dipole total and forward scattering up to $b_0 = 40$, so closely that the plotted curves are described as indistinguishable. At higher densities and small $N$, mean-field models fail because dipole-dipole correlations and diffraction perpendicular to the laser set in, and the Clausius-Mossotti susceptibility makes the disagreement worse rather than better.","pith_inferences":["A natural next test, not performed here, is to push the same comparison past $b_0 = 40$ at fixed $N = 2^{17}$; since the iterative solver already handles $N = 2^{18}$ to $2^{19}$, the onset density of mean-field breakdown could be mapped directly.","Because the eikonal regime depends only on OD, the paper's result implies that experiments comparing clouds of different $N$ at the same OD should see identical detuning curves, a clean observable check.","The failure of Clausius-Mossotti suggests that a more promising route to extend the continuum model is a susceptibility that absorbs short-range pair correlations explicitly, rather than a local-field factor.","The polarization inversion toward the back of the cloud at red detuning should be directly visible in spatially resolved imaging, not only in the scattering lineshape."],"forward_implications":["At fixed $b_0$ and $\\xi = 1$, the total scattering lineshape for $N = 2^{17}$ is essentially the same as for $N = 2^{11}$, so in the eikonal regime the system is controlled by optical depth OD rather than by $N$ or density separately.","At large $b_0$ the near-resonance total scattering narrows and develops a $|\\Delta|$-shaped cusp, while the forward scattering rate plateaus near resonance, both because light cannot penetrate the cloud; the continuum model reproduces both effects.","For elongated clouds with $\\xi = 2$ and small $N$, red-detuned light can focus back into the cloud and produce an extra scattering hump at $\\Delta < 0$, which the paraxial continuum calculation captures.","The Clausius-Mossotti susceptibility consistently gives worse agreement than the low-density susceptibility for these stationary-atom clouds, so local-field corrections are not the right way to extend the continuum model."],"supporting_citations":[{"why":"This reference gives the coherent forward broadening linewidth $\\Gamma' = (1+b_0\\xi/8)\\Gamma$ and the back-of-cloud refocusing effect that this paper tests against the continuum model.","marker":"[8]"},{"why":"This reference first showed that Clausius-Mossotti fails for stationary-atom scattering, the baseline this paper confirms and extends to Gaussian clouds.","marker":"[24]"},{"why":"This reference provides the timed-Dicke superradiance model whose broadened line is compared with the continuum result.","marker":"[2]"},{"why":"This reference supplies the collective-scattering context and the smaller-dense-cloud substitution that OD scaling makes quantitative.","marker":"[10]"},{"why":"This reference is the source of the point-dipole Green's function and the Clausius-Mossotti relation used in the continuum model.","marker":"[20]"},{"why":"This reference is the source of the low-density susceptibility formula used in Eq. (8).","marker":"[30]"},{"why":"This reference provides the paraxial expansion used to solve the continuum Maxwell equation beyond the eikonal limit.","marker":"[31]"}],"fun_headline_variants":["Cold atom cloud scattering shows three polarization regimes","Dielectric model reproduces cold-atom scattering up to b0=40","Mean-field models fail for dense cold-atom clouds","Polarization refocusing controls cold-atom lineshape","Clausius-Mossotti worsens cold-atom scattering fit"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The claimed agreement rests on averaging over only four random configurations for the largest clouds, and no error bars or convergence test are reported, so if run-to-run fluctuations are significant the whole quantitative match could shift.","fun_headline_variants_meta":{"raw":{"variants":["Cold atom cloud scattering shows three polarization regimes","Dielectric model reproduces cold-atom scattering up to b0=40","Mean-field models fail for dense cold-atom clouds","Polarization refocusing controls cold-atom lineshape","Clausius-Mossotti worsens cold-atom scattering fit"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000204,"raw_usage":{"total_tokens":1359,"prompt_tokens":885,"completion_tokens":474,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":501,"completion_tokens_details":{"reasoning_tokens":389}},"tokens_in":501,"tokens_out":474,"duration_ms":4585,"temperature":1.0,"reasoning_tokens":389,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:19:07.441351+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the total and forward scattering curves at $b_0 = 40$, $N = 2^{17}$ with tens of independent random configurations and plot the run-to-run spread; if the spread around the mean is comparable to or larger than the gap between the point-dipole and continuum curves, or if the mean shifts away from the continuum result, the central claim fails.","supporting_citations":[{"cited_title":"For all of the cases in Fig","cited_arxiv_id":null,"evidence_quote":"This reference gives the coherent forward broadening linewidth $\\Gamma' = (1+b_0\\xi/8)\\Gamma$ and the back-of-cloud refocusing effect that this paper tests against the continuum model."},{"cited_title":"Frequency shifts in emission and absorption by resonant systems of two-level atoms,","cited_arxiv_id":null,"evidence_quote":"This reference first showed that Clausius-Mossotti fails for stationary-atom scattering, the baseline this paper confirms and extends to Gaussian clouds."},{"cited_title":"Superradiance,","cited_arxiv_id":null,"evidence_quote":"This reference provides the timed-Dicke superradiance model whose broadened line is compared with the continuum result."},{"cited_title":"Coher- ent forward broadening in cold atom clouds,","cited_arxiv_id":null,"evidence_quote":"This reference supplies the collective-scattering context and the smaller-dense-cloud substitution that OD scaling makes quantitative."},{"cited_title":"Observation of the critical regime near anderson local- ization of light,","cited_arxiv_id":null,"evidence_quote":"This reference is the source of the point-dipole Green's function and the Clausius-Mossotti relation used in the continuum model."},{"cited_title":"Coherence effects in scattering order expansion of light by atomic clouds,","cited_arxiv_id":null,"evidence_quote":"This reference is the source of the low-density susceptibility formula used in Eq. (8)."},{"cited_title":"From Maxwell to paraxial wave optics,","cited_arxiv_id":null,"evidence_quote":"This reference provides the paraxial expansion used to solve the continuum Maxwell equation beyond the eikonal limit."}],"review_version":1}