R.T. Jones
August 10, 1997
Coherent bremsstrahlung is common technique used in high energy physics to produce intense beams of highly polarized and monochromatic photons. In the literature, the term ``coherent bremsstrahlung'' is used to refer to two different effects. The first takes place when a highly relativistic electron is travelling inside a solid in a direction parallel to a crystal axis. In an effect known as channeling, the electron settles into a transverse well formed by chains of atomic nuclei, and enters the strong field regime in which higher-order QED processes dominate over Born terms in the production of bremsstrahlung. This appears experimentally as an enhancement over the normal bremsstrahlung spectrum at low photon energies, characteristic of multi-photon radiation. Theoretical models for this process employ a semi-classical picture for the radiation field, in which coherent states are used in the place of states with definite numbers of photons, hence the use of the term ``coherent bremsstrahlung''. An example of such a calculation can be found here. In the more recent literature this is referred to as channeling radiation.
The second meaning for the term ``coherent bremsstrahlung,'' which is the how it is used in this article, is the enhancement that occurs in normal bremsstrahlung inside a crystal when the momentum transfer from the electron to the atom matches a reciprocal lattice vector. This process, which is analogous to Bragg scattering, can be well described by first-order QED. The underlying cross section is the same as for Bethe-Heitler radiation, with the form factor augmented by additional factors for the crystal structure. The characteristics of the radiation pattern, however, are quite different from Bethe-Heitler (incoherent) bremsstrahlung, and these differences are what makes it interesting for our purposes. The reason is that for a given photon angle and energy, the incoherent spectrum is a sum over a continuum of momentum-transfer vectors, whereas in the coherent case the sum is discrete. This is best illustrated by a general description of the kinematics of bremsstrahlung.
We are interested in the case where the electron is ultra-relativistic. The momentum transfer q to the atom is small compared to the atomic mass, so q = (0, qT, ql) and the energy transfer is negligible. The bulk of the cross section comes from momentum transfer values close to the peak of the function |1-F(q)|2/q2 which corresponds to the atomic scale, around 10KeV. This allows one to keep only the leading terms in q/me. Momentum conservation at each vertex leads to the following relations.
The small-angle approximation is used because we are interested in
forward angles within a few m/E radians. Apart from a
small transverse kick qT (order 1ur for an 8GeV
electron), the intermediate momentum vector points along the direction
of p. By contrast, the scattering angle at the second vertex
scales as the electron mass (order 100ur for an 8GeV electron) so that
one can consider the intermediate momentum to be along the incident
beam axis and the transverse components of k and p'
must balance each other. The initial-state radiation diagram looks
somewhat different, but the kinematics are the same. These relations
follow.
Except for bremsstrahlung very near the end-point or at large angle, where the cross section is going to zero, the longitudinal component of q must be of order me2/E, which is 32eV at 8GeV electron beam energy. The transverse component is not constrained by this condition and takes on values on the typical atomic scale of several KeV, leading to the conclusion that a momentum transfer q very nearly transverse to the incident direction is characteristic of bremsstrahlung scattering. The sum over all values of q which gives the incoherent bremsstrahlung spectrum of Bethe-Heitler leads to a distribution in bremsstrahlung angle which is independent of energy and forward-peaked with a characteristic opening angle of m/E, henceforth called 1. The integral over q conceals the angle-energy dispersion that is present at a particular value of q, derived above from kinematics. This dispersion reappears in the case of coherent bremsstrahlung, where only a few discrete vectors contribute to the sum.
While it is not possible in general to orient the crystal so that only one lattice coherent peak appears in the photon spectrum, in a particular range of the photon spectrum there is generally one that is dominant. For individual values of q with xmax well below the end-point there is a high degree of polarization of the outgoing photon in the plane of q and the incident beam. By rotating the crystal with respect to the beam it is possible to orient in any transverse plane the polarization of the coherent component of the photon beam in the vicinity of an isolated coherent peak.
As discussed above, the hardest part of the coherent photon spectrum in any peak is kinematically confined to a cone much smaller than the characteristic bremsstrahlung angle m/E, making it possible to increase the relative intensity of the photon spectrum at the peak by collimation. Ideally such a technique would permit the generation of essentially monochromatic photon beams, with the low-energy bremsstrahlung background reduced to a negligible level. In this study I explore the potential of this technique for use as an alternative to tagged-bremsstrahlung or backscattered-Compton beams at CEBAF in the 8GeV+ era.
An example of an actual measurement made with such a beam at SLAC is described by Kaune et.al. Their paper also contains a compact description of the physics of coherent bremsstrahlung, and formed the starting point for this study.
The scattering angle is assumed to be an implicit function of gi and x as given by Eq. 5. Consider the maximum of a given coherent peak, which occurs at xmax and =0. With crystal momentum g nearly transverse to the beam, the dependence on gi enters through factors |1-F(g)|2/g2 and exp(-Agi2). The first is the virtual photon flux factor. The second is the probability that a given photon can be absorbed or emitted by the crystal without absorbing or emitting additional phonons. The appearance of phonons in the process would introduce additional momenta to the problem, changing the kinematics and contributing to radiation that does not interfere with the phonon-free process. The following plot of these factors for carbon (diamond) is taken from Bilokon et.al.
The upper [lower] curve excludes [includes] the exponential factor. The abscissa is shown in units of me. Markers are placed on the x axis to show the magnitudes of the lowest-order reciprocal lattice vectors. In a same paper a survey of a large number of single-element crystals was performed, and the value of the maximum in the above figure calculated for each. It appears as column 7 in the following table.
The quantity M listed in last column, which the authors call the merit factor, is simply column 7 multiplied by the number of lattice sites of the crystal per unit volume. The lesson drawn from this table is that only elements with the lowest Z are interesting as crystal radiators. Apart from a density factor t/AV, the merit factor collects together almost all of the crystal-dependence of the cross section. The one thing left over is the factor (S/n)2 which has a maximum of 1. For a great many lattice vectors, symmetry in the unit cell gives rise to a cancellation that causes S to vanish. Therefore the relevant value of M is not the maximum, but rather the value at the lowest-order lattice vector for which S does not vanish. This the authors call the performance parameter P , tabulated in table 2 for the most relevant crystals from table 1, as a ratio of the performance parameter for diamond.
Considering its good mechanical properties and its availability as a high-quality crystal, diamond emerges the clear winner among the cases studied. The compound crystal Be2C is somewhat superior in its performance, although little is known about how easy these crystals would be fabricate and to handle. The performance increment is probably not sufficiently great to stimulate much interest in their development.
The vertical axis is intensity, a dimensionless quantity I = x dN/dx in which the incoherent cross section is approximately flat. The solid dots are the measurements. The dashed curve is the calculated ideal spectrum for a thin crystal, and the solid curve is the calculation including multiple-scattering and beam divergence effects. The dot-dashed curve is the incoherent contribution alone. The agreement with the solid curve supports the claim that experimental spectra, including the defects present in a real beam and target, can be reliably calculated within the simple framework described in this report.
The crystal structure of diamond is shown in the figure at the right. Consider first an incident electron beam directed along a1. The only reciprocal lattice vectors with small projections onto a1 would be those lying in the plane (0,1,0), (0,0,1). Ideally these would have vanishing projections onto the beam axis, but finite emittance and imperfect crystal alignment cause a certain spread in ql about zero. Since low ql means low x, this leads to enhanced bremsstrahlung production at the low end of the photon spectrum, such as appears in the above figure. To raise the energy of the coherent peak, one can tilt the crystal by a small angle about one of the symmetry axes perpendicular to the beam; either (0,1,1) or (0,-1,1). After this rotation only a few rows of lattice vectors contribute significantly to the sum, instead of a whole plane, with each row showing up as a discrete peak (see the above figure). There is still the problem, however, that there is a single row passing through the origin in q-space and parallel to the rotation axis. The lattice vectors in this row all have a vanishing projection onto the beam axis, and through beam divergence etc. give rise to a very large coherent peak at zero photon energy with a tail extending under the other coherent peaks. Note that in the above spectrum the vertical axis is Intensity, not photon flux, so that the normal bremsstrahlung spectrum is approximately flat (see dot-dashed line). The user of a bremsstrahlung photon beam would generally be ill- prepared to accept such a large surplus of background photons over what is already present in a conventional bremsstrahlung beam. The authors of the paper from which the above figure was taken were able to reproduce this peak in their calculation, and identify its cause. They cut it off in the above figure, otherwise it would suppress the scale.
The production of a coherent low energy peak can be eliminated completely by tilting around both axes (0,1,1) and (0,-1,1) simultaneously. The first tilt is by a large enough angle that all of the peaks in the above figure are pushed up to x=1 (where they disappear, leaving only the peak at zero) and then the second tilt is tuned to move the zero-peak up to the desired energy. The zero-peak turns out to be a whole sequence of peaks, but the lowest-energy one is the strongest. It is not possible to select only one lattice vector, and so obtain only one peak in the spectrum, but the remaining peaks are nicely separated from one another in energy and so do not really cause problems to the user of the coherent photon beam. One side effect of this dual-rotation procedure is that, by narrowing to just one the number of lattice vectors contributing to a peak, the photons within a peak are very highly polarized in the plane containing the lattice vector and the beam. A drawback of the dual-rotation procedure is that the coherent enhancement is less when a single lattice vector contributes than when a whole row of vectors contribute, a factor of 3 for diamond on the (0,2,2) vector. I take the point of view that a clean beam is more important than raw photon flux per electron.
Using the above procedure to fix the crystal orientation, I show below the total bremsstrahlung spectrum for a 8GeV electron beam. The angles were adjusted for peak photon flux at half the end-point. By plotting the total bremsstrahlung over incoherent, the degree of coherent enhancement is made clear. The wiggles in the vicinity of the sharp edges are smoothing artifacts of the plotting package. Note that the spectrum drops to unity (no coherent component) at the infrared limit.
Actually the distinction between coherent and incoherent components is somewhat artificial in these calculations. The incoherent radiation is not the same as that from an amorphous radiator of the same thickness. When one replaces the integral over q with a discrete sum, the sum is only carried out for relatively small integers (h,k,l), beyond which the terms are squelched by the phonon decoherence factor exp(-Aq2). Beyond that the crystal structure is irrelevant and the atoms contribute individually. The large-q part of the sum in these calculations is replaced by an integral, as for Bethe-Heitler, but here the integral must exclude the piece covered by the discrete sum in order to prevent double-counting. This large-q part, which has the same spectrum as ordinary bremsstrahlung, is called the incoherent component in these calculations. For diamond, this incoherent part is 15% less than would be obtained with an amorphous carbon radiator of the same thickness in g/cm2.
As discussed in the section above on kinematics, there is a tight correlation within the coherent peaks between the bremsstrahlung energy and angle of emission with respect to the incident electron direction. By collimating the photon beam, one can eliminate photons beyond a certain angle, and narrow the width of the peaks in energy. In the next three plots is shown the effects of progressively collimating the photon beam to smaller and smaller angles.
Plotting the spectra in this way displays in a clear way how collimation narrows the peaks and increases the coherent fraction of the beam. The polarization of the coherent component in these peaks is greater than 80%. The blue-shaded region of Fig. 5 gives the relative number of polarized photons. Outside the coherent peaks the polarization is zero.
The probability density p(u) differential in du2
represents the dispersion of impacts on the collimator plane from
photons generated at fixed bremsstrahlung angle, as illustrated in
the figure at the right below.
This function is calculated by numerical integration. It is plotted in Fig. 6 for three different thicknesses of diamond radiator, using a collimator of nominal radius 0.5 m/E. Note that at 100um thickness, the acceptance never rises above 50%, even at zero degrees, due to multiple scattering.
The effects of multiple scattering become more pronounced the more narrow is the collimation. For the 0.25 m/E collimator the shape of the energy peak is dominated by multiple scattering, even for the very thin radiator thickness of 10um. This is shown in Fig. 7. The colour codes have the same meaning as in Fig. 6.
Taking into account the effects of beam emittance on the collimator-plane angular resolution amounts to evaluating the constant 0 from a given beam specification. Referring to the above picture of photon impact on the collimator plane, consider the following for a very thin radiator so that multiple scattering can be neglected. In the limit of negligible angular divergence of the beam, the distribution p(u) of hits on the collimator plane will be essentially the image of the beam spot. If the beam is focused to a negligible spot at the radiator position then this distribution will be the momentum-space image of the beam in the plane transverse to the nominal incident direction. A realistic situation falls somewhere in between, so that 0 depends in some way both on coordinate-space and momentum-space r.m.s. widths of the beam at the radiator. The picture is much simpler if one considers the virtual image of the beam at the position of the collimator. The image is virtual because the real electron beam is broadened by multiple scattering and bent into the dump, so it looks quite different at the collimator plane. The virtual image is formed by free-motion transport of the electron beam phase-space distribution, as it is at the entrance to the radiator, downstream to the entrance plane of the collimator. Considered there, p(u) is simply the coordinate-space image of the beam in the transverse plane, converted to radians by dividing by the radiator-collimator distance.
The product of coordinate-space and momentum-space r.m.s. widths of the beam, called emittance, is a fixed property of the accelerator. It is only modified by radiation or interactions of the beam with matter, and not by passive beamline elements. For this study the emittance is taken to be 4 10-9rm, which I understand to be a conservative figure for CEBAF. With the emittance fixed, we are free to adjust the beam optics to project the desired spot size at the collimator. A lower limit on the spot size is given by the distance to the focusing elements, where the beam must be correspondingly large in order to focus down to the desired spot size, and yet small enough to avoid aberrations from to higher-order field components. If the quads-to-collimator distance is 40m and the beam is limited to 2mm r.m.s. at the quads, this gives a minimum spot size around 100um at the collimator. As will be seen from the next section, the divergence properties of such a strongly focused beam has other effects on the coherent bremsstrahlung beam, which favour relaxing the spot size requirement somewhat. In the following acceptance plot I show the trend in collimation resolution as the spot size is adjusted within what I consider to be reasonable limits.
Note that the spot size sets an effective lower limit on the distance from radiator to collimator. Simply reducing this distance by some factor in combination with a similar reduction in collimator diameter does not produce a beam of the same quality. This is because the beam spot sets a length scale for the smearing of the collimator cutoff which does not scale with the distance to the radiator, for a given beam optics. In general one should ensure that the radiator-collimator separation is sufficient that the smallest collimators of interest are much larger than the virtual beam spot size. In Fig. 9 is plotted the photon beam profiles corresponding to the three choices of beam spot size in Fig. 8. For these studies the effects of multiple scattering were suppressed in order to accent the effects of beam emittance.
This smearing takes place because each electron in the beam has a slightly different incidence angle to a given crystal plane, and hence a different value of ql which relates x to angle (see Eq. 5). Small non-planarities of the crystal also contribute to this smearing, an effect known in crystallography as mosaic spread. Multiple scattering also makes a contribution. To lowest order, a shift in ql with no change in the angular acceptance results in a shift of the peak in x.
Following the same treatment as outlined above under multiple scattering, the mean-square beam divergence enters the crystal with a value 02 and exits with 02 + ms2. This angular spread translates into a smearing in x given by the function f(x,x').
which is applied to the x-distribution as the last step in the calculation described in this report.
This smearing conserves the integral of counts in the spectrum, and has the effect of smoothing out the sharp edges of the peaks. To isolate the effects of crystal quality, I suppress beam emittance and multiple scattering in Fig. 10, and explore a range of values for mosaic spread.
The energy scale is expanded in Fig. 10 to highlight the differences between the curves. A diamond crystal with mosaic spread below 100ur is considered to be ``almost perfect'' by Bilokon et.al.. If this is the case then for electron beams available at CEBAF, the mosaic spread dominates over beam emittance and multiple scattering in determining the degree of energy/angle smearing. In an experimental spectrum it is possible to read the mosaic spread directly from the sharpness of the high-energy edge of the coherent peak.
The second effect which limits the use of thin radiators is the radiation damage suffered by the crystal during exposure to high currents. Kaune et.al. report that at SLAC they found it necessary the remove the diamond and treat it with an annealing process after it had seen 2 1019 electrons. The beam spot had a diameter of roughly 6mm at the radiator, which leads to a radiator lifetime of 7 1019 e-/cm2 for the SLAC 80um diamond. Bilokon et.al. estimate 1020 e-/cm2 for diamonds several mm thick. Others have estimated longer lifetimes by a factor of 2-3. In this report I take 1020 e-/cm2 to be a conservative estimate.
In order to translate this into duration for an experiment, I need to make some assumptions about the experimental rate requirements. For this study I take the design parameters that were used for the current RadPhi experiment. Detector and trigger were designed assuming a tagged bremsstrahlung beam of 4GeV end-point, with a total tagging rate of 5 107 photons/s in the energy range [3.0,3.8]GeV. I I furthermore assume that this represents a limit for rates in the detector elements, so that one would not run with higher currents even if tagging accidentals were not a problem. Since most of the rate in the detectors comes from the untagged flux at low energies, a coherent radiator could enhance the rate of good triggers without appreciably affecting the overall rates in the detectors. Under these conditions and without any collimation, the rate within the interval [3.0,3.8]GeV could be increased from 5 107 to 17 107 photons/s. Of course without collimation the coherent energy peaks are very broad, and use of the tagger would be required, limiting the usable rates to a few 107. In the absence of collimation, there is no reason to use ultra-thin crystals, and thicknesses up to 1mm are basically equivalent. The lifetime of the radiator under these circumstances is on the order of years, much longer than the duration of an experiment.
To bypass the flux limitations of tagging, one must give up the energy information from the tagger. Depending on the requirements of the experiment, the narrow peaks from collimated coherent bremsstrahlung may provide sufficient resolution for the data analysis. The primary use of the tagger information in RadPhi is to limit the online trigger rate by rejecting events from untagged flux. At a tagging rate of 5 107 and a coincidence window of 15ns the rejection power has dropped to 50%. At this point one might eliminate tagging in favour of collimated coherent bremsstrahlung. One could drop the incoherent bremsstrahlung flux by a factor of two to maintain the same low-level trigger rates, and still enhance the rate of good triggers by a considerable fraction.
I consider two cases: a 100um radiator (thick) and a 10um radiator (thin). In the case of the thick radiator, the multiple scattering limits useful collimation to roughly 0.5 m/E, or 2mm at 30m. The r.m.s. width of the peak is 10%, covering fully the original tagging range [3.0,3.8]GeV. The photon flux in this window is 2.2 108, a factor of 4 increase. The spectrum is shown by the blue curve in Fig. 11. In the case of the thin radiator, the collimation is reduced to 0.25 m/E, or 1mm at 30m. The r.m.s. width of the peak is 4%, centered at 3.6GeV. The total photon flux is 2.7 108 within the interval [3.0,3.8]GeV, an increase by a factor greater than 5. The spectrum is shown by the red curve in Fig. 11.
The beam currents and consequent crystal lifetimes are shown in the following table. The calculations assumed a virtual beam spot of 200um radius at the collimator. Depending on the flexibility of the beamline optics, one may be able to extend the crystal lifetime by making the beam spot size at the radiator somewhat larger than this. Three different possibilities for the spot size at the radiator are shown in the table. If a mechanism were put in place to raster the crystal across the beam then the spot size would be irrelevant; the only limitation on the duration of an experiment with a single crystal would then be given by the figure in the last column, proportional to the area of the crystal.
target thickness | beam current | beam r.m.s. at radiator | usable lifetime | ||
---|---|---|---|---|---|
200um | 500um | 2mm | |||
100um | 200nA | 28h | 175h | 120d | 930 days/cm2 |
10um | 4.5uA | 1.2h | 7.8h | 125h | 41 days/cm2 |
If beams of comparable quality could be produced at CEBAF, it would provide an attractive alternative to tagged bremsstrahlung for a number of experiments with real photons. The question arises, what are the requirements for a beamline that could deliver such a beam of real photons using electron energies of 8+GeV? In this section I wish to enumerate the key elements in the setup that was used to obtain the excellent results shown above, as a prototype for what might be done at CEBAF.
The figure above, taken from the SLAC paper, is set up for 20GeV electrons. The collimator at 91m from the radiator is 2mm in diameter, or 0.42 m/E. The beamline elements are discussed below under five headings.
The elements responsible for producing the focus are located before the exit from the accelerator, and so are not shown in the above layout. The two magnets which are shown just upstream of the radiator are vertical and horizontal dipoles. They perform the small corrections necessary to keep the photon beam well centered on the collimator. The dipole located between the radiator and the collimator bends the electron beam into the dump.
The position of the photon beam at the collimator is sensed by a four-quadrant tungsten-pin shower-emission detector located downstream of the first segment (a few radiation lengths) of the primary collimator. The current from these pins is proportional to the local shower intensity. The entire collimator assembly could be remotely rotated through 180 degrees so that the four quadrants could be gain-balanced. The four signals were added and subtracted in two separate circuits to produce horizontal vertical displacement signals. These are read out by the computer for each pulse, and any necessary corrections automatically applied to the corrector power supplies.
With no online monitor of the photon spectrum, one would have to rely on periodic checks with the pair spectrometer to monitor degradation of the crystal. If it were present, the tagger could provide this function. At reduced beam current the tagger could be run in coincidence with a total absorption counter to measured the photon spectrum instead of having to instrument the focal plane of a pair spectrometer. During normal running, even though the tagger rates might be too high for counting, reading out the current from the tagging counters at reduced voltage would provide a fast online monitor of the status of the crystal and the properties of the coherent photon beam.
Richard Jones
University of Connecticut / Jefferson Lab