动力学蒙特卡洛方法(KMC)及相关讨论

动力学蒙特卡洛方法(KMC)及相关讨论星期二,2010-05-1101:05—\o"浏览用户信息"satchel1979动态模拟在目前的计算科学中占据着非常重要的位置。随着计算能力和第一原理算法的发展，复杂的动态参数（扩散势垒、缺陷相互作用能等）均可利用第一原理计算得出。因此，部分复杂的体系动态变化，如表面形貌演化或辐射损伤中缺陷集团的聚合-分解演变等，已可以较为精确的予以研究。KMC——动力学蒙特卡洛方法（kineticMonteCarlo）原理简单，适应性强，因此在很多情况下都是研究人员的首选。此外，KMC在复杂体系或复杂过程中的算法发展也非常活跃。本文试图介绍KMC方法的基础理论和若干进展。KMC方法基本原理在原子模拟领域内，分子动力学(moleculardynamics,MD)具有突出的优势。它可以非常精确的描述体系演化的轨迹。一般情况下MD的时间步长在飞秒(s)量级，因此足以追踪原子振动的具体变化。但是这一优势同时限制了MD在大时间尺度模拟上的应用。现有的计算条件足以支持MD到10ns，运用特殊的算法可以达到10s的尺度。即便如此，很多动态过程，如表面生长或材料老化等，时间跨度均在s以上，大大超出了MD的应用范围。有什么方法可以克服这种局限呢？当体系处于稳定状态时，我们可以将其描述为处于维势能函数面的一个局域极小值（阱底）处。有限温度下，虽然体系内的原子不停的进行热运动，但是绝大部分时间内原子都是在势能阱底附近振动。偶然情况下体系会越过不同势阱间的势垒从而完成一次“演化”，这类小概率事件才是决定体系演化的重点。因此，如果我们将关注点从“原子”升格到“体系”，同时将“原子运动轨迹”粗化为“体系组态跃迁”，那么模拟的时间跨度就将从原子振动的尺度提高到组态跃迁的尺度。这是因为这种处理方法摈弃了与体系穿越势垒无关的微小振动，而只着眼于体系的组态变化。因此，虽然不能描绘原子的运动轨迹，但是作为体系演化，其“组态轨迹”仍然是正确的。此外，因为组态变化的时间间隔很长，体系完成的连续两次演化是独立的，无记忆的，所以这个过程是一种典型的马尔可夫过程(Markovprocess)，即体系从组态到组态，这一过程只与其跃迁速率有关。如果精确地知道，我们便可以构造一个随机过程，使得体系按照正确的轨迹演化。这里``正确''的意思是某条给定演化轨迹出现的几率与MD模拟结果完全一致(假设我们进行了大量的MD模拟，每次模拟中每个原子的初始动量随机给定)。这种通过构造随机过程研究体系演化的方法即为动力学蒙特卡洛方法(kineticMonteCarlo,KMC)[1]。指数分布与KMC的时间步长在KMC模拟中，构造呈指数分布的随机数是一个相当重要的步骤。这一节中我们对此进行讨论。因为体系在势能面上无记忆的随机行走，所以任意单位时间内，它找到跃迁途径的概率不变，设为。因此在区间内，体系不发生跃迁的概率为类似的，在区间内，体系不发生跃迁的概率为以此类推，当时，在区间内，体系不发生跃迁的概率为因此，当趋于时，体系不发生跃迁的概率为

(1)这一行为类似于原子核的衰变方程。从方程(1)我们可以得到单位时间内体系跃迁概率。从方程(1)的推导过程可以看出体系的跃迁概率是一个随时间积累的物理量，因此对时间积分到某一时刻必然等于，也即。因此我们立即可以得到[1]

(2)是体系处于态时所有可能的跃迁途径的速率之和，即

(3)对于每个具体的跃迁途径，上述讨论均成立。因此，我们可以定义单位时间内体系进行跃迁的概率为

(4)

单位时间内体系的跃迁概率呈指数分布这一事实说明KMC的时间步长也应是指数分布。因此我们需要产生一个指数分布的随机数序列。这一点可以非常容易的通过一个(0,1]平均分布的随机数序列转化得到：从而

(5)最后一步是因为和的分布相同。也可以通过上述步骤从方程(4)得到。计算跃迁速率过渡态理论(TST)决定了KMC模拟的精度甚至准确性。为避开通过原子轨迹来确定的做法(这样又回到了MD的情况)，一般情况下采用过渡态理论(transitionstatetheory,TST)进行计算[2]。在TST中，体系的跃迁速率决定于体系在鞍点处的行为，而平衡态(势阱)处的状态对其影响可以忽略不计。如果大量的相同的体系组成正则系综，则在平衡状态下体系在单位时间内越过某个垂直于跃迁途径的纵截面的流量即为。简单起见，假设有大量相同的一维双组态(势阱)体系，平衡状态下鞍点所在的假想面(对应于流量最小的纵截面)为，则TST给出该体系从组态A迁出到B的速率为[5,6]

(6)方程(6)中表示在组态A所属态空间里对正则系综的平均。表示只考虑体系从组态A迁出而不考虑迁入A的情况(后一种情况体系也对通过纵截面的流量有贡献)。根据普遍公式设体系的哈密顿量为，即可分解为动能和势能，同时设粒子坐标时体系处于组态A。则方程(6)可写为

