高等学校教材：计算材料科学

　　第1章计算物理学和计算材料科学1
　　10绪论1
　　11计算物理学的起源与发展1
　　12计算物理学的特点和方法2
　　13计算物理学的进展5
　　14计算材料科学6
　　141分子动力学的发展历程7
　　142蒙特卡洛方法及其发展历程8
　　参考文献9第2章误差理论10
　　21误差的基本概念10
　　211准确值，误差限，真值，期望值和平均值10
　　212误差的具体类型11
　　213有效数与可靠数11
　　22随机误差12
　　221误差分布函数12
　　222随机误差的特点13
　　223随机误差的估算13
　　23误差传递公式15
　　231单变量误差的传递16
　　232多变量误差的传递16
　　233标准误差的传递公式16
　　234关于测量数列算术平均值的标准误差17
　　24计算机中的数值运算18
　　241整数18
　　242浮点数表示及误差18
　　243计算函数值的坏条件判别法18
　　244保证计算数值稳定的原则19
　　25矩阵运算中的误差问题20
　　251矩阵的条件数20
　　252线性方程组解的相对误差20
　　参考文献22第3章数值计算概要23
　　31范数与谱半径23
　　311范数的定义23
　　312矩阵的谱半径24
　　32线性方程组的解法24
　　321矩阵求逆法24
　　322选主元的三角分解法25
　　33矩阵的特征值问题27
　　331矩阵折叠法求本征值28
　　332实对称矩阵的雅克比对角化变换法30
　　333特征值的敏感性33
　　34矩阵的镜象变换法34
　　341镜象变换法34
　　342实对称三对角矩阵的本征值36
　　35曲线插值法37
　　351抛物型插值法的拉格朗日公式37
　　352抛物型插值法的牛顿公式38
　　36最佳一致逼近39
　　361切比雪夫多项式39
　　362最佳平方逼近39
　　363最小二乘法42
　　364正交多项式逼近42
　　37差分和差商运算47
　　371差分算子47
　　372差商公式48
　　38积分运算49
　　381多项式逼近的数值积分公式49
　　382高斯数值积分法51
　　39迭代运算54
　　391直接迭代法54
　　392牛顿迭代法54
　　393弦截法55
　　310微分方程的数值解法56
　　3101Numerov法56
　　3102龙格库塔法56
　　3103高阶微分方程和一阶微分方程组的龙格库塔法求解57
　　参考文献58第4章分子动力学运动方程解法59
　　41理论力学原理60
　　411最小作用量原理60
　　412拉格朗日（Lagrange）方程61
　　413哈密顿(Hamilton)方程63
　　414哈密顿（Hamilton）量与半经典量子化计算64
　　42粒子与粒子系综67
　　421相互作用的多粒子系综67
　　422粒子互作用模型的平衡统计性质——热力学性质的时间平均68
　　43粒子系统运动方程的数值解法69
　　431边界条件69
　　432Verlet算法70
　　433Gear算法71
　　44Gear预测矫正方法71
　　441表象表示法71
　　442表象变换72
　　443预测矫正法（predictorcorrector）73
　　45几种算法的比较76
　　451Verlet算法76
　　452跳蛙式算法77
　　453Beeman算法77
　　454比较几种算法的精度78
　　46复合函数算法79
　　461第一型算法80
　　462第二型算法80
　　47高阶综合矫正复合函数算法82
　　471第一型推广算法82
　　472第二型推广算法84
　　473第三型推广算法85
　　474第四型推广算法86
　　475第五型推广算法89
　　476更高阶预测公式92
　　参考文献93第5章计算机模拟的粒子系综94
　　51微正则系综（E，V，N）95
　　511计算热力学量95
　　512运动方程组的建立97
　　513数值求解的具体实施98
　　52正则系综（N，V，T）100
　　521正则系综的原理100
　　522正则系综分子动力学算法102
　　523正则系综与微正则系综的区别104
　　53等温等压系综(N，P，T)105
　　531等温扩展法105
　　532等压扩展法106
　　54体积可变的等焓等压(H，P，N)系综107
　　541晶胞结构矩阵与矢量的建立108
　　542仅有均匀恒压的情况108
　　543系统受到外部胁强时的一般情况109
　　55新恒压分子动力学模型110
　　551分子系统的运动方程110
　　552分子系统的守恒定律112
　　553对称的晶胞矢量h约束下的运动方程113
　　554受约束动力学晶胞的运动方程115
　　56自由能的计算和其他系综117
　　561巨正则系综118
　　562Gibbs系综119
　　563半巨系综119
　　57非平衡系综动力学120
　　571输运过程120
　　572非平衡分子动力学模拟（NEMD）120
　　573非平衡分子动力学在线性谐振子上的应用122
　　574自扩散算法123
　　参考文献124第6章势函数理论与模型125
　　61势函数发展简介125
　　611原子间相互作用势的发展概况126
　　612构成分子的原子间势127
　　613固体中的对势模型127
　　614多体势129
　　62半经验势129
　　621长程势129
　　622短程势130
　　63经验势131
　　64赝势模型132
　　641等效势概念133
　　642电子密度和赝密度134
　　643离子赝势134
　　65共价和离子晶体中的原子间互作用势135
　　651二体加三体势（团簇势）136
　　652Abell势模型137
　　653Tersoff势模型138
　　654KD势模型139
　　655分子间势模型140
　　656电子气势模型141
　　66紧束缚势模型和密度泛函理论143
　　661密度泛函理论143
　　662态密度和带能147
　　663紧束缚近似148
　　664能矩和二次矩148
　　665过渡金属和合金的紧束缚势150
　　67分子动力学和密度泛函理论的统一方法151
　　671第一性原理分子动力学方法152
　　672多原子体系动力学——CP运算方法154
　　68固体、分子经验电子理论及相关势函数162
　　681固体、分子经验电子理论162
　　682过渡金属和合金的结合能公式164
　　683碱金属和碱土金属的Morse势函数模型167
　　684固体中LennardJones型势函数的谢氏模型168
　　685谢氏势函数模型的改进型172
　　686基于“经验电子理论”的势函数模型的特点及存在的困难174
　　参考文献174第7章金属中的嵌入势模型177
　　71嵌入原子法177
　　711嵌入原子法基本原理178
　　712晶体结合能公式178
　　713由结合能导出的力学性质公式178
　　714排斥势函数和电子密度函数179
　　715有效电荷及嵌入能函数180
　　716作用在每个原子上的合力181
　　717表面能计算公式181
　　718采用Rose普适方程的嵌入势函数拟合法181
　　719嵌入原子法的理论推导简介182
　　7110金属合金中的原子间相互作用势183
　　72分析型EAM模型183
　　721Johnson的分析型EAM模型184
　　722EAM模型的改进185
　　723改进的普适分析型EAM模型185
　　724分析型EAM模型的应用与展望187
　　73FinnisSinclair经验势函数188
　　731FinnisSinclair（FS）模型188
　　732FinnisSinclair模型的改进型189
　　733FinnisSinclair模型的样条函数结构190
　　734弹性形变和间隙原子能190
　　74HCP结构金属的嵌入势模型191
　　741金属Zr的分析型EAM模型192
　　742HCP结构金属的分析型FS模型194
　　743金属Zr的分析型FS模型196
　　744晶格振动问题196
　　745HCP结构金属的非球对称模型197
　　746考虑到内弹性常数的HCP结构金属的分析型模型200
　　75合金材料中的势函数207
　　751Johnson对势函数的解析形式208
　　752Ackland的FS合金势模型210
　　753FS合金势模型用于金属间化合物的研究实例212
　　754金属间化合物的EAM势模型217
　　76Rose普适理论219
　　761能量和原子间距之间的普适关系219
　　762普适关系函数221
　　参考文献224第8章Monte Carlo方法226
　　81MC方法的基本思想226
　　811Buffon试验226
　　812随机事件的概率定义227
　　813马尔可夫（Markov）过程228
　　814大数定理和中心极限定理228
　　82伪随机数的产生和检验229
　　821产生随机数的方法230
　　822伪随机数的随机性检验231
　　83给定分布的随机数抽样方法232
　　831离散型随机变量的直接抽样233
　　832连续分布的随机变量抽样234
　　833Metropolis抽样法241
　　84微正则系综MC方法244
　　85正则系综MC方法246
　　851跃迁概率246
　　852正则系综的MC模拟法247
　　86等温等压系综MC方法248
　　87巨正则系综MC方法249
　　871巨正则系综MC算法之一249
　　872巨正则系综MC算法之二250
　　参考文献251第9章Monte Carlo方法的应用252
　　91MC方法在数值计算中的应用252
　　911计算定积分252
　　912解非线性方程组256
　　913解微分方程257
　　92Ising模型258
　　921磁性材料259
　　922二元合金261
　　923Heisenberg连续自由度Ising模型261
　　93自由能的计算263
　　931正则系综自由能的计算263
　　932二元系金属模型自由能的计算265
　　933计算任意固体的自由能及其在刚球fcc和hcp相中的应用266
　　934有相互作用力的Einstein晶体自由能270
　　94MC在三维组织和模拟辐射屏蔽中的应用271
　　941用MC方法模拟材料三维组织及其演变271
　　942用MC方法求解辐射屏蔽问题273
　　参考文献279第10章计算程序设计方法280
　　101晶胞结构设计280
　　1011坐标变换方法281
　　1012HCP晶胞结构程序设计281
　　102周期性边界条件的程序设计284
　　103分子动力学计算程序设计286
　　1031预测矫正法程序设计286
　　1032原子系综的分子动力学程序设计290
　　104METROPOLIS抽样算法292
　　参考文献292附录293
　　附录1可以演示的几个计算程序293
　　例1双原子分子在相空间中的量子轨道和能级293
　　例2监测各次扫描之间可观察量的关联297
　　例3计算氢分子内聚能的程序303
　　例4利用局域自旋密度近似法计算氢原子有效势的程序305
　　例5利用四面体法计算状态密度306
　　例6用EAM法计算铜的原子内聚能和空位形成能309
　　例7sp3杂化轨道间交叠积分的计算312
　　附录2物理学常数315
　　附表1物理学常数315
　　附表2用于低能分子动力学的通用单位316

　　计算材料科学是计算物理学的一个分支，涉及范围甚广。本书则重点介绍固体材料（主要是金属材料）中微缺陷的分子动力学模拟计算和蒙特卡洛模拟计算的实用性方法和系统方面的概念。分别从误差理论、数值计算概要、分子动力学运动方程解法、计算机模拟的粒子系统、势函数理论与模型、金属中的嵌入势模型、Monte Carlo方法及应用、计算程序设计等方面展开论述。并配有大量程序、公式、图、表及附录。
　　可供从事材料科学的研究生及该专业高年级大学生和其他学习者参考。