Optimization of Semiempirical Quantum Chemistry Methods via Multiobjective Genetic Algorithms: Accurate Photodynamics for Larger Molecules and Longer Time Scales

Excited-state photodynamics is important in numerous varieties of important materials applications (e.g., liquid crystal display, light emitting diode), pharmaceuticals, and chemical manufacturing processing. We study the effectiveness of multiobjective genetic and evolutionary algorithms in multiscaling excited-state direct photodynamics via rapid reparameterization of semiempirical methods. Using a very limited set of ab initio and experimental data, semiempirical parameters are reoptimized to provide globally accurate potential energy surfaces, thereby eliminating the need for expensive ab initio dynamics simulations. Through reoptimization, excited-state energetics are predicted accurately via semiempirical methods, while retaining accurate ground-state predictions. In our initial study of small photo-excited molecules, our results show that the multiobjective evolutionary algorithm consistently yields solutions that are significantly better—up to 384% lower error in the energy and 87% lower error in the energy-gradient—than those reported previously. As verified with direct quantum dynamical calculations, multiple high-quality parameter sets obtained via genetic algorithms show near-ideal behavior on critical and untested excited-state geometries. The results demonstrate that the reparameterization via evolutionary algorithms is a promising way to extend direct dynamics simulations of photochemistry to multi-picosecond time scales and to larger molecules, with promise in more application beyond simple molecular chemistry.

[1]  Martin Pelikan,et al.  Scalable Optimization via Probabilistic Modeling: From Algorithms to Applications (Studies in Computational Intelligence) , 2006 .

[2]  Gary B. Lamont,et al.  Evolutionary Algorithms for Solving Multi-Objective Problems , 2002, Genetic Algorithms and Evolutionary Computation.

[3]  Ingolf V. Hertel,et al.  Internal conversion in highly excited benzene and benzene dimer: femtosecond time-resolved photoelectron spectroscopy , 1997 .

[4]  A. L. Thompson,et al.  Excited state direct dynamics of benzene with reparameterized multi-reference semiempirical configuration interaction methods , 2004 .

[5]  Samir W. Mahfoud Population Size and Genetic Drift in Fitness Sharing , 1994, FOGA.

[6]  Martin Pelikan,et al.  Hierarchical Bayesian optimization algorithm: toward a new generation of evolutionary algorithms , 2010, SICE 2003 Annual Conference (IEEE Cat. No.03TH8734).

[7]  John R. Koza,et al.  Genetic Programming IV: Routine Human-Competitive Machine Intelligence , 2003 .

[8]  Giovanni Granucci,et al.  Molecular gradients for semiempirical CI wavefunctions with floating occupation molecular orbitals , 2000 .

[9]  D. Goldberg,et al.  Modeling tournament selection with replacement using apparent added noise , 2001 .

[10]  Kalyanmoy Deb,et al.  Simulated Binary Crossover for Continuous Search Space , 1995, Complex Syst..

[11]  Kalyanmoy Deb,et al.  Messy Genetic Algorithms: Motivation, Analysis, and First Results , 1989, Complex Syst..

[12]  David E. Goldberg,et al.  Genetic Algorithms in Search Optimization and Machine Learning , 1988 .

[13]  Todd J. Martínez,et al.  Ab Initio Quantum Molecular Dynamics , 2002 .

[14]  David E. Goldberg,et al.  Genetic Programming for Multiscale Modeling , 2004 .

[15]  Kalyanmoy Deb,et al.  Real-coded Genetic Algorithms with Simulated Binary Crossover: Studies on Multimodal and Multiobjective Problems , 1995, Complex Syst..

[16]  Kalyanmoy Deb,et al.  Genetic Algorithms, Noise, and the Sizing of Populations , 1992, Complex Syst..

[17]  J. Stewart Optimization of parameters for semiempirical methods I. Method , 1989 .

[18]  M. Dewar,et al.  Ground States of Molecules. 38. The MNDO Method. Approximations and Parameters , 1977 .

[19]  Jane Michelle Owens Theoretical Studies of the Solvation, Dynamics, and Photochemistry of Ethylene, Retinal Protonated Schiff Base, Oligocellulose, and Gd(III) Clusters , 2004 .

[20]  David E. Goldberg,et al.  Genetic programming for multitimescale modeling , 2005 .

[21]  Todd J. Martínez,et al.  Photochemistry from first principles and direct dynamics , 2005 .

[22]  R. K. Ursem Multi-objective Optimization using Evolutionary Algorithms , 2009 .

[23]  Martin Pelikan,et al.  Scalable Optimization via Probabilistic Modeling , 2006, Studies in Computational Intelligence.

[24]  John R. Koza,et al.  Genetic programming - on the programming of computers by means of natural selection , 1993, Complex adaptive systems.

[25]  Todd J. Martínez,et al.  Photodynamics of ethylene: ab initio studies of conical intersections , 2000 .

[26]  Eamonn F. Healy,et al.  Development and use of quantum mechanical molecular models. 76. AM1: a new general purpose quantum mechanical molecular model , 1985 .

[27]  John H. Holland,et al.  Adaptation in Natural and Artificial Systems: An Introductory Analysis with Applications to Biology, Control, and Artificial Intelligence , 1992 .

[28]  T. Martínez,et al.  Ab Initio Multiple Spawning: Photochemistry from First Principles Quantum Molecular Dynamics , 2000 .

[29]  David E. Goldberg,et al.  The Design of Innovation: Lessons from and for Competent Genetic Algorithms , 2002 .

[30]  David E. Goldberg,et al.  The Gambler's Ruin Problem, Genetic Algorithms, and the Sizing of Populations , 1999, Evolutionary Computation.

[31]  Kalyanmoy Deb,et al.  A fast and elitist multiobjective genetic algorithm: NSGA-II , 2002, IEEE Trans. Evol. Comput..