各位版主、大神:

你们好。我想问几个关于LAMMPS建模的基本问题。我知道可能有点简单,但我真的是认真查看过Manual中相关的命令说明之后,仍有一些疑问。我自己想了很久,也问过很多人,但是还是没有解决。

命令来自一个很基础的LAMMPS例子,该例子的全部命令如下:

# Find minimum energy fcc configuration

# Mark Tschopp, 2010

# ---------- Initialize Simulation ---------------------

clear

units metal

dimension 3

boundary p p p

atom_style atomic

atom_modify map array

# ---------- Create Atoms ---------------------

lattice  fcc 4

region box block 0 1 0 1 0 1 units lattice

create_box 1 box

lattice fcc 4 orient x 1 0 0 orient y 0 1 0 orient z 0 0 1

create_atoms 1 box

replicate 1 1 1

# ---------- Define Interatomic Potential ---------------------

pair_style eam/alloy

pair_coeff * * Al99.eam.alloy Al

neighbor 2.0 bin

neigh_modify delay 10 check yes

# ---------- Define Settings ---------------------

compute eng all pe/atom

compute eatoms all reduce sum c_eng

# ---------- Run Minimization ---------------------

reset_timestep 0

fix 1 all box/relax iso 0.0 vmax 0.001

thermo 10

thermo_style custom step pe lx ly lz press pxx pyy pzz c_eatoms

min_style cg

minimize 1e-25 1e-25 5000 10000

variable natoms equal "count(all)"

variable teng equal "c_eatoms"

variable length equal "lx"

variable ecoh equal "v_teng/v_natoms"

print "Total energy (eV) = ${teng};"

print "Number of atoms = ${natoms};"

print "Lattice constant (Angstoms) = ${length};"

print "Cohesive energy (eV) = ${ecoh};"

print "All done!"

我的问题是:

(1)lattice fcc 4 和 lattice fcc 4 orient x 100 orient y 0 1 0 orient z 0 0 1这两句有什么区别?效果不是一样的么?为什么在lattice fcc 4建立格子之后,还要有lattice fcc 4 orient x 100 orient y 0 1 0 orient z 0 0 1这句?意义何在?

(2)replicate 1 1 1这句不相当于没有复制么?为什么要写这条命令呢?意义何在?

(3)如果我只想建模,不进行任何计算,然后想在LAMMPS中导出该模型,最后在VMD中查看该模型。那么在整个代码的最后一行需要添加一个怎样的dump语句?(我查过手册中的dump语句,也试过很多次自己写,但是都不对。)

跪求各路大神帮助。小弟在此谢过!

