On-chip high energy photon radiation source based on microwave-dielectric undulator

Fuming Jiang State Key Laboratory of High Field Laser Physics and CAS Center for Excellence in Ultra-intense Laser Science, Shanghai Institute of Optics and Fine Mechanics (SIOM), Chinese Academy of Sciences (CAS), Shanghai 201800, China University of Chinese Academy of Sciences, Beijing 100049, China    Xinyu Xie State Key Laboratory of High Field Laser Physics and CAS Center for Excellence in Ultra-intense Laser Science, Shanghai Institute of Optics and Fine Mechanics (SIOM), Chinese Academy of Sciences (CAS), Shanghai 201800, China University of Chinese Academy of Sciences, Beijing 100049, China    Chengpu Liu chpliu@siom.ac.cn State Key Laboratory of High Field Laser Physics and CAS Center for Excellence in Ultra-intense Laser Science, Shanghai Institute of Optics and Fine Mechanics (SIOM), Chinese Academy of Sciences (CAS), Shanghai 201800, China University of Chinese Academy of Sciences, Beijing 100049, China    Ye Tian tianye@siom.ac.cn State Key Laboratory of High Field Laser Physics and CAS Center for Excellence in Ultra-intense Laser Science, Shanghai Institute of Optics and Fine Mechanics (SIOM), Chinese Academy of Sciences (CAS), Shanghai 201800, China University of Chinese Academy of Sciences, Beijing 100049, China
Abstract

A new on-chip light source configuration has been proposed, which utilizes the interaction between microwave and a dielectric nanopillar array to generate a periodic electromagnetic near field, and applies periodic transverse acceleration to relativistic electrons to generate high-energy photon radiation. Here the dielectric nanopillar array interacting with microwave acts as the electron undulator, in which the near field drives electrons to oscillate. When an electron beam operates in this nanopillar array in this light source configuration, it is subjected to a periodic transverse near-field force, and will radiate X-ray or even γ𝛾\gammaitalic_γ-ray high energy photons after a relativistic frequency up-conversion. Compared with the laser-dielectric undulator based on the interaction between strong lasers and nanostructures to generate a plasmonic near field, this configuration is less prone to damage during operation.

I Introduction

Traditional solid and gas lasers cannot directly generate X-Ray or γ𝛾\gammaitalic_γ-Ray. In addition to decay radiation sources, high-energy electronic radiation can effectively generate such high-energy photons. For example, using electrons emitted by electron guns to bombard targets can generate X-ray through bremsstrahlung, or injecting relativistic electrons into undulators can also generate soft/hard X-ray. The former is the oldest X-ray generation plan in history, widely used in medicine, industrial production, and other fields, while the latter is free electron laser (FEL), currently one of the most cutting-edge topics. However, X-ray generated by bremsstrahlung have some defects, such as large divergence angle and poor coherence. And traditional FEL requires a large array of undulatory magnets, this large volume of strong magnetic system needs extremely short (micrometer level) oscillation period, so it is difficult to establish such a device. If we want to use FEL system to generate X-ray or even Gamma-ray, there are two direct ways: (1) we can accelerate the free electrons to extreme relativistic region, so the shortwave radiation can be generated by these high-energy electrons via relativistic frequency up-conversion effect; (2) We can also build a longer undulator system to generate high-order harmonics and then amplify the output. Both of these plans will make FEL systems larger and more expensive, and such complex precision systems will also reduce their reliability and increase the difficulty of construction and maintenance [1].

To solve the above problems, many scientists have taken different approaches and attempted non-magnetic array undulator configurations. Elias proposed the concept of "electromagnetic wave undulator" in 1979, pointing out that electromagnetic waves can interact with free electrons to generate a wave period much smaller than a magnetic array does [2]. The most intuitive principle is the inverse Compton scattering (ICS) gamma laser, where electrons interact with high-energy laser pulses and undergo transverse acceleration under the laser electromagnetic field, then stimulating radiation. At this point, the oscillation period of the electron is close to the oscillation period of the laser pulse, which is much shorter than the traditional undulator period. From a quantum perspective, it can be understood that a free electron undergoes ICS with a low-energy photon, generating a high-energy photon. MIT proposed an ICS gamma laser system based on a linear accelerator in 2009 [3]. In 2012, K. Ta Phuoc et al. proposed a "fully optical Compton gamma light source", where the generation, acceleration, and oscillation of free electrons are all achieved by a single light source [4]. With the development of Laser Wake Field Acceleration (LWFA), it has been found that when electrons are accelerated inside plasma bubbles, generated by ponderomotive force, they will produce cyclotron betatron radiation, with significantly shorter oscillation periods than solid undulators. Electrons, without extremely high energy, can also produce strong X-rays [5, 6].

In addition, the concept of "on-chip FEL light source" is also one of the cutting-edge types of research in recent years. There are two technical routes in this plan: (1) surface plasmon polaritons (SPPs) amplification, (2) and micro/nano-structure modulation of electrons. The former includes on-chip light sources that interact with free electrons and SPPs [7, 8], while the latter includes Smith Purcell radiation generated by the interaction between electrons and micro/nano gratings [9, 10, 11], and periodic near field generated by micro periodic structures effect electron beams [12, 13, 14, 15]. Although on-chip light sources based on Smith Purcell radiation can also produce X-rays [16], their directionality and monochromaticity are generally inferior to near-field undulators.

