Efficient evaluation of the Coulomb force in the Gaussian and finite-element Coulomb method

Kurashige, Yuki; Nakajima, Takahito; Sato, Takeshi; Hirao, Kimihiko
June 2010
Journal of Chemical Physics;6/28/2010, Vol. 132 Issue 24, p244107
Academic Journal
We propose an efficient method for evaluating the Coulomb force in the Gaussian and finite-element Coulomb (GFC) method, which is a linear-scaling approach for evaluating the Coulomb matrix and energy in large molecular systems. The efficient evaluation of the analytical gradient in the GFC is not straightforward as well as the evaluation of the energy because the SCF procedure with the Coulomb matrix does not give a variational solution for the Coulomb energy. Thus, an efficient approximate method is alternatively proposed, in which the Coulomb potential is expanded in the Gaussian and finite-element auxiliary functions as done in the GFC. To minimize the error in the gradient not just in the energy, the derived functions of the original auxiliary functions of the GFC are used additionally for the evaluation of the Coulomb gradient. In fact, the use of the derived functions significantly improves the accuracy of this approach. Although these additional auxiliary functions enlarge the size of the discretized Poisson equation and thereby increase the computational cost, it maintains the near linear scaling as the GFC and does not affects the overall efficiency of the GFC approach.


Related Articles

  • Linear-scaling multipole-accelerated Gaussian and finite-element Coulomb method. Watson, Mark A.; Kurashige, Yuki; Nakajima, Takahito; Hirao, Kimihiko // Journal of Chemical Physics;2/7/2008, Vol. 128 Issue 5, p054105 

    A linear-scaling implementation of the Gaussian and finite-element Coulomb (GFC) method is presented for the rapid computation of the electronic Coulomb potential. The current work utilizes the fast multipole method (FMM) for the evaluation of the Poisson equation boundary condition. The FMM...

  • Gaussian and finite-element Coulomb method for the fast evaluation of Coulomb integrals. Kurashige, Yuki; Nakajima, Takahito; Hirao, Kimihiko // Journal of Chemical Physics;4/14/2007, Vol. 126 Issue 14, p144106 

    The authors propose a new linear-scaling method for the fast evaluation of Coulomb integrals with Gaussian basis functions called the Gaussian and finite-element Coulomb (GFC) method. In this method, the Coulomb potential is expanded in a basis of mixed Gaussian and finite-element auxiliary...

  • Relativistic density functional calculations using two-spinor minimax finite-element method and linear combination of atomic orbitals for ZnO, CdO, HgO, UubO and Cu2, Ag2, Au2, Rg2. Kullie, O.; Zhang, H.; Kolb, J.; Kolb, D. // Journal of Chemical Physics;12/28/2006, Vol. 125 Issue 24, p244303 

    In previous work the authors have presented a highly accurate two-spinor fully relativistic solution of the two-center Coulomb problem utilizing the finite-element method (FEM) and furthermore developed a relativistic minimax two-spinor linear combination of atomic orbitals (LCAO). In the...

  • The accuracy of the Gaussian-and-finite-element-Coulomb (GFC) method for the calculation of Coulomb integrals. Przybytek, Michal; Helgaker, Trygve // Journal of Chemical Physics;Aug2013, Vol. 139 Issue 5, p054114 

    We analyze the accuracy of the Coulomb energy calculated using the Gaussian-and-finite-element-Coulomb (GFC) method. In this approach, the electrostatic potential associated with the molecular electronic density is obtained by solving the Poisson equation and then used to calculate matrix...

  • Numerical solution of Poisson’s equation in polyatomic molecules. Becke, A. D.; Dickson, R. M. // Journal of Chemical Physics;9/1/1988, Vol. 89 Issue 5, p2993 

    We have developed a completely numerical scheme for solution of Poisson’s equation in multicenter systems. We are thus able to numerically calculate the Coulomb potential of arbitrary charge distributions in polyatomic molecules. The method is based on a decomposition of the multicenter...

  • NUMERICAL SOLUTION OF KOHN-SHAM EQUATION FOR ATOM. Romanowski, Zbigniew // Acta Physica Polonica B;2007, Vol. 38 Issue 10, p3263 

    The numerical algorithm solving Kohn--Sham (KS) equation for the free or confined atom on uniform mesh is presented. Using the spherical symmetry of system, the three-dimensional KS equation is reduced to the one-dimensional functional eigenproblem. The functional eigenproblem is transformed to...

  • Fast and accurate Coulomb calculation with Gaussian functions. Füsti-Molnár, László; Kong, Jing // Journal of Chemical Physics;2/15/2005, Vol. 122 Issue 7, p074108 

    Coulomb interaction is one of the major time-consuming components in a density functional theory (DFT) calculation. In the last decade, dramatic progresses have been made to improve the efficiency of Coulomb calculation, including continuous fast multipole method (CFMM) and J-engine method, all...

  • Smoothed Gaussian molecular fields: an evaluation of molecular alignment problems. Leherte, Laurence; Vercauteren, Daniel // Theoretical Chemistry Accounts: Theory, Computation, & Modeling;Aug2012, Vol. 131 Issue 8, p1 

    Several smoothed Gaussian-based descriptors used in a molecular superposition algorithm are presented. One descriptor, as detailed in a previous work (Leherte in J Comput Chem 27:1800-1816, ), is the full electron density approximated through the promolecular atomic shell approximation (PASA)...

  • Direct calculation of the Coulomb matrix: Slater-type orbitals. Ema, Ignacio; López, Rafael; Ramírez, Guillermo; Rico, Jaime // Theoretical Chemistry Accounts: Theory, Computation, & Modeling;Jan2011, Vol. 128 Issue 1, p115 

    It is proved that the evaluation of the Coulomb potential and the calculation of its matrix elements can be carried out in separate steps whose costs formally increase as the square of the number of basis functions. The resulting method for computing the Coulomb matrix is reported, and its main...


Read the Article


Sorry, but this item is not currently available from your library.

Try another library?
Sign out of this library

Other Topics