Computer Physics Communications, Vol. 128 (1-2) (2000) pp. 399-411
© 2000 Elsevier Science B.V. All rights reserved.
PII: S0010-4655(99)00515-9

Efficient parallel algorithms in global optimization of potential energy functions for peptides, proteins, and crystals ¤

Jooyoung Lee a, Jaroslaw Pillardy a, Cezary Czaplewski a, Yelena Arnautova a, Daniel R. Ripoll b, Adam Liwo a,c, Kenneth D. Gibson a, Ryszard J. Wawak a and Harold A. Scheraga a * has5@cornell.edu

a Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, NY 14853-1301, USA
b Cornell Theory Center, Ithaca, NY 14853-3801, USA
c Faculty of Chemistry, University of Gdansk, Sobieskiego 18, 80-952 Gdansk, Poland

Abstract

Global optimization is playing an increasing role in physics, chemistry, and biophysical chemistry. One of the most important applications of global optimization is to find the global minima of the potential energy of molecules or molecular assemblies, such as crystals. The solution of this problem typically requires huge computational effort. Even the fastest processor available is not fast enough to carry out this kind of computation in real time for the problems of real interest, e.g., protein and crystal structure prediction. One way to circumvent this problem is to take advantage of massively parallel computing. In this paper, we provide several examples of parallel implementations of global optimization algorithms developed in our laboratory. All of these examples follow the master/worker approach. Most of the methods are parallelized on the algorithmic (coarse-grain) level and one example of fine-grain parallelism is given, in which the function evaluation itself is computationally expensive. All parallel algorithms were initially implemented on an IBM/SP2 (distributed-memory) machine. In all cases, however, message passing is handled through the standard Message Passing Interface (MPI); consequently the algorithms can also be implemented on any distributed- or shared-memory system that runs MPI. The efficiency of these implementations is discussed.

PACS: 87.15.Aa; 02.70.-c; 02.60.Pn

Keywords: Crystal structure prediction; Genetic algorithms; Global optimization; Monte Carlo methods; Parallel algorithms; Protein structure prediction

*Corresponding author.

 


¤This paper is published as part of a thematic issue on Parallel Computing in Chemical Physics.