Mejora en el modelado de las energías vibratorias de moléculas diatómicas mediante la derivada fraccionaria generalizada

Por qué importan las diminutas vibraciones de las moléculas

Cada respiración que das, cada llama que arde y cada nube en el espacio están llenas de moléculas sencillas de dos átomos que vibran constantemente. Estas vibraciones microscópicas determinan cómo las moléculas absorben la luz, almacenan energía y participan en reacciones químicas. Para entender y predecir ese comportamiento, los científicos construyen modelos matemáticos del movimiento molecular. Este artículo presenta una nueva forma de modelar esas vibraciones con mayor precisión, usando un giro moderno del cálculo aplicado a un modelo molecular antiguo pero poderoso.

Cómo suelen imaginar los científicos una molécula que tiembla

Una molécula diatómica —formada por solo dos átomos— puede imaginarse como dos bolas unidas por un muelle. Cuando los átomos se acercan o se separan, suben o bajan por un paisaje energético. Durante casi un siglo, una fórmula sencilla conocida como el potencial de Morse ha sido un caballo de batalla para describir ese paisaje. Captura cómo se estira el enlace, cómo se endurece y cómo finalmente se rompe. Combinado con la ecuación de Schrödinger de la mecánica cuántica, el potencial de Morse permite a los investigadores calcular las energías vibratorias permitidas, que aparecen como líneas nítidas en los espectros moleculares medidos en el laboratorio o en estrellas distantes.

Dónde flaquea la imagen clásica

El modelo estándar de Morse funciona mejor para moléculas que no están girando y que no están muy excitadas. Sin embargo, las moléculas reales rotan mientras vibran, y su comportamiento se vuelve más complejo a niveles de energía altos. Para manejar la rotación se añade una barrera energética extra—llamada término centrífugo—pero esto hace que las ecuaciones sean difíciles de resolver exactamente. A lo largo de los años se han desarrollado muchas aproximaciones ingeniosas y trucos numéricos, pero siguen existiendo pequeñas discrepancias entre la teoría y datos experimentales precisos, especialmente para niveles vibratorios altos y para moléculas con distintos estados electrónicos.

Un nuevo giro fraccionario al movimiento molecular

Los autores adoptan una idea reciente del cálculo "fraccionario", que generaliza las derivadas ordinarias de modo que su orden puede ser no entero. En este marco, la respuesta de un sistema puede llevar memoria de su pasado, o codificar efectivamente interacciones sutiles y no locales. Partiendo de una derivada fraccionaria generalizada que preserva reglas clave del cálculo estándar, el equipo extiende una técnica de solución conocida (el método de Nikiforov–Uvarov) a lo que denominan método NU fraccionario generalizado. Usan este enfoque para obtener fórmulas analíticas de los niveles de energía vibratoria de moléculas en cualquier número de dimensiones espaciales, partiendo del potencial de Morse y tratando la barrera rotacional con una aproximación refinada.

Poniendo a prueba el nuevo modelo

Para comprobar si el enfoque fraccionario realmente ayuda, los autores lo aplican a veintidós moléculas diatómicas diferentes, incluidas especies importantes en química, ciencia de materiales y astrofísica como CO, Na₂, AlH, SiO⁺, TaO y TaS. Para cada caso construyen curvas de energía potencial a partir de constantes moleculares conocidas y comparan los niveles vibracionales predichos con datos de referencia de alta calidad obtenidos mediante análisis Rydberg–Klein–Rees (RKR). Tratan el orden fraccionario como un pequeño parámetro ajustable que resume efectivamente interacciones complejas no captadas por el modelo más simple. Minimizando sistemáticamente el error porcentual medio para cada molécula, muestran que las fórmulas fraccionarias reproducen las energías vibratorias observadas con una precisión notable, a menudo mejorando el acuerdo por un factor de dos o más respecto al límite clásico y compitiendo con éxito con potenciales alternativos de última generación y métodos numéricos.

Por qué los errores más pequeños cuentan una historia mayor

Más allá de las medidas de error promedio, el estudio examina cómo cambia la discrepancia entre teoría y experimento nivel por nivel. Para muchas moléculas, el modelo tradicional se aleja más de los datos a medida que crece el número cuántico vibracional. En contraste, el modelo fraccionario mantiene esas desviaciones más pequeñas y uniformes, especialmente en la parte alta del espectro donde las predicciones precisas son más exigentes. Los autores también prueban su solución analítica para estados rotantes de la molécula de CO y encuentran un excelente acuerdo con puntos de referencia numéricos muy precisos, lo que muestra que sus aproximaciones siguen siendo fiables incluso para movimientos rotacionales considerables.

Qué significa esto para entender las moléculas

En términos cotidianos, el trabajo muestra que "flexionar" ligeramente las reglas habituales del cálculo —permitiendo derivadas de orden no entero— conduce a una descripción más flexible y fiel de cómo vibran y rotan las moléculas sencillas. Afinando un único parámetro fraccionario, el modelo captura efectos sutiles que antes requerían potenciales más complicados o cálculos numéricos pesados. Esto lo convierte en una herramienta potente y eficiente para interpretar espectros moleculares en el laboratorio, diseñar nuevos materiales y sondear la composición de entornos astronómicos donde abundan las moléculas diatómicas.

Cita: Khokha, E.M., Abu-Shady, M., Omugbe, E. et al. Improved modelling for vibrational energies of diatomic molecules using the generalized fractional derivative. Sci Rep 16, 12037 (2026). https://doi.org/10.1038/s41598-026-39091-5

Palabras clave: moléculas diatómicas, espectros vibracionales, potencial de Morse, cálculo fraccionario, espectroscopía molecular