Mapped discrete variable representations using the sine and cosine polynomials are purposed to accurately solve the stationary electronic Shrödinger equation of molecules. In this method, the grid points densely cluster around the positions of the nuclei but sparse enough in other regions. The electronic states of H2+ and H2 molecules are calculated by a preconditioned inexact spectral transform method in cylindrical coordinates for illustrating. A spectral convergence behaviour of eigenvalues of the electronic state of the H2+ molecule but slower convergence behaviour for the H2 molecule are observed, due to the repulsive singularity between electrons in the H2 molecule.
Read full abstract