From Atoms to Fullerene: Stochastic Surface Walking Solution for Automated Structure Prediction of Complex Material.

It is of general concern whether the automated structure prediction of unknown material without recourse to any knowledge from experiment is ever possible considering the daunting complexity of potential energy surface (PES) of material. Here we demonstrate that the stochastic surface walking (SSW) method can be a general and promising solution to this ultimate goal, which is applied to assemble carbon fullerenes containing up to 100 atoms (including 60, 70, 76, 78, 80, 84, 90, 96, and 100 atoms) from randomly distributed atoms, a long-standing challenge in global optimization. Combining the SSW method with a parallel replica exchange algorithm, we can locate the global minima (GM) of these large fullerenes efficiently without being trapped in numerous energy-nearly degenerate isomers. Detailed analyses on the SSW trajectories allow us to rationalize how and why the SSW method is able to explore the highly complex PES, which highlights the abilities of SSW method for surmounting the high barrier and the preference of SSW trajectories to the low energy pathways. The work demonstrates that the parallel SSW method is a practical tool for predicting unknown materials.

[1]  STRUCTURAL RELAXATION IN MORSE CLUSTERS : ENERGY LANDSCAPES , 1998, cond-mat/9808080.

[2]  S. Maruyama,et al.  A molecular dynamics demonstration of annealing to a perfect C60 structure , 1998 .

[3]  Zhi-Pan Liu,et al.  Constrained Broyden Dimer Method with Bias Potential for Exploring Potential Energy Surface of Multistep Reaction Process. , 2012, Journal of chemical theory and computation.

[4]  Tiffany R. Walsh,et al.  Relaxation dynamics of C60 , 1998 .

[5]  D. Wales,et al.  Theoretical studies of icosahedral C60 and some related species , 1986 .

[6]  J. Gale,et al.  The prediction of inorganic crystal structures using a genetic algorithm and energy minimisation , 1999 .

[7]  C. Shang,et al.  Constrained Broyden Minimization Combined with the Dimer Method for Locating Transition State of Complex Reactions , 2010 .

[8]  Zhipan Liu,et al.  Comprehensive mechanism and structure-sensitivity of ethanol oxidation on platinum: new transition-state searching method for resolving the complex reaction network. , 2008, Journal of the American Chemical Society.

[9]  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.

[10]  David J. Wales,et al.  Ab initio study of rearrangements between C60 fullerenes , 2003 .

[11]  Ho,et al.  Molecular geometry optimization with a genetic algorithm. , 1995, Physical review letters.

[12]  Che Ting Chan,et al.  The geometry of small fullerene cages: C20 to C70 , 1992 .

[13]  S. Maruyama,et al.  A molecular dynamics simulation of the fullerene formation process , 1998 .

[14]  Zhi-Pan Liu,et al.  Stochastic Surface Walking Method for Structure Prediction and Pathway Searching. , 2013, Journal of chemical theory and computation.

[15]  Stefan Goedecker,et al.  Energy landscape of fullerene materials: a comparison of boron to boron nitride and carbon. , 2010, Physical review letters.

[16]  Michael W Deem,et al.  Parallel tempering: theory, applications, and new perspectives. , 2005, Physical chemistry chemical physics : PCCP.

[17]  Wang,et al.  Replica Monte Carlo simulation of spin glasses. , 1986, Physical review letters.

[18]  C. Shang,et al.  Origin and activity of gold nanoparticles as aerobic oxidation catalysts in aqueous solution. , 2011, Journal of the American Chemical Society.

[19]  H. Scheraga,et al.  Global optimization of clusters, crystals, and biomolecules. , 1999, Science.

[20]  Jonathan Doye,et al.  Thermodynamics of Global Optimization , 1998 .

[21]  David J. Wales,et al.  Some further applications of discrete path sampling to cluster isomerization , 2004 .

[22]  Mark A. Miller,et al.  Archetypal energy landscapes , 1998, Nature.

[23]  E. Ōsawa,et al.  Computational modeling of the elemental catalysis in the Stone–Wales fullerene rearrangements , 2000 .

