| 查看: 2803 | 回复: 4 | |||
[交流]
【讨论】LAMMPS 晶界计算 输入脚本
|
|||
|
在一个网站上看到的一个研究晶界的能量随各个自由度变化的LAMMPS输入文件https://ccg.hpc.msstate.edu/medi ... boundary_generation，很强大的样子，可是有一些设置不是太明白，拿出来大家讨论一下。 我有这么几处不明白，为什么overlapinic 要设为86？是怎么得到86这个数的？ overlapdist 为什么要设为0.275+0.005*(v_c-1)？我知道它应该是测试晶界区域中不同的原子数带来的影响，可是这个0.275是怎么得到的？ if $d == 2 then "delete_atoms overlap ${overlapdist} upper lower" 这句看起来比较奇怪，好像d只能是1，为什么还有这个句子？ 望高手予以解答!! # LAMMPS Input File for Grain Boundaries # Mark Tschopp, Dec2009 # This file will generate numerous input files for LAMMPS # using a large number of grain boundaries # ---------- Setup Variables --------------------- variable etol equal 1.0e-25 variable ftol equal 1.0e-25 variable maxiter equal 5000 variable maxeval equal 10000 variable latparam equal 2.855312 variable minimumenergy equal -4.122435 variable overlapboth equal 1 variable gbname index Fe_100STGB1 variable counter equal 0 variable inc equal "v_latparam / 6" # Insert x,y,z sizes in LU and calculate in Angstroms variable xsize1 equal "sqrt(0^2 + 2^2 + 1^2)" variable zsize1 equal "sqrt(1^2 + 0^2 + 0^2)" variable xsize2 equal "sqrt(0^2 + 2^2 + -1^2)" variable zsize2 equal "sqrt(1^2 + 0^2 + 0^2)" if ${xsize1} <= ${xsize2} then "variable xsize equal ${xsize1}" else "variable xsize equal ${xsize2}" if ${zsize1} <= ${zsize2} then "variable zsize equal ${zsize1}" else "variable zsize equal ${zsize2}" variable xlen equal "v_xsize * v_latparam" variable zlen equal "v_zsize * v_latparam" # Determine number of increments for displacement grid in the in-plane GB directions variable xinc equal "floor(v_xlen / v_inc)" variable zinc equal "floor(v_zlen / v_inc)" # Implement overlap criterion variable overlapinc equal 86 # ---------- Define loops for simulation --------------------- label loopa variable a loop ${xinc} variable tx equal "(v_a-1) / v_xinc * v_xsize" label loopb variable b loop ${zinc} variable tz equal "(v_b-1) / v_zinc * v_zsize" label loopd variable d loop ${overlapboth} label loopc variable c loop ${overlapinc} variable overlapdist equal "0.275 + 0.005 * (v_c-1)" # ---------- Calculate counter and create data directory --------------------- variable ctemp equal ${counter}+1 variable counter equal ${ctemp} variable ctemp delete print "Counter: ${counter}" shell mkdir ${gbname} # ---------- Initialize Simulation --------------------- clear units metal dimension 3 boundary p p p atom_style atomic # ---------- Create Atomistic Structure --------------------- lattice bcc ${latparam} region whole block 0.000000 6.384672 -121.308763 121.308763 0.000000 2.855312 units box create_box 2 whole region upper block INF INF 0.000000 121.308763 INF INF units box lattice bcc ${latparam} orient x 0 2 1 orient y 0 -1 2 orient z 1 0 0 create_atoms 1 region upper region lower block INF INF -121.308763 0.000000 INF INF units box lattice bcc ${latparam} orient x 0 2 -1 orient y 0 1 2 orient z 1 0 0 create_atoms 2 region lower group upper type 1 group lower type 2 # ---------- Define Interatomic Potential --------------------- pair_style eam/fs pair_coeff * * /cavs/cmd/data1/users/mtschopp/LAMMPS/lammps-12Nov09/potentials/Fe_2.eam.fs Fe Fe neighbor 2.0 bin neigh_modify delay 10 check yes # ---------- Displace atoms and delete overlapping atoms --------------------- displace_atoms upper move ${tx} 0 ${tz} units lattice if $d == 1 then "delete_atoms overlap ${overlapdist} lower upper" if $d == 2 then "delete_atoms overlap ${overlapdist} upper lower" if $c == 1 then "variable atomprev equal 1" variable natoms equal "count(all)" print "Previous: ${atomprev}, Present: ${natoms}" if ${atomprev} == ${natoms} then "jump GB_Fe_100STGB1.in loopend" # ---------- Define Settings --------------------- compute csym all centro/atom compute eng all pe/atom compute eatoms all reduce sum c_eng # ---------- Run Minimization --------------------- reset_timestep 0 thermo 10 thermo_style custom step pe lx ly lz press pxx pyy pzz c_eatoms min_style cg minimize ${etol} ${ftol} ${maxiter} ${maxeval} # ---------- Run Minimization 2--------------------- # Now allow the box to expand/contract perpendicular to the grain boundary reset_timestep 0 thermo 10 thermo_style custom step pe lx ly lz press pxx pyy pzz c_eatoms fix 1 all box/relax aniso NULL 0.0 NULL vmax 0.001 min_style cg minimize ${etol} ${ftol} ${maxiter} ${maxeval} # ---------- Calculate GB Energy --------------------- variable esum equal "v_minimumenergy * count(all)" variable xseng equal "c_eatoms - (v_minimumenergy * count(all))" variable gbarea equal "lx * lz * 2" variable gbe equal "(c_eatoms - (v_minimumenergy * count(all)))/v_gbarea" variable gbemJm2 equal ${gbe}*16021.7733 variable gbernd equal round(${gbemJm2}) print "After third minimization:" print "GB energy is ${gbemJm2} mJ/m^2" # Store number of atoms for overlap criterion, i.e., do not rerun equivalent configurations variable atomprev equal "v_natoms" # ---------- Dump data into Data file ------------- shell cd Fe_100STGB1 reset_timestep 0 timestep 0.001 velocity all create 20 95812384 fix 2 all npt 1 1 100 xyz 0 0 100 drag 0.2 dump 1 all custom 1000 dump.${counter}_${gbernd} id type x y z c_csym c_eng run 0 shell cd .. # ---------- End of loop structure ------------- label loopend next c jump GB_Fe_100STGB1.in loopc variable c delete next d jump GB_Fe_100STGB1.in loopd variable d delete next b jump GB_Fe_100STGB1.in loopb variable b delete next a jump GB_Fe_100STGB1.in loopa print "All done" |
回复此楼» 猜你喜欢
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意！
已经有13人回复
-
同一个分子，一条来自文献，一条来自AI——不告诉你答案，你会选哪条？
已经有7人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 各位高手，怎么在windows下编译LAMMPS？
已经有4人回复
- lammps的时间步长问题
已经有16人回复
- 寻求lammps的应用 欢迎大家畅所欲言
已经有8人回复
- 【求助】LAMMPS-运行报错
已经有9人回复
- 【讨论】想学lammps是否应该先学习C++
已经有28人回复
- 【讨论】lammps里的boundary的位置能够随意更改？
已经有14人回复
- 【求助】将MS数据文件转化成lammps data
已经有25人回复
- 【求助】MS建模导入lammps问题求救
已经有14人回复
- 【求助】Lammps 多机并行的问题
已经有11人回复
» 抢金币啦！回帖就可以得到:
-
2026年材料结构与力学性能国际会议（MSMP 2026
+1/796
-
持续创业者，A8离异男征婚，事业遇到瓶颈，寻找生活伴侣
+1/479
-
专利代理人
+1/96
-
复旦大学生物与信息方向博士研究生等招聘
+1/95
-
急招哈尔滨工业大学（深圳）机械学科2027级"申请-考核"博士研究生1名
+1/87
-
招聘岗位：生物研究员（睡眠/脂代谢/抗衰老方向） 重庆 2人
+1/84
-
南科大薛亚辉课题组诚聘离子输运、二维材料、原子力显微镜等方向“快响行动”博士生
+1/66
-
南方科技大学力航系急招收纳米材料、力学、机械等方向推免硕士研究生两名
+1/62
-
华南师范大学胡勇军教授课题组招收2027年推荐面试攻读研究生（推免生）
+1/31
-
澳大利亚昆士兰大学博士生招聘/气体发酵、酵母代谢工程与合成生物学方向
+2/28
-
哈工大深圳（机器人方向）-姜欣课题组招收2027年博士研究生（10月7日截止）
+1/15
-
【西北工业大学-985西安】材料学院-余泓课题组-诚招2027年博士
+1/12
-
申请博后
+1/7
-
【2027硕士招生】招收材料、化工、环境、冶金类硕士研究生
+1/4
-
浙江大学Kemal Celebi教授课题组招收2027年秋季博士生-材料/化学/化工
+1/1
-
中科院东莞材料所｜新能源材料研究部招聘研究员 / 副研究员（可中央事业编制）
+1/1
-
ICTCPS（纺织、着色与高分子科学）Textiles, Coloration and Polymer Science
+1/1
-
2027 International Conference on Healthcare, Medicine and Medical Science
+1/1
-
2027 IFPAM（基础物理与先进材料）
+1/1
-
2027 CoASS（天文学与空间科学）
+1/1
2楼2011-03-29 18:34:29
3楼2012-06-05 21:54:28
5楼2015-05-11 10:54:00
简单回复
lcl1cxt4楼
2015-01-16 20:52
回复
顶一下，
