# arxiv: v1 [astro-ph.im] 6 Apr 2018

Save this PDF as:

Size: px
Start display at page:

## Transcription

4 A&A proofs: manuscript no. papier a gradient-descent optimization stage to the process, and reoptimize some of the best solutions found by the brute-force algorithm. To decide how many of the best grid values should be re-optimized, we ran the brute-force algorithm on a set of 25 images generated using our SPHERE/IRDIS simulator (see Sect. 2), in typical sky conditions (seeing of 0.8 ), in which no planet was introduced. This gave us a typical distribution of the noise values sampled by the grid, in which we found that the highest S/N value was 5.2, and that only 100 S/N values are higher than 4.5 and 1000 higher than 4.2. This means that if we re-optimize the best value found on the grid, a planet will be detected only if the corresponding S/N grid-maximum reaches a value higher than 5.4. However, if we re-optimize the 100 (resp. 1000) best values, then the planet will be detected if the corresponding grid value is greater than 5.1 (resp. 4.7). Compared to the total 10 8 points of the grid, re-optimizing 100 or 1000 values takes only little time, and allow to recover low S/N planets. In conclusion of this section, we built an algorithm which works in three steps: First, a brute-force algorithm is used to determine the value of the S/N function in each point of the grid. Then, the p = 100 highest values found are re-optimized by a gradient-descent algorithm or similar. We use the Broyden- Fletcher-Goldfarb-Shanno method, Fletcher Finally, we search for planets in all images corresponding to the re-optimized S/N values above Results 4.1. Blind test To test the K-Stacker algorithm, we performed a blind experiment. We used the SPHERE/IRDIS simulator presented in Section 2 to build 50 sets of 10 independent simulated observations, with a seeing value of 0.8. In each set, a planet was then injected on a random orbit, at a random S/N level. The orbit was selected by drawing the 6 parameters randomly using the distributions given in Table 2 allowing the algorithm to redraw a new set of parameters when the first one resulted in an orbit going outside of the AO corrected area in at least one image.the mass and the distance of the star are expected to be known with a good accuracy, as it is usually the case for the bright stars observed in coronagraphy. For each simulation, the planet was not injected at constant S/N, but rather at constant flux, to simulate what is expected from a constant exposure time serie of observations. For each serie of 10 observations, the noise level σ 1 at the position of the planet in the first image is computed, and the planet is injected at a flux value of F = snr 1 σ 1 on top of the background in each image, where snr 1 is randomly selected among four possible values : 0, 5/ 10, 7.5/ 10, or 10/ 10. Because the local value of the background behind the planet varies in each observation, as does the noise level if the planet is moving with respect to the central star, the exact S/N ratio differs from one image to the other, and the total expected S/N level of the simulation can only loosely be expected to be S/N 10 snr1 In all these simulations the star has a mass of 1M, is located at a distance of 10 pc, and has a magnitude of 8 in the R band (AO sensing band). The ten images of each set correspond to observations made at different times, selected arbitrarily to represent a plausible sequence of K-Stacker observations. The times are given in Table 3. The 50 sets of observations have been prepared using a dedicated computer program, which drew all the random parameters Parameter Range Distribution M star 1 M fixed value d star 10 pc fixed value a [0.2 A.U., 7.5 A.U.] uniform e [0, 0.5] uniform t 0 [-20 yr, 0 yr] uniform Ω [-180 deg, 180 deg] uniform i [0, 180 deg] uniform θ 0 [-180 deg, 180 deg] uniform 10 S/N {0, 5.0, 7.5, 10.0} uniform Table 2. Parameters used to inject the planet in the 50 simulation of our blind experiment. Observation number Date 1 January, 1, February, 6, March, 15, April, 20, Februray, 6, March, 15, February, 6, March, 15, Februray, 6, March, 14, 2020 Table 3. Dates used for each of the 10 observations constituting every set of the blind test and S/N values, and stored them in a file. All the simulations were then processed by the K-Stacker algorithm. For each simuation, K-Stacker produced a set of 100 best optimized recombinations, with their associated S/N values (before and after the re-optimization step) and orbital parameters. The observer has absolutely no idea of the S/N and orbital parameters of the planets injected in the data. The final images produced by K-Stacker can be checked by the observer, as usually done with high-contrast images. The observer can look at the final recombined images and, based on the total recombined S/N and the shape of the detection found by K-Stacker, assign a flag to each set of observation, validating or not the potential planet candidate. The file containing the parameters and S/N levels was only retrieved at the very end of the whole process in an effort to ensure that the observer never knew in advance what to expect from each simulation, and reduce overall bias Results obtained Among the 50 simulations prepared and analysed in our blind experiment, 15 were assigned by our algorithm to a group of S/N 0, 12 to S/N 5, 10 to S/N 7.5, and 13 to S/N 10. All the 15 simulations in which no planet was actually injected were correctly identified as containing no planet candidate by the observer. In the remaining 35 simulations, 25 planets were correctly identified as planet candidates by the observer, 9 planets were missed, and one false positive was found. In Figure 2, we show the distribution of the planets missed/found as a function of the total real S/N ratio of the simulation computed afterwards by combining the images using the true set of orbital parameters. Article number, page 4 of 12

