How Physicists Beamed Entangled Light Through Open City Air Without Any Fiber Cables

On August 19, at 12:26 a.m. Eastern Time, an optical transceiver atop Stony Brook University’s Health Sciences Center fired an unguided, continuous stream of entangled photons across 21.0 kilometers (13.05 miles) of open atmosphere, striking a 0.6-meter telescope at Brookhaven National Laboratory’s Quantum Lighthouse in Upton, New York. The received quantum states demonstrated a Bell state fidelity of 94.2% $\pm$ 1.1% and a Clauser-Horne-Shimony-Holt (CHSH) parameter $S = 2.684 \pm 0.038$—violating the classical realism limit of $S \le 2$ by 18 standard deviations without utilizing a single millimeter of underground optical fiber.

The test established the first permanent municipal-scale free-space quantum optical network leg in the United States. It recorded an average Quantum Bit Error Rate (QBER) of 2.38% under daytime background radiation exceeding $10^5$ lux and maintained a coincident photon pair detection rate of 4.22 kilobits per second across a channel exhibiting 28.5 decibels (dB) of atmospheric attenuation.

The demonstration provides empirical proof that unguided open air can serve as a transmission backbone for quantum information across metropolitan terrains, bypassing the physical, spectral, and financial constraints of dark fiber networks.

+----------------------------------------------------------------------------------------------------+
|               LONG ISLAND FREE-SPACE QUANTUM OPTICAL LINK (21.0 KM TRANSECT)                       |
+----------------------------------------------------------------------------------------------------+
|  TRANSMITTER NODE (Stony Brook "Watchtower")         RECEIVER NODE (Brookhaven "Lighthouse")       |
|  - Altitude: 110 m above sea level                   - Altitude: 85 m above sea level              |
|  - SPDC Source: PPKTP Crystal @ 810 nm               - Aperture: 0.60 m (23.6 in) Parabolic Mirror |
|  - Generation Rate: 1.2 x 10^7 pairs/s               - Adaptive Optics: 241-Actuator DM @ 2.0 kHz  |
|  - Transmitter Aperture: 250 mm Beam Expander        - Spatial Filter: 5.0 µm Single-Mode Core     |
|  - Temporal Synchronization: GPS/PTP < 50 ps         - Detector: SNSPD Array (Eff. 95%, DCR < 5 Hz)|
+--------------------------------------------------+-------------------------------------------------+
                                                   |
                             21.0 km Open Atmosphere Line-of-Sight
                             Beam Divergence: 1.65 µrad (Diffraction Limit)
                             Atmospheric Turbulence: Cn2 = 1.8 x 10^-14 m^-2/3
                             Total Channel Attenuation: 28.5 dB

Technical Metrics: Stony Brook–Brookhaven Free-Space Optical Quantum Link

Parameter	Measured Laboratory Value	Open-Air Field Metric (21 km)	Performance Tolerance / Threshold
Transmission Distance	0.00 km (Direct Bench)	21.0 km (13.05 miles)	Multi-node link criteria
Entangled State Fidelity ($F$)	$98.6\% \pm 0.3\%$	$94.2\% \pm 1.1\%$	$> 85.0\%$ for secure QKD
CHSH Bell Parameter ($S$)	$2.812 \pm 0.015$	$2.684 \pm 0.038$	$S > 2.000$ (Quantum Bound $\le 2\sqrt{2}$)
Violation Significance	$54.1\sigma$	$18.0\sigma$	$> 5\sigma$ (Standard discovery threshold)
Quantum Bit Error Rate (QBER)	$0.62\%$	$2.38\%$ (Avg Day/Night)	$< 11.0\%$ (Shor-Preskill security bound)
Secure Key Transmission Rate	$125.0\text{ kbps}$	$4.22\text{ kbps}$	Scalable commercial threshold
Total Channel Insertion Loss	$1.2\text{ dB}$	$28.5\text{ dB}$	Link budget limit: $45.0\text{ dB}$
Adaptive Optics Update Rate	N/A	$2,000\text{ Hz}$	$> 1,000\text{ Hz}$ (Turbulence Greenwood freq.)
Optical Temporal Coincidence Window	$250\text{ ps}$	$350\text{ ps}$	$< 500\text{ ps}$ to reject daylight solar flux

The Physics of Unguided Photonic Transport

Fiber-optic infrastructure has formed the backbone of early quantum networking prototypes, yet it enforces fundamental physical limits on metropolitan scalability. Standard Corning SMF-28 silica fiber exhibits minimum optical attenuation of $0.18\text{ dB/km}$ exclusively within the telecommunications C-band centered at $1550\text{ nm}$. Over a 100-kilometer fiber span, channel loss reaches $18\text{ dB}$, transmitting approximately 1.58% of input photons. Over 500 kilometers, attenuation climbs to $90\text{ dB}$, reducing transmission probability to $10^{-9}$ and rendering direct point-to-point quantum key distribution inoperable without trusted repeater nodes.

Photon Transmission Probability P(L) = 10^(-alpha * L / 10)
Where:
- alpha_fiber = 0.18 dB/km (at 1550 nm)
- alpha_atmosphere = 0.05 to 0.45 dB/km (clear air across 780-1550 nm windows)

Furthermore, atomic-vapor quantum memories, trapped-ion processors, and neutral-atom registers operate natively at shorter wavelengths: $780\text{ nm}$ and $795\text{ nm}$ for Rubidium ($^{87}\text{Rb}$), $852\text{ nm}$ for Cesium ($^{133}\text{Cs}$), and $369.5\text{ nm}$ for Ytterbium ions ($^{171}\text{Yb}^+$). In standard silica fiber, attenuation at $780\text{ nm}$ exceeds $3.50\text{ dB/km}$. A 20-kilometer run of fiber at $780\text{ nm}$ incurs $70\text{ dB}$ of loss, effectively extinguishing single-photon states.

FIBER VS. FREE-SPACE LOSS PROFILE AT ATOMIC RESONANCE (780 nm):
-----------------------------------------------------------------------------
Distance (km)     Silica Fiber Loss (3.5 dB/km)    Free-Space Geometric+Extinction
-----------------------------------------------------------------------------
 1 km              3.5 dB                           12.1 dB (Coupling Dominated)
 5 km             17.5 dB                           14.8 dB
