1
|
Jin WQ, Mao F, Li SM, Zuo WQ, Chen RD, Xiong GG, Mao H, Wang F, Zhang FS. The contribution of inner electron excitation to the electronic stopping power of palladium for protons. Phys Chem Chem Phys 2023; 25:9043-9050. [PMID: 36883618 DOI: 10.1039/d2cp05510d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/09/2023]
Abstract
The electronic stopping power of palladium (Pd) for protons is investigated based on time-dependent density functional theory combined with Ehrenfest molecular dynamics simulations. The electronic stopping power of Pd with explicitly considering inner electrons for protons is calculated and the excitation mechanism for the inner electrons of Pd is revealed. The velocity proportionality of the low-energy stopping power of Pd is reproduced. Our study verified that the inner electron excitation contributes significantly to the electronic stopping power of Pd in the high energy range, which is strongly dependent on the impact parameter. The electronic stopping power obtained from the off-channeling geometry is in quantitative agreement with the experimental data in a wide velocity range, and the discrepancy around the stopping maximum is further reduced by considering the relativistic correction on the binding energy of inner electrons. The velocity dependence of the mean steady-state charge of protons is quantified, and the results showed that the participation of 4p-electrons reduces the mean steady-state charge of protons, and consequently decreases the electronic stopping power of Pd in the low energy range.
Collapse
Affiliation(s)
- Wen-Qi Jin
- School of Nuclear Science and Technology, University of South China, Hengyang 421001, China.
| | - Fei Mao
- School of Nuclear Science and Technology, University of South China, Hengyang 421001, China.
| | - Shi-Ming Li
- School of Nuclear Science and Technology, University of South China, Hengyang 421001, China.
| | - Wen-Qi Zuo
- School of Nuclear Science and Technology, University of South China, Hengyang 421001, China.
| | - Rui-Da Chen
- School of Nuclear Science and Technology, University of South China, Hengyang 421001, China.
| | - Ge-Ge Xiong
- School of Nuclear Science and Technology, University of South China, Hengyang 421001, China.
| | - Hong Mao
- College of Mechanical Engineering, Hunan Institute of Science and Technology, Yueyang 414006, China
| | - Feng Wang
- School of Physics, Beijing Institute of Technology, Beijing 100081, China
| | - Feng-Shou Zhang
- The Key Laboratory of Beam Technology of Ministry of Education, College of Nuclear Science and Technology, Beijing Normal University, Beijing 100875, China
| |
Collapse
|
2
|
Andoh Y, Ichikawa SI, Sakashita T, Yoshii N, Okazaki S. Algorithm to minimize MPI communications in the parallelized fast multipole method combined with molecular dynamics calculations. J Comput Chem 2021; 42:1073-1087. [PMID: 33780021 DOI: 10.1002/jcc.26524] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/09/2020] [Revised: 03/09/2021] [Accepted: 03/11/2021] [Indexed: 11/07/2022]
Abstract
In the era of exascale supercomputers, large-scale, and long-time molecular dynamics (MD) calculations are expected to make breakthroughs in various fields of science and technology. Here, we propose a new algorithm to improve the parallelization performance of message passing interface (MPI)-communication in the MPI-parallelized fast multipole method (FMM) combined with MD calculations under three-dimensional periodic boundary conditions. Our approach enables a drastic reduction in the amount of communication data, including the atomic coordinates and multipole coefficients, both of which are required to calculate the electrostatic interaction by using the FMM. In communications of multipole coefficients, the reduction rate of communication data in the new algorithm relative to the amount of data in the conventional one increases as both the number of FMM levels and the number of MPI processes increase. The aforementioned rate increase could exceed 50% as the number of MPI processes becomes larger for very large systems. The proposed algorithm, named the minimum-transferred data (MTD) method, should enable large-scale and long-time MD calculations to be calculated efficiently, under the condition of massive MPI-parallelization on exascale supercomputers.
Collapse
Affiliation(s)
- Yoshimichi Andoh
- Center for Computational Science, Graduate School of Engineering, Nagoya University, Nagoya, Japan
| | - Shin-Ichi Ichikawa
- Computational Science Division, Technical Computing Business Unit, Fujitsu Limited, Chiba, Japan
| | - Tatsuya Sakashita
- Department of Information and Communication Technology, College of Engineering, Tamagawa University, Machida, Tokyo, Japan
| | - Noriyuki Yoshii
- Center for Computational Science, Graduate School of Engineering, Nagoya University, Nagoya, Japan
- Department of Materials Chemistry, Nagoya University, Nagoya, Japan
| | - Susumu Okazaki
- Department of Materials Chemistry, Nagoya University, Nagoya, Japan
- Department of Advanced Materials Science, The University of Tokyo, Kashiwa, Chiba, Japan
| |
Collapse
|
3
|
Urano R, Shinoda W, Yoshii N, Okazaki S. Exact long-range Coulombic energy calculation for net charged systems neutralized by uniformly distributed background charge using fast multipole method and its application to efficient free energy calculation. J Chem Phys 2020; 152:244115. [DOI: 10.1063/5.0007957] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Ryo Urano
- Department of Materials Chemistry, Nagoya University, Nagoya 464-8603, Japan
| | - Wataru Shinoda
- Department of Materials Chemistry, Nagoya University, Nagoya 464-8603, Japan
| | - Noriyuki Yoshii
- Department of Materials Chemistry, Nagoya University, Nagoya 464-8603, Japan
- Center for Computational Science, Graduate School of Engineering, Nagoya University, Nagoya 464-8603, Japan
| | - Susumu Okazaki
- Department of Materials Chemistry, Nagoya University, Nagoya 464-8603, Japan
- Center for Computational Science, Graduate School of Engineering, Nagoya University, Nagoya 464-8603, Japan
| |
Collapse
|
4
|
Andoh Y, Yoshii N, Okazaki S. Extension of the fast multipole method for the rectangular cells with an anisotropic partition tree structure. J Comput Chem 2020; 41:1353-1367. [PMID: 32100899 DOI: 10.1002/jcc.26180] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2019] [Revised: 02/04/2020] [Accepted: 02/16/2020] [Indexed: 02/06/2023]
Abstract
The fast multipole method (FMM) is an order N method for the numerically rigorous calculation of the electrostatic interactions among point charges in a system of interest. The FMM is utilized for massively parallelized software for molecular dynamics (MD) calculations. However, an inconvenient limitation is imposed on the implementation of the FMM: In three-dimensional case, a cubic MD unit cell is hierarchically divided by the octree partitioning under isotropic periodic boundary conditions along three axes. Here, we extended the FMM algorithm adaptive to a rectangular MD unit cell with different periodicity along the axes by applying an anisotropic hierarchical partitioning. The algorithm was implemented into the parallelized general-purpose MD calculation software designed for a system with uniform distribution of point charges in the unit cell. The partition tree can be a mixture of binary and ternary branches, the branches being chosen arbitrarily with respect to the coordinate axes at any levels. Errors in the calculated electrostatic interactions are discussed in detail for a selected partition tree structure. The extension enables us to execute MD calculations under more general conditions for the shape of the unit cell, partition tree, and boundary conditions, keeping the accuracy of the calculated electrostatic interactions as high as that with the conventional FMM. An extension of the present FMM algorithm to other prime number branches, such as 5 and 7, is straightforward.
Collapse
Affiliation(s)
- Yoshimichi Andoh
- Center of Computational Science, Graduate School of Engineering, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Japan
| | - Noriyuki Yoshii
- Center of Computational Science, Graduate School of Engineering, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Japan
| | - Susumu Okazaki
- Department of Materials Chemistry, Graduate School of Engineering, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Japan
| |
Collapse
|
5
|
Yoshii N, Andoh Y, Okazaki S. Fast multipole method for three-dimensional systems with periodic boundary condition in two directions. J Comput Chem 2020; 41:940-948. [PMID: 31930548 DOI: 10.1002/jcc.26141] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/14/2019] [Revised: 12/01/2019] [Accepted: 12/02/2019] [Indexed: 12/11/2022]
Abstract
We derived a new expression for the electrostatic interaction of three-dimensional charge-neutral systems with two-dimensional periodic boundary conditions (slab geometry) using a fast multipole method (FMM). Contributions from all the image cells are expressed as a sum of real and reciprocal space terms, and a self-interaction term. The reciprocal space contribution consists of two parts: zero and nonzero terms of the absolute value of the reciprocal lattice vector. To test the new expressions, electrostatic interactions were calculated for a randomly placed charge distribution in a cubic box and liquid water produced by molecular dynamics calculation. The accuracy could be controlled by the degree of expansion of the FMM. In the present expression, the computational complexity of the electrostatic interaction of N-particle systems is order N, which is superior to that of the conventional two-dimensional periodic Ewald method for a slab geometry and the particle mesh Ewald method with a large empty space at an interface of the unit cell. © 2020 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Noriyuki Yoshii
- Center for Computational Science, Graduate School of Engineering, Nagoya University, Nagoya, 464-8603, Japan.,Department of Materials Chemistry, Nagoya University, Nagoya, Japan
| | - Yoshimichi Andoh
- Center for Computational Science, Graduate School of Engineering, Nagoya University, Nagoya, 464-8603, Japan
| | - Susumu Okazaki
- Center for Computational Science, Graduate School of Engineering, Nagoya University, Nagoya, 464-8603, Japan.,Department of Materials Chemistry, Nagoya University, Nagoya, Japan
| |
Collapse
|
6
|
Yoshii N, Andoh Y, Okazaki S. Pressure tensor for electrostatic interaction calculated by fast multipole method with periodic boundary condition. J Comput Chem 2018; 39:1192-1199. [PMID: 29411407 DOI: 10.1002/jcc.25179] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2017] [Revised: 01/06/2018] [Accepted: 01/09/2018] [Indexed: 01/18/2023]
Abstract
A microscopic expression of the pressure tensor using the fast multipole method (FMM) with periodic boundary conditions has been derived. The pressure tensor calculated using this expression has been compared with that obtained using the Ewald method with high accuracy. The precision of the pressure tensor can be controlled as a function of expansion order p of FMM. Using the calculated pressure tensor, the constant pressure molecular dynamics calculation with fully fluctuating cell can be performed for anisotropic systems such as crystals, metals, liquid crystals, glasses, polymers, lipid bilayers, and interfacial regions between two phases. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Noriyuki Yoshii
- Graduate School of Engineering, Center for Computational Science, Nagoya University, Nagoya, 464-8603, Japan.,Department of Applied Chemistry, Nagoya University, Nagoya, 464-8603, Japan
| | - Yoshimichi Andoh
- Graduate School of Engineering, Center for Computational Science, Nagoya University, Nagoya, 464-8603, Japan
| | - Susumu Okazaki
- Graduate School of Engineering, Center for Computational Science, Nagoya University, Nagoya, 464-8603, Japan.,Department of Applied Chemistry, Nagoya University, Nagoya, 464-8603, Japan
| |
Collapse
|
7
|
Kawada S, Fujimoto K, Yoshii N, Okazaki S. Molecular dynamics study of the potential of mean force of SDS aggregates. J Chem Phys 2018; 147:084903. [PMID: 28863544 DOI: 10.1063/1.4998549] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In our previous study, all-atomistic molecular dynamics (MD) calculations have been carried out for the aggregation of ionic sodium dodecyl sulfate in water [S. Kawada et al., Chem. Phys. Lett. 646, 36 (2016)]. Aggregates of 20-30 dodecyl sulfate ions were formed within a short MD run for 10 ns. However, further aggregation did not occur despite a long MD calculation for more than 100 ns. This suggests that strong electrostatic repulsive interactions between the aggregates prevent the fusion of the aggregates. In the present study, mean force and potential of mean force acting between two aggregates with aggregation number N = 30 were evaluated as a function of their separation by MD calculations. The repulsive force becomes strong with decreasing distance between the two aggregates before they merge into one. An origin of the repulsive force is an electric double layer formed by the sulfate group and counter sodium ions. Strength of the repulsive force is in good agreement with the theoretical value given by the Derjaguin-Landau-Verwey-Overbeek (DLVO) theory. Once the aggregates establish contact, the force between them turns to be a large attractive force that can be explained by the interfacial tension. In order to form a single micelle from the two aggregates, it is necessary for them to climb over a free energy barrier of 23 kJ/mol. Once, the barrier is overcome, the micelle is stabilized by ∼200 kJ/mol. The time constant of aggregation evaluated from the calculated free energy barrier was about 28 μs at the concentration in our previous study.
Collapse
Affiliation(s)
- Shinji Kawada
- Department of Applied Chemistry, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8603, Japan
| | - Kazushi Fujimoto
- Department of Applied Chemistry, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8603, Japan
| | - Noriyuki Yoshii
- Department of Applied Chemistry, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8603, Japan
| | - Susumu Okazaki
- Department of Applied Chemistry, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8603, Japan
| |
Collapse
|
8
|
Akimov AV, Prezhdo OV. Large-Scale Computations in Chemistry: A Bird’s Eye View of a Vibrant Field. Chem Rev 2015; 115:5797-890. [DOI: 10.1021/cr500524c] [Citation(s) in RCA: 159] [Impact Index Per Article: 17.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Affiliation(s)
- Alexey V. Akimov
- Department
of Chemistry, University of South California, Los Angeles, California 90089, United States
| | - Oleg V. Prezhdo
- Department
of Chemistry, University of South California, Los Angeles, California 90089, United States
| |
Collapse
|
9
|
Lorenzen K, Schwörer M, Tröster P, Mates S, Tavan P. Optimizing the Accuracy and Efficiency of Fast Hierarchical Multipole Expansions for MD Simulations. J Chem Theory Comput 2012; 8:3628-36. [PMID: 26593008 DOI: 10.1021/ct300080n] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Based on p'th order Cartesian Taylor expansions of Coulomb interactions and on hierarchical decompositions of macromolecular simulation systems into hierarchies of nested, structure-adapted, and adaptively formed clusters of increasing size, fast multipole methods are constructed for rapid and accurate calculations of electrostatic interactions. These so-called SAMMp algorithms are formulated through totally symmetric and traceless tensors describing the multipole moments and the coefficients of local Taylor expansions. Simple recursions for the efficient evaluation and shifting of multipole moments are given. The required tensors are explicitly given up to order p = 4. The SAMMp algorithms are shown to guarantee the reaction principle. For systems with periodic boundaries, a reaction field (RF) correction is applied, which introduces at distances beyond the "minimum image convention" boundary a dielectric continuum surrounding each cluster at the top level of coarse graining. The correctness of the present SAMMp implementation is demonstrated by analyzing the scaling of the residuals and by checking the numerical accuracy of the reaction principle for a pair of distant molecular ions in vacuum. Molecular dynamics simulations of pure water and aqueous solutions containing artificial ions, which are enclosed by periodic boundaries, demonstrate the stability and low-noise behavior of SAMMp/RF.
Collapse
Affiliation(s)
- Konstantin Lorenzen
- Lehrstuhl für Biomolekulare Optik, Ludwig-Maximilians-Universität , Oettingenstrasse 67, 80538 München, Germany
| | - Magnus Schwörer
- Lehrstuhl für Biomolekulare Optik, Ludwig-Maximilians-Universität , Oettingenstrasse 67, 80538 München, Germany
| | - Philipp Tröster
- Lehrstuhl für Biomolekulare Optik, Ludwig-Maximilians-Universität , Oettingenstrasse 67, 80538 München, Germany
| | - Simon Mates
- Lehrstuhl für Biomolekulare Optik, Ludwig-Maximilians-Universität , Oettingenstrasse 67, 80538 München, Germany
| | - Paul Tavan
- Lehrstuhl für Biomolekulare Optik, Ludwig-Maximilians-Universität , Oettingenstrasse 67, 80538 München, Germany
| |
Collapse
|
10
|
Kurzak J, Pettitt BM. Fast multipole methods for particle dynamics. MOLECULAR SIMULATION 2006; 32:775-790. [PMID: 19194526 DOI: 10.1080/08927020600991161] [Citation(s) in RCA: 17] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
Abstract
The growth of simulations of particle systems has been aided by advances in computer speed and algorithms. The adoption of O(N) algorithms to solve N-body simulation problems has been less rapid due to the fact that such scaling was only competitive for relatively large N. Our work seeks to find algorithmic modifications and practical implementations for intermediate values of N in typical use for molecular simulations. This article reviews fast multipole techniques for calculation of electrostatic interactions in molecular systems. The basic mathematics behind fast summations applied to long ranged forces is presented along with advanced techniques for accelerating the solution, including our most recent developments. The computational efficiency of the new methods facilitates both simulations of large systems as well as longer and therefore more realistic simulations of smaller systems.
Collapse
Affiliation(s)
- J Kurzak
- Department of Computer Science, University of Houston, Houston, TX 77204, USA
| | | |
Collapse
|
11
|
Abstract
In the last five years we have witnessed a significant increase in the number publications describing accurate and reliable all-atom molecular dynamics simulations of nucleic acids. This increase has been facilitated by the development of fast and efficient methods for treating the long-range electrostatic interactions, the availability of faster parallel computers, and the development of well-validated empirical molecular mechanical force fields. With these technologies, it has been demonstrated that simulation is not only capable of consistently reproducing experimental observations of sequence specific fine structure of DNA, but also can give detailed insight into prevalent problems in nucleic acid structure, ion association and specific hydration of nucleic acids, polyadenine tract bending, and the subtle environmental dependence of the A-DNA-B-DNA duplex equilibrium. Despite the advances, there are still issues with the methods that need to be resolved through rigorous controlled testing. In general, these relate to deficiencies of the underlying molecular mechanical potentials or applied methods (such as the imposition of true periodicity in Ewald simulations and the need for energy conservation), and significant limits in effective conformational sampling. In this perspective, we provide an overview of our experiences, provide some cautionary notes, and provide recommendations for further study in molecular dynamics simulation of nucleic acids.
Collapse
Affiliation(s)
- T E Cheatham
- Department of Medicinal Chemistry, University of Utah, 30 South, 2000 East, Skaggs Hall 201, Salt Lake City, UT 84112-5820, USA.
| | | |
Collapse
|
12
|
|