6 A&A proofs: manuscript no. papier Fig. 3. One of the individual frame of simulation 48, which lead to the only fase positive of the serie. The planet that has not been detected is on the edge of the AO corrected area (indicated by the blue arrow). Simulation Grid search Gradient descent S/N 100th value Minimum S/N Maximum S/N Table 4. Values of the 100th best S/N ratio found by the grid search, and minimum/maximum values found by the reoptimization algorithm for each of the 15 simulations were no planet was injected. Fig. 4. Recombined K-Stacker image where the false planet candidate can be seen (indicated by the arrow). The observer noted the assymetry of the planet spot. the highest S/N value of the 100 best solutions found during the grid search as well as the minimum and maximum S/N values found by the gradient-descent algorithm. This table shows that if a simple S/N threshold were to be used to flag planet candidates, the number of false positives would be much higher. This result emphasizes the possibility of using a recombined spot shape criterion to disentangle false positives from true planets, and the importance of the observer s judgement for the reliability of this technique in its current implementation Can we get rid of the astronomer? In its current implemention, the K-Stacker algorithm produces 100 images for each time serie analyzed, and requires the intervention of a skillful observer to check these recombined images for the presence of a planet. This step, while also used in classical high-contrast imaging, is tedious, and may introduce uncon- trolled bias in the data reduction process. But for the time being, and as already discussed in Section 5.2, this step is absolutely necessary to avoid large number of false positives. In Figure 5 (resp. 6), we show what would have been the results of our blind test if we had used only a S/N threshold to flag planet candidates. For each simulation, the best solution found by K-Stacker is simply flagged as a planet candidate if the S/N is above 5 (resp. 7), without any intervention of an observer. As expected, the number of false positives (10, resp. 5 in total) is much higher, especially at low S/N. In our blind test, the observer was able to detect that these best solution found by K-Stacker were actually super-speckles, and he eliminated them (see also Fig. A.1). For most of them, he was also able to find the planet candidates. As discussed in Section 5.5, the shape of the recombined spot is important. To try to take that into account, we switched from our simple circular photometric aperture to a slightly more sophisticated gaussian weighted aperture. We found no real difference with our simple circular aperture, and this new method did not alleviate the need for the observer to carefully look at each image to assess the presence of the planet. A more complex algorithm, maybe one based on machine-learning trained to disentangle a planet PSF from recombined speckles, will be required in order to get rid of the observer s judgement. This would be a major improvement for K-Stacker, but left for future work. Other astrophysical informations (spectral, polarimetry) will also be used to remove the false alarms Orbital parameters determination K-Stacker can also provides an estimate of the orbital parameters, as a by-product of the optimization algorithm. The precision with which these parameters are estimated depends on many different factors, such as the actual orbit on which the planet is moving, the total time spanned by the observing sequence, etc. It is clear, for example, that if the planet does not move significantly during the sequence, one should not expect to have reliable information about the orbital parameters. In such cases, K-Stacker is able to recenter the images and detect planets not reachable with other methods, but several very different sets of orbital parameters can lead to a good recombination. Article number, page 6 of 12

7 M. Nowak, H. Le Coroller, et al.: Keplerian image recombination for the direct detection of exoplanets Fig. 5. Resulting distribution of the planet candidates found and missed among the 33 simulations were a planet was effectively injected, when using only a S/N > 5 treshold to flag planet candidates. Fig. 7. Mean distance between the true position of the planet and the position found by K-Stacker as a function of the total path length travelled by the planet, for the 25 simulations in which the planet has been detected. inated by the second effect, suggesting that the reoptimization stage could still be improved. In Figure 8, we show the error on each of the 6 parameters individually, as a function of the total path length travelled by the planet. On some parameters (e.g. a, i, t 0 ), there seems to be a clear advantage of having longer displacements of the planet to get better estimates. On some others (e. g. e, Ω, θ 0 ), it is not as clear, though. These parameters, which are not better fitted for longer travelled path, may be responsible for the overall increasing slope of Figure Possible improvements of the algorithm Fig. 6. Same as Figure 5, but for a S/N threshold of 7 In Figure 7, we give, for each of these 25 simulations in which the planet was found, the mean distance between the true position of the planet in each image and the one predicted according to the set of parameters optimized by the K-Stacker algorithm, as a function of the total path travelled, and of the S/N ratio. The two points with a mean error of more than 3 pixels correspond to two simulations in which the planet is very close to the edge of the AO corrected area in certain images of the sequence, reducing the ability of the algorithm to properly constrain the orbit. A linear fit on the other points reveal that the mean distance error is slightly increasing with the total path length travelled by the planet ( pixel per pixel travelled by the planet). When then planet travels along a large portion of its orbit in the sequence, we expect the orbit to be better constrained, and the algorithm to find a better fit to it. But at the same time, a fixed error on the orbital elements has a much bigger impact on the estimated positions if the planet moves along a large portion of its orbit. The result of Figure 7 shows that we are still dom- We find that the current version of the K-Stacker algorithm can achieve a 50 % detection rate on targets with a total S/N 5. Whereas this result already shows that it is possible to use K- Stacker to find very faint planets (S/N below 2 in individual frames for a serie of 10 observations), we strongly believe that there is still room for improvement in the algorithm. One sure way to improve the results of the K-Stacker algorithm would be to optimize the search-grid. For example, a larger steps in θ 0 could be used for smaller values of a. In our version of the algorithm, we are exploring the same range of values for t 0 no matter what the actual value of a is. But it is well known that the semi major-axis and the orbital period are directly related, and this could be used to reduce the range explored for smaller values of a. In fact, an optimal grid may be constructed for the K-Stacker algorithm by looking for a set of orbit which diverge by at least one FWHM of the instrument PSF in at least one image of the sequence. This could be done numerically, or maybe even analytically. This would ensure that the minimum possible number of orbits are used, and hence would help in reducing the total computing time, and the false alarm rate. Also, one area which has yet to be studied in details is the definition of the function to be optimized usually refered to as the gain function in optimization problems. As pointed out in Sect. 5.2, when the algorithm is trying to recombine speckles to create a super-speckle, the resulting pattern is blurred and Article number, page 7 of 12

