When I delete the surface atom of the 2D material to search for the vacancy energy, the atom will be lost. and my code are as follows: units metal dimension 3 boundary p p p atom_style atomic read_data data.origin change_box all boundary p p f remap units box pair_style sw pair_coeff * * wse2.sw W Se #计算势能 compute eng all pe/atom compute eatoms all…