Waves are naturally found in plasmas and have to be dealt with. This includes instabilities, fluctuations, waveinduced

Size: px
Start display at page:

Download "Waves are naturally found in plasmas and have to be dealt with. This includes instabilities, fluctuations, waveinduced"

Transcription

1 Lecure 1 Inroducion Why is i imporan o sudy waves in plasma? Waves are naurally found in plasmas and have o be deal wih. This includes insabiliies, flucuaions, waveinduced ranspor... Waves can deliver and/or manipulae energy-momenum in plasma: heaing, curren drive, paricle acceleraion mode sabilizaion, alpha channeling Raman and Brillouin scaering... Plasma diagnosics: inerferomery, Faraday roaion, Thomson scaering... Sudying plasma waves provides insigh ino wave physics in general. Basic equaions Maxwell s equaions A characerisic feaure of plasmas is ha hey conain free charges. Thus, any wave in a plasma is accompanied by oscillaing elecric currens, which unavoidably cause oscillaions of elecromagneic fields. Because of his, waves in plasmas are sudied using Maxwell s equaions. Overall, Maxwell s equaions include four equaions (or eigh equaions, if one couns componens independenly). However, Gauss s laws for elecric and magneic fields, aken ogeher wih charge conservaion, can be considered simply as iniial condiions for Ampere s and Faraday s laws. Thus, i will be enough for us o limi our consideraion o he laer wo, B 4 Π E 1 c c j 1 c B. E, A combinaion of hese equaions gives E 1 c or, equivalenly, 4 Π c j 1 c E 1 2 E c 2 4 Π 2 c 2 E j., Absen he curren densiy on he righ-hand side, his equaion describes ligh propagaion in vacuum. Thus, all feaures specific o plasma are, in fac, conained in he righ-hand side, i.e., he curren densiy. In general,

2 2 lecure-01.nb consiss of he self-consisen curren j ind induced by he field and i can also include an exernal curren j 0 independen of he wave field: j j ind j 0. Below, we consider j 0 prescribed and discuss how o model j ind. General linear media Linear media are defined as hose where j ind scales linearly wih he field E. In some cases, he relaion beween he wo can be as simple as j ind, x Σ, xe, x. This model is commonly known as Ohm s law. I is also called he local-response model, because i assumes ha he curren a, x is deermined only by he field a he same locaion, x. Tha can be a valid assumpion indeed, paricularly in srongly collisional plasma. Bu in general, he curren can be a nonlocal funcion of he field, i.e., exhibi dispersion. Nonlocaliy in ime is called emporal dispersion, and nonlocaliy in space is called spaial dispersion. The mos general expression for a dispersive linear curren can be wrien as follows. Suppose he dynamics is considered a 0, where 0 is some iniial momen of ime. Then, he curren induced by he field since ha ime can be wrien as j ind, x 0 ' 3 x' Σ, x, ', x'e ', x'. Commen: In lieraure, 0 is ofen aken o be o, bu he inegral does no always converge in his limi. Some subleies resuling from his fac will be discussed in his course laer. The ensor Σ, x, ', x', which is called (he ime-space represenaion of) conduciviy, can be defined formally as a funcional derivaive Σ j, x, x, ', x' E ', x' and serves as he weigh funcion ha deermines how much he field E ', x' conribues o he curren j ind, x. The fac ha he inegraion domain is limied o ' reflecs causaliy; namely, he curren a a given ime can be affeced by fields ha exised in he pas bu no in he fuure. In oher words, Σ, x, ', x' 0 for all '. Using j j 0 j ind and subsiuing he above expression for j ind, one ges he field equaion in he following inegro-differenial form: E 1 2 E c 2 4 Π j 0 4 Π 2 c 2 c 2 3 x' ' Σ, x, ', x'e ', x'. For he fuure references, le us represen his equaion also in wo oher forms as follows. Firs of all, one can wrie E 1 2 c 2 E, x ' 3 x' Χ, x, ', x'e ', x' 4 Π j 0 2 c 2, where Χ is called suscepibiliy (in he ime-space represenaion) defined via 0 ' 3 x' Χ, x, ', x'e ', x' 4 Π 0 ' 3 x' Σ, x, ', x'e ', x'. I is also convenien o inroduce he dielecric ensor (in he ime-space represenaion) as

3 lecure-01.nb 3 Ε, x, ', x' ' 0 x x' Χ, x, ', x', so he equaion for E (henceforh called wave equaion for breviy) becomes E 1 2 c 2 ' 3 x' Ε, x, ', x'e ', x' 4 Π j c 2. The added 0 in he ime argumen of he dela funcion denoes an infiniesimal posiive value. Adding his shif ensures ha he ime inegral over he domain 0, is well defined, because 0 ' 0 ' 0 0 Τ 0 Τ 0 Τ 0 Τ Τ Τ 1, whereas 0 Τ Τ is, sricly speaking, undefined. Saionary homogeneous linear media For now, le us limi our consideraion o saionary homogeneous linear media (SHLM). In such media, Σ is invarian wih respec o ime-space ranslaions, so i is allowed o depend only on he difference of ime-space coordinaes, x and ', x' raher han on heir absolue values. In oher words, one can express Σ as Σ, x, ', x' Σ ', x x'. Here, Σ is some ensor funcion, which we will also call conduciviy for breviy. Then, j ind, x 0 ' 3 x' Σ ', x x'e ', x'. Noe ha such Σ can be convenienly inerpreed as a Green s funcion. In oher words, Σ m,n, x can be undersood as he mh componen of he curren, j m, produced by he nh componen of he elecric field, E n, of he form E n ', x' ' 0 x'. The upper limi in he ime inegral can be replaced wih, since Σ is idenically zero a ' in any case. Le us also consider 0. Then, j is simply he convoluion of Σ and E, j ind, x ' 3 x' Σ ', x x'e ', x'. This expression is simplified significanly in he Fourier represenaion, j ind Ω, k ΣΩ, keω, k, so he specral represenaion of he conduciviy, ΣΩ, k, can be undersood as ΣΩ, k jω, k EΩ, k. (Noe ha his expression involves a regular parial derivaive, as opposed o a funcional derivaive ha emerges in he ime-space represenaion.) Explicily, i can be wrien as follows: ΣΩ, k Τ 3 Ω Τ ky y ΣΤ, y 0 Τ 3 y ΣΤ, y Ω Τ ky. The specral represenaions of he corresponding ensors Χ and Ε are as follows: ΧΩ, k 4 Π Ω ΣΩ, k,