(7)上式中无限小量是为了将函数全部包含进去。最后一项对于函数的系综平均可以直接通过MetropolisMonteCarlo方法计算出来：计算粒子落在范围内的次数相对于Metropolis行走总次数的比例。方程(7)最后等于

(8)将上述讨论扩展到3维情况非常直接，这里只给出结果，详细讨论请参阅文献[5]:

(9)其中是纵截面方程，代表3维情况中粒子流动方向与截面法向不平行对于计数的影响。简谐近似下的过渡态理论(hTST)虽然上一节已经给出了TST计算跃迁速率的方法，但是在具体工作中，更多地是利用简谐近似下的过渡态理论(harmonicTST,hTST)通过解析表达式给出。根据TST，跃迁速率为[3]

(10)其中为在跃迁中体系在鞍点和态处的自由能之差将上式代入方程(10)，可以得到

(11)hTST认为体系在稳态附近的振动可以用谐振子表示，因此其配分函数是经典谐振子体系的配分函数。分别写出体系在态和鞍点处的配分函数和：

根据Boltzmann公式，

(12)并将配分函数代入，则方程(11)得

(13)方程(13)在通常的文献上经常可以见到。声子谱可以通过Hessian矩阵对角化或者密度泛函微扰法(DFPT)求出，而就是的势垒，可以通过NEB或者drag方法求出。因此，方程(13)保证了可以通过原子模拟(MD或者DFT方法)解析地求出。事实上这个方程有两点需要注意。首先虽然方程(10)中出现了普朗克常数，但是在最终结果中被抵消了。这是因为TST本质上是一个经典理论，所以充分考虑了统计效应后不会出现[1]。其次，方程(13)表明对于每一个跃迁过程，鞍点处的声子谱应该单独计算。这样会大大增加计算量，因此在绝大部分计算中均设前置因子为常数，不随跃迁过程而变化。具体数值取决于体系，对于金属而言，一般取Hz。KMC几种不同的实现算法

点阵映射

到目前为止，进行KMC模拟的所有理论基础均已具备。但是前面所进行的讨论并没有联系到具体的模型。KMC在固体物理中的应用往往利用点阵映射将原子与格点联系起来。从而将跃迁(事件)具象化为原子格点关系的变化。比如空位(团)/吸附原子(岛)迁移等等。虽然与实际情况并不完全一致，但这样做在很多情况下可以简化建模的工作量，而且是非常合理的近似。很多情况下体系中的原子虽然对理想格点均有一定的偏离，但是并不太大()，因此这种原子点阵映射是有效的。这种做法的另一个好处是可以对跃迁进行局域化处理。每条跃迁途径只与其近邻的体系环境有关，这样可以极大的减少跃迁途径的数目，从而简化计算[1]。需要指出的是，这种映射对于KMC模拟并不是必须的。比如化学分子反应炉或者生物分子的生长等等，这些情况下根本不存在点阵。无拒绝方式

KMC的实现方法有很多种，这些算法大致可以分为拒绝(rejection)和无拒绝(rejection-free)两种范畴。每种范畴之下还有不同的实现方式。本文只选择几种最为常用的方法加以介绍。I.直接法直接法(directmethod)是最常用的一种KMC算法，其效率非常高。每一步只需要产生两个在之间平均分布的随机数和。其中被用来选定跃迁途径，确定模拟的前进时间。设体系处于态，将每条跃迁途径想象成长度与跃迁速率成正比的线段。将这些线段首尾相连。如果落在线段中，这个线段所代表的跃迁途径就被选中，体系移动到态，同时体系时间根据方程(5)前进。总结其算法如下：根据方程(4)计算体系处于态时的总跃迁速率；选择随机数；寻找途径，满足；体系移动到态，同时模拟时间前进；重复上述过程。需要指出的是，虽然一般步骤4中的根据方程(5)生成，但是如果将其换为并不会影响模拟结果。在文献[5]和[6]中均采用这种方式。II.第一反应法第一反应法(firstreactionmethod,FRM)在思路上比直接法更为自然。前面说过，对于处于稳态的体系而言，它可以有不同的跃迁途径可以选择。每条途径均可以根据方程(4)给出一个指数分布的"发生时间"，也即从当前算起第一次发生的时间。然后从中选出最小值(最先发生的"第一反应")，体系跃迁到相应的组态，模拟时间相应地前进。总结其算法如下：设共有条反应途径，生成个随机数；根据公式，给出每条路径的预计发生时间；找出的最小值；体系移动到态，同时模拟时间前进；重复上述过程。可以看出，这种算法的效率比直接法低下，因为每一步KMC模拟需要生成个随机数。通常情况下KMC模拟需要步来达到较好的统计性质，如果每一步都需要生成个随机数，则利用这种方法需要一个高质量的伪随机数发生器，这一点在比较大时尤为重要。III.次级反应法次级反应法(nextreactionmethod,NRM)是FRM方法的一种衍生方法，其核心思想是假设体系的一次跃迁并不会导致处于新态的体系对于其他跃迁途径的舍弃(比如充满可以发生种化学反应的分子，第一种反应发生并不会造成别的反应物的变化)，这样体系还可以选择中的次小值，从而跃迁到态，模拟时间前进。如果这次跃迁还可以满足上述假设条件，再重复上述过程。理想情况下，平均每一步KMC模拟只需要生成1个随机数。这无疑会大大提高效率以及时间跨度。但是实际上NRM的假设条件很难在体系每次跃迁之后都得到满足，在固体物理的模拟中尤其如此，因此其应用范围集中于研究复杂化学环境下的反应过程。试探-接受/拒绝方式

