Search Results (17 total)

In this paper we report on large-scale high resolution simulations of beam dynamics in electron linacs for the next-generation x-ray free electron lasers (FELs). We describe key features of a parallel macroparticle simulation code including three-dimensional (3D) space-charge effects, short-range structure wakefields, coherent synchrotron radiation (CSR) wakefields, and treatment of radio-frequency (rf) accelerating cavities using maps obtained from axial field profiles. We present a study of the microbunching instability causing severe electron beam fragmentation in the longitudinal phase space which is a critical issue for future FELs. Using parameters for a proposed FEL linac at Lawrence Berkeley National Laboratory (LBNL), we show that a large number of macroparticles (beyond 100 million) is generally needed to control the numerical macroparticle shot noise and avoid overestimating the microbunching instability. We explore the effect of the longitudinal grid on simulation results. We also study the effect of initial uncorrelated energy spread on the final uncorrelated energy spread of the beam for the FEL linac.

The transport properties, including electrical resistivity (ρ), thermopower (S), and thermal conductivity (κ) of bulk metallic glass alloys Cu₆₄Zr_28.5Ti_7.1, Cu_62.3Zr_23.7Ti₁₀, Cu_60.6Zr_26.9Ti_12.5 and Cu_58.8Zr_26.2Ti₁₅ in the temperature range 10–300 K are reported. The temperature variations of electrical resistivity in these alloys, with a negative temperature coefficient, are found to be rather weak. These findings are consistent with the metallic glass nature of these compounds. It is observed that the electrical resistivity increases with increasing Ti concentration, ascribed to the enhancement of disorder with such a composition change. The magnetoresistivity of Cu₆₄Zr_28.5Ti_7.5 alloy decreases with increasing temperature and increases with increasing magnetic field, suggesting that the weak localization effect dominates the electrical transport. The temperature-dependent thermopower and thermal conductivity characteristics are nearly identical and weakly independent of compositions. It is noted that the observed κ(T) of Cu-Zr-Ti metallic glass alloys show notable similarities with the quasicrystalline system. There are two main features in S(T), a knee around 120 K and a plateau below 50 K, which represent the deviation from the expected linear behavior. A detailed theoretical analysis has suggested that the appearance of the knee is due to the electron-phonon interaction at low temperatures and the plateau is associated with low-energy excitations in glasses.

Theoretical and numerical studies of the transport in vacuum of multi-nC, multi-MeV electron beams are performed using several methods, including envelope models, a novel semianalytic approach using ellipsoidal shell decomposition, a modified electrostatic particle-in-cell method, and a point-to-point interaction model. The effects of space-charge forces on the longitudinal and transverse bunch properties are evaluated for various bunch lengths, energies, energy spreads, and charges. An evaluation of the various methods for studying space-charge effects in large energy spread, high charge beams is summarized. Examples are given for beam distributions typical of those generated by plasma-based accelerators. It is found that, for the highly correlated beams produced in the self-modulated regime, the high energy portion of the beam can gain significant energy while propagating in vacuum due to space-charge effects.

In this paper, we present a three-dimensional quasistatic model for high brightness beam dynamics simulation in rf/dc photoinjectors, rf linacs, and similar devices on parallel computers. In this model, electrostatic space-charge forces within a charged particle beam are calculated self-consistently at each time step by solving the three-dimensional Poisson equation in the beam frame and then transforming back to the laboratory frame. When the beam has a large energy spread, it is divided into a number of energy bins or slices so that the space-charge forces are calculated from the contribution of each bin and summed together. Image-charge effects from conducting photocathode are also included efficiently using a shifted-Green function method. For a beam with large aspect ratio, e.g., during emission, an integrated Green function method is used to solve the three-dimensional Poisson equation. Using this model, we studied beam transport in one Linac Coherent Light Sources photoinjector design through the first traveling wave linac with initial misalignment with respect to the accelerating axis.

In recent high luminosity colliders, the finite crossing angle scheme becomes popular to gain the luminosity with multibunch or long bunch operation. Success of the KEKB factory showed that the finite crossing angle scheme has no problem achieving beam-beam parameters up to 0.05. We have studied the beam-beam interactions with and without crossing angle toward higher luminosity. We discuss how the crossing angle affects the beam-beam parameter and luminosity in the present KEKB using computer simulations. The simulations showed that crab cavities, which realize the head-on collision effectively, can be expected to double the luminosity.

