【LAMMPS学习】八、基础知识(2.8)长程色散设置

8. 基础知识

此部分描述了如何使用 LAMMPS 为用户和开发人员执行各种任务。术语表页面还列出了 MD 术语,以及相应 LAMMPS 手册页的链接。 LAMMPS 源代码分发的 examples 目录中包含的示例输入脚本以及示例脚本页面上突出显示的示例输入脚本还展示了如何设置和运行各种模拟。

8.1.通用基础知识

8.2. 设置入门

8.2.1.二维模拟

8.2.2.类型标签

8.2.3.三斜(非正交)模拟盒子

8.2.4.恒温器

8.2.5.恒压器 

8.2.6.Walls墙

8.2.7.NEMD 模拟

8.2.8.长程色散设置

PPPM 方法通过将对势分成两部分来计算相互作用,其中一部分以正常的成对方式计算,即所谓的实空间部分,而其中一部分则使用傅立叶变换计算,即所谓的倒数部分。 space 或 kspace 部分。对于这两个部分,电势并未精确计算,而是近似计算。因此,计算的两个部分都存在误差,即实空间误差和k空间误差。刚才提到的事实对于库仑的 PPPM 以及色散相互作用都是正确的。决定性的差异 - 也是必须更加谨慎选择 pppm/disp 参数的原因 - 是误差对结果的影响:库仑和色散相互作用的 PPPM 的 kspace 误差以及实空间误差库仑相互作用的 PPPM 具有噪声特征。相比之下,PPPM 的色散实空间误差有一个明确的物理解释:内聚力的低估。因此,实空间误差对 pppm/disp 模拟结果的影响比 k 空间误差强得多。因此,必须以该误差远小于 kspace 误差的方式选择参数。

当使用 pppm/disp 且未通过 kspacemodify 命令对 PPPM 参数进行任何指定时,将调整参数以使实际空间误差和 kspace 误差相等。这将导致模拟不准确或缓慢,这两者都是不可取的。为了选择可提供快速、准确模拟的 pppm/disp 参数,有两种方法,它们都有各自的优点和缺点。

第一种方法是通过 kspace_modifyforce/disp/real 和 kspace_modifyforce/disp/kspace 命令设置所需的实空间和 kspace 精度。请注意,精度必须以力单位指定,因此取决于所选的单位设置。对于实数单位,0.0001 和 0.002 似乎可以为实空间和 k 空间精度提供合理、准确且高效的计算。 0.002 和 0.05 对于大多数使用 lj 单位的系统来说效果很好。 PPPM 参数将根据所需的精度生成。这种方法的优点是它通常提供一组很好的参数,并且适用于 kspace_modify diff ad 和 kspace_modify diff ik 选项。该方法的缺点是在模拟初始化期间设置 PPPM 参数将花费一些时间。

第二种方法是使用 kspace_modify mesh/disp、kspace_modify order/disp 和 kspace_modify gewald/disp 命令显式设置 pppm/disp 的参数。这种方法需要经验丰富的用户,他们充分了解参数选择对模拟精度和性能的影响。这种方法提供了模拟的快速初始化。然而,它对错误很敏感:对于一个系统表现良好的参数组合可能会导致其他模拟的条件远非最佳。例如,为全原子力场提供准确和快速计算的参数可能会提供不够精确的精度或联合原子力场(这与后者通常具有较大的色散系数有关)。

为了避免不准确或低效的模拟,如果未采取任何措施来控制 PPPM 参数,pppm/disp 会停止模拟并显示错误消息。如果需要自动生成参数并且希望实空间和 kspace 精度相等,则可以使用 kspace_modify disp/auto yes 命令抑制此错误消息。

结合这两种方法优点的合理方法是使用 kspace_modify force/disp/real 和 kspace_modify force/disp/kspace 命令进行第一次运行,记下输出中的 PPPM 参数,并使用第二种方法指定这些参数在后续运行中(具有相同的成分、力场和大致相同的体积)。

关于 pppm/disp 的性能,还有两件事需要考虑。

首先是当使用 pppm/disp 时,cutoff 参数不再影响模拟的精度(以更改 cutoff 时调整 gewald/disp 为准)。因此,可以通过检查截止参数的不同值来提高性能。截止的下限仅由势对排斥项的截断误差设置。

第二个是当使用 pppm/disp 时,pair 风格的混合规则对计算时间有影响。使用几何混合规则时可以实现最快的计算。使用算术混合规则大大增加了计算成本。使用 kspace_modify mix/disp geom 和 kspace_modify splittol 命令可以减少计算开销。第一个命令只是在 kspace 计算中强制执行色散系数的几何混合。这会在计算中引入一些错误,但也会显着加快模拟速度。第二个关键字设置使用矩阵分解方法近似分散系数的精度。与使用第一个命令相比,这可能会带来更好的准确性,但通常也不会提供同样好的效率提高。

