Simulation of the interaction of macro-particles with the LHC proton beam

1. Simulation of the interaction of macro-particles with the LHC proton beam Zhao Yang, EPFL 14.6.2011

2. Equation of motion Introduction The purpose of simulation Study the influence of magnetic field/mass of particle/number of protons in the beam on the particle trajectory/beam loss rate/loss duration/total lost protons. -g Parameters

3. Aluminum Macro-particle Influence of magnetic field on particle trajectory. Np=0.9*1011*25 A equals to 1012~1015, B equals to 0 T, 8.33 T and 80 T (Blue, red and green) A=1012 A=1013 A=1014 A=1015 • The X-Y trajectory gets close to the Y-axis as the mass of particle decreases and the magnetic field intensity is increased. The strong magnetic field may result in a cyclotron oscillation of the particle. The period of this oscillation is reduced for increasing magnetic field and smaller particle mass.

4. Aluminous Macro-particle Influence of magnetic field on beam loss rate. Np=0.9*1011*25 A equals to 1012~1015, B equals to 0 T, 8.33 T and 80 T (Blue, red and green) 3.306 3.015 4.731 4.525 • The maximum of loss rate of the 1st crossing is independent of magnetic field; it increases with increasing mass A. • The maximum of loss rate of the 2nd crossing increases with increasing magnetic field strength. It decreases with increasing mass at normal magnetic field strength. A=1012 A=1013 -38.337 -63.448 -40.958 -66.453 7.470 6.123 5.582 -48.090 A=1014 A=1015 -123.319 -323.12 -126.780 -326.67 • The reason for a strange behavior of the loss rate in the case A=1014 B=80 T is unclear.

5. Aluminous Macro-particle Particle trajectory of A=1012, B=8.33 T Np=1.15*1011*2808 0.064946 s 0.060703 s 0.13530 s 0.14275 s At t=0.060703 s, the macro-particle gets close to the beam and there comes the 1st crossing, it is repelled upwards and hits the wall at t=0.064946 s. Then it falls down and approaches to the beam. During this process it hits the wall a 2nd time at t=0.13530 s. Finally it falls down at 0.14275 s. The time between two crossing is 0.082 s

6. Aluminous Macro-particle Plot the figures of beam loss rate for a varying total number of protons in the beam. We determine the maximum of each beam loss rate curve, and fit the curve around this maximum. Then we calculate the total number of lost protons. Np=1.15*1011*400 Np=1.15*1011*2808 Np=1.15*1011*1600 The maximum of each beam loss rate curve is fitted by the function log10y=a+bx+cx2 Original plots Fitting curves

7. Aluminous Macro-particle The number of lost protons calculated with sigma=0.0003 m and B=8.33T

8. Aluminous Macro-particle Plots of lost protons versus Np with different particle masses A=1013 A=1012 The solid lines are the original plots, the dotted line are the fitting ones. The number of lost protons is fitted by the function y=1/(a+bx+cx2) when the mass is relatively small. A=1014

9. Aluminous Macro-particle B=4.17 T is the present value of magnetic field in LHC at 3.5 TeV beam energy. We use sigma=0.0003 m, B=4.17 T and Np=1.15*1011*2808 to simulate the beam loss rate and calculate the total number of lost protons. The results are exactly the same as the results with B=8.33 T the magnetic field in LHC will not affect the total number of lost protons. We increase the value of transverse rms beam size by a factor 5, to sigma=0.0015 m, and keep the other parameters constant. The increase of the transverse beam size reduces the total number of lost protons.

10. Copper Macro-particle We change the material of macro-particle. Compared with aluminum, copper has larger density, cross section and atomic mass number. sigma=0.0003 m, x0=0.0001 m, Np=0.9*1011*25 The X-Y trajectories of Aluminum are closer to the Y-axis than the X-Y trajectories of Copper. The maximum of beam loss rate with aluminum macro-particle is larger than that with a copper macro-particle. A=1012 A=1013 A=1014 A=1013 A=1014 A=1012

11. Copper Macro-particle We plot the curves of beam loss rate versus time, then calculate the loss duration and the total number of lost protons for the 1st crossings x0=0.0001 m, B=8.33 T, A=1012, 1014 and 1016 Np=1.15*1011*2808 Np=1.15*1011*1600 Np=1.15*1011*400 The loss duration is almost independent of the material of the macro-particle. The total number of lost protons for a copper particle is smaller than that for an aluminum particle. Loss duration obtained for varying mass A and total proton intensity Np in unit of second (B=8.33 T). Total lost protons obtained for varying value of mass A and total protons intensity Np (B=8.33 T).

12. Final conclusions • The time separation between first and second crossing is consistent with some beam observations of multiple successive events. However, the losses at the second crossing are always much lower than for the first crossing, which is different from actual observations. • The computed peak losses should be compared with the quench threshold which corresponds to about 107 protons / second. • Effect of magnetic force on the macro-particle motion is weak and can be neglected, even for a field of 8.33 T. • The loss duration and the number of lost protons decrease with higher total beam intensity; the losses roughly in inverse proportion. • Increasing the beam size by a factor of 5 reduces the total proton loss by about a factor of 3. This might be part of the explanation why events have not been important at LHC injection. Future work may extend this discussion to some other materials and to other macro-particle shapes, e.g. to needle-like objects.

13. The end Thank you