大家好,我是小马老师。
如果lammps模拟体系中包含分子结构,最常见的流程是使用MS、packmol等软件建模,然后转换为lammps data文件,在lammps读入data文件。
其实,lammps也可以进行分子建模,对于一些比较简单的分子模型,使用lammps建模也比较方便。
本文以水球模型为例,介绍lammps分子类结构建模方法。

lammps创建分子的原理是把分子当做一个原子添加到体系中,和添加原子的命令一样,也是create_atoms。
命令语法规则为:
create_atoms mol template-ID seed
mol:添加分子的关键词
template-ID:分子模板,等同于ms中分子单体
seed:随机数种子,分子取向随机
具体的步骤为:
(1)创建分子模板(template-ID)文件

分子模板描述了一个分子单体内原子的类型、原子坐标、键、角等信息。
如下图为水分子的模板文件water.data部分内容,请加QQ群:754749935下载完整的water.data文件。

(2)编写in文件

units real
atom_style full
#创建box
region box block -20 20 -20 20 -20 20 units box
#box内原子数量、键、角等信息
create_box 2 box bond/types 1 angle/types 1 extra/bond/per/atom 2 &extra/special/per/atom 2 extra/angle/per/atom 1
#H2O的晶格,可调节晶格尺寸控制水分子数量或水的密度
lattice         bcc 3.92
#设定区域
region          H2O sphere 0 0 0 18 units box
#设定分子模板id为water,对应文件water.data
molecule        water water.data
#在球形区域内填充水分子
create_atoms    0 region H2O mol water 3453 units box
#设置原子摩尔质量
mass                1 16
mass                2 1
#保存模型文件
write_data          H2O.data


公众号案例代码以及学习交流,请加QQ群:754749935

扫描关注微信公众号:lammps加油站。

lammps自带命令create_atoms实现水分子建模相关推荐

  1. [lammps案例] region命令实现六面体核壳结构建模

    核壳结构在先进复合材料的发展中具有特殊的意义,因为它们可以有效地将不同的组分在纳米尺度上聚集在一起,这种结构的优点很大程度上依赖于核心和外壳的关键设计.     下图为本文建立的六面体核壳结构模型. ...

  2. 【转载】dos下通过wmic命令查看硬盘和内存/CPU信息(windows自带命令查看硬件信息)

    更新时间:2014年05月11日 00:26:16   投稿:mdxy-dxy 这篇文章主要介绍了windows系统自带命令查看硬件信息,怎样dos命令查看硬盘和内存/CPU信息的方法,需要的朋友可以 ...

  3. python编写请求参数带文件_转载:如何编写一个带命令行参数的Python文件

    看到别人执行一个支持命令行参数的python文件,瞬间觉得高大上起来.牛逼起来,那么如何编写一个带命令行参数的python脚本呢?不用紧张,下面将简单易懂地让你学会如何让自己的python脚本,支持命 ...

  4. c语言学习进阶-C语言带命令行参数的文件数据批量计算

    C语言带命令行参数的文件数据批量计算 (1)假设输入文件中的数据为多行整数,把数据文件命名为"InputData.txt",也可以命名为任意的文件名,尝试设计一个c 语言程序, 主 ...

  5. c语音学习进阶-C语言带命令行参数的程序设计

    C语言带命令行参数的程序设计 (1)深入理解并学习带命令行参数的程序设计方法,使用c 语言完成带命令行参数的科学计算程序,能从给定的任意文件名的文件当中输入数据,并能向给定的任意文件名的输出文件当中输 ...

  6. 转载《利用Windows系统自带命令手工搞定病毒》_原水_新浪博客

    转载 360软件百科的文章,其实WINDOWS有很多的命令我们不熟悉,你认真看下就会发现其他的操作系统其实就是这样做的,只是我们没有想到WINDOWS也可以这样做. 利用Windows系统自带命令手工 ...

  7. vs2013中带命令行参数的调试方法

    @VS2013带命令行参数的调试方法 在C++函数参数传递中,我们知道,main函数是演示C++程序如何向函数传递数组的好例子.一般来说,我们定义的大多main函数都只有空的形参列表: int mai ...

  8. 用C语言计算水分子的个数

    题目: 一个水分子的质量约为,1夸脱水大约有950g,编写一个程序,要求输入水的夸脱数,然后显示这么多水中含有多少水分子. 问题分析: 从题目中可以得到1夸脱水中水分子的含量 = 1夸脱水 * 950 ...

  9. 一个水分子的质量是 3.0 * 10^(-6) g,一升水的质量是 950g输入水的升数(整数),输出这些水中包含多少个水分子

    利用宏定义实现 注:科学计数法:3.0e-23 %e #include <stdio.h> #define water 3.0e-23 #define quality 950 int ma ...

  10. (win和win server)dos下通过wmic命令查看硬盘和内存/CPU信息(windows自带命令查看硬件信息)

    如何在windows系统自带命令查看硬件信息,怎样dos命令查看硬盘和内存/CPU信息?最直接的是:开始→运行→CMD打开命令提示符,在该窗口下输入systeminfo执行,即可看到几乎所有想知道的系 ...

最新文章

  1. 华华和月月种树(牛客)
  2. java将一个数字转换为数组_Java 数组
  3. 关系型数据库是如何运作的
  4. python机器学习---1. K近邻算法
  5. Merge Two Sorted Lists leetcode java
  6. XP访问Linux共享错误提示
  7. Wythoff's game
  8. 74AHC1G32GW 74系列逻辑芯片
  9. pyCharm最新2018.2激活 附激活文件/教育版免费使用
  10. 方便面又好卖了!康师傅天猫618狂增10倍,还说这只是起点
  11. 分布式的坑(队列、缓存、分库分表、事务)
  12. java网络爬虫如何控制爬取的深度_Java 动手写爬虫: 二、 深度爬取
  13. mat1 and mat2 shapes cannot be multiplied (3584x7 and 512x4096)
  14. 一键快速设置图层lisp程序_AE脚本-多图层一键快速排列 Staircase V1.03 + 使用教程...
  15. Java单元测试实践-09.Mockito的Stub参数条件
  16. 皕杰报表之小程序代码质量检测
  17. JAVA小项目之影院系统和员工管理系统
  18. ESP8266EX 串口WIFI无线模块
  19. 手变小是怎么回事_手小该怎么办呢?如何变大???
  20. 如何强制卸载软件,强制卸载的工具。

热门文章

  1. AI-语音处理理论和应用-语音识别、语音合成
  2. 计算机用户删除文件找回,电脑删除的文件怎么找回?失易得数据恢复帮助找回...
  3. HR告诉你,怎么回答“为什么从上家公司离职?”
  4. 基于Helm和Operator的K8S应用管理
  5. Ubuntu系统安装和运行n2n
  6. vue-video-player 断点续播
  7. gcd函数(C/C++)
  8. goodFeaturesToTrack——Shi-Tomasi角点检测
  9. springboot毕设项目网上商城系统8l44b(java+VUE+Mybatis+Maven+Mysql)
  10. win10--笔记本能连上手机热点,但不能连接WIFI