掌握自动批量计算平衡晶格常数:
- 平衡晶格常数(equilibrium lattice constant)对应的体系能量是最低的,因此只需要计算一些列不同晶格常数下体系的能量,那么体系能量最小时对应的晶格常数就是平衡晶格常数。
代码学习:
例如:
- 直接上in文件
- 提取log.lammps中的数据
学习内容:
#General
units metal
boundary p p p
atom_style atomic
# build the model
lattice fcc $x
region box block 0 10 0 10 0 10
create_box 3 box #设置原子种类
create_atoms 1 region box #向模拟区域添加原子
set type 1 type/fraction 2 0.40 873934
set type 1 type/fraction 3 0.20 123452
# specify the potential
pair_style eam/alloy
pair_coeff * * NiCoCr.lammps.eam Ni Co Cr
mass 1 58.69 # Ni
mass 2 58.9332 # Co
mass 3 51.996 # Cr
#variable v equal vol
timestep 0.005
thermo 10
# minimize the total energy
min_style cg
minimize 1.0e-12 1.0e-12 1000 1000
loop
jump in.eam lata
$ grep ^@ log.lammps > lat.vs.Ecoh.step1
提取log文件数据
#每次新窗口均需要先运行这两个moudle
module load fftw/3.3.10-openmpi411-gnu940
module load intel/parastudio/2017.8.262
#提取log。lammps中的数据
grep ^@ log.lammps > lat.vs.Ecoh.step1
学习总结:
要做自己舒适区之外的事,这些事才能带来最大的杠杆效益。
尽可能少做哪些假高回馈的重复性行动。
尽可能的做创新工作。