Некоторые вопросы поиска и формирования планетных систем

1. Некоторые вопросы поиска и формирования планетных систем С. И. Ипатов AESS, Катар; ИКИ, Москва

2. SIMULATOR FOR MICROLENS PLANET SURVEYSМодель яркости звездного неба и алгоритм оптимального выбора целей для наблюдений при поиске экзопланет методом микролинзирования. С. И. Ипатов • На основе данных наблюдений, полученных в 2011 с помощью четырех телескопов, использующихся для поиска планет методом микролинзирования, построена модель яркости звездного неба. Эта модель является частью алгоритма оптимального выбора целей для наблюдений при поиске экзопланет методом микролинзирования. Алгоритм позволяет определять какие цели доступны для наблюдений в различное время и какие цели следует выбирать для того, чтобы максимизировать вероятность обнаружения экзопланет. An older version of a presentation on this item can be found on http://star-www.st-and.ac.uk/~si8/skybrightnessiau.ppt

3. ---- Abstract We summarize the status of a computer simulator for microlens planet surveys. The simulator generates synthetic light curves of microlensing events observed with specified networks of telescopes over specified periods of time. The main purpose is to assess the impact on planet detection capabilities of different observing strategies, and different telescope resources, and to quantify the planet detection efficiency of our actual observing network, so that we can use the observations to constrain planet abundance distributions. At this stage we have developed models for sky brightness and seeing, calibrated by fitting to data from the OGLE survey and RoboNet observations in 2011. Time intervals during which events are observable are identified by accounting for positions of the Sun, the Moon and other restrictions on telescope pointing. Simulated observations are then generated for an algorithm that adjusts target priorities in real time with the aim of maximizing planet detection zone area summed over all the available events. 3

4. Studying planets bymicrolensing • There is no single technique known that can yield thecomplete planet demographics. Microlensing isunique in its sensitivity to wider-orbit (i.e. cool) planetary-mass bodies • 1) down to the mass of the • Moon, already with ground-based observations, • 2) in orbit around distant stars, even in neighbouringgalaxies, 3) in orbit around faint or dark stars or remnants, such as brown dwarfs, white dwarfs, neutronstars, or black holes, • 4) not bound at all. It thereby allows to explore a yetenigmatic region in massseparationspace, where even a small amount of data has the potential to make a large impact. • 4

5. If the source crosses a caustic, the deviations from a standard event can be large even for low mass planets. These deviations allow us to infer the existence and determine the mass and separation of the planet around the lens. Deviations typically last a few hours or a few days. Because the signal is strongest when the event itself is strongest, high-magnification events are the most promising candidates for detailed study.2009 Gliese 581 e is an exoplanet with 1.9 Earth masses.

6. Styding planets by microlensing The 15% blip lasting about 24 hrs thatrevealed5-Earth-mass planet OGLE-2005-BLG-390impressively demonstrated the sensitivity ofongoing microlensing efforts to Super-Earths.Had an Earth-mass planet been in the samespot, it would have been detectable from a 3%signal lasting 12 hrs. The detection of lessmassive planets requires photometry at the fewper cent level on Galactic bulge main-sequencestars, which, given the crowding levels,becomes possible with images of angularresolution below about 0.4″.6

7. Detection zone • Based on the approach presented in [3], at each time step for different events we calculate the detection zone area and the probability of detection of an exoplanet. The event with a maximum probability at a time step is chosen for observations. • We define the ‘detection zone’ (зона обнаружения) as the region on the lens plane (x,y) where the light curve (кривая блеска) anomaly δ(t,x,y,q) is large enough to be detected by the observations (q is the ratio of the planet to that of the star). [3] Horne K., Snodgrass C., Tsapras Y., MNRAS, 2009, v. 396, 2087-2102 Detection zones on the lens plane indicate the regions where a planet with mass ratio q=m/M=10-3 is detected with Δχ2>25. The light curve A(t) has maximum magnification Ao=5, and the accuracy of the measurements is σ=(5/A1/2) per cent. 7