8 A&A proofs: manuscript no. papier Fig. 8. The upper 6 panels show the errors made one the orbital parameters as a function of the total path travelled by the planet, for each of the simulation in which a planet has correctly been identified. In these panels, two cases, respectively corresponding to a good fit of the orbit (green dots) and to a bad fit (red dots) are highlighted. The corresponding orbits are shown in the two lower panels, with the real orbit of the planet in black, and the orbit found by K-Stakcer in green (good fit case), and red (bad fit case). For each orbit, the position of the planet in all the 10 images is also represented. does not ressemble a real PSF. The most efficient gain function might be one that takes into account the S/N but also the shape of the recombined spot. When searching for planets close to the star, small sample statistics should be taken into account to compute the S/N (Mawet et al. 2014). Along this line of thought, we also tried to use a noiseweighted averaging of the different images (close to what is suggested in Bottom et al., 2017), to take into account the fact that in our series, because of the motion of the planet with respect to the central star, some exposures might be better than others (e.g. when the planet is far away from the star). This modification did not drastically changed our results. The same blind test lead to 25 planets recovered out of 34 (instead of the 24 in the version presented here), with no false negative instead of one found here. A more sophisticated algorithm will be necessary to really improve the results. We also noted an unexpected behavior of the algorithm when using the noise-weighted average. This behavior is examplified in one particular simulation, in which the planet was found by the algorithm, but the orbit was not properly recovered (see Figure 9). An in-depth investigation revealed that this was due to some sort of interaction between the way the planets are injected in our images, and the noise-weighted averaging. We recall that in our images, the planet is injected at constant flux, rather than constant S/N, to reflect the reality of a constant exposure time. This means that from one image to the other, the actual S/N may vary. In this particular case, the planet was injected at a higher S/N in the first few images of the sequence than in the others Article number, page 8 of 12

9 M. Nowak, H. Le Coroller, et al.: Keplerian image recombination for the direct detection of exoplanets (S/N 1.9 vs 1.2). Because of the noise structure in high-contrast images, the noise-averaged combination also has a natural tendancy to favor images in which the computed position is far from the central star. The combination of these two factors made for a situation in which a highly eccentric orbit fitting correctly the first part of the real orbit and moving rapidly towards the central star in the last part (as in Figure 9, bottom panel) yields a good total S/N. The classical combination algorithm, which gives the same importance to each image of the sequence, is less subject to this effect. Finally, it has to be noted that the K-Stacking method does not necessarily implies the use of a brute-force optimization algorithm. Any type of optimization algorithm could be used to optimize the S/N function: simulated annealing, genetic algorithms, amoeba, etc. However, we believe that the brute-force method is one of the most appropriate for K-Stacker. Despite not being known for its efficiency, this method has a very interesting property: the sets of orbital parameters on which the S/N has to be calculated are known in advance and never change (these are the points of the search grid). This means that adding a new image to a set of n observations already processed by K- Stacker do not take much more computation time if the signal and noise terms computed in each image and for each point of the search grid are systematically stored. In this case, these terms can be combined to generate the S/N values using the equation S/N = k S k / k σ 2 k. Then to add a n + 1th image to a set of n images, one can just compute the values of S n+1 and σ n+1, and recompute the S/N. The new image can be processed separately and it is not necessary to reprocess the entire set of n + 1 images. On the opposite, when using Monte-Carlo methods, in which the optimization path is generated dynamically by the algorithm based on previous values found, adding a new image to a set of already processed observations does require a complete re-processing of the whole set, which can prove extremely timeconsuming. On the other hand, a Monte-Carlo type algorithm would explore more thoroughly the local search space around the best orbits, helping in determining error bars on the orbital parameters. The best way to proceed might be a combination of both a grid search to reduce the search space, and an MCMC (Monte-Carlo Markov Chain) based algorithm. Fig. 9. Illustration of an unexpected behavior of the K-Stacker algorithm when using a noise-weighted averaged combination of images. The top panel shows the best orbit found by the classical combination algorithm, and the bottom panel shows the best orbit found when using the noiseweighted averaged combination, for the same simulation. In this case, the planet has a S/N 1.9 in the first 4 images, and 1.2 in the last 6 images. The real orbit of the planet is in black. 6. Conclusion: what can we use K-Stacker for? In this study, we have shown how the K-Stacking method could be used to combine multiple high-contrast images obtained during different nights to detect exoplanets. We simulated SPHERE/IRDIS coronagraphic images taken over several months with planets at low S/N level (below the detection limit) in individual frames, and used the K-Stacker algorithm to combine these images and detect the planets. We find that when the total number n of available images is large enough to get n S/Nind 7 (where S/N ind is the S/N levels in individual frames), the K-Stacker algorithm is able to detect the planet with a very high level of reliability (> 90%). For cases where 5 n S/N ind 6, we find a recovery rate of about 50 %. The detection rate can probably be improved by refining and optimizing the search grid. This will be done in a future work, using real data. This first study was limited to simulated images in order to have complete control over all the important parameters (star magnitude, turbulence conditions, XAO performances, etc.). It has shown that K-Stacker could achieve good recovery rates (that is: simillar to what usual high-contrast imaging techiques are expected to provide) on simulated data, in a decent amount of time. Future studies on this subject should be done using real data, possibly with flase planets injected. In its current implementation, K-Stacker still heavily relies on the ability of the observer to identifiy false positives in the recombined images. Whereas our blind test has shown that this could be done reliably (only one false positive, which was identified as an uncertain planet candidate by the observer, among our 49 simulations), it is clear that this has to be tested using real data. In its current state, though, the K-Stacker technique could already prove useful to detect very faint planets. The ADI tech- Article number, page 9 of 12

11 M. Nowak, H. Le Coroller, et al.: Keplerian image recombination for the direct detection of exoplanets (a) (b) (c) (d) Article number, page 11 of 12

12 A&A proofs: manuscript no. papier (e) (f) (g) Fig. A.1. A set of best recombined images as found by K-Stacker, for different reoptimized S/N values, and for cases where the algorithm found a true/false planet candidate. In each panel, the blue arrow show the position of the plantet candidate found by K-Stacker. Panel 1a (resp. 1b): S/N = 5.7 with a correct (resp. false) planet candidate; Panel 1c (resp. 1d): S/N = 6.0 with a correct (resp. false) planet candidate; Panel 1e (resp. 1f): S/N = 6.5 with a correct (resp. false) planet candidate; Panel 1g: S/N = 8.3, with a true planet candidate. Note that the observer has been able to correctly identify false alarm (S/N > 6 without planets injected in simulations: 1b, 1d, and 1e) using the shape of the spot, without any a priori information on the planet injection. Article number, page 12 of 12

