On the simultaneous scattering of two photons by a single two-level atom

The interaction of light with a single two-level emitter is the most fundamental process in quantum optics, and is key to many quantum applications. As a distinctive feature, two photons are never detected simultaneously in the light scattered by the emitter. This is commonly interpreted by saying that a single two-level quantum emitter can only absorb and emit single photons. However, it has been theoretically proposed that the photon anticorrelations can be thought of as arising from quantum interference between two possible two-photon scattering amplitudes, which one refers to as coherent and incoherent. This picture is in stark contrast to the aforementioned one, in that it assumes that the atom has two different mechanisms at its disposal to scatter two photons at the same time. Here we experimentally validate the interference picture by showing that, when spectrally rejecting only the coherent component of the fluorescence light of a single two-level atom, the remaining light consists of photon pairs that have been simultaneously scattered by the atom. Our results offer fundamental insights into the quantum-mechanical interaction between light and matter and open up novel approaches for the generation of highly non-classical light fields enabling, for example, Fourier-limited photon-pair sources that approach the theoretical limit in brightness. The incoherent component of the fluorescence from a single two-level atom is investigated after rejecting the coherent component. Contrary to intuition, its photon statistics experimentally shows strong photon bunching. This result suggests that the atom does in fact simultaneously scatter two photons.

