A key component in calculations of exchange and correlation energies is the Coulomb operator, which requires the evaluation of two-electron integrals. For localized basis sets, these four-center integrals are most efficiently evaluated with the resolution of identity (RI) technique, which expands basisfunction products in an auxiliary basis. In this work we show the practical applicability of a localized RI-variant ('RI-LVL'), which expands products of basis functions only in the subset of those auxiliary basis functions which are located at the same atoms as the basis functions. We demonstrate the accuracy of RI-LVL for Hartree-Fock calculations, for the PBE0 hybrid density functional, as well as for RPA and MP2 perturbation theory. Molecular test sets used include the S22 set of weakly interacting molecules, the G3 test set, as well as the G2-1 and BH76 test sets, and heavy elements including titanium dioxide, copper and gold clusters. Our RI-LVL implementation paves the way for linear-scaling RI-based hybrid functional calculations for large systems and for all-electron manybody perturbation theory with significantly reduced computational and memory cost.

We present a series of capping-potentials designed as link atoms to saturate dangling bonds at the quantum/classical interface within density functional theory-based hybrid QM/MM calculations. We aim at imitating the properties of different carbon-carbon bonds by means of monovalent analytic pseudopotentials. These effective potentials are optimized such that the perturbations of the quantum electronic density are minimized. This optimization is based on a stochastic scheme, which helps to avoid local minima trapping. For a series of common biomolecular groups, we find capping-potentials that outperform the more common hydrogen-capping in view of structural and spectroscopic properties. To demonstrate the transferability to complex systems, we also benchmark our potentials with a hydrogen-bonded dimer, yielding systematic improvements in structural and spectroscopic parameters.

We present a first principles approach to compute the response of the molecular electronic charge distribution to a geometric distortion. The scheme is based on an explicit representation of the linear electronic susceptibility. The linear electronic susceptibility is a tensor quantity which directly links the first-order electronic response density to the perturbation potential, without requiring self-consistency. We first show that the electronic susceptibility is almost invariant to small changes in the molecular geometry. We then compute the dipole moments from the response density induced by the geometrical changes. We verify the accuracy by comparing the results to the corresponding values obtained from the self-consistent calculations of the ground-state densities in both geometries.