原子扰动python脚本

原子扰动python脚本

from dpdata import LabeledSystem,MultiSystems
from glob import glob
import dpdata
perturbed_system = dpdata.System('./POSCAR').perturb(pert_num=3, 
    cell_pert_fraction=0.05, 
    atom_pert_distance=0.6, 
    atom_pert_style='normal')
for i in range(3):
    perturbed_system[i].to_poscar('%s.vasp'%i)