13 março 2013

An Exercise of Astrometry with C/2011 L4 (Panstarrs)

Fig. 1 Intensity plot of C/2011 L4 (March 3 2013). Numbers indicate intensity levels in greyscale (0  - black ; 255 - white)
by Ademir Xavier

On last March 3rd I took some photos of C/2011 L4 (Panstars). Here I present an excercise of astrometry using one of those images to extract some meaningful data from this observation. What I do here can be applied to any other comet image or, in fact, celestial object in the sky. My aim is to show a simple example of practical applicaton of mathematics and geometry in the determination of sizes of celestial objects.

My interest is:
  • To estimate the apparent (in minutes of arc) and real dimension of the tail (in km or mi) of C/2011 L4;
  • To estimate the apparente (in minutes of arc) and real dimension of the coma (in km or mi);
A digital estimate of the comet brightness will not be attempted however, because this would involve a complex process of image calibration.

For that aim, we need:
  1. Comet image with date;
  2. Scale calibration;
  3. Software to extract brightness levels (the so called 'isointensity' curves);
  4. Distance of Earth to the comet at the date;
  5. Comet position angle in relation to the sun at the date;
I should also mention the need of a good sky simulator software (for all practical purposes I will use Stellarium, but other software could work as well). These are the 'inputs' of the work and the 'outputs' are described above, the physical dimension of the comet (at least an approximate value for this dimension).

Below I coment step-by-step all procedures that I used to find the final estimates. This exercise demonstrates a practical aspect of astronomical observation, something that is fully in agreement with the objectives of this blog.

1) Comet image with date

I use the image publish on last March 4 2013 reproduced below. The image was acquired on March 3rd at 22:15 UT and it is an important input for the determination of additional parameters as we will see. 

Fig. 2 Image used for the exercise. (click on the image to enlarge)
2) Scale calibration

Fig.2 shows not only Panstarrs but also a star named HIP 117488 of mag. 7.0. This star was easily identified with Stellarium using the date and time as input for the skymap display. The apparente distance between the star image centroid and the 'comet nucleus' is about 40' (forty minutes of arc) - Fig. 3.
Fig. 3 Estimate of apparent distance between HIP 117488 and comet C/2011 L4 at the date as given by Stellarium. 
Then we calculate the distance in pixels on the image between HIP 117488 and the comet and find 355.1. Therefore the scale factor (Sf) will be

Sf = 40'/355.1 ~ 0.113'/pixel.

This is nearly 6" per pixel and corresponds to the final resolution of the image. Note that this value is the overall resolution of both the combined optics and camera setup.

The resolution above is a practical scale for the determination of the comet dimension. Using Stellarium, we find that at the observation date, the comet-Earth (observer) distance (D) was:

D=1.09843629 AU.

Since 1 AU = 149,597,870,700 meters (92,955,807.273 mi) and Sf in radians is

Sf(rad) = 0.000032763 rad/pixel

then

Sf(km)= Sf(rad)*D(km) = 5383.74 km/pixel (=3345.3 mi/km)

Therefore, each pixel in the image at the comet position corresponds to about 5400 km. The smallest pixel in the image, in particular the one corresponding to the "comet nucleus", is a square of ~5400 x 5400 km, much larger than the expected physical size of that nucleus.

3) Detail analysis of the cometary image

The image is in fact a matrix of intensities on an arbitrary scale (in a grey scale 8 bit image, the intensity goes from 0 to 255). If we extract a small portion of the original image (after converting it to grey scale), say, a square of 35X35 centered at HIP 117488, we get Fig. 4.

Fig. 4 A small sample of the original image showing HIP 117488.
Fig. 5 Intensity surface of Fig. 3 of HIP 117488.

A 3D intensity plot of this image is shown in Fig. 5. This plot was made with Mathcad. What about the cometary image? I resampled the original image to 310 X 460 after converting it to greyscale and the resulting 3D intensity plot is:
Fig. 6 3D intensity plot of the original image showing the comet and HIP 117488. 
The intensity of the coma region is close to 255 (the maximum) and should be compared to the 'background' level between 10 and 15. Another possibility is to ask Mathcad to plot the intensity curves (Fig. 1). The tail length in pixels as estimated using Fig. 1 is about 90.5. Therefore the (aparent) tail length is
  • Apparent Tail(km) = 90.5 x Sf(km) ~ 490 000 km (= ~ 300 000 mi)
  • Apparent Tail(min of arc) = 10.2'

The apparent tail length was small (nearly 1/3 of the moon diameter).
Fig. 7 Intensity plot of the coma.
Fig. 7 is a zoomed version of Fig. 1. If we take the dimension of the coma region as nearly equal to que inner square (between 16 and 24 on the X-axis and 16 and 24 on the Y-axis), a good estimate of the coma condensation will span an area of 8 pixels X 8 pixels or almost 1'x 1'. This corresponds to sphere of ~50000km of diameter. However, if the outer intensity curves are regarded, the coma will have twice that size. Thus  the coma of Panstarrs was estimated to have an apparent size of 2'x2' or 90 000km of extension on last March 3rd.

A minor detail

