Transition radiation

The Relationship of Transition Radiation to X-ray Čerenkov Radiation

X-ray transition radiation (XTR) occurs when a relativistic charged particle passes from one medium to another of a different dielectric permittivity. In order to describe this process it is useful to begin with an explanation of X-ray Čerenkov radiation, which is closely related.

The mean number of X-ray Čerenkov radiation (XCR) photons of frequency ω emitted into an angle θ per unit distance along a particle trajectory is [Gri02b]:

(202)d3ˉNxcrdωdxdθ2=απcωcθ2Im{Z}.

Here the quantity Z is introduced as the complex formation zone of XCR in the medium:

(203)Z=L1iLl,L=cω[γ2+ω2pω2+θ2]1,γ2=1β2.

with l and ωp the photon absorption length and the plasma frequency, respectively, in the medium. For the case of a transparent medium, l and the complex formation zone reduces to the coherence length L of XCR. The coherence length roughly corresponds to that part of the trajectory in which an XCR photon can be created.

Introducing a complex quantity Z with its imaginary part proportional to the absorption cross-section (l1) is required in order to account for absorption in the medium. Usually, ω2p/ω2c/ωl. Then it can be seen from Eqs. (202) and (203) that the number of emitted XCR photons is considerably suppressed and disappears in the limit of a transparent medium. This is caused by the destructive interference between the photons emitted from different parts of the particle trajectory.

The destructive interference of X-ray Čerenkov radiation is removed if the particle crosses a boundary between two media with different dielectric permittivities, ϵ, where

ϵ=1ω2pω2+icωl.

Here the standard high-frequency approximation for the dielectric permittivity has been used. This is valid for energy transfers larger than the K-shell excitation potential.

If layers of media are alternated with spacings of order L, the X-ray radiation yield from a trajectory of unit length can be increased by roughly l/L times. The radiation produced in this case is called X-ray transition radiation (XTR).

Calculating the X-ray Transition Radiation Yield

Using the methods developed in Ref.[Gri02b] one can derive the relation describing the mean number of XTR photons generated per unit photon frequency and θ2 inside the radiator for a general XTR radiator consisting of n different absorbing media with fluctuating thicknesses:

d2ˉNindωdθ2=απc2ωθ2Re{n1i=1(ZiZi+1)2+2n1k=1k1i=1(ZiZi+1)[kj=i+1Fj](ZkZk+1)},
Fj=exp[tj2Zj].

In the case of gamma distributed gap thicknesses (foam or fiber radiators) the values Fj, (j=1,2) can be estimated as:

Fj=0dtj(νjˉtj)νjtνj1jΓ(νj)exp[νjtjˉtjitj2Zj]=[1+iˉtj2Zjνj]νj,

where Zj is the complex formation zone of XTR (similar to relation (203) for XCR) in the j-th medium [eal00, Gri02a]. Γ is the Euler gamma function, ˉtj is the mean thickness of the j-th medium in the radiator and νj>0 is the parameter roughly describing the relative fluctuations of tj. In fact, the relative fluctuation is δtj/ˉtj1/νj.

In the particular case of n foils of the first medium (Z1,F1) interspersed with gas gaps of the second medium (Z2,F2), one obtains:

(204)d2ˉNindωdθ2=2απc2ωθ2Re{R(n)},F=F1F2,
(205)R(n)=(Z1Z2)2{n(1F1)(1F2)1F+(1F1)2F2[1Fn](1F)2}.

Here R(n) is the stack factor reflecting the radiator geometry. The integration of ((204)) with respect to θ2 can be simplified for the case of a regular radiator (ν1,2), transparent in terms of XTR generation media, and n1 [Gar71]. The frequency spectrum of emitted XTR photons is given by:

(206)dˉNindω=10γ20dθ2d2ˉNindωdθ2=4αnπω(C1+C2)2kmaxk=kmin(kCmin)(kC1)2(k+C2)2sin2[πt1t1+t2(k+C2)],
C1,2=t1,2(ω21ω22)4πcω,Cmin=14πc[ω(t1+t2)γ2+t1ω21+t2ω22ω].

