Constant-pressure Monte Carlo: Difference between revisions

From SklogWiki
Jump to navigation Jump to search
(New page: In '''constant-pressure Monte Carlo''' a trial change in the volume becomes one of the Monte Carlo moves. The weighting function, ''w'', is given by (Ref 1 Eq. 2) :<math>w \simeq ...)
 
m (correction of first eq.)
Line 2: Line 2:
The [[weighting function]], ''w'', is given by (Ref 1 Eq. 2)
The [[weighting function]], ''w'', is given by (Ref 1 Eq. 2)


:<math>w \simeq c \frac{V^N}{N!} \exp \left[ -\beta ( U(q^N)) +pV \right]</math>
:<math>w \simeq c \frac{V^N}{N!} \exp \left[ -\beta ( U(q^N) +pV )\right]</math>


for a change in volume, <math>\Delta V</math>, one has
for a change in volume, <math>\Delta V</math>, one has

Revision as of 14:39, 20 February 2009

In constant-pressure Monte Carlo a trial change in the volume becomes one of the Monte Carlo moves. The weighting function, w, is given by (Ref 1 Eq. 2)

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle w \simeq c \frac{V^N}{N!} \exp \left[ -\beta ( U(q^N) +pV )\right]}

for a change in volume, Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle \Delta V} , one has

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle r = \frac{w_{\rm new}}{w_{\rm old}} = \exp \left[ -\beta \left( \Delta U + p\Delta V - Nk_BT \ln \frac{V+\Delta V}{V}\right)\right]}

If Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle r \ge 1} then the move is accepted, and if Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle 0 < r < 1} then r is compared with a random number Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle 0 < x < 1} . If Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle x\leq r} then the move is accepted.

References

  1. R. Eppenga and D. Frenkel "Monte Carlo study of the isotropic and nematic phases of infinitely thin hard platelets", Molecular Physics 52 pp. 1303-1334 (1984)