Order-N and embedded-cluster first-principles DFT calculations using SIESTA/Mosaico

Abstract

Order-N and embedded-cluster first-principles DFT calculations have been performed with the Mosaico method for energy optimization (Seijo and Barandiaran in J Chem Phys 121:6698, 2004) for the first time. The Hamiltonian matrix elements have been computed with the SIESTA code. The order-N behavior of the method in DFT calculations was shown in total energy calculations performed on bulk silicon using supercells up to Si8000. The sizes of the orbital-specific-basis-sets needed for precise calculations have been explored in demanding (bulk silicon) and favorable (water clusters) cases for a method based on the calculation of localized molecular orbitals. Embedded-cluster calculations, which are much faster than full-system calculations, have been performed on an Si-vacancy of bulk silicon and on a water cluster with a displacing water molecule. The feasiability of calculations of this type with Mosaico has been demonstrated. The sizes of the variationally free, active clusters which are needed for an agreement with full-system calculations have been explored and result to be reasonably small.