I am studying electronic structural features of 2D materials which contain almost more than 30-40 atoms in the unit cell. I think it is more tedious writing direct python scripts like for a 2 atom unit cell material. Can anyone please suggest me to succeed in this case?