It could be argued (with reason) that the estimated tail size must be corrected for the geometrical situation shown in Fig.8. The comet tail always points towards the sun (along the Sun-comet line), while we are observing the projection of this line on the perpendicular to the Earth-comet line.
Fig. 8
Fig. 8 depicts the geometrical situation: the apparent tail length is a function of the real tail length and the position angle (alpha) between the Sun-comet line and the perpendicular of the Earth-comet line. Therefore, the real tail length will be

Real Tail(km) = Apparent Tail(km)/cos(alpha).

It is not difficult to see that, for the geometry of Fig. 8:


Now, again using Stellarium we have:

Des = 0.99156172 AU;
phi = 18 deg 40' (elongation angle);

as the Earth-Sun distance and elongation angle for comet Panstarrs at the date, respectively. Therefore, using the above equation we find:

alpha = 26 deg 33',

so that cos (alpha) = 0.8944249989.

The real tail length in km (mi) will be

Real Tail(km) =  90.5 x Sf(km)/cos(alpha) = 545 000 km (~340 000 mi).

Compare this with the Earth-Moon distance (384 400km). There is no accurate definition of a comet tail (that depends on the density of particles such as dust, gas etc). What we can say here, however, is that, given the "definition" of tail as determined by the smallest intensity level on Fig. 1 (30.5), the real tail  extended itself for half a million kilometers in space on the date.

References

03 março 2013

Observação em 3 de Março 2013 do cometa C/2011 l4 (Panstarrs)

ATUALIZAÇÃO (12/3/2013): Haverá uma janela de 10 a 20 minutos apenas para observar esse cometa nas noites que seguem o dia 12/3. Ele não é um objeto visível facilmente à vista desarmada, exceto no caso de observadores muito bem treinados, que morem em lugares elevados e com horizonte ocidental livre de qualquer tipo de nuvem ou neblina. Dependendo da quantidade dessa neblina, nem mesmo binóculos poderão mostrar o cometa. A melhor opção é um telescópio com baixo aumento (grande campo de visão). O C/2011 L4 ou Panstarrs não será mais visível para observadores no hemisfério sul a partir do dia 16/3.


Imagem do cometa C/2011 l4 (Panstarrs) por entre fios de alta tensão em Campinas, SP às 7:15 do horário local. Esta foto foi tirada no alto da Cidade Universitária (UNICAMP).

O cometa pode ser visto à vista desarmada por entre as brumas do crepúsculo. A imagem abaixo é um 'zoom' de outra tomada onde aparece a estrela HIP 117488 de magnitude 7,0 da constelação de Aquário.



01 março 2013

Cometas em 2014: C/2013 A1 (colisão com marte?)

Imagem de Janeiro de 2013 do cometa Siding Spring tirada no Observatório do Vaticano (telescópio de 1,83 m). Este seria um cometa normal, se não fosse a chance de colisão com marte em Outubro de 2014...

O cometa C/2013 A1, chamado de 'Siding Spring', seria um cometa "normal" para 2014 que atingirá magnitude abaixo de 9.0 (visível por meio de binóculos em regiões de céu muito escuro). O problema é que, em estimativas mais recentes de sua órbita, ele tem a chance de colidir com marte!

Descoberto em 3 de Janeiro de 2013 por Robert McNaught no Observatório de Siding Spring com um telescópio de 50 cm, uma imagem recente desse cometa pode ser vista acima (imagem obtida no Observatório do Vaticano). 

Siding Spring: em rota de colisão com marte?
Se ocorrer colisão, será no dia 19 de Outubro de 2014. As estimativas de hoje (março de 2013) mostram que ele passará a menos de  0,0007 Unidades Astronômicas de marte, ou aproximadamente 100 mil quilômetros, o que é bem pouco para os padrões de distância planetários. Como há incertezas, essa estimativa pode representar uma colisão de fato. 

Estimativas ainda mais recentes (feitas por L. Elenin em 27/2/2013) mostram que ele poderá chegar a 41 mil quilômetros de marte. No caso das previsões mais catastrofistas (que estima um diâmetro de 50 km para o núcleo do cometa), a colisão seria realmente muito grande, sem nenhuma 'destruição' de marte, apenas para o cometa.

Hoje, este cometa brilha com mag. 18.1 (inacessível a telescópios amadores) e, em outubro de 2014, sua 'colisão' com marte, desde o Brasil, não poderá ser visível, pois se dará durante o dia (a hora de máxima aproximação prevista será por volta das 16:50 do tempo legal).  Nos dias anteriores, a magnitude estimada será pouco pouco abaixo de 9,0. Nessa data, o cometa e marte se encontrarão na constelação de Ofíuco.

De qualquer forma, na data de 19 de Outubro de 204, como visto desde marte, esse cometa poderá atingir mag. -8,0, o que o colocaria como um dos mais brilhantes objetos do céu. 

Posição no céu do cometa Siding Spring na noite de 18 de Outubro de 2014  (21:15 TL).  A 'colisão' não será visível  no  Brasil. A distância entre marte e o cometa nesta data é de aproximadamente 43 minutos de arco que serão consumidos nas horas seguintes.
Referências