Positron acceleration in a hollow plasma channel up to TeV regime.

Nowadays, human's understanding of the fundamental physics is somehow limited by the energy that our high energy accelerators can afford. Up to 4 TeV protons are realized in the Large Hadron Collider (LHC). Leptons, such as electrons and positrons, however gained energies of about 100 GeV or less. Multi-TeV lepton accelerators are still lacking due to the relatively low acceleration gradient of conventional methods, which may induce unbearable cost. On the other hand, plasmas have shown extraordinary potential in accelerating electrons and ions, providing orders of magnitude higher acceleration fields of 10-100 GV/m. In such context, we propose a plasma-based high-energy lepton accelerator, in which a weakly focusing plasma structure is formed near the beam axis. The structure preserves the emittance of the accelerated beam and produces low radiation losses. Moreover, the structure allows for a considerable decrease of the witness energy spread at the driver depletion stage.

During the past few decades, vital researches on plasma-based accelerators have made great progress in generating high-energy electron beams. Nowadays, electrons with energies beyond 1 GeV are obtained experimentally1234. In order to scale these approaches to the energy frontier of particle physics, i.e., to TeV energies, high energy proton beams, instead of laser pulses, were proposed to drive plasma wakefields and accelerate electrons5, since they can bear much more energy than lasers drivers. A lepton collider for particle physics also requires positron beams of comparable energy and a good beam quality. However, it is simply impossible to apply the proposed blowout-like scenario56 to positron acceleration. The main difficulty is that positrons, for their positive charges, are quickly expelled away from the bubble in the transverse direction789 resulting in no efficient acceleration. There is no such a problem in the linear wakefield regime1011, but the linear regime has lower efficiency12 and smaller accelerating fields, as compared to the nonlinear blowout regime.

To accelerate positrons by the plasma wakefield, the witness beam must be well confined in the accelerating region. We have found out that a simple hollow plasma channel can serve perfectly as the confiner1314. The proposed scheme is sketched in Fig. 1. A cylindrical vacuum channel is embedded in the plasma. When a TeV proton beam travels through, it sucks plenty of electrons from the channel walls and drives a highly nonlinear plasma wave. The drawn away electrons generate a strong charge separation field at the boundary, which is directed to the propagating axis. In the formed structure, the positrons are both accelerated and transversely confined by the proton-driven wakefield. As will be shown later, it makes possible acceleration in a single stage of the length of about 1 kilometer, leading to the TeV energy, high quality positron beams.

Figure 1

Conceptual map of plasma wakefield acceleration in a hollow plasma channel.

The quadrupole magnets are used to focus the proton driver. The accelerating field is shown by the colors.

The proposal shows several remarkable advantages. By optimizing the witness location, the total transverse force acting on the positrons almost vanishes. The beam then stays in a steady state without apparent transverse oscillations. Consequently, a very low transverse emittance can be maintained during the whole interaction length, and the radiation loss is tremendously suppressed.

Results

We demonstrate the overall accelerating process with two-dimensional (2D) simulations. A 2 TeV proton beam, which energy is available in nowadays conventional accelerators, is used as the driver. In order to overcome natural disperse during the 1 km-propagation in plasmas, external quadrupole magnets are used to provide the required focusing as shown in Fig. 1. The quadrupole focusing of the beam is not fully consistent with the axisymmetric geometry of the problem, so we use the axisymmetric analogue of quadrupole fields described in Ref. 6. An initially 1 GeV positron beam is injected behind the driver. The plasma has a constant initial density of n0 = 1015 cm−3 outside the plasma channel (r ≥ r0), where r0 = 0.75 mm is the radius of the vacuum channel. Note that the sharp channel boundaries are not the essential part of the concept, so the channel will also work with diffuse edges as well. The channel itself is wider than state-of-the art channels produced for laser driven wakefield acceleration. These factors make these channels conceivable, though challenging due to the long length required. Detailed simulation parameters can be found in Methods.

The hollow plasma channel enables a good radial confinement of the witness bunch

The accelerating length as long as 1.1 km is realized in a single stage. Figure 2(a–d) shows snapshots of the particle phase space (energy vs distance) at several moments. It is clearly seen that during the propagation though the plasma channel, the proton driver is gradually stretched, losing significant amounts of energy at the tail. Simultaneously, the witness positron beam is catching up with the driver and continuously picking up energies, until it reaches the dephasing region at around 1.1 km. The corresponding energy spectra are inlineed in Figure 2(e–h), respectively. At the end of acceleration, a positron beam with the peak energy of about 1.6 TeV is produced.

Figure 2

Phase space portraits of the driver and the witness bunch (energy W versus co-propagating coordinate ξ = z − ct) (a–d) and corresponding energy spectra (e–h).

The snapshots are taken at propagation distances L = 50, 400, 750, and 1100 m. The positrons and protons are shown in blue and red, respectively.

In Fig. 3(a), the structure of the proton-driven wakefield is shown. The accelerating field increases almost linearly in space along the upslope region (U region) and descends quickly in the downslope region (D region). As the positron beam enters the accelerating region, its peak energy grows proportional to the accelerating length, as is shown in Fig. 4(a). The energy gain saturates soon after the beam enters at the “D region” at around 900 meters.

