Comparison of algorithms for conical intersection optimisation using semiempirical methods

We present a comparison of three previously published algorithms for optimising the minimum energy crossing point between two Born–Oppenheimer electronic states. The algorithms are implemented in a development version of the MNDO electronic structure package for use with semiempirical configuration interaction methods. The penalty function method requires only the energies and gradients of the states involved, whereas the gradient projection and Lagrange–Newton methods also require the calculation of non-adiabatic coupling terms. The performance of the algorithms is measured against a set of well-known small molecule conical intersections. The Lagrange–Newton method is found to be the most efficient, with the projected gradient method also competitive. The penalty function method can only be recommended for situations where non-adiabatic coupling terms cannot be calculated.

[1]  Marco Garavelli,et al.  The C 5 H 6 NH 2 + Protonated Shiff Base: An ab Initio Minimal Model for Retinal Photoisomerization , 1997 .

[2]  W. Kabsch.  A solution for the best rotation to relate two sets of vectors , 1976 .

[3]  D. Yarkony,et al.  Nonadiabatic Interactions Between Potential Energy Surfaces: Theory and Applications , 2007 .

[4]  D. Yarkony.  Marching along ridges. Efficient location of energy-minimized conical intersections of two states using extrapolatable functions , 2004 .

[5]  T. Martínez,et al.  Optimization of Conical Intersections with Floating Occupation Semiempirical Configuration Interaction Wave Functions , 2002 .

[6]  M. Klessinger,et al.  Optimization of conical intersections using the semiempirical MNDOC–CI method with analytic gradients , 2000 .

[7]  Bruce A. Murtagh,et al.  Computational Experience with Quadratically Convergent Minimisation Methods , 1970, Comput. J..

[8]  Fernando Bernardi,et al.  Optimization and characterization of the lowest energy point on a conical intersection using an MC-SCF Lagrangian , 1992 .

[9]  J. Paier,et al.  Photochemistry of ethylene: a multireference configuration interaction investigation of the excited-state energy surfaces. , 2004, The Journal of chemical physics.

[10]  W. Kabsch.  A discussion of the solution for the best rotation to relate two sets of vectors , 1978 .

[11]  R. Fletcher.  Practical Methods of Optimization , 1988 .

[12]  Roland Lindh,et al.  New General Tools for Constrained Geometry Optimizations. , 2005, Journal of chemical theory and computation.

[13]  Todd J. Martínez,et al.  Conical intersections and double excitations in time-dependent density functional theory , 2006 .

[14]  D. Yarkony,et al.  Conical Intersections: Electronic Structure, Dynamics and Spectroscopy , 2004 .

[15]  K. Morokuma,et al.  Ab initio study of nonadiabatic interactions in the photodissociation of ketene , 1997 .

[16]  D. Yarkony,et al.  On the intersection of two potential energy surfaces of the same symmetry. Systematic characterization using a Lagrange multiplier constrained procedure , 1993 .

[17]  Michael J. Bearpark,et al.  A direct method for the location of the lowest energy point on a potential surface crossing , 1994 .

[18]  Keiji Morokuma,et al.  Determination of the lowest energy point on the crossing seam between two potential surfaces using the energy gradient , 1985 .

[19]  Thom Vreven,et al.  Photoisomerization Path for a Realistic Retinal Chromophore Model: The Nonatetraeniminium Cation , 1998 .

[20]  F. Bernardi,et al.  Mechanism of Carbene Formation from the Excited-States of Diazirine and Diazomethane - An MC-SCF Study , 1994 .

[21]  Walter Thiel,et al.  Implementation of a general multireference configuration interaction procedure with analytic gradients in a semiempirical context using the graphical unitary group approach , 2003, J. Comput. Chem..

[22]  Potential energy surface crossings in organic photochemistry , 1996 .

[23]  I. N. Ragazos,et al.  A conical intersection mechanism for the photochemistry of butadiene. A MC-SCF study , 1993 .

[24]  Walter Thiel,et al.  Generic implementation of semi-analytical CI gradients for NDDO-type methods , 2005 .

[25]  Michael J. Bearpark,et al.  Characterization of the indacene S0/S1 conical intersection: an MMVB and CASSCF study , 1999 .

[26]  Walter Thiel,et al.  Orthogonalization corrections for semiempirical methods , 2000 .

[27]  Thom Vreven,et al.  The Azulene S1 State Decays via a Conical Intersection: A CASSCF Study with MMVB Dynamics , 1996 .

[28]  Michel Dupuis,et al.  On the determination of the minimum on the crossing seam of two potential energy surfaces , 1991 .

[29]  G. Granucci,et al.  The photoisomerization mechanism of azobenzene: a semiclassical simulation of nonadiabatic dynamics. , 2004, Chemistry.

[30]  M. Barbatti,et al.  Ultrafast two-step process in the non-adiabatic relaxation of the CH2 molecule , 2006 .

[31]  R A Mathies,et al.  The first step in vision: femtosecond isomerization of rhodopsin. , 1991, Science.

[32]  Josep Maria Bofill,et al.  A reduced‐restricted‐quasi‐Newton–Raphson method for locating and optimizing energy crossing points between two potential energy surfaces , 1997 .

[33]  Ian J. Palmer,et al.  An MC-SCF Study of the S1 and S2 Photochemical Reactions of Benzene , 1993 .

[34]  Massimo Olivucci,et al.  Ground and excited state CASPT2 geometry optimizations of small organic molecules , 2003, J. Comput. Chem..