Este fin de semana toca, una vez más, práctica numérica de Fundamentos de Astrofísica. Una de las tareas que hay que hacer es, a partir de la curva de velocidad radial de una estrella binaria espectroscópica, hallar la velocidad del centro de masas.
Según el método de Lehmann-Filhes, lo que hay que hacer es trazar una linea horizontal en la gráfica de velocidad radial, de modo que las áreas encerradas por la curva por encima y por debajo de esta linea horizontal sean iguales. Parece trivial, pero en una licenciatura como la nuestra (CC. Físicas en la UCM) no nos enseñan a manejar programas de cálculo numérico lo suficiente como para resolver este tipo de problemas.
Yo he seguido un proceso iterativo. En primer lugar, he puesto los datos de la curva en una tabla, los he representado, y he usado la función de integración de Qtiplot (un programa clon de Origin) para realizar la integral de la curva por el método del trapecio.
Como se puede ver, el valle de la curva de velocidad radial es mucho más extenso que la cresta, y es por lo que al final las áreas suman negativo. La estrategia consiste en desplazar la curva hacia arriba hasta que la integral de cero, y la distancia que lo hayamos desplazado, cambiada de signo, equivaldrá a la velocidad del centro de masas en km/s. Es decir, si para que la integral se anule tenemos que desplazar la gráfica 8 km/s hacia arriba, la velocidad del centro de masas será -8 km/s.
En este caso, el desplazamiento adecuado fue de -13.3 km/s.
Son bienvenidos comentarios y sugerencias sobre otras formas de hacer el cálculo.
Recent Comments