## Abstract

High-gain parametric down-conversion (PDC) is a source of bright squeezed vacuum, which is a macroscopic nonclassical state of light and a promising candidate for quantum information applications. Here we study its properties, such as the intensity spectral width and the spectral width of pairwise correlations. In agreement with the theory, we observe an increase in the spectral width by 27% compared with the low-gain PDC. Frequency cross- and auto-correlations are registered by measuring the reduction of noise in the difference of PDC intensities at various pairs of wavelengths. The noise reduction plots also demonstrate super-bunching typical for collinear frequency-degenerate PDC.

© 2012 Optical Society of America

Bright squeezed vacuum (BSV), the state generated via high-gain parametric down-conversion (PDC) and having a large number of photons per mode, is currently in the focus of attention in quantum optics and quantum information. The reason is that BSV is considered as a state that can manifest macroscopic entanglement [1–4]. At the same time, because BSV can have huge numbers of photons per mode, it can be easily involved into nonlinear optical interactions [5] and interactions with material quantum objects such as atoms, molecules, quantum dots etc. In particular, its application to quantum memory is discussed [6]. Recently, BSV has been used for quantum imaging [7] and absolute calibration of photodetectors [8, 9]. Because of these applications, it is very important to know the spectral properties of BSV emitted via high-gain PDC.

It should be noted here that the spectra of PDC are essentially two-dimensional (see, for instance, Ref. [10]): the shape of the frequency spectrum depends on the selected angle of emission, and vice versa. In this connection, it is worth mentioning an interesting concept of two-dimensional coherence of PDC radiation [11]. Therefore, a general consideration should involve both frequency and angular spectra. However, in this work we only consider frequency spectra and all measurements are performed with an aperture selecting only collinear emission.

At high gain, PDC does not any more generate two-photon light, as is the case of low gain. The radiation of high-gain PDC contains not only photon pairs, but four-photon states, six-photon states, and so on, and the mean photon number per mode is large. The mean photon number per mode is in one-to-one correspondence with *brightness*; this is why the state can be called *bright squeezed vacuum*. Under these conditions, correlations in SV cannot be efficiently studied by measuring the number of coincidences. Instead, one should measure the noise (variance) of the intensity difference for the two selected beams [12–14].

The spectral and angular properties of high-gain PDC (known earlier as parametric super-fluorescence [15]) have been well studied theoretically [15–17]. In particular, it has been known for a long time [15] that the spectral width of both signal and idler beams should increase with the gain. Physically, this is clear from the fact that at high gain, there is exponential amplification of down-converted light along the crystal. Therefore, effectively it is not the whole crystal that contributes into the signal and idler intensity but only its exit part, where the signal/idler intensities (photon numbers) are most high. For the same reasons, the spectral width of correlations, usually determined by the inverse pump pulse width/coherence time, should also increase at high gain. Similar behavior should be observed for the angular spectra: both the total angular width and the far-field zone speckle size (correlation radius) should increase at high gain. Indirect experimental evidence for this was observed in Ref. [12]. The speckle size has been directly measured as a function of the parametric gain in Ref. [18]. However, no measurements on the frequency spectra of high-gain PDC have been performed, although there was an experiment on the study of spectral correlations for a similar effect, four-wave mixing [19].

In this work, we study the spectral properties of BSV generated via high-gain PDC from a picosecond-pulse pump. Due to a rather short coherence time of the pump, we are able to resolve the spectral width of correlations between signal and idler radiation. By scanning both the signal and idler frequencies separately, we can see the two-dimensional distribution of correlations, similar to the ‘joint two-photon intensity’ observed for two-photon light [20–22]. We also measure the spectral width versus the parametric gain and observe the spectral broadening in accordance with [15]. Finally, we observe some interesting effects caused by super-bunching in high-gain PDC.

The experimental setup is shown in Fig. 1. BSV was obtained in a 2 mm long BBO crystal cut for type-I collinear frequency-degenerate phasematching. As the pump, we used the 3rd harmonic of a Nd:YAG laser with the wavelength *λ _{p}* = 354.7 nm, pulse duration 18 ps, repetition rate 1 kHz and the energy per pulse up to 0.1 mJ. The coherence time of the pump was less than the pulse duration and amounted to 5 ps. To increase the parametric gain, the pump was focused into the crystal by means of a lens (L1) with the focal length 1 m. This resulted in the maximal peak intensity 10 GW/cm

