We extend the cleaving wall method to a nonpairwise additive potential. Using this method, we compute the anisotropy of crystal-melt interfacial free energy γ for Stillinger–Weber potential of silicon [F. H. Stillinger and T. A. Weber, Phys. Rev. B 31, 5262 (1985)]. The calculated γ for (100), (111), and (110) orientations are 0.42±0.02, 0.34±0.02, and 0.35±0.03J∕m2, respectively. The anisotropy in γ we found is consistent with the experimental observation that Si(100)-melt interface develops (111) facets and also helps in explaining a higher undercooling observed for Si(111)-melt interface in Czochralski method.

We present two methods suitable for controlling the translational degrees of freedom of a system when evaluating directly the free energy difference between the liquid and the solid phases by thermodynamic integration along a reversible path connecting these two phases. Such a constraint is crucial for an accurate prediction of the melting point by means of simulation. In one of the methods, the free energy difference was calculated by fixing one of the particles of the system at the center of the simulation box. In the second method, the free energy difference was calculated by constraining the center of mass of the system to a small region taken around the center of the simulation box. The correction to the free energy difference due to each constraint must be evaluated by a direct simulation. Both methods give consistent results when applied to a truncated and shifted Lennard-Jones system with cutoff radius of 2.5. However, the fixed particle constraint method is found to be more efficient computationally.

Stockmayer fluids are a prototype model system for dipolar fluids. We have computed the freezing temperatures of Stockmayer fluids at zero pressure using three different molecular-dynamics simulation methods, namely, the superheating-undercooling method, the constant-pressure and constant-temperature two-phase coexistence method, and the constant-pressure and constant-enthalpy two-phase coexistence method. The best estimate of the freezing temperature (in reduced unit) for the Stockmayer (SM) fluid with the dimensionless dipole moment μ*=1, √2, √3 is 0.656 ± 0.001, 0.726 ± 0.002, and 0.835 ± 0.005, respectively. The freezing temperature increases with the dipolar strength. Moreover, for the first time, the solid-liquid interfacial free energies γ of the fcc (111), (110), and (100) interfaces are computed using two independent methods, namely, the cleaving-wall method and the interfacial fluctuation method. Both methods predict that the interfacial free energy increases with the dipole moment. Although the interfacial fluctuation method suggests a weaker interfacial anisotropy, particularly for strongly dipolar SM fluids, both methods predicted the same trend of interfacial anisotropy, i.e., γ100 > γ110 > γ111.

Using a recently developed thermodynamic integration method, we compute the precise values of the excess Gibbs free energy (G e ) of the high density liquid (HDL) phase with respect to the crystalline phase at different temperatures (T ) in the supercooled region of the Stillinger-Weber (SW) silicon [F. H. Stillinger and T. A. Weber, Phys. Rev. B. 32, 5262 (1985)]. Based on the slope of G e with respect to T , we find that the absolute entropy of the HDL phase increases as its enthalpy changes from the equilibrium value at T ≥ 1065 K to the value corresponding to a non-equilibrium state at 1060 K.We find that the volume distribution in the equilibrium HDL phases become progressively broader as the temperature is reduced to 1060 K, exhibiting van-der-Waals (VDW) loop in the pressure-volume curves. Our results provides insight into the thermodynamic cause of the transition from the HDL phase to the low density phases in SW silicon, observed in earlier studies near 1060 K at zero pressure.Keywords amorphous silicon · liquid-liquid transition IntroductionThe liquid-amorphous transition in silicon, modeled by the Stillinger-Weber (SW) potential [23], has been intensely studied [9,19,20,2,22,7,24,15,18] with an aim of understanding the phase behavior of real silicon. In the initial molecular dynamics (MD) studies on SW silicon [20,2], it was found

The supercooled liquid silicon, modeled by Stillinger-Weber potential, shows anomalous increase in heat capacity C p , with a maximum C p value close to 1060 K at zero pressure. We study equilibration and relaxation of the supercooled SW Si, in the temperature range of 1060 K-1070 K at zero pressure. We find that as the relaxation of the metastable supercooled liquid phase initiates, a straight line region (SLR) is formed in cumulative potential energy distributions. The configurational temperature corresponding to the SLR is close to 1060 K, which was earlier identified as the freezing temperature of 4-coordinated amorphous network. The SLR is found to be tangential to the distribution of the metastable liquid phase and thus influences the broadness of the distribution. As the bath temperature is reduced from 1070 K to 1060 K, the effective temperature approaches the bath temperature which results in broadening of the metastable phase distribution.This, in turn, causes an increase in overall fluctuations of potential energy and hence an increase of heat capacity. We also find that during initial stages of relaxation, 4-coordinated atoms form 6membered rings with a chair-like structure and other structural units that indicate crystallization.Simultaneously a strong correlation is established between the number of chair-shaped 6-membered rings and the number of 4-coordinated atoms in the system. This shows that all properties related to 4-coordinated particles are highly correlated as the SLR is formed in potential energy distributions and this can be interpreted as a consequence of 'freezing' of amorphous network formed by 4-coordinated particles.