Figure 3

(a) On-axis acceleration field of the plasma wave. (b) The net charge density distribution and the radial field Er − Bφ at the cross section of witness positrons (black solid line).

Figure 4

(a) The mean energy (red) and energy spread (blue) of the witness positrons as a function of propagation distance. (b) Radiation loss rate (radiation loss/energy gain) for the positron beam as a function of witness energy (dots) and the γ2 scaling (red line). (c) Time dependence of the radial position (blue) and radiation power (red) for a test positron initially at rin = 0.3 mm, the green dashed line shows the critical radius rc. (d) Normalized emittance of the full witness bunch (red) and core particles (blue) in first 100 meters of the acceleration.

A novel compressing effect is seen in the behavior of the energy spread

In Fig. 4(a), the energy spread stays at about 10% all along the “U region” and then remarkably drops in the “D region” to mere 1.5%. Though very narrow compared to “U region”, the “D region” plays an important role in reducing the energy spread of the witness beam. More energetic positrons at the front experience weaker acceleration and vice versa, so that the longitudinal phase volume of the witness drastically shrinks, as presented in Fig. 2(d). We emphasize that this effect is specific to proton drivers, for which the driver-to-witness distance is continually reducing during the interaction. The positive charge of the witness is also important, since negatively charged electrons must reside on the left of the “U region” in Fig. 3(a), and the driver-to-witness distance must increase for the witness to cross the field extreme. Note also that the observed reduction of the energy spread does not rely on beam loading. Thus, the discovered “ballistic” reduction of the energy spread is independent and complimentary (to the beam loading) way of controlling the energy spread.

A weakly focusing structure is formed near the beam axis as shown in Fig. 3(b), which produces low radiation losses and preserves the emittance of the accelerated beam

Generation of the high quality TeV positron beam in 1 km implies not only good channeling and acceleration, but also suppression of the radiation loss. On the other hand, in order to build an electron-positron collider in future, a small beam emittance is required. The weakly focusing plasma structure in the proposed scheme gives a possible solution to these questions as shown in Fig. 4(b–d).

One can see that in the witness location, the transverse field Er − Bφ is almost-zero in the region close to the beam axis. The underlying physics for the formation of such structure is the introduction of hollow channel greatly weakens the “phase mixing” effect15, thus reducing the electron density in the “bubble” (especially near the beam axis) as shown in Fig. 3(b). This is due to the plasma electrons are sucked into the hollow channel with an initially pitch angel (with respect to the beam aixs) larger than a certain value (determined by the channel radius).

A distinctive feature shows up near the axis in Fig. 4(b), where all fields vanish within a critical radius r. The beauty of such a structure is that the radiation of paraxial positrons (r ≤ rc) is almost zero, resulting in a radiation-free region. To see this more clearly, the radiated power of a single test positron is plotted together with its trajectory in Fig. 4(c). The radiation power disappears when the positron enters the near axis region of r ≤ rc ≈ 0.15 mm, just as we predicted. The abaxial positrons, or the ones with large pitch angles, which try to escape from the confining potential well, would be reflected by the strong focusing force, lose transverse velocity by radiation and then fall into the radiation-free region.

The weakly focusing plasma structure also allows for the conservation of the normalized emittance during the acceleration as shown by Fig. 4(d). A much finer grid is used here to follow the emittance growth, so that we only simulate for the first 100 meters in the acceleration due to the limitations of the computing time. One can see that the normalized emittance of the full bunch grows slightly. A few halo particles are responsible for this emittance growth (as discussed in the Discussion). If those particles are excluded, the emittance is almost conserved for the core particles (containing 95% of the total charge). As we use a quasi-static code, there is no numerical Cherenkov radiation inherent to particle-in-cell codes1617, the unphysical emittance growth is fully controlled by grid and time steps.

Discussion

The radiation loss may be important for leptons accelerated to TeV range energies. Previous studies11 indicate the betatron oscillation produces significant energy damping for the electron accelerated in linear plasma wave. Positrons at TeV-level have relativistic factors of the order of 106. Hence the radiation loss in the channel mainly comes from their transverse oscillations.

Assuming the witness bunch propagates towards positive z-direction, the radiation power for a positron can be calculated as18 where r = e2/mc2 is the classic electron radius, e and m are the positron charge and mass, γ = (1 − v2/c2)−1/2 is the relativistic factor, and c is the light velocity, respectively; E and B are the radial electric field and azimuthal magnetic field in the channel. The term fq is the time-average force of the quadrupole magnets given by6 where S is the magnetic field gradient of the quadrupoles, Lq is the space period, and Ww is the energy of witness bunch. The term fr in Eq. (1) is the radial component of the radiation reaction (RR) force on the positron, which, in the relativistic case, is19 with EZ the accelerating field, θ = v/v the angle between the positron velocity and the axis, and β the transverse velocity of the positron in units of c. Since θ ≪ 1, the items related to Er − Bφ are the main source of the radiation loss. The radial Lorentz force Er − Bφ is shown in Fig. 3(b). Inside the channel (r ≤ r0), the focusing field is created by plasma electrons pulled out of the boundary.

