韦尔莱表

✍ dations ◷ 2025-11-06 22:16:22 #计算化学

韦尔莱表(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)} 。这两种方法的结合能进一步提高计算效率。

相关

  • 两生类两栖动物(学名:Amphibia)是两栖纲生物的通称,又名两生动物,包括所有生没有卵壳的卵,拥有四肢的脊椎动物(蚓螈的四肢已退化)。两栖动物的皮肤裸露,表面没有鳞片、毛发等覆盖,但是可以分
  • 瓦特瓦特(符号:W)是国际单位制的功率单位。瓦特的定义是1焦耳/秒(1 J/s),即每秒钟转换,使用或耗散的(以安培为量度的)能量的速率。日常生活中更常用千瓦作为单位,1千瓦=1000瓦特,千瓦又可合
  • EcoRVEcoRV(V是5的意思)是一种第二型限制内切酶,来自某些特定品系的大肠杆菌(Escherichia coli),此细菌的学名也是EcoRV前三个字母的由来。EcoRV是分子生物学与生物技术上常使用的限制
  • 食品加工业食品加工业是将食物透过物理或化学途径转化为其他形态的食物产业,旨于令食材更易包装及烹饪。典形的食物加工包括切碎、浸渍、液化、乳化、烹调(例如煮、烤、煎、烧)、腌制,甚至
  • 基奈半岛基奈半岛是美国的半岛,位于阿拉斯加州南岸,长245公里、宽190公里,最高点海拔高度2,130米,该地区蕴藏天然气、石油和煤等天然资源。坐标:60°10′N 150°15′W / 60.167°N 150.25
  • 奥克-罗曼语支奥克-罗曼语支 (加泰罗尼亚语:llengües occitanoromàniques,奥克语:lengas occitanoromanicas),又称东伊比利亚语支,是位于罗曼语族下西罗曼语支的一个分支,一般指奥克语和加泰罗
  • 阿尔达班四世阿尔达班四世(?-224年或226年或227年)是帕提亚帝国(安息)的最后一位国王(约213年-约224年在位)。也有一些文献把他列为阿尔达班五世。阿尔达班四世是沃洛吉斯五世的幼子,在击败了兄长
  • 人民激进党人民激进党(塞尔维亚语:Народна радикална странка / Narodna radikalna stranka,HPC/NRS),是塞尔维亚王国,及塞尔维亚人、克罗地亚人和斯洛文尼亚人王国(
  • 裕全裕全(?-1840年),满洲爱新觉罗氏。豫通亲王多铎后裔、豫良亲王修龄第五子、已革豫亲王裕丰、裕兴之弟。乾隆五十二年(1787年),父亲修龄逝世,其长兄裕丰袭豫亲王爵位。嘉庆四年(1799年),裕
  • 老君 (马来世界)老君本指印度尼西亚和马来西亚的一类土医或巫医,马、新与印尼的闽南方言(包括潮州话)中借用了这一词汇,成为“医生”的代名词。在马来世界里,他们的社会角色为传统的治疗师、灵媒