A simulation method for the estimation of anharmonic densities of stat
es of classical molecular models is described. The method is based on
the equilibrium energy distribution established in an uncoupled dimer
of the anharmonic molecule and a reference molecule whose density of s
tates is known analytically as a function of energy. Applications to o
ne-dimensional chain molecules and small clusters of atoms joined by M
orse bonds indicate that the method is simple and reliable.