Integración de Verlet
El algoritmo de Verlet y su forma de velocidad son los integradores estándar de la dinámica molecular, valorados porque son reversibles en el tiempo, simplécticos y conservan bien la energía a lo largo de los millones de pasos que requiere una simulación.
Definition
La integración de Verlet es un método simpléctico y reversible en el tiempo para integrar las ecuaciones de movimiento de Newton que actualiza las posiciones de las partículas utilizando posiciones actuales y previas y la aceleración, lo que produce trayectorias estables para la dinámica molecular.
Scope
Este tema cubre la familia de integradores de Verlet: el esquema original de Verlet de posición, las formulaciones equivalentes de leapfrog y Verlet de velocidad, su reversibilidad en el tiempo y estructura simpléctica, y la conservación de energía a largo plazo resultante. Sitúa estos métodos dentro de la teoría más amplia de la integración simpléctica de sistemas hamiltonianos.
Core questions
- ¿Cómo avanza el esquema de Verlet las posiciones y velocidades a partir de las fuerzas?
- ¿Por qué el algoritmo de Verlet es reversible en el tiempo y simpléctico?
- ¿Por qué la integración de Verlet conserva bien la energía en simulaciones muy largas?
- ¿Cómo se relacionan las formulaciones de Verlet de posición, leapfrog y velocidad entre sí?
Key theories
- Estructura simpléctica y reversible en el tiempo
- La integración de Verlet preserva la geometría simpléctica del espacio de fase y es invariante bajo la inversión del tiempo, lo que en conjunto previene la deriva energética sistemática que afecta a los integradores no simplécticos de sistemas conservativos.
- Conservación del hamiltoniano de sombra
- Aunque la trayectoria discreta de Verlet no conserva exactamente la energía verdadera, casi conserva un hamiltoniano de sombra estrechamente relacionado, manteniendo el error de energía acotado y oscilatorio en lugar de creciente.
- Formulaciones equivalentes
- Los esquemas de Verlet de posición, leapfrog y Verlet de velocidad producen la misma trayectoria, pero difieren en cómo y cuándo las velocidades están disponibles, siendo preferible el Verlet de velocidad cuando se necesitan posiciones y velocidades sincronizadas.
Clinical relevance
La integración de Verlet es el motor de avance temporal predeterminado en esencialmente todos los códigos de dinámica molecular, desde fluidos simples de Lennard-Jones hasta grandes simulaciones biomoleculares, y el mismo principio simpléctico se utiliza en la integración orbital a largo plazo en astronomía.
History
El esquema fue utilizado por el astrónomo Carl Stormer a principios del siglo XX y popularizado para la simulación molecular por Loup Verlet en su estudio de 1967 sobre fluidos de Lennard-Jones; análisis posteriores demostraron que era un integrador simpléctico, lo que explica su excelente estabilidad a largo plazo.
Key figures
- Loup Verlet
- Carl Stormer
- Ernst Hairer
Related topics
Seminal works
- verlet1967
- hairer1993
Frequently asked questions
- ¿Por qué se prefiere Verlet a un método de Runge-Kutta de orden superior en dinámica molecular?
- Aunque Runge-Kutta puede ser más preciso por paso, no es simpléctico y se desvía lentamente en energía en ejecuciones largas. La estructura simpléctica y reversible en el tiempo de Verlet mantiene la energía acotada durante millones de pasos, lo que importa mucho más que la precisión por paso para simulaciones de equilibrio.
- ¿Verlet conserva la energía exactamente?
- No. Conserva un hamiltoniano de sombra cercano en lugar de la energía exacta, por lo que la energía medida oscila dentro de una banda acotada en lugar de desviarse, lo cual es suficiente para calcular promedios termodinámicos estables.