分子建模與模擬導(dǎo)論

1、課程主頁：http:/ 期中口試：20% 期末作業(yè)：70%預(yù)修課程預(yù)修課程：1）四大力學(xué)；2）linux/unix 操作系統(tǒng)環(huán)境及命令；3）c/c+ 編程語言；4）數(shù)值計算基本知識。課本課本：1）frank jensen “introduction to computational chemistry” (1st ed.) john wiley & sons ltd. 1999.2）daan frenkel and berend smit “understanding molecular simulation from algorithms to applications” (2nd

2、ed.) academic press, 2002.參考書目參考書目：1) mpallen and d. j. tildesley “computer simulation of liquids” oxford science publications, 1987.2) david chandler “introduction to modern statistical mechanics” oxford university press, 1987.3) syed mansoor sarwar, robert koretsky, syed aqeel sarwar “unix the tex

3、t book” (2nd ed.) pearson education, 2005.i 第一性計算第一性計算 1 電子結(jié)構(gòu)方法 2 密度泛函理論 3 應(yīng)用舉例ii 分子建模分子建模 1 通用分子建模方法 2 lennard-jones 勢 3 金屬體系力場 4 化學(xué)體系力場 5 粗粒化方法iii 分子模擬方法分子模擬方法 1 系綜與邊界條件 2 蒙特卡羅方法 3 分子動力學(xué)方法 4 高級技術(shù) 5 數(shù)據(jù)分析方法iv 優(yōu)化技術(shù)優(yōu)化技術(shù) 1 尋找局部最小值 2 尋找全局最小值 3 尋找反應(yīng)路徑實際材料分子體系及其抽象模型解多體薛定諤方程：h(r) = e (r)方程中包含所有原子核和電子的哈密頓量和

4、波函數(shù) 原子間相互作用由主要由電子決定 計算系統(tǒng)最低能量構(gòu)型，原子光譜等 計算化學(xué)反應(yīng)勢壘，速率等 gaussian 軟件包one-unit kapton解多體薛定諤方程：h(r) = e (r)用電子密度代替電子坐標(biāo)，從而大大減少計算自由度 適于進(jìn)行較大體系的計算 vasp 軟件包a. s. barnard, x. m. lin, l. a. curtiss j. phys. chem. b 109, 24465 (2005). 以原子為質(zhì)點，原子間作用為經(jīng)驗勢 通用方法：給定原子間作用力函數(shù)形式，從實驗數(shù)據(jù)或第一性計算數(shù)據(jù)擬合參數(shù) lennard-jones 勢：惰性氣體氬原子之間的相互作

5、用，wikipedia 金屬體系力場：多體作用 化學(xué)體系力場：成鍵作用，靜電作用和范德華力作用 粗?；憾鄠€原子組合成一個粗?；cy. wang, g. a. voth to be submitted. iijijjttffr iiittmfa 1iiittttvva 11iiittttrrv empirical force fields are determined by fitting experimental results or data from first principles calculations quality of empirical force fields has big influence on simulation results capable of simulating up to millions of atoms (parallel computing) ijf solving newtons equations of motion./intro_simulation/node25.html 只計算勢能，不用計算力 基于 boltzmann 分布，只能模擬平衡態(tài)體系 模擬步長可以很大 總能運用標(biāo)準(zhǔn)方法找到局部最小值 尋找全局最小值具有很大的挑戰(zhàn)性http:/www.cryst.bbk