Partition Function at a Fixed Pressure

MathematicalPhysicist
Science Advisor
Gold Member
Messages
4,662
Reaction score
372

Homework Statement


I don't quite follow the solution to this problem (problem 2.11 in Bergersen's and Plischke's textbook), here are the quoted problem and its solution:

problem:
Consider a system of ##N## noninteracting molecules in a container of
cross-sectional area ##A##. The bottom of the container (at ##z = 0##) is rigid.
The top consists of an airtight piston of mass ##M## which slides without
friction.
(a) Construct the partition function ##Z## of the ##(N + 1)##-particle system
(##N## molecules of mass ##m##, one piston of mass ##M##, cross-sectional area
##A##). You may neglect the effect of gravity on the gas molecules.
(b) Show that the thermodynamic potential ##—k_BT\ln Z## is, in the ther-
modynamic limit, identical to the Gibbs potential of an ideal gas of
##N## molecules, subject to the pressure ##P = Mg/A##.

solution:
(a) The Hamiltonian for the system consisting of ##N## particles plus the frictionless piston of mass ##M## is $$H = \sum_{i=1}^N p_i^2/(2m)+P_z^2/(2M) +Mgz$$

The partition function for a single gas molecule in a volume ##V=Az## is ##Z_1 = \frac{Az}{\lambda^3}## where $$\lambda = \sqrt{\frac{h^2}{2\pi m k_B T}}$$
The partition function for the complete system is then:
$$ Z = \frac{A^N}{N!\lambda^{3N}} \int_{-\infty}^\infty \frac{dP_z}{h}e^{\frac{-\beta P_z^2}{2M}}\int_0^\infty dz z^N e^{-\beta Mgz}$$
or
$$Z = \frac{A^N(\beta Mgz)^{N+1}}{\lambda^{3N}} \sqrt{\frac{2\pi M}{\beta h^2}}$$
We find $$(2.10) \ \ \ \ -k_B T \ln Z = - Nk_B T \ln \bigg(\frac{Ak_B T}{\lambda^3 Mg}\bigg)$$

(b) It was shown in the text that the chemical potential is given by $$ (2.11) \ \ \ \ \mu = k_B T \ln \bigg( \frac{N\lambda^3}{V} \bigg) = G/N$$

Identifying the pressure as ##Mg/A## and using the ideal gas law ##PV=Nk_B T## we see that ##(2.10)## and ##(2.11)## are equivalent.

Homework Equations

The Attempt at a Solution


My problem is with the solution to (a), it seems they plugged into the LHS of (2.10) ##Z= Z_1^N## where ##k_B T = Mgz## and not the expression ##Z = \frac{A^N(\beta Mgz)^{N+1}}{\lambda^{3N}} \sqrt{\frac{2\pi M}{\beta h^2}}##, are they equivalent?

It doesn't look like that? what do you think?
 
Physics news on Phys.org
MathematicalPhysicist said:
My problem is with the solution to (a), it seems they plugged into the LHS of (2.10) ##Z= Z_1^N## where ##k_B T = Mgz## and not the expression ##Z = \frac{A^N(\beta Mgz)^{N+1}}{\lambda^{3N}} \sqrt{\frac{2\pi M}{\beta h^2}}##, are they equivalent?

I don't understand the reason for the ##z^N## factor in the integrand below
upload_2017-10-16_14-0-22.png

[EDIT: Never mind, I see where the ##z^N## is coming from.]

With this expression for ##Z##, then the next equation should read
upload_2017-10-16_14-8-41.png
where the exponent (N+1) should be -(N+1). Also, I don't think the ##z## in ##(\beta Mgz)## should be there. Of course, you should check this.

Then you get (2.10) if you assume N is very large and neglect certain small terms.
 
Last edited:
  • Like
Likes MathematicalPhysicist
Thread 'Need help understanding this figure on energy levels'
This figure is from "Introduction to Quantum Mechanics" by Griffiths (3rd edition). It is available to download. It is from page 142. I am hoping the usual people on this site will give me a hand understanding what is going on in the figure. After the equation (4.50) it says "It is customary to introduce the principal quantum number, ##n##, which simply orders the allowed energies, starting with 1 for the ground state. (see the figure)" I still don't understand the figure :( Here is...
Thread 'Understanding how to "tack on" the time wiggle factor'
The last problem I posted on QM made it into advanced homework help, that is why I am putting it here. I am sorry for any hassle imposed on the moderators by myself. Part (a) is quite easy. We get $$\sigma_1 = 2\lambda, \mathbf{v}_1 = \begin{pmatrix} 0 \\ 0 \\ 1 \end{pmatrix} \sigma_2 = \lambda, \mathbf{v}_2 = \begin{pmatrix} 1/\sqrt{2} \\ 1/\sqrt{2} \\ 0 \end{pmatrix} \sigma_3 = -\lambda, \mathbf{v}_3 = \begin{pmatrix} 1/\sqrt{2} \\ -1/\sqrt{2} \\ 0 \end{pmatrix} $$ There are two ways...
Back
Top