Abstract
Treecode algorithms efficiently approximate N-body interactions in O(N) or O(NlogN). In order to treat general 3D kernels, recent developments employ polynomial interpolation to approximate the kernels. The polynomials are a tensor product of 1-dimensional polynomials. Here, we develop an O(NlogN) tricubic interpolation based treecode method for 3D kernels. The tricubic interpolation is inherently three-dimensional and as such does not employ a tensor product. The form allows for easy evaluation of the derivatives of the kernel, required in dynamical simulations, which is not the case for the tensor product approach. We develop both a particle-cluster and cluster-particle variants and present results for the Coulomb, screened Coulomb and the real space Ewald kernels. We also present results of an MD simulation of a Lennard-Jones liquid using the tricubic treecode.
Sign-in/Register to access full text options 
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callMore From: Journal of Computational Mathematics and Data Science
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
