期刊
JOURNAL OF COMPUTATIONAL CHEMISTRY
卷 33, 期 25, 页码 1981-1996出版社
WILEY
DOI: 10.1002/jcc.23025
关键词
small angle scattering; Debye sums; fast multipole method; harmonic expansion; profile computation; Helmholtz equation; linear cost; SAXS; SANS; WAXS
资金
- National Institutes of Health [GM065334, GM095755]
- New Research Frontiers Award of the Institute of the Advanced Computer Studies of the University of Maryland
- Fantalgo, LLC
Debye summation, which involves the summation of sinc functions of distances between all pair of atoms in three-dimensional space, arises in computations performed in crystallography, small/wide angle X-ray scattering (SAXS/WAXS), and small angle neutron scattering (SANS). Direct evaluation of Debye summation has quadratic complexity, which results in computational bottleneck when determining crystal properties, or running structure refinement protocols that involve SAXS or SANS, even for moderately sized molecules. We present a fast approximation algorithm that efficiently computes the summation to any prescribed accuracy ? in linear time. The algorithm is similar to the fast multipole method (FMM), and is based on a hierarchical spatial decomposition of the molecule coupled with local harmonic expansions and translation of these expansions. An even more efficient implementation is possible when the scattering profile is all that is required, as in small angle scattering reconstruction (SAS) of macromolecules. We examine the relationship of the proposed algorithm to existing approximate methods for profile computations, and show that these methods may result in inaccurate profile computations, unless an error-bound derived in this article is used. Our theoretical and computational results show orders of magnitude improvement in computation complexity over existing methods, while maintaining prescribed accuracy. (c) 2012 Wiley Periodicals, Inc.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据