4 4 lecure-01.nb ΕΩ, k 1 ΧΩ, k 1 4 Π ΣΩ, k. Ω Then, he field equaion becomes kkeω, k Ω2 or, equivalenly, kkeω, k k 2 EΩ, k Ω2 4 Π Ω ΕΩ, keω, k 2 c c 2 j 0 Ω, k, 4 Π Ω ΕΩ, keω, k 2 c c 2 j 0 Ω, k, Unlike he equaion in he ime-space represenaion, his equaion is no differenial, so i can be sudied easily as follows. Waves wih real frequencies Dispersion relaions Le us consider free oscillaions, which correspond o j 0 0. (A more careful approach will be discussed on he nex lecure.) Then he field equaion can be expressed as i, j Ω, k E j 0, where we inroduced he following marix: i, j c2 Ω 2 k i k j k 2 i, j Ε i, j Ω, k. Equaion (1) saes ha he field specrum in Ω, k space can be nonrivial only on surfaces (or, in onedimensional problems, curves) saisfying de Ω, k 0. This can be considered as an equaion for Ωk, called dispersion relaion. Depending on he number of soluions of he dispersion relaion, called branches, he corresponding waves are aribued as scalar waves ( 1) or vecor waves ( 1). Commen: For example, a ypical second-order (in ) wave equaion for a scalar field (e.g., he Klein-Gordon equaion) acually describes vecor waves, because i allows for boh posiive and negaive frequencies, i.e., corresponds o 2. The oher definiion commonly used in lieraure, scalar waves are hose described by scalar funcions, is ambiguous. For now, le us assume ha all hese soluions are real. Hence he general field can be represened as EΩ, k Ω Ω b k 2 Π E b k, b where E b k are deermined by iniial condiions. Noe ha each of hese mus saisfy vecor equaions Ω b k, ke b k 0. One can hus express he field specral funcion as E b b e b, where b k is some scalar ampliude, and e b is a polarizaion vecor defined via Ω b k, ke b k 0, e b 1. The field in he, x represenaion is obained via he inverse Fourier ransform, (1) (2)

5 lecure-01.nb 5 E, x b 3 k Ω 2 Π 3 2 Π Ω Ω bk 2 Π E b k Ω kx E b, x, E b, x E b k Ω bk kx 3 k 2 Π. 3 In oher words, he field can be decomposed ino independen conribuions from differen branches of he dispersion relaion. If only one branch is excied, a vecor wave hen behaves exacly as a scalar wave. Commen: The laer saemen is rue only in SHLM, whereas in nonsaionary or inhomogeneous medium, differen branches can inerac. Such effecs will be discussed laer in his course. Envelope evoluion. Group and phase velociies Suppose ha only one branch is excied, so Ωk kx E, x Ek 3 k 2 Π, 3 where we omied he branch index b for breviy. Suppose also ha Ek is localized around some k k 0, so we can expand Ωk as Ω Ω 0 v g k, where Ω 0 Ωk 0, k k k 0, and v g k Ωk k k0. This leads o E, x k 0 xv p Ek 0 k kxv g 3 k 2 Π 3, where we inroduced v p k 0 k 0 Ω k 0. Equivalenly, his can be wrien as E, x k 0 xv p Ex v g, where we inroduced Ex Ek 0 k k x 3 k 2 Π 3. In paricular, his implies E0, x k 0x Ex, so Ex represens he iniial envelope. Since, in he expression for E, x, he funcion E depends on, x only hrough he combinaion x v g, he envelope urns ou o be saionary in he frame moving wih velociy v g. Then, in he laboraory frame, he envelope ravels wih velociy v g ; i.e., v g is he velociy of envelope propagaion. I is also called he group velociy. Commen: Laer, we will explicily show ha, for linear waves considered here, v g is also he velociy of he wave acion ( phoons ), energy, and momenum, as well as he velociy of informaion. In conras, nonlinear waves can have heir acion, energy, momenum, and informaion propagaing a differen velociies, so he b

6 6 lecure-01.nb group velociy mus be defined in a more suble manner. These issues are beyond he scope of he course, which will be resriced o linear waves. Similarly, he phase facor in he expression for E, x depends on, x only hrough he combinaion x v p. Thus, v p serves as he velociy of phase frons. I is also called phase velociy. The phase and group velociies can be very differen (e.g., in some cases hey are even direced opposiely o each oher). The geomerical meaning of v g and v p in 1D is also seen in he following figure illusraing he propagaion of a localized pulse (ligher regions correspond o larger E 2 ): I is seen now ha knowing he dispersion relaions Ωk gives informaion no only abou monochromaic oscillaions bu also abou he evoluion of quasimonochromaic waves. (We will also exend his saemen o weakly inhomogeneous and nonsaionary plasmas laer.) In his sense, dispersion relaions deermine all wave physics. Because of his, dispersion relaions will be he primary (albei no he only) focus of our course.

