Introduction to exoplanetology

Transcription

1 Introduction to exoplanetology Lesson #6: Direct exoplanet detection methods (1/2) Olivier Absil

2 Outline I. Direct detection: why and how? II. High contrast imaging from the ground III. High contrast imaging from space IV. Main results from high-contrast imaging V. Stellar interferometry

3 I. Direct detection Why and how?

4 Why direct detection? Characterization of planetary atmospheres Needs spectroscopy on actual planetary photons isolate planet from star

5 Why direct detection? Access non-transiting planets Opens up a wide range of separations (beyond a few 0.1 AU) Complementary with transit spectroscopy Transits mostly here

6 Why direct detection? Get full orbital solution 3D orbit > direct, model-independent access to the planet mass Study dynamics of planetary systems, interactions with dust disks, etc.

7 Challenge #1: contrast Visible: reflected light Infrared: thermal (blackbody-type) emission

8 Challenge #2: angular separation 1 10 pc = 0.1 arcsec Theoretically within reach of 10-m class telescope But watch out for turbulence!

9 Challenge summary A firefly close to a lighthouse 1000 miles away! (note: the star never turns off)

10 Technique #1: imaging Diffraction in circular aperture Airy pattern Angular resolution = size of Airy disk: θ = 1.22 λ/d Rayleigh criterion λ = 2 µm, D = 10 m θ = 0.05" (50 mas) = 1 AU at 20 pc Extended pattern noise!

11 The Airy pattern : aperture radius θ = 1.22 λ/d

12 Technique #2: interferometry Two separated telescopes fringes B Angular resolution set by baseline (B): θ = 0.5 λ/b Crests of 1st packet fall on troughs of 2nd packet λ = 2 µm, B = 100 m θ = 2 mas

13 The fringe pattern unresolved binary resolved binary (fringe packets offset for clarity)

14 II. High contrast imaging from the ground

15 Atmospheric windows V I J H K L M N

16 Imaging through the Earth atmosphere Temperature variations Distorted wavefront Short exposure: speckles Long exposure: wide PSF Short exposure Long exposure

17 Loss of angular resolution Fried parameter r 0 : diameter of a circular area over which the rms wavefront aberration due to passage through the atmosphere is equal to 1 radian La Palma r 0 ~ 10 cm at good astronomy site Same resolution as 10 cm telescope! D = 5 cm 10 cm 60 cm 1.2 m

18 Atmospheric turbulence Wavelength dependence: r 0 λ 6/ nm 50 2 µm 4 10 µm Seeing = FWHM of long exposure image Equal to 0.98 λ / r 0 1" seeing for r 0 = nm Varies slowly with wavelength (λ 1/5 )

19 Atmospheric turbulence Coherence time: t 0 = 0.31 r 0 / v wind Valid under Taylor s «frozen turbulence» hypothesis t 0 3 msec for r 0 = 10 cm and v wind = 10 m/s Isoplanatic angle: θ 0 = 0.31 r 0 / h θ 0 1.3" for r 0 = 10 cm and h = 5 km Stars separated by θ 0 have different short-exposure PSFs

20 Correction needed! Adaptive optics

21 Adaptive optics

22 Wavefront sensing Shack-Hartmann wave front sensor

23 Strehl ratio S = exp(iφ) 2 exp( σ φ2 ) φ = wavefront phase σ φ = rms phase on pupil Quantifies image quality peak intensity ratio wrt perfect image Perfect image S = 1 D = r 0 S = 0.36

24 AO correction MACAO on VLT 50% 20%

25 Detection in speckle noise? ESO-3.6m/ComeOn+ (1994) S=0.21 Speckles Planet

26 Getting rid of speckles 1. Observing strategies

27 Solution: PSF subtraction aka Reference-star Differential Imaging (RDI) β Pic Observe a reference star Similar to science target Scale its PSF (flux) α Pic Subtract it from science observation β Pic α Pic

28 Limitations of RDI No perfect reference star Should be of same magnitude, color and position as science target Time spent on reference star is «lost» (no planetary photon) Atmospheric conditions change with time Telescope/instrument aberrations also change with time Quasi-static speckles, can mimic planetary signal

