Researcher profile

Helen van Aggelen

Helen van Aggelen contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

6 published item(s)

preprint2013arXiv

Equivalence of Particle-Particle Random Phase Approximation Correlation Energy and Ladder-Coupled-Cluster-Double

We present an analytical proof and numerical demonstrations of the equivalence of the correlation energy from particle-particle random phase approximation (pp-RPA) and ladder-couple-cluster-doubles (ladder-CCD). These two theories reduce to the identical algebraic matrix equation and correlation energy expressions, under the assumption that the pp-RPA equation is stable. The numerical examples illustrate that the correlation energy missed by pp-RPA in comparison with couple-cluster single and double is largely canceled out when considering reaction energies. This theoretical connection will be beneficial to future pp-RPA studies based on the well established couple cluster theory.

preprint2013arXiv

Extensive v2DM study of the one-dimensional Hubbard model for large lattice sizes: Exploiting translational invariance and parity

Using variational density matrix optimization with two- and three-index conditions we study the one-dimensional Hubbard model with periodic boundary conditions at various filling factors. Special attention is directed to the full exploitation of the available symmetries, more specifically the combination of translational invariance and space-inversion parity, which allows for the study of large lattice sizes. We compare the computational scaling of three different semidefinite programming algorithms with increasing lattice size, and find the boundary point method to be the most suited for this type of problem. Several physical properties, such as the two-particle correlation functions, are extracted to check the physical content of the variationally determined density matrix. It is found that the three-index conditions are needed to correctly describe the full phase diagram of the Hubbard model. We also show that even in the case of half filling, where the ground-state energy is close to the exact value, other properties such as the spin-correlation function can be flawed.

preprint2012arXiv

Using tensor hypercontraction density fitting to achieve an O(L^4) CISD algorithm

Recently, Hohenstein et al[1] introduced tensor hypercontraction density fitting to decompose the rank-4 electron repulsion integral tensor as the product of five rank-2 tensors. In this paper, we use this methodology to construct an algorithm which calculates the approximate ground state energy in O(L^4) operations. We test our method using several small molecules and show that we quickly approach the CISD limit with a small number of auxiliary functions.

preprint2012arXiv

Variational two-particle density matrix calculation for the Hubbard model below half filling using spin-adapted lifting conditions

The variational determination of the two-particle density matrix is an interesting, but not yet fully explored technique that allows to obtain ground-state properties of a quantum many-body system without reference to an $N$-particle wave function. The one-dimensional fermionic Hubbard model has been studied before with this method, using standard two- and three-index conditions on the density matrix [J. R. Hammond {\it et al.}, Phys. Rev. A 73, 062505 (2006)], while a more recent study explored so-called subsystem constraints [N. Shenvi {\it et al.}, Phys. Rev. Lett. 105, 213003 (2010)]. These studies reported good results even with only standard two-index conditions, but have always been limited to the half-filled lattice. In this Letter we establish the fact that the two-index approach fails for other fillings. In this case, a subset of three-index conditions is absolutely needed to describe the correct physics in the strong-repulsion limit. We show that applying lifting conditions [J.R. Hammond {\it et al.}, Phys. Rev. A 71, 062503 (2005)] is the most economical way to achieve this, while still avoiding the computationally much heavier three-index conditions. A further extension to spin-adapted lifting conditions leads to increased accuracy in the intermediate repulsion regime. At the same time we establish the feasibility of such studies to the more complicated phase diagram in two-dimensional Hubbard models.

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.