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

Resumen Técnico: Un Modelo Geométrico Unificado de la Tensión de Anillo: Del Carbono al Nitrógeno y el Papel de la Repulsión Coulombiana de los Pares Solitarios

Planteamiento del Problema
La energía de tensión de anillo es un determinante fundamental de la estabilidad y reactividad de los compuestos cíclicos. A pesar de su centralidad en la química orgánica, una teoría cuantitativa unificada capaz de explicar la tensión en diferentes elementos del grupo principal ha permanecido esquiva. Los enfoques previos presentan limitaciones: los campos de fuerza de la mecánica molecular carecen de una conexión directa con observables de la mecánica cuántica; los modelos electrónicos tempranos (p. ej., "enlaces banana") proporcionaron una visión profunda pero carecieron de descriptores predictivos generalizables; y los modelos modernos de aprendizaje automático suelen funcionar como "cajas negras" sin interpretabilidad física. Este estudio aborda estas brechas buscando un conjunto mínimo de descriptores físicamente interpretables para predecir cuantitativamente la tensión de anillo para elementos estables y mapeando el panorama de estabilidad de los anillos homoatómicos del segundo periodo.

Metodología
El estudio emplea un flujo de trabajo computacional riguroso y multinivel que combina la teoría de la estructura electrónica de alto nivel con regresión estadística avanzada:

• Tamizaje de Estabilidad: Las geometrías iniciales para anillos de Be, B, C, N y O ($n=3–10$) se optimizaron al nivel B3LYP-D3(BJ)/def2-SVP. La estabilidad se definió por la convergencia a un mínimo local cíclico sin disociación.
• Muestreo Conformacional: Para anillos de Carbono (C3–C14) y Nitrógeno (N3–N14), se generaron conjuntos conformacionales utilizando CREST/GFN2-xTB. Los confórmeros se refinaron mediante un embudo de DFT multinivel (SVP $\to$ TZVPP $\to$ QZVPP) para obtener energías libres de Gibbs promediadas del conjunto.
• Cálculo de Energía de Alta Precisión: Las energías libres de Gibbs finales se derivaron mediante un esquema compuesto: energías de punto único DLPNO-CCSD(T)/def2-TZVPP corregidas por contribuciones térmicas/entrópicas de B3LYP-D3(BJ)/def2-QZVPP. Las energías de tensión de anillo ($E_{strain}$) se calcularon mediante reacciones de transferencia de hidrógeno isodesmicas.
• Extracción de Descriptores: Los descriptores geométricos ($\Delta\theta^2$, $\Delta R^2$) se derivaron de las geometrías optimizadas. Los descriptores electrónicos (cargas ADCH, momentos multipolares) se computaron con Multiwfn. Se formuló un término específico para la repulsión Coulombiana de los pares solitarios ($E_{coul}$) para los anillos de nitrógeno.
• Análisis Estadístico: Se realizó una regresión de Mínimos Cuadrados Ordinarios (OLS) de subconjunto completo contra las energías de referencia CCSD(T). La selección del modelo utilizó $R^2$, RMSE de Validación Cruzada Leave-One-Out (LOOCV), Factor de Inflación de la Varianza (VIF) y criterios AIC/BIC.
• Caracterización Electrónica: El comportamiento anómalo del N5 se investigó mediante cálculos CASSCF(6,5) y computaciones de Desplazamiento Isotópico Inducido por Núcleo (NICS) para sondear la $\sigma$-homoaromaticidad.

Contribuciones Clave y Resultados

1. Dicotomía de Estabilidad de los Anillos del Segundo Periodo:
   El tamizaje computacional revela un límite de estabilidad nítido. Solo el Carbono (C) y el Nitrógeno (N) soportan anillos homoatómicos estables de enlace simple a través de un amplio rango de tamaños (C3–C30+ y N3–N14+).

   • Berilio (Be): Los anillos colapsan en jaulas deficientes de electrones debido a la deficiencia electrónica.
   • Boro (B): Los anillos mayores a B4 son inestables debido a la desestabilización de enlaces multicéntricos.
   • Oxígeno (O): Los anillos más allá de O4 se descomponen mediante la formación de cadenas peróxido impulsada por una abrumadora repulsión de pares solitarios.
     Esta estabilidad se racionaliza mediante un principio de coincidencia de simetría: C y N poseen configuraciones electrónicas (tetraédrica/piramidal $sp^3$) compatibles con topologías cíclicas, mientras que Be, B y O favorecen estructuras alternativas.

