I'm trying to optimize a large number of structures( above thousand) with 60 atoms in each. I need a cheap method for optimization. I came across DFTB , which can be helpful for my calculation. Its there is Gaussian and Cp2K but I don't know if there is any module or keyword specific for such calculations. Can anyone please help me with the inputs ?