Field and intensity correlations: the Siegert relation from stars to quantum emitters

The Siegert relation relates field and intensity temporal correlations. After a historical review of the Siegert relation and the Hanbury Brown and Twiss effect, we discuss the validity of this relation in two different domains. We first show that this relation can be used in astrophysics to determine the fundamental parameters of stars, and that it is especially important for the observation with stellar emission lines. Second, we check the validity of this relation for moving quantum scatterers illuminated by a strong driving field.


Introduction
Light can be described using different tools, and in particular through the ones linked to its wave behavior, such as its coherence properties. The knowledge of these properties provides information about the light source itself, such as its angular intensity profile if one measures the spatial coherence, but also on the underlying light matter interaction processes when, for example, one measures the temporal coherence of light emitted or scattered by a medium. Temporal coherence properties are often characterized through the light spectrum S(ω), which corresponds to the light intensity distribution as a function of wavelength or frequency. More formally, for stationary processes, the spectrum is linked to the temporal field correlation function g (1) (τ ) through the Wiener-Khintchine theorem [1,2]: The temporal field correlation function g (1) (τ ), also called the first-order correlation function, is defined as: where . corresponds to the averaging over time t, and where the intensity of the field is given by I(t) = E (t)E(t).
To have a full description of the temporal coherence properties, one has to know the correlation functions g (n) (τ ) at all orders n [3]: g (n) (t 1 , t 2 , ..., t n , t n+1 , ..., t 2n ) However, there are cases where a full knowledge of correlation functions is not necessary. This is, for example, the case for chaotic light, such as the one emitted by a large number of independent scatterers, or more generally for light with with an electric field that can be considered as a random variable with a Gaussian (or normal) probability distribution. In this case, there exists a relation between the correlation functions at all orders. In particular, for spatially coherent polarized light one can relate the first orders of the correlation functions by the Siegert relation [4] g (2) (τ ) = 1 + |g (1) where g (2) (τ ) is the second-order temporal correlation function, or the temporal intensity correlation function, g (2) (τ ) = I(t)I(t + τ ) Physically, the g (1) (τ ) function quantifies the degree of mutual coherence between time-delayed electric fields: for chaotic light, it is equal to 1 at τ = 0 and decreases to zero at large time-delays over a characteristic time scale given by the coherence time of the field, inversely proportional to the spectrum width. From the Siegert relation, one can deduce that the g (2) (τ ) function decreases from 2 to 1, and that the intensity has a coherence time which is half the one of the electric field. The excess of intensity correlation at short delay (g (2) (τ = 0) > 1) is referred to the "bunching" of photons or "the Hanbury Brown and Twiss (HBT) effect", for historical reasons described in the next section.
But first, let us give a physical picture to understand this bunching effect in the simplest case of a spatially coherent chaotic source. Consider radiation with a finite optical spectrum of linewidth Δω. Since the source is chaotic, there is no phase relationship between the different spectral components. This is the case, for instance, if the light is generated from many independent emitters each with different velocities. In this configuration, the spectral phase φ(ω) can be considered as random and uniformly distributed between 0 and 2π. Let us now consider two frequency components from the optical spectrum. They induce a beat note at a frequency given by the difference between their optical frequencies. Since φ(ω) is a random variable, all the possible beat notes coming from all possible pairs sum up with random phases. This mechanism is responsible for intensity fluctuations. Note that this is a fully classical noise, due to the wave nature of the field and to its non-monochromaticity. This noise adds up to the photon noise (shot noise), which has a quantum origin. If the linewidth of the spectrum is infinite, one would get white noise and thus an intensity correlation function equal to unity, whatever the delay τ . On the other hand, a finite linewidth means that there is no beating at a frequency much larger than Δω. This cut-off in the power spectrum of the noise corresponds to a finite coherence time τ c ∼ 1/Δω and thus to correlated intensity fluctuations on this typical time scale.
The Siegert relation has been used in different domains. One can cite, for example, dynamic light scattering [5,6], where the intensity correlation is analyzed to determine the size of small scatterers. In this paper, we focus on two specific domains: astronomy and light scattered by quantum particles. After a brief introduction on the history of intensity correlations, we present how the Siegert relation can be used to determine fundamental parameters of stars. This relation is in particular interesting when one observes stars with strong emission lines. We then turn to the light scattered by quantum particles, namely cold atoms. We show that when those atoms are illuminated by a strong driving field, the incoherent scattered light still satisfies the Fig. 1 Simplified outline of an intensity interferometer for radio waves (a) and optical frequencies (b) taken from Ref. [13], with M1 and M2 the two mirrors, C1 and C2 the photocatodes of the two photomultipliers and A1 and A2 the two arials Siegert relation. Whereas this relation is usually derived in the classical domain [7,8], we give in Sect. 4.1 the detailed derivation for quantum emitters.

