LAMMPS计算Cu的空位形成能.pdf

LAMMPS计算Cu的空位形成能 金剑锋 材料学院 东北大学 Cu的空位形成能 N-1 个原子 N 个原子 + E1 Ec E E = E – [(N-1)/N]×E vac 1 E = (E +E )–N×E vac 1 c c 采用LAMMPS 内置的命令 建立脚本 units metal atom_style atomic lattice fcc 3.615 region box block 0 10 0 10 0 10 create_box 1 box E1 create_atoms 1 box region vacancy sphere 5.0 5.0 5.0 0.1 side in #晶格单位 delete_atoms region vacancy … fix 1 all box/relax aniso 0.0 run 100 minimize 1.0e-8 1.0e-10 50000 100000 Cu的空位形成能 + E1 Ec E E = -14156.06 eV E = -3.5402 eV E = -3.5402×4000 eV 1 c 3999个原子 1个原子 4000个原子 E = (E +E ) – E = 1.20 eV vac 1 0 = (-14156.06 - 3.5402) - (-3.5402×4000) Matlab 建立空位晶体结构 输入到 LAMMPS脚本 for i = 1 : 1 : nt if (abs(data(1:3,i)-[0, 0, 0]) 0.1) data(1:3,i) = data(1:3,nt); % 删除坐标(0,0,0)原子建立空位 nt=nt-1; i=i-1; end end ... fileID = fopen(‘Cu_vacancy.lam’, ‘w’); % 晶体结构文件名 fprintf(fileID, ‘%s\n’,‘Cu’); % Line 1 ：元素名称 fprintf(fileID, ‘\n’); % Line 2 ：空行 fprintf(fileID, ‘%9d%s\n’, num, ‘ atoms’); % Line 3 ：原子总数 fprintf(fileID, ‘\n’); % Line 4 ：空行 fprintf(fi