这一大类算法虽然在效率上不如直接法，但是它们所采用的试探-接受/拒绝在形式上更接近MetropolisMC方法，而且可以很方便的引入恒定步长，即固定。因此有必要进行详细的介绍。IV.选择直接法选择直接法在决定体系是否跃迁方面和MetropolisMC方法形式上非常相像，均是通过产生随机数和预定的阈值比较决定事件是否被采纳。具体算法如下：设共有条反应途径，选择反应速率最大值，设为。生成在均匀分布的随机数；设；如果<，则体系跃迁至新态，否则保持在态；模拟时间前进；重复上述过程。这种方法的长处在于每一步只需要生成一个随机数。但是缺点也很明显，对于反应速率相差太大，尤其是只有一个低势垒途径(与其他途径相比过大)的体系来讲，这种方法的效率会非常低下。某些情况下，这种低效率问题可以通过如下方法改进：将全部途径按照的大小分为几个亚组，每个亚组选定一个上限。但是这一步骤在整个KMC模拟过程中可能需要重复很多次，因此并不能完全解决问题。事实上低势垒在KMC中是个普遍的问题。这一点在后面还要简要提及。V.恒定步长法与上述四种方法不同，恒定步长法(constanttimestepmethod,CTSM)中体系的前进时间是个给定的参数\cite{dawnkaski}。在理想情况下，CTSM与直接法效率相同，每一步只需产生两个随机数。具体算法如下：给定恒定时间步长；将所有途径(共有个)设为长度恒为的线段，生成在均匀分布的随机数，选择途径；生成在均匀分布的随机数，如果<，则体系跃迁至新态，否则保持在态；模拟时间前进；重复上述过程。实际模拟中，需要满足(1)小于(见"第一反应法")，以及(2)对于最大的途径，接受率大致在0.5。其中第一个条件保证了所有的迁移途径发生概率都小于1，第二个条件则保证体系演化的效率不会过于低下。CTSM是非常行之有效的一类KMC算法，但是选择时需要特别的注意以保证效率。决定于具体体系以及模拟温度。这在一定程度上增加了CTSM的实现及使用难度。低势垒问题前面已经指出，低势垒的途径需要特别注意。如果体系在演化过程中一直存在着势垒较其他途径低很多的一个或几个途径，会对模拟过程产生不利的影响。这个问题被称之为低势垒问题。低势垒途径对于KMC模拟最直接的影响就是大大缩短了模拟过程所涵盖的时间跨度。这一点可以从方程(5)中看出。更为深刻的影响在于，这些由低势垒的途径联系起来的组态会组成一个近似于封闭的族。体系会频繁的访问这些态，而其他的对于体系演化更为重要的高势垒途径被选择的概率非常低，这显然会降低KMC的模拟效率。例如，吸附原子在高指数金属表面扩散，其沿台阶的迁移所对应的势垒要远低于与台阶分离的移动。这样，KMC模拟的绝大部分时间内吸附原子都在台阶处来回往复，而不会选择离开台阶在平台上扩散。这显然不是我们希望看到的情形。一种解决办法是人为地将这些低势垒加高以降低体系访问这些组态的几率，但是无法预测这种干扰是否会造成体系对于真实情况的严重偏离。另一种选择是利用NRM或者CTSM进行模拟，但是其效果如何尚待检测。如果考察体系的势能面，这类低势垒的途径一般处在一个"超势阱"之中。体系在这个超势阱中可以很快的达到热平衡，所需时间要短于从其中逸出的时间。如果可以明确的知道超势阱所包含的组态以及从超势阱逸出的所有途径，我们就可以按照Boltzmann分布合理的选择其中一条途径，使得体系向前演化。但是如何确定哪些组态包含在超势阱之中以及体系是否已在其中达到热平衡本身就是两个难题。对于第一点，Mason提出可以利用Zobrist密钥法标定访问过于频繁的组态[7]；Novotny则提出通过建立及对角化一个描述体系在这些组态间演化的传递矩阵来解决第二点[8]。对这个问题的详细讨论已超出了本文的讨论范围，请参阅文献[7]以及[8]。实体动力学蒙特卡洛方法OKMC上述的KMC都假设任何时候原子均处于其理想点阵格子上。但是很多情况下这种点阵映射是无效的，比如间隙原子或者位错。这类结构缺陷的运动在材料的辐射损伤和老化过程中扮演着非常重要的角色。而且与单个原子或者空位的运动相比，这类缺陷的运动时间跨度更长，也更为复杂，比如间隙原子团和空穴的湮没，间隙原子团的解构/融合，或者位错的攀移/交滑移等等。传统的KMC算法很难有效的处理这类问题，一方面是因为时间跨度太大，另一方面这类缺陷各自均可视为独立的实体(object)，其运动更近似于系统激发，因此单个或几个原子运动的积累效果很多情况下并不能有效地反应这些实体的整体运动。实体动力学蒙特卡洛方法(ObjectkineticMonteCarlo,OKMC)就是为了处理这类问题而被提出的。OKMC在算法上与普通的KMC完全一样。需要注意的地方是在OKMC中并不存在原子点阵。所有的实体在一个真空的箱子中按照其物理实质离散化运动，比如位错环的最小移动距离是其Burgers矢量大小，方向则为Burgers矢量方向；空位的移动距离为第一近邻或第二近邻的原子间距，等等。模拟过程中我们需要追踪该实体的形心，从而决定其位置、移动距离等等。此外，OKMC中对于跃迁速率的确定也和普通的KMC有所区别。本文前面已经指出，可以表达为的形式。普通的KMC假定为常数，不同途径的由决定。但是在OKMC的模拟中，的直接确定非常困难，因此一般的策略是对于特定的事件(包括实体自身的运动以及不同实体间的反应等)，跃迁势垒保持恒定，而将前置因子视为实体规模(所包含的原子/空位数目)的函数，通过MD模拟得出，一般而言可以表示为形如的表达式，其中和是拟合参量，是实体规模。最后需要注意的是在OKMC的模型中，实体有空间范围，因此需要一个额外的参数来表征其空间半径(假设为球形分布，否则的数目多于一个)。在模拟不同实体间的反应时，需要特别考虑其形心的间距，如果小于"反应距离"，即，反应一定进行，否则认为两个实体互相独立。Domain利用OKMC研究了Fe-Cu合金的辐射损伤[9]，在模拟中考虑了间隙原子(空位)的聚合、间隙原子(空位)团的发射、间隙原子-空位湮没、空位团对杂质的捕获、表面对于空位(团)的捕获、甚至辐射轰击引起的间隙原子(空位)萌生、增殖等等事件。从中可以看出，对于OKMC，一个棘手的问题是需要预先想到所有的事件。此外，OKMC所需要的所有参量基本上不可能通过原子模拟直接获得，人为的设定参数不可避免。这些参数会在多大程度上决定OKMC的准确程度无法预先得知。需要根据现有的实验数据进行修改、调试。这些困难都限制了OKMC的普及。但是如前所述，这种方法可以有效地进行大尺度的时间(天)和空间模拟(m以上)，而且对于缺陷的描述更为直接和符合直观，因此在材料研究中同样占有重要的地位。KMC的若干进展等时蛙跳算法(-leapKMC)