With the improvement of micro/nano-processing technology, micro/nano-structures with periodic of tens of nanometers can be prepared. By using appropriate external pumping, near field can be generated around these periodic micro/nano-structures. When an electron beam passes through the structure, it is subjected to the near field, generating periodic acceleration and radiation.

In this letter, we propose a model that utilizes the interaction between microwaves and nano-sized structures to generate quasi-static periodic near field, thereby inducing electron oscillations to generate X-Ray or even γ𝛾\gammaitalic_γ-Ray, and use simulation to obtain the radiation’s characters of the device.

II PERIODIC NEAR-FIELD CHANNEL

Refer to caption
FIG 1: High energy photon radiation source based on microwave dielectric undulator. (a) The interaction between a dielectric nanopillar array and a resonant cavity (purple square box) with a microwave standing wave (green wavefront) polarized along the y-axis with an amplitude of 80 MV/m is depicted in (b), which produces a periodic transverse near field as shown in (b). (c) is the detail view of (b), which shows that the nanopillar array generates a polarization field ,with the effect of microwaves, and the white arrows show the electric field lines around the periodic near-field channel. (d) The Feynman diagram of electrons oscillating and generating radiation in this structure.

One of the best ways to induce small periodic oscillations of free electrons under the influence of an external field is to use electromagnetic near-field modulation. For a dielectric medium in an electric field with no surface charge, if the medium is in a vacuum environment, it is easy to obtain (ε1E1ε0E2)en=0subscript𝜀1subscript𝐸1subscript𝜀0subscript𝐸2subscript𝑒𝑛0(\varepsilon_{1}\cdot\vec{E}_{1}-\varepsilon_{0}\cdot\vec{E}_{2})\cdot\vec{e}_% {n}=0( italic_ε start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_E end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - italic_ε start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_E end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) ⋅ over→ start_ARG italic_e end_ARG start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT = 0 and (E1E2)×en=0subscript𝐸1subscript𝐸2subscript𝑒𝑛0(\vec{E}_{1}-\vec{E}_{2})\times\vec{e}_{n}=0( over→ start_ARG italic_E end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - over→ start_ARG italic_E end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) × over→ start_ARG italic_e end_ARG start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT = 0 from Gaussian theorem, where ε0subscript𝜀0\varepsilon_{0}italic_ε start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and ε1subscript𝜀1\varepsilon_{1}italic_ε start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT are the vacuum’s and the dielectric’s dielectric constant, E1subscript𝐸1\vec{E}_{1}over→ start_ARG italic_E end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and E2subscript𝐸2\vec{E}_{2}over→ start_ARG italic_E end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT are the electric fields inside and outside the dielectric, and ensubscript𝑒𝑛\vec{e}_{n}over→ start_ARG italic_e end_ARG start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT is the normal vector of the interface. It can be seen that if a negative y-axis electric field is applied to the dielectric cylindrical array in vacuum [see FIG. 1 (a)], electric field distortion will occur near the cylinder [see FIG. 1 (b)]. This is because the dielectric will generate displacement polarization under the effect of an external electric field, and dipole enhancement will also occur between array elements. The most noteworthy aspect is that the induced electric field of the dielectric column contains an x𝑥xitalic_x-component electric field which is not included in the external electric field [red and blue blocks in FIGs. 1 (b) and (c)].

By using micro/nano-processing technology, dielectric nanorods with a diameter of 40 nm are arranged in two staggered columns along the y𝑦yitalic_y-axis with a center spacing of 50 nm, resulting in a x𝑥xitalic_x-component periodic staggered near-field channel along the y𝑦yitalic_y-axis [see FIGs. 1 (a) and (b)]. This near-field channel will impose a periodic transverse oscillation electric field force on the injected electrons, which can make them oscillate with extremely short periods equal to 50 nm, then radiating short-wavelength photons.

In order to apply an electric field in the y𝑦yitalic_y-axis on a dielectric nanopillar array, we cannot simply use a static electric field, because a high-voltage static electric field can easily breakdown the medium. So we use the radio frequency (RF) microwave standing wave to induce near field. By placing the dielectric nanopillars at specific positions in the RF cavity, the polarization of the microwave is parallel to the y𝑦yitalic_y-axis, the channel is at the antinode of the electric wave and the nodal position of the magnetic wave, in order to generate near field which will drive electrons lateral oscillate and minimize the influence of the background magnetic field as well [see FIG. 2 (b)]. When electromagnetic wave passes through an array of nanopillars which orthogonal to its polarization direction, if the nanopillars’ diameter and spacing period are far less than the wavelength, the electromagnetic wave will bypass the arry due to diffraction effects and do not hinder its propagation. Therefore, it can form a stable standing wave field in this structure. In addition, the nano dielectric materials’ low absorption rate of RF photons makes them less susceptible to be damaged under high-energy microwave heating.