The interaction of light with a single two-level emitter is the most fundamental process in quantum optics, and is key to many quantum applications.As a distinctive feature, two photons are never detected simultaneously in the light scattered by the emitter.This is commonly interpreted by saying that a single two-level quantum emitter can only absorb and emit single photons.However, it has been theoretically proposed that the photon anticorrelations can be thought of as arising from quantum interference between two possible two-photon scattering amplitudes, which one refers to as coherent and incoherent.This picture is in stark contrast to the aforementioned one, in that it assumes that the atom has two different mechanisms at its disposal to scatter two photons at the same time.Here we experimentally validate the interference picture by showing that, when spectrally rejecting only the coherent component of the fluorescence light of a single two-level atom, the remaining light consists of photon pairs that have been simultaneously scattered by the atom.Our results offer fundamental insights into the quantum-mechanical interaction between light and matter and open up novel approaches for the generation of highly non-classical light fields enabling, for example, Fourier-limited photon-pair sources that approach the theoretical limit in brightness.
The interaction of a single two-level quantum emitter with a nearresonant coherent light field is one of the cornerstones of quantum optics and lies at the heart of many modern experiments and applications in this field [1][2][3][4][5][6] .The quantum-mechanical description of this interaction 7,8 shows that the scattered light exhibits photon antibunching; that is, it never contains two photons at the same time and place 9,10 .This property can be seen as a consequence of the photon emission being associated with quantum jumps from the emitter's excited state to its ground state.It thus seems natural to state that a single two-level atom will never scatter two photons simultaneously.However, a more in-depth inspection of the optical Bloch equations reveals the existence of two distinct components of the scattered field, referred to as coherently and incoherently scattered light, reflecting their respective ability and inability to interfere with the driving field.Taken together, these components form the well-known Mollow structure in the fluorescence spectrum 11 .Interestingly, when considered individually, each of these components also contains higher photon-number components, that is, two or more photons at the same time and place.In view of this fact, it has been argued that the origin of antibunching in resonance fluorescence stems from the destructive interference between the higher photon-number components of the coherently and incoherently scattered light 12 .Cross-correlation measurements between the sidebands of the Mollow structure have revealed a time ordering between the detected photons, whereby the photon originating from the sideband furthest from the emitter's resonance typically arrives first [13][14][15] .This is commonly interpreted as a cascaded two-photon emission from the dressed states of the emitter [16][17][18] .Recently, motivated by theoretical works [19][20][21] , it has been experimentally shown that by spectrally rejecting individually, both the coherent and incoherent components have a non-zero probability for containing two simultaneous photons.Now, the two components exhibit different spectra, which for the coherently scattered component is given by a delta function at ω L , whereas for the incoherently scattered component, at low saturation, the spectrum is approximately given by a pair of Lorentzian lines, each of width γ, which are separated by ±Δ from ω L (ref.8; Fig. 1).Their distinct spectra allow us to separate the components by applying a frequency notch filter centred at the driving laser frequency.For this, we first collect the atomic fluorescence with a high-numerical-aperture (NA) objective before coupling it into a single-mode fibre, which guides the light to the filter (Fig. 2).The latter is realised by means of a fibre-ring resonator, where we use a low-loss variable-ratio fibre beamsplitter to control the coupling rate into the resonator, κ ext , and in turn the on-resonance transmission past the resonator (Fig. 2, right inset).The resonator is stabilized to resonance with the driving laser and has an extremely narrow linewidth that is adjustable between 1 and 15 MHz, which is much smaller than the separation of the coherent and incoherent components of |Δ|.As a consequence, the filter reduces the coherent component by the κ ext -dependent complex-valued field transmission factor, t coh (κ ext ), but does not substantially alter the incoherent component (Methods).To measure the modified photon statistics after filtering, we send the light to a Hanbury Brown and Twiss set-up consisting of a 50/50 beamsplitter and a single-photon counter in each of its outputs 29 , and record the second-order correlation function, , for different settings of the coupling parameter κ ext .Figure 3 showcases the second-order correlation functions as κ ext is increased.For κ ext = 0, the filter resonator has no effect on the collected atomic fluorescence.The resulting measured second-order correlation function (purple) is g (2) (0) = 0.43 ± 0.05, which is compatible with perfect photon antibunching when considering the background count rates of our photon counters.This illustrates the expected behaviour whereby, in the absence of spectral selection, the light scattered by a single atom never contains more than one photon at the same time the incoherently scattered component of the fluorescence light of a single quantum dot, one can modify the photon statistics of the scattered field in such a way that all correlations are lost, and the remaining light again recovers the spectral and temporal characteristics of the classical driving field 22,23 .In this Article we show that when performing the opposite, that is, rejecting the coherent component via spectral filtering, the remaining incoherent component of the fluorescence light of a single two-level atom consists of pairs of photons that appear to have been simultaneously scattered by the atom.Our results validate the picture that antibunching in resonance fluorescence arises from fully destructive interference between the two-photon components of the coherently and incoherently scattered light.
In more detail, we use a weak and detuned light field to excite the atom such that the coherent and incoherent components are spectrally separated while also ensuring that the latter consists purely of photon pairs (Fig. 1).To collect only the incoherently scattered light, we make use of a novel adjustable, narrowband optical notch filter to continuously reduce the amplitude of the coherently scattered light.Subsequently, we measure the second-order correlation function, g (2) (τ), of the residually transmitted light.By tuning the relative magnitude of the two scattered components, we observe an evolution from a clear photon antibunching of g (2) (0) = 0.43 ± 0.05 without filtering to a strong photon bunching of up to g (2) (0) = 7.65 ± 1.21 when maximally rejecting the coherently scattered light.The latter indicates that, counterintuitively, the incoherently scattered light only consists of pairs of simultaneously scattered photons.
In our experiment we prepare a single 85 Rb atom in an optical dipole trap that is loaded from a magneto optical trap (MOT).The MOT lasers, with frequency ω L , are red-detuned with respect to the Stark-shifted atomic resonance by Δ/2π = −57.9± 3.7 MHz.In this setting, the atom scatters photons from the MOT laser beams into free space.The quantum state of this fluorescence light can be separated into a coherent state, |α〉, and an incoherently scattered component, |ϕ〉, whereby the latter originates from the saturable nature of the emitter.We weakly drive the atom with a saturation parameter S = 0.025 ± 0.004.In this low-saturation regime, the incoherently scattered part consists solely of photon pairs [24][25][26][27] , with a probability amplitude for finding two photons within a time delay τ given by ϕ (2) This wavefunction is the temporal representation of an energy-time entangled state, whereby the photons jointly propagate within a time interval ~1/γ and exhibit perfect frequency correlations.The corresponding two-photon amplitude for the coherent component is given by the time-independent value Details are provided in the Methods.Here, n coh = γS/(S + 1) 2 is the photon scattering rate of the coherent component, and γ/2π = 3.0 MHz (ref.28) is the amplitude decay rate of the D 2 transition of 85 Rb.The total scattered field is given by the sum of the two components |ψ〉 = |α〉 + |ϕ〉 and for its two-photon amplitude one obtains Importantly, at zero time delay (τ = 0), the coherently and incoherently scattered components have equal amplitude but are π out of phase.This results in perfect destructive interference such that the light never contains two simultaneously scattered photons.Despite this fact, taken |α (2) (ω)| ~|ф (2) (ω)| g (2) (τ)  (2) (τ = 0) = 0.For |τ| > 0, Rabi oscillations are apparent.e, Temporal two-photon wavefunctions of the coherent and incoherent component, α (2) (τ) and ϕ (2) (τ), respectively, where ϕ (2) (0) = −α (2) (τ).f, Two-photon wavefunctions in the frequency domain.The coherent component, α (2) (ω), is a delta function centred at ω L , and the incoherent component, φ (2) (ω), is given by a double Lorentzian at ω L ± Δ.