A brief history of intensity correlations
The history of the Siegert relation is intimately linked to the controversy on equal-time intensity autocorrelations, also known as the Hanbury Brown and Twiss effect. The story starts during World War II, when radar technology drove a lot of research in the field of radio waves with, later, much repercussion on radio astronomy and optical sciences. The relation between electric field correlations and intensity correlations has been proposed in that context by A. J. F. Siegert [4] in a technical report. It was later named "Siegert relation", mainly in the field of mesoscopic physics [9,10].
The next important step has been achieved by Hanbury Brown and Twiss in the field of radio astronomy. In 1952, they proposed and demonstrated a novel type of radio interferometer. The intensities of radio waves collected by two telescopes were recorded and correlated (see Fig. 1), first for the Sun and then subsequently for two radio sources in Cygnus and Cassiopeia, without measuring the electromagnetic amplitude and phase information. The angular sizes of these sources were determined by measuring the spatial intensity correlations for several different configurations of telescope separations [11]. In their 1954 paper [12] they wrote: 'It is further shown that the correlator output, when suitably normalized, is equal to the square of the correlation coefficient measured by the Michelson interferometer '. This statement links the intensity correlation function with the field correlation function, establishing again the Siegert relation. Indeed, the original technical report by Siegert had remained largely unnoticed and the relation was independently rediscovered at that time.
It was then natural, for radio-astronomers, to extend this new concept to visible light. However, they faced a strong opposition from several physicists who preferred to think about light in terms of photons [14]. Indeed, a temporal (or spatial) intensity correlation measurement relies on the detection of at least two photons. The classical description given in the previous section only relies on interference. What was puzzling at that time is that physicists were convinced, following Dirac, that "interference between two different photons can never occur " [15].
HBT successfully tested their idea with a laboratory demonstration [13], and a few months later on the light from a star [16]. Their results were first disputed as other groups failed to reproduce the lab experiment, and as it was claimed that such results, if true, would call for a major revision of quantum mechanics [17]. Nevertheless, it was later shown that the other experiments were simply not sensitive enough [18]. These first experiments were performed with continuous wave detection, in which case the notion of photons is not needed and the classical explanation is perfectly appropriate. In the photon-counting regime, one can still assume that the quantization only occurs at the detection of an underlying continuous quantity, in which case the instantaneous value of the intensity I(t) gives the probability of detecting photons, and the classical picture is still valid.
However, if one insists on describing light in terms of photons, another physical description is needed. A first argument was given by Purcell in 1956 [19]: the bunching of photons is a consequence of the Bose-Einstein statistics to which they obey. This interpretation was further developed by HBT [20] and Kahn [21], and verified in an experiment done, again by HBT, in the photon-counting regime [22]. Finally, another interpretation in terms of two-photon interference was given by Fano a few years later [23].
The HBT experiment and its understanding in the framework of quantum mechanics can be considered as the birth of modern quantum optics (before the laser was invented!). In particular, and as acknowledged by Glauber in his Nobel prize speech [24], it triggered the development of the quantum theory of optical coherence [3,25,26], based in particular on correlation functions. In this context, the Siegert relation, which provides a relation between the first and second-order field correlations, is a particularly important tool to probe the quantum nature of light. Intensity correlations experiments have been used as a signature to distinguish laser light from classical fields [27]. This experiment illustrates that higher order photon statistics are a fundamental tool to identify non classical states of light, which could not be identified as such by field-field correlation functions or the optical spectrum. The use of intensity correlation functions have also allowed the pioneering experiments illustrating the violation of Bell inequalities and effects on single photon sources in the early 80s [28][29][30][31] which opened the path towards the area of modern quantum technologies.