10 km             35.0 dB                           18.2 dB
20 km             70.0 dB                           24.5 dB
50 km            175.0 dB                           36.0 dB
-----------------------------------------------------------------------------

Distributing quantum entanglement through air provides an alternative medium where low-altitude clear-atmosphere transmission displays attenuation windows below $0.45\text{ dB/km}$ across the near-infrared spectrum ($700\text{ nm}$ to $1600\text{ nm}$). This enables direct, unguided interfaces between atomic quantum processors and optical ground stations without requiring quantum frequency conversion stages, which typically introduce insertion losses between $40\%$ and $70\%$ ($2.2\text{ dB}$ to $5.2\text{ dB}$).

                  ATMOSPHERIC SPECTRAL ATTENUATION PROFILE (CLEAR SKY)
  Optical Loss (dB/km)
   5.0 |                                      
   4.0 |      H2O/O2 Absorption Bands         
   3.0 |        |                             
   2.0 |       |||                  H2O Band  
   1.0 |      |||||                   |||     
   0.5 |   |--|   |--|   |---|   |----| |---- (Low-loss transmission windows)
   0.0 +---+---------+---+---+---+-----------+--->
          780nm     810nm   852nm   1064nm  1550nm  Wavelength (nm)
         [Rb D2]   [SPDC]  [Cs D2]  [Nd:YAG][Telecom C]

Architectural Layout of the Long Island Experiment

The Stony Brook–Brookhaven system was constructed to establish an all-weather, high-bandwidth entanglement link across suburban Long Island topography. The transmitter node, designated the Quantum Watchtower, is located atop the 110-meter-high Health Sciences Center at Stony Brook University. The receiver node, designated the Quantum Lighthouse, sits 85 meters above ground on Building 510 at Brookhaven National Laboratory.

+------------------------------------------------------------------------------------+
|          DETAILED TRANSMITTER OPTICAL SCHEMATIC (STONY BROOK WATCHTOWER)           |
+------------------------------------------------------------------------------------+
|                                                                                    |
|  [405 nm CW Pump Laser] ---> [Half-Wave Plate] ---> [Focusing Lens]                |
|                                                            |                       |
|                                                            v                       |
|                                          +----------------------------------+      |
|                                          | Type-II PPKTP Nonlinear Crystal  |      |
|                                          | (Temperature: 28.40 °C +/- 0.01) |      |
|                                          +----------------------------------+      |
|                                                            |                       |
|                             +------------------------------+--------------------+  |
|                             | (Signal: 810 nm, H-Pol)      | (Idler: 810 nm, V) |  |
|                             v                              v                    |  |
|                   [Local Analysis Station]     [Active Fast-Steering Mirror]    |  |
|                   - Polarization Tomography                |                    |  |
|                   - SNSPD Detector Channel A               v                    |  |
|                                                [250 mm Primary Beam Expander]   |  |
|                                                            |                    |  |
|                                                            v                    |  |
|                                               === 21.0 km Atmospheric Link ===>  |
+------------------------------------------------------------------------------------+

1. Entangled Photon Source Generation

The source generates polarization-entangled photon pairs using spontaneous parametric down-conversion (SPDC) inside a 20-millimeter-long periodically poled potassium titanyl phosphate (PPKTP) crystal configured in a Sagnac interferometer loop.

Pump Wavelength: $\lambda_p = 405.0\text{ nm}$ (diode laser stabilized to within $\pm 0.1\text{ pm}$).
Emitted Pair Wavelength: Degenerate emission at $\lambda_s = \lambda_i = 810.0\text{ nm}$.
Phase-Matched Quantum State: Maximally entangled Bell singlet state:

$$\|\Psi^-\rangle = \frac{1}{\sqrt{2}} \left( \|H\rangle_s \|V\rangle_i - \|V\rangle_s \|H\rangle_i \right)$$

Brightness: $1.2 \times 10^7\text{ pairs/second per mW of pump power}$.
Source State Fidelity: $98.6\% \pm 0.3\%$ measured at the output collimator.

The signal photon is captured locally into a 5-micron core polarization-maintaining single-mode fiber (PM-SMF) and routed directly to a local single-photon detection module. The idler photon is coupled into a customized beam-expanding telescope with an output aperture diameter $D_{\text{tx}} = 250\text{ mm}$ to form the free-space collimated beam.

+------------------------------------------------------------------------------------+
|             DETAILED RECEIVER OPTICAL SCHEMATIC (BNL QUANTUM LIGHTHOUSE)           |
+------------------------------------------------------------------------------------+
|                                                                                    |
|   == Incoming 810 nm Wavefront (Expanded to ~1.4 m spot via turbulence) ==>        |
|                                      |                                             |
|                                      v                                             |
|                     [0.60 m Parabolic Primary Mirror]                              |
|                                      |                                             |
|                                      v                                             |
|                        [Dichroic Beam Splitter]                                    |
|                         /                    \                                     |
|   (850 nm Tracking Beacon)                  (810 nm Single Photons)                |
|             |                                         |                            |
|             v                                         v                            |
|    [Fast CCD Camera & FSM]                [241-Actuator Deformable Mirror]         |
|    - Bandwidth: 5 kHz                     - Bandwidth: 2.0 kHz                     |
|    - Spatial Tracking: <0.5 µrad                      |                            |
|                                                       v                            |
|                                           [Fabry-Perot Etalon (FWHM=30 pm)]        |
|                                                       |                            |
|                                                       v                            |
|                                           [Aspheric Coupling Lens]                 |
|                                                       |                            |
|                                                       v                            |
|                                           [5.0 µm Single-Mode Fiber]               |
|                                                       |                            |
|                                                       v                            |
|                                           [Polarization Analysis / SNSPD]          |
+------------------------------------------------------------------------------------+

2. Receiver and Coupling Topography

At the Brookhaven receiving node, the arriving photons are captured by a 0.60-meter astronomical-grade telescope:

Collecting Aperture: $D_{\text{rx}} = 600\text{ mm}$ ($23.6\text{ inches}$).
Effective Focal Length: $f_{\text{eff}} = 4,800\text{ mm}$ ($f/8\text{ optical configuration}$).
Deformable Mirror Architecture: 241 piezoelectric actuators driven by a closed-loop controller operating at $2.0\text{ kHz}$.
Spatial Coupling Target: Aspheric lens launching into a 5.0-micron core single-mode fiber coupled to a cryogenically cooled Superconducting Nanowire Single-Photon Detector (SNSPD) array operating at $0.85\text{ Kelvin}$ with a quantum detection efficiency $\eta_{\text{det}} = 95\%$ and dark count rate $< 5\text{ counts/second}$.

