π e ax2 dx = x 2 e ax2 dx or x 3 e ax2 dx = 1 x 4 e ax2 dx = 3 π 8a 5/2 (a) We are considering the Maxwell velocity distribution function: 2πτ/m

Similar documents
Probability Theory. The nth coefficient of the Taylor series of f(k), expanded around k = 0, gives the nth moment of x as ( ik) n n!

Week 9 Chapter 10 Section 1-5

10.40 Appendix Connection to Thermodynamics and Derivation of Boltzmann Distribution

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

Physics 607 Exam 1. ( ) = 1, Γ( z +1) = zγ( z) x n e x2 dx = 1. e x2

First Law: A body at rest remains at rest, a body in motion continues to move at constant velocity, unless acted upon by an external force.

Lecture 7: Boltzmann distribution & Thermodynamics of mixing

Physics 181. Particle Systems

Conservation of Angular Momentum = "Spin"

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

STATISTICAL MECHANICAL ENSEMBLES 1 MICROSCOPIC AND MACROSCOPIC VARIABLES PHASE SPACE ENSEMBLES. CHE 524 A. Panagiotopoulos 1

NUMERICAL DIFFERENTIATION

Structure and Drive Paul A. Jensen Copyright July 20, 2003

Lecture 4. Macrostates and Microstates (Ch. 2 )

1 Matrix representations of canonical matrices

Prof. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model

SUPPLEMENTARY INFORMATION

Georgia Tech PHYS 6124 Mathematical Methods of Physics I

8.592J: Solutions for Assignment 7 Spring 2005

Physics 207: Lecture 20. Today s Agenda Homework for Monday

Lecture 12: Discrete Laplacian

10. Canonical Transformations Michael Fowler

University of Washington Department of Chemistry Chemistry 453 Winter Quarter 2015

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential

EPR Paradox and the Physical Meaning of an Experiment in Quantum Mechanics. Vesselin C. Noninski

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

8.1 Arc Length. What is the length of a curve? How can we approximate it? We could do it following the pattern we ve used before

Module 1 : The equation of continuity. Lecture 1: Equation of Continuity

12. The Hamilton-Jacobi Equation Michael Fowler

PHYS 705: Classical Mechanics. Calculus of Variations II

Goodness of fit and Wilks theorem

Module 2. Random Processes. Version 2 ECE IIT, Kharagpur

CIVL 8/7117 Chapter 10 - Isoparametric Formulation 42/56

Section 8.3 Polar Form of Complex Numbers

Complex Numbers Alpha, Round 1 Test #123

STATISTICAL MECHANICS

Physics for Scientists and Engineers. Chapter 9 Impulse and Momentum

Poisson brackets and canonical transformations

1 (1 + ( )) = 1 8 ( ) = (c) Carrying out the Taylor expansion, in this case, the series truncates at second order:

Celestial Mechanics. Basic Orbits. Why circles? Tycho Brahe. PHY celestial-mechanics - J. Hedberg

MTH 263 Practice Test #1 Spring 1999

x = , so that calculated

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

Math1110 (Spring 2009) Prelim 3 - Solutions

Week 8: Chapter 9. Linear Momentum. Newton Law and Momentum. Linear Momentum, cont. Conservation of Linear Momentum. Conservation of Momentum, 2

PHYS 705: Classical Mechanics. Newtonian Mechanics

The Feynman path integral

Assignment 2. Tyler Shendruk February 19, 2010

Solution Thermodynamics

Workshop: Approximating energies and wave functions Quantum aspects of physical chemistry

10.34 Fall 2015 Metropolis Monte Carlo Algorithm

Physics 111: Mechanics Lecture 11

1 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations

Module 3: Element Properties Lecture 1: Natural Coordinates

So far: simple (planar) geometries

Lagrange Multipliers. A Somewhat Silly Example. Monday, 25 September 2013

Thermodynamics and Kinetics of Solids 33. III. Statistical Thermodynamics. Â N i = N (5.3) N i. i =0. Â e i = E (5.4) has a maximum.

U.C. Berkeley CS294: Spectral Methods and Expanders Handout 8 Luca Trevisan February 17, 2016

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

Strong Markov property: Same assertion holds for stopping times τ.

5.62 Physical Chemistry II Spring 2008

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

Thermodynamics and statistical mechanics in materials modelling II

Moments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is.

and Statistical Mechanics Material Properties

More metrics on cartesian products

Physics 240: Worksheet 30 Name:

TP A SOLUTION. For an ideal monatomic gas U=3/2nRT, Since the process is at constant pressure Q = C. giving ) =1000/(5/2*8.31*10)

