Observation of grating diffraction radiation at the KEK LUCX facility.

The development of linac-based narrow-band THz sources with sub-picosecond, [Formula: see text]-level radiation pulses is in demand from the scientific community. Intrinsically monochromatic emitters such as coherent Smith-Purcell radiation sources appear as natural candidates. However, the lack of broad spectral tunability continues to stimulate active research in this field. We hereby present the first experimental investigation of coherent grating diffraction radiation (GDR), for which comparable radiation intensity with central frequency fine-tuning in a much wider spectral range has been confirmed. Additionally, the approach allows for bandwidth selection at the same central frequency. The experimental validation of performance included the basic spectral, spatial and polarization properties. The discussion of the comparison between GDR intensity and other coherent radiation sources is also presented. These results further strengthen the foundation for the design of a tabletop wide-range tunable quasi-monochromatic or multi-colour radiation source in the GHz-THz frequency range.

Introduction

Various applications of THz radiation demand a high–brilliance monochromatic source with tunable characteristics such as spectral range, pulse duration, polarization, and directivity[1-4]. Semi–conventional linac–based THz sources can provide sub-picosecond radiation pulses at approximately the hundred -level with a continuous spectrum up to THz[5-9]. However, in many applied investigations, narrow–band sources at this energy level are desired[10-14]. Evidently, to achieve this objective, the usage of a monochromator is required. However, approaches offering more output energy are related to the utilization of sources based on radiation mechanisms which are intrinsically monochromatic; for instance, coherent Smith-Purcell radiation (SPR)[15] or coherent Cherenkov radiation from dielectric lined waveguides[16]. Typically, spectral line fine-tuning of SPR-based sources is performed by output photon angular selection, i.e., by detection aperture positioning[17] or by the rotation of an additional mirror upon parallel passage of an electron beam near a grating[18]. Superior spectral tuning can only be accomplished by the replacement of the grating with another having a different profile or period[19]. In the refs. [20,21] initial proposals to use non–parallel orientation of a grating relative to a beam to change the frequency of the SPR spectral lines for a given outgoing photon angle were presented. In this case, the position of the collimator or aperture, which provided the radiation beam spectral acceptance, was fixed. Moreover, in ref. [20] the generalized dispersion relation for the inclined grating was obtained as:where is the radiation wavelength, is the grating period, is the diffraction order, is the electron velocity in light units, is the grating orientation angle, and is the observation angle. For a charge trajectory parallel to the grating () Eq. 1 coincides with the well-known SPR dispersion relation . It is important to note that Eq. 1 is also valid for larger grating orientation angles that are typically unreachable for SPR generation experiments due to the longitudinal grating dimensions and the small distance to the electron beam. However, implementation of large angle radiation emission may be considered in a “diffraction radiation”–like arrangement. This should lead to different radiation polarization maps, intensities and directivity profiles, although the approach offers wide–range monochromatic spectral tunability which still follows Eq. 1. Radiation generated in this case, by analogy with the grating transition radiation[22-24], also has significant spectral tunability but without the drawback of the destructive interaction of the electron beam with the grating material. This radiation is referred to as the “grating diffraction radiation” (GDR).

In this report, we demonstrate the first experimental observation of coherent GDR which includes the investigation of the basic spectral, spatial and polarization properties in addition to a discussion on GDR intensity in comparison with coherent SPR and coherent diffraction radiation (CDR). The result improves the prospects of designing a tabletop wide–range tunable quasi–monochromatic or multi–colour radiation source in the GHz–THz frequency range.

Results

Details of the experimental geometry are illustrated in Fig. 1. The Coulomb field of the relativistic charge with the effective radius ( is the Lorentz-factor and is the radiation wavelength) interacts with a tilted periodic structure on the length for and on the length if , assuming the axis of rotation in the center of the grating and , where is the grating length and is the horizontal impact-parameter. If the relation is satisfied, one can expect that the resulting radiation would become quasi-monochromatic and the full width at half maximum (FWHM) bandwidth can be estimated as[25]:

