New Techniques in Numerical Integration: The Computation of Molecular Integrals over Exponential-Type Functions

  • Author / Creator
    Slevinsky, Richard M
  • The numerical evaluation of challenging integrals is a topic of interest in applied mathematics. We investigate molecular integrals in the B function basis, an exponentially decaying basis with a compact analytical Fourier transform. The Fourier property allows analytical expressions for molecular integrals to be formulated in terms of semi-infinite highly oscillatory integrals with limited exponential decay. The semi-infinite integral representations in terms of nonphysical variables stand as the bottleneck in the calculation. To begin our numerical experiments, we conduct a comparative study of the most popular numerical steepest descent methods, extrapolation methods and sequence transformations for computing semi-infinite integrals. It concludes that having asymptotic series representations for integrals and applying sequence transformations leads to the most efficient algorithms. For three-center nuclear attraction integrals, we find an analytical expression for the semi-infinite integrals. Numerical experiments show the resulting algorithm is approximately 300 times more efficient than the state-of-the-art. For the four-center two-electron Coulomb integrals, we take a different approach. The integrand has singularities in the complex plane that can be near the path of integration, making standard quadrature routines unreliable. The trapezoidal rule with double exponential variable transformations has been shown to have very promising properties as a general-purpose integrator. We investigate the use of conformal maps to maximize the convergence rate, resulting in a nonlinear program for the optimized variable transformation.

  • Subjects / Keywords
  • Graduation date
  • Type of Item
  • Degree
    Doctor of Philosophy
  • DOI
  • License
    This thesis is made available by the University of Alberta Libraries with permission of the copyright owner solely for non-commercial purposes. This thesis, or any portion thereof, may not otherwise be copied or reproduced without the written consent of the copyright owner, except to the extent permitted by Canadian copyright law.
  • Language
  • Institution
    University of Alberta
  • Degree level
  • Department
    • Department of Mathematical and Statistical Sciences
  • Specialization
    • Applied Mathematics
  • Supervisor / co-supervisor and their department(s)
    • Safouhi, Hassan (Mathematical & Statistical Sciences)
    • Lau, Anthony To-Ming (Mathematical & Statistical Sciences)
  • Examining committee members and their departments
    • Dai, Feng (Mathematical & Statistical Sciences)
    • Weniger, Ernst Joachim (External Reader, University of Regensburg)
    • Brown, Alex (Chemistry)
    • Li, Michael (Mathematical & Statistical Sciences)