Symmetry-adapted linear combinations of valence-bond (VB) diagrams are constructed for arbitrary point groups and total spin S using diagrammatic VB methods. VB diagrams are related uniquely to invariant subspaces whose size reflects the number of group elements; their nonorthogonality leads to sparser matrices and is fully incorporated into a binary integer representation. Symmetry-adapated linear combinations of VB diagrams are constructed for the 1764 singlets of a half-filled cube of eight sites, the 2.8 million π-electron singlets of anthracene, and for illustrative S≥0 systems.
Symmetry-adapted linear combinations of valence-bond (VB) diagrams are constructed for arbitrary point groups and total spin S using diagrammatic VB methods. VB diagrams are related uniquely to invariant subspaces whose size reflects the number of group elements; their nonorthogonality leads to sparser matrices and is fully incorporated into a binary integer representation. Symmetry-adapated linear combinations of VB diagrams are constructed for the 1764 singlets of a half-filled cube of eight sites, the 2.8 million π-electron singlets of anthracene, and for illustrative S≥0 systems.
Abstract Valence bond (VB) diagrams form a complete basis for model Hamiltonians that conserve total spin, S, and have one valence state, ϕp, per site. Hubbard and Pariser–Parr–Pople (PPP) models illustrate ionic problems, with zero, one, or two electrons in each ϕp, while isotropic Heisenberg models illustrate spin problems, with only purely covalent VB diagrams. The difficulty of nonorthogonal VB diagrams is by‐passed by exploiting the finite dimensionality of the complete basis and working with unsymmetric sparse matrices. We introduce efficient bit manipulations for generating, storing, and handling VB diagrams as integers and describe a new coordinate relaxation method for the ground and lowest excited states of unsymmetric sparse matrices. Antiferromagnetic spin‐½ Heisenberg rings and chains of N ⩽ 20 spins, or 2N spin functions, are solved in C2 symmetry as illustrative examples. The lowest S = 1 and 0 excitations are related to domain walls, or spin solitons, and studied for alternations corresponding to polyacetylene. VB diagrams with arbitrary S and nonneighbor interactions are constructed for both spin and ionic problems, thus extending diagrammatic VB theory to other topologies.
One of the long standing problems in quantum chemistry had been the inability to exploit full spatial and spin symmetry of an electronic Hamiltonian belonging to a non-Abelian point group. Here, we present a general technique which can utilize all the symmetries of an electronic (magnetic) Hamiltonian to obtain its full eigenvalue spectrum. This is a hybrid method based on Valence Bond basis and the basis of constant z-component of the total spin. This technique is applicable to systems with any point group symmetry and is easy to implement on a computer. We illustrate the power of the method by applying it to a model icosahedral half-filled electronic system. This model spans a huge Hilbert space (dimension 1,778,966) and in the largest non-Abelian point group. The C60 molecule has this symmetry and hence our calculation throw light on the higher energy excited states of the bucky ball. This method can also be utilized to study finite temperature properties of strongly correlated systems within an exact diagonalization approach. (C) 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012
A systematic procedure is described for calculating the reduced density matrices of first and second order for a resonating valence bond wavefunction on square planar nets with Born-von Karman boundary conditions of finite period. The density matrices explicitly display the electron correlation in such a wavefunction, which is a model ground state for a two-dimensional Mott insulator. Because of the extremely large number of contributing valence bond structures, it is essential to utilize their symmetry under the space group for minimum computation. For the totally symmetric singlet state described, density matrix eigenvectors form basis vectors for irreducible representations. The irreducible representations spanned by pair states on the square plane lattice are listed, and a set of irreducible basis vectors for the two-particle density matrix is defined which clearly displays the nature of electron pair correlations in such a system (correlated geminals). The apparatus has been employed in a calculation described in the paper accompanying this one (succeeding article).
One-electron wave functions of half-filled-band closed-shell systems have been decomposed into their orthogonal single-determinant valence-bond (VB) components (antisymmetrized products of orthogonal atomic orbitals). Despite the ordering effect of the Fermi hole among electrons of the same spin, favoring spin-density and charge-density waves, the relative independence between the \ensuremath{\alpha} and \ensuremath{\beta} parts of the wave function results in an overestimation of the weight of multiply ionic structures. Comparison is made with the internally correlated functions for different graphs and Hamiltonians (Hubbard, Pariser-Parr-Pople, and ab initio types). The scaling effect of the correlation does not appear as essentially dependent on the ionicity of the VB forms but exponentially linked to the mean energy of each VB component.