Display options
Share it on

J Chem Theory Comput. 2012 Oct 09;8(10):3565-73. doi: 10.1021/ct200897x. Epub 2012 Mar 14.

Linear Scaling Self-Consistent Field Calculations with Millions of Atoms in the Condensed Phase.

Journal of chemical theory and computation

Joost VandeVondele, Urban Borštnik, Jürg Hutter

Affiliations

  1. Department of Materials, ETH Zurich , Wolfgang-Pauli-Strasse 27, 8093 Zurich, Switzerland.
  2. Physical Chemistry Institute, University of Zurich , Winterthurerstrasse 190, CH-8057 Zurich, Switzerland.

PMID: 26593003 DOI: 10.1021/ct200897x

Abstract

In this work, the applicability and performance of a linear scaling algorithm is investigated for three-dimensional condensed phase systems. A simple but robust approach based on the matrix sign function is employed together with a thresholding matrix multiplication that does not require a prescribed sparsity pattern. Semiempirical methods and density functional theory have been tested. We demonstrate that self-consistent calculations with 1 million atoms are feasible for simple systems. With this approach, the computational cost of the calculation depends strongly on basis set quality. In the current implementation, high quality calculations for dense systems are limited to a few hundred thousand atoms. We report on the sparsities of the involved matrices as obtained at convergence and for intermediate iterations. We investigate how determining the chemical potential impacts the computational cost for very large systems.

Publication Types