Hostname: page-component-76d6cb85b7-mgxrv Total loading time: 0 Render date: 2026-07-16T21:09:23.944Z Has data issue: false hasContentIssue false
  • English
  • Français

Manipulating the Laser-Driven Proton Bunch with Plasma Wakefield

Published online by Cambridge University Press:  01 January 2024

Chao Jin Open the ORCID record for Chao Jin [Opens in a new window] ,
Xiao-ying Zhao Open the ORCID record for Xiao-ying Zhao [Opens in a new window] ,
Han-jie Cai Open the ORCID record for Han-jie Cai [Opens in a new window] ,
Xin Qi Open the ORCID record for Xin Qi [Opens in a new window] ,
Zhi-jun Wang Open the ORCID record for Zhi-jun Wang [Opens in a new window]  and
Yuan He Open the ORCID record for Yuan He [Opens in a new window]
Chao Jin
Affiliation:
Institute of Modern Physics, Chinese Academy of Science, Lanzhou 730000, China School of Nuclear Science and Technology, University of Chinese Academy of Science, Beijing 100049, China
Xiao-ying Zhao
Affiliation:
Institute of Modern Physics, Chinese Academy of Science, Lanzhou 730000, China School of Nuclear Science and Technology, University of Chinese Academy of Science, Beijing 100049, China
Han-jie Cai
Affiliation:
Institute of Modern Physics, Chinese Academy of Science, Lanzhou 730000, China School of Nuclear Science and Technology, University of Chinese Academy of Science, Beijing 100049, China
Xin Qi*
Affiliation:
Institute of Modern Physics, Chinese Academy of Science, Lanzhou 730000, China School of Nuclear Science and Technology, University of Chinese Academy of Science, Beijing 100049, China
Zhi-jun Wang
Affiliation:
Institute of Modern Physics, Chinese Academy of Science, Lanzhou 730000, China School of Nuclear Science and Technology, University of Chinese Academy of Science, Beijing 100049, China
Yuan He
Affiliation:
Institute of Modern Physics, Chinese Academy of Science, Lanzhou 730000, China School of Nuclear Science and Technology, University of Chinese Academy of Science, Beijing 100049, China
*
Correspondence should be addressed to Xin Qi; qixin2002@impcas.ac.cn

Abstract

With the advantages of short duration and extreme brightness, laser proton accelerators (LPAs) show great potential in many fields for industrial, medical, and research applications. However, the quality of current laser-driven proton beams, such as the broad energy spread and large divergence angle, is still a challenge. We use numerical simulations to study the propagation of such proton bunches in the plasma. Results show the bunch will excite the wakefield and modulate itself. Although a small number of particles at the head of the bunch cannot be manipulated by the wakefield, the total energy spread is reduced. Moreover, while reducing the longitudinal energy spread, the wakefield will also pinch the beam in the transverse direction. The space charge effect of the bunch is completely offset by the wakefield, and the transverse momentum of the bunch decreases as the bunch transports in the plasma. For laser-driven ion beams, our study provides a novel idea about the optimization of these beams.

Information

Type
Research Article
Information
Laser and Particle Beams , Volume 2022 , 2022 , e2
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright
Copyright © 2022 Chao Jin et al.

1. Introduction

During the last decade, plasma-based particle accelerators driven by high-intensity ultrashort laser pulses [Reference Ding, Schumacher, Jahn, Blazevic and Roth1–Reference Higginson, Gray and King3] or particle beams have shown great promise, primarily because of the extremely large accelerating electric fields that they can support, about a thousand times greater than conventional accelerators, enabling the realization of laboratory-scale applications ranging from high-energy physics to ultrabright light sources.

For proton accelerations, target normal sheath acceleration (TNSA) [Reference Hegelich, Albright and Cobble4, Reference Fuchs, Antici and D’Humières5], radiation pressure acceleration (RPA) [Reference Esirkepov, Borghesi, Bulanov, Mourou and Tajima6, Reference Macchi, Cattani, Liseykina and Cornolti7] and breakout afterburner (BOA) [Reference Yin, Albright, Bowers, Jung, Fernández and Hegelich8] are the most widely employed mechanisms. Despite the relatively low energy transfer efficiency, TNSA is considered the most robust and stable mechanism, where the ions (especially protons) from the surface contamination layer are accelerated by the charge separation sheath field. Since several experiments have demonstrated large accelerations, the resulting beam quality is still far from state-of-the-art conventional accelerators: due to the rapid diffusion of the electron layer, Coulomb explosion, and multiple instabilities, TNSA ion beams are characterized by an exponentially decaying energy spectrum and a large divergence angle [Reference Linz and Alonso9], Thus, one of the main goals is to control the energy divergence and shape the beam.