The Siegert relation in astronomy
The use of intensity correlations, also called intensity interferometry in astronomy, was pioneered by HBT, who measured the angular diameters of 32 stars and their fundamental characteristics from the spatial intensity correlation function [32]. After this series of impressive observations using the Narrabri Stellar Intensity Interferometer in Australia in the 1960s-1970s, this technique has been abandoned, mainly due to its poor signal-to-noise ratio compared to emerging techniques pioneered by A. Labeyrie in amplitude interferometry [33]. However, thanks to progress in photonics components, efficient single photon counting detectors, fast electronics and digital correlators, there is currently a strong effort from different groups towards the revival of intensity interferometry with modern photonic technologies. One strategy is to utilize large diameter (> 10 m) imaging air Cherenkov telescopes for intensity interferometry [34] and has recently resulted in successful on-sky measurements [35,36]. Our team is following an alternative approach by using traditional astronomical optical telescopes with photon-counting avalanche photodiodes (APDs) feeding a fast timetagger, which then computes the temporal correlations in real time [37][38][39][40].
Note that astronomers are mainly interested in the spatial intensity correlation function. While in the temporal domain, the bunching is related to the spectrum width of the source via the Siegert relation and the Wiener-Khintchine theorem, in the spatial domain, the bunching (which can be seen as the typical size of a speckle grain) is similarly related to the angular width of the source via the van Cittert-Zernike theorem [41]. In the following, however, we will focus on temporal intensity correlation measurements and show how this can be used to determine fundamental parameters of stars.

Experimental setup
The light coming from stars is generally assumed to be chaotic, so that it satisfies the Siegert relation. In order to test this assumption, one needs to measure both g (2) (τ ) and g (1) (τ ). S(ω) can also be measured instead of g (1) (τ ), since they are related through the Wiener-Khintchine theorem. In our specific case, the g (2) (τ ) function is measured thanks to the simplified setup depicted in Fig. 2 and described in detail in Ref. [42]. Very briefly, the light is first collected by a 1 m telescope at the C2PU facility on the Plateau de Calern site of Observatoire de la Côte d'Azur (OCA). It is then injected in a coupling assembly to perform spectral filtering with a narrow-band interference filter, polarization filtering with a polarizing beam splitter (PBS), and then injection in a 100 μm core multimode fiber (MMF). Since the stars we are looking at are not resolved by the telescope, we still inject only one spatial mode into the fiber. Finally, the output of the fiber is connected to a fibered beamsplitter (FBS) whose outputs illumi- nate two single-photon avalanche photodiodes (APDs). The photons detected by the APDs are time-tagged by a time-to-digital convertor (TDC) that also computes the temporal intensity correlation.

