PARALLEL ALGORITHM FOR EFFICIENT CALCULATION OF 2ND DERIVATIVES OF CONFORMATIONAL ENERGY FUNCTION IN INTERNAL COORDINATES

A parallel algorithm for efficient calculation of the second derivativ es (Hessian) of the conformational energy in internal coordinates is p roposed, This parallel algorithm is based on the master/slave model. A master processor distributes the calculations of components of the He ssian to one or more slave processors that, after finishing their calc ulations, send the results to the master processor that assembles all the components of the Hessian. Our previously developed molecular anal ysis system for conformational energy optimization, norm,al mode analy sis, and Monte Carlo simulation for internal coordinates is extended t o use this parallel algorithm for Hessian calculation on a massively p arallel computer. The implementation of our algorithm uses the message passing Interface and works effectively on both distributed-memory pa rallel computers and shared-memory parallel computers. We applied this system to the Newton-Raphson energy optimization of the structures of glutaminyl transfer RNA (Gln-tRNA) with 74 nucleotides and glutaminyl -tRNA synthetase (GlnRS) with 540 residues to analyze the performance of our system. The parallel speedups for the Hessian calculation were 6.8 for Gln-tRNA with 24 processors and 11.2 for GlnRS with 54 process ors. The parallel speedups for the Newton-Raphson optimization were 6. 3 for Gln-tRNA with 30 processors and 12.0 for GlnRS with 62 processor s. (C) 1998 John Wiley & Sons, Inc.