Search

This work proposes the variational determination of two-electron reduced density matrices corresponding to the ground state of N-electron systems within the doubly occupied-configuration-interaction methodology. The P, Q, and G two-index N-representability conditions have been extended to the T1 and T2 (T2') three-index ones and the resulting optimization problem has been addressed using a standard semidefinite program. We report results obtained from the doubly occupied-configuration-interaction method, from the two-index constraint variational procedure and from the two- and three-index constraint variational treatment. The discussion of these results along with a study of the computational cost demanded shows the usefulness of our proposal.

A class of polynomial scaling methods that approximate Doubly Occupied Configuration Interaction (DOCI) wave functions and improve the description of dynamic correlation is introduced. The accuracy of the resulting wave functions is analysed by comparing energies and studying the overlap between the newly developed methods and full configuration interaction wave functions, showing that a low energy does not necessarily entail a good approximation of the exact wave function. Due to the dependence of DOCI wave functions on the single-particle basis chosen, several orbital optimisation algorithms are introduced. An energy-based algorithm using the simulated annealing method is used as a benchmark. As a computationally more affordable alternative, a seniority number minimising algorithm is developed and compared to the energy based one revealing that the seniority minimising orbital set performs well. Given a well-chosen orbital basis, it is shown that the newly developed DOCI based wave functions are especially suitable for the computationally efficient description of static correlation and to lesser extent dynamic correlation.

We introduce new nonvariational orbital optimization schemes for the antisymmetric product of one-reference orbital geminal (AP1roG) wave function (also known as pair-coupled cluster doubles) that are extensions to our recently proposed projected seniority-two (PS2-AP1roG) orbital optimization method [ J. Chem. Phys. 2014 , 140 , 214114 )]. These approaches represent less stringent approximations to the PS2-AP1roG ansatz and prove to be more robust approximations to the variational orbital optimization scheme than PS2-AP1roG. The performance of the proposed orbital optimization techniques is illustrated for a number of well-known multireference problems: the insertion of Be into H2, the automerization process of cyclobutadiene, the stability of the monocyclic form of pyridyne, and the aromatic stability of benzene.

The geometric and electronic structure of the MIL-47(V) metal-organic framework (MOF) is investigated by using ab initio density functional theory (DFT) calculations. Special focus is placed on the relation between the spin configuration and the properties of the MOF. The ground state is found to be antiferromagnetic, with an equilibrium volume of 1554.70 Å(3). The transition pressure of the pressure-induced large-pore-to-narrow-pore phase transition is calculated to be 82 MPa and 124 MPa for systems with ferromagnetic and antiferromagnetic chains, respectively. For a mixed system, the transition pressure is found to be a weighted average of the ferromagnetic and antiferromagnetic transition pressures. Mapping DFT energies onto a simple-spin Hamiltonian shows both the intra- and inter-chain coupling to be antiferromagnetic, with the latter coupling constant being two orders of magnitude smaller than the former, suggesting the MIL-47(V) to present quasi-1D behavior. The electronic structure of the different spin configurations is investigated and it shows that the band gap position varies strongly with the spin configuration. The valence and conduction bands show a clear V d-character. In addition, these bands are flat in directions orthogonal to VO6 chains, while showing dispersion along the the direction of the VO6 chains, similar as for other quasi-1D materials.

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

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.