Test in the laboratory
The Siegert relation can be first demonstrated in the lab using an artificial unresolved star. The source is generated by injecting light from a halogen lamp, with a broad continuous spectrum, into a single mode fiber. The spectrum of the light on which g (2) (τ ) is computed thus corresponds to the spectrum of the light after the filter. This transmitted spectrum is plotted in Fig. 3a. The bandwidth of the filter is Δλ = 1 nm with a central frequency λ 0 = 656.3 nm.
From the spectrum transmitted by the filter, we can calculate its Fourier transform g (1) (τ ) and then g (2) (τ ) using the Siegert relation. The result is plotted in the inset of Fig. 3b. The coherence time τ c can be defined as the area of the bunching peak [43]. For technical reasons that will be explained in the following, we focus on the area, equal to 1.19 ps for this lab test. This coherence time is much lower than the electronic time resolution, mainly limited by the APD jitter, which is of the order of 500 ps, but which depends on different parameters such as the count rate and the beam size on the APD. This broader time response reduces the height of the measured bunching peak by approximately the Fig. 3 a Spectrum of the light coming from an halogen lamp injected in a monomode fiber (artificial star) and transmitted by a narrow filter, showing a spectral bandwidth of about 1 nm. b Gray line: temporal g (2) function measured on the artificial star. Dashed line: Theoretical g (2) (τ ) function, calculated from the spectrum assuming the Siegert relation, convolved by the electronic time response. Inset: theoretical g (2) (τ ) function calculated from the transmitted spectrum, assuming the Siegert relation and an infinite electronic response bandwidth ratio of the coherence time to the resolving time. On the other hand, the bunching width increases and its area remains constant.
The temporal jitter function can be assumed to be Gaussian at first order. Its full-width at half-maximum (FWHM) is estimated to be equal to √ 2×670 ps for this specific measurement, the √ 2 factor being due to the fact that we have two APDs to perform the measurement. The theoretical g (2) function now corresponds to the theoretical bunching peak with infinite electronic response bandwidth, shown in the inset of Fig. 3b, convolved by the Gaussian time response. This convolved g (2) function is shown in Fig. 3b with the dashed line. The area is still equal to 1.19 ps.
Finally, we plot the measured temporal g (2) function. This corresponds to the gray line in Fig. 3b. One can see that it is perfectly superimposed to the convolved theoretical g (2) function, as expected from the Siegert relation. The area measured with a Gaussian fit applied on the measured data gives 1.18 ± 0.08 ps, in agreement with the expected 1.19 ps. With a much longer integration time and a precise experimental characterization of the filtered spectrum, we have recently obtained an excellent agreement to the 1% level [40].

Measurements on stars
The next step is to check the Siegert relation for on-sky measurements. The first measurements were performed in the continuum of α Boo (Arcturus), α CMi (Procyon) and β Gem (Pollux), with a filter centered at 780 nm [44]. Since we are in the continuum, the spectra of the stars can be considered as constant over the transmission of the filter, and the bunching peak is thus simply limited by the bandwidth filter. In other words, the measured coherence time is limited by the filter and does not give any information on the star itself.
More interestingly, we can also focus on spectral emission lines. O, B and A supergiant stars have massive winds witnessed by those strong emission lines (e.g., hydrogen Balmer series, helium or carbon lines among others in the visible), and they mainly originate in a circumstellar envelope that extends from a few to tens and even hundreds of stellar radii. A well-known example in the northern hemisphere is P Cygni, which has a strong Hα emission line at 656.3 nm. To check the validity of the Siegert relation, one needs, as done in the previous section with the artificial star, to compute the spectrum transmitted by the filter thanks to the filter transmission and the measured star spectrum, and to measure the g (2) function of the light collected by one telescope. Since the star is not resolved by the telescope, we can consider that only one spatial mode is collected.
These lines exhibit time variability at different scales, from days to years and beyond. It is thus important to measure the stellar spectrum in the same period as the g (2) measurement. Figure 4a presents the spectrum emitted by P Cygni observed on August 8 th 2020 [45]. One can see that the linewidth, of the order of 0.5 nm, is smaller than the 1 nm width of the filter, with an emission line 15 times higher than the continuum. As already done in the previous section, we calculate the convolved g (2) from this spectrum, with a FWHM of the time response equal to 500 ps. This result is superimposed in Fig. 4b to the g (2) measurements performed on-sky between 3 August 2020 and 9 August 2020 [42]. Again, a good overlap is observed, illustrating the validity of the Siegert relation. The area of the theoretical g (2) function is equal to 2.35 ps, also in agreement with the area extracted from a Gaussian fit of the measurement, equal to 2.3 ± 0.3 ps.
We thus show here that the Siegert relation is valid for the light emitted by stars, in the continuum but also for specific cases such as strong emission lines. The measurement of simultaneous spectra and spatial intensity correlation functions g (2) (r) with two telescopes, the correlation at zero baseline being given by g (2) (τ ) measured with one telescope, is a powerful tool to estimate the angular diameter of the star, as pioneered by HBT. In association with spectral measurement and radiative transfer modeling, the distance of the star can be inferred, as shown in Refs. [42,46].  [45] on August 8th 2020. Black line: spectrum of P Cygni transmitted by the filter. b Grey line: temporal g (2) function measured on P Cygni between 3 August 2020 and 9 August 2020. Dashed line: Theoretical g (2) (τ ) function, calculated from the spectrum assuming the Siegert relation and convolved by the electronic time response

