next up previous contents index
Next: Laser Up: Physics and Numerical Methods Previous: Luminosity Integration

Beam Field

    One of the basic assumptions of CAIN is that the most particles in the beams have high energy and are almost either right- or left-going. This assumption leads to the following facts: In contrast to ABEL, CAIN does not assume that all the particles have the above property. Some particles may have low energies and large angles with respect to the s-axis. CAIN will work, unless the sum of their weight becomes a significant fraction of the beam. The equation of motion under the Lorentz force is integrated with possible low energies and large angles taken into account.

The calculation of the beam electric field is done in the following way. First, cut the right(left)-going particles into longitudinal slices (the width tex2html_wrap_inline8608 is defined by the parameter Smesh).    Within each slice the following Poisson equation  is solved.
where m is the electron mass in units of eV/ctex2html_wrap_inline7792, tex2html_wrap_inline9804 the classical electron radius in meters, tex2html_wrap_inline9806 is the charge (divided by the elementary charge) per unit transverse area, then tex2html_wrap_inline9784 is given in units of V/m.

For each slice and for each of right- and left going beams, a region tex2html_wrap_inline9810 is selected, where tex2html_wrap_inline9812 is the center-of-mass and tex2html_wrap_inline8274 is the width determined by the input parameters. The field created by the particles outside this region is ignored. Let us name this region [O].

  In the region [O], the Poisson equation  is solved using the FFT. Eq.(42) can formally be solved as
Divide this region by tex2html_wrap_inline9816 grid. Within each cell (i,j) (tex2html_wrap_inline9820, tex2html_wrap_inline9822), the the density tex2html_wrap_inline9806 is approximated by tex2html_wrap_inline9826, where tex2html_wrap_inline9828, tex2html_wrap_inline9830, and tex2html_wrap_inline9832 is the total charge in the cell:
where (tex2html_wrap_inline9834) is the cell center coordinate. Then, eq.(43) becomes a sum over the cells.
The kernel matrix tex2html_wrap_inline9836 has to be calculated by taking average over the source cell:
This averaging is important when tex2html_wrap_inline9838 is far from unity. The convolution in eq.(45) can be done efficiently by using FFT.

However, if we apply FFT for the finite region tex2html_wrap_inline8274 instead of the infinite region in eq.(43), we would be assuming a periodic charge distribution, i.e., the charge distribution in tex2html_wrap_inline8274 is infinitely repeated. To avoid this problem, we use the following trick. First double the region to tex2html_wrap_inline9844 by padding zero in the extended region and carry out FFT. This still means a periodic charge distribution as depicted in Fig.3. However, if we use the kernel matrix with zero padded   in the extended region (tex2html_wrap_inline9846 if tex2html_wrap_inline9848 or tex2html_wrap_inline9850), the field due to the ghost charges will never reach the real charge region because their horizontal(vertical) distance is larger than tex2html_wrap_inline8198 (tex2html_wrap_inline8204). Thus, the potential tex2html_wrap_inline9856 in the region tex2html_wrap_inline8274 is calculated correctly although incorrect in the extended region.

Figure 3: Doubled region for FFT. The solid frame indicates the doubled region for FFT and its left-bottom quadrant is the charge region tex2html_wrap_inline9860. The region hatched by solid lines is the real charge region and that by dotted lines the ghost charge due to the periodicity of Fourier transformation.

The obtained values of the potential are those at cell centers. They are interpolated by 2-dimensional cubic spline and differentiated to get tex2html_wrap_inline9862 and tex2html_wrap_inline9864.

When a charged particle gets out of the mesh region, the field created by it is ignored in CAIN. However, the force by the other beam is taken into account even if the particle is outside the mesh region of the other beam. To this end, CAIN  adopts three methods, namely, [A] direct Coulomb force by the charge distribution in the mesh, [B] harmonic expansion in polar coordinate, and [C] harmonic expansion in elliptic coordinate.

Let tex2html_wrap_inline8274 be the total width of the mesh region. If it is close to a square, or more precisely, if tex2html_wrap_inline9868, the whole region is divided into three regions [O],[A],[B], as depicted in Fig.4a. If the mesh region is far from square, the whole region is divided into four, [O],[A],[B],[C], as in Fig.4b. In the region [O] the mesh is used for calculating the field. In other regions, the methods mentioned above are used.

Figure 4: Regions for calculating the beam field

Since the sum is time consuming, this is used only in region [A], where two other methods fail to converge. The method is trivial and given by
However, this formula is not accurate when the bin size ratio tex2html_wrap_inline9838 is far from unity. It is needed to take average over a bin when the bin is close to the field point (x,y). CAIN makes a table for the Coulomb force by a bin (tex2html_wrap_inline9874) for faster computation.

  In the region [B] the following formula is used.

Here, tex2html_wrap_inline9876 is arbitrary (introduced for avoiding overflow/underflow). The formula is valid for tex2html_wrap_inline9878, where tex2html_wrap_inline9880 = tex2html_wrap_inline9882 is the maximum radius of the mesh region.

  When tex2html_wrap_inline9884 (otherwise, exchange x and y), the elliptic coordinate (u,v) defined by
is used. Here, f is chosen as
The maximum of the radial-like coordinate u in the mesh region is
which is taken at the four corners.

Then, the expansion of tex2html_wrap_inline9856 is

Actually, there is a finite relation between tex2html_wrap_inline9898 and tex2html_wrap_inline9900:
The formula converges if tex2html_wrap_inline9902, which corresponds to the region [C] (and [B]) in Fig.4b. The truncation of the series is defined by the operand NMOM of the command BBFIELD (common to the two types of expansions for simplicity).

next up previous contents index
Next: Laser Up: Physics and Numerical Methods Previous: Luminosity Integration

Toshiaki Tauchi
Thu Dec 3 17:27:26 JST 1998