Obrigado pela pergunta Carlos. Nós avaliamos o gap singleto-tripleto apenas em alguns casos de PAHs. Para o benzeno[2+] é mais grave, o estado singleto na estrutura "piramidal pentagonal" é o mínimo global mas é próximo em energia do tripleto na estrutura "benzenóide". Estamos mais curiosos em relação aos PAHs maiores no estado dicátion e tetracátion. Como existe uma tendência de abrir a camada quando o PAH crescer "linearmente", isso pode influenciar nos outros estados de carga. 

Olá Aline, muito obrigado pelo interesse! Devido ao pouco tempo eu optei por não falar dos métodos. A busca de mínimo global costuma ser feita em nível DFT utilizando o funcional PBE0 e uma base double zeta como a 6-31G(d). Depois reotimizamos com uma base maior (triple zeta). Em seguida fazemos cálculos single point em nível CCSD(T) ou CCSD-F12 para refinar. Em alguns casos quando suspeitamos de caráter multireferência, fazemos CASSCF/NEVPT2, se possível.

 No momento queremos estudar os caminhos de isomerização em função do estado de carga e do tipo de condensação dos PAHs.

Olá Eduardo, obrigado! 1) O efeito do N é similar ou fósforo porém menos pronunciado. O Cl levou a uma situação curiosa, foi caracterizada uma ligação C=Cl [10.1021/jacs.7b12749]. 2) Sim, queremos estudar diferentes estequiometrias de C/N/H e também queremos estudar comparativamente O e S.

Henrique Stark! hahaha