## The Magnetic Toroidal Sector as A broad-band Electron-Positron Pair Spectrometer I. Lepton Trajectories

S. Hagmann 1, 2 , P.M. Hillenbrand 1, 3 , Y. Litvinov , U. Spillmann , Th. Stöhlker 1 1 1, 4, 5

1 GSI Helmholtzzentrum Darmstadt, 64291 Darmstadt,  Inst. für Kernphysik, 2 Universität Frankfurt, 60438 Frankfurt,

3 Universität Giessen, 35392 Giessen,  Helmholtz Institut Jena, 07743 Jena, 4 5 Fakultät für Physik, Universität Jena, 07745 Jena

We  report  an  analysis  of  electron-optical  properties  of  a  toroidal  magnetic  sector spectrometer and examine parameters for its implementation in a relativistic heavy-ion storage ring like HESR. For studies of free-free pair production in heavy- ion atom collisions  this  spectrometer  exhibits  very  high  efficiencies  for  coincident  e + -e pair spectroscopy  over  a  wide  range  of  momenta  of  emitted  lepton  pairs.  The  high coincidence efficiency of the spectrometer is the key for stringent tests of theoretical predictions for the coincident positron- and electron emission characteristics and for the phase space correlation of lepton vector momenta in free-free pair production.

1. Introduction

2.

momentum calibration of intersections of trajectories on

Appendix 1

Appendix 2

8. References

57

59

| 1. Introduction                                                                                                | 1. Introduction                                                                                                | 3   |
|----------------------------------------------------------------------------------------------------------------|----------------------------------------------------------------------------------------------------------------|-----|
| 2.                                                                                                             | 2.                                                                                                             |     |
| 3. The high energy storage ring HESR                                                                           | 3. The high energy storage ring HESR                                                                           | 11  |
| 4. Some kinematical considerations for free-free pair production by ions in the relativistic storage ring HESR | 4. Some kinematical considerations for free-free pair production by ions in the relativistic storage ring HESR | 13  |
| 5. A magnetic toroidal sector for e + and e - in a storage ring                                                | 5. A magnetic toroidal sector for e + and e - in a storage ring                                                | 21  |
| 6. Invariants and other serendipitous features of lepton trajectories in                                       | a magnetic toroidal sector, towards 3DCS 𝑑 3 𝜎 𝑑𝐸 + 𝑑𝛺 + 𝑑𝐸 -                                                  | 32  |
| a)                                                                                                             |  lab =0; invariance of momentum arcs                                                                          | 32  |
| b)                                                                                                             |  lab >0                                                                                                       | 35  |
| c)                                                                                                             | momentum calibration of intersections of trajectories on detector planes                                       | 39  |
| d)                                                                                                             | geometric invariance of foci on arcs                                                                           | 48  |
| e)                                                                                                             | differential cross sections in emitter frame                                                                   | 50  |
| f) Impact parameter dependence of free-free pair production                                                    | f) Impact parameter dependence of free-free pair production                                                    | 51  |
| 7. Summary                                                                                                     | 7. Summary                                                                                                     | 54  |
| Appendix 1                                                                                                     | Appendix 1                                                                                                     | 54  |
| Appendix 2                                                                                                     | Appendix 2                                                                                                     | 57  |
| 8. References                                                                                                  | 8. References                                                                                                  | 59  |

## 1. Introduction

Electron -positron pair production, the excitation of an electron with negative energy from the completely filled Dirac sea  into an unoccupied state with either a discrete positive  energy  (resulting  in  bound-free  pairs)  or  positive  energy  in  a  continuum (resulting in free-free pairs) has evolved into a central topic of QED in extreme fields as the coupling between the lepton field and the electromagnetic field is close to one [1-3]. Already in the early times of the development of quantum mechanics and QED this process was assessed in the context of the Dirac-, Breit-Wheeler- and vacuumpolarization-in-strong-  electromagnetic-fields-  fundamental  processes,  and  also  the highly interesting and very profound relation of this particular channel in fast ion-atom collisions with other inelastic processes, e.g. electron nucleus bremsstrahlung (eNBS), has been thoroughly investigated early on [4-9]. These seminal calculations on eNBS have provided then in turn guidance for particular aspects in the evaluation of transition probabilities  in  pair  creation  calculations  [19].  It  is  within  this  context  that  pair production  in  ion-atom  collisions  with  projectile  energies  ranging  from  the  extreme adiabatic regime close to the Coulomb barrier ( i.e. below 6AMeV) up to energies in the extreme relativistic regime with collisions characterized by  few 1000 has since been subject of numerous theoretical treatises and experimental studies [10-51]. Most recently the surprisingly high cross sections even exceeding 100kb observed for pair production  in    highly  relativistic  very  heavy  ion  atom  collisions  have  aroused  new interest in this process even in the high-energy accelerator community as capture from pair production has been identified as a critical, potentially luminosity limiting process in relativistic heavy ion colliders like LHC and RHIC [13,19,21,55,56]. The enormous cross sections observed can be traced [1,3,13,19,21] to the large transverse electric fields E transv  (see fig. 1) which give rise to the unexpectedly high cross section  pair production for free-free pair production.

Fig.1. Dominance of transverse electric fields for a moving charge with increasing  . The  current  equivalent  collision  energy  range  of  interest  for  the  free-free  pair production  in  the  HESR  is  indicated.  The  strongly  adiabatic  region  for  positron production  at  6  AMeV  with  0.1,  i.e.  (  -1) « 1  and  as  well  pair  production  in  the extreme relativistic region with  » 100 are not subject of our present considerations.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01019v1-figure-89a74fbc-55ac-4cb0-8e73-17206208f251.jpg)

The relativistic  contraction  of  the  collision  time  coll =b/ vcoll  (  where  b  is  the  impact  parameter) combined with the high transverse electric field Etrans  Zproj/b 2 makes the electron positron pair-creation process in collisions of near relativistic heavy ions with atoms particularly attractive as a tool to investigate in high resolution correlated lepton dynamics for collision times t coll &lt;10 -18 s and at collision velocities vcoll larger than or large in comparison with the orbital velocity vK of the inner most bound electron v coll &gt;v K .

In  the  entire  discussion  following  this  introduction  we  will  only  be  concerned  with collisions with 1.2  6, i.e. covering the range of collisions energies for ions stored in the future HESR storage ring [57-59]. We are not considering electron positron pair creation in those other two collision regions which also have attracted recent attention, due the new windows opened for fundamental atomic processes by  intricate acceldecel modes in modern advanced accelerator facilities: these are very highly adiabatic collisions with vcoll&lt;&lt;vK ( this corresponds to collision energies  6AMeV, e.g for U 91+ projectiles), to be investigated at the future CRYRING/ESR, where during the collisions even a superheavy (i.e. Z1+Z2=ZUA  173) quasimolecule may be formed [1, 60-66], and

second,  ultraperipheral  collisions  with  1000-3000,  as  now  routinely  available  at CERN and RHIC [1, 14, 20, 21, 55].

At the future near relativistic storage ring HESR at FAIR the collision energy range up to  6 corresponding to Eproj  5AGeV will be available. For ion-induced pair-production in this collision velocity range theoretical calculations [5, 10, 11, 13-15, 18, 33, 36, 37, 41,  44,  46]  have  identified  various  major  channels  which  may  be  distinguished experimentally (see fig .2),

- i) free-free pair production,
- ii) bound-free pair production and
- iii) negative continuum dielectronic recombination:

Fig. 2. Energy diagram of the single particle Dirac equation and illustration of various open  channels  for  electron-positron  pair  production  besides  the  classical  electron capture, ~excitation and ~ionization active in swift ion-atom collisions [57, 58, 63].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01019v1-figure-1b2ce944-0187-4d4a-a758-5303ce049ff5.jpg)

The free-free pair creation described in eq.1a,1b is the main motivation for the current spectrometer design study strongly emphasizing a large combined detection efficiency (and solid angle) for simultaneous detection of electron and positron of a free-free pair: a  coincident  detection  of  both  outgoing  leptons  and  their  vector  momenta  -  as constitutive in the identification of free-free pair production - is not practically possible with conventional dispersive instruments covering only a small solid angle for each lepton because of the intrinsic very low pair coincidence efficiencies resulting from such configurations.

In free-free pair production both leptons are emitted into continuum states of Ekin &gt; 0:

$$X ^ { Z + } + A \rightarrow X ^ { Z + } + \{ A ^ { * } \} + e ^ { + } + e ^ { - } \quad f r e e { \text{-free} } \quad \text{ (1a)}$$

Theory predicts a scaling of the total cross section for free-free pair production [10, 13,

$$\sigma _ { f r e e - f r e e } \sim \left ( Z _ { p r o j } ^ { \quad 2 } \right ) ( Z _ { t a r } ^ { \quad 2 } ) l n ^ { 3 } \gamma \ \.$$

For this channel a total cross section of  3b has been measured for 940 AMeV Au +Au collisions by Belkacem et al. [34, 40, 54]. The complexity of two-lepton continua here presents a formidable impediment to theoretical differential cross sections and as well a  tremendous  obstacle  for  experimental  full  differential  cross  sections  (see  next section).

For positron production from bound-free pair creation channels according to eqs. 2a and 2b, the differential cross sections on the other hand can be well studied using traditional magnetic spectrometers at the HESR storage ring [57-60, 80] as here only one lepton, the positron, is found in the continuum, with a resulting satisfactory overall coincidence efficiency when employing traditional magnetic spectrometers.

In the bound- free pair production channel the electron of the lepton pair is captured into a bound state of the projectile while the positron is emitted into the continuum.

$$X ^ { Z + } + A \rightarrow X ^ { ( Z - 1 ) + } ( n l ) + \{ A ^ { * } \} + e ^ { + } \quad \ \ ( n l ) \text{mostly} \ ( 1 s ) \ ; \text{bound-free} \quad \ \ ( 2 a )$$

Here theory predicts a scaling of the total cross section following

$$\sigma _ { b o u n d - f r e e } \sim ( Z _ { p r o j } ^ { \quad 5 } ) ( Z _ { t a r } ^ { \quad 2 } ) l n \gamma$$

and the total cross section, again as measured by Belkacem et al for 940 AMeV Au +Au collisions, is  2b [34, 40, 54].

The  third  important  pair  creation  channel  is  the  negative  continuum  dielectronic recombination (NCDR) channel as described in eq. 3 and has not been experimentally observed yet; it was only recently proposed [ 67-69]

$$X ^ { Z + } + A & \to X ^ { ( Z - 2 ) + } + \{ A ^ { + * } \} + e ^ { + } & N C D R &$$

The NCDR channel has been predicted to exhibit a broad maximum as a function of  [60,  67-69] and not a monotonic dependence of the total collision cross section as established for free-free and bound-free pair production according to eqs. 1b and 2b. The NCDR is unambiguously characterized by a double capture into the  projectile coincident with positron emission. Its total cross section maximizes at approximately 3mb at around 1AGeV U 92+ [67], this is much lower than the two other channels (1) and (2).

All  three channels can be experimentally distinguished by their characteristics via  e /e + - X q+ coincidences, where the coincident projectile final charge state is q= Z, Z-1, Z-2, respectively.

The e + -e -pair production in the field of a nucleus induced by a high energy photon is -dominated by near back-to-back emission of both leptons and a near equal energy sharing 𝐸 𝑒 + ≈ 𝐸𝑒 of  electron  and  positron  [1].  For  free-free  pair  production  in  the transient strong electro-magnetic field created in fast ion-atom collisions, however, the resultant emission patterns exhibit significantly more complex features (18,28,29). In fig.3 we illustrate for the collision system 0.96AGeV U 92  + Au the total positron yield per  10 0   interval  as  function  of  the  laboratory  observation  angle  as  measured  by

Belkacem et al.[34]; the strong kinematic maximum at 1/  rad is clearly visible. 

Fig. 3. Relative positron production cross section in U 92  + Au collisions at 0.96 AGeV, data are from A. Belkacem et al., ref. 34.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01019v1-figure-f74ff4cf-4b38-4236-a1e0-60c3e5c3b8b1.jpg)

Few  theoretical  differential  cross  sections  for  free-free  pair  production  have  been reported  [18,  28,  29,  36],  however,  and,  as  follows  form  the  nature  of  the  full  6-

dimensional phase space and corresponding limitations for numerical computations, mostly for mutually not overlapping, and quite narrow subspaces. For the differential cross sections for pair creation in their fundamental relation to the high energy endpoint of enBS already Pratt et al. [50] and Fano et al. [51] in the early 60s have presented some general considerations pertaining to the single energy-differential cross sections with  strongly  different  characteristics  for  high  energy  end  of  positron  or  electron spectra;  these  far  reaching  predictions  have  never  been  subject  to  experimental scrutiny.

For e + -e  -pair production in relativistic ion-atom collisions the group of W. Scheid [28, -29] performed extensive calculations for effectively double differential cross sections yielding  some  characteristics  of  free-free  pair  production.  These  authors  used  a semiclassical  approach  and  solved  the  time  dependent  Dirac  equations  in  the semiclassical  approach  employing  free  wavepackets  in  a  continuum  discretized  in momentum space. They observe that there are no straightforward boundary conditions nor propensity rules on energy sharing and direction of emission, on the contrary, a very complex relation between the angular emission patterns of electron and positron emerged from Scheid's theoretical calculations for collision systems 10 AGeV Au 79+ + Au 79+ .  One  of  the  quite  interesting  findings  is  that  leptons  are  emitted  with  mean energies &lt;E lepton &gt;  &lt; 1MeV and that even at these high collision energies only a small fraction of the emission cross sections is seen with high energy leptons &gt;1MeV.

For further discussions of essential and indispensable spectrometer parameters we focus  on  two  surprising  and  highly  remarkable  results  reported  by  Scheid's  group, which are playing a decisive role in our further considerations: the angular distribution of positrons from the above reaction (integrated over all electron energies and electron emission  directions)  exhibits  two  peaks  in  the  forward  and  backward  direction  at around 50 and around 130 degree, respectively, superimposed over a near isotropic distribution (see fig. 4a) [ 28, 29].

Fig. 4a. Angular distribution of positrons in 10 AGeV Au 79+ + Au 79+ integrated over all energies and directions of the coincident electron [28, 29], discretization of momentum space results in angular bins of  =15 0  .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01019v1-figure-05c6f6d0-35eb-4d44-9375-0a0b037274c1.jpg)

An even more remarkable result was found for the distribution of difference angles 𝜃𝑒 + ,𝑒 -for the emitted electron and positron of a pair. For free-free electron- positron pairs  in  the  investigated  system  it  could  be  established  that  there  is  definitely  no evidence for a preferred back to back emission; their calculations rather exhibit a very broad  distribution.  Even  lepton  pairs  with  near  vanishing  difference  angles  are produced with a significant probability.  A  propensity  for difference  angles  in a  very broad distribution around 90 degree is observed in fig.4b; we also note a small but nonvanishing  share  of  pairs  to  appear  with  back  to  back  emission.  Tenzer  et  al.  also showed that only for transverse emission of the electron the corresponding positron is emitted near 160 0  with respect to the electron polar emission angle, whereas for other polar emission angles of the electron no such propensity is found [28, 29].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01019v1-figure-96d81cfc-56dd-421d-936a-aead3cac56a8.jpg)

Fig. 4b. Distribution of difference angles Δ𝜃𝑒 + ,𝑒 -in 10 AGeV Au 79+ + Au 79+ for arbitrary energy of leptons exhibiting that in ion-atom collisions there is no clear preference for back to back emission of the coincident leptons [28, 29]. Note, that a given difference angle Δ𝜃 𝑒 + ,𝑒 -for a lepton pair does not predict an absolute orientation angle 𝜃 𝑒 or 𝜃 𝑒 + for a lepton of the pair with respect to the beam axis.

Similarly, it had also been shown by H. Bhabha that electrons and positrons created in swift collisions of ions and atoms [25] emerge at angles independent of each other with a propensity for small pseudorapidity, i.e. perpendicular to the direction of the projectile beam, and a rather steep decrease with the lepton energy. Another differential cross section  for  a  near  symmetric  and  very  heavy  collision  system  was  studied  by  D. Ionescu in his Habilitation (18, 36 see fig 5 which only covers a small section of the entire available phase space, with a rather high kinetic energy of both leptons).

Fig. 5. Angular distribution of free electron positron pairs as function of positron angle with  electron  angle  e as  parameter,  for  lepton  energies  Eelectron=2.01mc 2   and Epositron=2mc 2 ,  respectively;  Sommerfeld-Maue  wavefunctions  for  the  free  electron positron states are used in this perturbative calculation [18].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01019v1-figure-b9a63892-e4a8-4224-a053-98f35217cf69.jpg)

