Size‐guided multi‐seed heuristic method for geometry optimization of clusters: Application to benzene clusters

Since searching for the global minimum on the potential energy surface of a cluster is very difficult, many geometry optimization methods have been proposed, in which initial geometries are randomly generated and subsequently improved with different algorithms. In this study, a size‐guided multi‐seed heuristic method is developed and applied to benzene clusters. It produces initial configurations of the cluster with n molecules from the lowest‐energy configurations of the cluster with n − 1 molecules (seeds). The initial geometries are further optimized with the geometrical perturbations previously used for molecular clusters. These steps are repeated until the size n satisfies a predefined one. The method locates putative global minima of benzene clusters with up to 65 molecules. The performance of the method is discussed using the computational cost, rates to locate the global minima, and energies of initial geometries. © 2018 Wiley Periodicals, Inc.

[1]  Hiroshi Takeuchi Two Perturbations for Geometry Optimization of Off‐lattice Bead Protein Models , 2017, Molecular informatics.

[2]  Faming Liang,et al.  Annealing contour Monte Carlo algorithm for structure optimization in an off-lattice protein model. , 2004, The Journal of chemical physics.

[3]  Johannes M. Dieterich,et al.  Composition‐induced structural transitions in mixed Lennard‐Jones clusters: Global reparametrization and optimization , 2011, J. Comput. Chem..

[4]  István Kolossváry,et al.  Global optimization of additive potential energy functions: predicting binary Lennard-Jones clusters. , 2010, Physical review. E, Statistical, nonlinear, and soft matter physics.

[5]  Jorge Nocedal,et al.  On the limited memory BFGS method for large scale optimization , 1989, Math. Program..

[6]  Bernd Hartke,et al.  Global Geometry Optimization of Molecular Clusters: TIP4P Water , 2000 .

[7]  Bernd Hartke,et al.  Conquering the hard cases of Lennard-Jones clusters with simple recipes , 2017 .

[8]  J. Straub,et al.  Statistical-temperature Monte Carlo and molecular dynamics algorithms. , 2006, Physical review letters.

[9]  Howard R. Mayne,et al.  A study of genetic algorithm approaches to global geometry optimization of aromatic hydrocarbon microclusters , 1998 .

[10]  Bernd Hartke,et al.  Larger water clusters with edges and corners on their way to ice: structural trends elucidated with an improved parallel evolutionary algorithm. , 2006, The journal of physical chemistry. A.

[11]  David J. Wales,et al.  Global minima of water clusters (H2O)n, n≤21, described by an empirical potential , 1998 .

[12]  Hiroshi Takeuchi,et al.  Novel Method for Geometry Optimization of Molecular Clusters: Application to Benzene Clusters , 2007, J. Chem. Inf. Model..

[13]  B. W. V. D. Waal,et al.  Computed structure of small benzene clusters , 1986 .

[14]  Wolfhard Janke,et al.  Multicanonical study of coarse-grained off-lattice models for folding heteropolymers. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

[15]  Hiroshi Takeuchi,et al.  Clever and Efficient Method for Searching Optimal Geometries of Lennard-Jones Clusters , 2006, J. Chem. Inf. Model..

[16]  Hiroshi Takeuchi,et al.  Ternary and quaternary Lennard-Jones atomic clusters: The effects of atomic sizes on the compositions, geometries, and relative stability , 2015 .

[17]  P. E. Abreu,et al.  A global optimization perspective on molecular clusters , 2017, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[18]  Thomas L. Starr,et al.  Calculation of the crystal structures of hydrocarbons by molecular packing analysis , 1977, Comput. Chem..

[19]  Mao Chen,et al.  Heuristic algorithm for off-lattice protein folding problem , 2006, Journal of Zhejiang University SCIENCE B.

[20]  Jooyoung Lee,et al.  Structure optimization by conformational space annealing in an off-lattice protein model. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

[21]  Hiroshi Takeuchi,et al.  Structural features of small benzene clusters (C6H6)n (n ≤ 30) as investigated with the all-atom OPLS potential. , 2012, The journal of physical chemistry. A.

[22]  K. Ho,et al.  Structural optimization of Lennard-Jones clusters by a genetic algorithm , 1996 .

[23]  Francisco B. Pereira,et al.  A study on diversity for cluster geometry optimization , 2009, Evol. Intell..

[24]  Stefan Goedecker,et al.  Efficient moves for global geometry optimization methods and their application to binary systems. , 2010, The Journal of chemical physics.

