Reach Us
+44-7482-864460

**Michael Obende ^{1*}, Okanigbuan OR^{2} and Kolebaje Olusola^{3}**

^{1}Department of Physics, AfeBabalola University, Ado Ekiti, Nigeria

^{2}Department of Physics, Ambrose Alli University Ekpoma, Edo State, Nigeria

^{3}Department of Physics, Adeyemi College of Education, Ondo, Ondo State, Nigeria

- *Corresponding Author:
- Michael Obende

Department of Physics, AfeBabalola University, Ado Ekiti, Nigeria

**E-mail:**[email protected]

**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 [1].

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 [2].

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.

(2.1.0)

Which satisfies

Now we can construct a new state that is orthogonal to the previous two as:

(2.1.1)

Where n = 0,1, 2,− − − − − − and the co − efficient are given by :

(2.1.2)

In this basis, it can be shown that the Hamiltonian matrix becomes

(2.1.3)

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 [16].

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 [17], 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

Recall equation1.1

(2.2.1)

For two electrons on two sites

(2.3.2)

Similary,

(2.3.3)

The procedure can be generalized by defining an orthogonal basis recursively as:

(2.3.4)

Where n = 0,1, 2,− − − − and the coefficients are given by

Supplemented by in this basis. We have the Hamiltonian matrix:

(2.3.5)

From equation 1.4

a is the normalizing factor

Also recall that,

In this basis, we generate the Hamiltonian matrix

(2.3.6)

Similarly,

In the basis

(2.3.7)

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.

Either:

(2.3.8)

Hence equation (2.3.9) is the ground state energy.

(2.3.9)

The Corresponding were function is given by:

Where:

(1)

(2)

(3)

From equation (i)

From equation (3)

Also,

**Ground state energy for two electrons on two sites (Exact Method)**

**Recall that:**

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:

Here,

(2.4.0)

For

Which gives

We now determine the eigenvalue, Where A = (Hij) is a square matrix, X is a column matrix and λ is a scalar quantity [18-22].

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.

U/4T | Lanczos | Exact |
---|---|---|

4 | -0.24 | -0.24 |

3.5 | -0.28 | -0.28 |

2.5 | -0.38 | -0.38 |

2 | -0.4 | -0.4 |

1.5 | -0.6 | -0.6 |

0.75 | -1 | -1 |

0.5 | -1.23 | -1.23 |

0.25 | -1.56 | -1.56 |

0.1 | -1.8 | -1.8 |

0.02 | -1.96 | -1.96 |

0.01 | -1.98 | -1.98 |

0 | -2 | -2 |

-0.01 | -2.02 | -2.02 |

-0.02 | -2.04 | -2.04 |

-0.1 | -2.2 | -2.2 |

-0.25 | -2.56 | -2.56 |

-1 | -4.8 | -4.8 |

-1.5 | -6.6 | -6.6 |

-2 | -8.4 | -8.4 |

-2.5 | -10.4 | -10.4 |

-3 | -12.32 | -12.32 |

-4 | -16.24 | -16.24 |

**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.

- Adrian cho “Second family of high temperature superconductors Discovered. science Now Daily New.
- Bardeen J,Cooper LN Schrieffer JR (1957) “Theory of superconductivity. Phys. Rev108: 1175-1205.
- Charles K (2004) Introduction to solid state Physics. New York.
- Chen L, Mei C (1989) Exact calculation of the two electron interaction in the ground state of the Hubbard. Model. Phys. Rev. B39:9006.
- Dagotto E (1994) correlated electrons in high temperature superconductors Rev. mod. Phys66:763-840.
- Dagotto E, Moreo A(1984) Hamiltonian variational study of SU(2) lattice gauge theory Phys. Rev. D 29:2350.
- Dagotto E, Adriana Moreo (1985) Exact Diagonalization of Small System Phys. Rev. D1:865.
- Dagotto E(1994) Correlated Electrons in High-Temperature superconductors. Reviews of Modern Physics 66:763.
- Doniach S, Sondhermer EH (1974) Green‘s Function for solid state physicists: 222-234.
- Enaible AE,Idiodi JOA (2000)The two-electron problem studied within the t-j model and the Hubbard model. Proceedings of Nig. Ass Math Phys 231-242.
- Enabible AE,Idiodi JOA(2003)The two electron interaction in the ground state of theHubbard-Hirsch Hamiltonian J. Nig Assoc. Math Phys275-280.
- Feiner LF, Jefferson JH, Raimonddi R (1996)Effective single-band models for the high-Tccuprates Coulomb interactions. Phys, Rev B53:8751.
- GinzburgVI, Landau ID (1950)“Microscopic derivation of the Ginzbury-Landau equations in the theory of superconductivity”.
- Hirsch JE (1989) Metallic ferromagnetism in a single-band model Phys. Rev40:2354.
- Hubbard, (1963).Electron Correlations in Narrow Energy Bands J. Proc. Roy. Soc. London.
- Hybertsen MS, Stechel EB, Schluter M JennisonsDR (1990) Phys. Rev. B41, P11068.
- Lanczos C(1950) J. Res. Nat. Bur. Strand45:255.
- Louis Eugene Felix (2001) Introduction of superconductivity.
- Pettifor DG,weaire DL (1985) The recursion method and its Applications, Springler series in solid-state sciences.
- Richard MB (1951) Ferromagnetis.
- Roomany, Wyld H (1980) Structural Elements in particle physics and statistical mechanics .Phys. Rev. D21:3341.
- Tinkhams M (2004) Introduction to Superconductivity.

Select your language of interest to view the total content in your interested language

- Aerodynamic Flow Control
- Carbon Nanotubes
- Computational Fluid Dynamics
- Cosmoparticle Physics
- Detectors and Optical Sensors
- Fusion Plasmas
- Luminescence
- Nuclear Physics
- Parallel Programming
- Quantum Gravity
- Soil Physics
- Solid State Plasmas
- Space Plasmas
- Stellar Spectroscopy
- Suprnova Remnants
- Surface Physics

- World Congress on Physics

MARCH 20-21, 2020 Paris, France - 13th International Conference on Optics, Photonics & Laser

April 22-23, 2020 Berlin, Germany - 4th International Conference on Astronomy and Space Technology

May 20-21, 2020 Berlin, Germany

- Total views:
**13110** - [From(publication date):

June-2015 - Jan 23, 2020] - Breakdown by view type
- HTML page views :
**9289** - PDF downloads :
**3821**

**Make the best use of Scientific Research and information from our 700 + peer reviewed, Open Access Journals**

International Conferences 2020-21