MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Statistical Physics I Spring Term Solutions to Problem Set #1

Similar documents
Name Solutions to Test 3 November 8, 2017

Quantum Physics I (8.04) Spring 2016 Assignment 8

(See Notes on Spontaneous Emission)

Quantum Physics II (8.05) Fall 2013 Assignment 2

221A Lecture Notes WKB Method

Physics 202H - Introductory Quantum Physics I Homework #08 - Solutions Fall 2004 Due 5:01 PM, Monday 2004/11/15

Improper Integrals, and Differential Equations

221B Lecture Notes WKB Method

Topic 1 Notes Jeremy Orloff

13.4 Work done by Constant Forces

Chapter 3 The Schrödinger Equation and a Particle in a Box

Simple Harmonic Motion I Sem

Continuous Quantum Systems

7 - Continuous random variables

1 Probability Density Functions

CHM Physical Chemistry I Chapter 1 - Supplementary Material

Summary: Method of Separation of Variables

Chapter 5. , r = r 1 r 2 (1) µ = m 1 m 2. r, r 2 = R µ m 2. R(m 1 + m 2 ) + m 2 r = r 1. m 2. r = r 1. R + µ m 1

Chapter 6 Notes, Larson/Hostetler 3e

Continuous Random Variables

SOLUTIONS FOR ADMISSIONS TEST IN MATHEMATICS, COMPUTER SCIENCE AND JOINT SCHOOLS WEDNESDAY 5 NOVEMBER 2014

Problem Set 3 Solutions

Massachusetts Institute of Technology Quantum Mechanics I (8.04) Spring 2005 Solutions to Problem Set 6

Math 8 Winter 2015 Applications of Integration

Physics Graduate Prelim exam

Jack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah

The Regulated and Riemann Integrals

7.2 The Definite Integral

Do the one-dimensional kinetic energy and momentum operators commute? If not, what operator does their commutator represent?

Quantum Mechanics Qualifying Exam - August 2016 Notes and Instructions

AP Calculus Multiple Choice: BC Edition Solutions

Practice Problems Solution

Jackson 2.26 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell

Energy Bands Energy Bands and Band Gap. Phys463.nb Phenomenon

University of Washington Department of Chemistry Chemistry 453 Winter Quarter 2010 Homework Assignment 4; Due at 5p.m. on 2/01/10

f(x) dx, If one of these two conditions is not met, we call the integral improper. Our usual definition for the value for the definite integral

Partial Derivatives. Limits. For a single variable function f (x), the limit lim

Improper Integrals. Type I Improper Integrals How do we evaluate an integral such as

Physics 116C Solution of inhomogeneous ordinary differential equations using Green s functions

Review of basic calculus

1B40 Practical Skills

A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007

f(a+h) f(a) x a h 0. This is the rate at which

CS667 Lecture 6: Monte Carlo Integration 02/10/05

Physics 215 Quantum Mechanics 1 Assignment 2

Describe in words how you interpret this quantity. Precisely what information do you get from x?

5.04 Principles of Inorganic Chemistry II

Consequently, the temperature must be the same at each point in the cross section at x. Let:

potentials A z, F z TE z Modes We use the e j z z =0 we can simply say that the x dependence of E y (1)

Overview of Calculus I

Math 231E, Lecture 33. Parametric Calculus

Calculus I-II Review Sheet

5.7 Improper Integrals

Quantum Physics III (8.06) Spring 2005 Solution Set 5

20 MATHEMATICS POLYNOMIALS

University of Washington Department of Chemistry Chemistry 453 Winter Quarter 2009

Best Approximation. Chapter The General Case

Math 1B, lecture 4: Error bounds for numerical methods

Physics 137A - Quantum Mechanics - Spring 2018 Midterm 1. Mathematical Formulas

Discrete Least-squares Approximations

MORE FUNCTION GRAPHING; OPTIMIZATION. (Last edited October 28, 2013 at 11:09pm.)

Math& 152 Section Integration by Parts

W. We shall do so one by one, starting with I 1, and we shall do it greedily, trying

