Smoothed particle hydrodynamics: a meshfree particle method, 英文以及中文译本，北大图书馆

这本书介绍了光滑粒子动力学方法（SPH）的起源、发展以及应用。我当时借看这本书，主要目的是看一看无网格方法有没有什么启示，另外，本书介绍了SPH与MD的一些耦合问题，还有一些并行的思想。我还下载了一些源程序，也许以后用得着。
浏览之后，这里仅记录下我感兴趣一些方面。

一、作为一种无网格粒子方法，和MD一样，需要使用最近粒子搜索法（NNPS，nearest neighbor search）。
1、全配对 all-pair（与MD中同，通过粒子之间距离判断是否发生作用，耗时最长，适合于小粒子数，O（N^2））
2、链表 linked-list（与MD中同，是用于光滑长度不变的情况，对应于MD为Rcut不变，O（N））
3、树形 tree O（NlgN）
4、Bucket Liu&Tu （2001， 2002）

二、与MD耦合
长度尺度耦合coupling length scale (CLS)
Q ：所谓多尺度，涉及空间尺度和时间尺度，分别都有哪些例子和应用，通常这两种尺度能分开吗？
各态遍历假设：时间平均=系综平均

三、分子间作用势
1、对于具有核子和电子相互作用的分子系统，需通过求解Schrödinger方程，可以通过Born-Oppenheimen近似大大降低复杂性。由于核子和电子的质量差异很大，需要分两步求解核子和电子的动力学性质：1）求解电子的Schrödinger方程得到对于固定核子的电子运动；2）求解核子的Schrödinger方程——ab initio MD
2、若选用经验参数简化ab initio MD，则可以减少计算时间——半经验方法
3、假设具有适当系数的解析势函数能够再生实验数据或者能够再生精确计算方法得到的精确解——MD中的经验势函数
Q：MD计算中一般也只是考虑两体势（pair potential），没有考虑多体势（multi-body potential）。是否说明也和DSMC一样，只考虑两体碰撞，因而只适用于稀薄气体？或者密度不能太大？
Morse势：适用于由化学键引起的相互吸引力的情况
Bunkingham势：由闭合壳层电子云的交叠而产生的强排斥力的情况
L-J：长距离的吸引项（-1/rij^6）表示了由于电子相关性而引起的范德华相互作用，以及深度为epsilon由电子云之间的非键交接而引起的强排斥项（1/rij^12）
建立原子势的方法 Xu&Liu （2003） MD+遗传算法，对实验数据进行逆向拟合