Source author record

Örs Legeza

Örs Legeza 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

29works
9topics
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

29 published item(s)

preprint2022arXiv

Compressing multireference character of wave functions via fermionic mode optimization

In this work, we present a brief overview of the fermionic mode optimization within the framework of tensor network state methods [PRL 117,210402(2016)], and demonstrate that it has the potential to compress the multireference character of the wave functions after finding optimal molecular orbitals (modes), based on entanglement minimization. Numerical simulations have been performed for the nitrogen dimer in the cc-pVDZ basis for the equilibrium and for stretched geometries.

preprint2022arXiv

Exact matrix product state representation and convergence of a fully correlated electronic wavefunction in the infinite basis limit

In this article we present the exact representation of a fully correlated electronic wavefunction as the single-particle basis approaches completeness. It consists of a half-infinite chain of matrices of exponentially increasing size. The complete basis limit is illustrated numerically using the density matrix renormalization group method by computing the core-valence entanglement in the C$_2$ ground state in increasing subsets of cc-pVTZ and pVQZ bases until convergence is reached.

preprint2022arXiv

Towards large-scale restricted active space calculations inspired by the Schmidt decomposition

We present an alternative, memory-efficient, Schmidt decomposition-based description of the inherently bipartite restricted active space (RAS) scheme, which can be implemented effortlessly within the density matrix renormalization group (DMRG) method via the dynamically extended active space procedure. Benchmark calculations are compared against state-of-the-art results of C$_2$ and Cr$_2$, which are notorious for their multi-reference character. Our results for ground and excited states together with spectroscopic constants demonstrate that the proposed novel approach, dubbed as DMRG-RAS, which is variational and free of uncontrolled method errors, has the potential to outperfom conventional methods for strongly correlated molecules.

preprint2021arXiv

How creating one additional well can generate Bose-Einstein condensation