Article
https://doi.org/10.1038/s41566-023-01260-7 and place.For |τ| > 0, we observe a damped oscillatory behaviour, originating from the driven atom undergoing Rabi oscillations at an effective frequency of Ω eff ≈ Δ, which are dampened on the timescale of the excited-state lifetime, 1/2γ.This oscillatory behaviour can also be considered to originate from interference between the coherently and incoherently scattered two-photon components, (equation ( 3)).
When increasing κ ext , the power transmission, |t coh (κ ext )| 2 , of the coherent component through the filter is reduced.Consequently, we upset the delicate balance of the two components, resulting in less destructive interference and hence increasing g (2) (τ = 0) (blue and green datasets).Around κ ext ≈ κ 0 , the coherent component is strongly attenuated.Consequently, the second-order correlation function resembles a double-exponential envelope centred at τ = 0, featuring strong photon bunching of up to g (2) (τ = 0) = 7.65 ± 1.21 (yellow and orange datasets).When further increasing κ ext , the transmission, t coh (κ ext ), of the coherent component will increase again, and the bunching decreases (red and violet datasets).For very large values of κ ext , antibunching would be recovered, although it is not observed here for technical reasons.
Figure 4 shows the evolution of the measured second-order correlation function at τ = 0, as a function of κ ext .For each setting, the value of g (2) (0, κ ext ) is obtained by averaging the data shown in Fig. 3 in a 3-ns window centred around τ = 0.The solid orange curve is the theoretical prediction of g (2) (0, κ ext ) (Methods).We note that our model takes into account residual drifts of the filter resonance, which we obtain from fits to the datasets in Fig. 3 (Methods).As a consequence of these drifts, maximum suppression of the coherently scattered component is reached for a setting of κ ext slightly larger than κ 0 .Across all settings of our filter, we observe a beating in the second-order correlation function resulting from the interference between the incoherently and remaining coherently scattered two-photon components, once again illustrating their coherence with respect to each other.For perfect suppression of the coherent component, these oscillations vanish and the g (2) (τ) would exhibit a pure double-exponential decay centred at τ = 0, with a decay time of 1/2γ.
From the two datasets exhibiting the best suppression of the coherently scattered light (κ ext = 1.00κ 0 and κ ext = 2.23κ 0 ), we obtain a rate of detected photon pairs of n (2)   meas.= (4.91 ± 0.47) × 10 −3 s −1 (Methods).Given the pair detection efficiency η 2 /2 of our set-up, this corresponds to a total rate of photon pairs incoherently scattered by the atom of n (2) total = 5.3 ± 0.5 kHz .This agrees well with the rate of incoherently scattered photon pairs, n inc /2 = 5.7 ± 1.4 kHz, that is expected for our saturation parameter (Methods).Note, that n inc is much smaller than the rate of coherently scattered photons, n coh = 455 ± 55 kHz, and, as the ratio of these two components is given by S, for a vanishing saturation, the relative contribution of the incoherent component tends to zero.Yet, intriguingly, when interfering with the coherent component, its small but finite value is still necessary for the observation of photon antibunching.
To conclude, we return to the initial discussion on whether a single two-level atom will simultaneously scatter two photons.As is often the case in quantum mechanics, the answer is both yes and no, depending on which picture is considered.Indeed, the unmodified light field emitted by the atom never contains two simultaneous photons.Here, we have highlighted the picture where this is a consequence of the atom continuously scattering the incident two-photon component in two different ways-coherently and incoherently-which perfectly destructively interfere to yield antibunched photon statistics.An alternative picture is to state that the atom only scatters photons one by one, and that the emitted stream of antibunched photons is transformed into bunched photons by the passive linear notch filter in our set-up.In this picture, the resonator acts as a storage medium for the scattered light, and Hong-Ou-Mandel quantum interference 30   incoupling beamsplitter between the incoming fluorescence light and the stored field results in photon bunching of the light transmitted past the resonator (Methods).Beyond these fundamental considerations, the demonstrated effect lends itself to realising spectrally narrowband photon-pair sources that are inherently compatible with optical transitions in quantum emitters 31,32 (Methods).In particular, for high photon collection efficiencies [33][34][35] , the achievable photon-pair rate would approach the theoretical limit in brightness of a Fourier-limited pair source.Consequently, we envision that photon-pair sources based on this effect will emerge as a key resource in optical quantum technologies and quantum information processing.

