Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum

Size: px
Start display at page:

Download "Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum"

Transcription

1 October 6, 4 ARDB Note Analytic Scaling Formulas for Crosse Laser Acceleration in Vacuum Robert J. Noble Stanfor Linear Accelerator Center, Stanfor University 575 San Hill Roa, Menlo Park, California 945 USA Abstract. Using an analytic result of Sprangle et al for the longituinal accelerating fiel experience by a relativistic electron from two crosse Gaussian lasers in vacuum, we erive a new approximation for the accelerating potential vali for any energy provie the electron is traveling near the spee of light. The expression is applicable when the laser crossing angle is small an the electron s istance from the focal point is small compare to the Rayleigh length of the focuse lasers. Both these conitions are the usual ones for high energy particle acceleration in a boune accelerating structure near a laser focus. We use our expression for the accelerating potential to obtain formulas for the maximum energy gain per slip length, the optimal crossing angle, an the corresponing optimal slip istance. In the high energy limit, the optimal crossing angle is times the iffraction angle an the graient is.65 times the prouct of the iniviual laser fiel strength an the laser iffraction angle. Near the optimal crossing angle, the graient increases approximately linearly with the ratio of crossing angle to iffraction angle.. Introuction The purpose of this note is to iscuss the optimiation of laser acceleration of electrons in vacuum with respect to key parameters etermine in experiments such as LEAP at Stanfor University an E-63 at SLAC []. In particular, we present new analytic scaling formulas for crosse laser acceleration vali for any energy provie the electron travels at nearly the spee of light. The prototypical vacuum laser acceleration configuration we consier was analye by Sprangle et al [,3] an is shown in Figure. In this scheme a pair of linearly polarie laser beams with Gaussian transverse profiles, having the same frequency an equal strength E, are focuse an intersecte in vacuum. The first laser propagates along the axis an the secon laser propagates along the axis, where the axis an the axis are rotate by the crossing angles of an, respectively, with respect to the axis. The phases of the lasers are such that the transverse electric fiels cancel along the axis while the axial fiels a. Properly phase electrons injecte along the axis can be accelerate by the net axial component of the laser fiel. Figure. Coorinate system an electric fiels for two intersecting laser beams.