Math 113 Exam 2 Practice

Physics 741 Graduate Quantum Mechanics 1 Solutions to Final Exam, Fall 2011

Unit 5. Integration techniques

Applied Physics Introduction to Vibrations and Waves (with a focus on elastic waves) Course Outline

Math Fall 2006 Sample problems for the final exam: Solutions

Physics 9 Fall 2011 Homework 2 - Solutions Friday September 2, 2011

Ph2b Quiz - 1. Instructions

Review of Gaussian Quadrature method

Math 31S. Rumbos Fall Solutions to Assignment #16

Orthogonal Polynomials and Least-Squares Approximations to Functions

Chapter 5 : Continuous Random Variables

ES.182A Topic 32 Notes Jeremy Orloff

Math 113 Fall Final Exam Review. 2. Applications of Integration Chapter 6 including sections and section 6.8

PHY4605 Introduction to Quantum Mechanics II Spring 2005 Final exam SOLUTIONS April 22, 2005

MTH 122 Fall 2008 Essex County College Division of Mathematics Handout Version 10 1 October 14, 2008

Math 426: Probability Final Exam Practice

Definition of Continuity: The function f(x) is continuous at x = a if f(a) exists and lim

Heat flux and total heat

Chapter 28. Fourier Series An Eigenvalue Problem.

Solutions to Problems in Merzbacher, Quantum Mechanics, Third Edition. Chapter 7

( ) 2. ( ) is the Fourier transform of! ( x). ( ) ( ) ( ) = Ae i kx"#t ( ) = 1 2" ( )"( x,t) PC 3101 Quantum Mechanics Section 1

This final is a three hour open book, open notes exam. Do all four problems.

Linear Differential Equations Physics 129a Solutions to Problems Frank Porter Revision F. Porter

Improper Integrals. The First Fundamental Theorem of Calculus, as we ve discussed in class, goes as follows:

Problem Set 9. Figure 1: Diagram. This picture is a rough sketch of the 4 parabolas that give us the area that we need to find. The equations are:

STEP FUNCTIONS, DELTA FUNCTIONS, AND THE VARIATION OF PARAMETERS FORMULA. 0 if t < 0, 1 if t > 0.

Orthogonal functions

Section 4.8. D v(t j 1 ) t. (4.8.1) j=1

Note 16. Stokes theorem Differential Geometry, 2005

Review of Calculus, cont d

Non-Linear & Logistic Regression

MATH 144: Business Calculus Final Review

4 The dynamical FRW universe

7.3 Problem 7.3. ~B(~x) = ~ k ~ E(~x)=! but we also have a reected wave. ~E(~x) = ~ E 2 e i~ k 2 ~x i!t. ~B R (~x) = ~ k R ~ E R (~x)=!

2.57/2.570 Midterm Exam No. 1 March 31, :00 am -12:30 pm

. Double-angle formulas. Your answer should involve trig functions of θ, and not of 2θ. sin 2 (θ) =

Designing Information Devices and Systems I Discussion 8B

Transcription:

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Deprtment 8.044 Sttisticl Physics I Spring Term 03 Problem : Doping Semiconductor Solutions to Problem Set # ) Mentlly integrte the function p(x) given in the figure. The result rises from zero t decresing rte, jumps discontinuously by 0. t x = d, then continues to rise symptoticlly towrd the vlue. This behvior is sketched below. b) < x > = 0.8 x p(x) dx = x exp( x/l) dx + 0. l } 0 0 {{}} l = 0.8 l + 0. d xδ(x d) dx {{} d c) < > = x 0.8 x p(x) dx = x exp( x/l) dx + 0. x δ(x d) dx l } 0 {{}} 0 {{} l 3 d =.6 l + 0. d Vr(x) < (x < x >) > = < x > < x > = (.6 l + 0. d ) (0.64 l + 0.3 ld + 0.04 d ) = 0.96 l 0.3 ld 0.6 d

