Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userOrganization of the Hilbert Space for Exact Diagonalization of Hubbard Model
142
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present an alternative scheme to the widely used method of representing the basis of one-band Hubbard model through the relation $I=I_{uparrow}+2^{M}I_{downarrow}$ given by H. Q. Lin and J. E. Gubernatis [Comput. Phys. 7, 400 (1993)], where $I_{uparrow}$, $I_{downarrow}$ and $I$ are the integer equivalents of binary representations of occupation patterns of spin up, spin down and both spin up and spin down electrons respectively, with $M$ being the number of sites. We compute and store only $I_{uparrow}$ or $I_{downarrow}$ at a time to generate the full Hamiltonian matrix. The non-diagonal part of the Hamiltonian matrix given as ${cal{I}}_{downarrow}otimes{bf{H}_{uparrow}} oplus {bf{H}_{downarrow}}otimes{cal{I}}_{uparrow}$ is generated using a bottom-up approach by computing the small matrices ${bf{H}_{uparrow}}$(spin up hopping Hamiltonian) and ${bf{H}_{downarrow}}$(spin down hopping Hamiltonian) and then forming the tensor product with respective identity matrices ${cal{I}}_{downarrow}$ and ${cal{I}}_{uparrow}$, thereby saving significant computation time and memory. We find that the total CPU time to generate the non-diagonal part of the Hamiltonian matrix using the new one spin configuration basis scheme is reduced by about an order of magnitude as compared to the two spin configuration basis scheme. The present scheme is shown to be inherently parallelizable. Its application to translationally invariant systems, computation of Greens functions and in impurity solver part of DMFT procedure is discussed and its extention to other models is also pointed out.
rate research
Read More
We have obtained the exact ground state wave functions of the Anderson-Hubbard model for different electron fillings on a 4x4 lattice with periodic boundary conditions - for 1/2 filling such ground states have roughly 166 million states. When compared to the uncorrelated ground states (Hubbard interaction set to zero) we have found strong evidence of the very effective screening of the charge homogeneities due to the Hubbard interaction. We have successfully modelled these local charge densities using a non-interacting model with a static screening of the impurity potentials. In addition, we have compared such wave functions to self-consistent real-space unrestricted Hartree-Fock solutions and have found that these approximate ground state wave functions are remarkably successful at reproducing the local charge densities, and may indicate the role of dipolar backflow in producing a novel metallic state in two dimensions.
We have used exact numerical diagonalization to study the excitation spectrum and the dynamic spin correlations in the $s=1/2$ next-next-nearest neighbor Heisenberg antiferromagnet on the square lattice, with additional 4-spin ring exchange from higher order terms in the Hubbard expansion. We have varied the ratio between Hubbard model parameters, $t/U$, to obtain different relative strengths of the exchange parameters, while keeping electrons localized. The Hubbard model parameters have been parametrized via an effective ring exchange coupling, $J_r$, which have been varied between 0$J$ and 1.5$J$. We find that ring exchange induces a quantum phase transition from the $(pi, pi)$ ordered Ne`el state to a $(pi/2, pi/2)$ ordered state. This quantum critical point is reduced by quantum fluctuations from its mean field value of $J_r/J = 2$ to a value of $sim 1.1$. At the quantum critical point, the dynamical correlation function shows a pseudo-continuum at $q$-values between the two competing ordering vectors.
We propose a distinct numerical approach to effectively solve the problem of partial diagonalization of the super-large-scale quantum electronic Hamiltonian matrices. The key ingredients of our scheme are the new method for arranging the basis vectors in the computers RAM and the algorithm allowing not to store a matrix in RAM, but to regenerate it on-the-fly during diagonalization procedure. This scheme was implemented in the program, solving the Anderson impurity model in the framework of dynamical mean-field theory (DMFT). The DMFT equations for electronic Hamiltonian with 18 effective orbitals that corresponds to the matrix with the dimension of 2.4 * 10^9 were solved on the distributed memory computational cluster.
We present an efficient exact diagonalization scheme for the extended dynamical mean-field theory and apply it to the extended Hubbard model on the square lattice with nonlocal charge-charge interactions. Our solver reproduces the phase diagram of this approximation with good accuracy. Details on the numerical treatment of the large Hilbert space of the auxiliary Holstein-Anderson impurity problem are provided. Benchmarks with a numerically exact strong-coupling continuous-time quantum-Monte Carlo solver show better convergence behavior of the exact diagonalization in the deep insulator. Special attention is given to possible effects due to the discretization of the bosonic bath. We discuss the quality of real axis spectra and address the question of screening in the Mott insulator within extended dynamical mean-field theory.
The tunable magnetism at graphene edges with lengths of up to 48 unit cells is analyzed by an exact diagonalization technique. For this we use a generalized interacting one-dimensional model which can be tuned continuously from a limit describing graphene zigzag edge states with a ferromagnetic phase, to a limit equivalent to a Hubbard chain, which does not allow ferromagnetism. This analysis sheds light onto the question why the edge states have a ferromagnetic ground state, while a usual one-dimensional metal does not. Essentially we find that there are two important features of edge states: (a) umklapp processes are completely forbidden for edge states; this allows a spin-polarized ground state. (b) the strong momentum dependence of the effective interaction vertex for edge states gives rise to a regime of partial spin-polarization and a second order phase transition between a standard paramagnetic Luttinger liquid and ferromagnetic Luttinger liquid.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
