Exact diagonalization

Exact diagonalization (ED) is a numerical technique used in physics to determine the eigenstates and energy eigenvalues of a quantum Hamiltonian. In this technique, a Hamiltonian for a discrete, finite system is expressed in matrix form and diagonalized using a computer. Exact diagonaliation is only feasible for systems with a few tens of particles, due to the exponential growth of the Hilbert space dimension with the size of the quantum system. It is frequently employed to study lattice models, including the Hubbard model, Ising model, Heisenberg model, t-J model, and SYK model.[1][2]

Expectation values from exact diagonalization

After determining the eigenstates and energies of a given Hamiltonian, exact diagonalization can be used to obtain expectation values of observables. For example, if is an observable, its thermal expectation value is

where is the partition function. If the observable can be written down in the initial basis for the problem, then this sum can be evaluated after transforming to the basis of eigenstates.

Green's functions may be evaluated similarly. For example, the retarded Green's function can be written

Exact diagonalization can also be used to determine the time evolution of a system after a quench. Suppose the system has been prepared in an initial state , and then for time evolves under a new Hamiltonian, . The state at time is

Memory requirements

The dimension of the Hilbert space describing a quantum system scales exponentially with system size. For example, consider a system of spins localized on fixed lattice sites. The dimension of the on-site basis is 2, because the state of each spin can be described as a superposition of spin-up and spin-down, denoted and . The full system has dimension , and the Hamiltonian represented as a matrix has size . This implies that computation time and memory requirements scale very unfavorably in exact diagonalization. In practice, the memory requirements can be reduced by taking advantage of symmetry of the problem, imposing conservation laws, working with sparse matrices, or using other techniques.

Number of sites Number of states Hamiltonian size in memory
4 16 2048 B
9 512 2 MB
16 65536 34 GB
25 33554432 9 PB
36 6.872e10 40 ZB
Naive estimates for memory requirements in exact diagonalization of a spin-½ system performed on a computer. It is assumed the Hamiltonian is stored as a matrix of double-precision floating point numbers.

Comparison with other techniques

Exact diagonalization is useful for extracting exact information about finite systems. However, often small systems are studied to gain insight into infinite lattice systems. If the diagonalized system is too small, its properties will not reflect the properties of the system in the thermodynamic limit, and the simulation is said to suffer from finite size effects.

Unlike some other exact theory techniques, such as Auxiliary-field Monte Carlo, exact diagonalization obtains Green's functions directly in real time, as opposed to imaginary time. Unlike in these other techniques, exact diagonalization results do not need to be numerically analytically continued. This is an advantage, because numerical analytic continuation is an ill-posed and difficult optimization problem.[3]

Applications

  • Studying various properties of the 2D Heisenberg model in a magnetic field, including antiferromagnetism and spin-wave velocity.[6]
  • Studying the Drude weight of the 2D Hubbard model.[7]
  • Studying out-of-time-order correlations (OTOCs) and scrambling in the SYK model.[8]
  • Simulating resonant x-ray spectra of strongly correlated materials.[9]

Implementations

Numerous software packages implementing exact diagonalization of quantum Hamiltonians exist. These include QuSpin, ALPS, DoQo, EdLib, edrixs, and many others.

Generalizations

Exact diagonalization results from many small clusters can be combined to obtain more accurate information about systems in the thermodynamic limit using the numerical linked cluster expansion.[10]

See also

References

  1. Bergeron, Dominic; Tremblay, A.-M. S. (5 August 2016). "Algorithms for optimized maximum entropy and diagnostic tools for analytic continuation". Physical Review E. 94 (2). arXiv:1507.01012. doi:10.1103/PhysRevE.94.023303.
  2. Medvedeva, Darya; Iskakov, Sergei; Krien, Friedrich; Mazurenko, Vladimir V.; Lichtenstein, Alexander I. (29 December 2017). "Exact diagonalization solver for extended dynamical mean-field theory". Physical Review B. 96 (23). arXiv:1709.09176. doi:10.1103/PhysRevB.96.235149.
  3. Hamer, C. J.; Barber, M. N. (1 January 1981). "Finite-lattice methods in quantum Hamiltonian field theory. I. The Ising model". Journal of Physics A: Mathematical and General. 14 (1): 241–257. doi:10.1088/0305-4470/14/1/024.
  4. Lüscher, Andreas; Läuchli, Andreas M. (5 May 2009). "Exact diagonalization study of the antiferromagnetic spin-1/2 Heisenberg model on the square lattice in a magnetic field". Physical Review B. 79 (19). arXiv:0812.3420. doi:10.1103/PhysRevB.79.195102.
  5. Nakano, Hiroki; Takahashi, Yoshinori; Imada, Masatoshi (15 March 2007). "Drude Weight of the Two-Dimensional Hubbard Model –Reexamination of Finite-Size Effect in Exact Diagonalization Study–". Journal of the Physical Society of Japan. 76 (3): 034705. arXiv:cond-mat/0701735. doi:10.1143/JPSJ.76.034705.
  6. Fu, Wenbo; Sachdev, Subir (15 July 2016). "Numerical study of fermion and boson models with infinite-range random interactions". Physical Review B. 94 (3). arXiv:1603.05246. doi:10.1103/PhysRevB.94.035135.
  7. Wang, Y.; Fabbris, G.; Dean, M.P.M; Kotliar, G. (2019). EDRIXS: An open source toolkit for simulating spectra of resonant inelastic x-ray scattering. 243. Computer Physics Communications. pp. 151–165. arXiv:1812.05735. doi:10.1016/j.cpc.2019.04.018.
  8. Tang, Baoming; Khatami, Ehsan; Rigol, Marcos (March 2013). "A short introduction to numerical linked-cluster expansions". Computer Physics Communications. 184 (3): 557–564. arXiv:1207.3366. doi:10.1016/j.cpc.2012.10.008.
This article is issued from Wikipedia. The text is licensed under Creative Commons - Attribution - Sharealike. Additional terms may apply for the media files.