ME 595M: Monte Carlo Simulation
|
|
- Vanessa Arnold
- 6 years ago
- Views:
Transcription
1 ME 595M: Monte Carlo Simulation T.S. Fisher Purdue University 1
2 Monte Carlo Simulation BTE solutions are often difficult to obtain Closure issues force the use of BTE moments Inherent approximations Ballistic transport is often not described well by perturbed equilibrium distributions Instead, seek to track representative carriers and use statistical averaging to predict transport behavior Monte Carlo methods are essentially averaging or integration tools In neither case can we easily incorporate quantum (wave) effects 2
3 Scattering Treatment of scattering underlies almost all MC methods Consider electron scattering mechanisms Acoustic deformation potential (ADP) 1 τ Intervalley scattering by phonon absorption 1 τ Others (τ 3-1, τ 4-1, ) 2 13 ( ) Combine with Mathiessen s rule 1 = 2 10 E( p)/ q ( 13) = E( p) / q Γ ( p) = τ i i 1 ( p) (Effective Scattering Rate) 3
4 Combining Scattering Events Lundstrom Fig
5 Equation of Motion The most complex transport equation employed in most MC simulations is Newton s law, F=ma F e dp = = qe dt Collision times are typically much smaller than free-flight times Collisions are treated as instantaneous events Carrier position can be expressed as an integral of velocity t r() t r(0) v(') t dt' = + 0 5
6 Free Flight Consider an electron under the action of an electric field directed along the z-axis px() t = px(0) p () t = p (0) Position relations become y y p () t = p (0) + ( q) E t z z z p (0) xt () = x(0) + x t * m py (0) y() t = y(0) + t * m Et () E(0) zt () = z(0) + ( q) E z E in numerator is energy last equality derives from the parabolic energy band assumption 6
7 Ensemble Scattering Consider a collection of carriers of concentration n CF that have not undergone scattering since time t=0 Assume a constant scattering rate Γ 0 (for now) The time evolution of n CF is expressed mathematically as dncf = Γ0nCF dt solution n () t = n (0) e Γ CF CF 0 t 7
8 Scattering Probability and Time Probability that an electron survives until time t without scattering is ncf () t Γ0t = e n (0) CF Probability that a carrier undergoes its first collision between t and t+dt is the scattering rate times the survival probability Ptdt () =Γ e 0 dt 0 Γ t 8
9 Random Selection Now we wish to choose a random number such that the probability of choosing a number between r and r+dr equals the probability of selecting a collision time between t and t+dt P() r dr = P() t dt let P(r) = 1 for a random number generator between 0 and 1 dr = Γ e 0 dt 0 Γ Γ t t r 1 0 c c = e 1 1 tc = ln(1 r ) = ln( r1 ) Γ c 0 Γ0 Note that r 1 is also a uniformly distributed random number between 0 and 1 Foregoing equation relates random numbers to collision times 9
10 Self-Scattering Lundstrom Fig
11 t c based on all scattering events (including self-scattering) Collision Distributions tc 1 Γ 0 t c based on only real scattering events (not including self-scattering) t c 1 Γ ( p) Lundstrom Fig
12 Identification of Scattering Events After selecting the duration of free flight (t c ), we must properly choose the type of scattering event that occurs Each scattering mechanism can alter the particle s momentum (direction and magnitude) differently Ultimately, the choices must be proportional to the relative likelihoods of occurrence 12
13 Cumulative Scattering Probabilities Consider k possible scattering mechanisms Imagine a cumulative bar chart that includes a sum of the fractional probabilities of each type Choose a uniformly distributed random number r 2 between 0 and 1 Select scattering event l if the following holds l 1 l 1 1 τ ( ) ( ) i 1 i p τ i 1 i p = r = 2 < Γ Γ
14 Graphical Interpretation Lundstrom, Fig
15 Selecting a Final State Generally the most expensive part of the computation Must select the final magnitude of momentum (energy) and its direction Consider particle states immediately before (t c- ) and after (t c+ ) scattering ( ) c * pt ( c + ) p' = 2m E t +ΔE assumes parabolic energy bands ΔE is a fundamental characteristic of the scattering event (e.g., ΔE=0 for elastic scattering, ΔE ~ ħω for events involving phonons) 15
16 Selecting a Final Direction Align local coordinate system with initial momentum vector p Assume azimuthal invariance 2π 0 1 P( β) dβ= 1 P( β) dβ= dβ 2π Thus, we can select azimuthal angle from a uniformly distributed random number r 3 as β=2πr 3 Lundstrom Fig
17 Determination of the Polar Angle Scattering often depends on the polar angle; hence, its treatment is more complicated π 0 P( α) dα = 1 P( α) dα = 2ππ π sin αdα S( p, p' ) dβp' dp' 0 0 S( p, p' )sin αdαdβp' dp If we assume a simple delta-function scattering mechanism [S~δ(E -E)], then 2 2 P( α) dα = sin αdα 2 17
18 Final Polar Angle Now we seek another uniformly distributed random number r 4 between 0 and 1 to satisfy sin αdα Prdr () = P( α) dα = Pr () = 1 2 r 4 α 1 1 dr = sin α' dα ' = ( 1 cos α) cos α = 1 2r Finally, we find p in the local coordinate system as p 'sinαcosβ p' = p 'sin sin α β p 'cosα Last step involves transformation from local to global coordinates 4 18
19 Other Monte Carlo Topics Ensemble vs. incident flux approaches Ensemble follows particles in parallel in a time-stepping procedure Uses superelectron approximation Incident flux follows each particle sequentially from beginning to end Treatment of Coulomb effects Charged particles alter the local field through Coulomb interactions Can be handled by summing individual contributions If foregoing is too onerous, then particle-in-cell-type methods can be applied Must ensure that the size of the time step does not exceed natural fluctuations at the plasma frequency 19
20 Relation between Monte Carlo Simulation at the BTE Consider a 1D slab on infinitesimal thickness that contains a single particle trajectory δ ni (,) r t = δ( r ri ), and δ f ( r, p,) t = δ( r r ()) t δ( p p ()) t Apply the chain rule i i δfi dr d = rδf i i i p + p δf i i i t dt dt p = r δf i i v+ p δf i i qe+ i t coll but we notice that δ f = δf i r, p r, p p δf and i p δf i i i = t coll t now sum over all trajectories f f + r f v qe p f = t t i i coll coll 20
21 Homework Problem Collector Emitter x 1. Two parallel infinite plates 2. Particles leave the same spot on the source plate at a fixed initial velocity 3. Assume no interactions among particles, i.e. neglect the force between particles 4. Particles scatter with a constant cross section 5. Post-collision velocity will be determined by the model provided 6. The spacing is so small that, for each particle, at most one collision would occur 7. Determine the radial and angular distributions on the collector, and the total number of collisions 21
22 Data structure: Code Breakdown What do we use to represent a particle? Something to bind together all the interesting kinetic properties of particles A particle class in C++ A particle structure in C A n by 7 array in C/C++, matlab, or Fortran 77 1 Vx V y Vz x y z 2 Vx Vy V z x y z 3 V x V y V z x y z 22
23 Algorithm 23
24 Algorithm Breakdown Generate the particles at the emitter Set the V x according to the initial energy, assuming y, z velocity components of particles are neligible Set X,Y,Z to zero V x = y z 2 E / m, V = 0, V = 0 X = 0, Y = 0, Z = 0 24
25 What time step to use? Move particles P( Δt) = 1 exp( nσυδt) < 0.1 Recommendation: ~ 1 femtosecond or less L >υδt 25
26 Position Updates x n+ 1 x n = υ x dt y n+ 1 y n = υ y dt z n+ 1 z n = υ z dt 26
27 Collision Calculate the probability: P( Δt) = 1 exp( nσυδt) r [0,1] Select a random number If r<p, then a collision occurs See notes about polar scattering angle and azimuthal angle See notes about post-collision velocity 27
28 Parameters Number of particles: N=50,000 Initial kinetic energy of all particles: E=5eV 1 2 E = m e υ 2 Gap: L = 10 nm Cross-section, σ =10-18 m 2 Number density of media (target particles), n = 2.5x10 25 m -3 Assume that all scattering is elastic 28
29 Record the particle information When the particle hits the collector or the emitter, we assume it is absorbed without other side effects The final particle positions should be stored only for particles that have collided 29
Decays and Scattering. Decay Rates Cross Sections Calculating Decays Scattering Lifetime of Particles
Decays and Scattering Decay Rates Cross Sections Calculating Decays Scattering Lifetime of Particles 1 Decay Rates There are THREE experimental probes of Elementary Particle Interactions - bound states
More informationElectron-phonon scattering (Finish Lundstrom Chapter 2)
Electron-phonon scattering (Finish Lundstrom Chapter ) Deformation potentials The mechanism of electron-phonon coupling is treated as a perturbation of the band energies due to the lattice vibration. Equilibrium
More informationMonte Carlo radiation transport codes
Monte Carlo radiation transport codes How do they work? Michel Maire (Lapp/Annecy) 23/05/2007 introduction to Monte Carlo radiation transport codes 1 Decay in flight (1) An unstable particle have a time
More informationMonte Carlo radiation transport codes
Monte Carlo radiation transport codes How do they work? Michel Maire (Lapp/Annecy) 16/09/2011 introduction to Monte Carlo radiation transport codes 1 Outline From simplest case to complete process : Decay
More informationLow field mobility in Si and GaAs
EE30 - Solid State Electronics Low field mobility in Si and GaAs In doed samles, at low T, ionized imurity scattering dominates: τ( E) ------ -------------- m N D πe 4 ln( + γ ) ------------- + γ γ E 3
More informationLecture 22: Ionized Impurity Scattering
ECE-656: Fall 20 Lecture 22: Ionized Impurity Scattering Mark Lundstrom Purdue University West Lafayette, IN USA 0/9/ scattering of plane waves ψ i = Ω ei p r U S ( r,t) incident plane wave ( ) = 2π H
More informationChapter V: Interactions of neutrons with matter
Chapter V: Interactions of neutrons with matter 1 Content of the chapter Introduction Interaction processes Interaction cross sections Moderation and neutrons path For more details see «Physique des Réacteurs
More informationStochastic Chemical Kinetics
Stochastic Chemical Kinetics Joseph K Scott November 10, 2011 1 Introduction to Stochastic Chemical Kinetics Consider the reaction I + I D The conventional kinetic model for the concentration of I in a
More informationELECTRONS AND PHONONS IN SEMICONDUCTOR MULTILAYERS
ELECTRONS AND PHONONS IN SEMICONDUCTOR MULTILAYERS В. К. RIDLEY University of Essex CAMBRIDGE UNIVERSITY PRESS Contents Introduction 1 Simple Models of the Electron-Phonon Interaction 1.1 General remarks
More informationAnisotropic fluid dynamics. Thomas Schaefer, North Carolina State University
Anisotropic fluid dynamics Thomas Schaefer, North Carolina State University Outline We wish to extract the properties of nearly perfect (low viscosity) fluids from experiments with trapped gases, colliding
More informationLecture 6 Scattering theory Partial Wave Analysis. SS2011: Introduction to Nuclear and Particle Physics, Part 2 2
Lecture 6 Scattering theory Partial Wave Analysis SS2011: Introduction to Nuclear and Particle Physics, Part 2 2 1 The Born approximation for the differential cross section is valid if the interaction
More informationCHAPTER 2 ELECTRON-PROTON COLLISION
CHAPTER ELECTRON-PROTON COLLISION.1 Electron-proton collision at HERA The collision between electron and proton at HERA is useful to obtain the kinematical values of particle diffraction and interaction
More informationenergy loss Ionization + excitation of atomic energy levels Mean energy loss rate de /dx proportional to (electric charge) 2 of incident particle
Lecture 4 Particle physics processes - particles are small, light, energetic à processes described by quantum mechanics and relativity à processes are probabilistic, i.e., we cannot know the outcome of
More informationPhysics Fall Mechanics, Thermodynamics, Waves, Fluids. Lecture 20: Rotational Motion. Slide 20-1
Physics 1501 Fall 2008 Mechanics, Thermodynamics, Waves, Fluids Lecture 20: Rotational Motion Slide 20-1 Recap: center of mass, linear momentum A composite system behaves as though its mass is concentrated
More informationParallel Ensemble Monte Carlo for Device Simulation
Workshop on High Performance Computing Activities in Singapore Dr Zhou Xing School of Electrical and Electronic Engineering Nanyang Technological University September 29, 1995 Outline Electronic transport
More informationClassical Statistical Mechanics: Part 1
Classical Statistical Mechanics: Part 1 January 16, 2013 Classical Mechanics 1-Dimensional system with 1 particle of mass m Newton s equations of motion for position x(t) and momentum p(t): ẋ(t) dx p =
More informationProblem Set 9: Momentum and Collision Theory. Nov 1 Hour One: Conservation Laws: Momentum and Collision Theory. Reading: YF
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Physics Physics 8.0T Fall Term 2004 Problem Set 9: Momentum and Collision Theory Available on-line October 29; Due: November 9 at 4:00 p.m. Please write
More informationConservation of Linear Momentum : If a force F is acting on particle of mass m, then according to Newton s second law of motion, we have F = dp /dt =
Conservation of Linear Momentum : If a force F is acting on particle of mass m, then according to Newton s second law of motion, we have F = dp /dt = d (mv) /dt where p =mv is linear momentum of particle
More informationProblem Set 9: Momentum and Collision Theory : Solutions
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Physics Physics 8.0T Fall Term 2004 Problem Set 9: Momentum and Collision Theory : Solutions Problem : Impulse and Momentum The compressive force per
More informationJ10M.1 - Rod on a Rail (M93M.2)
Part I - Mechanics J10M.1 - Rod on a Rail (M93M.2) J10M.1 - Rod on a Rail (M93M.2) s α l θ g z x A uniform rod of length l and mass m moves in the x-z plane. One end of the rod is suspended from a straight
More information14. Energy transport.
Phys780: Plasma Physics Lecture 14. Energy transport. 1 14. Energy transport. Chapman-Enskog theory. ([8], p.51-75) We derive macroscopic properties of plasma by calculating moments of the kinetic equation
More informationELECTRONS AND PHONONS IN SEMICONDUCTOR MULTILAYERS
ELECTRONS AND PHONONS IN SEMICONDUCTOR MULTILAYERS Second Edition B.K. RIDLEY University of Essex CAMBRIDGE UNIVERSITY PRESS Contents Preface Introduction 1 Simple Models of the Electron-Phonon Interaction
More informationDetectors in Nuclear Physics: Monte Carlo Methods. Dr. Andrea Mairani. Lectures I-II
Detectors in Nuclear Physics: Monte Carlo Methods Dr. Andrea Mairani Lectures I-II INTRODUCTION Sampling from a probability distribution Sampling from a probability distribution X λ Sampling from a probability
More informationKinetic Monte Carlo (KMC)
Kinetic Monte Carlo (KMC) Molecular Dynamics (MD): high-frequency motion dictate the time-step (e.g., vibrations). Time step is short: pico-seconds. Direct Monte Carlo (MC): stochastic (non-deterministic)
More informationThe dynamics of small particles whose size is roughly 1 µmt or. smaller, in a fluid at room temperature, is extremely erratic, and is
1 I. BROWNIAN MOTION The dynamics of small particles whose size is roughly 1 µmt or smaller, in a fluid at room temperature, is extremely erratic, and is called Brownian motion. The velocity of such particles
More informationFall Quarter 2010 UCSB Physics 225A & UCSD Physics 214 Homework 1
Fall Quarter 2010 UCSB Physics 225A & UCSD Physics 214 Homework 1 Problem 2 has nothing to do with what we have done in class. It introduces somewhat strange coordinates called rapidity and pseudorapidity
More informationDetectors of the Cryogenic Dark Matter Search: Charge Transport and Phonon Emission in Ge 100 Crystals at 40 mk
J Low Temp Phys (2008) 151: 443 447 DOI 10.1007/s10909-007-9666-5 Detectors of the Cryogenic Dark Matter Search: Charge Transport and Phonon Emission in Ge 100 Crystals at 40 mk K.M. Sundqvist B. Sadoulet
More informationSymmetries 2 - Rotations in Space
Symmetries 2 - Rotations in Space This symmetry is about the isotropy of space, i.e. space is the same in all orientations. Thus, if we continuously rotated an entire system in space, we expect the system
More informationThe Exchange Model. Lecture 2. Quantum Particles Experimental Signatures The Exchange Model Feynman Diagrams. Eram Rizvi
The Exchange Model Lecture 2 Quantum Particles Experimental Signatures The Exchange Model Feynman Diagrams Eram Rizvi Royal Institution - London 14 th February 2012 Outline A Century of Particle Scattering
More informationSUMMER SCHOOL - KINETIC EQUATIONS LECTURE II: CONFINED CLASSICAL TRANSPORT Shanghai, 2011.
SUMMER SCHOOL - KINETIC EQUATIONS LECTURE II: CONFINED CLASSICAL TRANSPORT Shanghai, 2011. C. Ringhofer ringhofer@asu.edu, math.la.asu.edu/ chris OVERVIEW Quantum and classical description of many particle
More informationElectrical Transport. Ref. Ihn Ch. 10 YC, Ch 5; BW, Chs 4 & 8
Electrical Transport Ref. Ihn Ch. 10 YC, Ch 5; BW, Chs 4 & 8 Electrical Transport The study of the transport of electrons & holes (in semiconductors) under various conditions. A broad & somewhat specialized
More informationKern- und Teilchenphysik I Lecture 2: Fermi s golden rule
Kern- und Teilchenphysik I Lecture 2: Fermi s golden rule (adapted from the Handout of Prof. Mark Thomson) Prof. Nico Serra Dr. Patrick Owen, Mr. Davide Lancierini http://www.physik.uzh.ch/de/lehre/phy211/hs2017.html
More informationλ φ φ = hc λ ev stop φ = λ φ and now ev stop λ ' = Physics 220 Homework #2 Spring 2016 Due Monday 4/11/16
Physics 0 Homework # Spring 06 Due Monday 4//6. Photons with a wavelength λ = 40nm are used to eject electrons from a metallic cathode (the emitter) by the photoelectric effect. The electrons are prevented
More informationCompound and heavy-ion reactions
Compound and heavy-ion reactions Introduction to Nuclear Science Simon Fraser University Spring 2011 NUCS 342 March 23, 2011 NUCS 342 (Lecture 24) March 23, 2011 1 / 32 Outline 1 Density of states in a
More informationMechanics Physics 151
Mechanics Physics 151 Lecture 4 Continuous Systems and Fields (Chapter 13) What We Did Last Time Built Lagrangian formalism for continuous system Lagrangian L Lagrange s equation = L dxdydz Derived simple
More informationLecture Models for heavy-ion collisions (Part III): transport models. SS2016: Dynamical models for relativistic heavy-ion collisions
Lecture Models for heavy-ion collisions (Part III: transport models SS06: Dynamical models for relativistic heavy-ion collisions Quantum mechanical description of the many-body system Dynamics of heavy-ion
More informationOlivier Bourgeois Institut Néel
Olivier Bourgeois Institut Néel Outline Introduction: necessary concepts: phonons in low dimension, characteristic length Part 1: Transport and heat storage via phonons Specific heat and kinetic equation
More informationPhysics Dec The Maxwell Velocity Distribution
Physics 301 7-Dec-2005 29-1 The Maxwell Velocity Distribution The beginning of chapter 14 covers some things we ve already discussed. Way back in lecture 6, we calculated the pressure for an ideal gas
More informationSelected Topics in the Theory of Heavy Ion Collisions Lecture 1
Selected Topics in the Theory of Heavy Ion Collisions Lecture 1 Urs chim Wiedemann CERN Physics Department TH Division Varenna, 19 July 2010 Based on http://cdsweb.cern.ch/record/1143387/files/p277.pdf
More informationPlasma Optimization in a Multicusp Ion Source by Using a Monte Carlo Simulation
Journal of the Korean Physical Society, Vol. 63, No. 7, October 2013, pp. 0 0 Plasma Optimization in a Multicusp Ion Source by Using a Monte Carlo Simulation M. Hosseinzadeh and H. Afarideh Nuclear Engineering
More information(Super) Fluid Dynamics. Thomas Schaefer, North Carolina State University
(Super) Fluid Dynamics Thomas Schaefer, North Carolina State University Hydrodynamics Hydrodynamics (undergraduate version): Newton s law for continuous, deformable media. Fluids: Gases, liquids, plasmas,...
More informationFermions in the unitary regime at finite temperatures from path integral auxiliary field Monte Carlo simulations
Fermions in the unitary regime at finite temperatures from path integral auxiliary field Monte Carlo simulations Aurel Bulgac,, Joaquin E. Drut and Piotr Magierski University of Washington, Seattle, WA
More informationMONTE CARLO SIMULATION FOR ELECTRON DYNAMICS IN SEMICONDUCTOR DEVICES
Mathematical and Computational Applications, Vol., No., pp. 9-26, 25. Association for Scientific Research MONTE CARLO SIMULATION FOR ELECTRON DYNAMICS IN SEMICONDUCTOR DEVICES Mustafa Akarsu, Ömer Özbaş
More informationA comparative study of no-time-counter and majorant collision frequency numerical schemes in DSMC
Purdue University Purdue e-pubs School of Aeronautics and Astronautics Faculty Publications School of Aeronautics and Astronautics 2012 A comparative study of no-time-counter and majorant collision frequency
More informationMonte Carlo Radiation Transfer I
Monte Carlo Radiation Transfer I Monte Carlo Photons and interactions Sampling from probability distributions Optical depths, isotropic emission, scattering Monte Carlo Basics Emit energy packet, hereafter
More informationLecture: Scattering theory
Lecture: Scattering theory 30.05.2012 SS2012: Introduction to Nuclear and Particle Physics, Part 2 2 1 Part I: Scattering theory: Classical trajectoriest and cross-sections Quantum Scattering 2 I. Scattering
More informationConservation of Momentum. Last modified: 08/05/2018
Conservation of Momentum Last modified: 08/05/2018 Links Momentum & Impulse Momentum Impulse Conservation of Momentum Example 1: 2 Blocks Initial Momentum is Not Enough Example 2: Blocks Sticking Together
More informationPHYS 5012 Radiation Physics and Dosimetry
Radiative PHYS 5012 Radiation Physics and Dosimetry Mean Tuesday 24 March 2009 Radiative Mean Radiative Mean Collisions between two particles involve a projectile and a target. Types of targets: whole
More informationThe Equipartition Theorem
Chapter 8 The Equipartition Theorem Topics Equipartition and kinetic energy. The one-dimensional harmonic oscillator. Degrees of freedom and the equipartition theorem. Rotating particles in thermal equilibrium.
More information22.54 Neutron Interactions and Applications (Spring 2004) Chapter 1 (2/3/04) Overview -- Interactions, Distributions, Cross Sections, Applications
.54 Neutron Interactions and Applications (Spring 004) Chapter 1 (/3/04) Overview -- Interactions, Distributions, Cross Sections, Applications There are many references in the vast literature on nuclear
More information5.62 Physical Chemistry II Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 5.6 Physical Chemistry II Spring 008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 5.6 Spring 008 Lecture #30
More informationBrownian motion and the Central Limit Theorem
Brownian motion and the Central Limit Theorem Amir Bar January 4, 3 Based on Shang-Keng Ma, Statistical Mechanics, sections.,.7 and the course s notes section 6. Introduction In this tutorial we shall
More informationExtensions of the TEP Neutral Transport Methodology. Dingkang Zhang, John Mandrekas, Weston M. Stacey
Extensions of the TEP Neutral Transport Methodology Dingkang Zhang, John Mandrekas, Weston M. Stacey Fusion Research Center, Georgia Institute of Technology, Atlanta, GA 30332-0425, USA Abstract Recent
More informationLow Bias Transport in Graphene: An Introduction
Lecture Notes on Low Bias Transport in Graphene: An Introduction Dionisis Berdebes, Tony Low, and Mark Lundstrom Network for Computational Nanotechnology Birck Nanotechnology Center Purdue University West
More informationLecture 15: Optoelectronic devices: Introduction
Lecture 15: Optoelectronic devices: Introduction Contents 1 Optical absorption 1 1.1 Absorption coefficient....................... 2 2 Optical recombination 5 3 Recombination and carrier lifetime 6 3.1
More information6.730 Physics for Solid State Applications
6.730 Physics for Solid State Applications Lecture 29: Electron-phonon Scattering Outline Bloch Electron Scattering Deformation Potential Scattering LCAO Estimation of Deformation Potential Matrix Element
More information2. Passage of Radiation Through Matter
2. Passage of Radiation Through Matter Passage of Radiation Through Matter: Contents Energy Loss of Heavy Charged Particles by Atomic Collision (addendum) Cherenkov Radiation Energy loss of Electrons and
More information2016 Lloyd G. Elliott University Prize Exam Compiled by the Department of Physics & Astronomy, University of Waterloo
Canadian Association of Physicists SUPPORTING PHYSICS RESEARCH AND EDUCATION IN CANADA 2016 Lloyd G. Elliott University Prize Exam Compiled by the Department of Physics & Astronomy, University of Waterloo
More informationPhysics 9e/Cutnell. correlated to the. College Board AP Physics 2 Course Objectives
correlated to the College Board AP Physics 2 Course Objectives Big Idea 1: Objects and systems have properties such as mass and charge. Systems may have internal structure. Enduring Understanding 1.A:
More informationPhysics 20 Homework 3 SIMS 2016
Physics 20 Homework 3 SIMS 2016 Due: Thursday, August 25 th Special thanks to Sebastian Fischetti for problems 1, 5, and 6. Edits in red made by Keith Fratus. 1. The ballistic pendulum is a device used
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 informationSimulations with MM Force Fields. Monte Carlo (MC) and Molecular Dynamics (MD) Video II.vi
Simulations with MM Force Fields Monte Carlo (MC) and Molecular Dynamics (MD) Video II.vi Some slides taken with permission from Howard R. Mayne Department of Chemistry University of New Hampshire Walking
More informationNotes for Special Relativity, Quantum Mechanics, and Nuclear Physics
Notes for Special Relativity, Quantum Mechanics, and Nuclear Physics 1. More on special relativity Normally, when two objects are moving with velocity v and u with respect to the stationary observer, the
More informationChapter V: Cavity theories
Chapter V: Cavity theories 1 Introduction Goal of radiation dosimetry: measure of the dose absorbed inside a medium (often assimilated to water in calculations) A detector (dosimeter) never measures directly
More informationz F 3 = = = m 1 F 1 m 2 F 2 m 3 - Linear Momentum dp dt F net = d P net = d p 1 dt d p n dt - Conservation of Linear Momentum Δ P = 0
F 1 m 2 F 2 x m 1 O z F 3 m 3 y Ma com = F net F F F net, x net, y net, z = = = Ma Ma Ma com, x com, y com, z p = mv - Linear Momentum F net = dp dt F net = d P dt = d p 1 dt +...+ d p n dt Δ P = 0 - Conservation
More informationPhysics 142 Energy in Mechanics Page 1. Energy in Mechanics
Physics 4 Energy in Mechanics Page Energy in Mechanics This set of notes contains a brief review of the laws and theorems of Newtonian mechanics, and a longer section on energy, because of its central
More informationMCRT: L4 A Monte Carlo Scattering Code
MCRT: L4 A Monte Carlo Scattering Code Plane parallel scattering slab Optical depths & physical distances Emergent flux & intensity Internal intensity moments Constant density slab, vertical optical depth
More informationA Quantum-Classical Approach for the Study of Cascade Processes in Exotic Hydrogen Atoms
PSAS 28 International Conference on Precision Physics of Simple Atomic Systems Windsor, July 21-26, 28 A Quantum-Classical Approach for the Study of Cascade Processes in Exotic Hydrogen Atoms M.P. Faifman
More informationCollision Processes. n n The solution is 0 exp x/ mfp
Collision Processes Collisions mediate the transfer of energy and momentum between various species in a plasma, and as we shall see later, allow a treatment of highly ionized plasma as a single conducting
More information8/31/2018. PHY 711 Classical Mechanics and Mathematical Methods 10-10:50 AM MWF Olin 103
PHY 7 Classical Mechanics and Mathematical Methods 0-0:50 AM MWF Olin 03 Plan for Lecture :. Brief comment on quiz. Particle interactions 3. Notion of center of mass reference fame 4. Introduction to scattering
More informationLaser Induced Control of Condensed Phase Electron Transfer
Laser Induced Control of Condensed Phase Electron Transfer Rob D. Coalson, Dept. of Chemistry, Univ. of Pittsburgh Yuri Dakhnovskii, Dept. of Physics, Univ. of Wyoming Deborah G. Evans, Dept. of Chemistry,
More informationLattice Boltzmann Method for Moving Boundaries
Lattice Boltzmann Method for Moving Boundaries Hans Groot March 18, 2009 Outline 1 Introduction 2 Moving Boundary Conditions 3 Cylinder in Transient Couette Flow 4 Collision-Advection Process for Moving
More informationPhysics 141. Lecture 15. No lecture on Tuesday 10/28. I will be across the Atlantic.
Physics 141. Lecture 15. No lecture on Tuesday 10/28. I will be across the Atlantic. Frank L. H. Wolfs Department of Physics and Astronomy, University of Rochester, Lecture 15, Page 1 Physics 141. Lecture
More informationPhysics 30 - Ballistic Pendulum Lab 2010, Science Kit All Rights Reserved
BACKGROUND Energy The maximum height achieved by the pendulum on the Ballistic Pendulum apparatus can be determined by using the angle it achieved. Figure S1 shows the pendulum in two different positions,
More informationBasic MonteCarlo concepts
Departamento de Física Atómica, Molecular y Nuclear Universidad de Sevilla Introductory tutorial to Geant4 ITA/IEav, São José dos Campos, SP, Brazil July 28 28 - August 1 st, 2014 Index 1 Probability distributions
More informationSupplement: Statistical Physics
Supplement: Statistical Physics Fitting in a Box. Counting momentum states with momentum q and de Broglie wavelength λ = h q = 2π h q In a discrete volume L 3 there is a discrete set of states that satisfy
More informationThe... of a particle is defined as its change in position in some time interval.
Distance is the. of a path followed by a particle. Distance is a quantity. The... of a particle is defined as its change in position in some time interval. Displacement is a.. quantity. The... of a particle
More informationPRESENTED BY: PROF. S. Y. MENSAH F.A.A.S; F.G.A.A.S UNIVERSITY OF CAPE COAST, GHANA.
SOLAR CELL AND ITS APPLICATION PRESENTED BY: PROF. S. Y. MENSAH F.A.A.S; F.G.A.A.S UNIVERSITY OF CAPE COAST, GHANA. OUTLINE OF THE PRESENTATION Objective of the work. A brief introduction to Solar Cell
More informationCHAPTER 9. Microscopic Approach: from Boltzmann to Navier-Stokes. In the previous chapter we derived the closed Boltzmann equation:
CHAPTER 9 Microscopic Approach: from Boltzmann to Navier-Stokes In the previous chapter we derived the closed Boltzmann equation: df dt = f +{f,h} = I[f] where I[f] is the collision integral, and we have
More informationSession 5: Solid State Physics. Charge Mobility Drift Diffusion Recombination-Generation
Session 5: Solid State Physics Charge Mobility Drift Diffusion Recombination-Generation 1 Outline A B C D E F G H I J 2 Mobile Charge Carriers in Semiconductors Three primary types of carrier action occur
More informationA theoretical study of the energy distribution function of the negative hydrogen ion H - in typical
Non equilibrium velocity distributions of H - ions in H 2 plasmas and photodetachment measurements P.Diomede 1,*, S.Longo 1,2 and M.Capitelli 1,2 1 Dipartimento di Chimica dell'università di Bari, Via
More informationPHY492: Nuclear & Particle Physics. Lecture 4 Nature of the nuclear force. Reminder: Investigate
PHY49: Nuclear & Particle Physics Lecture 4 Nature of the nuclear force Reminder: Investigate www.nndc.bnl.gov Topics to be covered size and shape mass and binding energy charge distribution angular momentum
More informationElectrodynamics of Radiation Processes
Electrodynamics of Radiation Processes 7. Emission from relativistic particles (contd) & Bremsstrahlung http://www.astro.rug.nl/~etolstoy/radproc/ Chapter 4: Rybicki&Lightman Sections 4.8, 4.9 Chapter
More informationIntro to Nuclear and Particle Physics (5110)
Intro to Nuclear and Particle Physics (5110) March 23, 2009 From Nuclear to Particle Physics 3/23/2009 1 Nuclear Physics Particle Physics Two fields divided by a common set of tools Theory: fundamental
More information9.3. Total number of phonon modes, total energy and heat capacity
Phys50.nb 6 E = n = n = exp - (9.9) 9... History of the Planck distribution or the Bose-Einstein distribution. his distribution was firstly discovered by Planck in the study of black-body radiation. here,
More informationAnnouncements. 1. Do not bring the yellow equation sheets to the miderm. Idential sheets will be attached to the problems.
Announcements 1. Do not bring the yellow equation sheets to the miderm. Idential sheets will be attached to the problems. 2. Some PRS transmitters are missing. Please, bring them back! 1 Kinematics Displacement
More informationDecay rates and Cross section. Ashfaq Ahmad National Centre for Physics
Decay rates and Cross section Ashfaq Ahmad National Centre for Physics 11/17/2014 Ashfaq Ahmad 2 Outlines Introduction Basics variables used in Exp. HEP Analysis Decay rates and Cross section calculations
More informationThursday Simulation & Unity
Rigid Bodies Simulation Homework Build a particle system based either on F=ma or procedural simulation Examples: Smoke, Fire, Water, Wind, Leaves, Cloth, Magnets, Flocks, Fish, Insects, Crowds, etc. Simulate
More informationNon-Continuum Energy Transfer: Phonons
Non-Continuum Energy Transfer: Phonons D. B. Go Slide 1 The Crystal Lattice The crystal lattice is the organization of atoms and/or molecules in a solid simple cubic body-centered cubic hexagonal a NaCl
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 informationDrude-Schwarzschild Metric and the Electrical Conductivity of Metals
Drude-Schwarzschild Metric and the Electrical Conductivity of Metals P. R. Silva - Retired associate professor Departamento de Física ICEx Universidade Federal de Minas Gerais email: prsilvafis@gmail.com
More informationSecond order Unit Impulse Response. 1. Effect of a Unit Impulse on a Second order System
Effect of a Unit Impulse on a Second order System We consider a second order system mx + bx + kx = f (t) () Our first task is to derive the following If the input f (t) is an impulse cδ(t a), then the
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationSolution of time-dependent Boltzmann equation for electrons in non-thermal plasma
Solution of time-dependent Boltzmann equation for electrons in non-thermal plasma Z. Bonaventura, D. Trunec Department of Physical Electronics Faculty of Science Masaryk University Kotlářská 2, Brno, CZ-61137,
More informationColumbia University Department of Physics QUALIFYING EXAMINATION
Columbia University Department of Physics QUALIFYING EXAMINATION Friday, January 17, 2014 1:00PM to 3:00PM General Physics (Part I) Section 5. Two hours are permitted for the completion of this section
More informationPart One Inelastic Collision:
Problem 3: Experiment 7: Collisions Analysis Part One Inelastic Collision: Analysis: Complete the analysis of your data table by following the two steps below, and answer Question below. You will analyze
More informationScattering in Cold- Cathode Discharges
Simulating Electron Scattering in Cold- Cathode Discharges Alexander Khrabrov, Igor Kaganovich*, Vladimir I. Demidov**, George Petrov*** *Princeton Plasma Physics Laboratory ** Wright-Patterson Air Force
More informationChapter 9: Statistical Mechanics
Chapter 9: Statistical Mechanics Chapter 9: Statistical Mechanics...111 9.1 Introduction...111 9.2 Statistical Mechanics...113 9.2.1 The Hamiltonian...113 9.2.2 Phase Space...114 9.2.3 Trajectories and
More informationCurriculum Map-- Kings School District Honors Physics
Curriculum Map-- Kings School District Honors Physics Big ideas Essential Questions Content Skills/Standards Assessment + Criteria Activities/Resources Motion of an object can be described by its position,
More informationLINEAR RESPONSE THEORY
MIT Department of Chemistry 5.74, Spring 5: Introductory Quantum Mechanics II Instructor: Professor Andrei Tokmakoff p. 8 LINEAR RESPONSE THEORY We have statistically described the time-dependent behavior
More information