d) < exp( x/s) > = exp( x/s) p(x) dx = 0.8 l exp( x/s) exp( x/l) dx + 0. exp( x/s)δ(x d) dx } 0 {{ } 0 } {{ } (/s + /l) exp( d/s) = ( ) 0.8 + 0. exp( d/s) + l/s Check to see tht this result is physiclly resonble. Note tht if the skin depth s is much less thn the distnce d, the impurities on the grin boundry do not contribute to the surfce impednce. Similrly, if the skin depth is much less thn the chrcteristic diffusion distnce l, the impurity contribution to the surfce impednce is gretly reduced.

Problem : A Peculir Probbility Density ) = p(x) dx = 0 dx b + x = dξ = (π/b) b } 0 + ξ {{ } π/ = (b/π) b) P (x) = x b x p(x ) dx = dx π b + x b = π [ x rctn(x /b) b = rctn(x/b) + π c) < x >= 0 by symmetry. p(x) is n even function nd x is odd. d) p(x) flls to hlf its vlue t x = ±b. e) b x < x >= dx π b + x However the limit of x /(b + x ) s x ± is unity, so this integrl diverges. Neither the men squre nor the Vrince of this distribution exist. 3

Problem 3: Visulizing the Probbility Density for Clssicl Hrmonic Oscilltor ) First find the velocity s function of time by tking the derivtive of the displcement with respect to time. ẋ(t) = d [x 0 sin(ωt + φ)] dt = ωx 0 cos(ωt + φ) But we don t wnt the velocity s function of t, we wnt it s function of the position x. And, we don t ctully need the velocity itself, we wnt the speed (the mgnitude of the velocity). Becuse of this we do not hve to worry bout losing the sign of the velocity when we work with its squre. ẋ (t) = (ωx 0 ) cos (ωt + φ) = (ωx ) [ sin 0 (ωt + φ)] = (ωx 0 ) [ (x(t)/x 0 ) ] Finlly, the speed is computed s the squre root of the squre of the velocity. ẋ(t) = ω(x 0 x (t)) / for x(t) x 0 b) We re told tht the probbility density for finding n oscilltor t x is proportionl to the the time given oscilltor spends ner x, nd tht this time is inversely proportionl to its speed t tht point. Expressed mthemticlly this becomes p(x) ẋ(t) = C(x 0 x ) / for x < x 0 where C is proportionlity constnt which we cn find by normlizing p(x). x0 p(x)dx = C (x 0 x ) / dx x 0 x0 dx/x 0 = C let x/x 0 y 0 (x/x 0) dy = C 0 y }} = πc {{ π/ = by normliztion The lst two lines imply tht C = /π. We cn now write (nd plot) the finl result. 4

p(x) = ( πx 0 (x/x0 ) ) x < x 0 = 0 x > x 0 As check of the result, note tht the re of the shded rectngle is equl to /π. The re is dimensionless, s it should be, nd is resonble frction of the nticipted totl re under p(x), tht is. c) The sketch of p(x) is shown bove. By inspection the most probble vlue of x is ±x 0 nd the lest probble ccessible vlue of x is zero. The men vlue of x is zero by symmetry. It is the divergence of p(x) t the turning points tht gives rise to the pprent imge of the pencil t these points in your experiment. COMMENTS If n oscilltor oscilltes bck nd forth with some fixed frequency, why is this p(x) independent of time? The reson is tht we did not know the strting time (or equivlently the phse φ) so we used n pproch which effectively verged over ll possible strting times. This wshed out the time dependence nd left time-independent probbility. If we hd known the phse, or equivlently the position nd velocity t some given time, then the process would hve been deterministic. In tht cse p(x) would be delt function centered t vlue of x which oscillted bck nd forth between x 0 nd +x 0. Those of you who hve lredy hd course in quntum mechnics my wnt to compre the clssicl result you found bove with the result for quntum hrmonic oscilltor in n energy eigenstte with high vlue of the quntum number n nd the sme totl energy. Will this probbility be time dependent? No. Recll why the energy eigensttes of potentil re lso clled sttionry sttes. 5

Problem 4: Quntized Angulr Momentum ) Using the expression for the normliztion of probbility density, long with expressions for the men nd the men squre, we cn write three seprte equtions relting the individul probbilities. p( ) + p(0) + p() = p( ) + 0 p(0) + p() = < L x > = 3 x p( ) + 0 p(0) + p() = < L > = 3 We now hve three simple liner equtions in three unknowns. The lst two cn be simplified nd solved for two of our unknowns. p( ) + p() = p() = 3 p( ) + p() = p( ) = 3 6 Substitute these results into the first eqution to find the lst unknown. + p(0) + = p(0) = 6 3 b) 6