Various beamlines have been proposed at several institutes, for instance, the light beam line at GSI Helmholtz Center [Reference Jahn, Schumacher and Brabetz10, Reference Busold, Schumacher and Brabetz11], the ELIMED beamline installed in Prague [Reference Romano, Schillaci and Cirrone12] and the CLAPA platform at Peking University [Reference Zhu, Wu and Zhu13, Reference Zhu, Wu and Liao14]. These beamlines are designed to realize the propagation of high current and dense ion beams with low energy spread and with reliability, availability, maintainability, and inspectability.

Compared with these beamlines, which are composed of conventional accelerator elements such as permanent magnet quadrupole lenses, solenoid magnets, and laser-triggered microlenses, recent experiments show that when electron beams propagate in a plasma element, the plasma wakefield will tune the longitudinal phase space of the electron beams and reduce the energy spread. For example, FLASHForward plasma-accelerator facility at DESY [Reference Pompili, Alesini and Anania15] showed that in beam-driven plasma acceleration, if an electron bunch acts as a driver, while an ultrashort witness bunch with positive energy chirp follows behind, the longitudinal phase space of the witness bunch will be rotated during acceleration, resulting in an ultralow energy spread that is even lower than the spread at the plasma entrance. Another experiment at the SPARC LAB test facility [Reference Lindstrøm, Garland and Schröder16] demonstrated that high efficiency and low energy spread can be achieved simultaneously by strong beam loading of plasma wakefields when accelerating bunches with carefully tailored current profiles. In both experiments above, there are two bunches: the driven bunch creates wakefields, while the energy spread of the witness bunch is reduced in these wakefields. Another experiment shows that even for a single electron bunch with negative chirp passing through plasmas, its wakefields will manipulate the longitudinal phase space (LPS) along the beam itself, and reduce the total energy spread [Reference Wu, Hua and Zhou17, Reference Shpakov, Anania and Bellaveglia18].

These experiments mainly concern about electron bunches. For proton accelerations, the proton beams have the typical exponentially decaying spectrum and positive energy chirp. In this paper, we employ the two-dimensional particle-in-cell (PIC) algorithm [Reference Zhang, Zhao, Qi, Duan and Yang19, Reference Kumar, Pukhov and Lotov20] to study the propagation of such beams in plasmas. The effect of plasma wakefield on longitudinal energy spread, as well as the transverse compression of a single proton bunch are investigated. The main contents of this paper are as follows: Section 2 introduces the physical model used in our simulations. In Section 3, the simulation results are presented and discussed. Finally, conclusions are made in Section 4.

2. Materials and Methods

The physical model is shown in Figure 1. The proton bunch propagates from left to right in the hydrogen plasma channel. Assuming the plasma is collision less, nonpolarized, and stable.

Figure 1:

Two-dimensional simulation model: proton bunch propagates in the plasma.

Normally, the proton beam produced by LPA has the characteristics of exponential decay energy distribution. The particle number spectrum dN/dE follows the formula: d N / d E = N 0 / 2 E k B T e − 2 E / k B T [Reference Fuchs, Antici and D’Humières5]. The bunch needs to be selected in advance by the sector magnet and a slit: when the beam with broad energy spread passing the sector magnet (x (horizontal) direction), protons with different energies are dispersed along the x axis. Then, a slit is used to choose the bunch with expected energy. In Figure 2(a), the blue line is the energy spectrum of the beam generated by the LPA and the red line is the one we chose to use as the initial bunch for simulations. The initial bunch has an exponentially decaying energy spectrum, and the energy of the particles in the bunch is 44.75 MeV to 45.25 MeV. In addition, since the duration of the bunch generated by the laser accelerator is usually very short (a few ps, and the longitudinal positions of protons with different energies are basically the same), the longitudinal phase space of the bunch after propagating a given distance depends mainly on the energy distribution. The longitudinal positions of 44.75 MeV protons and 45.25 MeV protons are basically the same when they are generated by the laser accelerator. After propagating 60.3 ns (propagation of about 5.4 m), they will be separated by 28 mm (the length of the beam used in the simulation) due to different velocities. The relationship between energy and velocity is not linear, but when the proton is between 44.75 MeV and 45.25 MeV, the relationship between energy and velocity is very close to linear, so the longitudinal phase space of the initial bunch is shown in Figure 2(b). Due to the space charge effect, disturbance is added to the initial bunch. The initial bunch used in the simulation is shown in Figure 2(c).

