23.05.2014 Views

Topics in Statistic Mechanics

Topics in Statistic Mechanics

Topics in Statistic Mechanics

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Adv. Sta. Phy. Homework 8 Markovian-cha<strong>in</strong> Monte Carlo Li,Zimeng PB06203182<br />

Markovian-cha<strong>in</strong> Monte Carlo<br />

Edited by Li, Zimeng<br />

Contents:<br />

1.Introduction<br />

1.1 Simple Sampl<strong>in</strong>g<br />

1.2 Importance Sampl<strong>in</strong>g<br />

2.Markov Process<br />

3.MCMC application on 2D Is<strong>in</strong>g model<br />

First, we have learned the MCMC method <strong>in</strong> Computational Physics, therefore I would<br />

only give a brief review of the method and, second, an extension to the 2D Is<strong>in</strong>g model.<br />

1.Introduction<br />

Partition functions and functional <strong>in</strong>tegrals (path <strong>in</strong>tegrals) reduce many-body problems<br />

to complicated multidimensional <strong>in</strong>tegrals or sums.<br />

The Monte Carlo technique has its orig<strong>in</strong>s <strong>in</strong> the numerical evaluation of <strong>in</strong>tegrals. The<br />

technique was later generalized to calculate the partition function and the mean value<br />

of observables <strong>in</strong> classical systems.<br />

The key of numerical <strong>in</strong>tegration is sampl<strong>in</strong>g. There are many sampl<strong>in</strong>g methods <strong>in</strong><br />

Monte-Carlo method. We will only <strong>in</strong>troduce two of them here.<br />

1.1 Simple Sampl<strong>in</strong>g<br />

We first generate random numbers<br />

probability density of these po<strong>in</strong>ts:<br />

uniformly distributed <strong>in</strong> [a,b], and get the<br />

We can calculate the expectation value of a function f with respect to the distribution P:<br />

E(f(x))=<br />

In order to calculate the above <strong>in</strong>tegration I <strong>in</strong> a numerical way, we def<strong>in</strong>e the Monte-<br />

Carlo Integral

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!