P. Bultinck

Polynomial scaling approximations and Dynamic Correlation Corrections to Doubly Occupied Configuration Interaction wave functions

M. Van Raemdonck, D. Alcoba, W. Poelmans, S. De Baerdemacker, A. Torre, L. Lain, G. Massaccesi, D. Van Neck, P. Bultinck
Journal of Chemical Physics
2015
A1

Abstract 

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 analyzed, by studying the overlap between DOCI and full configuration interaction (FCI) 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, a simulated annealing (SA) method is introduced as a global energy based optimisation procedure of the single-particle orbitals. Next, we
expand on the value of seniority number minimising orbitals from truncated configuration interaction (CI) wave functions for approximating the optimal single-particle orbitals of DOCI and related wave functions that contain predominantly seniority zero determinants. It is shown that the resulting seniority number and energy minimising orbitals are especially suitable for the description of static correlation.

Variational optimization of the second order density matrix corresponding to a seniority-zero configuration interaction wave function

W. Poelmans, M. Van Raemdonck, B. Verstichel, S. De Baerdemacker, A. Torre, L. Lain, G. Massaccesi, D. Alcoba, P. Bultinck, D. Van Neck
Journal of Chemical Theory and Computation (JCTC)
2015
A1

Abstract 

We perform a direct variational determination of the second-order (two-particle) density matrix corresponding to a many-electron system, under a restricted set of the two-index $N$-representability $\mathcal{P}$-, $\mathcal{Q}$-, and $\mathcal{G}$-conditions. In addition, we impose a set of necessary constraints that the two-particle density matrix must be derivable from a doubly-occupied many-electron wave function, i.e.\ a singlet wave function for which the Slater determinant decomposition only contains determinants in which spatial orbitals are doubly occupied. We rederive the two-index $N$-representability conditions first found by Weinhold and Wilson and apply them to various benchmark systems (linear hydrogen chains, He, $\text{N}_2$ and $\text{CN}^-$). This work is motivated by the fact that a doubly-occupied many-electron wave function captures in many cases the bulk of the static correlation. Compared to the general case, the structure of doubly-occupied two-particle density matrices causes the associate semidefinite program to have a very favorable scaling as $L^3$, where $L$ is the number of spatial orbitals. Since the doubly-occupied Hilbert space depends on the choice of the orbitals, variational calculation steps of the two-particle density matrix are interspersed with orbital-optimization steps (based on Jacobi rotations in the space of the spatial orbitals). We also point to the importance of symmetry breaking of the orbitals when performing calculations in a doubly-occupied framework.

A hybrid configuration interaction treatment based on seniority number and excitation schemes

D.R. Alcoba, A. Torre, L. Lain, O.B. Ona, P. Capuzzi, M. Van Raemdonck, P. Bultinck, D. Van Neck
Journal of Chemical Physics
141 (24), 244118
2014
A1

Abstract 

We present a configuration interaction method in which the Hamiltonian of an N-electron system is projected on Slater determinants selected according to the seniority-number criterion along with the traditional excitation-based procedure. This proposed method is especially useful to describe systems which exhibit dynamic (weak) correlation at determined geometric arrangements (where the excitation-based procedure is more suitable) but show static (strong) correlation at other arrangements (where the seniority-number technique is preferred). The hybrid method amends the shortcomings of both individual determinant selection procedures, yielding correct shapes of potential energy curves with results closer to those provided by the full configuration interaction method. (c) 2014 AIP Publishing LLC.

Open Access version available at UGent repository

Non-Variational Orbital Optimization Techniques for the AP1roG Wave Function

K. Boguslawski, P. Tecmer, P.W. Ayers, P. Bultinck, S. De Baerdemacker, D. Van Neck
Journal of Chemical Theory and Computation (JCTC)
10 (11), 4873-4882
2014
A1

Abstract 

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.

Efficient description of strongly correlated electrons with mean-field cost

K. Boguslawski, P. Tecmer, P.W. Ayers, P. Bultinck, S. De Baerdemacker, D. Van Neck
Physical Review B
89 (20), 201106
2014
A1

Abstract 

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-50 hydrogen chain.

Open Access version available at UGent repository

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

K. Boguslawski, P. Tecmer, P.A. Limacher, P.A. Johnson, P.W. Ayers, P. Bultinck, S. De Baerdemacker, D. Van Neck
Journal of Chemical Physics
140 (21), 214114
2014
A1

Abstract 

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 C-2 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. (C) 2014 AIP Publishing LLC.

Aliovalent doping of CeO2: DFT study of oxidation state and vacancy effects

D.E.P. Vanpoucke, P. Bultinck, S. Cottenier, V. Van Speybroeck, I. Van Driessche
Journal of Materials Chemistry A
2 (3), 13723-13737
2014
A1

Abstract 

The modification of CeO2 properties by means of aliovalent doping is investigated within the ab initio density functional theory framework. Lattice parameters, dopant atomic radii, bulk moduli and thermal expansion coefficients of fluorite type Ce1-xMxO2-y (with M = Mg, V, Co, Cu, Zn, Nb, Ba, La, Sm, Gd, Yb, and Bi) are presented for 0.00 < x < 0.25. The relative stability of the doped systems is discussed, and the influence of oxygen vacancies is investigated. It is shown that oxygen vacancies tend to increase the lattice parameter, and strongly decrease the bulk modulus. Defect formation energies are correlated with calculated crystal radii and covalent radii of the dopants, and are shown to present no simple trend. The previously observed inverse relationship between the thermal expansion coefficient and the bulk modulus in group IV doped CeO2 [J. Am. Ceram. Soc., 2014, 97(1), 258] is shown to persist independent of the inclusion of charge compensating vacancies.

Open Access version available at UGent repository

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

P.A. Limacher, T.D. Kim, P.W. Ayers, P.A. Johnson, S. De Baerdemacker, D. Van Neck, P. Bultinck
Molecular Physics
112 (5-6), 853-862
2014
A1

Abstract 

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.

Open Access version available at UGent repository

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

P.A. Limacher, P.W. Ayers, P.A. Johnson, S. De Baerdemacker, D. Van Neck, P. Bultinck
Physical Chemistry Chemical Physics (PCCP)
16 (11), 5061-5065
2014
A1

Abstract 

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 Surjan 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 N-2 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.

Open Access version available at UGent repository

Pages

Subscribe to RSS - P. Bultinck