【LAMMPS学习】八、基础知识(2.4)恒温器

8. 基础知识

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

8.1.通用基础知识

8.2. 设置入门

8.2.1.二维模拟

8.2.2.类型标签

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

8.2.4.恒温器

恒温是指在 MD 模拟中控制颗粒的温度。恒压调节 Barostatting 意味着控制压力。由于压力包括由粒子速度引起的动力学分量,因此这两种操作都需要计算温度。通常,目标温度 (T) 和/或压力 (P) 由用户指定,恒温器或恒压器尝试将系统平衡至所需的 T 和/或 P。

LAMMPS 中的恒温是通过 fixes来执行的,或者在一种情况下通过配对方式来执行。有几种恒温 fixes可用:Nose-Hoover (nvt)、Berendsen、CSVR、Langevin 和直接重新调整(temp/rescale)。耗散粒子动力学 (DPD) 恒温可以通过 dpd/tstat 对样式调用:

  • fix nvt

  • fix nvt/sphere

  • fix nvt/asphere

  • fix nvt/sllod

  • fix temp/berendsen

  • fix temp/csvr

  • fix langevin

  • fix temp/rescale

  • pair_style dpd/tstat

Fix nvt只有恒温器可以控制粒子的平移速度。 Fix nvt/sllod 也可以做到这一点,只不过它减去了由于变形框而产生的速度偏差并积分了 SLLOD 运动方程。有关更多详细信息,请参阅  Howto nemd 页面。Fix nvt/sphere and fix nvt/asphere  不仅可以调节球形和非球形粒子的平移速度,还可以调节旋转速度。

注:

Daivis 和 Todd 2017 年 出版的一本书讨论了 SLLOD 方法和非平衡 MD (NEMD) 恒温的一般用途,适用于简单流体和复杂流体,例如流体。分子系统。后者可能很难正确完成。

DPD 恒温以类似于fix langevin 的每粒子恒温的方式改变成对相互作用。

任何thermostatting fixes都可以指示使用自定义温度计算来消除偏差,这有两个影响:首先,与请求的目标温度进行比较的当前计算温度是在消除速度偏差的情况下计算的;其次,恒温器仅调整粒子速度的热温度分量,即去除偏置后的速度。然后将移除的偏差添加回调整后的速度。请参阅各个修复程序的文档页面,以及有关如何将温度计算分配给恒温修复程序的说明的  fix_modify 命令的说明。

例如,您可以将恒温器与 compute temp/region结合使用,仅将其应用于空间区域中的原子。或者,您可以将恒温器与 compute temp/partial结合使用,仅将其应用于速度的 x 和 z 分量。通过使用 compute temp/profile,您可以仅恒温颗粒流的热温度,而不影响流速度。

下面是可以这样使用的自定义温度计算列表:

  • compute temp/asphere command

  • compute temp/body command

  • compute temp/chunk command

  • compute temp/com command

  • compute temp/deform command

  • compute temp/partial command

  • compute temp/profile command

  • compute temp/ramp command

  • compute temp/region command

  • compute temp/rotate command

  • compute temp/sphere command

注:

只有nvt fixes 执行时间积分,这意味着它们分别更新由于力和速度而导致的粒子的速度和位置。另一个恒温器只修复调节速度;他们不执行时间积分更新。因此,它们应该与恒定的 NVE 集成修复结合使用,例如:

  • fix nve

  • fix nve/sphere

  • fix nve/asphere

热力学输出可以通过thermo_style 命令设置,通常包括温度值。正如thermo_style 命令页面上所解释的,默认温度是由thermo命令本身设置的。它不是与您定义的任何恒温修复相关的温度,也不是与您定义的计算温度的任何计算相关的温度。恒温修复的文档页面解释了他们创建的温度计算的 ID。因此,如果您想查看这些温度,则需要通过  thermo_style custom 命令显式指定它们。或者您可以使用 thermo_modify 命令重新定义默认热力学输出使用的温度计算。

Ref:

(Daivis and Todd) Daivis and Todd, Nonequilibrium Molecular Dynamics (book), Cambridge University Press, Nonequilibrium Molecular Dynamics, (2017).