Towards quantum astronomy
The fact that starlight fulfills the Siegert relation is obviously not a surprise. Nevertheless, one may find, in the future, some natural sources of light in space for which this would not be necessarily true. Indeed, it has been speculated that natural astrophysical "lasers" (emission lines amplified by stimulated emission) could exist [47].
In the microwave domain, space masers are very common [48,49], but as the wavelength decreases, obtaining a population inversion gets harder. Amplification by stimulated emission was indeed observed in the planetary atmospheres of Mars and Venus in the infrared (IR) (λ ∼ 10 μm) [50,51] and much later in stellar atmospheres in the far IR [52]. More recently, the discovery of astrophysical lasers in the near IR was claimed by Johansson and Letokhov, based on Fe ii [53] and O i [54], but this interpretation is disputed [55]. The most promising emission line seems to be the Fe ii line at 1.68μm, which should be present in the circumstellar envelope (the part of the star surrounding its core) of some spectral types of stars [55].
One could even imagine that scattering-induced feedback could enhanced the amplification, even reaching the oscillatory regime [56], a phenomenon that would be called a random laser today [57,58]. Measuring the g (2) (τ ) function of such exotic emission lines would thus open the way to the study of coherence and quantum effects in space and intensity correlation experiments performed in laboratory random lasing systems indicate this to be a convenient indicator [59,60].

The Siegert relation for quantum emitters
We now turn to the Siegert relation for light scattered by quantum scatterers. As will be discussed later, our scattering medium is made of cold Rb atoms, illuminated by a laser probe beam. Varying the intensity of this probe allows us to tune the saturation parameter s, and thus the ratio between elastic and inelastic scattering. For elastic scattering, the different moving scatterers each independently shift the frequency of the incident light by Doppler effect, so the scattered field has the properties of the incident one, yet broadened [61]. In this case, we can consider the scatterers as classical ones and they generate what is often called pseudothermal light [62]. To validate the Siegert relation, one needs a large number of independent scatterers and a process to randomize the phase between the fields emitted by the different scatterers. For elastic scattering, the phase is randomized thanks to nonzero temperature. The validity of the Siegert relation in this configuration is presented in Ref. [63].
When the saturation parameter is much larger than one, the scattered field is dominated by the inelastic one, which corresponds to the well-known Mollow triplet for the spectrum [64]. This effect is purely quantum, the fluorescence spectrum corresponding to the one of a two-level system driven by a strong and resonant incident field. The Mollow triplet is an effect for which strong correlations exist between the two sidebands and the carrier [65], as first observed by A. Aspect et al. in 1980 [66]. The existence of correlations in the photon radiation also leads to antibunching, usually observed for a single two-level emitter, but which can be observed in large systems with a phasematching experimental configuration [67]. Even if not obvious, the Siegert relation is still valid in our setup, as will be demonstrated both with its derivation and experimental results in the following.

