R l - U N I V E R S I T Y M E S O N F A C I L I T YTR IUMFMAGNET SYSTEMS FOR THE TRANSPORT OF ION BEAMS WITH NO CHANGE OF PROPERTIESbyMa urice F. T a u t zUniversity of Victoria Physics DepartmentUNIVERSITY OF ALBERTA TRI-68-9SIMON FRASER UNIVERSITY UNIVERSITY OF VICTORIA UNIVERSITY OF BRITISH COLUMBIAMAGNET SYSTEMS FOR THE TRANSPORT OF ION BEAMS WITH NO CHANGE OF PROPERTIESMaurice F. Tautz TRI-68-9Physics Department University of Victoria October 1968TABLE OF CONTENTS1. Introduction . . . . . . . . . . . .2. Thick Lens Quadrupole Identity Sectionsa. Four Quadrupole Systems . . . . . . . .b. Minimum Gradient System . . . . . . . .c. Infinite Periodic System . . . . . . . .d. Recovery of Thin Lens Formulae as L Approaches Zero3. Acceptance Calculations for Thick Lens Identity SectionsAcknowledgments . . . . . . . . . . . .References . . . . . . . . . . . . .Appendix I.- Non-Existence of Doublet or Triplet Thin Lens Identity Sectionsa. Doublet . . . . . . . . . . . .b. Triplet . . . . . . . . . . . .Appendix II.- Thin Lens Identity Sectionsa. Four Thin Lens Quadrupole Systems . . . . .b. Maximum Focal Length System . . . . . . .c. Infinite Periodic System . . . . . . . .d. Five Thin Lens Quadrupole System . . . . .- 1 -1. IntroductionParticle motion through a beam transport system may be described, to first order, by using the well known matrix formalism"''.We review this method, only considering aspects relevant to this report.We choose an orthogonal Cartesian co-ordinate system where x denotes the horizontal displacement, y denotes the vertical displacement, and z denotes the displacement in the beam direction.Beam handling elements (we consider only drift spaces andquadrupole magnets) are represented by 2 x 2 transformation matrices which operate on particle displacement-divergence vectors. Thus, for motion in the xz-plane we haveThe zero subscripts denote input values, and the resultantvector on the left gives the particle displacement and slope at the exit of the element. A similar equation holds for motion in the yz-plane.In a drift space the particles travel in straight lines and the transformation matrix iswhere D is the length of the field-free region. This matrix is the same in both the xz- and yz-planes.thwhere is the i element transformation matrix- 2 -A quadrupole magnet consists of four magnetic pole pieces mounted on a common yoke as shown below.If the co-ordinate system is oriented as shown (z is out of page) then the magnetic field B arising from such a pole configuration can be described, to good approximation, byBx - gy ,Hy - gX ,and B = 0z3B 3Bxwhere g = is the field gradientwhich is assumed to be constant over the effective length L of the magnet and zero outside of it (hard edge model).- 3 -By substituting this field into the equations of motiondtwe get, to first order, the differential equations^ - k^y - 0 ,with k2 = f - = .003 | ,and where the particle momentum P is to be expressed in Gev/c and the units of g are Gauss/cm.equations that the quadrupole is focussing in the xz-plane and de- focussing (with the same strength) in the yz-plane. If the magnet polarities are reversed, these focussing effects are interchanged.The solutions to the two equations of motion can be written in matrix form asIt is apparent from the signs in the above two differentialThese transformation matrices and may be factored^ into the product of three matrices describing two drift spaces separated by a thin lens; iMA, KL'1 K itan —0/l K 1tanh—^KL0i-K sin KLK sinhKLKLK 1 tan f -y_l .KL K 1 tanh—In the "thin lens approximation", we take sin KL s. sinhKL = KL ,KL ~ - iKL KL tan — = tanhy- == —Then both and can be written in the formm. = / 1 ^ Y 1 "V1 L/2'1 \ o i /Vi/f i/ Vo iwhere f = 1/K^L is the thin lens focal length which ispositive or negative depending on whether the effect of the lens is focussing or defocussing.The transformation matrix T representing a transport system is found by computing the product of the individual element matrices; i.e.NT = n M. , i-1 1where N is the number of elements in the system.We note that this product matrix always has unit determinant since this is true for the matrices representing drift spaces and quadrupoles.It is of interest to design "identity sections" for which the system transformation matrix in both xz- and yz-planes satisfiesT = ± I , where I is the unit matrix.These systems change particle displacement-divergence vectors according to- 5 -dxand thus any region in the input (x , o ) phase space which is° dzinvariant under the operation of reflection through the origin is reproduced at the exit of an identity section. In particular, if the input beam is represented by a central ellipse in displacement- divergence phase space, then the output beam is represented by the same ellipse. An identity section is therefore useful for transporting such a beam effectively unchanged between tw<? points on the beam line.The initial investigations of these identity sections were carried out using the thin lens approximation. It is shown in Appendices I and II that there exist no thin lens doublet or triplet systems but in the case of four thin lens quadrupoles there exist systems with two degrees of freedom. The corresponding four quadrupole thick lens systems are discussed in Section 2 using assumptions based on the thin lens analysis.The systems found have the following characteristics:- the system transformation matrix is equal to minus the unit matrix- the systems have a symmetry property; that is, they consist of two identical symmetric doublets in succession- the total system length D is taken as a design parameter and can easily be adjusted- a system of a given length has three degrees of freedom. These can be taken as:r, the ratio of the last drift length” co the first£, the ratio of the central drift space for the system to Dd, the ratio of the quadrupole magnet effective length to D.An expression is obtained which relates the quadrupole field gradient for the system to I and d. This gradient is a minimum when the drift spacesbetween the magnets are of equal length.Identity sections with 4n (n is a positive integer) quadrupoles can be formed by using n of the above systems in succession. If the r-, Z - and d- values are the same for each identity section in the series, then when n becomes infinite, we obtain an infinite periodic system where any segment with length equal to a multiple of D is an identity section. The infinite and four quadrupole systems are discussed in Section 2. The phase space acceptances for a number of four quadrupole identity sections were calculated numerically. It was found that for a system of a given length the acceptance is a minimum when the system field gradient is a minimum; i.e., when the drift spaces between magnets are equal in length. It was also found that a maximum acceptance occurred for particles of momentum greater than the design momentum for the system. The details of these computations are given in Section 3.* measured from the center of the magnet- 6 -- 7 -2. Thick Lens Identity Sectionsa. Four Quadrupole SystemsIt is demonstrated in Appendices I and II that there exist no doublet or triplet thin lens quadrupole identity sections but for the case of four thin lenses there exist a double infinity of solutions. Therefore the thick lens identity sections were investigated by studying systems with four thick quadrupole lenses.We consider first a theorem taken from elementary matrix theory. This is the Hamilton-Cayley5 theorem which states that a matrix satisfies its characteristic equation. For a 2 x 2 matrix M the characteristic equation isdet = X2 -(trace M)A + det M = 0 , ^Mii—A M22 \ M 2l M22~^and hence we can write M2 - (trace M)M + (det M)I = 0 . 2.1If M represents a system transformation matrix, thendet M = 1 ,and from equation 2.1 M z = - I , 2.2when we impose the condition trace M = 0 . 2.3From equation 2.2 we conclude that if we can find a thick lens doublet satisfying equation 2.3, then a four quadrupole identity section can be formed by using two such doublets in succession.To this end we consider the thick lens doublet shown below:Here and g2 are the quadrupole field gradients and L the quadrupole effective length. The drift lengths D 1# D2 , D 3 are measured between magnet edges. From our thin lens calculations (Appendix II), we suspect that the magnets must be of equal strength and alternately focussing and defocussing. Accordingly we assume8l = “82 - 8 .The horizontal plane transformation matrix is then obtained from 1 ® 3 \ / costj) sincj>\ j 1 D2H =, , sinh<f>\ ,,coshcj; — v \ ( 1 Dj\K.0 1 J V -K sin<{> cos<|,/ \ 0 1 /\K sinh<f> cosh<))/ \ 0 1where, for brevity, <j> = KL ,Multiplying out the above matrices we find that the diagonal components of H are given byHi 1 = -D2D3K2 sin<j)sinh(J> + D2Kcosc(>sinhcj) + D3K(cos<t>sinh<f> - sin<J>cosh<f>)+ coscf>cosh<f> + sinc|>sinh<j)andH22 = -DiD2K2sin<J>sinh<t> + DiK(cos<}>sinh<f) - sin<}>coshcj))- D2Ksin<j)cosh<}) + coscj>cosh<f) - sin<J>sinh<f> .2.42.52. 6- 9 -Now substituting equations 2,5 and 2.6 into equation 2.3 givesD2(Di + D3)K2tanh(f)tan(() - K(D^ + D2 + D 3) (tanh<Ji - tantfi) - 2 = 0* 2.7Also, setting trace V = 0, where V is the vertical plane transformation matrix for this doublet, leads again to equation 2.7.Since the doublet is to represent one half of the identity section the expression for the total system length D isD = 2 ( Dj + D2 + D3 + 2 L) . 2.8‘It is convenient to take D and L to be fixed parameters. Then we have two equations (2.7 and 2.8) in the four unknowns Dj, D2 , D 3, K. We choose to solve for D2 and D 3 + D3 in terms of K. Then the complete identity section, shown below, will be specified, except for the ratio of the first drift length to the last, when the K-value is given.We solve for D2 by rearranging equation 2.8 D2 = | - 2L - (Di + D3) , and substituting this into equation 2.7 yields2.9G 1+d3 = (T -L) ±/D tn2 2coth<j>cot<j) V. ~ " 1+K(-^ - -L) (tanh^ - tan<())! 2 .10- 10 -It is convenient to write these equations in terms of dimensionless parameters. We defined E L/D ,0 E KD ,Dl1 = <L + _ ,£ i = d_2 ' Dd 2*2 = d + — ,d °3 £3 = 2 + 3 ’2.11d2 ( da e d + d = £i + 5,3 ,r = £3/^1 •From the last two defining equations we get£1 = 1 + r ’£ 2.121 + r ■1Substituting from equations 2.11 into 2.9 gives£2 = 1/2 - £ ,and equation 2.10 becomes2.13£ - ± + 4 -,1 jXo 2cothd0cotd0 . „,1 ,w ,- d)z - ------------- j^ l + 0 - d) (tanhd0 - tand0) 2.14Using this notation, the total system length is equal to unity. To scale to a system of length D we multiply all drift lengths by D and choose the system field gradient according to g = .003 .003 p2 ‘2.15We note that a four quadrupole identity section of a given length is completely determined when the three dimensionless ratios r, £, and d are specified.- 11 -A typical system is shown below:£-d-> r d— ^ -d—5<--- SL--- f f f 1 „ V V 11+r 2 £A j - t —7 * - >1+r'1In Fig. 1 we have plotted 0 as a function of £ for values of d equal to .01, .02, .03 and .04. As would be expected, we see that if we use longer magnets (larger d-values) the values of 0 become smaller, i.e. we can use weaker magnet gradients.b. Minimum-Gradient-SystemIt is apparent that each curve shown in Fig. 1 goes througha minimum at £ = 1/4. This means that for a system (with fixed D, d)the minimum gradient occurs when the drift spaces between the magnetsare equal in length since from equation 2 .12, £ = 1/4 implies also thatSL2 = 1/4. These minimum gradient (maximum focal length) systems werealso predicted by the thin lens analysis of Appendix II. For thesesystems, we can find a relationship between the total system lengthand the quadrupole constant K. We set the quantity under the squareroot sign in equation 2.10 equal to zero and obtain/■D T\? _ 2 coth<j>cot<j> (, . „,D , x , . D(4 “ w jl + K(-^ - - L) (tanh<f> - tan<())> }which can be put in the formD = ^ |L + (cotcf) - coth<J>) + /coth2<f> + cotc)>2 ) . 2.15- 12 -We can get a feeling for the strength of magnets required for typical identity sections by plotting K as a function of D. In Figure 2 is shown such a curve with the parameter L fixed at 0.4 meters. The dashed lines represent the gradients required for a beam with average particle momentum of 1.09 Gev/c (corresponding to 500.0 Mev protons).We see that for systems less than 5.5 meters gradients of greater than 1000 Gauss/cm. will be required.c. Infinite Periodic SystemWe can use the theory just developed to form a periodic system consisting of an infinite number of symmetric thick lens doublets as shown below.is 9 9 0 _ is f N2 & A/ 2 X, X/ A ■2 76-d-s *d-SHere, d and £ are the dimensionless parameters defined previously and the magnet field gradients are obtained from equation 2.15 where 0 satisfies equation 2.14. Any segment of length 1 unit is an identity section independent of where the start of the unit length section is taken (since the curves in Fig. 1 are symmetric about £ = 1/4). In fact any section of length m (m is a positive integer) has a transformation matrix T given byT = M2™ = ( - 1 )m I , 2.14where M denotes the matrix representing each doublet.- 13 -It is known that the condition for overall focussing in a periodic system of quadrupoles is| trace M | < 2 , 2.15where M is the matrix for one period.From the derivation given for the identity sections (recall equation 2.3) we note that the infinite system under discussion satisfies this criterion.*d . Recovery of Thin Lens Formulae as L Approaches ZeroAs L -> 0 the thick lens relationships established in this section reduce to those found in Appendix I using the thin lens approximation.Thus equation 2.9 becomesD2 as | - (Dj. + D 3) , 2.16and taking tanh q - tan <f> s KL ,in equation 2.10 we getD 1 + °3 = ~ ±d4 !)2 K^L2Squaring to remove the radical givesorD1 + d3 “ Jf = 7 1 K^L2 ’ d 2 ( d 2 ( d^ (D! + D 3).>} 2.17where f is the thin lens focal length.In addition, the second and fourth drift lengths of the identity section are equal, due to the assumption of symmetry (M2 = -I). Thus equations 2.16 and 2.17 are equivalent to the thin lens formulae found in Appendix II.* Evidently, from equation 2.14, the trace of the periodic system oscillates between -2 and +2 with a period of length 2 units.- 14 -3. Acceptance Calculations for Thick Lens Identity SectionsThe acceptance of a beam transport system is defined as that area in displacement-divergence phase space which contains the initial conditions of all particles that will pass through the system. The UVic IBM 360/44 computer was used to calculate the acceptance, assuming a four inch diameter beam tube, for the four quadrupole thick lens identity sections described in Section 2. The computations were done using the following parameter values:D = 7.0, 10.0, 13.0, 16.0, 19.0 meters ,L = 0.4 meters (all magnet effective lengths were equal) tr = 1.0 (first and last drift lengths were equal) .The results for the xz-plane are shown in Fig. 3.* We see that in each case the acceptance is a slowly varying function of the independent variable £, with a minimum at £ = 1/4. These curves are a bit deceptive in that the actual acceptance area shapes change radically with changing £-values. This effect is seen in Fig. 4. The acceptance areas for both horizontal and vertical planes for a 10.0 meter identity section with £ = .05, .25, and .45 are shown. The acceptance areas in the two planesare related to each other via a reflection through the ordinate orabscissa.* The results for the yz-plane were the same. This can be predicted theoretically by noting that, for systems symmetrical about the midpoint, a particle traversing the system in one plane experiences the same focussing as one travelling in reverse order in the other plane.- 15 -It can be seen from Fig. 4 that large acceptance areas with small divergences are obtained by using £-values close to 1/2. Acceptance shapes for identity sections of length other than 10.0 meters were investigated and similar results were obtained.From Figs. 3 and 4 it is seen that there is no single value of £ which should be chosen in all design situations. If the phase space area representing the input beam is small (much less than the acceptance with £ = 1/4) then the minimum gradient system is the most economical.When the input beam has a large phase space area (due to multiple scattering from a target, say), then one may have to increase or decrease the £-value in order to accommodate the beam using reasonable quadrupole apertures.For the 10.0 meter identity section with £ = .25 a plot*was also made of acceptance as a function of particle momentum. It turns out that the maximum acceptance does not occur at the design momentum (1.09 Gev/c) but at a higher value. This indicates that identity sections would probably not be useful in designing a system (e.g. a muon channel) where maximum acceptance is required.* See Figure 5.- 16 -ACKNOWLEDGMENTSThe author wishes to thank Dr. R.M. Pearce for many interesting discussions regarding this report and in particular for pointing out the relationship of this work to that of Wolfe. The author also wishes to thank Dr. D. Lobb and Dr. L.P. Robertson for their valuable comments during the progress of the work. The author also wishes to thank Mr. B. Gray for checking some of the algebra for the thin lens analysis.The acceptance calculations for Fig. 3 and Fig. 4 were done on a computer program written by B. Harrison under the direction of Dr. Lobb.The calculations for Fig. 5 were done on the computer program "IPSO FACTO".REFERENCESPenner, S. "Calculations of Properties of Magnetic Deflection Systems" Rev. Sci. Inst. 32, 2. (1961)Banford, A.P. "The Transport of Charged Particle Beams" p. 40.SPON. (1966)Wolfe, W.A., Private CommunicationCitron, A., Morpurgo and H. Overas "The High Intensity Muon Beam with Low Pion Contamination at the CERN Synchro-Cyclotron"CERN 63-65. (1963)Ayres, F. "Theory and Problems of Matrices" p. 81. (1962)Tautz, M.F. M.Sc. Thesis. (1968)- 18 -Appendix I.Non-Existence of Doublet or Triplet Thin Lens Identity Sectionsa. DoubletWe consider the thin lens quadrupole doublet shown below.Here, F3 and F2 are the focal lengths and Lx, L2 , L3 are the drift lengths for the system. The transformation matrix components for this system (for the xz-plane) are obtained fromT =1 L3 \ / 1 00 1 J I 1/F2 1f 1+L2/F1+L3/Fi+L3/F2 +L2L 3/(F1F2)1/F1+1/F2+L2/(F1F2)1 L2\ ( 1 01/Fi 1 o 1N ETSYTSFTSOSYRHOTSOSFRHYTSYSFRHYTN EN FRP ETN EN YN FRIP EP YB1+Lx/Fi+Li/F2+L2/F2+L1L2/(F!F2)1-1The transformation matrix components for the yz-plane are the same except thatF1 - F 1 ,Fo — Fc- 19 -If the system is to be an identity section, the off diagonal components of the matrix T must be zero. Applying this condition to component WYO in both planes gives1 / F i + 1/F2 + L2/(F1F2) = 0 , 1-2- 1 / F i - 1/F2 + L2/(FiF2) = 0 . 1-3Adding equations 1-2 and 1-3 yieldsL2/(F1F2) = 0 . 1-4For finite F^ and F2 this equation has only the solutionL2 = 0 , 1-5and putting this in equation 1-2 we getF1 = - f2 . 1-6Substituting equations 1-5 and 1-6 into 1-1 shows that the transformation matrix for both the xz- and yz-planes becomesand we conclude that an identity section occurs only in the trivial case thatL-x = L3 * 0 .- 20 -b. TripletWe consider the thin lens quadrupole triplet shown below.Here, Iq , L2 , L 3 , are the drift lengths and , Fg , F3 are the focal lengths for the system. To make the calculations easier we write the system transformation matrix asT = SM ,where M represents a doublet and S represents the last thin lens and drift space. For the xz-plane the components of M are given by equation 1-1 and for S we have1 L h \ I 1 0 \ / i + L4/F3 L40 1 / \i/f3 1 \ ' 1/Fg 1S =1-71-8The transformation matrices for the yz-plane are the same except thatFi ->■ - Fi »f2 - F2 -f3 "*■ ~ F3 •The condition for the system to be an identity section isT = SM = ± I ,or M =where we have useds . . I s Z .V*® . Z 3 s \det S = 1- 21 -In particular, we haveM 21 = + S2i ,which for the horizontal plane is equivalent to (from equations 1-8 and 1-1)1/F! + 1/F2 + L2/(F1F2) = + 1/F3 , and similarily for the vertical plane- 1/Fi - 1/F2 + L2/(F1F2) = ± I/F3 .Adding equations 1-10 and 1-11 leads toL2/(F!F2) = 0 , and for finite F^ and F2 we haveL2 = 0 .Since L2 is zero, the first and second thin lenses are adjacent and we are effectively back to the doublet case which has already been treated. Hence there are no thin lens triplet identity sections.1-101-111-12- 22 -Appendix II.Thin Lens Identity Sectionsa. Four Thin Lens Quadrupole SystemsWe consider the four thin lens quadrupole system shown below.Here Fj , F2 , F3 , F^ are the focal lengths and Iq , L2 , L 3 , Li, , L 5 , Lg are the drift lengths for the system. D is the total system length. We note that this system has been divided at some arbitrary point P between the two middle thin lenses so that we have effectively two doublets. This division does not restrict the system in any way and will simplify the subsequent calculations.If we label the transformation matrices for the two doublets by M and D then the condition for an identity section isT = RM = ± I ,F3 pi+P< Iq L2 ----- SK— L 3 — — L5 -— Lg -----------Dorwhere we have used det R = 1 •- 23 -In terms of matrix components we getM 11 “ ± R22 > II-lMl2 ~ + R12 , II-2M21 = + r21 , M22 = ± R11 » II-4II-3where the upper sign is for T = + I.Using the matrix components for a doublet (given in equation 1-1) we get for the xz-plane the following four equations:- from equation II-l1+L2 / F 1+L3/F 2+L3/ F 1+L2L 3 / ( F 1F2) = ± (l+Li+/F 3+Lit/Fit+L 5 /Fit+LitL 5 / (F 3Fi+) ) , I I -5- from equation II-2L1+L2+L3+L1L2/F1+L1L3/F2+L2L3/F2+L i L 3 /F i +L i L2L 3 / ( F i F+ L i 3/F2(332-3232(2-32(32 ( 32-2(3- from equation II-41+L 1 /F 1+ L 1 /F 2+L 1L2 / ( F 1F2 ) = ± (I+L 5/F 3+L6/F^+Lg/F3+L 5L6/ (F3F4)) . I I - 8 Similarly, for the xz-plane we have+LltL6/F3+LltL5L6/(F3Flf)) , II-6- from equation II-31/F1+1/F2+L2/(F1F2) = ^(1/F3+1/F1++L5/(F3F1+)) , II-71-L2/F1-L3/F2-L3/F1+L2L3/(F1F2) = ±(l-Llt/F3-Lit/FIt-L5/F4+LltL5/(F3Flt)) , II-9L1+L2+L3—LiL2/F3 —L^L3/F2 -L 2L3 / F 2 - L i L3 / F i +L i L2L3 / ( F i F+ L i F2(3232/222 /22 /22 - L ltL6 /F 3+L ltL 5L6 / ( F 3 F l+) ) ,11-10- 24 --1/F1-1/F2+L2/(F1F2) = +(-l/F3-l/Flt+L5/(F3F£t)) , 11-111-L1/F1-L1/F2-L2/F2+L1L2/(F1F2) = ±(1-L5/F3-L6/F4-L6/F3+L5L6/(F3F4)) . 11-12It is useful to consider the total system length as a parameter.We writeLi + L2 + L3 + L^ + L5 + Lg = D . 11-13Since det M = det R = 1 only six of the eightequations II-5 to 11-12 are independent. We have in total seven independent equations in the nine unknowns L 1} L2, L3 + , L5, Lg,Fi, F2, F3 and F^. In the above, the central drift space for thesystem (given by L 3 + L^) is counted as a single unknown since itsdivision into lengths L3 and Li+ was made only for convenience and hasno physical significance.To solve this set of equations, we manipulate as follows: equation II-7 + 11-11 givesF3Ftt/FiF2 = ^ Lg/L;? , 11-14equation II-7 - 11-11 gives1/F]+ 1/F2 = + (1/F3 + 1/Fit) , 11-15equation II-5 + II-9 givesL2L 3/(F1F2) = (-1 ± 1 ± L1+L5/(F3Fit)) . 11-16If we choose the upper sign we getF3F4/(FiF2) = L4L5/L2L 3 = (L5/L2) L4/L3 ,and from equation 11-14L^/Ls = - 1 ■- 25 -This result indicates that there are no feasible solutions since the middle drift space of the system, L3 + , is zero and we are back tothe three lens case. Hence there are no four lens identity sections with T = + I.We will therefore consider only the case where T = - I.Then in equation 11—16 we use the lower sign to getL2L3/(FiF2) = - 2 - L^L5/(F3Ftt) .From equation II-8 + 11-12L 1L2/(F1F2) = - 2 - L5L6/(F3F J , and from equation 11-18 - 11-17F3Fit/(F1F2) = (Llt - L6)/(Ll - L3) .Substituting equation 11-14 into 11-19 givesLi* + L3 = L6 + 1^ .We have also equation II-5 - II-9(L2 + L3)/F! + L3/F2 = - L^/Fg - (L^ + L5)/Fit , and equation II-8 - 11-12Li/Fi + (Li + L2)/F2 = - (l5 + L6)/F3 - Lg/Fj* .Then equation 11-21 + 11-22 combined with 11-13 yields(1/Fl + 1/F2) = (1/F3 + l / F O U - D/dj+L^Lg)) . Taking equation 11-23 - 11-15 gives0 = (1/F3 + 1/F1+)(D/(L1+L2+L3)) .We ignore the trivial solution D = 0 and get1/F3 + 1 ^ = 0 ,and it follows from the above equation and 11-15 thatF3 = - Fit »Fi = - F2 .11-1711-1811-1911-2011-2111-2211-2311-24- 26 -We define F : F[ - - F2 )11-25and G = F3 = - F4. •Now from equation 11-14 G/F = /L5/L2 . 11-26Substituting equation 11-25 in 11-21 we findG/F = L5/L2 ,and from equation 11-26 we havewith solutions^Ls/Lz = L5/L2 , L5/L2 = 0, 1 .We discard the zero solution as this leads us back to the three lens case.It is apparent from equations 11-25 that the four focal lengths for the system are equal in absolute value and alternate in sign.We can determine L2 , L5 in terms of D and L3 + L^.From equations 11-13 and 11-27 we obtainL2 + L 5 = 2 L2 = 2 L 5 = D - (Li + L 3 + L4 + L 6) , and using equation 11-20Substitution of equations 11-27, 11-28, 11-25 and 11-29 into equation 11-16 gives the system focal length also in terms of D and L3 + L^.We are left with L5 = L2 , 11-27and from equation 11-26 we get F = G . 11-28L2 = L5 = D/2 - (L3 + Li*) • 11-29It is 11-30- 27 -We have now found the general solution for an identity section of length D. If we define the dimensionless quantitiesi = (L3 + l4)/d ,£. = L./D , i = 1, 2, 5, 6r = Lg/Li ,f = F/D ,then we can summarize the results obtained so far as follows:- from equation 11-30 f = - ^ 'fj ^ > 0 < SL < %^2 = ^-5 = H ~ SL ,i l + s . I Z® ,- from equation 11-29- from equation 11-20with £l -SL2 =1 + r ’I0 < r1 + r"1 'Using this notation, the total system length is equal to unity and to scale to a system of length D we multiply all focal lengths and drift lengths by D. We note that an identity section of a given length is completely determined when the two dimensionless ratios SL and r are specified. The system is as shown below.+f -f +f -f^ ^ ^ ^ P 'S SLV 1+r 1 ^ " x, - ■ C '2 ~ ~ f K 1+r"1 ' >We see that this system has a fundamental symmetry. It can be regarded as two identical symmetric doublets of length one-half a unit; i.e., if the doublet matrix is denoted by M, then the identity section transformation matrix is M2 .11-3111-3211-3311-3411-35- 28 -b . Maximum Focal Length SystemWe can fix the degree of freedom specified by £ by requiring that the system focal lengths be a maximum.We have df_ _ d_d£ d£ < k - SL) j = 0which yields £ = 1/4 .From equations 11-32 to 11-35 the rest of the system is then specified b]£2 = £5 = 1/4 ,£i + £6 = 1/4 ,andf = 47T -11-36- 29 -c. Infinite Periodic SystemWe now apply the results of the previous section to the periodic system consisting of an infinite number of thin lens quadrupole doublets all of length one-half unit. This system is shown below.Here, f and £ are the dimensionless parameters defined in equation 11-31.If we restrict f and £ by the relation0 < SL < h 11-36then from equations 11-32 to 11-34 any section of unit length is an identity section independent of where the start of the unit section occurs. In fact, any section of length n units (n is a positive integer) has the transformation matrix (-l)nI.- 30 -d. Five Thin Lens Quadrupole SystemA particular section of this periodic system has been discussed3by Wolfe . In this case, the end points of a unit length section lie on the "centers" of two thin lenses and a five quadrupole identity section is formed by "splitting" each of these end lenses into an adjacent pair, both of focal length 2f. The system obtained is shown below.Wolfe considered the case of maximum focal length and obtained the valueA - 1/4 .c- 32I TOTAL SYSTEM LENGTH D(m.)20.0 io o(OOI / snjpDJ uio) 3 0 N V l d 3 0 0 V6.0I- 34 -XoXocviXocviXo*DISPLACEMENT (cm.)- 35 -( JUJ UIO) 3 0 N V l d 3 0 0 VDESIGN MOMENTUM P (Gev/c)
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- TRIUMF: Canada's national laboratory for particle and nuclear physics /
- Magnet systems for the transport of ion beams with...
Open Collections
TRIUMF: Canada's national laboratory for particle and nuclear physics
Magnet systems for the transport of ion beams with no change of properties Tautz, M. F. Oct 31, 1968
pdf
Page Metadata
Item Metadata
Title | Magnet systems for the transport of ion beams with no change of properties |
Alternate Title | TRIUMF brown reports TRI-68-9 |
Creator |
Tautz, M. F. |
Publisher | TRIUMF |
Date Issued | 1968-10 |
Subject |
Design study |
Genre |
Report |
Type |
Text |
Language | eng |
Date Available | 2015-05-14 |
Provider | Vancouver : University of British Columbia Library |
Rights | Attribution-NonCommercial-NoDerivs 2.5 Canada |
DOI | 10.14288/1.0107810 |
URI | http://hdl.handle.net/2429/53160 |
Affiliation |
TRIUMF Non UBC |
Peer Review Status | Unreviewed |
Scholarly Level | Researcher |
Copyright Holder | TRIUMF |
Rights URI | http://creativecommons.org/licenses/by-nc-nd/2.5/ca/ |
Aggregated Source Repository | DSpace |
Download
- Media
- 51833-TRI-1968-09.pdf [ 15.93MB ]
- Metadata
- JSON: 51833-1.0107810.json
- JSON-LD: 51833-1.0107810-ld.json
- RDF/XML (Pretty): 51833-1.0107810-rdf.xml
- RDF/JSON: 51833-1.0107810-rdf.json
- Turtle: 51833-1.0107810-turtle.txt
- N-Triples: 51833-1.0107810-rdf-ntriples.txt
- Original Record: 51833-1.0107810-source.json
- Full Text
- 51833-1.0107810-fulltext.txt
- Citation
- 51833-1.0107810.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.51833.1-0107810/manifest