Multiscale Modelling, taking into account Collisions
|
|
- Teresa Higgins
- 5 years ago
- Views:
Transcription
1 Multiscale Modelling, taking into account Collisions Andreas Adelmann (Paul Scherrer Institut) October 28, 2017 Multiscale Modelling, taking into account Collisions October 28, 2017 Page 1 / 22
2 1 Motivation 2 Multiscale Nature of the Problem 3 The P 3 M algorithm 4 Disorder Induced Heating Simulation 5 Conclusion Multiscale Modelling, taking into account Collisions October 28, 2017 Page 2 / 22
3 Motivation I 1 model emission of ultra cold electrons 2 understand Coulomb scattering (Borsch effect) Qiang et al. [2007] 3 do we have to worry about it in next generation machines? 4 model non Gaussian tails (high intensity hadron machines) In Qiang et al. [2007] we wrote: Nano-tips with high acceleration gradient around the emission surface have been proposed to generate high brightness beams. However, due to the small size of the tip r = 10 nm, the charge density near the tip is very high even for a small number of electrons. The stochastic Coulomb scattering near the tip can degrade the beam quality and cause extra emittance growth and energy spread. Multiscale Modelling, taking into account Collisions October 28, 2017 Page 3 / 22
4 Motivation II Using a brute force N 2 summation we obtained the following observations: slice emittance over the bunch length 2 higher energy spread 100 higher Multiscale Modelling, taking into account Collisions October 28, 2017 Page 4 / 22
5 Multiscale Nature of the Problem Considering also high intensity high hadron machines (Cyclotrons & FFAGs) length scales: O(10 3 )... O(10 6 ) (m) time scales: (s) energy scales: ev... space charge or emittance dominated Space-Charge effects commonly treated in the mean-field or Vlasov- Poisson/Maxwell approximation. Valid mostly when timescales of interest are much smaller than relaxation timescales. In this approximation, particles move under the influence of external fields and a self-consistent force due to all the other particles. In the electrostatic case, this is the Vlasov-Poisson equation. Multiscale Modelling, taking into account Collisions October 28, 2017 Page 5 / 22
6 Vlasov-Poisson Equation Addressing the multi scale challange When neglecting collisions, and taking advantage of the electrostatic approximation, the Vlasov-Poisson equation describes the (time) evolution of the phase space f(x, v; t) > 0 when considering electromagnetic interaction with charged particles 1. df dt = f t + v xf + q m (E(x, t) + v B(x, t)) vf = 0. (1) Solving with ES-PIC Hockney and Eastwood, h x (t), h y (t), h z (t), M = M x M y M z SAAMG-PCG solver with geometry Adelmann et al. [2010] change M during simulation (many different field solver instances) change t adaptively Toggweiler et al. [2014] modern computational architectures Adelmann et al. [2015] 1 OPAL Multiscale Modelling, taking into account Collisions October 28, 2017 Page 6 / 22
7 Adding (back) Collisions The Boltzmann (Landau) equation Diffusion Damping {}}{{}}{ df dt = (1) = 1 2 (f v i ) + (f v i v j ), v i 2 v i v j the quantities v i and v i v j are the average change in velocity a particle experiences due to collisions with all other particles. large-angle collisions are hard small-angle collisions are soft. Because the Coulomb force is long-ranged, in many cases, the first fluctuation correction to the mean-field approximation comes from the collective contribution of a large number of small-angle collisions. Multiscale Modelling, taking into account Collisions October 28, 2017 Page 7 / 22
8 The P 3 M Algorithm The P 3 M = Particle-Particle + Particle-Mesh is a efficient way to accomplish this task. high resolution from PP part: O(K 2 ), K << N, 1/(x x + ε) good performance from PM part: Φ(x) = G(x, x )ρ(x )d 3 x adjustable influence of Coulomb collisions by fixing K in choosing r c Multiscale Modelling, taking into account Collisions October 28, 2017 Page 8 / 22
9 The P 3 M Algorithm Interaction Splitting I The main concept behind the P 3 M algorithm is a splitting of the interaction function G(r) into a short-range contribution G pp (r) and a long-range contribution G pm (r). This splitting can be done using a Gaussian screening charge distribution G(r) = 1 r = 1 erf(αr) r } {{ } G P P + erf(αr) r } {{ } G P M The long-range interactions are computed on a grid using a particle-mesh solver and the short range interactions are evaluated using an N-body solver with appropriate cutoff.. Multiscale Modelling, taking into account Collisions October 28, 2017 Page 9 / 22
10 The P 3 M Algorithm Interaction Splitting II f(r) r Greens function splitting for r c = 1 and α = 2: G(r) in green, G P M (r) in red and G P P (r) in blue Multiscale Modelling, taking into account Collisions October 28, 2017 Page 10 / 22
11 The P 3 M Algorithm Implementation following [Hockney and Eastwood, 1988] Particle-Mesh (PM): 1 interpolate charges to mesh (CIC, NGP,...) 2 solve for potential Φ using an FFT solver (convolution with greens function) 3 compute forces by F = Φ 4 interpolate forces to particles electric field Particle-Particle (PP): 1 compute linked lists for particles in interaction radius r e, fix K 2 compute short range forces 3 update electric field Multiscale Modelling, taking into account Collisions October 28, 2017 Page 11 / 22
12 The P 3 M Algorithm Implementation following [Hockney and Eastwood, 1988] Particle-Mesh (PM): 1 interpolate charges to mesh (CIC, NGP,...) 2 solve for potential Φ using an FFT solver (convolution with greens function) 3 compute forces by F = Φ 4 interpolate forces to particles electric field Particle-Particle (PP): 1 compute linked lists for particles in interaction radius r e, fix K 2 compute short range forces 3 update electric field Multiscale Modelling, taking into account Collisions October 28, 2017 Page 11 / 22
13 Disorder Induced Heating Simulation (DIH) I B fe acc kt E acc maximum achievable accelerating field kt the temperature of photoemitted electrons, kt < 25 mev Characteristic figures: conventional Cu laser cooled FEL ERL kt (mev) 25 1 N per bunch n 0 (m 3 ) Multiscale Modelling, taking into account Collisions October 28, 2017 Page 12 / 22
14 Disorder Induced Heating Simulation (DIH) II For a near-zero temperature bunch with such high density, we expect some contribution of individual stochastic Coulomb interactions from close encounters just after emission into vacuum to add to the total effective photoemission temperature. We need to determine this amount of stochastic heating as a function of beam density and initial temperature, as well as the nature of its evolution in time. We expect this effect to be most prevalent when the electrostatic potential energy of neighboring particles is comparable to their thermal energy Multiscale Modelling, taking into account Collisions October 28, 2017 Page 13 / 22
15 Disorder Induced Heating Simulation (DIH) III kt the Wigner-Seitz radius. e2 4πε 0 a with a = (3/4πn 0 ) 1/3 for a given density, we expect the heating O(e 2 /4πε 0 a) heating scales with cubic root of n 0 For a rough estimate of the importance of the effect, the plasma coupling parameter Γ can be used, defined as the ratio of kt to the average pair interaction potential. Γ = e2 1 4πε 0 a kt Γ = given an electron temperature of 5 mev, and for the density ranges described above. Multiscale Modelling, taking into account Collisions October 28, 2017 Page 14 / 22
16 Disorder Induced Heating Simulation (DIH) IV Consequence Thus, for applications requiring a large charge density, we expect this stochastic heating to serve as a hard limit to the lowest attainable electron temperature, and limiting the maximum attainable beam brightness. Multiscale Modelling, taking into account Collisions October 28, 2017 Page 15 / 22
17 Disorder Induced Heating Simulation (DIH) V The heating associated with the relaxation of a random, near-zero temperature distribution of charges is well known to the ultracold neutral plasma (UNP) community. In such systems, a cold gas is laser ionized, and after a time on the order of τ = 2πω 1 p e = 2π (n 2 ) 1/2 0 m e ε 0 as a consequence, the temperature T f rises to an equilibrium non zero value i.e. T f > T i. plasma physics terminology: this effect is referred to as DIH T f can be calculated analytically Multiscale Modelling, taking into account Collisions October 28, 2017 Page 16 / 22
18 Problem Setup Loosely connected to Filippetto et al. [2014], for an LBL UED project and Mitchell and Qiang [2015] spherical, cold beam of radius R = µm and charge Q = 25 fc constant focusing applied cubical domain with edge length L = 100 µm P 3 M simulation over 5 plasma periods Boundary conditions: open in x, y periodic in z M = r c varying from 0 µm to µm N = O(10 7 ) P core = 128 simulation over 1000 time-steps Goal compute T f Multiscale Modelling, taking into account Collisions October 28, 2017 Page 17 / 22
19 Simulation Results MSc thesis B. Ulmer Multiscale Modelling, taking into account Collisions October 28, 2017 Page 18 / 22
20 Simulation Results Coulomb collisions are responsible for beam heating! P3M can reproduce analytic results Multiscale Modelling, taking into account Collisions October 28, 2017 Page 19 / 22
21 Conclusion! With P 3 M collisions can be added within acceptable computing cost! The flexibility of turning on/off collisions resembles an other aspect of a multi-scale - or if you preferrer - reduced order model! P 3 M opens the door to S2E simulations with collisions! Low emittance beams at low energies can be accurately modelled Questions: how relevant is this test for: advanced photo cathodes IBS high intensity hardron machines need to better understand the r c as function of beam parameters Acknowledgements: B. Ulmer, C. Mitchell and J. Qiang Multiscale Modelling, taking into account Collisions October 28, 2017 Page 20 / 22
22 Bibliography I A. Adelmann, P. Arbenz, and Y. Ineichen. A fast parallel Poisson solver on irregular domains applied to beam dynamics simulations. 229(12): , A. Adelmann, U. Locans, and A. Suter. The Dynamical Kernel Scheduler - Part 1. ArXiv, Sept D. Filippetto, M. M. Munoz, H. Qian, F. Sannibale, R. Wells, W. Wan, and M. Zolotorev. High repetition rate ultrafast electron diffraction at lbnl. Proceedings of IPAC2014, Dresden, Germany, page 724, R. W. Hockney and J. W. Eastwood. Computer simulation using particles. CRC Press, doi: /phbl C. Mitchell and J. Qiang. A Parallel Particle-Particle, Particle-Mesh Solver for Studying Coulomb Collisions in the Code IMPACT-T. In Proceedings of IPAC2015, Multiscale Modelling, taking into account Collisions October 28, 2017 Page 21 / 22
23 Bibliography II J. Qiang, J. Corlett, S. Lidia, H. A. Padmore, W. Wan, A. Zholents, and A. Adelmann. Numerical study of coulomb scattering effects on electron beam from a nano-tip. In Proceedings of the IEEE Particle Accelerator Conference, M. Toggweiler, A. Adelmann, P. Arbenz, and J. Yang. A novel adaptive time stepping variant of the boris buneman integrator for the simulation of particle accelerators with space charge. JCP, 3(273): , Multiscale Modelling, taking into account Collisions October 28, 2017 Page 22 / 22
Granularity effects Bas van der Geer
Granularity effects Bas van der Geer Marieke de Loos Steinar Wouters Edgar Vredenbregt Jom Luiten Eindhoven University The Netherlands Globular cluster Messier 2 by Hubble Space Telescope.. Located in
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 informationOPAL Introduction Lecture 1
OPAL Introduction Lecture 1 Andreas Adelmann PSI OPAL lectures FFAG Workshop 2014 OPAL Introduction Lecture 1 PSI OPAL lectures FFAG Workshop 2014 Page 1 / 46 General Introduction Models in OPAL Architecture
More informationIntra Beam scattering at BESSY II, MLS and BESSY VSR
Intra Beam scattering at BESSY II, MLS and BESSY VSR T. Mertens 1, T. Atkinson 1, J. Feikes 1, P. Goslawski 1, A. Jankowiak 1 J. G. Hwang 1, J. Li 1, D. Malyutin 1, Y. Petenev 1, M. Ries 1, I. Seiler 1
More informationBasic Energy Science Workshop on Future Electron Sources
Basic Energy Science Workshop on Future Electron Sources Workshop Co-Organizers Pietro Musumeci (UCLA) & Xijie Wang (SLAC) DOE BES contact: Eliane Lessner Workshop charge Electron sources represent a critical
More informationFast Space Charge Calculations with a Multigrid Poisson Solver & Applications
Fast Space Charge Calculations with a Multigrid Poisson Solver & Applications Gisela Pöplau Ursula van Rienen Rostock University DESY, Hamburg, April 26, 2005 Bas van der Geer Marieke de Loos Pulsar Physics
More informationComputer simulation methods (2) Dr. Vania Calandrini
Computer simulation methods (2) Dr. Vania Calandrini in the previous lecture: time average versus ensemble average MC versus MD simulations equipartition theorem (=> computing T) virial theorem (=> computing
More information3D Space Charge Routines: The Software Package MOEVE and FFT Compared
3D Space Charge Routines: The Software Package MOEVE and FFT Compared Gisela Pöplau DESY, Hamburg, December 4, 2007 Overview Algorithms for 3D space charge calculations Properties of FFT and iterative
More informationDiagnostic Systems for Characterizing Electron Sources at the Photo Injector Test Facility at DESY, Zeuthen site
1 Diagnostic Systems for Characterizing Electron Sources at the Photo Injector Test Facility at DESY, Zeuthen site Sakhorn Rimjaem (on behalf of the PITZ team) Motivation Photo Injector Test Facility at
More informationA Fast and Scalable Low Dimensional Solver for Charged Particle Dynamics in Large Particle Accelerators
A Fast and Scalable Low Dimensional Solver for Charged Particle Dynamics in Large Particle Accelerators Yves Ineichen 1,2,3, Andreas Adelmann 2, Costas Bekas 3, Alessandro Curioni 3, Peter Arbenz 1 1 ETH
More informationSRF GUN CHARACTERIZATION - PHASE SPACE AND DARK CURRENT MEASUREMENTS AT ELBE*
SRF GUN CHARACTERIZATION - PHASE SPACE AND DARK CURRENT MEASUREMENTS AT ELBE* E. Panofski #, A. Jankowiak, T. Kamps, Helmholtz-Zentrum Berlin, Berlin, Germany P.N. Lu, J. Teichert, Helmholtz-Zentrum Dresden-Rossendorf,
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 informationComputer Simulation Using Particles
Computer Simulation Using Particles R W Hockney Emeritus Professor, University ofreading and J W Eastwood Culham Laboratory, Abingdon Institute of Physics Publishing Bristol and Philadelphia CONTENTS Foreword
More informationNeighbor Tables Long-Range Potentials
Neighbor Tables Long-Range Potentials Today we learn how we can handle long range potentials. Neighbor tables Long-range potential Ewald sums MSE485/PHY466/CSE485 1 Periodic distances Minimum Image Convention:
More informationLecture 4: The particle equations (1)
Lecture 4: The particle equations (1) Presenter: Mark Eric Dieckmann Department of Science and Technology (ITN), Linköping University, Sweden July 17, 2014 Overview We have previously discussed the leapfrog
More informationA particle-in-cell method with adaptive phase-space remapping for kinetic plasmas
A particle-in-cell method with adaptive phase-space remapping for kinetic plasmas Bei Wang 1 Greg Miller 2 Phil Colella 3 1 Princeton Institute of Computational Science and Engineering Princeton University
More informationNew Electron Source for Energy Recovery Linacs
New Electron Source for Energy Recovery Linacs Ivan Bazarov 20m Cornell s photoinjector: world s brightest electron source 1 Outline Uses of high brightness electron beams Physics of brightness High brightness
More informationBeam-Beam Simulations for e + e Colliders
Beam-Beam Simulations for e + e Colliders J. Rogers, Cornell University Motivation for simulations from observed effects Methods algorithms, simulation codes, and some results. Conclusions Brown Bag Seminar,
More informationExperimental Optimization of Electron Beams for Generating THz CTR and CDR with PITZ
Experimental Optimization of Electron Beams for Generating THz CTR and CDR with PITZ Introduction Outline Optimization of Electron Beams Calculations of CTR/CDR Pulse Energy Summary & Outlook Prach Boonpornprasert
More informationarxiv: v1 [physics.acc-ph] 21 Nov 2011
arxiv:1111.4971v1 [physics.acc-ph] 21 Nov 2011 On FFT based convolutions and correlations, with application to solving Poisson s equation in an open rectangular pipe Abstract Robert D. Ryne Lawrence Berkeley
More informationDark Current at Injector. Jang-Hui Han 27 November 2006 XFEL Beam Dynamics Meeting
Dark Current at Injector Jang-Hui Han 27 November 2006 XFEL Beam Dynamics Meeting Considerations for the guns Ultra-low slice emittance of electron beams higher gradient at the gun cavity solenoid field
More informationTime resolved transverse and longitudinal phase space measurements at the high brightness photo injector PITZ
Time resolved transverse and longitudinal phase space measurements at the high brightness photo injector PITZ 1. Motivation 2. Transverse deflecting structure 3. Longitudinal phase space tomography 4.
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 informationEXPERIMENTAL STUDIES WITH SPATIAL GAUSSIAN-CUT LASER FOR THE LCLS PHOTOCATHODE GUN*
SLAC-PUB-14571 EXPERIMETAL STUDIES WITH SPATIAL GAUSSIA-CUT LASER FOR THE LCLS PHOTOCATHODE GU* F. Zhou +, A. Brachmann, P. Emma, S. Gilevich, and Z. Huang SLAC ational Accelerator Laboratory, 575 Sand
More informationELECTRON DYNAMICS WITH SYNCHROTRON RADIATION
ELECTRON DYNAMICS WITH SYNCHROTRON RADIATION Lenny Rivkin Ecole Polythechnique Federale de Lausanne (EPFL) and Paul Scherrer Institute (PSI), Switzerland CERN Accelerator School: Introduction to Accelerator
More informationUltra-Cold Plasma: Ion Motion
Ultra-Cold Plasma: Ion Motion F. Robicheaux Physics Department, Auburn University Collaborator: James D. Hanson This work supported by the DOE. Discussion w/ experimentalists: Rolston, Roberts, Killian,
More informationExternal injection of electron bunches into plasma wakefields
External injection of electron bunches into plasma wakefields Studies on emittance growth and bunch compression p x x Timon Mehrling, Julia Grebenyuk and Jens Osterhoff FLA, Plasma Acceleration Group (http://plasma.desy.de)
More informationA facility for Femtosecond Soft X-Ray Imaging on the Nanoscale
A facility for Femtosecond Soft X-Ray Imaging on the Nanoscale Jan Lüning Outline Scientific motivation: Random magnetization processes Technique: Lensless imaging by Fourier Transform holography Feasibility:
More informationM2 TP. Low-Energy Electron Diffraction (LEED)
M2 TP Low-Energy Electron Diffraction (LEED) Guide for report preparation I. Introduction: Elastic scattering or diffraction of electrons is the standard technique in surface science for obtaining structural
More informationRecent Status of Polarized Electron Sources at Nagoya University
Recent Status of Polarized Electron Sources at Nagoya University M. Kuwahara, N. Yamamoto, F. Furuta, T. Nakanishi, S. Okumi, M. Yamamoto, M. Kuriki *, T. Ujihara ** and K. Takeda ** Graduate School of
More informationSimulations of the Microbunching Instability in FEL Beam Delivery Systems
Simulations of the Microbunching Instability in FEL Beam Delivery Systems Ilya Pogorelov Tech-X Corporation Workshop on High Average Power & High Brightness Beams UCLA, January 2009 Outline The setting:
More informationPHYSICS OF HOT DENSE PLASMAS
Chapter 6 PHYSICS OF HOT DENSE PLASMAS 10 26 10 24 Solar Center Electron density (e/cm 3 ) 10 22 10 20 10 18 10 16 10 14 10 12 High pressure arcs Chromosphere Discharge plasmas Solar interior Nd (nω) laserproduced
More informationLCLS-II SCRF start-to-end simulations and global optimization as of September Abstract
SLAC National Accelerator Lab LCLS-II TN-17-4 February 217 LCLS-II SCRF start-to-end simulations and global optimization as of September 216 G. Marcus SLAC, Menlo Park, CA 9425 J. Qiang LBNL, Berkeley,
More informationKinetic Solvers with Adaptive Mesh in Phase Space for Low- Temperature Plasmas
Kinetic Solvers with Adaptive Mesh in Phase Space for Low- Temperature Plasmas Vladimir Kolobov, a,b,1 Robert Arslanbekov a and Dmitry Levko a a CFD Research Corporation, Huntsville, AL 35806, USA b The
More informationPlasma Traps for Space-Charge Studies : Status and Perspectives. Hiromi Okamoto (Hiroshima Univ.)
Plasma Traps for Space-Charge Studies : Status and Perspectives Hiromi Okamoto (Hiroshima Univ.) Acknowlegements Past and Present Contributors to the S-POD Project Hiroshima University Hiroyuki Higaki,
More informationPlan of the lectures
Plan of the lectures 1. Introductory remarks on metallic nanostructures Relevant quantities and typical physical parameters Applications. Linear electron response: Mie theory and generalizations 3. Nonlinear
More informationModeling of the secondary electron emission in rf photocathode guns
Modeling of the secondary electron emission in rf photocathode guns J.-H. Han, DESY Zeuthen 8 June 2004 Joint Uni. Hamburg and DESY Accelerator Physics Seminar Contents 1. Necessity of secondary electron
More informationWorkshops on X-band and high gradients: collaboration and resource
Workshops on X-band and high gradients: collaboration and resource 25 October 2012 International workshop on breakdown science and high gradient technology 18-20 April 2012 in KEK 25 October 2012 International
More informationMODELING THE ELECTROMAGNETIC FIELD
MODELING THE ELECTROMAGNETIC FIELD IN ANISOTROPIC INHOMOGENEOUS MAGNETIZED PLASMA OF ECR ION SOURCES G. Torrisi (1), D. Mascali (1), A. Galatà (2), G. Castro (1), L. Celona (1), L. Neri (1), G. Sorbello
More informationAn Adventure in Marrying Laser Arts and Accelerator Technologies
An Adventure in Marrying Laser Arts and Accelerator Technologies Dao Xiang Beam Physics Dept, SLAC, Stanford University Feb-28-2012 An example sample Probe (electron) Pump (laser) Typical pump-probe experiment
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 informationChallenges in Simulating MW Beams in Cyclotrons
Challenges in Simulating MW Beams in Cyclotrons Y. J. Bi (CIAE & PSI & Tsinghua Univ.), A. Adelmann, R. Dölling, M. Humbel, W. Joho, M. Seidel (PSI), T. Zhang (CIAE) HB2010 Morschach, September 28, 2010
More informationReduction of thermal emittance of rf guns *
SLAC-PUB-884 LCLS TN 99-8 October 1999 Reduction of thermal emittance of rf guns * J. E. Clendenin, T. Kotseroglou, G. A. Mulhollan, D. T. Palmer, and J. F. Schmerge Stanford Linear Accelerator Center,
More informationUltracold Plasmas ( Zafar Yasin)
Ultracold Plasmas ( Zafar Yasin) Outline - Creation and why considered important? - Characterization. - Modeling. -My Past Research. - Current Research Trends. Why this topic? - Ultracold Plasmas of astrophysical
More informationBeam dynamics studies for PITZ using a 3D full-wave Lienard-Wiechert PP code
Beam dynamics studies for PITZ using a 3D full-wave Lienard-Wiechert PP code Y. Chen, E. Gjonaj, H. De Gersem, T. Weiland TEMF, Technische Universität Darmstadt, Germany DESY-TEMF Collaboration Meeting
More informationA Fast N-Body Solver for the Poisson(-Boltzmann) Equation
A Fast N-Body Solver for the Poisson(-Boltzmann) Equation Robert D. Skeel Departments of Computer Science (and Mathematics) Purdue University http://bionum.cs.purdue.edu/2008december.pdf 1 Thesis Poisson(-Boltzmann)
More informationExpansion and Equilibration of Ultracold Neutral Plasmas
Expansion and Equilibration of Ultracold Neutral Plasmas Thomas C. Killian Department of Physics and Astronomy 9 th NNP Conference, Columbia University Access Ultracold Temperatures with Laser Cooled Strontium
More informationComputational Challenges in High Intensity and High Brightness Beam Physics
Computational Challenges in High Intensity and High Brightness Beam Physics Robert D. Ryne Lawrence Berkeley National Laboratory Oct 10, 2005 Presented at PAHBEB2005 Erice, Sicily Acknowledgement Thanks
More informationNon-bonded interactions
speeding up the number-crunching continued Marcus Elstner and Tomáš Kubař December 3, 2013 why care? number of individual pair-wise interactions bonded interactions proportional to N: O(N) non-bonded interactions
More informationLecture 3: Optical Properties of Insulators, Semiconductors, and Metals. 5 nm
Metals Lecture 3: Optical Properties of Insulators, Semiconductors, and Metals 5 nm Course Info Next Week (Sept. 5 and 7) no classes First H/W is due Sept. 1 The Previous Lecture Origin frequency dependence
More informationMonte Carlo Collisions in Particle in Cell simulations
Monte Carlo Collisions in Particle in Cell simulations Konstantin Matyash, Ralf Schneider HGF-Junior research group COMAS : Study of effects on materials in contact with plasma, either with fusion or low-temperature
More informationxkcd.com It IS about physics. It ALL is.
xkcd.com It IS about physics. It ALL is. Introduction to Space Plasmas The Plasma State What is a plasma? Basic plasma properties: Qualitative & Quantitative Examples of plasmas Single particle motion
More informationTheory of electron cooling
Theory of electron cooling Daria Astapovych 03/12/2014 HSC Meeting Outline Motivation and idea of the particle beam cooling Cooler Low energy, high energy beam Electron beam Kinetics of electron cooling
More informationTowards a Low Emittance X-ray FEL at PSI
Towards a Low Emittance X-ray FEL at PSI A. Adelmann, A. Anghel, R.J. Bakker, M. Dehler, R. Ganter, C. Gough, S. Ivkovic, F. Jenni, C. Kraus, S.C. Leemann, A. Oppelt, F. Le Pimpec, K. Li, P. Ming, B. Oswald,
More informationSimulation of laser propagation in plasma chamber including nonlinearities by utilization of VirtualLab 5 software
Simulation of laser propagation in plasma chamber including nonlinearities by utilization of VirtualLab 5 software DESY Summer Student Programme, 2012 Anusorn Lueangaramwong Chiang Mai University, Thailand
More informationEmittance and Quantum Efficiency Measurements from a 1.6 cell S- Band Photocathode RF Gun with Mg Cathode *
LCLS-TN-4-3 SLAC PUB 763 September, 4 Emittance and Quantum Efficiency Measurements from a.6 cell S- Band Photocathode RF Gun with Mg Cathode * J.F. Schmerge, J.M. Castro, J.E. Clendenin, D.H. Dowell,
More informationShielded Scintillator for Neutron Characterization
Shielded Scintillator for Neutron Characterization A Thesis Submitted in Partial Fulfillment of the Requirements for Graduation with Research Distinction in Engineering Physics By Patrick X. Belancourt
More informationPlasma Astrophysics Chapter 1: Basic Concepts of Plasma. Yosuke Mizuno Institute of Astronomy National Tsing-Hua University
Plasma Astrophysics Chapter 1: Basic Concepts of Plasma Yosuke Mizuno Institute of Astronomy National Tsing-Hua University What is a Plasma? A plasma is a quasi-neutral gas consisting of positive and negative
More informationTight-Focusing of Short Intense Laser Pulses in Particle-in-Cell Simulations of Laser-Plasma Interaction
16/05/2017, CTU in Prague Tight-Focusing of Short Intense Laser Pulses in Particle-in-Cell Simulations of Laser-Plasma Interaction Bc. Petr Valenta (petr.valenta@eli-beams.eu) Supervisors: doc. Ing. Ondrej
More informationSimulating experiments for ultra-intense laser-vacuum interaction
Simulating experiments for ultra-intense laser-vacuum interaction Nina Elkina LMU München, Germany March 11, 2011 Simulating experiments for ultra-intense laser-vacuum interaction March 11, 2011 1 / 23
More informationPhotoelectric Effect Experiment
Experiment 1 Purpose The photoelectric effect is a key experiment in modern physics. In this experiment light is used to excite electrons that (given sufficient energy) can escape from a material producing
More informationBremsstrahlung. Rybicki & Lightman Chapter 5. Free-free Emission Braking Radiation
Bremsstrahlung Rybicki & Lightman Chapter 5 Bremsstrahlung Free-free Emission Braking Radiation Radiation due to acceleration of charged particle by the Coulomb field of another charge. Relevant for (i)
More informationBeam-Beam Simulation for PEP-II
Beam-Beam Simulation for PEP-II Yunhai Cai Accelerator Research Department-A, SLAC e + e Factories 2003, SLAC October 14, 2003 Outline Method of simulation: - Boundary condition - Parallel computing -
More informationFast Simulation of FEL Linacs with Collective Effects. M. Dohlus FLS 2018
Fast Simulation of FEL Linacs with Collective Effects M. Dohlus FLS 2018 A typical X-FEL gun environment photo cathode cavity, solenoid, drift straight cavity, quadrupole, drift dispersive bend, quadrupole,
More informationCluster fusion in a high magnetic field
Santa Fe July 28, 2009 Cluster fusion in a high magnetic field Roger Bengtson, Boris Breizman Institute for Fusion Studies, Fusion Research Center The University of Texas at Austin In collaboration with:
More informationSimulation of Coulomb Collisions in Plasma Accelerators for Space Applications
Simulation of Coulomb Collisions in Plasma Accelerators for Space Applications D. D Andrea 1, W.Maschek 1 and R. Schneider 2 Vienna, May 6 th 2009 1 Institut for Institute for Nuclear and Energy Technologies
More informationThe Gamma Factory proposal for CERN
The Gamma Factory proposal for CERN Photon-2017 Conference, May 2017 Mieczyslaw Witold Krasny LPNHE, CNRS and University Paris Sorbonne 1 The Gamma Factory in a nutshell Accelerate and store high energy
More informationIntroduction. Thermoionic gun vs RF photo gun Magnetic compression vs Velocity bunching. Probe beam design options
Introduction Following the 19/05/04 meeting at CERN about the "CTF3 accelerated programme", a possible french contribution has been envisaged to the 200 MeV Probe Beam Linac Two machine options were suggested,
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 informationApplication of Wavelets to N body Particle In Cell Simulations
Application of Wavelets to N body Particle In Cell Simulations Balša Terzić (NIU) Northern Illinois University, Department of Physics April 7, 2006 Motivation Gaining insight into the dynamics of multi
More informationParticle Accelerators
Experimental Methods of Particle Physics Particle Accelerators Andreas Streun, PSI andreas.streun@psi.ch https://ados.web.psi.ch/empp-streun Andreas Streun, PSI 1 Particle Accelerators 1. Introduction
More informationBrightness and Coherence of Synchrotron Radiation and Free Electron Lasers. Zhirong Huang SLAC, Stanford University May 13, 2013
Brightness and Coherence of Synchrotron Radiation and Free Electron Lasers Zhirong Huang SLAC, Stanford University May 13, 2013 Introduction GE synchrotron (1946) opened a new era of accelerator-based
More informationAnalysis of FEL Performance Using Brightness Scaled Variables
Analysis of FEL Performance Using Brightness Scaled Variables Michael Gullans with G. Penn, J. Wurtele, and M. Zolotorev Lawrence Berkeley National Laboratory, Berkeley, CA 94720 Outline Introduce brightness
More informationDESY/TEMF Meeting - Status 2011
Erion Gjonaj Technische Universität Darmstadt, Computational Electromagetics Laboratory (TEMF) Schlossgartenstr. 8 64289 Darmstadt, Germany DESY/TEMF Meeting - Status 211 DESY, Hamburg, 16.12.211 Contents
More informationX-ray Free-electron Lasers
X-ray Free-electron Lasers Ultra-fast Dynamic Imaging of Matter II Ischia, Italy, 4/30-5/3/ 2009 Claudio Pellegrini UCLA Department of Physics and Astronomy Outline 1. Present status of X-ray free-electron
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 informationarxiv: v1 [physics.acc-ph] 13 Apr 2016
Multi-objective Optimizations of a Normal Conducting RF Gun Based Ultra Fast Electron Diffraction Beamline Colwyn Gulliford, Adam Bartnik, and Ivan Bazarov CLASSE, Cornell University, 6 Synchrotron Drive
More informationEarly time dynamics of strongly coupled ultracold neutral Ca + and Ca 2+ plasmas
Early time dynamics of strongly coupled ultracold neutral Ca + and Ca 2+ plasmas M. Lyon and S. D. Bergeson Department of Physics and Astronomy, Brigham Young University, Provo, UT 84602, USA For interacting
More informationxkcd.com It IS about physics. It ALL is.
xkcd.com It IS about physics. It ALL is. Introduction to Space Plasmas! The Plasma State What is a plasma? Basic plasma properties: Qualitative & Quantitative Examples of plasmas! Single particle motion
More informationSingle-shot Ultrafast Electron Microscopy
Single-shot Ultrafast Electron Microscopy Renkai Li and Pietro Musumeci Department of Physics and Astronomy, UCLA 25 th North American Particle Accelerator Conference Sep 30 - Oct 4, 2013, Pasadena, CA,
More informationParticle In Cell Beam Dynamics Simulations With a Wavelet based Poisson Solver
Particle In Cell Beam Dynamics Simulations With a Wavelet based Poisson Solver Balša Terzić Beam Physics and Astrophysics Group, Department of Physics Northern Illinois University in collaboration with
More informationGeneration and characterization of ultra-short electron and x-ray x pulses
Generation and characterization of ultra-short electron and x-ray x pulses Zhirong Huang (SLAC) Compact XFEL workshop July 19-20, 2010, Shanghai, China Ultra-bright Promise of XFELs Ultra-fast LCLS Methods
More informationBeam-plasma Physics Working Group Summary
Beam-plasma Physics Working Group Summary P. Muggli, Ian Blumenfeld Wednesday: 10:55, Matt Thompson, LLNL, "Prospect for ultra-high gradient Cherenkov wakefield accelerator experiments at SABER 11:25,
More informationγmy =F=-2πn α e 2 y or y +ω β2 y=0 (1)
Relativistic Weibel Instability Notes from a tutorial at the UCLA Winter School, January 11, 2008 Tom Katsouleas USC Viterbi School of Engineering, LA, CA 90089-0271 Motivation: Weibel instability of relativistic
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 informationPerformance Metrics of Future Light Sources. Robert Hettel, SLAC ICFA FLS 2010 March 1, 2010
Performance Metrics of Future Light Sources Robert Hettel, SLAC ICFA FLS 2010 March 1, 2010 http://www-ssrl.slac.stanford.edu/aboutssrl/documents/future-x-rays-09.pdf special acknowledgment to John Corlett,
More informationLong-range interactions: P 3 M, MMMxD, ELC, MEMD and ICC
Long-range interactions: P 3 M, MMMxD, ELC, MEMD and ICC Axel Arnold Institute for Computational Physics Universität Stuttgart October 10, 2012 Long-range interactions gravity Coulomb interaction dipolar
More informationExternal Injection in Plasma Accelerators. R. Pompili, S. Li, F. Massimo, L. Volta, J. Yang
External Injection in Plasma Accelerators R. Pompili, S. Li, F. Massimo, L. Volta, J. Yang Why Plasma Accelerators? Conventional RF cavities: 50-100 MV/m due to electrical breakdown Plasma: E>100 GV/m
More informationWG2 on ERL light sources CHESS & LEPP
Charge: WG2 on ERL light sources Address and try to answer a list of critical questions for ERL light sources. Session leaders can approach each question by means of (a) (Very) short presentations (b)
More informationPIConGPU Bringing Large-Scale Laser Plasma Simulations to GPU Supercomputing
PIConGPU Bringing Large-Scale Laser Plasma Simulations to GPU Supercomputing Michael Bussmann 1, Heiko Burau 1, René Widera 1, Florian Berninger 1, Axel Hübl 1, Thomas Kluge 1, Alexander Debus 1, Ulrich
More informationOverview of Accelerated Simulation Methods for Plasma Kinetics
Overview of Accelerated Simulation Methods for Plasma Kinetics R.E. Caflisch 1 In collaboration with: J.L. Cambier 2, B.I. Cohen 3, A.M. Dimits 3, L.F. Ricketson 1,4, M.S. Rosin 1,5, B. Yann 1 1 UCLA Math
More informationNon-bonded interactions
speeding up the number-crunching Marcus Elstner and Tomáš Kubař May 8, 2015 why care? key to understand biomolecular structure and function binding of a ligand efficiency of a reaction color of a chromophore
More informationEP2Plus: a hybrid plasma. plume/spacecraft. interaction code. F. Cichocki, M. Merino, E. Ahedo
EP2Plus: a hybrid plasma plume/spacecraft interaction code F. Cichocki, M. Merino, E. Ahedo 24 th SPINE meeting ESTEC, Noordwijk, October 23 rd, 2017 Contents Introduction to EP2PLUS Overall structure
More informationTutorial: simulating a rod pinch diode for pulsed radiography with Trak and GamBet
Tutorial: simulating a rod pinch diode for pulsed radiography with Trak and GamBet Stanley Humphries, Copyright 2012 Field Precision PO Box 13595, Albuquerque, NM 87192 U.S.A. Telephone: +1-505-220-3975
More informationThe Luminosity Upgrade at RHIC. G. Robert-Demolaize, Brookhaven National Laboratory
The Luminosity Upgrade at RHIC G. Robert-Demolaize, Brookhaven National Laboratory RHIC accelerator complex: IPAC'15 - May 3-8, 2015 - Richmond, VA, USA 2 The Relativistic Heavy Ion Collider (RHIC) aims
More informationPIC solvers for Intense Beams in Synchrotrons
PIC solvers for Intense Beams in Synchrotrons Oliver Boine-Franenheim TU Darmstadt and GSI Contents o Application and challenges for beam simulations with space charge o D electrostatic Particle-In-Cell
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 informationStructure of Surfaces
Structure of Surfaces C Stepped surface Interference of two waves Bragg s law Path difference = AB+BC =2dsin ( =glancing angle) If, n =2dsin, constructive interference Ex) in a cubic lattice of unit cell
More informationThe VELA and CLARA Test Facilities at Daresbury Laboratory Peter McIntosh, STFC on behalf of the VELA and CLARA Development Teams
The VELA and CLARA Test Facilities at Daresbury Laboratory Peter McIntosh, STFC on behalf of the VELA and CLARA Development Teams Outline VELA & CLARA Accelerators VELA Commissioning VELA Exploitation
More informationProgress Report: New Particle-In-Cell Code For Numerical Simulations of Coherent Synchrotron Radiation
Progress Report: New Particle-In-Cell Code For Numerical Simulations of Coherent Synchrotron Radiation Balša Terzić and Rui Li CASA, Jefferson Lab CASA Group Meeting May 26, 2010 Notice: Authored by Jefferson
More information