Figure 2:

(a) The blue line is the energy spectrum of the bunch generated by the laser accelerator. The red line is the energy spectrum of the bunch after passing through the sector magnet and the slit, using this bunch as the initial bunch for the simulation. (b) Distribution of longitudinal positions of protons with different energies produced by laser accelerator after 60.3 ns of transmission. (c) Longitudinal phase space of initial bunch.

The wavelength of the wakefield is λ w = 2 π v b / ω e = 31.5 mm , where ω e = N e 0 e 2 / ε 0 m e , We choose the bunch length LB = 28mm, which is slightly smaller than the wakefield wavelength, and the bunch RMS radius = 0.36 mm. The bunch is initially a Gaussian distribution (Figure 3(a)). The transverse phase space of the initial bunch is shown in Figure 3(b) and y′ size (Rms) = 0.94 mrad. A laser accelerator produces a bunch with particle numbers between 108 and 1012 per MeV (the number of particles decreases exponentially with energy growth) [Reference Wagner, Deppert and Brabetz21, Reference Harres, Alber and Tauschwitz22]. To make the modulation of the beam in the plasma more obvious, we chose a higher current intensity than the experiment. The initial beam used in the simulation contains 1.92 × 10 12 protons, and the bunch density is N b 0 = 9.52 × 10 16 m − 3 . We select the plasma density Ne 0 according to the bunch density, N e 0 = 1 × 10 17 m − 3 . The plasma density distribution is shown in Figure 3(c), plasma temperature T e 0 = T i 0 = 1 e V .

Figure 3:

(a) Distribution of the initial bunch, (b) transverse phase space of the initial bunch, and (c) the distribution of plasma density in the longitudinal direction.

The software VORPAL [Reference Nieter and Cary23] code is used to simulate the propagation of proton bunch in the plasma. All simulations are performed by a 2D3V PIC algorithm. The length of the simulation area Lz is 0.338 m, while the width Ly 0.024 m. 7200 grids in the z direction and 512 grids in the y direction are used. Moreover, open boundary condition is employed in both z and y direction. Spatial step d z = d y = 2 ∗ λ e = 4.7 ∗ 10 − 5 m , where λe is Debye length of the plasma. The time step d t = 8.1475 × 10 − 14 s .

3. Results and Discussion

Figure 4 is the time evolution of the LPS as the proton bunch passes through the plasma. Figures 4(a)–4(d) are results when the bunch propagates in the plasma at distance of 0, 7.2, 14.4, and 21.6 cm, respectively. The solid line is the longitudinal electric field on the axis at the bunch position. Figures 4(e)–4(h) is the corresponding electron distribution in plasma. When passing through the plasma, the protons at the head of the bunch will disturb the electron distribution in the plasma and create a wakefield. This wakefield will affect the subsequent protons in the bunch. For a bunch with the positive energy chirper, the wakefield reduces the energy spread due to that the first half of the protons in the bunch lose their energy and create the wakefield, while the latter half of the protons gain energy from it. Thus, the total energy spread of the bunch is reduced. As shown in Figures 4(c) and 4(d), protons at the tail of the bunch are in a longitudinal accelerating electric field, while protons are in front under the effect of the deceleration field. It is noted that although the wakefield can modulate the energy spread, there is still a small accelerating electric field at the very head of the bunch. This accelerating electric field is caused by the space charge effect of the bunch: the wakefield only exists behind the bunch; thus, protons at the very head of the bunch cannot be affected. Therefore, these proton’s energies will continue to diverge by Coulomb’s repulsive force.

Figure 4:

(a–d) Beam longitudinal phase space (dot) and longitudinal plasma wakefield Ex (blue line). (e–h) The corresponding electron distribution in plasma. (a) Initial bunch energy distribution. (e) Initial electron distribution in plasma. (b, f) Propagation distance of bunch s = 7.2 cm. (c, g) Propagation distance of bunch s = 14.4 cm. (d, h) Propagation distance of bunch s = 21.6 cm.