Linear Response Theory: The connection between QFT and experiments

Linear Response Theory: The connection between QFT and experiments Phys540.nb 39 3 Linear Response Theory: The connecion beween QFT and experimens 3.1. Basic conceps and ideas Q: How do we measure he conduciviy of a meal? A: we firs inroduce a weak elecric field E, and

More information

Traveling Waves. Chapter Introduction

Traveling Waves. Chapter Introduction Chaper 4 Traveling Waves 4.1 Inroducion To dae, we have considered oscillaions, i.e., periodic, ofen harmonic, variaions of a physical characerisic of a sysem. The sysem a one ime is indisinguishable from

More information

Class Meeting # 10: Introduction to the Wave Equation

Class Meeting # 10: Introduction to the Wave Equation MATH 8.5 COURSE NOTES - CLASS MEETING # 0 8.5 Inroducion o PDEs, Fall 0 Professor: Jared Speck Class Meeing # 0: Inroducion o he Wave Equaion. Wha is he wave equaion? The sandard wave equaion for a funcion

More information

Two Coupled Oscillators / Normal Modes

Two Coupled Oscillators / Normal Modes Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own

More information

t 2 B F x,t n dsdt t u x,t dxdt

t 2 B F x,t n dsdt t u x,t dxdt Evoluion Equaions For 0, fixed, le U U0, where U denoes a bounded open se in R n.suppose ha U is filled wih a maerial in which a conaminan is being ranspored by various means including diffusion and convecion.

More information

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals

More information

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is UNIT IMPULSE RESPONSE, UNIT STEP RESPONSE, STABILITY. Uni impulse funcion (Dirac dela funcion, dela funcion) rigorously defined is no sricly a funcion, bu disribuion (or measure), precise reamen requires

More information

Some Basic Information about M-S-D Systems

Some Basic Information about M-S-D Systems Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,

More information

Physics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution

Physics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution Physics 7b: Saisical Mechanics Fokker-Planck Equaion The Langevin equaion approach o he evoluion of he velociy disribuion for he Brownian paricle migh leave you uncomforable. A more formal reamen of his

More information

Linear Surface Gravity Waves 3., Dispersion, Group Velocity, and Energy Propagation

Linear Surface Gravity Waves 3., Dispersion, Group Velocity, and Energy Propagation Chaper 4 Linear Surface Graviy Waves 3., Dispersion, Group Velociy, and Energy Propagaion 4. Descripion In many aspecs of wave evoluion, he concep of group velociy plays a cenral role. Mos people now i

More information

Lecture #8 Redfield theory of NMR relaxation

Lecture #8 Redfield theory of NMR relaxation Lecure #8 Redfield heory of NMR relaxaion Topics The ineracion frame of reference Perurbaion heory The Maser Equaion Handous and Reading assignmens van de Ven, Chapers 6.2. Kowalewski, Chaper 4. Abragam

More information

Math 334 Fall 2011 Homework 11 Solutions

Math 334 Fall 2011 Homework 11 Solutions Dec. 2, 2 Mah 334 Fall 2 Homework Soluions Basic Problem. Transform he following iniial value problem ino an iniial value problem for a sysem: u + p()u + q() u g(), u() u, u () v. () Soluion. Le v u. Then

More information

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004 ODEs II, Lecure : Homogeneous Linear Sysems - I Mike Raugh March 8, 4 Inroducion. In he firs lecure we discussed a sysem of linear ODEs for modeling he excreion of lead from he human body, saw how o ransform

More information

4.6 One Dimensional Kinematics and Integration

