The Monte Carlo method of estimating statistical mechanical averages in the petite canonical ensemble, described by Rosenbluthet al.,Wood and Parker (for fluid systems), and Salsburget al.(for lattice models), has been extended to a general multicomponent lattice model in a restricted grand canonical ensemble. The procedure is applied to the two‐dimensional traingular lattice gas with periodic boundary conditions at a supercritical temperature (&bgr;&egr;= —ln2), and numerical results are presented for the energy, specific heat, density, isothermal compressibility, thermal‐expansion coefficient, and grand partition function (pressure) at &Dgr;=0.1, 0.0, —0.1, —0.2, —0.3, —0.4, —0.6, —0.8 (withB=100); atB=16, 25, 36, 49, 64, 100, 196 (with &Dgr;=0.0); where &Dgr;=&bgr;&mgr;—3&bgr;&egr;; where &egr; is the nearest‐neighbor interaction; &bgr;=(kT)—1, where &mgr; is the chemical potential; andBis the number of sites. These are compared (where possible) with previous (Salsburget al.) and additional petite‐ensemble Monte Carlo results. This comparison emphasizes the differentBdependence of intensive properties in these two ensembles and is supplemented by an asymptotic analysis of this difference. Properties in both ensembles display the type of irregularBdependence predicted by Lebowitz and Percus for very small systems. For the larger lattices, properties in the petite ensemble show a strongerBdependence than in the grand ensemble, which is in quantitative agreement with the leading term in the asymptotic analysis. A comparison with the exact analytical results (B= ∞, &Dgr;=0) indicates that the accuracy of the Monte Carlo procedure for the grand ensemble can be reliably estimated by a statistical analysis of partial averages over the Markov chain.
展开▼
