Radial Expansion of a Low Energy Positron Beam ... - Science Direct

0 downloads 0 Views 845KB Size Report
The magnetic force on each particle is given by. (2). Here, is the Lorentz force constant ( in SI units), is the magnitude of a uniform magnetic field in the direction ...
Available online at www.sciencedirect.com

ScienceDirect Physics Procedia 66 (2015) 52 – 58

C 23rd Conference on Application of Accelerators in Research and Industry, CAARI 2014

Radial Expansion of a Low Energy Positron Beam Passing Through a Cold Electron Plasma within a Uniform Magnetic Field F. F. Aguirre* and C. A. Ordonez. Department of Physics, University of North Texas, Denton, Texas 76203

Abstract The radial expansion of a low energy positron beam passing through a cold electron plasma within a uniform magnetic field is studied using a classical trajectory Monte Carlo simulation. It was found that under certain conditions there is a range of impact parameters where electrons and positrons may form weakly bound states and drift perpendicular to the magnetic field before dissociating. The purpose of the simulation is to evaluate the diffusion coefficient due to the formation and dissociation of those states. In addition, the effects on the diffusion coefficient of varying the kinetic energy and the magnetic field are evaluated. © 2015 The Authors. Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license © 2014 The Authors. Published by Elsevier B.V. (http://creativecommons.org/licenses/by-nc-nd/4.0/). Selection and peer-review under responsibility of the Organizing Committee of CAARI 2014. Selection and peer-review under responsibility of the Organizing Committee of CAARI 2014 Keywords: Diffusion coefficient, magnetobound states of positronium, antihydrogen mixing experiment.

1. Introduction The growing success in the production and study of antihydrogen is associated with the sophistication of the techniques used for trapping, transferring, and mixing charged particles [1]. Penning traps are widely used for trapping of charged particles within strong magnetic fields. A better understanding of the collisional processes in these traps (e.g, Ref. [2]) can help in the design of experiments within them. An aspect that has been recently investigated is cross field transport. Important improvements were made to the theory of collisional transport , where is the cyclotron radius, is including collisions with impact parameters in the range off _________________ * Corresponding author. Tel.: +940-2318207; fax: +(940) 565-2515. E-mail address: [email protected]

1875-3892 © 2015 The Authors. Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/). Selection and peer-review under responsibility of the Organizing Committee of CAARI 2014 doi:10.1016/j.phpro.2015.05.009

53

F.F. Aguirre and C.A. Ordonez / Physics Procedia 66 (2015) 52 – 58

the impact parameter and is the Debye length [3-4]. References [5, 6] report studies on collisions involving antimatter particles. It has been predicted that formation of quasibound states may affect the cross field diffusion. Cross field diffusion has also been studied recently in Rydberg gases [7]. The present work is a study of the radial expansion that occurs during the interactions of two counterstreamming beams (one electron and one positron) within a uniform magnetic field using a classical trajectory Monte Carlo (CTMC) simulation. This study is conducted in the center of mass frame of reference, and the results can apply in a laboratory frame of reference where the electrons comprise a cold plasma. 2. Binary Collisions This section reports the phenomenon of giant cross-magnetic-field steps in collisions involving an electron and its antiparticle, a positron. Simulations results show that this phenomenon may occur at low energies in strong magnetic fields under conditions similar to the ones inside Penning traps in current antihydrogen experiments. The following describes the conditions and results of the simulation. The collision of a positron and an electron is treated classically, considering the electric interaction under the influence of a uniform magnetic field. According to Coulomb’s law, the electric force exerted on particle by is given by (1)

where

is the Coulomb force constant, and are the charges is the vector position of particle relative to particle . In SI units, permittivity of free space. The magnetic force on each particle is given by

is the relative distance, and where is the

(2)

