REDUÇÃO DE COMETAS
Objetivo: calcular os parâmetros fotométricos de um cometa usando dados observacionais de associações locais. O trabalho também consiste em comparar os dados obtidos com outras fontes de modo que os valores calculados se ajustem melhor com as observações.

1ª ETAPA:

A Home Page de Cometas/REA, possui registros observacionais desde 1988. No link "Arquivo de Observações" é possível obter dados para cada cometa observado. No nosso exercício iremos usar os dados do Cometa C/1999J3 (LINEAR).

Para calcular os parâmetros fotométricos usaremos principalmente as datas TU, as magnitudes estimadas e os instrumentos usados.

2ª ETAPA:

Devemos corrigir as magnitudes para o instrumento padrão. Normalmente durante o acompanhamento de um cometa, são usados diversos instrumentos: telescópios, binóculos, etc... Mas precisamos situar as magnitudes para um
instrumento padrão de abertura 6.8cm . Para isso usaremos as fórmulas abaixo:

refletores: m' = m - 0.019(A - 6.8)     &nbbsp;                 (1)

refratores: m' = m - 0.066(A - 6.8)     &nbbsp;                 (2)

Exemplo:

1999 Dez 6.32 S 7.7 6' 3 11x80B J.Aguiar, Campinas, Brasil

A magnitude estimada foi de 7.7 usando um binóculo 11x80.

Temos m = 7.7 e A = 8 cm

usa-se a formula (2) para refratores (o binóculo nada mais é do que 2 lunetas):

m' = m - 0.066(A - 6.8)
m' = 7.7 - 0.066 (8 - 6.8)
m' = 7.7 - 0.066 (1.2)
m' = 7.7 - 0.0792
m' = 7.6208 ou 7.62

7.62 seria a magnitude estimada por um instrumento padrão de 6.8 cm.

Com isso a tabela 1 mostra todas as observações do Cometa C/1999J3 com as magnitudes corrigidas.

 
Data
mag
mag. CORR.
OBS
INST
06/out/1999 7:45
7.7
7.62
JGA
11x80B
07/out/1999 7:45
7.6
7.52
JGA
11x80B
09/out/1999 8:00
7.6
7.52
JGA
11x80B
10/out/1999 7:45
7.7
7.62
JGA
11x80B
11/out/1999 8:00
7.8
7.72
JGA
11x80B
12/out/1999 6:00
7.7
7.62
PMR
11x80B
12/out/1999 8:00
7.9
7.82
JGA
11x80B
13/out/1999 8:00
8.0
7.92
JGA
11x80B
21/out/1999 5:20
8.0
8.12
AMO
7x50B
31/out/1999 2:10
8.3
8.16
AMO
14L
02/nov/1999 1:30
8.6
8.29
JGA
23L
06/nov/1999 2:00
8.7
8.56
AMO
14L
07/nov/1999 2:30
8.7
8.56
AMO
14L
07/nov/1999 23:20
9.0
8.86
AMO
14L
13/nov/1999 1:45
9.0
8.86
AMO
14L
 
Tabela 1

3ª ETAPA:

Devemos determinar os valores de Raio Vetor ( r ) e Delta ( D ) para cada Data TU

O Delta é a distância entre o Cometa e a Terra para a referida data. O Delta é definido em Unidades Astronômicas (UA). O Raio Vetor (r) é a distância entre o Cometa e o Sol para a referida data, também em UA.

Alimentaremos os programas SkyCharts, SkyMap, Guide ou qualquer outro software calculador de efemérides com os elementos orbitais recentes do Cometa a ser estudado, neste caso o C/1999J3 - MPC 36212.

Para cada data e hora TU da observação, haverá um valor de D e r calculados pelos programas.

Exemplo:

No dia 6 de outubro de 1999, o cometa C/1999J3 foi observado com magnitude 7.7
A fração de dia em TU é 0,32.

0,32 de dia equivale a 7:45 TU.

Usando o soft SkyMap 8, para o dia 6/10/1999 7:45 TU temos:

Distance: 0.6177 AU ( 92.4 million km)

Radius vector: 1.0166 AU (152.1 million km)

Agora vamos obter alguns cálculos:

1) Usando uma planilha eletrônica, colocaremos nas colunas subsequentes, os valores de raio (r) e delta (D).

2) Calcularemos individualmente os logaritmos de (r) e D .

3) O eixo das abcissas "X" do gráfico para redução possui os valores da fórmula abaixo:

X = 2,5 . log(r)                               [2 vezes e meia o valor de logaritmo do raio]                                                    (3)

4) O eixo das ordenadas "Y" do gráfico possui os valores da fórmula:

Y = m - 5.log(D)                            [ m menos 5 vezes o logaritmo de D]                                                              (4)

