CP2K is an open source electronic structure and molecular dynamics software package to perform atomistic simulations of solid-state, liquid, molecular, and biological systems. It is especially aimed at massively parallel and linear-scaling electronic structure methods and state-of-the-art ab initio molecular dynamics simulations. Excellent performance for electronic structure calculations is achieved using novel algorithms implemented for modern high-performance computing systems. This review revisits the main capabilities of CP2K to perform efficient and accurate electronic structure simulations. The emphasis is put on density functional theory and multiple post–Hartree–Fock methods using the Gaussian and plane wave approach and its augmented all-electron extension.

Common density-matrix functionals, the Müller and the power functional, have been benchmarked for the half-filled Hubbard dimer, which allows to model the bond dissociation problem and the transition from the weakly to the strongly correlated limit. Unbiased numerical calculations are combined with analytical results. Despite the well known successes of the Müller functional, the ground state is degenerate with a one-dimensional manifold of ferromagnetic solutions. The resulting infinite magnetic susceptibility indicates another qualitative flaw of the Müller functional. The derivative discontinuity with respect to particle number is not present indicating an incorrect metallike behavior. The power functional actually favors the ferromagnetic state for weak interaction. Analogous to the Hartree-Fock approximation, the power functional undergoes a transition beyond a critical interaction strength, in this case however, to a non-collinear antiferromagnetic state.

Abstract. In materials with strong electron correlation the proper treatment of local atomic physics described by orbital occupations is crucial. Reduced density-matrix functional theory is a natural extension of density functional theory for systems that are dominated by orbital physics. We review the current state of reduced density-matrix functional theory (RDMFT). For atomic structure relaxations or abinitio molecular dynamics the combination of density functional theory (DFT) and dynamical mean-field theory (DMFT) possesses a number of disadvantages, like the cumbersome evaluation of forces. We therefore describe a method, DFT+RDMFT, that combines many-particle effects based on reduced density-matrix functional theory with a density functional-like framework. A recent development is the construction of density-matrix functionals directly from many-particle theory such as methods from quantum chemistry or many-particle Green's functions. We present the underlying exact theorems and describe current progress towards quantitative functionals.

No abstract

A method, called the adaptive cluster approximation (ACA), for single-impurity Anderson models is proposed. It is based on reduced density-matrix functional theory, where the one-particle reduced density matrix is used as the basic variable. The adaptive cluster approximation introduces a unitary transformation of the bath states such that the effect of the bath is concentrated to a small cluster around the impurity. For this small effective system one can then either calculate the reduced density-matrix functional numerically exact from Levy's constrained-search formalism or approximate it by an implicit approximation of the reduced density-matrix functional. The method is evaluated for single-impurity Anderson models with finite baths. The method converges rapidly to the exact result with the size of the effective bath.