Theoretical description
Wavefunction description of the scattered field.In a standard picture, the atom-light interaction is described by optical Bloch equations, in which the atom is excited by a coherent laser beam and the scattered quantum field can be calculated from the time evolution of the atomic raising and lowering operators 8 .In the steady state, the scattered field exhibits two distinct components, usually referred to as coherently and incoherently scattered light, and their ratio depends on the saturation parameter of the excitation laser field (Ω/γ) Here, Δ = ω L − ω 0 is the laser-atom detuning with laser frequency ω L and atomic resonance frequency ω 0 .Ω is the Rabi frequency of the driving field and γ is the amplitude decay rate of the excited state.The coherent and incoherent scattering rates of the atom are given by and respectively, where the right-hand-side expressions are the low saturation expansions up to second order in S, which we apply throughout the manuscript to obtain a simple theoretical description of the scattered light.In the following, we consider the case where we drive the atom with an excitation laser beam for time duration T, where T is chosen such that it is much larger than the atomic decay time, T ≫ 1/γ (to neglect transient features), and small enough that the mean photon number is much smaller than one, that is, T ≪ 1/γS.We note that in the low saturation regime, S ≪ 1, this is always possible.The coherently scattered light can then be represented as a coherent state with mean photon number α 2 = n coh T, where, without loss of generality, we assume α to be real.In the time domain, in a frame rotating with the laser frequency, the two-photon component of the state can be written as where a † t (a t ) is the creation (annihilation) operator for a photon at time t, and α (2) = n coh /2 is the amplitude of the two-photon wavefunction.In the low saturation limit, the incoherently scattered light consists of energy-time entangled photon pairs [24][25][26][27] and the state in the time domain can be written as Here, the integrals extend over time duration T and the temporal envelope ϕ (2) (τ) is given by The total scattered field is thus given by |ψ⟩ = |α⟩ + |ϕ⟩ (11)   and for its two-photon component one obtains Modelling the photon statistics.The spectral distributions of the two components differ from each other; the spectrum of the coherent state |α〉 is approximately given by a delta function at laser frequency ω L , whereas the spectrum of the incoherently scattered component, at low saturation, is approximately given by a pair of Lorentzian lines centred at frequencies ω L ± Δ.This frequency separation between the coherent and incoherent components allows us to separate them by using a narrowband notch filter at the laser frequency, which in our experiment is realised as a fibre-ring resonator.The complex-valued field transmission factor 36 of the resonator is given by Here, κ ext is the external coupling rate into the filter resonator, κ 0 is the intrinsic loss rate of the filter resonator, and ω res is its resonance frequency.As the filter bandwidth is much smaller than the separation of the frequency components, in the following we make the simplifying assumption that the filter acts only on the coherently scattered light and reduces its field amplitude by the κ ext -dependent transmission factor t coh (κ ext ) ≡ t F (ω L , κ ext ), as presented in the main text.The wavefunction of the light after filtering is then given by and the corresponding envelope of the two-photon wavefunction is given by In our experiment, we analyse the photon statistics of the filtered atomic fluorescence by measuring its second-order correlation function where 〈…〉 t indicates the time-averaged expectation value and τ = t 2 − t 1 is the temporal separation between the two photons.
is the unnormalized second-order correlation function, for which we obtain from equation ( 15) The mean photon flux, n, after the filter is given by Here we introduce the total single-photon collection and detection efficiency of our set-up, η.In the experiment, in addition to the measured photon rate, we also have to take into account the total background count rate, d, that originates from stray light and detector dark counts.Including this background in the experimentally expected second-order correlation function, we obtain Description in the time domain.It has been shown that the spectral filtering of coherent light can alter its statistical properties, such that the phase noise of, for example, a laser light field 37 , is changed into intensity fluctuations by a narrow linear filter, resulting in chaotic light in its output.Although such a classical reasoning cannot be applied to our experiment, the transformation of the stream of antibunched photons emitted by the atom into a stream of bunched photons after the notch filter can still be explained by evoking Hong-Ou-Mandel quantum interference 30 .To this end, we make use of the fact that the resonator has a slow temporal response to the light field with a characteristic timescale of (κ 0 + κ ext ) −1 that is much slower than the amplitude fluctuations in the atomic fluorescence, which are on a timescale of ~Δ−1 .Consequently, the resonator averages over these fluctuations and, in the steady state, the light stored in the resonator can approximately be described by a coherent state with amplitude To calculate the probability of a photon coincidence detection with time delay τ after the filter resonator, we have to consider the different ways in which such a detection event can occur.In this picture, there are two fields, the resonator field and the incident fluorescence, which are overlapped on the incoupling beamsplitter of the resonator.Two photons can either both originate from the original fluorescence light (amplitude, ψ (2) f ), both from the resonator (amplitude, ψ (2) r ) or one photon from each (amplitude, ψ (2) fr ).Note that the incoupling beamsplitter is highly asymmetric with a transmission of approximately one for the incident fluorescence light.The three probability amplitudes are then given by The total two-photon detection probability amplitude is the coherent sum of the three components, and for the two-photon field after the resonator we obtain This is the same expression as equation (15), which demonstrates the equivalence of the two pictures.
Implementing a photon-pair source.Using the filtering technique presented in this Article to remove the coherently scattered light lends itself to realising spectrally narrowband energy-time entangled photon-pair sources that are inherently compatible with optical transitions in quantum emitters.As it relies on the fundamental mechanism of resonance fluorescence, this technique can be utilized with any two-level emitter.The crucial figure of merit for such a photon-pair source is the achievable photon pair rate.In general, the rate of photon pairs, n p , after the optical filtering is given by where η is the photon collection efficiency of the optical set-up.Note that, in this inequality, n p ≈ η 2 n inc /2 only holds for S < 1.For close-to-unity collection efficiencies that can be realised using high-NA optics 33 or solid-state quantum emitters with integrated structures 34,35 , the maximum entangled photon pair rate that can be realised is bounded by γ/2 for large saturation powers.This is the fundamental limit that can be achieved for any source of photon pairs with linewidth γ, as otherwise the temporal separation of subsequent photon pairs is smaller than their extension in time, ~1/γ.Consequently, one can conclude that for high photon-collection efficiencies, the filtered resonance fluorescence yields an ideal source of photon pairs.Comparison with alternative theoretical models 38 shows that our theoretical description should be valid for saturation parameters of S ≲ 0.25, otherwise the incoherent component also contains higher-order photon-scattering processes.
For S = 0.25, one would obtain a pair source with photon-pair rates of n p = γ/50 ≈ 4 × 10 5 pairs per second with a bandwidth of γ = 2π × 3 MHz, which is substantially better than state-of-the art photon-pair sources with a comparable bandwidth 39 .
Effect of non-perfect filtering.Imperfect filtering of the coherent light, that is, |t coh | > 0, would lead to a background of coherently scattered photons.The ratio of the coherent to incoherent components after optical filtering is given by For a low background of coherently scattered photons, this signal-to-noise ratio should be much larger than one.For saturation parameters around S ≈ 0.1, this can be readily achieved with a residual filter transmission of |t coh | 2 ≈ 1%.

