- Research Article
36
- 10.1063/1.453841
Monte Carlo simulation of the hard-sphere fluid with quantum correction and estimate of its free energy
- Mar 15, 1988
- The Journal of Chemical Physics
- Byoung-Jip Yoon + 1 more +1
A simple Monte Carlo (MC) method is introduced to treat the quantum-mechanical hard-sphere fluid. It uses the relative two-body Slater sum (instead of the classical Boltzmann factor) as a quantum-mechanical probability in the relative two-body configurational space. The MC calculations have been carried out for the hard-sphere fluid with quantum corrections at various densities and for various thermal wavelengths. The results are consistent with the interpretation that, in the quantum-mechanical picture, hard spheres repel each other at distances slightly greater (by an amount of the order of the thermal wavelength λT) than the contact collision distance of the classical treatment. The radial distribution functions of the classical and quantum-mechanical simulations differ considerably. The relative two-body free energy has been obtained from the relative two-body Slater sum, and the total free energy has been estimated (in preliminary calculations) by using the classical expression but replacing the classical hard-sphere diameter by the effective diameter. A good analytical fit to the MC results has been obtained for the equation of state of the quantum-mechanical hard-sphere fluid; the total free energy has been obtained (in extensive calculations) by integrating over the resulting pressures and densities. The pressure, internal energy, and free energy all increase with an increase in the thermal wavelength, i.e., with increasing quantum effect. From the equation of state, the second virial coefficients at various thermal wavelengths are in good agreement with theoretical values in the literature, and the third and fourth virial coefficients are predicted. Evidence has been found for a momentum condensation; the corresponding relative two-body entropy is negative. The heat capacity due to the relative two-particle kinetic energy has been found to be proportional to T−1/2 (approximately) whereas, classically, it is zero.
Read more