detección de Planetas Extrasolares

 

Objetivos específicos:

·         Determinar la masa de algunos compañeros oscuros de estrellas y discutir si se pueden clasificar como planetas extrapolares.

 

 

Conocimientos previos:

 

Leyes de Kepler.

Métodos de detección de planetas extrapolares basados en la observación de la variación de la velocidad radial de una estrella.

Dependencia de la curva de velocidad radial con los parámetros orbitales.

 

 

Descripción:

En base a datos reales de velocidad radial contra tiempo para dos sistemas extrasolares, se hallarán los períodos orbitales, distancias medias (semiejes mayores) y masas de dos objetos subestelares en torno a las estrellas 47UMa y 16CyB respectivamente.

De acuerdo a las masas halladas, se evaluará si pueden ser planetas o enanas marrones.

 

Materiales:

§         Tablas y gráficos con las velocidades radiales de las estrellas 16CyB y 47UMa obtenidas a partir de las observaciones (ver al final)

§         Software AVE para determinación de períodos (no imprescindible). Se puede bajar de: http://www.astrogea.org/soft/ave/introave.htm o obtener de aquí.

 

 

Procedimiento

1-      Determinación del período del objeto desconocido.

La determinación del período se podrá hacer por dos métodos:

i) determinación visual a a partir del estudio de la gráfica.

ii) determinación usando el programa AVE.

Las tablas antes mencionadas se proporcionarán en formato digital, bajo los nombres de 16cyb.djm y 46UM.djm respectivamente. La extensión djm es necesaria para ingresar los archivos al software AVE que calculará los períodos.

Para determinar el período de ambos planetas se deberán al cargar los datos, analizar la gráfica de velocidad radial contra tiempo y estimar visualmente un período.

Luego, utilizando la opción search period encontrar un valor numérico para el mismo.

Una vez hallado:

a)      compararlo con el valor tentativo hallado a partir del estudio de la gráfica.

b)      Discutir las posibles causas de error.

 

 

2. Determinar el semieje mayor de la órbita.

 

Previamente se deberá analizar si la órbita es circular o elíptica, lo cual se deduce de las gráficas. Se discutirá como la forma geométrica de la órbita influye en la simetría o asimetría de la gráfica.

Una vez obtenido este dato se calculará a partir de la tercer ley de Kepler el semieje de la órbita.

Según esta ley:

                       

T2  = k. a3                                                      (Ec. 1)

 

donde T es el período y a el semieje.

Si T está dado en años y a en Unidades Astronómicas, k es simplemente el inverso de la masa de la estrella (en masas solares) que es un dato conocido.

 

Masas estelares (en masas solares):

 

Tabla 3: masas estelares (en masas solares)

16 CyB 46 UMa
1.0 1.1

 

 

3. Determinar la masa mínima del compañero subestelar.

 

Hablamos de masa mínima porque tenemos la indeterminación del ángulo que forma el plano principal de la órbita respecto a nuestra línea de visual.

Las proyecciones de la velocidad no son las mismas si vemos el sistema ‘de canto’ que ‘de frente’ por lo cual hay un factor sin i (seno de la inclinación) que cuenta en la determinación de la masa.

De todas maneras calcularemos la masa mínima posible, sin tener en cuenta la inclinación.

Para ello consideramos la ecuación del baricentro del sistema:

Las coordenadas del baricentro B serán:

 

B = Mp xp + M*  x*                             (Ec. 2)

 

donde  Mp es la masa del objeto subestelar (nuestra incógnita)

            xp  es la posición del objeto subestelar respecto al baricentro

           M* la masa de la estrella (conocida)

            x* la posición de la estrella

 

Si derivamos respecto al tiempo:

 

0 = Mp  vp + M*  v*                       (Ec. 3)

 

donde vp es la velocidad del objeto subestelar, que debemos determinar y v*  es la velocidad de la estrella que obtenemos a partir de la semiamplitud de la gráfica de velocidad radial contra tiempo (el menos para un movimiento circular, donde el comportamiento de la velocidad radial será sinusoidal, esto es válido).

Los valores reales de la semiamplitud Doppler de las gráficas son (velocidad en m/s):

 

 

16 CyB

46 UMa

 

 

44

46

Tabla 4: semiamplitud en Doppler (en m/s)

                    

Estas velocidades deberán ser expresadas en UA/año.

¿Cómo calculamos la velocidad del objeto subestelar?

Existe una ecuación para el movimiento circular que plantea:

 

vp2 = m/a                                     (Ec. 4)

 

donde m= G.M*                          (Ec. 5)

 

El valor de G en caso de trabajar en UA y años es 4p2.

Por lo tanto:

 

vp=  (4.p2 . M* / a )                   (Ec. 6)

 

Conociendo este valor, podemos despejar en la Ec 3 el valor de la masa del objeto subestelar.

 

Mp = M*  v* / vp

 
 

 

 

 


Se deberá:

a)      calcular la masa del objeto subestelar con la semiamplitud que estimen de la gráfica.

b)      calcular la masa del objeto subestelar con los datos proporcionados en la Tabla 4.

c)      discutir las causas de error.

d)      evaluar de acuerdo a la masa hallada si el objeto subestelar es un planeta

 

 

Temas de discusión grupal:

Todos los datos obtenidos podrán compararse con los valores reales obtenidos por los investigadores, publicados en Sky & Telescope (Marzo de 1998).


 

Tabla 1.

Velocidad radial (en m/s) contra tiempo (en años) para 16CyB

 

Tabla 2.

Velocidad radial (en m/s) contra tiempo (en años) para 46UMa

 


 

 

 


16 Cy B

 

Año     Vel Rad.

 

1987.5     0.475

1989.8   24.240

1990.2   39.264

1990.5   45.259

1991.4  -28.429

1992.6   29.772

1992.7   40.823

1993.3  -56.980

1993.5  -19.744

1993.6  -41.922

1994.6     0.185

1994.7   29.379

1995.4  -56.353

1995.4  -35.216

1995.9  -36.313

1996.3  -15.225

1996.5    -6.195

1996.6     3.862

1997.2   35.974

1997.3   -1.304

1997.3   -29.512

1997.3      -53.690

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

46 UMa

 

Año     Vel Rad

 

1987.5  -30.959

1987.6  -18.659

1989.7   15.804

1989.9  -32.279

1990.9  -20.118

1991.7   49.172

1992.0   68.162

1992.7   25.572

1992.9   87.665

1993.0  -35.996

1993.1  -30.409

1993.2  -18.117

1993.7  -37.187

1993.8  -31.601

1995.0   51.077

1995.9   30.846

1996.0   17.404

1996.1   28.552

1996.2  -11.705

1996.3  -15.070         

1996.8  -41.980

1996.9  -31.914

1997.1  -27.461

1997.2  -11.808

1997.3    -62.218