Figure 1

Schematics of the experimental layout (a) and grating geometry (b,c). Abbreviations: - fixed interferometer mirror, - movable mirror, - beam splitter, - parabolic mirror with a focal length of mm, - grating length, - interaction length, - vertical impact-parameter, - horizontal impact-parameter, - grating inclination angle, - observation angle, - grating period, ☉ - grating rotation axis in XZ plane.

The monochromaticity estimation from the above–given equation is valid only for a small capture angle of the detection system , and if the chosen grating parameters and geometry satisfy the condition , where is the number of grating periods. The factor comes from the Fourier transform of the periodical function (grating periods), shaped by the rectangular window representing the finite grating length[26].

In the present investigation, we measured the spectral–angular distribution and polarization properties of the radiation produced by the MeV, pC single electron bunch with mm rms length passing below the mm period, mm2, , echelette profile grating (Fig. 1b,c). To reveal the actual full bandwidth of the spectrometer system GHz, the broadband coherent transition radiation (CTR) spectrum from the flat surface on the rear of the grating plate was measured. The CDR was generated from the same surface, but when the target’s vertical position was set to allow for the electron beam passage below the target. It is important to mention that two impact-parameters, the horizontal and the vertical (h < 0 for the CTR and h > 0 for the CDR), are considered as explained in the Methods section. Figure 2 shows a comparison of the typical charge–normalized auto–correlation dependencies (a-c), measured spectra (d) taken for the same observation angle and angular acceptance of the detection system , and GDR orientation dependencies (e). As can be seen in Fig. 2d, the Schottky Barrier Diode (SBD) detector has a full bandwidth of GHz[25,27] which limits the overall spectral sensitivity of the system. For a direct spectral intensity comparison, the spectrum of the coherent SPR, horizontal polarization was measured when the electron beam passed near the grating in SPR geometry: mm, mm, . Then, the spectra of GDR at for both horizontal polarization (HP) and vertical polarization (VP), were acquired. As expected, the yield of GDR at is smaller than that of SPR, since only part of the beam’s Coulomb field interacts with the grating. The integration over a given spectral range gives the following radiated power in arbitrary units: SPR-HP , GDR-HP and GDR-VP . The GDR emitted by the semi-plane grating has approximately times less intensity than the SPR at . Nevertheless, the observed spectral properties confirm that GDR is also monochromatic and its line widths for small angles and given practically coincide with the same characteristics of SPR. For different grating orientation angles , the GDR integral intensity of a few diffraction orders varies and reaches levels comparable to these of SPR as can be seen from the -scans acquired in the range from to while keeping mm, Fig. 2e.

Figure 2

Zoom-in of the measured auto–correlation curves of CTR (a), SPR (b) , mm, mm, and GDR horizontal polarization (c) for , . () Comparison of the CTR spectrum and SPR, GDR spectra taken for mm and mm. (e) GDR horizontal and vertical polarization angular distributions taken for mm.

Investigation of the GDR properties on the grating orientation angle shows wide–range spectral tunability. Different polarizations of the GDR have similar spectral content with minor variation in intensity (Fig. 3, for both polarizations) due to the –wedged sapphire vacuum window and different splitting efficiency of the m-thick silicon beam splitter used in the interferometer. The simulated dependence of , where and are the reflection and transmission coefficients for two linear polarizations[28,29], is shown in Fig. 4,b. In the case of ideal splitting .

Figure 3

Typical GDR horizontal polarization (a) and vertical polarization (b) spectra as a function of the grating orientation angle .

Figure 4

GDR spectral–angular distribution of diffraction orders calculated by Eqs. 1, 2 (a - colored bows) and measured GDR spectral peaks for both horizontal and vertical polarizations. Error-bars represent FWHM spectral line widths (a - markers). (b) Calculated efficiency of the m-thick silicon beam splitter, angle of incidence .