Figure 5 is the time evolution of the energy spectrum of the bunch. Figure 5(a) is the initial energy spectrum of the bunch with exponentially decaying energy distribution (The exponential decay of the particle number is not obvious due to the small energy spread). The initial energy spread is 0.5 MeV (44.75 MeV to 45.25 MeV). As the bunch passes through the plasma, although few protons at the very head of the bunch cannot be modulated, the energy of the rest protons begins to concentrate. When the bunch passes through the plasma, more than 50% of the proton energies are in the range from 44.85 MeV to 44.95 MeV. That is to say, the energy spread of half protons is concentrated to 0.2% by the plasma modulation, as shown in Figure 5(d).

Figure 5:

Energy spectrum of the bunch. (a) Initial bunch energy spectrum. Propagation distance s = 0 cm, (b) propagation distance s = 7.2 cm, (c) propagation distance s = 14.4 cm, and (d) propagation distance s = 21.6 cm.

The control of the transverse divergence of the bunch is another challenge. Figures 6(a) and 6(b) are the distribution and the transverse phase space of the bunch when the bunch passes through the plasma, while Figures 6(c) and 6(d) are the comparison of the bunch passing through a vacuum pipe at the same distance. Compared with these in a vacuum, the transverse distributions of the bunch in plasma are reduced for both position and momentum. In particular, the momentum distribution decreases significantly. This is due to that the protons are heavy. When passing through a plasma channel of only 21.6 cm, the wakefield significantly changes its momentum, but it has not significantly affected its position distribution yet.

Figure 6:

(a, b) The distribution and the transverse phase space of the bunch after the bunch transmits 21.6 cm in the plasma. (c, d) The distribution and the transverse phase space of the bunch after the bunch transmits 21.6 cm in the vacuum.

Another phenomenon in Figures 6(a) and 6(b) is that the transverse compression at the rear of the bunch is more obvious. To further investigate it, the corresponding longitudinal electric field E z , transverse electric field E y as well as the magnetic field B x are shown in Figures 7(a)–7(c), separately. Figure 7(a) is the distribution of the longitudinal electric field. As discussed above, it is periodic and goes with the bunch. The head of the bunch is in the deceleration field, while the rest is in the acceleration field. Figure 7(b) is the distribution of the transverse electric Field E y . The transverse defocuses and the focus regions alternately appear. The head of the bunch is in the defocus region while the rear of the bunch is in the focus region. Moreover, the total bunch is surrounded by the magnetic field, as shown in Figure 7(c). For protons in the first half of the bunch, although they are in the transversely divergent electric field, they are pinched by the magnetic field at the same time. Thus, the plasma still limits its transverse divergence. The rear protons are in the transverse focusing electric field and pinched magnetic field at the same time, so the transverse compression effect is more obvious.

Figure 7:

The electric and magnetic field (E x , E y , and B z ) in the wakefield region induced by the proton bunch traveling through the plasmas and propagation distance s = 21.6 cm. The blue ellipse in the figure represents the bunch. (a) Longitudinal electric field distribution in plasma (unit is MV/m), (b) transverse electric field distribution in plasma (unit is MV/m), and (c) magnetic field in the direction which perpendicular to the propagation plane (unit is mT).

The envelop curves of the transverse position and momentum versus the longitudinal position of the bunch are shown in Figure 8. The blue lines represent the result of bunch propagation in plasma, while the red lines are the result in a vacuum. Plasma can reduce the transverse divergence of the bunch. It is noted that the bunch stays in the plasma for only 2.4 ns (propagation distance 21.6 cm), therefore, the plasma cannot compress the transverse size of the bunch in such a short time. However, the wakefield of the plasma has significantly changed the transverse momentum of the bunch. As shown in Figure 8(b), the final transverse momentum of the bunch is lower than that before entering the plasma, which means that the wakefield not only completely neutralizes Coulomb’s repulsion force but also starts to pinch the bunch.

Figure 8:

Variation of the transverse size of the bunch when propagating in plasma.

4. Conclusions

In this paper, we use PIC simulations to study the propagation of a laser-driven proton bunch in the plasma. Our results show that the laser-driven proton beams, of which are characterized by large energy spread and divergence angles, can be modulated by the plasma wakefield. The bunch will excite the wakefield and modulate itself. Although a small number of particles at the head of the bunch cannot be modulated by the wakefield, the total energy spread is concentrated. Especially, for a proton bunch with energies between 44.75 MeV and 45.25 MeV, the energy of more than half of the protons is concentrated between 44.85 and 44.95 by plasma modulation. Moreover, while reducing the longitudinal energy spread, the wakefield will also pinch the beam in the transverse direction. The space charge effect of the bunch is completely offset by the wakefield, and the transverse momentum of the bunch decreases as the bunch transports in the plasma. For laser-driven proton beams, the bunches are characterized by large energy spread and divergence angles, which is difficult to optimize with traditional accelerator elements. The paper provides a novel idea for solving this problem.

