Método predictor-corrector siempre estable y reversible en el tiempo para la dinámica molecular de moléculas polarizables

Se propone un método mejorado para la dinámica molecular clásica de moléculas polarizables. El método utiliza un predictor, una evaluación del campo electrostático por paso de integración y relajación ( amortiguamiento). La solución autoconsistente se aproxima con error de segundo orden (respecto al paso de tiempo). El error de reversibilidad temporal (conservación de energía a largo plazo) es del orden (2n - 1), donde n es la longitud del predictor. El método es fácil de implementar, eficiente, preciso y adecuado para cualquier modelo de polarizabilidad.