An important GDR spectral feature is that not only central frequencies can be fine–tuned by selecting the angle, but also the spectral bandwidth for the same central frequency can be broadened by switching to lower diffraction orders at larger . Typical frequency shifts for sequentially small angles are shown for GDR VP in Fig. 3b and the bandwidth increase near similar central frequencies are shown for GDR HP in Fig. 3a. In the current example, the bandwidth was switched from to levels. According to splitter efficiency simulation Fig. 4b, drops to zero at GHz and GHz. Around these frequencies Michelson interferometer gradually loses performance, but still produces auto–correlation with reduced visibility which can be recalculated to a spectrum via Fourier transform. However, spectral amplitude in this case is decreasing and can not be directly compared with amplitudes from other regions of dependence. The red dashed line in Fig. 3a represents the re–normalized GDR HP spectra for which is valid to within an amplitude factor. The renormalization means addition of inversed spectral function multiplied by inversed splitter efficiency function, as can be derived from a generalized formula of Fourier transform spectroscopy[30,31]. This is done only for the GDR spectrum taken at in the range 430–470 GHz as other spectral lines lie in non-zero regions of dependence.

A summary of the GDR spectra measurements for both HP and VP is shown in Fig. 4a along with the GDR spectral–angular distribution for diffraction orders calculated using Eq. 1. The full width at half maximum (FWHM) spectral widths are calculated by Eq. 2. Measured data are in close agreement with the dispersion relation Eq. 1.

Discussion

It is expected that for the THz frequency range ( mm) electrons with energy MeV can generate monochromatic GDR from a grating with mm and . In this case the spectral line may be tuned over a broad interval, thereby changing the grating orientation angle in the range with an accuracy of can result in a spectral range of 0.1–1 THz and a spectral accuracy ranging from sub-GHz to a few GHz depending on and . Additional spectral selection can be performed by usage of bandpass filters[32] or different grating profiles[33,34].

The GDR radiated energy can be estimated through comparison of the GDR and CDR experimental data while CDR yield can be analytically found. Using the notation from ref. [35] one can write the following expression for CDR HP spectral density assuming a perfectly conducting semi-plane target and neglecting the terms lower than as:Here is the 3D Gaussian beam form–factor, is the characteristic diffraction radiation frequency, - radiation frequency, - vertical impact parameter, - Lorentz factor, - fine-structure constant, - speed of light, and are the projection angles between the specular reflection direction and the wave vector of the outgoing photon. Simulation of the CDR -scans by Eq. 3 is presented in Fig. 5a - red curves which shows a reasonable agreement with the experimental data for the real electron bunch sizes and the detector angular acceptance neglecting electron beam divergence. To find the CDR HP emitted energy per bunch Eq. 3 should be integrated over the detected radiation frequency band as:

Figure 5

(a) Measured and simulated CDR angular distributions. (b) Measured and simulated bremsstrahlung and CDR VP intensity as a function of the vertical impact parameter . (c) Measured and simulated coherent SPR HP intensity versus horizontal impact parameter .

Assuming an angular acceptance of the detection system , bunch population , measured SBD detector sensitivity bandwidth GHz ( GHz, GHz), electron bunch rms sizes in transverse and longitudinal directions   µm and µm we obtain pJ and, accordingly pJ/THz. The results of the GDR HP measurements (Fig. 2e) are presented in the same arbitrary units as the CDR yield (Fig. 5a). For the measured intensity of the GDR HP is about of the CDR HP taken at . For both cases mm. The spectral intensities ratio of GDR HP ( GHz) to ( GHz) in this case is or : , (Fig. 3a). Hence, the intensity of the single GDR spectral line GHz may be estimated as of the and taking into account the GDR line width ( GHz) one can obtain the estimation of the GDR spectral density as pJ/THz. This means, that the intensity of a single coherent GDR line with GHz from a train of short electron bunches with a total charge pC nC can achieve the level of nJ/THz. Such a source of monochromatic THz radiation based on a compact electron accelerator can be considered as a promising candidate for many practical applications. The errors include statistical errors due to the fitting or rms calculation in , , and systematic errors from the uncertainty in .

