Researcher profile

Dimitri Van Neck

Dimitri Van Neck contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
4topics
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

10 published item(s)

preprint2022arXiv

Bivariational Principle for an Antisymmetrized Product of Nonorthogonal Geminals Appropriate for Strong Electron Correlation

We develop a bivariational principle for an antisymmetric product of nonorthogonal geminals. Special cases reduce to the antisymmetric product of strongly-orthogonal geminals (APSG), the generalized valence bond-perfect pairing (GVB-PP), and the antisymmetrized geminal power (AGP) wavefunctions. The presented method employs wavefunctions of the same type as Richardson-Gaudin (RG) states, but which are not eigenvectors of a model Hamiltonian which would allow for more freedom in the mean-field. The general idea is to work with the same state in a primal picture in terms of pairs, and in a dual picture in terms of pair-holes. This leads to an asymmetric energy expression which may be optimized bivariationally, and is strictly variational when the two representations are consistent. The general approach may be useful in other contexts, such as for computationally feasible variational coupled-cluster methods.

preprint2021arXiv

The seniority quantum number in Tensor Network States

We employ tensor network methods for the study of the seniority quantum number - defined as the number of unpaired electrons in a many-body wave function - in molecular systems. Seniority-zero methods recently emerged as promising candidates to treat strong static correlations in molecular systems, but are prone to deficiencies related to dynamical correlation and dispersion. We systematically resolve these deficiencies by increasing the allowed seniority number using tensor network methods. In particular, we investigate the number of unpaired electrons needed to correctly describe the binding of the neon and nitrogen dimer and the $D_{6h}$ symmetry of benzene.

preprint2012arXiv

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

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

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

High-accuracy ab-initio quantum chemistry by means of an SU(2) x U(1) invariant matrix product state Ansatz: the static second hyperpolarizability

We have implemented the single-site density matrix renormalization group algorithm for the variational optimization of SU(2) \times U(1) (spin and particle number) invariant matrix product states for general spin and particle number symmetric fermionic Hamiltonians. This class also 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 obtained with the algorithm. A comparison with other methods is made.

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.

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.