Problem 5: A Coherent Stte of Quntum Hrmonic Oscilltor [ iωt i x αx Ψ( r, t) = (πx ) /4 cos 0 exp (αxx 0 sin t α ω ω x t x 0 ) 0 sin ω ) ( ( x 0 0 t ) ] ) First note tht the given wvefunction hs the form Ψ = exp[ib + c] = exp[ib] exp[c] where, b nd c re rel. Thus the squre of the mgnitude of the wvefunction is simply exp[c] nd finding the probbility density is not lgebriclly difficult. (x αx 0 cos ωt) p(x, t) = Ψ(x, t) = exp[ ] πx x 0 0 b) By inspection we see tht this is Gussin with time dependent men < x >= αx 0 cos ωt nd time independent stndrd devition σ = x 0. c) p(x, t) involves time independent pulse shpe, Gussin, whose center oscilltes hrmoniclly between αx 0 nd αx 0 with rdin frequency ω. t=/ T t=3/4 T t=/4 T t=0 -αx 0 0 αx 0 x Those lredy fmilir with quntum mechnics will recognize this s coherent stte of the hrmonic oscilltor, stte whose behvior is closest to the clssicl behvior. It is not n energy eigenstte since p(x) depends on t. It should be compred with clssicl hrmonic oscilltor with known phse φ nd the sme mximum excursion: x = αx 0 cos ωt. In this deterministic clssicl cse p(x, t) is given by p(x, t) = δ(x αx 0 cos ωt). The coherent stte is good representtion of the quntum behvior of the electromgnetic field of lser well bove the threshold for oscilltion. 7

Problem 6: Bose-Einstein Sttistics We re given the discrete probbility density p(n) = ( ) n n = 0,,, ) First we find the men of n. < n >= np(n) = ( ) n n } {{ S } The sum S cn be found by mnipulting the normliztion sum. p(n) = ( ) n = ( ) n must = Rerrnging the lst two terms gives the sum of geometric series: n =. But note wht hppens when we tke the derivtive of this result with respect to the prmeter. d n = n n = n n = S d d ( ) lso = = d ( ) Equting the two results gives the vlue of the sum we need, S = /( ), nd llows us to finish the computtion of the men of n: < n >=. c) To find the vrince we first need the men of the squre of n. < n >= n p(n) = ( ) n n }{{} S 8

Now try the sme trick used bove, but on the sum S. d d }{{} n n = n n = S n n = Then nd S d ( ) lso = = + d ( ) ( ) 3 ( ) < n > = [ ] ( ) ( ) ( ) + = + ( ) 3 ( ) = < n > + < n >, Vrince = < n > < n > =< n > + < n > = < n > (+ < n >). This is greter thn the vrince for Poisson, < n >, by fctor + < n >. c) p(x) = ( = f(x) δ(x ) n δ(x n) n) Try f(x) = Ce x/φ, then f(x = n) = Ce n/φ = C(e /φ ) n = ( ) n. This tells us tht C = nd exp( /φ) =. We cn invert the expression found bove for < n > to give s function of < n >: =< n > /(+ < n >). ( < n > ) /φ = ln = ln + < n > /φ = ( < n > + ) ( ) ln = ln + < n > < n > 9

Recll tht for smll x one hs the expnsion ln( + x) = x x / +.... Therefore in the limit < n > >>, /φ / < n > which implies φ < n >. 0

MIT OpenCourseWre http://ocw.mit.edu 8.044 Sttisticl Physics I Spring 03 For informtion bout citing these mterils or our Terms of Use, visit: http://ocw.mit.edu/terms.