Parallel molecular dynamics simulation code for dialkyl cationic surfactants

We have developed a new simulation code, COMFORT, for the study of assemblies of flexible surfactant molecules, structured for parallel execution and specialized to surfactants with dialkyl chain geometry. The approach is a hybrid domain-decomposition and systolic-loop algorithm which is suitable for systems composed of long chain molecules and with tens of thousands of atoms in total. The algorithm uses a modified Ewald technique for two dimensionally periodic systems which has been successfully parallelized. The code was designed to be highly portable between machines of different architectures. The code has been tested on a number of platforms including the Intel iPSC/860, the IBM SP1, the CRAY T3D, a SGI Power Challenge and a number of workstation clusters. We demonstrate that scalable parallel computing technology, combined with appropriate software, can provide a commercially viable simulation system for use in the exploration and development of surfactant assemblies.