The Diels–Alder reactivity of C59NH azafullerene has been explored computationally. The regioselectivity of the process and the factors controlling the reduced reactivity of this system with respect to the parent C60 fullerene have been analyzed in detail by using the activation strain model of reactivity and the energy decomposition analysis method. It is found that the presence of the nitrogen atom and the CH fragment in the fullerene reduces the interaction between the deformed reactants along the entire reaction coordinate

We are grateful for financial support from the Spanish MINECO-FEDER (Grants CTQ2013-44303-P, CTQ2016-78205-P, and CTQ2014-51912-REDC to I.F. and CTQ2014-54306-P to M.S.), Fundacion BBVA (Convocatoria 2015 de ́ Ayudas Fundación BBVA a Investigadores y Creadores Culturales), and Catalan DIUE (Projects 2014SGR931, ICREA Academia 2014 prize, and XRQTC to M.S.)