Modelli Clamfim Equazioni differenziali 7 ottobre 2013

Copyright 2017 by Taylor Enterprises, Inc., All Rights Reserved. Adjusted Control Limits for P Charts. Dr. Wayne A. Taylor

Chapter 1. Probability

Problem Set 9 Solutions

Inductance Calculation for Conductors of Arbitrary Shape

Math 217 Fall 2013 Homework 2 Solutions

One Dimension Again. Chapter Fourteen

PROBLEM SET 7 GENERAL EQUILIBRIUM

Lecture 5: Ideal monatomic gas

How Strong Are Weak Patents? Joseph Farrell and Carl Shapiro. Supplementary Material Licensing Probabilistic Patents to Cournot Oligopolists *

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family

Chapter 11 Angular Momentum

Hidden Markov Models & The Multivariate Gaussian (10/26/04)

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1

Homework & Solution. Contributors. Prof. Lee, Hyun Min. Particle Physics Winter School. Park, Ye

MATH 241B FUNCTIONAL ANALYSIS - NOTES EXAMPLES OF C ALGEBRAS

Open string operator quantization

Statistical mechanics handout 4

LECTURE 9 CANONICAL CORRELATION ANALYSIS

Conjugacy and the Exponential Family

Fermi-Dirac statistics

PHY688, Statistical Mechanics

Limited Dependent Variables

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)

EEE 241: Linear Systems

8.323 Relativistic Quantum Field Theory I

Bézier curves. Michael S. Floater. September 10, These notes provide an introduction to Bézier curves. i=0

2.3 Nilpotent endomorphisms

A particle in a state of uniform motion remain in that state of motion unless acted upon by external force.

CHAPTER 6. LAGRANGE S EQUATIONS (Analytical Mechanics)

8 Derivation of Network Rate Equations from Single- Cell Conductance Equations

Formal solvers of the RT equation

Transcription:

Homework Solutons Problem In solvng ths problem, we wll need to calculate some moments of the Gaussan dstrbuton. The brute-force method s to ntegrate by parts but there s a nce trck. The followng ntegrals should be known: e ax dx = π a xe ax dx = a Suppose we now want to calculate x e ax dx or x 3 e ax dx. We can evaluate these ntegrals by dfferentatng the above expressons. For example: d e ax dx = x e ax dx = d π π = da da a a 3/ d xe ax dx = x 3 e ax dx = d = da da a a d x e ax dx = x 4 e ax dx = d π π 3 = da da a 3/ 4 a 5/ Whch panlessly mples: π x e ax dx = a 3/ x 3 e ax dx = a x 4 e ax dx = 3 π 8a 5/ a We are consderng the Maxwell velocty dstrbuton functon: 3 pvdv = 4π e mv /τ v dv πτ/m where C = 4π πτ/m 3 = Ce av v dv and a = m τ. The most probable speed, v prob, s the speed whch maxmzes pv. dp = constant ve av + v e av av = dv τ v prob = a = m The average speed s calculated usng our Gaussan ntegrals: < v > = vpvdv = C v 3 e av dv = C a = 4π πτ/m 8τ < v >= πm 3 τ m

Smlarly: < v > = v pvdv = C v 4 e av dv = C 3 π 8a = 4π 5/ πτ/m v rms = 3τ < v > = m 3 3 π 8 b For N, the molecular weght s approxmately m N = 8m p where m p =.67 7 kg. Puttng n the numbers gves: τ m 5 = 3τ m v prob = τ m = [.38 3 J K 93K 8.67 7 kg ] / = 45 m s < v > = 468 m s v rms = 57 m s Problem V partcle wth velocty v,v,v x y z vdt θ A dv Fgure : The fgure shows the system at tme t. Partcles wth velocty v x, v y, v z must be n volume dv n order to ht area A by the tme t + dt. See fgure. We frst determne the contrbuton of those partcles havng velocty v = v x, v y, v z known henceforth as v-partcles to the flux. If our hole has area A, then n a tme nterval dt, only those v-partcles n the volume dv wll leave the hole. We may calculate dv : dv = Adtv cosθ where v = v. The number of v-partcles n the volume dv, whch we call nvdv, s gven by the number of v-partcles n the entre volume tmes the rato dv/v. The number of v-partcles actually the number of partcles wth veloctes between v and v + dv n the entre volume V s proportonal to the Maxwell velocty dstrbuton functon. nvdv = C dv V exp mv dv τ = C Adtexp mv v 3 dv cosθ snθdθdφ τ