[25]  Tao Ye,et al.  Global Optimization of Binary Lennard-Jones Clusters Using Three Perturbation Operators , 2011, J. Chem. Inf. Model..

[26]  Ajit J. Thakkar,et al.  Improved minima-hopping. TIP4P water clusters, (H2O)n with n⩽37 , 2009 .

[27]  Clifford E. Dykstra,et al.  A Model Study of Aggregation of Acetylene Molecules , 2000 .

[28]  Stefan Goedecker,et al.  The performance of minima hopping and evolutionary algorithms for cluster structure prediction. , 2009, The Journal of chemical physics.

[29]  Wayne J. Pullan,et al.  Structure Prediction of Benzene Clusters Using a Genetic Algorithm , 1997, J. Chem. Inf. Comput. Sci..

[30]  Yan Sun,et al.  Structural distribution in mixed ternary noble gas and Lennard-Jones clusters , 2013 .

[31]  Lucjan Piela,et al.  Theoretically predicted lowest-energy structures of water clusters , 1992 .

[32]  Yan Sun,et al.  Structural transitions in mixed ternary noble gas clusters , 2013, Journal of Molecular Modeling.

[33]  Kenneth D. Jordan,et al.  Theoretical study of small water clusters : low-energy fused cubic structures for (H2O)n, n = 8, 12, 16, and 20 , 1993 .

[34]  Bernd Hartke,et al.  Global optimization , 2011 .

[35]  Sergey Kazachenko,et al.  Water nanodroplets: predictions of five model potentials. , 2013, The Journal of chemical physics.

[36]  Fernando Pirani,et al.  Low‐energy structures of benzene clusters with a novel accurate potential surface , 2015, J. Comput. Chem..

[37]  Boon K. Teo,et al.  A Perspective on the Science of Clusters , 2013, Journal of Cluster Science.

[38]  S. Kazachenko,et al.  Are there any magic numbers for water nanodroplets, (H2O) n , in the range 36 ≤ n ≤ 50? , 2010 .

[39]  David J. Wales,et al.  Global minima for water clusters (H2O)n, n ⩽ 21, described by a five-site empirical potential , 2005 .

[40]  Jingfa Liu,et al.  Structure optimization in a three‐dimensional off‐lattice protein model , 2006, Biopolymers.

[41]  Jianpeng Ma,et al.  Simulation via direct computation of partition functions. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[42]  Lars Meyer,et al.  Mapping the magic numbers in binary Lennard-Jones clusters. , 2005, Physical review letters.

[43]  Francisco B. Pereira,et al.  An evolutionary algorithm for global minimum search of binary atomic clusters , 2010 .

[44]  Hiroshi Takeuchi,et al.  Geometry optimization of carbon dioxide clusters (CO2)n for 4 < or = n < or = 40. , 2008, The journal of physical chemistry. A.

[45]  John Arthur Niesse,et al.  Global optimization of atomic and molecular clusters using the space‐fixed modified genetic algorithm method , 1997 .

[46]  Francisco B. Pereira,et al.  An evolutionary algorithm for the global optimization of molecular clusters: application to water, benzene, and benzene cation. , 2011, The journal of physical chemistry. A.

[47]  J. Doye,et al.  Global Optimization by Basin-Hopping and the Lowest Energy Structures of Lennard-Jones Clusters Containing up to 110 Atoms , 1997, cond-mat/9803344.

[48]  V. Buch,et al.  Search for Low Energy Structures of Water Clusters (H2O)n, n = 20−22, 48, 123, and 293 , 2003 .

[49]  D. Vernon Inform , 1995, Encyclopedia of the UN Sustainable Development Goals.

[50]  Bernd Hartke,et al.  Structural transitions in clusters. , 2002, Angewandte Chemie.

[51]  Frank Potthast,et al.  Local interactions and protein folding: A three-dimensional off-lattice approach , 1997 .

[52]  Xiaoli Yang,et al.  A dynamic lattice searching method with constructed core for optimization of large lennard‐jones clusters , 2007, J. Comput. Chem..

[53]  Hiroshi Takeuchi,et al.  Binary Lennard-Jones atomic clusters: Structural features induced by large-sized atoms , 2014 .

[54]  Seung-Yeon Kim,et al.  Re‐examination of structure optimization of off‐lattice protein AB models by conformational space annealing , 2008, J. Comput. Chem..

[55]  Fabio Schoen,et al.  Global optimization of binary Lennard–Jones clusters , 2009, Optim. Methods Softw..