### arxiv: v1 [astro-ph.im] 10 Oct 2017

Experimental parametric study of the Self-Coherent Camera arxiv:1710.03520v1 [astro-ph.im] 10 Oct 2017 Johan Mazoyer a, Pierre Baudoz a, Marion Mas ab, Gerard Rousset a,raphaël Galicher cd a LESIA, Observatoire

### End to End Simulation of AO-assisted coronagraphic differential imaging: estimation of performance for SPHERE

End to End Simulation of AO-assisted coronagraphic differential imaging: estimation of performance for SPHERE Anthony Boccaletti a, Marcel Carbillet b, Thierry Fusco c, David Mouillet d, Maud Langlois

### Techniques for direct imaging of exoplanets

Techniques for direct imaging of exoplanets Aglaé Kellerer Institute for Astronomy, Hawaii 1. Where lies the challenge? 2. Contrasts required for ground observations? 3. Push the contrast limit Recycle!

### Suppressing Speckle Noise for Simultaneous Differential Extrasolar Planet Imaging (SDI) at the VLT and MMT

Direct Imaging of Exoplanets: Science and Techniques Proceedings IAU Colloquium No. 200, 2005 C. Aime and F. Vakili, Eds. c 2005 International Astronomical Union DOI: 00.0000/X000000000000000X Suppressing

### Date of delivery: 29 June 2011 Journal and vol/article ref: IAU Number of pages (not including this page): 5

Date of delivery: 29 June 2011 Journal and vol/article ref: IAU 1101538 Number of pages (not including this page): 5 Author queries: Typesetter queries: Non-printed material: The Physics of the Sun and

### Introduction to exoplanetology

Introduction to exoplanetology Lesson #6: Direct exoplanet detection methods (1/2) Olivier Absil absil@astro.ulg.ac.be Outline I. Direct detection: why and how? II. High contrast imaging from the ground

### Ground-Layer Adaptive Optics Christoph Baranec (IfA, U. Hawai`i)

