BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Dimitri Van Neck

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

ResearcherUnclaimed source record
physics.chem-phcond-mat.str-elphysics.comp-phquant-phmath-phmath.MPnlin.SIcond-mat.mes-hallcond-mat.stat-mechcond-mat.supr-connucl-thphysics.atm-clusphysics.atom-ph

Catalog footprint

What is connected

32works
13topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

DMRG-CASPT2 study of the longitudinal static second hyperpolarizability of all-trans polyenes

We have implemented internally contracted complete active space second order perturbation theory (CASPT2) with the density matrix renormalization group (DMRG) as active space solver [Y. Kurashige and T. Yanai, J. Chem. Phys. 135, 094104 (2011)]. Internally contracted CASPT2 requires to contract the generalized Fock matrix with the 4-particle reduced density matrix (4-RDM) of the reference wavefunction. The required 4-RDM elements can be obtained from 3-particle reduced density matrices (3-RDM) of different wavefunctions, formed by symmetry-conserving single-particle excitations op top of the reference wavefunction. In our spin-adapted DMRG code chemps2 [https://github.com/sebwouters/chemps2], we decompose these excited wavefunctions as spin-adapted matrix product states, and calculate their 3-RDM in order to obtain the required contraction of the generalized Fock matrix with the 4-RDM of the reference wavefunction. In this work, we study the longitudinal static second hyperpolarizability of all-trans polyenes C$_{2n}$H$_{2n+2}$ [n = 4 - 12] in the cc-pVDZ basis set. DMRG-SCF and DMRG-CASPT2 yield substantially lower values and scaling with system size compared to RHF and MP2, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Read-Green resonances in a topological superconductor coupled to a bath

We study a topological superconductor capable of exchanging particles with an environment. This additional interaction breaks particle-number symmetry and can be modelled by means of an integrable Hamiltonian, building on the class of Richardson-Gaudin pairing models. The isolated system supports zero-energy modes at a topological phase transition, which disappear when allowing for particle exchange with an environment. However, it is shown from the exact solution that these still play an important role in system-environment particle exchange, which can be observed through resonances in low-energy and -momentum level occupations. These fluctuations signal topologically protected Read-Green points and cannot be observed within traditional mean-field theory.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An eigenvalue-based method and determinant representations for general integrable XXZ Richardson-Gaudin models

We propose an extension of the numerical approach for integrable Richardson-Gaudin models based on a new set of eigenvalue-based variables. Starting solely from the Gaudin algebra, the approach is generalized towards the full class of XXZ Richardson-Gaudin models. This allows for a fast and robust numerical determination of the spectral properties of these models, avoiding the singularities usually arising at the so-called singular points. We also provide different determinant expressions for the normalization of the Bethe Ansatz states and form factors of local spin operators, opening up possibilities for the study of larger systems, both integrable and non-integrable. These expressions can be written in terms of the new set of variables and generalize the results previously obtained for rational Richardson-Gaudin models and Dicke-Jaynes-Cummings-Gaudin models. Remarkably, these results are independent of the explicit parametrization of the Gaudin algebra, exposing a universality in the properties of Richardson-Gaudin integrable systems deeply linked to the underlying algebraic structure.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Eigenvalue-based determinants for scalar products and form factors in Richardson-Gaudin integrable models coupled to a bosonic mode

Starting from integrable $su(2)$ (quasi-)spin Richardson-Gaudin XXZ models we derive several properties of integrable spin models coupled to a bosonic mode. We focus on the Dicke-Jaynes-Cummings-Gaudin models and the two-channel $(p+ip)$-wave pairing Hamiltonian. The pseudo-deformation of the underlying $su(2)$ algebra is here introduced as a way to obtain these models in the contraction limit of different Richardson-Gaudin models. This allows for the construction of the full set of conserved charges, the Bethe Ansatz state, and the resulting Richardson-Gaudin equations. For these models an alternative and simpler set of quadratic equations can be found in terms of the eigenvalues of the conserved charges. Furthermore, the recently proposed eigenvalue-based determinant expressions for the overlaps and form factors of local operators are extended to these models, linking the results previously presented for the Dicke-Jaynes-Cummings-Gaudin models with the general results for Richardson-Gaudin XXZ models.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

CheMPS2: a free open-source spin-adapted implementation of the density matrix renormalization group for ab initio quantum chemistry

The density matrix renormalization group (DMRG) has become an indispensable numerical tool to find exact eigenstates of finite-size quantum systems with strong correlation. In the fields of condensed matter, nuclear structure and molecular electronic structure, it has significantly extended the system sizes that can be handled compared to full configuration interaction, without losing numerical accuracy. For quantum chemistry (QC), the most efficient implementations of DMRG require the incorporation of particle number, spin and point group symmetries in the underlying matrix product state (MPS) ansatz, as well as the use of so-called complementary operators. The symmetries introduce a sparse block structure in the MPS ansatz and in the intermediary contracted tensors. If a symmetry is non-abelian, the Wigner-Eckart theorem allows to factorize a tensor into a Clebsch-Gordan coefficient and a reduced tensor. In addition, the fermion signs have to be carefully tracked. Because of these challenges, implementing DMRG efficiently for QC is not straightforward. Efficient and freely available implementations are therefore highly desired. In this work we present CheMPS2, our free open-source spin-adapted implementation of DMRG for ab initio QC. Around CheMPS2, we have implemented the augmented Hessian Newton-Raphson complete active space self-consistent field method, with exact Hessian. The bond dissociation curves of the 12 lowest states of the carbon dimer were obtained at the DMRG(28 orbitals, 12 electrons, D$_{\mathsf{SU(2)}}$=2500)/cc-pVDZ level of theory. The contribution of $1s$ core correlation to the $X^1Σ_g^+$ bond dissociation curve of the carbon dimer was estimated by comparing energies at the DMRG(36o, 12e, D$_{\mathsf{SU(2)}}$=2500)/cc-pCVDZ and DMRG-SCF(34o, 8e, D$_{\mathsf{SU(2)}}$=2500)/cc-pCVDZ levels of theory.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

DMRG-SCF study of the singlet, triplet, and quintet states of oxo-Mn(Salen)

We use CheMPS2, our free open-source spin-adapted implementation of the density matrix renormalization group (DMRG) [Wouters et al., Comput. Phys. Commun. 185, 1501 (2014)], to study the lowest singlet, triplet, and quintet states of the oxo-Mn(Salen) complex. We describe how an initial approximate DMRG calculation in a large active space around the Fermi level can be used to obtain a good set of starting orbitals for subsequent complete-active-space or DMRG self-consistent field (CASSCF or DMRG-SCF) calculations. This procedure mitigates the need for a localization procedure, followed by a manual selection of the active space. Per multiplicity, the same active space of 28 electrons in 22 orbitals (28e, 22o) is obtained with the 6-31G*, cc-pVDZ, and ANO-RCC-VDZP basis sets (the latter with DKH2 scalar relativistic corrections). Our calculations provide new insight into the electronic structure of the quintet.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Exact solution of the $p_x + ip_y$ pairing Hamiltonian by deforming the pairing algebra

The present paper makes a connection between collective bosonic states and the exact solutions of the $p_x + ip_y$ pairing Hamiltonian. This makes it possible to investigate the effects of the Pauli principle on the energy spectrum, by gradually reintroducing the Pauli principle. It also introduces an efficient and stable numerical method to probe all the eigenstates of this class of Hamiltonians.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Projector quantum Monte Carlo with matrix product states

We marry tensor network states (TNS) and projector quantum Monte Carlo (PMC) to overcome the high computational scaling of TNS and the sign problem of PMC. Using TNS as trial wavefunctions provides a route to systematically improve the sign structure and to eliminate the bias in fixed-node and constrained-path PMC. As a specific example, we describe phaseless auxiliary-field quantum Monte Carlo with matrix product states (MPS-AFQMC). MPS-AFQMC improves significantly on the DMRG ground-state energy. For the J1-J2 model on two-dimensional square lattices, we observe with MPS-AFQMC an order of magnitude reduction in the error for all couplings, compared to DMRG. The improvement is independent of walker bond dimension, and we therefore use bond dimension one for the walkers. The computational cost of MPS-AFQMC is then quadratic in the bond dimension of the trial wavefunction, which is lower than the cubic scaling of DMRG. The error due to the constrained-path bias is proportional to the variational error of the trial wavefunction. We show that for the J1-J2 model on two-dimensional square lattices, a linear extrapolation of the MPS-AFQMC energy with the discarded weight from the DMRG calculation allows to remove the constrained-path bias. Extensions to other tensor networks are briefly discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The density matrix renormalization group for ab initio quantum chemistry

During the past 15 years, the density matrix renormalization group (DMRG) has become increasingly important for ab initio quantum chemistry. Its underlying wavefunction ansatz, the matrix product state (MPS), is a low-rank decomposition of the full configuration interaction tensor. The virtual dimension of the MPS, the rank of the decomposition, controls the size of the corner of the many-body Hilbert space that can be reached with the ansatz. This parameter can be systematically increased until numerical convergence is reached. The MPS ansatz naturally captures exponentially decaying correlation functions. Therefore DMRG works extremely well for noncritical one-dimensional systems. The active orbital spaces in quantum chemistry are however often far from one-dimensional, and relatively large virtual dimensions are required to use DMRG for ab initio quantum chemistry (QC-DMRG). The QC-DMRG algorithm, its computational cost, and its properties are discussed. Two important aspects to reduce the computational cost are given special attention: the orbital choice and ordering, and the exploitation of the symmetry group of the Hamiltonian. With these considerations, the QC-DMRG algorithm allows to find numerically exact solutions in active spaces of up to 40 electrons in 40 orbitals.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Dicke model as the contraction limit of a pseudo-deformed Richardson-Gaudin model

The Dicke model is derived in the contraction limit of a pseudo-deformation of the quasispin algebra in the su(2)-based Richardson-Gaudin models. Likewise, the integrability of the Dicke model is established by constructing the full set of conserved charges, the form of the Bethe Ansatz state, and the associated Richardson-Gaudin equations. Thanks to the formulation in terms of the pseudo-deformation, the connection from the su(2)-based Richardson-Gaudin model towards the Dicke model can be performed adiabatically.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Variational optimization of the 2DM: approaching three-index accuracy using extended cluster constraints

The reduced density matrix is variationally optimized for the two-dimensional Hubbard model. Exploiting all symmetries present in the system, we have been able to study $6\times6$ lattices at various fillings and different values for the on-site repulsion, using the highly accurate but computationally expensive three-index conditions. To reduce the computational cost we study the performance of imposing the three-index constraints on local clusters of $2\times2$ and $3\times3$ sites. We subsequently derive new constraints which extend these cluster constraints to incorporate the open-system nature of a cluster on a larger lattice. The feasibility of implementing these new constraints is demonstrated by performing a proof-of-principle calculation on the $6\times6$ lattice. It is shown that a large portion of the three-index result can be recovered using these extended cluster constraints, at a fraction of the computational cost.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Linear Response Theory for the Density Matrix Renormalization Group: Efficient Algorithms for Strongly Correlated Excited States

Linear response theory for the density matrix renormalization group (DMRG-LRT) was first presented in terms of the DMRG renormalization projectors [Dorando et al., J. Chem. Phys. 130, 184111 (2009)]. Later, with an understanding of the manifold structure of the matrix product state (MPS) ansatz, which lies at the basis of the DMRG algorithm, a way was found to construct the linear response space for general choices of the MPS gauge in terms of the tangent space vectors [Haegeman et al., Phys. Rev. Lett. 107, 070601 (2011)]. These two developments led to the formulation of the Tamm-Dancoff and random phase approximations (TDA and RPA) for MPS. This work describes how these LRTs may be efficiently implemented through minor modifications of the DMRG sweep algorithm, at a computational cost which scales the same as the ground-state DMRG algorithm. In fact, the mixed canonical MPS form implicit to the DMRG sweep is essential for efficient implementation of the RPA, due to the structure of the second-order tangent space. We present ab initio DMRG-TDA results for excited states of polyenes, the water molecule, and a [2Fe-2S] iron-sulfur cluster.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Thouless theorem for matrix product states and subsequent post-density matrix renormalization group methods

The similarities between Hartree-Fock (HF) theory and the density-matrix renormalization group (DMRG) are explored. Both methods can be formulated as the variational optimization of a wave-function ansatz. Linearization of the time-dependent variational principle near a variational minimum allows to derive the random phase approximation (RPA). We show that the non-redundant parametrization of the matrix product state (MPS) tangent space [J. Haegeman et al., Phys. Rev. Lett. 107, 070601 (2011)] leads to the Thouless theorem for MPS, i.e. an explicit non-redundant parametrization of the entire MPS manifold, starting from a specific MPS reference. Excitation operators are identified, which extends the analogy between HF and DMRG to the Tamm-Dancoff approximation (TDA), the configuration interaction (CI) expansion, and coupled cluster theory. For a small one-dimensional Hubbard chain, we use a CI-MPS ansatz with single and double excitations to improve on the ground state and to calculate low-lying excitation energies. For a symmetry-broken ground state of this model, we show that RPA-MPS allows to retrieve the Goldstone mode. We also discuss calculations of the RPA-MPS correlation energy. With the long-range quantum chemical Pariser-Parr-Pople Hamiltonian, low-lying TDA-MPS and RPA-MPS excitation energies for polyenes are obtained.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Perturbations on the superconducting state of metallic nanoparticles: influence of geometry and impurities

The pair condensation energy of a finite-size superconducting particle is studied as a function of two control parameters. The first control parameter is the shape of the particle, and the second parameter is a position-dependent impurity introduced in the particle. Whereas the former parameter is known to induce strong fluctuations in the condensation energy, the latter control parameter is found to be a more gentle probe of the pairing correlations.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Faddeev Random Phase Approximation for Molecules

The Faddeev Random Phase Approximation is a Green's function technique that makes use of Faddeev-equations to couple the motion of a single electron to the two-particle--one-hole and two-hole--one-particle excitations. This method goes beyond the frequently used third-order Algebraic Diagrammatic Construction method: all diagrams involving the exchange of phonons in the particle-hole and particle-particle channel are retained, but the phonons are described at the level of the Random Phase Approximation. This paper presents the first results for diatomic molecules at equilibrium geometry. The behavior of the method in the dissociation limit is also investigated.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Electron penetration in the nucleus and its effect on the quadrupole interaction

A series expansion of the interaction between a nucleus and its surrounding electron distribution provides terms that are well-known in the study of hyperfine interactions: the familiar quadrupole interaction and the less familiar hexadecapole interaction. If the penetration of electrons into the nucleus is taken into account, various corrections to these multipole interactions appear. The best known one is a scalar correction related to the isotope shift and the isomer shift. This paper discusses a related tensor correction, which modifies the quadrupole interaction if electrons penetrate the nucleus: the quadrupole shift. We describe the mathematical formalism and provide first-principles calculations of the quadrupole shift for a large set of solids. Fully relativistic calculations that explicitly take a finite nucleus into account turn out to be mandatory. Our analysis shows that the quadrupole shift becomes appreciably large for heavy elements. Implications for experimental high-precision studies of quadrupole interactions and quadrupole moment ratios are discussed. A literature review of other small quadrupole-like effects is presented as well.

0 citations0 reviewsNo code linkedOpen access