Experimental methods
Trapping and detecting single atoms.We prepared a cloud of 85 Rb atoms inside an ultrahigh-vacuum chamber using a MOT.The MOT lasers, with frequency ω L , are red-detuned with respect to the unperturbed D 2 transition of 85 Rb (transition frequency ω 0 ) by δ MOT = ω L − ω 0 = −2π × 16.3 MHz.The MOT cloud contains several million 85 Rb atoms and is used as a reservoir of cold atoms for loading the optical dipole trap.To load a single atom, the MOT cloud is positioned at the focus region of an aspheric lens with NA of 0.55 (AS-AHL12-10, Asphericon), located inside the vacuum and with a focal length of f = 10 mm and a working distance of w d = 7.6 mm.The optical dipole trap is generated by focusing the incident trap laser beam (wavelength, λ trap = 784.65 nm) to a waist radius of w 0 = 1.8 ± 0.2 μm inside the MOT cloud.For a laser power of P trap = 2.5 mW, we obtain an optical trapping potential with a depth of U trap /k B = 1.66 mK, corresponding to trap frequencies of ω r = 2π × 96 kHz and ω z = 2π × 17 kHz in the radial and axial directions, respectively.The photon scattering rate of the trapping light field amounts to 0.86 kHz.The presence of an atom in the trap is registered by an increase in the photon detection rate from 120 s −1 to 1,050 s −1 for the case without spectral filtering (Extended Data Fig. 1).Here, we collect the fluorescence light with the same lens that we use for focusing the trap beam.Due to the microscopic trap volume, our trap operates in the collisional blockade regime 40 such that, at most, a single atom is present inside the trapping volume at any time (Extended Data Fig. 1).The collected fluorescence light is separated from the trapping light with a bandpass dichroic mirror centred at 780 nm that features https://doi.org/10.1038/s41566-023-01260-7linewidth of 3 nm (LL01-785-25, Semrock).We then couple the light into a single-mode fibre, which also acts as a spatial filter.The fibre-guided fluorescence is sent via the fibre-ring resonator filter to a fibre-based Hanbury Brown and Twiss set-up consisting of two fibre-coupled single-photon counting modules behind a 50/50 fibre-optic coupler.The two single-photon counting modules are connected to a field-programmable gate array-based time-tagger that records the arrival times of each detected photon.Spectral filtering.Our fibre-based optical notch filter consists of a variable ratio coupler (F-CPL-830-N-FA, Newport) that allows setting of the coupling rate, κ ext , into the resonator.A polarization controller is included in the fibre-ring section, such that the two polarization eigenmodes of the resonator are degenerate and the resonator acts as a polarization-independent filter.This is important for consistent spectral suppression of the collected atomic fluorescence.A section of the fibre ring is glued to a piezo-electric stack, which allows us to strain-tune the resonance frequency, ω res .The resonator is passively stabilized by placing its whole set-up inside a thermally and acoustically insulating box, isolating the set-up from external environmental fluctuations.To compensate for residual, slow drifts of the resonance frequency, every 5 s we launch a stabilization light field at frequency ω L into the resonator set-up, scan the resonator length, and set it such that the point of minimum transmission is reached.The resonator has a geometrical length of l = 2.25 ± 0.05 m, yielding a free spectral range of ν FSR = 89.0 ± 0.5 MHz.Extended Data Fig. 2 shows the onresonance transmission of the resonator for different coupling rates κ ext .From a fit to the data, we obtain an intrinsic resonator loss rate of κ 0 = 2π × 1.08 ± 0.02 MHz, which results in an unloaded resonator finesse of ℱ = πν FSR /κ 0 = 41.2 ± 0.8.

