Endre søk
RefereraExporteraLink to record
Permanent link

Direct link
Referera
Referensformat
  • apa
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Annet format
Fler format
Språk
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Annet språk
Fler språk
Utmatningsformat
  • html
  • text
  • asciidoc
  • rtf
The Spectral Ewald method for singly periodic domains
KTH, Skolan för teknikvetenskap (SCI), Matematik (Inst.).
KTH, Skolan för teknikvetenskap (SCI), Matematik (Inst.), Numerisk analys, NA.
2017 (engelsk)Inngår i: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 347, s. 341-366Artikkel i tidsskrift (Fagfellevurdert) Published
Abstract [en]

We present a fast and spectrally accurate method for efficient computation of the three dimensional Coulomb potential with periodicity in one direction. The algorithm is FFT-based and uses the so-called Ewald decomposition, which is naturally most efficient for the triply periodic case. In this paper, we show how to extend the triply periodic Spectral Ewald method to the singly periodic case, such that the cost of computing the singly periodic potential is only marginally larger than the cost of computing the potential for the corresponding triply periodic system. In the Fourier space contribution of the Ewald decomposition, a Fourier series is obtained in the periodic direction with a Fourier integral over the non-periodic directions for each discrete wave number. We show that upsampling to resolve the integral is only needed for modes with small wave numbers. For the zero wave number, this Fourier integral has a singularity. For this mode, we effectively need to solve a free-space Poisson equation in two dimensions. A very recent idea by Vico et al. makes it possible to use FFTs to solve this problem, allowing us to unify the treatment of all modes. An adaptive 3D FFT can be established to apply different upsampling rates locally. The computational cost for other parts of the algorithm is essentially unchanged as compared to the triply periodic case, in total yielding only a small increase in both computational cost and memory usage for this singly periodic case.

sted, utgiver, år, opplag, sider
Academic Press, 2017. Vol. 347, s. 341-366
Emneord [en]
Adaptive FFT, Coulomb potentials, Fast Ewald summation, Fast Fourier transform, Fourier integral, Single periodic, Spectral accuracy
HSV kategori
Identifikatorer
URN: urn:nbn:se:kth:diva-212215DOI: 10.1016/j.jcp.2017.07.001ISI: 000408045500017Scopus ID: 2-s2.0-85024927535OAI: oai:DiVA.org:kth-212215DiVA, id: diva2:1134039
Forskningsfinansiär
Swedish Research Council, 2011-3178Swedish e‐Science Research Center
Merknad

QC 20170817

Tilgjengelig fra: 2017-08-17 Laget: 2017-08-17 Sist oppdatert: 2017-12-12bibliografisk kontrollert
Inngår i avhandling
1. Fast methods for electrostatic calculations in molecular dynamics simulations
Åpne denne publikasjonen i ny fane eller vindu >>Fast methods for electrostatic calculations in molecular dynamics simulations
2018 (engelsk)Doktoravhandling, med artikler (Annet vitenskapelig)
Abstract [en]

This thesis deals with fast and efficient methods for electrostatic calculations with application in molecular dynamics simulations. The electrostatic calculations are often the most expensive part of MD simulations of charged particles. Therefore, fast and efficient algorithms are required to accelerate these calculations. In this thesis, two types of methods have been considered: FFT-based methods and fast multipole methods (FMM).

The major part of this thesis deals with fast N.log(N) and spectrally accurate methods for accelerating the computation of pairwise interactions with arbitrary periodicity. These methods are based on the Ewald decomposition and have been previously introduced for triply and doubly periodic problems under the name of Spectral Ewald (SE) method. We extend the method for problems with singly periodic boundary conditions, in which one of three dimensions is periodic. By introducing an adaptive fast Fourier transform, we reduce the cost of upsampling in the non periodic directions and show that the total cost of computation is comparable with the triply periodic counterpart. Using an FFT-based technique for solving free-space harmonic problems, we are able to unify the treatment of zero and nonzero Fourier modes for the doubly and singly periodic problems. Applying the same technique, we extend the SE method for cases with free-space boundary conditions, i.e. without any periodicity.

This thesis is also concerned with the fast multipole method (FMM) for electrostatic calculations. The FMM is very efficient for parallel processing but it introduces irregularities in the electrostatic potential and force, which can cause an energy drift in MD simulations. In this part of the thesis we introduce a regularized version of the FMM, useful for MD simulations, which approximately conserves energy over a long time period and even for low accuracy requirements. The method introduces a smooth transition over the boundary of boxes in the FMM tree and therefore it removes the discontinuity at the error level inherent in the FMM.

sted, utgiver, år, opplag, sider
Stockholm: KTH Royal Institute of Technology, 2018. s. 58
Serie
TRITA-MAT-A ; 2018:02
HSV kategori
Forskningsprogram
Tillämpad matematik och beräkningsmatematik
Identifikatorer
urn:nbn:se:kth:diva-219775 (URN)978-91-7729-640-9 (ISBN)
Disputas
2018-01-26, F3, Lindstedtsvägen 26, Stockholm, 10:00 (engelsk)
Opponent
Veileder
Forskningsfinansiär
Swedish e‐Science Research Center
Merknad

QC 20171213

Tilgjengelig fra: 2017-12-13 Laget: 2017-12-12 Sist oppdatert: 2017-12-13bibliografisk kontrollert

Open Access i DiVA

Fulltekst mangler i DiVA

Andre lenker

Forlagets fulltekstScopus

Søk i DiVA

Av forfatter/redaktør
Saffar Shamshirgar, DavoudTornberg, Anna-Karin
Av organisasjonen
I samme tidsskrift
Journal of Computational Physics

Søk utenfor DiVA

GoogleGoogle Scholar

doi
urn-nbn

Altmetric

doi
urn-nbn
Totalt: 53 treff
RefereraExporteraLink to record
Permanent link

Direct link
Referera
Referensformat
  • apa
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Annet format
Fler format
Språk
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Annet språk
Fler språk
Utmatningsformat
  • html
  • text
  • asciidoc
  • rtf