1. Introduction
Demands from applications promote the progress in preparing/generating ultrashort-duration, down to wavelength-level or shorter, laser pulse [
1,
2,
3,
4,
5,
6,
7,
8,
9,
10], which in turn promotes new applications areas appearing [
11,
12,
13,
14,
15,
16,
17,
18,
19,
20]. For applications emphasizing high energy density, spectral broadening technique, which relies on optical property of bulk material [
21,
22,
23,
24] and fine adjusting optical elements [
25,
26,
27], can supply single
m-level wavelength laser pulse with femtosecond(fs)-level duration and mJ-level energy. For applications emphasizing optics probing and optical manipulation of quantum system, the pulse is not required to be high-energy-density but to be of a shorter duration down to attosecond(as)-level. Current technique routes for achieving as-level duration are from (1) high-order harmonics generation (HHG) from laser-irradiated atoms in gaseous medium [
1,
2,
3,
4,
28,
29,
30,
31,
32,
33]. (2) HHG effect in laser–plasma interaction [
34,
35]. (3) The interaction between laser and electron beam free-electron-laser (FEL) facility [
36,
37,
38]. In these routes, Fourier components with shorter wavelength (below half of incident central wavelength) are produced from the interaction of the incident pulse with matter mentioned above. The ability of the pulse being shortened is determined by the efficiency of producing these Fourier components.
All these routes face the same difficulty, the efficiency. In the laser–atom route, the amplitude of higher-order harmonic component decreases exponentially with respect to the order. For example, the intensity of 50-order component is
of that of the fundamental component [
39]. Merely enhancing the intensity of the fundamental component is of little effect on this fact because the laser–atom route demands the survival of the involved quantum system, which does not favor the application of more intense fundamental components. The other two routes are based on the interaction of the fundamental component with classic electrons. Besides preparing plasmas, manipulating bulk plasma electrons to do a collective oscillation with a frequency multiple folds of the incident laser also demands the incident laser to be high-energy-density enough. The route based on FEL facility is also not optimistic in the efficiency because its background physics is associated with high-energy relativistic electrons and large-sized undulator magnetic field configuration, which are rather expensive for most university-level small research groups.
The commonality of these routes is to use easily-prepared electromagnetic (EM) Fourier components, such as visible and near-infrared laser with m-level wavelength (from m to m), as “raw material” to produce, through different technique routes, EM Fourier components with wavelength <0.1 m, which act as building components of desired pulse with as-level duration. Namely, producing the EM pulse with as-level duration is along an indirect way through the EM pulse with longer duration, or first producing a EM pulse with longer duration and then producing a EM pulse with as-level duration from this longer duration “raw material” supplied by a commercial available laser.
It is instructive and interesting to consider initially directly producing such a desired EM pulse with as-level duration. Retrospecting the history of laser, as well as its elder brother maser, we can find that the generation of EM pulse at visible and near-infrared band starts from flexible utilization of AC electric and magnetic field at power frequency (about tens Hz). A laser apparatus, due to electron structure of its kernel components/elements (i.e., gain medium), can convert power frequency input to optics-frequency output. The shortage of available material with electron structure suitable to 100 eV-level, or a level of tens nanometer (nm), wavelength optics output determines little hope, at least currently, of setting up a solid-state laser-like apparatus directly converting power frequency input to intense X-ray output. Of course, obtaining the as-level duration through plasma-based coherent X-ray source is obviously an expensive route.
However, if not being bound to suitable gain medium, we can try to consider the generation of radiation from low-energy electrons bunch in vacuum. The motive to this direction is for minimizing the dependence of generating radiations on kernel materials. Therefore, flexible usage of solid-state vacuum electron components/elements which are common in accelerator physics is the core of the technique route we consider.
Matured techniques in accelerator physics enables us to set up various configurations of DC magnetic field and DC electric field. These targeted designed field configurations can interact with incident electrons bunch to generate radiation. Many examples, such as FEL and synchronized radiation, have illustrated this point well. In FEL, a laser, or a EM optical-frequency, still need to be input and the incident electrons bunch is required to be sufficiently high-velocity. In contrast, the route presented in this work does not have these requirements.
The purpose is to seek for a cheap route of preparing EM pulse with as-level duration through the interaction of low-energy electrons bunch with targeted designed configuration of DC magnetic and electric fields. Based on previous works [
40,
41], it is feasible to achieve this goal because there is no severe requirement on components/elements. All components/elements are common in vacuum electron devices and working voltage and current are also low enough. The crucial part of this route is its idea, purposely preparing “defected” field configurations which can make electrons oscillate at desired frequency during a desired time scale. In the following sections, we will display details of this route.
2. Methods
The physics basis of our method can be summarized as “tailoring” Takeuchi-orbit. As shown in Takeuchi’s theory [
42], a single-electron orbit in an
configuration, where
and
are both DC, can be divided into two classes: elliptical and hyperbola, according to initial velocity of the electron entering into the configuration and values of
[
42,
43]. The elliptical orbit implies that the electron does a time-periodic motion, or oscillation, along a “large” orbit and hence could potentially act as a radiation source. Therefore, for practical purposes, we propose to compress/tailor the geometric size of such a Takeuchi-orbit to a small-enough level by flexible utilization of various
-field components/elements [
40].
-field components/elements we used are very common. For example, Helmholtz coil pair is a common basic component for screening an external magnetic field. We let its two coils be not co-axial on purpose and hence make the magnetic field between them to be of a “defect” or a discontinuity in energy density profile
[
41]. More exactly, the phrase “discontinuity” should be replaced by “slope” because of realistic finite-sized components [
40]. A pair of parallel metal-plates supplies a medium strength
vertical to
, and incident electrons, whose velocities are as slow as possible, are along the
direction and enter into the configuration from a position near the midpoint of the “slope” [
40]. The
-field can be designed to be of a finite thickness, denoted as
in
Figure 1, in the incident direction of electrons bunch by applying high-dielectric-constant insulator medium to “mask” the metal plates (see
Figure 1a). Thus, during the bunch crossing the
, the electron bunch will oscillate along a “spliced” Takeuchi-orbit and the time-duration of doing such an oscillation can be handily controlled by adjusting the
-value if initial incident electron velocity is given.
In such a configuration, the transverse aspect of the single-electron dynamics can be strictly controlled, through values of several parameters, to be an oscillation over the finite-sized “spliced” orbit [
40,
41]. The time cycle, or the frequency, of the oscillation can be controlled by values of those parameters and is, in principle, able to be arbitrary. Thus, relative to a far-field observer on the way of the electrons bunch, there is a mono-color light source. The magnitudes of
and
are not required to be high because one can adjust
D-value to control the time cycle of the tailored Takeuchi-orbit [
40,
41]. To some extent, the tailoring Takeuchi-orbit scheme above-described is alike to well-known synchronization radiation apparatus being stretched to a straight-line by applying a
.
The time-duration of electrons travelling in such a configuration, which depends on the
-value and incident longitudinal velocity of the bunch, determines the time-duration of electrons doing transverse oscillation. Besides adjusting the area of the parallel metal-plates, adjusting the relative position, (represented by
D in
Figure 1) and “gesture” of the plates to the midpoint of the slope, or the
contour, is also effective. As shown in
Figure 1b, infinitely short time-duration is in principle able to be achieved.
However, realistic factors, such as available mechanics processing skill, determine the minimum time-duration to be large than a threshold. Currently available mechanics process skill can warrant
m-level positioning accuracy. For achieving
-level time-duration, nm-level positioning accuracy should be warranted. This can be achieved either through corresponding mechanics device in scanning electron microscope (SEM) or through simple component/element which can minify, according to “lever” principle (see
Figure 2), the space step.
A practical radiation source should have an obvious output power. Therefore, the charge of the bunch, or the number of electrons in the bunch, cannot be too low. The charge is dependent on the density of the parent material which supplies electrons. Considering the fact that the density of solid-state cathode is usually at cmm-level, we should apply more realistic multiple-body dynamics method to this question.
Simulating up to
or more electrons is a challenge. If each electron is described by 6 double-precision format (real-valued digit) floating-point numbers for its 3 space position coordinate components and 3 velocity components,
double-precision format (real-valued digit) floating-point numbers needs up to 6 × 24 TB storage mount (i.e., a double-precision format (real-valued digit) floating-point numbers needs 24 Byte storage mount). Currently popular multiple-body dynamics methods, such as Particle-in-Cell simulation [
44,
45], are mostly in an exhaustive manner to simulate electrons. Because
TB storage mount is beyond the volume of faster semiconductor storage medium (i.e., memory) in computer groups and hence the simulation inevitably involves in time-consuming Input/Output (I/O) operation, which reads/writes data on slower magnetic harddisk (whose access time is usually at
GB per second-level), people have to invent approximation for speeding up the simulation at the cost of scientific reliability. For example, the PIC method usually contains two approximations, one is rigid-macroparticle approximation (RMPA) [
46], and the other is alternative-updation approximation (AUA) [
44,
45]. Even though the number of electrons is so low, such as <10, that the RMPA is not hired, the AUA still threats the scientific reliability of the simulation.
This is because the PIC method is to solve a coupled set of
N relativistic Newton equations (RNEs) and four Maxwell equations (MEs). As pointed out elsewhere [
46], at any time
t, although self-consistent electric field
E can be known from all particles’ information at the time
t through a ME
, self-consistent magnetic field
B cannot. Due to the term
in the ME
, the information of
B at
t has dependence on all particles’ information at
. This implies that each difference version of RNE involves indeed all particles’ coordinates at
, i.e.,
or
Thus, N difference versions form a linear equation set of whose solutions involve inevitably in a matrix. The larger N is, the more hopeless exact solution is. It is well-known that, even for a matrix algebra equation with tens, solving it via computer, whose kernel is an integrated circuit with GHz-level clock frequency, is a formidable task. Likewise, the time-cost of computer solving a matrix algebra equation with a astronomical-figure is also an astronomical-figure in units of seconds.
The AUA avoids time-consuming solving matrix equation in every time step and represents a relation differing from Equation (
1)
which corresponds to a diagonal matrix linking
and
. Clearly, the AUA is a “non-synchronous” trick, and disagrees with the basic principle of calculation mathematics, which requires all coupled difference equations to be solved simultaneously or synchronously. Therefore, our multiple-body dynamics simulation is on a stricter route [
46,
47].
Our simulation is to treat
N RNEs and four MEs in a non-exhaustive manner. No matter how large
N is,
N particles’ RNEs are all Lagrangian versions of the fluid RNE of the combination of a mono-valued field
and a multiple-valued field
[
46,
47]
The phrase “mono-valued” refers to that at each space-time point , the field-value is certain. The phrase “multiple-valued” refers to otherwise cases and hence means that at each space-time point , the field-value is allowed to be over a range with a breadth. Moreover, at each space-time point , the range of allowed values of the -field is from the negative to the positive, and the summation of these allowed values is 0. Due to these properties, it is unnecessary to solve N RNEs and four MEs step-by-step. Instead, N particles’ RNEs leads to a RNE of the mono-valued u-field, and what wait to be solved is this RNE and four MEs.
This result can also be automatically derived from the Vlasov–Maxwell description of plasmas because the solution of microscopic distribution function
f should meet a physical constraint
[
46]. A power series expression
where
,
and
are functions of space-time coordinates
and
, can also lead to a feasible scheme of calculating
f in terms of those
. The constraint
determines an important property of the coefficient function set
. That is,
must be fulfilled for avoiding possible presence of
-region around
-surface within the 6D phase space.
Although an electron bunch has its own self-fields, its global motion is only affected by external fields and such self-fields affect relative motions of parts of the bunch to its “mass-center” which is defined as the zero-value point of the self-fields [
47]. Namely, the self-field is responsible for the variation in the shape of the density profile of a bunch. For a bunch with considerable total charge
, obvious global oscillation of its “mass-center” means a bright light source relative to a far-field observer.
Due to the self-fields of bunch electrons, the radiation spectrum cannot be monocolor and instead of a breadth. However, if the center of a spectrum of a breadth is at the X-ray band or higher and the breadth of the spectrum is finite, such a not-too-much-broad spectrum is of application value. After all, it is very desirable than a broad spectrum starting from an optics-frequency (OF) fundamental frequency component to a component with tens-fold frequency and relative strengths of these components decreases with respect to the harmonics order.
If viewing the behavior of each electron to be governed by the combination of the external fields and the bunch self-fields, its acceleration rate’s Fourier spectrum is not a Dirac function, the spectrum should be if no self-fields exist, but instead has a finite breadth and around the frequency of the external field. Thus, the superposition of those spectrums from different electrons will be still centered/around at the frequency of the external field. Of course, for ensuring EM energy to be mainly at the center frequency, the density of the bunch cannot be overly high. Namely, although high electronic density favors high-intensity at the center frequency, it also favors those at other frequencies. It needs a compromise between high-intensity at the center frequency and the breadth of the spectrum, they are both dependent on, in different manner, the electronic density.
The self-fields cause EM energy to distribute over components differing from the center frequency determined by the external fields. If the quantity “pulse maximum strength” is defined as the maximum of superposition of the square of vector potential components at different frequency components , it might be not, or weakly, affected by the self-fields. Definite answer should be resorted to self-consistent calculation, which will be presented elsewhere. However, the qualitative discussions above have outlined that higher values of electronic density at the zero-value point of the self-fields are favorable to yield higher intensity at the center frequency.
The dependence of the center frequency,
, of the generated pulse on
D, the initial electron position in the configuration, is described by a formula [
41]
where detailed expressions of
,
,
M and
can be found elsewhere [
40,
41]. Note that
is linear dependent on
D.
The polarization state of the generated EM pulse is usually elliptically-polarized because of the shape of the tailored Takeuchi-orbit. Because the shape can be controlled by choosing suitable values of those parameters, the polarization state is in principle able to be controlled. Some authors who are along the route 1 [
48] have achieved polarization control by controlling the ellipticity of the driving laser.