# High density terahertz frequency comb produced by coherent synchrotron radiation

S Tammaro, Olivier Pirali, P Roy, Jean-Francois Lampin, Guillaume Ducournau, A. Cuisset, Francis Hindle, Gaël Mouret

► To cite this version:

S Tammaro, Olivier Pirali, P Roy, Jean-Francois Lampin, Guillaume Ducournau, et al. High density terahertz frequency comb produced by coherent synchrotron radiation. Nature Communications, 2015, 6 (1), pp.1-6 / Article n° 7733. 10.1038/ncomms8733. hal-04830388

HAL Id: hal-04830388

https://hal.science/hal-04830388v1

Submitted on 11 Dec 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
ARTICLE

Received 22 Mar 2015 | Accepted 5 Jun 2015 | Published 20 Jul 2015

DOI: 10.1038/ncomms8733

OPEN

# High density terahertz frequency comb produced by coherent synchrotron radiation

S. Tammaro<fnref n="1" />,², O. Pirali<fnref n="1" />,³, P. Roy<fnref n="1" />, J.-F. Lampin⁴, G. Ducournau⁴, A. Cuisset², F. Hindle² & G. Mouret²

Frequency combs have enabled significant progress in frequency metrology and high-resolution spectroscopy extending the achievable resolution while increasing the signal-to-noise ratio. In its coherent mode, synchrotron radiation is accepted to provide an intense terahertz continuum covering a wide spectral range from about 0.1 to 1 THz. Using a dedicated heterodyne receiver, we reveal the purely discrete nature of this emission. A phase relationship between the light pulses leads to a powerful frequency comb spanning over one decade in frequency. The comb has a mode spacing of 846 kHz, a linewidth of about 200 Hz, a fractional precision of about $2 \times 10^{-10}$ and no frequency offset. The unprecedented potential of the comb for high-resolution spectroscopy is demonstrated by the accurate determination of pure rotation transitions of acetonitrile.

<footnote n="1">¹ AILES Beamline, Synchrotron SOLEIL, l'Orme des Merisiers, Saint-Aubin, Gif-sur-Yvette 91192, France. ² Laboratoire de Physico-Chimie de l'Atmosphere, Université du Littoral Côte d'Opale, EA4493 CNRS, 189A Avenue Maurice Schumann, Dunkerque 59140, France. ³ Institut des Sciences Moléculaires d'Orsay, UMR8214 CNRS—Université Paris-Sud, Bâtiment 210, Orsay 91405, France. ⁴ Institut d'Electronique de Microélectronique et de Nanotechnologie, UMR8520 CNRS—Université de Lille 1, Avenue Poincaré-Cité Scientifique CS 60069, Villeneuve d'Ascq 59652, France. Correspondence and requests for materials should be addressed to G.M. (email: mouret@univ-littoral.fr).</footnote>
Frequency combs (FC) have radically changed the landscape of optical instrumentation containing a vast number of narrow optical modes at a series of well-defined frequencies, which may be used as a highly accurate ruler over its spectral range. Initially developed in the visible and near-infrared spectral regions¹, the use of FC has been expanded to the mid-infrared², extreme ultra-violet³, X-ray⁴ regions and significant effort is currently being dedicated to the generation of FC at THz frequencies. One solution based on converting a stabilized optical FC using a photoconductive terahertz emitter remains hampered by the low available THz power⁵,⁶. Another approach is based on active mode locked THz quantum-cascade-lasers providing intense FC over a relatively limited spectral range⁷⁻⁹.

As an alternative to laser based experiments and electronic generation, broadband THz radiation emitted from short electron bunches is now used for high-resolution spectroscopy and spectro-microscopy applications. In the so-called low-α mode, intense THz power is obtained by reducing the electron bunch length to produce phase coherent synchrotron radiation (CSR)¹⁰. Following the first demonstration, such radiation has been developed in several synchrotron facilities over the world. It has provided unprecedented power in the region 0.1–1 THz, up to 4 orders of magnitude brighter than standard incoherent synchrotron radiation¹⁰,¹¹. Super-Radiance effects obtained from bunch-to-bunch interference¹² further enhance specific frequencies (1/T<sub>bunch</sub> where T<sub>bunch</sub> is the time interval between successive bunches) of the THz CSR.

In this work, using a dedicated high-performance sub-THz heterodyne receiver (see Methods), we reveal that THz CSR is an intense, stable (in both frequency and amplitude), dense zero offset FC covering a wide spectral range 0.1–1 THz with no continuum background.