The Channel Loss Budget: Decibel Breakdown Over 21 Kilometers

Transmitting single-photon quantum states over long-distance terrestrial links requires accounting for geometric beam expansion, atmospheric absorption, particulate scattering, wavefront phase distortion, and receiver collection mechanics.

Total Channel Loss (dB) = L_geo + L_atm + L_turb + L_opt + L_coup

                     LINK ATTENUATION BUDGET DISTRIBUTION (28.5 dB TOTAL)
 +------------------------------------------------------------------------------------+
 | Component                                            Loss Value (dB)  % of Total   |
 +------------------------------------------------------------------------------------+
 | Geometric Divergence (L_geo)                         13.4 dB          [47.0%]      |
 | Fiber Re-Coupling Mismatch (L_coup)                   6.8 dB          [23.9%]      |
 | Adaptive Optics Residual Aberrations (L_turb)         3.8 dB          [13.3%]      |
 | Atmospheric Extinction (L_atm: Rayleigh/Mie)          2.7 dB          [9.5%]       |
 | Receiver/Transmitter Optics Transmittance (L_opt)     1.8 dB          [6.3%]       |
 +------------------------------------------------------------------------------------+

===================================================================================
DECIBEL LOSS BUDGET: STONY BROOK TO BROOKHAVEN NATIONAL LAB (21.0 KM FREE-SPACE)
===================================================================================
Loss Parameter                      Mechanism / Physical Cause         Loss (dB)
-----------------------------------------------------------------------------------
1. Geometric Expansion (L_geo)      Free-space diffraction over 21 km    13.4 dB
2. Atmospheric Extinction (L_atm)    Rayleigh (0.9 dB) + Mie (1.8 dB)      2.7 dB
3. Scintillation/Turbulence (L_turb) Residual uncorrected wavefront phase  3.8 dB
4. Optical Component Losses (L_opt) 14 optical surfaces + bandpass filter  1.8 dB
5. Fiber Re-Coupling Losses (L_coup) SMF mode-field mismatch and jitter    6.8 dB
-----------------------------------------------------------------------------------
TOTAL ACCUMULATED CHANNEL LOSS                                            28.5 dB
===================================================================================

Mathematical Breakdown of Loss Components:

1. Geometric Loss ($L_{\text{geo}}$):

The theoretical diffraction-limited divergence angle $\theta_{\text{diff}}$ of the emitted beam is dictated by the transmitter aperture $D_{\text{tx}} = 0.25\text{ m}$:

$$\theta_{\text{diff}} = 1.22 \frac{\lambda}{D_{\text{tx}}} = 1.22 \frac{810 \times 10^{-9}\text{ m}}{0.25\text{ m}} = 3.95 \times 10^{-6}\text{ rad} = 3.95\text{ }\mu\text{rad}$$

Over $L = 21,000\text{ meters}$, the diffraction spot diameter $D_{\text{spot}}$ arriving at the receiver plane is:

$$D_{\text{spot}} = D_{\text{tx}} + L \cdot \theta_{\text{diff}} = 0.25\text{ m} + (21,000\text{ m} \times 3.95 \times 10^{-6}\text{ rad}) \approx 0.333\text{ m} = 33.3\text{ cm}$$

Because the receiver primary mirror diameter is $D_{\text{rx}} = 0.60\text{ m}$, the geometric collection efficiency in a vacuum would approach $100\%$. However, uncorrected atmospheric turbulence induces beam broadening and centroid displacement (beam wander), expanding the time-averaged beam footprint $\langle W_{\text{turb}} \rangle$ to approximately $1.42\text{ meters}$ in diameter:

$$L_{\text{geo}} = -10 \log_{10} \left( \frac{D_{\text{rx}}^2}{\langle W_{\text{turb}} \rangle^2} \right) = -10 \log_{10} \left( \frac{0.60^2}{1.42^2} \right) = -10 \log_{10}(0.1785) \approx 7.48\text{ dB}$$

Adding clipping margins and secondary mirror central obscuration ($15\%$ area loss, $0.71\text{ dB}$) brings the physical collection penalty to $13.4\text{ dB}$.

2. Atmospheric Extinction ($L_{\text{atm}}$):

Computed via the Beer-Lambert law:

$$T_{\text{atm}} = \exp\left( - \int_0^L \left[ \alpha_{\text{Rayleigh}}(z) + \alpha_{\text{Mie}}(z) \right] dz \right)$$

At $\lambda = 810\text{ nm}$:

$\alpha_{\text{Rayleigh}} \approx 0.043\text{ dB/km}$ ($\approx 0.90\text{ dB}$ across $21\text{ km}$).
$\alpha_{\text{Mie}} \approx 0.086\text{ dB/km}$ under clear visibility conditions $>20\text{ km}$ ($\approx 1.80\text{ dB}$).
Accumulated $L_{\text{atm}} = 2.70\text{ dB}$.

3. Fiber Re-Coupling Loss ($L_{\text{coup}}$):

Coupling the collected, aberrated wavefront from a 0.6-meter aperture into a single-mode core of $d_{\text{core}} = 5.0\text{ }\mu\text{m}$ requires overlap matching between the incoming spatial mode $E_{\text{rx}}(r, \theta)$ and the fundamental fiber mode $\psi_0(r)$:

$$\eta_{\text{coup}} = \frac{\left| \int E_{\text{rx}}(r, \theta) \psi_0^*(r) r \, dr \, d\theta \right|^2}{\int |E_{\text{rx}}(r, \theta)|^2 r \, dr \, d\theta \cdot \int |\psi_0(r)|^2 r \, dr \, d\theta}$$

Under static lab conditions, $\eta_{\text{coup}} \approx 82\%$ ($-0.86\text{ dB}$). Under field turbulence corrected by the 2.0 kHz adaptive optics system, the dynamic mean coupling efficiency settles at $\eta_{\text{coup}} = 20.9\%$, yielding $L_{\text{coup}} = 6.80\text{ dB}$.

Overcoming Atmospheric Turbulence: The Adaptive Optics Architecture