The realization of Bose-Einstein condensation in ultracold trapped gases has led to a revival of interest in that fascinating quantum phenomenon. This experimental achievement necessitated both extremely low temperatures and sufficiently weak interactions. Particularly in reduced spatial dimensionality even an infinitesimal interaction immediately leads to a departure to quasi-condensation. We propose a system of strongly interacting bosons which overcomes those obstacles by exhibiting a number of intriguing related features: (i) The tuning of just a single control parameter drives a transition from quasi-condensation to complete condensation, (ii) the destructive influence of strong interactions is compensated by the respective increased mobility, (iii) topology plays a crucial role since a crossover from one- to `infinite'-dimensionality is simulated, (iv) a ground state gap opens which makes the condensation robust to thermal noise. Remarkably, all these features can be derived by analytical and exact numerical means despite the non-perturbative character of the system.

preprint2021arXiv

Simulating Lindbladian evolution with non-abelian symmetries: Ballistic front propagation in the $SU(2)$ Hubbard model with a localized loss

We develop a non-Abelian time evolving block decimation (NA-TEBD) approach to study of open systems governed by Lindbladian time evolution, while exploiting an arbitrary number of abelian or non-abelian symmetries. We illustrate this method in a one-dimensional fermionic $SU(2)$ Hubbard model on a semi-infinite lattice with localized particle loss at one end. We observe a ballistic front propagation with strongly renormalized front velocity, and a hydrodynamic current density profile. For large loss rates, a suppression of the particle current is observed, as a result of the quantum Zeno effect. Operator entanglement is found to propagate faster than the depletion profile, preceding the latter.

preprint2021arXiv

Topologically protected, correlated end spin formation in carbon nanotubes

For most chiralities, semiconducting nanotubes display topologically protected end states of multiple degeneracies. We demonstrate using density matrix renormalization group based quantum chemistry tools that the presence of Coulomb interactions induces the formation of robust end spins. These are the close analogues of ferromagnetic edge states emerging in graphene nanoribbons. The interaction between the two ends is sensitive to the length of the nanotube, its dielectric constant, as well as the size of the end spins: for $S=1/2$ end spins their interaction is antiferromagnetic, while for $S>1/2$ it changes from antiferromagnetic to ferromagnetic as the nanotube length increases. The interaction between end spins can be controlled by changing the dielectric constant of the environment, thereby providing a possible platform for two-spin quantum manipulations.

preprint2020arXiv

DMRG on top of plane-wave Kohn-Sham orbitals: case study of defected boron nitride

In this paper, we analyze the numerical aspects of the inherently multi-reference density matrix renormalization group (DMRG) calculations on top of the periodic Kohn-Sham density functional theory (DFT) using the complete active space (CAS) approach. Following the technical outline related to the computation of the Hamiltonian matrix elements and to the construction of the active space, we illustrate the potential of the framework by studying the vertical many-body energy spectrum of hexagonal boron nitride (hBN) nano-flakes embedding a single boron vacancy point defect with prominent multi-reference character. We investigate the consistency of the DMRG energy spectrum from the perspective of sample size, basis size, and active space selection protocol. Results obtained from standard quantum chemical atom-centered basis calculations and plane-wave based counterparts show excellent agreement. Furthermore, we also discuss the spectrum of the periodic sheet which is in good agreement with extrapolated data of finite clusters. These results pave the way toward applying DMRG method in extended correlated solid state systems, such as point qubit in wide band gap semiconductors.

preprint2020arXiv

Ground State of the Fe(II)-porphyrin Model System Corresponds to the Quintet State: A DFT and DMRG-based Tailored CC Study

Fe(II)-porphyrins (FeP) play an important role in many reactions relevant to material science and biological processes, due to their closely lying spin states. However, this small energetic separation also makes it challenging to establish the correct spin state ordering. Although the prevalent opinion is that these systems posses the triplet ground state, the recent experiment on Fe(II)-phthalocyanine under conditions matching those of an isolated molecule points toward the quintet ground state. We present a thorough study of FeP model by means of the density functional theory and density matrix renormalization group based tailored coupled clusters, in which we address all previously discussed correlation effects. We examine the importance of geometrical parameters, the Fe-N distances in particular, and conclude that the system possesses the quintet ground state, which is in our calculations well-separated from the triplet state.

preprint2020arXiv

Massively parallel quantum chemical density matrix renormalization group method

We present, to the best of our knowlegde, the first attempt to exploit the supercomputer platform for quantum chemical density matrix renormalization group (QC-DMRG) calculations. We have developed the parallel scheme based on the in-house MPI global memory library, which combines operator and symmetry sector parallelisms, and tested its performance on three different molecules, all typical candidates for QC-DMRG calculations. In case of the largest calculation, which is the nitrogenase FeMo cofactor cluster with the active space comprising 113 electrons in 76 orbitals and bond dimension equal to 6000, our parallel approach scales up to approximately 2000 CPU cores.

preprint2020arXiv

Near-linear Scaling in DMRG-based Tailored Coupled Clusters: An Implementation of DLPNO-TCCSD and DLPNO-TCCSD(T)

We present a new implementation of DMRG-based tailored coupled clusters method (TCCSD), which employs the domain-based local pair natural orbital approach (DLPNO-TCCSD). Compared to the previous LPNO version of the method, the new implementation is more accurate, offers more favorable scaling and provides more consistent behavior across the variety of systems. On top of the singles and doubles, we include the perturbative triples correction (T), which is able to retrieve even more dynamic correlation. The methods were tested on three systems: tetramethyleneethane, oxo-Mn(Salen) and Iron(II)-porphyrin model. The first two were revisited to assess the performance with respect to LPNO-TCCSD. For oxo-Mn(Salen), we retrieved between 99.8-99.9% of the total canonical correlation energy which is the improvement of 0.2% over the LPNO version in less than 63% of the total LPNO runtime. Similar results were obtained for Iron(II)-porphyrin. When the perturbative triples correction was employed, irrespective of the active space size or system, the obtained energy differences between two spin states were within the chemical accuracy of 1 kcal/mol using the default DLPNO settings.

preprint2020arXiv

Towards DMRG-tailored coupled cluster method in the 4c-relativistic domain

There are three essential problems in computational relativistic chemistry: electrons moving at relativistic speeds, close lying states and dynamical correlation. Currently available quantum-chemical methods are capable of solving systems with one or two of these issues. However, there is a significant class of molecules, in which all the three effects are present. These are the heavier transition metal compounds, lanthanides and actinides with open d or f shells. For such systems, sufficiently accurate numerical methods are not available, which hinders the application of theoretical chemistry in this field. In this paper, we combine two numerical methods in order to address this challenging class of molecules. These are the relativistic versions of coupled cluster methods and density matrix renormalization group (DMRG) method. To the best of our knowledge, this is the first relativistic implementation of the coupled cluster method externally corrected by DMRG. The method brings a significant reduction of computational costs, as we demonstrate on the system of TlH, AsH and SbH.

preprint2019arXiv

Quantum information-based analysis of electron-deficient bonds

Recently, the correlation theory of the chemical bond was developed, which applies concepts of quantum information theory for the characterization of chemical bonds, based on the multiorbital correlations within the molecule. Here for the first time, we extend the use of this mathematical toolbox for the description of electron-deficient bonds. We start by verifying the theory on the textbook example of a molecule with three-center two-electron bonds, namely the diborane(6). We then show that the correlation theory of the chemical bond is able to properly describe bonding situation in more exotic molecules which have been synthetized and characterized only recently, in particular the diborane molecule with four hydrogen atoms [diborane(4)] and neutral zerovalent s-block beryllium complex, whose surprising stability was attributed to a strong three-center two-electron $π$ bond stretching across the C-Be-C core. Our approach is of a high importance especially in the light of a constant chase after novel compounds with extraordinary properties where the bonding is expected to be unusual.

preprint2019arXiv

Symmetric single-impurity Kondo model on a tight-binding chain: a comparison of analytical and numerical ground-state approaches

We analyze the ground-state energy, local spin correlation, impurity spin polarization, impurity-induced magnetization, and corresponding zero-field susceptibilities of the symmetric single-impurity Kondo model on a tight-binding chain with bandwidth $W=2{\cal D}$ and coupling strength $J_{\rm K}$. We compare perturbative results and variational upper bounds from Yosida, Gutzwiller, and first-order Lanczos wave functions to the numerically exact data obtained from the Density-Matrix Renormalization Group (DMRG) and from the Numerical Renormalization Group (NRG) methods. The Gutzwiller variational approach becomes exact in the strong-coupling limit and reproduces the ground-state properties from DMRG and NRG for large couplings. We calculate the impurity spin polarization and its susceptibility in the presence of magnetic fields that are applied globally/locally to the impurity spin. The Yosida wave function provides qualitatively correct results in the weak-coupling limit. In DMRG, chains with about $10^3$ sites are large enough to describe the susceptibilities down to $J_{\rm K}/{\cal D}\approx 0.5$. For smaller Kondo couplings, only the NRG provides reliable results for a general host-electron density of states $ρ_0(ε)$. To compare with results from Bethe Ansatz, we study the impurity-induced magnetization and zero-field susceptibility. For small Kondo couplings, the zero-field susceptibilities at zero temperature approach $χ_0(J_{\rm K}\ll {\cal D})/(gμ_{\rm B})^2\approx \exp[1/(ρ_0(0)J_{\rm K})]/(2C{\cal D}\sqrt{πe ρ_0(0)J_{\rm K}})$, where $\ln(C)$ is the regularized first inverse moment of the density of states. Using NRG, we determine the universal sub-leading corrections up to second order in $ρ_0(0)J_{\rm K}$.

preprint2019arXiv

Towards the Efficient Local Tailored Coupled Cluster Approximation and the Peculiar Case of Oxo-Mn(Salen)

We introduce a new implementation of the coupled cluster method tailored by matrix product states wave functions (DMRG-TCCSD), which employs the local pair natural orbital approach (LPNO). By exploiting locality in the coupled cluster stage of the calculation, we were able to remove some of the limitations that hindered the application of the canonical version of the method to larger systems and/or with larger basis sets. We assessed the accuracy of the approximation using two systems: tetramethyleneethane (TME) and oxo-Mn(Salen). Using the default cut-off parameters, we were able to recover over 99.7% and 99.8% of canonical correlation energy for the triplet and singlet state of TME respectively. In case of oxo-Mn(Salen), we found out that the amount of retrieved canonical correlation energy depends on the size of the active space (CAS) - we retrieved over 99.6% for the larger 27 orbital CAS and over 99.8% for the smaller 22 orbital CAS. The use of LPNO-TCCSD allowed us to perform these calculations up to quadruple-$ζ$ basis set amounting to 1178 basis functions. Moreover, we examined dependance of the ground state of oxo-Mn(Salen) on CAS composition. We found out that the inclusion of 4d$_{xy}$ orbital plays an important role in stabilizing the singlet state at the DMRG-CASSCF level via double-shell effect. However, by including dynamic correlation the ground state was found to be triplet regardless of the size of the basis set or composition of CAS, which is in agreement with previous findings by canonical DMRG-TCCSD in smaller basis.

preprint2016arXiv

Analysis of two-orbital correlations in wavefunctions restricted to electron-pair states

Wavefunctions constructed from electron-pair states can accurately model strong electron correlation effects and are promising approaches especially for larger many-body systems. In this article, we analyze the nature and the type of electron correlation effects that can be captured by wavefunctions restricted to electron-pair states. We focus on the Antisymmetric Product of 1-reference orbital Geminal (AP1roG) method combined with an orbital optimization protocol presented in [Phys. Rev. B, 89, 201106(R), 2014] whose performance is assessed against electronic structures obtained form DMRG reference data. Our numerical analysis covers model systems for strong correlation: the one-dimensional Hubbard model with periodic boundary condition as well as metallic and molecular hydrogen rings. Specifically, the accuracy of AP1roG is benchmarked using the single-orbital entropy, the orbital-pair mutual information as well as the eigenvalue spectrum of the one-orbital and two-orbital reduced density matrices. Our study indicates that contributions from singly occupied states become important in the strong correlation regime which highlights the limitations of the AP1roG method. Furthermore, we examine the effect of orbital rotations within the AP1roG model on correlations between orbital pairs.

preprint2016arXiv

Coupled cluster method with single and double excitations tailored by matrix product state wave functions

In the last decade, the quantum chemical version of the density matrix renormalization group (DMRG) method has established itself as the method of choice for calculations of strongly correlated molecular systems. Despite its favourable scaling, it is in practice not suitable for computations of dynamic correlation. We present a novel method for accurate "post-DMRG" treatment of dynamic correlation based on the tailored coupled cluster (CC) theory in which the DMRG method is responsible for the proper description of non-dynamic correlation, whereas dynamic correlation is incorporated through the framework of the CC theory. We illustrate the potential of this method on prominent multireference systems, in particular N$_2$, Cr$_2$ molecules and also oxo-Mn(Salen) for which we have performed the first "post-DMRG" computations in order to shed light on the energy ordering of the lowest spin states.

preprint2016arXiv

On the Multi-Reference Nature of Plutonium Oxides: PuO$_2^{2+}$, PuO$_2$, PuO$_3$ and PuO$_2$(OH)$_2$

Actinide-containing complexes present formidable challenges for electronic structure methods due to the large number of degenerate or quasi-degenerate electronic states arising from partially occupied 5f and 6d shells. Conventional multi-reference methods can treat active spaces that are often at the upper limit of what is required for a proper treatment of species with complex electronic structures, leaving no room for verifying their suitability. In this work we address the issue of properly defining the active spaces in such calculations, and introduce a protocol to determine optimal active spaces based on the use of the Density Matrix Renormalization Group algorithm and concepts of quantum information theory. We apply the protocol to elucidate the electronic structure and bonding mechanism of volatile plutonium oxides (PuO$_3$ and PuO$_2$(OH)$_2$), species associated with nuclear safety issues for which little is known about the electronic structure and energetics. We show how, within a scalar relativistic framework, orbital-pair correlations can be used to guide the definition of optimal active spaces which provide an accurate description of static/non-dynamic electron correlation, as well as to analyse the chemical bonding beyond a simple orbital model. From this bonding analysis we are able to show that the addition of oxo- or hydroxo-groups to the plutonium dioxide species considerably changes the pi-bonding mechanism with respect to the bare triatomics, resulting in bent structures with considerable multi-reference character.

preprint2015arXiv

Dissecting the Bond Formation Process of $d^{10}$-Metal-Ethene Complexes with Multireference Approaches

The bonding mechanism of ethene to a nickel or palladium center is studied by the density matrix renormalization group algorithm, the complete active space self consistent field method, coupled cluster theory, and density functional theory. Specifically, we focus on the interaction between the metal atom and bis-ethene ligands in perpendicular and parallel orientations. The bonding situation in these structural isomers is further scrutinized using energy decomposition analysis and quantum information theory. Our study highlights the fact that when two ethene ligands are oriented perpendicular to each other, the complex is stabilized by the metal-to-ligand double-back-bonding mechanism. Moreover, we demonstrate that nickel-ethene complexes feature a stronger and more covalent interaction between the ligands and the metal center than palladium-ethene compounds with similar coordination spheres.

preprint2015arXiv

Hückel--Hubbard-Ohno modeling of $\boldsymbolπ$-bonds in ethene and ethyne with application to trans-polyacetylene

Quantum chemistry calculations provide the potential energy between two carbon atoms in ethane (H$_3$C$-$CH$_3$), ethene (H$_2$C$=$CH$_2$), and ethyne (HC$\equiv$CH) as a function of the atomic distance. Based on the energy function for the $σ$-bond in ethane, $V_σ(r)$, we use the Hückel model with Hubbard--Ohno interaction for the $π$~electrons to describe the energies $V_{σπ}(r)$ and $V_{σππ}(r)$ for the $σπ$ double bond in ethene and the $σππ$ triple bond in ethyne, respectively. The fit of the force functions shows that the Peierls coupling can be estimated with some precision whereas the Hubbard-Ohno parameters are insignificant at the distances under consideration. We apply the Hückel-Hubbard-Ohno model to describe the bond lengths and the energies of elementary electronic excitations of trans-polyacetylene, (CH)$_n$, and adjust the $σ$-bond potential for conjugated polymers.

preprint2015arXiv

On the calculation of complete dissociation curves of closed-shell pseudo-onedimensional systems through the multireference method of increments

The Method of Increments (MoI) has been employed using a multireference approach to calculate the dissociation curve of beryllium ring-shaped clusters Be$_n$ of different sizes. Benchmarks obtained through different single and multireference methods including the ab initio Density Matrix Renormalization Group (DMRG) were used to verify the validity of the MoI truncation which showed a reliable behavior for the whole dissociation curve. Moreover we investigated the size dependence of the correlation energy at different distances in order to extrapolate the values for the periodic chain and to discuss the transition from a metal-like to a insulating-like behavior of the wave function through quantum chemical considerations.

preprint2015arXiv

Second-order Peierls transition in the spin-orbital Kumar-Heisenberg model

We add a Heisenberg interaction term $\proptoλ$ in the one-dimensional SU(2)$\otimes$XY spin-orbital model introduced by B. Kumar. At $λ=0$ the spin and orbital degrees of freedom can be separated by a unitary transformation leading to an exact solution of the model. We show that a finite $λ>0$ leads to spontaneous dimerization of the system which in the thermodynamic limit becomes a smooth phase transition at $λ\to 0$, whereas it remains discontinuous within the first order perturbation approach. We present the behavior of the entanglement entropy, energy gap and dimerization order parameter in the limit of $λ\to 0$ confirming the critical behavior. Finally, we show the evidence of another phase transition in the Heisenberg limit, $λ\to\infty$, and give a qualitative analytical explanation of the observed dimerized states both in the limit of small and large $λ$.

preprint2014arXiv

A quantum informational approach for dissecting chemical reactions

We present a conceptionally different approach to dissect bond-formation processes in metal-driven catalysis using concepts from quantum information theory. Our method uses the entanglement and correlation among molecular orbitals to analyze changes in electronic structure that accompany chemical processes. As a proof-of-principle example, the evolution of nickel-ethene bond-formation is dissected which allows us to monitor the interplay of back-bonding and $π$-donation along the reaction coordinate. Furthermore, the reaction pathway of nickel-ethene complexation is analyzed using quantum chemistry methods revealing the presence of a transition state. Our study supports the crucial role of metal-to-ligand back-donation in the bond-forming process of nickel-ethene.

preprint2014arXiv

Tensor product methods and entanglement optimization for ab initio quantum chemistry

The treatment of high-dimensional problems such as the Schrödinger equation can be approached by concepts of tensor product approximation. We present general techniques that can be used for the treatment of high-dimensional optimization tasks and time-dependent equations, and connect them to concepts already used in many-body quantum physics. Based on achievements from the past decade, entanglement-based methods, -- developed from different perspectives for different purposes in distinct communities already matured to provide a variety of tools -- can be combined to attack highly challenging problems in quantum chemistry. The aim of the present paper is to give a pedagogical introduction to the theoretical background of this novel field and demonstrate the underlying benefits through numerical applications on a text book example. Among the various optimization tasks we will discuss only those which are connected to a controlled manipulation of the entanglement which is in fact the key ingredient of the methods considered in the paper. The selected topics will be covered according to a series of lectures given on the topic "New wavefunction methods and entanglement optimizations in quantum chemistry" at the Workshop on Theoretical Chemistry, 18 - 21 February 2014, Mariapfarr, Austria.

preprint2013arXiv

Excited states in poly-diacetylene chains: A Density-matrix-renormalization-group study

We study theoretically poly-diacetylene chains diluted in their monomer matrix. We employ the density-matrix renormalization group method (DMRG) on finite chains to calculate the ground state and low-lying excitations of the corresponding Peierls--Hubbard-Ohno Hamiltonian which is characterized by the electron transfer amplitude t0 between nearest neighbors, by the electron-phonon coupling constant α, by the Hubbard interaction U, and by the long-range interaction V. We treat the lattice relaxation in the adiabatic limit, i.e., we calculate the polaronic lattice distortions for each excited state. Using chains with up to 102 lattice sites, we can safely perform the extrapolation to the thermodynamic limit for the ground-state energy and conformation, the single-particle gap, and the energies of the singlet exciton, the triplet ground state, and the optical excitation of the triplet ground state. The corresponding gaps are known with high precision from experiment. We determine a coherent parameter set (t0*=2.4 eV, α*=3.4 eV/Å, U*=6 eV, V*=3 eV) from a fit of the experimental gap energies to the theoretical values which we obtain for 81 parameter points in the four dimensional search space (t0, α, U, V). We identify dark in-gap states in the singlet and triplet sectors as seen in experiment. Using a fairly stiff spring constant, the length of our unit cell is about one percent larger than its experimental value.

preprint2013arXiv

The density matrix renormalization group algorithm on kilo-processor architectures: implementation and trade-offs

In the numerical analysis of strongly correlated quantum lattice models one of the leading algorithms developed to balance the size of the effective Hilbert space and the accuracy of the simulation is the density matrix renormalization group (DMRG) algorithm, in which the run-time is dominated by the iterative diagonalization of the Hamilton operator. As the most time-dominant step of the diagonalization can be expressed as a list of dense matrix operations, the DMRG is an appealing candidate to fully utilize the computing power residing in novel kilo-processor architectures. In the paper a smart hybrid CPU-GPU implementation is presented, which exploits the power of both CPU and GPU and tolerates problems exceeding the GPU memory size. Furthermore, a new CUDA kernel has been designed for asymmetric matrix-vector multiplication to accelerate the rest of the diagonalization. Besides the evaluation of the GPU implementation, the practical limits of an FPGA implementation are also discussed.

preprint2012arXiv

Accurate ab initio spin densities

We present an approach for the calculation of spin density distributions for molecules that require very large active spaces for a qualitatively correct description of their electronic structure. Our approach is based on the density-matrix renormalization group (DMRG) algorithm to calculate the spin density matrix elements as basic quantity for the spatially resolved spin density distribution. The spin density matrix elements are directly determined from the second-quantized elementary operators optimized by the DMRG algorithm. As an analytic convergence criterion for the spin density distribution, we employ our recently developed sampling-reconstruction scheme [J. Chem. Phys. 2011, 134, 224101] to build an accurate complete-active-space configuration-interaction (CASCI) wave function from the optimized matrix product states. The spin density matrix elements can then also be determined as an expectation value employing the reconstructed wave function expansion. Furthermore, the explicit reconstruction of a CASCI-type wave function provides insights into chemically interesting features of the molecule under study such as the distribution of $α$- and $β$-electrons in terms of Slater determinants, CI coefficients, and natural orbitals. The methodology is applied to an iron nitrosyl complex which we have identified as a challenging system for standard approaches [J. Chem. Theory Comput. 2011, 7, 2740].

preprint2012arXiv

Entanglement Measures for Single- and Multi-Reference Correlation Effects

Electron correlation effects are essential for an accurate ab initio description of molecules. A quantitative a priori knowledge of the single- or multi-reference nature of electronic structures as well as of the dominant contributions to the correlation energy can facilitate the decision regarding the optimum quantum chemical method of choice. We propose concepts from quantum information theory as orbital entanglement measures that allow us to evaluate the single- and multi-reference character of any molecular structure in a given orbital basis set. By studying these measures we can detect possible artifacts of small active spaces.

preprint2010arXiv

Simulating Strongly Correlated Quantum Systems with Tree Tensor Networks

We present a tree-tensor-network-based method to study strongly correlated systems with nonlocal interactions in higher dimensions. Although the momentum-space and quantum-chemistry versions of the density matrix renormalization group (DMRG) method have long been applied to such systems, the spatial topology of DMRG-based methods allows efficient optimizations to be carried out with respect to one spatial dimension only. Extending the matrix-product-state picture, we formulate a more general approach by allowing the local sites to be coupled to more than two neighboring auxiliary subspaces. Following Shi. et. al. [Phys. Rev. A, 74, 022320 (2006)], we treat a tree-like network ansatz with arbitrary coordination number z, where the z=2 case corresponds to the one-dimensional scheme. For this ansatz, the long-range correlation deviates from the mean-field value polynomially with distance, in contrast to the matrix-product ansatz, which deviates exponentially. The computational cost of the tree-tensor-network method is significantly smaller than that of previous DMRG-based attempts, which renormalize several blocks into a single block. In addition, we investigate the effect of unitary transformations on the local basis states and present a method for optimizing such transformations. For the 1-d interacting spinless fermion model, the optimized transformation interpolates smoothly between real space and momentum space. Calculations carried out on small quantum chemical systems support our approach.

preprint2009arXiv

Density-matrix-renormalization-group study of excitons in poly-diacetylene chains

We study the elementary excitations of a model Hamiltonian for the $π$-electrons in poly-diacetylene chains. In these materials, the bare band gap is only half the size of the observed single-particle gap and the binding energy of the exciton of 0.5 eV amounts to 20% of the single-particle gap. Therefore, exchange and correlations due to the long-range Coulomb interaction require a numerically exact treatment which we carry out using the density-matrix renormalization group (DMRG) method. Employing both the Hubbard--Ohno potential and the screened potential in one dimension, we reproduce the experimental results for the binding energy of the singlet exciton and its polarizability. Our results indicate that there are optically dark states below the singlet exciton, in agreement with experiment. In addition, we find a weakly bound second exciton with a binding energy of 0.1 eV. The energies in the triplet sector do not match the experimental data quantitatively, probably because we do not include polaronic relaxation effects.