韦尔莱表

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

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

相关

  • 义位义位(英语:Sememe)是语义学中能够独立存在、运用的最小意义单位。义位是形成语义场的核心。义位可以是一个语素所表达的语义,比如英语的复数语素-s ,对应的义位即其语义为;也可以
  • 拉斐尔房间四间拉斐尔客房(Stanze di Raffaello)是梵蒂冈宫内的一组客房,教宗住所的公共部分。它们以拉斐尔及其工作室创作的壁画著称。它们连同西斯廷礼拜堂内米开朗基罗的天顶画,构成标
  • 百立建设百立建设股份有限公司(简称百立建设)成立于1982年11月17日,每个建案皆能获得建筑业界的奖项,是高雄知名的豪宅品牌。负责人为邱祐勋。曾获ISO9002品质认证及内政部优良企业识别
  • 大专院校高等教育(法语:Études Supérieures;德语:Höhere Bildung;英语:Higher Education;西班牙语:Educación Superior),又称专上教育(英语:Post-secondary Education)是一个教育层级的概念,广
  • 台北府台北府为清廷于光绪元年12月20日(公元1876年1月16日),沈葆祯奏请在台湾岛北部设置之官署及行政区划。于官署方面,台北府设置台北府知府等官员,衙门地点则位于台北城内,与台湾巡抚
  • 20112011年伊朗人口普查是伊朗自1956年以来的第7次人口普查,普查标准日在2011年(伊朗历1390年)10月24日。
  • 赫伯特·弗勒利希赫伯特·弗勒利希,FRS(德语:Herbert Fröhlich,1905年12月9日-1991年1月23日),德国出生的英国物理学家。
  • 台北市政府消防局防灾科学教育馆台北市政府消防局防灾科学教育馆位于台湾台北市内湖区,为台湾第一座防灾科学教育馆,隶属台北市政府消防局。台北市政府消防局在1995年改制成立后,对于容易在台北造成的灾害加以
  • 1,5-二氨基四唑1,5-二氨基四唑是一种有机化合物,可以二氨基胍盐酸盐为原料制备。1,5-二氨基四唑是一种配位能力较强的多齿配体,可以和金属形成配合物,如:它也可以和强酸反应并成盐,如盐酸盐、高
  • 克里斯琴·戈特弗里德·丹尼尔·尼斯·冯埃森贝克克里斯琴·戈特弗里德·丹尼尔·尼斯·冯埃森贝克(Christian Gottfried Daniel Nees von Esenbeck,1776年2月14日-1858年3月16日)为德国植物学家。他描述了约7000种植物(几乎和林