pith. sign in

arxiv: 2505.15751 · v2 · pith:PXPU2ZPZnew · submitted 2025-05-21 · 🪐 quant-ph · cond-mat.mes-hall

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

classification 🪐 quant-ph cond-mat.mes-hall
keywords quantum emittersbound states in the continuummetasurfaceslong-range interactionsentanglementbeta factorsnanophotonicsqubit arrays
0
0 comments X

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.

The paper establishes that non-local metasurfaces supporting bound states in the continuum enable strong interactions between quantum emitters placed far apart. These interactions depend almost entirely on the efficiency with which emitters couple to the bound states, which reaches over 80 percent in the proposed geometry without extra engineering. This performance matches one-dimensional waveguides while fitting large two-dimensional arrays naturally. The setup produces entanglement between remote emitters that forms faster than in free space, grows stronger, and holds across distances of several wavelengths. Optimal results occur with large coupling efficiencies but only moderately small Purcell factors, both within experimental reach.

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

Figures reproduced from arXiv: 2505.15751 by Adam Stokes, Ahsan Nazir, Diego Romero Abujetas, Emmanuel Lassalle, Hannah Riley, Ramon Paniagua-Dominguez.

Figure 1
Figure 1. Figure 1: FIG. 1: Long-range metasurface-mediated interactions between two qubits. (a) Initially uncorrelated QEs (qubits) sponta [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: LDOS Γ [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Lateral spatial dependence of the CDOS Γ [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: The CDOS Γ [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: (a) and (b): Decay rate Γ [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: Concurrence [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7: a) and b) Collective decay rate Γ [PITH_FULL_IMAGE:figures/full_fig_p013_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8: Absolute error between numerics and analytics for the collective decay rate, in the case of (a) the ED-BIC and (b) the [PITH_FULL_IMAGE:figures/full_fig_p018_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9: Concurrence in free space as a function of the normalized time Γ [PITH_FULL_IMAGE:figures/full_fig_p019_9.png] view at source ↗
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.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

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)
  1. [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.
  2. [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)
  1. Define acronyms such as BIC, QE, and Purcell factor on first use in the main text for clarity.
  2. Clarify the numerical methods used to compute β-factors and entanglement dynamics, including any discretization or truncation parameters, to aid reproducibility.

Simulated Author's Rebuttal

2 responses · 0 unresolved

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
  1. 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

  2. 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

0 steps flagged

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

0 free parameters · 1 axioms · 0 invented entities

The central claims rest on standard electromagnetic modeling of metasurfaces and quantum optics master equations for emitter dynamics. No new particles or forces are introduced. The key unverified elements are the numerical values of β-factors and Purcell factors obtained from the specific geometry.

axioms (1)
  • domain assumption The metasurface supports bound states in the continuum that couple to emitters placed in a 2D lattice.
    Invoked when stating that β-factors exceed 80% without additional mode engineering.

pith-pipeline@v0.9.0 · 5717 in / 1341 out tokens · 32623 ms · 2026-05-22T13:42:47.522795+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

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

59 extracted references · 59 canonical work pages · 2 internal anchors

  1. [1]

    (35) 12 We see from Eq

    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. [2]

    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

    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. [3]

    (49) The zz component of the renormalized polarizability reads: Gb,zz(k∥) = lim r→0 Gl,zz(k∥, r) − Gzz(r, 0)

    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. [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. [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. [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. [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. [8]

    Lodahl, S

    P. Lodahl, S. Mahmoodian, and S. Stobbe, Rev. Mod. Phys. 87, 347 (2015), URL https://link.aps.org/doi/10.1103/ RevModPhys.87.347

  9. [9]

    J. B. Trebbia, Q. Deplano, P. Tamarat, and B. Lounis, Nature Commun. 13, 2962 (2022), 2109.10584

  10. [10]

    Aharonovich, J.-P

    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. [11]

    Zundel, A

    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. [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. [13]

    F. Shah, T. L. Patti, O. Rubies-Bigorda, and S. F. Yelin, Phys. Rev. A 109, 012613 (2024), 2306.08555

  14. [14]

    Gonzalez-Tudela, D

    A. Gonzalez-Tudela, D. Martin-Cano, E. Moreno, L. Martin-Moreno, C. Tejedor, and F. J. Garcia-Vidal, Physical review letters 106, 020501 (2011)

  15. [15]

    Asenjo-Garcia, M

    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. [16]

    Albrecht, L

    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. [17]

    A. S. Sheremet, M. I. Petrov, I. V. Iorsh, A. V. Poshakinskiy, and A. N. Poddubny, Reviews of Modern Physics 95, 015002 (2023)

  18. [18]

    Tiranov , author V

    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. [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

  20. [20]

    Zheng, D

    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. [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. [22]

    Yu and F

    N. Yu and F. Capasso, Nature materials 13 2, 139 (2014), URL https://api.semanticscholar.org/CorpusID:25771706

  23. [23]

    Soukoulis and M

    C. Soukoulis and M. Wegener, Nature Photonics 5, 523 (2011)

  24. [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. [25]

    Nat. Photon. 17 (2023)

  26. [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

  27. [27]

    Zundel, A

    L. Zundel, A. Cuartero-Gonz´ alez, S. Sanders, A. I. Fern´ andez-Dom´ ınguez, and A. Manjavacas, ACS Photonics 9, 540 (2022)

  28. [28]

    C. W. Hsu, B. Zhen, A. Stone, J. Joannopoulos, and M. Soljaˇ ci´ c, Nature Reviews Materials1, 16048 (2016)

  29. [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. [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/

  31. [31]

    von Neuman and E

    J. von Neuman and E. Wigner, Physikalische Zeitschrift 30, 467 (1929)

  32. [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. [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)

  34. [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)

  35. [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)

  36. [36]

    Canaguier-Durand and R

    A. Canaguier-Durand and R. Carminati, Physical Review A 93, 033836 (2016)

  37. [37]

    Lecamp, P

    G. Lecamp, P. Lalanne, and J. Hugonin, Physical review letters 99, 023902 (2007)

  38. [38]

    Manga Rao and S

    V. Manga Rao and S. Hughes, Physical Review B—Condensed Matter and Materials Physics 75, 205437 (2007)

  39. [39]

    Arcari, I

    M. Arcari, I. S¨ ollner, A. Javadi, S. Lindskov Hansen, S. Mahmoodian, J. Liu, H. Thyrrestrup, E. H. Lee, J. D. Song, S. Stobbe, et al., Physical review letters 113, 093603 (2014)

  40. [40]

    Tana´ s and Z

    R. Tana´ s and Z. Ficek, Journal of Optics B: Quantum and Semiclassical Optics 6, S90 (2004)

  41. [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

  42. [42]

    G. S. Agarwal, Quantum Optics (Cambridge University Press, 2012)

  43. [43]

    Stokes and A

    A. Stokes and A. Nazir, New Journal of Physics 20, 043022 (2018)

  44. [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)

  45. [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. [46]

    W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998), URL https://link.aps.org/doi/10.1103/PhysRevLett.80.2245

  47. [47]

    Wootters, Quantum Information & Computation 1, 27 (2001)

    W. Wootters, Quantum Information & Computation 1, 27 (2001)

  48. [48]

    Johansson, P

    J. Johansson, P. Nation, and F. Nori, Computer Physics Communications 183, 1760 (2012), URL https://www. sciencedirect.com/science/article/pii/S0010465512000835

  49. [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

  50. [50]

    K¨ astel and M

    J. K¨ astel and M. Fleischhauer, Physical Review A—Atomic, Molecular, and Optical Physics 71, 011804 (2005)

  51. [51]

    Barberis-Blostein, J

    Alvarez-Giron and P. Barberis-Blostein, J. Phys. A: Math. Theor. 53 (2020)

  52. [52]

    Available: https://link.aps.org/doi/10.1103/PhysRevA

    S. Ribeiro and S. A. Gardiner, Phys. Rev. A 105, L021701 (2022), URL https://link.aps.org/doi/10.1103/PhysRevA. 105.L021701

  53. [53]

    T¨ orm¨ a and W

    P. T¨ orm¨ a and W. L. Barnes, Reports on Progress in Physics78, 013901 (2014)

  54. [54]

    Stokes and A

    A. Stokes and A. Nazir, Rev. Mod. Phys. 94, 045003 (2022), URL https://link.aps.org/doi/10.1103/RevModPhys.94. 045003

  55. [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

  56. [56]

    Liesen, Z

    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

  57. [57]

    Cramer and J

    A. Stokes and A. Nazir, New Journal of Physics 20, 043022 (2018), URL https://dx.doi.org/10.1088/1367-2630/ aab29d

  58. [58]

    Dzsotjan, J

    D. Dzsotjan, J. K¨ astel, and M. Fleischhauer, Phys. Rev. B 84, 075419 (2011), URL 10.1103/PhysRevB.84.075419

  59. [59]

    π/2 out-of-phase, see Eqs

    Indeed, in free space Γ 12 and Ω12 are in quadrature, i.e. π/2 out-of-phase, see Eqs. (41) and (42)