引入这类算法前，我们先简要介绍两个常用的离散分布：泊松分布(PoissonDistribution,PD)以及二项式分布(BinomialDistribution,BD)。泊松随机数定义为给定事件发生率以及观测时间下事件发生的数目。如果用代表给定的发生数目，则恰好等于的概率是一个泊松分布：

(14)也即如果产生一个泊松随机数序列，则这个序列符合泊松分布PD。需要指出，是无界的，范围是任意非负整数。与其类似，二项式随机数定义为重复次独立的成功率均为的伯努利实验的成功数。如果给定成功数，则恰好等于的概率是一个二项式分布：

(15)为了和本文中的标号一致，我们将跃迁的成功率表示为，将方程(15)重新写为

(16)与PD不同，BD中的是有界的，为0到之间的任意整数。可以看出，如果将这两种随机数理解为给定跃迁路径(发生率为)在一定的时间步长()内发生的次数，则可以立即运用于粒子数空间内的KMC中，其时间范围可以得到很大提高。这就是等时蛙跳算法-leapKMC[10,11]。-leapKMC方法最早由Gillespie提出，通过PD[方程(14)]，在给定时间步长下决定每个跃迁途径发生的次数，然后将体系移到这些跃迁累计发生后产生的新态。因为每一步模拟体系不止发生一次跃迁，所以模拟的速度可以大大加快。我们以多种反应物在化学反应炉中的演化为例加以详细说明。设在炉内共有种分子，在时刻各自的个数为，则在粒子数空间中构成一个矢量，或称为一个组态。总共有种反应路径。对于给定的，反应速率是占据态的函数。此外，我们单独定义一个矢量

，其中由通过反应而得，即。因此的元素代表反应所引起的种分

子的数目变化。由此建立算法如下：VI.PD--leapKMC[10]给定恒定时间步长；对于每条反应途径按照方程(14)生成泊松随机数序列，按照模拟步数从序列中找出每种反应发生的次数；按照更新体系；模拟时间前进；重复上述过程。Gillespie仔细考虑了的选择条件，称为蛙跳条件(leapcondition)：

(17)其中

如前所述，没有上限，因此即使满足方程(17)，在模拟过程中也可能会出现某种分子总数为负数的情况，这显然不符合实际，也是PD--leapKMC的一个弱点。Tian和Burrage提出可以用二项式分布BD取代PD，因为有上限，所以可以有效的解决这个问题。此外，他们对于某种分子参与多种反应的情况也进行了考虑，从而提高了-leapKMC的稳定性和普适性。其算法如下：VII.BD--leapKMC[11]

