| 查看: 3264 | 回复: 9 | ||
[求助]
【求助】关于LAMMPS建模
|
|
各位版主、大神： 你们好。我想问几个关于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语句，也试过很多次自己写，但是都不对。） 跪求各路大神帮助。小弟在此谢过！ |
回复此楼» 猜你喜欢
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
-
材料工程专硕274一志愿211求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助lammps晶界建模如何保持box的周期性~~
已经有11人回复
- 求助LAMMPS关于MSD计算的问题
已经有17人回复
- 请求高手关于ATK/VNL中分子建模的问题
已经有11人回复
- MS建模问题求助
已经有13人回复
- 用MS建模导入LAMMPS计算，是不是不能用compass力场？
已经有9人回复
- lammps求助：计算含有不同原子体系的热导率的in文件该怎么修改
已经有7人回复
- 水分子建模方法
已经有6人回复
- 关于用Lammps模拟Cu熔点的例子问题
已经有7人回复
- lammps里面的pair_coeff设置问题
已经有3人回复
- lammps的时间步长问题
已经有16人回复
- 请教关于materials studio建模和OPLS力场
已经有7人回复
- 在LAMMPS中怎么建立图示的模型
已经有6人回复
- 【讨论】想学lammps是否应该先学习C++
已经有28人回复
- 【求助】fedora中LAMMPS编译问题
已经有4人回复
- 【讨论】编译lammps时Makefile的修改问题
已经有8人回复
- 【求助】lammps如何测量模拟以后材料的晶格参数？
已经有7人回复
- 【讨论】lammps中关于fix ave/spatial命令的输出问题
已经有4人回复
- 【求助】lammps中多元化合物的晶格创建问题
已经有14人回复
- 【求助】将MS数据文件转化成lammps data
已经有25人回复
- 【求助】MS建模导入lammps问题求救
已经有14人回复
- 【求助】Lammps 多机并行的问题
已经有11人回复
2楼2013-10-14 19:15:00
天道and酬勤
铁杆木虫 (著名写手)
- 应助: 18 (小学生)
- 金币: 14371.8
- 散金: 71
- 红花: 11
- 帖子: 1469
- 在线: 111.3小时
- 虫号: 1097985
- 注册: 2010-09-13
- 专业: 力学
【答案】应助回帖
感谢参与，应助指数 +1
|
对于问题1：好好看看手册，说的蛮详细。The other 3 arguments are the crystallographic direction in the lattice that you want to orient along that axis, specified as integers. E.g. "orient x 2 1 0" means the x-axis in the simulation box will be the [210] lattice direction. 两条命令有差别； 对于问题2：我觉得可以去掉； 对于问题3：dump 1 all atom 1 dump.lammpstrj 就可以了，其中的1不可以和其他的dump重复；1000为每多少步输出一次，可以自己设定，看手册就知道啥意思了。关键在于文件名后缀为lammpstrj。 PS：我没仔细看你的程序，只是针对你的问题。你可以试下！ |
3楼2013-10-16 17:31:04
4楼2013-10-27 09:20:31
Ergon
金虫 (著名写手)
- 应助: 9 (幼儿园)
- 金币: 1332.5
- 散金: 220
- 红花: 5
- 沙发: 2
- 帖子: 1553
- 在线: 105.5小时
- 虫号: 3295777
- 注册: 2014-06-27
- 性别: GG
- 专业: 理论和计算化学
5楼2015-05-23 16:05:18
FunWithFan
禁虫 (初入文坛)
|
本帖内容被屏蔽 |
6楼2015-11-20 22:07:24
7楼2016-01-04 09:29:47
8楼2016-01-10 22:18:28
ruibingqilin
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 682.7
- 帖子: 285
- 在线: 152.6小时
- 虫号: 1943683
- 注册: 2012-08-17
- 性别: GG
- 专业: 高分子材料结构与性能
9楼2019-02-10 15:26:06
raolixiang
新虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 509.5
- 散金: 1540
- 红花: 6
- 帖子: 739
- 在线: 115.8小时
- 虫号: 3721641
- 注册: 2015-03-09
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
10楼2019-09-15 16:05:05
20