These first theoretical benchmark calculations show that a detailed experimental study of the correlated vector momenta of coincidently emitted electrons and positrons in free-free pair production is urgently needed in order to execute a comprehensive test of the current most advanced theories.

Due to the formidable experimental obstacles there are up to this time no experimental full fourfold differential cross sections (4DCS) for free-free pair production

𝒅 𝟒 𝝈 𝒅𝑬 𝒆 +𝒅𝜴 𝒆 +𝒅𝑬 𝒆 -𝒅𝜴 𝒆 -, or 𝒅 𝟒 𝝈 𝒅𝒑 𝒆 +𝒅𝜴 𝒆 +𝒅𝒑 𝒆 -𝒅𝜴 𝒆 respectively, in the literature, fully differential in energy/momentum and emission direction of both coincident leptons. Such 4DCSs evidently would present the most complete information on the collision dynamics and would simultaneously constitute the most stringent test of theories. Similarly, on the theoretical side it is the immense computing power required for a true unrestricted twolepton continuum  process  which  has  purportedly kept theorists from actually calculating full 4DCS for free-free pair production but focus rather on some selected benchmarks  and  a  rather  coarse  momentum  and  angular  resolution,  as  shown  in Scheid's work.

Under these circumstances it appears attractive to first address for an experimental study threefold differential cross sections

$$3 \text{DCS} ( \mathrm e ^ { + } ) = \frac { d ^ { 3 } \sigma } { d E _ { e ^ { + } } d \Omega _ { e ^ { + } } d E _ { e ^ { - } } } \quad \text{and} \ 3 \text{DCS} ( \mathrm e ^ { + } ) = \frac { d ^ { 3 } \sigma } { d E _ { e ^ { - } } d \Omega _ { e ^ { - } } d E _ { e ^ { + } } }$$

or the corresponding 3DCS in momentum space, where one has integrated over the emission direction of one member of the lepton pair; these 3DCSs may serve as a first, moderately restricted though urgently needed experimental test of advanced theories.

Our present study serves to show the feasibility of a spectrometer for an experiment detecting  with  a  large  solid  angle  for  both  leptons  in  coincidence  electrons  and positrons from free-free pairs produced in relativistic ion-atom collisions and derive selected  3DCSs  and  partially  4DCSs  for  heavy  ion  atom  collisions  systems.  For experimental investigations using magnetic spectrometers the momentum space is the natural  frame  and  thus  in  the  next  sections  all  differential  cross  sections  will  be provided and discussed in momentum space, e.g. the 4DCS 𝒅 𝟒 𝝈 𝒅𝒑 𝒆 + 𝒅𝜴 𝒆 + 𝒅𝒑 𝒆 -𝒅𝜴𝒆 and

$$\text{corresponding } 3 \text{DCSs} \ 3 \text{DCSp} ( e ^ { + } ) \ = \ \frac { d ^ { 3 } \sigma } { d p _ { e ^ { + } } d \Omega _ { e ^ { + } } d p _ { e ^ { - } } } \quad \text{ and } \ 3 \text{DCSp} ( e ^ { - } ) \ =$$

𝒅 𝟑 𝝈 𝒅𝑬𝒑𝒆 -𝒅𝜴𝒆 -𝒅𝒑 𝒆 + ., for details see section 6d and for the corresponding transformations appendix 2.

## 3. The high energy storage ring HESR

The relativistic high energy storage ring for ions and antiprotons HESR was originally conceived to store and cool antiprotons for the PANDA collaboration at FAIR [58]; but very soon additional extended ion-optical calculations [57, 59, 60] proved the HESR to

be a highly promising facility for studies of atomic collision dynamics and few-electron high-Z atomic spectroscopy in the realm of near relativistic heavy-ion atom collisions. The  collision  energy  range  offered  by  the  HESR  is  unique  in  the  world  and  is particularly interesting for Coulomb dominated free-free pair production (see fig 6a). In the following we will only list a few essential parameters of this unique storage ring, which are useful for our further discussions, but we refer the inclined reader to the literature for more detailed information [57-60].

The HESR with a magnetic rigidity ranging from 5 to 50 Tm has a circumference of 550 m and can store highly charged ions up to U 92+  with energies between 170 AMeV and 5AGeV, corresponding to  max  6.

Fig 6a. Comparison of collision energy range of the future HESR storage ring with current Storage rings ESR and Cryring at GSI. It is nicely visible that the HESR extends from  near  relativistic  regime  into  the  pure  relativistic  regime  where  γ  is  clearly dominated by 𝜸 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01019v1-figure-aa50e7b5-a121-4521-9f14-93600695107f.jpg)

Beams from synchrotrons, the current SIS18, or the future SIS100, will be injected at 740AMeV energy and can be decelerated and/or accelerated to the energy desired by the respective experiments. Stochastic cooling and electron cooling will be possible between 740AMeV and 5AGeV in order to provide a beam with a momentum spread  p/p below 10 -4 . The storage ring is equipped with two long straight sections, one for the PANDA detector, and the other for a variety of experimental installations including

the electron cooler. At two alternative locations indicated in fig. 6 a supersonic jet may be implemented [60] for experimental installations ranging from high resolution x ray crystal spectrometers and calorimeters to magnetic pair spectrometers.

Fig. 6b. The relativistic high energy storage ring HESR with selected installations, the two options for the installation of the magnetic pair spectrometer in the arc section immediately in front of and following the electron cooler, respectively are indicated by red arrows, for details see [57, 59, 60] .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01019v1-figure-a8331fdd-d5eb-4525-91cb-5f485436c93e.jpg)

From  extensive  electron-optical  and  ion-optical  calculations  which  were  performed when implementing our magnetic 0 0  electron spectrometer [124] in the supersonic jet target zone of the current ESR storage ring we know that the effect of the magnetic field of the spectrometer on the coasting ion beam could completely be compensated by two small magnetic steerers in the ring lattice. Corresponding calculations for the toroidal spectrometer in the HESR are in progress.

## 4. Kinematics of free-free pair production in an emitter of near relativistic collisions velocity

The kinematic maximum in the positron spectra reported by Belkacem et al. [34, 40] and the calculations of Scheid et al.[28, 29] demonstrate that the continuum of the two emitted  leptons  for  free-free  pair  production  is  in  the  emitter  frame  dominated    by leptons  of  low  kinetic  energies.  It  thus  appears  indicated  to  revisit  for  relativistic collision systems the spectral features which are expected for the 2-lepton continua in the laboratory and which can be generated by the highly diverse theoretically predicted

emission patterns of a free-free pair in the projectile emitter frame. For near-relativistic collision  systems  with  (  -1)&gt;&gt;0  the  Lorentz  transformation  from  the  emitter-  to  the laboratory frame complicates the easy graphic derivation of kinematic relations so well established  from  Newton  diagrams  in  non-relativistic  kinematics.  A  relativistically correct graphic representation[71 - 74] does nevertheless also exhibit all the important material features, e.g. a maximum laboratory observation angle  max, which is suitably calculated following 𝑠𝑖𝑛𝜃𝑚𝑎𝑥 = 𝛽′𝛾′ 𝛽𝛾 (  ',  ' for the lepton in the emitter frame,  ,  for the emitter in the laboratory frame) for a slow lepton emitted by a fast projectile and also the momenta of leptons in the laboratory as function of laboratory observation angle;  this  instructive  and  nearly  indispensable  tool  [74]  had  a  high  practical  value before the advent of computers, today it retains a very desirable didactic value (see fig. 7), though.

Fig. 7. Illustration of useful kinematic relations for momenta in emitter and laboratory frame in relativistic kinematics

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01019v1-figure-49353668-6119-450d-b4c9-cf155673e95a.jpg)

In order to illustrate the features adequate for a useful pair spectrometer we consider as a first step the kinetic energy and the angular distribution in the laboratory expected for a hypothetical scenario at the future HESR: we begin with leptons emitted with fixed energy Ekin=500 keV in the projectile frame for 5.1AGeV projectile energy.

Using the relations

$$& E _ { e, k i n } ^ { l a b } = \gamma \left ( E _ { e, k i n } ^ { p r o j } + \beta \sqrt { \left ( E _ { e, k i n } ^ { p r o j } \right ) ^ { 2 } + 2 E _ { e, k i n } ^ { p r o j } * m c ^ { 2 } } \cos \theta _ { p r o j } \quad \right ) + \\ & ( \gamma - 1 ) m c ^ { 2 }$$

or in a more compact form using momenta

$$E _ { e, k i n } ^ { l a b } = ( \gamma \gamma _ { e } ^ { p r o j } ( 1 + \beta \beta _ { e } ^ { p r o j } \cos \theta _ { p r o j } ) - 1 ) \ m c ^ { 2 }$$

𝑙𝑎𝑏

for the dependence of the kinetic energy of the lepton in the laboratory frame 𝐸𝑒,𝑘𝑖𝑛 on energy and emission angle in the emitter frame, 𝐸 𝑒,𝑘𝑖𝑛 𝑝𝑟𝑜𝑗 and 𝜃𝑝𝑟𝑜𝑗 ,respectively, and

$$\sin \theta = \frac { \sin \theta _ { p r o j } } { \sqrt { \gamma ^ { 2 } ( c o s \theta _ { p r o j } + \frac { \beta } { \beta _ { e } ^ { p r o j } } ) ^ { 2 } + s i n ^ { 2 } \theta _ { p r o j } ) } } \quad \text{ or }$$

$$t a n \theta = \frac { \beta _ { e } ^ { p r o j } s i n \theta _ { p r o j } } { \gamma ( \beta _ { e } ^ { p r o j } c o s \theta _ { p r o j } + \beta _ { 0 } ) } \quad$$

respectively, for the relation between emitter frame emission angle 𝜃𝑝𝑟𝑜𝑗 and laboratory observation angle  we arrive at the instructive graphical relation in fig. 8.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01019v1-figure-b8f0ea28-8a06-469d-bbe6-3d24111f14e3.jpg)

Fig.  8.  Kinematical  relations  for  leptons  emitted  by  a  fast  emitter,  e.g  during  pair creation: a hypothetical e  - e  pair with + -𝐸 𝑒,𝑘𝑖𝑛 𝑝𝑟𝑜𝑗 = 500 keV for each leptons is emitted in a moving frame of 5.1 AGeV (  6). If such a pair of leptons is assumed to be emitted back to back e.g. at 60 0  and 120 0  in the emitter frame (dashed lines) it will appear in the laboratory as coincident pair with 9 MeV at 4.5 0  and 3.3 MeV at 11.5 0 , respectively.

We observe that the lepton emission in the fast moving emitter produces as a general feature the well-known narrow forward cone of leptons in the laboratory with a wide range of kinetic energies. It is apparent that even  very strict kinematical conditions assumed for the emission characteristics of a hypothetical lepton pair in the emitter frame will still encompass a wide range of kinetic energies of leptons, e.g. between 400keV  and  20MeV  which  will  have  to  be  covered  by  actual  lepton  detectors,  as illustrated  in  the  example  in  fig.  8.  However,  only  when  assuming  back-to-back emission of both leptons (and with a further restriction of both leptons emitted with the same  kinetic  energy  in  the  emitter  frame),  it  is  straightforward  to  determine  from coincident measurement of the laboratory energies of the leptons the emitter frame energy and emission angle, as is illustrated in fig. 9. On the other hand, a hypothetical assumption of equal energies in the emitter frame with no assumption on emission angles, though, will not result in a unique complete identification of primordial emission characteristics.

For ion-atom collisions, however, the open phase space spanned by the two leptons of a free-free pair in the continuum is very large and is practically not constrained by any kinematic conservation rules, such as back to back emission of leptons in the emitter  frame.  Theory  [24-25,  28-29,  36]  has  predicted  that  the  two  leptons  of  a coincident pair may be indeed emitted with a nearly unrestricted, very wide range of difference angles (in the emitter frame) and with no strict rules applying, neither for respective emission directions, nor for energy sharing. We illustrate for typical lepton energies in the emitter frame for 5.1AGeV Au + Au, i.e. 200 keV, 500 keV, 1MeV and 2 MeV, the kinematic relations resulting in the laboratory for emission of leptons over a wide range of kinetic energies, in fig. 9; this illustrates the very large phase space which needs to be covered by a spectrometer when both leptons of the pair are to be detected in coincidence.

We wish to draw attention to two interesting features:

i)  for  lepton emission up to around 90 0  in the projectile frame the relation between emitter frame emission angle  proj and laboratory observation angle  exhibits only a very  weak  dependence  on  the  primordial  kinetic  energy  of  the  lepton.  For  a hypothetically preferred transverse emission at  90 0 of the free-free pair in the emitter frame, and equal energy sharing for the leptons and kinetic energies 𝐸 𝑒,𝑘𝑖𝑛 𝑝𝑟𝑜𝑗 of .5, 1 and  2  MeV,  respectively,  both  leptons  would  appear  under  laboratory  observation angles quite close to  =8 0  for all corresponding laboratory kinetic energies, 6, 10 and 17 MeV.

ii)  near  the  maximum  laboratory  scattering  angle,  i.e.  at  backward  emission  in  the emitter frame, the laboratory energy of the leptons is rather low as a consequence of kinematic  focussing  and  the  relation  Elab  Eemitterframe  for  the  leptons  is  inverted. Leptons from (e - ,e + )-free-free pair production with a low energy in the laboratory are all originating from high energy leptons emitted in the backward direction in the emitter frame  due  to  their  production  in  a  fast  moving  emitter.  This  has  a  considerable implication for the observation of leptons and the spectrometer design: the transverse momentum of  the  leptons  which  determines  the  magnitude  of  the  gyroradius  in  a magnetic field is rather small. This will be discussed in more detail in the context of the lepton trajectories.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01019v1-figure-a5da6550-244e-4f97-b3df-8016b6308402.jpg)

Fig. 9. Kinematical relations for leptons emitted by a fast projectile, for lepton energies 200keV,  500keV,  1MeV  and  2MeV  in  the  emitter  frame.  Note  the  wide  range  of laboratory  kinetic  energies  of  the  leptons  extending  beyond  30MeV  in  the  forward direction. The diagram shows that the large laboratory observation angles correspond to small laboratory energies for backward emission in the emitter frame  em  150 0 .

The electron branch of the toroid spectrometer will see high-energy electrons mainly from three processes:

- a)  -electrons from target ionization processes (i.e. bound state electrons of the target  ionized  by  the  time  dependent  perturbation  by  the  projectile  potential).  The momentum  transfer 〈𝑞〉 = 𝐸 𝑏 +𝐸 𝑐𝑜𝑛𝑡 ℏ𝑣 (E b initial  state  binding  energy,  E cont final  state continuum energy of electron, v collision velocity) necessary to generate high energy (  MeV)  -electrons  implies  that  these  electrons  originate  mostly  from  very  high momentum components of the Compton profile of the initially bound target electron [75] if they appear at MeV laboratory energies. Corresponding cross sections are very small.
- b) As well ionization of the projectile via electron loss to continuum (ELC) may contribute,  when  a  non-bare  projectile  carries  electrons  into  the  collision;  these electrons are easily identified: they appear in coincidence with a projectile whose final charge state is increased by n for n-fold projectile ionization and exhibit a narrow cusp at  0 0 with ve  vproj [124].
- c) Electrons from free-free pair creation mainly under consideration here will see some overlap in the same spectral region as the two above ionization channels a) and b). For free-free pair creation according to eq (1a) the projectile does not experience a charge exchange [36].

For the positron branch in the spectrometer positrons from weak interaction nuclear processes  may  arise  besides  those  of  atomic  origin;  however,  the  lifetime  of  the corresponding emitter states will lead to negligible decay before the ion passes the spectrometer entrance.

We illustrate  with  this  example,  that  due  to  the  lack  of  strong  kinematic  boundary conditions  in  lepton  pairs  produced  in  ion  atom  collisions  true  electron-positron coincidences may be located anywhere in the very extended 2-lepton phase space given by the very large range of 𝐸𝑒,𝑘𝑖𝑛 𝑙𝑎𝑏 and  lab (see figs. 9 and 10a,b). In order to

achieve reasonable electron-positron coincidence count rates in an experiment any new spectrometers will have to cover simultaneously a large range of kinetic energy ∆𝐸𝑒,𝑘𝑖𝑛 𝑙𝑎𝑏 , or electron momenta, respectively, and solid angle d  lab for the observation angles  lab without sacrificing the desired energy and angular resolution. For a lepton pair originating in a 5.1AGeV collision each lepton's phase space is depicted in fig. 10a.

The  largest  angles  lab  are  associated  with  low  lepton  momenta  and  will  thus  be mapped onto small perpendicular momenta/radii. It is the electron and positron from near forward emission at small  lab which will dominate the events at perpendicular momenta.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01019v1-figure-89f9fdf9-23ed-43d0-875e-3889353ab8b5.jpg)