## Results

**Spectral characterization of CSR.** In this work, the CSR spectral structure has been probed by a Bruker Fourier Transform Infrared (FTIR) instrument with an ultimate resolution of 30 MHz, as well as by the proposed heterodyne receiver shown in Fig. 1. The receiver, capable of operation in a large spectral region with a bandwidth of about 10 GHz, probes the synchrotron radiation with unprecedented spectral resolution, high signal-to-noise ratio and excellent frequency accuracy. Using both spectrometers, we are able to detect the fine spectral features of the CSR. In particular, we identified two main ultra-stable FC sequences in both regions studied (200 ± 10 GHz and 400 ± 10 GHz). The FC sequences are illustrated in Fig. 2 with a progressively finer frequency scale. The upper panel of this figure shows the CSR spectrum as recorded by the FTIR spectrometer at the highest possible spectral resolution (maximum optical path difference of 8.82 m, corresponding to a resolution of 30 MHz). The full CSR range convoluted by FTIR response is shown in Fig. 2a. A detailed view ~600 GHz reveals a 352 MHz FC superimposed on the broad THz emission, Fig. 2b, as previously observed in other facilities¹². The lower panel shows the heterodyne analysis of the CSR. Figure 2c shows 3 GHz of heterodyne spectrum at 200 GHz, in which the 352 MHz FC originating from the bunch-to-bunch repetition is also clearly resolved. This FC is the signature of a degree of coherence among the electron bunches¹³ and is designated as super-radiance emission. Expanding the frequency scale further reveals a second FC composed of sharp teeth regularly spaced by 846 kHz, Fig. 2d. This comb is indeed related to the very stable revolution period of the electron bunches in the storage ring (1.18 µs). It produces a spectrally dense THz FC with >10⁶ components covering the THz range from 0.1 to 1 THz. The observed structure as shown in

<figure><img src="image_2.png" /></figure>

**Figure 1 | CSR heterodyne detection schematic.** Short bunches of relativistic electrons circulating in the SOLEIL storage ring emit CSR in the 0.1–1 THz range. The bunch-to-bunch frequency 1/T<sub>bunch</sub> is 352 MHz, while the storage ring revolution frequency is 846 kHz (see Methods). The THz CSR pulses are mixed with the monochromatic radiation (f<sub>LO</sub>) from a local oscillator (microwave synthesizer and multiplier) creating an electronic signal at the intermediate frequency (IF): f<sub>IF</sub> = |f<sub>LO</sub> - f<sub>CSR</sub>| within the mixer bandwidth. The resulting signal is amplified by the IF amplifier A and analysed using a microwave spectrum analyser. The mixing produces a frequency downconversion of f<sub>CSR</sub> into the microwave range, and an aliasing of the spectrum due to the superposition of the upper (f<sub>CSR</sub> > f<sub>LO</sub>) and the lower (f<sub>CSR</sub> < f<sub>LO</sub>) frequency bands.
<figure><img src="image_3.png" /><figcaption>**Figure 2** | CSR FC structure recorded using the FTIR interferometer and heterodyne receiver. All amplitude scales are linear unless stated otherwise. Heterodyne frequency scales are established as a function of the intermediate frequency. (a) Full range of the CSR measured by the FTIR interferometer with a resolution of 30 MHz; the low-frequency response is limited by the optical elements of the instrument. (b) Zoom of 3 GHz revealing the 352 MHz FC corresponding to $1/T_{\text{bunch}}$. (c) Heterodyne spectrum obtained by mixing CSR with $f_{\text{LO}} = 200$ GHz. The FC at 352 MHz is also easily distinguished. (d) The second FC composed of sharp teeth regularly spaced by 846 kHz is revealed by an expanded frequency scale. This second comb is produced by the very stable revolution period of the electron bunches in the storage ring (1.18 µs). (e) A zoom of a given dual comb structure caused by the aliasing of the upper and lower frequency bands. (e), inset: the excellent resolving power of the frequency scale allows the detailed examination of a single tooth, highlighting the ultrafine splitting due to low-frequency instabilities of the electron bunches.</figcaption></figure>

