The problem of global geometry optimization of clusters is addressed with a phenotype variant of the method of genetic algorithms, with several novel performance enhancements. The resulting algorithm is applied to Lennard–Jones clusters as benchmark system, with up to 150 atoms. The well‐known, difficult cases involving nonicosahedral global minima can be treated reliably using the concept of niches. The scaling of computer time with cluster size is approximately cubic, which is crucial for future applications to much larger clusters. © 1999 John Wiley & Sons, Inc. J Comput Chem 20: 1752–1759, 1999