Fig. 10a. Laboratory phase space for leptons emitted by a fast projectile: kinetic energy -observation angle range for free-free pair production for the case of projectiles with a specific kinetic energy 5.1 AGeV and a sample of four emitter frame energies of a lepton. For a lepton (e.g. e + ) with a given energy the laboratory observation angles of its coincidently emitted lepton (corresponding e - ) can cover a wide range according to Scheid [ 28-29] and the electron principally may appear anywhere in this graph with no a priori reasonable limit for the range of kinetic energy and emission direction provided by theory. It  is  apparent  that  a  combination of  two  independent  position dispersive instruments, even when equipped with 2D-PSD lepton detectors will have prohibitively small combined coincidence efficiency. The ticks for emitter frame angles 60 0 , 120 0 and 150 0  illustrate a favourable strong kinematic forward emission in the laboratory frame. Virtually the entire forward hemisphere is mapped into  lab  10 0 ; even at Eemitter= 2MeV the emitter frame emission up to 120 0  is mapped into lab angles  15 0 . For every emitter frame energy of a lepton the cross indicates on the abscissa the laboratory angle where the transverse momentum has its maximum. It is apparent that for the

leptons  the  large  laboratory  emission  angles  correspond  to  a  small  laboratory momenta and corresponding laboratory kinetic energy .

Even though the measured laboratory kinetic energy combined with the laboratory observation angle of a detected lepton permits to uniquely identify its primordial vector momentum in the emitter frame it is apparent that  traditional  spectrometers  with  a meaningful energy  and angular resolution will have a prohibitively small coincidence efficiency, due to their type of position dispersion and small solid angle.

Fig.  10b.  Relation  of  ptransverse  and  ptot= 𝛽𝛾 in  the  laboratory  of  leptons  emitted  by 5.1AGeV projectiles. A lepton is uniquely and  completely characterized by ptransverse and ptot and thus its vector momentum in the emitter frame can be reconstructed. For emitter  frame  kinetic  energies  of  2MeV  and  1MeV,  respectively,  leptons  with  low laboratory momenta p= 𝛽𝛾 &lt;30 and 𝛽𝛾 &lt;20, respectively, present small fractions of the total  cross  section  with  emitter  frame  backward  emission  cones  of  10 0   and  15 0 , respectively around the backward direction (see also fig 9).Note that the apex of each curve, the largest ptransverse, does not correspond to the largest assumed  lab.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01019v1-figure-98016066-3f48-46af-b1fa-d6695dcad53f.jpg)

The laboratory emitter angle  lab determines for a lepton's momentum 𝛽𝛾 its transverse component (𝛽𝛾)𝑡𝑟𝑎𝑛𝑠𝑣𝑒𝑟𝑠𝑒 , which will in the small angle paraxial approximation permit to  derive  an  approximate  radius  of  gyration 𝜚 in  a  magnetic  field  B  via 𝐵𝜚[𝑇𝑚] = 1.7045 ∗ 10 -3 (𝛽𝛾)𝑡𝑟𝑎𝑛𝑠𝑣𝑒𝑟𝑠𝑒 This results for perpendicular momenta (𝛽𝛾)𝑡𝑟𝑎𝑛𝑠𝑣𝑒𝑟𝑠𝑒 =4.8 (  Ekin= 2 MeV) in a radius of gyration of 𝜚  8 cm for a magnetic field B=1000G [and gives  equiv. (𝛽𝛾) = 0.48, 𝑖. 𝑒. 55 𝑘𝑒𝑉 in a 100𝐺 field ]. All  leptons  whose  momentum (𝛽𝛾)𝑡𝑟𝑎𝑛𝑠𝑣𝑒𝑟𝑠𝑒 is  smaller  than 586.7𝐵𝜚[𝑇𝑚] are  detected  within  a  circular  area  with

diameter 2 𝜚 in the detector plane. (Note for electron momenta pe[MeV/c]=0.511  and in atomic units pe[a.u.]=   / )

Even  position  dispersive  spectrometers  (e.g.  magnetic  sector  spectrometers  [7578,137])  with  2D  position  sensitive  lepton  detectors  can  only  partially  improve  the overall  coincidence  detection  efficiency  as  the  energy  or  momentum  window  both instruments see at any given time is small compared to the very large range of lepton momentum/energy  and  emission  angle  which  needs  to  be  covered  for  acquiring differential cross sections over a meaningful dynamic range.

## 5. The Magnetic Toroidal Sector

A major stimulus and valuable guide in pursuing a study of electron-optical properties of a magnetic toroidal sector for lepton pairs originating from relativistic heavy ion atom collisions  has been the TORI spectrometer conceived for experiments on adiabatic positron  production  in  superheavy  transient  quasimolecules  [78-80],  e.g.  collision energies  5.9AMeV  for  U q+ +Pb,  close  to  the  Coulomb  barrier;  TORI  exploits  a fundamental  property  of  a  toroidal  magnetic  field,  separating  particles  of  opposite charge  perpendicular  to  the  bend  plane  of  the  toroid;  TORI  was  developed  and operated in a series of experiments by the Kankeleit group at the UNILAC. In the TORI spectrometer electrons and positrons emitted by transient superheavy quasimolecules [64, 76, 77] created in a solid high Z target during adiabatic collisions with incident high Z heavy ions are guided in the first section by a toroidal +90 0  magnetic sector from the target to an intermediate plane where a baffle intercepts the electrons, which follow at this location trajectories geometrically well separated from those of the positrons. The positrons pass this midplane and enter a second section, point mirror symmetric to the first sector with respect to the midplane; then, in this second 90 0  toroidal sector the positrons are guided back onto the optical axis at the end of this section and there onto an array of PSD detectors [84-87,91].

The  excellent  suppression  of  background-electrons  arising  from  the  solid  target combined with a very high overall efficiency for positron detection are reason for the success  of  these  experiments  studying  positron  production  (and  delta-electron production  as  well,  in  a  corresponding  configuration  [75])  in  transient  superheavy quasimolecules created in adiabatic collisions of very heavy ions and heavy atoms [76, 77].

The fundamental property of a toroidal magnetic field, to disperse particles of opposite charge perpendicular to B and to  B into opposite directions and thus allowing to guide them to spatially isolated PSDs, appears to be also particularly suitable for experiments investigating  free-free-pair  production  planned  at  the  future  HESR  storage  ring  at GSI/FAIR; here, for free-free pair production, in a significant step beyond positron-only spectroscopy as used in TORI, vector momenta of the positron and the corresponding electron must  be  determined  simultaneously.  For  accomplishing  this  class  of experiments we require, however, just one toroidal sector, because such a toroidal Bfield configuration intrinsically permits to simultaneously analyse the vector momenta of both leptons of the pair in coincidence at a location of optimal separation (see fig 11), when using position sensitive detectors for both, electrons and positrons in the detector plane.

## B fields in midplane of solenoid-toroid

Fig. 11.  Cut through the mid-plane of the coil assembly of a toroidal magnetic sector spectrometer for simultaneous electron and positron spectroscopy at the supersonic gas jet target in the HESR storage ring. The centre line of the toroid is depicted as a solid blue circular arc- sector extending to coil 5 (location of 2D-PSD detectors). Two dashed lines indicate locations of  300mm from the centreline in the mid plane. The coil spacing in the solenoidal part reflects the option to facilitate access to the target zone for other detectors, e.g. x-ray and other photon detectors and as well 2D-PSD recoil  ion  detectors.  We  assume  for  illustration  of  the  geometry  a  local  cylindrical coordinate system (   , , z) with its origin at the centre of curvature of the toroid, here

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01019v1-figure-8a1e2cf6-5c28-4f0a-b443-23c7806ac128.jpg)

we have  =1000mm and  =0 at the entrance of the toroid. The corresponding z-axis is perpendicular to the centre plane, +z pointing out of plane. In this coordinate system the  magnetic  induction  B  in  the  toroidal  section  -  when  generated  by  an  effective meridional current in the coils - is ideally tangential 𝐵 ⃗ ⃗ ⃗ = 𝐴(𝑖) 𝑅 𝑒 𝜑 , where R=distance from the origin to the location of the field measurement.

For  implementation  of  a  magnetic  toroidal  spectrometer  inside  a  storage  ring  the boundary conditions of  operation  - a minimum  perturbation to  the  coasting  particle beam passing through several coils of the set (and subsequent compensation of this perturbation to the coasting beam using additional beam optical elements in the ring lattice downstream of the toroid) - impose significant limitations on the geometry for any  configuration  of  the  magnetic  coils  which  shall  produce  the  desired  toroidal magnetic field; obviously a single contiguous toroidal coil, although desirable, is not compatible with boundary conditions characteristic for  operation in a storage ring. (In storage rings the injection orbit of ions is in general quite distinct  [59, 60] from the orbit chosen for the data taking cycle, which may only commence after cooling and optional cycles for deceleration to the desired collision energy - thus imposing strict conditions on all in-ring spectrometer designs.) The transverse opening required by the storage ring  beam  optics  of  min  200mm  H  and  100mm  V  combined  with  the  desire  of minimizing the distortion of the homogeneity of the toroidal magnetic field configuration due  to  separated  coils  producing  the  toroidal  B  field  entails  a  toroidal  radius  R  of approximately 1000mm and corresponding coil sizes.

In order to maintain the flexibility for complementing the lepton pair spectrometer with a reaction microscope [88] containing a full 2D position sensitive detector in the target zone in the straight section [94] our design of up to 10 coils on an arc radius of 1000mm for generation of the toroidal field is complemented by up to 7 coils for a solenoidal section upstream of the supersonic jet target zone.

In an ideal magnetic toroid the magnetic induction B is tangential 𝐵 ⃗ = 𝐴(𝑖) 𝑅 𝑒 𝜑 assuming that the current in the coil is purely meridional, i.e. has no component along 𝑒 𝜑 . This normally  implies  a  contiguous  toroidal  coil  for  generating  this  field.  The  boundary conditions for in-ring spectrometers in the HESR storage ring, however, prohibit this solution.  We  have  thus  studied  the  magnetic  fields  produced  by  various  coil configurations. In fig. 12 and 13 we show - using calculations of a corresponding coil configuration shown in fig. 11 developed with TOSCA/OPERA-3D [81] - that the current configuration of 10+7 coils permits to generate a toroidal B-field with approximately

 B/B  1*10 -3 deviation from the ideal field over the centre arc and a region of a cross section of more than  200 x 200 mm around the centre arc perpendicular to the optical centre line from the target zone all the way to a detector plane. We investigate detector planes  positioned  between  =45 0   and  =83 0 ,  depending  on  the  relative  spatial separations of the leptons of different energies desired; in fig. 11 a detector array in position at  =66 0  (i.e. coil 5 of the toroid section) is schematically depicted.

An  array  of  multiple  2D  PSD-scintillators  [85-87]  with  energy  resolution  down  to  E/E&lt;0.01  provides  a  position  resolution  down  to  0.5  mm;  this  good  energy resolution for incident leptons adds for incident leptons highly useful redundancies in the correct lepton identification.

For experiment configurations requiring detection of leptons over a wide range of low kinetic  energies  in  the  laboratory  frame,  including  energies  close  to  threshold,  and simultaneously of very high energy leptons close to the beam direction, e.g. Cusp or binary  encounter  electrons,  an  option  for  a  detector  plane  at  =45 0 will  cover  all energies for leptons from threshold to very high energies at forward emission, however with a slightly decreased momentum dispersion for leptons.

Additional  simulation  runs  verified  that  optional  jet  target  positions  (desirable  for operation in conjunction with a reaction microscope) extending up to -500mm upstream (with  respect  to  the  position  (1000,0,0))  in  the  solenoidal  section  do  not  visibly deteriorate lepton momentum dispersion in the detector planes.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01019v1-figure-5dd8dbec-6009-4ae1-9599-113d33891994.jpg)

Fig. 12. Magnetic B-field uniformity: B-Field along centre line of solenoid (coil 1 to 3) and centre arc of toroid calculated using OPERA-3D for the coils shown in fig. 10. The abscissa 0 corresponds to the entrance of the first coil of the toroidal sector (=coil 4 counted from left), the detector position depicted in fig 10 (=coil 5 of toroidal section) is located at 1187mm on the arc. The deviation from the mean value does not exceed  B/B  1*10 -3 in the region of interest.

All magnetic field calculations were executed using the TOSCA section of the code OPERA\_3D\_10.5 which besides  calculating  the  magnetic  field  from  solving  partial differential equations with boundary conditions

$$\nabla \cdot \mu \nabla A - \nabla \cdot \mu \left ( \int _ { \Omega _ { J } \, | R | ^ { 3 } } \frac { J \times R } { | R | ^ { 3 } } \, d \Omega _ { J } \right ) = 0,$$

with  the  permeability  tensor  and A the  magnetic  potential,  also  calculates  and displays trajectories of charged particles in the calculated magnetic fields[81].

Fig. 13. The B-field of the coil configuration exhibits a toroidal behaviour over a large volume: The magnetic field on various arc sections for an inner radius R=800mm, a centre radius R=1000mm and an outer radius R=1200mm are shown. This shows that over  an  area  200mm  x  200mm  throughout  the  toroid  the  B  field  exhibits  a satisfactory toroidal behaviour; the relative deviations are particularly small not only on the centre line( as seen in fig. 12) but also in the important low field region outside of the centre line, where Reff is larger than the centre line radius R.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01019v1-figure-9dfe854f-af9e-46d9-9d6e-efebc61f08e2.jpg)

We have previously studied in conjunction with kinematic complete experiments in ionatom collisions in the few AMeV collision energy range the electron-optical properties of  a  toroidal  branch  combined  with  a  reaction  microscope  [82],  but  with  emphasis exclusively on very low energy electron spectroscopy.

## a) Toroidal coordinates

Before  analyzing  the  lepton  trajectories  in  solenoidal  and  toroidal  fields  it  may  be interesting to consider the choice of a suitable coordinate system. Given the toroidal geometry of the planned spectrometer we examine some useful properties of toroidal coordinates. For a toroidal geometry of the coil configuration with a purely meridional current  the  cylindrical  coordinate  system  frequently  used  in  calculations  for  a magnetically  confined  plasma  is  not  optimal  [95-97];  it  is  compelling  to  treat  the magnetic fields and the lepton dynamics in toroidal coordinates [98,122,123] (we use the nomenclature chosen by Moon and Spencer [98]):

$$x = \frac { a S i n h \eta \, C o s \psi } { C o s h \eta - C o s \theta } \quad y = \frac { a S i n h \eta \, S i n \psi } { C o s h \eta - C o s \theta } \quad z = \frac { a S i n \theta } { C o s h \eta - C o s \theta } \ \, \, 0 \leq \eta \leq \infty, \, - \pi \leq \theta \leq \pi, \, 0 \leq \psi \leq 2 \pi, \, \text{and}$$

a the scaling parameter for the foci (see figures 14a, b); the toroidal coordinates are a natural generalization of the bipolar coordinates 𝑥𝑝 = 𝑎𝑆𝑖𝑛ℎ𝜂 𝐶𝑜𝑠ℎ𝜂-𝐶𝑜𝑠𝜃 and 𝑦𝑝 = 𝑎𝑆𝑖𝑛𝜃 𝐶𝑜𝑠ℎ𝜂-𝐶𝑜𝑠𝜃 . For illustrative purposes we add two figures representing the corresponding coordinate surfaces.

Fig. 14a. Bipolar coordinates  and  , - ,  - . The curves of constant  (blue in fig. 14a) are non-intersecting, non-concentric circles around the foci -a and +a. For increasing |  | the radius of the circles decreases. The curves of constant  are similarly non-concentric circles intersecting the foci -a and +a; for increasing |  | the radius of

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01019v1-figure-b87ad821-0ddb-459d-859c-a0be2e37b16f.jpg)

the circle decreases. The curves of constant  and  ,  respectively, intersect at right angles;  they  turn,  upon  generalization,  into  toroidal  coordinates  with  toroids  and spheroids as coordinate surfaces.

The  fact  that  surfaces  of  constant  are  toroids  makes  toroidal  coordinates  very appropriate for investigation of magnetic fields generated by coils configured in toroidal geometry [99-102].

In fig 14a the parameters are selected such that the largest blue circles correspond to the actual coil geometry currently investigated with radius of curvature 𝜚 =1000mm and coil-radius  r=800mm and the foci calculated according to 𝑎 2 = (𝜚 - 𝑟)( 𝜚 + 𝑟) [100102] to be at a=  600mm. The set of blue circles are the geometric location of points with a fixed ratio of distance to the foci and, respectively, the set of red circles are the geometric  location  of  points  under  which  the  foci  appear  under  a  fixed  angle;  this system of orthogonally intersecting circles was elaborated on in antiquity as the circles of Apollonius [103].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01019v1-figure-fd446fd6-b19b-41f8-8d14-3a65666c0ade.jpg)

