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:
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?
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:
La transformación de O a O com es simplemente traducción por [-X com , -Y com , -Z com ]:
Ahora, tomemos como ejemplo la fórmula para I xx en el marco O:
(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
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.
frijol
Adán Zalcman