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

where (1.1.1)<br />

Due to the Central Limit Theorem, we have<br />

E[ ]=<br />

Therefore we can use (1.1.1) to calculate the Integral<br />

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

In order to dismiss the dramatic change <strong>in</strong> functions which might lead to great error <strong>in</strong><br />

the numerical <strong>in</strong>tegration <strong>in</strong> the simple sampl<strong>in</strong>g method, we <strong>in</strong>troduce <strong>in</strong> important<br />

sampl<strong>in</strong>g to select po<strong>in</strong>ts accord<strong>in</strong>g to the trend of function, and this selection is<br />

previously set.<br />

Referr<strong>in</strong>g to (1.1.1), we change it to<br />

In importance sampl<strong>in</strong>g, P is relevant with , and therefore we have<br />

We can also demonstrate that E[<br />

]=I <strong>in</strong> the sense of central limit theorem.<br />

An application of importance sampl<strong>in</strong>g will be illustrated here. In statistic physics, an<br />

observable is calculated <strong>in</strong> the follow<strong>in</strong>g means:<br />

This is easily done by simple sampl<strong>in</strong>g method, but the exponent function is such a<br />

dramatic function that we would use important sampl<strong>in</strong>g <strong>in</strong>stead to <strong>in</strong>crease efficiency.<br />

Therefore we write:<br />

(1.2.1)<br />

if P is carefully selected so that it is proportional to the distribution e , we will<br />

get <strong>in</strong>stantly that

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

Saved successfully!

Ooh no, something went wrong!