Fig.  14b.  Surfaces  of  constant  ,  and  in  toroidal  coordinates  following  the nomenclature  of  Moon-Spencer.  Note,  that  for  increasing  and  the  radii  of  the corresponding toroid and spheroid decreases (note corresponding error for toroidal coordinate directions in fig. 4.04 in Moon-Spencer [98]). For the actual magnetic coil configuration discussed here we have  0  =0.7.  In  order  to  find  the  vector  potential

inside the toroid, e.g. for the future calculation of the lepton trajectories in classical Hamiltonian equations of motion, one has accordingly to evaluate the vector potential for  0.7  (see  appendix  1).  Note,  that  in  toroidal  coordinates  the  magnetic  vector potential has only one non-vanishing component A  ≠ 0, whereas A =A =0 [101].  

## b) Helical Trajectories in a Solenoidal B-field

In order to better comprehend the character of lepton trajectories in a toroidal field we shall  recall  the  main  properties  of  the  corresponding  trajectories  in  a  solenoidal magnetic B-field.

In  a  homogeneous  solenoidal  magnetic  B-field  a  lepton  with  kinetic  energy 𝐸𝑘𝑖𝑛 = (𝛾 - 1)𝑚0 executes  a  helical  trajectory  around  the  B-field  lines  with  a  relativistic cyclotron  frequency 𝑓 = 𝑞𝐵 2𝜋𝛾𝑚0 as  is  illustrated  in  fig.  15.  This  corresponds  to  a revolution  time 𝑇 = 𝑓 -1 = 2𝜋𝛾𝑚0 𝑞𝐵 ,  a  gyroradius 𝑟 𝑔 = 𝛽⊥𝑚0𝑐 𝑞𝐵 ,  where  the  perpendicular component  of  is  used  in  calculating  rg  and  the  wavelength  of  motion 𝜆𝐿 = 𝛽∥𝑐 𝑓 = 2𝜋𝛾𝑚0𝛽∥𝑐 𝑞𝐵 . In the non-relativistic approximation  =1.0 cyclotron frequency  𝑓 𝑛𝑟 = 𝑞𝐵 2𝜋𝑚0 and revolution time 𝑇𝑛𝑟 = 2𝜋𝑚0 𝑞𝐵 are independent of the particular lepton energy. Due to its opposite  charge,  an  electron  experiences  in  a  B-field  a  Lorentz-force  equal  in magnitude but opposite in direction to the force experienced by a positron.

Fig. 15. Helical trajectory of an electron around a solenoidal B field. The helix for a positron with identical initial vector momentum has opposite sense of rotation and a common tangent with the electron's helix in the starting point.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01019v1-figure-c2dd4955-e5a4-402f-8f19-2598ef10782f.jpg)

Electrons and positrons, emitted from a single point of origin in the solenoidal B-field, execute left- and right helical motions, but around a common axis, which contains the joint  point  of origin  and is oriented parallel to the B-field lines. However, this is not sufficient  for  the  positive  identification  of  the  respective  leptons  and  subsequent determination of their primordial vector momenta when both, electrons and positrons

are present. In fig. 16 we compare the helices electrons and positrons emitted with momenta and azimuth angles (p-, 0 0  120 0 ) and (p+, 0 0  120 0 ), respectively, are executing with respect to the B-field in a homogeneous solenoidal B-field. For leptons emitted uniformly over azimuthal angles  from the target zone in a solenoidal B field, the helices of both, electrons and positrons, fill a common cylinder with a respective gyroradius  rg  centred  on  the  common  solenoid  axis;  this  prohibits  to  differentiate electrons from positrons simply by exploiting the location on the detector.

Fig.  16.  Positron  trajectories  with  momentum  p+  and  electron  trajectories  with momentum p- launched from a common origin on the axis in the jet target zone in a solenoidal B-field oriented parallel to the beam axis. Only those leptons launched with emission directions 0 0  120 0  are displayed for reasons of visibility, for 0 0  360 0 helices  would  exhibit  the  complete  overlap  in  space.  The  electrons  and  positrons execute left and right, respectively, helical trajectories and always intersect the axis after  a  completion  of  one  full  helical  turn,  i.e.  the  common  revolution  time  T.  This intersection for positrons and electrons occurs, however, at locations on the axis given by the wavelength of motion  L, according to the different velocity components of e + and e -  parallel to the beam axis.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01019v1-figure-c9b4faae-03ee-40dc-9db5-d8e9612886d5.jpg)

It  is  apparent  and  important  to  note  that  all  trajectories  of  electrons  and  positrons always return to a point intersecting the axis containing their origin after time of flight given by T, (which in its non-relativistic expression is of course independent of the energy).  The  location,  where  the  particle  returns  to  the  axis  after  the  common revolution time T, depends on the component of lepton velocity parallel to the B field

and is given by the wavelength of motion  L. For a given hit on the detector plane the distance to the axis even in combination with T is not sufficient to determine the vector momentum  of  a  lepton,  but  most  importantly,  not  even  whether  it's  a  positron  or electron.

A  solenoidal  B-field  has  seen  a  powerful  application  in  the  longitudinal  reaction microscope  [88  and  refs  therein,  104],  where  in  the  analysis  of  the  kinematically complete ionization and transfer cross sections of low energy electrons emitted in ionatom collisions the solenoidal structure of the field permits to reconstruct the entire primordial vector momentum of electrons from their TOF and the radial displacement mapping their gyroradius in the B-field onto a 2D-PSD.

## c) Electron and Positron Trajectories in a Toroidal Field

Whereas the calculation of orbits of leptons in a pure solenoidal field is straightforward  and  the  resulting  helix  a  firmly  established  powerful  tool  in  deriving kinematically complete ionization cross sections using reaction microscopes [88,104], the 1/R-dependence of the toroidal magnetic field

$$\vec { B } = \frac { A ( i ) } { R } \vec { e } _ { \varphi }$$

(e.g. as also employed as one of the confining fields in a fusion plasma) presents a formidable obstacle for calculations of charged particle dynamics and does not allow to find an analytic solution in closed form for the Newton classical equation of motion [79,99b,120,123]. The 1/R variation (see also fig. 11) of the magnitude of the B-field across the presumed helix implies that the radius of curvature on the 'inside' of the helix, i.e. closer to the center of the toroid, is smaller due to the larger B-field and is larger on the 'outside' of the helix, i.e. locations further from the centre of the toroid. This entails that the orbit of the lepton, or any charged particle, will never intersect the centre line again resulting in an increasing displacement perpendicular to the B-field lines, whose  direction  is  opposite  for  positive  and  negative  charged  particles, superimposed onto the helical motion.

Fig.  17.  Illustration  of  opposite  drift  perpendicular  to  magnetic  field 𝐵 ⃗ (with  B-field direction out of the plane) and to the gradient of the magnetic field ∇𝐵 ⃗ for positively charged positrons and for negatively charged electrons

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01019v1-figure-c5465445-ee75-449f-b933-a1cf5688c2f3.jpg)

This property of toroidal fields is the basis for its favorable disposition for lepton pair spectroscopy.

Historically AlfvØn, analysing particle trajectories in magnetically stabilized plasmas, was the first to introduce an approximation [105] based on observations of the physical properties of a collisionless magnetized plasma: here the Larmor radius characterizing the gyration is much smaller than the typical variation length scale of the generating magnetic B field. This permits to decouple the fast gyromotion from the slow motion of the centre of the confining gyromotion of the electron (called the guiding centre). This technique  has  been  used  numerous  times  in  studying  the  dynamics  of  plasma  in magnetic confinement [97,110,111], but notably also for positron spectroscopy using a toroidal spectrometer, e.g. TORI, as in transient superheavy quasimolecules created in adiabatic heavy ion -atom collisions [76,77].

However,  in  the  current  context  we  will  not  impose  such  seriously  restraining restrictions  as  for  TORI  on  digression  from  the  center  line  as  they  would  confine resulting  trajectories  to  unpractically  small  amplitudes  of  the  helical  motion  of  the lepton.

For the equations for the gyromotion and for the (slow) motion of the guiding centre across magnetic field lines in the AlfvØn approximation more powerful technique have been developed in the Hamiltonian formalism [97,106-116, 123,141 ]. The key idea is to simplify the 6D (or 8D, respectively, in relativistic context) problem and to derive after some suitable canonical transformations e.g. simple Hamiltonian equations for

the guiding centre and a gyrating 3D-motion of the charged particle presented as a Hamiltonian  in  an  effective  1D  potential  well  along  the  radial  coordinate  with  slow variation along  [96,97,110,117,118,123].

The corresponding vector potential A may either be approximated by analytical models or may be derived directly from given current density distributions, preferably in toroidal coordinates [100-102,141],but see also Appenix 1.

It  was shown that corresponding numerical evaluations of expansion coefficients in toroidal coordinates may be performed significantly faster than standard Biot-Savart techniques [99,100-102,119-122,141] which is of relevance as guidance  for further analytical analysis, e.g. for Hamilton equations of motion, is desirable.

In the present first part of our study, however, both, the B-field distribution in the toroid and solenoidal sections of the spectrometer and all corresponding lepton trajectories from production in the jet-target zone to the detector plane are calculated with OPERA3D such that also trajectories with large gyroradii can be analysed ( for results see next section, see also Appendix 1).

## 6. Invariants and other serendipitous features of lepton trajectories in a toroidal B-field

Toroidal  magnetic  spectrometers  permit  to  simultaneously  investigate  e +   and  e continua from free-free pair production emitted in heavy ion atom collisions due to the separation of leptons of opposite charge perpendicular to the centre plane of the toroid. In  the  following  we  will  illustrate  the  highly  useful  and  expedient  characteristics  of electron  and  positron  trajectories  in  toroidal  B-fields and  how  the  momentum dispersion of lepton trajectories perpendicular to the centre plane of the spectrometer permits to tune the B-field for any wanted lepton momentum window and determine the corresponding differential cross sections (DCS).

- a) 𝜃 𝑙𝑎𝑏 𝑙𝑒𝑝𝑡𝑜𝑛 = 0 0   p (lab)=0 ; invariance of momentum arcs

Fig 18a

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01019v1-figure-a264e6ec-8713-4687-acc7-031e1d54e92f.jpg)

Fig 18b

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01019v1-figure-711d058b-9873-463b-b61d-e9ba4eb223aa.jpg)

Fig.  18  a,  b.  Invariance  of  lepton  momentum  arcs  under  variation  of  the  toroidal magnetic  field.  The  momentum  dispersion  of  leptons  is  displayed  for  leptons  with laboratory momentum p =(px=0,py=0,pz) in the local coordinate system with origin at the launch  point  in  the  jet  target,  (i.e.  leptons  only  possessing  an  initial  momentum

component in the beam direction are considered). Leptons are dispersed perpendicular to the x-z-plane = bend-plane, which is a symmetry plane of the toroidal sector.  For  clarity  of  the  illustration  only  a  few  coils  of  the  toroid  including  coil  7 containing the detector plane - here positioned at  =45 0  - are shown (see fig.11 for details of coordinate frame and another option with detector plane coinciding with coil 8). For very low momenta p electrons and positrons with momentum p  0 will follow the B-field lines and both intersect the detector plane at the midpoint (indicated as O) of the coil. For a given magnetic field B, with increasing momentum of the leptons the intersection points will move away from the centre plane, the x-z- plane, and map out characteristic arcs C+ for positrons and C- for electrons. For 108.7G in the upper figure the trajectories for lepton laboratory energies up to 3 MeV (  =6.8) are shown. In the lower figure one sees that geometrically identical arcs are traced out for 1087 G and lepton energies up to around 30MeV (  =59.7). For the grid shown in the respective detector planes the separation of both horizontal and vertical lines is 50 mm for the low field and for the high field case. For the low B-field the location of the electron cusp for 5.1AGeV projectile beam energy at Elab=2.8MeV (  =6.40) is indicated on the electron arc C-; the cusp's location on the arc constitutes a tool of considerable usefulness for spectrometer calibration. B-fields given are understood as field strength on the centre ray of the toroidal bending radius  .

In  this  subsection  we  focus  on  trajectories  of  leptons  with  no  initial  transverse momentum,  ptransverse=p  =0,  and  the  formation  of  invariant  momentum  arcs  in  the detector  plane.  In  fig.  18  the  momentum  dispersion  of  electrons  and  positrons perpendicular  to  the  toroidal  bend  plane  is  illustrated  for  leptons  with  an  initial transverse  momentum  p =0.  The  intersection  of  their  trajectories  with  the  detector  plane (in figs. 18 selected to match the plane of coil 7 at  =45 0 ) traces out, starting from the midpoint O , arcs C+ and C- in the detector plane for pz increasing, for positrons and electrons, respectively; for a given B-field and the detector plane positioned at a given  the location of this intersection on the arcs uniquely determines the momentum of  a  lepton.  It  is  important  to  note  that  for  fixed  the  arcs  are  invariant  and  their geometry  does  not  change  with  a  varying  B-field;  the  geometry of  the  arcs  only depends on the toroidal radius  and the angle  ; we emphasize that w ith increasing B-field the locus of leptons of a fixed momentum will slide on the respective arc toward the midpoint O of the coil. By increasing/decreasing the B-field the intersection/locus of leptons of a fixed momentum may thus be moved to any desired location on the invariant arc toward the midpoint O of the coil or away from O . The actual geometric locations in the detector plane, i.e. the arcs C- and C+, traced out by these electrons and positrons, respectively, starting with p  =0 and p =0 at the origin  O in the detector plane, for all values of laboratory momentum p =   in the upper and lower half plane

are established by optics calculations using OPERA3D [81]. This invariance of the arcs for fixed  is an essential property of the spectrometer and will considerably facilitate the configuration of the 2D PSD detector arrays and calibration.

Apparently,  the  spatial  resolution  for  low  lepton  momenta  may  be  conveniently enhanced on the expense of total range of lepton momenta covered, either by properly tuning  the  B-field  in  the  toroid  to  an  adequate  low  field  strength  or  by  placing  the detector plane at suitable larger angles  . Low energy electrons with energies as low as 1 keV and up to around 400 keV may be conveniently studied for a magnetic field as  low  as  30G.  The  resulting  increase  of  spatial  dispersion  entails,  however,  a reduction  of  the  useful  lepton  momentum  range  covered  in  the  detector  plane.  By moving  in  a  1087G  field  the  detector  plane  from  =45 0   to  =73 0 ,  e.g.  the  vertical intercept for 13 MeV leptons emitted at 0 0  with respect to the beam axis moves from 250mm above/below the midplane of the coils at  =45 0  to a location outside the coils, theoretically near 700mm above/below the midplane, at  =73 0  .

For the desired experimental momentum calibration of the electron-arc C- one has a supremely powerful  and  very  practical  tool:  the  prominent  electron  cusp  [124,  and references therein], which contains electrons emitted with very low Ekin in the emitter frame and under near 0 0  with respect to the projectile direction and thus with velocity ve  vpar  vproj (i.e. v  =0). The location of the cusp on the arc C- for varying B-field or for different projectile energies will map out and calibrate locations corresponding to cusp lepton  momentum  (  )lepton=(  )proj  on  this  arc  C-.  The  symmetry  of  trajectories  for electrons  and  positrons  in  the  toroidal  magnetic  field  facilitates  construction  of  the homologue arc C+ for positrons as mirror image of C- with respect to the xz bend plane of the toroidal configuration.

- b) 𝜃 𝑙𝑎𝑏 𝑙𝑒𝑝𝑡𝑜𝑛 ≥ 0 0   p (lab)  0;

In the following we also consider leptons with initial transverse laboratory momentum p &gt;0 and show how the imaging properties of a toroidal field in combination with a set  of  2D  position-  and  energy  sensitive  detectors  permit  to  retrieve  all  laboratory momentum components, ptransverse and pparallel, of both leptons and thus obtain then the desired  angular  emission  characteristic  of  the  leptons  in  the  emitter  frame.  In  the remainder of this section we restrict discussions to laboratory emission angles 𝜃 𝑙𝑎𝑏 𝑙𝑒𝑝𝑡𝑜𝑛 around 20 0  and below, which in view of fig 9 amounts to a tolerable restriction of the emitter frame phase space of the leptons.

Some characteristic features of lepton trajectories described in the previous section for solenoidal fields (see also fig. 16) will be approximately conserved and will persist with some adequate modifications also in toroidal fields [see ref. 82 for case of low energy electrons in a toroidal reaction microscope]: trajectories for leptons of momenta p=  /n, n 1 and not too large forward emission angles mostly not extending beyond 20 0  with  respect to the beam axis, i.e. p  (lab)  0, will intersect at or very near a focus Xn when their  time  of  flight  (TOF)  T=Sj/  parc  (Sj  =arc  travelled  by  lepton  with  perpendicular momentum pperp=0 to Xn in detector plane j) from the origin at the jet target to the focus Xn (determined by their longitudinal momentum ppar) is a multiple of the inverse of the effective cyclotron frequency f= qB/2  m 0 for the B field given:

$$T = \frac { S } { \beta _ { p a r ^ { c } } } = n f ^ { - 1 } = n \left ( \frac { q B } { 2 \pi \gamma m _ { 0 } } \right ) ^ { - 1 } = n \frac { 2 \pi \gamma m _ { 0 } } { q B }$$

It is then straightforward - by tuning the B-field - to position the anticipated focus Xn for a well-defined lepton momentum pn on the respective arc in a chosen detector plane (recall that the geometry of an arc is only depending on the detector plane location and the radius of curvature  of the toroidal field, but is invariant under variation of the Bfield).

In fig 19 we show trajectories in a 108.7G field for leptons with two different momenta in the laboratory,  =1.2/n, and 1  n 2, respectively, for laboratory emission angles up  to  L=20 0  which intersect for all azimuthal angles in four distinct foci on the arcs C- and C+.

The usefulness of the relation given in eq. 9 is, of course, contingent on the assumption that only trajectories in the confidence volume of the magnetic field B is considered (recall that following eq. 7 the B-field inside the toroid strictly varies with 1/R). This implies, for the foci to be useful (i.e. the foci exhibiting a negligible astigmatism and thus reasonably small extension) a restriction to excursions perpendicular to the centre of curvature  of the toroid (see fig 11) not extending into the fringes of the coils, i.e. a moderate  gyro-radius  rg  of  the  lepton,  is  imposed.  This  assumption  is  met  in  the example illustrated in fig 20.

Fig.  19.  Electrons  and  positrons  launched  at  the  supersonic  jet  target  (Cartesian coordinates (1000, 0, -50), see fig. 11) with polar angles  L=20 0  and azimuthal angles 0 0  360 0   with  respect  to  the  projectile  axis  =  local  z-axis  at  the  jet.  In  the  given configuration  with  B=108.7  G  and  the  detector  plane  located  at  66 0   all  lepton trajectories with  =1.2/n, n=1, 2, are intersecting at well distinguishable locations Xn in the detector plane, where an array of 2D PSDs may be appropriately configured, centred  on  the  arcs  C+  and  C-.  The  lepton  kinetic  energies  corresponding  to  the trajectories given here are 85 keV (  =1.2/2) and 287 keV (  =1.2/1), respectively. The detector plane has grid-tiles of 100mm x 100mm. Note, that trajectories for different azimuth  (in coordinates with origin at jet) exhibit the completion of one full circle upon intersection at the first focus following eq. 9 for leptons with 287 keV. For electrons with  =1.2/2 leptons have completed two full circle upon intersection at the detector plane; the option of an intercept of electrons at a suitable intermediate focus by a small disc or even a 2D PSD is illustrated for the lower energy electrons with Elab=85keV, i.e.  =0.6.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01019v1-figure-1028c23a-eb16-416a-8dbe-93781d88cc58.jpg)

For any B-field selected these foci Xn , n  1, in the detector plane are located along the arcs  C+  and  C-  according  to  their  momentum  pn=(  )/n  shown  in  fig  18.  Lepton trajectories for leptons with pn=(  )/n, n&gt;1, contain (n-1) intermediate foci before hitting the nth focus in the detector plane at Xn. With increasing n, i.e. decreasing momenta, the spatial separation of foci Xn located on the arcs C+ and C- decreases rapidly and foci  Xn  converge  towards O ;  this  implies  that  the  entire  backward  emission  of  fast leptons (as seen in the projectile frame, see fig. 9) out of the fast projectile will thus virtually collapse to one-dimensional arc sectors in the detector plane either side of O .

The corresponding single differential cross section 𝑑𝜎 𝑑𝑝 for  lepton production mapped onto a 2D PSD may thus be conveniently recovered. It follows that the coincident DDCSs 𝑑 2 𝜎 𝑑𝑝𝑑Ω for  free-free  pair  production  where  both  leptons  are  emitted  into  the backward  direction  (in  the  emitter  frame)  is  similarly  acquired  in  a  straightforward fashion  via  coincident  detection  of  electrons  and  positrons  in  the  corresponding detectors covering the arcs close to O . (In passing we note that from fig 9a one sees that under 𝜃𝑙𝑎𝑏 𝑖𝑛 𝑡ℎ𝑒 𝑣𝑖𝑐𝑖𝑛𝑖𝑡𝑦 𝑜𝑓 0 0 for decreasing laboratory lepton energy/momentum only leptons with increasing emitter frame energy are contributing).

For free-free lepton pairs with backward-forward emission pattern in the emitter frame correspondingly the coincident laboratory frame DDCS are derived from coincidences between a lepton on the arc section close to O and the coincident second lepton of the pair at individual foci Xn, n  1.

Leptons  with  momenta  p&gt;p1  will  not  exhibit  an  interception  focus  independent  of perpendicular momenta in  the detector  plane  (see  also  fig  26);  this  enables  direct determination of full laboratory DCS for p&gt;p1.

The apparent near convergence onto a focus Xn after the time of flight TOF = n f -1 , with f the cyclotron frequency, is obviously facilitated and its usefulness augmented by the considerable  kinematic  restriction  to  small  laboratory  angles  lab arising  at  large collision energies as provided by the HESR, for which ppar and thus the TOF only exhibit a weak dependence upon  lab (see also fig. 9). However, for large laboratory transverse momenta, i.e. for gyro-radii 𝑟 𝑔 = 𝑚𝛾𝛽⊥𝑐 |𝑒|𝐵 corresponding to lepton trajectories sampling regions with B-fields noticeably larger (at R&lt;R0) and smaller (at R&gt;R0) than the B-field at the centre radius R0 of the toroid, the inverse cyclotron frequency varies sufficiently as to generate azimuth-dependent interceptions [82] instead of single nearly azimuthindependent foci Xn as seen for the small gyro-radii rg.

This salient feature of lepton trajectories in a toroidal magnetic field - appearance of periodic  intermediate  foci  along  trajectories  to  the  detector  plane  for  not  too  large gyroradii - may also be employed to intercept electrons and positrons of any desired momentum p&lt;p  at suitable locations before the detector plane with high efficiency by 1 small  diameter  detectors  or  beam  stops,  with  minor  effects  on  leptons  of  other momenta due to the strong spatial dispersion. This permits to introduce for particular

combinations  of  lepton  energies  highly  efficient  multiple  detector  planes  (see  e.g. single intercept in fig. 19), as also the location of the intermediate foci are space fixed and do not depend on the magnetic field B ( for details see section d).

Figures 20 show for a high B-field of 1087 G the intersection in the detector plane of lepton trajectories over a wide range of laboratory energies and two sets of emission angles. This corresponds to covering emitter frame energies up to 500keV and 1MeV (which is the mean lepton energy Eav predicted by Scheid et al. for a similar collision velocity) and even above Eav for the highest laboratory energies (see also figs. 9 and 10).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01019v1-figure-d470fccd-9de1-4b92-bff9-734545cae3c6.jpg)

Fig. 20. Lepton trajectories in the toroidal B-field =1087G. For the laboratory lepton energies  selected  here  the  range  of  projectile  frame  angles  and  laboratory  frame angles are given. The arc C+ and C- (see fig. 18) along which leptons are dispersed are indicated as dotted lines. The tile dimensions are 100mm x100mm. For illustration only lepton trajectories with pn=  /n for  =12.7 and [n=1, 4] are displayed besides two cases with higher lepton momenta. Intercepts in the detector plane for pn with n  2 then easily produce laboratory differential cross sections for backwards emission (as seen in the emitter frame, and summed over a wide range of emitter frame energies). The spatial  overlap  of  lepton  intercepts  for  the  larger  lab will  not  confuse  the  correct identification of lepton energy as the 2D-PSD detector may identify the lepton energy due to energy resolution down to  E/E 1%. 

For completeness we add that the entire very-low energy electron emission out of the projectile (Ekin  1-2keV, again as seen in the projectile frame, in contrast to high Ekin around 500keV to 2000keV as in fig. 9) is well known to be entirely contained in the electron  cusp  [124  and  refs.  therein],  which  for  5.1AGeV  specific  projectile  energy

Ebeam  is  found  at  an  electron  laboratory  kinetic  energy  of  Ecusp  =Ebeam/1822.8878= 2.8MeV  at  lab  0 0 (see  also  the  cusp  in  fig  9  as  the  point  of  convergence  in  the laboratory for emitter frame energy Eemitter →0). The low energy electrons in the cusp have been shown to mostly arise over a range of close to more peripheral collisions (i.e.  with  larger  internuclear  separations)  from  radiative  and  non-radiative  electron capture into the projectile continuum, RECC and ECC, respectively and electron loss to continuum, ELC, for non-bare projectiles [124].

In the domain of very low energy (emitter frame) leptons besides above mentioned cusp electrons also positrons may appear, originating in free-free as well as in boundfree pair creation with very-low kinetic energy in the emitter frame; they will then appear in  a  corresponding  mirror  location  on  the  C+  arc.  Positrons  from  bound-free  pair creation (eq. 2a) will be unambiguously identified via coincidences with a projectile with an outgoing charge reduced by one unit from the incoming charge state and thus are clearly distinguished from free-free lepton pair creation (eq. 1a) and their identification via electron-positron coincidences.

Any  electron  emission  from  direct  target  ionization  is,  of  course,  not  subject  to kinematic restrictions concerning transverse momenta discussed above; the corresponding low energy electron distribution will be rather isotropic and their energy distribution is a continuum with a near exponential fall-off above zero.

- c) General momentum/energy calibration of lepton trajectory intersections in the detector plane, i.e. including p  (lab)  0.

In the general case, i.e. the absence of any kinematic restrictions for coincident lepton pairs (as e.g. back to back emission in the emitter frame), a unique attribution of the primordial  vector  momentum  of  the  lepton  in  the  emitter  frame  (the  emitter  frame characterized by 𝛽0 and 𝛾0 , a lepton in the emitter frame by p',  ',  ')  using

$$( p ^ { \prime } c ) ^ { 2 } = ( p c ) ^ { 2 } \sin ^ { 2 } \theta _ { l a b } + \gamma _ { 0 } ^ { 2 } \Big ( \text{pc} \cos \theta _ { l a b } - \beta _ { 0 } \sqrt { ( p c ) ^ { 2 } + m ^ { 2 } c ^ { 4 } } \Big ) ^ { 2 } \, \Big ( 1 0 )$$

can only be provided when both, the lepton's laboratory momentum p or energy Elab and the laboratory observation angle  lab, are simultaneously determined, as follows from fig. 9. The corresponding equation relating laboratory frame and emitter frame angles is

$$c o s \theta ^ { \prime } = \frac { - \frac { \beta _ { 0 } } { \beta } \gamma _ { 0 } ^ { 2 } s i n ^ { 2 } \theta _ { l a b } \pm \sqrt { \gamma _ { 0 } ^ { 2 } ( 1 - ( \frac { \beta _ { 0 } } { \beta } ) ^ { 2 } ) s i n ^ { 2 } \theta _ { l a b } + c o s ^ { 2 } \theta _ { l a b } } \ c o s \theta _ { l a b } } { \gamma _ { 0 } ^ { 2 } s i n ^ { 2 } \theta _ { l a b } + c o s ^ { 2 } \theta _ { l a b } }$$

where 𝛽 ′ corresponds to the momentum of the lepton in the moving frame.

Instead of directly measuring  lab  of  a  lepton,  however, it is far more convenient to experimentally  determine  its  gyro-radius  rg  in  the  toroidal  B-field  and  subsequently calculate its laboratory transverse momentum.  The  relation  between  the  true gyroradius rg and the effective rg-eff obtained from the location of the intercept in the detector  plane  is,  however,  best  established  in  the  spectrometer  calibration  using conversion electron lines, because only for locations outside the foci Xn on the arcs C+ and  C-  in  the  detector  plane  one  may  eventually  obtain  rg  and  the  transverse momentum and its azimuth from rg-eff (see below).

The first step of the experimental momentum/energy calibration of the arcs C+ and Cuses the narrow electron cusp with (𝛽𝛾) 𝐶𝑢𝑠𝑝 = (𝛽𝛾) 𝑏𝑒𝑎𝑚 and energy Ecusp =Ebeam/1822.8878 correspondingly ; ( for cusp electrons  lab=0 0   p (lab)=0 to a very good approximation) at varying beam energies (as indicated in the previous sections). Additionally  conversion  electron  lines  from  radioactive  sources  equipped  with  an aperture  for  lab=0 0 will be  employed.  We  may  thus  attribute  for  any  suitable combination of B-field and projectile energy to every point XE on the arc C- an electron energy E and electron momentum p, as illustrated in figs 18.

For electron transverse momenta p (lab)   0 intersections on the detector plane off the arc C- appear; here suitable conversion electron lines from radioactive sources (now optionally equipped with selected angle defining apertures for distinct  lab˃0 0 resulting in well-defined p  (lab)  0 for every angle  lab ) are employed  (see figs. 21, 22).

Fig. 21. Momentum calibration with a  207 Bi source in 108.7G toroidal magnetic field. Energy/momentum calibration  of  the  toroidal  magnetic  spectrometer  using  Pb  KLL Auger and conversion electron lines emitted by a  207 Bi source and detected at 66 0  for a B-field setting of 108.7G. Besides the Pb KLL Auger group centred at 56.7 keV we show in this illustration only the strongest components of the conversion electron line groups; besides the strong 481.7 KeV line the 553.9 keV line with  25% intensity of the former and the 975.7 keV line with the 1047.8 keV line which also has around a quarter the intensity of the former. The emitter source situated beneath the jet in the target zone is furnished with an angle defining aperture, which can be suitably adjusted to  angles  between  0  and  20  degree,  here  an  angle  of  5 0   is  used  in  the  trajectory calculations. The use of various angles on the source serves to calibrate the effective gyro-radius seen in the detector plane. The combination of  1% energy resolution and 2D  position  sensitive  detection  of  the  electron  accomplishes  in  an  experiment  the unambiguous identification  where  electrons  of  different  energies  may  hit  the  same location on the detector.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01019v1-figure-69f6b8dd-cc36-41ef-9251-a3abe044deef.jpg)

For the general gyro-radius/momentum calibration, i.e. intercepts in the detector plane off the arc C-, it is convenient to use angle defining apertures, e.g. 0 0 , 5 0  and 15 0  with respect  to  the  beam  axis,  for  a 207 Bi  conversion  electron  source  over  a  range  of different  magnetic  B-fields  in  the  toroid.  This  permits  to  calibrate  the  apparent gyroradius  rg-eff for the respective transverse  momentum  component  (and  the corresponding true gyro-radius rg) selected by the angle setting in the aperture, e.g. of the prevailing 975 keV conversion electron line over a range of magnetic fields.  In fig. 22 we show the mapping of the 975keV conversion electron line for a range of B-fields.

For  165  G  we  illustrate  in  fig.  22  how  using  apertures  for  two  different  angles  the procedure for momentum calibration of the locations of the intercept in the detector plane is accomplished. It can also be nicely seen that for 234G all electrons from the 975.7 keV line emerging in the forward direction with  5 0 exhibit a focus after 66 0 , i.e. at coil 8 of the toroid, in one location on the arc C-.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01019v1-figure-413382f3-840e-4b7d-ac24-75a6a1316a0c.jpg)

Fig. 22. Mapping of the 975.7 keV conversion electron line from  207 Bi onto the detector plane at 66 0  and toroidal magnetic B-field for fields ranging from 125 G top 935 G. The 207 Bi  source  located  in  the  target  zone  beneath  the  supersonic  jet  target  is  in  this OPERA  simulation  equipped  with  angle  defining  apertures  to  better  illustrate  the mapping used for calibration. It is apparent that for the electron energy selected, 975 keV, the B-field may either be tuned to focus all electrons onto one focal spot in the detector plane, as here shown for 234G, or to any lower value, e.g. here 165 G, where the effective gyro-radius on the detector plane can be related to the transverse electron momentum via the known laboratory emission angle in the calibration source.

On the path to obtain differential cross sections in the emitter frame, it is instructive to look at leptons emitted with a fixed emitter frame energy, but over a wide range of emitter frame angles; they will appear over a wide range of energies and angles in the laboratory frame; this results in very useful allowed 'bands' on either side of the arcs as is illustrated in fig. 23a,b.

Figs. 23a,b exhibit for 2 magnetic field settings, 272G and 1087G, for an emitter frame kinetic energy of 0.2MeV (  =0.96) and 1MeV (  =2.8), and 0.1, 0.2, 0.5, 1 and 2 MeV, respectively, for emitted leptons the regions covered in the detector plane around the