[56]  Reinhard Hentschke,et al.  Global Potential Energy Minima of SPC/E Water Clusters without and with Induced Polarization Using a Genetic Algorithm , 1999 .

[57]  R. Hentschke,et al.  Global Minima of Water Clusters (H2O)N, N ≤ 25, Described by Three Empirical Potentials , 2003 .

[58]  Francisco B. Pereira,et al.  Colloidal clusters from a global optimization perspective , 2015 .

[59]  D. E. Williams,et al.  Calculated energy and conformation of clusters of benzene molecules and their relationship to crystalline benzene , 1980 .

[60]  J. Marques,et al.  On the use of different potential energy functions in rare-gas cluster optimization by genetic algorithms: application to argon clusters. , 2008, The journal of physical chemistry. A.

[61]  Xueguang Shao,et al.  A dynamic lattice searching method for fast optimization of Lennard–Jones clusters , 2004, J. Comput. Chem..

[62]  Hiroshi Takeuchi Global minimum geometries of acetylene clusters (HCCH)n with n ≤ 55 obtained by a heuristic method combined with geometrical perturbations , 2010, J. Comput. Chem..

[63]  Wolfhard Janke,et al.  Two-state folding, folding through intermediates, and metastability in a minimalistic hydrophobic-polar model for proteins. , 2007, Physical review letters.

[64]  David C. Easter Toward the Experimental Structure Determination of Larger Molecular Clusters: Application and Limitations of the Weak Interaction Model to the (C6H6)13B2u ← A1g000 Spectrum , 2004 .

[65]  John E Straub,et al.  Structure optimization and folding mechanisms of off-lattice protein models using statistical temperature molecular dynamics simulation: Statistical temperature annealing. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[66]  Bernd Hartke,et al.  Size-dependent transition from all-surface to interior-molecule structures in pure neutral water clusters , 2003 .

[67]  Haiyan Jiang,et al.  An Efficient Method Based on Lattice Construction and the Genetic Algorithm for Optimization of Large Lennard-Jones Clusters. , 2004, The journal of physical chemistry. A.

[68]  Gustavo G. Rondina,et al.  Revised Basin-Hopping Monte Carlo Algorithm for Structure Optimization of Clusters and Nanoparticles , 2013, J. Chem. Inf. Model..

[69]  Xueguang Shao,et al.  Structural Distribution of Lennard-Jones Clusters Containing 562 to 1000 Atoms. , 2004, The journal of physical chemistry. A.

[70]  Tirtharaj Dash,et al.  Gradient gravitational search: An efficient metaheuristic algorithm for global optimization , 2015, J. Comput. Chem..

[71]  H. Kabrede,et al.  Using vibrational modes in the search for global minima of atomic and molecular clusters , 2006 .

[72]  Yan Sun,et al.  Parametric effects of the potential energy function on the geometrical features of ternary Lennard-Jones clusters. , 2012, The journal of physical chemistry. A.

[73]  Hsiao-Ping Hsu,et al.  Structure optimization in an off-lattice protein model. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[74]  Yanchao Wang,et al.  Particle-swarm structure prediction on clusters. , 2012, The Journal of chemical physics.

[75]  C. Zhang,et al.  Enhanced sampling in generalized ensemble with large gap of sampling parameter: case study in temperature space random walk. , 2009, The Journal of chemical physics.

[76]  R. Johnston Evolving better nanoparticles: Genetic algorithms for optimising cluster geometries , 2003 .

[77]  Adam Liwo,et al.  An Efficient Deformation-Based Global Optimization Method (Self-Consistent Basin-to-Deformed-Basin Mapping (SCBDBM)). Application to Lennard-Jones Atomic Clusters , 1999 .

[78]  Roy L. Johnston,et al.  Global optimization analysis of water clusters (H2O)n (11⩽n⩽13) through a genetic evolutionary approach , 2002 .

[79]  Xiaoli Yang,et al.  A dynamic lattice searching method with interior operation for unbiased optimization of large Lennard‐Jones clusters , 2008, J. Comput. Chem..

[80]  William L. Jorgensen,et al.  Aromatic-aromatic interactions: free energy profiles for the benzene dimer in water, chloroform, and liquid benzene , 1990 .

[81]  Handan Arkin Determination of the structure of the energy landscape for coarse-grained off-lattice models of folding heteropolymers. , 2008, Physical review. E, Statistical, nonlinear, and soft matter physics.

[82]  Max L. Berkowitz,et al.  Cube to cage transitions in (H2O)n (n=12, 16, and 20) , 1996 .