Monte Carlo simulation of crystal-fluid coexistence states in the hard-sphere system under gravity with stepwise control

Monte Carlo (MC) simulations were performed for hard spheres (with diameter sigma and mass m) placed between well-separated upper and lower hard walls. A periodic boundary condition was imposed in the horizontal direction. The system was exposed to the gravitational field with the acceleration due to gravity g. After preparing a melt as the initial state, g was increased stepwise up to mgsigma/k(B)T(identical with g(*))=1.5 or 2.0 with an increment Deltag(*) = 0.1; k(B)T is the temperature multiplied by Boltzmann's constant. We maintained g(*) at each value for 2.0 x 10(5) MC cycles. The transition of the system into a metastable state such as a polycrystalline state due to trapping phenomena was successfully avoided. A monotonic increase and subsequent saturation were observed for the development of the crystalline region formed at the bottom of the system. The development of this region accompanied a shrinkage of the defective (or less ordered) crystalline region that was formed between the bottom region and the fluid phase. As the development of the bottom region almost saturated, the defective region grew upward again.