韦尔莱表

✍ dations ◷ 2025-08-05 19:26:06 #计算化学

韦尔莱表(Verlet table 或 Verlet list)是分子模拟中常用的一种减少粒子间距离计算量的方法,由法国物理学家卢普·韦尔莱(英语:Loup Verlet)首先提出。

分子模拟中,为减少计算量,通常为体系中每一个粒子规定一个“截断半径”,对于一个粒子,只有当某个其他粒子与其距离处于截断半径以内时才计算它们之间的相互作用。由于粒子间作用力通常都是短程力,这种近似广泛用于蒙特卡洛方法和分子动力学模拟中。然而,当模拟的体系进一步增大时,计算每两个粒子间的距离变得非常耗时,韦尔莱表应运而生。韦尔莱提出为每一个粒子建立一个列表,用来保存在它截断半径之内的其他粒子的编号,这个列表就称为韦尔莱表。为使韦尔莱表不必每个模拟步长都需要更新,韦尔莱表的构建引入“第二截断半径”'Rv'大于粒子的截断半径'Rc'。例如,对于蒙特卡洛方法,此值为 R v = R c + 2 n d {\displaystyle R_{v}=R_{c}+2nd} ,其中 n {\displaystyle n} 为韦尔莱表更新步长间隔, d {\displaystyle d} 为一步中粒子的最大移动距离,以此保证所有应当计算的粒子都得到统计。更新韦尔莱表的时间复杂度为 O ( N 2 ) {\displaystyle O(N^{2})} (N为粒子总数),对于蒙特卡洛方法经优化可达到 O ( N 5 / 3 ) {\displaystyle O(N^{5/3})}

以下是以Fortran描述的构建韦尔莱表的算法。

subroutine new_listdo i = 1 , npart ! 初始化列表,npart为体系中粒子总数    nlist(i) = 0    xv(i) = x(i)end dodo i = 1 , npart - 1    do j = i + 1 , npart ! 遍历所有粒子对        xr = x(i) - x(j) ! 计算两粒子距离        call period_condition(xr) ! 依周期性边界条件校正粒子距离        if(abs(xr) .lt. rv) then ! 找到符合条件的粒子对            ! 往韦尔莱表中添加信息            nlist(i) = nlist(i) + 1            nlist(j) = nlist(j) + 1 ! MC模拟中每个粒子独自考虑,故ij均保留完全的列表。而MD中可只保留粒子i的列表,粒子j的作用力由牛顿第三定律求算。            list(i,nlist(i)) = j            list(j,nlist(j)) = i        end if    end doend do

不足与改进

韦尔莱表的 O ( N 2 ) {\displaystyle O(N^{2})} 复杂度使其在体系增大时耗时骤增,直至成为整个模拟中最耗时的步骤。在更大的体系时,通常采用“元胞列表”(Cell lists)的方法,其复杂度为 O ( N ) {\displaystyle O(N)} 。这两种方法的结合能进一步提高计算效率。

相关

  • 共价化合物共价键(英语:covalent bond),是化学键的一种。两个或多个非金属原子共同使用它们的外层电子(砷化镓为例外),在理想情况下达到电子饱和的状态,由此组成比较稳定和坚固的化学结构叫做
  • 乳酸钠乳酸钠(英语:Sodium lactate)是乳酸的钠盐,常温下为白色粉末,具有较淡的咸味。它是通过含糖类的植物(如玉米和甜菜)发酵之后,将产生的乳酸中和而制得的。早在1836年,人们已经知道乳酸
  • 招待局中国人民解放军军徽中央军委机关事务管理总局招待局,位于北京市,是中央军事委员会机关事务管理总局下属局,负责中央军委机关的招待工作。2016年改革前,中国人民解放军总参谋部管
  • 苗栗span class=inline-unihan style=border-bottom: 1px dotted; font-variant: norm苗栗
  • 徐小明徐小明可能指以下其中一人:
  • 石龙子亚目石龙子亚目(Scincomorpha)是有鳞目的一个亚目。Polyglyphanodon sternbergiEutropis macularia
  • 新加坡图书馆列表这是一个位于新加坡的图书馆列表:这些图书馆由国家图书馆管理局负责管理。
  • 克莱蒙费朗克莱蒙费朗(法语:Clermont-Ferrand,法语发音:.mw-parser-output .IPA{font-family:"Charis SIL","Doulos SIL","Linux Libertine","Segoe UI","Lucida Sans Unicode","Code2000"
  • 基督教赉恩堂基督教赉恩堂,位于中国广东省佛山市禅城区莲花路27号,为佛山市市级文物保护单位,类型为近现代重要史迹及代表性建筑,公布时间为1998年4月30日。基督教赉恩堂的历史年代为民国。
  • 我的智能情人《我的智能情人》(韩语:나 홀로 그대,英语:My Holo Love),2月7日正式上线,由《认识的妻子》的李尚烨导演与《精神病患者日记》的刘勇在作家合作打造,为2020年Netflix韩国原创剧集。