Abstract
本文以正則系綜分子模擬法, 包括蒙地卡羅分子模擬法和分子動態法來配合Biosym分子模擬的商用軟體所提供的氬分子力場函數的參數和型式來模擬氬分子系統, 計算溫度在120K 時的體積-壓力等溫線, 發現不論使用蒙地卡羅分子模擬法或使用分子動態法來模擬氬分子系統, 都可以得到一個相近的壓力-體積-溫度關係, 並與 Petal和Teja 的狀態方程式所預測氬分子體積-壓力等溫線的結果做一個比較, 發現在高密度區和低密度區所得的結果差異不大, 但在相轉移區則有明顯的不同.我們以Maxwell 等面積定律來得到氬分子的相平衡數據, 發現分子模擬的結果不比三次方程式遜色. 我們另外以Cerius分子模擬的商用軟體所提供的氮分子力場函數的參數和型式來模擬氮分子系統, 並討論在分子模擬時的多原子分子系統的壓力計算方法, 發現以分子為觀點的壓力計算方法比以原子為觀點壓力計算方法較容易得到一個收斂的結果. 所以我們以蒙地卡羅分子模擬法,使用以分子為觀點的壓力計算法, 來計算氮分子系統在溫度100k 時的體積-壓力等溫線, 同樣與 Patel和Teja狀態方程式所預測氮分子體積-壓力等溫線的結果做一個比較, 發現在高密度區和低密度區所得的結果差異不大, 但在相轉移區則有明顯的不同. 我們仍以Maxwell等面積定律來得到氬分子的相平衡數據, 發現分子模擬的結果並不良好, 無法得到氮分子系統良好的相平衡數據, 此原因可能是由於相轉移區液氣界面的影響,因此我們建議使用吉伯氏系綜蒙地卡羅分子模擬法來忽略界面的影響, 以期得到一個正確的相平衡數據。