Dear Developers, I think it's better if we can use more than one basis set for one element such as: C SZV C_surface DZVP It can save a lot of time, especially for very large system. I'm trying to modiffy some parts of the code but it might be not sufficient. Best regards, Nguyen, Ngoc Ha Hanoi National University of Education