Dimitri Van Neck

A New Mean-Field Method Suitable for Strongly Correlated Electrons: Computationally Facile Antisymmetric Products of Nonorthogonal Geminals.

We propose an approach to the electronic structure problem based on noninteracting electron pairs that has similar computational cost to conventional methods based on noninteracting electrons. In stark contrast to other approaches, the wave function is an antisymmetric product of nonorthogonal geminals, but the geminals are structured so the projected Schrödinger equation can be solved very efficiently. We focus on an approach where, in each geminal, only one of the orbitals in a reference Slater determinant is occupied. The resulting method gives good results for atoms and small molecules. It also performs well for a prototypical example of strongly correlated electronic systems, the hydrogen atom chain.

The density matrix renormalization group for ab initio quantum chemistry

During the past 15 years, the density matrix renormalization group (DMRG) has become increasingly important for ab initio quantum chemistry. Its underlying wavefunction ansatz, the matrix product state (MPS), is a low-rank decomposition of the full configuration interaction tensor. The virtual dimension of the MPS, the rank of the decomposition, controls the size of the corner of the many-body Hilbert space that can be reached with the ansatz. This parameter can be systematically increased until numerical convergence is reached. The MPS ansatz naturally captures exponentially decaying correlation functions. Therefore DMRG works extremely well for noncritical one-dimensional systems. The active orbital spaces in quantum chemistry are however often far from one-dimensional, and relatively large virtual dimensions are required to use DMRG for ab initio quantum chemistry (QC-DMRG). The QC-DMRG algorithm, its computational cost, and its properties are discussed. Two important aspects to reduce the computational cost are given special attention: the orbital choice and ordering, and the exploitation of the symmetry group of the Hamiltonian. With these considerations, the QC-DMRG algorithm allows to find numerically exact solutions in active spaces of up to 40 electrons in 40 orbitals.

Communication: DMRG-SCF study of the singlet, triplet, and quintet states of oxo-Mn(Salen)

We use CheMPS2, our free open-source spin-adapted implementation of the density matrix renormalization group (DMRG) [S. Wouters, W. Poelmans, P. W. Ayers, and D. Van Neck, Comput. Phys. Commun. 185, 1501 (2014)], to study the lowest singlet, triplet, and quintet states of the oxo-Mn(Salen) complex. We describe how an initial approximate DMRG calculation in a large active space around the Fermi level can be used to obtain a good set of starting orbitals for subsequent complete-active-space or DMRG self-consistent field calculations. This procedure mitigates the need for a localization procedure, followed by a manual selection of the active space. Per multiplicity, the same active space of 28 electrons in 22 orbitals (28e, 22o) is obtained with the 6-31G(*), cc-pVDZ, and ANO-RCC-VDZP basis sets (the latter with DKH2 scalar relativistic corrections). Our calculations provide new insight into the electronic structure of the quintet.

Comparison of the Hirshfeld-I and iterated stockholder atoms in molecules schemes

Two recently introduced self-consistent Hirshfeld procedures for obtaining atoms in molecules are compared in detail. The Hirshfeld-I scheme introduces self consistency by requiring that the atomic population of the promolecular atom is equal to that of the atom-in-the-molecule. In the iterated stockholder atoms (ISA) approach, self consistency is obtained by requiring that for every value of the radius of a sphere around every nucleus, the average electron density on the surface of this sphere is the same in the promolecular atom and in the atom in the molecule. The relationships between the two schemes are examined, and common backgrounds and differences are discussed. Whereas it can be argued that the Hirshfeld-I approach has a stronger physical background, the ISA scheme avoids having to define what states of the atoms are to be considered when constructing the promolecule.

The influence of orbital rotation on the energy of closed-shell wavefunctions

The orbital dependence of closed-shell wavefunction energies is investigated by performing doubly-occupied configuration interaction (DOCI) calculations, representing the most general class of these wavefunctions. Different local minima are examined for planar hydrogen clusters containing two, four, and six electrons applying (spin) symmetry-broken restricted, unrestricted, and generalised orbitals with real and complex coefficients. Contrary to Hartree–Fock (HF), restricted DOCI is found to properly break bonds and thus unrestricted orbitals, while providing a quantitative improvement of the energy, are not needed to enforce a qualitatively correct bond dissociation. For the beryllium atom and the BH diatomic, the lowest possible HF energy requests symmetry-broken generalised orbitals, whereas accurate results for DOCI can be obtained within a restricted formalism. Complex orbital coefficients are shown to increase the accuracy of HF and DOCI results in certain cases. The computationally inexpensive AP1roG geminal wavefunction is proven to agree very well with all DOCI results of this study.

Comparative Study of Various Normal Mode Analysis Techniques Based on Partial Hessians

Standard normal mode analysis becomes problematic for complex molecular systems, as a result of both the high computational cost and the excessive amount of information when the full Hessian matrix is used. Several partial Hessian methods have been proposed in the literature, yielding approximate normal modes. These methods aim at reducing the computational load and/or calculating only the relevant normal modes of interest in a specific application. Each method has its own (dis)advantages and application field but guidelines for the most suitable choice are lacking. We have investigated several partial Hessian methods, including the Partial Hessian Vibrational Analysis (PHVA), the Mobile Block Hessian (MBH), and the Vibrational Subsystem Analysis (VSA). In this article, we focus on the benefits and drawbacks of these methods, in terms of the reproduction of localized modes, collective modes, and the performance in partially optimized structures. We find that the PHVA is suitable for describing localized modes, that the MBH not only reproduces localized and global modes but also serves as an analysis tool of the spectrum, and that the VSA is mostly useful for the reproduction of the low frequency spectrum. These guidelines are illustrated with the reproduction of the localized amine‐stretch, the spectrum of quinine and a bis‐cinchona derivative, and the low frequency modes of the LAO binding protein.

Efficient description of strongly correlated electrons with mean-field cost

We present an efficient approach to the electron correlation problem that is well suited for strongly interacting many-body systems, but requires only mean-field-like computational cost. The performance of our approach is illustrated for one-dimensional Hubbard rings with different numbers of sites, and for the nonrelativistic quantum-chemical Hamiltonian exploring the symmetric dissociation of the H

Longitudinal static optical properties of hydrogen chains: finite field extrapolations of matrix product state calculations.

{}_{50}

Subsystem constraints in variational second order density matrix optimization: curing the dissociative behavior.

hydrogen chain.

Projected seniority-two orbital optimization of the antisymmetric product of one- reference orbital geminal

We have implemented the sweep algorithm for the variational optimization of SU(2) U(1) (spin and particle number) invariant matrix product states (MPS) for general spin and particle number invariant fermionic Hamiltonians. This class includes non-relativistic quantum chemical systems within the Born-Oppenheimer approximation. High-accuracy ab initio finite field results of the longitudinal static polarizabilities and second hyperpolarizabilities of one-dimensional hydrogen chains are presented. This allows to assess the performance of other quantum chemical methods. For small basis sets, MPS calculations in the saturation regime of the optical response properties can be performed. These results are extrapolated to the thermodynamic limit.