分子體積的計算

1、談?wù)劮肿芋w積的計算文Sobereva2011Sep20我曾在網(wǎng)上數(shù)次看到有人問怎么算分子體積，以及關(guān)于Gaussian的volume關(guān)鍵詞使用問題，從提問以及很多人的回答上看有不少人對這問題存在錯誤的認(rèn)識。本文就來談?wù)勥@個問題。首先要認(rèn)識到，分子體積不是一個可觀測量，在計算方法上也不可能有唯一的定義，因此“怎么計算分子體積”這個問題本身就是不嚴(yán)格的。分子體積就是分子表面內(nèi)部空間的體積，由于分子表面也沒有唯一的定義，所以不同的分子表面的定

2、義就會給出不同的分子體積定義。首先看一個簡單分子的例子作為范德華表面，這種表面通常能夠囊括分子98%以上的電子密度，這種范德華體積通常比范德華球疊加得到的范德華體積要大。而對于處于凝聚態(tài)的分子，考慮到分子間擠壓，建議用電子密度為0.002的等值面作為范德華表面（顯然，其體積小于0.001等值面對應(yīng)的體積）。計算范德華體積最常用的辦法是蒙特卡羅方法。首先，設(shè)立一個矩形盒子，將整個分子擴住，并且各個方向上都預(yù)留一定空間以避免將范德華表面截斷

3、，記這個盒子的體積為V_box。然后，在盒子里隨機分布m個點，依次檢驗這些點是否符合條件。對于范德華球疊加方式的定義，如果當(dāng)前點與任何一個原子核的距離小于相應(yīng)原子的范德華半徑，則認(rèn)為此點符合條件；而對于Bader的定義，若當(dāng)前點的電子密度大于閾值（0.001或0.002），就認(rèn)為符合條件。假設(shè)最后有n個點符合條件，那么分子的范德華體積就是nmV_box。如果m較小，那么算出來的體積是不精確的，想要增加精度，就必須增加m。這就像人口普查，