反應場算法
2022-10-05 07:53 作者:邱新龍分子模擬小破屋 | 我要投稿
在反應場(reaction field)算法中,作用在一個偶極矩i上的力被分為近程 相互作用和長程相互作用兩個部分。近程相互作用部分包括了反應球內(即截斷半徑rc內)所有偶極矩對偶極矩i的相互作用,長程相互作用部分是反應球外(即截斷半徑rc以外)的所有偶極矩對偶極矩i的相互作用。其中,反應球外的偶極矩被近似為介電常數(shù)ε的連續(xù)介質,并在反應球內形成一個反應場。由連 續(xù)介質電動力學可以知道,作用在偶極矩i上的反應場的大小正比于反應球的矩。

式中,μj為反應球內的偶極矩,求和遍及反應球內所有偶極矩,也包括偶極矩i。偶極矩i與反應場的作用能的貢獻為Ei? μj,被加到反應球內的相互作用上。
反應場算法的缺點是反應球內的偶極子數(shù)變化時,能量和受力的變化是不連 續(xù)的,但可以通過一個開關函數(shù)來解決不連續(xù)問題。

標簽: