No Title next up previous
Next: About this document

G25.2651: Statistical Mechanics

Problem set #8

Due: May 5, 2003

1.
The quantum time correlation function corresponding to an absorption process is

displaymath15

where Q is the canonical partition function. Consider a more ``symmetric'' correlation function

displaymath23

where tex2html_wrap_inline186 is a complex time variable. Show that the Fourier transforms tex2html_wrap_inline188 and tex2html_wrap_inline190 and give a relationship between tex2html_wrap_inline188 and tex2html_wrap_inline190 . Finally, give an expression for the absorption probability spectrum tex2html_wrap_inline196 in terms of tex2html_wrap_inline190

2.
Consider the free rotational motion of a rigid heteronuclear diatomic molecule of (fixed) bond length R and moment of inertia tex2html_wrap_inline202 , where tex2html_wrap_inline204 is the reduced mass, about an axis through its center of mass perpendicular to the internuclear bond axis. The molecule is constrained to rotate in the x-y plane only. One of the atoms carries a charge q and the other a charge -q.

a.
Ignoring center of mass motion, write down the Hamiltonian, tex2html_wrap_inline214 for the molecule.

b.
Find the eigenvalues and eigenvectors of tex2html_wrap_inline214 .

c.
If the molecule is exposed to spatially homogeneous, monochromatic radiation with an electric field tex2html_wrap_inline218 given by

displaymath46

write down the perturbation Hamiltonian tex2html_wrap_inline220 .

d.
Calculate the energy spectrum tex2html_wrap_inline222 for tex2html_wrap_inline224 . Interpret your results, and in particular, explain how the allowed absorptions and emissions are manifest in your final expression. Plot the absorption part of your spectrum. Where do you expect the peak intensity to occur?

Hint: Consider using a convergence factor, tex2html_wrap_inline226 , and let tex2html_wrap_inline228 to go 0 at the end of the calculation.

e.
Based on your results from parts (a)-(d), plot the spectrum 3-dimensional rigid rotor, for which the energy eigenvalues are tex2html_wrap_inline230 and m=-l,...,l is the quantum number for the z-component of angular momentum. Where do you expect the peak intensity to occur in the 3-dimensional case?

3.
In order to describe a system in a non-cubic box, such as might occur under the action of a directional shearing force, one often introduces a box matrix tex2html_wrap_inline236 , whose columns contain the three cell vectors tex2html_wrap_inline238 , tex2html_wrap_inline240 and tex2html_wrap_inline242 that describe a parallelpiped:

   figure64
Figure 1:

a.
We know from elementary geometry that the volume of the box is given by tex2html_wrap_inline244 . Show that the volume is also given by

displaymath76

b.
In the canonical ensemble, tex2html_wrap_inline246 , with a Cartesian Hamiltonian

displaymath81

In order to make the dependence on the box matrix explicit, introduce the transformation

displaymath87

where tex2html_wrap_inline248 is a momentum conjugate to the scaled coordinate tex2html_wrap_inline250 . Here, tex2html_wrap_inline252 and tex2html_wrap_inline254 index the three spatial directions, x, y, and z. Show that this transformation preserves the phase space measure:

displaymath99

Show, therefore, that the pressure tensor estimator discussed in class can be deduced from the thermodynamic relation

displaymath109

where tex2html_wrap_inline262 are indices that run over the three cartesian directions x,y,z.

c.
The isothermal-isobaric (NPT) ensemble is particularly useful for determining the bulk viscosity of a substance. Recall the partition function of this ensemble is

displaymath119

where tex2html_wrap_inline268 is the external applied pressure.

d.
Show that the linear response formula is preserved is tex2html_wrap_inline270 is chosen to the the isothermal-isobaric distribution function tex2html_wrap_inline272 .
e.
Next, consider coupling a system to an external compression field described by the equations of motion

eqnarray131

where tex2html_wrap_inline252 and tex2html_wrap_inline254 index the three spatial directions, x, y, and z. Show that the equations of motion satisfy the incompressibility condition.

f.
Consider the specific choice:

displaymath146

where tex2html_wrap_inline284 is the compression rate. The coefficient of bulk viscosity, tex2html_wrap_inline286 is given by a generalization of Newton's law of viscosity:

displaymath152

where tex2html_wrap_inline288 means average over the equilibrium NPT distribution function and tex2html_wrap_inline292 is the full nonequilibrium average. Using the linear response formula to evaluate tex2html_wrap_inline294 , derive the appropriate Green-Kubo expression for tex2html_wrap_inline286 .




next up previous
Next: About this document

Mark Tuckerman
Mon Apr 28 14:37:27 EDT 2003