Fig. 2d is composed of the lower and upper frequency band contributions as expected for heterodyne detection technique. Aliasing of the FC results in the observation of a doubled FC structure whose frequency separation is dependent on the local oscillator (LO) frequency. Shifting the LO frequency allowed the two components to be distinguished, here a separation of 100 kHz was selected, Fig. 2e. The detailed examination of any individual tooth, Fig. 2e inset, indicates an even narrower structure: typically, two components separated by 1.3 kHz, each with a full width at half maximum (FWHM) systematically <500 Hz. The substructure is interpreted as a consequence of the kHz range amplitude modulation caused by the low-frequency instabilities of the electron bunches¹⁴.

**Power measurement.** The heterodyne analysis allowed the power of an individual CSR line ~400 GHz to be estimated at ~77 dBm (that is, 20 pW) (see Methods). The CSR spectral power density is 3 orders of magnitude greater than can be achieved from a classical dipolar antenna used in a context of THz time domain spectroscopy⁵. The dynamic range of the amplitude of the FC teeth is systematically in excess of 35 dB with the strongest lines exhibiting over 70 dB.

**Offset determination.** In general, FC modes are described as multiples of a reference frequency ($f_R$) shifted by an offset ($f_0$): $f_{\text{FC}} = n \times f_R + f_0$, where n is an integer. The offset frequency was assessed by measuring the revolution frequency with a counter (CNT 90, Pendulum). Over a period of 20 s, the revolution frequency is determined with a typical uncertainty of 30 µHz. The comb mode centre frequencies were determined by the spectral analysis of the intermediate frequency (IF) signal. In excess of 50 modes were examined, each being found to be an exact multiple of the revolution frequency with a confidence interval of 1 Hz. Unlike optical FC, the CSR FC is therefore considered to be offset free.

**Pure rotational spectroscopy of acetonitrile.** We recorded the absorption spectrum of acetonitrile (CH₃CN) as the K series of this symmetric top molecule is particularly well suited to test the frequency accuracy and the dynamic range of the CSR FC. A 1.3-m-long absorption gas cell was inserted into the CSR path. Figure 3 shows a small section of the acetonitrile absorption spectrum obtained in 65 s using the CSR FC with the LO at 202 GHz. The FC gives rise to a discrete spectrum built from the harmonics of the repetition rate of 846 kHz (violet vertical lines).
<figure><img src="image_4.png" /><figcaption>**Figure 3** | Normalized absorbance of acetonitrile. Normalized absorbance versus frequency showing the K structure of the pure rotational R(10) transition of acetonitrile recorded with the heterodyne receiver, *f*<sub>LO</sub> = 202.0 GHz (violet vertical lines). The gas pressure was set to 50 µbar, the measurement and background IF spectra are processed identically to isolate the upper frequency band FC only (*f*<sub>CSR</sub> > *f*<sub>LO</sub>). The violet solid line is the result of the fit of the 5 individual lines using Voigt profiles (see Methods), while the red curve corresponds to a simulation at the maximum possible resolution (30 MHz) of commercially available Fourier Transform (FT) spectrometers. The instrument configurations are schematically represented in the insets. Inset (a): heterodyne measurement, the synchrotron radiation transmitted through a cell filled with gas is incident on a mixer simultaneously with the radiation from a local oscillator. (b), inset: absorbance measurement based on FT method, the interferometre-modulated radiation passes through a gas-containing cell before arriving at the detector.</figcaption></figure>

The complete *K* structure of the R(10) transition of CH₃CN is clearly resolved (violet profile, as explained in the associated content). This is in strong contrast to the simulated spectrum corresponding to the 30 MHz resolution of the IFS125 interferometre (red profile), only displaying a broad unresolved feature. In particular, the THz FC allows the doublet *K*=0 and *K*=1 separated by 3.9 MHz (1.3 × 10⁻⁴ cm⁻¹) to be resolved. The fit of the five individual absorption lines with Voigt profiles provides their central frequencies, which are in excellent agreement with reference data (r.m.s<80 kHz) (see Methods). The accuracy obtained is comparable with frequency multiplication techniques currently employed in this frequency range.

Discussion

The measurements described in the previous section shows clearly the discrete nature of the CSR. Moreover, the relatively high power together with the frequency stability of this source allows for room temperature heterodyne fast detection. Such results demonstrate that the phase relation between the emissions from bunches exists and is preserved over a long time to form a dense FC. In this sense, in the low alpha configuration, the operational mode of the synchrotron can be likened to ultrashort laser pulses where the first repetition rate is determined by the length of the ring resonator, and a second being defined by the bunch-to-bunch spacing. It is noteworthy that by using the same