The primary obstacle to maintaining unguided quantum entanglement through air across low altitudes is refractive index turbulence. Solar heating of asphalt, concrete, and varied canopy structures generates turbulent air cells (eddies) characterized by varying temperatures ($T$) and pressures ($P$).

                                 ATMOSPHERIC TURBULENCE CASCADE
 +------------------------------------------------------------------------------------+
 | Outer Scale (L0 ~ 10-100 m)  ==> Energy injected via thermal convection/wind       |
 |                                                                                    |
 | Inertial Subrange            ==> Kolmogorov cascade transfers turbulent kinetic   |
 |                                  energy to smaller eddies without dissipation      |
 |                                                                                    |
 | Inner Scale (l0 ~ 1-10 mm)   ==> Viscous dissipation into heat                     |
 +------------------------------------------------------------------------------------+

These eddies range from the outer scale ($L_0 \approx 10\text{ to }100\text{ m}$) down to the inner scale ($l_0 \approx 1\text{ to }10\text{ mm}$), altering the local refractive index $n$ according to:

$$n \approx 1 + 77.6 \times 10^{-6} \left( 1 + 7.52 \times 10^{-3} \lambda^{-2} \right) \frac{P}{T}$$

                TEMPORAL VARIATION OF TURBULENCE COEFFICIENT (Cn2) OVER 24 HOURS
  Cn2 (m^-2/3)
  10^-13 |                             /---------\  (Intense midday thermal turbulence)
  10^-14 |      /---\                 /           \                 /---\
  10^-15 | ----/     \---------------/             \---------------/     \----
  10^-16 |____________________________________________________________________
          00:00 03:00 06:00 09:00   12:00   15:00 18:00 21:00 24:00 (Local Time)

Across the 21.0 km Stony Brook–Brookhaven corridor, the refractive index structure parameter $C_n^2$ was continuously logged:

Nighttime Minimum (Quiescent): $C_n^2 = 1.2 \times 10^{-15}\text{ m}^{-2/3}$
Midday Maximum (High Thermal Convection): $C_n^2 = 3.8 \times 10^{-13}\text{ m}^{-2/3}$
Mean Operational Value: $C_n^2 = 1.8 \times 10^{-14}\text{ m}^{-2/3}$

The Fried Parameter ($r_0$) and Scintillation Index ($\sigma_I^2$)

The spatial coherence diameter of the atmosphere, known as the Fried parameter $r_0$, dictates the maximum effective telescope aperture usable before atmospheric turbulence degrades spatial resolution:

$$r_0 = \left( 0.423 \, k^2 \, C_n^2 \, L \right)^{-3/5}$$

Where wavenumber $k = \frac{2\pi}{\lambda} = \frac{2\pi}{810 \times 10^{-9}\text{ m}} \approx 7.757 \times 10^6\text{ m}^{-1}$.

Calculated Fried Parameter over 21 km path:
- At Night (Cn2 = 1.2 x 10^-15 m^-2/3):  r0 = 8.84 cm
- At Mean (Cn2 = 1.8 x 10^-14 m^-2/3):    r0 = 1.74 cm
- At Midday (Cn2 = 3.8 x 10^-13 m^-2/3):  r0 = 0.28 cm

Because the receiving aperture $D_{\text{rx}} = 60\text{ cm} \gg r_0$, the incoming optical wavefront breaks into dozens of distinct spatial speckles, preventing direct coupling into the $5.0\text{ }\mu\text{m}$ single-mode fiber core without active phase correction.

The intensity fluctuations (scintillation) are modeled via the Rytov variance $\sigma_R^2$ for a spherical wave:

$$\sigma_R^2 = 0.563 \, C_n^2 \, k^{7/6} \, L^{11/6}$$

At $L = 21\text{ km}$ and $C_n^2 = 1.8 \times 10^{-14}\text{ m}^{-2/3}$:

$$\sigma_R^2 = 0.563 \times (1.8 \times 10^{-14}) \times (7.757 \times 10^6)^{7/6} \times (21,000)^{11/6} \approx 14.2$$

This indicates a strong scintillation regime ($\sigma_R^2 > 1$), characterized by deep signal fading that requires sub-millisecond adaptive wavefront stabilization.

+------------------------------------------------------------------------------------+
|                CLOSED-LOOP ADAPTIVE OPTICS SERVO CONTROL DIAGRAM                   |
+------------------------------------------------------------------------------------+
|                                                                                    |
|  [Incoming Aberrated Wavefront]                                                    |
|                |                                                                   |
|                v                                                                   |
|   +----------------------------+                                                   |
|   | Fast Steering Mirror (FSM) | <---------- [Tip-Tilt Correction Signal @ 5 kHz]  |
|   +----------------------------+                               |                   |
|                |                                               |                   |
|                v                                               |                   |
|   +----------------------------+                               |                   |
|   | 241-Actuator Deformable M. | <---- [Wavefront Phase Correction @ 2 kHz]        |
|   +----------------------------+             |                 |                   |
|                |                             |                 |                   |
|                +-----------------------------+                 |                   |
|                |                                               |                   |
|                v                                               |                   |
|     [Dichroic Beam Splitter]                                   |                   |
|            |              |                                    |                   |
|            | (Beacon)     +--> [Shack-Hartmann Sensor] --------+                   |
|            |                   (128 Sub-Apertures)             |                   |
|            v                                                   |                   |
|    [Position-Sensitive Diode] ---------------------------------+                   |
|                                                                                    |
+------------------------------------------------------------------------------------+

Adaptive Optics Correction Performance

To cancel these distortions, the Brookhaven receiver incorporates an adaptive optics pipeline:

Tip-Tilt Tracking: An auxiliary continuous-wave beacon laser at $\lambda_{\text{beacon}} = 850\text{ nm}$ is co-propagated with the entangled photons. A position-sensitive diode reads beam centroid wander, driving a piezo-actuated Fast Steering Mirror (FSM) at $5.0\text{ kHz}$ to limit angular jitter below $0.45\text{ }\mu\text{rad}$.
High-Order Wavefront Correction: A Shack-Hartmann wavefront sensor comprising 128 microlenses samples the residual phase error. A 241-actuator deformable mirror applies inverse phase conjugations at $2.0\text{ kHz}$, exceeding the atmospheric Greenwood frequency:

