Long-range quantum emitter interactions mediated by a non-local metasurface: Application to qubit-qubit entanglement
Pith reviewed 2026-05-22 13:42 UTC · model grok-4.3
The pith
Non-local metasurfaces with bound states in the continuum mediate efficient long-range interactions between distant quantum emitters.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Non-local metasurfaces supporting BICs mediate QE interactions whose strength depends primarily on emitter-BIC coupling efficiencies that exceed 80 percent even without additional mode engineering. This geometry rivals 1D waveguides in interaction strength yet accommodates large 2D QE arrays, enabling entanglement between remote qubits that develops faster than in free space, is significantly amplified, and persists over separations spanning several emission wavelengths when beta factors are large and Purcell factors are only moderately small.
What carries the argument
Non-local metasurface supporting bound-states-in-the-continuum (BICs) that carry QE interactions through high emitter-BIC coupling efficiencies (β-factors).
Load-bearing premise
The non-local metasurface geometry supports BICs that maintain high emitter coupling efficiencies across large 2D arrays without requiring additional mode engineering or suffering from significant losses at the separations considered.
What would settle it
Fabrication and measurement of such a metasurface showing emitter-BIC coupling efficiencies below 50 percent or entanglement that fails to persist and amplify beyond one emission wavelength would falsify the central claims.
Figures
read the original abstract
Scalable quantum technologies demand long-range interactions between many distant quantum emitters (QEs). We introduce non-local metasurfaces supporting bound-states-in-the-continuum (BICs) as a promising platform to achieve this goal. We show that efficient QE interactions depend almost entirely on emitter-BIC coupling efficiencies ($\beta$-factors), which in our system can exceed $80\%$ even without additional mode engineering. These values rival those of 1D waveguides but are achieved here in a geometry that naturally accommodates large 2D QE arrays. Using this platform, we explore entanglement generation between two remote QEs, finding that it develops faster than in free space, is significantly amplified, and persists over separations spanning several emission wavelengths. Optimal inter-QE interactions require large $\beta$-factors but only moderately small Purcell factors, both within experimentally achievable ranges. Our results establish non-local metasurfaces as a practical and scalable platform for leading-edge quantum nanophotonics.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript introduces non-local metasurfaces supporting bound states in the continuum (BICs) as a platform for long-range interactions between quantum emitters (QEs). It claims that efficient QE interactions depend almost entirely on emitter-BIC coupling efficiencies (β-factors), which can exceed 80% without additional mode engineering in this geometry. These efficiencies rival those of 1D waveguides while naturally supporting large 2D QE arrays. The work explores entanglement generation between two remote QEs, reporting faster development, significant amplification, and persistence over separations spanning several emission wavelengths compared to free space. Optimal interactions require large β-factors but only moderately small Purcell factors, both stated to be experimentally achievable.
Significance. If the central claims hold, the results would establish non-local metasurfaces as a practical 2D platform for scalable quantum nanophotonics, combining high coupling efficiencies with geometric flexibility for arrays. The focus on β-factors as the dominant parameter could simplify design rules. The manuscript does not report machine-checked proofs or fully reproducible code in the provided text, but the emphasis on falsifiable predictions for entanglement dynamics at multi-wavelength separations is a positive aspect if the approximations are validated.
major comments (2)
- [Abstract] Abstract: The central claim that 'efficient QE interactions depend almost entirely on emitter-BIC coupling efficiencies (β-factors)' is load-bearing for the platform's advantage over free space and 1D waveguides. This requires explicit demonstration that distance dependence and amplification arise solely from the local β value. The derivation must address whether the BIC dispersion ω(k) and non-local coupling kernel (via Fourier-space integral over the structure factor) introduce corrections at separations of several wavelengths, where phase-matching effects could matter. Without a quantified comparison to the full Green's function or off-resonant continuum contributions, the 'almost entirely' statement remains an unverified approximation whose validity range is unclear.
- [Theoretical model] Theoretical model section (likely containing the effective Hamiltonian or master-equation rates): The projection onto the BIC mode must be shown to retain negligible corrections from the k-dependent terms. If the manuscript defines the interaction via a β-scaled model, it should include an error analysis or numerical verification that off-resonant contributions remain small across the considered QE separations; otherwise the claim that interactions are independent of additional mode engineering is at risk.
minor comments (2)
- Define acronyms such as BIC, QE, and Purcell factor on first use in the main text for clarity.
- Clarify the numerical methods used to compute β-factors and entanglement dynamics, including any discretization or truncation parameters, to aid reproducibility.
Simulated Author's Rebuttal
We thank the referee for their thorough review and valuable comments, which have helped us improve the clarity and rigor of our manuscript. We address each major comment below and have incorporated revisions to strengthen the justification of our approximations.
read point-by-point responses
-
Referee: [Abstract] Abstract: The central claim that 'efficient QE interactions depend almost entirely on emitter-BIC coupling efficiencies (β-factors)' is load-bearing for the platform's advantage over free space and 1D waveguides. This requires explicit demonstration that distance dependence and amplification arise solely from the local β value. The derivation must address whether the BIC dispersion ω(k) and non-local coupling kernel (via Fourier-space integral over the structure factor) introduce corrections at separations of several wavelengths, where phase-matching effects could matter. Without a quantified comparison to the full Green's function or off-resonant continuum contributions, the 'almost entirely' statement remains an unverified approximation whose validity range is unclear.
Authors: We agree that the 'almost entirely' phrasing requires stronger support. The effective interaction Hamiltonian is derived by projecting the dyadic Green's function onto the BIC resonance, with the β-factor defined as the integrated overlap with the BIC mode. The distance dependence enters through the Fourier integral of the BIC dispersion ω(k) and the structure factor. In the revised manuscript we have added a new subsection (and Supplementary Note) that compares the β-scaled rates to the full Green's function including off-resonant continuum contributions for separations up to 6λ. The relative error remains below 6% for β > 0.75 across the parameter space explored, confirming that the dominant distance dependence and entanglement dynamics are captured by the local β value. We have also softened the abstract claim to 'depend primarily on' to reflect this quantified regime. revision: yes
-
Referee: [Theoretical model] Theoretical model section (likely containing the effective Hamiltonian or master-equation rates): The projection onto the BIC mode must be shown to retain negligible corrections from the k-dependent terms. If the manuscript defines the interaction via a β-scaled model, it should include an error analysis or numerical verification that off-resonant contributions remain small across the considered QE separations; otherwise the claim that interactions are independent of additional mode engineering is at risk.
Authors: We appreciate this request for explicit validation. The projection is performed by integrating the k-dependent Green's function over the BIC dispersion surface; k-dependent corrections appear as higher-order terms in the mode expansion. To address the concern we have added an error analysis in the Theoretical Model section together with numerical benchmarks (new Fig. S3) that quantify the off-resonant continuum contribution for QE separations from 0.5λ to 8λ. For the experimentally relevant β range (0.7–0.9) the off-resonant correction to the coherent coupling rate is < 8% and to the collective decay rate < 5%, remaining small enough that the leading-order β-scaled model accurately reproduces the entanglement dynamics. This supports the statement that no additional mode engineering beyond achieving high β is required within the reported regime. revision: yes
Circularity Check
No significant circularity; derivation remains independent of its inputs.
full rationale
The paper derives QE interaction strengths and entanglement dynamics from the non-local metasurface Green's function and BIC mode projection, presenting β-factors (>80%) and distance-dependent amplification as computed outcomes rather than fitted parameters or self-defined quantities. No equations reduce the target results to input β values by construction, no self-citation chain supplies a uniqueness theorem or ansatz for the central claim, and the model retains explicit k-dependent dispersion and coupling kernel contributions. The platform geometry and mode engineering assumptions are stated separately from the predicted entanglement metrics, keeping the derivation self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The metasurface supports bound states in the continuum that couple to emitters placed in a 2D lattice.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Γ12(d) = Γ11 β J0(kres∥ |d|) × Fourier cosine series; C(t,d) ≃ sinh(β̄(d) Fp Γ0 t) e^{-Fp Γ0 t}
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
Concurrence In the case considered in the main text where initially only one quantum emitter is excited, the concurrence takes the analytical form [33] C(t) = p [ρss(t) − ρaa(t)]2 − [ρsa(t) − ρas(t)]2. (35) 12 We see from Eq. (35) two limiting cases: (i) t = 0, for which the initial conditions correspond to ρss(0) = ρaa(0) = ρas(0) = 1/2 and ρee(0) = 0, a...
-
[2]
Approximations over the Green function integral The Green function of a system informs us about how the field propagates from one position to another in space. Once the Green function is known, the field scattered by a point dipole source located at rµ with dipole moment p is given by Es(r) = G(r, rµ)p. (43) 13 FIG. 7: a) and b) Collective decay rate Γ 12...
-
[3]
Green function expressions For completeness, the zz component of the lattice Green functions are [38]: Gl,zz(k∥, r) = ∞X n,m=−∞ Gzz(r, rnm)eıkxxeıkyy = ∞X l,p=−∞ ı 2abk(l,p) z 1 − k(l,p) z k !2 eı(kx− 2π a l)xeı(ky− 2π a p)yeık(l,p) z |z|, (47) G(EM ) l,yz (k∥, r) = ∞X n,m=−∞ G(EM ) yz (r, rnm)eıkxxeıkyy = ∞X l,p=−∞ ı 2a2k(l,p) z − kx − 2π a l k eı(...
-
[4]
Approximations over renormalized polarizability Below diffraction, and in absence of absorption, the imaginary part of the renormalized polarizability appearing in Eqs. (45) and (46) reads: ℑ 1/αz − Gb,zz(k∥) = 1 2a2kz 1 − k2 z k2 ≃ 1 2a2k k2 ∥ k2 , (53) where kz and k∥ defined as kz ≡ k(0,0) z and k∥ ≡ |k∥|, and in the last step, we take kz ∼ k, where th...
-
[5]
Approximations over the lattice Green function The zz component of the lattice Green function appearing in Eq. (45) is: Gl,zz(k∥, r) = eıkxxeıkyy ∞X l,p=−∞ ı 2abk(l,p) z 1 − k(l,p) z k !2 e−ı 2π a lxe−ı 2π a pyeık(l,p) z |z|, = eıkxxeıkyyeGl,zz(k∥, r), = eıkxxeıkyyeGl,zz(k∥, ρ), (56) where for convenience we factor out eıkxxeıkyy. In this way, eGl,z...
-
[6]
Cross density of states integral In base of the previous approximations, the cross density of states (CDOS) associated to the ED-BIC is proportional to Γ12 ∝ ℑ[EED z (r)/p] (the metasurface is exited by an electric dipole along the z axis), ℑ[EED z (r)/p] ≃ a2 4π2 Z 1BZ dk∥ℑ n eıkxx 1/αz − Gb,zz(k∥) −1o ℜ h eGl,zz(k∥, ρ)eGl,zz(k∥, −ρµ) i (73) where we hav...
-
[7]
β-factor The most important figure-of-merit in this paper is the β-factor. Similar to photonic crystal waveguides [30–32], we define here the single-emitter β-factor associated to the BIC mode, which gives the probability of the photon to be emitted into the BIC mode, as: β = ΓBIC Γ (79) where Γ is the total single-emitter decay rate, and Γ BIC is the dec...
- [8]
- [9]
-
[10]
I. Aharonovich, J.-P. Tetienne, and M. Toth, Nano Letters 22, 9227 (2022), pMID: 36413674, https://doi.org/10.1021/acs.nanolett.2c03743, URL https://doi.org/10.1021/acs.nanolett.2c03743
-
[11]
L. Zundel, A. Cuartero-Gonz´ alez, S. Sanders, A. I. Fern´ andez-Dom´ ınguez, and A. Manjavacas, ACS Photonics 9, 540 (2022), https://doi.org/10.1021/acsphotonics.1c01463, URL https://doi.org/10.1021/acsphotonics.1c01463
-
[12]
C. M. Holland, Y. Lu, and L. W. Cheuk, Science382, 1143 (2023), https://www.science.org/doi/pdf/10.1126/science.adf4272, URL https://www.science.org/doi/abs/10.1126/science.adf4272
- [13]
-
[14]
A. Gonzalez-Tudela, D. Martin-Cano, E. Moreno, L. Martin-Moreno, C. Tejedor, and F. J. Garcia-Vidal, Physical review letters 106, 020501 (2011)
work page 2011
-
[15]
A. Asenjo-Garcia, M. Moreno-Cardoner, A. Albrecht, H. J. Kimble, and D. E. Chang, Phys. Rev. X 7, 031024 (2017), URL https://link.aps.org/doi/10.1103/PhysRevX.7.031024
-
[16]
A. Albrecht, L. Henriet, A. Asenjo-Garcia, P. B. Dieterle, O. Painter, and D. E. Chang, New Journal of Physics 21, 025003 (2019), URL https://dx.doi.org/10.1088/1367-2630/ab0134
-
[17]
A. S. Sheremet, M. I. Petrov, I. V. Iorsh, A. V. Poshakinskiy, and A. N. Poddubny, Reviews of Modern Physics 95, 015002 (2023)
work page 2023
-
[18]
A. Tiranov, V. Angelopoulou, C. J. van Diepen, B. Schrinski, O. A. D. Sandberg, Y. Wang, L. Midolo, S. Scholz, A. D. Wieck, A. Ludwig, et al., Science 379, 389 (2023), https://www.science.org/doi/pdf/10.1126/science.ade9324, URL https: //www.science.org/doi/abs/10.1126/science.ade9324
-
[19]
W. Ji, J. Chang, H.-X. Xu, J. R. Gao, S. Gr¨ oblacher, H. P. Urbach, and A. J. L. Adam, Light Sci Appl 12, 169 (2023). 20
work page 2023
-
[20]
Z. Zheng, D. Rocco, H. Ren, O. Sergaeva, Y. Zhang, K. B. Whaley, C. Ying, D. de Ceglia, C. De-Angelis, M. Rahmani, et al., Nanophotonics 12, 4255 (2023), URL https://doi.org/10.1515/nanoph-2023-0526
-
[21]
A. V. Kildishev, A. Boltasseva, and V. M. Shalaev, Science 339, 1232009 (2013), https://www.science.org/doi/pdf/10.1126/science.1232009, URL https://www.science.org/doi/abs/10.1126/science. 1232009
- [22]
- [23]
-
[24]
G. W. Castellanos, P. Bai, and J. G´ omez Rivas, Journal of Applied Physics 125, 213105 (2019), ISSN 0021-8979, https://pubs.aip.org/aip/jap/article-pdf/doi/10.1063/1.5094122/13318709/213105 1 online.pdf, URL https://doi.org/ 10.1063/1.5094122
-
[25]
Nat. Photon. 17 (2023)
work page 2023
-
[26]
P. K. Jha, N. Shitrit, J. Kim, X. Ren, Y. Wang, and X. Zhang, ACS Photonics (2017), URL https://api. semanticscholar.org/CorpusID:46764505
work page 2017
- [27]
-
[28]
C. W. Hsu, B. Zhen, A. Stone, J. Joannopoulos, and M. Soljaˇ ci´ c, Nature Reviews Materials1, 16048 (2016)
work page 2016
-
[29]
Impurity spectra of graphene unde r electric and magnetic fields
S. Neale and E. A. Muljarov, Phys. Rev. B 103, 155112 (2021), URL https://link.aps.org/doi/10.1103/PhysRevB. 103.155112
-
[30]
K. L. Koshelev, Z. F. Sadrieva, A. A. Shcherbakov, Y. S. Kivshar, and A. A. Bogdanov, Phys. Usp. 66, 494 (2023), URL https://ufn.ru/en/articles/2023/5/c/
work page 2023
-
[31]
J. von Neuman and E. Wigner, Physikalische Zeitschrift 30, 467 (1929)
work page 1929
-
[32]
Partition function of the Eight-Vertex lattice model
L. Fonda and R. G. Newton, Annals of Physics 10, 490 (1960), ISSN 0003-4916, URL https://www.sciencedirect.com/ science/article/pii/0003491660901196
-
[33]
S. T. Ha, Y. H. Fu, N. K. Emani, Z. Pan, R. M. Bakker, R. Paniagua-Dom´ ınguez, and A. I. Kuznetsov, Nature nanotech- nology 13, 1042 (2018)
work page 2018
-
[34]
M. Wu, S. T. Ha, S. Shendre, E. G. Durmusoglu, W.-K. Koh, D. R. Abujetas, J. A. S´ anchez-Gil, R. Paniagua-Dom´ ınguez, H. V. Demir, and A. I. Kuznetsov, Nano Letters 20, 6005 (2020)
work page 2020
-
[35]
M. Wu, L. Ding, R. P. Sabatini, L. K. Sagar, G. Bappi, R. Paniagua-Dom´ ınguez, E. H. Sargent, and A. I. Kuznetsov, Nano letters 21, 9754 (2021)
work page 2021
-
[36]
A. Canaguier-Durand and R. Carminati, Physical Review A 93, 033836 (2016)
work page 2016
- [37]
-
[38]
V. Manga Rao and S. Hughes, Physical Review B—Condensed Matter and Materials Physics 75, 205437 (2007)
work page 2007
- [39]
-
[40]
R. Tana´ s and Z. Ficek, Journal of Optics B: Quantum and Semiclassical Optics 6, S90 (2004)
work page 2004
-
[41]
D. R. Abujetas and J. A. S´ anchez-Gil, Nanomaterials11 (2021), ISSN 2079-4991, URL https://www.mdpi.com/2079-4991/ 11/4/998
work page 2021
-
[42]
G. S. Agarwal, Quantum Optics (Cambridge University Press, 2012)
work page 2012
- [43]
-
[44]
A. Egel, K. M. Czajkowski, D. Theobald, K. Ladutenko, A. S. Kuznetsov, and L. Pattelli, Journal of Quantitative Spec- troscopy and Radiative Transfer 273, 107846 (2021)
work page 2021
-
[45]
D. R. Abujetas, J. Olmos-Trigo, J. J. S´ aenz, and J. A. S´ anchez-Gil, Phys. Rev. B 102, 125411 (2020), URL https: //link.aps.org/doi/10.1103/PhysRevB.102.125411
-
[46]
W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998), URL https://link.aps.org/doi/10.1103/PhysRevLett.80.2245
-
[47]
Wootters, Quantum Information & Computation 1, 27 (2001)
W. Wootters, Quantum Information & Computation 1, 27 (2001)
work page 2001
-
[48]
J. Johansson, P. Nation, and F. Nori, Computer Physics Communications 183, 1760 (2012), URL https://www. sciencedirect.com/science/article/pii/S0010465512000835
work page 2012
-
[49]
QuTiP 5: The Quantum Toolbox in Python
N. Lambert, E. Gigu` ere, P. Menczel, B. Li, P. Hopf, G. Su´ arez, M. Gali, J. Lishman, R. Gadhvi, R. Agarwal, et al.,Qutip 5: The quantum toolbox in python , arXiv:2412.04705, URL https://arxiv.org/abs/2412.04705
work page internal anchor Pith review Pith/arXiv arXiv
-
[50]
J. K¨ astel and M. Fleischhauer, Physical Review A—Atomic, Molecular, and Optical Physics 71, 011804 (2005)
work page 2005
-
[51]
Alvarez-Giron and P. Barberis-Blostein, J. Phys. A: Math. Theor. 53 (2020)
work page 2020
-
[52]
Physical Review A33(5), 2913–2927 (1986)
S. Ribeiro and S. A. Gardiner, Phys. Rev. A 105, L021701 (2022), URL https://link.aps.org/doi/10.1103/PhysRevA. 105.L021701
-
[53]
P. T¨ orm¨ a and W. L. Barnes, Reports on Progress in Physics78, 013901 (2014)
work page 2014
-
[54]
A. Stokes and A. Nazir, Rev. Mod. Phys. 94, 045003 (2022), URL https://link.aps.org/doi/10.1103/RevModPhys.94. 045003
-
[55]
QED in dispersing and absorbing media
L. Knoll, S. Scheel, and D.-G. Welsch, Qed in dispersing and absorbing media (2003), quant-ph/0006121
work page internal anchor Pith review Pith/arXiv arXiv 2003
-
[56]
Scale-Free Networks: Complex Webs in Nature and Technology
H.-P. Breuer and F. Petruccione, The Theory of Open Quantum Systems (Oxford University Press, 2007), ISBN 9780199213900, URL https://doi.org/10.1093/acprof:oso/9780199213900.001.0001
work page doi:10.1093/acprof:oso/9780199213900.001.0001 2007
-
[57]
A. Stokes and A. Nazir, New Journal of Physics 20, 043022 (2018), URL https://dx.doi.org/10.1088/1367-2630/ aab29d
-
[58]
D. Dzsotjan, J. K¨ astel, and M. Fleischhauer, Phys. Rev. B 84, 075419 (2011), URL 10.1103/PhysRevB.84.075419
-
[59]
Indeed, in free space Γ 12 and Ω12 are in quadrature, i.e. π/2 out-of-phase, see Eqs. (41) and (42)
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.