Estrellas Dobles y Múltiples:

Resultados del Proyecto Conjunto AAT-AAGC.

Orlando José Rodríguez Santana


Desde el inicio de la Sección de Cielo Profundo y Variables de la AAGC (Agrupación Astronómica de Gran Canaria), año 1991, hubo un gran interés en la observación de estrellas dobles y múltiples. Después de algún año de experiencia en su observación, en el cual, la mayoría de nosotros conoció esta rama tan interesante de la astronomía. Gracias al II Encuentro Astronómico Canario pudimos saber que el incansable observador "Kiko" Hernández, ahora presidente, de la AAT (Agrupación Astronómica de Tenerife), también se había dedicado a observar estrellas dobles, no hacía demasiado tiempo atrás. En este encuentro surgió la idea de realizar el proyecto que estábamos llevando a cabo en la AAGC sobre el "Cálculo de las órbitas de las estrellas dobles y múltiples" de forma conjunta con la AAT.

En la observación de las estrellas dobles y múltiples lo que nos interesa es la obtención del ángulo de posición y la separación de la estrella secundaria (como tal se suele elegir a la menos brillante) frente a la primaria. (Ver fig.1)

Existen diferentes tipos de estrellas dobles :

1. Visuales:

- ópticas: no están unidas por la gravedad, sino que están en nuestra línea de visión y debido a la perspectiva desde la Tierra aparecen juntas en la visual. (Ver fig.2)

- orbitales (físicas): son aquellas que están unidas por la gravedad, tal como los planetas al Sol en el caso de nuestro sistema solar. (Ver fig.3)

2. Espectroscópicas: son parecidas a las anteriores, pero en este caso sólo se descubren debido al corrimiento Doppler de sus líneas espectrales.

 

 

 

 

 

 

 

 

Para el caso orbital, tendremos tres tipo diferentes de órbitas que se describen:

1.a Aparente: es la que nos dibuja una estrella cuando representamos a la otra en el centro de nuestro sistema de coordenadas. Es la proyección de la órbita real sobre la esfera celeste.

1.b Real: en ella se considera a una de las estrellas (a la primaria normalmente) en el foco de la elipse que describe una alrededor de la otra, según las leyes de Kepler.

- física: en ella se consideran los movimientos de ambas estrellas en torno a su centro de masas común.

A partir de las observaciones del ángulo de posición y de la separación, podremos obtener la órbita aparente. Si sólo poseemos un par de observaciones, sólo obtendremos una parte, pero podremos ya hacernos una idea del tipo que es: óptica, física,... Según podemos ver en la figura 4 de aquí a la órbita real sólo hay un paso de "perspectiva", en el cual el ángulo que nos define la distancia entre la esfera celeste y el plano de la órbita real es la inclinación i. Este ángulo poseerá dos valores posibles, debido a que la perspectiva que nos da una estrella situada con un ángulo i y otra con un ángulo -i , que poseerá la misma visión terrestre. Por tanto su valor será ±i . Esta duda podríamos resolverla si viajáramos al sistema, pero por ahora eso no es posible...

Este ángulo junto con 6 magnitudes más, son los que nos definen la posición de la órbita real tanto en el tiempo como en el espacio. Son los elementos orbitales de la estrella doble orbital.

- P, periodo; nos indica el periodo de revolución: tiempo en dar una vuelta completa. Se mide en años.

