Hi,
The attached script was written to determine the melting point of silicon. I am using the NPT ensemble but am not seeing any change in volume at a temperature below or above the melting point or at any temperature. I am trying to find the Gibbs free energy for the solid and liquid phases at various temperatures to ultimately calculate the melting point. Could someone let me know what I am doing wrong?
Thanks,
Mihir.
in.SiliconMP (514 Bytes)