The displacement polarization of bound electrons in a dielectric requires a certain response time, which is generally in the range of femtoseconds to sub-nanoseconds, much shorter than the oscillation period of RF microwave. Therefore, this near field is relatively stable. The RF cavity of an electron induction accelerator can not sustain an electromagnetic field which strength is higher than the GV/m level [17, 18]. We hope that electrons will experience a near-field electric field strength of 10 MV/m on the x𝑥xitalic_x-axis in the oscillation channel. If the nanopillar array is made of silicon, the standing wave amplitude should be about 80 MV/m. This is feasible in engineering and there is no risk of thermal damage in the short term. If dielectric ceramic materials with higher dielectric coefficients and higher thermal damage threshold are used, the strength of the standing wave can be further reduced while the near-field electric field strength remains unchanged, or a stronger near field can be generated by using a stronger driving field to enhance the power of the light source.

It is worth noting that when the direction of the microwave electric field is reversed, the direction of the induced near field will also be reversed, which will cause the oscillation of electrons a sudden change with the π𝜋\piitalic_π phase. This is quite unfavorable for the monochromaticity of electron radiation. Therefore, considering that the microwave resonant cavity and microwave source of RF accelerators widely use S-band microwave (with a wavelength of about 10 cm), the length of the near-field channel of the device should be less than 5 cm, so that the direction of the electric field remains unchanged when electrons travel through the channel. In fact, in order to ensure that the driving near field intensity is sufficiently high, the near-field channel’s length is preferably less than 1 cm, allowing electrons to be injected into the device with the highest microwave electric field.

Electrons are subjected to both transverse near-field oscillations and longitudinal acceleration driven by the background microwave electric field. But, for relativistic electrons, the total increment of their longitudinal velocity when passing through the oscillation channel does not exceed the order of 105 m/s, which can be ignored compared to the longitudinal velocity close to the speed of light and the corresponding energy gain is not higher than the energy dissipation of the high-quality electron beam [19]. Moreover, there is no necessary precise phase-locked condition between electrons and electric field in this model, such small changes in longitudinal velocity will not have a significant impact on the radiation output.

In addition, due to the long RF wavelength, comparing with the length of the near-field channel, the near-field can be regarded as an approximate electrostatic field when the relativistic electrons undergo several oscillation cycles, which is used in the mathematical derivation in the following text.

III ELECTRONIC OSCILLATION RADIATION MODEL

Refer to caption
FIG 2: Dependence of electron radiation on azimuthal angle θ𝜃\thetaitalic_θ for electron energies of 1 MeV and 10 MeV.

According to the Lienart-Wiechart potential, the instantaneous electric field strength of a moving electron at the observation point with distance R𝑅Ritalic_R can be expressed as (Gaussian unit) [20],

E(x,t)=e[nβk3γ2R2]ret+ec[n×(nβ)×β˙k3R]ret.𝐸𝑥𝑡𝑒subscriptdelimited-[]𝑛𝛽superscript𝑘3superscript𝛾2superscript𝑅2ret𝑒𝑐subscriptdelimited-[]𝑛𝑛𝛽˙𝛽superscript𝑘3𝑅ret\vec{E}\left(\vec{x},t\right)=e\left[\frac{\vec{n}-\vec{\beta}}{k^{3}\gamma^{2% }R^{2}}\right]_{\text{ret}}+\frac{e}{c}\left[\frac{\vec{n}\times\left(\vec{n}-% \vec{\beta}\right)\times\dot{\vec{\beta}}}{k^{3}R}\right]_{\text{ret}}.over→ start_ARG italic_E end_ARG ( over→ start_ARG italic_x end_ARG , italic_t ) = italic_e [ divide start_ARG over→ start_ARG italic_n end_ARG - over→ start_ARG italic_β end_ARG end_ARG start_ARG italic_k start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_R start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ] start_POSTSUBSCRIPT ret end_POSTSUBSCRIPT + divide start_ARG italic_e end_ARG start_ARG italic_c end_ARG [ divide start_ARG over→ start_ARG italic_n end_ARG × ( over→ start_ARG italic_n end_ARG - over→ start_ARG italic_β end_ARG ) × over˙ start_ARG over→ start_ARG italic_β end_ARG end_ARG end_ARG start_ARG italic_k start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_R end_ARG ] start_POSTSUBSCRIPT ret end_POSTSUBSCRIPT . (1)

Where e𝑒eitalic_e is elementary charge, n𝑛\vec{n}over→ start_ARG italic_n end_ARG is the unit vector from the electron to the observation point, β𝛽\vec{\beta}over→ start_ARG italic_β end_ARG is the normalized velocity of the electron, k=1nβ(t)=dt/dt𝑘1𝑛𝛽superscript𝑡d𝑡dsuperscript𝑡k=1-\vec{n}\cdot\vec{\beta}\left(t^{\prime}\right)=\text{d}t/\text{d}t^{\prime}italic_k = 1 - over→ start_ARG italic_n end_ARG ⋅ over→ start_ARG italic_β end_ARG ( italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) = d italic_t / d italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT is the wave vector, γ=1β2𝛾1superscript𝛽2\gamma=\sqrt{1-\beta^{2}}italic_γ = square-root start_ARG 1 - italic_β start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG is the Lorentz factor, []retsubscriptret[]_{\textbf{ret}}[ ] start_POSTSUBSCRIPT ret end_POSTSUBSCRIPT is the retarded bracket, and c𝑐citalic_c is the light speed in the vacuum. Only the second item can radiate to a distance. It can be seen that when electrons are subjected to external forces, they will emit electromagnetic waves transmitting perpendicular to the direction of the force.

