韦尔莱表

✍ dations ◷ 2025-04-02 17:15:20 #计算化学

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

相关

  • 屈折变化在语法学中,词形变化(又译作形态变化、屈折变化)(Inflection or inflexion)指单词(或词根)的变化,以导致语法功能改变,进而使其代表的意义也有所改变。印欧语屈折变化又可以分为变位
  • 迦太基迦太基(腓尼基语:
  • 俄罗斯宪法俄罗斯联邦宪法(俄语:Конституция Российской Федерации, Konstitutsiya Rossiyskoy Federatsii; 发音.mw-parser-output .IPA{font-family:"C
  • 原发性肝细胞癌肝细胞癌 (Hepatocellular carcinoma,简称HCC)是成年人中最常见类型的慢性肝癌,并且在肝硬化的患者中也是最常见的死亡原因。它发生在慢性肝脏炎症的环境中,并且与慢性病毒性肝
  • 克伦特罗克伦特罗(Clenbuterol),是一种β2-肾上腺素受体促效剂(β2-adrenergic agonist),类似麻黄素(Ephedrine)作用,临床上经常用来治疗慢性阻塞性肺疾(COPD),亦被作为缓和气喘急性发作时的支气
  • 大稻埕圆环防空蓄水池坐标:25°03′14.11″N 121°30′52.39″E / 25.0539194°N 121.5145528°E / 25.0539194; 121.5145528大稻埕圆环防空蓄水池是第二次世界大战中太平洋战争(日本称之为大东亚
  • 交换式多兆位数据服务交换式多兆位数据服务(英语:Switched Multi-megabit Data Service,缩写SMDS)是一项1990年代早期用于连接LAN、MAN和WAN以交换数据的无连接服务。在欧洲,该服务被称为无连接宽带数
  • 费迪南·艾森斯坦费迪南·哥德霍尔特·马克斯·艾森斯坦(Ferdinand Gotthold Max Eisenstein,1823年4月16日-1852年10月11日),德国数学家。艾森斯坦的数学才华很早便显露了,他在一篇自传式的文章写
  • 伍德维尔 (阿拉巴马州)伍德维尔(英文:Woodville),是美国阿拉巴马州下属的一座城市。面积约为6.62平方英里(约合 17.15平方公里)。根据2010年美国人口普查,该市有人口746人,人口密度为112.65/平方英里(约合4
  • 峇株巴辖县峇株巴辖县(马来语:Daerah Batu Pahat),是马来西亚柔佛州西部的一个县。其面积为1872.58平方公里,人口于2010年为417,458。该县北临昔加末县,西北临麻坡县,东临居銮县,东南临笨珍县,