Source author record

Patrick Bultinck

Patrick Bultinck appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

13works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2023arXiv

Spin-constrained Hartree-Fock and the generator coordinate method for the 2-site Hubbard model

We present a mathematical analysis of the spin-constrained Hartree-Fock solutions (CHF) of the 2-site Hubbard model. The analysis sheds light on the spin symmetry breaking process around the Coulson-Fischer point. CHF states are useful as input states for the Generator Coordinate Method (GCM) in which CHF states can be used as a basis for multiconfigurational calculations in the Hill-Wheeler equations. The spin degree of freedom provides the generator coordinate in the GCM related to static electron correlation.

preprint2022arXiv

Near-exact treatment of seniority-zero ground and excited states with a Richardson-Gaudin mean-field

Eigenvectors of the reduced Bardeen-Cooper-Schrieffer Hamiltonian, Richardson-Gaudin (RG) states, are used as a variational wavefunction Ansatz for strongly-correlated electronic systems. These states are geminal products whose coefficients are solutions of non-linear equations. Previous results showed un-physical behaviour but in this contribution it is shown that with only the variational solution for the ground state, all the seniority-zero states are quite well approximated. The difficulty is in choosing the correct RG state. The systems studied showed a clear choice and we expect it should always be possible to reason physically which state to choose.

preprint2015arXiv

Maximum Probability Domains for Hubbard Models

The theory of Maximum Probability Domains (MPDs) is formulated for the Hubbard model in terms of projection operators and generating functions for both exact eigenstates as well as Slater determinants. A fast MPD analysis procedure is proposed, which is subsequently used to analyse numerical results for the Hubbard model. It is shown that the essential physics behind the considered Hubbard models can be exposed using MPDs. Furthermore, the MPDs appear to be in line with what is expected from Valence Bond Theory-based knowledge.

preprint2014arXiv

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. %which is based on orbital optimization of electron pairs. The performance of our approach is illustrated for the one-dimensional Hubbard model with periodic boundary conditions for different chain lengths, and for the non-relativistic quantum chemical Hamiltonian exploring the symmetric dissociation of the H$_{50}$ hydrogen chain.

preprint2014arXiv

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 [J. Chem. Phys., 135, 044119 (2011)]. A numerical analysis is performed for the C$_2$, LiF and CH$_2$ molecules 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.

preprint2013arXiv

Aliovalent Doping of CeO$_2$: DFT-study of Oxidation State and Vacancy Effects

The modification of the properties of CeO$_2$ through aliovalent doping are investigated within the \emph{ab-initio} density functional theory framework. Lattice parameters, dopant atomic radii, bulk moduli and thermal expansion coefficients of fluorite type Ce$_{1-x}$M$_{x}$O$_{2-y}$ (with M$ = $ Mg, V, Co, Cu, Zn, Nb, Ba, La, Sm, Gd, Yb, and Bi)are presented for dopant concentrations in the range $0.00 \leq x \leq 0.25$. The stability of the dopants is compared and 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, but are shown to present no simple trend. The previously observed inverse relation between the thermal expansion coefficient and the bulk modulus is shown to persist independent of the inclusion of charge compensating vacancies.

preprint2013arXiv

Tuning of CeO$_2$ buffer layers for coated superconductors through doping

The appearance of microcracks in CeO$_2$ buffer layers, as used in buffer layer architectures for coated superconductors, indicates the presence of stress between this buffer layer and the substrate. This stress can originate from the differences in thermal expansion or differences in lattice parameters between the CeO$_2$ buffer layer and the substrate. In this article, we study, by means of \textit{ab initio} density functional theory calculations, the influence of group IV doping elements on the lattice parameter and bulk modulus of CeO$_2$. Vegard's law behavior is found for the lattice parameter in systems without oxygen vacancies, and the Shannon crystal radii for the doping elements are retrieved from the lattice expansions. We show that the lattice parameter of the doped CeO$_2$ can be matched to that of the La$_2$Zr$_2$O$_7$ coated NiW substrate substrate for dopant concentrations of about $5\%$, and that bulk modulus matching is either not possible or would require extreme doping concentrations.

preprint2011arXiv

Exact Ionization Potentials from Wavefunction Asymptotics: The Extended Koopmans' Theorem, Revisited