$$f_G = 0.426 \left( \frac{k^2 C_n^2 v_{\text{wind}}^{5/3} L}{\cos(\zeta)} \right)^{3/5} \approx 620\text{ Hz}$$

(calculated at transverse wind velocity $v_{\text{wind}} = 8.5\text{ m/s}$).

By operating at more than three times the Greenwood frequency, the system maintained a Strehl ratio $S_r \ge 0.48$ throughout the trial, preserving fiber coupling efficiency and keeping link loss within operational limits.

Mitigating Daytime Noise: Spatial, Spectral, and Temporal Filtering

Achieving stable quantum entanglement through air during daylight hours requires isolating single-photon signals from solar background flux. Midday solar spectral radiance at sea level reaches $H_{\text{sun}} \approx 0.10\text{ W}/(\text{m}^2 \cdot \text{sr} \cdot \text{nm})$ at $810\text{ nm}$. Without filtering, a 0.6-meter collection aperture captures over $1.5 \times 10^9$ background photons per second, overwhelming single-photon detectors.

                    SOLAR BACKGROUND FILTERING ATTENUATION STAGES
 Incoming Solar Photons: ~1.5 x 10^9 counts/sec
   |
   +---> [Spatial Filter: 5.0 µm Core Pin-Hole Mode Selection]
   |     Remaining Noise: 3.2 x 10^6 counts/sec (99.78% Rejection)
   |
   +---> [Spectral Filter: 30 pm Fabry-Perot Etalon Cavity]
   |     Remaining Noise: 4.8 x 10^3 counts/sec (99.85% Rejection)
   |
   +---> [Temporal Filter: 350 ps Coincidence Gating Window]
         Remaining Noise: 0.082 counts/sec Effective Coincidence Noise

===================================================================================
DAYLIGHT NOISE REJECTION MATRIX (810 NM CHANNEL, 100,000 LUX AMBIENT ILLUMINANCE)
===================================================================================
Filtering Stage            Mechanism Applied                   Noise Suppression
-----------------------------------------------------------------------------------
1. Spatial Filtering       $5.0\text{ }\mu\text{m}$ SMF spatial mode selection       $26.7\text{ dB}$ ($470\times$)
2. Spectral Filtering      Fabry-Perot Cavity ($\Delta \lambda = 30\text{ pm}$)     $28.2\text{ dB}$ ($660\times$)
3. Temporal Gating         Coincidence gating ($\Delta t = 350\text{ ps}$)          $34.5\text{ dB}$ ($2,850\times$)
4. Polarization Filtering  Cross-polarization rejection prisms                      $12.0\text{ dB}$ ($16\times$)
-----------------------------------------------------------------------------------
TOTAL BACKGROUND PHOTON REJECTION EFFICIENCY                               101.4 dB
===================================================================================

