《分子建模与模拟导论.ppt》由会员分享，可在线阅读，更多相关《分子建模与模拟导论.ppt（13页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、分子建模与模拟导论分子建模与模拟导论课课 程程 作业与考作业与考 核核课程主页：http:/ 操作系统环境及命令；3）C/C+编程语言；4）数值计算基本知识。课本课本：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 Ed.)Academic Press,2002.
2、参考书目参考书目：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 Text Book”(2nd Ed.)Pearson Education,2005.课
3、课 程程 大大 纲纲I 第一性计算第一性计算 1 电子结构方法 2 密度泛函理论 3 应用举例II 分子建模分子建模 1 通用分子建模方法 2 Lennard-Jones 势 3 金属体系力场 4 化学体系力场 5 粗粒化方法III 分子模拟方法分子模拟方法 1 系综与边界条件 2 蒙特卡罗方法 3 分子动力学方法 4 高级技术 5 数据分析方法IV 优化技术优化技术 1 寻找局部最小值 2 寻找全局最小值 3 寻找反应路径研研 究究 对对 象象实际材料分子体系及其抽象模型I.1.I.1.第第 一一 性性 计计 算算 电子结构方法电子结构方法解多体薛定谔方程：H(r)=E (r)方程中包含所有
4、原子核和电子的哈密顿量和波函数 原子间相互作用由主要由电子决定 计算系统最低能量构型，原子光谱等 计算化学反应势垒，速率等 Gaussian 软件包One-unit Kapton电子结构方法举例电子结构方法举例 最小能量态最小能量态I.2.I.2.第第 一一 性性 计计 算算 密密 度度 泛泛 函函 理理 论论解多体薛定谔方程：H(r)=E (r)用电子密度代替电子坐标，从而大大减少计算自由度 适于进行较大体系的计算 VASP 软件包密度泛函方法举例密度泛函方法举例 纳米金的晶体结构纳米金的晶体结构A.S.Barnard,X.M.Lin,L.A.Curtiss J.Phys.Chem.B 10
5、9,24465(2005).II.II.分分 子子 建建 模模 以原子为质点，原子间作用为经验势 通用方法：给定原子间作用力函数形式，从实验数据或第一性计算数据拟合参数 Lennard-Jones 势：惰性气体氩原子之间的相互作用，wikipedia 金属体系力场：多体作用 化学体系力场：成键作用，静电作用和范德华力作用 粗粒化：多个原子组合成一个粗粒化点Y.Wang,G.A.Voth to be submitted.Empirical force fields are determined by fitting experimental results or data from first
6、principles calculations Quality of empirical force fields has big influence on simulation results Capable of simulating up to millions of atoms(parallel computing)Solving Newtons Equations of Motion.III.1.III.1.分分 子子 动动 力力 学学 模模 拟拟III.2.III.2.蒙蒙 特特 卡卡 罗罗 模模 拟拟http:/cmm.cit.nih.gov/intro_simulation/node25.html 只计算势能，不用计算力 基于 Boltzmann 分布，只能模拟平衡态体系 模拟步长可以很大IV.1.IV.1.优优 化化 技技 术术 寻寻 找找 最最 小小 值值 总能运用标准方法找到局部最小值 寻找全局最小值具有很大的挑战性http:/www.cryst.bbk.ac.uk/PPS2/projects/day/TDayDiss/KineticStability.htmlIV.2.IV.2.优优 化化 技技 术术 寻寻 找找 反反 应应 路路 径径By Christoph Dellago体系从初态到终态的过程路径及发生的几率