29 Four differential solutions Keep observing the same target (no reference star) Tune an observing parameter to discriminate between stellar PSF and planet Spectral differential imaging (SDI) Spectral deconvolution (SD) Angular differential imaging (ADI) Polarimetric differential imaging (PDI)

30 Spectral differential imaging

31 SDI in practice λ1 λ2 λ3 λ1 λ2 (λ1 λ2) k(λ1 λ3)

32 Pros and cons of SDI + Detection down to diffraction limit (1 λ/d) Companion needs strong molecular absorption (e.g., methane) Differential aberrations between the three wavelength optical paths Narrow filters lower sensitivity

33 Spectral deconvolution Integral Field Spectrograph (IFS) Provides field image as function of wavelength («image cube») Diffraction and speckle pattern scale as function of wavelength Pattern moves out from star with increasing wavelength Exoplanet position is fixed Can be distinguished from speckles

34 SD in practice Similar principle as SDI, on much wider wavelength band Observed (x,λ) slice Rescaled (x,λ) slice

35 Pros and cons of SD + Works with any type of planet (no specific feature) + No differential aberrations / simultaneous observations + Does not rely on specific feature in planet spectrum + End product = spectrum of the planet! Detect and characterize planet at the same time Speckle pattern not perfectly constant over wavelength Limited inner and outer working angles (depend on wavelength range and spectral resolution)

36 Angular differential imaging Use field rotation while keeping telescope fixed Usually done by switching off derotator Planet moves around star as a function of parallactic angle Quasi-static speckles stay at fixed position

37 ADI in practice

38 Pros and cons of ADI + Works with any type of planet (no specific feature) + Does not require specific hardware Does not work well for stars far from zenith (small variation of the parallactic angle) Limited inner working angle (planet must move by more than 1 λ/d in the field for ADI to work) Speckle pattern evolves with time

39 Speckle decorrelation t0 t min (t min) t0 t min (t min) t0

40 Polarimetric differential imaging Reflected light from planet is partially polarized Typically 10% polarisation Star produces unpolarized light Can be exploited by polarimetric imager

41 PDI in practice Fig. 1. Processed images of HD (top row) and HD (bottom row) in H band. From left to right: intensity image (I) in logarithmic

42 Pros and cons of PDI + Speckle subtraction can be very good + No limitation in inner or outer working angle Small fraction of planet light is polarized low sensitivity Works only in reflected light (best in visible range) Requires specific, non-standard hardware

43 Getting rid of speckles 2. Image processing

44 The LOCI algorithm Locally Optimized Combination of Images Goal: make best use of a set of reference images Can be used with various observing techniques Solve a linear system to minimize residuals to be minimized pixels target PSF coefficients to be determined ref PSF individual frames

45 LOCI working with ADI

46 Practical use of LOCI To be applied locally because correlation between frames expected to depend on position Optimize LOCI coefficient on «optimization» zone Perform PSF subtraction on «subtraction» zone Optim. zone > subtr. zone to avoid planet signal subtraction

47 Pros and cons of LOCI Several free parameters + Can be fine tuned to any particular data set Not very straightforward to use Strong self-subtraction of planetary signal LOCI tries to remove everything bias on photometry Can be evaluated with fake companions, or mitigated by introducing even more parameters (masking) CPU intensive

48 Principal Component Analysis (PCA) Method to reduce the number of variables needed to describe a data set Build an orthogonal basis on which to decompose the images Truncated basis to perform PSF subtraction

49 How PCA works Subtract mean (M) from data set (D): A = D M Calculate covariance matrix: C = AA T Calculate eigenvectors (Q) and eigenvalues CQ = QΛ, with Λ a diagonal matrix of eigenvalues Choose first K components in Q Q K Project data onto these components F = Q K A new data set of reduced dimensions

50 PCA with a stack of images Use the women s face to compute orthogonal basis Project man s face on first 50 coefficients and add Describe man s face with 50 coefficient instead of thousands of pixels

51 Practical use of PCA

