熔化和凝固氩,铜,铝.PDFVIP

Project #5 熔化与凝固：氩，铜，铝 铜和铝的熔化转变: 对于铜和铝，LAMMPS 建立 8×8×5 的FCC 晶格体系;充分弛豫后利用 Nose-Hover 方法，保 持压强为零，使体系从 T=2.5K 开始加热，直至发生熔化转变。 下面是铜熔化的输入文件： # LAMMPS Melt_Cu or Al units metal # 单位，指定为 lammps 里的金属类的单位，长度为 Å ，能量为 eV 。 boundary p p p # 周期性边界条件 atom_style atomic # 原子模式 variable x equal 2.5 # 定义变量 x 为初始温度 lattice fcc 3.61 # Cu 的晶格常数3.61 #lattice fcc 4.05 region box block 0 8 0 8 0 5 # x,y,z 各方向上的晶胞重复单元数，也即区域大小 create_box 1 box # 将上述区域指定为模拟的盒子 create_atoms 1 box # 将原子按晶格填满盒子 timestep 0.01 # 步长 0.005fs thermo 1000 #每隔1000 步输出热力学结果 pair_style eam/alloy # 选取 Cu 的EAM 势作为模型 pair_coeff * * jin_copper_lammps.setfl Cu # EAM 势文件名称 #pair_style eam/fs #pair_coeff * * Al_FM.eam.fs Al neighbor 0.5 bin neigh_modify every 5 delay 0 check yes #velocity all create $x 825577 dist gaussian fix 1 all nvt $x $x 1.0 drag 0.2 #保持初始温度，在 NVT 下弛豫 #compute 3 all pe/atom #compute 4 all ke/atom #compute 5 all coord/atom 3.0 #dump 1 all custom 1 dump.atom id xs ys zs c_3 c_4 c_5 run 10000 # 运行10000 步 unfix 1 fix 1 all npt $x 2000 4.00 xyz 0.0 0.0 6.0 drag 0.2 #在 NPT 下加热至 2000K #fix 1 all npt $x 1500 4.00 xyz 0.0 0.0 6.0 drag 0.2 run 1200000 #运行 1200000 步 对于铜，其熔点为 1357.77K，但在我们的模拟中其在 1609K 附近发生一级相变，比其平 衡时熔点增大了 18.7%。 对于铝，其熔点为 934.477K, 在我们的模拟中，其在 1110.0K 附近发生一级相变，单位体 积发生突变。其相对于平衡熔点增大了了 18.9%。 从铜和铝的熔化过热，以及后面所涉及的氩的熔化过热与凝固过冷，我们可以看出，在利 用分子动力学方法模拟熔化与凝固时往往会发生过热与过冷，其值基本在 10%-30%之间。 产生过热与过冷的因素基本上可以从热力学与动力学的方面阐述，例如均匀形核而导致。 铜在特定温度下的性质: 利用 LAMMPS 建立 8×8×8 的FCC 格子，分别在 10K，500K，1000K， 1800K 和 2000