水,作為地球上最重要的物質之一,其獨特的物理和化學特性對生命的存在至關重要。自從計算化學興起以來,科學家們便一直致力於用數學模型來模擬水的行為。這些模型不僅能夠預測水的物理特性,還有助於深入了解水在不同環境下的反應以及在生物體內的角色。
「水的特性與其分子結構密切相關,而這些結構又可以透過計算化學方法來模擬和預測。」
計算化學中的水模型主要用來模擬水分子聚集體、液態水及水溶液等。這些模型是根據量子力學、分子力學、實驗數據以及這些方法的組合來建立的。為了模仿水分子的特定特性,研究人員開發了多種類型的模型,這些模型通常可以從以下三個方面進行分類:(i) 交互作用點的數量稱為“位點”,(ii) 模型是剛性還是柔性,(iii) 模型是否包含極化效應。
水的模擬中,一個常見的做法是採用明確的溶劑模型,即基於具體分子的模型。相對於這些明確模型,還可以使用隱式溶劑模型,這類模型使用連續體模型來處理水的行為。這一領域的一個例子包括COSMO溶劑模型或極化連續體模型(PCM),甚至一些混合溶劑模型。
剛性模型被認為是最簡單的水模型,依賴於非鍵合的相互作用。在這些模型中,鍵合相互作用通過全局約束來隱式處理。靜電相互作用是基於庫侖定律而建模,排斥和色散力則使用Lennard-Jones潛力來描述。這些潛力模型,如TIP3P(可轉移的三點分子潛力)和TIP4P等,表示為:
E = ∑(kC * qi * qj / rij) + (A / rOO^12) - (B / rOO^6)
其中,kC是靜電常數,qi和qj是相對於電子電荷的部分電荷,rij是兩個原子之間的距離。許多水模型中,Lennard-Jones項僅適用於氧原子之間的相互作用。各種水模型的幾何參數(如OH距離和HOH角度)根據模型而有所不同。
二位模型基於熟悉的三位SPC模型,能夠使用位點重整化的分子流體理論來預測水的介電特性。而三位模型則對應於水分子的三個原子,每個位點都有電荷。三位模型計算效率高,廣泛用於許多分子動力學模擬中。
「目前常見的三位模型,如TIP3P,對於計算特定熱量的性能表現良好。」
例如,SPC/E模型為勢能函數添加了極化修正,使得其所得到的水的密度和擴散常數更優於SPC模型。TIP3P模型則在CHARMM力場中得到廣泛應用,並在原始模型上進行稍微的修改,使其更適合生物分子的模擬。
靈活的SPC水模型是一種重新參數化的三位水模型,與剛性的SPC模型不同,靈活模型在分子動力學模擬中能正確描述水的密度和介電常數。這種模型在多個計算程序中得以實現,例如MDynaMix和Abalone。
四位模型通過在三位模型的氧原子附近添加一個虛擬原子來改善水分子的電荷分佈,最早的此類模型可追溯至1933年的Bernal–Fowler模型。雖然該模型在歷史上有重要意義,但並未能夠很好地再現水的主要性質。
TIP4P模型被廣泛應用於計算化學軟件中,並在針對生物分子系統的模擬中發揮關鍵作用,而如OPC模型等新型水模型則能夠更準確地描述水的電特性。
五位模型,儘管計算成本較高,近年來在TIP5P模型的引入下逐漸取得進展。五位模型更好地復現了水二聚體的幾何結構,並且能夠精確捕捉實驗數據。六位模型則在數據中融合了過去的所有模型特點,專為研究水和冰系統而設計。
「計算化學中,水的模擬不僅是一個技術挑戰,更是理解生命運作的一把鑰匙。」
水模型的計算成本隨著位點數量的增加而上升。對於分子動力學模擬,隨著位點數增加,所需計算的原子間距數量也跟著增長。然而,這些模型的發展不僅僅是數學敘述,更是自然界中水的真實表現的縮影。隨著技術的進步,我們能否在不久的將來找到顯示水更多奧秘的模型呢?