52 Pros and cons of PCA + Can be applied to whole image at once + CPU time can be drastically reduced! In practice, annulus-wise version is generally better Self-subtraction reduced, but still present + Reduced bias wrt LOCI, generally more linear Fake companion injection still used in most cases + Forward modeling also possible

53 A «real-life» example HR 8799

54 2008: ADI detection Contrast ~ 14 mag (~10 6 )

55 4 th planet with ADI+LOCI

56 2010: long slit spectroscopy Figure 2. Image with the HR 8799 c spectrum before extraction. The spatial direction along the star planet axis is vertical. The spectral direction is horizontal, with wavelength increasing from left to right. ANSON ET AL. Vol % 400 s, taken nt. A plied 799 c re 1). re an on of c and his is tions. ound ation were l sky. dures ected ringe e flat Since on on to be aging speco this ectrorame on of Figure 4. into three (circles, d dash-dotte No. 1, 2010 SPECTROSCOPY OF HR 8799 c Figure 2. Image with the HR 8799 c spectrum before extraction. The spatial direction along the star planet axis is vertical. The spectral direction is horizontal, with wavelength increasing from left to right. L37 Figure 3. Upper panel: spectrum of HR 8799 c. The dashed lines and faintly shaded area (light blue in the online version) denote the errors. Lower panel: log g = detailed that non Non-equ large dif particula presence 4 µm fo Fortney (2008) t being pr behavio of whet thereby The r planet de ing purp narrowe since it models. true also in the sp 8.66 ma 9.50 ma and is ex In the

57 2011: HK band IFS

58 2013: full IFS study with SD

59 Getting rid of speckles 3. Hardware solutions

60 Coronagraphy

61 Coronagraphy

62 With central obscuration

63 Phase mask coronagraph Proposed by Roddier (1997) Apply 180 phase shift to PSF center Ideal mask size = 43% of 1st Airy ring Very chromatic design

64 Four Quadrant Phase Mask

65 Beyond the FQPM Quadrants octants continuous phase ramp

66 Vortex phase mask Continuous phase ramp from 0 to 4π Main problem: chromaticity of the phase ramp

67 Achromatic phase mask? Sub-wavelength gratings achromatic half wave plate Fig. 3. 4QZOG implementation: s and p are the vectorial complex Annular Groove Phase Mask (AGPM)

68 The AGPM performance 1000" Peak%rejec)on% 100" 10" 1" août(10" nov.(10" févr.(11" juin(11" sept.(11" déc.(11" avr.(12" juil.(12" oct.(12"

69 On-sky example Raw image Post-processed (ADI+PCA)

70 Sensitivity to planets β Pic b

71 Apodization «Pupil plane coronagraphy» Act on amplitude or phase of wavefront in pupil Redistribute the intensity in the focal plane to make it more compact or create a dark region

72 Amplitude apodization Goal: reduce PSF side lobes Degraded angular resolution

73 Apodized coronagraph Apodization makes PSF more compact on the focal plane mask

74 Shaped pupils

75 Microdot mask How it s done

76 Phase apodization Act on wavefront phase instead of amplitude Can produce asymmetric PSF APP: Apodizing Phase Plate

77 On-sky example

78 Phase-induced amplitude apodization

79 Other hardware solutions Low-order wavefront sensing Performed at instrument level (science camera, unexploited stellar light, etc) Corrects for non-common path errors between AO and instrument Speckle nulling Measure the complex amplitude of quasi static speckles and correct with deformable mirror

80 Tree of life

81 All steps in one picture Telescope pointed (i) (i) (i) (i)(i) AO turned on (ii) (ii) (ii) Speckle calibration Centered on mask (ii)(ii) (iii) (iii) (iii) (iii) (iii) ADI + postprocessing b" b" b"b"b" c" c" Project Project Project Project Project c"c"c" Toward Toward Exploration Exploration of Other of Other Worlds Worlds Toward Exploration of Other Worlds Toward Toward Exploration Exploration of of Other Other Worlds Worlds e" e" d" d" e"e"e" d"d"d" N" N" N" N" N" E" (iv) (iv) (iv) (iv) (iv) (v) (v) (v) (v)(v) E"E"E" 0.5 " 0.5 " 0.5 " E" 0.5 "0.5 "