detection scheme, no discrete spectral components and no FC have been observed in standard synchrotron mode proving that for longer bunches of electrons there is no phase coherence between subsequent pulses, preventing the formation of a dense FC.

As shown by the quality of the THz acetonitrile absorption spectrum, the CSR FC is particularly suited to high-resolution gas phase THz spectroscopy as the mode spacing of 846 kHz is the same order of magnitude of the Doppler-limited linewidths of a typical rotational transition of a light polar compound at room temperature. This is not the case for many laser-generated THz combs with mode spacing typically > 100 MHz, requiring complex interleaving strategies to be adopted to obtain Doppler-limited rotational spectra⁶,¹⁵. Indeed the CSR FC composed of powerful, fine modes over an extended spectral range (0.1–1 THz) should outperform FTIR instruments in terms of resolution and analysis times. In comparison to electronic sources that offer an excellent resolution, the advantage of CSR FC is the spectral coverage of one decade, whereas a frequency multiplier may offer a tunability of ~20%, numerous multiplier combinations are required to cover an octave of frequency. The excellent resolving power, frequency metrology and high dynamic range of the CSR FC are demonstrated by the high-resolution spectroscopy measurements undertaken using the heterodyne receiver initially constructed for the spectral analysis.
To harness the full potential of this FC, a solution must be sought to simultaneously access its entire frequency band (that is, over more than one decade of frequency). Mixing the CSR FC with an optical comb in a suitable photoconductor should yield the required information via a dual comb spectroscopy approach unlocking the advantages of the CSR¹⁵.

The CSR FC characteristics revealed here clearly prompt the development of receivers for ultra-high resolution over broader frequency ranges. This is an enabling approach for new high-resolution THz time-resolved experiments. Among highly desirable new applications are ultrafast monitoring of reaction kinetics, ultra-high resolution atomic and molecular spectroscopy¹⁶, and spatio-temporal dynamics of relativistic electron bunches.

## Methods

**Coherent synchrotron radiation.** The AILES beamline collecting the edge emission¹⁷ from the 2.75 GeV, third generation Synchrotron Radiation source SOLEIL, has an extended programme of research with synchrotron radiation in the THz spectral range. In particular, a special attention is paid into increasing the flux while maintaining source stability through optimization of CSR¹¹. The SOLEIL magnet lattice allows for the storage ring low alpha operation appropriate for the generation of the CSR in the THz spectral range. For this operation, relativistic electrons entering the dipole magnetic field are known to emit coherently at wavelengths comparable with or larger than the electron bunch length. In presence of a storage ring vacuum chamber, the coherent spontaneous emission can be extracted for wavelengths smaller than the cut-off wavelength value of the vacuum chamber, acting as a waveguide. The principle consists in providing conditions for stable electron dynamics with the smallest possible bunch length. In storage rings, this typically requires a low value of the momentum compaction factor and a smaller electron bunch length. The complete sequences of intense THz FC combs were observed in the low α coherent mode characterized by a 4.8 ps r.m.s bunch length and corresponding to a compaction factor $α = α_0/25$, with $α_0 = 1.7 \times 10^{-5}$ being the nominal value¹⁴. The ring was injected with 208 electron bunches distributed on one half of the SOLEIL ring and by one isolated bunch in the centre of the other half. In this operation mode the total current was about 16 mA with periodic refilling in top-up mode. Two principal periods characterize the temporal radiation distribution: (i) the 208 bunches produce a 2.8 ns periodic pulsed radiation train and (ii) the revolution period in the ring is about 1.18 µs.

**Heterodyne detection.** To investigate the spectral structure of the THz CSR, we constructed a high sensitivity coherent electronic receiver¹⁸. Heterodyne detection is very convenient for high-resolution spectral analysis with a high dynamic range. The operation is based on a frequency translation of the analysed spectra by electrical mixing between a known reference frequency (LO) and the signal to be analysed. The mixing element is a Schottky-based sub-harmonic-mixer (SHM), fed by the LO at mW level (continuous wave mode), this power being required to operate the Schottky device in a non-linear regime. Two different SHM mixers WR5.1 (140-220 GHz) and WR2.2 (325-500 GHz) are used for operation at 200 GHz and 400 GHz, respectively. The LO is generated by an electronic frequency multiplication chain driven by a Rohde and Schwarz SMF100A synthesizer. This provides equivalent frequency multiplication factors of × 12 and × 24, respectively. The IF signal is amplified with two base-band amplifiers, providing a total gain of 45 dB over a bandwidth of 500 MHz. The downconverted CSR signal is analysed by a 3.6 GHz Agilent MXA Signal Analyzer N9020A.

