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
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
This thesis consists of six chapters. The first chapter gives an introduction to the field of low-dimensional magnetic and electronic systems and relevant numerical techniques. The recent developments in molecular magnets are highlighted. The numerical techniques are reviewed along with their advantages and disadvantages from the present perspective. Study of entanglement of a system can give a great insight into the system. At the last part of this chapter a general overview is given regarding entanglement, its measures and its significance in studying many-body systems.
Chapter 2 deals with the technique that has been developed by us for the full symmetry adaptation of non-relativistic Hamiltonians. It is advantageous both computationally and physically/chemically to exploit both spin and spatial symmetries of a system. It has been a long-standing problem to target a state which has definite total spin and also belongs to a definite irreducible representation of a point group, particularly for non-Abelian point groups. A very general technique is discussed in this chapter which is a hybrid method based on valence-bond basis and the basis of the z-component of the total spin. This technique is not only applicable to a system with arbitrary site spins and belonging to any point group symmetry, it is also quite easy to implement computationally. To demonstrate the power of the method, it is applied to the molecular magnetic system, Cu6Fe8, with cubic symmetry.
In chapter 3, the extension of the previous hybrid technique to electronic systems is discussed. The power of the method is illustrated by applying it to a model icosahedral half-filled electronic system. This model spans a huge Hilbert space (dimension 1,778,966) and is in the largest non-Abelian point group. All the eigenstates of the model are obtained using our technique.
Chapter 4 deals with the thermodynamic properties of an important class of single-chain magnets (SCMs). This class of SCMs has alternate isotropic spin-1/2 units and anisotropic high spin units with the anisotropy axes being non-collinear. Here anisotropy is assumed to be large and negative, as a result, anisotropic units behave like canted spins at low temperatures; but even then simple Ising-type model does not capture the essential physics of the system due to quantum mechanical nature of the isotropic units. A transfer matrix (TM) method is developed to study statistical behavior of this class of SCMs. For the first time, it is also discussed in detail that how weak inter-chain interactions can be treated by a TM method. The finite size effect is also discussed which becomes important for low temperature dynamics. This technique is applied to a real helical chain magnet, which has been studied experimentally.
In the fifth chapter a bipartite entanglement entropy of finite systems is studied using exact diagonalization techniques to examine how the entanglement changes in the presence of long-range interactions. The PariserParrPople model with long-range interactions is used for this…
We extend our previous work [S. Sharma and G. K.-L. Chan, J. Chem. Phys. 136, 124121 (2012)], which described a spin-adapted (SU(2) symmetry) density matrix renormalization group algorithm, to additionally utilize general non-Abelian point group symmetries. A key strength of the present formulation is that the requisite tensor operators are not hard-coded for each symmetry group, but are instead generated on the fly using the appropriate Clebsch-Gordan coefficients. This allows our single implementation to easily enable (or disable) any non-Abelian point group symmetry (including SU(2) spin symmetry). We use our implementation to compute the ground state potential energy curve of the C2 dimer in the cc-pVQZ basis set (with a frozen-core), corresponding to a Hilbert space dimension of 10(12) many-body states. While our calculated energy lies within the 0.3 mEh error bound of previous initiator full configuration interaction quantum Monte Carlo and correlation energy extrapolation by intrinsic scaling calculations, our estimated residual error is only 0.01 mEh, much more accurate than these previous estimates. Due to the additional efficiency afforded by the algorithm, the excitation energies (Te) of eight lowest lying excited states: a(3)Πu, b(3)Σg (-), A(1)Πu, c(3)Σu (+), B(1)Δg, B(') (1)Σg (+), d(3)Πg, and C(1)Πg are calculated, which agree with experimentally derived values to better than 0.06 eV. In addition, we also compute the potential energy curves of twelve states: the three lowest levels for each of the irreducible representations (1)Σg (+), (1)Σu (+), (1)Σg (-), and (1)Σu (-), to an estimated accuracy of 0.1 mEh of the exact result in this basis.
We present two methods to study model Hamiltonians defined on clusters. One is the symmetrization of the Hilbert space basis according to irreducible representations of the symmetry group of the Hamiltonian. The other one adapts the Lanczos diagonalization method directly to the second quantized form of the Hamiltonian studied. By the symmetrization of the basis, one takes advantage of the symmetries of the Hamiltonian, and uses the group theory to find the invariant subspaces of the Hamiltonian in the complete Hilbert space. Then one needs only to work within the invariant subspaces. This method cuts down substantially the size of the matrix representation of the Hamiltonian to be diagonalized. Therefore, one could study relatively larger systems. As a result, all the eigenvectors and eigenvalues are obtained. Consequently, one could calculate thermodynamical quantities and various correlation functions, such as, the spin correlation, the charge correlation, the superconductivity pairing correlation, and the photoemission spectra. One the other hand, the Lanczos method, mentioned here in particular, is a specially developed computer program for the study of the model Hamiltonians. This program uses the key feature of the Lanczos method in diagonalization, allows one to use the second quantized form of the Hamiltonian directly without pregeneration or storing the matrix elements of the Hamiltonian. Two methods could be combined to take the full advantage of the symmetries of the Hamiltonian by input of the symmetrized wavefunction as the starting trial wavefunction in the Lanczos method. Here we illustrate the above two methods by applying to the Hubbard, and extended Hubbard models on a simple cubic cluster, and the lattice Anderson model on the two dimensional square lattice. Other interesting results for different systems are also briefly discussed here.
The essence of quantum magnetism lies in correlated many-body states. They determine the physical properties beyond paramagnetism. Even constructions as simple as singlets or triplets cannot be understood without the mathematical concept of superpositions. This is even more evident if time evolution of quantum states and phenomena such as quantum coherenceQuantum coherence are considered. The major approach to assess correlated many-body states in quantum magnetism of molecules is exact diagonalizationDiagonalization of the many-body Hamiltonian describing the physical system. Symmetries are of great help when tackling the numerical task of matrix diagonalization since they allow to decompose a problem of large Hilbert spaceHilbert space dimension into a number of smaller problems. This approach will be presented first. It rests on the identification of the available symmetries—e.g., U(1), SU(2), as well as point groupsPoint group—and their application. However, the applicability of exact diagonalizationDiagonalization is limited by the exponential growth of the Hilbert spaceHilbert space of the spin system with the number of spins: $$\dim \left( {\mathcal{H}} \right) = \mathop \prod \nolimits_{i} \left( {2 s_{i} + 1} \right)$$ . Nowadays, complex Hermitian matricesHermitian matrices of a linear dimension of about 100,000 can be completely diagonalized on supercomputers which limits the number of spins depending on their size si to just a few. I will therefore present a very accurate approximation to exact diagonalization that rests on Lanczos procedures and the concept of typicality. The most popular version, the finite-temperature Lanczos methodFinite-Temperature Lanczos Method (FTLM), allows to evaluate low-energy spectroscopic data, thermal averages as well as time evolutions for systems with Hilbert space dimensions of up to about 1011. Hands-on experience and examples will be provided.
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.
While either a spin or point-group adaptation is straightforward when considered independently, the standard technique for factoring isotropic spin Hamiltonians by the total spin S and the irreducible representation Γ of the point group is limited by the complexity of the transformations between different coupling schemes that are related in terms of their site permutations. To overcome these challenges, we apply projection operators directly to uncoupled basis states, enabling the simultaneous treatment of spin and point-group symmetry without the need for recoupling transformations. This provides a simple and efficient approach for the exact diagonalization of isotropic spin models, which we illustrate, with applications in Heisenberg spin rings and polyhedra, including systems that are computationally inaccessible with conventional coupling techniques.