Source author record

Stijn De Baerdemacker

Stijn De Baerdemacker 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

22works
11topics
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

22 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

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.

preprint2020arXiv

The decomposition of an arbitrary $2^w\times 2^w$ unitary matrix into signed permutation matrices

Birkhoff's theorem tells that any doubly stochastic matrix can be decomposed as a weighted sum of permutation matrices. A similar theorem reveals that any unitary matrix can be decomposed as a weighted sum of complex permutation matrices. Unitary matrices of dimension equal to a power of~2 (say $2^w$) deserve special attention, as they represent quantum qubit circuits. We investigate which subgroup of the signed permutation matrices suffices to decompose an arbitrary such matrix. It turns out to be a matrix group isomorphic to the extraspecial group {\bf E}$_{2^{2w+1}}^+$ of order $2^{2w+1}$. An associated projective group of order $2^{2w}$ equally suffices.

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.

preprint2016arXiv

The block-ZXZ synthesis of an arbitrary quantum circuit

Given an arbitrary $2^w \times 2^w$ unitary matrix $U$, a powerful matrix decomposition can be applied, leading to four different syntheses of a $w$-qubit quantum circuit performing the unitary transformation. The demonstration is based on a recent theorem by Führ and Rzeszotnik, generalizing the scaling of single-bit unitary gates ($w=1$) to gates with arbitrary value of~$w$. The synthesized circuit consists of controlled 1-qubit gates, such as NEGATOR gates and PHASOR gates. Interestingly, the approach reduces to a known synthesis method for classical logic circuits consisting of controlled NOT gates, in the case that $U$ is a permutation matrix.

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.

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.

preprint2015arXiv

From reversible computation to quantum computation by Lagrange interpolation

Classical reversible circuits, acting on $w$~bits, are represented by permutation matrices of size $2^w \times 2^w$. Those matrices form the group P($2^w$), isomorphic to the symmetric group {\bf S}$_{2^w}$. The permutation group P($n$), isomorphic to {\bf S}$_n$, contains cycles with length~$p$, ranging from~1 to $L(n)$, where $L(n)$ is the so-called Landau function. By Lagrange interpolation between the $p$~matrices of the cycle, we step from a finite cyclic group of order~$p$ to a 1-dimensional Lie group, subgroup of the unitary group U($n$). As U($2^w$) is the group of all possible quantum circuits, acting on $w$~qubits, such interpolation is a natural way to step from classical computation to quantum computation.

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.

preprint2015arXiv

Scaling a unitary matrix

The iterative method of Sinkhorn allows, starting from an arbitrary real matrix with non-negative entries, to find a so-called 'scaled matrix' which is doubly stochastic, i.e. a matrix with all entries in the interval (0, 1) and with all line sums equal to 1. We conjecture that a similar procedure exists, which allows, starting from an arbitrary unitary matrix, to find a scaled matrix which is unitary and has all line sums equal to 1. The existence of such algorithm guarantees a powerful decomposition of an arbitrary quantum circuit.

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

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.

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.

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.

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.

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.

preprint2012arXiv

Richardson-Gaudin integrability in the contraction limit of the quasispin

Background: The reduced, level-independent, Bardeen-Cooper-Schrieffer Hamiltonian is exactly diagonalizable by means of a Bethe Ansatz wavefunction, provided the free variables in the Ansatz are the solutions of the set of Richardson-Gaudin equations. On the one side, the Bethe Ansatz is a simple product state of generalised pair operators. On the other hand, the Richardson-Gaudin equations are strongly coupled in a non-linear way, making them prone to singularities. Unfortunately, it is non-trivial to give a clear physical interpretation to the Richardson-Gaudin variables because no physical operator is directly related to the individual variables. Purpose: The purpose of this paper is to shed more light on the critical behavior of the Richardson-Gaudin equations, and how this is related to the product wave structure of the Bethe Ansatz. Method: A pseudo-deformation of the quasi-spin algebra is introduced, leading towards a Heisenberg-Weyl algebra in the contraction limit of the deformation parameter. This enables an adiabatic connection of the exact Bethe Ansatz eigenstates with pure bosonic multiphonon states. The physical interpretation of this approach is an adiabatic suppression of the Pauli exclusion principle. Results: The method is applied to a so-called "picket-fence" model for the BCS Hamiltonian, displaying a typical critical behavior in the Richardson-Gaudin variables. It was observed that the associated bosonic multiphonon states change collective nature at the critical interaction strengths of the Richardson-Gaudin equations. Conclusions: The Pauli exclusion principle is the main responsible for the critical behavior of the Richardson-Gaudin equations, which can be suppressed by means of a pseudo deformation of the quasispin algebra.

preprint2010arXiv

The Tamm-Dancoff Approximation as the boson limit of the Richardson-Gaudin equations for pairing

A connection is made between the exact eigen states of the BCS Hamiltonian and the predictions made by the Tamm-Dancoff Approximation. This connection is made by means of a parametrised algebra, which gives the exact quasi-spin algebra in one limit of the parameter and the Heisenberg-Weyl algebra in the other. Using this algebra to construct the Bethe Ansatz solution of the BCS Hamiltonian, we obtain parametrised Richardson-Gaudin equations, leading to the secular equation of the Tamm-Dancoff Approximation in the bosonic limit. An example is discussed in depth.