4.6 One Dimensional Kinematics and Integration 4.6 One Dimensional Kinemaics and Inegraion When he acceleraion a( of an objec is a non-consan funcion of ime, we would like o deermine he ime dependence of he posiion funcion x( and he x -componen of

More information

Maxwell s Equations and Electromagnetic Waves

Maxwell s Equations and Electromagnetic Waves Phsics 36: Waves Lecure 3 /9/8 Maxwell s quaions and lecromagneic Waves Four Laws of lecromagneism. Gauss Law qenc all da ρdv Inegral From From he vecor ideni da dv Therefore, we ma wrie Gauss Law as ρ

More information

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,

More information

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems. Mah 2250-004 Week 4 April 6-20 secions 7.-7.3 firs order sysems of linear differenial equaions; 7.4 mass-spring sysems. Mon Apr 6 7.-7.2 Sysems of differenial equaions (7.), and he vecor Calculus we need

More information

Final Spring 2007

Final Spring 2007 .615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o

More information

Chapter 2. First Order Scalar Equations

Chapter 2. First Order Scalar Equations Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.

More information

10. State Space Methods

10. State Space Methods . Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he

More information

Numerical Dispersion

Numerical Dispersion eview of Linear Numerical Sabiliy Numerical Dispersion n he previous lecure, we considered he linear numerical sabiliy of boh advecion and diffusion erms when approimaed wih several spaial and emporal

More information

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still. Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in

More information

Guest Lectures for Dr. MacFarlane s EE3350 Part Deux

Guest Lectures for Dr. MacFarlane s EE3350 Part Deux Gues Lecures for Dr. MacFarlane s EE3350 Par Deux Michael Plane Mon., 08-30-2010 Wrie name in corner. Poin ou his is a review, so I will go faser. Remind hem o go lisen o online lecure abou geing an A

More information

LAPLACE TRANSFORM AND TRANSFER FUNCTION

LAPLACE TRANSFORM AND TRANSFER FUNCTION CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION Professor Dae Ryook Yang Spring 2018 Dep. of Chemical and Biological Engineering 5-1 Road Map of he Lecure V Laplace Transform and Transfer funcions

More information

Physics for Scientists & Engineers 2

Physics for Scientists & Engineers 2 Direc Curren Physics for Scieniss & Engineers 2 Spring Semeser 2005 Lecure 16 This week we will sudy charges in moion Elecric charge moving from one region o anoher is called elecric curren Curren is all

More information

The Maxwell Equations, the Lorentz Field and the Electromagnetic Nanofield with Regard to the Question of Relativity

The Maxwell Equations, the Lorentz Field and the Electromagnetic Nanofield with Regard to the Question of Relativity The Maxwell Equaions, he Lorenz Field and he Elecromagneic Nanofield wih Regard o he Quesion of Relaiviy Daniele Sasso * Absrac We discuss he Elecromagneic Theory in some main respecs and specifically

More information

( ) = b n ( t) n " (2.111) or a system with many states to be considered, solving these equations isn t. = k U I ( t,t 0 )! ( t 0 ) (2.

( ) = b n ( t) n  (2.111) or a system with many states to be considered, solving these equations isn t. = k U I ( t,t 0 )! ( t 0 ) (2. Andrei Tokmakoff, MIT Deparmen of Chemisry, 3/14/007-6.4 PERTURBATION THEORY Given a Hamilonian H = H 0 + V where we know he eigenkes for H 0 : H 0 n = E n n, we can calculae he evoluion of he wavefuncion

More information

Continuous Time. Time-Domain System Analysis. Impulse Response. Impulse Response. Impulse Response. Impulse Response. ( t) + b 0.

Continuous Time. Time-Domain System Analysis. Impulse Response. Impulse Response. Impulse Response. Impulse Response. ( t) + b 0. Time-Domain Sysem Analysis Coninuous Time. J. Robers - All Righs Reserved. Edied by Dr. Rober Akl 1. J. Robers - All Righs Reserved. Edied by Dr. Rober Akl 2 Le a sysem be described by a 2 y ( ) + a 1

More information

Week 1 Lecture 2 Problems 2, 5. What if something oscillates with no obvious spring? What is ω? (problem set problem)

Week 1 Lecture 2 Problems 2, 5. What if something oscillates with no obvious spring? What is ω? (problem set problem) Week 1 Lecure Problems, 5 Wha if somehing oscillaes wih no obvious spring? Wha is ω? (problem se problem) Sar wih Try and ge o SHM form E. Full beer can in lake, oscillaing F = m & = ge rearrange: F =

More information

Basic Circuit Elements Professor J R Lucas November 2001

Basic Circuit Elements Professor J R Lucas November 2001 Basic Circui Elemens - J ucas An elecrical circui is an inerconnecion of circui elemens. These circui elemens can be caegorised ino wo ypes, namely acive and passive elemens. Some Definiions/explanaions

More information

Unsteady Flow Problems

Unsteady Flow Problems School of Mechanical Aerospace and Civil Engineering Unseady Flow Problems T. J. Craf George Begg Building, C41 TPFE MSc CFD-1 Reading: J. Ferziger, M. Peric, Compuaional Mehods for Fluid Dynamics H.K.

More information

k 1 k 2 x (1) x 2 = k 1 x 1 = k 2 k 1 +k 2 x (2) x k series x (3) k 2 x 2 = k 1 k 2 = k 1+k 2 = 1 k k 2 k series

k 1 k 2 x (1) x 2 = k 1 x 1 = k 2 k 1 +k 2 x (2) x k series x (3) k 2 x 2 = k 1 k 2 = k 1+k 2 = 1 k k 2 k series Final Review A Puzzle... Consider wo massless springs wih spring consans k 1 and k and he same equilibrium lengh. 1. If hese springs ac on a mass m in parallel, hey would be equivalen o a single spring

More information

Ground Rules. PC1221 Fundamentals of Physics I. Kinematics. Position. Lectures 3 and 4 Motion in One Dimension. A/Prof Tay Seng Chuan

Ground Rules. PC1221 Fundamentals of Physics I. Kinematics. Position. Lectures 3 and 4 Motion in One Dimension. A/Prof Tay Seng Chuan Ground Rules PC11 Fundamenals of Physics I Lecures 3 and 4 Moion in One Dimension A/Prof Tay Seng Chuan 1 Swich off your handphone and pager Swich off your lapop compuer and keep i No alking while lecure

More information

The expectation value of the field operator.

The expectation value of the field operator. The expecaion value of he field operaor. Dan Solomon Universiy of Illinois Chicago, IL dsolom@uic.edu June, 04 Absrac. Much of he mahemaical developmen of quanum field heory has been in suppor of deermining

More information

A Note on Fractional Electrodynamics. Abstract

A Note on Fractional Electrodynamics. Abstract Commun Nonlinear Sci Numer Simula 8 (3 589 593 A Noe on Fracional lecrodynamics Hosein Nasrolahpour Absrac We invesigae he ime evoluion o he racional elecromagneic waves by using he ime racional Maxwell's

More information

Theory of! Partial Differential Equations!

Theory of! Partial Differential Equations! hp://www.nd.edu/~gryggva/cfd-course/! Ouline! Theory o! Parial Dierenial Equaions! Gréar Tryggvason! Spring 011! Basic Properies o PDE!! Quasi-linear Firs Order Equaions! - Characerisics! - Linear and

More information

EXERCISES FOR SECTION 1.5

EXERCISES FOR SECTION 1.5 1.5 Exisence and Uniqueness of Soluions 43 20. 1 v c 21. 1 v c 1 2 4 6 8 10 1 2 2 4 6 8 10 Graph of approximae soluion obained using Euler s mehod wih = 0.1. Graph of approximae soluion obained using Euler

More information

Second quantization and gauge invariance.

Second quantization and gauge invariance. 1 Second quanizaion and gauge invariance. Dan Solomon Rauland-Borg Corporaion Moun Prospec, IL Email: dsolom@uic.edu June, 1. Absrac. I is well known ha he single paricle Dirac equaion is gauge invarian.

More information

8. Basic RL and RC Circuits

8. Basic RL and RC Circuits 8. Basic L and C Circuis This chaper deals wih he soluions of he responses of L and C circuis The analysis of C and L circuis leads o a linear differenial equaion This chaper covers he following opics

More information

System of Linear Differential Equations

System of Linear Differential Equations Sysem of Linear Differenial Equaions In "Ordinary Differenial Equaions" we've learned how o solve a differenial equaion for a variable, such as: y'k5$e K2$x =0 solve DE yx = K 5 2 ek2 x C_C1 2$y''C7$y

More information

Lecture Notes 2. The Hilbert Space Approach to Time Series

Lecture Notes 2. The Hilbert Space Approach to Time Series Time Series Seven N. Durlauf Universiy of Wisconsin. Basic ideas Lecure Noes. The Hilber Space Approach o Time Series The Hilber space framework provides a very powerful language for discussing he relaionship

More information

Solution: b All the terms must have the dimension of acceleration. We see that, indeed, each term has the units of acceleration

Solution: b All the terms must have the dimension of acceleration. We see that, indeed, each term has the units of acceleration PHYS 54 Tes Pracice Soluions Spring 8 Q: [4] Knowing ha in he ne epression a is acceleraion, v is speed, is posiion and is ime, from a dimensional v poin of view, he equaion a is a) incorrec b) correc

More information

Mon Apr 9 EP 7.6 Convolutions and Laplace transforms. Announcements: Warm-up Exercise:

Mon Apr 9 EP 7.6 Convolutions and Laplace transforms. Announcements: Warm-up Exercise: Mah 225-4 Week 3 April 9-3 EP 7.6 - convoluions; 6.-6.2 - eigenvalues, eigenvecors and diagonalizabiliy; 7. - sysems of differenial equaions. Mon Apr 9 EP 7.6 Convoluions and Laplace ransforms. Announcemens:

More information

Echocardiography Project and Finite Fourier Series

Echocardiography Project and Finite Fourier Series Echocardiography Projec and Finie Fourier Series 1 U M An echocardiagram is a plo of how a porion of he hear moves as he funcion of ime over he one or more hearbea cycles If he hearbea repeas iself every

More information

Differential Equations

Differential Equations Mah 21 (Fall 29) Differenial Equaions Soluion #3 1. Find he paricular soluion of he following differenial equaion by variaion of parameer (a) y + y = csc (b) 2 y + y y = ln, > Soluion: (a) The corresponding

More information

15. Vector Valued Functions

15. Vector Valued Functions 1. Vecor Valued Funcions Up o his poin, we have presened vecors wih consan componens, for example, 1, and,,4. However, we can allow he componens of a vecor o be funcions of a common variable. For example,

More information

Challenge Problems. DIS 203 and 210. March 6, (e 2) k. k(k + 2). k=1. f(x) = k(k + 2) = 1 x k

Challenge Problems. DIS 203 and 210. March 6, (e 2) k. k(k + 2). k=1. f(x) = k(k + 2) = 1 x k Challenge Problems DIS 03 and 0 March 6, 05 Choose one of he following problems, and work on i in your group. Your goal is o convince me ha your answer is correc. Even if your answer isn compleely correc,

More information

Solution of Integro-Differential Equations by Using ELzaki Transform

Solution of Integro-Differential Equations by Using ELzaki Transform Global Journal of Mahemaical Sciences: Theory and Pracical. Volume, Number (), pp. - Inernaional Research Publicaion House hp://www.irphouse.com Soluion of Inegro-Differenial Equaions by Using ELzaki Transform

More information

where the coordinate X (t) describes the system motion. X has its origin at the system static equilibrium position (SEP).

where the coordinate X (t) describes the system motion. X has its origin at the system static equilibrium position (SEP). Appendix A: Conservaion of Mechanical Energy = Conservaion of Linear Momenum Consider he moion of a nd order mechanical sysem comprised of he fundamenal mechanical elemens: ineria or mass (M), siffness

More information

Physics 1402: Lecture 22 Today s Agenda

Physics 1402: Lecture 22 Today s Agenda Physics 142: ecure 22 Today s Agenda Announcemens: R - RV - R circuis Homework 6: due nex Wednesday Inducion / A curren Inducion Self-Inducance, R ircuis X X X X X X X X X long solenoid Energy and energy

More information

2. Nonlinear Conservation Law Equations

2. Nonlinear Conservation Law Equations . Nonlinear Conservaion Law Equaions One of he clear lessons learned over recen years in sudying nonlinear parial differenial equaions is ha i is generally no wise o ry o aack a general class of nonlinear

More information

20. Applications of the Genetic-Drift Model

20. Applications of the Genetic-Drift Model 0. Applicaions of he Geneic-Drif Model 1) Deermining he probabiliy of forming any paricular combinaion of genoypes in he nex generaion: Example: If he parenal allele frequencies are p 0 = 0.35 and q 0

More information

Lecture 4 Kinetics of a particle Part 3: Impulse and Momentum

Lecture 4 Kinetics of a particle Part 3: Impulse and Momentum MEE Engineering Mechanics II Lecure 4 Lecure 4 Kineics of a paricle Par 3: Impulse and Momenum Linear impulse and momenum Saring from he equaion of moion for a paricle of mass m which is subjeced o an

More information

ψ(t) = V x (0)V x (t)

ψ(t) = V x (0)V x (t) .93 Home Work Se No. (Professor Sow-Hsin Chen Spring Term 5. Due March 7, 5. This problem concerns calculaions of analyical expressions for he self-inermediae scaering funcion (ISF of he es paricle in

More information

Structural Dynamics and Earthquake Engineering

Structural Dynamics and Earthquake Engineering Srucural Dynamics and Earhquae Engineering Course 1 Inroducion. Single degree of freedom sysems: Equaions of moion, problem saemen, soluion mehods. Course noes are available for download a hp://www.c.up.ro/users/aurelsraan/

More information

SOLUTIONS TO ECE 3084

SOLUTIONS TO ECE 3084 SOLUTIONS TO ECE 384 PROBLEM 2.. For each sysem below, specify wheher or no i is: (i) memoryless; (ii) causal; (iii) inverible; (iv) linear; (v) ime invarian; Explain your reasoning. If he propery is no

More information

Electrical Circuits. 1. Circuit Laws. Tools Used in Lab 13 Series Circuits Damped Vibrations: Energy Van der Pol Circuit

Electrical Circuits. 1. Circuit Laws. Tools Used in Lab 13 Series Circuits Damped Vibrations: Energy Van der Pol Circuit V() R L C 513 Elecrical Circuis Tools Used in Lab 13 Series Circuis Damped Vibraions: Energy Van der Pol Circui A series circui wih an inducor, resisor, and capacior can be represened by Lq + Rq + 1, a

More information

Laplace Transform and its Relation to Fourier Transform

Laplace Transform and its Relation to Fourier Transform Chaper 6 Laplace Transform and is Relaion o Fourier Transform (A Brief Summary) Gis of he Maer 2 Domains of Represenaion Represenaion of signals and sysems Time Domain Coninuous Discree Time Time () [n]

More information

EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE

EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE Version April 30, 2004.Submied o CTU Repors. EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE Per Krysl Universiy of California, San Diego La Jolla, California 92093-0085,

More information

Hamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation:

Hamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation: M ah 5 7 Fall 9 L ecure O c. 4, 9 ) Hamilon- J acobi Equaion: Weak S oluion We coninue he sudy of he Hamilon-Jacobi equaion: We have shown ha u + H D u) = R n, ) ; u = g R n { = }. ). In general we canno

