⚗️ Superficie de Energía Potencial

Es una superficie de energía potencial analítica bidimensional muy utilizada, construida a partir de una suma de cuatro términos exponenciales de tipo gaussiano, cada uno con su propia amplitud, posición y curvatura. Fue diseñada para tener exactamente tres mínimos locales y dos puntos de silla de primer orden dispuestos de modo que un camino de un mínimo a otro deba pasar por ambos puntos de silla, convirtiéndola en un caso de prueba estándar para algoritmos que buscan caminos de reacción y estados de transición.

La partícula sigue una dinámica de tipo Langevin: en cada paso temporal su velocidad se actualiza usando el gradiente negativo del potencial (la fuerza que la atrae cuesta abajo), un término de fricción proporcional a su velocidad actual y un impulso térmico aleatorio escalado por el deslizador de temperatura. El gradiente se calcula con una aproximación de diferencias finitas centrada, y la partícula rebota con velocidad reducida si alcanza el borde del dominio graficado.

La fricción amortigua la velocidad de la partícula en cada paso, por lo que una fricción más alta hace que el movimiento sea más lento y sobreamortiguado, asentándose en un mínimo sin mucha oscilación. La temperatura añade una fuerza aleatoria escalada a su valor en cada paso, dando a la partícula suficiente energía cinética para ocasionalmente volver a subir sobre una barrera que ya ha cruzado — exactamente análogo a cómo la energía térmica en una molécula real le permite cruzar una barrera de activación aunque la reacción vaya cuesta arriba en energía potencial.

La energía de activación Ea es la diferencia de energía entre un estado de transición y el pozo de reactivo del que parte — en esta simulación, Ea = V(TS1) − V(R), leída en vivo desde el panel de lectura. Una energía de activación más alta significa que la barrera entre reactivo y producto es más alta, por lo que a una temperatura fija la partícula necesita una fluctuación térmica mayor y más rara para cruzarla, reflejando exactamente por qué las reacciones reales con altas energías de activación transcurren más lentamente a una temperatura dada.

Muchas reacciones químicas reales no van directamente de reactivos a productos en un solo paso; en cambio, pasan por uno o más intermedios — especies estables pero de mayor energía — cada uno separado de sus vecinos por su propio estado de transición. Las reacciones catalizadas por enzimas, las sustituciones nucleofílicas SN1 y los mecanismos orgánicos de varios pasos son ejemplos clásicos. Trazar la superficie de energía potencial completa, como hace esta simulación en una forma bidimensional simplificada, es exactamente cómo los químicos computacionales identifican estos intermedios y el estado de transición limitante de la velocidad que controla la velocidad global de la reacción.