Ola professor Luciano, obrigada!!!! Essa formação da ligação ainda é pra nós um caso a ser investigado. Por dinâmica molecular pensamos numa investigação da primeira camada de solvatação pela função de distribuição radial como forma de completar a segunda parte do trabalho que é a interação do CO2. Fazendo de maneira "manual", ao inserir o CO2, notamos que a partir de mais ou menos 18 moléculas é quase que constante a influência. Então essa opção surgiu numa tentativa de analisar a primeira camada de solvatação, mas ainda é so uma ideia. Entretanto estou considerando novas contribuições. 

Ola, 

Os cálculos foram todos realizados no software gaussian, sendo o DFT na versão 09 e o DFTBA no 16. A distância inicial de 5A foi escolhida com base na literatura em relação a alguns trabalhos sobre os planos aromáticos como as estruturas de grafeno. Essa distância inicial é entre o centro de massa dos monômeros que compõem o dímero, na posição Face to Face, isto é, em paralelo.  Portanto todas as configurações partem da mesma distância.

A energia de agregação na presença do CO2 é definida pela equação abaixo:

Eagg = EAB – (EA + EB)

Onde, EAB refere-se a energia do dímero com todas as moléculas de CO2 e, EA e EB aos respectivos monômeros com as moléculas de CO2 sobre a sua superfície. As energias foram retiradas do cálculo counterpoise.

Ola Rodolfo, obrigada pela contribuição.

Ainda não conseguimos racionalizar o real motivo da influencia da posição na mudança do gap entre o Ncentral e o Nperiférico. Somente avaliamos a relação entre a aromaticidade e o gap, como mostrado a posição Nperiferico por ter o maior gap é o mais aromático.

Quanto a segunda pergunta, neste momento só avaliamos de maneira geral a influência do CO2 nas interações entre os monômeros. O segundo passo dessa contribuição é avaliar a ação dele como um todo no dímero. 

O trabalho ainda esta em fase inicial, e temos como metas futuras essas avaliações, além das fragmentações dessas energias para cada aglomerado.

Oi Rodolfo, então para quantificar a aromaticidade neste trabalho eu usei somente a diferença entre os orbitais de fronteira. Entretanto para tese, estou realizando os cálculos HOMA e I6. Obrigada pelas contribuições.

Ola Ricardo, muito Obrigada!!! 1) Sim, na minha tese que está em andamento eu realizei cálculos na outra posição interna. O N esta substituindo o C acima do que neste trabalho é a substituição Ncentral. 2) Na formação do dímero, por aproximação e comparação com as demais, mantivemos o singleto ao invés do tripleto. Entretanto já esta como meta realizar este cálculo para avaliar essa diferença. Ainda não conseguimos chegar a um resultado específico sobre essa ligação entre os C na posição Ncentral. O que eu tenho até agora é que essa ligação não acontece na outra posição do centro a qual você perguntou. 3) Quanto a usar dinâmica, p este caso em particular nao tinha pensado e qualquer contribuição será bem vinda. Estamos pensando em dinâmica para completar a segunda parte no estudo da influência do CO2 mas ainda é algo recente. 

Obrigada Ricardo pelas contribuições!