Efficient computation of the 2D periodic Green’s function using the Ewald method

• Published in: Journal of computational physics Vol. 219; no. 2; pp. 899 - 911
• Main Authors: Oroskar, Siddharth, Jackson, David R., Wilton, Donald R.
• Format: Journal Article
• Language: English
• Published: Amsterdam Elsevier Inc 10.12.2006; Elsevier
• Subjects: Accuracy; Asymptotic properties; Computational efficiency; Computational techniques; Convergence; Ewald method; Exact sciences and technology; Green's functions; Lattice; Mathematical methods in physics; Mathematical models; Optimization; Periodic Green’s function; Physics; Two dimensional; Ewald method; Periodic Green’s function; Lattice; Point sources; Periodic Green's function; Digital simulation; Periodic function; Calculation; Asymptotic convergence; Green function; Calculation methods
• ISSN: 0021-9991; 1090-2716
• DOI: 10.1016/j.jcp.2006.06.050

An efficient computation of the periodic Helmholtz Green’s function for a 2D array of point sources using the Ewald method is presented. Limitations on the numerical accuracy when using the “optimum” E parameter (which gives optimum asymptotic convergence) at high frequency are discussed. A “best” E parameter is then derived to overcome these limitations, which allows for the fastest convergence while maintaining a specific level of accuracy (loss of significant figures) in the final result. The actual loss of significant figures has been verified through numerical simulations. Formulas for the number of terms needed for convergence have also been derived for both the spectral and the spatial series that appear in the Ewald method and are found to be accurate in almost all cases.