Derivation in the quantum formalism
For the derivation of the Siegert relation, we consider a large number of uncorrelated scatterers. A quantum treatment of the correlation functions g (1) and g (2) requires the use of electric field operatorsÊ, carefully accounting for their ordering, also known as "normal ordering" [3]. This also requires a quantum treatment of the emitters, here modeled as two-level atoms. Furthermore, the derivation of the Siegert relation requires that the mean value, which corresponds to the expected value in the quantum case, of the single-atom electric field vanishes. For two-level atoms, the electric field operatorÊ + j of atom j is proportional to the atomic lowering operatorσ j , which is characterized by a fluctuation operator δσ j =σ j − σ j t→∞ , see Ref. [68]. By definition, the expectation value of the fluctuations is equal to zero (yet not its variance). The scattered electric field is given by the sum of the contribution of each scatterer δσ j , so the (non-normalized) correlation function for the inelastically scattered intensity can be written as: We consider the scatterers to be independent from each other, so the expectation value may be split for different indexes ( Â jBm = Â j B m for j = m). Since the expectation value of fluctuation operators is zero, most terms vanish. When one expands this equation, only the terms in which the operator acts on its complex conjugate remain: Since there are only terms of the formÂ †Â , the specific phase of the driving field on each atom is unimportant, and all scatterers can be treated equally (we neglect shadow and multiple scattering effects since the saturation parameter is very large and the optical thickness is low). Thus, we can write: For a large number of scatterers, the term scaling as N 2 dominates, and the other terms become negligible. Normalizing Eq. (8) by the squared intensity, see Eq. (3) where the first (n = 1) and second (n = 2) order coherence are computed using the operator electric field, we obtain the Siegert relation for the inelastically scattered light: We note that the derivation is very similar to the classical case, although the operator nature of the electric field requires a careful ordering of the operators involved. Finally, the zero average value of the electric field scattered by single atoms here stems from the  (1) (τ ) and g (2) (τ ) of the light scattered by quantum scatterers. A cold atomic cloud is illuminated by a laser beam circularly polarized thanks to λ/2 and λ/4 plates. The scattered light is collected by a polarization-maintaining (PM) single-mode fiber after polarization selection with a λ/2 and a polarizing beam splitter (PBS). The light is split with a fibered beam splitter (FBS) and its two outputs illuminate two avalanche photodiodes (APDs). Each photon arrival is time-tagged by a time-to-digital converter (TDC) and the correlations are computed by a computer. Finally, a local oscillator (LO), derived from the same laser that illuminates the scattering medium and frequency-shifted by an acoustooptical modulator (AOM), is injected in the second input of the FBS nature of the inelastic scattering, rather than macroscopic features such as collisions or temperature [7].

