Warning: foreach() argument must be of type array|object, bool given in /var/www/html/web/app/themes/studypress-core-theme/template-parts/header/mobile-offcanvas.php on line 20

In this section I've formulated the cluster expansion for a gas with a fixed number of particles, using the "canonical" formalism of Chapter 6. A somewhat cleaner approach, however, is to use the "grand canonical" formalism introduced in Section 7.1, in which we allow the system to exchange particles with a much larger reservoir.

(a) Write down a formula for the grand partition function (Z) of a weakly interacting gas in thermal and diffusive equilibrium with a reservoir at fixed T andµ. Express Z as a sum over all possible particle numbers N, with each term involving the ordinary partition function Z(N).

(b) Use equations 8.6 and 8.20 to express Z(N) as a sum of diagrams, then carry out the sum over N, diagram by diagram. Express the result as a sum of similar diagrams, but with a new rule 1 that associates the expression (>./vQ) J d3ri with each dot, where >. = e13µ,. Now, with the awkward factors of N(N - 1) · · · taken care of, you should find that the sum of all diagrams organizes itself into exponential form, resulting in the formula

Note that the exponent contains all connected diagrams, including those that can be disconnected by removal of a single line.

(c) Using the properties of the grand partition function (see Problem 7.7), find diagrammatic expressions for the average number of particles and the pressure of this gas.

(d) Keeping only the first diagram in each sum, express N(µ) and P(µ) in terms of an integral of the Mayer /-function. Eliminate µ to obtain the same result for the pressure (and the second virial coefficient) as derived in the text.

(e) Repeat part (d) keeping the three-dot diagrams as well, to obtain an expression for the third virial coefficient in terms of an integral of /-functions. You should find that the A-shaped diagram cancels, leaving only the triangle diagram to contribute to C(T).

Short Answer

Expert verified

The required equation isZ=1N!·(2πmkT)3·d3r1d3rN·e-βUpot.

Step by step solution

01

Given information

We have been given that Z1=1h3·d3rd3p·e-βU

02

Simplify

For same molecules

Z1=1h3·d3rd3p·e-βU

Energy can be written as:

U=p122m+p222m++pN22m

Z=1N!·1h3N·d3r1d3rNd3p1d3pN·e-β·p122me-β·pN22m·e-βUppt

We can write virial equation

=2π·0r2e-βU(r)-1dr

Unlock Step-by-Step Solutions & Ace Your Exams!

  • Full Textbook Solutions

    Get detailed explanations and key concepts

  • Unlimited Al creation

    Al flashcards, explanations, exams and more...

  • Ads-free access

    To over 500 millions flashcards

  • Money-back guarantee

    We refund you if you fail your exam.

Over 30 million students worldwide already upgrade their learning with Vaia!

One App. One Place for Learning.

All the tools & learning materials you need for study success - in one app.

Get started for free

Most popular questions from this chapter

By changing variables as in the text, express the diagram in equation 8.18 in terms of the same integral as in the equation8.31. Do the same for the last two diagrams in the first line of the equation8.20. Which diagrams cannot be written in terms of this basic integral?

Modify the ising program to simulate a one-dimensional Ising model.

(a) For a lattice size of 100, observe the sequence of states generated at various temperatures and discuss the results. According to the exact solution (for an infinite lattice), we expect this system to magnetise only as the temperature goes to zero; is the behaviour of your program consistent with this prediction? How does the typical cluster size depend on temperature?

(b) Modify your program to compute the average energy as in Problem 8.27. Plot the energy and heat capacity vs. temperature and compare to the exact result for an infinite lattice.

(c) Modify your program to compute the magnetisation as in Problem 8.28. Determine the most likely magnetisation for various temperatures and sketch a graph of this quantity. Discuss.

In this problem you will use the mean field approximation to analyse the behaviour of the Ising model near the critical point.

(a) Prove that, when x1,tanhxx-13x3

(b) Use the result of part (a) to find an expression for the magnetisation of the Ising model, in the mean field approximation, when T is very close to the critical temperature. You should find MTc-Tβ¯,whereβ(not to be confused with 1/kT) is a critical exponent, analogous to the f defined for a fluid in Problem 5.55. Onsager's exact solution shows that β=1/8in two dimensions, while experiments and more sophisticated approximations show that β1/3in three dimensions. The mean field approximation, however, predicts a larger value.

(c) The magnetic susceptibility χis defined as χ(M/B)T. The behaviour of this quantity near the critical point is conventionally written as χT-Tc-γ , where y is another critical exponent. Find the value of in the mean field approximation, and show that it does not depend on whether T is slightly above or slightly below Te. (The exact value of y in two dimensions turns out to be 7/4, while in three dimensions γ1.24.)

Modify the ising program to compute the average energy of the system over all iterations. To do this, first add code to the initialise subroutine compute the initial energy of the lattice; then, whenever a dipole is flipped, change the energy variable by the appropriate amount. When computing the average energy, be sure to average over all iterations, not just those iterations in which a dipole is actually flipped (why?). Run the program for a 5 x 5 lattice for T values from 4 down to l in reasonably small intervals, then plot the average energy as a function of T. Also plot the heat capacity. Use at least 1000 iterations per dipole for each run, preferably more. If your computer is fast enough, repeat for a 10x 10 lattice and for a 20 x 20 lattice. Discuss the results. (Hint: Rather than starting over at each temperature with a random initial state, you can save time by starting with the final state generated at the previous, nearby temperature. For the larger lattices you may wish to save time by considering only a smaller temperature interval, perhaps from 3 down to 1.5.)

To quantify the clustering of alignments within an Ising magnet, we define a quantity called the correlation function, c(r). Take any two dipoles i and j, separated by a distance r, and compute the product of their states: sisj. This product is 1 if the dipoles are parallel and -1 if the dipoles are antiparallel. Now average this quantity over all pairs that are separated by a fixed distance r, to |obtain a measure of the tendency of dipoles to be "correlated" over this distance. Finally, to remove the effect of any overall magnetisation of the system, subtract off the square of the average s. Written as an equation, then, the correlation function is

c(r)=sisj¯-si¯2

where it is understood that the first term averages over all pairs at the fixed distance r. Technically, the averages should also be taken over all possible states of the system, but don't do this yet.

(a) Add a routine to the ising program to compute the correlation function for the current state of the lattice, averaging over all pairs separated either vertically or horizontally (but not diagonally) by r units of distance, where r varies from 1 to half the lattice size. Have the program execute this routine periodically and plot the results as a bar graph.

(b) Run this program at a variety of temperatures, above, below, and near the critical point. Use a lattice size of at least 20, preferably larger (especially near the critical point). Describe the behaviour of the correlation function at each temperature.

(c) Now add code to compute the average correlation function over the duration of a run. (However, it's best to let the system "equilibrate" to a typical state before you begin accumulating averages.) The correlation length is defined as the distance over which the correlation function decreases by a factor of e. Estimate the correlation length at each temperature, and plot graph of the correlation length vs.

See all solutions

Recommended explanations on Physics Textbooks

View all explanations

What do you think about this solution?

We value your feedback to improve our textbook solutions.

Study anywhere. Anytime. Across all devices.

Sign-up for free