A Tabela 2 mostra os resultados obtidos.

 
- - - - - -
Abcissas X
ordenadas Y
Data
mag
r
D
log ( r )
log (D)
2,5 log ( r )
m - 5log (D)
6/10/1999 7:45
7.62
1.0166
0.6177
0.0072
-0.2092
0.0179
8.6656
07/10/1999 7:45
7.52
1.0215
0.6012
0.0092
-0.2210
0.0231
8.6244
09/10/1999 8:00
7.52
1.0321
0.5714
0.0137
-0.2431
0.0343
8.7348
10/10/1999 7:45
7.62
1.0377
0.5584
0.0161
-0.2531
0.0402
8.8848
11/10/1999 8:00
7.72
1.0435
0.5465
0.0185
-0.2624
0.0462
9.0315
12/10/1999 6:00
7.62
1.0492
0.537
0.0209
-0.2700
0.0521
8.9696
12/10/1999 8:00
7.82
1.0498
0.5362
0.0211
-0.2707
0.0528
9.1728
13/10/1999 8:00
7.92
1.0561
0.5274
0.0237
-0.2779
0.0593
9.3088
21/10/1999 5:20
8.12
1.1137
0.5186
0.0468
-0.2852
0.1169
9.5433
31/10/1999 2:10
8.16
1.2012
0.6483
0.0796
-0.1882
0.1990
9.1039
02/11/1999 1:30
8.29
1.2204
0.6869
0.0865
-0.1631
0.2163
9.1074
06/11/1999 2:00
8.56
1.2609
0.7732
0.1007
-0.1117
0.2517
9.1214
07/11/1999 2:30
8.56
1.2713
0.7962
0.1042
-0.0990
0.2606
9.0577
07/11/1999 23:20
8.86
1.2804
0.8162
0.1073
-0.0882
0.2684
9.3038
13/11/1999 1:45
8.86
1.3348
0.9383
0.1254
-0.0277
0.3135
9.0011
 Tabela 2

Depois criamos o gráfico com os valores das duas últimas colunas. A Figura 1 mostra os pontos e duas regressões aplicadas.
 

Figura 1

4ª ETAPA:

Agora vem a parte mais interessante que é a interpretação dos dados.

No caso do Cometa C/1999J3 temos alguns fatos curiosos decorrentes do cálculo de regressão linear e polinomial.

    1. A regressão linear foi aplicada nas 9 observações iniciais feitas por meio de binóculos. O MS-Excel calculou a função pertencente a esses pontos (valores arredondados):     y = 15.5x  +  8.3
    2. Uma regressão polinomial de 2ºgrau foi aplicada a todas as observações e obtivemos a seguinte função:  y = -36x2 + 12.4x + 8.5
    3. Os dois pontos vazados não foram incluídos no cálculo da regressão.
    4. A construção deste gráfico permite que obtenhamos os valores de Ho e n da função linear ax + b em que o parâmetro n corresponde ao coeficiente angular a ; e Ho é o coeficiente linear b.
Deste modo teríamos Ho = 8.3 ou Ho = 8.5 ;

O valor de n poderia ser 15.5 ou 12.4 .

Mas qual valor utilizar?

Bem, devemos lembrar que só utilizamos 15 observações realizadas num intervalo de tempo de 35 dias. Isso pode não ser suficiente na obtenção de valores mais precisos dos parâmetros fotométricos, apesar de usarmos o procedimento comum.

Então é interessante comparar os resultados com as fontes internacionais.

Neste caso, as efemérides do MPC/IAU forneciam os valores de Ho = 12 e n = 4 (MPC 36212). Ora , vemos que o valor de Ho está muito alto em relação às observações. Por outro lado, n está no valor usual para novos cometas. Normalmente o valor de n situa-se no intervalo 2 £ n £ 4.

Já a análise de 197 observações feitas pela BAA Comet Section revelou Ho = 9 , valor bem mais próximo do que foi obtido com as observações locais. Já o valor de n ficou em 4.7 . Mesmo assim, o coordenador desta secção indicou que tais valores ainda seriam "incertos".

Dessa forma, não exigiríamos muito das observações locais se usarmos o valor Ho = 8.5 e n = 4 , sendo que o valor de n foi tomado como maior valor do intervalo comum aos cometas novos.

Na Figura 2 temos diversas curvas de luz. Aquela formada por pontos corresponde aos valores visuais estimados. A linha tracejada corresponde a curva baseada nos valores do MPC e foi obtida pela fórmula

m = 12 + 5logD + 10log r                                                        (5)

A linha cinza foi obtida com os valores da regressão linear e segue a fórmula abaixo:

m = 8.3 + 5logD + 38.8log r                                                    (6)

Já a linha contínua em preto foi obtida da fórmula

m = 8.5 + 5logD + 10log r                                                       (7)

CONCLUSÃO

Vemos que a fórmula (7) se aproxima melhor das observações, de modo que daremos preferência a tais valores. Por essa razão nem sempre os resultados de uma regressão devem ser seguidos a risca. Caberá ao calculista interpretar seus resultados e aplicá-los às observações da melhor maneira possível.

Figura 2

 

BIBLIOGRAFIA:

-KEITCH,G.S. - Visual Comet Photometry, JBAA, vol 93-5, pp 200-3 (1983)
-BRASIL Jr,C.L. e NAPOLEÃO,T.A. - Determinação dos Parâmetros Fotométricos de Cometas Recentes, REA, vol 5, pp 41-5 (1992).
-The Comet's Tale (BAA Comet Section), vol 7-1, p 15 (2000)
-AMORIM,A. - Construindo Curvas de Luz com o MS-Excel (2002)
-Secção de Cometas/REA - Observações do Cometa C/1999J3 (2002) em
http://www.geocities.com/costeira1/cometa/observ99j3.htm

 
©2002 Costeira1 - Alexandre Amorim
(volta)
1