Calibración Fotométrica en Huancayo
Calibración Fotométrica en Huancayo
EN EL OBSERVATORIO DE HUANCAYO
RESUMEN
ABSTRACT
8” de apertura (Celestron Celestar 8), de razón focal dark current fueron tomadas al comienzo y al final
f/10, instalado sobre una montura ecuatorial tipo de la noche.
horquilla dentro de la cúpula automática del
observatorio. La cúpula posee aproximadamente 2.5 Para el cálculo de los coeficientes de extinción
m de diámetro por 3 m de alto y alberga fueron observadas un total de 7 estrellas de tipo A0
confortablemente hasta tres observadores. Una (o aproximado) en varias masas de aire y para la
cámara CCD (SBIG ST-7E)1 fue adaptada al foco del determinación de los coeficientes de transformación
telescopio donde son formadas las imágenes fue escogido el aglomerado abierto M41 (ver Figura
estelares. La cámara posee un chip Kodak KAF- 1) localizado en la constelación del Can Mayor. En
0401E (765 x 510 píxeles) con un tamaño de píxel de cada filtro, fueron tomadas imágenes individuales
9µ x 9µ y que para la distancia focal del telescopio para cada estrella A0 y una imagen respectiva para el
(F = 2032 mm) proporciona una escala de placa de aglomerado M41. El campo del CCD cubre
0.91”/píxel. El campo total de cielo cubierto por el aproximadamente una tercera parte del área del
CCD es de 12” x 8”. El sistema de enfriamiento aglomerado, que fue suficiente para tener algunas
termoeléctrico de la cámara trabaja a 25° por debajo estrellas brillantes (indicadas por círculos en la
de la temperatura ambiente, la cual fue de aprox. 10° figura) para hacer el análisis.
C al comienzo de la noche.
La Tabla 1 lista las estrellas de tipo A0 observadas.
Las medidas fueron realizadas usando los filtros U, B En la columna (1) está el número de observaciones
y V en las bandas ultravioleta, azul, y visual, realizadas. Todas las estrellas fueron observadas en
respectivamente, del sistema fotométrico Johnson- una única secuencia de filtros a una determinada
Morgan. Los filtros fueron instalados en una rueda masa de aire, excepto las estrellas SAO 144150
automática (CFW8) que permite un rápido (observaciones 3 y 4) y SAO 146044 (observaciones
intercambio de filtros durante las observaciones. La 6 y 7) que fueron observadas en dos masas de aire
cámara CCD como la rueda de filtros son diferentes. En la columna (2) está el nombre de la
controlados a través del software de control estrella respectiva según el catálogo SAO2; en las
CCDOPS (ver. 4.12-NT) instalado en una columnas (3) y (4) están indicadas la coordenadas
computadora portátil y donde la conexión es hecha ecuatoriales (2000) para cada objeto; la columna (5)
usando el puerto paralelo. Las imágenes digitales indica el tipo espectral de cada estrella según la base
generadas son almacenadas en el disco duro de la de datos astronómica SIMBAD3; la columna (6)
portátil para posteriormente ser transferidas a una muestra el valor de seeing medido para cada objeto
computadora donde son procesadas. en segundos de arco. La calidad media del cielo
(seeing) para la noche de observación calculada en
Se presentaron algunos inconvenientes en la base al ancho total a media altura (fwhm, full with at
observación de estrellas en masas de aire mayores. half-maximum) de los perfiles estelares fue 3.48”,
Por un lado el área efectiva disponible de la abertura con un valor mínimo de 2.26” y un valor máximo de
de la cúpula limita la observación a alturas pequeñas 5.06”. Las columnas (7) y (8) indican la fecha y hora
sobre el horizonte y, por otro, el poco espacio local de la observación; y finalmente la columna (9)
disponible en la montura tipo horquilla para la muestra la masa de aire (X) calculada, para la
colocación de la cámara CCD más la rueda de filtros posición del objeto al momento de ser observado
también limita el mismo tipo de observación. (ver siguiente sección).
El sistema telescopio, cámara y rueda de filtros fue La Tabla 2 lista las estrellas observadas en el
contrapesado convenientemente para permitir una aglomerado M41. La información mostrada en
fácil localización de las estrellas a estudiar. Un semejante a la indicada en la Tabla 1. Los nombres
nuevo sistema de contrapesos bastante simple en su de las estrellas indicadas en la columna (2) se
diseño fue exitosamente probado. La localización de refieren al catálogo HD4 [2,3].
los objetos estelares fue hecha manualmente
utilizando un buscador de 7 x 50 (pequeño telescopio Los tiempos de integración (TI) para todas los
de 7 aumentos y 50 cm de apertura) y por estrellas listadas en la Tabla 1 fueron de 1 segundo
comparación con cartas estelares electrónicas.
Imágenes de calibración de flat field en cada filtro y
2
Smithsonian Astrophysical Observatory Star Catalog
1 3
[Link] (Santa Barbara Instrument Group, [Link]
4
Inc.) Henry Draper Catalogue and Extension
5
4
1 2
6
Figura 1. Aglomerado Abierto M41. El área mostrada corresponde a una región de 24’ x 24’. Las estrellas utilizadas para el análisis de
los coeficientes de transformación están indicadas con círculos y números. La numeración se corresponde a la indicada en la columna 1
de la Tabla 2.
en los tres filtros. Las estrellas de tipo A0 escogidas REDUCCIÓN Y ANÁLISIS DE DATOS
eran más brillantes que 6ta. magnitud para que
integraciones rápidas puedan realizarse. Problemas El proceso de reducción y análisis de las imágenes
en el motor de acompañamiento sideral y en el digitales fue realizado en el sistema operativo Linux
alineamiento polar actual no permiten hacer RedHat 6.2, usando el ambiente IRAF5 [4,5].
integraciones largas para objetos a distancias
cenitales grandes. El TI para la imagen de M41 en U 5
IRAF (Image Reduction and Analysis Facility), es un
fue de 5 s y en B y V fue de 1 s. Un tiempo mayor (5 sistema de software de propósito general para la reducción
s) fue necesario para la imagen de M41 en U para y análisis de datos astronómicos. IRAF es escrito por el
asegurar que un número suficiente de fotones sean IRAF programming group ([Link] en NOAO
captados. Esto a pesar de la pérdida de cualidad en la (National Optical Astronomy Observatories) en Tucson,
imagen por los problemas de seguimiento ya Arizona. NOAO is operado por AURA (Association of
mencionados. Universities for Reserach in Astronomy, Inc.) bajo
acuerdo cooperativo con la NSF (National Science
Foundation).
Tabla 1. Lista de estrellas observadas de tipo A0.
observación estrella A. R.2000 Dec.2000 T. E. filtro seeing día hora X
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
1 SAO 110543 02 28 09.54 +08 27 36.2 B9III U 4.16 09/18/2001 [Link] 1.0226
B 2.88 09/18/2001 [Link] 1.0240
V 2.85 09/18/2001 [Link] 1.0271
2 SAO 122754 17 47 53.56 +02 42 26.2 A0V U 5.06 09/17/2001 [Link] 2.0302
B 3.95 09/17/2001 [Link] 2.0004
V 3.78 09/17/2001 [Link] 1.9693
3 SAO 144150 20 11 18.29 -00 49 17.3 B9.5III U 4.10 09/17/2001 [Link] 1.0736
B 3.91 09/17/2001 [Link] 1.0802
V 4.06 09/17/2001 [Link] 1.0832
4 U 3.33 09/18/2001 [Link] 1.9111
B 3.39 09/18/2001 [Link] 1.8903
V 3.29 09/18/2001 [Link] 1.8561
5 SAO 144810 20 47 40.55 -09 29 44.8 A1V U 2.83 09/17/2001 [Link] 1.0783
B 4.50 09/17/2001 [Link] 1.0778
V 3.29 09/17/2001 [Link] 1.0772
6 SAO 146044 22 21 39.38 -01 23 14.4 A0V U 3.69 09/17/2001 [Link] 1.0284
B 3.90 09/17/2001 [Link] 1.0282
V 3.34 09/17/2001 [Link] 1.0280
7 U 3.57 09/18/2001 [Link] 1.1870
B 2.64 09/18/2001 [Link] 1.1826
V 2.71 09/18/2001 [Link] 1.1779
8 SAO 55397 02 15 56.28 +33 21 32.0 AOV U 3.79 09/18/2001 [Link] 1.1186
B 2.53 09/18/2001 [Link] 1.1207
V 2.26 09/18/2001 [Link] 1.1236
9 SAO 55427 02 17 18.87 +33 50 49.9 A1Vnn U 4.15 09/18/2001 [Link] 1.1382
B 3.07 09/18/2001 [Link] 1.1411
V 2.96 09/18/2001 [Link] 1.1441
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14)
SAO 110543 4.26 -0.045 -0.125 1.0246 ± 0.0023 -11.570 ± 0.004 12 -13.613 ± 0.001 14 -13.917 ± 0.001 15 18.177 ± 0.001 -0.349 ± 0.001 -2.168 ± 0.004
SAO 122754 3.74 0.045 0.045 2.0000 ± 0.0305 -11.349 ± 0.005 14 -13.752 ± 0.001 15 -14.242 ± 0.001 15 17.982 ± 0.001 -0.445 ± 0.001 -2.358 ± 0.005
SAO 144150 3.29 -0.08 -0.16 1.0790 ± 0.0049 -12.666 ± 0.002 15 -14.585 ± 0.001 15 -14.934 ± 0.001 15 18.224 ± 0.001 -0.429 ± 0.001 -2.079 ± 0.002
1.8858 ± 0.0278 -12.260 ± 0.003 15 -14.305 ± 0.001 15 -14.653 ± 0.001 15 17.943 ± 0.001 -0.428 ± 0.001 -2.205 ± 0.003
SAO 144810 3.77 0.01 0.04 1.0778 ± 0.0006 -11.975 ± 0.003 11 -14.095 ± 0.001 15 -14.411 ± 0.001 15 18.181 ± 0.001 -0.306 ± 0.001 -2.08 ± 0.003
SAO 146044 3.88 -0.055 -0.095 1.0282 ± 0.0002 -11.971 ± 0.003 12 -14.020 ± 0.001 15 -14.313 ± 0.001 15 18.193 ± 0.001 -0.348 ± 0.001 -2.144 ± 0.003
1.1825 ± 0.0045 -11.935 ± 0.003 11 -13.924 ± 0.001 13 -14.225 ± 0.001 14 18.105 ± 0.001 -0.356 ± 0.001 -2.084 ± 0.003
SAO 55397 5.28 0 -0.03 1.1210 ± 0.0026 -10.256 ± 0.008 7 -12.41 ± 0.002 7 -12.879 ± 0.002 13 18.159 ± 0.002 -0.469 ± 0.003 -2.184 ± 0.008
SAO 55427 4 0.02 0.02 1.1411 ± 0.0029 -11.403 ± 0.004 9 -13.637 ± 0.001 15 -14.072 ± 0.001 14 18.072 ± 0.001 -0.415 ± 0.001 -2.214 ± 0.004
Después de la corrección rutinaria de darks y flats (4). La masa de aire escogida para una secuencia de
sobre las imágenes, se procedió al cálculo de las medidas UBV para un determinado objeto es el
masas de aire para cada imagen. promedio de las masas de aire de cada medida en
cada filtro y está indicada en la columna (5) con su
Para el correcto cálculo de las masas de aire, es respectivo desvío patrón. Las magnitudes
necesario actualizar los valores correctos de la fecha instrumentales calculadas (u, b y v) con sus
de observación, hora local (local time, LC), tiempo respectivos errores y la apertura escogida en cada
sideral (sideral time, ST) y coordenadas ecuatoriales medida están indicados en las columnas (6) a (11).
con su época respectiva del objeto observado. Con Las aperturas están indicadas en píxeles. Las
estos valores y conociendo las coordenadas cantidades de interés (V – v), (B – V) – (b – v) y (U –
geográficas del sitio de observación, la masa de aire B) – (u – b) con sus respectivos errores son
puede ser calculada, por ejemplo, utilizando la rutina calculadas en las columnas (12) a (14).
setairmass6 de IRAF. Los valores de la columna (10)
de las Tabla 1 y 2 fueron calculados de esta manera. La Figura 2 muestra un gráfico (V – v) vs. X para las
Esta rutina está basada en la ec. (5). estrellas tipo A0 de nuestra muestra. El ajuste lineal
según la ec. (13) proporciona un valor para el
El siguiente paso es realizar fotometría de apertura coeficiente primario de extinción K′v = 0.2356 ±
sobre cada objeto estelar para obtener las magnitudes 0.0410. Claramente podemos notar una falta de
instrumentales (u, b y v). Cuando el analizador es un valores para masas de aire intermediarias (X = 1.2 –
CCD, como es nuestro caso, la medida fotométrica 1.8). Básicamente esto se debió al poco tiempo
es la suma de ADUs (unidades analógico digitales) disponible para hacer el seguimiento de una misma
de cada uno de los píxeles dentro de la apertura estrellas en varias masas de aire, así como a las
circular definida para realizar la fotometría. Un paso complicaciones en el proceso manual de apuntar a un
importante es la definición de cual es la apertura objeto dado que se ven disminuidas cuando se
correcta para el cálculo de las magnitudes estelares. observa hacia el cenit con el equipo disponible.
Usualmente, fotometría de apertura es realizada para Especial cuidado debe tenerse en futuras medidas
una serie de aperturas concéntricas de radios cada para tratar de cubrir la mayor cantidad de masas de
vez mayores centradas en el pico de intensidad del aire posible. De manera análoga las Figuras 3 y 4
perfil estelar. En nuestro caso la fotometría fue grafican (B – V) – (b – v) vs. X y (U – B) – (u – b)
realizada para radios de apertura desde 2 (= 1.8”) vs. X, según las ecs. (14) y (15) respectivamente. Los
hasta 15 (= 13.7”) píxeles, donde la magnitud valores calculados por los ajustes lineales para los
instrumental estelar es calculada en cada apertura. coeficientes de extinción K′bv y K′ub son 0.0696 ±
Una curva de crecimiento de magnitudes es útil para 0.0484 y 0.1742 ± 0.0593, respectivamente. Los
investigar el comportamiento de la medida en valores encontrados son bastante razonables
función del radio de apertura y analizar también el indicando la tendencia general de considerar una
error esperado [6,7]. La apertura escogida estuvo en mayor corrección atmosférica en el visual que para
función del menor error obtenido para la magnitud los índices de color. En particular, K′ub parece no ser
que sea coherente con la mayor razón señal-ruido despreciable.
medida. La fotometría de apertura fue realizada
usando la rutina phot7 de IRAF. 18.25
18.20
18.10
La Tabla 3 muestra la información obtenida del
V-v
-2.05
-0.28
-0.30
-2.10
-0.32
-0.34 -2.15
(B - V) - (b - v)
(U - B) - (u - b)
-0.36
-2.20
-0.38
-0.40 -2.25
-0.42
-2.30
-0.44
-0.46 -2.35
-0.48
-2.40
1.0 1.2 1.4 1.6 1.8 2.0 1.0 1.2 1.4 1.6 1.8 2.0
X X
Figura 3. Ajuste lineal (B – V) – (b – v) vs. X para la obtención Figura 4. Ajuste lineal (U – B) – (u – b) vs. X para la obtención
de K′bv. El valor calculado es K′bv = 0.0696 ± 0.0484. de K′ub. El valor calculado es K′ub = 0.1742 ± 0.0593.
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13)
HD 49150 8.5 -0.07 -0.15 -6.765 ± 0.023 5 -8.932 ± 0.016 5 -9.515 ± 0.012 6 18.364 ± 0.066 -0.497 ± 0.061 -2.009 ± 0.070
HD 49126 7.28 0.58 0.2 -7.137 ± 0.017 5 -9.663 ± 0.01 5 -10.755 ± 0.006 8 18.384 ± 0.065 -0.356 ± 0.059 -2.018 ± 0.067
HD 49106 8.95 -0.04 -0.16 -6.268 ± 0.027 4 -8.409 ± 0.021 4 -8.999 ± 0.016 5 18.298 ± 0.067 -0.474 ± 0.064 -1.993 ± 0.072
HD 49049 8.61 -0.06 -0.21 -6.728 ± 0.02 4 -8.844 ± 0.018 5 -9.431 ± 0.013 6 18.390 ± 0.066 -0.491 ± 0.062 -2.018 ± 0.069
HD 49023 8.4 -0.05 -0.23 -6.974 ± 0.019 5 -9.036 ± 0.015 5 -9.528 ± 0.011 5 18.277 ± 0.066 -0.386 ± 0.061 -1.984 ± 0.068
HD 49091 6.9 1.53 1.74 -4.898 ± 0.068 3 -9.256 ± 0.013 5 -11.12 ± 0.005 9 18.369 ± 0.065 -0.178 ± 0.059 -2.310 ± 0.094
18.32
coeficientes de extinción, así como para los
18.28 coeficientes de transformación al sistema patrón de
magnitudes son bastantes razonables. Se hace
18.24
necesario continuar con este tipo de medidas para
18.20 hacer una caracterización profunda de la extinción
-0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 atmosférica en el sitio de observación y a su vez
B-V analizar sus variaciones estacionales. Nuevas
medidas también podrán discriminar la estabilidad de
Figura 5. Ajuste lineal (V – v0) vs. (B – V) para la obtención de las constantes de transformación halladas para el
ε y ζv. Los valores calculados son ε = 0.0268 ± 0.0337 y ζv =
instrumento utilizado.
18.3388 ± 0.0229.
-0.1
AGRADECIMIENTOS
-0.2
Los autores están especialmente agradecidos al
personal del Observatorio de Huancayo por su
(B - V) - (b - v)0
-0.3
colaboración para la realización de este trabajo. Este
trabajo forma parte del proyecto “Observaciones
-0.4
Astronómicas con CCD” financiado por la Facultad
de Ciencias – UNI e “Implementación del
Laboratorio de Tratamiento de Imágenes
-0.5
Astronómicas CCD, FC-UNI”, financiado por el
CONCYTEC (C.S. No. 018-2002-CONCYTEC-
-0.6
-0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 OAJ). Este trabajo ha hecho uso de la base de datos
B-V SIMBAD, operada por el CDS (Centre de Données
astronomiques de Strasbourg), Estrasburgo, Francia.
Figura 6. Ajuste lineal (B – V) – (b – v)0 vs. (B – V) para la
obtención de µ y ζbv. Los valores calculados son µ = 1.2169 ±
0.0440 y ζbv = -0.5510 ± 0.0320.
REFERENCIAS
-1.9
1. Johnson, H.L. and Morgan, W.W.,
“Fundamental stellar photometry for standards of
-2.0
spectral type on the revised system of the Yerkes
spectral atlas” “Astrophys. J.”, 117, 313-352
(U - B) - (u - b)0
-2.1 (1953)
2. Cannon, A.J., and Pickering, E. C., “The Henry
-2.2 Draper Catalogue” “Ann. Astron. Obs. Harvard.
College”, 91-99 (1918-1924)
-2.3 3. Cannon, A.J., “The Henry Draper Extension”
“Ann. Astron. Obs. Harvard College”, 100
-2.4 (1925-1936)
-0.5 0.0 0.5 1.0 1.5 2.0 4. Tody, D., “The IRAF Data Reduction and
U-B Analysis System”, Proc. SPIE Instrumentation in
Astronomy VI, ed. D.L. Crawford, 627, 733
Figura 7. Ajuste lineal (U – B) – (u – b)0 vs. (U – B) para la
(1986)
obtención de ψ y ζub. Los valores calculados son ψ = 0.8645 ±
0.0137 y ζub = -1.7487 ± 0.0292.
9