/*! 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.26 Implement the ising program on y... [FREE SOLUTION] | 91Ó°ÊÓ

91Ó°ÊÓ

Implement the ising program on your favourite computer, using your favourite programming language. Run it for various lattice sizes and temperatures and observe the results. In particular:

(a) Run the program with a 20 x 20 lattice at T = 10, 5, 4, 3, and 2.5, for at least 100 iterations per dipole per run. At each temperature make a rough estimate of the size of the largest clusters.

(b) Repeat part (a) for a 40 x 40 lattice. Are the cluster sizes any different? Explain. (c) Run the program with a 20 x 20 lattice at T = 2, 1.5, and 1. Estimate the average magnetisation (as a percentage of total saturation) at each of these temperatures. Disregard runs in which the system gets stuck in a metastable state with two domains.

(d) Run the program with a 10x 10 lattice at T = 2.5. Watch it run for 100,000 iterations or so. Describe and explain the behaviour.

(e) Use successively larger lattices to estimate the typical cluster size at temperatures from 2.5 down to 2.27 (the critical temperature). The closer you are to the critical temperature, the larger a lattice you'll need and the longer the program will have to run. Quit when you realise that there are better ways to spend your time. Is it plausible that the cluster size goes to infinity as the temperature approaches the critical temperature?

Short Answer

Expert verified

The pictures and codes are given.

Step by step solution

01

Given information and explanation

Implement the ising program on your favourite computer, using your favourite programming language. Run it for various lattice sizes and temperatures and observe the results.

A Python application is written, and this is how it looks:

02

Explanation

a)We can observe that when the temperature rises, the size of the domain shrinks.

b)We can see that the clusters are smaller for the larger lattice.

03

Explanation

c) We have the following scenarios for these three temperatures. The average magnetisation (number of black squares) is around 50% or higher, as can be shown.

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

Use the cluster expansion to write the total energy of a monatomic nonideal gas in terms of a sum of diagrams. Keeping only the first diagram, show that the energy is approximatelyU≈32NkT+N2V·2π∫0∞r2u(r)e-βu(r)drUse a computer to evaluate this integral numerically, as a function of T, for the Lennard-Jones potential. Plot the temperature-dependent part of the correction term, and explain the shape of the graph physically. Discuss the correction to the heat capacity at constant volume, and compute this correction numerically for argon at room temperature and atmospheric pressure.

Consider a gas of molecules whose interaction energy u(r)u is infinite for r<r0and negative for r>r0, with a minimum value of -u0. Suppose further that kT≫u0, so you can approximate the Boltzmann factor forr>r0using ex≈1+x. Show that under these conditions the second virial coefficient has the form B(T)=b-(a/kT), the same as what you found for a van der Waals gas in Problem 1.17. Write the van der Waals constants aand b in terms of r0and u(r), and discuss the results 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?

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?

Problem 8.13. Use the cluster expansion to write the total energy of a monatomic nonideal gas in terms of a sum of diagrams. Keeping only the first diagram, show that the energy is approximately
U≈32NkT+N2V·2π∫0∞r2u(r)e-βu(r)dr
Use a computer to evaluate this integral numerically, as a function of T, for the Lennard-Jones potential. Plot the temperature-dependent part of the correction term, and explain the shape of the graph physically. Discuss the correction to the heat capacity at constant volume, and compute this correction numerically for argon at room temperature and atmospheric pressure.

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.