**Power measurements.** From IF signal spectral analysis, the individual CSR lines located around 400 GHz display a typical signal strength of −40 dBm. Taking into account the 45 dB IF gain and a 8 dB conversion loss of SHM, the corresponding power of CSR lines should be around −77 dBm (20 pW). The integrated THz power of the CSR was measured to be 60.5 µW by a Gentec THz 5B-BL-DZ pyroelectric detector working at room temperature. The beam was filtered with a black polyethylene film to avoid any effect of visible or mid-infrared radiation. The ratio of the integrated power to the power of one line is then equal to about $3 \times 10^6$. This gives a rough estimation of the number of lines assuming that all lines of the FC have the same intensity. This value is compatible with the fact that in a 1 THz range there should be about $1.2 \times 10^6$ lines separated by 846 kHz. The spectral power density is hence calculated to be $23.6 \text{ nW GHz}^{-1}$.

**Spectroscopic analysis on acetonitrile.** The IF spectra are processed to isolate the lower and upper frequency band contributions. Once all the peaks in the spectra have been identified, the FC spacing of 846 kHz is used to distinguish each FC component. In Fig. 3, the normalized absorbance of CH₃CN measured with the THz FC was simulated with Voigt profiles. The Doppler linewidth (Gaussian contributions) of the five components has been fixed to a calculated value at room

<table><caption>Table 1 | fitted frequencies compared with the frequencies listed in the JPL database¹⁹.</caption><thead><tr><th>K</th><th>f<sub>fitted</sub><br>(MHz)</th><th>f<sub>calc</sub><br>(MHz)</th><th>f<sub>fitted</sub> - f<sub>calc</sub><br>(MHz)</th><th>100 × (f<sub>fitted</sub> - f<sub>calc</sub>)/<br>f<sub>calc</sub></th></tr></thead><tbody><tr><td>0</td><td>202,355.39</td><td>202,355.51</td><td>-0.12</td><td>5.93E-05</td></tr><tr><td>1</td><td>202,351.73</td><td>202,351.61</td><td>0.12</td><td>5.93E-05</td></tr><tr><td>2</td><td>202,339.88</td><td>202,339.92</td><td>-0.04</td><td>1.97E-05</td></tr><tr><td>3</td><td>202,320.42</td><td>202,320.44</td><td>-0.02</td><td>9.89E-06</td></tr><tr><td>4</td><td>202,293.19</td><td>202,293.18</td><td>0.01</td><td>4.94E-06</td></tr></tbody></table>

temperature of 0.4 MHz (FWHM). The centre frequencies and the collisional linewidths (Lorentzian contributions) were fitted for $K=0-5$ components of the R(10) transitions. In Table 1, the fitted frequencies $f_{fitted}$ are compared with the frequencies $f_{calc}$ listed in the JPL database¹⁹ and calculated from the latest combined fit reported by Müller *et al.*²⁰. The r.m.s determined from the ($f_{fitted} - f_{calc}$) values is evaluated to be 79 kHz. The collisional linewidths were fitted to 3.4 MHz (FWHM) for the five components (no K dependence of the R(10) self-broadening was observed). At a pressure of 50 µbar a collisional linewidth of ~4.5 MHz (FWHM) may be estimated for the R(10) transitions from CH₃CN self-broadening measurements²¹,²².

## References

1. Udem, T. H., Holzwarth, R. & Hänsch, T. W. Optical frequency metrology. *Nature* **416**, 233–237 (2002).

2. Schliesser, A., Picqué, N. & Hänsch, T. W. Mid-infrared frequency combs. *Nat. Photonics* **6**, 440–449 (2012).

3. Zinkstok, R. T. h., Witte, S., Ubachs, W., Hogervorst, W. & Eikema, K. S. E. Frequency comb laser spectroscopy in the vacuum-ultraviolet region. *Phys. Rev. A* **73**, 061801 (2006).

4. Cavaletto, S. M. *et al.* Broadband high-resolution X-ray frequency combs. *Nat. Photonics* **8**, 520–523 (2014).

5. Tani, M., Matsuura, S., Sakai, K. & Nakashima, S. I. Emission characteristics of photoconductive antennas based on low-temperature-grown GaAs and semi-insulating GaAs. *Appl. Opt.* **36**, 7853–7859 (1997).

