¿Cuál es la forma correcta de calcular el momento principal de inercia?

Estoy escribiendo un programa que incorpora el cálculo del momento principal de inercia para un residuo de proteína basado en las coordenadas XYZ de su átomo componente. Estoy extremadamente confundido acerca de qué fórmulas usar para calcular el momento principal de inercia para mi situación.

Hasta ahora mi programa hace lo siguiente:

  1. Calcular elementos (Ixx, Iyy, Izz, etc.) de la matriz simétrica A
  2. Encuentre los valores propios (y por lo tanto los principales momentos de inercia)

Las fórmulas que usé originalmente para calcular los elementos en el paso 1 eran de la pág. 5 de este documento . Sin embargo, mirando Wikipedia , las fórmulas para los elementos de la matriz simétrica son diferentes.

También estuve buscando en Internet y encontré un código fuente abierto donde el centro de masa se resta de las coordenadas x, y y z antes de comenzar el cálculo de los elementos Ixx, Iyy, Izz, etc.

¿Qué fórmulas/algoritmo utilizo para calcular el momento principal de inercia para mi caso? No estoy seguro de por dónde empezar a elegir qué fórmulas usar. ¿Hay alguna fuente accesible para aquellos con poca experiencia en física que me ayude a comprender qué fórmula/algoritmo usar para calcular el momento principal de inercia?

Respuestas (1)

Las fórmulas del documento de la NASA y Wikipedia están simplemente en marcos diferentes: las de Wikipedia asumen un marco de coordenadas cuyos ejes pasan por el centro de masa, mientras que el documento de la NASA asume un marco de coordenadas cartesiano arbitrario.

Sean mayúsculas X i , Y i y Z i las coordenadas del i-ésimo átomo en un marco arbitrario de coordenadas cartesianas O y x i , y i y z i denoten las respectivas coordenadas en el marco O com trasladadas al centro de masa de la molécula .

Las coordenadas del centro de masa de la molécula en el marco O son:

X C o metro = 1 METRO i metro i X i Y C o metro = 1 METRO i metro i Y i Z C o metro = 1 METRO i metro i Z i

La transformación de O a O com es simplemente traducción por [-X com , -Y com , -Z com ]:

X i = X i X C o metro = X i 1 METRO i metro i X i y i = Y i Y C o metro = Y i 1 METRO i metro i Y i z i = Z i Z C o metro = Z i 1 METRO i metro i Z i

Ahora, tomemos como ejemplo la fórmula para I xx en el marco O:

I X X = i metro i ( y i 2 + z i 2 )

(Tenga en cuenta que la suma de dos términos entre paréntesis es simplemente la distancia euclidiana entre el i-ésimo átomo y el eje x, exactamente como cabría esperar).

Sustituyendo las fórmulas para y i y z i arriba tenemos

I X X = i metro i [ ( Y i Y C o metro ) 2 + ( Z i Z C o metro ) 2 ] I X X = i metro i ( Y i 2 + Z i 2 ) + i metro i ( Y C o metro 2 + Z C o metro 2 ) 2 i metro i ( Y i Y C o metro + Z i Z C o metro ) I X X = i metro i ( Y i 2 + Z i 2 ) + METRO ( Y C o metro 2 + Z C o metro 2 ) 2 METRO ( Y C o metro 2 + Z C o metro 2 ) I X X = i metro i ( Y i 2 + Z i 2 ) 1 METRO ( i metro i Y i ) 2 1 METRO ( i metro i Z i ) 2

que es exactamente la fórmula para I xx en el documento de la NASA. Análogamente para las otras fórmulas.

Curiosamente, la suma de los dos últimos términos en la fórmula final es exactamente lo que cabría esperar del teorema de los ejes paralelos .

La conclusión es que debe usar las fórmulas que se ajusten a su marco de coordenadas. Si el centro de su marco de coordenadas coincide con el centro de masa de la molécula, entonces puede usar las fórmulas más simples de Wikipedia. De lo contrario, debe usar los del documento de la NASA.

Mi programa devuelve números mucho más grandes de lo esperado (obtengo 30000 para los valores propios cuando espero algo entre 0 y 40). ¿Qué unidades debo usar para la masa atómica? Estoy usando el peso atómico estándar dado en Wikipedia, así que, por ejemplo, el nitrógeno es solo 14.0067.
Como siempre, puedes usar las unidades que prefieras siempre y cuando mantengas la coherencia. Es más fácil si anota/comenta las posibles declaraciones de variables y el código y con las unidades que usa en cualquier etapa de entrada, cálculo y salida. En particular, cuando comparas los valores propios que obtienes con lo que esperas, las unidades también deben ser las mismas.