Estoy tratando de escribir un código para calcular el potencial y las fuerzas, para lo mismo usando la suma de Ewald. Para este propósito, la fórmula para el potencial y la fuerza que he usado es:
donde la contribución del potencial en el espacio k viene dada por
y las ecuaciones de fuerza -
Actualmente estoy usando este método para resolver la simulación de dinámica molecular de moléculas de agua. El problema al que me enfrento es que el cálculo del potencial depende en gran medida de , que puede ser, pero hay un gran cambio en el potencial si cambio levemente. En segundo lugar, ¿no debería el cambio con respecto disminuir más allá de cierto valor? No veo que eso también suceda.
Aquí deseo saber si hay una forma precisa de averiguar si mi código basado en las fórmulas anteriores funciona correctamente.
EDICIÓN 1: Como nota, he simulado el modelo de agua SPC/E usando esto ahora. El problema al que me enfrento en este momento es que la energía potencial que se calcula es de 5 a 6 veces mayor que la que se informó anteriormente en las revistas, pero sin embargo, la conservación de energía durante la ejecución está intacta. No entiendo la razón aquí !!
Solo he usado la suma de Ewald, nunca la he implementado yo mismo.
Sin embargo, mencionas que no estás convergiendo como aumenta ni está convergiendo al valor correcto. Parecería que, independientemente del problema, si su implementación es correcta, debería converger en algún punto.
Si alcanza la convergencia wrt ; en cuanto al punto de que no está obteniendo los valores informados. Recuerda que el tamaño del cuadro de simulación también es importante. ¿Está seguro de que está utilizando un espacio de simulación lo suficientemente grande con suficientes moléculas de agua?
EDITAR:
Basado en sus comentarios a continuación:
Bueno, la energía debe conservarse, :P. Si su campo de fuerza, sea lo que sea que implementó, es conservador, entonces la energía debe conservarse (o su integrador de movimiento está mal).
¿Qué conjunto estás usando? ¿El papel? Si es NVE, ¿está seguro de que la V es exactamente correcta? recuerde que incluso un cambio del 0,1% en el tamaño de la celda de la caja provocará un gran cambio en la energía de un material condensado. (¡intenta comprimir agua al 0,1 %!)
¿Estás simulando agua fluida o agua cristalina? De cualquier manera, ¿estás seguro de que no estás atrapado en un mínimo local? Intenta pulsar el volumen de la caja del sim para salir.
usuario35952
usuario35952