More information

14 Autoregressive Moving Average Models

14 Autoregressive Moving Average Models 14 Auoregressive Moving Average Models In his chaper an imporan parameric family of saionary ime series is inroduced, he family of he auoregressive moving average, or ARMA, processes. For a large class

More information

Theory of! Partial Differential Equations-I!

Theory of! Partial Differential Equations-I! hp://users.wpi.edu/~grear/me61.hml! Ouline! Theory o! Parial Dierenial Equaions-I! Gréar Tryggvason! Spring 010! Basic Properies o PDE!! Quasi-linear Firs Order Equaions! - Characerisics! - Linear and

More information

A Special Hour with Relativity

A Special Hour with Relativity A Special Hour wih Relaiviy Kenneh Chu The Graduae Colloquium Deparmen of Mahemaics Universiy of Uah Oc 29, 2002 Absrac Wha promped Einsen: Incompaibiliies beween Newonian Mechanics and Maxwell s Elecromagneism.

More information

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n Module Fick s laws of diffusion Fick s laws of diffusion and hin film soluion Adolf Fick (1855) proposed: d J α d d d J (mole/m s) flu (m /s) diffusion coefficien and (mole/m 3 ) concenraion of ions, aoms

More information

A Bayesian Approach to Spectral Analysis

A Bayesian Approach to Spectral Analysis Chirped Signals A Bayesian Approach o Specral Analysis Chirped signals are oscillaing signals wih ime variable frequencies, usually wih a linear variaion of frequency wih ime. E.g. f() = A cos(ω + α 2

More information

Reading from Young & Freedman: For this topic, read sections 25.4 & 25.5, the introduction to chapter 26 and sections 26.1 to 26.2 & 26.4.

Reading from Young & Freedman: For this topic, read sections 25.4 & 25.5, the introduction to chapter 26 and sections 26.1 to 26.2 & 26.4. PHY1 Elecriciy Topic 7 (Lecures 1 & 11) Elecric Circuis n his opic, we will cover: 1) Elecromoive Force (EMF) ) Series and parallel resisor combinaions 3) Kirchhoff s rules for circuis 4) Time dependence

