Introduction
With the EM polarization extension it is possible to track polarized particles (leptons and photons). Special emphasis will be put in the proper treatment of polarized matter and its interaction with longitudinal polarized electrons/positrons or circularly polarized photons, which is for instance essential for the simulation of positron polarimetry. The implementation is base on Stokes vectors [McM61]. Further details can be found in [Lai08].
In its current state, the following polarization dependent processes are considered:
Bhabha/Møller scattering,
Positron Annihilation,
Compton scattering,
Pair creation,
Bremsstrahlung.
Several simulation packages for the realistic description of the development of electromagnetic showers in matter have been developed. A prominent example of such codes is EGS (Electron Gamma Shower) [NHR85]. For this simulation framework extensions with the treatment of polarized particles exist [Flottmann93a, LKNS, NBH93]; the most complete has been developed by K. Flöttmann [Flottmann93a]. It is based on the matrix formalism [McM61], which enables a very general treatment of polarization. However, the Flöttmann extension concentrates on evaluation of polarization transfer, i.e. the effects of polarization induced asymmetries are neglected, and interactions with polarized media are not considered.
Another important simulation tool for detector studies is Geant3 [eal85]. Here also some effort has been made to include polarization [aal, Hoo97], but these extensions are not publicly available.
In general the implementation of polarization in this EM polarization library follows very closely the approach by K. Flöttmann [Flottmann93a]. The basic principle is to associate a Stokes vector to each particle and track the mean polarization from one interaction to another. The basics for this approach is the matrix formalism as introduced in [McM61].
Stokes vector
The Stokes vector [McM61, Sto52] is a rather simple object (in comparison to e.g. the spin density matrix), three real numbers are sufficient for the characterization of the polarization state of any single electron, positron or photon. Using Stokes vectors all possible polarization states can be described, i.e. circular and linear polarized photons can be handled with the same formalism as longitudinal and transverse polarized electron/positrons.
The Stokes vector can be used also for beams, in the sense that it defines a mean polarization.
In the EM polarization library the Stokes vector is defined as follows:
Photons |
Electrons |
|
---|---|---|
\(\xi_1\) |
linear polarization |
polarization in x direction |
\(\xi_2\) |
linear polarization but \(\pi/4\) to right |
polarization in y direction |
\(\xi_3\) |
circular polarization |
polarization in z direction |
This definition is assumed in the particle reference frame, i.e. with the momentum of the particle pointing to the z direction, cf. also next section about coordinate transformations. Correspondingly a 100% longitudinally polarized electron or positron is characterized by
where \(\pm1\) corresponds to spin parallel (anti parallel) to particle’s momentum. Note that this definition is similar, but not identical to the definition used in McMaster [McM61].
Many scattering cross sections of polarized processes using Stokes vectors for the characterization of initial and final states are available in [McM61]. In general a differential cross section has the form
i.e. it is a function of the polarization states of the initial particles \({{\boldsymbol{\zeta}}}^{(1)}\) and \({{\boldsymbol{\zeta}}}^{(2)}\), as well as of the polarization states of the final state particles \({{\boldsymbol{\xi}}}^{(1)}\) and \({{\boldsymbol{\xi}}}^{(2)}\) (in addition to the kinematic variables \(E\), \(\theta\), and \(\phi\)).
Consequently, in a simulation we have to account for
Asymmetries:
Polarization of beam (\({{\boldsymbol{\zeta}}}^{(1)}\)) and target (\({{\boldsymbol{\zeta}}}^{(2)}\)) can induce azimuthal and polar asymmetries, and may also influence on the total cross section (Geant4: GetMeanFreePath()).
Polarization transfer / depolarization effects
The dependence on the final state polarizations defines a possible transfer from initial polarization to final state particles.
Transfer matrix
Using the formalism of McMaster, differential cross section and polarization transfer from the initial state (\({{\boldsymbol{\zeta}}}^{(1)}\)) to one final state particle (\({{\boldsymbol{\xi}}}^{(1)}\)) are combined in an interaction matrix \(T\):
where \(I\) and \(O\) are the incoming and outgoing currents, respectively. In general the \(4\times4\) matrix \(T\) depends on the target polarization \({{\boldsymbol{\zeta}}}^{(2)}\) (and of course on the kinematic variables \(E\), \(\theta\), \(\phi\)). Similarly one can define a matrix defining the polarization transfer to second final state particle like
In this framework the transfer matrix \(T\) is of the form
The matrix elements \(T_{ij}\) can be identified as (unpolarized) differential cross section (\(S\)), polarized differential cross section (\(A_j\)), polarization transfer (\(M_{ij}\)), and (de)polarization (\(P_i\)). In the Flöttmann extension the elements \(A_j\) and \(P_i\) have been neglected, thus concentrating on polarization transfer only. Using the full matrix takes now all polarization effects into account.
The transformation matrix, i.e. the dependence of the mean polarization of final state particles, can be derived from the asymmetry of the differential cross section w.r.t. this particular polarization. Where the asymmetry is defined as usual by
The mean final state polarizations can be determined coefficient by coefficient. In general, the differential cross section is a linear function of the polarizations, i.e.
In this form, the mean polarization of the final state can be read off easily, and one obtains
Note that the mean polarization states do not depend on the correlation matrix \(M_{(\zeta^{(1)},\zeta^{(2)})}\). In order to account for correlation one has to generate single particle Stokes vector explicitly, i.e. on an event by event basis. However, this implementation generates mean polarization states, and neglects correlation effects.
Coordinate transformations
Three different coordinate systems are used in the evaluation of polarization states:
World frame The geometry of the target, and the momenta of all particles in Geant4 are noted in the world frame \(X\), \(Y\), \(Z\) (the global reference frame, GRF). It is the basis of the calculation of any other coordinate system.
Particle frame Each particle is carrying its own coordinate system. In this system the direction of motion coincides with the \(z\)-direction. Geant4 provides a transformation from any particle frame to the World frame by the method
G4ThreeMomemtum::rotateUz()
. Thus, the \(y\)-axis of the particle reference frame (PRF) lies in the \(X\)-\(Y\)-plane of the world frame.The Stokes vector of any moving particle is defined w.r.t. the corresponding particle frame. Particles at rest (e.g. electrons of a media) use the world frame as particle frame.
Interaction frame For the evaluation of the polarization transfer another coordinate system is used, defined by the scattering plane, cf. Fig. 34. There the \(z\)-axis is defined by the direction of motion of the incoming particle. The scattering plane is spanned by the \(z\)-axis and the \(x\)-axis, in a way, that the direction of particle 1 has a positive \(x\) component. The definition of particle 1 depends on the process, for instance in Compton scattering, the outgoing photon is referred as particle 11.
All frames are right handed.
Polarized beam and material
Polarization of beam particles is well established. It can be used for simulating low-energy Compton scattering of linear polarized photons. The interpretation as Stokes vector allows now the usage in a more general framework. The polarization state of a (initial) beam particle can be fixed using the standard ParticleGunMessenger class. For example, the class G4ParticleGun provides the method SetParticlePolarization(), which is usually accessible via:
/gun/polarization <Sx> <Sy> <Sz>
in a macro file.
In addition for the simulation of polarized media, a possibility to assign Stokes vectors to physical volumes is provided by a new class, the so-called G4PolarizationManager. The procedure to assign a polarization vector to a media, is done during the detector construction. There the logical volumes with certain polarization are made known to polarization manager. One example DetectorConstruction might look like follows:
G4double Targetthickness = .010*mm;
G4double Targetradius = 2.5*mm;
G4Tubs* solidTarget =
new G4Tubs("solidTarget",
0.0,
Targetradius,
Targetthickness/2,
0.0*deg,
360.0*deg );
G4LogicalVolume * logicalTarget =
new G4LogicalVolume(solidTarget,
iron,
"logicalTarget",
0,0,0);
G4VPhysicalVolume * physicalTarget =
new G4PVPlacement(0,G4ThreeVector(0.*mm, 0.*mm, 0.*mm),
logicalTarget,
"physicalTarget",
worldLogical,
false,
0);
G4PolarizationManager * polMgr = G4PolarizationManager::GetInstance();
polMgr->SetVolumePolarization(logicalTarget,G4ThreeVector(0.,0.,0.08));
Once a logical volume is known to the G4PolarizationManager, its polarization vector can be accessed from a macro file by its name, e.g. the polarization of the logical volume called “logicalTarget” can be changed via:
/polarization/volume/set logicalTarget 0. 0. -0.08
Note, the polarization of a material is stated in the world frame.
Ionisation
Method
The class G4ePolarizedIonization provides continuous and discrete energy losses of polarized electrons and positrons in a material. It evaluates polarization transfer and – if the material is polarized – asymmetries in the explicit delta rays production. The implementation baseline follows the approach derived for the class G4eIonization described in Mean Energy Loss and Ionisation. For continuous energy losses the effects of a polarized beam or target are negligible provided the separation cut \(T_{\rm cut}\) is small, and are therefore not considered separately. On the other hand, in the explicit production of delta rays by Møller or Bhabha scattering, the effects of polarization on total cross section and mean free path, on distribution of final state particles and the average polarization of final state particles are taken into account.
Total cross section and mean free path
Kinematics of Bhabha and Møller scattering is fixed by initial energy
and variable
which is the part of kinetic energy of initial particle carried out by scatter. Lower kinematic limit for \(\epsilon\) is \(0\), but in order to avoid divergences in both total and differential cross sections one sets
where \(T_{min}\) has meaning of minimal kinetic energy of secondary electron. And, \(\epsilon_{\rm max}=1(1/2)\) for Bhabha(Møller) scatterings.
Total Møller cross section
The total cross section of the polarized Møller scattering can be expressed as follows
where the \(r_e\) is classical electron radius, and
Total Bhabha cross section
The total cross section of the polarized Bhabha scattering can be expressed as follows
where
Mean free path
With the help of the total polarized Møller cross section one can define a longitudinal asymmetry \(A^M_L\) and the transverse asymmetry \(A^M_T\), by
and
Similarly, using the polarized Bhabha cross section one can introduce a longitudinal asymmetry \(A^B_L\) and the transverse asymmetry \(A^B_T\) via
and
These asymmetries are depicted in Fig. 35, Fig. 36 for Møller and Fig. 37, Fig. 38 for Bhabha.
If both beam and target are polarized the mean free path as defined in Ionisation has to be modified. In the class G4ePolarizedIonization the polarized mean free path \(\lambda^{\rm pol}\) is derived from the unpolarized mean free path \(\lambda^{\rm unpol}\) via
Sampling the final state
Differential cross section
The polarized differential cross section is rather complicated. The full result can be found in [eal, FM57, Ste58]. In G4PolarizedMollerCrossSection the complete result is available taking all mass effects into account, with only binding effects neglected. Here we state only the ultra-relativistic approximation (URA), to show the general dependencies.
The corresponding cross section for Bhabha cross section is implemented in G4PolarizedBhabhaCrossSection. In the ultra-relativistic approximation it reads
where
\(r_e\) |
classical electron radius |
\(\gamma\) |
\(E_{k_1}/m_e c^2\) |
\(\epsilon\) |
(\(E_{p_1}-m_e c^2)/(E_{k_1}-m_e c^2)\) |
\(E_{k_1}\) |
energy of the incident electron/positron |
\(E_{p_1}\) |
energy of the scattered electron/positron |
\(m_e c^2\) |
electron mass |
\({{\boldsymbol{\zeta}}}^{(1)}\) |
Stokes vector of the incoming electron/positron |
\({{\boldsymbol{\zeta}}}^{(2)}\) |
Stokes vector of the target electron |
\({{\boldsymbol{\xi}}}^{(1)}\) |
Stokes vector of the outgoing electron/positron |
\({{\boldsymbol{\xi}}}^{(2)}\) |
Stokes vector of the outgoing (2nd) electron . |
Sampling
The delta ray is sampled according to methods discussed in Monte Carlo Methods. After exploitation of the symmetry in the Møller cross section under exchanging \(\epsilon\) versus \((1-\epsilon)\), the differential cross section can be approximated by a simple function \(f^M(\epsilon)\):
with the kinematic limits given by
A similar function \(f^B(\epsilon)\) can be found for Bhabha scattering:
with the kinematic limits given by
The kinematic of the delta ray production is constructed by the following steps:
\(\epsilon\) is sampled from \(f(\epsilon)\)
calculate the differential cross section, depending on the initial polarizations \({{\boldsymbol{\zeta}}}^{(1)}\) and \({{\boldsymbol{\zeta}}}^{(2)}\).
\(\epsilon\) is accepted with the probability defined by ratio of the differential cross section over the approximation function.
The \(\varphi\) is diced uniformly.
\(\varphi\) is determined from the differential cross section, depending on the initial polarizations \({{\boldsymbol{\zeta}}}^{(1)}\) and \({{\boldsymbol{\zeta}}}^{(2)}\)
Note, for initial states without transverse polarization components, the \(\varphi\) distribution is always uniform. In Fig. 39 and Fig. 40 the asymmetries indicate the influence of polarization. In general the effect is largest around \(\epsilon=1/2\).
After both \(\phi\) and \(\epsilon\) are known, the kinematic can be constructed fully. Using momentum conservation the momenta of the scattered incident particle and the ejected electron are constructed in global coordinate system.
Polarization transfer
After the kinematics is fixed the polarization properties of the outgoing particles are determined. Using the dependence of the differential cross section on the final state polarization a mean polarization is calculated according to method described in Introduction.
The resulting polarization transfer functions \(\xi^{(1,2)}_3(\epsilon)\) are depicted in Fig. 41, Fig. 42, and Fig. 43, Fig. 44.
Positron - Electron Annihilation
Method
The class G4eplusPolarizedAnnihilation simulates annihilation of polarized positrons with electrons in a material. The implementation baseline follows the approach derived for the class G4eplusAnnihilation described in Positron - Electron Annihilation. It evaluates polarization transfer and – if the material is polarized – asymmetries in the produced photons. Thus, it takes the effects of polarization on total cross section and mean free path, on distribution of final state photons into account. And calculates the average polarization of these generated photons. The material electrons are assumed to be free and at rest.
Total cross section and mean free path
Kinematics of annihilation process is fixed by initial energy
and variable
which is the part of total energy available in initial state carried out by first photon. This variable has the following kinematical limits
Total Cross Section
The total cross section of the annihilation of a polarized \(e^+e^-\) pair into two photons could be expressed as follows
where
Mean free path
With the help of the total polarized annihilation cross section one can define a longitudinal asymmetry \(A^A_L\) and the transverse asymmetry \(A^A_T\), by
and
These asymmetries are depicted in Fig. 45, Fig. 46.
If both incident positron and target electron are polarized the mean free path as defined in section Positron - Electron Annihilation has to be modified. The polarized mean free path \(\lambda^{\rm pol}\) is derived from the unpolarized mean free path \(\lambda^{\rm unpol}\) via
Sampling the final state
Differential Cross Section
The fully polarized differential cross section is implemented in the class G4PolarizedAnnihilationCrossSection, which takes all mass effects into account, but binding effects are neglected [eal, Pag57]. In the ultra-relativistic approximation (URA) and concentrating on longitudinal polarization states only the cross section is rather simple:
where
\(r_e\) |
classical electron radius |
\(\gamma\) |
\(E_{k_1}/m_e c^2\) |
\(E_{k_1}\) |
energy of the incident positron |
\(m_e c^2\) |
electron mass |
\({{\boldsymbol{\zeta}}}^{(1)}\) |
Stokes vector of the incoming positron |
\({{\boldsymbol{\zeta}}}^{(2)}\) |
Stokes vector of the target electron |
\({{\boldsymbol{\xi}}}^{(1)}\) |
Stokes vector of the 1st photon |
\({{\boldsymbol{\xi}}}^{(2)}\) |
Stokes vector of the 2nd photon . |
Sampling
The photon energy is sampled according to methods discussed in Monte Carlo Methods. After exploitation of the symmetry in the Annihilation cross section under exchanging \(\epsilon\) versus \((1-\epsilon)\), the differential cross section can be approximated by a simple function \(f(\epsilon)\):
with the kinematic limits given by
The kinematic of the two photon final state is constructed by the following steps:
\(\epsilon\) is sampled from \(f(\epsilon)\)
calculate the differential cross section, depending on the initial polarizations \({{\boldsymbol{\zeta}}}^{(1)}\) and \({{\boldsymbol{\zeta}}}^{(2)}\).
\(\epsilon\) is accepted with the probability defined by the ratio of the differential cross section over the approximation function \(f(\epsilon)\).
The \(\varphi\) is diced uniformly.
\(\varphi\) is determined from the differential cross section, depending on the initial polarizations \({{\boldsymbol{\zeta}}}^{(1)}\) and \({{\boldsymbol{\zeta}}}^{(2)}\).
A short overview over the sampling method is given in Monte Carlo Methods. In Fig. 48 the asymmetries indicate the influence of polarization for an 10MeV incoming positron. The actual behavior is very sensitive to the energy of the incoming positron.
Polarization transfer
After the kinematics is fixed the polarization of the outgoing photon is determined. Using the dependence of the differential cross section on the final state polarizations a mean polarization is calculated for each photon according to method described in section Introduction.
The resulting polarization transfer functions \(\xi^{(1,2)}(\epsilon)\) are depicted in Fig. 49, Fig. 50.
Annihilation at Rest
The method AtRestDoIt
treats the special case where a positron comes
to rest before annihilating. It generates two photons, each with energy
\(E_{p_{1/2}}=m c^2\) and an isotropic angular distribution.
Starting with the differential cross section for annihilation with
positron and electron spins opposed and parallel,
respectively, [Pag57]
In the limit \(\beta\to0\) the cross section \(d\sigma_1\) becomes one, and the cross section \(d\sigma_2\) vanishes. For the opposed spin state, the total angular momentum is zero and we have a uniform photon distribution. For the parallel case the total angular momentum is 1. Here the two photon final state is forbidden by angular momentum conservation, and it can be assumed that higher order processes (e.g. three photon final state) play a dominant role. However, in reality 100% polarized electron targets do not exist, consequently there are always electrons with opposite spin, where the positron can annihilate with. Final state polarization does not play a role for the decay products of a spin zero state, and can be safely neglected (is set to zero).
Polarized Compton scattering
Method
The class G4PolarizedCompton simulates Compton scattering of polarized photons with (possibly polarized) electrons in a material. The implementation follows the approach described for the class G4ComptonScattering introduced in Compton scattering. Here the explicit production of a Compton scattered photon and the ejected electron is considered taking the effects of polarization on total cross section and mean free path as well as on the distribution of final state particles into account. Further the average polarizations of the scattered photon and electron are calculated. The material electrons are assumed to be free and at rest.
Total cross section and mean free path
Kinematics of the Compton process is fixed by the initial energy
and the variable
which is the part of total energy available in initial state carried out by scattered photon, and the scattering angle
The variable \(\epsilon\) has the following limits:
Total Cross Section
The total cross section of Compton scattering reads
where
and
Mean free path
When simulating the Compton scattering of a photon with an atomic electron, an empirical cross section formula is used, which reproduces the cross section data down to 10 keV (see Compton scattering). If both beam and target are polarized this mean free path has to be corrected.
In the class G4ComptonScattering the polarized mean free path \(\lambda^{\rm pol}\) is defined on the basis of the the unpolarized mean free path \(\lambda^{\rm unpol}\) via
where
is the expected asymmetry from the the total polarized Compton cross section given above. This asymmetry is depicted in Fig. 51, Fig. 52.
Sampling the final state
Differential Compton Cross Section
In the ultra-relativistic approximation the dependence of the differential cross section on the longitudinal/circular degree of polarization is very simple. It reads
where
\(r_e\) |
classical electron radius |
\(X\) |
\(E_{k_1}/m_e c^2\) |
\(E_{k_1}\) |
energy of the incident photon |
\(m_e c^2\) |
electron mass |
The fully polarized differential cross section is available in the class G4PolarizedComptonCrossSection. It takes all mass effects into account, but binding effects are neglected [eal, LT54a, LT54b]. The cross section dependence on \(\epsilon\) for right handed circularly polarized photons and longitudinally polarized electrons is plotted in Fig. 53, Fig. 54.
Sampling
The photon energy is sampled according to methods discussed in Monte Carlo Methods. The differential cross section can be approximated by a simple function \(\Phi(\epsilon)\):
with the kinematic limits given by
The kinematic of the scattered photon is constructed by the following steps:
\(\epsilon\) is sampled from \(\Phi(\epsilon)\)
calculate the differential cross section, depending on the initial polarizations \({{\boldsymbol{\zeta}}}^{(1)}\) and \({{\boldsymbol{\zeta}}}^{(2)}\), which the correct normalization.
\(\epsilon\) is accepted with the probability defined by ratio of the differential cross section over the approximation function.
The \(\varphi\) is diced uniformly.
\(\varphi\) is determined from the differential cross section, depending on the initial polarizations \({{\boldsymbol{\zeta}}}^{(1)}\) and \({{\boldsymbol{\zeta}}}^{(2)}\).
In Fig. 53, Fig. 54 the asymmetries indicate the influence of polarization for an 10 MeV incoming positron. The actual behavior is very sensitive to energy of the incoming positron.
Polarization transfer
After the kinematics is fixed the polarization of the outgoing photon is determined. Using the dependence of the differential cross section on the final state polarizations a mean polarization is calculated for each photon according to the method described in section Introduction.
The resulting polarization transfer functions \(\xi^{(1,2)}(\epsilon)\) are depicted in Fig. 55, Fig. 56.
Polarized Bremsstrahlung for electron and positron
Method
The polarized version of Bremsstrahlung is based on the unpolarized cross section. Energy loss, mean free path, and distribution of explicitly generated final state particles are treated by the unpolarized version G4eBremsstrahlung. For details consult Bremsstrahlung.
The remaining task is to attribute polarization vectors to the generated final state particles, which is discussed in the following.
Polarization in gamma conversion and bremsstrahlung
Gamma conversion and bremsstrahlung are cross-symmetric processes (i.e. the Feynman diagram for electron bremsstrahlung can be obtained from the gamma conversion diagram by flipping the incoming photon and outgoing positron lines) and their cross sections closely related. For both processes, the interaction occurs in the field of the nucleus and the total and differential cross section are polarization independent. Therefore, only the polarization transfer from the polarized incoming particle to the outgoing particles is taken into account.
For both processes, the scattering can be formulated by:
Where \(\mathcal{N}_{1}(k_{\mathcal {N}_{1}}, {{\boldsymbol{\zeta}}}^{(\mathcal {N}_{1})})\) and \(\mathcal{N}_{2}(p_{\mathcal{N}_{2}}, {{\boldsymbol{\xi}}}^{(\mathcal{N}_{2})})\) are the initial and final state of the field of the nucleus respectively assumed to be unchanged, at rest and unpolarized. This leads to \(k_{\mathcal {N}_{1}} = k_{\mathcal {N}_{2}} = 0\) and \({{\boldsymbol{\zeta}}}^{(\mathcal {N}_{1})} = {{\boldsymbol{\xi}}}^{(\mathcal{N}_{2})} = 0\).
In the case of gamma conversion process: \(\mathcal{K}_{1}(k_{1},{{\boldsymbol{\zeta}}}^{(1)})\) is the incoming photon initial state with momentum \(k_{1}\) and polarization state \({{\boldsymbol{\zeta}}}^{(1)}\). \(\mathcal{P}_{1}(p_{1},{{\boldsymbol{\xi}}}^{(1)})\) and \(\mathcal{P}_{2}(p_{2},{{\boldsymbol{\xi}}}^{(2)})\) are the two photons final states with momenta \(p_{1}\) and \(p_{2}\) and polarization states \({{\boldsymbol{\xi}}}^{(1)}\) and \({{\boldsymbol{\xi}}}^{(2)}\).
In the case of bremsstrahlung process: \(\mathcal{K}_{1}(k_{1},{{\boldsymbol{\zeta}}}^{(1)})\) is the incoming lepton \(e^{-}(e^{+})\) initial state with momentum \(k_{1}\) and polarization state \({{\boldsymbol{\zeta}}}^{(1)}\). \(\mathcal{P}_{1}(p_{1},{{\boldsymbol{\xi}}}^{(1)})\) is the lepton \(e^{-}(e^{+})\) final state with momentum \(p_{1}\) and polarization state \({{\boldsymbol{\xi}}}^{(1)}\). \(\mathcal{P}_{2}(p_{2},{{\boldsymbol{\xi}}}^{(2)})\) is the bremsstrahlung photon in final state with momentum \(p_{2}\) and polarization state \({{\boldsymbol{\xi}}}^{(2)}\).
Polarization transfer from the lepton e+e- to a photon
The polarization transfer from an electron (positron) to a photon in a bremsstrahlung process was first calculated by Olsen and Maximon [OM59] taking into account both Coulomb and screening effects. In the Stokes vector formalism, the \(e^{-}(e^{+})\) polarization state can be transformed to a photon polarization finale state by means of interaction matrix \(T_{\gamma}^{b}\). It defined via
and
where
and
\(\epsilon_{1}\) |
Total energy of the incoming lepton \(e^{+}(e^{-})\) in units \(mc^{2}\) |
\(\epsilon_{2}\) |
Total energy of the outgoing lepton \(e^{+}(e^{-})\) in units \(mc^{2}\) |
\(k\) |
\(=(\epsilon_{1}-\epsilon_{2})\), the energy of the bremsstrahlung photon in units of \(mc^{2}\) |
\({{\boldsymbol{p}}}\) |
Electron (positron) initial momentum in units \(mc\) |
\({{\boldsymbol{k}}}\) |
Bremsstrahlung photon momentum in units \(mc\) |
\({{\boldsymbol{u}}}\) |
Component of \({{\boldsymbol{p}}}\) perpendicular to \({{\boldsymbol{k}}}\) in units \(mc\) and \(u=\vert {{\boldsymbol{u}}} \vert\) |
\(\hat\xi\) |
\(= 1/(1+u^{2})\) |
Coulomb and screening effects are contained in \(\Gamma\), defined as follows
with
\(f(Z)\) is the Coulomb correction term derived by Davies, Bethe and Maximon [eal54]. \(\mathcal{F}({\hat\xi}/{\delta})\) contains the screening effects and is zero for \(\Delta \le 0.5\) (No screening effects). For \(0.5 \le \Delta \le 120\) (intermediate screening) it is a slowly decreasing function. The \(\mathcal{F}({\hat\xi}/{\delta})\) values versus \(\Delta\) are given in Table 28 [KM59] and used with a linear interpolation in between.
The polarization vector of the incoming \(e^{-}(e^{+})\) must be rotated into the frame defined by the scattering plane (x-z-plane) and the direction of the outgoing photon (z-axis). The resulting polarization vector of the bremsstrahlung photon is also given in this frame.
\(\Delta\) |
\(-\mathcal{F}\left({\hat\xi}/{\delta}\right)\) | \(\Delta\) |
\(-\mathcal{F}\left({\hat\xi}/{\delta}\right)\) |
|
---|---|---|---|
0.5 |
0.0145 |
40.0 |
2.00 |
1.0 |
0.0490 |
45.0 |
2.114 |
2.0 |
0.1400 |
50.0 |
2.216 |
4.0 |
0.3312 |
60.0 |
2.393 |
8.0 |
0.6758 |
70.0 |
2.545 |
15.0 |
1.126 |
80.0 |
2.676 |
20.0 |
1.367 |
90.0 |
2.793 |
25.0 |
1.564 |
100.0 |
2.897 |
30.0 |
1.731 |
120.0 |
3.078 |
35.0 |
1.875 |
Using Eq.(193) and the transfer matrix given by Eq.(194) the bremsstrahlung photon polarization state in the Stokes formalism [McM61, McM54] is given by
Remaining polarization of the lepton after emitting a bremsstrahlung photon
The \(e^{-}(e^{+})\) polarization final state after emitting a bremsstrahlung photon can be calculated using the interaction matrix \(T_{l}^{b}\) which describes the lepton depolarization. The polarization vector for the outgoing \(e^{-}(e^{+})\) is not given by Olsen and Maximon. However, their results can be used to calculate the following transfer matrix [Flottmann93b, Hoo97].
where
and
\(\epsilon_{1}\) |
Total energy of the incoming \(e^{+}/e^{-}\) in units \(mc^{2}\) |
\(\epsilon_{2}\) |
Total energy of the outgoing \(e^{+}/e^{-}\) in units \(mc^{2}\) |
\(k\) |
\(=(\epsilon_{1}-\epsilon_{2})\), energy of the photon in units of \(mc^{2}\) |
\({{\boldsymbol{p}}}\) |
Electron (positron) initial momentum in units \(mc\) |
\({{\boldsymbol{k}}}\) |
Photon momentum in units \(mc\) |
\({{\boldsymbol{u}}}\) |
Component of \({{\boldsymbol{p}}}\) perpendicular to \({{\boldsymbol{k}}}\) in units \(mc\) and \(u=\vert {{\boldsymbol{u}}} \vert\) |
Using Eq.(196) and the transfer matrix given by Eq.(197) the \(e^{-}(e^{+})\) polarization state after emitting a bremsstrahlung photon is given in the Stokes formalism by
Polarized Gamma conversion into an electron–positron pair
Method
The polarized version of gamma conversion is based on the EM standard process G4GammaConversion. Mean free path and the distribution of explicitly generated final state particles are treated by this version. For details consult Gamma Conversion into e+e- Pair.
The remaining task is to attribute polarization vectors to the generated final state leptons, which is discussed in the following.
Polarization transfer from the photon to the two leptons
Gamma conversion process is essentially the inverse process of bremsstrahlung and the interaction matrix is obtained by inverting the rows and columns of the bremsstrahlung matrix and changing the sign of \(\epsilon_{2}\), cf. Polarized Bremsstrahlung for electron and positron. It follows from the work by Olsen and Maximon [OM59] that the polarization state \(\xi^{(1)}\) of an electron or positron after pair production is obtained by
and
where
and
\(\epsilon_{1}\) |
total energy of the first lepton \(e^{+}(e^{-})\) in units \(mc^{2}\) |
\(\epsilon_{2}\) |
total energy of the second lepton \(e^{-}(e^{+})\) in units \(mc^{2}\) |
\(k=(\epsilon_{1}+\epsilon_{2})\) |
energy of the incoming photon in units of \(mc^{2}\) |
\({{\boldsymbol{p}}}\) |
electron=positron initial momentum in units \(mc\) |
\({{\boldsymbol{k}}}\) |
photon momentum in units \(mc\) |
\({{\boldsymbol{u}}}\) |
electron/positron initial momentum in units \(mc\) |
\(u\) |
\(=\vert {{\boldsymbol{u}}} \vert\) |
Coulomb and screening effects are contained in \(\Gamma\), defined in section Polarized Bremsstrahlung for electron and positron.
Using Eq.(198) and the transfer matrix given by Eq.(199) the polarization state of the produced \(e^{-}(e^{+})\) is given in the Stokes formalism by:
Polarized Photoelectric Effect
Method
This section describes the basic formulas of polarization transfer in the photoelectric effect class (G4PolarizedPhotoElectricEffect). The photoelectric effect is the emission of electrons from matter upon the absorption of electromagnetic radiation, such as ultraviolet radiation or x-rays. The energy of the photon is completely absorbed by the electron and, if sufficient, the electron can escape from the material with a finite kinetic energy. A single photon can only eject a single electron, as the energy of one photon is only absorbed by one electron. The electrons that are emitted are often called photoelectrons. If the photon energy is higher than the binding energy the remaining energy is transferred to the electron as a kinetic energy
In Geant4 the photoelectric effect process is taken into account if:
Where \(k\) is the incoming photon energy and \(B_{shell}\) the electron binding energy provided by the class G4AtomicShells.
The polarized version of the photoelectric effect is based on the EM standard process G4PhotoElectricEffect. Mean free path and the distribution of explicitly generated final state particles are treated by this version. For details consult section PhotoElectric Effect.
The remaining task is to attribute polarization vectors to the generated final state electron, which is discussed in the following.
Polarization transfer
The polarization state of an incoming polarized photon is described by the Stokes vector \(\vec{\zeta}^{(1)}\). The polarization transfer to the photoelectron can be described in the Stokes formalism using the same approach as for the bremsstrahlung and gamma conversion processes, cf. Polarized Bremsstrahlung for electron and positron and Polarized Gamma conversion into an electron–positron pair. The relation between the photoelectron’s Stokes parameters and the incoming photon’s Stokes parameters is described by the interaction matrix \(T_{l}^{p}\) derived from H. Olsen [OV58] and reviewed by H.W McMaster [McM61]:
In general, for the photoelectric effect as a two-body scattering, the cross section should be correlated with the spin states of the incoming photon and the target electron. In our implementation the target electron is not polarized and only the polarization transfer from the photon to the photoelectron is taken into account. In this case the cross section of the process remains polarization independent. To compute the matrix elements we take advantage of the available kinematic variables provided by the generic G4PhotoelectricEffect class. To compute the photoelectron spin state (Stokes parameters), four main parameters are needed:
The incoming photon Stokes vector \(\vec{\zeta}^{(1)}\)
The incoming photon’s energy \(k\).
the photoelectron’s kinetic energy \(E_{kin}^{e^-}\) or the Lorentz factors \(\beta\) and \(\gamma\).
The photoelectron’s polar angle \(\theta\) or \(\cos\theta\).
The interaction matrix derived by H. Olsen [OV58] is given by:
where
Using Eq.(201) and the transfer matrix given by Eq.(202) the polarization state of the produced \(e^{-}\) is given in the Stokes formalism by:
From equation (203) one can see that a longitudinally (transversally) polarized photoelectron can only be produced if the incoming photon is circularly polarized.
Bibliography
- aal
G. Alexander at al. Undulator-based production of polarized positrons: a proposal for the 50-gev beam in the fftb. Technical Report SLAC-TN-04-018, SLAC-PROPOSAL-E-166, SLAC.
- eal54(1,2)
H. Davies et al. Theory of bremsstrahlung and pair production. ii. integral cross section for pair production. Phys. Rev., 93(4):788–795, 1954.
- eal(1,2,3)
P. Starovoitov et al. in preparation.
- eal85
R. Brun et al. The geant3 electromagnetic shower program and a comparison with the egs3 code. Technical Report CERN-DD/85/1, CERN, 1985.
- Flottmann93a(1,2,3)
K. Flöttmann. DESY-93-161, DESY Hamburg, phd thesis edition, 1993.
- Flottmann93b
K. Flöttmann. Investigations toward the development of polarized and unpolarized high intensity positron sources for linear colliders. Universitat Hamburg, phd thesis edition, 1993.
- FM57
G. W. Ford and C. J. Mullin. Phys. Rev., 108:477, 1957.
- Hoo97(1,2)
Johannes Marinus Hoogduin. Electron, positron and photon polarimetry. Rijksuniversiteit Groningen, phd thesis edition, 1997.
- KM59(1,2)
H. W. Koch and J. W. Motz. Bremsstrahlung cross-section formulas and related data. Rev. Mod. Phys., 31(4):920–955, 1959.
- Lai08
K. Laihem. Measurement of the positron polarization at an helical undulator based positron source for the International Linear Collider ILC. Humboldt University Berlin, Germany, phd thesis edition, 2008.
- LT54a
F.W. Lipps and H.A. Tolhoek. Polarization phenomena of electrons and photons i. Physica, 20:85, 1954.
- LT54b
F.W. Lipps and H.A. Tolhoek. Polarization phenomena of electrons and photons ii. Physica, 20:395, 1954.
- LKNS
J. C. Liu, T. Kotseroglou, W. R. Nelson, and D. C. Schultz. Polarization study for nlc positron source using egs4. Technical Report SLAC-PUB-8477, SLAC.
- McM61(1,2,3,4,5,6,7,8)
W. H. McMaster. Rev. Mod. Phys., 33():8, 1961.
- McM54
W.H. McMaster. Polarization and the stokes parameters. American Journal of Physics, 22(6):351–362, 1954.
- NBH93
Y. Namito, S. Ban, and H. Hirayama. Implementation of linearly polarized photon scattering into the egs4 code. Nucl. Instrum. Meth. A, 332:277, 1993.
- NHR85
W.R. Nelson, H. Hirayama, and D.W.O. Rogers. EGS4 code system. SLAC, Dec 1985. SLAC-265, UC-32.
- OM59(1,2)
H. Olsen and L.C. Maximon. Photon and electron polarization in high- energy bremsstrahlung and pair production with screening. Physical Review, 114:887–904, 1959.
- OV58(1,2)
H. Olsen and Kgl. N. Videnskab. Selskabs Forh., 1958.
- Pag57(1,2)
L. A. Page. Polarization effects in the two-quantum annihilation of positrons. Phys. Rev., 106:394–398, 1957.
- Ste58
P. Stehle. Phys. Rev., 110:1458, 1958.
- Sto52
G. Stokes. Trans. Cambridge Phil. Soc., 9 edition, 1852.
- 1
Note, for an incoming particle travelling on the \(Z\)-axis (of GRF), the \(y\)-axis of the PRF of both outgoing particles is parallel to the \(y\)-axis of the interaction frame.