给定恒定时间步长，满足；对于每条反应途径按照方程(16)生成二项式随机数序列，按照模拟步数从序列中找出每种反应发生的次数；如果有某种分子同时参与了和，则首先生成然后通过确定的发生次数；按照更新体系；模拟时间前进；重复上述过程。步骤1、2中出现的是参与反应的各类分子的个数的最小值，即此外Gillespie，Tian和Burrage还考虑用预测时刻体系状态的方法来进一步提高精度。具体请参阅文献[10,11]。如果-leap算法和OKMC结合起来可以进一步加大模拟的时间尺度，但是目前还没有这方面工作的介绍。基于即时动态分析的KMC方法(on-the-flyKMC)到目前为止，所有的KMC都是在模拟之前建立好所有可能的跃迁途径。但是实际上"所有"是很难达到的目标。因为很多途径远离一般的直觉，而且在演化过程中体系有可能寻找到新的途径。因此，跃迁途径应该随着体系的演化而不断更新，是动态的过程。Henkelman和Jönsson将途径搜索和KMC结合起来，提出了即时动态的KMC方法on-the-flyKMC[12]：在每一个稳态(势阱)处，选定一个激活原子(一般是近邻不饱和的原子)，在以其为中心的局部区域内引入呈高斯分布的随机位移，即加入扰动，然后利用dimer方法[13]寻找所有可能的跃迁途径。建立起即时的途径库之后再通过普通KMC算法进行模拟。显然，这种方法的计算量非常大，需要一个有效的标识方法来识别所有已经遇到过的途径以避免重复计算。Trushin提出可以利用包括至激活原子第三壳层的所有格点(顺时针排列)的占据与否(分别标记为1和0)来构建二进制数，从而根据始态和终态的标号来唯一地标识某条途径[14]，例如，激活原子标为"1"，其第一壳层的原子标记为"2","3",,""，依此类推，然后将原子的标号""作为二进制的数位，这样，每一个稳态都有唯一的一个二进制数与之对应。虽然仍不完善，但是这种方法具有非常清晰的逻辑结构，具有良好的扩展性。和KMC方法一般情况下KMC的大部分时间花费在选择途径上。如果采用普通的方法，即循环叠加直至从而选择，这种情况下计算用时与途径数目呈线性增长，即算法。按照二叉树安排不同数目的之和可以改进到[15]：

将所有作为树叶(不足2整数次幂的叶子由0填补)，每两片叶子之和作为父节点，依次类推直至树根。一株二叉树构建完毕后，生成一个随机数，由树根开始寻找，若不大于左子节点，沿左分支向下寻找；否则设，沿右分支向下寻找，直至树叶，体系按途径演化。Slepoy和Thompson等进一步提出分流-拒绝(composition-rejection,CR)方法以实现搜索用时与途径总数无关的算法[16]：(1)先找出和，按照()将条途径分为个组，，(2)然后生成随机数，按照上述二叉树寻找所落入的组别，(3)再生成两个随机数和，设，其中为该组中包含的途径数，，如果，则选择途径，否则重复步骤(3)，直至有一条途径被选中为止。可以看出，CR算法虽然搜索速度很快，但是每一步KMC需要产生至少4个随机数(用于确定前进时间)，因此需要高质量的随机数发生器。不过对于跃迁途径复杂的体系演化而言，CR的效率无疑是很有吸引力的。

[1]A.F.Voter,{\itRadiationEffectsinSolids}(Springer2006)p.1-24.

[2]H.Eyring,J.Chem.Phys.3,107(1935).

[3]P.Kratzer,MultiscaleSimulationMethodinMolecularScience(NICSerices,Vol.42,Forschungszentrum,Jülich2009)p.51-76.

[4]E.J.Dawnkaski,D.SrivastavaandB.J.Gamson,J.Chem.Phys.102,9401(1995).

[5]A.F.VoterandJ.D.Doll,J.Chem.Phys.80,5832(1984).

[6]A.F.Voter,Phys.Rev.B34,6819(1986).

[7]D.R.Mason,T.S.HudsonandA.P.Sutton,Comp.Phys.Comm.165,37(2005).

[8]M.A.Novotny,Phys.Rev.Lett.74,1(1994);Erratum75,1424(1995).

[9]C.Domain,C.S.BecquartandL.Malerba,J.Nucl.Mater.335,121(2004).

[10]D.T.Gillespie,J.Chem.Phys.115,1716(2001).

[11]T.TianandK.Burrage,J.Chem.Phys.121,10356(2004).

[12]G.HenkelmanandH.J\'{o}nsson,J.Chem.Phys.115,9657(2001).

[13]G.HenkelmanandH.J\'{o}nsson,J.Chem.Phys.111,7010(1999).

[14]O.Trushin,A.Karim,A.KaraandT.S.Rahman,Phys.Rev.B72,115401(2005).

[15]M.A.GibsonandJ.Bruck,J.Phys.Chem.A104,1876(2000).