最后,当不适用混合规则时,也可以使用 pppm/disp。这可以使用 kspace_modify mix/disp none 命令来实现。请注意,代码不会自动检查是否满足任何混合规则。如果混合规则不适用,用户将必须明确指定此命令。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.hqwc.cn/news/618980.html

如若内容造成侵权/违法违规/事实不符,请联系编程知识网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

【opencv】示例-train_svmsgd.cpp 随机梯度下降支持向量机(SVMSGD)对二维数据进行分类的UI...

#include "opencv2/core.hpp" // 引入opencv2的核心头文件 #include "opencv2/video/tracking.hpp" // 引入opencv2视频跟踪相关功能的头文件 #include "opencv2/imgproc.hpp" // 引入opencv2的…

什么是 MVVM、mvc 模型

mvc模型 MVC: MVC 即 model-view-controller(模型-视图-控制器)是项目的一种分层架构思想,它把复杂的业务逻辑, 抽离为职能单一的小模块,每个模块看似相互独立,其实又各自有相互依赖关系。它的好处是:保证了…

微服务相关

1. 微服务主要七个模块 中央管理平台:生产者、消费者注册,服务发现,服务治理,调用关系生产者消费者权限管理流量管理自定义传输协议序列化反序列化 2. 中央管理平台 生产者A在中央管理平台注册后,中央管理平台会给他…

移动Web学习06-移动端适配Less预处理器项目案例

项目目标:实现在不同宽度设备中等比缩放的网页效果 Less代码 import ./base; import ./normalize;// 变量: 存储37.5 rootSize: 37.5rem; *{margin: 0;padding: 0; } body {background-color: #F0F0F0; }// 主体内容 .main {// padding-bottom: (50 / 37.5rem);pa…

【uniapp踩坑记】——微信小程序转发保存图片

关于微信小程序转发&保存图片 微信小程序图片转发保存简单说明网络图片的转发保存base64流形式图片转发保存 已经好多年没写博客了,最近使用在用uniapp开发一个移动版管理后台,记录下自己踩过的一些坑 吃相别太难看,搞一堆下头僵尸号来点…

Spring Cloud学习笔记:Eureka集群搭建样例

这是本人学习的总结,主要学习资料如下 - 马士兵教育 1、项目架构2、Dependency3、项目启动类4、application.yml5、启动项目 1、项目架构 因为这是单机模拟集群搭建,为了方便管理就都放在了一个项目中。这次准备搭建三个项目server1, server2, server3 …

Linux网络基础 (二) ——(IP、MAC、端口号、TCPUDP协议、网络字节序)

文章目录 IP 地址基本概念源IP地址 & 目的IP地址 MAC 地址基本概念源MAC地址 & 目的MAC地址 端口号基本概念源端口号 & 目的端口号 TCP & UDP 协议基本概念TCP 与 UDP 的抉择 网络字节序大端、小端字节序 🎖 博主的CSDN主页:Ryan.Alask…

攻防世界---Web_php_include

1.题目链接 2.补充知识&#xff1a; 3.构造&#xff1a;执行成功 /?pagedata://text/plain,<?php phpinfo()?> 4.构造下面url&#xff0c;得到目录路径 /?pagedata://text/plain,<?php echo $_SERVER[DOCUMENT_ROOT]?> 5构造下面url&#xff0c;读取该路径的…

Linux的学习之路:10、进程(2)

摘要 本章主要是说一下fork的一些用法、进程状态、优先级和环境变量。 目录 摘要 一、fork 1、fork的基本用法 2、分流 二、进程状态 三、优先级 四、环境变量 1、常见环境变量 2、和环境变量相关的命令 3、通过代码如何获取环境变量 五、导图 一、fork 1、fork…

微信小程序实现预约生成二维码

业务需求&#xff1a;点击预约按钮即可生成二维码凭码入校参观~ 一.创建页面 如下是博主自己写的wxml&#xff1a; <swiper indicator-dots indicator-color"white" indicator-active-color"blue" autoplay interval"2000" circular > &…

Zookeeper和Kafka的部署

目录 一、Zookeeper的基本概念 1. Zookeeper定义 2. Zookeeper工作机制 3. Zookeeper特点 4. Zookeeper数据结构 5. Zookeeper应用场景 5.1 统一命名服务 5.2 统一配置管理 5.3 统一集群管理 5.4 服务器动态上下线 5.5 软负载均衡 6. Zookeeper 选举机制 6.1 第一…

Cortex-M3/M4处理器的bit-band(位带)技术

ARM Cortex-M3/M4的位带&#xff08;Bit-Band&#xff09;技术是一种内存映射技术&#xff0c;它允许对单个位进行直接操作&#xff0c;而不需要对整个字&#xff08;通常是32位&#xff09;进行操作。这项技术主要用于对特定的位进行高效的读写&#xff0c;特别是在需要对GPIO…