The sum in (206) is defined by terms with kkmin corresponding to the region of θ0. Therefore kmin should be the nearest to Cmin integer kminCmin. The value of kmax is defined by the maximum emission angle θ2max10γ2. It can be evaluated as the integer part of

Cmax=Cmin+ω(t1+t2)4πc10γ2,kmaxkmin1021031.

Numerically, however, only a few tens of terms contribute substantially to the sum, that is, one can choose kmaxkmin+20. Eq.(206) corresponds to the spectrum of the total number of photons emitted inside a regular transparent radiator. Therefore the mean interaction length, λXTR, of the XTR process in this kind of radiator can be introduced as:

λXTR=n(t1+t2)[ωmaxωmindωdˉNindω]1,

where ωmin1 keV, and ωmax100 keV for the majority of high energy physics experiments. Its value is constant along the particle trajectory in the approximation of a transparent regular radiator. The spectrum of the total number of XTR photons after regular transparent radiator is defined by (206) with:

nneff=n1k=0exp[k(σ1t1+σ2t2)]=1exp[n(σ1t1+σ2t2)]1exp[(σ1t1+σ2t2)],

where σ1 and σ2 are the photo-absorption cross-sections corresponding to the photon frequency ω in the first and the second medium, respectively. With this correction taken into account the XTR absorption in the radiator ((206)) corresponds to the results of [FS75]. In the more general case of the flux of XTR photons after a radiator, the XTR absorption can be taken into account with a calculation based on the stack factor derived in [GMGY75]:

(207)R(n)flux=(L1L2)2{1Qn1Q(1+Q1)(1+F)2F12Q1F22(1F)(1F1)(Q1F1)F2(QnFn)(1F)(QF)},
Q=Q1Q2,Qj=exp[tj/lj]=exp[σjtj],j=1,2.

Both XTR energy loss (205) and flux (207) models can be implemented as a discrete electromagnetic process (see below).

Simulating X-ray Transition Radiation Production

A typical XTR radiator consists of many (100) boundaries between different materials. To improve the tracking performance in such a volume one can introduce an artificial material [eal00], which is the geometrical mixture of foil and gas contents. Here is an example:

// In DetectorConstruction of an application
// Preparation of mixed radiator material
foilGasRatio  = fRadThickness/(fRadThickness+fGasGap);
foilDensity  = 1.39*g/cm3;     // Mylar
gasDensity   = 1.2928*mg/cm3 ; // Air
totDensity   = foilDensity*foilGasRatio +
               gasDensity*(1.0-foilGasRatio);
fractionFoil =  foilDensity*foilGasRatio/totDensity;
fractionGas  =  gasDensity*(1.0-foilGasRatio)/totDensity;
G4Material* radiatorMat = new G4Material("radiatorMat",
                                          totDensity,
                                          ncomponents = 2 );
radiatorMat->AddMaterial( Mylar, fractionFoil );
radiatorMat->AddMaterial( Air,   fractionGas  );
G4cout << *(G4Material::GetMaterialTable()) << G4endl;
// materials of the TR radiator
fRadiatorMat = radiatorMat;   // artificial for geometry
fFoilMat     = Mylar;
fGasMat      = Air;

This artificial material will be assigned to the logical volume in which XTR will be generated:

solidRadiator = new G4Box("Radiator",
                           1.1*AbsorberRadius ,
                           1.1*AbsorberRadius,
                           0.5*radThick        );
logicRadiator = new G4LogicalVolume( solidRadiator,
                                     fRadiatorMat,  // !!!
                                    "Radiator");
physiRadiator = new G4PVPlacement(0,
                                   G4ThreeVector(0,0,zRad),
                                   "Radiator", logicRadiator,
                                   physiWorld, false, 0       );

XTR photons generated by a relativistic charged particle intersecting a radiator with 2n interfaces between different media can be simulated by using the following algorithm. First the total number of XTR photons is estimated using a Poisson distribution about the mean number of photons given by the following expression:

ˉN(n)=ω2ω1dωθ2max0dθ2d2ˉN(n)dωdθ2=2απc2ω2ω1ωdωθ2max0θ2dθ2Re{R(n)}.

