| 查看: 1618 | 回复: 4 | ||
| 【悬赏金币】回答本帖问题，作者木棉和我将赠送您 100 个金币 | ||
木棉和我木虫 (职业作家)
|
[求助]
求助：建立的初始简单晶体体系(NaCl结构)，npt弛豫时即发生体系“爆炸” 已有1人参与
|
|
|
如题，建立了一个NbC完美晶体，尺寸与取向分别为：10a*10a*10a; X-[100], Y-[010], and Z-[001]。 准备简单计算一下相关势函数下的晶格常数，即在5K下进行npt保温弛豫，模型开始计算则发生“爆炸”。具体Lammps运行脚本如下： units metal boundary p p p atom_style atomic read_data nbc_c10-10-10.lmp #perfect NbC crystal (NaCl structure)；由Atomsk建立. pair_style meam/c pair_coeff * * library.meam Nb Fe C NbFeC.meam Nb C timestep 0.001 thermo 100 thermo_style custom step temp pe ke etotal epair press vol lx ly lz restart 100000 fcc-Fe.restart neighbor 0.5 bin neigh_modify every 2 delay 10 check yes dump 1 all custom 100000 dump.lattice-nbc-5K.* id type x y z velocity all create 5 123456 dist gaussian fix 1 all npt temp 5 5 0.1 iso 0 0 0.1 run 500000 unfix 1 undump 1 相关势函数、in文件及初始晶体请见附件。 另外，MD构建初始晶体时，是否需要考虑晶体结构的化学键类型？还是只是需要把各类原子的初始位置确定，而化学键类型主要影响作用力，作用力通过势函数体现？或者说，建立MD结构模型时，是不是仅有位置信息，没有其他信息，其他信息由势函数决定。 请相关老师们帮忙解答一下，十分感谢。 |
回复此楼
valenhou001
至尊木虫 (职业作家)
- 应助: 241 (大学生)
- 金币: 25701.7
- 散金: 602
- 红花: 166
- 帖子: 3782
- 在线: 873.8小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
2楼2021-12-11 06:34:54
3楼2021-12-12 00:26:14
木棉和我
木虫 (职业作家)
- 应助: 5 (幼儿园)
- 金币: 9175.1
- 散金: 2362
- 红花: 3
- 帖子: 3337
- 在线: 681.3小时
- 虫号: 1852823
- 注册: 2012-06-08
- 性别: GG
- 专业: 冶金物理化学与冶金原理
4楼2021-12-13 08:34:26
木棉和我
木虫 (职业作家)
- 应助: 5 (幼儿园)
- 金币: 9175.1
- 散金: 2362
- 红花: 3
- 帖子: 3337
- 在线: 681.3小时
- 虫号: 1852823
- 注册: 2012-06-08
- 性别: GG
- 专业: 冶金物理化学与冶金原理
|
谢谢老师您的解答！ 1. 也就说MD晶体构型中，只有位置信息(结果中，可以根据原子位置与速度，提取相关参量)；然后原子间作用力的影响因素，如化学键类型、电荷等，是直接由势函数体现，对吗？ 2. 这个模型是这样的，就是根据文献中的势函数计算出的NbC晶格常数，我建立了一个完美NbC晶体(NaCl结构)。我想确定一下，这个势函数下NbC的晶格常数是不是如文献中所述，于是简单的对完美NbC晶体做了NPT保温，结果体系“爆炸”了。我也尝试过对这个NbC晶体进行cg法的能量最小化，但最小化好像对整个晶体原子位置毫无影响，因为本身建立的NbC晶体就是完美的Nb：C=1:1的NaCl结构。 |
5楼2021-12-13 08:47:32
5