Analysis of measured correlation functions.
To fit the theory to the measured correlation functions in Fig. 3, we perform the following procedure.
First, the measured correlation functions exhibit a weak bunching effect on a microsecond timescale, which originates from the diffusive motion of the trapped atoms through the complex intensity and polarization pattern of the MOT laser field 41 (Extended Data Fig. 3).To account for this, we fit the function 1 + Ae |τ|/t b to the data for large time delays |τ| ≫ 1/2γ, from which we obtain typical values of its amplitude A = 1.38 and decay time t b = 3.48 μs for the measurement without filtering (κ ext = 0).For each g (2) measurement, our model is then adjusted such that the value 1 + A serves as the new baseline.Second, the finite temperature of the atom in the dipole trap gives rise to a temperaturedependent distribution of atomic positions in the trap and, consequently, of the a.c.Stark shifts.To account for this, we assume a thermal position distribution and fit the second-order correlation function of the unfiltered case (κ ext = 0) with the atom temperature and the maximum a.c.Stark shift as free fit parameters.The latter takes into account that the atom experiences a substantial tensor light shift, which results in a Zeeman state-dependent transition frequency.From the fit, we obtain a temperature of 144 ± 47 μK and a mean detuning of the atom from the MOT lasers of Δ/2π = −57.9± 3.7 MHz.These values are used to fit the data obtained for all other settings of the filter resonator.
Finally, for the measurements that incorporate the filter resonator, an additional experimental uncertainty occurs due to drifts of the narrow resonator absorption line with respect to the MOT laser frequency.To account for these residual resonator drifts in our model, we assume that the filter exhibits a Gaussian probability distribution of resonance frequencies with a width of σ centred around an average resonator-laser detuning δ = ω res − ω L .Using this model, we fit all measured correlation functions in Fig. 3 with k ext ≠ 0. From the results we obtain an average laser-resonator detuning of δ = −0.26± 0.58 MHz with a width of σ = 1.92 ± 0.36 MHz.
Detection efficiency.From the atomic temperature obtained from the previous section, we calculate the distribution of saturation parameters S from which we obtain the expected scattering rates of the incoherent and coherent components.Together with the laserfilter detuning, we can then calculate the expected photon detection rate.We fit this model with the total photon detection efficiency, η, as the only free parameter to the data shown in the inset of Fig. 4. The fit yields η = 0.136 ± 0.002%, which agrees with the value expected from our collection efficiency of ~1.3%, together with the propagation losses through the filter set-up and the limited detector efficiency of ~0.55.
Correlation function at zero time delay versus filter setting.To evaluate the theory prediction of the zero-time-delay value of the second-order correlation as a function of κ ext shown in Fig. 4, we use the approximate expression for g (2) (τ) given in equation (19).Here, we use the temperature-induced distribution of S, the residual laserresonator detunings and the photon detection efficiency η, which we obtained from previous measurements, as well as the mean detector background rates as parameters.
Photon-pair rate.To determine the number of detected photon pairs, we start with the histogram of measured coincidences and subtract the base level found for |τ| ≫ 1/2γ, which stems from accidental coincidences.By summing up all remaining coincidences within a time interval of ±100 ns around τ = 0 and dividing by the effective measurement time, we obtain the rate of detected photon pairs.Given our photon-pair detection efficiency η 2 /2, we then determine the total rate of photon pairs scattered by the atom, 2n (2) meas./η 2 .For the case of strong suppression of the coherent component, that is, t F ≈ 0, these detected pairs all originate from the incoherently scattered light, so the rate of scattered photon pairs should be given by n inc /2 (equation ( 6)). when one atom is trapped.The absence of higher count rates indicates the sub-poissonian occupation statistics due to the collisional blockade effect present in microscopic optical traps 40 .

κFig. 3 |Fig. 4 |
Fig. 3 | Measured second-order correlation functions g (2) (τ) for different filter settings κ ext .Each dataset has the same vertical scale.The solid lines are fits of our theoretical model to the data points (open circles), with the temperature of the trapped atom and a residual resonator-laser detuning as fit parameters (Methods).For the unfiltered dataset (κ ext = 0), the purple horizontal line indicates the noise floor due to background counts.

Extended Data Fig. 1 |
Histogram of the detected photon count rate originating from the trap volume.Without spectral filtering, we observe an increase in the most likely count rate from a background of 120 s −1 , to 1050 s −1