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 />

However, P cannot not be easily solved on many conditions. We would use Markov<br />

Process <strong>in</strong>stead to solve this problem.<br />

2.Markov Process<br />

Markov process is an example of 'random walk' through phase space <strong>in</strong> a series of .<br />

We would first give the three properties of Markov cha<strong>in</strong> below:<br />

[1]<br />

[2]<br />

[3]f(x)<br />

means the probability that a po<strong>in</strong>t transfer from site x to site x'. [1] is<br />

obvious. [2] says the random walk is not restrict to history, the sites can be visited many<br />

times, and we will f<strong>in</strong>ally reach to the border or equilibrium state. We call this<br />

"ergodicity hypothesis". [3] is just the learned "detailed balance" where the equilibrium<br />

distribution f(x) is reached.<br />

How to efficiently converge <strong>in</strong> the Marcov Cha<strong>in</strong> depends on the selection of ,<br />

which has a large freedom to choose. A simple choice is seen below, also known as the<br />

Metropolis method.<br />

How to pursue the random walk? We use (1.2.1) as an example and use Metropolis<br />

method. We def<strong>in</strong>e the transition probability as<br />

Follow<strong>in</strong>g Metropolis, we choose another transition probability p, which reads<br />

or<br />

Then follow the follow<strong>in</strong>g steps:<br />

[1] Select a site i, and its first random step to the nearest neighbour i'<br />

[2] Compute<br />

[3] Calculate , if , then p=1, accept the sp<strong>in</strong> value of this site; otherwise,<br />

go to [4]

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

Saved successfully!

Ooh no, something went wrong!