Project #4
表面与界面能
铜的表面能
当物体形成表面时，表面上的原子键发生断裂，接近表面的几层原子不再如之前处于平衡状
态，从而导致能量的升高，升高的温度便是物体的表面能。
利用 LAMMPS 做出 20*20*40 fcc 的盒子，删去边缘的原子制造出一段真空层；算出此时
体系的总能量
0
E
，
然后从中间把盒子切成两半并移至足够远的距离，此时的体系总能量为
E
f
inal
，
从而表面能：
0
2
final
surface
EE
A
γ
=
−
A
为表面的面积
(100) 面与 (111) 面
如下是输入文件 in.surface_Cu_100
# LAMMPS Cu _Surface_100
units metal
boundary p p p
atom_style atomic
lattice fcc 3.61
region box block 0 20 0 20 0 40
create_box 1 box
create_atoms 1 box
timestep 0.005
thermo 5
pair_style eam/alloy
pair_coeff * * jin_copper_lammps.setfl Cu
region boundary1 block INF INF INF INF 29.9 INF
region boundary2 block INF INF INF INF INF 9.9
group boundary1 region boundary1
group boundary2 region boundary2
group boundary union boundary1 boundary2
评论0