CellListMap.jl: Efficient and customizable cell list implementation for calculation of pairwise particle properties within a cutoff

Autor(es)

Martinez, Leandro

• Tipo de documento: Artigo Científico
• Fonte: COMPUTER PHYSICS COMMUNICATIONS; v. 279, p. 12-pg., 2022-10-01.

Resumo

N-body simulations and trajectory analysis rely on the calculation of attributes that depend on pairwise particle distances within a cutoff. Interparticle potential energies, forces, distribution functions, neighbor lists, and distance-dependent distributions, for example, must be calculated. Cell lists are widely used to avoid computing distances outside the cutoff. However, efficient cell list implementations are difficult to customize. Here, we provide a fast and parallel implementation of cell lists in Julia that allows the mapping of custom functions dependent on particle positions in 2 or 3 dimensions. Arbitrary periodic boundary conditions are supported. Automatic differentiation and unit propagation can be used. The implementation provides a framework for the development of new analysis tools and simulations with custom potentials. The performance of resulting computations is comparable to state-of-the-art implementations of neighbor list algorithms and cell lists, available in specialized software. Examples are provided for the computation of potential energies, forces, distribution of pairwise velocities, neighbor lists and other typical calculations in molecular and astrophysical simulations. The Julia package is freely available at http://m3g .github .io /CellListMap .jl. Interfacing with Python and R with minimal overhead is possible. (C) 2022 Elsevier B.V. All rights reserved. (AU)

• Processo FAPESP: 10/16947-9 - Estrutura, dinâmica e função em proteínas: simulação computacional e desenvolvimento de algoritmos
• Beneficiário: Leandro Martinez
• Modalidade de apoio: Auxílio à Pesquisa - Regular
• Processo FAPESP: 13/08293-7 - CECC - Centro de Engenharia e Ciências Computacionais
• Beneficiário: Munir Salomao Skaf
• Modalidade de apoio: Auxílio à Pesquisa - Centros de Pesquisa, Inovação e Difusão - CEPIDs
• Processo FAPESP: 18/14274-9 - Determinação de estruturas de proteínas usando restrições de distância obtidas de experimentos de ligação cruzada: métodos computacionais e aplicações
• Beneficiário: Leandro Martinez
• Modalidade de apoio: Auxílio à Pesquisa - Regular
• Processo FAPESP: 18/24293-0 - Métodos computacionais de otimização
• Beneficiário: Sandra Augusta Santos
• Modalidade de apoio: Auxílio à Pesquisa - Temático