Starting from a $\ce{BaTiO3}$ structure, I want to dope this with Strontium to obtain $\ce{Ba_{0.67}Sr_{0.33}TiO3}$. Is there a way to do this using LAMMPS?
I began a discussion about this in Chinese on the Computational Chemistry Commune. I know this can be done using a scripting language (e.g. Bash/Python) to edit my input file, but I was curious if LAMMPS or some other existing tool can already automate this process.