A comparison of neighbor search algorithms for large rigid molecules.

Fast determination of neighboring atoms is an essential step in molecular dynamics simulations or Monte Carlo computations, and there exists a variety of algorithms to efficiently compute neighbor lists. However, most of these algorithms are general, and not specifically designed for a given type of application. As a result, although their average performance is satisfactory, they might be inappropriate in some specific application domains. In this article, we study the case of detecting neighbors between large rigid molecules, which has applications in, e.g., rigid body molecular docking, Monte Carlo simulations of molecular self-assembly or diffusion, and rigid body molecular dynamics simulations. More precisely, we compare the traditional grid-based algorithm to a series of hierarchy-based algorithms that use bounding volumes to rapidly eliminate large groups of irrelevant pairs of atoms during the neighbor search. We compare the performance of these algorithms based on several parameters: the size of the molecules, the average distance between them, the cutoff distance, as well as the type of bounding volume used in the culling hierarchy (AABB, OBB, wrapped, or layered spheres). We demonstrate that for relatively large systems (> 100,000 atoms) the algorithm based on the hierarchy of wrapped spheres shows the best results and the traditional grid-based algorithm gives the worst timings. For small systems, however, the grid-based algorithm and the one based on the wrapped sphere hierarchy are beneficial.

A broad assessment of factors determining Culicoides imicola abundance: modelling the present and forecasting its future in climate change scenarios.

Bluetongue (BT) is still present in Europe and the introduction of new serotypes from endemic areas in the African continent is a possible threat. Culicoides imicola remains one of the most relevant BT vectors in Spain and research on the environmental determinants driving its life cycle is key to preventing and controlling BT. Our aim was to improve our understanding of the biotic and abiotic determinants of C. imicola by modelling its present abundance, studying the spatial pattern of predicted abundance in relation to BT outbreaks, and investigating how the predicted current distribution and abundance patterns might change under future (2011-2040) scenarios of climate change according to the Intergovernmental Panel on Climate Change. C. imicola abundance data from the bluetongue national surveillance programme were modelled with spatial, topoclimatic, host and soil factors. The influence of these factors was further assessed by variation partitioning procedures. The predicted abundance of C. imicola was also projected to a future period. Variation partitioning demonstrated that the pure effect of host and topoclimate factors explained a high percentage (>80%) of the variation. The pure effect of soil followed in importance in explaining the abundance of C. imicola. A close link was confirmed between C. imicola abundance and BT outbreaks. To the best of our knowledge, this study is the first to consider wild and domestic hosts in predictive modelling for an arthropod vector. The main findings regarding the near future show that there is no evidence to suggest that there will be an important increase in the distribution range of C. imicola; this contrasts with an expected increase in abundance in the areas where it is already present in mainland Spain. What may be expected regarding the future scenario for orbiviruses in mainland Spain, is that higher predicted C. imicola abundance may significantly change the rate of transmission of orbiviruses.