A simple explanation is given for the exactness of the extended Koopmans' theorem for computing the removal energy of any many-electron system to the lowest-energy ground state ion of a given symmetry. In particular, by removing the electron from a "removal orbital" of appropriate symmetry that is concentrated in the asymptotic region, one obtains the exact ionization potential and the exact Dyson orbital for the corresponding state of the ion. It is argued that the EKT is not restricted to many-electron systems but holds for any finite many-body system, provided the interaction vanishes for increasing interparticle distance. A necessary and sufficient condition for the validity of the extended Koopmans' theorem for any state (not just the lowest-energy states of a given symmetry) in terms of the third-order reduced density matrix is stated and derived.

preprint2011arXiv

Fast density-matrix based partitioning of the energy over the atoms in a molecule consistent with the Hirshfeld-I partitioning of the electron density

For the Hirshfeld-I atom-in-molecule model, associated single-atom energies and interaction energies at the Hartree-Fock level are determined efficiently in one-electron Hilbert space. In contrast to most other approaches, the energy terms are fully consistent with the partitioning of the underlying one-electron density matrix. Starting from the Hirshfeld-I atom-in-molecule model for the electron density, the molecular one-electron density matrix is partitioned with a previously introduced double-atom scheme [Vanfleteren D. et al., J Chem Phys 2010, 132, 164111]. Single-atom density matrices are constructed from the atomic and bond contributions of the double-atom scheme. Since the Hartree-Fock energy can be expressed solely in terms of the one-electron density matrix, the partitioning of the latter over the atoms in the molecule leads naturally to a corresponding partitioning of the Hartree-Fock energy. When the size of the molecule or the molecular basis set does not grow too large, the method shows considerable computational advantages compared to other approaches that require cumbersome numerical integration of the molecular energy integrals weighted by atomic weight functions.

preprint2011arXiv

Hilbert-space partitioning of the molecular one-electron density matrix with orthogonal projectors

A double-atom partitioning of the molecular one-electron density matrix is used to describe atoms and bonds. All calculations are performed in Hilbert space. The concept of atomic weight functions (familiar from Hirshfeld analysis of the electron density) is extended to atomic weight matrices. These are constructed to be orthogonal projection operators on atomic subspaces, which has significant advantages in the interpretation of the bond contributions. In close analogy to the iterative Hirshfeld procedure, self-consistency is built in at the level of atomic charges and occupancies. The method is applied to a test set of about 67 molecules, representing various types of chemical binding. A close correlation is observed between the atomic charges and the Hirshfeld-I atomic charges.

preprint2011arXiv

Partitioning of the molecular density matrix over atoms and bonds

A double-index atomic partitioning of the molecular first-order density matrix is proposed. Contributions diagonal in the atomic indices correspond to atomic density matrices, whereas off-diagonal contributions carry information about the bonds. The resulting matrices have good localization properties, in contrast to single-index atomic partitioning schemes of the molecular density matrix. It is shown that the electron density assigned to individual atoms, when derived from the density matrix partitioning, can be made con- sistent with well-known partitions of the electron density over AIM basins, either with sharp or with fuzzy boundaries. The method is applied to a test set of about 50 molecules, representative for various types of chemical binding. A close correlation is observed between the trace of the bond matrices and the SEDI (shared electron density index) bond index.

preprint2011arXiv

Variational density matrix optimization using semidefinite programming

We discuss how semidefinite programming can be used to determine the second-order density matrix directly through a variational optimization. We show how the problem of characterizing a physical or N -representable density matrix leads to matrix-positivity constraints on the density matrix. We then formulate this in a standard semidefinite programming form, after which two interior point methods are discussed to solve the SDP. As an example we show the results of an application of the method on the isoelectronic series of Beryllium.

preprint2010arXiv

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

A previous study of diatomic molecules revealed that variational second-order density matrix theory has serious problems in the dissociation limit when the N-representability is imposed at the level of the usual two-index (P, Q, G) or even three-index (T1, T2) conditions [H. van Aggelen et al., Phys. Chem. Chem. Phys. 11, 5558 (2009)]. Heteronuclear molecules tend to dissociate into fractionally charged atoms. In this paper we introduce a general class of N-representability conditions, called subsystem constraints, and show that they cure the dissociation problem at little additional computational cost. As a numerical example the singlet potential energy surface of BeB+ is studied. The extension to polyatomic molecules, where more subsystem choices can be identified, is also discussed.