- T, tiempo de paso por el periastro, se mide en segundos de arco (").

- a, semieje mayor de la órbita real, se mide en segundos de arco (").

- e, excentricidad de la órbita real,

0<e<1

- i, inclinación de la órbita real, se mide en grados (°). Se define:

=90° para la línea de visión

<90° ángulo de posición decrece (sentido horario)

>90° ángulo de posición crece (sentido antihorario)

0<i<180°

-W , posición angular del nodo ascendente medido en el plano del cielo en el sentido de las agujas del reloj. Se mide en grados (°)

0<W <180°

-w , longitud del periastro. Se mide en grados (°). Angulo que mide en el plano de la órbita real, distancia entre el nodo ascendente hasta el punto del periastro y en la dirección del movimiento.

Debemos hacer notar que en la realidad el "plano celeste" no es un plano sino la esfera celeste, pero debido a que la distancia es muy grande se puede considerar como un plano, pues esta esfera tiene el tamaño de la distancia hasta la estrella.

Resultados de las Observaciones

A finales del año 1996 ya habíamos obtenido más de 150 observaciones (del ángulo de posición y la separación) de determinadas estrellas dobles, que debido a la falta de "tiempo observacional" los observadores implicados: "Kiko" Hernández Cabrera (AAT), Máximo Suárez Tejera (AAGC) y yo mismo; no se pudieron elegir estrellas orbitales con periodos muy cortos para así comprobar el sistema a usar para el cálculo de las órbitas. Así que entre cometa y cometa, y entre otros eventos astronómicos observábamos las estrellas dobles de la constelación que mejor se podía observar esa noche. A pesar de eso dos observaciones de "Kiko" poseían dos estrellas orbitales cuyos elementos orbitales pude encontrar en el "Sky Catalogue 2000.0". Eran S 422 en Tauro y V Aquari. Precisamente estrellas que ya se sabían sus elementos orbitales eran lo que buscábamos para comprobar nuestros cálculos con los elementos orbitales antiguos e intentar calcular los elementos orbitales nuevos. Es precisamente en este punto donde hemos tenido mayores problemas y ha hecho que nos retrasemos tanto en la obtención de los primeros resultados.

Las observaciones de "Kiko" Hernández Cabrera fueron realizadas con un Meade de 255 mm de diámetro a f6.3, en donde usábamos un ocular de 12.5 mm (128,52 aumentos) reticulado con círculo graduado adaptado realizado por el propio "Kiko". Cada división del retículo iluminado tenía 4" por tanto nuestra incertidumbre es ± 2" y la incertidumbre en el ángulo de posición es de 0.5°. Si calculáramos la resolución teórica del telescopio mediante el criterio de Rayleigh sería de 0.55".

Gracias a la Sección de Astroinformática y Multimedia de la Agrupación Astronómica de Gran Canaria (SECAM-AAGC), representada por Victor R. Ruíz realizamos un programa de ordenador (después de haber hecho los cálculos "a lápiz y papel" para comprobar las ecuaciones usadas) mediante el cual, gracias a los elementos orbitales obtenemos la posición de la estrella orbital en función del tiempo. Esto nos podría indicar para la época deseada (a pesar de los grandes periodos usados) la posición exacta de la estrella secundaria.

A continuación expongo la tabla de los datos de las estrellas tiene los números de catálogo siguientes y los datos observacionales disponibles:

La primera recibe los siguientes nombres 55, V Aquari, ADS (15971) (Aitken Double Stars) y S 2909 Aqur (Struve).

Año

Ang Pos (° )

Sep (")

Fuente

1991

211

2

Observación Kiko

1987

211

1.9

Jose Luis Comellas

1935

289

2.6

Aster

1940

263

2.5

Sky Catalogue 2000.0

(Calculado)

1960

227

2.0

SC 2000.0

1980

192

1.8

SC 2000.0

2000

180

2.1

SC 2000.0

2020

203.7

2.7

SC 2000.0

La siguiente estrella fue, como ya hemos dicho: 422 en Tauro, que resultó ser la variable V711 de esta constelación, según pude comprobar en el Sky Catalogue 2000.0. Era una estrella variable del tipo RS-Canum Venaticorum, que son estrellas binarias de rápida rotación tipo F o K, subgigantes o gigantes. Sus periodos varían desde fracciones de días hasta bastantes meses. Su variación es producida por la actividad de manchas y en algunos casos (como este) por eclipses. Muchos tipos RS-variables son fuentes de radio y de rayos X (como este caso) debido a perturbaciones coronales. La magnitud conjunta era de 5.81, siendo en el máximo 5.70 y el mínimo 5.81 (visual). Su periodo es de 2.83782 días. Sus tipos espectrales son G5 IV y K1 IV.

Año

Ang Pos (° )

Sep (² )

Fuente

1993

265

6

Observaciones Kiko

1987

265

6.3

J.L.Comellas

1940

257

6.5

Sky Catalogue 2000.0

1950

261

6.6

S.C.2000

1980

266

6.6

S.C.2000

2000

274

6.7

S.C.2000

Los elementos orbitales de ambas estrellas según el Sky Catalogue 2000.0 se pueden ver en la siguiente tabla:

Estrella

P(años)

T(años)

a(² )

e

i

w

W

S 422 Tau

2101

1900.0

8.023

0.18

32.1

152.0

92.3

V Aqu

856.0

1957.6

5.055

0.49

131.2

55.12

310.2

Utilizando el programa de ordenador que realizamos los datos para S 422 Tau calculados a partir de los elementos orbitales del Sky Catalogue 2000.0, en este calculo usamos un intervalo de tiempo de 50 años desde el año 100 a.J. Para V Aquari usamos un intervalo de tiempo de 25 años desde el año de nuestra era, de estos datos obtenemos las siguientes gráficas:

Como se puede ver en las siguientes dos gráficas las observaciones concuerdan bien con lo predicho por el cálculo mediante la órbita.

Estos han sido los cálculos para las estrellas que ya poseían la órbita calculada, pero además obtuvimos interesantes datos para algunas de las estrellas de las que obtuvimos alguna observación, como ejemplo tenemos las siguientes observaciones seleccionadas de las cerca de dos centenares observadas.

Listado de las estrellas seleccionadas de las 171 observadas, en donde se puede apreciar movimientos que podrían indicar que son orbitales
                 

a

d

Nombre

m1

m2

Años

Sep"

Ang Pos(º)

Fuente

hhmm.s

ºº ''

Constelación

           
                 

159.5

2101

S196 AB Ari

9.5

12

1832

2.4

56

 
         

1958

2.2

48

 
   

AC

 

10.6

1832

39.5

167

 
         

1918

30.9

163

 
         

1993

28

162

Kiko

   

AD

 

6.1

1862

183.7

1

 
         

1914

188.1

4

 
         

1993

198

8

Kiko

255.6

2652

S 326 AB Ari

7.6

9.8

1831

216

9

 
         

1976

220

5.9

 
         

1993

222

7

Kiko y Orlando

310.9

-1108

Bpm? AB Eri

7.2

8.8

1894

144.2

57

 
         

1907

143.7

56

 
         

1993

139

52

Kiko

319.4

322

Bpm?AB Cet

4.8

9.3

1853

271.6

153

 
         

1909

268.7

157

 
         

1993

247

163

Kiko

340.9

-1237

S436 Eri

7.5

9.1

1832

30.2

232

 
         

1935

40.1

236

 
         

1993

45

236

Kiko

400.8

1812

Bpm? Tau

5.9

9.2

1909

161.3

277

 
         

1921

164.2

278

 
         

1993

175

279

Kiko

404.7

2205

Bpm? Tau

4.4

9.4

1853

138.7

187

 
         

1921

137.2

190

 
         

1993

135

191

Kiko

409.5

-1729

Bpm? AB Eri

8.2

8.7

1894

126.3

139

 
         

1908

125.2

139

 
         

1993

117

135

Kiko y Orlando

414.4

-1015

S516 AB Eri

5

8

1836

6.3

154

 
         

1967

6.4

146

 
         

1993

9

145

Kiko

   

AC

 

9.5

1911

148.9

159

 
         

1993

135

157

Kiko

427

1908

S546 Tau

7.9

9.7

1836

6.6

190

 
         

1954

6.8

184

 
         

1993

7

179

Kiko

439.1

752

Bpm? AB Tau

5.5

10.4

1884

69.2

247

 
         

1931

71.6

249

 
         

1993

72

251

Kiko

   

AC

 

9

1884

299.4

317

 
         

1921

297.8

316

 
         

1993

292

316

Kiko

451.4

1850

Bpm? Tau

5.1

10.4

1909

174.5

303

 
         

1922

175.8

303

 
         

1993

180

305

Kiko

456.4

-510

S,h 48? Eri

5.5

9.1

1821

65.9

75

 
         

1913

67.3

75

 
         

1993

63

73

Kiko

507

928

Bpm? AB Ori

6.2

10

1852

128.2

255

 
         

1909

124

264

 
         

1993

127

280

Kiko

   

AC

 

?

1904

402.7

268

 
         

1993

405

274

Kiko

1859.1

1337

S 2424 Aql

5.2

8.7

1831

118.7

242

 
         

1957

18.5

286

 
         

1991

18

294

Kiko

1902.2

846

S 2436 AB Aql

8.3

9.1

1830

34.6

309

 
         

1925

31.9

313

 
         

1991

28

310

Kiko

1902.7

-42

S2434 AB Aql

8.5

8.5

1831

25.6

147

 
         

1955

24.3

106

 
         

1991

24

95

Kiko

2207.1

34

S2862 Aqr

8.2

8.6

1828

2.3

104

 
         

1959

2.5

98

 
         

1991

3

97

Kiko

2226.6

-1645

S,h 345 AB Aqr

6.4

6.6

1823

10.1

303

 
         

1976

3.1

334

 
         

1991

2

352

Kiko

2243.1

-819

S2935 AB Aqr

6.9

7.9

1831

2.6

313

 
         

1959

2.3

308

 
         

1991

2

305

Kiko

2247.7

-1403

S2943 Aqr

5.8

9

1783

35.6

110

 
         

1959

23.7

121

 
         

1991

22

126

Kiko

2331.4

1613

S2021 AB Peg

7.9

9.1

1830

8.6

308

 
         

1951

8.6

308

 
         

1991

8

309

Kiko

   

AC

 

10.1

1913

118

24

 
         

1931

118.7

24

 
         

1991

112

26

Kiko

2347.9

1703

S3041 AB Peg

8

7.8

1900

64.4

351

 
         

1936

61

352

 
         

1991

52

355

Kiko

   

AC

 

9.2

1900

68.1

351

 
         

1921

66.4

351

 
         

1991

54

356

Kiko

Quizás en estos casos otro tipo de efectos (como la precisión del eje de la tierra, fallos observacionales,....) puedan explicar el movimiento que parece ser orbital. Hace falta una informatización de los datos para poder utilizar toda la información posible para ir descartando las diferentes hipótesis, para cada caso en particular.

Un ejemplo de la información que podríamos obtener es el siguiente caso en donde si suponemos que el movimiento de las dos estrellas compañeras es debido al movimiento propio , de ellas obtendríamos los siguientes resultados:

Como se puede ver con no demasiados medios se pueden llevar a cabo interesantes trabajos, solo hay que armarse de paciencia a la hora de obtener información. Como ejemplo la mayor parte del tiempo lo he perdido buscando información para calcular la órbita de una estrella a partir de 4 datos observacionales, y eso a pesar de que disponía de las mejores bibliotecas para la astronomía y astrofísica en Canarias, pues he estado estudiando en La Laguna, donde esta la sede del IAC. Pero como se sabe hoy día la observación de las estrellas dobles y múltiples pertenece al pasado, hoy tan solo tiene interés científico la observación de planetas extrasolares. Creo que es importante saber el número de estrellas que son sistemas múltiples, para así tener una idea de la posibilidad de que se puedan formar sistemas planetarios como el nuestro, pues en un sistema múltiple es más dificil la formación de planetas, según estudios gravitacionales.

No ha sido necesario el uso de CCDs, pero podrían ayudar y simplificar quizás las tomas de observaciones. Como se demuestra en el articulo, llegamos a la conclusión de que es necesario la informatización de los datos que se poseen y de las futuras observaciones. En esto es en lo que estoy trabajando ahora con ayuda de Victor R. Ruiz. Esperamos tener el CREM (Catálogo renovado de estrellas múltiples para antes del fin de 1998 ó 1999). A esto hay que añadir, sin ninguna duda, el trabajo realizado por el satélite Hipparcos (High Precision Parallax Collecting Satellite) el cual, él solito y gracias a que ha estado fuera de la atmósfera, ha estado investigando aquellas sistemas estelares que poseían diferencias de 4 magnitudes o 0.1 segundos de arco entre una y otra. Hipparcos ha resuetlo 23900 sistemas múltiples (la tercera parte de todos los que se conocían hasta ahora), de los cuales cerca de 300 son nuevos, 235 son sistemas orbitales, 2900 sospechosas de ser binarias astrométricas (nuevas) y 8500 sospechosas de ser múltiples. Pero Hipparcos no ha realizado todo el trabajo que hace falta, y es aquí donde los no profesionales podremos aportar nuestro granito de arena , .... o debemos decir estrellita??

Cálculo del movimiento interanual suponiendo que este es debido al movimiento del grupo de estrellas
               

Años

Sep(")

AP(grados)

x

y

     

1884

69.2

247

-27.0379519

230.9021898

     

1921

71.6

249

-25.6584659

236.2703566

     

1993

72

251

-23.4402097

238.7154135

     

1884

299.4

317

218.96994

-276.172877

     

1921

297.8

316

214.2220605

-279.525797

     

1993

292

316

210.0498377

-292.988687

     
               
               
Calculo de la recta hacia 1993, para la primera estrella        
Pendiente x=  

0.037283407

         
Ordenada en el origen x =

-97.2798912

         
Calculo teórico x=  

-22.97406

         
Pendiente y=  

0.145085589

         
Ordenada en el origen y =

-42.4390595

         
Calculo teórico y=  

246.71652

         
Diferencia en x =  

-0.46614917

         
Diferencia en y =  

-8.00110546

         
Diferencia total=  

8.014673

         
Esto equivaldría a un movimiento interanual de :