8c04c0a551554cabb26f1a7098a8d179.png

 

 

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

如若内容造成侵权/违法违规/事实不符,请联系我们进行投诉反馈qq邮箱809451989@qq.com,一经查实,立即删除!

相关文章

【JavaWeb】Day38.MySQL概述——数据库设计-DQL(一)

数据库设计——DQL 介绍 DQL英文全称是Data Query Language(数据查询语言),用来查询数据库表中的记录。 查询关键字:SELECT 查询操作是所有SQL语句当中最为常见,也是最为重要的操作。在一个正常的业务系统中,查询操作的使用频次…

2024妈妈杯Mathorcup数学建模竞赛选题建议

关于十四届妈妈杯题目点击链接加入群聊【2024年第十四届MathorCup数学建模】:http://qm.qq.com/cgi-bin/qm/qr?_wv1027&kNol416eFZyg4AFPHCZsYfGkLnEnZ8H36&authKeyo9lIFgRfnk2U%2FfYYg1AlFJd5I456CSA2M6nlvJX2qcUockIKuMHj%2BhKMP6RnPeo1&noverify0&a…

在Ubuntu上搭建Prometheus + Grafana监控系统

1.Prometheus 部署 从官网下载页面找到最新的二进制文件下载 cd ~ curl -LO https://github.com/prometheus/prometheus/releases/download/v2.51.1/prometheus-2.51.1.linux-amd64.tar.gz将文件解压到指定目录 tar xf prometheus-2.51.1.linux-amd64.tar.gz -C /usr/local为…

每天学点儿Python(5) -- 序列索引和切片

Python中,序列是指一块可存放多个值的连续内存空间,这些值按一定顺序排列,可通过每个值所在位置的编号(称为索引)访问它们。它类似于C/C中的数组或字符串,但又比数组或字符串强大很多 序列类型包括字符串、…

挖掘未来:私有LTE/5G网络驱动智慧矿山的自动化

私有LTE/5G网络为世界上一些最偏远的角落提供无线连接。如果没有无线通信网络,各行业就无法满足增加产量、降低运营成本和减少环境破坏的需求。 在本案例研究中,我们着眼于自动化如何改变无线网络的动态。智慧矿山要求运营商无缝集成多个系统和应用程序…

Keil #include “stm32f10x.h“报错

给我的提示 我看到懵逼了,我就一直在网上找各种方法,发现都不行,我就想没有找到,我挨个挨个去找它想要的,发现是有的 我就想这不是有吗,怎么会找不到呢 我就又想是我路径写错了?我就看了一下路…

【深度学习】Fine-Grained Face Swapping via Regional GAN Inversion高保真换脸范式

文章目录 代码介绍实践效果 帮助、问询 代码 https://github.com/e4s2022/e4s 介绍 Fine-Grained Face Swapping via Regional GAN Inversion 提出一种新的高保真换脸范式,能够保留期望的微妙几何和纹理细节。从微观面部编辑的角度重新思考换脸任务,基…

npm创建Vue3项目

npm创建Vue3项目 1 创建Vue项目说明 2 安装3 运行 1 创建Vue项目 创建最新版的Vue项目,已经不推荐使用CLI构建方式了。参考如下即可。 npm create vuelatest如果发现一直动不了,切换网络试一下,个人热点尝试一下。 按下图的选项按需引入自…

[STM32+HAL]DengFOC移植之闭环位置控制

一、源码来源 DengFOC官方文档 二、HAL库配置 1、开启硬件IIC低速模式 低速更稳定 2、PWM波开启 三、keil填写代码 1、AS5600读取编码器数值 #include "AS5600.h" #include "math.h"float angle_prev0; int full_rotations0; // full rotation trac…

OSCP靶场--Hawat

OSCP靶场–Hawat 考点(目录扫描zip文件下载java代码审计web1的sql注入写到web2的webshel) 1.nmap扫描 ## ┌──(root㉿kali)-[~/Desktop] └─# nmap 192.168.158.147 -sV -sC -Pn --min-rate 2500 -p- Starting Nmap 7.92 ( https://nmap.org ) at 2024-04-11 07:18 EDT …