^{2}. The gain value could be varied by changing the pump power, which was done by rotating a half-wave plate (HWP) in front of a Glan prism (GP). After the crystal, the pump radiation was reflected by a dichroic mirror (DM), which transmitted nearly all the down-converted light (97%), and its power was measured with a power-meter. The PDC radiation was filtered in the transverse wavevector (angle of emission) by means of an aperture (A) placed in the far-field zone. The emission angle selected this way was 0.17° (inside the crystal). All radiation passing through the aperture was focused, using a lens (L2), on the input slit of the HORIBA Jobin Yvon Micro HR monochromator equipped with a SYG-NATURE PDA CCD array. The spectral measurements were performed with the best possible resolution, which was 0.2 nm.

In all our measurements, the spectral properties of PDC were studied depending on the parametric gain values. In its turn, the parametric gain was found by fitting the dependence of the PDC output intensity *I* on the input pump power *P* (inset in Fig. 1) by the function

*I*

_{0},

*κ*were the fitting parameters [23]. The gain was measured for collinear frequency-degenerate regime. It is important that for such a measurement, the PDC radiation should be very well filtered both in the angle and in frequency, so that less than a single mode is registered. After the fitting, each pump power is put into correspondence to a certain value of the parametric gain $G\equiv \kappa \sqrt{P}$.

The dependence of the PDC spectral width (FWHM) on the parametric gain was measured for the range of gain values 3.9 < *G* < 6.5. At lower gain values, the PDC signal was too low to make spectral measurements while higher gain values could not be achieved with the available pump powers. The measured dependence is shown in Fig. 2. One can see that in the chosen range of gain values, the spectral width changes by 12%. Compared to the spectral width of low-gain PDC, the increase is 27%. Note that the effect of the gain on the spectral width is even stronger in the case of type-II or frequency-nondegenerate phase matching [15]. The insets show the shapes of the spectra recorded at the lowest and highest gain values. The shapes are well fitted with the theoretical dependencies calculated using the approach of Ref. [17]. The same approach was used in calculating the theoretical dependence of the FWHM on the parametric gain (dashed line). One can observe perfect agreement between the theory and experiment. The only fitting parameter used in the theoretical dependence was the orientation of the crystal, which turned out to differ by 0.0025° from the one for exact collinear degenerate phase matching. This led to about 2 nm shift for the whole dependence, which was difficult to notice in the experiment.

Here, it is worth mentioning that the ‘raw’ spectra obtained at the output of the spectrometer were always asymmetric, the long-wavelength part having smaller intensity. A typical distribution is shown in Fig. 3 by red circles. It is interesting that similar distributions have been reported in many papers on degenerate (broadband) PDC - see, for instance, Ref. [22].

