Constant-force approach to discontinuous potentials

Abstract

Aiming to approach the thermodynamical properties of hard-core systems by standard molecular dynamics simulation, we propose setting a repulsive constant-force for overlapping particles. That is, the discontinuity of the pair potential is replaced by a linear function with a large negative slope. Hence, the core-core repulsion, usually modeled with a power function of distance, yields a large force as soon as the cores slightly overlap. This leads to a quasi-hardcore behavior. The idea is tested for a triangle potential of short range. The results obtained by replica exchange molecular dynamics for several repulsive forces are contrasted with the ones obtained for the discontinuous potential and by means of replica exchange Monte Carlo. We found remarkable agreements for the vapor-liquid coexistence densities as well as for the surface tension.

11 pages, 3 figures

References in corpus (6)

• Canonical sampling through velocity-rescaling
• Exchange Monte Carlo Method and Application to Spin Glass Simulations
• Simulated Tempering: A New Monte Carlo Scheme
• Calculation of surface tension via area sampling
• Direct calculation of interfacial tensions from computer simulation: Results for freely jointed tangent hard sphere chains
• Vapor-liquid surface tension of strong short-range Yukawa fluid

Cited by in corpus (1)

• Self-diffusion coefficient of the square-well fluid from molecular dynamics within the constant force approach