build a serials models of different point defect concentration and then optimize the models by nvt and npt. As the point defect concentration, the system will become unstable. And you can detect that concentration by analysing the optimized structure and energy.