【LAMMPS学习】八、基础知识(2.2)类型标签

8. 基础知识

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

8.1.通用基础知识

8.2. 设置入门

8.2.1.二维模拟

8.2.2.类型标签

2022 年 9 月 15 日版本中的新增功能。

LAMMPS 中的每个原子都有一个关联的数字原子类型。类似地,每个 bond, angle, dihedral, 和 improper都被指定为键类型、角类型等。

这些类型的主要用途是将势(力场)参数映射到原子、 bond, angle, dihedral, 和 improper的相互作用。

默认情况下,类型值以 1 到 Ntypes 之间的整数形式输入,无论它们出现在 LAMMPS 输入或输出文件中。每次交互的 Ntype 总数在创建模拟框时被“锁定”。

LAMMPS 最近添加的一项功能是使用字符串(称为类型标签)作为替代方案。使用类型标签代替数字类型在各种场景中都具有优势。例如,类型标签可以使输入更具可读性和通用性(即,可以通过为类型分配不同数值的不同系统使用 include 命令。这种通用性也适用于其他输入,例如由 read_data 读取的数据文件或由 read_data 读取的分子模板文件分子命令。请参阅下面的列表,了解可以以不同方式使用类型标签的其他命令。

LAMMPS 将在内部继续使用数字类型,这意味着许多以前的限制仍然适用。例如,在创建模拟框时,类型的总数被锁定,并且即使没有被任何交互使用,也必须提供每种类型的潜在参数。

所有类型种类(原子类型、键类型等)的类型标签集合存储为“标签映射”,它只是数字类型及其关联类型标签的列表。在类型种类中,每个类型标签必须是唯一的。它只能分配给一种数字类型。要在给定类型种类的数据文件中读取和写入类型标签,所有关联的数值类型都需要分配一个类型标签。可以使用 labelmap write 命令保存部分映射,并使用 include 命令读回。

有效的类型标签可以包含大多数 ASCII 字符,但不能以数字、“#”或“*”开头。此外,标签不得包含空格字符。在 LAMMPS 输入中使用 labelmap 命令时,如果类型标签中出现某些字符,例如单引号 (') 或双引号 (") 或 '#' 字符,则标签必须放入双引号、单引号、或三引号(“””)。三引号允许使用最通用的类​​型标签字符串,但它们需要有前导和尾随空格。定义类型标签时,空格将被忽略。例子:

labelmap angle 1 """ C1'-C2"-C3# """

此命令会将字符串 `C1'-C2"-C3#` 映射到角度类型 1。

有两种定义标签映射的方法。一种是通过 labelmap 命令。另一种是通过 read_data 命令。数据文件可以包含原子类型标签、键类型标签等部分,这些部分将类型标签分配给数字类型。可以通过  write_data 命令将标签映射写入数据文件。该映射还可以通过 write_restart 和 read_restart 命令写入重新启动文件或从重新启动文件中读取。

在 LAMMPS 输入或输出中使用类型标签 

许多采用数字类型作为参数的 LAMMPS 输入脚本命令可以改用关联的类型标签。如果没有为特定数字类型定义类型标签,则只能使用其数字类型。

此示例为原子类型分配标签,然后使用类型标签重新定义对系数。

pair_coeff 1 2 1.0 1.0              # numeric types
labelmap atom 1 C 2 H
pair_coeff C H 1.0 1.0              # type labels

向各种命令添加对类型标签的支持是一个正在进行的项目。如果输入脚本命令(或命令读取的文件中的部分)允许用类型标签替换数字类型参数,则该命令的文档页面中将明确提及。

作为临时措施,输入脚本命令可以利用变量以及如何在输入处理过程中扩展变量。变量可以使用将类型标签字符串转换为当前标签映射中定义的各自数字的函数。详细信息请参见变量命令。

例如,如果尚不支持类型标签,则pair_coeff命令可以与类型标签一起使用,无论是使用显式变量命令还是pair_coeff命令中使用的隐式变量。

labelmap atom 1 C 2 H
variable atom1 equal label2type(atom,C)
variable atom2 equal label2type(atom,H)
pair_coeff ${atom1} ${atom2} 1.0 1.0
labelmap atom 1 C 2 H
pair_coeff $(label2type(atom,C)) $(label2type(atom,H)) 80.0 1.2

可以使用标签类型的命令 

任何涉及读取多个数据文件、分子模板或两者组合的工作流程都可以通过使用类型标签而不是数字类型来简化,因为类型会在文件之间自动同步。使用类型标签时,为fix bond/react 创建模拟就绪反应模板要简单得多,并且生成的模板可以在多个模拟或不同系统中无需修改即可使用。

fb0790c10a094df09aad28dac816ac7a.png

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

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

相关文章

arhtas idea plugin 使用手册

arthas idea plugin 使用文档 语雀

GridView 多次更新数据后,导致notifyDataSetChanged,只更新了UI,而数据不对,点击某个Item时位置不对。

最近在做作项目的时候使用到了GridView,需要显示的数据分为多类,例如:全部、工作、学习、社交和电商等等,当用户点击某个一级分类时,需要把这个分类下的所有数据重新设置到GridView,我每次通过adpter自定义…

基于ssm的校园短期闲置资源置换平台(java项目+文档+源码)

风定落花生,歌声逐流水,大家好我是风歌,混迹在java圈的辛苦码农。今天要和大家聊的是一款基于ssm的校园短期闲置资源置换平台。项目源码以及部署相关请联系风歌,文末附上联系信息 。 项目简介: 校园短期闲置资源置换…

大话设计模式——9.单例模式(Singleton Pattern)

简介 确保一个类只有一个实例,并提供全局访问点来获取该实例,是最简单的设计模式。 UML图: 单例模式共有两种创建方式: 饿汉式(线程安全) 提前创建实例,好处在于该实例全局唯一,不…

C程序的编译

经过预处理后的源文件,退去一切包装,注释被删除,各种预处理命令也基本上被处理掉,剩下的就是原汁原味的C代码了。接下来的第二步,就开始进入编译阶段。编译阶段主要分两步:第一步,编译器调用一系列解析工具,去分析这些C代码,将C源文件编译为汇编文件;第二步,通过汇编…

anaconda创建了虚拟python环境,且安装了pytorch,但是pycharm中import torch运行时报错

报错如下: C:\Users\tashi\.conda\envs\test1\python.exe D:\project\python\transformer\main.py C:\Users\tashi\.conda\envs\test1\lib\site-packages\numpy\__init__.py:127: UserWarning: mkl-service package failed to import, therefore Intel(R) MKL init…

电子邮件协议学习

电子邮件协议学习 电子邮件服务器的核心功能是发送、接收和存储电子邮件;依赖于一系列协议,SMTP,IAMP,POP3,SMTPS ,IAMPS等协议。 SMTP(Simple Mail Transfer Protocol) 原理 SMTP基于文本的…

AI预测体彩排3第2弹【2024年4月13日预测--第1套算法开始计算第2次测试】

各位小伙伴,今天实在抱歉,周末回了趟老家,回来比较晚了,数据今天上午跑完后就回老家了,晚上8点多才回来,赶紧把预测结果发出来吧,虽然有点晚了,但是咱们前面说过了,目前的…

将Visio绘图导出PDF文件,使其自适应大小,并去掉导入Latex的边框显示

问题描述 将Visio绘图导成pdf文件,首先在Visio绘图如下: 如果直接导出或者另存为pdf文件,则会发现pdf文件是整个页面大小,而不是图片大小。而且在导入latex等排版工具现实时,会显示边框。 问题解决 1.调整Visio中的页…

mongodb中的多表查询aggregate中排序不是按全表排序,而是当前页排序问题如何解决?

在上篇文章中讲到了多表查询https://blog.csdn.net/guige8888811/article/details/133880362 其中还少一个知识点则是排序。而起初写了排序条件之后发现排序并不是按全表排的。代码如下: await mongodb.getConnection(attendanceRecord).aggregate([{$lookup: {from: userInf…

AI-漫画推文

🍨🍨项目源码 私信 🍨🍨技术选型 前端:Vue Vite Electron TS后端:Spring Boot MybatisPlus Redis Mysql

vox2vec论文速读

vox2vec: A Framework for Self-supervised Contrastive Learning of Voxel-Level Representations in Medical Images 摘要 本文介绍了 vox2vec——一种体素级表示的自监督学习 (SSL) 对比方法 vox2vec 表示由特征金字塔网络 (FPN&#xf…

Cascader 级联选择器 - 选择器最后一级数据为空

原因:将扁平数据转化为树形数据时,给每个项都添加了 children export const transList2Tree (list, rootPid) > {const result []list.forEach(item > {if (item.pid rootPid) {const children transList2Tree(list, item.id)item.children …

深圳风控建模岗薪资水平如何?

不管是学生还是工作的小伙伴,估计都对不同岗位工作几年的薪酬水平比较感兴趣。之前的文章提供爬取招聘网站,获取某类工作招聘信息的实现逻辑和代码。感兴趣可翻看批量爬取招聘网站上工作岗位的招聘信息。本文对爬取下来的信息进行清洗处理,以深圳风控建模岗为例,分析不同工…

Linux C++ 034-STL之谓词

Linux C 034-STL之谓词 本节关键字:Linux、C、谓词 相关库函数: 谓词概念 概念: 返回值为bool类型的仿函数称为谓词 如果operator()接受一个参数,那么叫做一元谓词 如果operator()接受两个参数,那么叫做二元谓词 一…

Unity绘制地图

首先在项目/Assets文件夹下创建一个Tiles文件夹 在层级下点击鼠标右键选择2D对象选择瓦片地图创建Tilemap。 选择地图素材 如果素材需要裁剪,在检查器Sprite模式选择多个,点击Sprite Editor,选择切 ,选择类型Grid By Cell Count,…

BoostCompass(建立正排索引和倒排索引模块)

阅读导航 一、模块概述二、编写正排索引和倒排索引模块✅安装 jsoncpp✅Jieba分词库的安装1. 代码基本框架2. 正排索引的建立3. 倒排索引的建立 三、整体代码⭕index.hpp 一、模块概述 这个模块我们定义了一个名为Index的C类,用于构建和维护一个文档索引系统。该系…

3D视觉技术

1. 简介 3D视觉是一个多学科相融合的技术,可以总结为:计算图形学计算机视觉人工智能3D视觉。3D视觉技术是通过3D摄像头采集视野空间内每个点位的三维座标信息,通过算法复原获取三维立体成像,不会轻易受到外界环境、复杂光线的影响…

Fiddler的安装和使用

Fiddler是一款强大的网络调试工具,可以帮助开发者进行网络请求和响应的调试和分析。以下是关于Fiddler的安装和使用的简要说明: 安装: 访问Fiddler的官方网站,下载对应版本的安装程序(一般选择经典版下载 - Fiddler Classic&…

【计算机毕业设计】停车场管理系统——后附源码

🎉**欢迎来到琛哥的技术世界!**🎉 📘 博主小档案: 琛哥,一名来自世界500强的资深程序猿,毕业于国内知名985高校。 🔧 技术专长: 琛哥在深度学习任务中展现出卓越的能力&a…