A Unified Geometric Model of Ring Strain: From Carbon to Nitrogen and the Role of Lone-Pair Coulomb Repulsion

Resumo Técnico: Um Modelo Geométrico Unificado da Tensão de Anel: Do Carbono ao Nitrogênio e o Papel da Repulsão de Coulomb de Pares Isolados

Definição do Problema
A energia de tensão de anel é um determinante fundamental da estabilidade e reatividade de compostos cíclicos. Apesar de sua centralidade na química orgânica, uma teoria quantitativa unificada capaz de explicar a tensão em diferentes elementos do grupo principal permaneceu elusiva. Abordagens anteriores apresentam limitações: campos de força de mecânica molecular carecem de conexão direta com observáveis da mecânica quântica; modelos eletrônicos iniciais (ex: "ligações banana") forneceram insights, mas careceram de descritores preditivos generalizáveis; e modelos modernos de aprendizado de máquina frequentemente funcionam como "caixas pretas" sem interpretabilidade física. Este estudo aborda essas lacunas ao buscar um conjunto mínimo de descritores fisicamente interpretáveis para prever quantitativamente a tensão de anel para elementos estáveis e ao mapear o panorama de estabilidade de anéis homoatômicos do segundo período.

Metodologia
O estudo emprega um fluxo de trabalho computacional rigoroso de múltiplos níveis, combinando teoria de estrutura eletrônica de alto nível com regressão estatística avançada:

• Triagem de Estabilidade: Geometrias iniciais para anéis de Be, B, C, N e O ($n=3–10$) foram otimizadas no nível B3LYP-D3(BJ)/def2-SVP. A estabilidade foi definida pela convergência para um mínimo local cíclico sem dissociação.
• Amostragem Conformacional: Ensembles conformacionais para anéis de Carbono (C3–C14) e Nitrogênio (N3–N14) foram gerados usando CREST/GFN2-xTB. Os conformeros foram refinados através de um funil de DFT de múltiplos níveis (SVP $\to$ TZVPP $\to$ QZVPP) para obter energias livres de Gibbs médias do ensemble.
• Cálculo de Alta Precisão de Energia: As energias livres de Gibbs finais foram derivadas de um esquema composto: energias de ponto único DLPNO-CCSD(T)/def2-TZVPP corrigidas por contribuições térmicas/entrópicas de B3LYP-D3(BJ)/def2-QZVPP. As energias de tensão de anel ($E_{strain}$) foram calculadas via reações isodesmicas de transferência de hidrogênio.
• Extração de Descritores: Descritores geométricos ($\Delta\theta^2$, $\Delta R^2$) foram derivados das geometrias otimizadas. Descritores eletrônicos (cargas ADCH, momentos multipolares) foram computados usando Multiwfn. Um termo específico para a repulsão de Coulomb de pares isolados ($E_{coul}$) foi formulado para anéis de nitrogênio.
• Análise Estatística: Regressão de Mínimos Quadrados Ordinários (OLS) de subconjunto total foi realizada contra energias de referência CCSD(T). A seleção de modelos utilizou $R^2$, RMSE de Validação Cruzada de Um Único Elemento Deixado Fora (LOOCV), Fator de Inflação de Variância (VIF) e critérios AIC/BIC.
• Caracterização Eletrônica: O comportamento anômalo do N5 foi investigado através de cálculos CASSCF(6,5) e computações de Deslocamento Químico Independente de Núcleo (NICS) para sondar a $\sigma$-homoaromaticidade.

Principais Contribuições e Resultados