Data Availability

The source code and data that support the findings of this study are available upon request from the corresponding authors.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant nos. 11775282 and 11805253), the Youth Innovation Promotion Association CAS (Grant no. 2018452), and Large Research Infrastructures “China Initiative Accelerator Driven System” (Grant no. 2017-000052-75-01-000590).

References

Ding, J., Schumacher, D., Jahn, D., Blazevic, A., and Roth, M., “Simulation studies on generation, handling and transport of laser-accelerated carbon ions,” Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, vol. 909, 2018.CrossRefGoogle Scholar
Kim, I. J., Pae, K. H., Choi, I. W. et al., “Radiation pressure acceleration of protons to 93 MeV with circularly polarized petawatt laser pulses,” Physics of Plasmas, vol. 23, no. 7, Article ID 070701, 2016.CrossRefGoogle Scholar
Higginson, A., Gray, R. J., King, M. et al., “Near-100 MeV protons via a laser-driven transparency-enhanced hybrid acceleration scheme,” Nature Communications, vol. 9, no. 1, p. 724, 2018.CrossRefGoogle Scholar
Hegelich, B. M., Albright, B. J., Cobble, J. et al., “Laser acceleration of quasi-monoenergetic MeV ion beams,” Nature, vol. 439, no. 7075, pp. 441–444, 2006.CrossRefGoogle ScholarPubMed
Fuchs, J., Antici, P., D’Humières, E. et al., “Laser-driven proton scaling laws and new paths towards energy increase,” Nature Physics, vol. 2, no. 1, pp. 48–54, 2006.CrossRefGoogle Scholar
Esirkepov, T., Borghesi, M., Bulanov, S. V., Mourou, G., and Tajima, T., “Highly efficient relativistic-ion generation in the laser-piston regime,” Physical Review Letters, vol. 92, no. 17, Article ID 175003, 2004.CrossRefGoogle ScholarPubMed
Macchi, A., Cattani, F., Liseykina, T. V., and Cornolti, F., “Laser acceleration of ion bunches at the front surface of overdense plasmas,” Physical Review Letters, vol. 94, no. 16, Article ID 165003, 2005.CrossRefGoogle ScholarPubMed
Yin, L., Albright, B. J., Bowers, K. J., Jung, D., Fernández, J. C., and Hegelich, B. M., “Three-dimensional dynamics of breakout afterburner ion acceleration using high-contrast short-pulse laser and nanoscale targets,” Physical Review Letters, vol. 107, no. 4, Article ID 045003, 2011.CrossRefGoogle ScholarPubMed
Linz, U. and Alonso, J., “Laser-driven ion accelerators for tumor therapy revisited,” Physical Review Accelerators and Beams, vol. 19, no. 12, p. 124802, 2016.CrossRefGoogle Scholar
Jahn, D., Schumacher, D., Brabetz, C. et al., “First application studies at the laser-driven LIGHT beamline: improving proton beam homogeneity and imaging of a solid target,” Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, vol. 909, pp. 173–176, 2018.CrossRefGoogle Scholar
Busold, S., Schumacher, D., Brabetz, C. et al., “Towards highest peak intensities for ultra-short MeV-range ion bunches,” Scientific Reports, vol. 5, no. 1, Article ID 12459, 2015.CrossRefGoogle ScholarPubMed
Romano, F., Schillaci, F., Cirrone, G. et al., “The ELIMED transport and dosimetry beamline for laser-driven ion beams,” Nuclear Instruments and Methods in Physics Research Section A Accelerators Spectrometers Detectors and Associated Equipment, vol. 829, pp. 153–158, 2016.CrossRefGoogle Scholar
Zhu, J., Wu, M., Zhu, K. et al., “Demonstration of tailored energy deposition in a laser proton accelerator,” Physical Review Accelerators and Beams, vol. 23, no. 12, Article ID 121304, 2020.CrossRefGoogle Scholar
Zhu, J. G., Wu, M., Liao, Q. et al., “Experimental demonstration of a laser proton accelerator with accurate beam control through image-relaying transport,” Physical Review Accelerators and Beams, vol. 22, no. 6, Article ID 061302, 2019.CrossRefGoogle Scholar
Pompili, R., Alesini, D., Anania, M. P. et al., “Energy spread minimization in a beam-driven plasma wakefield accelerator,” Nature Physics, vol. 17, no. 4, pp. 499–503, 2021.CrossRefGoogle Scholar
Lindstrøm, C., Garland, J. M., Schröder, S. et al., “Energy-spread preservation and high efficiency in a plasma-wakefield accelerator,” Physical Review Letters, vol. 126, no. 1, p. 014801, 2021.CrossRefGoogle Scholar
Wu, Y., Hua, J., Zhou, Z. et al., “Phase space dynamics of a plasma wakefield dechirper for energy spread reduction,” Physical review letters, vol. 122, 2019.CrossRefGoogle ScholarPubMed
Shpakov, V., Anania, M. P., Bellaveglia, M. et al., “Longitudinal phase-space manipulation with beam-driven plasma wakefields,” Physical Review Letters, vol. 122, no. 11, p. 114801, 2019.CrossRefGoogle ScholarPubMed
Zhang, L. Y., Zhao, X. Y., Qi, X., Duan, W. S., and Yang, L., “Particle-in-cell simulations of the long proton beam focusing in background plasmas,” Laser and Particle Beams, vol. 34, no. 2, pp. 299–305, 2016.CrossRefGoogle Scholar
Kumar, N., Pukhov, A., and Lotov, K., “Self-modulation instability of a long proton bunch in plasmas,” Physical Review Letters, vol. 104, no. 25, Article ID 255003, 2010.CrossRefGoogle ScholarPubMed
Wagner, F., Deppert, O., Brabetz, C. et al., “Maximum proton energy above 85 MeV from the relativistic interaction of laser pulses with micrometer ThickCH2Targets,” Physical Review Letters, vol. 116, no. 20, Article ID 205002, 2016.CrossRefGoogle Scholar
Harres, K., Alber, I., Tauschwitz, A. et al., “Beam collimation and transport of quasineutral laser-accelerated protons by a solenoid field,” Physics of Plasmas, vol. 17, no. 2, Article ID 023107, 2010.CrossRefGoogle Scholar
Nieter, C. and Cary, J. R., “VORPAL: a versatile plasma simulation code,” Journal of Computational Physics, vol. 196, no. 2, pp. 448–473, 2004.CrossRefGoogle Scholar
Figure 0