One can numerically calculate the radiation loss rate from Eqs. (1–3). Let G = eEzvz be the energy gain of the witness positron per unit time. We define the radiation loss rate by R = P/G. A bunch of test positrons is tracked in the electro-magnetic fields taken from simulations. As predicted by Eq. (1), the quadratic dependence of the radiation loss rate on the witness energy holds up to multi-TeV regime [Fig. 4(b)]. In our case, R is only 3.2 × 10−7 for 1.6 TeV energy. The radiation loss is extremely low so that more precise calculations are unnecessary.

Since the radial force below the critical radius is almost zero, the positrons accelerated there are free from emittance growth2021. Nonetheless, there would always be some positrons escaping from this area and being reflected by the boundary field. As the threshold r varies slightly with time (acceleration distance), such reflection is usually not a “mirror reflection”, resulting in increase of the beam emittance. It is thus desirable to avoid touching the critical radius by the witness body. Let us estimate the witness emittance required for that. Considering a test positron initially located at radius σwr, having pitch angle αi and the relativistic factor γi. This positron can be thought as a marginal one for the positron beam of the normalized emittance εcn ≈ γiαiσwr. If uniformly accelerated by the field Ez to the relativistic factor γf, this positron will change its radial position to Assuming rf = rc, we find the threshold witness emittance If the emittance is below this value, the most part of the witness propagates in vacuum and sees no particles or transversely deflecting fields on its path.

For the present parameters, Eq. (5) gives the threshold of approximately 2 mm-mrad. We then chose the initial witness emittance of 1 mm-mrad in our simulation. The normalized emittance almost conserves in the simulation as inlineed in Fig. 4(d), which confirms our discussion above.

In conclusion, we have suggested a PWFA scheme capable of accelerating positrons to energies beyond 1 TeV. The novel feature of this scheme is a combination of the proton driver and the vacuum channel created inside a uniform plasma. An initially 1 GeV positron beam is not only accelerated efficiently by the proton-driven wakefield, but also confined radially in the almost rectangular potential well. In this scheme, the radiation damping is negligibly small, and the normalized emittance is almost preserved for the core particles. The relatively low velocity of the proton driver is turned to advantage and used to reduce the final energy spread of accelerated positrons. In the illustrative simulation, a quasi-monoenergetic positron beam is obtained in a single acceleration stage, with the final energy of 1.6 TeV and the energy spread of only 1.5%.

Methods

The simulations are performed with computationally efficient code LCODE22 capable of simulating long distance beam propagation. The simulation is carried out in the axisymmetric geometry. In the code, the simulation window moves with the light velocity. The beams and plasma are modeled by fully relativistic macro-particles. Quasi-static approximation is used for calculating the plasma response so that the radiation reaction is not included in this code and was estimated separately. Two different simulation grid sizes are used in this work. The finer one of size 0.01c/ωp ≈ 1.68 μm is used to study the detailed shape of the plasma wave and the emittance conservation (ωp is the plasma frequency of the surrounding plasma), and the coarse one of the size 0.05c/ωp ≈ 8.4 μm is used to simulate acceleration process over the long distance, The time steps for these two cases are and , respectively.

The detailed simulation parameters are listed in Table 1. The parameters of the driver are similar to those used in Refs. 5 and 13. The driver energy is larger (2 TeV) so that it can accelerate positrons beyond 1 TeV, and a typical 1 GeV positron beam serves as the witness. We have scanned over a range of plasma densities (near the matching condition 59) and channel radii and found out that n0 = 1015 cm−3 and r0 = 0.75 mm provides the best performance.

Table 1

Parameters for the simulation

Parameters	Symbols	Values
Number of protons in the drive beam	Np	1011
Initial energy of protons	WP	2 TeV
Initial longitudinal size of the driver	σz	100 μm
Initial spot size of the driver	σr	430 μm
Initial longitudinal momentum spread of the driver	ΔPz/Pz	0.1
Normalized emittance of the driver	εn	1 mm-mrad
Initial witness energy	Ww	1 GeV
Total witness charge	Qw	800 pC
Initial bunch length of the witness beam	σwz	25 μm
Initial spot size of the witness beam	σwr	50 μm
Initial normalized emittance of the witness beam	εwn	1 mm-mrad
Plasma density	n0	1015 cm−3
Plasma channel radius	r0	0.75 mm
Magnetic field gradient	S	0.5 T/mm
Space period of the quadrupoles	Lq	3 m

Author Contributions

L.-Q.Y., B.-F.S. and L.-L.J. contributed to all aspects of this work; K.L. and A.S. contributed to code development; X.-M.Z., W.-P.W., J.-C.X. and K.L. provided inspiring ideas and helped L.-Q.Y. to write the paper; Y.S. and L.-G.Z. helped to analyze the radiation loss in the acceleration process; Z.-Z.X. gave some useful comments and suggestions to this work.