We have used state-of-the-art ab initio restricted active RASPT2 computations using a 16 orbitals, 18 electrons active space to produce an extended three-dimensional map of the potential energy surfaces (PESs) of the ground and first nπ∗ excited states of azobenzene along CNNC torsion and the two CNN bending angles, which are the most relevant coordinates for the trans-cis photoisomerization process. Through comparison with fully unconstrained optimizations performed at the same level of theory, we show that the three selected coordinates suffice to correctly describe the photoisomerization mechanism and the S1-S0 crossing seam. We also provide a map of the nonadiabatic coupling between the two states in the region where they get closer in energy. Eventually, we show that treating the two CNN bending angles as independent coordinates is fundamental to break the symmetry and couple the two electronic states. The accuracy of the S0 and S1 PESs and couplings was validated with semiclassical dynamics simulations in the reduced space of the scanned coordinates, showing results in good agreement with published full-coordinate dynamics.

Multidimensional Potential Energy Surfaces Resolved at the RASPT2 Level for Accurate Photoinduced Isomerization Dynamics of Azobenzene

Aleotti F.;Soprani L.;Nenov A.;Berardi R.;Arcioni A.;Zannoni C.;Garavelli M.
2019

Abstract

We have used state-of-the-art ab initio restricted active RASPT2 computations using a 16 orbitals, 18 electrons active space to produce an extended three-dimensional map of the potential energy surfaces (PESs) of the ground and first nπ∗ excited states of azobenzene along CNNC torsion and the two CNN bending angles, which are the most relevant coordinates for the trans-cis photoisomerization process. Through comparison with fully unconstrained optimizations performed at the same level of theory, we show that the three selected coordinates suffice to correctly describe the photoisomerization mechanism and the S1-S0 crossing seam. We also provide a map of the nonadiabatic coupling between the two states in the region where they get closer in energy. Eventually, we show that treating the two CNN bending angles as independent coordinates is fundamental to break the symmetry and couple the two electronic states. The accuracy of the S0 and S1 PESs and couplings was validated with semiclassical dynamics simulations in the reduced space of the scanned coordinates, showing results in good agreement with published full-coordinate dynamics.
Aleotti F.; Soprani L.; Nenov A.; Berardi R.; Arcioni A.; Zannoni C.; Garavelli M.
File in questo prodotto:
File Dimensione Formato  
JCTC_AB_revision.pdf

embargo fino al 24/10/2020

Descrizione: Main text of the paper
Tipo: Postprint
Licenza: Licenza per Accesso Aperto. Creative Commons Attribuzione (CCBY)
Dimensione 8.04 MB
Formato Adobe PDF
8.04 MB Adobe PDF Visualizza/Apri
JCTC_SI.pdf

embargo fino al 24/10/2020

Descrizione: Supplementary file of the final revised version
Tipo: File Supplementare
Licenza: Licenza per Accesso Aperto. Creative Commons Attribuzione (CCBY)
Dimensione 19.71 MB
Formato Adobe PDF
19.71 MB Adobe PDF Visualizza/Apri
JCTC_AB_revision_cop.pdf

embargo fino al 24/10/2020

Descrizione: postprint
Tipo: Postprint
Licenza: Licenza per Accesso Aperto. Creative Commons Attribuzione (CCBY)
Dimensione 8.18 MB
Formato Adobe PDF
8.18 MB Adobe PDF Visualizza/Apri

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: http://hdl.handle.net/11585/731503
Citazioni
  • ???jsp.display-item.citation.pmc??? 7
  • Scopus 26
  • ???jsp.display-item.citation.isi??? 25
social impact