2 The analysis by Sprangle et al was one assuming that the focuse Gaussian laser propagates as a nearly uniirectional wave of finite cross section escribe by a fiel E t = E(x,y, exp(ik-iωt, an E t >> E. Uner this assumption the usual vector wave equation can be replace by the scalar Helmholt equation for the transverse fiel E(x,y,, an the longituinal fiel is etermine later by emaning iv E =. The Helmholt equation is solve in the so-calle paraxial approximation in which the fiel E(x,y, is assume to vary slowly along the irection on the scale of a wavelength. The paraxial approximation is vali when the laser spot sie (waist w is much greater than the wavelength λ. For small crossing angles <<, the axial accelerating fiel, as seen by an electron traveling with velocity v c along the axis, was calculate to be E ˆ ( / E = exp cosψ, t (a + ˆ ˆ + where ˆ = / R, R = π / λ is the Rayleigh length, = w / is the laser iffraction angle, w ψ ˆ ˆ + = tan ˆ φ ( γ t, (b + ˆ t = /v has been assume, v / ( γ, = ( / c / γ v, an φ = R φ is the arrival phase of the electron relative to the laser maximum fiel value (. In Fig. the normalie electric fiel E /E as a function of the position / R is plotte for the parameters of the E-63 experiment [] using nominal 6 MeV kinetic energy electrons an for the same parameters in the limit of high energy electrons (γ. E /E =.5 mra = 6.4 mra / R Figure. Normalie electric fiel for the E-63 experiment using nominal 6 MeV electrons (upper curve; re an in the high energy limit (lower curve; blue.

3 The phase velocity of the accelerating fiel is greater than the spee of light, an therefore the fiel slips ahea of the electron. An electron interacting with the laser fiels over the entire region = - to + will experience no net energy gain. But if the laser fiels are terminate, say by bounaries, then net energy gain can occur. An electron injecte π/ in phase relative the fiel crest will gain energy until the phase slip changes by π, at which point it begins to be ecelerate. From Eq. (b near the focal point (<< R, the istance require for the electron to slip by π is λγ + γ / γ, ( s c /( c / ( + where γ efines a critical energy. Whenγ < γ, slippage is ominate by the low c = velocity of the electron, an whenγ > γ, slippage is ominate by the phase shift ue to laser iffraction. c Experimentally s etermines the length of the accelerator cell, which is marke by bounaries to terminate the laser fiels. In Fig. these points correspon to the locations of the fiel ero crossings on either sie of the origin. The axial electric fiel can in principle be written as the graient of an effective potential, E = U / ˆ, provie one can ientify the potential from the integral U = E ζ. R As far as the author knows, the general integral of Eq. ( is not known in close form. However in the special limit where γ tan ˆ = tan ˆ / ˆ to c ˆ R, one can use the trigonometric ientity ( ( rewrite cosψ t as a prouct of sine an cosine functions an absorb the leaing ( ˆ + foun in E. This observation allowe Sprangle et al to ientify the potential in the high energy limit as 4E ˆ ˆ U ( ˆ = exp sin φ sinφ, (3 + + ˆ ˆ where k = π/λ, an the potential is efine to be ero at the origin. Note that high energy here actually means γ >>. The normalie potential U/(E /k is shown in Fig.3 along with the numerical integral of the electric fiel ( for the parameters of E-63. U/(E /k =.5 mra = 6.4 mra Figure 3. Normalie potential from the numerical integration of Eq. ( for 6 MeV electrons (smaller curve; re an the potential in the high energy limit from Eq. (3 (larger curve; blue. / R 3

4 . Effective Accelerating Potential for Arbitrary Energy Although useful as a limit, Eq. (3 oes not allow one to unerstan the optimiation of parameters for a real experiment at finite energy or to erive simple scaling results for arbitrary energy. The analytic integration of Eq. ( to obtain the potential for arbitrary γ is hampere by the fact that the phase ψ t is nonlinear in ẑ. To make progress we note that one is only intereste in E an the potential U over a slip length aroun the origin, i.e. in the range [- s /, + s /]. Diviing the slip length s in Eq. ( by R one fins ˆ s π / =, (4 + + γ which is always less than. This suggests the following approximation for E vali for a relativistic electron when ˆ <<, which is essentially all of the accelerating region, ( ˆ E E exp ˆ cos ˆ ˆ + φ. (5 γ Here we have replace + ˆ by an linearie the arctangent function in Eq. ( within this approximation. In Fig. 4 the normalie electric fiel (5 for the E-63 parameters is plotte along with the fiels shown earlier in Fig.. The approximation is in excellent agreement with Eq. ( in the central region of π phase slip where the accelerator cell is locate. E /E =.5 mra = 6.4 mra Figure 4. Normalie electric fiel from the approximation in Eq. (5 for the E-63 parameters (ashe; brown compare with the exact fiel from Eq. ( (mile curve; re an high energy limit (lower; blue. The integral of E in Eq. (5 over the range [, ˆ] can be expresse in terms of the real part of the error function erf(w with complex argument w = u+iv [4]. The approximate effective potential is then / R 4

5 π E i ( = φ U ˆ exp( v R[ e ( erf ( u + iv erf ( iv], (6a where the quantities u an v are efine by u( ˆ = ˆ an v = + +. (6b γ Note that u ( ẑ is proportional to the longituinal coorinate, an v is proportional to the phase shift per unit length in Eq. (5. From Eq. (4, we can write v = / π / ˆ. ( s The error function is convenient to work with because so many approximations exist for it epening on the magnitues of u, v, an w. We anticipate that the optimal value of for the maximum energy gain will be of orer because if the crossing angle is much greater than the iffraction angle the slip istance is too short for the electron to gain much energy, while if the crossing angle is too small the projecte longituinal electric fiel will not prouce much energy gain. The analytic results in the next section will be consistent with this assumption. The error function in Eq. (6 can be greatly simplifie uner the conitions that u << an w >>. The first conition is always true within the slip region since ˆ << an is of orer, an the secon conition is a statement that v >>, which is vali in practice since v 3/ for all. Performing an asymptotic approximation of the error function uner these conitions yiels E ( ( u U ˆ exp u + sin( uv + φ + cos( uv + φ + sinφ, (7 v v v v where terms up to orer /v 3 an to first orer in u/v have been retaine. The /v 3 terms are essential for a goo approximation. The error function potential (6 an the asymptotic approximation (7 are plotte in Fig. (5 along with the two potentials shown earlier in Fig. 3. These potentials are excellent approximations to the numerical integral in the central region of π phase slip where the accelerator cell is locate. U/(E /k =.5 mra = 6.4 mra Figure 5. Approximate potentials from Eqs. (6 (ashe; brown an (7 (ash-ot; magenta compare with the potentials plotte in Fig. 3. / R 5

6 3. Analytic Scaling Formulas Because there are only exponential an trigonometric functions in the potential of Eq. (7, it is very useful for eriving simple scaling formulas for crosse laser acceleration at arbitrary energy. The energy gain per unit charge is efine by W = U ˆ U ˆ. (8 ( ( ( final initial For energy gain over a slip istance, we have ˆ = ˆ /, an ˆ = ˆ /. At these points, uv = ± π /. Using Eq. (7, the energy gain per unit charge in a slip istance is 8 E 3 W = cosφ ( exp( r + r π r / 4, (9a where / r =. (9b + + γ The optimal crossing angle which maximies the energy gain is etermine by the conition W / ( / =, or using the chain rule ( W / r ( r / ( / =. Solving ( W / r = yiels one real root for r, but when substitute into Eq. (9b, the corresponing solution for / is complex. The other partial erivative yiels ( / r / ( / = =, ( + ( / + ( / γ ( + ( / + ( / γ which yiels the optimal crossing angle to iffraction angle ratio ( / = +. ( opt γ The maximum energy gain per unit charge over a slip length is 4E π Wmax cosφ + exp 3 / opt / 6 / opt initial = ( ( ( opt s final s. ( In the limit γ, the optimum crossing angle becomes, an the maximum energy gain is 5 E cosφ E cosφ Wmax ( γ = exp( π / 3.6. (3 The optimal slip length (= accelerator cell length corresponing to the optimal crossing angle is λ λ 4 s ( opt = =. (4 + ( / opt γ The optimal slip length ivie by the Rayleigh length ( π R = λ is π / 4 ˆ s ( opt =, (5 + γ which is less than one for all γ. 6

7 From Eqs. ( an (4, one obtains the average graient when the energy gain is maximie, 4E π G = Wmax / s( opt = cosφ ( / + exp opt π ( / ( opt 6 / opt. (6 In the limit γ, this graient becomes 5 Ginf = exp( π / 3 E cosφ.65e cosφ. (7 π It shoul be stresse that G in Eq. (6 is not the maximum graient that can be achieve in the accelerator cell, but rather is the particular value when the energy gain is maximie. In fact W tens to be nearly constant for / near the optimal value, while the slip length s tens to fall off rapily with increasing /. Hence one can increase the graient by increasing / somewhat above the optimum an still achieve about the same energy gain in a shorter accelerator cell. For a given laser strength E, which may be limite by surface amage threshols, the accelerating graient can be increase in this way. One fins that the graient increases approximately linearly with / within the approximations of our analysis. In Fig. 6 the normalie energy gain from Eq. (9, the slip length an graient are plotte as a function of / to illustrate this relation for the E-63 parameters. We have not plotte the exact energy gain calculate from the numerical integration of Eq. ( for comparison in Fig. 6 since the curve woul lie irectly over the approximation. Normalie graient s / R W/(-4Ecosφ/ / Figure 6. Normalie energy gain W/(-4E cosφ /, slip length s / R, an normalie graient (ratio of these quantities as a function of the laser crossing angle ratio / for the E-63 parameters. The optimal value of / for maximal energy gain in a cell is.94, but larger graients are possible by increasing the crossing angle with little reuction in the energy gain. In Table we summarie the E-63 experimental parameters from the proposal [] an the optimal values for the key parameters calculate from our analysis above. It shoul be note that the esign crossing angle of.5 mra originally chosen for E-63 was etermine by using the Sprangle high energy limit for the potential (3 an the slip length ( for a 6 MeV electron, an then varying the crossing angle to achieve a maximum energy gain. This was the only metho available without an analytic formula vali for arbitrary energy. The preicte energy gain using the previous metho overestimates the energy gain by % 7

8 accoring to our analysis. This can be seen graphically in Fig. 5 where the potential in the high energy limit is significantly larger than the finite energy potential near the turning points. Table : Parameters for the E-63 laser acceleration experiment. 4. Conclusion Parameter Value Laser E 5.9x 9 V/m λ.8 µm w 4 µm R Design 6.8 mm 6.4 mra.5 mra γ c 68.5 Electron KE 6 MeV γ 8.4 -ew (using Eq.(3 9 kev (/ opt.94 opt.5 mra -ew max 34 kev s(opt.63 mm W max / s(opt 89 MV/m 6 MV/m G inf Using an analytic result of Sprangle et al for the longituinal accelerating fiel experience by a relativistic electron from two crosse Gaussian lasers in vacuum, we have erive a new approximation for the accelerating potential vali for any energy provie the electron is traveling near the spee of light. The optimal laser crossing angle for maximum energy gain in a slip length, an the corresponing slip length were etermine. It was shown that in the high energy limit, the optimal crossing angle is times the iffraction angle, an the graient is approximately.65 times the prouct of the iniviual laser fiel strength an the laser iffraction angle. The energy gain function has a broa maximum aroun the optimal crossing angle, while the slip length falls off rapily with angle. As a result in the neighborhoo of the optimal value, the graient increases almost linearly with crossing angle. The calculations escribe in this paper have been written into a Mathca program, available from the author, which allows the user to vary the key parameters an plot the corresponing accelerating fiels an potentials. ACKNOWLEDGMENTS This work was supporte by the U.S. Department of Energy uner contract DE-AC-76SF55. REFERENCES. E. Colby et al, Laser Acceleration of Electrons in Vacuum, SLAC Proposal E-63, August 4,.. E. Esarey, P. Sprangle, an J. Krall, Laser acceleration of electrons in vacuum, Phys. Rev. E, Vol. 5, 5443 ( P. Sprangle, E. Esarey, J. Krall, an A. Ting, Vacuum laser acceleration, Optics Comm., Vol. 4, 69 ( M. Abramowit an I.A. Stegun, Hanbook of Mathematical Functions, p. 95, Ninth eition, November 97. 8

Table of Common Derivatives By David Abraham

Table of Common Derivatives By David Abraham Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec

More information

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+

More information

Math 342 Partial Differential Equations «Viktor Grigoryan

Math 342 Partial Differential Equations «Viktor Grigoryan Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite

More information

Conservation Laws. Chapter Conservation of Energy

Conservation Laws. Chapter Conservation of Energy 20 Chapter 3 Conservation Laws In orer to check the physical consistency of the above set of equations governing Maxwell-Lorentz electroynamics [(2.10) an (2.12) or (1.65) an (1.68)], we examine the action

More information

3.7 Implicit Differentiation -- A Brief Introduction -- Student Notes

3.7 Implicit Differentiation -- A Brief Introduction -- Student Notes Fin these erivatives of these functions: y.7 Implicit Differentiation -- A Brief Introuction -- Stuent Notes tan y sin tan = sin y e = e = Write the inverses of these functions: y tan y sin How woul we

More information

Physics 505 Electricity and Magnetism Fall 2003 Prof. G. Raithel. Problem Set 3. 2 (x x ) 2 + (y y ) 2 + (z + z ) 2

Physics 505 Electricity and Magnetism Fall 2003 Prof. G. Raithel. Problem Set 3. 2 (x x ) 2 + (y y ) 2 + (z + z ) 2 Physics 505 Electricity an Magnetism Fall 003 Prof. G. Raithel Problem Set 3 Problem.7 5 Points a): Green s function: Using cartesian coorinates x = (x, y, z), it is G(x, x ) = 1 (x x ) + (y y ) + (z z

More information

Implicit Differentiation

Implicit Differentiation Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,

More information

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University

More information

Linear First-Order Equations

Linear First-Order Equations 5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)

More information

Angles-Only Orbit Determination Copyright 2006 Michel Santos Page 1

Angles-Only Orbit Determination Copyright 2006 Michel Santos Page 1 Angles-Only Orbit Determination Copyright 6 Michel Santos Page 1 Abstract This ocument presents a re-erivation of the Gauss an Laplace Angles-Only Methos for Initial Orbit Determination. It keeps close

More information

12.5. Differentiation of vectors. Introduction. Prerequisites. Learning Outcomes

12.5. Differentiation of vectors. Introduction. Prerequisites. Learning Outcomes Differentiation of vectors 12.5 Introuction The area known as vector calculus is use to moel mathematically a vast range of engineering phenomena incluing electrostatics, electromagnetic fiels, air flow

More information

18 EVEN MORE CALCULUS

18 EVEN MORE CALCULUS 8 EVEN MORE CALCULUS Chapter 8 Even More Calculus Objectives After stuing this chapter you shoul be able to ifferentiate an integrate basic trigonometric functions; unerstan how to calculate rates of change;

More information

Lecture 1b. Differential operators and orthogonal coordinates. Partial derivatives. Divergence and divergence theorem. Gradient. A y. + A y y dy. 1b.

Lecture 1b. Differential operators and orthogonal coordinates. Partial derivatives. Divergence and divergence theorem. Gradient. A y. + A y y dy. 1b. b. Partial erivatives Lecture b Differential operators an orthogonal coorinates Recall from our calculus courses that the erivative of a function can be efine as f ()=lim 0 or using the central ifference

More information

Qubit channels that achieve capacity with two states

Qubit channels that achieve capacity with two states Qubit channels that achieve capacity with two states Dominic W. Berry Department of Physics, The University of Queenslan, Brisbane, Queenslan 4072, Australia Receive 22 December 2004; publishe 22 March

More information

Chapter 2 Derivatives

Chapter 2 Derivatives Chapter Derivatives Section. An Intuitive Introuction to Derivatives Consier a function: Slope function: Derivative, f ' For each, the slope of f is the height of f ' Where f has a horizontal tangent line,

More information

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x) Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)

More information

Tutorial 1 Differentiation

Tutorial 1 Differentiation Tutorial 1 Differentiation What is Calculus? Calculus 微積分 Differential calculus Differentiation 微分 y lim 0 f f The relation of very small changes of ifferent quantities f f y y Integral calculus Integration

More information

The Principle of Least Action

The Principle of Least Action Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of

More information

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1 Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition

More information

Chapter 2 Lagrangian Modeling

Chapter 2 Lagrangian Modeling Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie

More information

Lecture 2 Lagrangian formulation of classical mechanics Mechanics

Lecture 2 Lagrangian formulation of classical mechanics Mechanics Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,

More information

Differentiation ( , 9.5)

Differentiation ( , 9.5) Chapter 2 Differentiation (8.1 8.3, 9.5) 2.1 Rate of Change (8.2.1 5) Recall that the equation of a straight line can be written as y = mx + c, where m is the slope or graient of the line, an c is the

More information

Math 1272 Solutions for Spring 2005 Final Exam. asked to find the limit of the sequence. This is equivalent to evaluating lim. lim.

Math 1272 Solutions for Spring 2005 Final Exam. asked to find the limit of the sequence. This is equivalent to evaluating lim. lim. Math 7 Solutions for Spring 5 Final Exam ) We are gien an infinite sequence for which the general term is a n 3 + 5n n + n an are 3 + 5n aske to fin the limit of the sequence. This is equialent to ealuating

More information

The Kepler Problem. 1 Features of the Ellipse: Geometry and Analysis

The Kepler Problem. 1 Features of the Ellipse: Geometry and Analysis The Kepler Problem For the Newtonian 1/r force law, a miracle occurs all of the solutions are perioic instea of just quasiperioic. To put it another way, the two-imensional tori are further ecompose into

More information

MATHEMATICS BONUS FILES for faculty and students

MATHEMATICS BONUS FILES for faculty and students MATHMATI BONU FIL for faculty an stuents http://www.onu.eu/~mcaragiu1/bonus_files.html RIVD: May 15, 9 PUBLIHD: May 5, 9 toffel 1 Maxwell s quations through the Major Vector Theorems Joshua toffel Department

More information

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent

More information

1 Heisenberg Representation

1 Heisenberg Representation 1 Heisenberg Representation What we have been ealing with so far is calle the Schröinger representation. In this representation, operators are constants an all the time epenence is carrie by the states.

More information

How the potentials in different gauges yield the same retarded electric and magnetic fields

How the potentials in different gauges yield the same retarded electric and magnetic fields How the potentials in ifferent gauges yiel the same retare electric an magnetic fiels José A. Heras a Departamento e Física, E. S. F. M., Instituto Politécnico Nacional, México D. F. México an Department

More information

23 Implicit differentiation

23 Implicit differentiation 23 Implicit ifferentiation 23.1 Statement The equation y = x 2 + 3x + 1 expresses a relationship between the quantities x an y. If a value of x is given, then a corresponing value of y is etermine. For

More information

Math 1271 Solutions for Fall 2005 Final Exam

Math 1271 Solutions for Fall 2005 Final Exam Math 7 Solutions for Fall 5 Final Eam ) Since the equation + y = e y cannot be rearrange algebraically in orer to write y as an eplicit function of, we must instea ifferentiate this relation implicitly

More information

Lagrangian and Hamiltonian Mechanics

Lagrangian and Hamiltonian Mechanics Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical

More information

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson JUST THE MATHS UNIT NUMBER 10.2 DIFFERENTIATION 2 (Rates of change) by A.J.Hobson 10.2.1 Introuction 10.2.2 Average rates of change 10.2.3 Instantaneous rates of change 10.2.4 Derivatives 10.2.5 Exercises

More information

TAYLOR S POLYNOMIAL APPROXIMATION FOR FUNCTIONS

TAYLOR S POLYNOMIAL APPROXIMATION FOR FUNCTIONS MISN-0-4 TAYLOR S POLYNOMIAL APPROXIMATION FOR FUNCTIONS f(x ± ) = f(x) ± f ' (x) + f '' (x) 2 ±... 1! 2! = 1.000 ± 0.100 + 0.005 ±... TAYLOR S POLYNOMIAL APPROXIMATION FOR FUNCTIONS by Peter Signell 1.

More information

Solutions to Math 41 Second Exam November 4, 2010

Solutions to Math 41 Second Exam November 4, 2010 Solutions to Math 41 Secon Exam November 4, 2010 1. (13 points) Differentiate, using the metho of your choice. (a) p(t) = ln(sec t + tan t) + log 2 (2 + t) (4 points) Using the rule for the erivative of

More information

Day 4: Motion Along a Curve Vectors

Day 4: Motion Along a Curve Vectors Day 4: Motion Along a Curve Vectors I give my stuents the following list of terms an formulas to know. Parametric Equations, Vectors, an Calculus Terms an Formulas to Know: If a smooth curve C is given

More information

Strength Analysis of CFRP Composite Material Considering Multiple Fracture Modes

Strength Analysis of CFRP Composite Material Considering Multiple Fracture Modes 5--XXXX Strength Analysis of CFRP Composite Material Consiering Multiple Fracture Moes Author, co-author (Do NOT enter this information. It will be pulle from participant tab in MyTechZone) Affiliation

More information

and from it produce the action integral whose variation we set to zero:

and from it produce the action integral whose variation we set to zero: Lagrange Multipliers Monay, 6 September 01 Sometimes it is convenient to use reunant coorinates, an to effect the variation of the action consistent with the constraints via the metho of Lagrange unetermine

More information

inflow outflow Part I. Regular tasks for MAE598/494 Task 1

inflow outflow Part I. Regular tasks for MAE598/494 Task 1 MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the

More information

1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7.

1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7. Lectures Nine an Ten The WKB Approximation The WKB metho is a powerful tool to obtain solutions for many physical problems It is generally applicable to problems of wave propagation in which the frequency

More information

G4003 Advanced Mechanics 1. We already saw that if q is a cyclic variable, the associated conjugate momentum is conserved, L = const.

G4003 Advanced Mechanics 1. We already saw that if q is a cyclic variable, the associated conjugate momentum is conserved, L = const. G4003 Avance Mechanics 1 The Noether theorem We alreay saw that if q is a cyclic variable, the associate conjugate momentum is conserve, q = 0 p q = const. (1) This is the simplest incarnation of Noether

More information

6. Friction and viscosity in gasses

6. Friction and viscosity in gasses IR2 6. Friction an viscosity in gasses 6.1 Introuction Similar to fluis, also for laminar flowing gases Newtons s friction law hols true (see experiment IR1). Using Newton s law the viscosity of air uner

More information

The Exact Form and General Integrating Factors

The Exact Form and General Integrating Factors 7 The Exact Form an General Integrating Factors In the previous chapters, we ve seen how separable an linear ifferential equations can be solve using methos for converting them to forms that can be easily

More information

Homework 7 Due 18 November at 6:00 pm

Homework 7 Due 18 November at 6:00 pm Homework 7 Due 18 November at 6:00 pm 1. Maxwell s Equations Quasi-statics o a An air core, N turn, cylinrical solenoi of length an raius a, carries a current I Io cos t. a. Using Ampere s Law, etermine

More information

1 Lecture 20: Implicit differentiation

1 Lecture 20: Implicit differentiation Lecture 20: Implicit ifferentiation. Outline The technique of implicit ifferentiation Tangent lines to a circle Derivatives of inverse functions by implicit ifferentiation Examples.2 Implicit ifferentiation

More information

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION The Annals of Statistics 1997, Vol. 25, No. 6, 2313 2327 LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION By Eva Riccomagno, 1 Rainer Schwabe 2 an Henry P. Wynn 1 University of Warwick, Technische

More information

Calculus of Variations

Calculus of Variations Calculus of Variations Lagrangian formalism is the main tool of theoretical classical mechanics. Calculus of Variations is a part of Mathematics which Lagrangian formalism is base on. In this section,

More information

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y Ph195a lecture notes, 1/3/01 Density operators for spin- 1 ensembles So far in our iscussion of spin- 1 systems, we have restricte our attention to the case of pure states an Hamiltonian evolution. Toay

More information

Problem Set 2: Solutions

Problem Set 2: Solutions UNIVERSITY OF ALABAMA Department of Physics an Astronomy PH 102 / LeClair Summer II 2010 Problem Set 2: Solutions 1. The en of a charge rubber ro will attract small pellets of Styrofoam that, having mae

More information

Vectors in two dimensions

Vectors in two dimensions Vectors in two imensions Until now, we have been working in one imension only The main reason for this is to become familiar with the main physical ieas like Newton s secon law, without the aitional complication

More information

ELECTRON DIFFRACTION

ELECTRON DIFFRACTION ELECTRON DIFFRACTION Electrons : wave or quanta? Measurement of wavelength an momentum of electrons. Introuction Electrons isplay both wave an particle properties. What is the relationship between the

More information

In the usual geometric derivation of Bragg s Law one assumes that crystalline

In the usual geometric derivation of Bragg s Law one assumes that crystalline Diffraction Principles In the usual geometric erivation of ragg s Law one assumes that crystalline arrays of atoms iffract X-rays just as the regularly etche lines of a grating iffract light. While this

More information

MATH2231-Differentiation (2)

MATH2231-Differentiation (2) -Differentiation () The Beginnings of Calculus The prime occasion from which arose my iscovery of the metho of the Characteristic Triangle, an other things of the same sort, happene at a time when I ha

More information

Mathematics. Circles. hsn.uk.net. Higher. Contents. Circles 1. CfE Edition

Mathematics. Circles. hsn.uk.net. Higher. Contents. Circles 1. CfE Edition Higher Mathematics Contents 1 1 Representing a Circle A 1 Testing a Point A 3 The General Equation of a Circle A 4 Intersection of a Line an a Circle A 4 5 Tangents to A 5 6 Equations of Tangents to A

More information

Separation of Variables

Separation of Variables Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical

More information

Gravitation as the result of the reintegration of migrated electrons and positrons to their atomic nuclei. Osvaldo Domann

Gravitation as the result of the reintegration of migrated electrons and positrons to their atomic nuclei. Osvaldo Domann Gravitation as the result of the reintegration of migrate electrons an positrons to their atomic nuclei. Osvalo Domann oomann@yahoo.com (This paper is an extract of [6] liste in section Bibliography.)

More information

Tutorial Test 5 2D welding robot

Tutorial Test 5 2D welding robot Tutorial Test 5 D weling robot Phys 70: Planar rigi boy ynamics The problem statement is appene at the en of the reference solution. June 19, 015 Begin: 10:00 am En: 11:30 am Duration: 90 min Solution.

More information

Schrödinger s equation.

Schrödinger s equation. Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of

More information

Quantum Mechanics in Three Dimensions

Quantum Mechanics in Three Dimensions Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.

More information

This section outlines the methodology used to calculate the wave load and wave wind load values.

This section outlines the methodology used to calculate the wave load and wave wind load values. COMPUTERS AND STRUCTURES, INC., JUNE 2014 AUTOMATIC WAVE LOADS TECHNICAL NOTE CALCULATION O WAVE LOAD VALUES This section outlines the methoology use to calculate the wave loa an wave win loa values. Overview

More information

Tensors, Fields Pt. 1 and the Lie Bracket Pt. 1

Tensors, Fields Pt. 1 and the Lie Bracket Pt. 1 Tensors, Fiels Pt. 1 an the Lie Bracket Pt. 1 PHYS 500 - Southern Illinois University September 8, 2016 PHYS 500 - Southern Illinois University Tensors, Fiels Pt. 1 an the Lie Bracket Pt. 1 September 8,

More information

Diagonalization of Matrices Dr. E. Jacobs

Diagonalization of Matrices Dr. E. Jacobs Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is

More information

De Broglie s Pilot Waves

De Broglie s Pilot Waves De Broglie s Pilot Waves Bohr s Moel of the Hyrogen tom: One way to arrive at Bohr s hypothesis is to think of the electron not as a particle but as a staning wave at raius r aroun the proton. Thus, nλ

More information

Calculus and optimization

Calculus and optimization Calculus an optimization These notes essentially correspon to mathematical appenix 2 in the text. 1 Functions of a single variable Now that we have e ne functions we turn our attention to calculus. A function

More information

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation

More information

Further Differentiation and Applications

Further Differentiation and Applications Avance Higher Notes (Unit ) Prerequisites: Inverse function property; prouct, quotient an chain rules; inflexion points. Maths Applications: Concavity; ifferentiability. Real-Worl Applications: Particle

More information

Chapter-2. Steady Stokes flow around deformed sphere. class of oblate axi-symmetric bodies

Chapter-2. Steady Stokes flow around deformed sphere. class of oblate axi-symmetric bodies hapter- Steay Stoes flow aroun eforme sphere. class of oblate axi-symmetric boies. General In physical an biological sciences, an in engineering, there is a wie range of problems of interest lie seimentation

More information

Applications of the Wronskian to ordinary linear differential equations

Applications of the Wronskian to ordinary linear differential equations Physics 116C Fall 2011 Applications of the Wronskian to orinary linear ifferential equations Consier a of n continuous functions y i (x) [i = 1,2,3,...,n], each of which is ifferentiable at least n times.

More information

ector ition 1. Scalar: physical quantity having only magnitue but no irection is calle a scalar. eg: Time, mass, istance, spee, electric charge, etc.. ector: physical quantity having both magnitue an irection

More information

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5

More information

PHY 114 Summer 2009 Final Exam Solutions

PHY 114 Summer 2009 Final Exam Solutions PHY 4 Summer 009 Final Exam Solutions Conceptual Question : A spherical rubber balloon has a charge uniformly istribute over its surface As the balloon is inflate, how oes the electric fiel E vary (a)

More information

Gravitation as the result of the reintegration of migrated electrons and positrons to their atomic nuclei. Osvaldo Domann

Gravitation as the result of the reintegration of migrated electrons and positrons to their atomic nuclei. Osvaldo Domann Gravitation as the result of the reintegration of migrate electrons an positrons to their atomic nuclei. Osvalo Domann oomann@yahoo.com (This paper is an extract of [6] liste in section Bibliography.)

More information

5.4 Fundamental Theorem of Calculus Calculus. Do you remember the Fundamental Theorem of Algebra? Just thought I'd ask

5.4 Fundamental Theorem of Calculus Calculus. Do you remember the Fundamental Theorem of Algebra? Just thought I'd ask 5.4 FUNDAMENTAL THEOREM OF CALCULUS Do you remember the Funamental Theorem of Algebra? Just thought I' ask The Funamental Theorem of Calculus has two parts. These two parts tie together the concept of

More information

A. Incorrect! The letter t does not appear in the expression of the given integral

A. Incorrect! The letter t does not appear in the expression of the given integral AP Physics C - Problem Drill 1: The Funamental Theorem of Calculus Question No. 1 of 1 Instruction: (1) Rea the problem statement an answer choices carefully () Work the problems on paper as neee (3) Question

More information

Math 1B, lecture 8: Integration by parts

Math 1B, lecture 8: Integration by parts Math B, lecture 8: Integration by parts Nathan Pflueger 23 September 2 Introuction Integration by parts, similarly to integration by substitution, reverses a well-known technique of ifferentiation an explores

More information

Calculus in the AP Physics C Course The Derivative

Calculus in the AP Physics C Course The Derivative Limits an Derivatives Calculus in the AP Physics C Course The Derivative In physics, the ieas of the rate change of a quantity (along with the slope of a tangent line) an the area uner a curve are essential.

More information

ECE 6310 Spring 2012 Exam 1 Solutions. Balanis The electric fields are given by. E r = ˆxe jβ 0 z

ECE 6310 Spring 2012 Exam 1 Solutions. Balanis The electric fields are given by. E r = ˆxe jβ 0 z ECE 6310 Spring 2012 Exam 1 Solutions Balanis 1.30 The electric fiels are given by E i ˆxe jβ 0 z E r ˆxe jβ 0 z The curl of the electric fiels are the usual cross prouct E i jβ 0 ẑ ˆxe jβ 0 z jβ 0 ŷe

More information

Free rotation of a rigid body 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012

Free rotation of a rigid body 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012 Free rotation of a rigi boy 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012 1 Introuction In this section, we escribe the motion of a rigi boy that is free to rotate

More information

Research Article Relativistic Propagation of Linearly/Circularly Polarized Laser Radiation in Plasmas

Research Article Relativistic Propagation of Linearly/Circularly Polarized Laser Radiation in Plasmas ISRN Optics Volume 213, Article ID 642617, 8 pages http://x.oi.org/1.1155/213/642617 Research Article Relativistic Propagation of Linearly/Circularly Polarize Laser Raiation in Plasmas Sonu Sen, 1 Meenu

More information

A note on the Mooney-Rivlin material model

A note on the Mooney-Rivlin material model A note on the Mooney-Rivlin material moel I-Shih Liu Instituto e Matemática Universiae Feeral o Rio e Janeiro 2945-97, Rio e Janeiro, Brasil Abstract In finite elasticity, the Mooney-Rivlin material moel

More information

d dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1

d dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1 Lecture 5 Some ifferentiation rules Trigonometric functions (Relevant section from Stewart, Seventh Eition: Section 3.3) You all know that sin = cos cos = sin. () But have you ever seen a erivation of

More information

Static Equilibrium. Theory: The conditions for the mechanical equilibrium of a rigid body are (a) (b)

Static Equilibrium. Theory: The conditions for the mechanical equilibrium of a rigid body are (a) (b) LPC Physics A 00 Las Positas College, Physics Department Staff Purpose: To etermine that, for a boy in equilibrium, the following are true: The sum of the torques about any point is zero The sum of forces

More information

Convective heat transfer

Convective heat transfer CHAPTER VIII Convective heat transfer The previous two chapters on issipative fluis were evote to flows ominate either by viscous effects (Chap. VI) or by convective motion (Chap. VII). In either case,

More information

A Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential

A Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential Avances in Applie Mathematics an Mechanics Av. Appl. Math. Mech. Vol. 1 No. 4 pp. 573-580 DOI: 10.4208/aamm.09-m0946 August 2009 A Note on Exact Solutions to Linear Differential Equations by the Matrix

More information

Chapter 6: Energy-Momentum Tensors

Chapter 6: Energy-Momentum Tensors 49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.

More information

New Zealand Institute of Physics

New Zealand Institute of Physics New Zealan Institute of Physics ASSESSMENT SCHEDULE Physics Level 2 90258 v2 Demonstrate unerstaning of physics in an integrate context Note: Minor computational errors will not be penalise. A wrong answer

More information

Hyperbolic Systems of Equations Posed on Erroneous Curved Domains

Hyperbolic Systems of Equations Posed on Erroneous Curved Domains Hyperbolic Systems of Equations Pose on Erroneous Curve Domains Jan Norström a, Samira Nikkar b a Department of Mathematics, Computational Mathematics, Linköping University, SE-58 83 Linköping, Sween (

More information

Generalization of the persistent random walk to dimensions greater than 1

Generalization of the persistent random walk to dimensions greater than 1 PHYSICAL REVIEW E VOLUME 58, NUMBER 6 DECEMBER 1998 Generalization of the persistent ranom walk to imensions greater than 1 Marián Boguñá, Josep M. Porrà, an Jaume Masoliver Departament e Física Fonamental,

More information

Chapter 9 Method of Weighted Residuals

Chapter 9 Method of Weighted Residuals Chapter 9 Metho of Weighte Resiuals 9- Introuction Metho of Weighte Resiuals (MWR) is an approimate technique for solving bounary value problems. It utilizes a trial functions satisfying the prescribe

More information

Basic Differentiation Rules and Rates of Change. The Constant Rule

Basic Differentiation Rules and Rates of Change. The Constant Rule 460_00.q //04 4:04 PM Page 07 SECTION. Basic Differentiation Rules an Rates of Change 07 Section. The slope of a horizontal line is 0. Basic Differentiation Rules an Rates of Change Fin the erivative of

More information

A Second Time Dimension, Hidden in Plain Sight

A Second Time Dimension, Hidden in Plain Sight A Secon Time Dimension, Hien in Plain Sight Brett A Collins. In this paper I postulate the existence of a secon time imension, making five imensions, three space imensions an two time imensions. I will

More information

Radial Basis-Function Networks

Radial Basis-Function Networks Raial Basis-Function Networks Back-Propagation Stochastic Back-Propagation Algorithm Step by Step Example Raial Basis-Function Networks Gaussian response function Location of center u Determining sigma

More information

Optimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations

Optimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations Optimize Schwarz Methos with the Yin-Yang Gri for Shallow Water Equations Abessama Qaouri Recherche en prévision numérique, Atmospheric Science an Technology Directorate, Environment Canaa, Dorval, Québec,

More information

A simple model for the small-strain behaviour of soils

A simple model for the small-strain behaviour of soils A simple moel for the small-strain behaviour of soils José Jorge Naer Department of Structural an Geotechnical ngineering, Polytechnic School, University of São Paulo 05508-900, São Paulo, Brazil, e-mail:

More information

Experimental Robustness Study of a Second-Order Sliding Mode Controller

Experimental Robustness Study of a Second-Order Sliding Mode Controller Experimental Robustness Stuy of a Secon-Orer Sliing Moe Controller Anré Blom, Bram e Jager Einhoven University of Technology Department of Mechanical Engineering P.O. Box 513, 5600 MB Einhoven, The Netherlans

More information

Witten s Proof of Morse Inequalities

Witten s Proof of Morse Inequalities Witten s Proof of Morse Inequalities by Igor Prokhorenkov Let M be a smooth, compact, oriente manifol with imension n. A Morse function is a smooth function f : M R such that all of its critical points

More information

5-4 Electrostatic Boundary Value Problems

5-4 Electrostatic Boundary Value Problems 11/8/4 Section 54 Electrostatic Bounary Value Problems blank 1/ 5-4 Electrostatic Bounary Value Problems Reaing Assignment: pp. 149-157 Q: A: We must solve ifferential equations, an apply bounary conitions

More information

Track Initialization from Incomplete Measurements

Track Initialization from Incomplete Measurements Track Initialiation from Incomplete Measurements Christian R. Berger, Martina Daun an Wolfgang Koch Department of Electrical an Computer Engineering, University of Connecticut, Storrs, Connecticut 6269,

More information

Final Exam Study Guide and Practice Problems Solutions

Final Exam Study Guide and Practice Problems Solutions Final Exam Stuy Guie an Practice Problems Solutions Note: These problems are just some of the types of problems that might appear on the exam. However, to fully prepare for the exam, in aition to making

More information

Unit #4 - Inverse Trig, Interpreting Derivatives, Newton s Method

Unit #4 - Inverse Trig, Interpreting Derivatives, Newton s Method Unit #4 - Inverse Trig, Interpreting Derivatives, Newton s Metho Some problems an solutions selecte or aapte from Hughes-Hallett Calculus. Computing Inverse Trig Derivatives. Starting with the inverse

More information

Physics 2212 GJ Quiz #4 Solutions Fall 2015

Physics 2212 GJ Quiz #4 Solutions Fall 2015 Physics 2212 GJ Quiz #4 Solutions Fall 215 I. (17 points) The magnetic fiel at point P ue to a current through the wire is 5. µt into the page. The curve portion of the wire is a semicircle of raius 2.

More information