性能优化-01

当看到性能指标时,你会首先想到什么呢?我相信 “高并发” 和 “响应快” 一定是最先出现在你脑海里的两个词,而它们也正对应着性能优化的两个核心指标—— “吞吐” 和 “延时” 。这两个指标是从应用负载的视角来考察性能,直接影…

python爬虫----BeautifulSoup(第二十天)

🎈🎈作者主页: 喔的嘛呀🎈🎈 🎈🎈所属专栏:python爬虫学习🎈🎈 ✨✨谢谢大家捧场,祝屏幕前的小伙伴们每天都有好运相伴左右,一定要天天…

组合数学<1>——组合数学基础

今天我们聊聊组合数学。(本期是给刚刚学习组合数学的同学看的&#xff0c;dalao们可以自行忽略) 建议:不会求逆元的出门左转数论<2>&#xff0c;不会数论的出门右转数论<1>。 加乘原理 加乘原理小学奥数就有。 总的来说:加法原理:分类;乘法原理:分步 比如说&a…

景芯2.5GHz A72训练营dummy添加(一)

景芯A72做完布局布线之后导出GDS&#xff0c;然后进行GDS merge&#xff0c;然后用Calibre对Layout添加Dummy。在28nm以及之前的工艺中&#xff0c;Dummy metal对Timing的影响不是很大&#xff0c;当然Star RC也提供了相应的解决方案&#xff0c;可以考虑Dummy metal来抽取RC。…

【Vector-Map-路径规划(0)】卷首语

因为城市NOA 的开发过程中&#xff0c;十字路口这类场景非常不好处理&#xff0c;个人对路径规划没有什么基础&#xff0c;只知道深度优先&#xff0c;广度优先&#xff0c;A*&#xff0c;Dijkstra等算法&#xff0c;不知道在矢量地图中如何使用&#xff1f;因此花几天时间读几…

【LangChain系列】2. 一文全览LangChain数据连接模块:从文档加载到向量检索RAG,理论+实战+细节

本文学习 LangChain 中的 数据连接&#xff08;Retrieval&#xff09; 模块。该模块提供文档加载、切分&#xff0c;向量存储、检索等操作的封装。最后&#xff0c;结合RAG基本流程&#xff0c;我们将利用LangChain实现RAG的基本流程。 0. 模块介绍 在前面文章中我们已经讲了…

ssm042在线云音乐系统的设计与实现+jsp

在线云音乐系统的设计与实现 摘 要 随着移动互联网时代的发展&#xff0c;网络的使用越来越普及&#xff0c;用户在获取和存储信息方面也会有激动人心的时刻。音乐也将慢慢融入人们的生活中。影响和改变我们的生活。随着当今各种流行音乐的流行&#xff0c;人们在日常生活中经…

MySQL 连接查询

目录 连接查询 命令格式&#xff1a; 内连接&#xff1a; 等值连接&#xff1a; 格式&#xff1a; 非等值连接&#xff1a; 格式&#xff1a; 外连接&#xff1a; 左连接&#xff1a; 格式: 结果&#xff1a; 右连接&#xff1a; 格式: 结果&#xff1a; 全外连…

D-LinkNAS 远程命令执行漏洞(CVE-2024-3273)RCE漏

声明&#xff1a; 本文仅用于技术交流&#xff0c;请勿用于非法用途 由于传播、利用此文所提供的信息而造成的任何直接或者间接的后果及损失&#xff0c;均由使用者本人负责&#xff0c;文章作者不为此承担任何责任。 简介 D-LinkNAS是由D-Link公司制造的网络附加存储设备。…

01 SQL基础 -- 初识数据库与安装

一、初识数据库 数据库是将大量数据保存起来,通过计算机加工而成的可以进行高效访问的数据集合。该数据集合称为数据库(Database, DB)。用来管理数据库的计算机系统称为数据库管理系统(Database Management System, DBMS) 1.1 DBMS 的种类 DBMS 主要通过数据的保存格式…