6. Finneran, I. A. *et al.* Decade-spanning high-precision terahertz frequency comb. *Phys. Rev. Lett.* **114**, 163902 (2015).

7. Barbieri, S. *et al.* Coherent sampling of active mode-locked terahertz quantum cascade lasers and frequency synthesis. *Nat. Photonics* **5**, 306–313 (2011).

8. Maysonnave, J. *et al.* Mode-locking of a terahertz laser by direct phase synchronization. *Opt. Exp.* **20**, 20855–20862 (2012).

9. Burghoff, D. *et al.* Terahertz laser frequency combs. *Nat. Photonics* **8**, 462–467 (2014).

10. Abo-Bakr, M. *et al.* Brilliant, coherent far-infrared (THz) synchrotron radiation. *Phys. Rev. Lett.* **90**, 094801 (2003).

11. Barros, J. *et al.* Coherent synchrotron radiation for broadband terahertz spectroscopy. *Rev. Sci. Instrum.* **84**, 033102 (2013).

12. Billingham, B. E. *et al.* Observation of superradiant synchrotron radiation in the terahertz region. *Phys. Rev. Spec. Top. Accel. Beams* **16**, 060702 (2013).

13. Shibata, Y. *et al.* Observation of interference between coherent synchrotron radiation from periodic bunches. *Phys. Rev. A* **44**, R3445 (1991).

14. Evain, C. *et al.* Spatio-temporal dynamics of relativistic electron bunches during the micro-bunching instability in storage rings. *Europhys. Lett.* **98**, 40006 (2012).

15. Hsieh, Y.-D. *et al.* Spectrally interleaved, comb-mode-resolved spectroscopy using swept dual terahertz combs. *Sci. Rep.* **4**, 3816 (2014).

16. Bartalini, S. *et al.* Frequency-comb-assisted terahertz quantum cascade laser spectroscopy. *Phys. Rev. X* **4**, 021006 (2014).

17. Roy, P. *et al.* Spectral distribution of infrared synchrotron radiation by an insertion device and its edges: a comparison between experimental and simulated spectra. *Phys. Rev. Lett.* **84**, 483 (2000).

18. Maas, S. A. Microwave Mixers (Artech House, 1986).

19. Pickett, H. M. *et al.* Submillimeter, millimeter, and microwave spectral line catalog. *J. Quant. Spectrosc. Rad. Transfer* **60**, 883–890 (1998).

20. Müller, H. S. P., Drouin, B. J. & Pearson, J. C. Rotational spectra of isotopic species of methyl cyanide, CH₃CN, in their ground vibrational states up to terahertz frequencies. *Astron. Astrophys.* **506**, 1487–1499 (2009).

21. Rinsland, C. P. *et al.* Multispectrum analysis of the v₄ band of CH₃CN: Positions, intensities, self- and N₂-broadening, and pressure-induced shifts. *J. Quant. Spectrosc. Rad. Transfer* **109**, 974–994 (2008).

22. Buffa, G. *et al.* Far-infrared self-broadening in methylcyanide: absorber-perturber resonance. *Phys. Rev. A* **45**, 6443–6450 (1992).
## Acknowledgements

This work was partly supported by the Université du Littoral Côte d'Opale in the framework of the 'Bonus Qualité Recherche' and by the PALM Labex (project FC THz, theme Emergence), by the synchrotron SOLEIL under the proposal 20131346, as well as by the Agence Nationale pour la Recherche (project ANR-12-NANO-0014). We would like to acknowledge the support of Jean-Blaise Brubach (SOLEIL), Laurent Manceron (SOLEIL) and Marc Fourmentin (LPCA). We are indebted to the SOLEIL accelerators department (in particular M.A. Tordeux, C. Evain and M.E. Couprie) and to the scientific director P. Morin for their continuous support.

## Author contributions

O.P. and G.M. wrote the synchrotron beamtime proposal and supervised the project. Experiments were performed by S.T., O.P., J.-F.L., G.D., F.H. and G.M. All the authors contributed equally to the analysis of the results and to the preparation of the manuscript.

## Additional information

**Competing financial interests:** The authors declare no competing financial interests.

**Reprints and permission information:** is available online at http://npg.nature.com/reprintsandpermissions/

**How to cite this article:** Tammaro, S. et al. High density THz frequency comb produced by coherent synchrotron radiation. *Nat. Commun.* 6:7733 doi: 10.1038/ncomms8733 (2015).

<figure><img src="image_5.png" /></figure>

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/