Configurational temperature: Verification of Monte Carlo simulations

被引:127
作者
Butler, BD
Ayton, O
Jepps, OG
Evans, DJ
机构
[1] Natl Inst Stand & Technol, Phys & Chem Properties Div, Boulder, CO 80303 USA
[2] Australian Natl Univ, Res Sch Chem, Canberra, ACT 0200, Australia
关键词
D O I
10.1063/1.477301
中图分类号
O64 [物理化学(理论化学)、化学物理学];
学科分类号
070304 ; 081704 ;
摘要
A new diagnostic that is useful for checking the algorithmic correctness of Monte Carlo computer programs is presented. The check is made by comparing the Boltzmann temperature, which is input to the program and used to accept or reject moves, with a configurational temperature k(B)T(config) = \del(q)Phi\(2)/del(q)(2)Phi. Here, Phi is the potential energy of the system and del(q) represents the dimensionless gradient operator with respect to the particle positions q. We show, using a simulation of Lennard-Jones particles, that the configurational temperature rapidly and accurately tracks changes made to the input temperature even when the system is not in global thermodynamic equilibrium. Coding and/or algorithm errors can be detected by checking that the input temperature and T-config agree. The effects of system size and continuity of Phi and its first derivative on T-config are also discussed. (C) 1998 American Institute of Physics. [S0021-9606(98)52640-2].
引用
收藏
页码:6519 / 6522
页数:4
相关论文
共 3 条