How to add randomly distributed vacancies or inter.. in the crystal??

Thanks a lot?
I want to add 2% vacancies in a bulk randomly.
How to get to through MS modelling?