Finite element calculations for systems with multiple Coulomb centers

• Autores: M. Braun
• Localización: Journal of computational and applied mathematics, ISSN 0377-0427, Vol. 236, Nº 18, 2012 (Ejemplar dedicado a: FEMTEC 2011: 3rd International Conference on Computational Methods in Engineering and Science, May 9-13, 2011), págs. 4840-4846
• Idioma: inglés
• DOI: 10.1016/j.cam.2012.02.005
• Texto completo no disponible (Saber más ...)
• Resumen

  • The presence of multiple Coulomb centers in molecules or solids poses a challenge when solving the effective Schrödinger equation, required as a crucial ingredient in density functional or Hartree�Fock calculations. This is primarily because Kato�s cusp condition needs to be satisfied close to each nucleus and the matrix elements of the Coulomb potential at the nuclei are rather difficult to evaluate when using global basis functions.

    A novel method for dealing with these challenges is introduced, rewriting the wavefunction as a product of a function satisfying the nuclear cusp conditions and a smooth function, resulting in a transformed variational principle and a regularized potential. Results of three-dimensional finite element calculations based on this ansatz for the ground state of the molecule H+ 2 in the Born�Oppenheimer approximation are presented, which were obtained using custom written Python/Fortran code.