Beam-beam effects limit the luminosity of circular colliders. Once the bunch population exceeds a threshold, the luminosity increases at a slower rate. This phenomenon is called the beam-beam limit. Onset of the beam-beam limit has been analyzed with various simulation methods based on the weak-strong and strong-strong models. We have observed that an incoherent phenomenon is mainly concerned in the beam-beam limit. The simulation have shown that equilibrium distributions of the two colliding beams are distorted from Gaussians when the luminosity is limited. The beam-beam limit is estimated to be ξ∼0.1 for a B factory with damping time of several thousand turns.

We investigate the phenomenon of space-charge driven emittance growth in a three-dimensional mismatched anisotropic charged particle beam with relevance to high-intensity linear accelerators. The final emittance growth can be understood as a superposition of the contributions from the mismatch-induced halo formation and from the anisotropy-induced energy exchange. The averaged emittance growth per degree of freedom is bounded from above by the so-called “free energy limit” extended by the contributions from energy exchange. The partition of the growth into longitudinal or transverse is, however, a strong function of the tune ratio including the possibility that an initially equipartitioned beam is even driven substantially away from equipartition. The growth of the beam halo extent is dominated by the effect of mismatch, whereas anisotropy itself generates practically no halo.

In this paper, we report on the halo formation and emittance growth driven by a parametric resonance during beam-beam collisions with optically mismatched injection. In the regime of the weak-strong beam-beam interaction, if two beams have the same machine tunes, on-axis head-on collisions between a mismatched strong beam and a weak beam will not cause the formation of halo. However, if the two beams collide with an initial offset, the beam-beam force from the mismatched strong beam can cause halo formation and emittance growth in the weak beam. Meanwhile, if two beams have different machine tunes, for opposite charged colliding beams, when the machine tune of the weak beam is smaller than that of the strong beam, there is emittance growth in the weak beam. When the machine tune of the weak beam is larger than that of the strong beam, there is little emittance growth. In the regime of the strong-strong beam-beam interaction, halo is formed in both beams even when the two beams collide head-on on the axis with equal machine tunes. This puts a strong requirement for a good beam match during the injection to colliders in order to avoid emittance growth if there are no other mechanisms after the injection to decohere the mismatched envelope oscillation.

The theory and simulation of coherent resonant coupling due to space charge in coasting or bunched anisotropic equilibrium beams is presented. Our work confirms that analytical results on coherent oscillations and instabilities of anisotropic KV (Kapchinskij-Vladimirskij) distributions are a valid tool to interpret the findings from 2D and 3D self-consistent particle-in-cell simulations for both KV and waterbag distributions. With reference to rings we discuss space charge coherent tune shifts up to fourth order and introduce a coherent coupled mode coefficient, which enables us to resolve the issue of KV anomalies by relating them to negative energy modes. The second emphasis of this study is with reference to linacs and a detailed discussion of “stability charts” describing resonant regions where approach to equipartition may occur.

We report macroparticle simulations for comparison with measured results from a proton beam halo experiment in a 52-quadrupole periodic-focusing channel. An important issue is that the input phase-space distribution is not experimentally known. Three different initial distributions with different shapes predict different beam profiles in the transport system. Simulations have been fairly successful in reproducing the core of the measured matched-beam profiles and the trend of emittance growth as a function of the mismatch factor, but underestimate the growth rate of halo and emittance for mismatched beams. In this study, we find that knowledge of the Courant-Snyder parameters and emittances of the input beam is not sufficient for reliable prediction of the halo. Input distributions with greater population in the tails produce larger rates of emittance growth, a result that is qualitatively consistent with the particle-core model of halo formation in mismatched beams.

We present results from an experimental study of the beam halo in a high-current 6.7-MeV proton beam propagating through a 52-quadrupole periodic-focusing channel. The gradients of the first four quadrupoles were independently adjusted to match or mismatch the injected beam. Emittances and beamwidths were obtained from measured profiles for comparisons with maximum emittance-growth predictions of a free-energy model and maximum halo-amplitude predictions of a particle-core model. The experimental results support both models and the present theoretical picture of halo formation.

