| 查看: 2527 | 回复: 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" |
回复此楼» 猜你喜欢
268求调剂
已经有9人回复
-
301求调剂
已经有4人回复
-
0703化学336分求调剂
已经有5人回复
-
301求调剂
已经有4人回复
-
材料专硕306英一数二
已经有9人回复
-
304求调剂
已经有8人回复
-
材料工程专硕调剂
已经有5人回复
-
344求调剂
已经有4人回复
-
326求调剂
已经有6人回复
-
296求调剂
已经有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人回复
» 抢金币啦！回帖就可以得到:
-
坐标上海，诚征女友，非常 着急，私信必回
+1/471
-
物理学 调剂
+1/85
-
物理学 调剂
+1/84
-
物理学 调剂
+1/80
-
福建师范大学招收2026年化学、材料硕士3-4名
+1/78
-
何时使用 CODA™ 科里奥利质量流量仪表- 艾里卡特（Alicat）
+2/52
-
湖北工业大学 生命科学与健康学院 课题组招收2026级食品/生物方向硕士
+1/42
-
华南师范大学（211）申请审核制博士招生- 电子、自动化、机械、生物学、物理相关专业
+1/25
-
国家级领军人才团队高分子材料方向2026年博士研究生招生
+1/16
-
【2026博士招生/博后招聘】北京航空航天大学潘彪课题组——AI芯片设计方向
+1/11
-
重庆大学诚招2026年生物材料方向博士生
+1/8
-
哈尔滨理工大学材料与化学工程学院 纳米功能材料与光电催化团队 招收硕士研究生
+1/7
-
英国埃克塞特大学 & 法国巴黎萨克雷大学联合培养博士
+1/7
-
东北林业大学（211，双一流）宋老师课题组招收高分子、化学、复合材料等方向研究生
+1/6
-
层流压差式流量测控产品在辐射防护及环保分析领域的应用与优势
+1/3
-
加拿大阿尔伯塔大学招收电磁、无线通信、机器学习方向全奖硕士/博士/博士后/访问学者
+2/2
-
招收理论凝聚态物理/纳米光学/量子计算方向博士、硕士研究生/博士后
+1/1
-
085404 270求调剂，b区求导师收留接受跨专业
+1/1
-
北京工业大学材料工程调剂03
+1/1
-
【2026年春季补录】厦门大学田中群院士/易骏教授纳米光子学团队招收2026年入学博士生
+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
回复
顶一下，