arcs C+ and C- for leptons of the corresponding emitter frame energies. It is apparent that 2D PSD detectors do not have to cover the entire cross section of the coil in the detector plane but that the leptons are restricted to a reasonably confined area either side of the arcs in the detector plane, e.g. as shown at  =45 0  .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01019v1-figure-97c1fb04-470d-4070-8147-693d87a0f230.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01019v1-figure-c89fd795-f2aa-4e21-bb47-8eb970fba251.jpg)

Fig. 23a,b. Kinematically allowed regions in the detector plane for leptons emitted by 5.1AGeV projectiles for two B-field settings of the toroidal magnetic spectrometer with the detector plane at  =45 0 . The top figure displays for B=272G the regions adjacent to  the  arcs  C+ and C- where leptons with 0.2MeV and 1 MeV kinetic energy in the emitter frame can be detected; in the lepton trajectory calculation for the low B field only laboratory energies up to 7 MeV were used, for EEmitter =1MeV this corresponds to

p (lab)   0 . In the lower figure the corresponding regions for leptons with 0.1, 0.2, 0.5, 1 and 2MeV are displayed for a magnetic field of 1087 G.

The areas either side of the arcs exhibit how possible locations of the intersections of lepton trajectories with the detector plane are confined. In these 'allowed regions' now a fine structure evolves superimposed on these areas in accordance with the foci in the mapping of lepton trajectories in the toroidal B-field, following eq. 9 (also figs. 1922). This turns out to be extremely useful for obtaining differential cross sections and will be discussed in the following.

For electrons detected at a location P in the lower half plane on the 2D PSD detector its total kinetic energy El measured attributes uniquely a corresponding location XE on the arc C-, respectively, of an electron with kinetic energy El, but vanishing transverse momentum. The distance 𝑃𝑋𝐸 ̅ ̅ ̅ ̅ ̅ ̅ of  the  actual  location  P  of  the  hit  to  the  location  XE determines  the apparent gyro-radius 𝑟 𝑔-𝑒𝑓𝑓 = 𝑃𝑋𝐸. ̅ ̅ ̅ ̅ ̅ ̅ For leptons with momenta approximately fulfilling eq. 9, i.e. the intersection of their trajectories in the detector plane converging on a focus, obviously the distance 𝑃𝑋𝐸 ̅ ̅ ̅ ̅ ̅ ̅ will be vanishingly small. For other locations 𝑃𝑋𝐸 ̅ ̅ ̅ ̅ ̅ ̅ will underestimate the true gyro-radius rg of a lepton by a location dependent factor (see fig 24). For this reason for a given magnetic field B the true gyroradius rg has to be derived always from calibrations with conversion electron lines from the  measured  apparent  gyro-radius  rg-eff.  The  momentum  calibration  found  for electrons can then be applied in analogy for positions of positrons in the upper half plane.

The  true  gyro-radius  rg=(   ) m 0 c/qB  then  determines  p =(    ) of  the  lepton,  as illustrated  in  fig.  24.  El  and  p  now uniquely  determine  for  the  detected  lepton  the emitter frame energy (see fig. 10b) and also its emission direction in the emitter frame (fig.9).

Due to the invariance of arcs C+ and C- at a given toroidal angle  , all leptons with pperp &gt;0 will intercept the detector plane in well-defined areas on either side of the arcs C+  and  C-  as  illustrated  in  fig  23.  This  will  considerably  facilitate  design  of  the appropriate 2D PSDs detector frames for lepton detection.

In figs 24 we illustrate the relation between apparent gyro-radius and the transverse laboratory lepton momentum for two B-field settings, B=1087G and B=272 G.

Fig. 24 a b. Transverse lepton momentum ptransverse  and true and apparent gyro-radii  vs. laboratory momentum  for four emitter frame lepton energies 200, 500, 1000 and 2000 keV at emitter specific energy of 5.1 AGeV for two magnetic fields, 1087G and 272 G. The red curves represent the relation between a lepton's transverse momentum

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01019v1-figure-5741b987-d40b-4cb3-9e99-0d9305a0bc20.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01019v1-figure-b1f150bc-839a-4876-8e7a-65872a5a3156.jpg)

and total laboratory momentum for the four selected emitter frame lepton energies. The  vertical  blue  lines  indicate  those  laboratory  lepton  momenta  (  ),  for  which according to eq. 9 at the two magnetic fields B=272G and B=1087G foci are generated at the detector plane, i.e. all trajectories of this momentum will intercept at one locus in  the  detector  plane  (here  at  =66 0 )  with  an  apparent  gyro-radius=0,  thus  not permitting determination of  lab or transverse momenta; numbers above the blue lines give the number of nodes completed when the trajectory crosses the detector plane: n=1 means first focus in detector plane, n=2 means one intermediate focus: (n-1) gives the number of intermediate foci (for momentum pn) before hitting the detector plane (see also fig. 19). The ordinate on the right in the two figures gives the particular gyroradius for the leptons in the two respective B-fields using  [m]= 1.7045 10 -3   (  )perp /B[T], e.g.  [cm]=6.27 (  )perp for B=272G. Vertical lines for lepton momenta pn, n&gt;6 will appear with increasing density and have thus been omitted for simplicity of the figure. Analogously for a field B=272G the group of lines is shifted to lower momenta on the abscissa, the first node appears at  =3.15 and all subsequent nodes correspondingly at lower momenta  /n, n  2. Red lines represent the true gyro-radii and the green lines between  the  nodes  represent  the  apparent  gyro-radii  rg-eff  to  be  determined  by calibrations. Importantly, the geometric location of the nodes in the detector plane is invariant under change of B-field ( see below).

In figs 24 we illustrate for a collision energy of 5.1AGeV and four emitter frame lepton energies between 0.2 and 2MeV the laboratory transverse momentum of the lepton as function of laboratory momentum  and corresponding gyro-radii in the two B-fields which  is  in  a  satisfactory  approximation  calculated  using  here  the  equation  for  a homogenous B-field: the true gyro-radius rg for a lepton with momentum perpendicular to the field follows from 𝑟 [𝑚] = 1.7045 10 -3 (𝛽𝛾) 𝑔 𝑡𝑟𝑎𝑛𝑠𝑣 /𝐵[𝑇] . The toroidal nature of the field,  however,  produces  intermediate  foci  along  the  path  to  the  detector  plane  as outlined above. For a 1087G toroidal B-field and the detector plane at 66 0  these foci appear at momenta indicated by six vertical lines in fig. 24 (similarly for the 272G field the  lines  appear  at  corresponding  lower  momenta).  At  these  momenta  leptons  will intersect almost independent of their transverse momentum ptransverse at one focal point on the detector plane. For these momenta the transverse momentum is integrated over and no determination of ptransverse is possible, as the apparent gyroradius rg-eff is 0. For momenta between the momenta pn-1, pn of the foci a determination of true ptransverse and true gyroradii via measured apparent gyro-radii is possible as outlined before; this is indicated by the green arcs centred on the foci at the blue vertical lines in figs. 24a,b For lepton momenta situated approximately between two foci a measured apparent gyro-radius (to be imagined as situated on a virtual green arc) needs to be corrected accordingly  using  the  calibration  to  derive  a  true  gyro-radius  rg  and  the  correct

transverse momentum, on the corresponding red line; - this now allows to determine the emitter frame energy and emission direction of the lepton.

In the calibration with a conversion electron source the angle defining aperture on the source determines for every electron line (e.g. Pb KLL-Auger and 500keV or 1000 keV conversion  electron  lines  from 207 Bi)  of  momentum  (  )  the  transverse  momentum component (  )transverse, and thus a corresponding 'red' point in a diagram analogue to those in figs 24 for a suitably chosen toroidal B field. The actually measured 'green' point in the calibration determines for the momentum (  ) on the abscissa an apparent gyro-radius rg-eff and thus a correction factor. This is to be executed for various angles of the source and all suitable conversion lines, i.e. lepton energies on the abscissa. Thus correction factors are obtained to derive for measured hits in the detector plane the true (  )transverse at this abscissa.

For lower B-fields leptons with correspondingly lower momenta will be focussed, but onto exactly the same geometric locations Xn, n  1 (see below). (This is of considerable usefulness  for  designing  the  2D-PSD  detector  configuration  and  the  calibration procedures, as well as for the subsequent derivation of emitter frame DCS).

A full identification of a lepton's primordial vector momentum in the emitter frame just requires the determination of the lepton's laboratory energy E and of the hit location P in the detector plane in order to determine ptransverse. This is the basis for the next step, determining the differential production cross sections in the emitter frame.

## d) geometric invariance of location of foci on arcs C+ and C-

As outlined already above, for a given B-field in the toroid the collision velocity and emitter frame lepton energy and -emission direction ( i.e. lepton vector momentum) determine the kinematically allowed regions adjacent to the arcs C+ and C-; they are highlighted for selected emitter frame lepton energies in figs 23a,b. These regions will, however, not be uniformly illuminated by the leptons due to the characteristics of lepton trajectories in a magnetic toroidal field.

The  near  helical  character  of  the  lepton's  orbits  in  the  toroidal  field  and  the  axial intercepts at well determined momenta pn=  /n, n  1,resulting in a series of foci Xn on the arcs C+ and C-, will modify the mapping of each lepton's phase space onto the 2D-

PSDs, as illustrated in figs. 24 and 25. When comparing the location of the foci on the arcs C+ and C- for 3 different B-fields, we recognize some highly significant features:

- i) as exhibited above( fig 18a,b) for leptons with pperp=0 the geometry of the arcs C+ and C- for a given toroidal radius  and angle  is invariant under the change of the toroidal B-field,
- ii) with  increasing  B-field  leptons  of  correspondingly  increased  momenta  are focussed,
- iii) and most importantly, also the geometric locations of the foci Xn on the arcs C+ and C- are invariant under a change of toroidal B-field. The momentum of these leptons at location Xi increases proportional to the B-field.

B=109 G

B=272 G

B=1090 G

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01019v1-figure-3646b4d3-353d-4270-b4dc-5a8d3583fcd8.jpg)

Fig. 25. Invariance of geometric location of foci under variation of the toroidal magnetic field, illustrated for foci Xi, 1   i 3, on the arc C+ and C- at fixed position of the detector plane  at  =66 0 .  For  three  B-fields  the  lepton  energies  at  these  invariant  locations, which exhibit a focus at Xi, are shown. The momentum of these leptons at location Xi increases proportional to the B-field, i.e. leptons at X1 for B=1090G have a momentum  ten times that for leptons at X1 and B=109 G.

It is seen immediately that the invariance of both, the geometric location of arcs and of the foci on these arcs under variation of B, introduces significant simplifications in the configuration of  the 2D-PSD lepton detectors in the detector plane; this apparently entails simplifications in the calibration procedures as well.

For  momenta well  above  p1  the  determination  of  the  DDCS  is  not  affected  by  the characteristics of lepton trajectories as long as the azimuthal twist of the trajectory from origin  to  the  detector  plane,  as  illustrated  in  fig.  19  and  fig.  20,  is  taken  into consideration. This azimuthal twist will only affect the location of the plane of emission but not the determination of the gyro radius rg.

From figs 24 and 9 we note that for large laboratory momenta of leptons below the kinematic maximum only small transverse momenta contribute. E.g., for 1087G the largest perpendicular excursion a lepton with emitter frame energy 1MeV will have to the arcs C + /C in the detector plane is 46mm corresponding to near 90 0  emission in the emitter  frame.  For  locations  between  the  momentum  foci  on  the  arcs  C + /C one samples for increasing n of the foci more and more  backward emission in the emitter frame with decreasing fraction of the DDCS.

It is interesting to note that due to the kinematics of a fast moving emitter frame large transverse momenta can only appear for small laboratory emission angles  lab; large angels  lab are related via the kinematics with very small laboratory momenta. Lepton momenta will thus trace out a narrow region left and right of C+ or C- and only at large laboratory lepton momenta a broad region on either side of C+ and C- is kinematically possible.

e) Differential cross sections in emitter frame When the 2D PSD detector arrays are suitably distributed along the arcs in the detector plane they selectively cover foci and intermediate regions as needed and so permit to determine unconditional double differential cross sections (DDCS) 𝑑 2 𝜎 𝑑𝑝𝑑Ω for electron and  positron  production.  When  one  applies  suitable  coincidence  conditions  for detection of electron-positron pairs the triple and partially the full quadruple differential cross sections 3DCS and 4DCS for free-free pair production may then be derived ( e.g. 3𝐷𝐶𝑆(𝑒 + ) = 𝑑 3 𝜎 𝑑𝑝+𝑑Ω+𝑑𝑝-)  .

As electrons and positrons emitted as lepton pairs are not kinematically coupled, i.e. the relation between emitter frame coordinates and laboratory frame coordinates of one  lepton  of  the  pair  is  entirely  independent  of  the  emitter  and  laboratory  frame coordinates  of  the  other  lepton,  the  Jacobi  matrix  for  transformation  of  the  cross sections form  the  emitter  into  the  laboratory  frame  is  diagonal  in  the  leptons.  This means that for calculating the essential Jacobian of the Lorentz transformation we may focus on the DDCS.

For the double differential lepton production cross sections the desired emitter frame cross section (primed variables) is derived from the measured laboratory frame cross section by evaluation the Jacobian of the corresponding Lorentz transformation.  K. Dedrick  reports  [132]  that  Wolfgang  Panofsky  was  the  first  to  give  in  high  energy

physics the transformation for a laboratory DDCS 𝑑 2 𝜎 𝑑𝐸𝑑Ω

to provide the corresponding emitter frame DDCS 𝑑 2 𝜎 ′ 𝑑𝐸 ′ 𝑑Ω ′ ; he showed that evaluation of the Jacobian

$$\begin{vmatrix} \frac { \partial E ^ { \prime } } { \partial E } \Big | _ { \cos \theta } & \frac { \partial \cos \theta ^ { \prime } } { \partial E } \Big | _ { \cos \theta } \\ \frac { \partial E ^ { \prime } } { \partial \cos \theta } \Big | _ { E } & \frac { \partial \cos \theta ^ { \prime } } { \partial \cos \theta } \Big | _ { E } \end{vmatrix} \text{leads to the equation (see also appendix 2)}$$

$$\frac { d ^ { 2 } \sigma ^ { \prime } } { d E ^ { \prime } d \Omega ^ { \prime } } = \frac { p } { p ^ { \prime } } \frac { d ^ { 2 } \sigma } { d E d \Omega } \quad,$$

a relation of astonishing simplicity.

In the current study of pair production in relativistic collision system using a magnetic toroidal spectrometer the momentum space provides the more suitable frame for the differential  cross  sections.  Once  more,  after  some  algebra  the  determinant  for  the corresponding Jacobian matrix

$$\begin{vmatrix} \cdot \\ \left | \frac { \partial p ^ { \prime } } { \partial p } \right | _ { c o s \theta } \\ \frac { \partial p ^ { \prime } } { \partial c o s \theta } \right | _ { p } \\ \cdot \end{vmatrix}$$

$$\begin{array} { c c c } \cdot & \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
 \cdot & \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
 \cdot & \cdot \\ \cdot & \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \\ \cdot & \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \\ \cdot & \quad \\ \cdot & \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \

\cdot & \quad \\ \cdot & \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \\ \cdot & \quad \\ \cdot & & \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \\ \cdot & & \quad \\ \cdot & & \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \\ \cdot & & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cd.. & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad\\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \\ \cd. & & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \ \ \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdocks & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad
 \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cd. & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad \\ \cdots & \quad } \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

transformation the DDCS in the emitter frame as

$$\frac { d ^ { 2 } \sigma ^ { \prime } ( p ^ { \prime }, \Omega ^ { \prime } ) } { d p ^ { \prime } d \Omega ^ { \prime } } = \frac { \gamma } { \gamma ^ { \prime } } \frac { ( p ^ { \prime } ) ^ { 2 } } { p ^ { 2 } } \frac { d ^ { 2 } \sigma } { d p d \Omega } \quad,$$

where again, as in the derivation of eq. 12, only the well-known relations for the Lorentz transformation of energy and momentum are used. This surprisingly simple relation was first given by Hagedorn [142] and Dedrick [132]. In the appendix we rederive the explicit  expressions  for  all  individual  matrix  elements  of  both  Jacobians  used  in derivation of eqs. 12 and 13, because these individual terms are considerably more elaborate than the resulting products appearing in eqs 12 and 13. For the 4DCS then one has the product of the Jacobian determinants for electron and positron.

- f) Recoil ions and  electrons  in a magnetic  toroidal spectrometer,  impact parameter of free-free pair production

