韦尔莱表

✍ dations ◷ 2025-04-26 11:57:58 #计算化学

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

相关

  • 双及物动词在语法学上,双及物动词是需要支配一个主词及两个受词的动词。此与支配单个受词的单及物动词相对。1.双及物结构的语义表达研究,外语教学与研究,2009.1
  • 中石器时代中国中石器时代的考古工作远没有旧石器时代与新石器时代丰硕。中石器时代处于旧石器时代晚期与新石器时代早期之间,其特征是细石器的制造,所以又称作“细石器时代”。中国在约
  • Wacker法瓦克法(Wacker process),又称Hoechst-Wacker法,最早是指乙烯在含有四氯钯酸盐催化剂的水中,被空气中的氧气氧化为乙醛的反应。这是第一个工业化的有机金属(有机钯)反应,亦是均相催化
  • 美国食品药物管理署美国食品药品监督管理局(英语:U.S. Food and Drug Administration,缩写为FDA)为美国卫生与公众服务部直辖的联邦政府机构,其主要职能为负责对美国国内生产及进口的食品、膳食补充
  • 菲律宾比索菲律宾比索(菲律宾语:Piso,ISO 4217货币编码:PHP)是菲律宾的法定货币,1菲律宾比索相等于100 centavos(菲律宾语:sentimo)。菲律宾是世界上其中一个曾被西班牙殖民统治而使用比索作为
  • 面部 (部首)面部,为汉字索引中的部首之一,康熙字典214个部首中的第一百七十六个(九划的则为第一个)。就繁体和简体中文中,面部归于九划部首。面部只以左方、下方为部字。且无其他部首可用者
  • 轨道根数轨道根数(或称轨道要素或轨道参数)是描述在牛顿运动定律和牛顿万有引力定律的作用下的天体或航天器,在其开普勒轨道上运动时,确定其轨道所必要的六个参数。由于运动的方式有许多
  • 社会主义者和民主人士进步联盟S&D(2009年6月23日起) PES(1993年4月21日-2009年6月22日) SOC(1958年-1993年4月21日) S&D(2009年6月23日起) PSE(1993年4月21日-2009年6月22日) SOC(1958年-1993年4月21日) Group of the Pro
  • 2014年国际足联世界杯外围赛 – 欧洲区G组2014年世界杯外围赛欧洲区G组(英语:2014 FIFA World Cup qualification UEFA Group G)是欧洲足协辖下2014年世界杯欧洲赛区外围赛的一个分组,包括希腊、斯洛伐克、波斯尼亚和黑
  • 新加坡环境及水源部新加坡环境及水源部(马来语:Kementerian Alam Sekitar dan Sumber Air;英语:Ministry of the Environment and Water Resources (MEWR))是新加坡政府辖下的部门,原名为新加坡环境