where C and C are constants. The prevous expresson gves the number of v-partcles leavng the hole of area A n tme dt. The flux of v-partcles, whch we call Fv, s defned as ths expresson dvded by Adt: Fvdv = C exp mv v 3 dv cosθ sn θdθdφ τ So we have calculated the velocty dstrbuton of the flux. In order to get a speed dstrbuton, we must ntegrate over sold angle. Fvdv = π dφ π/ dθfv = D exp mv v 3 dv τ where D s a constant. Note that the θ ntegraton was only over [, π ] because the range [ π, π] corresponds to partcles movng away from the hole. We can get D by requrng our functon to be normalzed. Fvdv = D v 3 exp mv dv = D τ = τ m m D = τ Problem 3 A sngle spn can be n one of two energy states: +E and E. The sngle partcle partton functon s gven by: Zτ = e E/τ + e E/τ The partton functon for N partcles s a product of the sngle partcle functons: Z N τ = Zτ N = e E/τ + e E/τ N In partcular, there s no factor of N! n the denomnator because the spns are dstngushable, beng labeled by ther lattce postons. If the spns were not fxed n space but nstead could move from ste to ste, then we would need the N!. Whle ths does not effect the energy, t wll effect the entropy. The free energy s obtaned n the usual way: and the entropy: Sτ and the energy: Fτ = τ lnzτ = Nτ lne E/τ + e E/τ F = τ V = N lne E/τ + e E/τ τ E τ e E/τ + E τ e E/τ e E/τ + e E/τ = N lne E/τ + e E/τ E τ tanh E τ U = F + τs = NE tanh E τ 3

Problem 4 Suppose we have two systems whch nteract suffcently weakly that ther nteracton energy may be gnored. The fact that they nteract at all mples they should be at the same temperature. A state of the combned system s specfed by gvng the sngle system states of ndvdual systems. The combned partton functon s: Z tot = e E+E/τ states of,states of = e E/τ e E/τ states of states of = Z Z Problem 5 a The states of our system are unquely specfed by the number of open lnks. The state havng n open lnks corresponds to an energy of nɛ. Thus: N Z = e ɛ/τn = n= exp N + ɛ/τ exp ɛ/τ b Let x = ɛ/τ. In the lmt x >> : N n= ne ɛ/τn = d dx Z = e N+x N + e x e N+x e x + e x e x e x e x Problem 6 a Ths problem s the same as havng a one-dmensonal chan of spns. The spn on each lattce ste can pont left or rght. If s s the rght excess, then: N gn, s = N + s s the number of confguratons wth that rght excess. Because the energy of the system only depends on the magntude of s, we can fnd the total number of states wth a gven s : gn, s + gn, s = N! N + s! N s! b To fnd the entropy, we make use of the Strlng approxmaton and also the fact that s << N here we wrte s but mean s : σl = lngn, s + gn, s = ln N! ln N + s! ln N s! 4

N lnn N = gn, N = gn, N = gn, N gn, s N = gn, l Nρ s + N ln N s + N N s N ln N s + N s s + ln + N N N s s ln N N s + N s N s N N s N s N s N s + N s N s N + N s N s N + s N where n the last lne, we used the defnton l = sρ. c σ f = τ l = lτ U Nρ Problem 7 We model a sngle partcle as a free partcle n a one-dmensonal box. The energy levels of ths system are gven by: E n = h π ML n where n s a postve nteger. The partton functon s gven by: Z = e En/τ = dne α n π Mτ = n= α = π h L h where α = π MτL. The N partcle partton functon s gven by: Z N = ZN N! Here the N! s needed because the partcles are ndstngushable. We get the free energy by the usual method: Mτ F = τ lnz N = τ N ln π h L N lnn + N = τ N ln And fnally, the entropy: σ F = τ = N ln Mτ π h L N + N Mτ = N ln L L π h N + N + N Mτ L π h N + 3 5

Problem 8 The probabltes whch maxmzed the entropy n the canoncal ensemble were gven by p = e λ λe. We note that the probabltes must sum to by defnton: p = e λ e λe = e λ = = e λe Zλ p = e λe Zλ where the notaton Zλ dentfes the sum wth a partton functon. The average energy and entropy of ths system are gven by: U = p E σ = p lnp The parameter λ may be vared, as t s somethng mposed on the system from the outsde to mantan the average energy constrant. Varyng λ mples the followng varatons n the energy and entropy: δu = E δp δσ = δp lnp + δp Because the p s are probabltes, they must always sum to meanng that the sum of ther varatons must gve : δp =. Ths mples: δσ = = δp lnp + δp δp lnp = δp λ E lnzλ = λ E δp Fnally, we use the statstcal defnton of temperature: τ = σ U V to conclude that λ = τ. 6