Abstract
Abstract
The phase transition temperature of a two-dimensional (2D) lattice-gas model for the Ih-ice prism surface ((1010) surface on a wurtzite crystal structure) is calculated using the density-matrix renormalization-group (DMRG) method. Since the unit cell on the prism surface contains four atoms without inversion symmetry, the asymmetric version of the DMRG is required to obtain a precise value of the phase transition temperature. The effective bond energy on the surface is obtained by comparing the phase transition temperature with the faceting transition temperature observed for a real ice crystal. Using the effective bond energy, the temperature dependence of the coverage on sub-lattice sites is calculated. The surface tensions on the (0001) surface (the basal plane) and the prism plane are also calculated. The surface tension of the basal plane for the ice/water interface agrees well with the recent calculated values obtained by large scale molecular dynamics.
Subject
General Physics and Astronomy
Reference33 articles.
1. Growth kinetics of ice from the vapour phase and its growth forms. 1982;Kuroda;J. Cryst. Growth
2. Glancing-angle X-ray scattering studies of the premelting of ice surfaces. 1995;Dosch;Surf. Sci.
3. Ellipsometric study of the transition layer on the surface of an ice crystal. 1987;Furukawa;J. Cryst. Growth
4. X-ray diffraction pattern of quasi-liquid layer on ice crystal surface. 1987;Kouchi;J. de Physique