Development of 2D particle-in-cell code to simulate high current, low energy beam in a beam transport system
|
|
- Brenda Armstrong
- 5 years ago
- Views:
Transcription
1 PRAMANA c Indian Academy of Sciences Vol. 9, No. journal of October 7 physics pp Development of D particle-in-cell code to simulate high current, low energy beam in a beam transport system S C L SRIVASTAVA, S V L S RAO and P SINGH Nuclear Physics Division, Bhabha Atomic Research Centre, Mumbai 85, India shashi uv@yahoo.com MS received 1 March 7; revised 1 May 7; accepted 5 July 7 Abstract. A code for D space-charge dominated beam dynamics study in beam transport lines is developed. The code is used for particle-in-cell (PIC) simulation of z-uniform beam in a channel containing solenoids and drift space. It can also simulate a transport line where quadrupoles are used for focusing the beam. Numerical techniques as well as the results of beam dynamics studies are presented in the paper. Keywords. Space-charge; particle-in-cell; beam dynamics; Poisson s equation; solenoids; quadrupole magnets. PACS Nos 9.7.Bd; 5.5.Rr; 1..cv 1. Introduction Newly proposed accelerators with applications to nuclear waste transmutation and spallation neutron sources require high intensity Linacs. To develop the technology of high intensity Linacs for our ADS programme, a MeV, 3 ma CW proton accelerator is being built at the Bhabha Atomic Research Centre, Mumbai. It consists of 5 kev ECR ion-source, low energy beam transport line (LEBT), 3 MeV radio frequency quadrupole (RFQ), medium energy beam transport line (MEBT) and MeV drift tube Linac (DTL) [1]. In the low energy section of such accelerators, beams are strongly subjected to the Coulomb repulsion and understanding the behaviour of such space-charge dominated beams is a challenging task. This requires a careful control of particle dynamics by suitably incorporating space-charge forces in the beam. Several codes (TRACED [], TRANSPORT [3] etc.) are available to study the beam dynamics and to match the beam parameters between two structures (for example, ion-source and RFQ). As these codes are formulated by linearizing the space-charge forces and equations of motion, they are expected to be accurate in the linear regime. These codes cannot therefore simulate the nonlinear space-charge effects in beam dynamics. There are numerous particle-in-cell (PIC) codes (PARMILA [], PARMTEQ [5], WARP3D [], BEAMPATH [7] etc.) 551
2 S C L Srivastava, S V L S Rao and P Singh developed in the field of accelerator physics to study dynamics of space-charge dominated beams. To our knowledge, only BEAMPATH [7] simulates the continuous beam in transport lines using PIC method. In this paper, we present an algorithm that forms the basis of the code, which has been developed to study various aspects of space-charge dominated continuous beam in a LEBT. In addition to tracking the particles, generation of different kinds of beam distributions is an integral part of this code. In future, we plan to incorporate, in order to study space-charge compensation, collisions of protons with neutral gas atoms and molecules using Monte Carlo technique in the code.. Transport of space-charge dominated beam Consider the propagation of an intense charged particle beam through a beam transport system consisting of different focusing elements. For the beam dynamics study, the PIC method [8] is used. The beam is represented as a combination of a large number of macroparticles with the same charge-to-mass ratio as that of the real beam. The simulation is performed in D phase space of particle transverse positions x, y and transverse velocities v x, v y. The macroparticles following certain kind of distribution will be tracked in the presence of external field and the field exerted by the beam. To calculate field exerted by the beam (self-field), we need to solve the Poisson equation and from that electric fields are calculated. Since motion of particles is influenced by the self-fields which depend on their spatial distribution inside the beam, we need to solve the Poisson equation at each and every step and this problem has to be solved self-consistently. During simulations the macroparticle is lost if it touches the boundaries. From the single particle Hamiltonian, equations of motion can be derived and in general they can be written as follows: dx dt = v x, dy dt = v y, dv x dt = q m (E + v B) x, dv y dt = q m (E + v B) y. (1) Integration is performed with fixed time step δt. In (1), the electric field is a combination of external and space-charge fields, while magnetic field is external only (at low energy, self-magnetic field is very small as compared to electric field). In the following sections, the numerical algorithms of the code and results are discussed. 3. Beam distribution generator Initially, the particles will be distributed in the D transverse phase space according to a distribution function. For the generation of particle distribution in transverse 55 Pramana J. Phys., Vol. 9, No., October 7
3 Development of D PIC code D phase space x, x, y, y, consider a class of distributions with elliptical symmetry. The distribution function n(x, x, y, y ) depends on total emittance ɛ which has a meaning of radius-vector in D phase space: n(x, x, y, y ) = dn dxdx = n(ɛ), () dydy ɛ = A x + ca y, (3) where c is the ratio of beam emittances [9]. For any value of c except 1 the area of the D projection (ellipse) of (3) will be unequal. Parameters A x, A y describe a family of ellipses. A x = A y = ( βx x + α xx βx ) + ( x βx ), () ( ) ( ) βy y + α yy y +, (5) βy βy where α x,y, β x,y and γ x,y are Courant Snyder parameters in the x and y planes. The equation ɛ = constant describes a hyperellipsoid surface in the phase space x, x, y, y. As the distribution function depends on ɛ, the phase space density, n, will be constant on one hyperellipsoid surface while it will vary from one surface to another. The distribution function is normalized under the following condition: n(x, x, y, y )dxdx dydy = 1. () For generating the phase space distributions, we have to generate the N values of x, x, y, y which correspond to a given distribution function n(ɛ). To generate various distributions, we need to know the distribution of ɛ, g(ɛ) = dn(ɛ). (7) dɛ Let us transform x, x, y, y to new coordinates r, θ, φ, ψ ( θ π, φ π, ψ π): βx x + α xx βx = r sin θ cos φ sin ψ, (8) x βx = r sin θ sin φ sin ψ, (9) c ( βy y + α yy βy ) = r cos θ sin ψ, (1) ( ) y c = r cos ψ, r = ɛ. (11) βy Pramana J. Phys., Vol. 9, No., October 7 553
4 S C L Srivastava, S V L S Rao and P Singh The phase space element transforms as dxdx dydy = r3 c sin ψ sin θdrdθdφdψ. (1) Then the number of particles in the phase space element is dn(r, θ, φ, ψ) = n(r) r3 c sin ψ sin θdrdθdφdψ. (13) Integration of (13) over angle variables gives the number of particles dn as a function of r, g(r) = dn(r) dr = π r3 n(r). (1) c The algorithm for generation of different distributions is based on the above equations and it is given below: 1. To simulate the distribution g(r) we will use the inverse transform method in which we take the integral distribution defined by (15) to find r = r(g) under the assumption that the values of G are uniformly distributed in the interval [,1]: G(r) = r g(r )dr. (15). For each value of r, two random numbers A x and A y are chosen such that they satisfy (3) and then the points x, x, y, y are calculated by the following equations: x = A x βx cos p, (1) ( x = A x α x cos p + sin p ), βx βx (17) y = A y βy cos q, (18) ( ) y = A y α y cos q + sin q, βy βy (19) where p and q are two random numbers which are uniformly distributed in [, π]. The steps 1 and will be repeated N times for generating N particles. The generated particles will follow the distribution n(ɛ). The above method is implemented in the program, which presently generates D Kapchinskij Vladimirskij (KV), waterbag and parabolic distributions. This needs emittance, and Courant Snyder parameters of the beam to generate the specified distributions. The definition of different distributions is given in table 1. The (x, x ) projections of all the three distributions are shown in figure 1. The representation of the emittance was found to be better than 1% in all the cases for 1 and more particles. 55 Pramana J. Phys., Vol. 9, No., October 7
5 Development of D PIC code Table 1. Definition of different phase space distributions. Distributions Definition c rδ(r π r KV c π r Waterbag c (1 π r Parabolic 5 3 r ) r Waterbag Distribution 3 5 KV Distribution r ) 5 (a) (b) 5 3 Parabolic Distribution (c) Figure 1. Projection of (a) KV, (b) waterbag and (c) parabolic distribution on x x plane for Twiss parameters βx = βy =.78 cm/rad, γx = γy =.171 rad/cm, αx = αy = 1.8, emittance ²x = ²y =.π cm mrad and c = 1.. Poisson solver The space-charge potential of the beam, U, for an instantaneous space-charge density distribution, ρ, is calculated from the solution of Poisson s equation U ρ(x, y) U + =. x y ² () For a z-uniform beam, this problem reduces to a D problem in x y coordinates. First, we distribute the space-charge of macroparticles among grid nodes, then solve Pramana J. Phys., Vol. 9, No., October 7 555
6 S C L Srivastava, S V L S Rao and P Singh the Poisson s equation on grid and finally after calculation of field on grid nodes, scale it at macroparticle positions. The simulation region is divided into uniform rectangular meshes. Charge of every particle with coordinates (x n, y n ) is distributed among the nearest four nodes utilizing area weighting method [8]. The charge density at node point, ρ ij, is given by ρ ij = N n=1 ( ρ xy 1 x ) ( n x i 1 y ) n y i, (1) h x h y where ρ xy is the space-charge density of an individual particle and h x, h y are the mesh sizes. ρ xy is related to the beam current I and to the speed of the beam in the z-direction ρ xy = I v z Nh x h y. () The Poisson s solver is implemented such that it can handle either Dirichlet or Neumann boundary condition in one direction and mixed boundary condition in the other direction..1 For Dirichlet and mixed boundary condition Suppose that U(x, y) satisfies mixed boundary conditions in the x-direction, i.e., a l U(x, y) + b l U(x, y) x = g l (y) (3) at x = x l and a similar equation for x = x h with new constants. Here, a l, b l, etc. are known constants, whereas g l is a known function of y. Furthermore, suppose that U(x, y) satisfies the following simple Dirichlet boundary conditions in the y- direction: U(x, ) = U(x, L) =. () Let us write U(x, y) as a Fourier series in the y-direction: U(x, y) = u j (x) sin(jπy/l). (5) j= The functions sin(jπy/l) are orthogonal, and form a complete set in the interval y [, L]. In fact, L L sin(jπy/l) sin(kπy/l)dy = δ jk. () Thus, we can write the source term as 55 Pramana J. Phys., Vol. 9, No., October 7
7 Development of D PIC code ρ(x, y) = ϱ j (x) sin(jπy/l), (7) j= where ϱ j (x) = L L ρ(x, y) sin(jπy/l) dy. (8) Furthermore, the boundary conditions in the x-direction become a l u j (x) + b l du j (x) dx = Γ l j (9) at x = x l and a similar equation for x = x h with new constants, where Γ l j = L L g l (y) sin(jπy/l) dy. (3) Using (5) and (7) in () and equating the coefficients of sin(jπy/l) (since these functions are orthogonal), we obtain d u j (x) dx j π L u j(x) = ϱ j (x), (31) for j =. Now, we can discretize the problem in the y-direction by truncating our Fourier expansion, i.e., by only solving the above equations for j = J, rather than j =. This is essentially equivalent to discretization in the y- direction on the equally-spaced grid-points y j = jl/j. The problem is discretized in the x-direction by dividing the domain into equal segments. Thus, we obtain u i 1,j ( + j κ ) u i,j + u i+1,j = ϱ i,j (h x ), (3) for i = 1 N x and j = J. Here, u i,j u j (x i ), ϱ i,j ϱ j (x i ) and κ = πh x /L. The boundary conditions (9) discretize to give u,j = Γ ljh x b l u 1,j a l h x b l, (33) u Nx+1,j = Γ hjh x + b h u Nx,j a h h x + b h, (3) for j = J. Equations (3) (3) constitute a set of uncoupled tridiagonal matrix equations (with one equation for each j value). These equations can be inverted, to give u i,j. Finally, U(x i, y j ) values can be reconstructed from (5).. For Neumann and mixed boundary condition If there is Neumann boundary condition in the y-direction as Pramana J. Phys., Vol. 9, No., October 7 557
8 S C L Srivastava, S V L S Rao and P Singh U(x, y = ) y = U(x, y = L) y then we can express U(x, y) in the form U(x, y) = =, (35) u j (x) cos(jπy/l) (3) j= so that it automatically satisfies the boundary conditions in the y-direction. Likewise, we can write the source term ρ(x, y) as where ρ(x, y) = ϱ j (x) cos(jπy/l), (37) j= since ϱ j (x) = L L ρ(x, y) cos(jπy/l) dy, (38) L L cos(jπy/l) cos(kπy/l) dy = δ jk. (39) Finally, the boundary conditions in the x-direction become a l u j (x) + b l du j (x) dx = Γ l j, () at x = x l and a similar equation for x = x h with new constants, where Γ l j = L L g l (y) cos(jπy/l) dy. (1) Note, however, that the factor in front of the integrals in (38) and (1) takes the special value 1/L for the j = harmonic. As before, we truncate the Fourier expansion in the y-direction, and discretize in the x-direction, to obtain the set of tridiagonal matrix equations and solve as before. Fast Fourier transform is implemented using FFTW libraries [1]. Electric field on the grid nodes are calculated using E i,j = U i,j ; () electric field at the particle position is calculated from field values at grid nodes using the same area weight method (1), as for charge. Since these fields are in a moving frame of reference, the perpendicular component of the field in laboratory frame will be modified by the relativistic factor γ. The parallel component, however, will remain unaffected. (For low energies, the correction is insignificant because γ is close to unity.) 558 Pramana J. Phys., Vol. 9, No., October 7
9 Development of D PIC code 5. Integration of particle trajectories The particle trajectories are integrated using leap-frog method to preserve the kinematic time reversal symmetry in the simulation. In the most general case, integration of particle trajectories in this code is carried out through the following steps: 1. The particle performs a half-step acceleration in an electric field: v n = v n + qδt m E n. (3). The particle velocity undergoes rotation in the magnetic field. This is implemented by employing Boris scheme [11]: v n = v n + v n T, () T = q B δt m, (5) v n+1 = v n + v n s, () s = T 1 + T. (7) 3. The particle performs again half-step acceleration in an electric field: v n+1 = v n+1 + qδt m E n. (8). Finally, particles are advanced with a velocity, v n+1 : x n+3/ = x n+1/ + v n+1 δt. (9) In all the above equations v and x possess only x- and y-components (not z) while B has its usual meaning.. Focusing fields In this section, we describe the implementation of two focusing devices in the code, namely, solenoids and quadrupoles..1 Solenoid In this code the hard edge model of solenoid is implemented, i.e. the effect of fringe field is considered at a point. To see the effect of fringe field on the motion of charged particle, we take as Gaussian surface a truncated square circular cylinder of radius R coaxial with the solenoid. One end is well inside the solenoid, where Pramana J. Phys., Vol. 9, No., October 7 559
10 S C L Srivastava, S V L S Rao and P Singh Gaussian Surface(Truncated Cylinder) B Circular Area Cylindrical Area Figure. Side view of Gaussian surface in the solenoid. the field is B in the axial direction, and the other is well outside, where the field is zero as in figure. Then over this surface, B d S =. (5) So πr B = πr B r dz. (51) Simply put, all the flux leaving the cylinder cap inside the solenoid must have come in through the sides of the cylinder. This means that there is an irreducible amount of integrated radial magnetic field through the fringe field region: B r dz = B R/. (5) The effect of the fringe field on a particle of charge q traveling with axial velocity component v z = dz/dt is to kick it sideways with a force or F θ = qv z B r = dp θ /dt (53) dp θ = qb r dz, (5) p θ = q B r dz = qb R/. (55) The negative (positive) denotes the entrance (exit) of the beam in the solenoid field. There is no change in p r and p z. Let us transform these changes in Cartesian coordinates. As p x p y p z = cos φ sin φ sin φ cos φ 1 p θ, (5) 5 Pramana J. Phys., Vol. 9, No., October 7
11 Development of D PIC code so p x = ± qb y, (57) p y = qb x. (58) At the beginning and at the end of the solenoid, the transverse velocity changes according to (57) and (58) but there will be no change in position as we have represented the field as a Dirac-delta distribution. Inside the solenoid, magnetic field is constant and so the equations of motion are dv x dt = q m (E x + v y B ), (59) dv y dt = q m (E y v x B ). (). Quadrupole The field of quadrupoles [9] is given by the following equations: B x = B y a and B y = B x a, (1) where a is the aperture of the quadrupole. The equations of motion (1) are now augmented with the magnetic field values given by (1). 7. Code validation and low energy beam transport simulation To validate the code, we have chosen to study the evolution of round uniform beam following KV distribution in the drift space. As the space-charge forces are linear, the envelope equation can be written as [9] r m = ɛ rm 3 + I I β 3 γ 3, () r m where β and γ are relativistic parameters. This equation is numerically evaluated for a 5 kev, 3 ma proton beam and then evolution is compared with that obtained by the PIC code. The difference in the two solutions (defined as difference of two envelope radius divided by envelope radius calculated from ()) is less than 1. The two solutions are plotted in figure 3 and are found to be in good agreement. We have used our code to simulate the LEBT for the MeV, 3 ma proton accelerator being built in BARC [1]. A comparison between the parameters obtained by TRACED and our code is presented. The LEBT is designed using TRACED with the following input Twiss parameters, ɛ =.π cm mrad, β x = β y =.78 Pramana J. Phys., Vol. 9, No., October 7 51
12 S C L Srivastava, S V L S Rao and P Singh..5 Envelope Radius (m)..3. PIC KV Eq z (m) Figure 3. Envelope radius as a function of distance traveled in drift space. cm/rad, γ x = γ y =.171 rad/cm and α x = α y = 1.8 and the output parameters are found to be ɛ =.π cm mrad, α x = α y = 1.8, β x = β y =.3 cm/rad and γ x = γ y =.59 rad/cm. The lattice parameters of the LEBT is given in table. As TRACED is an envelope tracking code, a constant emittance is obvious. With the same input parameters and lattice, we obtain the output parameters as ɛ x =.31π cm mrad, ɛ y =.π cm mrad, α x = 1.83, α y = 1.9, β x =.1 cm/rad, β y =.87 cm/rad and γ x =.9 rad/cm, γ y =.7 rad/cm. The particle trajectories calculated using these parameters are shown in figure. As the numerical simulations are sensitive to mesh size, we found that for this problem there is no significant change in the output for mesh size greater than All the calculations are therefore performed on a mesh size of Effect of different distributions on the beam envelope To bring out the capabilities of simulating the non-linear space-charge fields with this code, we have simulated the behaviour of beam with different particle distributions in the designed LEBT. Since waterbag and parabolic distributions are non-uniform in D (x y) projections, they will give rise to non-linear space-charge field. The effect of non-linearity is clearly seen in the form of increase in the output emittance values. They are found to be ɛ x =.31,.7 and.33π cm mrad, ɛ y =.,.5 and.31π cm mrad for KV, waterbag and parabolic distribution respectively in each plane. The phase space projection on x x at the input and output are shown in figure 5. The initial Twiss parameters and emittance values are taken to be the same as mentioned in 7. Furthermore, the effect of various distributions is clearly visible on the maximum beam radius as shown in figure. 5 Pramana J. Phys., Vol. 9, No., October 7
13 Development of D PIC code x (cm) Distance (cm) Figure. Trajectories of particles in LEBT in x z plane. Table. LEBT (using solenoids and drift) parameters. Element Length (cm) Drift Solenoid Drift Solenoid Drift Strength (T) Table 3. LEBT (using quadrupole and drift) parameters. Element Drift Quadrupole Drift Quadrupole Drift Quadrupole Drift Quadrupole Drift Length (cm) Strength (T/m) For comparison, we designed a beam transport system using magnetic quadrupole lenses instead of solenoids. The input/output parameters are chosen the same as in case of LEBT containing solenoids and the designed LEBT (beam) optical parameters using TRACED are given in table 3. On simulation with our code with these parameters for KV distribution, the output turns out to be: ²x =.1π cm mrad, ²y =.π cm mrad, αx =.93, Pramana J. Phys., Vol. 9, No., October 7 53
14 S C L Srivastava, S V L S Rao and P Singh Input KV Distribution Output (KV) 1 15 (a) Input Waterbag Distribution 5 Output (Waterbag) (c) (d) Input Parabolic Distribution 1 Output (parabolic) (b) (e) (f) Figure 5. Phase space projections of different distributions: (a) KV, (c) waterbag and (e) parabolic in x x plane at the input while (b), (d) and (f ) are at the output of LEBT in the corresponding order. αy = 1.1, βx = cm/rad, βy = 3.57 cm/rad, γx =.593 rad/cm and γy =.581 rad/cm. The emittances were poorer in cases of waterbag (²x =.59, ²y =.98), and parabolic (²x =.3, ²y =.39π cm mrad) distributions. The poorer emittances were due to non-linearity in the space-charge field. The maximum beam radius (1.1 cm) and beam radius at the end of the LEBT (.3 cm) are found to be in agreement with the values obtained by TRACED (1 cm,.3 cm respectively). The beam radius as a function of z is shown in figure 7. 5 Pramana J. Phys., Vol. 9, No., October 7
15 Development of D PIC code 7 5 x (cm) 3 1 KV distribution Waterbag distribution Parabolic distribution z (cm) Figure. Effect of different distributions on maximum envelope radius. 1 Beam radius (cm) 1 8 x envelope y envelope z (m) Figure 7. Variation of beam radius in the two transverse planes as a function of longitudinal distance in the LEBT. 9. Conclusion The D PIC program for studying the beam dynamics in LEBT system has been developed. The results are found to be in good agreement with theoretical analysis. The algorithm implemented has good accuracy and takes a reasonable amount of computer time. It is planned to incorporate collisions of protons with neutral gas atoms and molecules for charge compensation estimation using Monte Carlo technique. Acknowledgements We thank Drs V C Sahni, S Kailas and R K Choudhury for their keen interest in this work. We also thank Dr Kartik Patel, Rajni Pande, Tushima Basak and Shweta Roy for many useful discussions. Pramana J. Phys., Vol. 9, No., October 7 55
16 References S C L Srivastava, S V L S Rao and P Singh [1] P Singh et al, Pramana J. Phys. 8, 331 (7) [] K R Crandall and D R Rusthoi, Computer program TRACED, LANL Report, LA- 135 (198) [3] K L Brown and S K Howry, TRANSPORT/3, a computer program for design charged particle beam transport system, SLAC report No. 91 (197) [] B Austin, T W Edwards, J E O Meara, M L Palmer, D A Swenson and D E Young, MURA Report No. 713 (195) D A Swenson and J Stovall, LANL Internal Report MP-3-19 (198) J P Boicourt and Charles R Eminhizer (eds), AIP Conference Proceedings (New York, 1988) vol. 177, p. 1 [5] K R Crandall, T P Wangler and Charles R Eminhizer (eds), AIP Conference Proceedings (New York, 1988) vol. 177, p. [] A Friedman, D P Grote and I Haber, Phys. Fluids B, 3 (199) [7] Y K Batygin, Nucl. Instrum. Methods in Phys. Res. A539, 55 (5) [8] C K Birdsall and A B Langdon, Plasma physics via computer simulation (IOP, Bristol 1991) [9] M Reiser, Theory and design of charged particle beams (John Wiley & Sons, New York, 199) [1] Matteo Frigo and Steven G Johnson, The design and implementation of FFTW3, Proc. IEEE 93(), 1 (5) [11] J P Boris, Relativistic plasma simulation-optimization of a hybrid code, Proc. th Conf. Num. Sim. Plasmas, Naval Res. Lab., November Pramana J. Phys., Vol. 9, No., October 7
Developmental Studies of High Current Proton Linac for ADS Program
1 Developmental Studies of High Current Proton Linac for ADS Program P. Singh, S.V.L.S. Rao, Rajni Pande, Shweta Roy, Rajesh Kumar, Piyush Jain, P.K. Nema, S. Krishnagopal, R.K. Choudhury, S. Kailas, V.C.
More informationPhysics design of a CW high-power proton Linac for accelerator-driven system
PRAMANA c Indian Academy of Sciences Vol. 78, No. 2 journal of February 2012 physics pp. 247 255 for accelerator-driven system RAJNI PANDE, SHWETA ROY, S V L S RAO, P SINGH and S KAILAS Physics Group,
More informationA Low Energy Beam Transport Design with high SCC for TAC Proton Accelerator
A Low Energy Beam Transport Design with high SCC for TAC Proton Accelerator * A. Caliskan 1, H. F. Kisoglu 2, S. Sultansoy 3,4, M. Yilmaz 5 1 Department of Engineering Physics, Gumushane University, Gumushane,
More informationSC forces RMS Envelope SC Child Langmuir law SCC Simulation Codes LEBT Simulation
Space Charge Effects N. Chauvin Commissariat à l Energie Atomique et aux Energies Alternatives, DSM/Irfu; F-91191 Gif-sur-Yvette, France. Nicolas.Chauvin@cea.fr May 31, 2012 N. Chauvin Space Charge Effects
More informationPreliminary design studies of a 100 MeV H /H + LINAC as injector for SNS synchrotron/ads LINAC
PRAMANA cfl Indian Academy of Sciences Vol. 59, No. 5 journal of November 2002 physics pp. 859 869 Preliminary design studies of a 100 MeV H /H + LINAC as injector for SNS synchrotron/ads LINAC S A PANDE,
More informationHALO SIMULATION IN A REALISTIC PROTON LINAC DESIGN
HALO SIMULATION IN A REALISTIC PROTON LINAC DESIGN M. Pabst and K. Bongardt, Forschungszentrum Jülich GmbH, 52425 Jülich, Germany A.P. Letchford, Rutherford Appleton Laboratory, Chilton, Didcot, UK Abstract
More informationBEAM INSTABILITIES IN LINEAR MACHINES 1.
BEAM INSTABILITIES IN LINEAR MACHINES 1 Massimo.Ferrario@LNF.INFN.IT CERN 4 November 2015 SELF FIELDS AND WAKE FIELDS The realm of collecdve effects Direct self fields Image self fields Space Charge Wake
More informationCOMPARISON OF TRACKING SIMULATION WITH EXPERIMENT ON THE GSI UNILAC
COMPARISON OF TRACKING SIMULATION WITH EXPERIMENT ON THE GSI UNILAC X.Yin 1,2, L. Groening 2, I. Hofmann 2, W. Bayer 2, W. Barth 2,S.Richter 2, S. Yaramishev 2, A. Franchi 3, A. Sauer 4 1 Institute of
More informationParticle-in-Cell Code BEAMPATH for Beam Dynamics Simulations in Linear Accelerators and Beamlines*
SLAC-PUB-1081 5 October 004 Particle-in-Cell Code BEAMPATH for Beam Dynamics Simulations in Linear Accelerators and Beamlines* Yuri K. Batygin Stanford Linear Accelerator Center, Stanford University, Stanford,
More informationSpace Charge Mi-ga-on
Space Charge Mi-ga-on Massimo.Ferrario@LNF.INFN.IT Hamburg June nd 016 OUTLINE The rms emicance concept rms envelope equa-on Space charge forces Space charge induced emicance oscilla-ons Matching condi-ons
More informationSTUDY ON SPACE CHARGE COMPENSATION OF LOW ENERGY HIGH INTENSITY ION BEAM IN PEKING UNIVERSITY*
STUDY ON SPACE CHARGE COMPENSATION OF LOW ENERGY HIGH INTENSITY ION BEAM IN PEKING UNIVERSITY* S. X. Peng 1,, A. L. Zhang 1, 2, H. T. Ren 1, T. Zhang 1, J. F. Zhang 1, Y. Xu 1, J. M. Wen 1, W. B. Wu 1,
More informationAnalysis of Slice Transverse Emittance Evolution in a Photocathode RF Gun. Abstract
SLAC PUB 868 October 7 Analysis of Slice Transverse Emittance Evolution in a Photocathode RF Gun Z. Huang, Y. Ding Stanford Linear Accelerator Center, Stanford, CA 9439 J. Qiang Lawrence Berkeley National
More informationLinac Beam Dynamics Code Benchmarking
Journal of Modern Physics, 2015, 6, 1044-1050 Published Online July 2015 in SciRes. http://www.scirp.org/journal/jmp http://dx.doi.org/10.4236/jmp.2015.68108 Linac Beam Dynamics Code Benchmarking Xuejun
More informationCERN Accelerator School. Intermediate Accelerator Physics Course Chios, Greece, September Low Emittance Rings
CERN Accelerator School Intermediate Accelerator Physics Course Chios, Greece, September 2011 Low Emittance Rings Part 1: Beam Dynamics with Synchrotron Radiation Andy Wolski The Cockcroft Institute, and
More informationSpace Charge in Linear Machines
Space Charge in Linear Machines Massimo.Ferrario@LNF.INFN.IT Egham September 6 th 017 Relativistic equation of motion dp dt = F p = γm o v γm o dv dt + m ov dγ dt = F β = v c dγ dt = d dt " a v % m o γ
More informationS2E: Solenoidal Focusing
S2E: Solenoidal Focusing Writing out explicitly the terms of this expansion: The field of an ideal magnetic solenoid is invariant under transverse rotations about it's axis of symmetry (z) can be expanded
More informationS2E: Solenoidal Focusing
S2E: Solenoidal Focusing The field of an ideal magnetic solenoid is invariant under transverse rotations about it's axis of symmetry (z) can be expanded in terms of the on axis field as as: solenoid.png
More informationSelf-consistent analysis of three dimensional uniformly charged ellipsoid with zero emittance
1 SLAC-PUB-883 17 May 1 Self-consistent analysis of three dimensional uniformly charged ellipsoid with zero emittance Yuri K. Batygin Stanford Linear Accelerator Center, Stanford University, Stanford,
More informationThree Loose Ends: Edge Focusing; Chromaticity; Beam Rigidity.
Linear Dynamics, Lecture 5 Three Loose Ends: Edge Focusing; Chromaticity; Beam Rigidity. Andy Wolski University of Liverpool, and the Cockcroft Institute, Daresbury, UK. November, 2012 What we Learned
More informationOptimization of radial matching section 1
Optimization of radial matching section Aleander Ovsannikov Dmitri Ovsannikov Saint-Petersburg State Universit (SPbSU) St. Petersburg Russia Sheng-Luen Chung National Taiwan Universit of Science and Technolog
More informationLinac JUAS lecture summary
Linac JUAS lecture summary Part1: Introduction to Linacs Linac is the acronym for Linear accelerator, a device where charged particles acquire energy moving on a linear path. There are more than 20 000
More informationBeam halo formation in high-intensity beams
Beam halo formation in high-intensity beams Alexei V. Fedotov,1,2 Brookhaven National Laboratory, Upton, NY 11973, USA Abstract Studies of beam halo became an unavoidable feature of high-intensity machines
More informationUPGRADE OF THE HIT INJECTOR LINAC-FRONTEND
UPGRADE OF THE HIT INJECTOR LINAC-FRONTEND S. Yaramyshev, W. Barth, M. Maier, A. Orzhekhovskaya, B. Schlitt, H. Vormann, GSI, Darmstadt R. Cee, A. Peters, HIT, Heidelberg Abstract The Therapy Linac in
More informationNotes on the HIE-ISOLDE HEBT
EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH HIE-ISOLDE-PROJECT-Note-13 Notes on the HIE-ISOLDE HEBT M.A. Fraser Abstract The HEBT will need to transfer the beam from the HIE-ISOLDE linac to up to four experimental
More informationSPACE CHARGE NEUTRALIZATION AND ITS DYNAMIC EFFECTS
SPACE CHARGE NEUTRALIZATION AND ITS DYNAMIC EFFECTS Ahmed Ben Ismail, Romuald Duperrier, Didier Uriot, CEN Saclay, 99 Gif sur Yvette cedex Nicolas Pichoff, CEN Bruyères le Châtel, 9680 Bruyères le Châtel,
More informationIFMIF High energy beam line design and beam expansion using non-linear multipole lenses and "step-like" magnet
IFMIF High energy beam line design and beam expansion using non-linear multipole lenses and "step-like" magnet N. Chauvin 1 R. Duperrier 1 P.A.P. Nghiem 1 J. Y. Tang 2 D. Uriot 1 Z.Yang 2 1 Commissariat
More informationA Compact Magnetic Focusing System for Electron Beams Suitable with Metamaterial Structures
A Compact Magnetic Focusing System for Electron Beams Suitable with Metamaterial Structures Ms. Kimberley Nichols University of New Mexico Advised by Dr. Edl Schamiloglu work performed in collaboration
More informationLund/Barnard USPAS Problem 1
Problem Lund/Barnard USPAS 207 Consider a round uniform ion beam with a current of ampere, composed of Hg + ions (atomic mass A =200), a kinetic energy of 2 MeV, a beam radius of 2 cm and normalized emittance
More informationDesign of a 10 MeV normal conducting CW proton linac based on equidistant multi-gap CH cavities *
Design of a 10 MeV normal conducting CW proton linac based on equidistant multi-gap CH cavities * LI Zhi-Hui( ) 1) The Key Labratory of Radiation Physics and Technology of Ministy of Eduction, Institute
More informationUPGRADE OF THE HIT INJECTOR LINAC-FRONTEND
UPGRADE OF THE HIT INJECTOR LINAC-FRONTEND S. Yaramyshev, W. Barth, M. Maier, A. Orzhekhovskaya, B. Schlitt, H. Vormann, GSI, Darmstadt R. Cee, A. Peters, HIT, Heidelberg Abstract The Therapy Linac in
More informationDEVELOPMENT OF LARGE SCALE OPTIMIZATION TOOLS FOR BEAM TRACKING CODES*
Proceedings of Hadron Beam 8, Nashville, Tennessee, USA DEVELOPMENT OF LARGE SCALE OPTIMIZATION TOOLS FOR BEAM TRACKING CODES* B. Mustapha # and P. N. Ostroumov Argonne National Laboratory, 97 S. Cass
More informationBEAM DYNAMICS ISSUES IN THE SNS LINAC
BEAM DYNAMICS ISSUES IN THE SNS LINAC A. Shishlo # on behalf of the SNS Accelerator Group, ORNL, Oak Ridge, TN 37831, U.S.A. Abstract A review of the Spallation Neutron Source (SNS) linac beam dynamics
More informationEFFICIENT 3D POISSON SOLVERS FOR SPACE-CHARGE SIMULATION
Abstract EFFICIENT 3D POISSON SOLVERS FOR SPACE-CHARGE SIMULATION Three-dimensional Poisson solver plays an important role in the self-consistent space-charge simulation. In this paper, we present several
More informationLecture 3: Modeling Accelerators Fringe fields and Insertion devices. X. Huang USPAS, January 2015 Hampton, Virginia
Lecture 3: Modeling Accelerators Fringe fields and Insertion devices X. Huang USPAS, January 05 Hampton, Virginia Fringe field effects Dipole Quadrupole Outline Modeling of insertion devices Radiation
More informationEvaluating the Emittance Increase Due to the RF Coupler Fields
Evaluating the Emittance Increase Due to the RF Coupler Fields David H. Dowell May 2014 Revised June 2014 Final Revision November 11, 2014 Abstract This technical note proposes a method for evaluating
More informationInvestigation of the Effect of Space Charge in the compact-energy Recovery Linac
Investigation of the Effect of Space Charge in the compact-energy Recovery Linac Ji-Gwang Hwang and Eun-San Kim, Kyungpook National University. 1370 Sankyok-dong, Buk-ku, Daegu, 702-701, Korea Tsukasa
More informationIntroduction to electron and photon beam physics. Zhirong Huang SLAC and Stanford University
Introduction to electron and photon beam physics Zhirong Huang SLAC and Stanford University August 03, 2015 Lecture Plan Electron beams (1.5 hrs) Photon or radiation beams (1 hr) References: 1. J. D. Jackson,
More informationLow Emittance Machines
CERN Accelerator School Advanced Accelerator Physics Course Trondheim, Norway, August 2013 Low Emittance Machines Part 1: Beam Dynamics with Synchrotron Radiation Andy Wolski The Cockcroft Institute, and
More information07. The Courant Snyder Invariant and the Betatron Formulation *
07. The Courant Snyder Invariant and the Betatron Formulation * Prof. Steven M. Lund Physics and Astronomy Department Facility for Rare Isotope Beams (FRIB) Michigan State University (MSU) US Particle
More informationBeam Dynamics with Space- Charge
Beam Dynamics with Space- Charge Chris Prior, ASTeC Intense Beams Group, RAL and Trinity College, Oxford 1 1. Linear Transverse Review of particle equations of motion in 2D without space-charge - Courant-Snyder
More informationPreliminary Simulation of Beam Extraction for the 28 GHz ECR Ion Source
Preliminary Simulation of Beam Extraction for the 28 GHz ECR Ion Source Bum-Sik Park*, Yonghwan Kim and Seokjin Choi RISP, Institute for Basic Science, Daejeon 305-811, Korea The 28 GHz ECR(Electron Cyclotron
More informationLongitudinal Beam Dynamics
Longitudinal Beam Dynamics Shahin Sanaye Hajari School of Particles and Accelerators, Institute For Research in Fundamental Science (IPM), Tehran, Iran IPM Linac workshop, Bahman 28-30, 1396 Contents 1.
More informationd2rx = F2 2 June 912 IEEE TRANSACTIONS ON NUCLEAR SCIENCE DESIGN OF A TRANSPORT SYSTEM FOR HIGH INTENSITY BEAMS
196 EEE. Personal use of this material is permitted. However, permission to reprint/republish this material for advertising or promotional purposes or for creating new collective works for resale or redistribution
More informationEnvelope and Matrix Codes (TRACE 3-D & TRANSPORT Introduction)
Envelope and Matrix Codes (TRACE 3-D & TRANSPORT Introduction) George H. Gillespie G. H. Gillespie Associates, Inc. P. O. Box 2961 Del Mar, California 92014, U.S.A. Presented at U. S. Particle Accelerator
More informationRF LINACS. Alessandra Lombardi BE/ ABP CERN
1 RF LINACS Alessandra Lombardi BE/ ABP CERN Contents PART 1 (yesterday) : Introduction : why?,what?, how?, when? Building bloc I (1/) : Radio Frequency cavity From an RF cavity to an accelerator PART
More informationOptimization of the electron-beam transport in the Israeli tandem FEL
Nuclear Instruments and Methods in Physics Research A 407 (1998) 350 355 Optimization of the electron-beam transport in the Israeli tandem FEL I. Merhasin, A. Abramovich *, Y. Pinhasi, A. Gover, J.S. Sokolowski
More informationIonenstrahlsimulationen außerhalb. Martin Droba
Ionenstrahlsimulationen außerhalb verfügbarer Standardprogramme Martin Droba Contents Motivation LORASR Magnetic codes TNSA LASIN Conclusion & Outlook Standard Tracking Programs TRACE3D, PARMILA, PARMTEQ,
More informationPIP-II Injector Test s Low Energy Beam Transport: Commissioning and Selected Measurements
PIP-II Injector Test s Low Energy Beam Transport: Commissioning and Selected Measurements A. Shemyakin 1, M. Alvarez 1, R. Andrews 1, J.-P. Carneiro 1, A. Chen 1, R. D Arcy 2, B. Hanna 1, L. Prost 1, V.
More informationAccelerator Physics Final Exam pts.
Accelerator Physics Final Exam - 170 pts. S. M. Lund and Y. Hao Graders: C. Richard and C. Y. Wong June 14, 2018 Problem 1 P052 Emittance Evolution 40 pts. a) 5 pts: Consider a coasting beam composed of
More informationTheoretical Task 3 (T-3) : Solutions 1 of 9
Theoretical Task 3 (T-3) : Solutions of 9 The Design of a Nuclear Reactor Uranium occurs in nature as UO with only 0.70% of the uranium atoms being 35 U. Neutron induced fission occurs readily in 35 U
More informationStudy of Analyzing and Matching of Mixed High Intensity Highly Charged Heavy Ion Beams
Study of Analyzing and Matching of Mixed High Intensity Highly Charged Heavy Ion Beams Youjin Yuan Institute of Modern Physics (IMP) Chinese Academy of Sciences 2016-7-6 HB2016, Malmö, Sweden New HIRFL
More informationHIGH CURRENT PROTON BEAM INVESTIGATIONS AT THE SILHI-LEBT AT CEA/SACLAY
TU31 Proceedings of LINAC 26, Knoxville, Tennessee USA HIGH CURRENT PROTON BEAM INVESTIGATIONS AT THE SILHI-LEBT AT CEA/SACLAY R. Hollinger, W. Barth, L. Dahl, M. Galonska, L. Groening, P. Spaedtke, GSI,
More informationLow Emittance Machines
Advanced Accelerator Physics Course RHUL, Egham, UK September 2017 Low Emittance Machines Part 1: Beam Dynamics with Synchrotron Radiation Andy Wolski The Cockcroft Institute, and the University of Liverpool,
More information4. Statistical description of particle beams
4. Statistical description of particle beams 4.1. Beam moments 4. Statistical description of particle beams 4.1. Beam moments In charged particle beam dynamics, we are commonly not particularly interested
More informationContents Motivation Particle In Cell Method Projects Plasma and Ion Beam Simulations
PIC Method for Numerical Simulation Ninad Joshi NNP Group 1 Contents Motivation Particle In Cell Method Projects Plasma and Ion Beam Simulations Motivation 3 Particle simulation Ion beams and Plasmas Accelerators
More informationIntro. Lecture 06: Initial Distributions
Intro. Lecture 06: Initial Distributions * Prof. Steven M. Lund Physics and Astronomy Department Facility for Rare Isotope Beams (FRIB) Michigan State University (MSU) US Particle Accelerator School (USPAS)
More informationUltra-Relativistic Heavy Ion Physics (FYSH551), May 31, 2013 Jan Rak and Thorsten Renk
Ultra-Relativistic Heavy Ion Physics (FYSH551), May 31, 2013 Jan Rak and Thorsten Renk Final Exam Instructions: Please write clearly. Do not just answer the questions, but document the thoughts leading
More informationBeam Dynamics and Electromagnetic Design Studies of 3 MeV RFQ for SNS Programme
J. Electromagnetic Analysis & Applications, 2010, 2, 519-528 doi:10.4236/jemaa.2010.29068 Published Online September 2010 (http://www.scirp.org/journal/jemaa) 519 Beam Dynamics and Electromagnetic Design
More informationMultiparameter optimization of an ERL. injector
Multiparameter optimization of an ERL injector R. Hajima a, R. Nagai a a Japan Atomic Energy Research Institute, Tokai-mura, Ibaraki 319 1195 Japan Abstract We present multiparameter optimization of an
More informationPhase Space Gymnastics
Phase Space Gymnastics As accelerator technology advances, the requirements on accelerator beam quality become increasingly demanding. Phase space gymnastics becomes a new focus of accelerator physics
More informationPhysics 663. Particle Physics Phenomenology. April 9, Physics 663, lecture 2 1
Physics 663 Particle Physics Phenomenology April 9, 2002 Physics 663, lecture 2 1 History Two Principles Electrostatic Cockcroft-Walton Accelerators Van de Graaff and tandem Van de Graaff Transformers
More information04.sup Equations of Motion and Applied Fields *
04.sup Equations of Motion and Applied Fields * Prof. Steven M. Lund Physics and Astronomy Department Facility for Rare Isotope Beams (FRIB) Michigan State University (MSU) S2: Transverse Particle Equations
More informationRFQ BEAM DYNAMICS DESIGN FOR LARGE SCIENCE FACILITIES AND ACCELERATOR-DRIVEN SYSTEMS
RFQ BEAM DYNAMICS DESIGN FOR LARGE SCIENCE FACILITIES AND ACCELERATOR-DRIVEN SYSTEMS Chuan Zhang # Institut für Angewandte Physik, Goethe-Universität, D-60438 Frankfurt a. M., Germany Abstract Serving
More information03.lec Solenoid Focusing*
03.lec Solenoid Focusing* Prof. Steven M. Lund Physics and Astronomy Department Facility for Rare Isotope Beams (FRIB) Michigan State University (MSU) PHY 905 Lectures Steven M. Lund and Yue Hao Michigan
More informationCompact Linac for Deuterons (Based on IH Structures with PMQ Focusing)
Compact Linac for Deuterons (Based on IH Structures with PMQ Focusing) S. Kurennoy, J. O Hara, L. Rybarcyk LANL, Los Alamos, NM Thanks to D. Barlow, F. Neri, and T. Wangler We are developing a compact
More informationPractical Lattice Design
Practical Lattice Design S. Alex Bogacz (JLab) and Dario Pellegrini (CERN) dario.pellegrini@cern.ch USPAS January, 15-19, 2018 1/48 D. Pellegrini - Practical Lattice Design Purpose of the Course Gain a
More informationCharged particle motion in external fields
Chapter 2 Charged particle motion in external fields A (fully ionized) plasma contains a very large number of particles. In general, their motion can only be studied statistically, taking appropriate averages.
More informationBuncher-System for FRANZ
Buncher-System for FRANZ Concept Beam-Dynamics Chopper-System Buncher-System for FRANZ 150 kv Terminal Wb = 120 kev Pb = 2.4 x 104 W Wb = 1 MeV Pb,max = 1 x 104 W Wb = 1.87-2.1 MeV Pb,max = 2.1 x 104 W
More informationComparison of simulated and observed beam profile broadening in the Proton Storage Ring and the role of space charge
PHYSICAL REVIEW SPECIAL TOPICS - ACCELERATORS AND BEAMS, VOLUME 3, 3421 (2) Comparison of simulated and observed beam profile broadening in the Proton Storage Ring and the role of space charge J. D. Galambos,
More informationTransverse emittance measurements on an S-band photocathode rf electron gun * Abstract
SLAC PUB 8963 LCLS-01-06 October 2001 Transverse emittance measurements on an S-band photocathode rf electron gun * J.F. Schmerge, P.R. Bolton, J.E. Clendenin, F.-J. Decker, D.H. Dowell, S.M. Gierman,
More informationTransverse Dynamics II
Transverse Dynamics II JAI Accelerator Physics Course Michaelmas Term 217 Dr. Suzie Sheehy Royal Society University Research Fellow University of Oxford Acknowledgements These lectures have been produced
More informationAccelerator Option for Neutron Sources & the
Accelerator Option for Neutron Sources & the Front end Injector of Proton Driver P.K. Nema P. Singh Bhabha Atomic Research Centre Mumbai, India US India workshop on ADS & thorium utilization Virginia Tech.
More informationLecture 2: Modeling Accelerators Calculation of lattice functions and parameters. X. Huang USPAS, January 2015 Hampton, Virginia
Lecture 2: Modeling Accelerators Calculation of lattice functions and parameters X. Huang USPAS, January 2015 Hampton, Virginia 1 Outline Closed orbit Transfer matrix, tunes, Optics functions Chromatic
More informationNovel Features of Computational EM and Particle-in-Cell Simulations. Shahid Ahmed. Illinois Institute of Technology
Novel Features of Computational EM and Particle-in-Cell Simulations Shahid Ahmed Illinois Institute of Technology Outline Part-I EM Structure Motivation Method Modes, Radiation Leakage and HOM Damper Conclusions
More informationComputations on Gabor lens having two different field distributions
IOSR Journal of Applied Physics (IOSR-JAP) e-issn: 2278-4861.Volume 6, Issue 6 Ver. II (Nov.-Dec. 2014), PP 06-11 Computations on Gabor lens having two different field distributions Saif KamilShnain Department
More informationAccelerator Physics Homework #3 P470 (Problems: 1-5)
Accelerator Physics Homework #3 P470 (Problems: -5). Particle motion in the presence of magnetic field errors is (Sect. II.2) y + K(s)y = B Bρ, where y stands for either x or z. Here B = B z for x motion,
More informationSimulation of transverse emittance measurements using the single slit method
Simulation of transverse emittance measurements using the single slit method Rudolf Höfler Vienna University of Technology DESY Zeuthen Summer Student Program 007 Abstract Emittance measurements using
More informationTransverse Beam Optics of the FLASH Facility
Transverse Beam Optics of the FLASH Facility ( current status and possible updates ) Nina Golubeva and Vladimir Balandin XFEL Beam Dynamics Group Meeting, 18 June 2007 Outline Different optics solutions
More informationMethod of Perturbative-PIC Simulation for CSR Effect
Method of Perturbative-PIC Simulation for CSR Effect Jack J. Shi Department of Physics & Astronomy, University of Kansas OUTLINE Why Do We Want to Do This? Perturbation Expansion of the Retardation of
More informationNumerical Methods II
Numerical Methods II Kevin Li With acknowledgements to: H. Bartosik, X. Buffat, L.R. Carver, S. Hegglin, G. Iadarola, L. Mether, E. Metral, N. Mounet, A. Oeftiger, A. Romano, G. Rumolo, B. Salvant, M.
More informationTransverse Equilibrium Distributions
* Transverse Equilibrium Distribution Functions: Outline Vlasov Model Vlasov Equilibria The KV Equilibrium Distribution Continuous Focusing Limit of the KV Equilibrium Distribution Equilibrium Distributions
More informationSynchrotron radiation
Synchrotron radiation When a particle with velocity v is deflected it emits radiation : the synchrotron radiation. Relativistic particles emits in a characteristic cone 1/g The emitted power is strongly
More informationTHEORETICAL COMPETITION. 1A. SPRING CYLINDER WITH MASSIVE PISTON (5 points)
Question 1 1A. SPRING CYLINDER WITH MASSIVE PISTON (5 points) Consider n=2 moles of ideal Helium gas at a pressure P 0, volume V 0 and temperature T 0 = 300 K placed in a vertical cylindrical container
More informationInstitute of Modern Physics, China Academy of Science, Lanzhou , China
Submitted to Chinese Physics C KONUS Beam Dynamics Design of Uranium IH-DTL for HIAF Dou Wei-Ping( 窦为平 ) a He Yuan( 何源 ) a Lu Yuan-Rong( 陆元荣 ) b, a Institute of Modern Physics, China Academy of Science,
More informationA NEW PARADIGM FOR MODELING, SIMULATIONS AND ANALYSIS OF INTENSE BEAMS
WEO2C6 Proceedings of HB21, Morschach, Switzerland A NEW PARADIGM FOR MODELING, SIMULATIONS AND ANALYSIS OF INTENSE BEAMS E. Nissen, B. Erdelyi, Department of Physics, Northern Illinois University, Dekalb,
More informationMeasurement of transverse Emittance
Measurement of transverse Emittance The emittance characterizes the whole beam quality, assuming linear behavior as described by second order differential equation. 1 It is defined within the phase space
More informationTransverse dynamics Selected topics. Erik Adli, University of Oslo, August 2016, v2.21
Transverse dynamics Selected topics Erik Adli, University of Oslo, August 2016, Erik.Adli@fys.uio.no, v2.21 Dispersion So far, we have studied particles with reference momentum p = p 0. A dipole field
More informationLecture 4: Emittance Compensation. J.B. Rosenzweig USPAS, UW-Madision 6/30/04
Lecture 4: Emittance Compensation J.B. Rosenzweig USPAS, UW-Madision 6/30/04 Emittance minimization in the RF photoinjector Thermal emittance limit Small transverse beam size Avoid metal cathodes? n,th
More informationChopping High-Intensity Ion Beams at FRANZ
Chopping High-Intensity Ion Beams at FRANZ C. Wiesner, M. Droba, O. Meusel, D. Noll, O. Payir, U. Ratzinger, P. Schneider IAP, Goethe-Universität Frankfurt am Main Outline 1) Introduction: The FRANZ facility
More informationEmittance Compensation. J.B. Rosenzweig ERL Workshop, Jefferson Lab 3/20/05
Emittance Compensation J.B. Rosenzweig ERL Workshop, Jefferson Lab 3//5 Emittance minimization in the RF photoinjector Thermal emittance limit Small transverse beam size Avoid metal cathodes? " n,th #
More informationS5: Linear Transverse Particle Equations of Motion without Space Charge, Acceleration, and Momentum Spread S5A: Hill's Equation
S5: Linear Transverse Particle Equations of Motion without Space Charge, Acceleration, and Momentum Spread S5A: Hill's Equation For a periodic lattice: Neglect: Space charge effects: Nonlinear applied
More informationS5: Linear Transverse Particle Equations of Motion without Space Charge, Acceleration, and Momentum Spread S5A: Hill's Equation
S5: Linear Transverse Particle Equations of Motion without Space Charge, Acceleration, and Momentum Spread S5A: Hill's Equation Neglect: Space charge effects: Nonlinear applied focusing and bends: Acceleration:
More informationSimulation of H ion source extraction systems for the Spallation Neutron Source with IBSimu
Simulation of H ion source extraction systems for the Spallation Neutron Source with IBSimu Taneli Kalvas 1, R. F. Welton 2, O. Tarvainen 1, B. X. Han 2 and M. P. Stockli 2 1. Department of Physics, University
More informationSimulations for photoinjectors C.Limborg
Simulations for photoinjectors C.Limborg 1- GTF Simulations Parmela modeling improvements Comparison to experimental results: 2ps & 4ps Sensitivity study Plans for future simulations 2- LCLS Injector Simulations
More informationTRIUMF Document: 10049
TRIUMF Document: 10049 DESIGN NOTE TRI-DN-07-18 Beam Line 2A Optics: Calculations and Measurements Legacy Document Document Type: Design Note Name: Signature: Date: Approved By: N/A Note: Before using
More informationSingle Particle Motion
Single Particle Motion C ontents Uniform E and B E = - guiding centers Definition of guiding center E gravitation Non Uniform B 'grad B' drift, B B Curvature drift Grad -B drift, B B invariance of µ. Magnetic
More informationSelected electrostatic problems (Lecture 14)
Selected electrostatic problems (Lecture 14) February 1, 2016 236/441 Lecture outline In this lecture, we demonstrate how to solve several electrostatic problems related to calculation of the fields generated
More informationSupporting Information
Supporting Information A: Calculation of radial distribution functions To get an effective propagator in one dimension, we first transform 1) into spherical coordinates: x a = ρ sin θ cos φ, y = ρ sin
More informationComparison of Different Simulation Codes with UNILAC Measurements for High Beam Currents
Comparison of Different Simulation Codes with UNILAC Measurements for High Beam Currents L. Groening, W. Barth, W. Bayer, G. Clemente, L. Dahl, P. Forck, P. Gerhard, I. Hofmann, M.S. Kaiser, M..Maier,
More informationBeam Dynamics and Emittance Growth
2010/03/09 Beam Dynamics and Emittance Growth Christoph Wiesner 1. Solenoidal Focusing and Emittance Growth 2. Beam Deflection and Emittance Growth 3. Time-dependent Kicker Fields, Electron Effects and
More information