Wide–range tunability and the possibility of bandwidth selection, along with further grating and electron beam parameter optimization will lead to a much higher peak radiation power. Due to the unperturbed interaction of the electron beam with the grating, one can consider GDR applications in time-resolved THz spectroscopy experiments[36], to non–invasive beam diagnostics[37-39] and multi–colour radiation generation[40] by an array of gratings. To increase the GDR radiated power further, a slit–grating (i.e. when the electron beam is passing through a horizontal slit made in a grating) with an optimized profile should be considered.

Methods

The experiment was performed at the KEK LUCX facility[41]. The detailed description of the accelerator, grating, and THz interferometer can be found in ref. [25]. As it was discussed in ref. [42], the Michelson interferometer is optimized for a large bandwidth of both HP (polarization in the diffraction plane) and VP of the incoming radiation excluding only GHz and GHz spectral lines where silicon splitter efficiency drops to zero. The SBD detector and a mm diameter wire-grid polarizer consisted of tungsten wires with diameters of µm and spacings of µm (installed in front of the SBD) were mounted to the rotation stage to allow for polarization selection during spectral measurements. All spectra (Figs. 2d and 3) were obtained by Fourier transform[43,44] of the charge–normalized auto–correlation curves which were measured with an rms resolution determined by the interferometer’s movable arm travel range mm (zero–path difference point mm) as and so are plotted with identical scales. The Michelson interferometer used in the experiment has a resolution higher than that of the natural GDR VP spectral line width determined for a point-like aperture , since , . However, the angular acceptance of the detection system limited the measured spectral line widths to level. The observed anomaly of the curve around GHz in Fig. 3a suggests to perform a higher resolution study to check for a more complex spectral structure.

Initially, the interferometer was set to a zero-path difference point and a number of preliminary scans were performed in order to verify the grating angle and the position with respect to the electron beam. The CDR yield versus the angle between the electron beam and the flat side on the rear of the grating (-scan for a large angles around ) was measured for impact parameter mm, Fig. 5a. This allows the grating orientation angles to be determined with an accuracy better than .

It was noted that the current experimental geometry for the investigation of the GDR properties has two impact parameters: vertical - the distance between the electron beam and the edge of the grating; and horizontal which can be defined as the distance between the electron beam and the grating main plane at , Fig. 1. For corresponding to the maximum of CDR VP, (specular reflection direction from the flat rear side of the grating) the radiation yield dependence on the vertical impact parameter was obtained (Fig. 5b - green markers). The peak of the curve coincides with the target edge, the right slope corresponds to the CDR, and the left part is connected with CTR. A similar dependence for non-coherent radiation was measured in the optical wavelength range at the KEK ATF accelerator, and was discussed in ref. [45]. In order to avoid direct beam interaction with the grating and to maintain a significant CDR/GDR yield, the vertical impact-parameter for further measurements was chosen as mm.

Electron beam parameter simulations[46,47] show that the transverse rms bunch size at the grating location was equal to µm. The measurements of the bremsstrahlung produced by electrons versus the vertical impact parameter confirm the simulation results (Fig. 5b - black markers). If the assumption is made that the transverse profile of the beam can be described by a Gaussian function (confirmed by measurements with luminescent screen installed mm downstream of the experimental station) the measured dependence can be described by the function:with free parameters mm and µm.

The horizontal impact parameter mm was set in accordance with the radiation yield dependence on the horizontal distance between the electron beam and the grating in SPR geometry: , mm, Fig. 5c. If is known, this allows for a comparison of yields at for GDR and SPR to be made.