1
|
Nakai H, Kobayashi M, Yoshikawa T, Seino J, Ikabata Y, Nishimura Y. Divide-and-Conquer Linear-Scaling Quantum Chemical Computations. J Phys Chem A 2023; 127:589-618. [PMID: 36630608 DOI: 10.1021/acs.jpca.2c06965] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
Fragmentation and embedding schemes are of great importance when applying quantum-chemical calculations to more complex and attractive targets. The divide-and-conquer (DC)-based quantum-chemical model is a fragmentation scheme that can be connected to embedding schemes. This feature article explains several DC-based schemes developed by the authors over the last two decades, which was inspired by the pioneering study of DC self-consistent field (SCF) method by Yang and Lee (J. Chem. Phys. 1995, 103, 5674-5678). First, the theoretical aspects of the DC-based SCF, electron correlation, excited-state, and nuclear orbital methods are described, followed by the two-component relativistic theory, quantum-mechanical molecular dynamics simulation, and the introduction of three programs, including DC-based schemes. Illustrative applications confirmed the accuracy and feasibility of the DC-based schemes.
Collapse
Affiliation(s)
- Hiromi Nakai
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo169-8555, Japan.,Waseda Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo169-8555, Japan
| | - Masato Kobayashi
- Department of Chemistry, Faculty of Science, Hokkaido University, Kita 10 Nishi 8, Kita-ku, Sapporo, Hokkaido060-0810, Japan.,Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Kita 21 Nishi 10, Kita-ku, Sapporo, Hokkaido001-0021, Japan
| | - Takeshi Yoshikawa
- Faculty of Pharmaceutical Sciences, Toho University, 2-2-1 Miyama, Funabashi, Chiba274-8510, Japan
| | - Junji Seino
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo169-8555, Japan.,Waseda Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo169-8555, Japan
| | - Yasuhiro Ikabata
- Information and Media Center, Toyohashi University of Technology, 1-1 Hibarigaoka, Tempaku-cho, Toyohashi, Aichi441-8580, Japan.,Department of Computer Science and Engineering, Toyohashi University of Technology, 1-1 Hibarigaoka, Tempaku-cho, Toyohashi, Aichi441-8580, Japan
| | - Yoshifumi Nishimura
- Waseda Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo169-8555, Japan
| |
Collapse
|
2
|
Nakai H. Development of Linear-Scaling Relativistic Quantum Chemistry Covering the Periodic Table. BULLETIN OF THE CHEMICAL SOCIETY OF JAPAN 2021. [DOI: 10.1246/bcsj.20210091] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Affiliation(s)
- Hiromi Nakai
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
- Waseda Research Institute for Science and Engineering (WISE), Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
- Elements Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, Katsura, Kyoto 615-8520, Japan
| |
Collapse
|
3
|
Le HA, Shiozaki T. Occupied-Orbital Fast Multipole Method for Efficient Exact Exchange Evaluation. J Chem Theory Comput 2018; 14:1228-1234. [DOI: 10.1021/acs.jctc.7b00880] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Hai-Anh Le
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States
| | - Toru Shiozaki
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States
| |
Collapse
|
4
|
Orimoto Y, Liu K, Aoki Y. Elongation method for electronic structure calculations of random DNA sequences. J Comput Chem 2015; 36:2103-13. [DOI: 10.1002/jcc.24047] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2015] [Revised: 07/20/2015] [Accepted: 07/23/2015] [Indexed: 01/01/2023]
Affiliation(s)
- Yuuichi Orimoto
- Department of Material Sciences, Faculty of Engineering Sciences; Kyushu University; 6-1 Kasuga-Park Fukuoka 816-8580 Japan
| | - Kai Liu
- Department of Material Sciences, Faculty of Engineering Sciences; Kyushu University; 6-1 Kasuga-Park Fukuoka 816-8580 Japan
| | - Yuriko Aoki
- Department of Material Sciences, Faculty of Engineering Sciences; Kyushu University; 6-1 Kasuga-Park Fukuoka 816-8580 Japan
- Japan Science and Technology Agency, CREST; 4-1-8 Hon-Chou Kawaguchi Saitama 332-0012 Japan
| |
Collapse
|
5
|
Orimoto Y, Yamamoto R, Xie P, Liu K, Imamura A, Aoki Y. Ab initio O(N) elongation-counterpoise method for BSSE-corrected interaction energy analyses in biosystems. J Chem Phys 2015; 142:104111. [DOI: 10.1063/1.4913931] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yuuichi Orimoto
- Department of Material Sciences, Faculty of Engineering Sciences, Kyushu University, 6-1 Kasuga-Park, Fukuoka 816-8580, Japan
| | - Ryohei Yamamoto
- Department of Molecular and Material Sciences, Interdisciplinary Graduate School of Engineering Sciences, Kyushu University, 6-1 Kasuga-Park, Fukuoka 816-8580, Japan
| | - Peng Xie
- Department of Material Sciences, Faculty of Engineering Sciences, Kyushu University, 6-1 Kasuga-Park, Fukuoka 816-8580, Japan
| | - Kai Liu
- Department of Material Sciences, Faculty of Engineering Sciences, Kyushu University, 6-1 Kasuga-Park, Fukuoka 816-8580, Japan
| | - Akira Imamura
- Hiroshima Kokusai Gakuin University, 6-20-1 Nakano, Aki-ku, Hiroshima 739-0321, Japan
| | - Yuriko Aoki
- Department of Material Sciences, Faculty of Engineering Sciences, Kyushu University, 6-1 Kasuga-Park, Fukuoka 816-8580, Japan
- Japan Science and Technology Agency, CREST, 4-1-8 Hon-chou, Kawaguchi, Saitama 332-0012, Japan
| |
Collapse
|
6
|
Toivanen EA, Losilla SA, Sundholm D. The grid-based fast multipole method – a massively parallel numerical scheme for calculating two-electron interaction energies. Phys Chem Chem Phys 2015; 17:31480-90. [DOI: 10.1039/c5cp01173f] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A grid-based fast multipole method has been developed for calculating two-electron interaction energies for non-overlapping charge densities.
Collapse
Affiliation(s)
| | | | - Dage Sundholm
- Department of Chemistry
- University of Helsinki
- Finland
| |
Collapse
|
7
|
|
8
|
SEINO J, NAKAI H. Large-Scale and Highly Accurate Relativistic Quantum-Chemical Scheme:toward Establishment ofTheoretical Foundation for Element Strategy. JOURNAL OF COMPUTER CHEMISTRY-JAPAN 2014. [DOI: 10.2477/jccj.2013-0010] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
|
9
|
Seino J, Nakai H. Local unitary transformation method toward practical electron correlation calculations with scalar relativistic effect in large-scale molecules. J Chem Phys 2013; 139:034109. [DOI: 10.1063/1.4813595] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
10
|
A modified localization scheme for the three-dimensional elongation method applied to large systems. Chem Phys Lett 2013. [DOI: 10.1016/j.cplett.2013.02.039] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
11
|
Liu K, Peng L, Gu FL, Aoki Y. Three dimensional elongation method for large molecular calculations. Chem Phys Lett 2013. [DOI: 10.1016/j.cplett.2012.12.046] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
12
|
Liu K, Inerbaev T, Korchowiec J, Gu FL, Aoki Y. Geometry optimization for large systems by the elongation method. Theor Chem Acc 2012. [DOI: 10.1007/s00214-012-1277-9] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
13
|
Choi CH, Fedorov DG. Reducing the scaling of the fragment molecular orbital method using the multipole method. Chem Phys Lett 2012. [DOI: 10.1016/j.cplett.2012.06.018] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
14
|
Kobayashi M, Nakai H. How does it become possible to treat delocalized and/or open-shell systems in fragmentation-based linear-scaling electronic structure calculations? The case of the divide-and-conquer method. Phys Chem Chem Phys 2012; 14:7629-39. [DOI: 10.1039/c2cp40153c] [Citation(s) in RCA: 66] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
|
15
|
Aoki Y, Gu FL. An elongation method for large systems toward bio-systems. Phys Chem Chem Phys 2012; 14:7640-68. [DOI: 10.1039/c2cp24033e] [Citation(s) in RCA: 50] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
16
|
Rubensson EH, Rudberg E, Salek P. Methods for Hartree-Fock and Density Functional Theory Electronic Structure Calculations with Linearly Scaling Processor Time and Memory Usage. CHALLENGES AND ADVANCES IN COMPUTATIONAL CHEMISTRY AND PHYSICS 2011. [DOI: 10.1007/978-90-481-2853-2_12] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
17
|
Divide-and-Conquer Approaches to Quantum Chemistry: Theory and Implementation. CHALLENGES AND ADVANCES IN COMPUTATIONAL CHEMISTRY AND PHYSICS 2011. [DOI: 10.1007/978-90-481-2853-2_5] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
|
18
|
Kobayashi M, Yoshikawa T, Nakai H. Divide-and-conquer self-consistent field calculation for open-shell systems: Implementation and application. Chem Phys Lett 2010. [DOI: 10.1016/j.cplett.2010.10.005] [Citation(s) in RCA: 33] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
19
|
Dachsel H. Corrected Article: “An error-controlled fast multipole method” [J. Chem. Phys.131, 244102 (2009)]. J Chem Phys 2010; 132:119901. [DOI: 10.1063/1.3358272] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
20
|
Mullin JM, Roskop LB, Pruitt SR, Collins MA, Gordon MS. Systematic fragmentation method and the effective fragment potential: an efficient method for capturing molecular energies. J Phys Chem A 2010; 113:10040-9. [PMID: 19739681 DOI: 10.1021/jp9036183] [Citation(s) in RCA: 61] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The systematic fragmentation method fragments a large molecular system into smaller pieces, in such a way as to greatly reduce the computational cost while retaining nearly the accuracy of the parent ab initio electronic structure method. In order to attain the desired (sub-kcal/mol) accuracy, one must properly account for the nonbonded interactions between the separated fragments. Since, for a large molecular species, there can be a great many fragments and therefore a great many nonbonded interactions, computations of the nonbonded interactions can be very time-consuming. The present work explores the efficacy of employing the effective fragment potential (EFP) method to obtain the nonbonded interactions since the EFP method has been shown previously to capture nonbonded interactions with an accuracy that is often comparable to that of second-order perturbation theory. It is demonstrated that for nonbonded interactions that are not high on the repulsive wall (generally >2.7 A), the EFP method appears to be a viable approach for evaluating the nonbonded interactions. The efficacy of the EFP method for this purpose is illustrated by comparing the method to ab initio methods for small water clusters, the ZOVGAS molecule, retinal, and the alpha-helix. Using SFM with EFP for nonbonded interactions yields an error of 0.2 kcal/mol for the retinal cis-trans isomerization and a mean error of 1.0 kcal/mol for the isomerization energies of five small (120-170 atoms) alpha-helices.
Collapse
|
21
|
Time-dependent Hartree–Fock frequency-dependent polarizability calculation applied to divide-and-conquer electronic structure method. Chem Phys Lett 2010. [DOI: 10.1016/j.cplett.2009.12.043] [Citation(s) in RCA: 33] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
22
|
Abstract
We present a two-stage error estimation scheme for the fast multipole method (FMM). This scheme can be applied to any particle system. It incorporates homogeneous as well as inhomogeneous distributions. The FMM error as a consequence of the finite representation of the multipole expansions and the operator error is correlated with an absolute or relative user-requested energy threshold. Such a reliable error control is the basis for making reliable simulations in computational physics. Our FMM program on the basis of the two-stage error estimation scheme is available on request.
Collapse
Affiliation(s)
- Holger Dachsel
- Institute for Advanced Simulation, Jülich Supercomputing Centre, Forschungszentrum Jülich, 52425 Jülich, Germany.
| |
Collapse
|
23
|
Korchowiec J, Lewandowski J, Makowski M, Gu FL, Aoki Y. Elongation cutoff technique armed with quantum fast multipole method for linear scaling. J Comput Chem 2009; 30:2515-25. [DOI: 10.1002/jcc.21252] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
|
24
|
Neese F, Wennmohs F, Hansen A, Becker U. Efficient, approximate and parallel Hartree–Fock and hybrid DFT calculations. A ‘chain-of-spheres’ algorithm for the Hartree–Fock exchange. Chem Phys 2009. [DOI: 10.1016/j.chemphys.2008.10.036] [Citation(s) in RCA: 1551] [Impact Index Per Article: 103.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
25
|
KOBAYASHI M, AKAMA T, NAKAI H. Implementation of Divide-and-Conquer (DC) Electronic Structure Code to GAMESS Program Package. ACTA ACUST UNITED AC 2009. [DOI: 10.2477/jccj.h2027] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
|
26
|
Rudberg E, Rubensson EH, Sałek P. Hartree–Fock calculations with linearly scaling memory usage. J Chem Phys 2008; 128:184106. [DOI: 10.1063/1.2918357] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
27
|
Rubensson EH, Rudberg E, Sałek P. Density matrix purification with rigorous error control. J Chem Phys 2008; 128:074106. [DOI: 10.1063/1.2826343] [Citation(s) in RCA: 37] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
28
|
Kobayashi M, Akama T, Nakai H. Second-order Møller-Plesset perturbation energy obtained from divide-and-conquer Hartree-Fock density matrix. J Chem Phys 2006; 125:204106. [PMID: 17144689 DOI: 10.1063/1.2388261] [Citation(s) in RCA: 66] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The density matrix (DM) obtained from Yang's [Phys. Rev. Lett. 66, 1438 (1991)] divide-and-conquer (DC) Hartree-Fock (HF) calculation is applied to the explicit second-order Møller-Plesset perturbation (MP2) energy functional of the HF DM, which was firstly mentioned by Ayala and Scuseria [J. Chem. Phys. 110, 3660 (1999)] and was improved by Surján [Chem. Phys. Lett. 406, 318 (2005)] as DM-Laplace MP2. This procedure, termed DC-DM MP2, requires the HF DM of holes, for which we propose two evaluation schemes in DC manner. Numerical studies reveal that the DC-DM MP2 energy deviation from canonical MP2 is the same order of magnitude as DC-HF energy deviation from conventional HF whichever type of hole DM is adopted. It is also confirmed that the central processing unit time of DC-DM MP2 is less than that of DM-Laplace MP2 because the DC-HF DM is sparser than conventional DM.
Collapse
Affiliation(s)
- Masato Kobayashi
- Department of Chemistry, School of Science and Engineering, Waseda University, Tokyo 169-8555, Japan
| | | | | |
Collapse
|
29
|
Abstract
A number of computational techniques are described that reduce the effort related to the continuous fast multipole method, used for the evaluation of Coulomb matrix elements as needed in Hartree-Fock and density functional theories. A new extent definition for Gaussian charge distributions is proposed, as well as a new way of dividing distributions into branches. Also, a new approach for estimating the error caused by truncation of multipole expansions is presented. It is found that the use of dynamically truncated multipole expansions gives a speedup of a factor of 10 in the work required for multipole interactions, compared to the case when all interactions are computed using a fixed multipole expansion order. Results of benchmark calculations on three-dimensional systems are reported, demonstrating the usefulness of our present implementation of the fast multipole method.
Collapse
Affiliation(s)
- Elias Rudberg
- Department of Chemistry, University of Warwick, Coventry CV4 7AL, United Kingdom
| | | |
Collapse
|
30
|
Gan CK, Tymczak CJ, Challacombe M. Linear scaling computation of the Fock matrix. VII. Parallel computation of the Coulomb matrix. J Chem Phys 2004; 121:6608-14. [PMID: 15473715 DOI: 10.1063/1.1790891] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present parallelization of a quantum-chemical tree-code for linear scaling computation of the Coulomb matrix. Equal time partition is used to load balance computation of the Coulomb matrix. Equal time partition is a measurement based algorithm for domain decomposition that exploits small variation of the density between self-consistent-field cycles to achieve load balance. Efficiency of the equal time partition is illustrated by several tests involving both finite and periodic systems. It is found that equal time partition is able to deliver 91%-98% efficiency with 128 processors in the most time consuming part of the Coulomb matrix calculation. The current parallel quantum chemical tree code is able to deliver 63%-81% overall efficiency on 128 processors with fine grained parallelism (less than two heavy atoms per processor).
Collapse
Affiliation(s)
- Chee Kwan Gan
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA.
| | | | | |
Collapse
|
31
|
Kudin KN, Scuseria GE. Revisiting infinite lattice sums with the periodic fast multipole method. J Chem Phys 2004; 121:2886-90. [PMID: 15291598 DOI: 10.1063/1.1771634] [Citation(s) in RCA: 36] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The evaluation of lattice sums as well as stress lattice sums encountered in the periodic fast multipole method is reinvestigated. Simple, accurate, and efficient recurrence expressions for such sums following the ideas of the renormalization method are derived. The first few nonzero lattice sum terms in a three-dimensional cubic lattice are computed and given in Tables. The practical considerations accompanying the computation of the sums such as convergence and accuracy are discussed.
Collapse
|
32
|
Choi CH. Direct determination of multipole moments of Cartesian Gaussian functions in spherical polar coordinates. J Chem Phys 2004; 120:3535-43. [PMID: 15268515 DOI: 10.1063/1.1642597] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A new way of generating the multipole moments of Cartesian Gaussian functions in spherical polar coordinates has been established, bypassing the intermediary of Cartesian moment tensors. A new set of recurrence relations have also been derived for the resulting analytic integral values. The new method furnishes a conceptually simple and numerically efficient evaluation procedure for the multipole moments. The advantages over existing methods are documented. The results are relevant for the linear scaling quantum theories based on the fast multipole method.
Collapse
Affiliation(s)
- Cheol Ho Choi
- Department of Chemistry, College of Natural Sciences, Kyungpook National University, Taegu 702-701, South Korea.
| |
Collapse
|
33
|
|
34
|
|
35
|
Amisaki T, Toyoda S, Miyagawa H, Kitamura K. Development of hardware accelerator for molecular dynamics simulations: a computation board that calculates nonbonded interactions in cooperation with fast multipole method. J Comput Chem 2003; 24:582-92. [PMID: 12632472 DOI: 10.1002/jcc.10193] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
Evaluation of long-range Coulombic interactions still represents a bottleneck in the molecular dynamics (MD) simulations of biological macromolecules. Despite the advent of sophisticated fast algorithms, such as the fast multipole method (FMM), accurate simulations still demand a great amount of computation time due to the accuracy/speed trade-off inherently involved in these algorithms. Unless higher order multipole expansions, which are extremely expensive to evaluate, are employed, a large amount of the execution time is still spent in directly calculating particle-particle interactions within the nearby region of each particle. To reduce this execution time for pair interactions, we developed a computation unit (board), called MD-Engine II, that calculates nonbonded pairwise interactions using a specially designed hardware. Four custom arithmetic-processors and a processor for memory manipulation ("particle processor") are mounted on the computation board. The arithmetic processors are responsible for calculation of the pair interactions. The particle processor plays a central role in realizing efficient cooperation with the FMM. The results of a series of 50-ps MD simulations of a protein-water system (50,764 atoms) indicated that a more stringent setting of accuracy in FMM computation, compared with those previously reported, was required for accurate simulations over long time periods. Such a level of accuracy was efficiently achieved using the cooperative calculations of the FMM and MD-Engine II. On an Alpha 21264 PC, the FMM computation at a moderate but tolerable level of accuracy was accelerated by a factor of 16.0 using three boards. At a high level of accuracy, the cooperative calculation achieved a 22.7-fold acceleration over the corresponding conventional FMM calculation. In the cooperative calculations of the FMM and MD-Engine II, it was possible to achieve more accurate computation at a comparable execution time by incorporating larger nearby regions.
Collapse
Affiliation(s)
- Takashi Amisaki
- Department of Biological Regulation, Faculty of Medicine, Tottori University, 86 Nishi-machi, Yonago, Tottori 683-8503, Japan
| | | | | | | |
Collapse
|