Ground-Layer Adaptive Optics Christoph Baranec (IfA, U. Hawai`i) Photo credit: T. Stalcup What is Ground-layer Adaptive Optics (GLAO)? Benefits of GLAO to astronomy. MMT multiple-laser AO system. Ground-layer

### THE SPHERE VIEW OF BETELGEUSE

Title : will be set by the publisher Editors : will be set by the publisher EAS Publications Series, Vol.?, 2013 THE SPHERE VIEW OF BETELGEUSE O. Chesneau 1, H.-M. Schmid 2, M. Carbillet 3, A. Chiavassa

### Infrared Differential Imager and Spectrograph for SPHERE: Performance Status with Extreme Adaptive Optics before shipment to ESO/VLT

Infrared Differential Imager and Spectrograph for SPHERE: Performance Status with Extreme Adaptive Optics before shipment to ESO/VLT M. Langlois a, A. Vigan b,, C. Moutou b, K. Dohlen b, A. Costille c,

### Exoplanetary transits as seen by Gaia

Highlights of Spanish Astrophysics VII, Proceedings of the X Scientific Meeting of the Spanish Astronomical Society held on July 9-13, 2012, in Valencia, Spain. J. C. Guirado, L.M. Lara, V. Quilis, and

### An Introduction to. Adaptive Optics. Presented by. Julian C. Christou Gemini Observatory

An Introduction to Adaptive Optics Presented by Julian C. Christou Gemini Observatory Gemini North in action Turbulence An AO Outline Atmospheric turbulence distorts plane wave from distant object. How

### Response of DIMM turbulence sensor

Response of DIMM turbulence sensor A. Tokovinin Version 1. December 20, 2006 [tdimm/doc/dimmsensor.tex] 1 Introduction Differential Image Motion Monitor (DIMM) is an instrument destined to measure optical

### A Fast Algorithm for Cosmic Rays Removal from Single Images

A Fast Algorithm for Cosmic Rays Removal from Single Images Wojtek Pych David Dunlap Observatory, University of Toronto P.O. Box 360, Richmond Hill, Ontario, Canada L4C 4Y6 and Copernicus Astronomical

### Image requirements for the detection and caracterisation of Earth-like planets and hot Jupiter s

Image requirements for the detection and caracterisation of Earth-like planets and hot Jupiter s Alain Chelli Laboratoire d Astrophysique de Grenoble & Jean-Marie Mariotti Center Content The Optical system

### SPHERE IFS: the spectro differential imager of the VLT for Exoplanets search

SPHERE IFS: the spectro differential imager of the VLT for Exoplanets search R.U. Claudi *a, M. Turatto b, E. Giro a, D. Mesa a, U. Anselmi a, P. Bruno b, E. Cascone c, V. De Caprio d, S. Desidera a, R.

### Analyzing Spiral Galaxies Observed in Near-Infrared

Analyzing Spiral Galaxies Observed in Near-Infrared Preben Grosbøl European Southern Observatory Karl-Schwarzschild-Str. 2, D-85748 Garching, Germany Abstract A sample of 54 spiral galaxies was observed

### Planet Detection. AST 105 Intro Astronomy The Solar System

Review AST 105 Intro Astronomy The Solar System MIDTERM III this THURSDAY 04/8 covering LECT. 17 through We ve talked about the Terrestrial Planets and the Jovian Planets - What about planets around other

### Sky Projected Shack-Hartmann Laser Guide Star

Sky Projected Shack-Hartmann Laser Guide Star T. Butterley a, D.F. Buscher b, G. D. Love a, T.J. Morris a, R. M. Myers a and R. W. Wilson a a University of Durham, Dept. of Physics, Rochester Building,

### The Wolf-Rayet + O binary WR 140 in Cygnus

The Wolf-Rayet + O binary WR 140 in Cygnus http://spektroskopie.fg-vds.de Fachgruppe SPEKTROSKOPIE 1. The system The archetype of colliding-wind binary (CWB) systems is the 7.9-year period WR+O binary

### Detection of Artificial Satellites in Images Acquired in Track Rate Mode.

Detection of Artificial Satellites in Images Acquired in Track Rate Mode. Martin P. Lévesque Defence R&D Canada- Valcartier, 2459 Boul. Pie-XI North, Québec, QC, G3J 1X5 Canada, martin.levesque@drdc-rddc.gc.ca

### AstroBITS: Open Cluster Project

AstroBITS: Open Cluster Project I. Introduction The observational data that astronomers have gathered over many years indicate that all stars form in clusters. In a cloud of hydrogen gas, laced with helium

### Brown Dwarfs and Planets around Nearby Stars. A Coronagraphic Search for. B. R. Oppenheimer (AMNH) Morino,, H. Suto,, M. Ishii, M. K.

A Coronagraphic Search for Brown Dwarfs and Planets around Nearby Stars T. Nakajima, J.-I. Morino,, H. Suto,, M. Ishii, M. Tamura, N. Kaifu,, S. Miyama, H. Takami,, N. Takato, S. Oya, S. Hayashi, M. Hayashi

### Astronomie et astrophysique pour physiciens CUSO 2015

Astronomie et astrophysique pour physiciens CUSO 2015 Instruments and observational techniques Adaptive Optics F. Pepe Observatoire de l Université Genève F. Courbin and P. Jablonka, EPFL Page 1 Adaptive

### Modern Image Processing Techniques in Astronomical Sky Surveys

Modern Image Processing Techniques in Astronomical Sky Surveys Items of the PhD thesis József Varga Astronomy MSc Eötvös Loránd University, Faculty of Science PhD School of Physics, Programme of Particle

### Planets Around M-dwarfs Astrometric Detection and Orbit Characterization

Planets Around M-dwarfs Page of 7 Planets Around M-dwarfs Astrometric Detection and Orbit Characterization N. M. Law (nlaw@astro.caltech.edu), S. R. Kulkarni, R. G. Dekany, C. Baranec California Institute

### EXONEST The Exoplanetary Explorer. Kevin H. Knuth and Ben Placek Department of Physics University at Albany (SUNY) Albany NY

EXONEST The Exoplanetary Explorer Kevin H. Knuth and Ben Placek Department of Physics University at Albany (SUNY) Albany NY Kepler Mission The Kepler mission, launched in 2009, aims to explore the structure

### arxiv: v1 [astro-ph.ep] 18 Jun 2012

Astronomy & Astrophysics manuscript no. paper c ESO 2012 June 20, 2012 The International Deep Planet Survey I. The frequency of wide-orbit massive planets around A-stars A. Vigan 1, J. Patience 1,2, C.

### The Potential of Ground Based Telescopes. Jerry Nelson UC Santa Cruz 5 April 2002

The Potential of Ground Based Telescopes Jerry Nelson UC Santa Cruz 5 April 2002 Contents Present and Future Telescopes Looking through the atmosphere Adaptive optics Extragalactic astronomy Planet searches

### Open Cluster Research Project

Open Cluster Research Project I. Introduction The observational data indicate that all stars form in clusters. In a cloud of hydrogen gas, laced with helium and a trace of other elements, something triggers

### SPIE The Differential Tip-Tilt Sensor of SPHERE

SPIE The Differential Tip-Tilt Sensor of SPHERE Pierre Baudoz a,b, Reinhold Dorn c, Jean-Louis Lizon c, Thierry Fusco d, Kjetil Dohlen e, Julien Charton f,, Jean-Luc Beuzit f, Pascal Puget f, David Mouillet

### DATA LAB. Data Lab Page 1

NOTE: This DataLab Activity Guide will be updated soon to reflect April 2015 changes DATA LAB PURPOSE. In this lab, students analyze and interpret quantitative features of their brightness graph to determine

### A Fast Algorithm for Cosmic-Ray Removal from Single Images

Publications of the Astronomical Society of the Pacific, 116:148 153, 2004 February 2004. The Astronomical Society of the Pacific. All rights reserved. Printed in U.S.A. A Fast Algorithm for Cosmic-Ray

### Reading list... Fischer et al. chapter in PPVI, Exoplanet Detection Techniques. Chapter 2 in Exoplanets (Fischer & Lovis) pg

Reading list... Fischer et al. chapter in PPVI, Exoplanet Detection Techniques Chapter 2 in Exoplanets (Fischer & Lovis) pg 27-53. (RV method) Radial Velocity (RV) Detection: Early 20th century RV precision

### The table summarises some of the properties of Vesta, one of the largest objects in the asteroid belt between Mars and Jupiter.

Q1.(a) The table summarises some of the properties of Vesta, one of the largest objects in the asteroid belt between Mars and Jupiter. Diameter / m Distance from the Sun / AU smallest largest 5.4 10 5

### Hubble's Law and the Age of the Universe

Hubble's Law and the Age of the Universe Procedure: Name: 1. Login into the network using your user ID and your password. 2. Double click on the Astronomy shortcuts folder on the desktop. 3. Double click

### Introduction to Telescopes Pre-lab

AST 114 Spring 2005 Introduction to Telescopes Introduction to Telescopes Pre-lab 1. Read through the "Parts of a Telescope" section of the lab. You may be asked to point out the different parts of the

### NA LASER GUIDE STAR AO WITH DYNAMICAL REFOCUS

Florence, Italy. Adaptive May 2013 Optics for Extremely Large Telescopes III ISBN: 978-88-908876-0-4 DOI: 10.12839/AO4ELT3.13893 NA LASER GUIDE STAR AO WITH DYNAMICAL REFOCUS Sebastian Rabien 1,a, Fernando

### arxiv: v1 [astro-ph.im] 17 Mar 2014 F. Millour 1, R. Petrov 1, S. Lagarde 1, P. Berio 1, Y. Bresson 1, L. Abe 1

Improving the performances of current optical interferometers & future designs Proceedings of Haute Provence Observatory Colloquium (23-27 September 2013) Edited by L. Arnold, H. Le Coroller & J. Surdej

### arxiv: v1 [astro-ph.ep] 18 Feb 2009

Astronomy & Astrophysics manuscript no. 1531ms c ESO 2016 November 20, 2016 arxiv:0902.2997v1 [astro-ph.ep] 18 Feb 2009 Letter to the Editor Bayesian analysis of the radial velocities of HD 11506 reveals

### Exoplanet Microlensing Surveys with WFIRST and Euclid. David Bennett University of Notre Dame

Exoplanet Microlensing Surveys with WFIRST and Euclid David Bennett University of Notre Dame Why Space-based Microlensing? Space-based microlensing is critical for our understanding of exoplanet demographics

### Designing a Space Telescope to Image Earth-like Planets

Designing a Space Telescope to Image Earth-like Planets Robert J. Vanderbei Rutgers University December 4, 2002 Page 1 of 28 Member: Princeton University/Ball Aerospace TPF Team http://www.princeton.edu/

### Fundamentals of Exoplanet Observing

Fundamentals of Exoplanet Observing Dennis M. Conti Chair, AAVSO Exoplanet Section Copyright Dennis M. Conti 2017 1 The Strange World of Exoplanets Most exoplanets we have discovered are close-in, large

### Fundamentals of Exoplanet Observing

Fundamentals of Exoplanet Observing Dennis M. Conti Chair, AAVSO Exoplanet Section Copyright Dennis M. Conti 2017 1 The Strange World of Exoplanets Most exoplanets we have discovered are close-in, large

### Lecture 12: Distances to stars. Astronomy 111

Lecture 12: Distances to stars Astronomy 111 Why are distances important? Distances are necessary for estimating: Total energy released by an object (Luminosity) Masses of objects from orbital motions

### MASS DETERMINATIONS OF POPULATION II BINARY STARS

MASS DETERMINATIONS OF POPULATION II BINARY STARS Kathryn E. Williamson Department of Physics and Astronomy, The University of Georgia, Athens, GA 30602-2451 James N. Heasley Institute for Astronomy, University

### Astronomy. Astrophysics. Deriving structural parameters of semi-resolved star clusters. FitClust: a program for crowded fields

DOI: 10.1051/0004-6361/201322577 c ESO 2014 Astronomy & Astrophysics Deriving structural parameters of semi-resolved star clusters FitClust: a program for crowded fields D. Narbutis 1, D. Semionov 1, R.

### Extrasolar planets. Lecture 23, 4/22/14

Extrasolar planets Lecture 23, 4/22/14 Extrasolar planets Extrasolar planets: planets around other stars Also called exoplanets 1783 exoplanets discovered as of 4/21/14 Orbitting 1105 different stars Number

### HIGH-CONTRAST IMAGING OF YOUNG PLANETS WITH JWST

ExSoCal Conference September 18-19th 2017, Pasadena HIGH-CONTRAST IMAGING OF YOUNG PLANETS WITH JWST Marie Ygouf HR 8799 planetary system Jason Wang / Christian Marois Keck data Giants planets and brown

### The Mass of Jupiter Student Guide

The Mass of Jupiter Student Guide Introduction: In this lab, you will use astronomical observations of Jupiter and its satellites to measure the mass of Jupiter. We will use the program Stellarium to simulate

### Exhaustive strategy for optical survey of geosynchronous region using TAROT telescopes

Exhaustive strategy for optical survey of geosynchronous region using TAROT telescopes Pascal Richard, Carlos Yanez, Vincent Morand CNES, Toulouse, France Agnès Verzeni CAP GEMINI, Toulouse, France Michel

### Virgo Ego Scientific Forum VESF school 2011 Optical Image Analysis Tutorial by Marica Branchesi

Virgo Ego Scientific Forum VESF school 2011 Optical Image Analysis Tutorial by Marica Branchesi The main goal of present laboratory is to introduce the students to Optical Image Analysis processing for

### Field Rotation in Altitude over Azimuth Mounts and Its Effect on CCD Imaging What is the Maximum Exposure? Bill Keicher.

Field Rotation in Altitude over Azimuth Mounts and Its Effect on CCD Imaging What is the Maximum Exposure? Bill Keicher 12 March 2005 When planning astronomical imaging sessions with a telescope equipped

### Astronomy 421. Lecture 8: Binary stars

Astronomy 421 Lecture 8: Binary stars 1 Key concepts: Binary types How to use binaries to determine stellar parameters The mass-luminosity relation 2 Binary stars So far, we ve looked at the basic physics

### Report on the new EFOSC2 VPH grisms

Report on the new EFOSC2 VPH grisms Ivo Saviane Lorenzo Monaco v 1.0 March 01, 2008 1 Introduction In January 2008 the ULTRASPEC project delivered two volume-phased holographic grisms (VPHG) to be used

### Source Detection in the Image Domain

Source Detection in the Image Domain Astro 6525, Fall 2015 Shami Chatterjee Oct 2015 Source Identification Identify a source in the presence of noise. Localize it as best as possible given instrumental

### DOME C AS A SETTING FOR THE PERMANENT ALL SKY SURVEY (PASS)

Title : will be set by the publisher Editors : will be set by the publisher EAS Publications Series, Vol.?, 2005 DOME C AS A SETTING FOR THE PERMANENT ALL SKY SURVEY (PASS) H.J. Deeg, J.A. Belmonte, R.

### MIRI, METIS and the exoplanets. P.O. Lagage CEA Saclay

MIRI, METIS and the exoplanets P.O. Lagage CEA Saclay French Co-PI of JWST-MIRI and Coordinator of European MIRI GTO on exoplanets Member of the ELT-METIS science team Why MIRI and METIS? Because of my

### ABSTRACT 1. INTRODUCTION

FASS: The full aperture seeing sensor A.Guesalaga *,1, J.Osborn 2, M.Sarazin 3, B.Neichel 4, S. Perera 2, R.Wilson 2 1 Pontificia Universidad Católica de Chile, 4860 Vicuña Mackenna, 7820436 Santiago,

### Reduced data products in the ESO Phase 3 archive (Status: 02 August 2017)

Reduced data products in the ESO Phase 3 archive (Status: 02 August 2017) The ESO Phase 3 archive provides access to reduced and calibrated data products. All these data are stored in standard formats

### Infra-red imaging of perpendicular nested bars in spiral galaxies with the Infra-red Camera at the Carlos Sanchez Telescope

Infra-red imaging of perpendicular nested bars in spiral galaxies with the Infra-red Camera at the Carlos Sanchez Telescope S.N. Kemp (skemp@ll.iac.es) Instituto de Astrofísica de Canarias, E-38200 La

### NACO-SDI DIRECT IMAGING SEARCH FOR THE EXOPLANET ERI b 1

The Astronomical Journal, 133:2442 2456, 2007 June # 2007. The American Astronomical Society. All rights reserved. Printed in U.S.A. NACO-SDI DIRECT IMAGING SEARCH FOR THE EXOPLANET ERI b 1 Markus Janson,

### SPHERE ZIMPOL: Overview and performance simulation

SPHERE ZIMPOL: Overview and performance simulation Christian Thalmann a, Hans M. Schmid a, Anthony Boccaletti b, David Mouillet c, Kjetil Dohlen d, Ronald Roelfsema e, Marcel Carbillet f, Daniel Gisler

### Measuring The Mass of Uranus by Observing the Revolution of the Moons Using Kepler s Third Law

Measuring The Mass of Uranus by Observing the Revolution of the Moons Using Kepler s Third Law Ryan Rubenzahl AST Final Project December 7th, 204 Abstract In this project we determined the mass of Uranus

### DETECTING TRANSITING PLANETS WITH COROT. Stefania Carpano ESAC (25th of November 2009)

DETECTING TRANSITING PLANETS WITH COROT Stefania Carpano ESAC (25th of November 2009) Outline 1) Introduction: Exoplanet science and detection methods Detection of exoplanets using the transit method The

### Scott Gaudi The Ohio State University Sagan Summer Workshop August 9, 2017

Single Lens Lightcurve Modeling Scott Gaudi The Ohio State University Sagan Summer Workshop August 9, 2017 Outline. Basic Equations. Limits. Derivatives. Degeneracies. Fitting. Linear Fits. Multiple Observatories.

### How to measure star formation rates in galaxies?

How to measure star formation rates in galaxies? Véronique Buat 1 1. Aix-Marseille Université, CNRS, LAM (Laboratoire d Astrophysique de Marseille) UMR7326, 13388, Marseille, France Our knowledge of cosmic

### Lab 7: The H-R Diagram of an Open Cluster

Lab 7: The H-R Diagram of an Open Cluster Due Date: 2007 Nov 27 (after thanksgiving) 1 Introduction: The HR Diagram In this two week project you will do absolute (not differential) photometry with a CCD

### Detecting Cosmic Rays in Infrared Data

The 2010 STScI Calibration Workshop Space Telescope Science Institute, 2010 Susana Deustua and Cristina Oliveira, eds. Detecting Cosmic Rays in Infrared Data Rachel E. Anderson and Karl Gordon Space Telescope

### High accuracy transit photometry of the planet OGLE-TR-113b with a new deconvolution-based method ABSTRACT

A&A 459, 249 255 (2006) DOI: 0.05/0004-636:20065844 c ESO 2006 Astronomy & Astrophysics High accuracy transit photometry of the planet OGLE-TR-3b with a new deconvolution-based method M. Gillon,2,F.Pont,

### Extraction of Point Source Spectra from STIS Long Slit Data

1997 HST Calibration Workshop Space Telescope Science Institute, 1997 S. Casertano, et al., eds. Extraction of Point Source Spectra from STIS Long Slit Data J. R. Walsh Spect Telescope European Coordinating

### EXOPLANET DISCOVERY. Daniel Steigerwald

EXOPLANET DISCOVERY Daniel Steigerwald WHAT IS AN EXOPLANET? An exoplanet is a planet outside of our solar system Extrastellar Rogue 1853 Planets 1162 planetary systems 473 Multiple planetary systems HISTORY

### THE GALACTIC BULGE AS SEEN BY GAIA

143 THE GALACTIC BULGE AS SEEN BY GAIA C. Reylé 1, A.C. Robin 1, M. Schultheis 1, S. Picaud 2 1 Observatoire de Besançon, CNRS UMR 6091, BP 1615, 25010 Besançon cedex, France 2 IAG/USP Departamento de

### Searching for Other Worlds: The Methods

Searching for Other Worlds: The Methods John Bally 1 1 Center for Astrophysics and Space Astronomy Department of Astrophysical and Planetary Sciences University of Colorado, Boulder The Search Extra-Solar

### Pole searching algorithm for Wide-field all-sky image analyzing monitoring system

Contrib. Astron. Obs. Skalnaté Pleso 47, 220 225, (2017) Pole searching algorithm for Wide-field all-sky image analyzing monitoring system J. Bednář, P. Skala and P. Páta Czech Technical University in

### PHYS133 Lab 4 The Revolution of the Moons of Jupiter

PHYS133 Lab 4 Goals: Use a simulated remotely controlled telescope to observe iter and the position of its four largest moons. Plot their positions relative to the planet vs. time and fit a curve to them

### Image Processing in Astronomy: Current Practice & Challenges Going Forward

Image Processing in Astronomy: Current Practice & Challenges Going Forward Mario Juric University of Washington With thanks to Andy Connolly, Robert Lupton, Ian Sullivan, David Reiss, and the LSST DM Team

### Ay 1 Lecture 2. Starting the Exploration

Ay 1 Lecture 2 Starting the Exploration 2.1 Distances and Scales Some Commonly Used Units Distance: Astronomical unit: the distance from the Earth to the Sun, 1 au = 1.496 10 13 cm ~ 1.5 10 13 cm Light

### Astr As ome tr tr ome y I M. Shao

Astrometry I M. Shao Outline Relative astrometry vs Global Astrometry What s the science objective? What s possible, what are fundamental limits? Instrument Description Error/noise sources Photon noise

### Keck Adaptive Optics Note 1069

Keck Adaptive Optics Note 1069 Tip-Tilt Sensing with Keck I Laser Guide Star Adaptive Optics: Sensor Selection and Performance Predictions DRAFT to be updated as more performance data becomes available

### Expected Performance From WIYN Tip-Tilt Imaging

Expected Performance From WIYN Tip-Tilt Imaging C. F. Claver 3 September 1997 Overview Image motion studies done at WIYN show that a significant improvement to delivered image quality can be obtained from

### Extrasolar Planets = Exoplanets III.

Extrasolar Planets = Exoplanets III http://www.astro.keele.ac.uk/~rdj/planets/images/taugruishydra2.jpg Outline Gravitational microlensing Direct detection Exoplanet atmospheres Detecting planets by microlensing:

### How Giovanni s Balloon Borne Telescope Contributed to Today s Search for Life on Exoplanets

How Giovanni s Balloon Borne Telescope Contributed to Today s Search for Life on Exoplanets Wesley A. Jet Propulsion Laboratory, California Ins:tute of Technology Symposium for Giovanni Fazio Harvard Smithsonian

### C 2 n. profilometry from Shack-Hartmann data: model and experiment. 1 Introduction. 2 Problem statement with SH slope and scintillation correlations

C 2 n profilometry from Shack-Hartmann data: model and experiment J. Voyez 1a, C. Robert 1, B. Fleury 1, N. Védrenne 1, V. Michau 1, T. Fusco 1, and E. Samain 2 1 Office National d Etudes et de Recherches

### PIONIER. JB Le Bouquin, Lazareff, Berger, Zins, Traub, Millan-Gabet... for the full PIONIER team

PIONIER JB Le Bouquin, Lazareff, Berger, Zins, Traub, Millan-Gabet... for the full PIONIER team The Very Large Telescope Interferometer UT1 UT2 UT3 UT4 VLTI ELT The context: milliarcsecond angular resolution

### Physical Parameters of KID , a Low-Mass, Double- Lined Eclipsing Binary

Physical Parameters of KID 6131659, a Low-Mass, Double- Lined Eclipsing Binary Gideon Bass Jerome Orosz, William Welsh, Gur Windmiller, Trevor Gregg, Tara Fetherolf, Richard Wade, Samuel Quinn Paper available

### Geostationary Satellites and Astrophotography

Geostationary Satellites and Astrophotography By Clay S. Turner 11/14/08 Introduction When I first started doing astrophotography last January, it was only natural to take pictures of some of Orion s deep

### MAORY (Multi conjugate Adaptive Optics RelaY) for E-ELT. Paolo Ciliegi. On behalf of the MAORY Consortium

MAORY (Multi conjugate Adaptive Optics RelaY) for E-ELT Paolo Ciliegi INAF Osservatorio Astronomico di Bologna On behalf of the MAORY Consortium Science ELT Workshop Team Meeting ESO Garching, MPE Garching,

### The asteroids. Example for the usage of the Virtual Observatory

Example for the usage of the Virtual Observatory The asteroids Florian Freistetter, ZAH, Heidelberg florian@ari.uni-heidelberg.de Asteroids in the solar system There are not only planets in our solar system.

### Statistical validation of PLATO 2.0 planet candidates

Statistical validation of PLATO 2.0 planet candidates Rodrigo F. Díaz Laboratoire d Astrophysique de Marseille José Manuel Almenara, Alexandre Santerne, Claire Moutou, Anthony Lehtuillier, Magali Deleuil

### The magnetic properties of Am stars

Contrib. Astron. Obs. Skalnaté Pleso 48, 48 52, (2018) The magnetic properties of Am stars A.Blazère 1,2, P.Petit 3,4 and C.Neiner 2 1 Institut d Astrophysique et de Géophysique, Université de Liège, Quartier

### Basic Theory of Speckle Imaging

1 Basic Theory of Speckle Imaging Elliott Horch, Southern Connecticut State University 1 arcsec BU 151AB 2 Speckle Often Means Binary Stars Stellar Masses. Mass-Luminosity Relation (MLR) Initial Mass Function

### HD Transits HST/STIS First Transiting Exo-Planet. Exoplanet Discovery Methods. Paper Due Tue, Feb 23. (4) Transits. Transits.

Paper Due Tue, Feb 23 Exoplanet Discovery Methods (1) Direct imaging (2) Astrometry position (3) Radial velocity velocity Seager & Mallen-Ornelas 2003 ApJ 585, 1038. "A Unique Solution of Planet and Star

### Polarimetry of exoplanets

Polarimetry of exoplanets Dr. Daphne Stam Aerospace Engineering Technical University Delft The Netherlands This lecture is about Polarimetry of Exoplanets. Exoplanets or extrasolar planets are planets

### UltraPhot, a photometric mode for FLAMES. Françoise Roques LESIA Observatoire de Paris

UltraPhot, a photometric mode for FLAMES Françoise Roques LESIA Observatoire de Paris UltraPhot Science : Variable Objects, Transits, Occultations Temporal resolution is a relatively unexplored region

### Web sites and Worksheets

Introduction Web sites and Worksheets Sandra Woodward Oakhill College Castle Hill (swoodward@oakhill.nsw.edu.au) Astronomy is a practical subject by nature. Unfortunately for teachers, many of the practical

### Gravitational Fields

Gravitational Fields Examples 00 Currently, the space probe, Cassini, is between Jupiter and Saturn. Cassini s mission is to deliver a probe to one of Saturn s moons, Titan, and then orbit Saturn collecting

### Astronomical image reduction using the Tractor

the Tractor DECaLS Fin Astronomical image reduction using the Tractor Dustin Lang McWilliams Postdoc Fellow Carnegie Mellon University visiting University of Waterloo UW / 2015-03-31 1 Astronomical image

### On the FPA infrared camera transfer function calculation

On the FPA infrared camera transfer function calculation (1) CERTES, Université Paris XII Val de Marne, Créteil, France (2) LTM, Université de Bourgogne, Le Creusot, France by S. Datcu 1, L. Ibos 1,Y.