In this paper we present a new approach, based on a shifted Green function, to evaluate the electromagnetic field in a simulation of colliding beams. Unlike a conventional particle-mesh code, we use a method in which the computational mesh covers only the largest of the two colliding beams. This allows us to study long-range parasitic collisions accurately and efficiently. We have implemented this algorithm in a new parallel strong-strong beam-beam simulation code. As an application, we present a study of a beam sweeping scheme for the LBNL luminosity monitor of the Large Hadron Collider.

Macroparticle simulation plays an important role in modern accelerator design and operation. Most linear rf accelerators have been designed based on macroparticle simulations using longitudinal position as the independent variable. In this paper, we have done a systematic comparison between using longitudinal position as the independent variable and using time as the independent variable in macroparticle simulations. We have found that, for an rms-matched beam, the maximum relative moment difference for second, fourth moments and beam maximum amplitudes between these two types of simulations is 0.25% in a 10 m reference transport system with physical parameters similar to the Spallation Neutron Source linac design. The maximum z-to- t transform error in the space-charge force calculation of the position dependent simulation is about 0.1% in such a system. This might cause a several percent error in a complete simulation of a linac with a length of hundreds of meters. Furthermore, the error may be several times larger in simulations of mismatched beams. However, if such errors are acceptable to the linac designer, then one is justified in using position dependent macroparticle simulations in this type of linac design application.

Energy exchange between the longitudinal and transverse degrees of freedom of nonequipartitioned bunched beams (non-neutral plasmas) is investigated by means of 3D simulation. It is found that collective instability may lead to energy transfer in the direction of equipartition, without full progression to it, in certain bounded regions of parameter space where internal resonance conditions are satisfied, in good agreement with stability charts from an earlier derived 2D Vlasov analysis. Nonequipartitioned stable equilibria, however, exist in relatively wide regimes of parameter space. This provides evidence that such regimes may be safely used in the design of future high-intensity linacs.

A stochastic leapfrog algorithm for the numerical integration of Brownian motion stochastic differential equations with multiplicative noise is proposed and tested. The algorithm has a second-order convergence of moments in a finite time interval and requires the sampling of only one uniformly distributed random variable per time step. The noise may be white or colored. We apply the algorithm to a study of the approach towards equilibrium of an oscillator coupled nonlinearly to a heat bath and investigate the effect of the multiplicative noise (arising from the nonlinear coupling) on the relaxation time. This allows us to test the regime of validity of the energy-envelope approximation method.

In this paper we present a study of beam halo based on a three-dimensional particle-core model of an ellipsoidal bunched beam in a constant focusing channel including the effects of nonlinear rf focusing. For an initially mismatched beam, three linear envelope modes—a high frequency mode, a low frequency mode, and a quadrupole mode—are identified for an azimuthally symmetric bunched beam. The high frequency mode has three components all in phase; the low frequency mode has the transverse components in phase and the longitudinal component 180° out of phase; the quadrupole mode has no longitudinal component, and the two transverse components in the mode are 180° out of phase. We also study the case of an ellipsoidal bunched beam without azimuthal symmetry and find that the high frequency mode and the low frequency mode are still present but the quadrupole mode is replaced by a new mode with transverse components 180° out of phase and a nonzero longitudinal component. Previous studies, which generally addressed the situation where the longitudinal-to-transverse focusing strength is roughly 0.6 or less, conclude that the oscillation of the high frequency mode is predominantly transverse, and that of the low frequency mode is predominantly longitudinal. In this paper we present a systematic study of the features of the modes as a function of the longitudinal-to-transverse focusing strength ratio. We find that, when the ratio is greater than unity, the high frequency mode may contain a significant longitudinal component. Thus, excitation of the high frequency mode in this situation can be responsible for the formation of longitudinal beam halo. Furthermore, while previous studies have observed halo amplitudes roughly 2–3 times the matched beam edge, for the present parameters we observe much larger amplitudes (5 times or more). This is due to the fact that the longitudinal-to-transverse focusing ratio used here is greater than that of previous studies. The finding of large transverse halo amplitude can have significant impact on the design of high-intensity ion accelerators where the longitudinal-to-transverse focusing ratio is slightly greater than unity in some parts of the linac.