Here, is the Lorentz force constant ( in SI units), is the magnitude of a uniform magnetic field in the ) are the unit vectors of a Cartesian coordinate system, and are velocity components of direction, ( particle . , where is the The motion of the particles is governed by Newton’s second law, is its acceleration. Therefore, the acceleration of particle is given by mass of particle and (3)

The velocity and position of each particle are functions of time. If the position is written as and the time derivative is written as then the equations of motion are (4)

(5)

(6)

where

,

applies for the positron and

,

applies for the electron.

54

F.F. Aguirre and C.A. Ordonez / Physics Procedia 66 (2015) 52 – 58

The electron and positron are treated (but not simulated) as having traveled in opposite directions from two sources located at infinite distances from each other and at infinite distances from the coordinate origin. At infinite separation, the electric potential energy is defined to be zero. When the particles are separated by a finite distance at the start of a simulation, conservation of energy requires (7)

Here are the components of the initial velocity of the positron, are the components of the initial velocity of the electron, the kinetic energies for the two particles at infinite separation are denoted , respectively, and is their initial distance of separation. and The positron and the electron are considered to be initially positioned at Cartesian coordinates and respectively, and . The , , following values are used for the present study: and . The initial conditions in the present study represent oppositely directed particle is referred to as the initial axial separation. beams. Hereafter, is referred to as the impact parameter, and Equation (7) is written as

(8)

Rearrangement provides an expression for the nonzero component of the initial velocity of the positron (with ) and the electron (with ),

(9)

Each simulation is stopped at a time , when the axial separation becomes larger than the initial axial separation. The condition for the simulation to stop is when the inequality (10)

T, = = 6κ, , where has is first satisfied. The parameter values used in the simulation are the value of Boltzmann’s constant but the unit of energy. An output of the simulation is the cross-magnetic-field drift distance experienced by each particle. The crossmagnetic-field drift distance is (11)

The cyclotron radius is useful for normalizing lengths. The cyclotron radius is defined for the present study as the radius of a circular trajectory when one charged particle is under the influence of the magnetic field and has a kinetic energy that is only associated with motion perpendicular to the magnetic field. The cyclotron radius is , where and are the mass and kinetic energy of the particle. The equations of motion and initial conditions have been written in such a way as to also be applicable for simulating a collision between two positrons. Each simulation is run for both electron-positron and positron-positron binary collisions. Figure 1 contrasts drift distances for electron-positron and positron-positron (or electron-electron) collisions. The in increments of . For electron-positron collisions, giant drift impact parameter was varied from to to . Outside of this range, the distances occur at values of the impact parameter ranging from less than

F.F. Aguirre and C.A. Ordonez / Physics Procedia 66 (2015) 52 – 58

electron-positron drift distances drop to values closer to the corresponding drift distances for positron-positron (or electron-electron) collisions. Although the range of impact parameters within which giant cross-magnetic-field steps occur is well defined, the drift distances are sensitive to small changes in the impact parameter.

Figure 1: Normalized cross-magnetic-field drift distance versus normalized impact parameter.

Figure 2. Electron and positron trajectories rajectories pprojected onto the y-z plane when .

55

56

F.F. Aguirre and C.A. Ordonez / Physics Procedia 66 (2015) 52 – 58

Figure 2 shows the electron and positron drifting together as a correlated pair across the magnetic field. Each plane. Here, a drift distance of is particle oscillates several times with variable amplitude across the observed for the electron-positron collision. The electron-positron collision results in a drift distance that is two orders of magnitude larger than that for a positron-positron (or electron-electron) collision having the same impact parameter. To summarize, these results show that giant cross-magnetic-field steps can occur as a result of positronelectron collisions. Simulations predict the existence of a continuous range of impact parameters within which the system experiences giant cross-magnetic-field drift distances. The drift distances can be to two orders of magnitude larger than that associated with like-charge collisions. 3. Cross Field Diffusion Coefficient The diffusion coefficient due to formation and dissociation of temporarily bound states of positronium characterizes the rate at which each species expands across the magnetic field as a result of these binary interactions. The drift is modeled as a random walk process and the positron diffusion coefficient is evaluated as (12)

Here, is the electron density, is equal to , where is the maximum impact parameter at which , is the initial relative speed, is the azimuthal temporarily bound states are formed for a given value of and represents an average over the impact parameter, azimuthal angle and velocity distributions. angle and Considering the monoenergeticity of the beams, Eq. (12) reduces to (13)

From symmetry in the azimuthal direction, (14)

The diffusion coefficient reduces to .

(15)

For the evaluation of the ratio , around trajectories were produced for each set of parameters . A monitoring of a constant of the motion, the total energy, was done for some of the trajectories involving and to . long drifts. The relative deviations were typically within the orders for an initial energy equal to 6κ. For this case, the magnetic Figure 3 shows the variation of the ratio field was varied from 1 to 5 T in units of 0.5 T. Figure 4 shows that the mean-square of the steps decreases almost linearly with the magnetic field.

Figure 3. Ratio

versus .

F.F. Aguirre and C.A. Ordonez / Physics Procedia 66 (2015) 52 – 58

Figure 4.

versus .

for a fixed magnetic field (1 T) by varying the energy from Figure 5 shows the variation of the ratio 3 to 7 . In this case, it was observed that the mean-square of the steps approximately decrease quadratically with magnetic field (Fig. 6).

Figure 5. Ratio

Figure 6.

versus

versus

.

57

58

F.F. Aguirre and C.A. Ordonez / Physics Procedia 66 (2015) 52 – 58

As an example application of the simulation results, the diffusion coefficient for radial expansion is for an electron density of , an initial kinetic energy of and a magnetic field of 1 T. A possible implication of the simulation results is related to radial expansion when preparing an antihydrogen mixing experiment (for example, when antiproton clouds are driven into the positron plasma in a Penning trap). Since a cascade of long drifts may expel any remaining electrons from antiproton plasmas when overlap with a positron cloud commences, if the number of electrons is sufficiently small compared to the number of positrons, the electrons could be quickly ejected. Acknowledgment This material is based upon work supported by the Department of Energy under Grant No. DE-FG02-06ER54883 and by the National Science Foundation under Grant No. PHY-1202428. References [1]G. B. Andresen, M. D.Ashkezari, M. Baquero-Ruiz, W. Bertsche, P. D. Bowe, E. Butler, C. L. Cesar, S. Chapman, M. Charlton, A. Deller, et al., Nature 468, 673 (2010). [2]M. E. Glinsky, T. M. O`Neil, M. N. Rosenbluth, K. Tsuruta, and S. Ichimaru, Physics of Fluids B: Plasma Physics 4, 1156 (1992). [3]C. F. Driscoll, F. Anderegg, D. H. E. Dubin, D.-Z. Jin, J. M. Kriesel, E. M. Hollmann, and T. M. O`Neil, Physics of Plasmas 9, 1905 (2002). [4]Y. Ahat, C. E. Correa, and C. A. Ordonez, Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms 261, 252 (2007). [5]C. R. Correa, J. R. Correa, and C. A. Ordonez, Physical Review E 72, 046406 (2005). [6]J. R. Correa and C. A. Ordonez, Physics of Plasmas 21, 082115 (2014). [7]Y. V. Dumin, Physical Review Letters 110, 033004 (2013).