Editing Gibbs ensemble

Jump to navigation Jump to search
Warning: You are not logged in. Your IP address will be publicly visible if you make any edits. If you log in or create an account, your edits will be attributed to your username, along with other benefits.

The edit can be undone. Please check the comparison below to verify that this is what you want to do, and then publish the changes below to finish undoing the edit.

Latest revision Your text
Line 1: Line 1:
Here we have the ''N-particle distribution function''
Phase separation is one of the topics for which simulation techniques have preferentially been focused in the recent past. Different procedures have been used for this purpose. Thus, for the particular case of chain systems, we can perform simulations in the semi-grand canonical ensemble, histogram reweighting, or characterization of the spinodal curve from the study of computed collective scattering function.
(Ref. 1 Eq. 2.2)
The Gibbs ensemble Monte Carlo method mainly developed by Panagiotopoulos5 avoids the problem of finite size interfacial effects. In this method, a NVT ensemble containing two species is divided into two boxes. In addition to the usual particle moves in each one of the boxes, the algorithm includes moves steps to change the volume and composition of the boxes (at mechanical and chemical equilibrium). Transferring a chain molecule from a box to the other requires the use of an efficient method to insert chains. The configurational bias Monte Carlo method6,7 is specially recommended for this purpose.
 
:<math>\mathcal{G}_{(N)} ({\mathbf X}_{(N)},t)= \frac{\Gamma_{(N)}^{(0)}}{\mathcal{N}} \frac{{\rm d}\mathcal{N}}{{\rm d}\Gamma_{(N)}}</math>
 
where <math>\Gamma_{(N)}^{(0)}</math> is a normalized constant with the dimensions
of the [[phase space]] <math>\left. \Gamma_{(N)} \right.</math>.
 
:<math>{\mathbf X}_{(N)} = \{ {\mathbf r}_1 , ...,  {\mathbf r}_N ; {\mathbf p}_1 , ..., {\mathbf p}_N \}</math>
 
Normalization condition (Ref. 1 Eq. 2.3):
 
:<math>\frac{1}{\Gamma_{(N)}^{(0)}} \int_{\Gamma_{(N)}} \mathcal{G}_{(N)} {\rm d}\mathcal{N} =1</math>
 
it is convenient to set (Ref. 1 Eq. 2.4)
 
:<math>\Gamma_{(N)}^{(0)} = V^N \mathcal{P}^{3N}</math>
 
where <math>V</math> is the volume of the system and <math>\mathcal{P}</math> is the characteristic momentum
of the particles (Ref. 1 Eq. 3.26),
 
:<math>\mathcal{P} = \sqrt{2 \pi m \Theta}</math>
 
Macroscopic mean values are given by (Ref. 1 Eq. 2.5)
 
:<math>\langle \psi ({\mathbf r},t)\rangle= \frac{1}{\Gamma_{(N)}^{(0)}}
\int_{\Gamma_{(N)}}  \psi  ({\mathbf X}_{(N)}) \mathcal{G}_{(N)} ({\mathbf X}_{(N)},t) {\rm d}\Gamma_{(N)}
</math>
 
===[[Ergodic hypothesis |Ergodic theory]]===
Ref. 1 Eq. 2.6
 
:<math>\langle \psi \rangle = \overline \psi</math>
 
===[[Entropy]]===
Ref. 1 Eq. 2.70
 
:<math>S_{(N)}= - \frac{k_B}{ V^N \mathcal{P}^{3N}} \int_\Gamma  \Omega_1,... _N  \mathcal{G}_1,... _N {\rm d}\Gamma_{(N)}</math>
 
where <math>\Omega</math> is the ''N''-particle [[thermal potential]] (Ref. 1 Eq. 2.12)
 
:<math>\Omega_{(N)} ({\mathbf X}_{(N)},t)= \ln \mathcal{G}_{(N)} ({\mathbf X}_{(N)},t)</math>
 
==References==
# G. A. Martynov  "Fundamental Theory of Liquids. Method of Distribution Functions", Adam Hilger (out of print)
[[category: statistical mechanics]]
Please note that all contributions to SklogWiki are considered to be released under the Creative Commons Attribution Non-Commercial Share Alike (see SklogWiki:Copyrights for details). If you do not want your writing to be edited mercilessly and redistributed at will, then do not submit it here.
You are also promising us that you wrote this yourself, or copied it from a public domain or similar free resource. Do not submit copyrighted work without permission!

To edit this page, please answer the question that appears below (more info):

Cancel Editing help (opens in new window)