2010
DOI: 10.1063/1.3357981
|View full text |Cite
|
Sign up to set email alerts
|

An efficient algorithm for classical density functional theory in three dimensions: Ionic solutions

Abstract: Classical density functional theory ͑DFT͒ of fluids is a valuable tool to analyze inhomogeneous fluids. However, few numerical solution algorithms for three-dimensional systems exist. Here we present an efficient numerical scheme for fluids of charged, hard spheres that uses O͑N log N͒ operations and O͑N͒ memory, where N is the number of grid points. This system-size scaling is significant because of the very large N required for three-dimensional systems. The algorithm uses fast Fourier transforms ͑FFTs͒ to e… Show more

Help me understand this report
View preprint versions

Search citation statements

Order By: Relevance

Paper Sections

Select...
1
1
1

Citation Types

0
51
0

Year Published

2011
2011
2018
2018

Publication Types

Select...
7
2

Relationship

1
8

Authors

Journals

citations
Cited by 67 publications
(51 citation statements)
references
References 28 publications
0
51
0
Order By: Relevance
“…In general, the DFT equations must be solved numerically. This is more involved than for the PB or Bikerman models, but several straightforward algorithms have been described [19,49,52]. One important point to note is that the DFT Eqs.…”
Section: Dft Of Charged Hard Spheresmentioning
confidence: 98%
“…In general, the DFT equations must be solved numerically. This is more involved than for the PB or Bikerman models, but several straightforward algorithms have been described [19,49,52]. One important point to note is that the DFT Eqs.…”
Section: Dft Of Charged Hard Spheresmentioning
confidence: 98%
“…Over the years, a number of numerical approaches have been proposed for solving DFT equations. 9,27,[36][37][38] The crucial step in the numerical implementation is to accurately and efficiently compute the non-local integral terms corresponding to the fluid-fluid interactions, i.e., the fourth term in Eq. (14).…”
Section: Numerical Approachmentioning
confidence: 99%
“…This integration is in the form of convolution. Gillespie et al have used fast Fourier transform to deal with the convolution [28], while Meng et al have used the definition of Dirac delta function and change of variables to transform these 3D integrals into 2D integrals on spheres and remove the singularity in the integrands [40]. Though both algorithms can solve the integro-differential equations, they cost lots of computer memory and time during calculation.…”
mentioning
confidence: 99%