Experimental setup
In our experiment, the scattering medium corresponds to a cold atomic cloud. This medium is produced by loading a magneto-optical trap (MOT) from a vapor of 85 Rb, as described in Refs. [61,63,69,70], and then released for a given time of flight. The number of atoms N is typically of the order of 10 8 − 10 9 . The number of scatterers is thus large, fulfilling one of the conditions needed to mimic chaotic light and thus to validate the Siegert relation. The temperature T is of the order of 100 μK. The atomic cloud is finally characterized by its on-resonance optical thickness b 0 , measured by recording the transmission of a small probe beam going through the cloud as a function of the detuning on the |3 → |4 D 2 hyperfine transition. In this paper, this parameter is lower than one to avoid any collective effects.
The experimental setup is depicted in Fig. 5. The cloud is illuminated by a flat-top intensity profile provided by a fibered laser (EYLSA from Quantel laser) injected in a beam shaper from Asphericon (TopShape model). Its frequency is locked on the |3 → |4 hyperfine transition of the 85 Rb D 2 line. We use waveplates to get a circularly polarized light. The beam diameter, equal to 14.7 mm, is adjusted to be larger than the cloud radius in order to get a uniform intensity on the atoms, and thus a constant Rabi frequency. Its intensity can be changed to tune the saturation parameter between less than 0.01, for which scattering is mainly elastic, to more than 60 with mainly inelastic scattering. The duration of the probe is adjusted to limit the number of photons scattered per atoms to a few hundreds and thus limit heating effect. Furthermore, for all parameters this duration is much larger than the coherence time of the speckle grains induced by the interference of all scattered fields, which guarantees that a statistical average is realized.
The scattered light is collected at typically 90 • from the probe beam, by a polarization-maintaining (PM) single-mode fiber. This ensures the selection of a single spatial mode. The polarization is also selected before the fiber with a λ/2 plate and a polarization beam splitter, to maximize the amount of collected photons as well as to adjust the incident polarization parallel to the PM fiber axis. To measure g (1) (τ ), we superimpose the collected scattered light to a local oscillator (LO) with a FBS. This local oscillator is derived from the laser which delivers the probe beam, frequency shifted by ω BN with an acousto-optical modulator. Its polarization is also adjusted before the entrance of the fiber to correspond to the PM fiber axis. The polarizations of the collected scattered light and the LO are thus parallel. The output of the FBS illuminates two APDs which are connected to a TDC. This last device allows to time-tag the arrival of each photon from which the intensity correlation function is finally calculated.
The intensity I BN of the beat note between the LO and the collected scattered light is used to compute the intensity correlation function g BN (2) (τ ) in our setup. For a LO coming from a laser and if one assumes either that the temporal coherence of the LO is much larger than that of the scattered light or that the LO is derived from the laser illuminating the scattering medium, one obtains: where I BN , I sc and I LO correspond to the intensity of the beat note, the collected scattered light and the LO, respectively. The g sc (1) (τ ) function corresponds to the temporal electric field correlation function of the scattered light, and g sc (2) (τ ) to its temporal intensity correlation function. The detailed derivation of this equation can be found in Ref. [63]. As can be seen in Eq. (11), this setup allows measuring the electric field and the intensity correlation functions of the light under study at the same time. This is particularly suited to check the validity of the Siegert relation, avoiding in particular any drift effect as when the g (1) and g (2) functions are measured separately.

Results
An example of the temportal g BN (2) function obtain ed for s = 60 is plotted in the inset of Fig. 6. One observes a low frequency oscillation corresponding to the beating between the carrier and the sidebands of the Mollow triplet, as well as a fast oscillation corresponding to the beating between the scattered light and the LO. To extract the contribution of g sc (1) (τ ) and g sc (2) (τ ), we take the Fourier transform of g BN (2) (τ ). Indeed, while the Fourier transform of g sc (2) (τ ), i.e., g sc (2) (ω), is centered around the DC value, the Fourier transform of g sc (1) (τ ),g sc (1) (ω), is frequency shifted by ω BN . Figure 6 presents the corresponding signal in the Fourier space. As soon asg sc (2) (ω) does not overlap withg sc (1) (ω), the two quantities can be extracted separately.
The Siegert relation given by Eq. (4) in the temporal domain can be written in the Fourier space: g (2) (ω) = δ(0) +g (1) (ω) * g (1) (ω), with δ(0) the Dirac function, * the convolution and the complex conjugate. To check its validity on our experimental setup, we use the following procedure. We first correct the Fourier transform to take into account the limited temporal response of the APDs, which can be approximated by a first-order low pass filter with a bandwidth limited by the inverse of the APDs jitter. We then extractg sc (2) (ω) andg sc (1) (ω). We shift g sc (1) (ω) back to zero frequency and we calculate its self-convolution. We finally renormalize this last quantity in order to match the height ofg sc (2) (ω) to the one ofg sc (1) (ω) * g sc (1) (ω) 1 . Note that this normaliza- tion depends on the intensity imbalance between I sc and I LO , as can be seen in Eq. (11). For this experiment, this normalization has been done by measuring the mean value ofg sc (2) (ω) andg sc (1) (ω) * g sc (1) (ω) close to the zero frequency.
The two previous extracted quantities are plotted in Fig. 7a for s = 60, thus with scattered light mainly composed of inelastic scattering. The open circles corresponds tog sc (2) (ω), while the plain curve corresponds tog sc (1) (ω) * g sc (1) (ω). One can see a very good overlap, with five peaks visible due to the convolution of the Mollow triplet. We also performed the same measurement for s = 1, for which we have equal contribution from elastic and inelastic scattering in terms of power. Figure 7b presents the corresponding data, showing again a good overlap.
A good overlap betweeng sc (2) (ω) andg sc (1) (ω) * g sc (1) (ω) indicates that |g sc (1) (τ )| 2 and g sc (2) (τ ) have the same temporal shape, as expected from Eq. (4). However, to fully validate the Siegert relation, we also have to check the quantitative overlap, meaning the contrast of the bunching peak g sc (2) (0) − 1 should be equal to 1, the value of g sc (1) (0) being equal to 1 by definition. We thus take different measurements of g sc (2) (τ ) of the scattered light without the LO and for different values of s between 1 and 60. The bunching peak contrast is plotted in Fig. 8. For most of the experimental outputs of the FBS. This negative sign is removed during the normalization.

