Density Functional Theory Investigation of Electronic, Spectroscopic, Reactivity, and Nonlinear Optical Properties of Uracil in Neutral, Cation, and Anion States

Resumen Técnico: Investigación mediante la Teoría del Funcional de la Densidad de las Propiedades Electrónicas, Espectroscópicas, de Reactividad y Ópticas No Lineales del Uracilo en los Estados Neutro, Catiónico y Aniónico

Planteamiento del Problema
El uracilo es una nucleobase fundamental en el ARN, desempeñando un papel crítico en la regulación y expresión genética. Aunque su relevancia biológica ha impulsado extensas investigaciones experimentales y teóricas, los estudios computacionales previos se han centrado predominantemente en la forma neutra del uracilo o en derivados específicos. Persiste una brecha significativa en la literatura con respecto a una comparación sistemática de múltiples propiedades del uracilo a través de sus diferentes estados de ionización (neutro, catiónico y aniónico). Específicamente, un análisis exhaustivo que integre optimización de geometría, orbitales moleculares frontera, propiedades termodinámicas, espectros UV-Vis, descriptores de reactividad (Fukui/dual), distribución de huecos-electrones, características topológicas (ELF, LOL, STM) y propiedades ópticas no lineales (NLO) para los tres estados de carga sigue estando poco explorado. Este estudio aborda esta brece para comprender mejor los efectos de la ionización relevantes para la química de la radiación, la electrónica biomolecular y los materiales NLO.

Metodología
La investigación emplea la Teoría del Funcional de la Densidad (DFT) utilizando el funcional de intercambio-correlación B3LYP y el conjunto de bases 6-311++G(d,p). Todos los cálculos se realizaron utilizando el software Gaussian 16. El estudio abarca los siguientes marcos analíticos:

• Optimización de Geometría: Realizada para los estados neutro, catiónico y aniónico en fase gaseosa.
• Termodinámica: La entropía (S), la capacidad calorífica (Cp) y la entalpía (H) se analizaron como funciones de la temperatura utilizando scripts de thermo.pl.
• Descriptores de Reactividad: Los índices de reactividad global (brechas HOMO-LUMO, dureza química, suavidad, electrofilicidad) y los descriptores locales (funciones de Fukui condensadas y descriptores duales) se calcularon utilizando el software UCI FUKUI.
• Espectroscopía: Se utilizó la DFT dependiente del tiempo (TD-DFT) con el modelo de solvatación C-PCM (agua) para simular los espectros de absorción UV-Vis.
• Análisis de Transferencia de Carga: El análisis de la distribución de huecos-electrones se llevó a cabo utilizando Multiwfn 3.8 para evaluar el carácter de transferencia de carga (CT) frente al de excitación local (LE).
• Análisis Topológico: Se computaron simulaciones de Microscopía de Túnel de Barrido (STM), Densidad Local de Estados (LDOS), Función de Localización de Electrones (ELF) y el Localizador de Orbitales Localizados (LOL) para visualizar la localización electrónica y el comportamiento de tunelización.
• Propiedades NLO: Las polarizabilidades estáticas y las hiperpolarizabilidades se evaluaron para valorar el potencial óptico no lineal, comparándolas con la urea.
• Estructura Electrónica: El análisis de Orbitales de Enlace Natural (NBO) y el análisis de población de Mulliken se utilizaron para estudiar las interacciones intramoleculares y la distribución de carga.

Resultados Clave

• Cambios Estructurales: La ionización altera significativamente las longitudes y los ángulos de enlace. Se observaron cambios notables en los enlaces N4-C5 y N4-C7. El ángulo de enlace C6-C8-C7 varió de 119.90° (neutro) a 120.27° (catión) y 122.31° (anión), confirmando una geometría plana con hibridación sp².
• Estabilidad Electrónica y Reactividad: Las brechas de energía HOMO-LUMO fueron de 5.622 eV (neutro) y 5.618 eV (catión α), indicando una alta estabilidad cinética. En contraste, los estados aniónico α (1.422 eV) y catiónico β (2.148 eV) exhibieron brechas significativamente menores, lo que sugiere una mayor reactividad química. La suavidad global fue mayor para el anión (1.406 eV), lo que indica una mayor polarizabilidad.
• Termodinámica: La entropía, la capacidad calorífica y la entalpía aumentaron con la temperatura para todos los estados. El estado aniónico mostró consistentemente los valores más altos para estas propiedades, lo que indica un sistema más desordenado con mayor contenido de calor en comparación con los estados neutro y catiónico.
• Espectroscopía y Transferencia de Carga: La TD-DFT reveló una absorción principal a 242 nm (f = 0.1954) en el segundo estado excitado, impulsada principalmente por una transición HOMO→LUMO (68.33%). El análisis de huecos-electrones indicó que la mayoría de las excitaciones poseen un carácter de transferencia de carga (CT), excepto la segunda excitación en el estado neutro, que mostró un carácter de excitación local (LE). El estado aniónico exhibió la mayor longitud de transferencia de carga (índice D) y la menor energía de atracción de Coulomb entre el hueco y el electrón.
• Sitios Reactivos: Las funciones de Fukui condensadas y los descriptores duales identificaron al C8 como el sitio primario para el ataque electrofílico y al C7 para el ataque nucleofílico. El orden del descriptor dual para valores positivos fue C7 > C6 > N3 > H11, mientras que los valores negativos siguieron el orden C8 > N4 > O1 > O2 > C5.
• Características Topológicas: Las simulaciones de STM mostraron corrientes de tunelización fuertes en el estado catiónico en los átomos C7, C8 y N3. Los mapas de ELF y LOL confirmaron nubes de deslocalización alrededor de los átomos de oxígeno, nitrógeno y carbono, identificándose regiones covalentes entre los átomos C-C y C-N.
• Propiedades Ópticas No Lineales (NLO): El estado aniónico demostró una respuesta NLO notablemente mejorada. La hiperpolarizabilidad estática (β) del anión de uracilo se calculó aproximadamente 172 veces mayor que la de la urea (el umbral estándar). Los estados catiónico y neutro también mostraron propiedades Nológicas mejoradas en comparación con la urea (8.66 y 1.55 veces, respectivamente).
• Interacciones Intramoleculares: El análisis NBO destacó energías de estabilización fuertes derivadas de las interacciones LP(N) → π*(C=O) y LP(O) → σ*. La mayor energía de estabilización (58.96 kcal/mol) se observó para la interacción LP(1)(N3) → π*(O1-C5).

Significancia y Reivindicaciones
Los autores afirman que este estudio sistemático de múltiples estados proporciona una comprensión más profunda de cómo la ionización afecta las propiedades estructurales, electrónicas y ópticas del uracilo. El trabajo establece que el estado aniónico del uracilo posee brechas de energía significativamente reducidas y una transferencia de carga mejorada, lo que conduce a un aumento dramático en la respuesta óptica no lineal. El estudio sugiere que estos hallazgos son relevantes para comprender el daño inducido por la radiación en el ARN, los procesos de transferencia de carga en sistemas biomoleculares y el potencial de aplicación de nucleobases ionizadas en la electrónica molecular y los materiales NLO. El artículo se posiciona como una solución que llena un vacío de investigación específico al ofrecer una comparación exhaustiva y comparativa de los tres estados de carga, un alcance que anteriormente carecía de la literatura.