Introduction to Astrophysics Tutorial 2: Polytropic Models
|
|
- Jean Dean
- 5 years ago
- Views:
Transcription
1 Itroductio to Astrophysics Tutorial : Polytropic Models Iair Arcavi 1 Summary of the Equatios of Stellar Structure We have arrived at a set of dieretial equatios which ca be used to describe the structure of a star i equilibrium. The equatios ca be writte either as fuctio of the radius coordiate r or the mass coordiate m: Hydrostatic Equilibrium: dr Gm = ρ r dm = Gm 4πr 4 Cotiuity: dm dr = 4πr ρ dr dm = 1 4πr ρ Radiative Trasfer: (if eergy is trasferred oly by radiative diusio) dt dr = 3 4ac κρ T 3 F 4πr dt dm = 3 4ac κ T 3 F (4πr ) Thermal Equilibrium: df dr = 4πr df ρq dm = q These structure equatios are supplemeted by: P = P I + P e + P rad = R µ I ρt + P e at 4 κ = κ 0 ρ a T b q = q 0 ρ m T Solvig these equatios is very dicult, sice they are o-liear, coupled ad have two-poit boudary coditios. However, we ca gai isight ito the structure of stars by aalyzig the equatios without solvig them, ad by usig simple models. Oe such model is called the Polytropic Model. 1
2 The Polytropic Model Note that the rst pair of equatios is coected to the secod pair by the equatio of state. However, if the pressure is idepedet of temperature, the the rst pair of equatios is separated from the secod pair. Multiplyig the equatio of hydrostatic equilibrium by r /ρ ad dieretiatig with respect to r, we d: ( ) d r = G dm dr ρ dr dr Usig the cotiuity equatio o the right had side gives: ( ) 1 d r r = 4πGρ dr ρ dr Now eters the model. We cosider equatios of state of the form: P = Kρ γ where K ad γ are costats. This is kow as a polytropic equatio of state, ad we have see that it pops up with degeerate gases ad adiabatic processes. We dee the polytropic idex as: γ = so a degeerate gas has idex = 1.5 for the o-relativistic case ad = 3 for the ultra relativistic case. Usig this model i the structure equatio we costructed gives: ( ) ( + 1) K 1 d r dρ 4πG r = ρ dr dr With the boudary coditios: ρ 1 ρ (R) = 0 dρ dr = 0 r=0 (where the secod coditio comes from the equatio of hydrostatic equilibrium). The solutio ρ (r) is called a polytrope, ad is uiquely deed by K, ad R. It is coveiet to dee dimesioless variables θ ad ξ by: ρ = ρ c θ r = αξ where ρ c is the cetral desity ad α = ( + 1) K 4πGρ 1 c
3 Our equatio the becomes the Lae-Emde equatio: with boudary coditios: ( 1 d ξ ξ ) = θ θ ξ=0 = 1 = 0 ξ=0 For < 5, θ (ξ) decreases mootoically, ad the radius of the star ca be foud by lookig for the rst zero of θ (usually umerically), called ξ 1 (the R = αξ 1 ), sice at the surface ρ = 0 so θ = 0. There are oly three cases where a aalytical solutio to the Lae-Emde equatio exists: 1. = 0, θ 0 = 1 ξ 6, ξ 1 = 6, γ =. = 1, θ 1 = si ξ ξ, ξ 1 = π, γ = 1 3. = 5, θ 5 =, ξ 1 =, γ = 6 1+ξ /3 5 I case 1, we have P = P c θ, ad this is applicable for a icompressible uid (sice = 0 meas costat desity). Case 3 is ot physical because the desity ever reaches zero. The total mass of a star ca be foud by: M = ˆR 0 4πr ρdr = 4πα 3 ρ c Substitutig from the Lae-Emde equatios we have: which gives: The mea desity ρ = M = 4πα 3 ρ c ˆξ 1 0 ˆξ 1 0 d M = 4πα 3 ρ c ξ 1 ξ θ ( ξ ) M 4 3 πr3 ca be show to be related to the cetral desity by: ρ c = D ρ ξ1 3
4 with: [ 3 D = ξ 1 ξ1 ] 1 This eables us to d a relatio betwee the mass ad the radius: ( ) 1 ( ) 3 GM R = M R [( + 1) K] 4πG with: M = ξ 1 R = ξ 1 ξ1 It is iterestig to see what happes i the special cases = 1 ad = 3. For = 3, the mass is idepedet of the radius: M = 4πM 3 ( K πg while for = 1 we ca d a radius idepedet of mass: ) 3/ ( ) 1/ K R = R 1 πg Geerally, for betwee these values, we ote that: R 3 1 M 1 The more massive the star, the smaller (ad hece deser) it gets! Fially, we would like to d a expressio for the cetral pressure, which ca be doe by substitutig the cetral desity ito the equatio of state P = Kρ 1+ 1 : which ca be writte as: P c = (4πG) ( GM M ) 1 ( R R P c = (4π) 1/3 B GM /3 ρ 4/3 c ) 3 ρ +1 c where B icludes all the -depedet coeciets. It turs out that B varies quite slowly with, implyig that the above expressio is almost uiversal, ad therefore will be useful later. 3 The Chadrasekhar Mass Stars which are domiated by degeeracy pressure, ca be described by a polytropic equatio of state with idex = 1.5 ad K = K 1 (from the previous tutorial). White dwarfs are a example of such 4
5 stars (havig roughly the mass of the su but the radius of the earth, givig them a mea desity of ρ 10 5 g/cm 3 ). For such stars the mass-radius relatio is: R M 1/3 ad so the average desity icreases with the mass as: ρ MR 3 M For icreasig mass, the, the star becomes deser, ad so ultimately its electros become relativistic. I this case, = 3, for which the mass-radius relatio gives a mass that is idepedet of the radius: M = 4πM 3 ( K πg ) 3/ Pluggig i K = K from the relativistic degeerate equatio of state, gives the Chadrasekhar Mass: M Ch = M ( π Substitutig all the costats, this ca be writte as: hc Gm 4/3 H M Ch = 5.83 µ M e ) 3/ µ e This is the maximal mass a stable coguratio ca have. For a white dwarf composed mostly of He ad heavier elemets, for which we ca take µ e =, we d M Ch = 1.46M. 5
6 6
Lecture 5-2: Polytropes. Literature: MWW chapter 19
Lecture 5-2: Polytropes Literature: MWW chapter 9!" Preamble The 4 equatios of stellar structure divide ito two groups: Mass ad mometum describig the mechaical structure ad thermal equilibrium ad eergy
More informationpoint, requiring all 4 curves to be continuous Discontinuity in P
. Solutio Methods a Classical approach Basic equatios of stellar structure + boudary coditios i Classical (historical method of solutio: Outward itegratio from ceter, where l = m = ward itegratio from
More informationdp dr = ρgm r 2 (1) dm dt dr = 3
Week 5: Simple equilibrium configurations for stars Follows on closely from chapter 5 of Prialnik. Having understood gas, radiation and nuclear physics at the level necessary to have an understanding of
More informationECE Spring Prof. David R. Jackson ECE Dept. Notes 20
ECE 6341 Sprig 016 Prof. David R. Jackso ECE Dept. Notes 0 1 Spherical Wave Fuctios Cosider solvig ψ + k ψ = 0 i spherical coordiates z φ θ r y x Spherical Wave Fuctios (cot.) I spherical coordiates we
More informationDETERMINATION OF MECHANICAL PROPERTIES OF A NON- UNIFORM BEAM USING THE MEASUREMENT OF THE EXCITED LONGITUDINAL ELASTIC VIBRATIONS.
ICSV4 Cairs Australia 9- July 7 DTRMINATION OF MCHANICAL PROPRTIS OF A NON- UNIFORM BAM USING TH MASURMNT OF TH XCITD LONGITUDINAL LASTIC VIBRATIONS Pavel Aokhi ad Vladimir Gordo Departmet of the mathematics
More informationStreamfunction-Vorticity Formulation
Streamfuctio-Vorticity Formulatio A. Salih Departmet of Aerospace Egieerig Idia Istitute of Space Sciece ad Techology, Thiruvaathapuram March 2013 The streamfuctio-vorticity formulatio was amog the first
More information( ) (( ) ) ANSWERS TO EXERCISES IN APPENDIX B. Section B.1 VECTORS AND SETS. Exercise B.1-1: Convex sets. are convex, , hence. and. (a) Let.
Joh Riley 8 Jue 03 ANSWERS TO EXERCISES IN APPENDIX B Sectio B VECTORS AND SETS Exercise B-: Covex sets (a) Let 0 x, x X, X, hece 0 x, x X ad 0 x, x X Sice X ad X are covex, x X ad x X The x X X, which
More informationThe axial dispersion model for tubular reactors at steady state can be described by the following equations: dc dz R n cn = 0 (1) (2) 1 d 2 c.
5.4 Applicatio of Perturbatio Methods to the Dispersio Model for Tubular Reactors The axial dispersio model for tubular reactors at steady state ca be described by the followig equatios: d c Pe dz z =
More informationFundamental Concepts: Surfaces and Curves
UNDAMENTAL CONCEPTS: SURACES AND CURVES CHAPTER udametal Cocepts: Surfaces ad Curves. INTRODUCTION This chapter describes two geometrical objects, vi., surfaces ad curves because the pla a ver importat
More informationPHY4905: Nearly-Free Electron Model (NFE)
PHY4905: Nearly-Free Electro Model (NFE) D. L. Maslov Departmet of Physics, Uiversity of Florida (Dated: Jauary 12, 2011) 1 I. REMINDER: QUANTUM MECHANICAL PERTURBATION THEORY A. No-degeerate eigestates
More informationCastiel, Supernatural, Season 6, Episode 18
13 Differetial Equatios the aswer to your questio ca best be epressed as a series of partial differetial equatios... Castiel, Superatural, Seaso 6, Episode 18 A differetial equatio is a mathematical equatio
More informationChapter 4. Fourier Series
Chapter 4. Fourier Series At this poit we are ready to ow cosider the caoical equatios. Cosider, for eample the heat equatio u t = u, < (4.) subject to u(, ) = si, u(, t) = u(, t) =. (4.) Here,
More informationWhere do eigenvalues/eigenvectors/eigenfunctions come from, and why are they important anyway?
Where do eigevalues/eigevectors/eigeuctios come rom, ad why are they importat ayway? I. Bacgroud (rom Ordiary Dieretial Equatios} Cosider the simplest example o a harmoic oscillator (thi o a vibratig strig)
More informationx = Pr ( X (n) βx ) =
Exercise 93 / page 45 The desity of a variable X i i 1 is fx α α a For α kow let say equal to α α > fx α α x α Pr X i x < x < Usig a Pivotal Quatity: x α 1 < x < α > x α 1 ad We solve i a similar way as
More informationLinear Regression Demystified
Liear Regressio Demystified Liear regressio is a importat subject i statistics. I elemetary statistics courses, formulae related to liear regressio are ofte stated without derivatio. This ote iteds to
More informationFinally, we show how to determine the moments of an impulse response based on the example of the dispersion model.
5.3 Determiatio of Momets Fially, we show how to determie the momets of a impulse respose based o the example of the dispersio model. For the dispersio model we have that E θ (θ ) curve is give by eq (4).
More informationPolynomials with Rational Roots that Differ by a Non-zero Constant. Generalities
Polyomials with Ratioal Roots that Differ by a No-zero Costat Philip Gibbs The problem of fidig two polyomials P(x) ad Q(x) of a give degree i a sigle variable x that have all ratioal roots ad differ by
More informationPHYC - 505: Statistical Mechanics Homework Assignment 4 Solutions
PHYC - 55: Statistical Mechaics Homewor Assigmet 4 Solutios Due February 5, 14 1. Cosider a ifiite classical chai of idetical masses coupled by earest eighbor sprigs with idetical sprig costats. a Write
More informationLecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting
Lecture 6 Chi Square Distributio (χ ) ad Least Squares Fittig Chi Square Distributio (χ ) Suppose: We have a set of measuremets {x 1, x, x }. We kow the true value of each x i (x t1, x t, x t ). We would
More informationLecture 25 (Dec. 6, 2017)
Lecture 5 8.31 Quatum Theory I, Fall 017 106 Lecture 5 (Dec. 6, 017) 5.1 Degeerate Perturbatio Theory Previously, whe discussig perturbatio theory, we restricted ourselves to the case where the uperturbed
More informationStatistical Inference Based on Extremum Estimators
T. Rotheberg Fall, 2007 Statistical Iferece Based o Extremum Estimators Itroductio Suppose 0, the true value of a p-dimesioal parameter, is kow to lie i some subset S R p : Ofte we choose to estimate 0
More informationCS284A: Representations and Algorithms in Molecular Biology
CS284A: Represetatios ad Algorithms i Molecular Biology Scribe Notes o Lectures 3 & 4: Motif Discovery via Eumeratio & Motif Represetatio Usig Positio Weight Matrix Joshua Gervi Based o presetatios by
More information3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,
3. Z Trasform Referece: Etire Chapter 3 of text. Recall that the Fourier trasform (FT) of a DT sigal x [ ] is ω ( ) [ ] X e = j jω k = xe I order for the FT to exist i the fiite magitude sese, S = x [
More informationLecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting
Lecture 6 Chi Square Distributio (χ ) ad Least Squares Fittig Chi Square Distributio (χ ) Suppose: We have a set of measuremets {x 1, x, x }. We kow the true value of each x i (x t1, x t, x t ). We would
More informationSeptember 2012 C1 Note. C1 Notes (Edexcel) Copyright - For AS, A2 notes and IGCSE / GCSE worksheets 1
September 0 s (Edecel) Copyright www.pgmaths.co.uk - For AS, A otes ad IGCSE / GCSE worksheets September 0 Copyright www.pgmaths.co.uk - For AS, A otes ad IGCSE / GCSE worksheets September 0 Copyright
More informationCEE 522 Autumn Uncertainty Concepts for Geotechnical Engineering
CEE 5 Autum 005 Ucertaity Cocepts for Geotechical Egieerig Basic Termiology Set A set is a collectio of (mutually exclusive) objects or evets. The sample space is the (collectively exhaustive) collectio
More informationM06/5/MATHL/HP2/ENG/TZ0/XX MATHEMATICS HIGHER LEVEL PAPER 2. Thursday 4 May 2006 (morning) 2 hours INSTRUCTIONS TO CANDIDATES
IB MATHEMATICS HIGHER LEVEL PAPER DIPLOMA PROGRAMME PROGRAMME DU DIPLÔME DU BI PROGRAMA DEL DIPLOMA DEL BI 06705 Thursday 4 May 006 (morig) hours INSTRUCTIONS TO CANDIDATES Do ot ope this examiatio paper
More informationA widely used display of protein shapes is based on the coordinates of the alpha carbons - - C α
Nice plottig of proteis: I A widely used display of protei shapes is based o the coordiates of the alpha carbos - - C α -s. The coordiates of the C α -s are coected by a cotiuous curve that roughly follows
More informationUNIVERSITY OF CALIFORNIA - SANTA CRUZ DEPARTMENT OF PHYSICS PHYS 116C. Problem Set 4. Benjamin Stahl. November 6, 2014
UNIVERSITY OF CALIFORNIA - SANTA CRUZ DEPARTMENT OF PHYSICS PHYS 6C Problem Set 4 Bejami Stahl November 6, 4 BOAS, P. 63, PROBLEM.-5 The Laguerre differetial equatio, x y + ( xy + py =, will be solved
More informationPOWER SERIES SOLUTION OF FIRST ORDER MATRIX DIFFERENTIAL EQUATIONS
Joural of Applied Mathematics ad Computatioal Mechaics 4 3(3) 3-8 POWER SERIES SOLUION OF FIRS ORDER MARIX DIFFERENIAL EQUAIONS Staisław Kukla Izabela Zamorska Istitute of Mathematics Czestochowa Uiversity
More informationPhysics 232 Gauge invariance of the magnetic susceptibilty
Physics 232 Gauge ivariace of the magetic susceptibilty Peter Youg (Dated: Jauary 16, 2006) I. INTRODUCTION We have see i class that the followig additioal terms appear i the Hamiltoia o addig a magetic
More informationRay Optics Theory and Mode Theory. Dr. Mohammad Faisal Dept. of EEE, BUET
Ray Optics Theory ad Mode Theory Dr. Mohammad Faisal Dept. of, BUT Optical Fiber WG For light to be trasmitted through fiber core, i.e., for total iteral reflectio i medium, > Ray Theory Trasmissio Ray
More informationReview Problems 1. ICME and MS&E Refresher Course September 19, 2011 B = C = AB = A = A 2 = A 3... C 2 = C 3 = =
Review Problems ICME ad MS&E Refresher Course September 9, 0 Warm-up problems. For the followig matrices A = 0 B = C = AB = 0 fid all powers A,A 3,(which is A times A),... ad B,B 3,... ad C,C 3,... Solutio:
More informationChE 471 Lecture 10 Fall 2005 SAFE OPERATION OF TUBULAR (PFR) ADIABATIC REACTORS
SAFE OPERATION OF TUBULAR (PFR) ADIABATIC REACTORS I a exothermic reactio the temperature will cotiue to rise as oe moves alog a plug flow reactor util all of the limitig reactat is exhausted. Schematically
More informationCHAPTER 5. Theory and Solution Using Matrix Techniques
A SERIES OF CLASS NOTES FOR 2005-2006 TO INTRODUCE LINEAR AND NONLINEAR PROBLEMS TO ENGINEERS, SCIENTISTS, AND APPLIED MATHEMATICIANS DE CLASS NOTES 3 A COLLECTION OF HANDOUTS ON SYSTEMS OF ORDINARY DIFFERENTIAL
More information1 Adiabatic and diabatic representations
1 Adiabatic ad diabatic represetatios 1.1 Bor-Oppeheimer approximatio The time-idepedet Schrödiger equatio for both electroic ad uclear degrees of freedom is Ĥ Ψ(r, R) = E Ψ(r, R), (1) where the full molecular
More informationNICK DUFRESNE. 1 1 p(x). To determine some formulas for the generating function of the Schröder numbers, r(x) = a(x) =
AN INTRODUCTION TO SCHRÖDER AND UNKNOWN NUMBERS NICK DUFRESNE Abstract. I this article we will itroduce two types of lattice paths, Schröder paths ad Ukow paths. We will examie differet properties of each,
More informationLinear Differential Equations of Higher Order Basic Theory: Initial-Value Problems d y d y dy
Liear Differetial Equatios of Higher Order Basic Theory: Iitial-Value Problems d y d y dy Solve: a( ) + a ( )... a ( ) a0( ) y g( ) + + + = d d d ( ) Subject to: y( 0) = y0, y ( 0) = y,..., y ( 0) = y
More informationTR/46 OCTOBER THE ZEROS OF PARTIAL SUMS OF A MACLAURIN EXPANSION A. TALBOT
TR/46 OCTOBER 974 THE ZEROS OF PARTIAL SUMS OF A MACLAURIN EXPANSION by A. TALBOT .. Itroductio. A problem i approximatio theory o which I have recetly worked [] required for its solutio a proof that the
More informationTopic 9: Sampling Distributions of Estimators
Topic 9: Samplig Distributios of Estimators Course 003, 2016 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be
More informationMost text will write ordinary derivatives using either Leibniz notation 2 3. y + 5y= e and y y. xx tt t
Itroductio to Differetial Equatios Defiitios ad Termiolog Differetial Equatio: A equatio cotaiig the derivatives of oe or more depedet variables, with respect to oe or more idepedet variables, is said
More informationNEW FAST CONVERGENT SEQUENCES OF EULER-MASCHERONI TYPE
UPB Sci Bull, Series A, Vol 79, Iss, 207 ISSN 22-7027 NEW FAST CONVERGENT SEQUENCES OF EULER-MASCHERONI TYPE Gabriel Bercu We itroduce two ew sequeces of Euler-Mascheroi type which have fast covergece
More informationSOLUTION SET VI FOR FALL [(n + 2)(n + 1)a n+2 a n 1 ]x n = 0,
4. Series Solutios of Differetial Equatios:Special Fuctios 4.. Illustrative examples.. 5. Obtai the geeral solutio of each of the followig differetial equatios i terms of Maclauri series: d y (a dx = xy,
More informationThe z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j
The -Trasform 7. Itroductio Geeralie the complex siusoidal represetatio offered by DTFT to a represetatio of complex expoetial sigals. Obtai more geeral characteristics for discrete-time LTI systems. 7.
More informationRecurrence Relations
Recurrece Relatios Aalysis of recursive algorithms, such as: it factorial (it ) { if (==0) retur ; else retur ( * factorial(-)); } Let t be the umber of multiplicatios eeded to calculate factorial(). The
More informationLinear Elliptic PDE s Elliptic partial differential equations frequently arise out of conservation statements of the form
Liear Elliptic PDE s Elliptic partial differetial equatios frequetly arise out of coservatio statemets of the form B F d B Sdx B cotaied i bouded ope set U R. Here F, S deote respectively, the flux desity
More informationThe time evolution of the state of a quantum system is described by the time-dependent Schrödinger equation (TDSE): ( ) ( ) 2m "2 + V ( r,t) (1.
Adrei Tokmakoff, MIT Departmet of Chemistry, 2/13/2007 1-1 574 TIME-DEPENDENT QUANTUM MECHANICS 1 INTRODUCTION 11 Time-evolutio for time-idepedet Hamiltoias The time evolutio of the state of a quatum system
More informationMathematical Methods for Physics and Engineering
Mathematical Methods for Physics ad Egieerig Lecture otes Sergei V. Shabaov Departmet of Mathematics, Uiversity of Florida, Gaiesville, FL 326 USA CHAPTER The theory of covergece. Numerical sequeces..
More informationMath 155 (Lecture 3)
Math 55 (Lecture 3) September 8, I this lecture, we ll cosider the aswer to oe of the most basic coutig problems i combiatorics Questio How may ways are there to choose a -elemet subset of the set {,,,
More informationDirection: This test is worth 250 points. You are required to complete this test within 50 minutes.
Term Test October 3, 003 Name Math 56 Studet Number Directio: This test is worth 50 poits. You are required to complete this test withi 50 miutes. I order to receive full credit, aswer each problem completely
More informationNonequilibrium Excess Carriers in Semiconductors
Lecture 8 Semicoductor Physics VI Noequilibrium Excess Carriers i Semicoductors Noequilibrium coditios. Excess electros i the coductio bad ad excess holes i the valece bad Ambiolar trasort : Excess electros
More informationAll Excuses must be taken to 233 Loomis before 4:15, Monday, April 30.
Miscellaeous Notes The ed is ear do t get behid. All Excuses must be take to 233 Loomis before 4:15, Moday, April 30. The PYS 213 fial exam times are * 8-10 AM, Moday, May 7 * 8-10 AM, Tuesday, May 8 ad
More informationMA541 : Real Analysis. Tutorial and Practice Problems - 1 Hints and Solutions
MA54 : Real Aalysis Tutorial ad Practice Problems - Hits ad Solutios. Suppose that S is a oempty subset of real umbers that is bouded (i.e. bouded above as well as below). Prove that if S sup S. What ca
More information17 Phonons and conduction electrons in solids (Hiroshi Matsuoka)
7 Phoos ad coductio electros i solids Hiroshi Matsuoa I this chapter we will discuss a miimal microscopic model for phoos i a solid ad a miimal microscopic model for coductio electros i a simple metal.
More informationMath 475, Problem Set #12: Answers
Math 475, Problem Set #12: Aswers A. Chapter 8, problem 12, parts (b) ad (d). (b) S # (, 2) = 2 2, sice, from amog the 2 ways of puttig elemets ito 2 distiguishable boxes, exactly 2 of them result i oe
More informationMiscellaneous Notes. Lecture 19, p 1
Miscellaeous Notes The ed is ear do t get behid. All Excuses must be take to 233 Loomis before oo, Thur, Apr. 25. The PHYS 213 fial exam times are * 8-10 AM, Moday, May 6 * 1:30-3:30 PM, Wed, May 8 The
More informationPhys 6303 Final Exam Solutions December 19, 2012
Phys 633 Fial Exam s December 19, 212 You may NOT use ay book or otes other tha supplied with this test. You will have 3 hours to fiish. DO YOUR OWN WORK. Express your aswers clearly ad cocisely so that
More informationKinetics of Complex Reactions
Kietics of Complex Reactios by Flick Colema Departmet of Chemistry Wellesley College Wellesley MA 28 wcolema@wellesley.edu Copyright Flick Colema 996. All rights reserved. You are welcome to use this documet
More informationFirst, note that the LS residuals are orthogonal to the regressors. X Xb X y = 0 ( normal equations ; (k 1) ) So,
0 2. OLS Part II The OLS residuals are orthogoal to the regressors. If the model icludes a itercept, the orthogoality of the residuals ad regressors gives rise to three results, which have limited practical
More information6a Time change b Quadratic variation c Planar Brownian motion d Conformal local martingales e Hints to exercises...
Tel Aviv Uiversity, 28 Browia motio 59 6 Time chage 6a Time chage..................... 59 6b Quadratic variatio................. 61 6c Plaar Browia motio.............. 64 6d Coformal local martigales............
More informationCS / MCS 401 Homework 3 grader solutions
CS / MCS 401 Homework 3 grader solutios assigmet due July 6, 016 writte by Jāis Lazovskis maximum poits: 33 Some questios from CLRS. Questios marked with a asterisk were ot graded. 1 Use the defiitio of
More informationLECTURE 14. Non-linear transverse motion. Non-linear transverse motion
LETURE 4 No-liear trasverse motio Floquet trasformatio Harmoic aalysis-oe dimesioal resoaces Two-dimesioal resoaces No-liear trasverse motio No-liear field terms i the trajectory equatio: Trajectory equatio
More informationSECTION 2 Electrostatics
SECTION Electrostatics This sectio, based o Chapter of Griffiths, covers effects of electric fields ad forces i static (timeidepedet) situatios. The topics are: Electric field Gauss s Law Electric potetial
More informationAnalytic Continuation
Aalytic Cotiuatio The stadard example of this is give by Example Let h (z) = 1 + z + z 2 + z 3 +... kow to coverge oly for z < 1. I fact h (z) = 1/ (1 z) for such z. Yet H (z) = 1/ (1 z) is defied for
More informationTHE NUMERICAL SOLUTION OF THE NEWTONIAN FLUIDS FLOW DUE TO A STRETCHING CYLINDER BY SOR ITERATIVE PROCEDURE ABSTRACT
Europea Joural of Egieerig ad Techology Vol. 3 No., 5 ISSN 56-586 THE NUMERICAL SOLUTION OF THE NEWTONIAN FLUIDS FLOW DUE TO A STRETCHING CYLINDER BY SOR ITERATIVE PROCEDURE Atif Nazir, Tahir Mahmood ad
More informationLarge holes in quasi-random graphs
Large holes i quasi-radom graphs Joaa Polcy Departmet of Discrete Mathematics Adam Mickiewicz Uiversity Pozań, Polad joaska@amuedupl Submitted: Nov 23, 2006; Accepted: Apr 10, 2008; Published: Apr 18,
More informationChapter 9: Numerical Differentiation
178 Chapter 9: Numerical Differetiatio Numerical Differetiatio Formulatio of equatios for physical problems ofte ivolve derivatives (rate-of-chage quatities, such as velocity ad acceleratio). Numerical
More informationDEGENERACY AND ALL THAT
DEGENERACY AND ALL THAT Te Nature of Termodyamics, Statistical Mecaics ad Classical Mecaics Termodyamics Te study of te equilibrium bulk properties of matter witi te cotext of four laws or facts of experiece
More informationCS321. Numerical Analysis and Computing
CS Numerical Aalysis ad Computig Lecture Locatig Roots o Equatios Proessor Ju Zhag Departmet o Computer Sciece Uiversity o Ketucky Leigto KY 456-6 September 8 5 What is the Root May physical system ca
More informationREAL ANALYSIS II: PROBLEM SET 1 - SOLUTIONS
REAL ANALYSIS II: PROBLEM SET 1 - SOLUTIONS 18th Feb, 016 Defiitio (Lipschitz fuctio). A fuctio f : R R is said to be Lipschitz if there exists a positive real umber c such that for ay x, y i the domai
More informationLet us give one more example of MLE. Example 3. The uniform distribution U[0, θ] on the interval [0, θ] has p.d.f.
Lecture 5 Let us give oe more example of MLE. Example 3. The uiform distributio U[0, ] o the iterval [0, ] has p.d.f. { 1 f(x =, 0 x, 0, otherwise The likelihood fuctio ϕ( = f(x i = 1 I(X 1,..., X [0,
More informationSemiconductor Statistical Mechanics (Read Kittel Ch. 8)
EE30 - Solid State Electroics Semicoductor Statistical Mechaics (Read Kittel Ch. 8) Coductio bad occupatio desity: f( E)gE ( ) de f(e) - occupatio probability - Fermi-Dirac fuctio: g(e) - desity of states
More informationAIT. Blackbody Radiation IAAT
3 1 Blackbody Radiatio Itroductio 3 2 First radiatio process to look at: radiatio i thermal equilibrium with itself: blackbody radiatio Assumptios: 1. Photos are Bosos, i.e., more tha oe photo per phase
More informationCS537. Numerical Analysis and Computing
CS57 Numerical Aalysis ad Computig Lecture Locatig Roots o Equatios Proessor Ju Zhag Departmet o Computer Sciece Uiversity o Ketucky Leigto KY 456-6 Jauary 9 9 What is the Root May physical system ca be
More informationUniversity of Colorado Denver Dept. Math. & Stat. Sciences Applied Analysis Preliminary Exam 13 January 2012, 10:00 am 2:00 pm. Good luck!
Uiversity of Colorado Dever Dept. Math. & Stat. Scieces Applied Aalysis Prelimiary Exam 13 Jauary 01, 10:00 am :00 pm Name: The proctor will let you read the followig coditios before the exam begis, ad
More informationSimilarity Solutions to Unsteady Pseudoplastic. Flow Near a Moving Wall
Iteratioal Mathematical Forum, Vol. 9, 04, o. 3, 465-475 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.988/imf.04.48 Similarity Solutios to Usteady Pseudoplastic Flow Near a Movig Wall W. Robi Egieerig
More informationInfinite Sequences and Series
Chapter 6 Ifiite Sequeces ad Series 6.1 Ifiite Sequeces 6.1.1 Elemetary Cocepts Simply speakig, a sequece is a ordered list of umbers writte: {a 1, a 2, a 3,...a, a +1,...} where the elemets a i represet
More informationNumerical Methods in Fourier Series Applications
Numerical Methods i Fourier Series Applicatios Recall that the basic relatios i usig the Trigoometric Fourier Series represetatio were give by f ( x) a o ( a x cos b x si ) () where the Fourier coefficiets
More informationMath 61CM - Solutions to homework 3
Math 6CM - Solutios to homework 3 Cédric De Groote October 2 th, 208 Problem : Let F be a field, m 0 a fixed oegative iteger ad let V = {a 0 + a x + + a m x m a 0,, a m F} be the vector space cosistig
More informationHydrogen (atoms, molecules) in external fields. Static electric and magnetic fields Oscyllating electromagnetic fields
Hydroge (atoms, molecules) i exteral fields Static electric ad magetic fields Oscyllatig electromagetic fields Everythig said up to ow has to be modified more or less strogly if we cosider atoms (ad ios)
More information62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +
62. Power series Defiitio 16. (Power series) Give a sequece {c }, the series c x = c 0 + c 1 x + c 2 x 2 + c 3 x 3 + is called a power series i the variable x. The umbers c are called the coefficiets of
More informationMATH 205 HOMEWORK #2 OFFICIAL SOLUTION. (f + g)(x) = f(x) + g(x) = f( x) g( x) = (f + g)( x)
MATH 205 HOMEWORK #2 OFFICIAL SOLUTION Problem 2: Do problems 7-9 o page 40 of Hoffma & Kuze. (7) We will prove this by cotradictio. Suppose that W 1 is ot cotaied i W 2 ad W 2 is ot cotaied i W 1. The
More informationLemma Let f(x) K[x] be a separable polynomial of degree n. Then the Galois group is a subgroup of S n, the permutations of the roots.
15 Cubics, Quartics ad Polygos It is iterestig to chase through the argumets of 14 ad see how this affects solvig polyomial equatios i specific examples We make a global assumptio that the characteristic
More informationEigenvalues and Eigenvectors
5 Eigevalues ad Eigevectors 5.3 DIAGONALIZATION DIAGONALIZATION Example 1: Let. Fid a formula for A k, give that P 1 1 = 1 2 ad, where Solutio: The stadard formula for the iverse of a 2 2 matrix yields
More information1. Hydrogen Atom: 3p State
7633A QUANTUM MECHANICS I - solutio set - autum. Hydroge Atom: 3p State Let us assume that a hydroge atom is i a 3p state. Show that the radial part of its wave fuctio is r u 3(r) = 4 8 6 e r 3 r(6 r).
More informationAssignment 2 Solutions SOLUTION. ϕ 1 Â = 3 ϕ 1 4i ϕ 2. The other case can be dealt with in a similar way. { ϕ 2 Â} χ = { 4i ϕ 1 3 ϕ 2 } χ.
PHYSICS 34 QUANTUM PHYSICS II (25) Assigmet 2 Solutios 1. With respect to a pair of orthoormal vectors ϕ 1 ad ϕ 2 that spa the Hilbert space H of a certai system, the operator  is defied by its actio
More informationCHAPTER 10 INFINITE SEQUENCES AND SERIES
CHAPTER 10 INFINITE SEQUENCES AND SERIES 10.1 Sequeces 10.2 Ifiite Series 10.3 The Itegral Tests 10.4 Compariso Tests 10.5 The Ratio ad Root Tests 10.6 Alteratig Series: Absolute ad Coditioal Covergece
More informationThe Born-Oppenheimer approximation
The Bor-Oppeheimer approximatio 1 Re-writig the Schrödiger equatio We will begi from the full time-idepedet Schrödiger equatio for the eigestates of a molecular system: [ P 2 + ( Pm 2 + e2 1 1 2m 2m m
More informationFor use only in [the name of your school] 2014 FP2 Note. FP2 Notes (Edexcel)
For use oly i [the ame of your school] 04 FP Note FP Notes (Edexcel) Copyright wwwpgmathscouk - For AS, A otes ad IGCSE / GCSE worksheets For use oly i [the ame of your school] 04 FP Note BLANK PAGE Copyright
More informationEconomics 241B Relation to Method of Moments and Maximum Likelihood OLSE as a Maximum Likelihood Estimator
Ecoomics 24B Relatio to Method of Momets ad Maximum Likelihood OLSE as a Maximum Likelihood Estimator Uder Assumptio 5 we have speci ed the distributio of the error, so we ca estimate the model parameters
More informationChapter Vectors
Chapter 4. Vectors fter readig this chapter you should be able to:. defie a vector. add ad subtract vectors. fid liear combiatios of vectors ad their relatioship to a set of equatios 4. explai what it
More informationECE-S352 Introduction to Digital Signal Processing Lecture 3A Direct Solution of Difference Equations
ECE-S352 Itroductio to Digital Sigal Processig Lecture 3A Direct Solutio of Differece Equatios Discrete Time Systems Described by Differece Equatios Uit impulse (sample) respose h() of a DT system allows
More information6 Integers Modulo n. integer k can be written as k = qn + r, with q,r, 0 r b. So any integer.
6 Itegers Modulo I Example 2.3(e), we have defied the cogruece of two itegers a,b with respect to a modulus. Let us recall that a b (mod ) meas a b. We have proved that cogruece is a equivalece relatio
More informationNumerical Methods for Ordinary Differential Equations
Numerical Methods for Ordiary Differetial Equatios Braislav K. Nikolić Departmet of Physics ad Astroomy, Uiversity of Delaware, U.S.A. PHYS 460/660: Computatioal Methods of Physics http://www.physics.udel.edu/~bikolic/teachig/phys660/phys660.html
More informationSimple Polygons of Maximum Perimeter Contained in a Unit Disk
Discrete Comput Geom (009) 1: 08 15 DOI 10.1007/s005-008-9093-7 Simple Polygos of Maximum Perimeter Cotaied i a Uit Disk Charles Audet Pierre Hase Frédéric Messie Received: 18 September 007 / Revised:
More informationSeunghee Ye Ma 8: Week 5 Oct 28
Week 5 Summary I Sectio, we go over the Mea Value Theorem ad its applicatios. I Sectio 2, we will recap what we have covered so far this term. Topics Page Mea Value Theorem. Applicatios of the Mea Value
More information4. Partial Sums and the Central Limit Theorem
1 of 10 7/16/2009 6:05 AM Virtual Laboratories > 6. Radom Samples > 1 2 3 4 5 6 7 4. Partial Sums ad the Cetral Limit Theorem The cetral limit theorem ad the law of large umbers are the two fudametal theorems
More informationAstronomy 112: The Physics of Stars. Class 9 Notes: Polytropes
Astronomy 112: The Physics of Stars Class 9 Notes: Polytropes With our discussion of nuclear reaction rates last time, we have mostly completed our survey of the microphysical properties of stellar matter
More informationPhysics 324, Fall Dirac Notation. These notes were produced by David Kaplan for Phys. 324 in Autumn 2001.
Physics 324, Fall 2002 Dirac Notatio These otes were produced by David Kapla for Phys. 324 i Autum 2001. 1 Vectors 1.1 Ier product Recall from liear algebra: we ca represet a vector V as a colum vector;
More informationE. Incorrect! Plug n = 1, 2, 3, & 4 into the general term formula. n =, then the first four terms are found by
Calculus II - Problem Solvig Drill 8: Sequeces, Series, ad Covergece Questio No. of 0. Fid the first four terms of the sequece whose geeral term is give by a ( ) : Questio #0 (A) (B) (C) (D) (E) 8,,, 4
More information