/*! This file is auto-generated */ .wp-block-button__link{color:#fff;background-color:#32373c;border-radius:9999px;box-shadow:none;text-decoration:none;padding:calc(.667em + 2px) calc(1.333em + 2px);font-size:1.125em}.wp-block-file__button{background:#32373c;color:#fff;text-decoration:none} Q. 8.9 Show that the Lennard-Jones pote... [FREE SOLUTION] | 91Ó°ÊÓ

91Ó°ÊÓ

Show that the Lennard-Jones potential reaches its minimum value at r=r0, and that its value at this minimum is -u0. At what value of rdoes the potential equal zero?

Short Answer

Expert verified

At, r=r0216the value of Lennard-Jones potential become 0.

Step by step solution

01

Step 1. Given information

Lennard- Jones potential=

u(r)=u0r0r12-2r0r6

02

Step 2. Differentiating both sides with respect to r

We get,

dudr=u0r012×-121r13-2r06×-61r7

=u0-12r012r13+12r06r7

For extreme points,

dudr=0

∴u0-12r012r13+12r06r7=0

12r06u0r13r6-r06=0

r6-r06=0

r=r0

Thus the Lennard-Jones potential has extremum point at r=r0

03

Step 3. To check that the potential at r=r0 is maximum or minimum take double derivative

We get,

d2udr2=u0-12×-13r012r14+12×-7r06r8

=12u013r012r14-7r06r8

d2udr2r=s0=12u013r012r014-7r06r08

=12u013r02-7r02

=12u0r02[6]

=72u0r02>0

∴d2udr2r=r0>0

So the Lennard- Jones potential has minimum value at r'=r0

u(r)r=r0=u0r0r12-2r0r6r=r0

=u0r0r12-2r0r6

=u0[1-2]

u(r)r=r0=-u0

04

Step 4. If the potential becomes 0.

u(r)=0

u0r0*r12-2r0r6=0

u0r12r012-2r06r6=0

r012-2r06r6=0

r6=r0122r06=r062

r6=r0122r06

r=1216r0

Thus the Lennard-Jones potential becomes zero atr=1216r0

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 91Ó°ÊÓ!

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

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.)

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).

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.

Consider a gas of "hard spheres," which do not interact at all unless their separation distance is less than r0, in which case their interaction energy is infinite. Sketch the Mayer f-function for this gas, and compute the second virial coefficient. Discuss the result briefly.

Consider an Ising model of just two elementary dipoles, whose mutual interaction energy is ±ϵ. Enumerate the states of this system and write down their Boltzmann factors. Calculate the partition function. Find the probabilities of finding the dipoles parallel and antiparallel, and plot these probabilities as a function of kT/ϵ. Also calculate and plot the average energy of the system. At what temperatures are you more likely to find both dipoles pointing up than to find one up and one down?

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.