[24]  Alessandro Laio,et al.  Efficient exploration of reactive potential energy surfaces using Car-Parrinello molecular dynamics. , 2003, Physical review letters.

[25]  Craig E. Larson,et al.  Graph-Theoretic Independence as a Predictor of Fullerene Stability , 2003 .

[26]  J. Doye,et al.  Tetrahedral global minimum for the 98-atom Lennard-Jones cluster. , 1999, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[27]  Donald W. Brenner,et al.  A second-generation reactive empirical bond order (REBO) potential energy expression for hydrocarbons , 2002 .

[28]  Y. Sugita,et al.  Replica-exchange molecular dynamics method for protein folding , 1999 .

[29]  Stephan Irle,et al.  The C60 formation puzzle "solved": QM/MD simulations reveal the shrinking hot giant road of the dynamic fullerene self-assembly mechanism. , 2006, The journal of physical chemistry. B.

[30]  Julian D. Gale,et al.  The General Utility Lattice Program (GULP) , 2003 .

[31]  Hanchul Kim,et al.  Dynamic pathway model for the formation of C(60). , 2004, The Journal of chemical physics.

[32]  Fabio Schoen,et al.  Global Optimization of Morse Clusters by Potential Energy Transformations , 2004, INFORMS J. Comput..

[33]  Steven Hobday,et al.  Optimisation of carbon cluster geometry using a genetic algorithm , 1997 .

[34]  David J Wales,et al.  Equilibrium thermodynamics from basin-sampling. , 2006, The Journal of chemical physics.

[35]  C. D. Gelatt,et al.  Optimization by Simulated Annealing , 1983, Science.

[36]  Milani,et al.  Simulated annealing of carbon clusters. , 1990, Physical review. B, Condensed matter.

[37]  A. Oganov,et al.  Crystal structure prediction using ab initio evolutionary techniques: principles and applications. , 2006, The Journal of chemical physics.

[38]  Reduction of activation energy barrier of Stone-Wales transformation in endohedral metallofullerenes , 2005, cond-mat/0507402.

[39]  Julian D. Gale,et al.  Empirical potential derivation for ionic materials , 1996 .

[40]  H. W. Kroto,et al.  The stability of the fullerenes Cn, with n = 24, 28, 32, 36, 50, 60 and 70 , 1987, Nature.

[41]  J. Doye,et al.  The effect of the range of the potential on the structures of clusters , 1995 .

[42]  Sándor Suhai,et al.  Self-consistent-charge density-functional tight-binding method for simulations of complex materials properties , 1998 .

[43]  S. C. O'brien,et al.  C60: Buckminsterfullerene , 1985, Nature.

[44]  J. Grossman,et al.  Structure and stability of molecular carbon: Importance of electron correlation. , 1995, Physical review letters.

[45]  Zhipan Liu,et al.  Mechanism and Tafel lines of electro-oxidation of water to oxygen on RuO2(110). , 2010, Journal of the American Chemical Society.

[46]  Jonathan P. K. Doye,et al.  Stationary points and dynamics in high-dimensional systems , 2003 .

[47]  T. Frauenheim,et al.  DFTB+, a sparse matrix-based implementation of the DFTB method. , 2007, The journal of physical chemistry. A.

[48]  F. Stillinger,et al.  Packing Structures and Transitions in Liquids and Solids , 1984, Science.

[49]  Xu,et al.  Tight-binding molecular dynamics simulations of fullerene annealing and fragmentation. , 1994, Physical review letters.

[50]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[51]  Modelling diamond-like carbon with the environment-dependent interaction potential , 2002 .

[52]  A. V. Duin,et al.  ReaxFF: A Reactive Force Field for Hydrocarbons , 2001 .

[53]  Jinlong Yang,et al.  Global minimum structures of morse clusters as a function of the range of the potential: 81, 2007, The journal of physical chemistry. A.

[54]  Douglas J. Klein,et al.  GENERALIZATIONS OF THE STONE-WALES REARRANGEMENT FOR CAGE COMPOUNDS, INCLUDING FULLERENES , 1996 .