¿Qué significa f=−∇uf=−∇u\textbf{f} = -\boldsymbol{\nabla} u en la práctica y cómo se calcula?

En las simulaciones por computadora clásicas, como las simulaciones de dinámica molecular (MD), se integran las ecuaciones de movimiento de Newton para determinar las trayectorias de las partículas. Si pensamos en la Segunda Ley de Newton como simplemente

F i = d pag i d t
dónde F i es la fuerza neta que actúa sobre una partícula i y pag i es el momento de la partícula, entonces está claro que para integrarlo en última instancia y propagar la trayectoria, necesitaremos calcular la fuerza en cada partícula en el sistema.

Sin embargo, en la práctica, creo que la energía potencial tu , en lugar de la fuerza, se calcula inicialmente (un campo de fuerza especifica la forma de la energía potencial V ). Ahora, sé que existe una relación entre la energía potencial tu y fuerza F :

F = tu

En otras palabras, la fuerza es el gradiente negativo de la energía potencial.

Estoy leyendo Comprensión de la simulación molecular de Frenkel y Smit (segunda edición) , y en la página 69 (Google Books tiene algunas páginas aquí ), veo este párrafo:

Supongamos que deseamos calcular el X -componente de la fuerza:

F X ( r ) = tu ( r ) X = ( X r ) ( tu ( r ) r )

Me temo un poco que esta es una pregunta tonta, pero ¿cómo pasan los autores de tu ( r ) X a ( X r ) ( tu ( r ) r ) ? ¿Es esta la regla de la cadena? ¿Podría ayudarme a ver el paso que están dando los autores?

Respuestas (2)

Tienes razón en que el resultado que ves se debe a la regla de la cadena . El autor usa coordenadas esféricas o cilíndricas , por lo que

r = X 2 + y 2 + z 2

o

r = X 2 + y 2

que puedes diferenciar para obtener

r X = X r

Por eso

F X ( r ) = tu ( r ) X = tu ( r ) r r X = X r tu ( r ) r

La respuesta anterior está bien, pero no debes perder el tiempo haciendo la regla de la cadena. Proviene de una relación diferencial que se aprende una vez y se convierte para siempre en una segunda naturaleza:

r d r = X d X + y d y + z d z

Esta es la forma diferencial de la ley de Pitágoras:

r 2 = X 2 + y 2 + z 2

y te dice el tamaño del tamaño relativo de los incrementos x,y,z, cuando haces un incremento de r. El incremento x es X r d r , por lo que el cambio en la función es

tu ( X + d X , y , z ) = tu ( r + d r ) = tu ( r ) + tu ( r ) d r = tu ( r ) + tu ( r ) X r d X

es importante hacerlo así, para que entre en tu cabeza y se quede ahí.