Conclusion
In this paper, we have checked the validity of the Siegert relation in two different domains: in astronomy, for light coming from stars, and in cold-atom physics, when resonant laser light is scattered by a large number of quantum two-level systems. In astronomy, we have shown that this relation can be used to extract astrophysical information, such as the coherence time due to emission lines. Combined to spatial interferometry, this can be used to determine the angular diameter of the source. This is especially interesting to perform such measurements on emission lines to characterize the extended atmosphere of the star, which can give access to other fundamental parameters such as its distance [42,46].
To validate the Siegert relation, one needs to satisfy different assumptions. The first one is that the scatterers or emitters should be independent and their number should be large. The second assumption is that the emitted or scattered phase should be random and uncorrelated. While the process of phase randomization is obvious in stars with thermal radiation, it is a bit more complex for light scattered by quantum scatterers. The process that randomizes the phase actually depends on the scattering regime. On one hand, for low saturation parameter, scattering is mainly elastic and the phase randomization is due to temperature. The coherence time is thus linked to the atom velocities as well as the optical thickness which roughly corresponds to the number of scattering events that occurs inside the medium [61]. On the other hand, when the saturation parameter is large, light is mainly inelastically scattered. The origin of the phase randomization is completely different since it comes from the finite lifetime of the two-level excited state. The coherence time is thus of the order of a few nanoseconds, much lower than the few hundreds of nanoseconds when due to temperature. We have shown in this paper that the Siegert relation is also valid in this configuration. As far as we know, this is the first time that the Siegert relation is checked with the simultaneous measurement of g (1) (τ ) and g (2) (τ ) for large s where the particles behave as quantum emitters. Note that the Siegert relation is valid for any saturation parameter and any laser detuning, since these two parameters will simply change the ratio between elastically and inelastically scattered light [71].
The future of intensity correlations in astronomy can address two complementary goals. The first one is to increase the angular resolution in the visible, and more particularly for blue wavelength for which stellar amplitude interferometry [33] is still hardly feasible. The second one is to look for deviation from the Siegert relation, with the ultimate goal of being able to have a direct signature of astrophysical lasing. On the quantum emitters side, one may search for deviation from the Siegert relation as a signature of correlations between the scatterers, like, e.g., in the specific experimental configurations of Ref. [72], where antibunching was observed with many-atoms trapped in the vicinity of a nanofiber. One could also look for signatures of random lasing in such samples, based on high optical thickness to induce multiple scattering and radiation trapping [73], and gain mechanisms such as Mollow gain [74][75][76], Raman gain [58,75,[77][78][79][80][81] or parametric gain with four-wave mixing [75,82,83]. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/.