2. Modelo Geométrico Unificado para Anillos de Carbono:
   Para los anillos de carbono (C3–C14), un modelo de dos parámetros puramente geométrico basado en la desviación del ángulo de enlace ($\Delta\theta^2$) y la distorsión de la longitud de enlace ($\Delta R^2$) captura el 98.2% de la varianza de la tensión ($R^2 = 0.982$, RMSE = 0.36 kcal mol$^{-1}$ átomo$^{-1}$). La inclusión de descriptores electrónicos produce una mejora insignificante, confirmando que la tensión de los anillos de carbono está dominada por factores geométricos.

3. Papel de la Repulsión de Pares Solitarios en Anillos de Nitrógeno:
   Para los anillos de nitrógeno, el esqueleto geométrico por sí solo explica el 92.5% de la varianza. Sin embargo, la incorporación de un término de repulsión Coulombiana de pares solitarios físicamente derivado ($E_{coul}$) mejora significativamente el modelo a $R^2 = 0.968$ (RMSE = 0.64 kcal mol$^{-1}$ átomo$^{-1}$). Esto confirma que la tensión residual persistente en los anillos de nitrógeno (que permanece elevada en comparación con el carbono) es impulsada por las interacciones repulsivas entre los pares solitarios del nitrógeno.

4. La Anomalía del N5 y la $\sigma$-Homoaromaticidad:
   El N5 es identificado como un valor atípico estadístico con un residuo de $-4.0$ kcal mol$^{-1}$, lo que indica una estabilización inesperada. Esto se atribuye a la $\sigma$-homoaromaticidad: los cinco pares solitarios del nitrógeno en el plano se deslocalizan a través del marco $\sigma$, creando una corriente de anillo diatrópica (NICS(0) = $-15.18$ ppm). Este efecto contrarresta la repulsión de los pares solitarios, un fenómeno distinto de la $\pi$-aromaticidad del benceno, como lo evidencia la rápida atenuación de la corriente de anillo por encima del plano molecular (NICS(1)).

Significado y Marco Propuesto
El artículo propone una teoría estratificada de la tensión de anillo que unifica efectos geométricos, electrostáticos y cuánticos de aromaticidad en un único marco cuantitativo:
$$E_{strain} = E_{geom}(\Delta\theta^2, \Delta R^2) + E_{LP}(E_{coul}) + E_{\sigma-aromatic}(NICS)$$

• $E_{geom}$: Un esqueleto geométrico universal común a todos los anillos homoatómicos de los grupos principales, recuperando la tensión angular de Baeyer y añadiendo la distorsión de la longitud de enlace.
• $E_{LP}$: Una corrección Coulombiana de pares solitarios requerida para elementos con pares solitarios (pnictógenos), que se convierte en el factor no geométrico dominante para el nitrógeno.
• $E_{\sigma-aromatic}$: Una corrección cuántica no aditiva para casos específicos (como el N5) donde las condiciones geométricas y electrónicas permiten la deslocalización de los pares solitarios a través de los enlaces.

Este trabajo establece que, si bien la tensión de anillo es fundamentalmente geométrica para el carbono, requiere un tratamiento explícito de la repulsión de los pares solitarios para el nitrógeno, siendo la $\sigma$-aromaticidad una corrección crítica para anillos pequeños. El estudio cierra la brecha entre los campos de fuerza empíricos y la realidad de la mecánica cuántica al proporcionar descriptores directamente vinculados a las energías CCSD(T) y a la estructura electrónica, ofreciendo un modelo transparente y físicamente fundamentado para predecir la estabilidad cíclica.