Peter A. Limacher

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 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.

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

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.

Assessing the Accuracy of New Geminal-Based Approaches

We present a systematic theoretical study on the dissociation of diatomic molecules and their spectroscopic constants using our recently presented geminal-based wave function ansätze. Specifically, the performance of the antisymmetric product of rank two geminals (APr2G), the antisymmetric product of 1-reference-orbital geminals (AP1roG) and its orbital-optimized variant (OO-AP1roG) are assessed against standard quantum chemistry methods. Our study indicates that these new geminal-based approaches provide a cheap, robust, and accurate alternative for the description of bond-breaking processes in closed-shell systems requiring only mean-field-like computational cost. In particular, the spectroscopic constants obtained from OO-AP1roG are in very good agreement with reference theoretical and experimental data.

Cross-conjugation: Cross-Conjugation

Cross‐conjugation, even though very common in chemistry, has only sparsely been visited by experiment or by theory. Until very recently, this form of conjugation, in which two separate conjugated branches are linked to the same carbon atom which is part of a (shared) double bond, was considered to be much less able to promote electron delocalization than through‐conjugation. Therefore, exploiting cross‐conjugation to design new materials was considered difficult and not very promising. In this article, we show that this view, at least for neutral molecules in their electronic ground state, is essentially correct. For cross‐conjugated radical ions, the situation is different: There is appreciable charge delocalization across the bifurcation point, i.e., the point where the two conjugated branches meet. Simple Hückel molecular orbital considerations show that the connectivity pattern encountered in cross‐conjugation will lead to enhanced electron delocalization effects. This observation, confirmed by density functional calculations, also applies for electronically excited species. Therefore, cross‐conjugation may be exploited to build molecular switches or may be used to design devices such as molecular transistors.

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

We present a new, non-variational orbital-optimization scheme for the antisymmetric product of one-reference orbital geminal wave function. Our approach is motivated by the observation that an orbital-optimized seniority-zero configuration interaction (CI) expansion yields similar results to an orbital-optimized seniority-zero-plus-two CI expansion [L. Bytautas, T. M. Henderson, C. A. Jimenez-Hoyos, J. K. Ellis, and G. E. Scuseria, J. Chem. Phys. 135, 044119 (2011)]. A numerical analysis is performed for the C2 and LiF molecules, for the CH2 singlet diradical as well as for the symmetric stretching of hypothetical (linear) hydrogen chains. For these test cases, the proposed orbital-optimization protocol yields similar results to its variational orbital optimization counterpart, but prevents symmetry-breaking of molecular orbitals in most cases.

A systematic analysis of the structure and (hyper)polarizability of donor-acceptor substituted polyacetylenes using a Coulomb-attenuating density functional

In this paper we perform a systematic investigation on all-trans polyacetylene chains of different lengths, end-capped with moieties of different donor or acceptor natures and different strengths, to infer useful structure/property relationship rules and behavioral patterns. The values for bond length alternation (BLA), longitudinal polarizability, and first and second hyperpolarizabilities have been computed with the Coulomb-attenuating density function (CAM-B3LYP), using response theory. A comparison of the relative effect that each end-capping combination contributes to BLA, linear, and nonlinear optical coefficients has been performed. This results in useful insights and general rules to ad hoc tailoring the molecular response for a specific characteristic.

On the effect of electron correlation on the static second hyperpolarizability of π conjugated oligomer chains

In this article, we report on the ab initio calculation of the static longitudinal second hyperpolarizability (γ) of π conjugated unsaturated oligomer chains using polyacetylene and polyyne as model compounds. The common observation is that the electron correlation enhances γ in these systems. The present study reveals that for extended chain lengths the opposite appears to be true: Electron correlation may have a damping effect on this property. For double-zeta basis sets, a negative contribution from electron correlation to γ is found within the range of chain lengths investigated. For triple-zeta basis sets, the same behavior must be anticipated at larger chain lengths based on extrapolation schemes. The analysis of the excitation energies and transition moments shows that transition moments between excited states as predicted by the Hartree-Fock and coupled cluster methods have a different response to chain length extension. There also are indications that higher order correlation effects will enhance γ.

Simple and inexpensive perturbative correction schemes for antisymmetric products of nonorthogonal geminals

A new multireference perturbation approach has been developed for the recently proposed AP1roG scheme, a computationally facile parametrization of an antisymmetric product of nonorthogonal geminals. This perturbation theory of second-order closely follows the biorthogonal treatment from multiconfiguration perturbation theory as introduced by Surján et al., but makes use of the additional feature of AP1roG that the expansion coefficients within the space of closed-shell determinants are essentially correct already, which further increases the predictive power of the method. Building upon the ability of AP1roG to model static correlation, the perturbation correction accounts for dynamical electron correlation, leading to absolute energies close to full configuration interaction results. Potential surfaces for multiple bond dissociation in H2O and N2 are predicted with high accuracy up to bond breaking. The computational cost of the method is the same as that of conventional single-reference MP2.

A scheme for the evaluation of electron delocalization and conjugation efficiency in linearly π-conjugated systems

In this study, we present a scheme for the evaluation of electron delocalization and conjugation efficiency in lineraly π-conjugated systems. The scheme, based on the natural bond orbital theory, allows monitoring the evolution of electron delocalization along an extended conjugation path as well as its response to chemical modification. The scheme presented is evaluated and illustrated by means of a computational investigation of π-conjugation in all-trans polyacetylene [PA; H(-CH═CH)n-H], polydiacetylene [PDA, H(-C≡C-CH═CH)n-H], and polytriacetylene [PTA, H(-C≡C-CH═CH-C≡C)n-H] with up to 180 carbon atoms, all related by the number of ethynyl units incorporated in the chain. We are able to show that for short oligomers the incorporation of ethynyl spacers into the PA chain increases the π-delocalization energy, but, on the other hand, reduces the efficiency with which π-electron delocalization is promoted along the backbone. This explains the generally shorter effective conjugation lengths observed for the properties of the polyeneynes (PDA and PTA) relative to the polyenes (PA). It will also be shown that the reduced conjugation efficiency, within the NBO-based model presented in this work, can be related to the orbital interaction pattern along the π-conjugated chain. We will show that the orbital interaction energy pattern is characteristic for the type and the length of the backbone and may therefore serve as a descriptor for linearly π-conjugated chains.