Abstract:With the development of computing technology, numerical exact diagonalization method plays a vital role in modern computational condensed matter physics, especially in the research area of strongly correlated electron systems: it becomes a benchmark for other numerical computational techniques, such as quantum Monte Carlo, numerical renormalization group, density matrix renormalization group, and dynamic mean field theory. In this paper, we first numerically exactly diagonalize the three-dimensional hydrogen atom with the combination of finite-difference method, and find that the numerical wave function of ground state is in good agreement with the analytical calculations. We then turn to discuss the space dimension confinement hydrogen system, two-dimensional hydrogen atom, and notice that the numerical wave function is no longer in agreement with the analytical calculation, where the ground state wave function has a numerical singularity as radius approaches to zero. Compared with the case of the three-dimensional hydrogen atom, this issue mainly comes from the nature of space dimension confinement. To resolve such an issue of numerical singularity in two-dimensional hydrogen atom, we need to construct a new discrete and normalized Bessel function as a basis to study the ground state behavior of dimension confinement system based on the framework of Lanczos-type numerical exact diagonalization. The constructed normalized Bessel basis is orthogonal and discrete, and thus becomes suitable for practical calculation. Besides, these prominent properties of such a Bessel basis greatly reduce the complexity and difficulty in practical calculation, and thus makes computing work efficient. In addition, Lanczos-type numerical exact diagonalization method can extremely speed up the process of solving the eigenvalue equation. As a result, such a high efficient calculation of our method demonstrates the consistence between numerical and analytical ground state energy value, and the corresponding wave function with enough truncated basis number. Since this kind of numerical singularity occurs in many space dimension confinement systems, our finding for constructing a new discrete Bessel basis function may be helpful in studying the quantum systems with numerical singularity behaviors in wavefunctions in future. On the other hand, it should be pointed out that the Bessel basis is incorporated into the linear augment plane wave method in the density functional theory to study the electronic band structure of the condensed material and obtain high accurate results, especially in the theoretical prediction of topological insulators and in experimental realization as well.

High-accuracy solution of the Dirac equation for three-dimensional hydrogen atom using a variationally optimized diagonalization method

High-precision solution of the Dirac Equation for the hydrogen molecular ion using a basis-set expansion

Self-Consistent Solution Of Kohn-Sham Equation By Real-Space Finite-Difference Method

Calculating the many-potential vacuum polarization density of the Dirac equation in the finite-basis approximation

Numerical Exact Diagonalization of Singularity in the Ground State of Two-Dimensional Hydrogen Atom

Solving Dirac Equations on A 3d Lattice with Inverse Hamiltonian and Spectral Methods

Efficient Solution for the Dirac Equation in 3D Lattice Space with the Conjugate Gradient Method

4-component Relativistic Calculations in a Multiwavelet Basis with Improved Convergence

Self-consistent Density Functional Calculations without Diagonalization

All-electron fully relativistic Kohn-Sham theory for solids based on the Dirac-Coulomb Hamiltonian and Gaussian-type functions

Solving the relativistic Hartree-Bogoliubov equation with the finite-difference method

Optimized Dirac Woods-Saxon basis for covariant density functional theory

The non-relativistic expansion of Dirac-Coulomb energy and the non-retarded Breit interaction correction up to $α^8$order

An accurate Fortran code for computing hydrogenic continuum wave functions at a wide range of parameters

Accurate and Efficient Calculations of Hellmann-Feynman Forces for Quantum Computation.

Analytical Structure Matching and Very Precise Approach to the Coulombic Quantum Three-Body Problem

Gaussian Beam Methods for the Dirac Equation in the Semi-classical Regime

Electronic structure calculations with interpolating tensor product wavelet basis

Solution of the Dirac–Coulomb equation using the Rayleigh–Ritz method: Results for He-like atoms

High-precision evaluation of Wigner's d-matrix by exact diagonalization

Efficient iterative method for solving the Dirac–Kohn–Sham density functional theory