Here θ2max10γ2, ω11 keV, ω2100 keV, and R(n) correspond to the geometry of the experiment. For events in which the number of XTR photons is not equal to zero, the energy and angle of each XTR quantum is sampled from the integral distributions obtained by the numerical integration of expression (204). For example, the integral energy spectrum of emitted XTR photons, ˉN(n)>ω, is defined from the following integral distribution:

ˉN(n)>ω=2απc2ω2ωωdωθ2max0θ2dθ2Re{R(n)}.

In Geant4 XTR generation inside or after radiators is described as a discrete electromagnetic process. It is convenient for the description of tracks in magnetic fields and can be used for the cases when the radiating charge experiences a scattering inside the radiator. The base class G4VXTRenergyLoss is responsible for the creation of tables with integral energy and angular distributions of XTR photons. It also contains the PostDoIt function providing XTR photon generation and motion (if fExitFlux=true) through a XTR radiator to its boundary. Particular models like G4RegularXTRadiator implement the pure virtual function GetStackFactor, which calculates the response of the XTR radiator reflecting its geometry. Included below are some comments for the declaration of XTR in a user application.

In the physics list one should pass to the XTR process additional details of the XTR radiator involved:

// In PhysicsList of an application
else if (particleName == "e-")  // Construct processes for electron with XTR
{
   pmanager->AddProcess(new G4MultipleScattering, -1, 1,1 );
   pmanager->AddProcess(new G4eBremsstrahlung(),  -1,-1,1 );
   pmanager->AddProcess(new Em10StepCut(),        -1,-1,1 );
// in regular radiators:
   pmanager->AddDiscreteProcess(
   new G4RegularXTRadiator        // XTR dEdx in general regular radiator
// new G4XTRRegularRadModel        - XTR flux after general regular radiator
// new G4TransparentRegXTRadiator  - XTR dEdx in transparent
//                                   regular radiator
// new G4XTRTransparentRegRadModel - XTR flux after transparent
//                                   regular radiator
                         (pDet->GetLogicalRadiator(), // XTR radiator

                          pDet->GetFoilMaterial(), // real foil
                          pDet->GetGasMaterial(),  // real gas
                          pDet->GetFoilThick(),    // real geometry
                          pDet->GetGasThick(),
                          pDet->GetFoilNumber(),
                          "RegularXTRadiator"));
// or for foam/fiber radiators:
   pmanager->AddDiscreteProcess(
   new G4GammaXTRadiator    // - XTR dEdx in general foam/fiber radiator
// new G4XTRGammaRadModel   - XTR flux after general foam/fiber radiator
                          ( pDet->GetLogicalRadiator(),
                            1000.,
                            100.,
                            pDet->GetFoilMaterial(),
                            pDet->GetGasMaterial(),
                            pDet->GetFoilThick(),
                            pDet->GetGasThick(),
                            pDet->GetFoilNumber(),
                            "GammaXTRadiator"));
}

Here for the foam/fiber radiators the values 1000 and 100 are the ν parameters (which can be varied) of the Gamma distribution for the foil and gas gaps, respectively. Classes G4TransparentRegXTRadiator and G4XTRTransparentRegRadModel correspond (206) to n and neff, respectively.

Bibliography

eal00(1,2)

J. Apostolakis et al. Comput. Phys. Commun., 132():241, 2000.

FS75

C.W. Fabian and W. Struczinski. Physics Letters, B57():483, 1975.

GMGY75

G.M. Garibian, L.A. Gevorgian and C. Yang. Sov. Phys.- JETP, 39():265, 1975.

Gar71

G.M. Garibyan. Sov. Phys. JETP, 32():23, 1971.

Gri02a

V.M. Grichine. Generation of x-ray transition radiation inside complex radiators. Physics Letters B, 525(3-4):225–239, jan 2002. URL: https://doi.org/10.1016/S0370-2693(01)01443-5.

Gri02b(1,2)

V.M. Grichine. On the energy-angle distribution of cherenkov radiation in an absorbing medium. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 482(3):629–633, apr 2002. URL: https://doi.org/10.1016/S0168-9002(01)01927-1.