IOAA 2012, Анализ на данни, английски — Задача 1. Question 1

6th International Olympiad on Astronomy and Astrophysics, 04 to 14 August, 2012 - Rio de Janeiro - Brazil — Data analysis · 4 август 2012 г. · 13,5 т.

Условие

A few facts about photometry of asteroids:

  • Asteroids are small, irregularly shaped objects of our Solar System that orbit the Sun in approximately elliptical orbits.

  • Their brightness observed at a given instant from Earth depend on the surface area illuminated by the Sun and the part of the asteroid which is visible to the observer. Both vary as the asteroid moves.

  • The way the sunlight is reflected by the surface of the asteroid depends on its texture and on the angle between the Sun, the asteroid and the observer (phase angle), which varies as the Earth and the asteroid move along their orbits. In particular, asteroids with surfaces covered by fine dust (called regoliths) exhibit a sharp increase in brightness at phase angles φ\varphi close to zero (i.e., when they are close to the opposition).

  • Since the observed flux of any source decreases with the square of distance, the observed magnitude of an asteroid also depends on its distance from the Sun and from the observer at the time of the observation. Their apparent magnitude mm outside the atmosphere is then

    M(t)=Mr(t)+5log⁡(RD)M(t) = \mathrm{M_r}(t)+5\log(RD)

    where mrm_r is usually called reduced magnitude (meaning the magnitude the asteroid would have if its distances from the Sun and the Earth were reduced both to 1 AUAU) and depends only on the visible illuminated surface area and on phase angle effects. RR and DD are, respectivelly, the heliocentric and geocentric distances.

Consider now the following scenario. Light curves of a given asteroid were obtained at three different nights at different points of its orbit, and at each time a photometric standard was observed in the same frame of the asteroid. Table 1 shows the geometric configuration of the asteroid at each night (phase angle φ\varphi, in degrees, RR and DD in AUAUs), and the calibrated magnitude of the photometric standard star that was observed along with the asteroid. Consider the calibrated magnitude as the final apparent magnitude after all the effects are taken into account.

Tables 2, 3 and 4 contain, for each night, the time interval of each observation with respect to the first one (in hours) the air mass, the instrumental magnitude of the asteroid and the instrumental magnitude of the star.

Air mass is the dimensionless thickness of the atmosphere along the line of sight, and is equal to 1 when looking to zenith.

1. Plot the star magnitude versus air mass for each set [1,5 т.]

2. Calculate the extinction coefficient for each night (see Appendix A at the end of the text) [3 т.]

3. Were the observations affected by clouds in one night? Express your answer as: a) Night A, b) Night B, c) Night C, d) None of the nights. [0,5 т.]

4. Plot the calibrated magnitude versus time for each set of observations of the asteroid (see Appendix B). [3 т.]

5. Determine the rotation period for each night. Consider that the light curve for this asteroid has two minima and two maxima, and that the semi-period is the average of the intervals between the two maxima and the two minima. [1,5 т.]

6. Determine the amplitude (difference from maximum to minimum) of the light curve for each night [1,5 т.]

7. Plot the calibrated reduced magnitude Mr_{\mathrm{r}} versus phase angle φ\varphi (use the mean value of each light curve) [1,5 т.]

8. Calculate the angular coefficient of the phase curve (the plot of the calibrated reduced magnitude versus the phase angle) considering only the points away from the opposition (See 3rd^{\mathrm{rd}} bullet of facts about photometry of asteroids above) [0,5 т.]

9. Is there any reason to assume a surface covered by fine dust (regolith)? Answer YES/NO. [0,5 т.]

Table 1

NightDRφMstar_{\mathrm{star}}
A0.361.350.08.2
B1.152.138.68.0
C2.701.8915.68.1

Table 2: Night A

ΔtAir massmast_{\mathrm{ast}}mstar_{\mathrm{star}}
0.001.287.448.67
0.441.187.388.62
0.891.117.348.59
1.331.067.288.58
1.771.027.328.58
2.211.007.338.56
2.661.007.338.56
3.101.017.308.56
3.541.037.278.58
3.991.077.278.58
4.431.137.318.61
4.871.217.378.63
5.311.327.428.67
5.761.487.498.73
6.201.717.598.81
6.642.067.698.92
7.092.627.879.14
7.533.678.219.49

Table 3: Night B

ΔtAir massmast_{\mathrm{ast}}mstar_{\mathrm{star}}
0.001.2813.248.38
0.441.1813.218.36
0.891.1113.138.34
1.331.0613.118.33
1.771.0213.118.32
2.211.0013.158.32
2.661.0013.178.32
3.101.0113.178.32
3.541.0313.138.33
3.991.0713.158.34
4.431.1313.148.34
4.871.2113.148.37
5.311.3213.218.38
5.761.4813.308.43
6.201.7113.348.47
6.642.0613.398.54
7.092.6213.448.65
7.533.6713.678.87

Table 4: Night C

ΔtAir massmast_{\mathrm{ast}}mstar_{\mathrm{star}}
0.001.2811.648.58
0.441.1811.538.54
0.891.1111.568.60
1.331.0611.498.52
1.771.0211.588.48
2.211.0011.798.63
2.661.0011.678.53
3.101.0111.538.46
3.541.0311.478.48
3.991.0711.638.67
4.431.1311.518.51
4.871.2111.658.55
5.311.3211.778.61
5.761.4811.888.75
6.201.7111.868.78
6.642.0612.039.03
7.092.6212.149.19
7.533.6712.639.65

Appendix A: extinction coefficients

The magnitude outside the atmosphere is given by:

M=m−A.X+B\mathrm{M = m - A.X + B}

where AA is the extinction coefficient and BB is the zero point coefficient for the night, XX is the airmass of the observation and mm is the instrumental magnitude (that is, the magnitude obtained directly from the image)

Appendix B: calibrated differential magnitude

If the standard star is the same frame as the object, the calibrated magnitude of the object can be obtained by:

Mast=mast−mstar+MstarM_{ast} = m_{ast} - m_{star} + M_{star}

Appendix C: least squares estimation of angular coefficients

Given straight lines described by

yi=α+βxi, i=1..ny_i = \alpha + \beta x_i,\ i = 1..n

the least squares estimator of β\beta is given by

β=∑inxiyi−1n∑inxi∑inyi∑inxi2−1n(∑inxi)2\beta = \frac{\sum_i^n x_i y_i - \frac{1}{n}\sum_i^n x_i \sum_i^n y_i}{\sum_i^n x_i^2 - \frac{1}{n}\left(\sum_i^n x_i\right)^2}

Отговори

Решение


Оригинал в Архива: DataAnalysisFINAL.doc · официални решения: DataAnalysis_Solutions.pdf