When a relativistic electron traveling in a straight channel line under an extremely short period of transverse reciprocating electric field, it can be seen from Eq. (1) that this kind of electric field, with extremely short period and not very high intensity, will cause the electron to apply a periodic oscillating electric field (i.e. radiate electromagnetic waves) to the forward observation point, and the electron’s trajectory will hardly produce significant deflection.

In this case, β𝛽\vec{\beta}over→ start_ARG italic_β end_ARG approaches 1 and has a direction similar to n𝑛\vec{n}over→ start_ARG italic_n end_ARG, i.e βnsimilar-to𝛽𝑛\vec{\beta}\sim\vec{n}over→ start_ARG italic_β end_ARG ∼ over→ start_ARG italic_n end_ARG, and k=1nβ(t)1/(2γ2)𝑘1𝑛𝛽superscript𝑡12superscript𝛾2k=1-\vec{n}\cdot\vec{\beta}(t^{\prime})\approx 1/(2\gamma^{2})italic_k = 1 - over→ start_ARG italic_n end_ARG ⋅ over→ start_ARG italic_β end_ARG ( italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) ≈ 1 / ( 2 italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ), nβn(11+12γ2)=12γ2n𝑛𝛽𝑛1112superscript𝛾212superscript𝛾2𝑛\vec{n}-\vec{\beta}\approx\vec{n}(1-1+{1\over{2{\gamma^{2}}}})={1\over{2{% \gamma^{2}}}}\vec{n}over→ start_ARG italic_n end_ARG - over→ start_ARG italic_β end_ARG ≈ over→ start_ARG italic_n end_ARG ( 1 - 1 + divide start_ARG 1 end_ARG start_ARG 2 italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) = divide start_ARG 1 end_ARG start_ARG 2 italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG over→ start_ARG italic_n end_ARG. By substituting the second term of Eq. (1), the electric field of the radiation field can be obtained,

