{"id":"138ed7c3-999b-481b-9a34-089347436e65","arxiv_id":"2411.14098","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A weakly driven chiral atom-waveguide array with one undriven atom reaches a steady state in which up to two-thirds of the excitation localizes on that atom as the interatomic spacing approaches zero.","lead":"This paper predicts that removing the laser drive on one atom in a waveguide-coupled atomic array makes steady-state atomic excitation concentrate strongly on that single undriven atom, up to a saturation fraction of 2/3. The result suggests a simple way to create and position single-site excitations, potentially useful for quantum memory and quantum simulation.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Equation (16) silently assumes γ_ng=0; including non-guided decay changes the ξ→2π saturation from 2/3 to ≈(N-1)/N, so the universal 2/3 limit is not robust.","rationale":"The reader's CONDITIONAL verdict is appropriate, but the specific weakest assumption they identify (near-field corrections at ξ < 0.1π) is not the most load-bearing issue. The paper explicitly states its validity for ξ ≥ 0.1π, and the saturation at ξ = 2π corresponds to a physical separation of one guided wavelength, where near-field corrections are suppressed. A more direct, internal problem is that Eq. (16) omits γ_ng, which appears in the model's matrix M in Eq. (7). At the critical point ξ = 2π and γ_ng = 0, M has a zero eigenvalue, so the steady state is singular; the paper even excludes this point in Sec. III A. When γ_ng is nonzero, the exact steady state at ξ = 2π is regular, and the limiting population at the undriven atom approaches (N-1)/[(1+γ_ng/γ)^2 + N-1], which for large N is close to 1 rather than 2/3. Thus the claimed universal saturation value 2/3 depends on taking γ_ng→0 before ξ→2π; for any finite non-guided decay, the value is different. This does not overturn the paper's qualitative conclusion of strong single-site localization — indeed the localization is even stronger with γ_ng > 0 — but it means the headline quantitative claim is not robust. The paper's statement of robustness against γ_ng up to 0.5γ appears to refer to finite ξ (e.g., ξ = 0.1π) and does not rescue the exact 2/3 at the singular limit. The manuscript should state the γ_ng = 0 assumption behind Eq. (16) and clarify the order-of-limits issue, or provide the γ_ng-dependent generalization. With that condition, the work remains a valid contribution, so the reader's CONDITIONAL verdict stands.","tokens_in":15521,"tokens_out":28763,"duration_ms":267122,"concrete_test":"Numerically solve Eq. (8) for the defect-driving scheme with D = 0, normal incidence, N = 100, ξ = 2π, and γ_ng/γ = 0.1 and 0.5, computing P_m from M^{-1} b. If the results match (N-1)/[(1+γ_ng/γ)^2 + N-1] (approximately 0.988 and 0.978) rather than 2/3, the saturation claim fails for nonzero non-guided decay; if they reproduce 2/3, the concern is refuted.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The analytical solution in Eq. (16), and hence the saturation limit lim_{ξ→0(2π)} P_m = 2/3 in Eq. (17), implicitly sets γ_ng = 0 even though γ_ng appears in the model matrix (7). This matters because the limit is singular: at D = 0 and ξ = 2π with γ_ng = 0, M = -(γ/2)J has a zero eigenvalue and no steady state exists in the single-excitation subspace; the paper itself excludes such parameter points in Sec. III A. If the physical non-guided decay γ_ng > 0 is retained, then at ξ = 2π, M = -(γ/2)J - (γ_ng/2)I, and solving M p = b for b_j = -iΩ (j≠m), b_m = 0 yields |p_m|/|p_j| = γ(N-1)/(γ+γ_ng) for all driven j. Consequently, P_m → (N-1)/[(1+γ_ng/γ)^2 + N-1] as ξ→2π, which for N = 100 and γ_ng = 0.5γ is approximately 0.978, not 2/3. The limits ξ→2π and γ_ng→0 therefore do not commute; the universal 2/3 is an artifact of the γ_ng = 0 assumption and is not protected by the claimed robustness to γ_ng up to 0.5γ.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies a one-dimensional array of two-level atoms coupled to a chiral waveguide, driven weakly by lasers from both sides. It introduces two driving schemes: an asymmetric scheme in which the array is split into two angular zones, and a defect scheme in which one or more atoms are left undriven. The paper characterizes steady-state single-excitation populations through numerical solution of the weak-excitation master equation, and reports strongly confined interfaced or single-site localization, an empirical formula for localization-favoring parameters, and an analytic closed-form solution for the defect scheme under reciprocal coupling that yields a purported universal saturation P_m → 2/3 as the interatomic distance tends to 0 or 2π. It closes with proposals for precise multiple-site localization.","tokens_in":15810,"tokens_out":11789,"duration_ms":103504,"significance":"If the central claims hold, the paper offers an analytically tractable route to strong single-site excitation localization in waveguide QED, with the 2/3 saturation as a clean parameter-independent signature and a set of control protocols for multi-site localization. The strengths are the explicit closed-form amplitude solution in Eq. (16), which I verified by substitution for the γ_ng=0, θ1=θ2=π/2, D=0 case, the systematic use of IPR/IIPR/RIEL diagnostics, and the clearly stated multi-site localization proposals. However, the robustness of the 2/3 limit to non-guided decay, the overlap of the optimal parameter with the model's own validity cutoff, and the unsupported generalization to arbitrary angles need to be addressed before the central claim is established.","major_comments":[{"comment":"The analytic solution in Eq. (16) is obtained only for γ_ng=0, but the text does not state this condition. Retaining the full diagonal of Eq. (7), which is -(γ+γ_ng)/2 at D=0, and taking ξ→2π (equivalently ξ→0) gives |p_m|/|p_j| = γ(N-1)/(γ+γ_ng) for a driven atom j≠m, hence P_m → (N-1)/[(1+γ_ng/γ)^2 + N-1]. For N=100 and γ_ng=0.5γ, this is ≈0.978, not 2/3. The limits ξ→2π and γ_ng→0 do not commute; indeed, for γ_ng=0 and ξ=2π the matrix M=-(γ/2)J is singular and no steady state exists in the single-excitation subspace. Therefore the claimed universal 2/3 saturation and its stated robustness to non-guided decay up to 0.5γ cannot both hold as stated. The paper should either explicitly present Eqs. (16)–(17) as the γ_ng=0 result and re-examine the robustness claim at finite ξ, or provide the corresponding closed-form solution with γ_ng>0.","section":"§IV.B, Eqs. (16)–(17); §II, Eq. (7)"},{"comment":"The optimal interatomic distance ξ_max = 2 arctan(1/√(4N-13)) at which P_m is maximal lies below 0.1π for N≳30, while §III.B states that the model is valid only for ξ≥0.1π. For N=100, ξ_max≈0.032π; for N=500, ξ_max≈0.014π. Thus the parameter region where the 2/3 saturation is approached is precisely the region where near-field corrections are neglected. The paper should either extend the model to include near-field terms, or show numerically that a substantial fraction of the 2/3 saturation survives for all ξ≥0.1π.","section":"§III.B and §IV.B, Eq. (18)"},{"comment":"The statement that 'these results hold true for any θ∈(0,π)' is not substantiated. Equation (16) is derived only for the normal-incidence configuration θ1=θ2=π/2, and Figs. 7(a)–(b) show a maximum at θ≈0.86π for ξ=0.1π but do not establish the ξ→0 limit for all θ. Please provide the generalized solution or a systematic (θ,ξ) scan showing that the saturation value is approached as ξ→0 for every θ.","section":"§IV.B, after Eq. (17)"}],"minor_comments":[{"comment":"The factor e^{ikϕ_j} is dimensionally inconsistent because ϕ_j is a phase; it should likely read e^{iϕ_j}, or the phase should be defined as kx_j cosθ_j inside the exponential.","section":"§II, Eq. (6)"},{"comment":"References [16] and [22] cite the same paper by Söllner et al.; the duplicate should be consolidated.","section":"References [16] and [22]"},{"comment":"The range n∈[1-(mξ/π),(mξ/π)-1] is only well-defined when mξ/π is an integer; for general ξ the paper should specify how the integer values of n are selected.","section":"§III.C, Eq. (14)"},{"comment":"The IIPR definition in Eq. (12) is indeterminate when the denominator vanishes (e.g., for a perfectly uniform distribution); the paper should state how such cases are handled.","section":"§III.C, Eq. (13)"},{"comment":"There is a typo in 'interpartical spacing'; it should be 'interparticle spacing'.","section":"§V, Conclusion"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within scope for a quantum-optics journal. The first two major concerns are technical and fixable with additional calculations; the core analytic solution appears correct in its stated simplified regime, but the paper's key claims about universality and robustness need to be narrowed or corrected. I saw no indication of any integrity issue; the empirical formula is disclosed as empirical."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper is a real contribution to the narrow problem of excitation localization in driven waveguide QED. The genuinely new piece is the defect-driving protocol: leave one atom undriven and the steady state localizes on it, with a closed-form solution, Eq. (16), and a clean 2/3 saturation limit. I checked the algebra to Eqs. (17)–(18); it is consistent. The numerical phase diagrams are extensive, and the empirical formula Eq. (14) is a helpful guide, though it is fit to the same data it predicts.\n\nThe main soft spot is the handling of non-guided decay. Eq. (16) comes from setting γ_ng = 0, but the paper doesn't say that. If you keep γ_ng > 0, the ξ→2π limit of the undriven-atom population is not 2/3; it is (N−1)/[(N−1)+(1+γ_ng/γ)^2], which is larger, approaching 1 for large N. So the exact 2/3 is a property of the lossless model, and the limit does not commute with γ_ng→0. The authors claim robustness against γ_ng up to 0.5γ; that is true for the localization phenomenon, but not for the saturation value. They need to state the assumption and report the finite-γ_ng limit explicitly.\n\nA second, smaller caveat: the ξ→0 limit sits at the singular point they exclude in Sec. III A. The limit can be approached with Ω scaled down, so it's not fatal, but the validity boundary should be stated. Eq. (14) being empirical is fine, but it is not derived; readers should treat it as a curve fit.\n\nNone of this kills the central result. The defect-driving scheme works; the localization is strong and analytically explained. The 2/3 saturation should be reframed as the γ_ng=0 limit, with the robust finite-loss limit given. That's a revision, not a rejection. I would send it to peer review. A serious referee will want the γ_ng issue clarified, but the paper deserves the time.","headline":"Defect-driving scheme gives a real single-site localization and a closed-form solution, but the headline 2/3 saturation limit is a lossless artifact; the finite-γ_ng limit is N-dependent and larger.","tokens_in":16347,"tokens_out":8674,"would_cite":true,"duration_ms":80085,"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":"This paper claims that in a weakly driven atomic array coupled to a waveguide, removing the drive on a single atom pins exactly two-thirds of the steady-state excitation onto that atom as the interatomic spacing approaches an integer…","keywords":["atom-waveguide interface","chiral coupling","driven-dissipative system","steady-state localization","defect-driving scheme","inverse participation ratio","dipole-dipole interactions","quantum memory"],"falsifier":"Directly measure or simulate the steady-state population of the undriven atom in a defect-driven reciprocal array at small spacing (e.g., $N=100$, $\\xi=0.05\\pi$) with a full near-field dipole-dipole tensor; if $\\tilde P_m$ deviates from $2/3$, the idealized phase-only coupling is the failing assumption. A less expensive version: compute the first near-field correction to Eq. (17) analytically and check whether the limiting ratio persists.","tokens_in":15318,"feed_emoji":"🎯","tokens_out":6465,"duration_ms":54154,"temperature":0.7,"pith_summary":"The paper considers a one-dimensional array of two-level atoms coupled to a photonic-crystal waveguide and driven weakly by lasers from both sides at different incident angles. It claims that the competition between the photon-mediated dipole-dipole interaction and the phase differences of the drives can concentrate the steady-state excitation onto the atoms at the interface between the two driving zones, sometimes into a soliton-like two-site peak. In the defect-driving variant, where one atom is left undriven, the paper derives an exact analytical steady state showing that, in the limit of vanishing interatomic spacing, the undriven atom carries exactly $2/3$ of the population regardless of the array size. A second empirical formula predicts, as a function of system size and angle difference, which angular configurations favor interfaced localization. The paper argues these effects give a practical way to write, move, and multiply local excitations, with implications for quantum simulation and quantum memory.","feed_headline":"One undriven atom traps two-thirds of the excitation","feed_subtitle":"In a waveguide-coupled atomic array, a single dark site pins the steady-state population, offering a controllable quantum memory.","key_machinery":"The argument runs on the steady-state solution of the linear weak-excitation equations (6)-(8), in which the inverse coupling matrix $M^{-1}$ encodes both the chiral decay rates $\\gamma_L,\\gamma_R$ and the infinite-range photon-mediated dipole-dipole phases $e^{ik|x_j-x_\\ell|}$. For the defect-driving scheme at reciprocal coupling $D=0$ and normal incidence, the paper obtains the exact amplitude vector of Eq. (16), whose special entries at the undriven site and its neighbours produce the saturation formula (17); the limit $\\xi\\to 0$ is controlled by the divergent csc and cot terms, leaving the $2/3$ ratio. For the asymmetric-driving scheme, the empirical rule (14) is read off the RIEL (ratio of interface-edge localization), and the measures IPR and IIPR (inverse participation ratio and its interface-restricted version) are used to distinguish interfaced localization from one-edge excitation.","core_discovery":"The central result is that a single undriven atom in a weakly driven, waveguide-coupled atomic array with reciprocal coupling acts as a defect that pins the steady state. For normal-incidence driving with symmetric phases, the steady-state amplitudes have the closed form of Eq. (16), and the population of the undriven atom obeys Eq. (17), with $\\lim_{\\xi\\to 0(2\\pi)} \\tilde P_m = 2/3$ for any $N\\ge 5$; the neighboring atoms each take $1/6$ in the same limit. The same saturation is robust to the incident angle, so for any $\\theta\\in(0,\\pi)$ the population collapses onto the interface as the spacing approaches an integer multiple of the wavelength. In the asymmetric-driving scheme without a defect, the paper identifies strong interfaced localization at small directionality $D$ and small $\\xi$, and derives the condition $\\cos\\theta_1-\\cos\\theta_2 = 2n\\pi/(m\\xi)$ that traces the local minima of the interface-edge ratio and thereby predicts the parameter regimes of interfaced localization. The paper closes by showing that removing drives on several non-adjacent atoms yields multiple, spatially separated single-site localizations, which it proposes as a route to quantum memory.","pith_inferences":["The $2/3$ saturation is independent of $N$ and of the incident angle, which suggests a universal ratio that may also appear in other one-dimensional driven-dissipative lattices with reciprocal infinite-range interactions; a test would be to look for the same fixed point in a chain with power-law instead of phase-only interactions.","Because the paper's model excludes spacings $\\xi < 0.1\\pi$ while the saturation is reached at $\\xi\\to 0$, the exact $2/3$ could be modified by near-field corrections; computing the first-order correction in the near-field coupling would show whether the ratio shifts or stays.","The multiple-site localization by removing drives on non-adjacent atoms is, in effect, a classical bit pattern written into the drive; one could test whether the pattern fidelity degrades with $N$ or with moderate nonreciprocity $D\\neq 0$.","Eq. (14) is presented as an empirical formula; a first-principles derivation from the phase-matching condition of the two subchains would likely turn it into a proof and may reveal an underlying Talbot-like interference structure."],"forward_implications":["With one undriven atom and spacing near an integer multiple of the wavelength, the steady-state population of that atom saturates at exactly $2/3$ (and its neighbors at $1/6$ each), independent of the number of atoms in the array.","Under asymmetric driving, interfaced two-site localization appears predominantly for reciprocal or weak directionality and small spacing; the condition $\\cos\\theta_1 - \\cos\\theta_2 = 2n\\pi/(m\\xi)$ predicts the angular combinations that sustain it and shows that localization becomes harder as the array grows.","Removing the drive from several non-adjacent atoms produces multiple separated single-site localizations, providing a deterministic way to place excitations at chosen positions along the array.","The single-site localization survives non-guided decay up to $\\gamma_{ng}\\approx 0.5\\gamma$, so the effect is not an artifact of perfect coupling into the waveguide.","Shifting the interface position (or choosing different undriven atoms) moves the localized peak, offering a control knob for routing excitations without changing the interatomic spacing."],"supporting_citations":[{"why":"Provides the atom-waveguide QED platform and the coupling-efficiency numbers used to justify the model.","marker":"[1]"},{"why":"Defines the directionality parameter $D$ used throughout the localization analysis.","marker":"[10]"},{"why":"Supplies the chiral Lindblad master equation and coupling matrix underlying Eq. (1).","marker":"[12]"},{"why":"Gives the dissimilar-array steady-state phases that the asymmetric-driving scheme is compared against.","marker":"[45]"},{"why":"Provides the steady-state phase diagram and five-phase classification that organize the results.","marker":"[60]"},{"why":"Cited as the near-field correction that restricts the validity of the model to $\\xi\\gtrsim 0.1\\pi$.","marker":"[63]"},{"why":"Also cited for near-field and polarization effects limiting the small-$\\xi$ regime.","marker":"[64]"}],"fun_headline_variants":["One dark atom pins two-thirds of the field","Single defect atom traps excitation in a waveguide","Atomic defect localizes 2/3 of steady-state population","Weak drive, strong pin: one atom holds the state","Two-thirds on one atom: a quantum memory trick"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The model treats the waveguide-mediated dipole-dipole coupling through the simple phase factors $e^{ik|x_j-x_\\ell|}$ and neglects near-field corrections; the paper itself notes this is valid only for spacings $\\xi\\gtrsim 0.1\\pi$, yet the central $2/3$ saturation is derived at $\\xi\\to 0$.","fun_headline_variants_meta":{"raw":{"variants":["One dark atom pins two-thirds of the field","Single defect atom traps excitation in a waveguide","Atomic defect localizes 2/3 of steady-state population","Weak drive, strong pin: one atom holds the state","Two-thirds on one atom: a quantum memory trick"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000297,"raw_usage":{"total_tokens":1772,"prompt_tokens":1047,"completion_tokens":725,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":663,"completion_tokens_details":{"reasoning_tokens":649}},"tokens_in":663,"tokens_out":725,"duration_ms":7537,"temperature":1.0,"reasoning_tokens":649,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T15:32:53.111254+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Directly measure or simulate the steady-state population of the undriven atom in a defect-driven reciprocal array at small spacing (e.g., $N=100$, $\\xi=0.05\\pi$) with a full near-field dipole-dipole tensor; if $\\tilde P_m$ deviates from $2/3$, the idealized phase-only coupling is the failing assumption. A less expensive version: compute the first near-field correction to Eq. (17) analytically and check whether the limiting ratio persists.","supporting_citations":[{"cited_title":"Chung, I","cited_arxiv_id":null,"evidence_quote":"Gives the dissimilar-array steady-state phases that the asymmetric-driving scheme is compared against."},{"cited_title":"Fiscelli, L","cited_arxiv_id":null,"evidence_quote":"Cited as the near-field correction that restricts the validity of the model to $\\xi\\gtrsim 0.1\\pi$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Also cited for near-field and polarization effects limiting the small-$\\xi$ regime."}],"review_version":1}