SAST Lecture - III. Linear Accelerators. P N Prakash. IUAC, New Delhi
|
|
- Benedict Dennis
- 5 years ago
- Views:
Transcription
1 Lecture - III Linear Accelerators P N Prakash IUAC, New Delhi School on Accelerator Science and Technology (SAST ) Inter-University Accelerator Centre, New Delhi May 16-27,
2 Outline: Lecture - III Figures of Merit (continued) TEM Class Resonator Design P N Prakash - SAST-
3 we continue with Figures of Merit from where we left it off. P N Prakash - SAST-
4 Peak Electric Field In superconducting cavities the issue is not so much with electric breakdown or sparking as in DC accelerators, but rather it is with electron field emission. Field Emission (FE) arises at sufficiently high surface electric fields from micron or sub-micron size contaminant particles. When field emission begins, the cavity quality factor falls sharply. If the electron current from field emission becomes sufficiently intense, the heating caused by it can result in thermal breakdown (i.e. cavity going from Superconducting to normal state). So, it is important to reduce the Peak Surface Electric Field E Peak. In fact, the quantity to reduce in a cavity design is E Peak /E a. P N Prakash - SAST-
5 Peak Electric Field When field emission begins: quality factor Q falls sharply, heating caused by the current can result in thermal breakdown. So, it is important to reduce the Peak Surface Electric Field E Peak. In fact, The quantity to reduce in a cavity design is E Peak /E a. onset of field emission Q drop due to field emission Components which support high electric fields, such as the Drift Tube, must have sufficiently P N Prakash large - SAST- radii (no sharp corners).
6 Peak Electric Field E peak In order to properly compare different resonator geometries / models / designs, the peak electric field per unit accelerating gradient, E peak, is commonly used as a design parameter. E peak E peak _ total E a U 1J The peak surface electric field occurs at the open-end of the coaxial line, i.e. on the drift tube. The radius of curvature (or radii, if there are more than one) primarily decides (decide) the value of the peak electric field. P N Prakash - SAST-
7 Field Emission Fowler Nordheim model: The tunneling current density J in A/m 2 is given by: J E exp E where: is work function of the metal in ev E is the electric field in MV/m P N Prakash - SAST- The Fowler Nordheim theory has been experimentally verified for DC electric fields.
8 Field Emission For the tunneling current I (=JA e ) in RF cavities the formula has to be modified, and becomes: I e A E exp E whiskers or sharp projections enhance the electric field where β is the enhancement factor. Semi-ellipsoidal micro protrusion. Calculated geometric enhancement factor for the semi-ellipsoidal micro protrusion. P N Prakash - SAST-
9 Field Emission How to identify if a Cavity is suffering from FE? Field emission is always accompanied by X-Ray radiation!! onset of field emission onset of X-Ray radiation P N Prakash - SAST-
10 Field Emission Field emission is different from resonant electron multipacting, where it becomes impossible to increase the accelerating gradient even after pumping more power into the cavity, until the multipacting barrier is broken. Field emission is also different from high field Q- slope, where Q-drop is observed in the absence of any X-Ray radiation. Many field emitters in field emitting cavities can be processed through high power pulsed conditioning or through helium pulse conditioning. P N Prakash - SAST-
11 Peak Magnetic Field For superconducting niobium cavities the value of the critical magnetic field beyond which the material ceases to remain in the Meissner superconducting state, i.e. H C1, dictates the value of the peak surface magnetic field. Like the peak electric field, the quantity to be reduced in a cavity design is B Peak /E a. Recall that: TEM class resonators, which are usually designed around MHz frequency, typically operate at 4.2 K temperature (boiling point of LHe), so H C1 ~1350 G. The TM class elliptical cavities, designed for much higher frequencies, usually operate at 2 K, so H C1 ~1700 G. P N Prakash - SAST-
12 Peak Magnetic Field Here are some pointers: How to reduce B Peak? The Central Loading Element, also called Central Conductor, is usually tapered rather than keeping it cylindrical (it appears like the frustum of a cone). In elliptical cavities the cell is shaped at the equator to reduce the flux density. The diameter of the Outer Housing is increased in the high magnetic field region, which further reduces the value of the peak magnetic field. In addition, it also improves the Shunt Impedance and Geometry Factor of the structure. more on this later!!! P N Prakash - SAST-
13 Efficiency figures of merit P N Prakash - SAST-
14 Quality Factor Apart from the accelerating electric field (gradient) the Quality Factor is perhaps the most important figure of merit in terms of the operational performance of the resonator. The efficiency with which an accelerating structure stores (or retains) the energy is given by the figure of merit Quality Factor, usually denoted by Q (or Q 0 ). Energy stored Energy stored U Q=2 2 f Energy dissipated per cycle Power dissipated P where ω is the angular frequency, U is the stored energy in the cavity and P is the power dissipated in the cavity walls. P N Prakash - SAST-
15 Quality Factor The total energy in the cavity oscillates between the electric and magnetic fields (remember, the two standing waves have a phase difference of π/2). The time averaged stored energy in the electric field is equal to that in the magnetic field: U= 0 H dv= v 0 E dv 2 2 v The power loss due to Joule heating is defined through the power dissipated per unit area as: dp 1 R s H ds 2 2 P N Prakash - SAST- where R s is the surface resistance.
16 Quality Factor Assuming R s to be constant over the entire surface, we get: 1 2 P= Rs H ds 2 s The quality factor Q then becomes: U 0 H dv Q= v P 2 R H ds s s 2 or QR = 0 v s 2 s 2 H dv H ds P N Prakash - SAST-
17 Geometry Factor The Geometry Factor, G (=QRs), of the resonator is given by: Geometry Factor G: G= 0 v 2 s does not depend on surface resistance, so, it doesn t matter whether the cavity is in the normal state or in superconducting state. 2 H dv H ds Any figure of merit which does not depend on the surface resistance is going to be extremely useful for comparing cavities of different designs, or different models / versions of the same design. P N Prakash - SAST-
18 Shunt Impedance The Shunt Impedance is a figure of merit that measures how efficiently the cavity transfers its stored energy to the beam passing through it. The Shunt Impedance R sh is defined as: V 2 acc R sh= P The Shunt Impedance must be made as large as possible in order to maximize the accelerating voltage V acc for a given power P (or conversely minimize P for a given V acc ). P N Prakash - SAST-
19 Geometric Shunt Impedance The quality factor Q is given by: Q U P but U U E 2 0 a where U 0 is the stored energy at unit accelerating electric field and E a is the accelerating electric field (gradient). (the stored energy at unit accelerating electric field U 0 comes from bead pull measurement) but E a V L acc eff where L eff is the effective acceleration length. P N Prakash - SAST-
20 Geometric Shunt Impedance Using these equations, we get the Geometric Shunt Impedance R Q: sh R Q sh = 2 Leff U 0 The Geometric Shunt Impedance depends only on the geometric design of the resonator and is independent of the surface resistance of the material or its state (normal or superconducting). P N Prakash - SAST-
21 A word on Effective Length, L eff P N Prakash - SAST-
22 Effective Length L eff Often there is considerable confusion on what L eff really is!! L eff may not be the length over which the accelerating electric field extends. Beam Port Remember that the field penetrates into the beam ports and tapers off. But the value of the field inside the ports is generally not high. The geometry of the resonator along with the port size decides the extent of field penetration. L eff?? L eff?? An approximate thumb rule is that P N Prakash the - field SAST- falls 20 db every r BP.
23 Effective Length L eff Few important points that designers may remember are: - The definition of L eff is quite clear in the TM class cavities; it is the length of the cell. - In TEM class structures many definitions have been used by different groups. - Currently the most accepted definition of L eff in TEM 2 class resonators is 3. - When designing accelerating structures, however, it is important to retain the same definition of L eff in all the design models. Otherwise comparing them can become tricky, difficult and may be wrong. P N Prakash - SAST-
24 Effective Length L eff - You may, for example, choose L eff for the design calculations to be from the inside of the beam port at one end to the inside of the beam port at the other end. This can help in uniformly comparing the various design models. - Remember that L eff affects E a - This can in turn affect U 0 acc 2 1 E. - So, the definition of L eff must be consistent. a V L. eff P N Prakash - SAST-
25 Other important parameters P N Prakash - SAST-
26 Accelerator Cost Accelerators are expensive instruments to build. A good portion of the project cost and effort goes into civil construction. Acquiring land can also be difficult. It is always the endevour of accelerator designers to obtain maximum beam energy from the system for a given amount of power over a given length of the accelerator. Much of the R&D in accelerating structures like cavities, resonators, DTLs etc. is aimed at maximizing the energy gain over the shortest length. P N Prakash - SAST-
27 Real Estate Gradient It is important to realize that merely achieving large voltage gains alone from the accelerating structures is not sufficient. It is also important to consider over what physical length such energy gain is achieved by a group of accelerating structures (e.g. in a cryomodule), since it has direct bearing on the building size, and therefore the cost of the project. The Real Estate Gradient E REG is given by: E REG V L TOTAL ACC_LEN P N Prakash - SAST-
28 Real Estate Gradient What does E REG depend on? E REG V L TOTAL ACC_LEN To some extent E REG depends upon the resonator design (remember that the design can be optimized to produce large energy gains over short lengths). To a slightly larger extent E REG depends on the packing of the resonators. P N Prakash - SAST- pack em up!!!
29 Electron Multipacting The word Multipacting comes from the words multiple and impact. Electron Multipacting is a resonant process in which a large number of electrons build up in the resonator, absorbing all the power, thus making it impossible to increase the accelerating electric field in it by increasing the input power. The electron multipacting process starts when an electron present in the resonator gets accelerated by the RF field and impacts on the walls thereby producing secondary electrons. The secondary electrons in turn get accelerated by the fields and produce more electrons by again impacting on the cavity walls. P N Prakash - SAST-
30 Electron Multipacting If the number of secondary electrons produced per impact is more than the number of primary electrons, the electron current grows exponentially. The secondary electron emission coefficient (SEC), δ, of a material gives the number of secondary electrons produced per primary electron impacting the surface. Material δ max K max (ev) K 1 (ev) K 2 (ev) Copper Niobium The upper and lower crossovers for which δ>1 are given by K 1 and K 2. K max is the electron energy for which δ is maximum. SEC also depends on the state of the material. P N Prakash - SAST-
31 Electron Multipacting Although the secondary electron emission coefficient is primarily a function of the kinetic energy of the impacting electron, it also depends on the state of the material. R. Calder et al., CERN Secondary electron emission coefficient as a function of electron kinetic energy for niobium material P N following Prakash - SAST- different treatments.
32 Electron Multipacting Two kinds of electron multipacting are possible in RF cavities: One point multipacting Two point multipacting One point (or sometimes called single point) multipacting is predominant in TM class elliptical type cavities. Cyclotron frequency: H c g n n 0He m m n c e e 0 n (order) is an integer g e P N Prakash - SAST- K ee m 2 2 e 2 g
33 Electron Multipacting Two kinds of electron multipacting are possible in RF cavities: One point multipacting Two point multipacting Two point multipacting is predominant in TEM class resonators where the combination of RF frequency and gap length results in the involvement of two impact sites. The time period between the impacts must be half multiples of the RF time period T. (2n 1) T 2 where n is an integer. Drift Tube voltage V n : Kinetic energy of electron K n : V n 2 2 d me (2n 1) e K 2e V m d 2 2 n n 2 2 where d is the gap between the impacting sites. P N Prakash - SAST- K f d (2n 1) n where f is in MHz and d in cm
34 Tackling Multipacting In TM class cavities, electron multipacting is tackled by appropriately shaping those regions where it dominates (mostly near the equator). BEAM H E E Pill Box with Beam Ports Single Cell Cavity Use multipacting simulation codes for detailed study. Some of the codes available are: MultiPac, MUPAC, MULTIP, XING, TRACK-3D etc. P N Prakash - SAST-
35 Tackling Multipacting In TEM class resonators, electron multipacting is tackled by avoiding combinations of frequency and gap lengths which can result in electron energies close to δ max. For n=1: K f d Avoid cylindrical symmetry. Often the central loading element in TEM class resonator is tapered (albeit not to reduce multipacting), but it breaks the symmetry and reduces electron multipacting. If combinations of f and d which result in K being in the multipacting range cannot be avoided, then attempts must be made to reduce the area over which multipacting is likely to occur. P N Prakash - SAST-
36 Operating Temperature Critical Temperature T C of Niobium Material is ~9.2K. Resonators made out of niobium material must be cooled below 9.2K to make them superconducting. The boiling point of liquid helium (LHe) is 4.2K, which is ideally suited for cooling down niobium resonators. Helium liquefiers are commercially available for producing liquid helium. The surface resistance, R S, of superconducting niobium is: Rs RBCS Rres where R BCS is the BCS resistance and R res is the residual resistance. P N Prakash - SAST-
37 Operating Temperature T For T<T C /2, a convenient expression for R BCS (nω) is: R BCS f T 4 2 e T (where f in GHz & T in K) Frequency f R BCS (T=4.2 K) R BCS (T=2K) 100 MHz 3.2 n 0.06 n 300 MHz 28.4 n 0.6 n 1300 MHz 533 n 11 n TEM class resonators are usually designed around MHz frequency (essentially to keep its size neither too small nor too big). At this frequency R BCS 3 nω at 4K. P N Prakash - SAST-
38 Operating Temperature T A reasonable value for the Geometry Factor G (=QR S ) for a TEM class resonator ~15. With R BCS 3 nω at 4.2K, the limiting quality factor will be Q (=G/R S ) , which is quite high. For Q (factor of 5 lower than what we have estimated above), the power, P, required to operate the resonator at E a =4 MV/m, given that f=100 MHz, U 0 =0.1 J/(MV/m) 2, will be ~1 W!! Note that I have not taken into account the Residual Resistance, R res. So, resonators designed around MHz can be very conveniently operated at 4.2K temperature (boiling point P N Prakash of - SAST- liquid helium).
39 Operating Temperature T Frequency f R BCS (T=4.2 K) R BCS (T=2K) 100 MHz 3.2 n 0.06 n 300 MHz 28.4 n 0.6 n 1300 MHz 533 n 11 n TM class elliptical cavities are generally designed to operate at ~ MHz frequencies. At these frequencies R BCS > 100 nω at 4.2K. Clearly R BCS has to be reduced to a more tolerable value. At 2K, R BCS ~5-10 nω; which is tolerable. TM class Elliptical Cavities are operated at 2K. P N Prakash - SAST-
40 Residual Resistance R res The residual resistance R res arises from inclusion of foreign materials, impurities, residues from chemical processing, formation of oxides and hydrides on the niobium surface and trapped magnetic flux during cool down. By taking precautions during resonator fabrication, processing and in handling the material, a carefully prepared niobium surface can reach R res <5 nω, which is acceptable. P N Prakash - SAST-
41 Some more Parameters There are some more parameters which are very important. These are: Frequency Tunability Frequency Sensitivity Structural Stability Mechanical Frequency (microphonics). Due to shortage of time, I can t discuss them. P N Prakash - SAST-
42 What are the Key Parameters for designing a Niobium Resonator? P N Prakash - SAST-
43 Key Parameters Key Parameters for designing a Superconducting Niobium Resonator: 1. Synchronous Velocity β 0 2. Resonator Frequency f 3. Energy Gain E Gain 4. Transit Time Factor TTF 5. Stored Energy U 0 6. Geometry Factor QR S 7. Peak Magnetic Field B peak 8. Peak Electric Field E peak 9. Geometric Shunt Impedance R sh /Q 10. Operating Temperature T 11. Real Estate Gradient 12. Electron Multipacting 13. Frequency Tunability Δf ST 14. Frequency Sensitivity 15. Mechanical Stability Remember, out of these, some parameters have to be selected while others need to be optimized. P N Prakash - SAST-
44 Key Parameters It is simply impossible to optimize all the key parameters in a resonator design since they are coupled to each other. Remember the Car example that I had given in the second lecture. It is impossible to optimize all the parameters simultaneously. Therefore only a Global Optimization in a resonator design is possible to optimize its parameters. P N Prakash - SAST-
45 Let us look at a couple of Key Parameters in detail for designing TEM Class Resonators. Due to shortage of time, I won t be able to go over all of them. P N Prakash - SAST-
46 Resonator Frequency f There are two competing requirements for frequency selection. Lower f 1. Fewer Resonators are needed to reach the final energy goal. For example, in a 2-gap QWR the center of the first gap to center of the second gap d is given by: d 2 Leff d g where g is the gap between the drift tube and the end cap. P N Prakash - SAST- Higher f 1. Length of the Loading element (also called the Central Conductor), L λ/2 or λ/4, reduces. This is a major advantage, especially in the Top Loading type cryomodule where the height of the beam hall (more so in existing beam halls) is a major constraint.
47 Mechanical Frequency Lowest mechanical vibration frequency of a hollow rod anchored at one end (approximate formula neglecting the inner radius r i ): where: mech ro Y L 2 L is the length of the rod, r i & r o are the inner and outer radii of the rod, Y is the Young s Modulus of the material ρ is the density of the material Response to a particular frequency of disturbance falls as So, it is logical to increase ν mech as high as possible. For the Tapered Conductor, to a good approx.: R R 2 S O r o, So, in addition to reducing the peak magnetic field B peak (as we shall see), tapering also increases ν mech r o tapered r o untapered r o R O 1 2 mech Note: The mechanical frequency also depends on the mass at the open-end, which should be kept small. P N Prakash - SAST-
48 Resonator Frequency f There are two competing requirements for frequency selection. Lower f Higher f 2. Energy & Time spreads of the beam are better controlled. E t q Eaf sin s A E A E t Phase Acceptance: =5 (±2.5º about s ); at low frequencies this is more. P N Prakash - SAST- 2. Increases ν mech : mech ro Y L 2 this reduces microphonic coupling, and hence the RF frequency jitter Δω shake 2 PC shakeu shakeu 0Ea
49 Summary - Resonator Frequency There is no sharp cutoff where one of the (conflicting) requirements begins to dominate the other. In the frequency range MHz the conflicts are largely reconciled and reasonably taken care of. The final choice of frequency is made keeping β 0 in mind to obtain reasonably long effective acceleration length (L eff ) while keeping other constraints in mind. TEM class resonators are usually designed around MHz frequency. P N Prakash - SAST-
50 Peak Magnetic Field B peak In order to properly compare different resonator geometries / models / designs, the peak magnetic field per unit accelerating gradient, B peak, is commonly used as a design parameter. B peak B E total a U 1J Since the peak magnetic field produced as a result of the applied RF must remain below H C1 of superconducting niobium, one of the main efforts in the resonator design is directed towards reducing B peak, so as to maximize the achievable accelerating gradient. P N Prakash - SAST-
51 Reducing B peak One way to reduce B peak in the resonator is by tapering the central conductor (loading element). Consider a quarter wave coaxial line (this is the simplest TEM class structure). Note that for simplicity, we are considering the coaxial line without a drift tube. B peak 0V0 2 R Z open 0 P N Prakash - SAST- where Z 0 is the characteristic impedance of the coaxial line and V 0 is the peak voltage at the openend of the coaxial line.
52 Reducing B peak Consider now a quarter wave coaxial line with a linearly tapering central conductor (loading element). The tapered central conductor may be approximated to a stepped coaxial line of two sections of radii R short & R open, each of length l/2. P N Prakash - SAST-
53 Reducing B peak B For such a stepped (tapered) coaxial line, it can be shown that: 0 0 peak 1 2 Rshort sin(2 ) 2 Z1 Z2 for / 4, It can be shown that: for R short >R open : V 1 sin(2 ) 1 ½ R short [Z 1 +Z 2 ]>R open Z 0 B P N Prakash - SAST- peak where Z 1 & Z 2 are the characteristic impedances of the two sections of the coaxial line, and V 0 is the peak voltage at the open-end. Therefore for the same V 0 : stepped B peak untapered
54 Reducing B peak Another way to reduce B peak is by increasing the diameter of the outer cylinder (housing). Characteristic Impedance Z 0 of a coaxial line is given by: Z 0 1 R 0 ln 2 R hsg 0 open where R CC is the radius of the loading element / central conductor, (which can be R short or R open in our case). For the straight (un-tapered) line: B peak 0V0 2 R Z open 0 Z 0, B peak Physically, there is more volume P N Prakash - for SAST- the magnetic flux to expand.
55 Reducing B peak New shapes of tapered outer housing further reduce B peak. This in turn increases the achievable accelerating gradient E a. P N Prakash - SAST- Peter Ostroumov & Michael Kelly Argonne National Lab
56 What I could not cover? Perhaps more than what I covered in the lectures!!! TM class cavity design Complete design aspects of TEM class resonators Simulation codes Issues connected with niobium material Niobium cavity fabrication Processing of superconducting niobium cavities Bead pull measurements Cavity testing (dressing, loading, cooldown, MP conditioning, high power & helium pulse conditioning etc.) Issues related to electronic control of cavities (phase & amplitude locking, microphonics, etc.) Operational issues (in actual Linacs) New developments and R&D in this field R&D on new materials Last, but not the least, P the N Prakash work - SAST- done at IUAC.
57 Thank you P N Prakash - SAST-
Superconducting RF Accelerators: Why all the interest?
Superconducting RF Accelerators: Why all the interest? William A. Barletta Director, United States Particle Accelerator School Dept. of Physics, MIT The HEP prespective ILC PROJECT X Why do we need RF
More informationSpoke and other TEM-class superconducting cavities. J.L. Muñoz, ESS-Bilbao Academy-Industry Matching Event CIEMAT, Madrid, 27.May.
Spoke and other TEM-class superconducting cavities J.L. Muñoz, ESS-Bilbao Academy-Industry Matching Event CIEMAT, Madrid, 27.May.2013 Outline Introduction Basic design of TEM cavities Cavity design issues
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 informationFrom the Wideröe gap to the linac cell
Module 3 Coupled resonator chains Stability and stabilization Acceleration in periodic structures Special accelerating structures Superconducting linac structures From the Wideröe gap to the linac cell
More informationSimulation of RF Cavity Dark Current in Presence of Helical Magnetic Field
Preprint FERMILAB-TM-2467-TP. Simulation of RF Cavity Dark Current in Presence of Helical Magnetic Field GennadyRomanov, Vladimir Kashikhin Abstract. In order to produce muon beam of high enough quality
More informationCERN Accelerator School. RF Cavities. Erk Jensen CERN BE-RF
CERN Accelerator School RF Cavities Erk Jensen CERN BE-RF CERN Accelerator School, Varna 010 - "Introduction to Accelerator Physics" What is a cavity? 3-Sept-010 CAS Varna/Bulgaria 010- RF Cavities Lorentz
More informationDe Prisco, Renato; Comunian, Michele; Grespan, Franscesco; Piscent, Andrea; Karlsson, Anders
ESS DTL RF MODELIZATION: FIELD TUNING AND STABILIZATION De Prisco, Renato; Comunian, Michele; Grespan, Franscesco; Piscent, Andrea; Karlsson, Anders 2013 Link to publication Citation for published version
More informationReview of proposals of ERL injector cryomodules. S. Belomestnykh
Review of proposals of ERL injector cryomodules S. Belomestnykh ERL 2005 JLab, March 22, 2005 Introduction In this presentation we will review injector cryomodule designs either already existing or under
More informationDesign of an RF Photo-Gun (PHIN)
Design of an RF Photo-Gun (PHIN) R. Roux 1, G. Bienvenu 1, C. Prevost 1, B. Mercier 1 1) CNRS-IN2P3-LAL, Orsay, France Abstract In this note we show the results of the RF simulations performed with a 2-D
More informationLecture 1 Introduction to RF for Accelerators. Dr G Burt Lancaster University Engineering
Lecture 1 Introduction to RF for Accelerators Dr G Burt Lancaster University Engineering Electrostatic Acceleration + - - - - - - - + + + + + + Van-de Graaff - 1930s A standard electrostatic accelerator
More informationLecture 5: Photoinjector Technology. J. Rosenzweig UCLA Dept. of Physics & Astronomy USPAS, 7/1/04
Lecture 5: Photoinjector Technology J. Rosenzweig UCLA Dept. of Physics & Astronomy USPAS, 7/1/04 Technologies Magnetostatic devices Computational modeling Map generation RF cavities 2 cell devices Multicell
More informationFundamentals of RF Cavities
USPAS June, 2016 Fundamentals of RF Cavities Suba De Silva, S. A. Bogacz, G. A. Krafft, and R. Gamage Old Dominion University / Jefferson Lab Colorado State University Lecture 10 Outline RF Cavities Cavity
More informationMicrowave multipactor breakdown in open two-wire transmission systems
Microwave multipactor breakdown in open two-wire transmission systems Joel Rasch*, D. Anderson*, V. E. Semenov, E. Rakova, J. F. Johansson, M. Lisak* and J. Puech *Chalmers University of Technology Outline
More informationRadiation field inside the tunnel of the Linear Collider TESLA
Laboratory Note DESY D3 113 April 2000 Radiation field inside the tunnel of the Linear Collider TESLA Dark current, first attempt A. Leuschner, S. Simrock Deutsches Elektronen-Synchrotron DESY Abstract
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 informationRF cavities (Lecture 25)
RF cavities (Lecture 25 February 2, 2016 319/441 Lecture outline A good conductor has a property to guide and trap electromagnetic field in a confined region. In this lecture we will consider an example
More informationFundamental Concepts of Particle Accelerators III : High-Energy Beam Dynamics (2) Koji TAKATA KEK. Accelerator Course, Sokendai. Second Term, JFY2012
.... Fundamental Concepts of Particle Accelerators III : High-Energy Beam Dynamics (2) Koji TAKATA KEK koji.takata@kek.jp http://research.kek.jp/people/takata/home.html Accelerator Course, Sokendai Second
More informationPH2200 Practice Final Exam Summer 2003
INSTRUCTIONS 1. Write your name and student identification number on the answer sheet. 2. Please cover your answer sheet at all times. 3. This is a closed book exam. You may use the PH2200 formula sheet
More informationLongitudinal Dynamics
Longitudinal Dynamics F = e (E + v x B) CAS Bruges 16-25 June 2009 Beam Dynamics D. Brandt 1 Acceleration The accelerator has to provide kinetic energy to the charged particles, i.e. increase the momentum
More informationMonday July 14. Capacitance demo slide 19 Capacitors in series and parallel slide 33 Elmo example
Monday July 14 Lecture 5 Capacitance demo slide 19 Capacitors in series and parallel slide 33 Elmo example Lecture 6 Currents and esistance Lecture 9 Circuits Wear Microphone 1 3 Lecture 6 Current and
More informationEstimates of local heating due to trapped modes in vacuum chamber
Estimates of local heating due to trapped modes in vacuum chamber Gennady Stupakov SLAC National Accelerator Laboratory, Menlo Park, CA 94025 CERN, April 29, 2016 2 Motivation The motivation for this analysis
More informationChemistry Instrumental Analysis Lecture 35. Chem 4631
Chemistry 4631 Instrumental Analysis Lecture 35 Principle components: Inlet Ion source Mass analyzer Ion transducer Pumps Signal processor Mass analyzers Quadrupole Time of Flight Double Focusing Ion
More informationDesign of a Double Spoke Cavity for C-ADS
Design of a Double Spoke Cavity for C-ADS JIANG Tian-Cai ( 蒋天才 ) 1,2 Huang Yu-Lu ( 黄玉璐 ) 1,2 HE Yuan( 何源 ) 1 LU Xiang-Yang ( 鲁向阳 ) 1,3 Zhang Sheng-Hu ( 张生虎 ) 1 1 Institute of Modern Physics, Chinese Academy
More informationPhysics 610. Adv Particle Physics. April 7, 2014
Physics 610 Adv Particle Physics April 7, 2014 Accelerators History Two Principles Electrostatic Cockcroft-Walton Van de Graaff and tandem Van de Graaff Transformers Cyclotron Betatron Linear Induction
More informationElliptical Cavities: Proven SRF Option
Elliptical Cavities: Proven SRF Option Terry L. Grimm Michigan State University July 2005 Outline Multi-cell elliptical cavities Demonstrated for β > 0.4 Future trends Multi-spoke cavities Limited experimental
More informationElectricity & Magnetism Study Questions for the Spring 2018 Department Exam December 4, 2017
Electricity & Magnetism Study Questions for the Spring 2018 Department Exam December 4, 2017 1. a. Find the capacitance of a spherical capacitor with inner radius l i and outer radius l 0 filled with dielectric
More informationTHz Electron Gun Development. Emilio Nanni 3/30/2016
THz Electron Gun Development Emilio Nanni 3/30/2016 Outline Motivation Experimental Demonstration of THz Acceleration THz Generation Accelerating Structure and Results Moving Forward Parametric THz Amplifiers
More informationSemiconductor Physics and Devices
Introduction to Quantum Mechanics In order to understand the current-voltage characteristics, we need some knowledge of electron behavior in semiconductor when the electron is subjected to various potential
More informationand another with a peak frequency ω 2
Physics Qualifying Examination Part I 7-Minute Questions September 13, 2014 1. A sealed container is divided into two volumes by a moveable piston. There are N A molecules on one side and N B molecules
More informationOptimization of the End Cells in SRF Cavities
Optimization of the End Cells in SRF Cavities Jonathan Luk Universtity of California, San Diego, La Jolla, CA 92093 (Dated: August 12, 2005) The shape of the end cells of the RF superconducting cavities
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 informationIIT JEE PAPER The question paper consists of 3 parts (Chemistry, Mathematics and Physics). Each part has 4 sections.
IIT JEE - 2009 PAPER 2 A. Question paper format: 1. The question paper consists of 3 parts (Chemistry, Mathematics and Physics). Each part has 4 sections. 2. Section I contains 4 multiple choice questions.
More informationHigh Gradient RF Studies
High Gradient RF Studies J. Norem Argonne Jefferson Lab Aug. 4, 05 Collaborators Experiments in Fermilab Muon Test Area (MTA) J. Norem, Argonne A. Moretti, A. Bross, Z. Qian FNAL Y. Torun, IIT D. Li, M.
More information(a) zero. B 2 l 2. (c) (b)
1. Two identical co-axial circular loops carry equal currents circulating in the same direction: (a) The current in each coil decrease as the coils approach each other. (b) The current in each coil increase
More informationHomework 2: Forces on Charged Particles
Homework 2: Forces on Charged Particles 1. In the arrangement shown below, 2 C of positive charge is moved from plate S, which is at a potential of 250 V, to plate T, which is at a potential of 750 V.
More informationBasic Concepts of Measurements Made on Superconducting RF Cavities
Basic Concepts of Measurements Made on Superconducting RF Cavities Jens Knobloch, Laboratory of Nuclear Studies, Cornell University, August 1991 Presented here is an introduction to some of the concepts
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 informationDownloaded from
Question 4.1: A circular coil of wire consisting of 100 turns, each of radius 8.0 cm carries a current of 0.40 A. What is the magnitude of the magnetic field B at the centre of the coil? Number of turns
More informationANALYSIS OF HIGH ORDER MODES IN 1.3 GHZ CW SRF ELECTRON LINAC FOR A LIGHT SOURCE
ANALYSIS OF HIGH ORDER MODES IN 1.3 GHZ CW SRF ELECTRON LINAC FOR A LIGHT SOURCE A. Sukhanov, A. Vostrikov, V. Yakovlev, Fermilab, Batavia, IL 60510, USA Abstract Design of a Light Source (LS) based on
More informationParticle physics experiments
Particle physics experiments Particle physics experiments: collide particles to produce new particles reveal their internal structure and laws of their interactions by observing regularities, measuring
More informationJoint Entrance Examination for Postgraduate Courses in Physics EUF
Joint Entrance Examination for Postgraduate Courses in Physics EUF First Semester/01 Part 1 4 Oct 011 Instructions: DO NOT WRITE YOUR NAME ON THE TEST. It should be identified only by your candidate number
More information$)ODW%HDP(OHFWURQ6RXUFHIRU/LQHDU&ROOLGHUV
$)ODW%HDP(OHFWURQ6RXUFHIRU/LQHDU&ROOLGHUV R. Brinkmann, Ya. Derbenev and K. Flöttmann, DESY April 1999 $EVWUDFW We discuss the possibility of generating a low-emittance flat (ε y
More informationDemonstration Some simple theoretical models Materials How to make superconductors Some applications
Superconductivity Demonstration Some simple theoretical models Materials How to make superconductors Some applications How do we show superconductivity? Superconductors 1. have an electrical resistivity
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 informationSection 4 : Accelerators
Section 4 : Accelerators In addition to their critical role in the evolution of nuclear science, nuclear particle accelerators have become an essential tool in both industry and medicine. Table 4.1 summarizes
More informationHigh-Gradient, Millimeter Wave Accelerating Structure
DPF2015-337 August 31, 2015 High-Gradient, Millimeter Wave Accelerating Structure S.V. Kuzikov 1,2, A.A. Vikharev 1,3, N.Yu. Peskov 1 1 Institute of Applied Physics, 46 Ulyanova Str., Nizhny Novgorod,
More informationStatus of linear collider designs:
Status of linear collider designs: Main linacs Design overview, principal open issues G. Dugan March 11, 2002 Linear colliders: main linacs The main linac is the heart of the linear collider TESLA, NLC/JLC,
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 informationGood Luck! Mlanie LaRoche-Boisvert - Electromagnetism Electromagnetism and Optics - Winter PH. Electromagnetism and Optics - Winter PH
1 Notes: 1. To submit a problem, just click the Submit button under it. The Submit All button is not necessary. 2. A problem accepted as correct by CAPA will be highlighted in green. Once you see this,
More informationChapter 4 - Moving Charges and Magnetism. Magnitude of the magnetic field at the centre of the coil is given by the relation,
Question 4.1: A circular coil of wire consisting of 100 turns, each of radius 8.0 cm carries a current of 0.40 A. What is the magnitude of the magnetic field B at the centre of the coil? Number of turns
More informationSchematic for resistivity measurement
Module 9 : Experimental probes of Superconductivity Lecture 1 : Experimental probes of Superconductivity - I Among the various experimental methods used to probe the properties of superconductors, there
More informationIntroduction to accelerators for teachers (Korean program) Mariusz Sapiński CERN, Beams Department August 9 th, 2012
Introduction to accelerators for teachers (Korean program) Mariusz Sapiński (mariusz.sapinski@cern.ch) CERN, Beams Department August 9 th, 2012 Definition (Britannica) Particle accelerator: A device producing
More informationNovel, Hybrid RF Injector as a High-average. Dinh Nguyen. Lloyd Young
Novel, Hybrid RF Injector as a High-average average-current Electron Source Dinh Nguyen Los Alamos National Laboratory Lloyd Young TechSource Energy Recovery Linac Workshop Thomas Jefferson National Accelerator
More information1P22/1P92 Exam Review Problems 2013 Friday, January 14, :03 AM. Chapter 20
Exam Review Problems 2011 Page 1 1P22/1P92 Exam Review Problems 2013 Friday, January 14, 2011 10:03 AM Chapter 20 True or false? 1 It's impossible to place a charge on an insulator, because no current
More informationVarying accelerating fields
Varying accelerating fields Two approaches for accelerating with time-varying fields Linear Accelerators Circular Accelerators Use many accelerating cavities through which the particle beam passes once.
More informationAcceleration to higher energies
Acceleration to higher energies While terminal voltages of 20 MV provide sufficient beam energy for nuclear structure research, most applications nowadays require beam energies > 1 GeV How do we attain
More informationPhysics (
Question 2.12: A charge of 8 mc is located at the origin. Calculate the work done in taking a small charge of 2 10 9 C from a point P (0, 0, 3 cm) to a point Q (0, 4 cm, 0), via a point R (0, 6 cm, 9 cm).
More informationAnalysis of HIE-ISOLDE cavity results
Analysis of HIE-ISOLDE cavity results A. Miyazaki 1,2, Y. Kadi 1, K. Schirm 1, A. Sublet 1, S. Teixeira 1, M. Therasse 1, W. Venturini Delsolaro 1 1 Organisation européenne pour la recherche nucléaire
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 informationTo be published in the Proceedings of ICEC-22, Seoul Korea, July 2008 MICE Note 232 1
To be published in the Proceedings of ICEC-22, Seoul Korea, 21-25 July 2008 MICE Note 232 1 AC Loss Analysis on the Superconducting Coupling in MICE H. Wu, L. Wang, M. A. Green*, L. K. Li, F. Y. Xu, X.
More informationExpressions for f r (T ) and Q i (T ) from Mattis-Bardeen theory
8 Appendix A Expressions for f r (T ) and Q i (T ) from Mattis-Bardeen theory The Mattis-Bardeen theory of the anomalous skin effect in superconductors [0] may be used to derive the behavior of the resonance
More informationBernhard Holzer, CERN-LHC
Bernhard Holzer, CERN-LHC * Bernhard Holzer, CERN CAS Prague 2014 x Liouville: in reasonable storage rings area in phase space is constant. A = π*ε=const x ε beam emittance = woozilycity of the particle
More informationTools of Particle Physics I Accelerators
Tools of Particle Physics I Accelerators W.S. Graves July, 2011 MIT W.S. Graves July, 2011 1.Introduction to Accelerator Physics 2.Three Big Machines Large Hadron Collider (LHC) International Linear Collider
More informationShort Introduction to CLIC and CTF3, Technologies for Future Linear Colliders
Short Introduction to CLIC and CTF3, Technologies for Future Linear Colliders Explanation of the Basic Principles and Goals Visit to the CTF3 Installation Roger Ruber Collider History p p hadron collider
More informationThe TESLA Dogbone Damping Ring
The TESLA Dogbone Damping Ring Winfried Decking for the TESLA Collaboration April 6 th 2004 Outline The Dogbone Issues: Kicker Design Dynamic Aperture Emittance Dilution due to Stray-Fields Collective
More informationIntroduction to Elementary Particle Physics I
Physics 56400 Introduction to Elementary Particle Physics I Lecture 9 Fall 2018 Semester Prof. Matthew Jones Particle Accelerators In general, we only need classical electrodynamics to discuss particle
More informationIIT-JEE 2012 PAPER - 1 PART - I : PHYSICS. Section I : Single Correct Answer Type
IIT-JEE PAPER - PART - I : PHYSICS Section I : Single Correct Answer Type This section contains multiple choice questions. Each question has four choices (A), (B), (C) and (D) out of which ONLY ONE is
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 informationSUPERCONDUCTING HEAVY ION CYCLOTRON
Atomic Energy of Canada Limited SUPERCONDUCTING HEAVY ION CYCLOTRON by C.B. BIGHAM, J.S. FRASER and H.R. SCHNEIDER Chalk River Nuclear Laboratories Chalk River, Ontario November 1973 AECL-4654 - 1 - SUPERCONDUCTING
More informationTemperature Waves in SRF Research
4th International Particle Accelerator Conference Shanghai China May 12 17, 2013 Temperature Waves in SRF Research Nicholas Valles, Andrei Ganshin, Don Hartill, Georg Hoffstaetter, Xiao Mi and Eric Smith
More informationLecture 6 Current and Resistance Ch. 26
Lecture 6 Current and esistance Ch. 6 Cartoon -nvention of the battery and Voltaic Cell Warm-up problem Topics What is current? Current density Conservation of Current esistance Temperature dependence
More informationPHYS 241 EXAM #2 November 9, 2006
1. ( 5 points) A resistance R and a 3.9 H inductance are in series across a 60 Hz AC voltage. The voltage across the resistor is 23 V and the voltage across the inductor is 35 V. Assume that all voltages
More informationCLASS XII WB SET A PHYSICS
PHYSICS 1. Two cylindrical straight and very long non magnetic conductors A and B, insulated from each other, carry a current I in the positive and the negative z-direction respectively. The direction
More informationConsider a point P on the line joining the two charges, as shown in the given figure.
Question 2.1: Two charges 5 10 8 C and 3 10 8 C are located 16 cm apart. At what point(s) on the line joining the two charges is the electric potential zero? Take the potential at infinity to be zero.
More informationExtraction from cyclotrons. P. Heikkinen
Extraction from cyclotrons P. Heikkinen Classification of extraction schemes Linear accelerators Circular accelerators No extraction problem Constant orbit radius (sychrotrons, betatrons) Increasing orbit
More informationCHAPTER 7 ELECTRODYNAMICS
CHAPTER 7 ELECTRODYNAMICS Outlines 1. Electromotive Force 2. Electromagnetic Induction 3. Maxwell s Equations Michael Faraday James C. Maxwell 2 Summary of Electrostatics and Magnetostatics ρ/ε This semester,
More information18 - ELECTROMAGNETIC INDUCTION AND ALTERNATING CURRENTS ( Answers at the end of all questions ) Page 1
( Answers at the end of all questions ) Page ) The self inductance of the motor of an electric fan is 0 H. In order to impart maximum power at 50 Hz, it should be connected to a capacitance of 8 µ F (
More informationPractice Final C. 1. The diagram below shows a worker using a rope to pull a cart.
1. The diagram below shows a worker using a rope to pull a cart. 6. The graph below represents the relationship between gravitational force and mass for objects near the surface of Earth. The worker s
More informationAcceleration of magnetic dipoles by the sequence of current turns
Acceleration of magnetic dipoles by the sequence of current turns S. N. Dolya Joint Institute for Nuclear Research, Joliot Curie str. 6, Dubna, Russia, 141980 Abstract Acceleration of magnetic dipoles
More informationSecond Year Electromagnetism Summer 2018 Caroline Terquem. Vacation work: Problem set 0. Revisions
Second Year Electromagnetism Summer 2018 Caroline Terquem Vacation work: Problem set 0 Revisions At the start of the second year, you will receive the second part of the Electromagnetism course. This vacation
More informationThe CIS project and the design of other low energy proton synchrotrons
The CIS project and the design of other low energy proton synchrotrons 1. Introduction 2. The CIS project 3. Possible CMS 4. Conclusion S.Y. Lee IU Ref. X. Kang, Ph.D. thesis, Indiana University (1998).
More informationChapter 22. Induction
Chapter 22 Induction Induced emf A current can be produced by a changing magnetic field First shown in an experiment by Michael Faraday A primary coil is connected to a battery A secondary coil is connected
More informationKey-Holes Magnetron Design and Multiphysics Simulation
Key-Holes Magnetron Design and Multiphysics Simulation A. Leggieri* 1, F. Di Paolo 1, and D. Passi 1 1 Department of Electronic Engineering, University of Rome Tor Vergata, Italy *A. Leggieri: Department
More informationChapter 27. Current and Resistance
Chapter 27 Current and Resistance Electric Current Most practical applications of electricity deal with electric currents. The electric charges move through some region of space. The resistor is a new
More informationCLIC ACCELERATING STRUCTURE DEVELOPMENT
CLIC ACCELERATING STRUCTURE DEVELOPMENT W. Wuensch, CERN, Geneva, Switzerland. Abstract One of the most important objectives of the CLIC (Compact Linear Collider) study is to demonstrate the design accelerating
More informationLongitudinal dynamics Yannis PAPAPHILIPPOU CERN
Longitudinal dynamics Yannis PAPAPHILIPPOU CERN United States Particle Accelerator School, University of California - Santa-Cruz, Santa Rosa, CA 14 th 18 th January 2008 1 Outline Methods of acceleration
More informationChem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Spring Semester 2006 Christopher J. Cramer. Lecture 10, February 10, / 4
Chem 350/450 Physical Chemistry II (Quantum Mechanics 3 Credits Spring Semester 006 Christopher J. Cramer Lecture 10, February 10, 006 Solved Homework We are asked to find and for the first two
More informationCavity basics. 1 Introduction. 2 From plane waves to cavities. E. Jensen CERN, Geneva, Switzerland
Cavity basics E. Jensen CERN, Geneva, Switerland Abstract The fields in rectangular and circular waveguides are derived from Maxwell s equations by superposition of plane waves. Subsequently the results
More informationPhysics 1308 Exam 2 Summer 2015
Physics 1308 Exam 2 Summer 2015 E2-01 2. The direction of the magnetic field in a certain region of space is determined by firing a test charge into the region with its velocity in various directions in
More informationModules 5&6. Beam dynamics overview and Linac architecture
Modules 5&6 Beam dynamics overview and Linac architecture 1 Introduction General problem: Need accurate beam optics design to minimize beam loss and emittance growth, to Reduce the requirements on the
More informationDielectric Accelerators at CLARA. G. Burt, Lancaster University On behalf of ASTeC, Lancaster U., Liverpool U., U. Manchester, and Oxford U.
Dielectric Accelerators at CLARA G. Burt, Lancaster University On behalf of ASTeC, Lancaster U., Liverpool U., U. Manchester, and Oxford U. Dielectric Accelerators Types Photonic structures Dielectric
More informationMAGNETIC PROBLEMS. (d) Sketch B as a function of d clearly showing the value for maximum value of B.
PHYS2012/2912 MAGNETC PROBLEMS M014 You can investigate the behaviour of a toroidal (dough nut shape) electromagnet by changing the core material (magnetic susceptibility m ) and the length d of the air
More information1. Consider the biconvex thick lens shown in the figure below, made from transparent material with index n and thickness L.
Optical Science and Engineering 2013 Advanced Optics Exam Answer all questions. Begin each question on a new blank page. Put your banner ID at the top of each page. Please staple all pages for each individual
More informationLecture 5 Notes, Electromagnetic Theory II Dr. Christopher S. Baird, faculty.uml.edu/cbaird University of Massachusetts Lowell
Lecture 5 Notes, Electromagnetic Theory II Dr. Christopher S. Baird, faculty.uml.edu/cbaird University of Massachusetts Lowell 1. Waveguides Continued - In the previous lecture we made the assumption that
More informationHandout 5: Current and resistance. Electric current and current density
1 Handout 5: Current and resistance Electric current and current density Figure 1 shows a flow of positive charge. Electric current is caused by the flow of electric charge and is defined to be equal to
More informationRadiation Heat Transfer Prof. J. Srinivasan Centre for Atmospheric and Oceanic Sciences Indian Institute of Science, Bangalore
Radiation Heat Transfer Prof. J. Srinivasan Centre for Atmospheric and Oceanic Sciences Indian Institute of Science, Bangalore Lecture - 10 Applications In the last lecture, we looked at radiative transfer
More information2010 Physics GA 3: Examination 2
2010 Physics GA 3: Examination 2 GENERAL COMMENTS The number of students who sat for the 2010 Physics examination 2 was 6839. The mean score was 63 per cent; this indicated that students generally found
More informationAccelerators and radiation spectra
X-ray tube Accelerators and radiation spectra Electrons are released from the cathode (negative electrode) by thermionic emission accelerated in an evacuated tube hit the anode (target, positive electrode)
More informationPRELIMINARY EXAMINATION 2018 H2 PHYSICS 9749/01. Paper 1 SEP 2018
PRELIMINARY EXAMINATION 2018 H2 PHYSICS 9749/01 Paper 1 SEP 2018 Additional Materials: Multiple Choice Answer Sheet Duration: 1 hour DO NOT OPEN THIS BOOKLET UNTIL YOU ARE TOLD TO DO SO READ THESE INSTRUCTIONS
More informationPhase Space Study of the Synchrotron Oscillation and Radiation Damping of the Longitudinal and Transverse Oscillations
ScienceAsia 28 (2002 : 393-400 Phase Space Study of the Synchrotron Oscillation and Radiation Damping of the Longitudinal and Transverse Oscillations Balabhadrapatruni Harita*, Masumi Sugawara, Takehiko
More information