8. Maximising planet detection zone area The photometric S/N (signal to noise) ratio and hence the area w of an isolated planet detection zone scales as the square root of the exposure time : S/N = (Δt /τ) 1/2, w = g Δt 1/2 . Here τ is the exposure time required to reach S/N=1. The 'goodness' gi of an available target depends on the target's brightness and magnification, the telescope and detector characteristics, and observing conditions (airmass, sky brightness, seeing). The simulator evaluates 'goodness' of available targets in real time, and observes the one offering the greatest increase in w with exposure time. Moves to a new target occur when the increase in w for the new target is better than the current target, accounting for the slew time required to move to the new target. As the CCD camera takes a finite time tread to read out, and the telescope takes a finite time tslew to slew from one target and settle into position on the next, the on-target exposure time accumulated during an observation time t is Δt=t- tslew –n tread (tslew ~ 1-3 min, tread ~10-20 s). At a time step Δt the detection area of i-th event increases by gi[(Δt+tdone) 1/2 - tdone1/2 ], where tdone is an exposure time already has been done. For a new target, the area is gi (Δt-tslew) 1/2 . See [3] for details. [3] Horne K., Snodgrass C., Tsapras Y., MNRAS, 2009, v. 396, 2087-2102. 8

9. Example of comparison of gi[(Δt+tdone) 1/2 - tdone1/2 ] for a choice of the best event to be observed (OGLE observations of 10 events: 110251 – 110260) at a time step of 200 s and tslew=100 s. All but one of the targets require slew time before the exposure can begin. 9

10. Target observability • The observability of a target is limited by its own position on the sky, as well as that of the Sun and the Moon, and telescopes moreover have pointing restrictions. Taking the LT (from http://telescope.livjm.ac.uk/) as example, we particularly require not to make observations at: • Air mass (≈1/cos[zenith angle of a target; зенитное расстояние]) > 3 orCos (zenith of the Sun) < sin (-8.8o) or • Altitude of a target (высота цели): alt<altmin=25o or alt>altmax=87o orHour angle: ha<hamin or ha>hamax. For LT there are no limits on ha: hamin=-12 h and hamax=12 h. • Telescopes considered: • 1.3m OGLE - The Optical Gravitational Lensing Experiment - Las Campanas, Chile. • 2m FTS - Faulkes Telescope South - Siding Springs, Australia. • 2m FTN - Faulkes Telescope North - Haleakela, Hawaii. • 2m LT - Liverpool Telescope - La Palma, Canary Islands. • 10

12. Observations analyzed for construction of sky model: • For studies of sky brightness for FTS, FTN, and LT, we considered those events observed in 2011 for which .dat files are greater than 1 kbt: FTS - 39 events; FTN - 19 events, LT – 20 events. For OGLE we considered 20 events (110251-110270). The observations were infrared. The used sky model was based mainly on K. Krisciunas & B. Schaefer, 1991, PASP, v. 103, 1033-1039. • Calculations of vskyzen (яркости звездного неба в зените) and the coefficients (k1 and ko)presented in the tables and on the plots were based on χ2 optimization of the straight line fit (y=k1·x+ko, χ2=∑[(yi-k1·xi-ko)/σi]2, σi2 is variance - дисперсия). The value of vskyzen (sky brightness at zenith) was chosen in such a way that the sum of squares of differences between observational and model sky brightness magnitudes (разности между наблюдаемыми и модельными звездными величинами яркости неба) were minimum in the case when the Moon is below the horizon. 12

13. Dependences of seeing on air mass and values of sky brightness at zenith obtained based on analysis of observations Values of sky brightness at zenith (I magnitude per square arcsec) for an extinction coefficient extmag=0.05 (for extmag equal to 0 and 0.1, values of vskyzen differed by less than 0.3 %): 13 Seeing (FWHM in arcsec) vs airmass (χ2optimization): seeing=ko+k1×(airmass-1). Seeing is the blurring effects of air turbulence in the atmosphere (угловой диаметр кружка, в виде которого изображение звезды предстает в телескопе; эффект расплывчатости изображения из-за турбулентности в атмосфере). σ– среднеквадратичное отклонение

14. Seeing vs. air mass. Seeing (in arcsec) vs. air mass. FTS observations of 39 events. A thick straight line is based on χ2 optimization (y=a+b(x-1), a=1.334, b=0.519). Thinner straight lines differ from this line by +/- Ϭ (Ϭ=0.367). Non-straight lines show mean and median values (the line for the mean value is thicker). 14

15. ) Sky brightness (mag) (звездная величина яркости неба) vs. air mass. Different points are for OGLE observations of 20 different events (110251 – 110270) for the Moon below the horizon and solar elevation < -18o. The lines are for the χ2 optimisation with different ko (different values for different events) and the same k1. The most solid line is for the model for which ko is the same for all events. 15