lammps建模_【求助】关于LAMMPS建模 - 分子模拟 - 小木虫 - 学术 科研 互动社区相关推荐

  1. vasp计算压电系数_求助DFTP算出来的压电系数 - 第一原理 - 小木虫 - 学术 科研 互动社区...

    我用的VASP5.2 下的LEPSILON-tag算的压电系数和BEC等东西,之前发过类似的贴,不过这次想一起贴上来在好好求助一下: 我的INCAR SYSTEM=PLPK ENCUT=400 PRE ...

  2. arcgis 同名图层合并_【求助】ArcGIS中怎样合并图层 - 地学 - 小木虫 - 学术 科研 互动社区...

    这个问题 最简单的方式是用editor中的切割多边形功能  可以把一个矢量文件按照你定义的线分割  可直接在属性表中算出面积  要是想直接用一个线文件去切割面文件的话 操作有点复杂  下面是用edit ...

  3. 差分电荷密度怎么画_【求助】如何做电荷密度图和差分电荷密度图? - 量子化学 - 小木虫 - 学术 科研 互动社区...

    人们从化学实践出发,对分子的电子结构一直非常关注.多年来,人们总是想探讨分子中原子的性质,特别是分子中原子的净电荷,提出了各种计算方法,并进行了不断地改进,然而至今都还不能令人满意.特别是原子净电荷划 ...

  4. matlab 摄氏度符号怎么打,【求助】弱弱的问下,Matlab如何在公式中输入希腊字母? - 程序语言 - 小木虫 - 学术 科研 互动社区...

    很多时候都要在matlab画图的时候添加一些公式符号之类的,有一些特殊的字符并不能直接从键盘上输入,比如希腊字母等等.但是有想用,因为这样使图看起来漂亮而且容易理解. 例如:我想输入摄氏度的符号,怎么 ...

  5. matlab decomposition filters,MATLAB小波去噪求助(附算法和显示图片)!不知自己哪个地方出了问题,求指点! - 信息科学 - 小木虫 - 学术 科研 互动社区...

    star2883 不是非常明白楼主的意思.我的理解是: 问题1:为什么中间两幅图,显示的都是空白? 答:你是指Figure2和Figure3的显示中,全是白色吗?这个是因为你用了colormap(ma ...

  6. lisp 计算三点的夹角_已知三点坐标求夹角 - 程序语言 - 小木虫 - 学术 科研 互动社区...

    最近需要处理大量数据,需要对化学键键角批量输出. 我已经把原子坐标以xyz的形式给出. 现在搞不定如何把夹角以degree(度数)的方式给求出来. 哪位知道怎么弄? 本人很菜,别笑话. #It is ...

  7. matlab 水文频率曲线,【求助】如何计算水文频率,外行求教 - 地学 - 小木虫 - 学术 科研 互动社区...

    根据某水文现象的统计特性,利用现有水文资料,分析水文变量设计值与出现频率(或重现期)之间的定量关系的工作过程称为水文频率计算. 自然界的现象按发生情况可分成:必然事件,即在一定条件下必然会发生的事情, ...

  8. matlab自家自减,【求助】VC环境下自增(++)自减(--) - 信息科学 - 小木虫 - 学术 科研 互动社区...

    VC++环境下连续自增(自减)运算规则 1.        在一般表达式中和在函数参数中计算的结果相同(这一点与Turbo C环境不同): 2.        后缀:同表达式中Turbo C的情况相同 ...

  9. qpsk的matlab平方根升余,【求助】求助:【QPSK功率谱密度的matlab的代码】 - 仿真模拟 - 小木虫 - 学术 科研 互动社区...

    CODE: clear; %随机产生长度大于1000的'0'.'1'信号序列,对其进行QPSK调制 %定义待仿真序列的维数 N global N N=2000; %定义产生'1'的概率为 p glob ...

最新文章

  1. VI连接esx主机时报错“client.xml”错误
  2. mysql 环形复制_mysql复制(Replication)
  3. plot与legend画图与图例
  4. 第二届「星斗奖」获奖名单正式公布
  5. 利用计算机或图形计算器在,图形计算器在函数教学中的应用
  6. 如何用50行代码构建情感分类器
  7. 作者:秦兵(1968-),女,哈尔滨工业大学计算机科学与技术学院教授、博士生导师。...
  8. ThreadPoolExecutor解析
  9. python选取列中数据的范围_Python Pandas中根据列的值选取多行数据
  10. python re span_Python的re模块与正则表达式小结
  11. linux centos安装wine qq,ubuntu安装wine QQ
  12. 一步步学习微软InfoPath2010和SP2010--第五章节--添加逻辑和规则到表单(1)--InfoPath中初级类型的表单逻辑
  13. HTML5 标签audio添加网页背景音乐代码
  14. Python基础知识详解(十)小结,用python实现教师信息管理系统
  15. MATLAB符号计算相关函数
  16. 003_java三大版本
  17. Java base64转inputStream
  18. 数据库管理员、系统分析员、数据库设计人员,应用程序员的职责是什么?
  19. 三天晚上看了24集 央视版《神雕侠侣》 还不错
  20. 大数据下机械智能诊断的机遇与挑战,阅读文献系列(一)

热门文章

  1. 王齐 。《后科技时代的思考》
  2. 宠物成长记录小程序开发,开辟宠物市场新流量
  3. 多级分销系统数据库模型(猜想)
  4. ORACLE---函数
  5. InnoDB索引数据结构
  6. es查询所欲,检索ElasticSearch中所有_id的高效方法
  7. NX10实例教程——草图初识
  8. 各种增删改查(员工管理)
  9. Reversed Numbers
  10. C语言标准 —— C89(C90)、C99、C11、C17、C2X