[16]A.Slepoy,A.P.ThompsonandS.J.Plimpton,J.Chem.Phys.128,205101(2008).资料来源量子化学网/content/%E5%8A%A8%E5%8A%9B%E5%AD%A6%E8%92%99%E7%89%B9%E5%8D%A1%E6%B4%9B%E6%96%B9%E6%B3%95kmc%E5%8F%8A%E7%9B%B8%E5%85%B3%E8%AE%A8%E8%AE%BA附录资料：不需要的可以自行删除教你如何保护笔记本HYPERLINK""一、每天关机前要做的清洗:双击“我的电脑”——右键点C盘——点“属性”——点“磁盘清理”——点“确定”——再点“是”——再点“确定”。清理过程中，您可看得到未经您许可（您可点“查看文件”看，就知道了）进来的“临时文件”被清除了，盘的空间多了。对D，E，F盘也要用这法进行。二、随时要进行的清理:打开网页——点最上面一排里的“工具”——点“Internet选项”——再点中间的“Internet临时文件”中的“删除文件”——再在“删除所有脱机内容”前的方框里打上勾——再点“确定”——清完后又点“确定”。这样，可为打开网和空间提高速度。三、一星期进行的盘的垃圾清理:点“开始”——用鼠标指着“所有程序”，再指着“附件”，再指着“系统工具”，再点“磁盘粹片整理程序”——点C盘，再点“碎片整理”（这需要很长时间，最好在您去吃饭和没用电脑时进行。清理中您可看到您的盘里的状况，可将清理前后对比一下）——在跳出“清理完成”后点“关闭”。按上述，对D，E，F盘分别进行清理。电脑系统越来越慢，怎么删除临时文件啊1.关闭"休眠"方法:打开[控制面板]→[电源选项]→[休眠],把"启用休眠"前面的勾去掉说明:休眠是系统长时间一种待机状态,使您在长时间离开电脑时保存操作状态,如果您不是经常开着电脑到别处去的话,那就把它关了吧!☆立即节省:256M2.关闭"系统还原"方法:打开[控制面板]→[系统]→[系统还原],把"在所有驱动器上关闭系统还原'勾上说明:系统还原是便于用户误操作或产生软件问题时的一种挽救手段,可以回复到误操作以前的状态.不建议初级用户使用.当然,它采用的是跟踪手段,需要记录大量信息,所消耗的资源也要很大的.☆立即节省:数百M(根据还原点的多少而不同)您也可以在不关闭系统还原的前提下,相应的减少系统还原所占的磁盘空间,这只会减少可用还原点的数目,一般还原点有一两个就够了吧.方法:...[系统还原]-选择一个"可用驱动器"-[设置]-调整"要使用的磁盘空间"3.关闭"远程管理"方法:打开[控制面板]→[系统]→[远程],把"允许从这台计算机发送远程协助邀请"前面的勾去掉.说明:谁会经常用到这种功能呢?它占用的不是磁盘空间,但是会影响系统运行速度.☆提高系统性能4.关闭"自动更新"方法:打开[控制面板]→[自动更新]→选"关闭自动更新"说明:系统更新占用的也不是磁盘空间,您可以在有可用更新时到微软的更新页面上更新,而不是总需要一个进程监视那是不是有可用更新.☆提高系统性能5.关闭"索引服务"方法:[开始]→[运行]→输入"msconfig"→[服务]去掉indexingservise前面的勾说明:索引服务是提高系统搜索速的的,谁没事总在系统里搜来搜去啊☆提高系统性能6.减少开机加载的启动项方法:[开始]→[运行]→输入"msconfig"→[启动]去掉不必要的启动项说明:必要的启动项一般"系统盘\windows\system32下的最好不要去掉;而肯定能去掉的有TINTSETPIMJPMIGIMSCMIGQQ等☆加快开机速度7.合理设置虚拟内存方法:打开[控制面板]→[系统]→[高级]-[性能]设置-[高级]-[虚拟内存]更改-在"自定义大小"中填入数值说明:一般的最小值是物理内存的1.5倍,最大值是物理内存的3倍;如果您经常运行大型软件的话最小值应该设为(物理内存*1.5*0.2+物理内存*1.5),最大值应为(物理内存*3*0.2+物理内存*3)☆提高系统性能8.取消系统失败时的日志文件和询问打开[控制面板]→[系统]→[高级]→[启动和故障恢复]→[系统失败]去掉"将事件写入系统日志""发送错误警报"☆提高系统性能9.删除无用的文件方法:打开[控制面板]→[文件夹选项]→[查看]→[显示所有文件和文件夹]C:\DocumentsandSettings\用户名\Cookies\除index文件外的所有文件(浏览网页产生的记录文件)C:\DocumentsandSettings\用户名\LocalSettings\Temp\下的所有文件（用户临时文件）C:\DocumentsandSettings\用户名\LocalSettings\TemporaryInternetFiles\下的所有文件(浏览网页记录的多媒体信息,加速以后浏览)C:\DocumentsandSettings\用户名\LocalSettings\History\下的所有文件（访问历史纪录）C:\DocumentsandSettings\用户名\Recent\下的所有文件（最近浏览文件的快捷方式）C:\WINDOWS\Temp\下的所有文件（临时文件）C:\WINDOWS\ServicePackFiles下的所有文件（升级sp1或sp2后的备份文件）C:\WINDOWS\DriverCache\i386下的压缩文件(驱动程序的备份文件)C:\WINDOWS\SoftwareDistribution\download下的所有文件(未完成的软件下载)C:\Windows\下以$u...开头的隐藏文件(微软更新安装文件备份)☆立即节省:数百M10.磁盘碎片整理下面该做的就是磁盘碎片整理了,只有整理之后才能获得更多的空间哦^_^方法:[开始]→[程序]→[附件]→[系统工具]→[磁盘碎片整理]...(记得先"分析"后"整理")二。发信人:liushafeng(终结者),信区:Hardware标题:减少WinXP资源占用的八大技巧发信站:紫金飞鸿(SunJun416:25:482006)WindowsXP被微软称为其历史上最优秀的操作系统，有让你眼花缭乱的各种功能、更快的速度，当然这一切都对计算机的硬件提出了更高的要求，如果你希望WindowsXP能够尽可能少地占用你有限的虽然WindowsXP被微软自称为有史以来最稳定、功能最强大的Windows操作系统，并且运行速度飞快——启动速度甚至比同样配置的Win2000还要快许多，你可能依然不满足于此，希望一、使用朴素界面XP安装后默认的界面包括任务栏、开始选单、桌面背景、窗口、按钮等都采用的是XP的豪华、炫目的风格，但缺点显而易见，它们将消耗掉不少系统资源，但实用意义不大。[方法]鼠标右键单击桌面空白处，在弹出选单点击“属性”进入显示属性设置窗口，将“主题、外观”都设置为“Windows经典”，将桌面背景设置为“无”，按确定保存退出。二、减少启动时加载项目许多应用程序在安装时都会自作主张添加至系统启动组，每次启动系统都会自动运行，这不仅延长了启动时间，而且启动完成后系统资源已经被吃掉不少![方法]选择“开始”选单的“运行”，键入“msconfig”启动“系统配置实用程序”，进入“启动”标，在此窗口列出了系统启动时加载的项目及来源，仔细查看你是否需要它自动加载，否则清除项目前的复选框，加载的项目愈少，启动的速度自然愈快。此项需要重新启动方能生效。三、优化视觉效果[方法]选择“系统属性”中的“高级”标签进入“性能选项”界面，其中“视觉效果”中可供选择的包括:自动设置为最佳、最佳外观、最佳性能、自定义。选中的效果越多则占用的系统资源越多，选定“最佳性能”项将关闭列表中列出诸如淡入淡出、平滑滚动、滑动打开等所有视觉效果。四、关闭系统还原默认情况下系统还原功能处于启用状态，每个驱动器约被占用高达4%～12%的硬盘空间，并且系统还原的监视系统会自动创建还原点，这样在后台运行就会占用较多的系统资源。[方法]鼠标右键点击桌面“我的电脑”中的“属性”进入“系统属性”设置窗口，选择“系统还原”标签，将“在所有驱动器上关闭系统还原”置为选中状态五、加快选单显示速度[方法]运行注册表编辑器，进入“HKEY_CURRENT_USERControlPanelDesktop”，将名称为MenuShowDelay的数据值由原来默认的400修改为0，修改后XP的开始选单、甚至应用软件的选单显示速度都会明显加快。六、启用DMA传输模式所谓DMA，即直接存储器存储模式，指计算机周边设备(主要指硬盘)可直接与内存交换数据，这样可加快硬盘读写速度，提高速据传输速率。[方法]选择“系统属性”中的“硬件”标签，打开“设备管理器”，其中“IDE控制器”有两项“PrimaryIDEChannel”及“SecondaryIDEChannel”，依次进入“属性→高级设置”，该对话框会列出目前IDE接口所连接设备的传输模式，点击列表按钮将“传输模式”设置为“DMA(若可用七、移动临时文件储存路径多数应用软件在运行时都会产生临时文件，而且这些临时文件都默认保存于启动分区C盘，长时间频繁读写C盘极易产生大量文件碎片，从而影响C盘性能，而C盘又是储存系统启动核心文件的分区，C盘的性能直接影响到系统的稳定性与运行效率。应尽量将应用软件安装于启动盘以外的分区并定期对硬盘进行整理，此举可最大程度避免产生磁盘碎片，将启动或读写速度保持在最佳状态。InternetExplorer临时文件夹[方法]在IE主窗口中，依次进入“工具→Internet选项→常规”标签，打开“Internet临时文件”设置界面，点击“移动文件夹”按钮将原来保存于C盘的临时目录移动至C盘以外的驱动器中，如果你使用的是宽带，可将“临时文件夹”使用空间设置为最小值1M刻录时产生的临时文件[方法]文件在刻录之前都会保存于C盘的刻录临时文件夹中，进入资源管理器，选择刻录机盘符并单击鼠标右键选单的“属性”项，在“录制”标签下可将此临时文件夹安置于其它驱动器。我的文档[方法]鼠标右键点击“我的文档”，在属性设置项中可将“我的文档”默认的保存路径修改至其它盘符。八、增加虚拟内存[方法]进入“性能选项”的“高级”设置窗口，首先将“处理器计划”及“内存使用”都调整为“程序”优化模式。点击“更改”按钮进入虚拟内存设置窗口，若你的内存大于256M，建议你禁用分页文件。默认的分页文件为物理内存的1.5倍。禁用系统缓存需重新启动系统。如果你的内存低于256M，请勿禁用分页文件，否则会导致系统崩溃或无法再启动XP!三windows使用时间长了，自然就会产生这样那样的临时文件，影响系统速度，让人烦心。对于老鸟来说，可以通过修改注册表、手动优化系统性能来提高速度，而对新手来说比较困难。这里我们可以下载一些专门的系统优化软件，进行一些简单的设计，就能达到我们想要的效果。这里我推荐一些常用又好使的软件：WINDOWS优化王、优化大师、超级兔子等。其中WINDOWS优化王非常好用，一看就会、功能全面、省时省心。完全清理[毒]垃圾藏在XP中的一个秘密武器，可以完整清除垃圾文件藏在XP中的一个秘密武器，可以完整清除垃圾文件你有用过Windows内置的清理磁盘功能吗？它并不能完全地清洗Windows内不需要的档案，因为它的功能被隐藏了，本篇将会把它被封印了的功能完全打开。适用的窗口板本除了Win95及Win98外这个方法阶适用于Win98se、Win2000、WinME、WinXP现在介绍两个「清理磁盘」工具的指令：SAGESET及SAGERUN首先在「开始」>「执行」然后输入cleanmgr/sageset:99设定:特别模式「清理磁盘」工具会执行，你会发觉多了很多清理选择，选择你想要清理的档案，通常全部都可以删除，完成你的选择后再按「确定」。然后再打开「开始」>「运行」输入：cleanmgr/SAGERUN:99杀毒1这种情况往往表现在打开IE时，在IE界面的左下框里提示：正在打开网页，但老半天没响应。在任务管理器里查看进程，（进入方法，把鼠标放在任务栏上，按右键—任务管理器—进程）看看CPU的占用率如何，如果是100%，可以肯定，是感染了病毒，这时你想运行其他程序简直就是受罪。这就要查查是哪个进程贪婪地占用了CPU资源．找到后，最好把名称记录下来，然后点击结束，如果不能结束，则要启动到安全模式下把该东东删除，还要进入注册表里，（方法：开始—运行，输入regedit）在注册表对话框里，点编辑—查找，输入那个程序名，找到后，点鼠标右键删除，然后再进行几次的搜索，往往能彻底删除干净。杀毒2今天在这里为大家提供两则小技巧，以便帮你强行杀死顽固病毒进程。根据进程名查杀这种方法是通过WinXP系统下的taskkill命令来实现的，在使用该方法之前，首先需要打开系统的进程列表界面，找到病毒进程所对应的具体进程名。接着依次单击“开始→运行”命令，在弹出的系统运行框中，运行“cmd”命令；再在DOS命令行中输入“taskkill/imaaa”格式的字符串命令，单击回车键后，顽固的病毒进程“aaa”就被强行杀死了。比方说，要强行杀死“conime。exe”病毒进程，只要在命令提示符下执行“taskkill/imconime。exe”命令，要不了多久，系统就会自动返回结果。根据进程号查杀上面的方法，只对部分病毒进程有效，遇到一些更“顽固”的病毒进程，可能就无济于事了。此时你可以通过Win2000以上系统的内置命令——ntsd，来强行杀死一切病毒进程，因为该命令除System进程、SMSS。EXE进程、CSRSS。EXE进程不能“对付”外，基本可以对付其它一切进程。但是在使用该命令杀死病毒进程之前，需要先查找到对应病毒进程的具体进程号。考虑到系统进程列表界面在默认状态下，是不显示具体进程号的，因此你可以首先打开系统任务管理器窗口，再单击“查看”菜单项下面的“选择列”命令，在弹出的设置框中，将“PID(进程标志符)”选项选中，单击“确定”按钮。返回到系统进程列表页面中后，你就能查看到对应病毒进程的具体PID了。接着打开系统运行对话框，在其中运行“cmd”命令，在命令提示符状态下输入“ntsd-cq-pPID”命令，就可以强行将指定PID的病毒进程杀死了。例如，发现某个病毒进程的PID为“444”，那么可以执行“ntsd-cq-p444”命令，来杀死这个病毒进程。杀毒3Windows任务管理器是大家对进程进行管理的主要工具，在它的“进程”选项卡中能查看当前系统进程信息。在默认设置下，一般只能看到映像名称、用户名、CPU占用、内存使用等几项，而更多如I/O读写、虚拟内存大小等信息却被隐藏了起来。可别小看了这些被隐藏的信息，当系统出现莫名其妙的故障时，没准就能从它们中间找出突破口。1.查杀会自动消失的双进程木马前段时间朋友的电脑中了某木马，通过任务管理器查出该木马进程为“system.exe”，终止它后再刷新，它又会复活。进入安全模式把c:＼windows＼system32＼system.exe删除，重启后它又会重新加载，怎么也无法彻底清除它。从此现象来看，朋友中的应该是双进程木马。这种木马有监护进程，会定时进行扫描，一旦发现被监护的进程遭到查杀就会复活它。而且现在很多双进程木马互为监视，互相复活。因此查杀的关键是找到这“互相依靠”的两个木马文件。借助任务管理器的PID标识可以找到木马进程。调出Windows任务管理器，首先在“查看→选择列”中勾选“PID（进程标识符）”，这样返回任务管理器窗口后可以看到每一个进程的PID标识。这样当我们终止一个进程，它再生后通过PID标识就可以找到再生它的父进程。启动命令提示符窗口，执行“taskkill/imsystem.exe/f”命令。刷新一下电脑后重新输入上述命令，可以看到这次终止的system.exe进程的PID为1536，它属于PID为676的某个进程。也就是说PID为1536的system.exe进程是由PID为676的进程创建的。返回任务管理器，通过查询进程PID得知它就是“internet.exe”进程。找到了元凶就好办了，现在重新启动系统进入安全模式，使用搜索功能找到木马文件c:＼windows＼internet.exe，然后将它们删除即可。前面无法删除system.exe，主要是由于没有找到internet.exe（且没有删除其启动键值），导致重新进入系统后internet.exe复活木马。2.揪出狂写硬盘的P2P程序单位一电