1. Dicotomia de Estabilidade de Anéis do Segundo Período:
   A triagem computacional revela uma fronteira de estabilidade nítida. Apenas o Carbono (C) e o Nitrogênio (N) suportam anéis homoatômicos estáveis com ligações simples em uma ampla faixa de tamanho (C3–C30+ e N3–N14+).

   • Berílio (Be): Os anéis colapsam em gaiolas deficientes de elétrons devido à deficiência eletrônica.
   • Boro (B): Anéis maiores que B4 são instáveis devido à desestabilização de ligações multicêntricas.
   • Oxigênio (O): Anéis além de O4 se decompõem via formação de cadeias peróxidas, impulsionados pela esmagadora repulsão de pares isolados.
     Esta estabilidade é racionalizada por um princípio de correspondência de simetria: C e N possuem configurações eletrônicas (tetraédrica/piramidal $sp^3$) compatíveis com topologias cíclicas, enquanto Be, B e O favorecem estruturas alternativas.

2. Modelo Geométrico Unificado para Anéis de Carbono:
   Para anéis de carbono (C3–C14), um modelo puramente geométrico de dois parâmetros baseado no desvio do ângulo de ligação ($\Delta\theta^2$) e na distorção do comprimento de ligação ($\Delta R^2$) captura 98,2% da variância da tensão ($R^2 = 0,982$, RMSE = 0,36 kcal mol$^{-1}$ átomo$^{-1}$). A inclusão de descritores eletrônicos produz melhoria negligenciável, confirmando que a tensão do anel de carbono é dominada por fatores geométricos.

3. Papel da Repulsão de Pares Isolados em Anéis de Nitrogênio:
   Para anéis de nitrogênio, a espinha dorsal geométrica sozinha explica 92,5% da variância. No entanto, a incorporação de um termo de repulsão de Coulomb de pares isolados derivado fisicamente ($E_{coul}$) melhora significamente o modelo para $R^2 = 0,968$ (RMSE = 0,64 kcal mol$^{-1}$ átomo$^{-1}$). Isso confirma que a tensão residual persistente nos anéis de nitrogênio (que permanece elevada em comparação ao carbono) é impulsionada pelas interações repulsivas entre os pares isolados do nitrogênio.

4. A Anomalia do N5 e $\sigma$-Homoaromaticidade:
   O N5 é identificado como um outlier estatístico com um resíduo de $-4,0$ kcal mol$^{-1}$, indicando uma estabilização inesperada. Isso é atribuído à $\sigma$-homoaromaticidade: os cinco pares isolados do nitrogênio no plano se deslocalizam através da estrutura $\sigma$, criando uma corrente de anel diatrópica (NICS(0) = $-15,18$ ppm). Este efeito neutraliza a repulsão de pares isolados, um fenômeno distinto da $\pi$-aromaticidade do benzeno, como evidenciado pela rápida atenuação da corrente de anel acima do plano molecular (NICS(1)).

Significância e Estrutura Proposta
O artigo propõe uma teoria em camadas da tensão de anel que unifica efeitos geométricos, eletrostáticos e quânticos aromáticos em um único framework quantitativo:
$$E_{strain} = E_{geom}(\Delta\theta^2, \Delta R^2) + E_{LP}(E_{coul}) + E_{\sigma-aromatic}(NICS)$$

• $E_{geom}$: Uma espinha dorsal geométrica universal comum a todos os anéis homoatômicos de grupos principais, recuperando a tensão de ângulo de Baeyer e adicionando a distorção do comprimento de ligação.
• $E_{LP}$: Uma correção de Coulomb de pares isolados necessária para elementos com pares isolados (pnictogênios), que se torna o fator não geométrico dominante para o nitrogênio.
• $E_{\sigma-aromatic}$: Uma correção quântica não aditiva para casos específicos (como o N5) onde as condições geométricas e eletrônicas permitem a deslocalização de pares isolados através da ligação.

Este trabalho estabelece que, embora a tensão de anel seja fundamentalmente geométrica para o carbono, ela requer o tratamento explícito da repulsão de pares isolados para o nitrogênio, com a $\sigma$-aromaticidade servindo como uma correção crítica para anéis pequenos. O estudo preenche a lacuna entre campos de força empíricos e a realidade da mecânica quântica, fornecendo descritores diretamente ligados às energias CCSD(T) e à estrutura eletrônica, oferecendo um modelo transparente e fisicamente fundamentado para prever a estabilidade cíclica.