Item


Conical intersection optimization based on a double Newton-Raphson algorithm using composed steps

An algorithm for conical intersection optimization based on a double Newton-Raphson step (DNR) has been implemented and tested in 11 cases using CASSCF as the electronic structure method. The optimization is carried out in redundant coordinates, and the steps are the sum of two independent Newton-Raphson steps. The first step is carried out to reach the energy degeneracy and uses the gradient of the energy difference between the crossing states and the so-called branching space Hessian. The second step minimizes the energy in the intersection space and uses the projected excited state gradient and the intersection space Hessian. The branching and intersection space Hessians are obtained with a Broyden-Fletcher-Goldfarb-Shanno update from the gradient difference and projected excited state gradients, respectively. In some cases, mixing of the quasi-degenerate states near the seam causes changes in the direction of the gradient difference vector and induces a loss of the degeneracy. This behavior is avoided switching to a composed step (CS) algorithm [Sicilia et al. J. Chem. Theory Comput.2008, 4, 27], i.e., a hybrid DNR-CS implementation. Compared to the composed gradient (CG) [Bearpark et al. Chem. Phys. Lett.1994, 223, 269] and hybrid CG-CS algorithms, the DNR-CS algorithm reaches the MECI in 30% and 15% less steps, respectively. The improvement occurs mostly because the approach to the seam is more efficient, and a degeneracy threshold of 0.001 hartree is reached at lower energies than in the CG and CG-CS cases

This work has been supported by grants CTQ2008-06696 and CTQ2011-26573 from the Spanish Ministerio de Ciencia e Innovacion (MICINN) and Ministerio de Economia y Competividad (MINECO), respectively, SGR0528 from the Catalan Agencia de Gestio d’Ajuts Universitaris i de Recerca (AGAUR), UNGI08-4E-003 from MICINN and the European Fund for Regional Development, and the Xarxa de Referencia en Quimica Teorica i Computacional de Catalunya from AGAUR S.R.-B. thanks the MICINN for grant BES-2009-029177

© Journal of Chemical Theory and Computation, 2013, vol. 9, núm. 3, p. 1433-1442

American Chemical Society (ACS)

Author: Ruiz-Barragán, Sergi
Robb, Michael A.
Blancafort San José, Lluís
Date: 2013
Abstract: An algorithm for conical intersection optimization based on a double Newton-Raphson step (DNR) has been implemented and tested in 11 cases using CASSCF as the electronic structure method. The optimization is carried out in redundant coordinates, and the steps are the sum of two independent Newton-Raphson steps. The first step is carried out to reach the energy degeneracy and uses the gradient of the energy difference between the crossing states and the so-called branching space Hessian. The second step minimizes the energy in the intersection space and uses the projected excited state gradient and the intersection space Hessian. The branching and intersection space Hessians are obtained with a Broyden-Fletcher-Goldfarb-Shanno update from the gradient difference and projected excited state gradients, respectively. In some cases, mixing of the quasi-degenerate states near the seam causes changes in the direction of the gradient difference vector and induces a loss of the degeneracy. This behavior is avoided switching to a composed step (CS) algorithm [Sicilia et al. J. Chem. Theory Comput.2008, 4, 27], i.e., a hybrid DNR-CS implementation. Compared to the composed gradient (CG) [Bearpark et al. Chem. Phys. Lett.1994, 223, 269] and hybrid CG-CS algorithms, the DNR-CS algorithm reaches the MECI in 30% and 15% less steps, respectively. The improvement occurs mostly because the approach to the seam is more efficient, and a degeneracy threshold of 0.001 hartree is reached at lower energies than in the CG and CG-CS cases
This work has been supported by grants CTQ2008-06696 and CTQ2011-26573 from the Spanish Ministerio de Ciencia e Innovacion (MICINN) and Ministerio de Economia y Competividad (MINECO), respectively, SGR0528 from the Catalan Agencia de Gestio d’Ajuts Universitaris i de Recerca (AGAUR), UNGI08-4E-003 from MICINN and the European Fund for Regional Development, and the Xarxa de Referencia en Quimica Teorica i Computacional de Catalunya from AGAUR S.R.-B. thanks the MICINN for grant BES-2009-029177
Format: application/pdf
ISSN: 1549-9618 (versió paper)
1549-9626 (versió electrònica)
Document access: http://hdl.handle.net/10256/11474
Language: eng
Publisher: American Chemical Society (ACS)
Collection: MICINN/PN 2012-2014/CTQ2011-26573
MEC/2008/UNGI08-4E-003
MEC/PN 2009-2011/CTQ2008-06696
AGAUR/2009-2014/2009 SGR-528
Reproducció digital del document publicat a: http://dx.doi.org/10.1021/ct301059t
Articles publicats (D-Q)
Is part of: © Journal of Chemical Theory and Computation, 2013, vol. 9, núm. 3, p. 1433-1442
Rights: Tots els drets reservats
Subject: Optimització matemàtica
Mathematical optimization
Química quàntica
Quantum chemistry
Química de l’estat excitat
Excited state chemistry
Title: Conical intersection optimization based on a double Newton-Raphson algorithm using composed steps
Type: info:eu-repo/semantics/article
Repository: DUGiDocs

Subjects

Authors