Received Date: March 24, 2015; Accepted Date: April 29, 2015; Published Date: May 20, 2015
Citation: Obende M, Okanigbuan OR, Olusola K (2015) Correlation between Lanczos Method and Exact Diagonalisation Method in the Study of Highly Correlated Electrons System. J Astrophys Aerospace Technol 3: 113. doi:10.4172/2329-6542.1000113
Copyright: © 2015 Obende M, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Visit for more related articles at Journal of Astrophysics & Aerospace Technology
The ground state energy and wave function of the single band hubbard model on the one dimensional lattice is computed using the Lanczos methods. It is shown that the ground state energies obtained for different values of the interaction strength it, compare nicely with that obtained using the method of exact calculation.
Hubbard model; Lanczos method; Exact method
The Hubbard model is an approximate model used to describe the transition between conducting and insulating systems. The Hubbard model, named after John Hubbard 1963, is the simplest model of interacting particles in lattices, with the only two tennis in the Hamiltonian; a kinetic term allowing for hopping of particles between soles of the lattice and a potential term consisting of an on sites interaction. The particles can either be fermions or bosons .
The Hubbard model is a good approximate for particles in a periodic potential at sufficiently low temperature that all the particles are in the lowest Bloch band, as long as any long-range interactions between the particles can be ignored. If interaction between particles on different sites of the lattice are included, the model is often referred to as the extended Hubbard model .
The model was originally proposed to describe electrons in solids and has since been the focus of particular interest as a model for high temperature super conductivity. For electron in a solid, the Hubbard model can be considered as an improvement on the tight-binding model, which includes only the hopping term for strong interaction, it can give qualitatively different behavior from the tight. Binding model, and correctly predicts the existence of so called mott insulators, which are prevented from becoming conducting by the strong repulsion between the particles [3,4].
The Hubbard model is based on the tight –binding approximation.
In the tight-binding approximation, electrons are viewed as occupying the standard orbital of their hopping between atoms, and then hopping between atoms, and then hopping between atoms during conduction. Mathematically, this is represented as hopping integral or transfer integral between neighboring atoms, which can be viewed as the physical principle that creates electrons bonds in crystalline materials [5-7].
The Hubbard Hamiltonian takes the form:
Where are the creation, annihilation, and number operations, respectively, for an electron of spin σ in the wannier state on the on the ith lattice Fi, j means that only interest neighbor site hopping are allowed [8,9].
The Hubbard model has been extensively studied the single band Hubbard model in both one dimensional ring and two dimensional torus. By the exact calculation of the pair correlation function Fi, j defined as the possibility of finding an electron at site when there has been an electron of opposite at site I, it is shown that for two electrons, the interaction is always repulsive in the ground state for any positive value of the on–site coulomb interaction u .
Enable and Idiodi studied the single –band tight-binding model with on-site repulsion u and nearest neighbor-exchange interaction j (the so called Hubbard Hirsch Hamiltonian) with the help of a correlated variational approach [10,11]. Two finite sized lattices with periodic boundary conditions were considered and the criteria for the occurrence of a transition from an anti-ferromagnetic phase to a ferromagnetic phase were discussed.
The Lanczos Tridiagonalization Method
The basic idea of the Lanczos method is that a special basis can be constructed where the Hamiltonian has a tridiagonal representation. This is carried out iteratively as shown below. First, it is necessary to select an arbitrary vector in the Hilbert space of the model being studied. If the Lanczos method is used to obtain the ground state energy of the model, then it is necessary that the overlap between the actual ground stat and the initial state be non-zero [12,13]. If there is no priori information about the ground state is know, this requirement is usually satisfied by selecting an initial state with randomly satisfied by selecting an initial state randomly chosen co-efficient in the working basis that is being used. If some other information about the ground state is known, like it is convenient to initiate the iterations with a state already belonging to the subspace having those quantum numbers (and still with random coefficient within this subspace) [14,15].
After is selected, we can define a new vector by applying Hamiltonian H to the initial state subtraction the projection over , we obtain.
Now we can construct a new state that is orthogonal to the previous two as:
Where n = 0,1, 2,− − − − − − and the co − efficient are given by :
In this basis, it can be shown that the Hamiltonian matrix becomes
i.e it is tridiagonal, as expected. Once in this form the matrix can be diagonalized easily using standard library subroutine. However, to diagonalise completely the model being studied on a smote cluster, a number of iterations to the size of the Hilbert space (or of the subspace under consideration) are needed .
However, one of the advantages of this technique is that accurate enough information about the ground state of the problem can be obtained after a small number of iterations /typically of the order of ≈100 or less). Thus the method is suitable for the analysis of low temperature properties of the models of correlated electrons described.
Ground state energy of the two electrons on 2 sites of aid lattice using the lanczos method
In the Lanczos method, a special basis is constructed such that the Hamiltonian has a tridigonal representation , Firstly, we select an arbitrary vector in the Hilbert space of the model being studied for a two site system containing 2 electrons, we have the following many body states;
Let We now define a new vector by applying the Hubbard Hamiltonian to the state subtracting the projection over we obtain
For two electrons on two sites
The procedure can be generalized by defining an orthogonal basis recursively as:
Where n = 0,1, 2,− − − − and the coefficients are given by
Supplemented by in this basis. We have the Hamiltonian matrix:
From equation 1.4
a is the normalizing factor
Also recall that,
In this basis, we generate the Hamiltonian matrix
In the basis
Diagonalising equation (2.3.8), we solve the equation.
By expansions, we have that:
Hence, by using polynomial division, we have that:
To determine the eigenvalue (x), we solve the above equation.
Hence equation (2.3.9) is the ground state energy.
The Corresponding were function is given by:
From equation (i)
From equation (3)
Ground state energy for two electrons on two sites (Exact Method)
For n electron on N sites, the number of many body states is given by
It has six states which are
With these Hamiltonian set up:
Ground state energy for two electron system on two sites were analytically solved using 4× 4 matrix and the results obtained are presented in the figures and table below.
Figure 1 shows ground state energy as a function of positive U energy for 2 electrons on 2 sites for Lanczos method while Figure 2 shows ground state energy as a function of positive U for 2 electrons on 2 sites for exact method. Thus, the correlation between the two graphs are in perfect agreement.
Table 1 below shows the numerical results of Lanczos method and exact method which is in a good agreement with experimental values.
Table 1: Numerical Result Correlation Between Lanczos and Exact Method.
In this study, we have shown that Lanczos technique is a reliable numerical method in determining the ground state energy of a system described by the Hubbard Hamiltonian. Ground state energies obtained using the Lanczos method compare nicely with that obtained using exact method. The size of the Hamiltonian matrix to diagonalise is reduced from (4×4) matrix to a (3×3) matrix when Lanczos technique is applied.
Make the best use of Scientific Research and information from our 700 + peer reviewed, Open Access Journals