The observation that a toroidal magnetic field separates particles of opposite charge transverse to the bend plane of the spectrometer as illustrated above for electrons and positrons, immediately encourages  curiosity as to the suitability of this field configuration for coincident detection of electrons and recoil ions produced in a heavy-

ion atom collision. It is evident that besides the different mass the emission direction and momentum have to be considered. In the following we use a symmetric very heavy collision system at adiabatic and near relativistic collision velocities to demonstrate the imaging properties for simultaneous mapping of lepton and recoil ions.

Fig 26. Simultaneous emission of electrons (blue and green tracks extending from the jet towards the left) and selected recoil ions with emission angels 89.97 0 /89.98 0  with azimuth  0 0   and180  energy  energies  20  eV  (yellow  trajectories).  Note  reflection  of recoils in the B=272G field for 20 eV and azimuth 0 0 , but not for azimuth 180 0 . Red trajectories correspond to recoil energies 1000eV with the same emission angles as for 20 eV.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01019v1-figure-f93e6f27-131d-4d94-97f7-3b98a1314317.jpg)

Theory (28, 29) using coupled channel calculations predicts &lt;b&gt;  100fm, i.e. E kin = 0.1AkeV for recoiling target ions ( see also figs. 27 and 28). For B=270 G we find  = 35 cm for Au 40+ recoils extracted onto 2D PSDs antiparallel to v proj .  Bhabha [24, 25] has  arrived  at  a  similar  result  for  the  relevant  impact  parameters  for  lepton  pair creation.

Fig 27. Polar plot for impact parameter and recoil ion energy as function of laboratory recoil  emission  angle  for  5  AGeV  Au 79+  +  Au.  Even  for  collisions  predicted  to  be dominantly contributing  to  free-free  pair  production  the  recoil  ions  are  emitted  with angles very close to 90 0  with respect to the beam axis.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01019v1-figure-8f144efc-83f7-453b-8d6b-9419d88ae79e.jpg)

More than the high mass, it is the extreme polar emission angle (corresponding to the relevant impact parameters) with respect to the magnetic field direction, which affects the mapping of target recoiling from the collisions. From fig. 26 it is apparent that the position dependent detection of the recoil ions moving in the toroidal field opens the door to extreme sensitive determination of the azimuth of the recoil emission angle. We will  analyse  the  potential  for  very  high  resolution  spectroscopy  employing  the extreme azimuthal sensitivity of the spectrometer in a forthcoming paper.

Fig. 28 Recoil ion energy as function of impact parameter for 5AGeV Au 79+  + Au. Note the specific recoil energy of around 50AeV for the predicted most probable impact parameter 100fm facilitating the determination of impact parameters via recoil TOF.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01019v1-figure-2474e848-b049-48ec-94d1-e287de208492.jpg)

## 7. Summary

Free-free  pair  production  cross  sections  in  near  relativistic  ion-atom  collisions  are predicted to produce lepton pairs over a wide range of momenta and relative emission directions-  contrary  to  photon  induced  case.  Due  to  considerable  computational requirements theory shows only partial DDCS, but no fully 4DCS 𝑑 4 𝜎 𝑑𝐸 𝑒 +𝑑𝛺 .

𝑒 +𝑑𝐸 𝑒 -𝑑𝛺 𝑒 -

We present an investigation of lepton trajectories in toroidal magnetic fields and we analyze  a  large  acceptance  magnetic  toroidal  sector  spectrometer  for  coincident studies of electrons and positrons from ion-atom collisions .

The toroidal magnetic configuration offers the advantage to image simultaneously the vector momenta of electron and positron of a free-free lepton pair with near 4  solid 𝑑 3 and angle.  This  opens  the  door  to  selected  benchmarks  for  3DCS 𝜎 𝑑𝐸 𝑒 +𝑑𝛺 𝑒 +𝑑𝐸𝑒 -

$$\frac { d ^ { 3 } \sigma } { d E _ { e ^ { - } d \Omega _ { - } d E _ { e ^ { + } } } }, \, \text{at this time there are no experimental 4DCS } \frac { d ^ { 4 } \sigma } { d E _ { e ^ { + } d \Omega _ { e ^ { + } } d E _ { e ^ { - } } d \Omega _ { e ^ { - } } } } \, \text{for stringent}$$

test of theories.

A preliminary explorative analysis of recoil trajectories in the toroidal field indicates that the  dynamics  of  free-free  pair  production  may  be  accessible  via  kinematic  coincidences of recoiling target ions with lepton pairs.

## Appendix 1: The vector potential in toroidal coordinates

Interest  in  the  behavior  of  the  vector  potential A inside  a  toroid  did  not  arise  with modern toroidal devices in plasma physics and corresponding Hamiltonian treatment of charged particle dynamics in plasma, but can be traced to Niven's 2 nd  edition of Maxwell's treatise [131 ]; this first derivation of A was, however, shown to be erroneous by Schenkel et al. [100,101].

For a quick visualization concerning the overall configuration of the magnetic vector potential A for a current distribution j, Maxwell's equations of magnetostatics provide a nice observation, which to our knowledge was first given by Schlomka et al. [ 143 ], but has since been independently reported by other authors [133] as well, because of its great usefulness, e.g. in toroidal configurations. The vector potential A1 resulting from a current density distribution j 1 is in direction and magnitude equal to the magnetic field B2 generated by a current density distribution j 2 , when j 2 equals up to a factor c/4  (in Gauß units) the magnetic field B1 originating from j 1. In other words, the vector potential A may be constructed form the field B in the same way as B is calculated

from j. For the direction of A inside the toroid with a uniform and strictly meridional current distribution they showed that A has approximately the same direction as the generating current.

This  has  been  exploited  to  visualize  magnetic  vector  potentials  for  various  current density distributions deviating from circular shapes [100, 101, 133].

To our knowledge Schenkel et al. [100,101]  were then the first to calculate the vector potential A inside  and  outside  a  toroid  using  toroidal  coordinates,  implementing  a series expansion of 1/r in toroidal functions given explicitly already in the 19 th  century by Hübschmann [102], for various current density distributions on a toroidal surface. The toroidal (or ring-) functions used in the expansion are the associated Legendre functions with odd-half-integer degree and integer order. With increasing interest in magnetic fields generated by currents in toroidal geometry in plasma physics, eg. of toroidal stellarators, a large number of investigations, besides Schenkel et al. [100],e.g. by Carron et al. [ 133] Mirin [ 141], Haas [ 99a] and Cohn et al. [ 120-122]  began to focus  on  the  vector  potential  inside  a  toroidal  configuration,  apparently  without knowledge of previous treatises.

Methods for quasi-analytic solutions in terms of a vector potential A were investigated by  [141,  99a,135].  Extensive  applications  of  Hamiltonian  theory  and  equations  of motion assuring a much broader applicability beyond non-Hamiltonian guiding -center theory has been introduced by [ 95,96,99,106-115,118, 123].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01019v1-figure-8413699c-fdc4-4978-aa2c-9850346faf8d.jpg)

Fig. 29 Cut through a toroidal coil, e.g. coil 7 in fig 11, with orientation of the magnetic vector potential A inside the toroidal coil configuration with a strictly meridional current distribution  on  the  surface  of  the  coils,  calculated  using  expansions  for A given by Schenkel [100, 101]. The coil center lies on the toroidal arc as in fig 11; at the focus the toroidal coordinate  goes to infinity (see also fig 14a). In the proximity of the inner coil  surface  the  quality  of  the  convergence  of  the  series  expansion  is  decreasing. Lepton  trajectories  originating  at  the  origin  below  the  paper  plane  fly  towards  the observer in the detector plane (see fig. 11); in this figure all detectors are to be imagined to the right of the coil center on the invariant arcs.

Useful for application in the Hamilton equations, one derives in toroidal coordinates (   , , 𝜓) from the observation of the pure azimuthal B-field from 𝐵 ⃗ = ∇   × 𝐴 ⃗ ⃗ that the magnetic vector potential has only a  -component 𝐴 = (𝐴𝜂 = 0, 𝐴 [𝜂, 𝜃], 𝐴  𝜓 = 0) .

Fig  29  gives  for  the  cross  section  of  a  toroid  (in  a  geometry  as  considered  in  our trajectory  calculations with  OPERA3D) the vector potential A as  following  from  the series expansions of Schenkel et al. [100,101] and Mirin et.[141]; we notice a small disagreement with respect to expectations from current density on the inside of the toroidal coil in the proximity to the coil, probably due to slower series convergence.

We have embarked to find suitably  parametrized  approximations for the  numerical vector potential A (as displayed in fig 30) to be used in the toroidal Hamiltonian 𝐻 =

$$\frac { \left ( \cosh \eta - \cos \theta \right ) ^ { 2 } } { 2 m a ^ { 2 } } \left [ p _ { \eta } ^ { 2 } + p _ { \theta } ^ { 2 } + \frac { \left ( p _ { \psi } - A _ { \psi } \right ) ^ { 2 } } { s i n h ^ { 2 } \eta } \right ] \text{ whi}$$

$$+ \frac { ( p _ { \psi } - A _ { \psi } ) } { \circledast i m h ^ { 2 } m } \Big | \text{ which may be amenable to describe trajectories via}$$

the  Hamilton equations  related  to  the  observed  invariants  in  the  lepton  trajectories which we found using the numerical techniques in OPERA-3D. Details will be  subject of a forthcoming paper.

Besides analysis of particle transport in thermonuclear plasma in tokamaks or other fusion devices also experiments related to the Aharonov-Bohm effect have  stimulated a search for expressions for the magnetic vector potential A in toroidal configurations [133].

On the mathematical side it is important to note a few errors in the formulae given in the  literature  in  the  context  of  toroidal  coordinates  which  unfortunately  even  have propagated into Mathematica.

The  defining  equations  for  toroidal  functions  given  in  Wolfram  World[125]  are  the erroneous functions from the 9th printing of Abramowitz-Stegun 1 st  ed [126] which  in turn was only later corrected in their 2 nd  edition; it, was unfortunately propagated into

Mathematica  10  and  MathWorld.  The  equations  given  in  the  corrected  version  of Abramowitz-  Stegun  [127]  are  in  complete  agreement  with  the  equations  given  by Hübschmann[ 102] already in 1890; Mirin [141] has used a slightly different formulation for the toroidal functions.

We  also  point  out,  that  in  Fig,  4.04  in  Moon/Spencer  [98]  illustrating  toroidal coordinates the arrow on their  -coordinate is pointing outward  instead of inward. The direction of increasing coordinate theta is displayed correctly in fig. 2.09, which treats the corresponding bipolar coordinates.

## Appendix 2: transformation of cross sections

For the problem of relating measured  energy-angle and/or momentum/angle distributions in the laboratory frame to the corresponding distributions in the emitter frame it is necessary to consider the Jacobians of the Lorentz transformations of the distributions. As the individual terms of the matrix are significantly more elaborate than the final determinant, we give here the explicit expressions for both transformations, for the cross sections differential in the energy and also for those differential in the lepton momentum, as they were first given by Panofsky [132] and e.g. Hagedorn[143].

For evaluating the determinant of the Jacobian

$$\begin{smallmatrix} & & \bullet & & \bullet \\ & & \bullet & & \bullet \\ & & \overline { \partial E } \Big | _ { \cos \theta } & & \frac { \partial \cos \theta } { \partial E } \Big | _ { \cos \theta } \\ & & \frac { \partial E ^ { ^ { \prime } } } { \partial \cos \theta } \Big | _ { E } & & \frac { \partial \cos \theta ^ { ^ { \prime } } } { \partial \cos \theta } \Big | _ { E } \\ & & & & & & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bul. & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bulle & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bul! & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bul!. & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullets & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bulots & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bul; & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bul.: & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bul: & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bulclots & &
 \end{smallmatrix} & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bullet & & \bul. & & \bul. & &
 \end{smallmatrix} \end{smallmatrix}$$

it is straightforward to derive the four matrix elements

$$\text{it is straightforward to derive the four matrix elements} \\ \\ J _ { 1 1 } = \frac { \partial E ^ { \prime } } { \partial E } \Big | _ { \cos \theta } = \gamma _ { 0 } ( 1 - \beta _ { 0 } \cos \theta \ \frac { E } { p } ) \\ \\ J _ { 1 2 } = \frac { \partial \cos \theta ^ { \prime } } { \partial E } \Big | _ { \cos \theta } \\ \\ = \frac { \gamma _ { 0 } } { p ^ { \prime } c } ( \frac { E \cos \theta } { p c } - \beta _ { 0 } - \frac { E ^ { \prime } } { p ^ { \prime } c } \cos \theta ^ { \prime } + \frac { E ^ { \prime } E } { p ^ { \prime } p } \beta _ { 0 } \cos \theta \cos \theta ^ { \prime } ) \\ \\ J _ { 2 1 } = \frac { \partial E ^ { \prime } } { \partial \cos \theta } \Big | _ { E } = - \gamma _ { 0 } \beta _ { 0 } p c$$

$$J _ { 2 2 } = \frac { \partial \cos \theta ^ { \prime } } { \partial \cos \theta } \Big | _ { E } = \gamma _ { 0 } \, \frac { p } { p ^ { \prime } } \, ( 1 + \frac { E ^ { \prime } \beta _ { 0 } } { p ^ { \prime } c } \cos \theta ^ { \prime } )$$

One finds 𝐽 11 22 𝐽 -𝐽12𝐽21 = 𝑝 𝑝 ′ . From this directly follows eq.12:

$$\frac { d ^ { 2 } \sigma ^ { \prime } } { d E ^ { \prime } d \Omega ^ { \prime } } = \frac { p } { p ^ { \prime } } \frac { d ^ { 2 } \sigma } { d E d \Omega }.$$

For transformation of the momentum differential cross sections 𝑑 2 𝜎 𝑑𝑝𝑑Ω

one has to evaluate elements for the determinant of the corresponding Jacobian

$$\begin{vmatrix} \frac { \partial p ^ { ^ { \prime } } } { \partial p } \Big | _ { c o s \theta } & \frac { \partial c o s \theta ^ { ^ { \prime } } } { \partial p } \Big | _ { c o s \theta } \\ \frac { \partial p ^ { ^ { \prime } } } { \partial c o s \theta } \Big | _ { p } & \frac { \partial c o s \theta ^ { ^ { \prime } } } { \partial c o s \theta } \Big | _ { p } \end{vmatrix}, \\. \quad \. \quad. \quad. \quad. \quad.$$

here one derives the following  matrix elements

$$J _ { 1 1 } ^ { p } = \frac { \partial p \prime } { \partial p } = \frac { E ^ { \prime } } { p ^ { \prime } c } \frac { p c } { E } \gamma _ { 0 } ( 1 - \beta _ { 0 } \cos \theta \frac { E } { p c } )$$

$$J _ { 1 2 } ^ { p } = \frac { \partial \cos \theta ^ { \prime } } { \partial p } =$$

$$\frac { p } { p ^ { \prime } } \frac { \gamma _ { 0 } } { E } \left ( \frac { E } { p c } \cos \theta - \beta _ { 0 } - \frac { E ^ { \prime } } { p ^ { \prime } c } \cos \theta ^ { \prime } + \frac { \beta _ { 0 } E ^ { \prime } E } { p ^ { \prime } c p c } \cos \theta ^ { \prime } \cos \theta \right )$$

$$J _ { 2 1 } ^ { p } = \frac { \partial p ^ { \prime } } { \partial \cos \theta } = - \gamma _ { 0 } \beta _ { 0 } ( \frac { p E ^ { \prime } } { p ^ { \prime } } )$$

$$J _ { 2 2 } ^ { p } = \frac { \partial \cos \theta ^ { \prime } } { \partial \cos \theta } = \gamma _ { 0 } \, \frac { p } { p ^ { \prime } } \, ( 1 + \frac { E ^ { \prime } \beta _ { 0 } } { p ^ { \prime } c } \cos \theta ^ { \prime } )$$

$$\text{Now one finds} \quad J _ { 1 1 } ^ { \nu }$$

$$\ n d s \quad J _ { 1 1 } ^ { p } J _ { 2 2 } ^ { p } - J _ { 1 2 } ^ { p } J _ { 2 1 } ^ { p } = \frac { \left ( p ^ { \prime } \right ) ^ { 2 } } { p ^ { 2 } } \frac { E } { E ^ { \prime } } = \frac { \left ( p ^ { \prime } \right ) ^ { 2 } } { p ^ { 2 } } \frac { \gamma } { \gamma ^ { \prime } }.$$

From this now expression eq. 13 for obtaining the DDCS in the emitter frame from the DDCS in the laboratory frame

$$\frac { d ^ { 2 } \sigma ^ { \prime } ( p ^ { \prime }, \Omega ^ { \prime } ) } { d p ^ { \prime } d \Omega ^ { \prime } } = \frac { \gamma } { \gamma ^ { \prime } } \frac { ( p ^ { \prime } ) ^ { 2 } } { p ^ { 2 } } \frac { d ^ { 2 } \sigma } { d p d \Omega }$$