More information

Chapter 12: Velocity, acceleration, and forces

Chapter 12: Velocity, acceleration, and forces To Feel a Force Chaper Spring, Chaper : A. Saes of moion For moion on or near he surface of he earh, i is naural o measure moion wih respec o objecs fixed o he earh. The 4 hr. roaion of he earh has a measurable

More information

Effects of Coordinate Curvature on Integration

Effects of Coordinate Curvature on Integration Effecs of Coordinae Curvaure on Inegraion Chrisopher A. Lafore clafore@gmail.com Absrac In his paper, he inegraion of a funcion over a curved manifold is examined in he case where he curvaure of he manifold

More information

Concourse Math Spring 2012 Worked Examples: Matrix Methods for Solving Systems of 1st Order Linear Differential Equations

Concourse Math Spring 2012 Worked Examples: Matrix Methods for Solving Systems of 1st Order Linear Differential Equations Concourse Mah 80 Spring 0 Worked Examples: Marix Mehods for Solving Sysems of s Order Linear Differenial Equaions The Main Idea: Given a sysem of s order linear differenial equaions d x d Ax wih iniial

More information

A New Perturbative Approach in Nonlinear Singularity Analysis

A New Perturbative Approach in Nonlinear Singularity Analysis Journal of Mahemaics and Saisics 7 (: 49-54, ISSN 549-644 Science Publicaions A New Perurbaive Approach in Nonlinear Singulariy Analysis Ta-Leung Yee Deparmen of Mahemaics and Informaion Technology, The

More information

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x WEEK-3 Reciaion PHYS 131 Ch. 3: FOC 1, 3, 4, 6, 14. Problems 9, 37, 41 & 71 and Ch. 4: FOC 1, 3, 5, 8. Problems 3, 5 & 16. Feb 8, 018 Ch. 3: FOC 1, 3, 4, 6, 14. 1. (a) The horizonal componen of he projecile

More information

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3 and d = c b - b c c d = c b - b c c This process is coninued unil he nh row has been compleed. The complee array of coefficiens is riangular. Noe ha in developing he array an enire row may be divided or

More information

k B 2 Radiofrequency pulses and hardware

k B 2 Radiofrequency pulses and hardware 1 Exra MR Problems DC Medical Imaging course April, 214 he problems below are harder, more ime-consuming, and inended for hose wih a more mahemaical background. hey are enirely opional, bu hopefully will

More information

Robotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1.

Robotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1. Roboics I April 11, 017 Exercise 1 he kinemaics of a 3R spaial robo is specified by he Denavi-Harenberg parameers in ab 1 i α i d i a i θ i 1 π/ L 1 0 1 0 0 L 3 0 0 L 3 3 able 1: able of DH parameers of

More information

( ) a system of differential equations with continuous parametrization ( T = R + These look like, respectively:

( ) a system of differential equations with continuous parametrization ( T = R + These look like, respectively: XIII. DIFFERENCE AND DIFFERENTIAL EQUATIONS Ofen funcions, or a sysem of funcion, are paramerized in erms of some variable, usually denoed as and inerpreed as ime. The variable is wrien as a funcion of

More information

Finish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details!

Finish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details! MAT 257, Handou 6: Ocober 7-2, 20. I. Assignmen. Finish reading Chaper 2 of Spiva, rereading earlier secions as necessary. handou and fill in some missing deails! II. Higher derivaives. Also, read his

More information

( ) = Q 0. ( ) R = R dq. ( t) = I t

( ) = Q 0. ( ) R = R dq. ( t) = I t ircuis onceps The addiion of a simple capacior o a circui of resisors allows wo relaed phenomena o occur The observaion ha he ime-dependence of a complex waveform is alered by he circui is referred o as

More information

Notes on Kalman Filtering

Notes on Kalman Filtering Noes on Kalman Filering Brian Borchers and Rick Aser November 7, Inroducion Daa Assimilaion is he problem of merging model predicions wih acual measuremens of a sysem o produce an opimal esimae of he curren

More information

6.2 Transforms of Derivatives and Integrals.

6.2 Transforms of Derivatives and Integrals. SEC. 6.2 Transforms of Derivaives and Inegrals. ODEs 2 3 33 39 23. Change of scale. If l( f ()) F(s) and c is any 33 45 APPLICATION OF s-shifting posiive consan, show ha l( f (c)) F(s>c)>c (Hin: In Probs.

More information

Solutions from Chapter 9.1 and 9.2

Solutions from Chapter 9.1 and 9.2 Soluions from Chaper 9 and 92 Secion 9 Problem # This basically boils down o an exercise in he chain rule from calculus We are looking for soluions of he form: u( x) = f( k x c) where k x R 3 and k is

More information

Announcements: Warm-up Exercise:

Announcements: Warm-up Exercise: Fri Apr 13 7.1 Sysems of differenial equaions - o model muli-componen sysems via comparmenal analysis hp//en.wikipedia.org/wiki/muli-comparmen_model Announcemens Warm-up Exercise Here's a relaively simple

More information

- If one knows that a magnetic field has a symmetry, one may calculate the magnitude of B by use of Ampere s law: The integral of scalar product

- If one knows that a magnetic field has a symmetry, one may calculate the magnitude of B by use of Ampere s law: The integral of scalar product 11.1 APPCATON OF AMPEE S AW N SYMMETC MAGNETC FEDS - f one knows ha a magneic field has a symmery, one may calculae he magniude of by use of Ampere s law: The inegral of scalar produc Closed _ pah * d

More information

In this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should

In this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should Cambridge Universiy Press 978--36-60033-7 Cambridge Inernaional AS and A Level Mahemaics: Mechanics Coursebook Excerp More Informaion Chaper The moion of projeciles In his chaper he model of free moion

More information

Institute for Mathematical Methods in Economics. University of Technology Vienna. Singapore, May Manfred Deistler

Institute for Mathematical Methods in Economics. University of Technology Vienna. Singapore, May Manfred Deistler MULTIVARIATE TIME SERIES ANALYSIS AND FORECASTING Manfred Deisler E O S Economerics and Sysems Theory Insiue for Mahemaical Mehods in Economics Universiy of Technology Vienna Singapore, May 2004 Inroducion

More information

Vanishing Viscosity Method. There are another instructive and perhaps more natural discontinuous solutions of the conservation law

Vanishing Viscosity Method. There are another instructive and perhaps more natural discontinuous solutions of the conservation law Vanishing Viscosiy Mehod. There are anoher insrucive and perhaps more naural disconinuous soluions of he conservaion law (1 u +(q(u x 0, he so called vanishing viscosiy mehod. This mehod consiss in viewing

More information

From Particles to Rigid Bodies

From Particles to Rigid Bodies Rigid Body Dynamics From Paricles o Rigid Bodies Paricles No roaions Linear velociy v only Rigid bodies Body roaions Linear velociy v Angular velociy ω Rigid Bodies Rigid bodies have boh a posiion and

More information

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities: Mah 4 Eam Review Problems Problem. Calculae he 3rd Taylor polynomial for arcsin a =. Soluion. Le f() = arcsin. For his problem, we use he formula f() + f () + f ()! + f () 3! for he 3rd Taylor polynomial

More information

= ( ) ) or a system of differential equations with continuous parametrization (T = R

= ( ) ) or a system of differential equations with continuous parametrization (T = R XIII. DIFFERENCE AND DIFFERENTIAL EQUATIONS Ofen funcions, or a sysem of funcion, are paramerized in erms of some variable, usually denoed as and inerpreed as ime. The variable is wrien as a funcion of

More information

Nature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time.

Nature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time. Supplemenary Figure 1 Spike-coun auocorrelaions in ime. Normalized auocorrelaion marices are shown for each area in a daase. The marix shows he mean correlaion of he spike coun in each ime bin wih he spike

More information

Lecture Outline. Introduction Transmission Line Equations Transmission Line Wave Equations 8/10/2018. EE 4347 Applied Electromagnetics.

Lecture Outline. Introduction Transmission Line Equations Transmission Line Wave Equations 8/10/2018. EE 4347 Applied Electromagnetics. 8/10/018 Course Insrucor Dr. Raymond C. Rumpf Office: A 337 Phone: (915) 747 6958 E Mail: rcrumpf@uep.edu EE 4347 Applied Elecromagneics Topic 4a Transmission Line Equaions Transmission These Line noes

More information

Wall. x(t) f(t) x(t = 0) = x 0, t=0. which describes the motion of the mass in absence of any external forcing.

Wall. x(t) f(t) x(t = 0) = x 0, t=0. which describes the motion of the mass in absence of any external forcing. MECHANICS APPLICATIONS OF SECOND-ORDER ODES 7 Mechanics applicaions of second-order ODEs Second-order linear ODEs wih consan coefficiens arise in many physical applicaions. One physical sysems whose behaviour

More information

Problem Set 5. Graduate Macro II, Spring 2017 The University of Notre Dame Professor Sims

Problem Set 5. Graduate Macro II, Spring 2017 The University of Notre Dame Professor Sims Problem Se 5 Graduae Macro II, Spring 2017 The Universiy of Nore Dame Professor Sims Insrucions: You may consul wih oher members of he class, bu please make sure o urn in your own work. Where applicable,

More information

Ordinary Differential Equations

Ordinary Differential Equations Lecure 22 Ordinary Differenial Equaions Course Coordinaor: Dr. Suresh A. Karha, Associae Professor, Deparmen of Civil Engineering, IIT Guwahai. In naure, mos of he phenomena ha can be mahemaically described

More information

!!"#"$%&#'()!"#&'(*%)+,&',-)./0)1-*23)

!!#$%&#'()!#&'(*%)+,&',-)./0)1-*23) "#"$%&#'()"#&'(*%)+,&',-)./)1-*) #$%&'()*+,&',-.%,/)*+,-&1*#$)()5*6$+$%*,7&*-'-&1*(,-&*6&,7.$%$+*&%'(*8$&',-,%'-&1*(,-&*6&,79*(&,%: ;..,*&1$&$.$%&'()*1$$.,'&',-9*(&,%)?%*,('&5

More information

Module II, Part C. More Insight into Fiber Dispersion

Module II, Part C. More Insight into Fiber Dispersion Moule II Par C More Insigh ino Fiber Dispersion . Polariaion Moe Dispersion Fiber Birefringence: Imperfec cylinrical symmery leas o wha is known as birefringence. Recall he HE moe an is E x componen which

More information

Let us start with a two dimensional case. We consider a vector ( x,

Let us start with a two dimensional case. We consider a vector ( x, Roaion marices We consider now roaion marices in wo and hree dimensions. We sar wih wo dimensions since wo dimensions are easier han hree o undersand, and one dimension is a lile oo simple. However, our

More information

IB Physics Kinematics Worksheet

IB Physics Kinematics Worksheet IB Physics Kinemaics Workshee Wrie full soluions and noes for muliple choice answers. Do no use a calculaor for muliple choice answers. 1. Which of he following is a correc definiion of average acceleraion?

More information

Section 3.8, Mechanical and Electrical Vibrations

Section 3.8, Mechanical and Electrical Vibrations Secion 3.8, Mechanical and Elecrical Vibraions Mechanical Unis in he U.S. Cusomary and Meric Sysems Disance Mass Time Force g (Earh) Uni U.S. Cusomary MKS Sysem CGS Sysem fee f slugs seconds sec pounds

More information

2.3 SCHRÖDINGER AND HEISENBERG REPRESENTATIONS

2.3 SCHRÖDINGER AND HEISENBERG REPRESENTATIONS Andrei Tokmakoff, MIT Deparmen of Chemisry, 2/22/2007 2-17 2.3 SCHRÖDINGER AND HEISENBERG REPRESENTATIONS The mahemaical formulaion of he dynamics of a quanum sysem is no unique. So far we have described

More information