16. The difference in sky brightness near different events.Coefficients for sky brightness=bo+b1×(airmass-1). One value of b1; values of bo are different for different events. The range of boi for Moon below the horizon and solar elevation < -18o: • min max max-min telescope • 19.49 20.41 0.92 FTS • 19.05 20.15 1.10 FTN • 19.89 20.60 0.71 LT • 19.64 20.35 0.70 OGLE The values of max-min were about 0.7-1.1 mag. The difference max-min characterizes the difference in sky brightness due to surrounding stars near different events. The range ofboi for FTS for different positions of the Moon and the Sun: min max max-min data considered • 18.2 20.2 2.0 - all observations • 19.3 20.4 1.1 - moon below the horizon • 19.5 20.4 0.9 - moon below and solar elevation < -18o The upper limit (less bright observations) in the above table does not vary much; the difference in the lower limit (more bright sky) is greater and is up to 1.5 mag. When the Moon is below the horizon and solar elevation <-18o, then σ is smaller than that for all observations by a factor of 3 (0.11 instead of 0.36). The data of the table show the influence of positions of the Moon and the Sun on a typical sky brightness near an event. 16

17. Residuals of sky brightness Most of residuals (разности) (observations minus χ2 optimization which is different for different events) of sky brightness are in a small range (-0.4 to 0.4 mag) even for all Moon and Sun positions; for the Moon below the horizon there are many of values in the range [-0.2, 0.2]; greater values of residuals are for a small number of observations. The range of sky brightness residuals for the Moon below the horizon and solar elevation <-18o (for example, [-0.42, 0.87] for FTS) is smaller by a factor of several than for all positions of the Moon and Sun. For considered observations, there was no bright sky (i.e. the lower limit of residuals was >-1 mag) only in the case when both the Moon was below the horizon and the solar altitude was <-18o. 17

18. Sky brightness residuals (mag) vs. air mass for the model with different ko for FTS observations of 39 events. Left plot is for all positions of the Moon and the Sun. Right plot is for the Moon below the horizon and solar elevation < - 18o. 18

19. Sky brightness residuals vs. solar elevation Analysis of the plots shows that the influence of solar elevation (высоты солнца) on sky brightness began to play a role at solar elevation > -14o, and was considerable at solar elevation >-7o. For example, if we consider only observations with the Moon below the horizon, then for FTS: sbr >-0.4 mag at se<-14o, sbr>-1 at se<-8o, sbr can be up to -3 mag at se in the range (-8o,-7o), where se is the value of solar elevation, and sbr is the value of sky brightness residual (in mag). Sky brightness residuals (in mag.) vs. solar elevation for FTS observations of 39 events. More dense signs are for the Moon below the horizon. 19

20. Time intervals when it is better observe events Time intervals for events selected for observations with OGLE. Considered events: 110250-110400 (numbered from 1 to 151). 0 event corresponds to ‘no observations’. Initial Julian time=2455728.5 (June 15, 2011). 20

21. Light curves Light curves for events selected for observations with OGLE. Considered events: 110250-110400. Initial Julian time=2455728.5 (June 15, 2011). 21

22. CONCLUSIONS At this stage we have developed models for sky brightness and seeing, calibrated by fitting to data from the OGLE survey and RoboNet observations in 2011. Time intervals during which events are observable are identified by accounting for positions of the Sun, the Moon and other restrictions on telescope pointing. Simulated observations are then generated for an algorithm that adjusts target priorities in real time with the aim of maximising planet detection zone area summed over all the available events. REFERENCES • [1] M. Dominik, 2010, Gen. Relat. Gravit. 42, 2075 • [2] Y. Tsapras, R. Street, K. Horne, et al., 2009, AN 330, 4 • [3] K. Horne, C. Snodgrass, Y. Tsapras, 2009, MNRAS 396, 2087 • [4]K.Krisciunas&B.Schaefer,1991,PASP, 103,1033 • 22