is directly derived.

For convenience we add the relation for the double differential cross sections in energy and momentum space:

$$\frac { d ^ { 2 } \sigma } { d E d \Omega } = \frac { d ^ { 2 } \sigma } { d p d \Omega } \frac { E } { p c }$$

## Acknowledgments

Within this project Yu. A. Litvinov has received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (grant agreement No 682841 'ASTRUm'')

P.M. Hillenbrand acknowledges support by DFG project HI 2009/1-1.

## 8. References

- 1. R. Ruffini et al. Physics Report 487 1-140 (2010)
- 2. J. Eichler, W. E. Meyerhof, Relativist. Atomic Coll., Academic Press (1995)
- 3. A. Voitkiv, J. Ullrich, Relativistic Collisions, Springer Series AOPP, Vol. 49 (2008)
- 4. H. Bethe et al., Proc Roy Soc A146 83-112 (1934)
- 5. H. A. Bethe et al., Phys. Rev. 93 , 768 (1954)
- 6. F. Sauter, Ann. Physik 412, 404 (1934)
- 7. A. Sommerfeld, Ann Physik 403 , 257 (1931)
- 8. J. Oppenheimer et al., Phys. Rev. 44, 53 (1933)
- 9. W. Heitler, Z. f. Phys. 84 , 145 (1933)
- 10. L. D. Landau, E. M. Lifshitz, Phys. Z. Sowjet Union 6 , 244 (1934)
- 11. G. Racah, Il Nuovo Cim 11 477-481 (1934) pair prod
- 12. G. Racah, Nuovo Cimento 14 (1937) 93
- 13. G. Baur et al. Phys. Rep. 453 1-27 (2007)
- 14. G. Baur et al. Phys. Rep. 364 359 (2002)
- 15. G. Baur et al. J. Phys G24 , 1657 (1998) pair
- 16.   K. Hencken et al. Phys. Rev. A59 841 (1999)
- 17.    B. Segev et al. Phys. Rev. C59 2753 (1999)
- 18. D. Ionescu, Habilitationsschrift, Freie Universität Berlin  (1997)
- 19. C. A. Bertulani et al. Phys. Rep. 163 299 (1988)
- 20. C. A. Bertulani et al., Ann. Rev. Nuc. Part. Sci. 55 , 271 (2005) UPC
- 21. A. Baltz et al. Phys. Rep. 458 1 (2008) UPC
- 22. A. Baltz et al.  Phys. Rev. A47 3444 (1993), Phys. Rev. A50 4842 (1994)

- 23. M. Strayer et al. Phys. Rev. A41 1399 (1990),
- 24. H. J. Bhabha, Proc. Roy Soc A152, 559 (1935)
- 25. H. J. Bhabha, Proc. Cambridge Phil. Soc. 31 , 394 (1935)
- 26. U. Becker et al. J. Phys. B19 , 1347 (1986)   pair
- 27. U. Becker et al. J Phys. B20 , 2075 and 6563 (1987)   bound free pair
- 28.   R.W. Tenzer, Dissertation, Univ. Giessen (1999)
- 29.  R. W. Tenzer et al. Eur. Phys. J. D 11 347-353 (2000 )
- 30.  M. Gail et al. 2003 J. Phys. B: At. Mol. Opt. Phys. 36 1397
- 31.  M. Gail et al. Phys. Scr. T92 317-318 (2001)   non-pert-th
- 32.    C. Agger   Dissertation Univ Aarhus (1996)   ph-induced pair
- 33.   A. Artemyev et al. Phys. Rev. A79 032713 (2009) neg cont
- 34.  A. Belkacem et al. Phys. Rev. A56 2806 (1997)
- 35.  D. Ionescu et al. Phys. Rev. A54 4960-4967 (1996) bound pair,
- 36.   D. Ionescu et al. Phys. Rev. A48 1176 (1993) free-free
- 37.   K. Rumrich et al. Phys. Rev. Lett. 66 2613 (1991) Phys. Rev. A47 215 (1993)
- 38.   J. Thiel et al. Phys. Rev. A46 2607 (1992)
- 39.   C. Vane et al Phys. Rev. A56 3682-3686 (1997)
- 40.   A. Belkacem et al. Phys. Rev. A58 1253-1255 (1998)
- 41.   H. Meier et al. Phys. Rev. A63 032713 (2001)   bound-free basel
- 42.   R. Lee et al   et al. Phys. Rev. A63 062712 (2001) bound-pair  theory
- 42a.  D. Condren 2007 J. Phys.: Conf. Ser. 58 315  pair 2 center
- 43. O.Busic et al .  Phys. Rev. A70 062707 (2004) bound pair th Scheid
- 44. B. Najjari et al. Phys. Rev. A80 012701 (2009)  neg cont capt
- 45. A. Baltz Phys. Rev. C74 054903 (2006) free pair th
- 46. A. Baltz et al. Phys. Rep. 458 1-171 (2008)   UPc-cern-pair
- 47. H. Olsson et al.   Phys. Rev. 114 887-904 (1959)  enBS-pair
- 48. H. Olsson et al.   Phys. Rev. 106 27 (1957)  enBS-pair
- 49. R. Pratt  Phys. Rev. 117 1017 (1960)  enBS-pair
- 49a. R. Pratt  Phys. Rev. 119 1619 (1961)  enBS-pair
- 50. R.  Jabbur et al. Rev. 129 1619 (1961)  enBS-pair
- 51. U. Fano Phys. Rev. 116 1156-1158 (1963)  enbs-pair ang distr

- 52. H. Meier et al Eu. Phys. J. C2 741-746 (1998)             pair enBS
- 53. R. Lee et al. Phys. Rev. A65 022102  (2002)      UPC
- 54. N. Claytor et al.  Phys. Rev. A55 R842 (1997) e-capture/loss
- 55. R. Bruce et al., Phys. Rev. ST Accel Beams 11 , 071002(2009)
- 56. R. Anholt et al. Adv. At. Mol. Phys. 22, 315 (1986)  pair-lumin-accel
- 57. T. Stoehlker et al. SPARC@HESR , A feasibility study , Darmstadt (2012)

## https://www.gsi.de/fileadmin/SPARC/documents/SPARC@HESR\_FS\_V26.pdf

- 58. W. Henning et al. (ed) FAIR, CDR , http://www.gsi.de/GSI-future/cdr 2001
- 59. R. Maier et al. HESR TDR V3.1.2 Darmstadt (2008)
- 60. Z. Anjelkovic et al. Technical Design Report SPARC@HESR: Instrumentation https://www.gsi.de/work/forschung/appamml/atomphysik/ap\_und\_fair/sparc/dokumen te.htm
- 61. V. Popov JETP 65 , 35 (1973)
- 62. R. Gerstein et al., Sov. J. Exp Theor. Phys. 30 , 358 (1969)
- 63. R. Gerstein et al., Lettere al  Nuovo Cimen. 1 ,835 (1969)
- 64. J. Reinhardt et al. Rep. Prog. Phys. 40 ,219 (1977)
- 65. W. Greiner et al. QED of Strong Fields Oxford Univ. Press, Oxford (1995)
- 66. O. Graf et al. Phys. Rev. Lett 61 2831-2834 (1988) pair ang corr
- 67. I. Tupitsyn et al., Phys. Rev. A 82 , 042701 (2010)
- 68. A. Artemyev et al. Phys. Rev. A 79 , 032713 (2009);
- A. Artemyev et al. Phys. Rev. A 67 , 052711 (2003);  Neg cont DR
- 69. B. Najjari et al. Phys. Rev. A 80 , 012701 (2009)
- 70. Y Yerokhin et al. Phys. Rev. A 92 , 042708 (2015)neg cont DR
- 71. R. Hagedorn Relativistic Kinematics, W. A. Benjamin, NY, (1963)
- 72. J. Blaton Mat.-Fys. Medd 24 1-31 (1950)  kinematik
- 73. V. Goldanski et al. Kinematic methods in Hgh Energy Physics (1989)
- 74. W. Trower Kinematics of particles UCRL-2426 U. of California (1965)
- 75. M.A. Herath Banda et al. , Phys. Rev. A 33 , 861 (1986)
- 76. P. Kienle et al. in Proc. Conf. NATO Adv. Study Inst. on QED of Strong Fields,
- Lahnstein, NATO Adv. Study Inst. Series B80 ( 1983)
- 77. H. Backe et al. Phys. Rev. Lett 40 1443 (1978)
- 78.  R. Krieg et al.  Nuc. Inst. Meth Phys B9 762-766  (1985)   spectrom

- 79.   E.Kankeleit et al. Nuc. Inst. Meth Phys A234 81-90  (1985)
- 80. G. Ter-Akopian et al. Int Journal Mod Phys. E24 1550016 (2015) positron spec
- 81.   OPERA-3D User Guide, Vector Fields Limited, 24 Bankside,
- Kidington, Oxford OX5 1JE, UK
- 82.   X. Wang et al. Meas. Sci. Technol. 18 161-166   (2007)
- 83. H. Wollnik, Optics of charged particles, (1987) Academic Press, N.Y.
- 84.  A. Kawachi et al.     Nuc. Inst. Meth Phys A416 253-262  (1998)  tori
- 85.  S. Wüstling et al.   Nuc. Inst. Meth Phys A568 382-390  (2006)   64pinSili
- 86.  S. Wüstling et al.    Nuc. Inst. Meth Phys A624 295-302  (2010)  64pin
- 87.  A. Efremov et al.  Nuc. Inst. Meth Phys A585 182-200  (2008)    panda-solen
- 88. J. Ullrich, V. Shevelko, Many-Particle Quantum Dynamics in Atomic and Molecular Fragmentation, Springer Series AOPP, Vol. 35 (2003)
- 89. G. Deco et al. 1989 J. Phys. B: At. Mol. Opt. Phys. 22 1043
- 91.  J. Picone  Nuc. Inst. Meth Phys 105 377-379  (1972) beta ICspec

;UCLEAR

- 92.  S. Dörscher, Dissertation    Univ. Hamburg (2013)    Yb-Mot
- 93. D. Robin et al Nuc. Inst. Meth Phys A659 484-493  (2011) toroid inclined coils
- 94. N. Petridis et al. Technical Report: Design of the TARGET@HESR (2013)
- 95. S. Abdullaev and K. Finken, Phys Plasma 9 , 4193 (2002)
- 96. S. Abdullaev, Lecture notes in physics 691, Berlin Springer (2006)
- 97. S. Abdullaev, Springer Series on Atomic, Optical and Plasma Physics 78,1 Berlin, Springer (2014)
- 98. P. Moon, D. Spencer, Field theory handbook, 1961 Springer Heidelberg

Figure 4.04 in Moon-Spencer depicting the coordinate surfaces contains an error in the coordinate-direction of  : the radius of the spheroid decreases with increasing  .

- 99. G. Carter et al., Quart. J. of Mechanics and Appl. Math 18, 87-106(1965)
- 99a. H. Haas, Archiv. f. Elektrotechnik 58 (1976) 197ff   Springer Berlin
- 99b. A.H. Boozer, Physics of Fluids 25 , 520 (1982) toroidal coord
- 100. G. Schenkel et al., Annalen d. Physik 426 541 (1939) Toroid
- 101. G. Schenkel, 'Über das Vektorpotential und elektrische Feld stromumflossener Toroide', Dissertation, Techn. Hochschule Hannover   1938
- 102. H. Hübschmann, 'Die Ringfunktion und ihre Anwendung auf die elektrostatischen Probleme des Ringes' ,Dissertation Univ. Leipzig (1890)
- 103. https://en.wikipedia.org/wiki/Circles\_of\_Apollonius
- 104. T. Jahnke et al. J. Electron Spectroscopy, Related Phenomena 141, 229(2004)
- 105. H. AlfvØn, Ark. Mat. Astron. Fys. 27A-22, 1 (1940)
- 106. H. E. Mynick, Physics of Fluids 23 , 1888 (1980) guiding Hamiltonian

- 107. H. D. Littlejohn, The Physics of Fluids 24 , 1730 (1981) guide Hamiltonian
- 108. H. Weitzner, Physics of Fluids 24 , 2280 (1981), axisymm field
- 109. R. B. White et al., Physics of Fluids 25 , 575 (1982) drift Hamiltonian
- 110. A. Boozer, Physics of Fluids 27 , 2441 (1984); drift Hamiltonian
- 111. R.B. White et. al. Physics of Fluids 27 , 2455 (1984) drift Hamilton
- 112. J. D. Meiss et al., Physics of Fluids B 2 , 2563 (1990); guiding center
- 113. Q. Yao et al., Physics of Fluids B 5 , 3852 (1993); guiding center-toroid
- 114. S. Zenitani et al., Physics of Plasmas 23 , 102102 (2016); guiding  in plasma
- 115. Q. Yao et al., Physics of Fluids B 4, 2429 (1992); guiding center-toroid
- 116. J. M. Greene et al., The Physics of Fluids 4 , 875 (1961) cylind-toroid-perturb
- 117. A. Salat, Z. f. Naturforschung 40a, 959 (1985)
- 118. J. Cary et al. Rev. Mod. Phys. 81 , 693 (2009)
- 119.  E. Thomas et al. Comp. Phys. Comm. 100, 31 (1997)
- 120. H. Cohl et al. Phys. Rev. A64 , 052509 (2001)
- 121. H. Cohl et al. Astrophys. J. 527 , 86 (1999)
- 122. H. Cohl et al. Astron. Nachr. 321 363 (2000)
- 123. A. Boozer, Princeton Plasma Phys. Rep. PPPL-1989 and PPPL-2094 (1984)
- 124. P. M. Hillenbrand et al. Phys. Rev. A90 ,  042713 (2014), A90 ,  022707 (2014), A91 , 022705 (2015)
- 125. http://mathworld.wolfram.com/ToroidalFunction.html
- 126. M. Abramowitz, and I. Stegun, Handbook of Math. Function., 9 th  printing, NY, Dover, p336, 1972
- 127. NIST Digital Library of Mathematical Functions, http://dlmf.nist.gov/ Sect 14.19 Toroidal and Ringfunctions
- 128. http://people.math.sfu.ca/~cbm/aands/ corrected version of Abramowitz-Stegun, Handbook of Math Functions; Sect 8.11 Toroidal Functions
- 129. G. Bonheure et al. Czech. J. Phys. 56 B61 (2006)
- 130. V. Mukhovatov et al. Nucl. Fusion 11 , 605 (1971)
- 131. J. Cl. Maxwell, A treatise on electricity and magnetism II, 2 nd  ed (1881)
- 132. K. Dedrick, Rev.  Mod.  Phys. 34 429  (1962)  and  K.  Dedrick, Stanford University,Report. M-226(1960)
- 133. N. J. Carron, Am J. Phys 63, 717(1965) toroid
- 134. H. Essen et al. Progress In Electromagnetics Research B. 46 357 (2013)
- 135. D. Bhadra, Review of Scientific Instruments 39 1536 1968
- 136. L. Brouwer et al. Proc. PAC2013 WEPBA12 p 907(2013) scal mag. Potential
- 137. H. Wollnik, Optics of Charged Particles, Academic Press, Orlando (1987)
- 138. B. van Milligen et al. Computer Physics Communications 81, 74 (1994)

- 139. E. Heine, Handbuch der Kugelfunktionen, 2 nd  ed. Berlin 1881
- 140. D. Stern, GSFC-NASA Report X-641-71-373, (1971) Greenbelt, Md, USA
- 141. A. A. Mirin et al.  UCRL-52069 Report, Berkeley (1976).
- 142. R. Hagedorn, Relativistic Kinematics, Lect. Notes in Physics, J.D. Jackson. ed., W. Benjamin Inc, New York (1963)
- 143. T: Schlomka , Ztschft f. techn. Physik 16 480 (1935)
Siegbert Hagmann Pierre-Michel Hillenbrand Yurij Litvinov Uwe Spillmann Thomas Stöhlker
The Magnetic Toroidal Sector as A broad-band Electron-Positron Pair Spectrometer I. Lepton Trajectories
We report an analysis of electron-optical properties of a toroidal magnetic
sector spectrometer and examine parameters for its implementation in a
relativistic heavy-ion storage ring like HESR. For studies of free-free pair
production in heavy-ion atom collisions this spectrometer exhibits very high
efficiencies for coincident $e^+e^-$ pair spectroscopy over a wide range of
momenta of emitted lepton pairs. The high coincidence efficiency of the
spectrometer is the key for stringent tests of theoretical predictions for the
coincident positron- and electron- emission characteristics and for the phase
space correlation of lepton vector momenta in free-free pair production.