+------------------------------------------------------------------------------------+
|            TEMPORAL COINCIDENCE DISCRIMINATION (PICOSECOND GATING)                 |
+------------------------------------------------------------------------------------+
|                                                                                    |
|   Channel A (Local Photon Arrival Time)                                            |
|   ---|--------|-----------------|--------|-------------------|--------|--------->  |
|      t1       t2                t3       t4                  t5       t6           |
|                                                                                    |
|   Channel B (21 km Remote Arrival Time, Synchronized via GPS-Disciplined Rubidium) |
|   ---|--------------------------|----------------------------|------------------>  |
|      t1'                        t3'                          t5'                   |
|                                                                                    |
|   Coincidence Window: Delta_t = |t_A - t_B| <= 350 ps                              |
|   - Real Pair Event (t1-t1'): Correlated Entangled Detection  ==> REGISTERED       |
|   - Ambient Solar Photon (t2): Outside Gating Window         ==> DISCARDED        |
|   - Real Pair Event (t3-t3'): Correlated Entangled Detection  ==> REGISTERED       |
|   - Dark Count Noise (t4): Outside Gating Window             ==> DISCARDED        |
+------------------------------------------------------------------------------------+

Quantitative Signal-to-Noise Ratio (SNR) Derivation

The background noise count rate $R_{\text{noise}}$ passing through the system is determined by:

$$R_{\text{noise}} = H_{\text{sun}} \cdot \Omega_{\text{fov}} \cdot A_{\text{rx}} \cdot \Delta \lambda \cdot \eta_{\text{rx}} \cdot \eta_{\text{det}}$$

Where:

Receiver field of view solid angle: $\Omega_{\text{fov}} = \pi \left( \frac{\theta_{\text{fov}}}{2} \right)^2 \approx \pi (1.95 \times 10^{-6}\text{ rad})^2 \approx 1.19 \times 10^{-11}\text{ sr}$.
Collecting area: $A_{\text{rx}} = \pi (0.30\text{ m})^2 \approx 0.283\text{ m}^2$.
Optical bandwidth: $\Delta \lambda = 0.030\text{ nm}$ ($30\text{ pm}$).
Transmittance: $\eta_{\text{rx}} = 0.25$, detector efficiency: $\eta_{\text{det}} = 0.95$.

$$R_{\text{noise}} = 0.10 \times (1.19 \times 10^{-11}) \times 0.283 \times 0.030 \times 0.25 \times 0.95 \times \left( \frac{1\text{ photon}}{2.45 \times 10^{-19}\text{ J}} \right) \approx 4.88 \times 10^3\text{ noise counts/sec}$$

The coincidence accidental rate $R_{\text{acc}}$ inside a time window $\Delta t = 350\text{ ps}$ is:

$$R_{\text{acc}} = R_{\text{local}} \cdot R_{\text{remote}} \cdot \Delta t$$

With local detection rate $R_{\text{local}} = 1.8 \times 10^5\text{ counts/s}$ and remote rate $R_{\text{remote}} = (4.22 \times 10^3\text{ signal}) + (4.88 \times 10^3\text{ noise}) = 9.10 \times 10^3\text{ counts/s}$:

$$R_{\text{acc}} = (1.8 \times 10^5) \times (9.10 \times 10^3) \times (350 \times 10^{-12}) \approx 0.573\text{ accidental coincidences/sec}$$

Given the detected entangled coincidence rate $R_{\text{coinc}} = 4.22 \times 10^3\text{ pairs/sec} \times 10^{-2.85} \approx 5.96\text{ true pairs/sec}$, the coincidence-to-accidental ratio (CAR) is:

$$\text{CAR} = \frac{R_{\text{coinc}}}{R_{\text{acc}}} = \frac{5.96}{0.573} \approx 10.4$$

A CAR exceeding 10 confirms that quantum correlations remain distinguishable from ambient daylight noise, allowing continuous daytime operation without blinding the single-photon array.

Experimental Validation: Density Matrix and Bell Inequality Violations

To verify the integrity of the distributed quantum states, the research team performed complete Two-Qubit Quantum State Tomography (QST) across the 21.0 km link.

+------------------------------------------------------------------------------------+
|             RECONSTRUCTED DENSITY MATRIX REAL COMPONENTS Re(rho)                   |
+------------------------------------------------------------------------------------+
|                                                                                    |
|             |HH>             |HV>             |VH>             |VV>                |
|      +----------------+----------------+----------------+----------------+         |
| <HH| |    +0.024      |    -0.011      |    +0.008      |    -0.015      |         |
|      +----------------+----------------+----------------+----------------+         |
| <HV| |    -0.011      |    +0.478      |    -0.456      |    +0.012      |         |
|      +----------------+----------------+----------------+----------------+         |
| <VH| |    +0.008      |    -0.456      |    +0.474      |    -0.009      |         |
|      +----------------+----------------+----------------+----------------+         |
| <VV| |    -0.015      |    +0.012      |    -0.009      |    +0.024      |         |
|      +----------------+----------------+----------------+----------------+         |
|                                                                                    |
| Target Singlet State |Psi-> = (|HV> - |VH>) / sqrt(2)                              |
| Measured Fidelity: F = <Psi-| rho |Psi-> = 94.2% +/- 1.1%                          |
| Purity: Tr(rho^2) = 89.8% | Concurrence: C(rho) = 0.884                            |
+------------------------------------------------------------------------------------+

Using 16 independent projective measurement combinations across horizontal ($H$), vertical ($V$), diagonal ($D = \frac{H+V}{\sqrt{2}}$), anti-diagonal ($A = \frac{H-V}{\sqrt{2}}$), right-circular ($R = \frac{H-iV}{\sqrt{2}}$), and left-circular ($L = \frac{H+iV}{\sqrt{2}}$) polarization bases, the density matrix $\rho$ was reconstructed via maximum-likelihood estimation:

Reconstructed Density Matrix Parameters:
- State Fidelity: F = <Psi-| rho |Psi-> = 0.942 +/- 0.011
- State Purity: gamma = Tr(rho^2) = 0.898 +/- 0.014
- Concurrence: C(rho) = 0.884 +/- 0.018
- Entanglement of Formation: E_F = 0.841 e-bits

               DIURNAL QUANTUM BIT ERROR RATE (QBER) STABILITY OVER 24 HOURS
  QBER (%)
  12.0 | - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - (Shor-Preskill Limit: 11.0%)
  10.0 |
   8.0 |
   6.0 |
   4.0 |                  /-----\  (Solar Noon Peak: 3.24%)
   2.0 |  /---\          /       \          /---\
   0.0 +--+---+----------+-------+----------+---+--->
         00:00 04:00   08:00   12:00      16:00 20:00 24:00 (Time of Day)
         [Night: 1.12%] [Morning: 1.85%] [Dusk: 1.45%]

Clauser-Horne-Shimony-Holt (CHSH) Test

The Bell test correlation function was evaluated across measurement angles $a = 0^\circ, a' = 45^\circ$ for Alice (Stony Brook) and $b = 22.5^\circ, b' = 67.5^\circ$ for Bob (Brookhaven):

$$S = |E(a, b) - E(a, b') + E(a', b) + E(a', b')|$$

Correlation Outcomes:
- E(a, b)   = +0.678 +/- 0.019
- E(a, b')  = -0.669 +/- 0.020
- E(a', b)  = +0.672 +/- 0.018
- E(a', b') = +0.665 +/- 0.019

Calculated CHSH Parameter:
S = |0.678 - (-0.669) + 0.672 + 0.665| = 2.684 +/- 0.038

The observed value exceeds the local realism threshold ($S \le 2$) by 18 standard deviations, confirming the preservation of non-local quantum states despite atmospheric turbulence, building sway, and daylight background radiation.

Comparative Global Free-Space Entanglement Deployments

The Stony Brook–Brookhaven transmission builds upon a lineage of urban and long-distance free-space experiments, highlighting the technological trajectory of terrestrial quantum links.

                     TRANSMISSION DISTANCE VS. INSERTION LOSS BENCHMARKS
  Channel Loss (dB)
   60 |                                                  [Hefei-Huailai (51.4 dB, 7.0 km)]
   50 |
   40 |                                     [Rome Urban (32.0 dB, 0.27 km)]
   30 |                 [Stony Brook-BNL (28.5 dB, 21.0 km)]
   20 |    [Vienna Danube (18.2 dB, 1.4 km)]
   10 |
    0 +----+------------+-------------------+------------+----------------------------+--->
          0.1 km       1.0 km              10.0 km      50.0 km                     100.0 km

========================================================================================================================
COMPARATIVE MATRIX OF FREE-SPACE QUANTUM OPTICAL NETWORK EXPERIMENTS
========================================================================================================================
Project / Location           Distance   Encoding Scheme        Wavelength   Loss (dB)  Fidelity / Visibility  Key Rate / Coinc
------------------------------------------------------------------------------------------------------------------------
BNL–Stony Brook (2026)     21.0 km    Polarization (Bell)    810 nm       28.5 dB    F = 94.2%              4.22 kbps
Sapienza Rome (2026)       0.27 km    Time-Evolving QD       850 nm       32.0 dB    V = 91.5% (<50ps sync)  1.80 kbps
Ottawa Intra-City (2025)   5.4 km     High-Dim OAM (d=4)     780 nm       26.0 dB    V = 88.4%              12.4 kbps
Jinan Terrestrial (2025)       7.0 km     Time-Bin BB84          1550 nm      51.4 dB    QBER = 0.87%           134 bps
Vienna Intra-City (Classic) 7.8 km     Polarization SPDC      810 nm       35.0 dB    S = 2.27 (14-sigma)    25 Hz coinc
Munich Urban Link            1.4 km     Phase-Encoded QKD      1550 nm      24.0 dB    V = 99.07%             4.22 kbps
Canary Islands (Inter-Island) 143.0 km   Polarization SPDC      810 nm       55.0 dB    S = 2.45 (Loopholes)   0.3 Hz coinc
========================================================================================================================

High-Dimensional State Encoding (OAM and Time-Bins)

While the Long Island link utilized two-dimensional polarization qubits ($d = 2$), parallel urban programs are transitioning to high-dimensional quantum systems (qudits, $d > 2$). At the University of Ottawa, a 5.4 km urban link employed Orbital Angular Momentum (OAM) states carrying helical phase fronts $\exp(i \ell \phi)$ with topological charges $\ell = \pm 1, \pm 2$ ($d = 4$):

Information Density: Each detected photon carries $\log_2(d) = \log_2(4) = 2\text{ classical bits of information}$.
Noise Resilience: The theoretical error threshold required to compromise security increases from $11.0\%$ for two-dimensional states up to $18.9\%$ for $d=4$ and $25.0\%$ for $d=8$, offering enhanced resilience against urban atmospheric turbulence.

Economic and Infrastructure Analysis: Free-Space vs. Dark Fiber

Deploying quantum networks across dense metropolitan centers faces substantial economic and physical obstacles when relying entirely on subterranean fiber-optic cables.

       CAPITAL EXPENDITURE (CAPEX) COMPARISON PER MILE: URBAN FIBER VS. FREE-SPACE
 Cost ($ USD / Mile)
 $1,200,000 |  +--------------------------------------------------------+
            |  | Dense Metro Trenching (Manhattan/London/Tokyo Core)   | = $1,200,000/mile
   $800,000 |  +--------------------------------------------------------+
            |  | Suburban Fiber Conduit Construction                   | = $350,000/mile
   $400,000 |  +--------------------------------------------------------+
            |  | Leased Dark Fiber (10-Year Amortized CapEx Equivalent)| = $180,000/mile
    $15,000 |  | Free-Space Optical Rooftop Terminal Installation      | = $15,300/mile (Over 13 mi)
         $0 +--+--------------------------------------------------------+------------------->

===================================================================================
FINANCIAL AND DEPLOYMENT COST MODEL: METROPOLITAN QUANTUM NETWORK (20-MILE SPAN)
===================================================================================
Cost Element                        Subterranean Dark Fiber        Free-Space Optical Node Pair
-----------------------------------------------------------------------------------
1. Initial Capital Deployment       $7,000,000 ($350k/mi trenching) $280,000 (Two optical terminals)
2. Right-of-Way / Permitting Fees   $450,000 (Municipal permits)   $25,000 (Rooftop lease fees)
3. Specialized Optical Hardware     $120,000 (Standard interfaces) $180,000 (Adaptive optics + AO telescope)
4. Annual Maintenance / Leasing     $72,000/year (Fiber lease/rep) $14,000/year (Window cleaning & calib)
5. Physical Deployment Timeline     14 to 26 Months                3 to 6 Weeks
-----------------------------------------------------------------------------------
TOTAL 10-YEAR EXPENDITURE TCO       $8,290,000                     $625,000
===================================================================================

Deploying quantum entanglement through air reduces the 10-year Total Cost of Ownership (TCO) by roughly 92% compared to new civil engineering trenching projects, while shortening deployment timelines from years to weeks.

+------------------------------------------------------------------------------------+
|                  HYBRID QUANTUM METROPOLITAN ROUTING TOPOLOGY                      |
+------------------------------------------------------------------------------------+
|                                                                                    |
|                     [Central Quantum Server / Mainframe Node]                      |
|                                         |                                          |
|                    +--------------------+--------------------+                     |
|                    |                                         |                     |
|         (Underground Fiber Channel)               (Free-Space Line-of-Sight)       |
|         - 1550 nm Telecom C-Band                  - 780 nm / 810 nm Native Atoms   |
|         - Fixed Point-to-Point                    - Dynamic Beam Steering          |
|                    |                                         |                     |
|                    v                                         v                     |
|         [Subterranean Data Center]                [Rooftop Quantum Watchtower]     |
|                                                              |                     |
|                                                 ========================           |
|                                                 21.0 km Free-Space Link            |
|                                                 ========================           |
|                                                              |                     |
|                                                              v                     |
|                                                   [Island / Marine Node]           |
|                                                   (Yale / Long Island Sound)       |
|                                                                                    |
+------------------------------------------------------------------------------------+

Integration with Quantum Repeaters and Atomic Memories

A permanent free-space link provides a foundational testbed for integrating quantum memories and repeater nodes. Direct optical transmission without repeaters remains constrained by link losses exceeding $45\text{ to }50\text{ dB}$, where coincidence rates drop below single-photon detector dark count thresholds.

+------------------------------------------------------------------------------------+
|             QUANTUM REPEATER NODE WITH ROOM-TEMPERATURE ATOMIC MEMORY              |
+------------------------------------------------------------------------------------+
|                                                                                    |
|  Incoming 810 nm Free-Space Photon                                                 |
|               |                                                                    |
|               v                                                                    |
|  +-------------------------------------+                                           |
|  | Electromagnetically Induced         |                                           |
|  | Transparency (EIT) Vapor Cell       | <--- [795 nm Control Laser Pulse]         |
|  | (Rubidium-87 Vapor @ 45.0 °C)        |                                           |
|  +-------------------------------------+                                           |
|               |                                                                    |
|               | (Photon Stored as Collective Atomic Spin Wave)                     |
|               | Storage Time: T1 = 1.25 ms | Retrieval Fidelity: F = 91.4%         |
|               |                                                                    |
|               v                                                                    |
|  +-------------------------------------+                                           |
|  | Retrieved 810 nm Single Photon      | ---> [Bell State Measurement (BSM)]       |
|  +-------------------------------------+               ^                           |
|                                                        |                           |
|  Entangled Idler from Adjacent Node -------------------+                           |
|                                                                                    |
+------------------------------------------------------------------------------------+

At Stony Brook University, researchers have integrated the free-space transceiver with room-temperature rubidium vapor cells utilizing Electromagnetically Induced Transparency (EIT):

Storage Medium: Caesium ($^{133}\text{Cs}$) or Rubidium ($^{87}\text{Rb}$) atomic vapor operating at $45.0^\circ\text{C}$ (eliminating cryogenic cooling systems).
Optical Storage Coherence Time ($T_1$): $1.25\text{ milliseconds}$, sufficient to buffer photonic qubits during round-trip classical acknowledgment signaling over distances up to $375\text{ kilometers}$.
Storage and Retrieval Fidelity: $91.4\% \pm 0.8\%$ for single-photon polarization qubits.
Entanglement Swapping Feasibility: By buffering photons arriving from the 21 km free-space link and performing a Bell State Measurement (BSM) against an adjacent node, entanglement can be extended across cascaded free-space segments without direct end-to-end photon transmission.

Scaling to Maritime and Spaceborne Channels: 2026–2035 Roadmap

The validation of the Stony Brook–Brookhaven optical channel initiates a multi-phase deployment roadmap extending across regional water barriers and into orbit.

+------------------------------------------------------------------------------------+
|                 TRI-NODE REGIONAL FREE-SPACE QUANTUM TOPOGRAPHY                    |
+------------------------------------------------------------------------------------+
|                                                                                    |
|                                [Yale University]                                   |
|                           (New Haven, CT / Node 3)                                 |
|                                      ^                                             |
|                                     /                                              |
|                                    /                                               |
|                    48.0 km (30 mi)/ Cross-Sound Maritime Channel                   |
|                    Marine Boundary Layer (Cn2 = 4.5 x 10^-15 m^-2/3)               |
|                                  /                                                 |
|                                 v                                                  |
|   [Stony Brook Watchtower] <========================> [BNL Quantum Lighthouse]     |
|   (Stony Brook, NY / Node 1)  21.0 km Terrestrial     (Upton, NY / Node 2)         |
|                               Suburban Link                                        |
|                                                                                    |
+------------------------------------------------------------------------------------+

The Trans-Sound Leg: Stony Brook to Yale University

The completed next phase extends the network from Stony Brook's Quantum Watchtower across the Long Island Sound to a newly completed optical station at Yale University in New Haven, Connecticut:

Span Distance: $48.0\text{ kilometers}$ ($29.8\text{ miles}$) across open marine water.
Atmospheric Boundary Layer: Marine boundary layers exhibit lower thermal convection than overland paths ($C_n^2 \approx 4.5 \times 10^{-15}\text{ m}^{-2/3}$), increasing the Fried parameter to $r_0 \approx 5.2\text{ cm}$ despite the doubled link length.
Projected Total Channel Loss: $34.2\text{ dB}$, targeting an operational key generation rate $> 850\text{ bits/second}$.

========================================================================================================================
STRATEGIC ROADMAP: FREE-SPACE AND HYBRID QUANTUM NETWORK INFRASTRUCTURE (2026–2035)
========================================================================================================================
Phase / Milestone   Target Timeline   Channel Configuration   Target Distance   Projected Key Rate  Primary Hardware Milestone
------------------------------------------------------------------------------------------------------------------------
**Phase 1: Metro**   2026–2027         Terrestrial FSO Link    21 km to 48 km    1.0 to 5.0 kbps     Adaptive Optics + 0.6m Telescopes
**Phase 2: Regional**2027–2029         Over-Water + Fiber Mesh 150 km Regional   10.0 to 50.0 kbps   Warm Vapor Quantum Memory Repeater
**Phase 3: Space LEO**2029–2032         LEO Satellite-to-Ground 500 km to 1200 km 50.0 to 200.0 kbps 100 GHz Integrated Modulators
**Phase 4: Global**  2032–2035         Intercontinental Const. > 5,000 km Global 1.0 to 10.0 Mbps    Spaceborne Quantum Repeaters
========================================================================================================================

           PROJECTED GLOBAL FREE-SPACE QUANTUM NETWORK THROUGHPUT SCALING
  Distilled Key Rate (bps)
  10^7 |                                                      [2035: Global Satellite-Repeater Mesh]
  10^6 |
  10^5 |                                      [2030: LEO Satellite Downlinks (50-200 kbps)]
  10^4 |                  [2028: Regional Memory Repeaters (10-50 kbps)]
  10^3 | [2026: Terrestrial FSO (4.22 kbps)]
  10^2 |______________________________________________________________________________
        2026                 2028                 2030                 2032                 2035

Orbital Interconnects and Global Reach

Because the effective thickness of Earth's atmospheric boundary layer is roughly 6 to 8 kilometers when pointing vertically, the 21.0-kilometer horizontal ground-level link spans an atmospheric mass equivalent to more than 2.5 vertical atmospheric thicknesses.

Operating through 21.0 kilometers of ground turbulence confirms that optical tracking and adaptive optics systems can maintain quantum channels with Low Earth Orbit (LEO) satellites at 500-kilometer altitudes, where over 95% of the propagation path occurs in the near-vacuum of space.

+------------------------------------------------------------------------------------+
|              ATMOSPHERIC DENSITY COMPARISON: HORIZONTAL VS. SATELLITE              |
+------------------------------------------------------------------------------------+
|                                                                                    |
|  1. Horizontal Terrestrial Path (Stony Brook to BNL):                              |
|     [Transmitter] ================== 21.0 km ==================> [Receiver]        |
|     - Constant Pressure: P = 1.0 atm throughout entire 21.0 km span                |
|     - Total Integrated Air Mass = 21.0 km-atm equivalent                           |
|                                                                                    |
|  2. Vertical / Slant Satellite Uplink/Downlink (LEO at 500 km, 45° Elevation):     |
|     [Ground Station] -------- 8 km Atmosphere --------> [Vacuum of Space] -> [LEO] |
|     - Exponentially decaying pressure: P(z) = P0 * exp(-z / 8.4 km)                |
|     - Total Integrated Air Mass = ~8.4 / sin(45°) = 11.8 km-atm equivalent         |
|                                                                                    |
|  CONCLUSION: Transmitting quantum states across 21 km of ground air imposes        |
|  nearly DOUBLE the optical turbulence and scattering load of a satellite downlink. |
+------------------------------------------------------------------------------------+

The successful transmission between Stony Brook and Brookhaven demonstrates that unguided optical channels can support high-fidelity quantum entanglement across metropolitan landscapes.

Operating with a Bell parameter $S = 2.684$, an average QBER of 2.38%, and sustained resilience against solar and atmospheric noise, this link provides a verified blueprint for building scalable, multi-node quantum networks without trenching fiber infrastructure.