Figure 1: Two-dimensional simulation model: proton bunch propagates in the plasma.

Figure 1

Figure 2: (a) The blue line is the energy spectrum of the bunch generated by the laser accelerator. The red line is the energy spectrum of the bunch after passing through the sector magnet and the slit, using this bunch as the initial bunch for the simulation. (b) Distribution of longitudinal positions of protons with different energies produced by laser accelerator after 60.3 ns of transmission. (c) Longitudinal phase space of initial bunch.

Figure 2

Figure 3: (a) Distribution of the initial bunch, (b) transverse phase space of the initial bunch, and (c) the distribution of plasma density in the longitudinal direction.

Figure 3

Figure 4: (a–d) Beam longitudinal phase space (dot) and longitudinal plasma wakefield Ex (blue line). (e–h) The corresponding electron distribution in plasma. (a) Initial bunch energy distribution. (e) Initial electron distribution in plasma. (b, f) Propagation distance of bunch s = 7.2 cm. (c, g) Propagation distance of bunch s = 14.4 cm. (d, h) Propagation distance of bunch s = 21.6 cm.

Figure 4

Figure 5: Energy spectrum of the bunch. (a) Initial bunch energy spectrum. Propagation distance s = 0 cm, (b) propagation distance s = 7.2 cm, (c) propagation distance s = 14.4 cm, and (d) propagation distance s = 21.6 cm.

Figure 5

Figure 6: (a, b) The distribution and the transverse phase space of the bunch after the bunch transmits 21.6 cm in the plasma. (c, d) The distribution and the transverse phase space of the bunch after the bunch transmits 21.6 cm in the vacuum.

Figure 6

Figure 7: The electric and magnetic field (Ex, Ey, and Bz) in the wakefield region induced by the proton bunch traveling through the plasmas and propagation distance s = 21.6 cm. The blue ellipse in the figure represents the bunch. (a) Longitudinal electric field distribution in plasma (unit is MV/m), (b) transverse electric field distribution in plasma (unit is MV/m), and (c) magnetic field in the direction which perpendicular to the propagation plane (unit is mT).

Figure 7

Figure 8: Variation of the transverse size of the bunch when propagating in plasma.