By analyzing the detection scheme (standard for such measurements), which includes the selection of a certain small solid angle and a certain wavelength interval, we come to the conclusion that the spectra should be multiplied by a factor scaling as the fourth degree of the wavelength, *λ*^{4}. The reason is that the total number of photons is the product of the number of photons per mode (wavelength-independent in the vicinity of the degenerate wavelength) and the number of modes collected. At the same time, the number of transverse and longitudinal modes within a given solid angle interval ΔΩ and a given wavelength interval Δ*λ* scales as 1/*λ*^{4}. This can be understood from the following considerations. The number of modes is given by the ratio of the detection volume and the coherence volume [24],

*V*=

_{coh}*S*, where

_{coh}ct_{coh}*S*is the transverse coherence area,

_{coh}*c*is the speed of light, and

*t*is the coherence time. The detection volume does not depend on the wavelength, while the coherence area and time are wavelength-dependent,

_{coh}*k*

_{⊥}is the interval of transverse wavevector variation and Δ

*ω*is the frequency bandwidth corresponding to Δ

*λ*. As a result, the number of collected modes scales as

*λ*

^{−4}, and this factor should be taken into account in the spectra. Figure 3 shows the shape of the spectrum before (red circles) and after (black squares) the correction. The correction also took into account the fact that in our CCD camera, the wavelength interval selected by the one pixel slightly depended on the wavelength. This lead to a 10% variation of the wavelength resolution along the CCD camera.

It should be pointed out that the *λ*^{4} factor must be taken into account not only for the spectral measurements on PDC, but for all spectral measurements involving wavelength and angular selection. It is the flatness of type-I PDC spectrum that makes the asymmetry very noticeable without this correction.

In order to study the spectral correlations in high-gain PDC, we used the technique of noise reduction measurement [12–14], in which the observable under study is the variance of the difference of intensities for two beams. Correlation between two beams leads to a decrease in the noise of the difference intensity. We have studied the difference-intensity noise between beams at various pairs of wavelengths in the PDC spectrum. For different gain values, 10000 – 20000 frames were made, each one collecting the intensity from 500 – 1000 pulses, and then processed. First, one wavelength was fixed (by fixing the pixel of the CCD camera) and the other one was scanned. The variance was calculated for the difference between the signals from the chosen pixel and all remaining ones. The results are shown in Fig. 4. In all spectra, the fixed wavelength (‘signal’) was chosen to be *λ _{s}* = 702.4 nm, while the other one (‘idler’),

*λ*, was scanned within the range from 700 nm to 720 nm.

_{i}One can see that each spectral distribution of the variance contains two ‘dips’ and one peak. The left-hand ‘dip’ is caused by the classical correlations existing within a single mode of PDC radiation. The dip goes down to zero because the variance is identically zero for the case where the reading of a pixel is subtracted from itself. This point is an artifact and should not be considered. Another artifact is the oscillating noise seen on some (but not all) spectra in the vicinity of the left-hand ’dip’. The right-hand dip is caused by the nonclassical correlation between light beams scattered into parametrically conjugated modes, with the wavelengths related as ${\lambda}_{s}^{-1}+{\lambda}_{i}^{-1}={\lambda}_{p}^{-1}$.

The peak in Fig. 4 occurs at exactly the degenerate wavelength 709.3 nm. It appears because at the degenerate wavelength, PDC radiation is single-mode squeezed vacuum, sometimes referred to as ‘light with even photon numbers’, and this kind of radiation has peculiar statistical properties. Namely, single-mode squeezed vacuum manifests super-bunching even at strong parametric gain [25–27]. Its second-order normalized Glauber’s correlation function is *g*^{(2)} = 3 + 1/*N*, where *N* is the mean photon number per mode, and in the high-gain regime *g*^{(2)} = 3. Therefore, it has increased fluctuations, which causes the peak in Fig. 4.

The dips and the peak in Fig. 4 can be explained in more detail by recalling that the variance of the photon-number difference between the beams labeled *s* and *i* is [28]

*m*independent modes are registered: the variance in the left-hand side, as well as the mean values in the right-hand side, will acquire a factor of

*m*, and the factors containing the correlation functions in the right-hand side will change as

*g*

^{(2)}– 1 → (

*g*

^{(2)}– 1)/

*m*[23].

From the second term we see that the variance of the difference signal should have a peak if *λ _{i}* corresponds to the degeneracy. Indeed, if 〈

*N*〉 = 〈

_{s}*N*〉 ≡

_{i}*N*, then outside of the peak, where ${g}_{\mathit{ii}}^{(2)}={g}_{\mathit{ss}}^{(2)}=2$, the variance of the difference photon number is 2

*N*

^{2}while in the peak, ${g}_{\mathit{ii}}^{(2)}=3$, and the variance is 3

*N*

^{2}. In the ‘autocorrelation dip’,

*λ*≈

_{i}*λ*, ${g}_{\mathit{ii}}^{(2)}={g}_{\mathit{ss}}^{(2)}=2$, and Var(

_{s}*N*–

_{s}*N*) ≈ 2

_{i}*N*, i.e., the variance of the difference photon number should be suppressed down to the shot-noise level. In the ‘cross-correlation dip’, where

*λ*and

_{s}*λ*are conjugated, ${g}_{\mathit{ij}}^{(2)}=2+1/N$, and in theory, Var(

_{i}*N*–

_{s}*N*) should be zero. However, this complete suppression of noise can be observed only provided that three conditions are satisfied [14]: (1) the quantum efficiency is unity; (2) a large number of transverse and longitudinal modes are collected and (3) the signal and idler modes are properly matched. In our experiment, the quantum efficiency is not more than 20%, only a few transverse modes are selected, and no proper mode matching is performed.

_{i}According to Eq. 4, the super-bunching peak in Fig. 4 should be 50% higher than the background value, independently of the gain. In all dependencies, we actually see only a 30 ± 5% increase of the variance, which is caused by the insufficient frequency selection. As to the cross-correlation dip, ideally it should go down to zero but, due to the imperfect mode matching, the depth of the dip should be a certain fraction of the background. In all plots of Fig. 4, the value at the ‘dip’ is 40 ± 5% of the background value. At smaller gain values, the measurement is not accurate enough to reveal the correlations.

Another observation one can make from Fig. 4 is that the peak width is considerably less than the width of the ‘dips’. Accurate fitting of all of them by Gaussian functions yields, in the case of the middle dependence in Fig. 4, the values of 0.26 ± 0.06 nm, 0.45 ± 0.06 nm and 0.52 ± 0.06 nm for the peak, the ‘autocorrelation dip’ and the ‘cross-correlation dip’, respectively. The widths of the dips coincide with the widths of the second-order correlation functions and indicate the spectra of frequency correlations. Theoretically, they should be given by the pump spectral width, which is Δ*ω _{p}* = 1.1 THz. This would result in the low-gain correlation width of Δ

*λ*= 0.15 nm at wavelengths around 709 nm, but at high gain should be larger due to the nonlinear dependence of the PDC intensity on the pump intensity.

As to the peak, its width, theoretically, should be narrower than the correlation width by a factor of 2. The observed widths of 0.26 ± 0.06 nm for the peak and 0.52 ± 0.06 nm for the dip are in agreement with this figure.

More complete information about the frequency correlations in the spectrum of high-gain PDC can be obtained by measuring a two-dimensional distribution of some parameter characterizing the correlations, as a function of both wavelengths *λ _{s}*,

*λ*. As such a parameter, one can choose the difference-signal variance given by Eq. (4) but it is more convenient to measure the covariance of the signals at two various wavelengths,

_{i}The covariance is equal, up to a constant factor, to the third term in (4), since
$({g}_{\mathit{si}}^{(2)}-1)\u3008{N}_{s}\u3009\u3008{N}_{i}\u3009=\text{Cov}({N}_{s},{N}_{i})$. It gives a full account of the pairwise correlations, which now manifest themselves as peaks (Fig. 5). The ‘autocorrelation’ peaks occur on the straight line *λ _{s}* =

*λ*and correspond to the classical correlations within a single PDC mode. The ‘cross-correlation’ peaks correspond to the signal-idler correlations and are on the line given by ${\lambda}_{s}^{-1}+{\lambda}_{i}^{-1}={\lambda}_{p}^{-1}$. The distribution of these peaks is similar to the ’joint intensity’ observed for two-photon light [20–22]: its cross-section (conditional distribution) shows the spectral range of correlations while its projections on the axes

_{i}*λ*,

_{s}*λ*(marginal distributions) show the PDC spectrum. The increased value of the correlation function for the collinear frequency-degenerate PDC also reveals itself in the vicinity of the point

_{i}*λ*=

_{s}*λ*= 709.3 nm. While the shape of the distribution in Fig. 5 is the same as for the joint intensity (see, for instance, Ref [22]), the widths are different. Namely, in the low-gain regime the full spectral width is 70.2 nm and the correlation width is 0.15 nm, while at gain 6.5 the corresponding numbers are 89.5 nm and 0.52 nm, respectively.

_{i}As a possible extension of our work, it would be interesting to combine the spectral expansion with the angular expansion (‘crossed dispersion’ method of obtaining two-dimensional spectra, see, for instance, Ref. [10]). Then, by putting a 2D CCD camera into the focal plane of a lens, one could study both angular and frequency correlations in high-gain PDC spectra. Such measurements could be done, for instance, with a setup of Ref. [11]. In particular, the dependence of the spectral widths on the parametric gain could be studied for non-collinear PDC. However, we would expect a similar dependence as observed here.

In conclusion, we have performed an experiment aimed at the study of the spectral properties of high-gain PDC. In agreement with the theory, we observed the broadening of the spectrum with the increase in the parametric gain. We have shown that the dependence of the coherence volume on the wavelength causes the asymmetry of the measured PDC spectra, which can be removed by multiplying the spectra by a factor of *λ*^{4}. By measuring the variance of the intensity difference for two wavelengths, we saw the correlations in the PDC spectrum: both auto-correlations, due to the thermal statistics of each mode of the PDC radiation, and cross-correlations, due to the pairwise statistics of signal and idler modes. The technique of noise reduction measurement allows us to resolve the spectral width of the intensity correlations. The measured variance distributions also show the super-bunching of PDC radiation at exactly degenerate wavelength. Finally, we have measured the distribution of the covariance for the PDC intensities at two wavelengths, both scanned within the range 650 – 760 nm. The obtained two-dimensional distributions show both the intensity correlations within a single mode of PDC radiation and the cross-correlation between the intensities at the signal and idler wavelengths.

## Acknowledgments

We would like to thank Gerd Leuchs for illuminating discussions. This work has been supported in part by the Russian Foundation for Basic Research, grants 10-02-00202 and 11-02-01074. T.Sh.I. acknowledges funding from Alexander von Humboldt Foundation and K.Yu.S., from the ‘Dynasty’ foundation.

## References and links

**1. **W. J. Munro and M. D. Reid, “Violation of Bell’s inequality by macroscopic states generated via parametric down-conversion,” Phys. Rev. A **47**, 4412–4421 (1993). [CrossRef] [PubMed]

**2. **Ch. Simon and D. Bouwmeester, “Theory of an entanglement laser,” Phys. Rev. Lett. **91**, 053601 (2003). [CrossRef] [PubMed]

**3. **A. Gatti, R. Zambrini, M. San Miguel, and L. A. Lugiato, “Multiphoton multimode polarization entanglement in parametric down-conversion,” Phys. Rev. A **68**, 053807 (2003). [CrossRef]

**4. **C. Vitelli, N. Spagnolo, L. Toffoli, F. Sciarrino, and F. De Martini, “Quantum-to-classical transition via fuzzy measurements on high-gain spontaneous parametric down-conversion,” Phys. Rev. A **81**, 032123 (2010). [CrossRef]

**5. **O. Jedrkiewicz, J.-L. Blanchet, A. Gatti, E. Brambilla, and P. Di Trapani, “High visibility pump reconstruction via ultra broadband sum frequency mixing of intense phase-conjugated twin beams,” Opt. Express **19**, 12903–12912 (2011). [CrossRef] [PubMed]

**6. **L. V. Gerasimov, I. M. Sokolov, D. V. Kupriyanov, and M. D. Havey, “Quantum hologram of macroscopically entangled light via the mechanism of diffuse light storage,” arXiv:1111.6669v1 [quant-ph] 29 Nov 2011.

**7. **G. Brida, M. Genovese, and I. Ruo Berchera, “Experimental realization of sub-shot-noise quantum imaging,” Nat. Photonics **4**, 227–230 (2010). [CrossRef]

**8. **G. Brida, I. P. Degiovanni, M. Genovese, M. L. Rastello, and I. Ruo-Berchera, “Detection of multimode spatial correlation in PDC and application to the absolute calibration of a CCD camera,” Opt. Express **18**, 20572–20584 (2010). [CrossRef] [PubMed]

**9. **I. N. Agafonov, M. V. Chekhova, T. S. Iskhakov, A. N. Penin, G. O. Rytikov, and O. A. Shcherbina, “Absolute calibration of photodetectors: photocurrent multiplication versus photocurrent subtraction,” Opt. Lett. **36**, 1329–1331 (2011). [CrossRef] [PubMed]

**10. **A. V. Burlakov, M. V. Chekhova, D. N. Klyshko, S. P. Kulik, A. N. Penin, Y. H. Shih, and D. V. Strekalov, “Interference effects in spontaneous two-photon parametric scattering from two macroscopic regions,” Phys. Rev. A **56**, 3214–3225 (1997). [CrossRef]

**11. **O. Jedrkiewicz, A. Picozzi, M. Clerici, D. Faccio, and P. Di Trapani, “Emergence of x-shaped spatiotemporal coherence in optical waves,” Phys. Rev. Lett. **97**, 243903 (2006). [CrossRef]

**12. **O. Jedrkiewicz, Y.-K Jiang, E. Brambilla, A. Gatti, M. Bache, L. A. Lugiato, and P. Di Trapani, “Detection of sub-shot-noise spatial correlation in high-gain parametric down conversion,” Phys. Rev. Lett. **93**, 243601 (2004). [CrossRef]

**13. **M. Bondani, A. Allevi, G. Zambra, M. G. A. Paris, and A. Andreoni, “Sub-shot-noise photon-number correlation in a mesoscopic twin beam of light,” Phys. Rev. A **76**, 013833 (2007). [CrossRef]

**14. **T. Sh. Iskhakov, E. D. Lopaeva, A. N. Penin, G. O. Rytikov, and M. V. Chekhova “Two methods for detecting nonclassical correlations in parametric scattering of light,” JETP Lett. **88**, 660–664 (2008). [CrossRef]

**15. **D. N. Klyshko*Photons and Nonlinear Optics* (Gordon and Breach, 1988).

**16. **E. Brambilla, A. Gatti, M. Bache, and L. A. Lugiato, “Simultaneous near-field and far-field spatial quantum correlations in the high-gain regime of parametric downconversion,” Phys. Rev. A **69**, 023802 (2004). [CrossRef]

**17. **B. Dayan, “Theory of two-photon interactions with broadband down-converted light and entangled photons,” Phys. Rev. A **76**, 043813 (2007). [CrossRef]

**18. **G. Brida, M. Genovese, A. Meda, I. Ruo-Berchera, and E. Predazzi, “Tailoring PDC speckle structure,” J. Mod. Opt. **56**, 201–211 (2009). [CrossRef]

**19. **S. Spälter, N. Korolkova, F. König, A. Sizmann, and G. Leuchs, “Observation of multimode quantum correlations in fiber optical solitons,” Phys. Rev. Lett. **81**, 786–789 (1998). [CrossRef]

**20. **Y.-H. Kim and W. P. Grice, “Measurement of the spectral properties of the two-photon state generated via type II spontaneous parametric downconversion,” Opt. Lett. **30**, 908–910 (2005). [CrossRef] [PubMed]

**21. **H. S. Poh, C. Y. Lum, I. Marcikic, A. Lamas-Linares, and C. Kurtsiefer, “Joint spectrum mapping of polarization entanglement in spontaneous parametric down-conversion,” Phys. Rev. A **75**, 043816 (2007). [CrossRef]

**22. **S.-Y. Baek and Y.-H. Kim, “Spectral properties of entangled photon pairs generated via frequency-degenerate type-I spontaneous parametric down-conversion,” Phys. Rev. A **77**, 043807 (2008). [CrossRef]

**23. **O. A. Ivanova, T. Sh. Iskhakov, A. N. Penin, and M. V. Chekhova, “Multiphoton correlations in parametric down-conversion and their measurement in the pulsed regime,” Quantum Electron. **36**, 951–956 (2006). [CrossRef]

**24. **D. N. Klyshko, *Physical Foundations of Quantum Electronics* (World Scientific, 2011).

**25. **R. Loudon and P. L. Knight, “Squeezed states,” JMO **34**, 709–759 (1987). [CrossRef]

**26. **J. Janszky and Y. Yushin, “Many-photon processes with the participation of squeezed light,” Phys. Rev. A **36**, 1288–1292 (1987). [CrossRef] [PubMed]

**27. **T. Sh. Iskhakov, K. Yu. Spasibko, A. Perez, M. V. Chekhova, and G. Leuchs, “Super-bunched bright squeezed vacuum state,” (submitted) (2012).

**28. **T. Iskhakov, A. Allevi, D. A. Kalashnikov, V. G. Sala, M. Takeuchi, M. Bondani, and M. Chekhova, “Intensity correlations of thermal light,” Eur. Phys. J. Special Top. **199**, 127–138 (2011). [CrossRef]