E(x,t)=ec[n×(12γ2n×β˙)(12γ2)3R]ret=ec[n×(n×β˙)R4γ4]ret,𝐸𝑥𝑡𝑒𝑐subscriptdelimited-[]𝑛12superscript𝛾2𝑛˙𝛽superscript12superscript𝛾23𝑅ret𝑒𝑐subscriptdelimited-[]𝑛𝑛˙𝛽𝑅4superscript𝛾4ret\vec{E}(\vec{x},t)={e\over c}{\left[{{{\vec{n}\times({1\over{2{\gamma^{2}}}}% \vec{n}\times\dot{\vec{\beta}})}\over{{{({1\over{2{\gamma^{2}}}})}^{3}}R}}}% \right]_{{\rm{ret}}}}={e\over c}{\left[{{{\vec{n}\times(\vec{n}\times\dot{\vec% {\beta}})}\over R}4{\gamma^{4}}}\right]_{{\rm{ret}}}},over→ start_ARG italic_E end_ARG ( over→ start_ARG italic_x end_ARG , italic_t ) = divide start_ARG italic_e end_ARG start_ARG italic_c end_ARG [ divide start_ARG over→ start_ARG italic_n end_ARG × ( divide start_ARG 1 end_ARG start_ARG 2 italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG over→ start_ARG italic_n end_ARG × over˙ start_ARG over→ start_ARG italic_β end_ARG end_ARG ) end_ARG start_ARG ( divide start_ARG 1 end_ARG start_ARG 2 italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_R end_ARG ] start_POSTSUBSCRIPT roman_ret end_POSTSUBSCRIPT = divide start_ARG italic_e end_ARG start_ARG italic_c end_ARG [ divide start_ARG over→ start_ARG italic_n end_ARG × ( over→ start_ARG italic_n end_ARG × over˙ start_ARG over→ start_ARG italic_β end_ARG end_ARG ) end_ARG start_ARG italic_R end_ARG 4 italic_γ start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT ] start_POSTSUBSCRIPT roman_ret end_POSTSUBSCRIPT , (2)

the corresponding magnetic field as,

B(x,t)=[n×E(x,t)]ret,𝐵𝑥𝑡subscriptdelimited-[]𝑛𝐸𝑥𝑡ret\vec{B}(\vec{x},t)={\left[{\vec{n}\times\vec{E}(\vec{x},t)}\right]_{{\rm{ret}}% }},over→ start_ARG italic_B end_ARG ( over→ start_ARG italic_x end_ARG , italic_t ) = [ over→ start_ARG italic_n end_ARG × over→ start_ARG italic_E end_ARG ( over→ start_ARG italic_x end_ARG , italic_t ) ] start_POSTSUBSCRIPT roman_ret end_POSTSUBSCRIPT , (3)

and the Poynting vector as,

S(x,t)=c4π[E(x,t)×B(x,t)]ret=c4π[|E|2n]ret.𝑆𝑥𝑡𝑐4𝜋subscriptdelimited-[]𝐸𝑥𝑡𝐵𝑥𝑡ret𝑐4𝜋subscriptdelimited-[]superscript𝐸2𝑛ret\vec{S}(\vec{x},t)={c\over{4\pi}}{\left[{\vec{E}(\vec{x},t)\times\vec{B}(\vec{% x},t)}\right]_{{\rm{ret}}}}={c\over{4\pi}}{\left[{{{\left|{\vec{E}}\right|}^{2% }}\vec{n}}\right]_{{\rm{ret}}}}.over→ start_ARG italic_S end_ARG ( over→ start_ARG italic_x end_ARG , italic_t ) = divide start_ARG italic_c end_ARG start_ARG 4 italic_π end_ARG [ over→ start_ARG italic_E end_ARG ( over→ start_ARG italic_x end_ARG , italic_t ) × over→ start_ARG italic_B end_ARG ( over→ start_ARG italic_x end_ARG , italic_t ) ] start_POSTSUBSCRIPT roman_ret end_POSTSUBSCRIPT = divide start_ARG italic_c end_ARG start_ARG 4 italic_π end_ARG [ | over→ start_ARG italic_E end_ARG | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT over→ start_ARG italic_n end_ARG ] start_POSTSUBSCRIPT roman_ret end_POSTSUBSCRIPT . (4)

Observing at a small angle near the axis of electronic motion, the radiation power per unit solid angle is:

dP(t)dΩ=c4π[R2|E|2]ret=e24πc4γ4|n×(n×β˙(t))|2=e2γ4πc|n×(n×β˙(t))|2.𝑑𝑃𝑡𝑑Ω𝑐4𝜋subscriptdelimited-[]superscript𝑅2superscript𝐸2retsuperscript𝑒24𝜋𝑐4superscript𝛾4superscript𝑛𝑛˙𝛽𝑡2superscript𝑒2superscript𝛾4𝜋𝑐superscript𝑛𝑛˙𝛽𝑡2\begin{split}{{dP\left(t\right)}\over{d\Omega}}&={c\over{4\pi}}{{\rm{[}}{R^{2}% }{\left|{\vec{E}}\right|^{2}}{\rm{]}}_{{\rm{ret}}}}={{{e^{2}}}\over{4\pi c}}% \cdot 4{\gamma^{4}}\cdot{\left|{\vec{n}\times\left({\vec{n}\times\dot{\vec{% \beta}}\left(t\right)}\right)}\right|^{2}}\\ &={{{e^{2}}{\gamma^{4}}}\over{\pi c}}{\left|{\vec{n}\times\left({\vec{n}\times% \dot{\vec{\beta}}\left(t\right)}\right)}\right|^{2}}.\end{split}start_ROW start_CELL divide start_ARG italic_d italic_P ( italic_t ) end_ARG start_ARG italic_d roman_Ω end_ARG end_CELL start_CELL = divide start_ARG italic_c end_ARG start_ARG 4 italic_π end_ARG [ italic_R start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT | over→ start_ARG italic_E end_ARG | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] start_POSTSUBSCRIPT roman_ret end_POSTSUBSCRIPT = divide start_ARG italic_e start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 4 italic_π italic_c end_ARG ⋅ 4 italic_γ start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT ⋅ | over→ start_ARG italic_n end_ARG × ( over→ start_ARG italic_n end_ARG × over˙ start_ARG over→ start_ARG italic_β end_ARG end_ARG ( italic_t ) ) | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL = divide start_ARG italic_e start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_γ start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG start_ARG italic_π italic_c end_ARG | over→ start_ARG italic_n end_ARG × ( over→ start_ARG italic_n end_ARG × over˙ start_ARG over→ start_ARG italic_β end_ARG end_ARG ( italic_t ) ) | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT . end_CELL end_ROW (5)

For relativistic electrons, there is a conversion relationship between the electromagnetic wave frequency ν(θ)𝜈𝜃\nu\left(\theta\right)italic_ν ( italic_θ ) radiated by them and the electron oscillation frequency ν0subscript𝜈0{\nu_{0}}italic_ν start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT on their oscillation plane (x𝑥xitalic_x-y𝑦yitalic_y plane) [21]: ν(θ)2γ2ν0/(1+γ2θ2)𝜈𝜃2superscript𝛾2subscript𝜈01superscript𝛾2superscript𝜃2\nu\left(\theta\right)\approx 2{\gamma^{2}}{\nu_{0}}/(1+{\gamma^{2}}{\theta^{2% }})italic_ν ( italic_θ ) ≈ 2 italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ν start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT / ( 1 + italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_θ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ). That is, the relativistic frequency up-conversion relationship [see FIG. 2]. Among them, θ𝜃\thetaitalic_θ is the angle between the radiation direction and the direction of electron motion, which can be approximated as the angle between the radiation direction and the y-axis in this model [see FIG. 1 (a)]. For relativistic electrons with higher energy, their radiation frequency has a weak dependence on the radiation direction and the angle between the z𝑧zitalic_z-axis, so can be approximated as independent.

If electric polarization near field and photon number are used to describe the power, according to the up-conversion relationship, Eq. (5) will be rewritten as:

dN(t)dΩ=e2γ4πc2|n×(n×EPγme)|212γ2hν0=e22πc2hν0|n×(n×EPme)|2.𝑑𝑁𝑡𝑑Ωsuperscript𝑒2superscript𝛾4𝜋superscript𝑐2superscript𝑛𝑛subscript𝐸𝑃𝛾subscript𝑚𝑒212superscript𝛾2subscript𝜈0superscript𝑒22𝜋superscript𝑐2subscript𝜈0superscript𝑛𝑛subscript𝐸𝑃subscript𝑚𝑒2\begin{split}{{dN\left(t\right)}\over{d\Omega}}&={{{e^{2}}{\gamma^{4}}}\over{% \pi{c^{2}}}}{\left|{\vec{n}\times\left({\vec{n}\times{{{{\vec{E}}_{P}}}\over{% \gamma{m_{e}}}}}\right)}\right|^{2}}{1\over{2{\gamma^{2}}h{\nu_{0}}}}\\ &={{{e^{2}}}\over{2\pi{c^{2}}h{\nu_{0}}}}{\left|{\vec{n}\times\left({\vec{n}% \times{{{{\vec{E}}_{P}}}\over{{m_{e}}}}}\right)}\right|^{2}}.\end{split}start_ROW start_CELL divide start_ARG italic_d italic_N ( italic_t ) end_ARG start_ARG italic_d roman_Ω end_ARG end_CELL start_CELL = divide start_ARG italic_e start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_γ start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG start_ARG italic_π italic_c start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG | over→ start_ARG italic_n end_ARG × ( over→ start_ARG italic_n end_ARG × divide start_ARG over→ start_ARG italic_E end_ARG start_POSTSUBSCRIPT italic_P end_POSTSUBSCRIPT end_ARG start_ARG italic_γ italic_m start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT end_ARG ) | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT divide start_ARG 1 end_ARG start_ARG 2 italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_h italic_ν start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL = divide start_ARG italic_e start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_π italic_c start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_h italic_ν start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG | over→ start_ARG italic_n end_ARG × ( over→ start_ARG italic_n end_ARG × divide start_ARG over→ start_ARG italic_E end_ARG start_POSTSUBSCRIPT italic_P end_POSTSUBSCRIPT end_ARG start_ARG italic_m start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT end_ARG ) | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT . end_CELL end_ROW (6)

Where EPsubscript𝐸𝑃{\vec{E}_{P}}over→ start_ARG italic_E end_ARG start_POSTSUBSCRIPT italic_P end_POSTSUBSCRIPT is the polarization field, and mesubscript𝑚𝑒m_{e}italic_m start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT is the electron mass. It can be seen that for relativistic electrons with γ1much-greater-than𝛾1\gamma\gg 1italic_γ ≫ 1, the number of photons radiated per unit time and unit solid angle is independent of the electron energy [FIG. 1 (d) is a schematic diagram of the radiation generated by this model].

IV RADIATION SPECTRA OF ELECTRONS WITH DIFFERENT ENERGIES IN THE UNDULATOR

In order to generate electromagnetic waves with shorter wavelengths, the energy of electrons should be high enough to generate strong relativistic frequency up-conversion effects. However, the emissivity of high-energy electron beams is high, making it difficult to pass through channels with a width of only about 30 nm without colliding with nanorods, which will cause ionization. Therefore, we will mainly discuss low-energy relativistic electrons with energies of 1similar-to\sim6 MeV and high-energy electrons with energies of 10similar-to\sim20 MeV. The lower energy electrons can be supplied from compact accelerators, RF electron guns, even on-chip electron sources, which can be easily obtained in laboratory.

Refer to caption
FIG 3: (a) The radiation spectra of single 1 MeV, 3 MeV, and 6 MeV electron in the device. (b) The amplification of the 1 MeV electron radiation spectrum in (a).

Due to the fact that the photon momentum of low-energy electron’s radiation is much lower than the electron’s own momentum, radiation damping can be disregarded. Use the particle tracking module of the finite element analysis (FEA) program to simulate the motion of electrons and inputs their acceleration into the Lienart-Wiechart potential calculation to obtain the electron radiation intensity. In this simulation, the device is composed of silicon columns with a diameter of 40 nm, and generates a near field of 108superscript10810^{8}10 start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT V/m under the drive of a 10 cm wavelength RF microwave. In common experimental environments, electron beams have a certain degree of energy dissipation, usually 1% or even lower, which can cause spectral broadening. However, for this model, the oscillation period of relativistic electrons is fixed, and their radiation photon frequency γ2proportional-toabsentsuperscript𝛾2\propto{\gamma^{2}}∝ italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT, while the energy of relativistic electrons γproportional-toabsent𝛾\propto\gamma∝ italic_γ. Therefore, the spectral broadening effect caused by electron dispersion should be less than 0.01% and can be ignored.

The photon energy of lower energy relativistic electron radiation is in the range of 0.1similar-to\sim10 keV [see FIG. 3]. The photon energy of 1 MeV electron radiation is about 0.2 keV, belonging to the extreme ultraviolet region [see FIG. 3 (b)], while the photons of 3 MeV and 6 MeV electron radiation are in the soft X-ray and X-ray region. From FIG. 3 (a), it can be observed that the number of photons radiated by 1 MeV electron per unit of time per unit of solid angle is much lower than that of 3, 6 MeV electrons, which deviates from the results of Eq. (6) to a certain extent. This is because the approximation mentioned earlier was made when the Lorentz factor γ1much-greater-than𝛾1\gamma\gg 1italic_γ ≫ 1, which is not accurate for lower energy electrons.

In order to obtain radiation with shorter wavelengths, we can use higher energy electron beams, but correspondingly, higher energy electron beams have higher electron emissivity and it is more difficult to compress the beam diameters. Due to the fact that nanorods can extend along the z𝑧zitalic_z-axis, elliptical electron beams compressed along the x𝑥xitalic_x direction can be used, but this process is also very difficult for extreme relativistic electrons. In order to reduce the possibility of ionization caused by electron beams colliding with nanorods, the current intensity must also significantly decrease, which means that the number of photons generated per unit time will also decrease. Therefore, an electron beam of 10similar-to\sim 20 MeV is an appropriate choice for generating hard X-band.

Refer to caption
FIG 4: (a) The radiation spectra of single 10 MeV, 15 MeV, and 20 MeV electrons in the device. (b) The amplification of the 10 MeV electron radiation spectrum in (a).

FIG. 4 shows the single electron radiation intensity of 10, 15, 20 MeV electrons. It can be seen that for electrons exceeding 10 MeV, their radiation spectrum is in the hard X-Ray region. FIG. 4 (a) shows that the number of photons per unit time and per unit solid angle emitted by strongly relativistic electrons with different energies is approximately equal, consistent with that given by Eq. (6).

Refer to caption
FIG 5: Transverse displacement of 1 MeV electron within 50 cycles. The figure shows the electron’s trajectory without initial lateral velocity modulation, the trajectory with initial lateral velocity modulation, and the theoretical motion trajectory. It can be seen that the lateral displacement of electrons is very small, and the initial lateral velocity modulation can greatly reduce the electron’s walk off.

Combining FIGs. 3 and 4, it can be seen that the electromagnetic wave spectrum radiated by a single energy electron beam in this device is a quasi-monochromatic spectrum with excellent monochromaticity. As mentioned earlier, compared with traditional FEL systems, the transverse displacement and the changes of velocity direction of electrons in this system are very small. FIG. 5 shows the transverse displacement of a 1 MeV electron during approximately 50 oscillation cycles, only at the femtometer level. For higher energy electrons, due to the increase in dynamic mass, their transverse displacement comes to the sub-femtometer level. Due to the relaxation process of the dipole polarization field and dipole enhancement generated by dielectric materials, the oscillation trajectory of electrons under the polarization field is not a standard sine curve [purple dashed line in FIG. 5], and often moves away to one side. If an appropriate transverse velocity is given to the electrons before they are injected into the undulator (such as through a parallel charged electrode plate), the walk-off effect can be greatly reduced, further compressing the transverse motion range of the electron.

The extremely small transverse oscillation amplitude causes the undulator factor of this device Ksimilar-to\sim0 [22], resulting in a very small divergence angle of radiation, which can be approximated as 1/γ1𝛾1/\gamma1 / italic_γ. Therefore, one can say that the higher energy the electron has, the smaller the divergence angle of radiation is. The output area of the device is only on the order of μ𝜇\muitalic_μm2, so the radiation brightness is quite high.

V CONCLUSION

We propose for the first time a mechanism based on the interaction between microwave and nanostructures to generate periodic oscillations of relativistic electrons in the near field, promoting the generation of extremely short wavelength electromagnetic radiation. The melting point of dielectric materials such as monocrystalline silicon and dielectric ceramics is mostly above 1000 C and have low absorption rate of RF photons, which make this configuration is less prone to thermal damage during operation.

Because the oscillation period of electrons in this model is entirely determined by the period of nanopillar array, high brightness X-ray radiation with excellent monochromaticity and minimal divergence angle can be generated by using an RF photocathode electron gun or a desktop accelerator system. Moreover, due to the many similarities between this structure and the structure of the Laser Dielectric Acceleration (DLA) device, both using an external driving electromagnetic wave field to irradiate periodic nanostructures and generate near-field modulation of electron motion. Therefore, after adjusting the period and arrangement of the nanostructure, higher energy laser irradiation can also be used to apply for higher acceleration to electrons in this device, thereby increasing the power of electron oscillation radiation. Furthermore, this structure can be etched on the same substrate with the DLA system, allowing the acceleration, regulation, and oscillation radiation processes of the electron beam to be completed on the same chip, within a distance of a few centimeters.

If laser is used instead of microwave as the driving source in the future, the electromagnetic damage threshold of the dielectric material itself is much higher than that of the metal resonant cavity. Therefore, lasers with higher electric field intensity can be used to generate stronger near field. This means that the power of electron radiation can be increased by several orders, or the width of the oscillation channel can be greatly expanded while maintaining the force intensity of electrons on the central axis. This means that extreme relativistic electrons with higher emissivity can also be injected into it to generate γ𝛾\gammaitalic_γ-ray or even hard γ𝛾\gammaitalic_γ-ray.

Acknowledgments: This work has been supported by National Natural Science Foundation of China (No. 12388102, 12074398, 12325409, U226720057); CAS Project for Young Scientists in Basic Research (YSBR-060) and Shanghai Pilot Program for Basic Research, Chinese Academy of Sciences, Shanghai Branch.

References

  • [1] D. Zhang, Y. Zeng, Y. Tian, and R. Li, Coherent free-electron light sources, Photonics Insights 2, R07 (2023).
  • [2] L. R. Elias, High-power, cw, efficient, tunable (uv through ir) free-electron laser using low-energy electron beams, Phys. Rev. Lett. 42, 977 (1979).
  • [3] W. S. Graves, W. Brown, F.X. Kaertner, and D.E. Moncton, MIT inverse Compton source concept, Nucl. Instrum. Methods Phys. Res. Sect. A-Accel. Spectrom. Dect. Assoc. Equ 608, S103 (2009).
  • [4] K. Ta Phuoc, S. Corde, C. Thaury, V. Malka, A. Tafzi, J.-P. Goddet, R.C. Shah, S. Sebban, and A. Rousse, All-optical Compton gamma-ray source, Nat. Photonics 6, 308 (2012).
  • [5] A. Rousse, K.T. Phuoc, R. Shah, A. Pukhov, E. Lefebvre, V. Malka, S. Kiselev, F. Burgy, J.-P. Rousseau, and D. Umstadter, Production of a keV X-ray beam from synchrotron radiation in relativistic laser-plasma interaction, Phys. Rev. Lett. 93, 135005 (2004).
  • [6] S. Kiselev, A. Pukhov, I. Kostyukov, X-ray generation in strongly nonlinear plasma waves, Phys. Rev. Lett. 93, 135004 (2004).
  • [7] D. Zhang, Y. Zeng, Y. Bai, Z. Li, Y. Tian, and R. Li, Coherent surface plasmon polariton amplification via free-electron pumping, Nature 611, 55 (2022).
  • [8] L. J. Wong, I. Kaminer, O. Ilic, J. D. Joannopoulos, and M. Soljačić, Towards graphene plasmon-based free-electron infrared to X-ray sources, Nat. Photonics 10, 46 (2016).
  • [9] H. Ishizuka, Y. Kawamura, K. Yokoo, H. Shimawaki, and A. Hosono, Smith–Purcell experiment utilizing a field-emitter array cathode: measurements of radiation, Nucl. Instrum. Methods Phys. Res. Sect. A-Accel. Spectrom. Dect. Assoc. Equ 475, 593 (2001).
  • [10] Y. Ye, F. Liu, M. Wang, L. Tai, K. Cui, X. Feng, W. Zhang, and Y. Huang, Deep-ultraviolet smith–purcell radiation, Optica 6, 592 (2019).
  • [11] W. Liu, M. Cao, Y. Wang, and K. Li, Start current of dielectric-loaded grating in Smith-Purcell radiation, Phys. Plasmas 23, (2016).
  • [12] I. A. Andriyash, R. Lehe, A. Lifschitz, C. Thaury, J.M. Rax, K. Krushelnick, and V. Malka, An ultracompact X-ray source based on a laser-plasma undulator, Nat. Commun. 5, 4736 (2014).
  • [13] T. Plettner, and R.L. Byer, Proposed dielectric-based microstructure laser-driven undulator, Phys. Rev. Spec. Top.-Accel. Beams 11, 030704 (2008).
  • [14] D. Rohrbach, C. B. Schroeder, A. Pizzi, R. Tarkeshian, M. Hayati, W.P. Leemans, and T. Feurer, THz-driven surface plasmon undulator as a compact highly directional narrow band incoherent x-ray source, Phys. Rev. Spec. Top.-Accel. Beams 22, 090702 (2019).
  • [15] Y. Tian, J. Liu, Y. Bai, S. Zhou, H. Sun, W. Liu, J. Zhao, R. Li, and Z. Xu, Femtosecond-laser-driven wire-guided helical undulator for intense terahertz radiation, Nat. Photonics 11, 242 (2017).
  • [16] M. Shentcis, A. K. Budniak, X. Shi, R. Dahan, Y. Kurman, M. Kalina, H. Herzig Sheinfux, M. Blei, M. K. Svendsen, and Y. Amouyal, Tunable free-electron X-ray radiation from van der Waals materials, Nat. Photonics 14, 686 (2020).
  • [17] T. Abe, T. Kageyama, H. Sakai, Y. Takeuchi, and K. Yoshino, Direct observation of breakdown trigger seeds in a normal-conducting rf accelerating cavity, Phys. Rev. Accel. Beams 21, 122002 (2018).
  • [18] R. J. England, R. J. Noble, K. Bane, D.H. Dowell, C.-K. Ng, J. E. Spencer, S. Tantawi, Z. Wu, R. L. Byer, and E. Peralta, Dielectric laser accelerators, Rev. Mod. Phys. 86, 1337 (2014).
  • [19] J. Wang, M. Zeng, D. Li, X. Wang, and J. Gao, High quality beam produced by tightly focused laser driven wakefield accelerators, Phys. Rev. Accel. Beams 26, 091303 (2023).
  • [20] J. D. Jackson, and R. F. Fox, Classical electrodynamics, 3rd ed., Am. J. Phys. 67, 841 (1999).
  • [21] P. Schmüser, M. Dohlus, J. Rossbach, and C. Behrens, Free-electron lasers in the ultraviolet and X-ray regime, Springer Tracts in Modern Physics 258, (2014).
  • [22] G. Dattoli, M. Del Franco, M. Labat, P.L. Ottaviani, and S. Pagnutti, Introduction to the physics of free electron laser and comparison with conventional laser sources, Free Electron Lasers 1, 20 (2012).