The Impact of EMCCDs on Speckle Imaging

Nic Scott

Research Scientist

BAERI/NASA ARC

nicscott.org

https://www.gemini.edu/node/21236 Credit: Gemini Observatory/NSF/AURA/Artwork by Joy Pollard

Exoplanets can't hide their secrets from innovative new instrument

https://phys.org/news/2019-08-exoplanets-secrets-instrument.html

Space science research highlight: High-resolution imaging and exoplanets

https://www.nasa.gov/feature/space-science-research-highlight-high-resolution-imaging-and-exoplanets

Hidden Secrets of Elusive Exoplanet Revealed by Innovative New Instrument

https://scitechdaily.com/hidden-secrets-of-elusive-exoplanet-revealed-by-innovative-new-instrument/

 

1''

speckle

aperture photometry

blue

red

High-resolution Imaging Transit Photometry of Kepler-13AB -Howell et al.  

2019AJ....158..113H

Both speckle analysis and high speed aperture photometry performed with `Alopeke confirmed that the highly irratiated gas giant, Kepler-13b orbits Kepler-13A.

a single star, a circle representing the isoplanatic patch, and the small star shapes are "speckles"

a single star, a circle representing the isoplanatic patch, and the small star shapes are "speckles"

a single star, the smaller circles show worse seeing and smaller isoplanatic patches, each producing "speckles"

a single star, a circle representing the isoplanatic patch, and the small star shapes are "speckles"

a single star, the smaller circles show worse seeing and smaller isoplanatic patches, each producing "speckles"

a binary pair, close enough that they share an isoplanatic patch, producing "speckles" that correspond to their separation & position angle

a single star, a circle representing the isoplanatic patch, and the small star shapes are "speckles"

a single star, the smaller circles show worse seeing and smaller isoplanatic patches, each producing "speckles"

a binary pair, close enough that they share an isoplanatic patch, producing "speckles" that correspond to their separation & position angle

a binary pair, wherein the ratio of their separation to the isoplanatic patch size is such that their "speckles" are not correlated

Kepler-13AB

Kepler/K2

TESS

Multiplicity of stars

single :  54%

double : 33%

triple : 8%

higher: 4%

 

 

 

40-50% of Kepler/K2/TESS objects of interest may have companions

 

  • biases occurence rates & underestimates exoplanet radii by 50%

The Rayleigh criterion and the tyranny of the atmosphere

\Theta = 1.22 \lambda/d
\approx \lambda / r_0

but in atmosphere, turbulent cells limit resolution to 

Fizeau & Michelson

Brown & Twiss

1850s

1890s

1920 - measured Betelgeuse (with Pease)

1956 - HBT effect, correlation b/t coherent photons

detail lost to the atmosphere can be regained through interferometric analysis 

Intensity interferometry

Labeyrie

1970

Text

a method to obtain diffraction-limited resolution across the full aperture of a large telescope

 long exposure

speckles blur

produce Airy pattern

true images are impossible, only centrosymmetric objects can be reconstructed

speckle pattern is the Fourier transform of telescope pupil

autocorrelation of speckles

(in Fourier space)

modulus

(time-averaged intensity)

Knox & Thompson

cross-spectrum, a 2nd order correlation

1974

non-symmetric input

long exposure avg

Labeyrie technique

Knox & Thompson method

mean square of the image transform

modulus of the object transform

autocorrelation of the image transform

phase of the object transform

diffraction-limited image of the object

unambiguous reconstruction of arbitrary shapes but has            ambiguity

+

180\degree

The importance of phase in image reconstruction

1) we're limited to measuring the Fourier magnitude, but phase carries the majority of the information

2) phase is needed for a unique solution for reconstruction

In 1980’s bispectral analysis was found to have higher S/N and be less susceptible to systematic error.

Weigelt & Lohmann

1977-1983

double star simulation

bispectrum modulus

triple correlation

record PSF of object

produce synthetic reference star by shifting the speckle pattern

phase is preserved

Speckle masking/triple correlation theory/bispectral analysis, a 3rd order correlation

deconvolve

  true images

Fringe spacing

Fringe orientation

Fringe Depth

Binary separation

Binary position angle

Binary magnitude difference

1s

40s

20min

An aside on Lucky Imaging

  • reaches higher resolution than typical seeing
  • does not utilize Fourier analysis
  • not capable of reaching the diffraction limit
  • requires a very large number of  images

http://inspirehep.net/record/823349/plots

Fried

1978

record a large series of images

discard instances of poor resolution

combine the remainder through shift-and-add techniques

Can reach high angular resolution and is often confused with speckle imaging.

Changes in technology

ICCDs

EMCCDs

  • CCD + image intensifier tube
  • amplify signal prior to ADC
  • tube gating and amplification
  • up to 500kHz
  • lower resolution and low QE
  • interline detector      50% QE
  • CCD + series of HV gain registers
  • + photoelectric events w/ each register
  • greatly amplified signal w/ single readout of the array
  • 2008: Tokovinin published the first speckle imaging results
  • frame transfer, QE > 90%, low readout noise
  • higher S/N than ICCDs
  • 10-100 times faster than a CCD
  • active pixel sensor: amplification and  conversion happens directly on each pixel
  • cheap, less power and lower gain needed
  • higher S/N, 
  • fast readout
  • no clocking the charge across the array
  • further on-chip processing can be added

CMOS (complementary metal–oxide–semiconductor)

EMCMOS

  • each pixel has a series of HV registers
  • + photoelectric events w/ each register

2015: first EMCMOS results reported by Brugière

Title Text

Title Text

Text

DSSI/NESSI/`Alopeke/Zorro

Speckle

  • 10mas/pxl
  • mag limit ~17
  • contrast limit ~8

Wide Field

  • 73mas/pxl

Speckle

  • 18mas/pxl
  • mag limit ~14
  • contrast limit ~6

Wide Field

  • 81mas/pxl

0.011'' @u

0.026'' @832nm

0.025'' @u

0.060'' @832nm

6.7''

60''

19''

56''

Dual Andor iXon Ultra 888:

  • 1,000x0.06s exp
  • 256x256 ROI frames 
  • sensitivity ~5.5 e-/A/D
  • pixel noise 150e-

EMCCD factors

  • Readnoise is 100200 e (speckle mode)
  • At high EM gain settings,
    • effective readnoise =1 ADU / 1 e 
    • CIC creates spurious counts in individual pixels
    • minimized by using fastest possible readout rates and vertical shift speeds
    • At slow readout speeds, the CIC is 50,00070,000 pixels per image with 6% of the pixels affected.
    • At the highest speeds this drops to 300500 pixels per image (<0.05%).

constant source

non-linear

dome flat PTC

no EM

~5e-/ADU

non-linearity

5.4

65259

Science Programs

58 total proposals

138.3 total nights

(NOAO 2016A to 2019B)

  • ~100n/yr at Gemini

  • ~25n/yr at WIYN

  • 20A has 17 proposals asking for >400 hours of total time at both Gemini's

  • 2-3 :1 subscription

  • Constrain NEA diameters. Image SS objects
  • Determine multiplicity of nearby K and M-dwarfs, does it vary across spectral type?
  • Imaging of brown dwarfs and distant large planets, particularly around M dwarfs
  • Investigate differences in planetary system architectures between multiple vs not (known) multiple host stars
  • Examine long-term RV trends/determine binarity of RV planet hosts.  
  • K2, TESS follow-up
  • Provide an unbiased sample for TESS, so statistical determinations of planet occurrence rates can be made
  • Occultations, transit photometry, pulsar time scales, observe pulsating WDs at high cadence
  • Interacting binary systems, jets, CVs, R Aquarii
  • WR stars - jets, outflows

Some proposals:

HST R Aquarii

Science Results Highlights

  • optical, NIR imaging of host stars of KOIs
  • ~ 90% of the confirmed & candidate exoplanet hosts
  • separations, PA, and dm for all detected, bound and LoS companion stars.
  • 2297 companions around 1903 primary stars
    • ~ 10% obs stars 1+ companions w/i 1”
    • ~ 30% w/i 4
  • correction factors for exoplanet radii caused by the dilution of the transit depth
    • decreases the number of KOI planets with radii smaller than 2 Earth radii by 2% - 23%

Furlan's previous results applied to planets w/ known masses & radii, analyze the effects of a close stellar companion on planetary density.

  • 50 planets orbiting 26 stars in the Kepler field
  • transit dillution         planet radii      , planet density 
  • with faint companion star, density may decrease by ~ 3x

~ 0.5” companions of Kepler/K2 planet candidate hosts and the exoplanet radius distribution.

  • Fulton gap is robust regarding undetected stellar companions
  • gap became broader & shallower when accounting for possible undetected stellar companions
  • core composition of super-Earth and sub-Neptune exoplanets may not have so strong a divide as is suggested initially
  • w/o high-resolution imaging of Kepler and TESS host stars, the exoplanet radius distribution will be incorrectly inferred.

Fulton

mini-neptunes

super-Earths

  • Large radius errors originally hid distribution features
  • Fulton gap revealed after CKS (10% stellar radius errors)
  • Accounting for binarity shifts gap in the distribution

brightest companion 1''

brightest companion 2''

Shift from 1.8 to 2.2 

  • increased water/ice vs pure Si rock
R_E
  • 170 KOI companions < 2”
    • AO, speckle, lucky, or HST
  • constrained stellar properties
  • assessed the probability that the companions are physically bound
  • 60 - 80% of companions < 1”  
  • > 90% of companions < 0.5” were found to be bound
  • assuming the planet is equally likely to be orbiting the primary or secondary, unless they are vetted, nearly half of all Kepler planets may have radii underestimated by an average of 65%.
  • DSSI @ G-S
  • highest-resolution images to date
  • 27 mas
  • imaging from 0.32 to 14.5 au
  • excludes all possible stellar and brown dwarf companions
  • Detection % vs. projected angular separation for TRAPPIST-1.
  • 5σ detection limit convolved with the detection likelihood.
  • eliminates all companions in the blue region
    • separations of 0.32–17 au
  • Numbers are the dmag & spectral type limits at the corresponding points.
    • ex. all companions earlier than T7 are eliminated < ∼8 au.

observed

bound (simulated)

line-of-sight (simulated)

K2 binaries detectable with DSSI-

most companions within 1'' are bound, while only ~50% within 2'' are bound

  • Speckle imaging can detect potential companions to the left of the dashed lines
    • eliminating companions that may cause false positives or dilute transits.
  • Solid gray lines show the contrast ratios at which planet radii will have 1% and 0.1% correction factors due to the presence of a stellar companion.

non-speckle uses - time series photometry

  • binary WD+WD eclipsing system 
  • Usually one is a DA, 1.2M-sun, earth size, the other is a HE WD, 0.2 M-sun, much larger in radius.
  • 320 frames sequence, each at 3 sec
  • period < 1hr
  • EMCCD sensitivity + speed means sub-second time series are possible

Future

  • Quad-channel
  • Wave-front sensing
  • Possible near-IR channel?
  • DSSI heritage at the DCT+NPOI
  • novel construction

QWSSI

PI: van Belle

NEAs

  • sizes
    • radar typically quotes 40% uncertainties
    • albedo/radar size mismatch
  • shapes
    • adapt stellar surface modeling tools
    • light curve inversion/illumination model 

a~2.8AU

d~270km x 80km

(neck~50-65km)

model: Franck Marchis

Phaethon

Dec 2017 ~ 0.07AU

d~6km

Point source PS

Phaethon power spectrum (resolved)

Seeing-limited

Reconstructed

42''

Two-color wide-field speckle reconstruction

from NESSI

  • 0.25'' resolution from 500 frames (20s)
  • compromise b/t angres and contrast
  • Seeing ~ 0.85''

3-4x improvement over native seeing

Text

Text

M13

  • astrometry of clusters
  • very early results and little calibration or modelling over the FoV

Reconstructed

Seeing-limited

Abs astrometry residuals

~6-7 mas

4.2 mas

Further improvement possible to obtain ~0.02 pix (~0.4mas)

  • optics, dither, deeper obs

0.05 mag accuracy on aperture photometry

FWHM ~7pix (0.57'')

FWHM ~4pix (0.3'')

Pluto opposition

SOFIA occultation target duplicity

WF and speckle optics

Conv and EMCCD imaging

466/832 and g/i filters

TNO Varda ⌀~700km, d~50au

multiplicity is the single largest source of error to occultation path prediction

next steps: grisms & h-alpha filters

  • observe jets of Wolf-Rayet stars in h-alpha compared to "continuum" SDSS

HST WR 124

  • leverage the 0.01"/pixel plate scale of the instruments
  • extremely sensitive EMCCD, large aperture (Gemini), low-dispersion grism 
  • sci case1: high-z transients, GRBs 
  • sci case2: exoplanet host flares

Open to the community

  • NOAO proposal process
  • NESSI - WIYN@KPNO
  • `Alopeke - Gemini-N
  • Zorro - Gemini-S
  • Transitioning from block queue operation to pseudo-facility instrument

Summary


  • a brief history and introduction to speckle interferometry
  • the relevance of high angular resolution imaging to exoplanet and stellar research
  • how EMCCDs have revitalized speckle techniques
  • current speckle projects and science
  • development of future community instruments
  • other benefits and uses of EMCCDs with speckle imagers

  • small space telescopes, optimized for large surveys of the sky, like TESS and GAIA, perform relatively poorly at high-contrast, sub-arcsecond resolution where speckle imaging excels

  • a cost-effective, efficient means to reach the diffraction limit with single aperture telescopes

Additional References 

Brown, R. H. and Twiss, R. Q., Nature 177, 27{29 (Jan. 1956).
Hanbury Brown, R., Davis, J., and Allen, L. R., MNRAS 137, 375 (1967).
Labeyrie, A., A&A 6, 85 (May 1970).
Knox, K. T. and Thompson, B. J., ApJ 193, L45{L48 (Oct. 1974).
Weigelt, G., Scientifc Importance of High Angular Resolution at Infrared and Optical Wavelengths, Ulrich, M. H. and Kjaer, K., eds., 95{114 (1981).
Weigelt, G., Lowell Observatory Bulletin 9, 144{152 (1983).
Weigelt, G. P., Optics Communications 21, 55{59 (Apr. 1977).
Lohmann, A. W., Weigelt, G., and Wirnitzer, B., Appl. Opt. 22, 4028{4037 (Dec. 1983).

Fried, D. L., Journal of the Optical Society of America (1917-1983) 68, 1651–1658 (Dec. 1978).
Horch, E. P., Veillette, D. R., Baena Galle, R., Shah, S. C.,O'Rielly, G. V., and van Altena, W. F., AJ 137, 5057{5067 (June 2009).

Brugière, T., Mayer, F., Fereyre, P., Gu´erin, C., Dominjon, A., and Barbier, R.,  Nuclear Instruments and Methods in Physics Research A 787, 336–339 (July 2015).
Scott, N. J., Howell, S. B., Horch, E. P., and Everett, M. E., PASP 130, 054502 (May 2018).
Hofmann, K.-H. and Weigelt, G., A&A 278, 328-339 (Oct. 1993).
Law, N. M., Mackay, C. D., and Baldwin, J. E.,  A&A 446, 739{745 (Feb. 2006).
Tokovinin, A. and Cantarutti, R., PASP 120, 170 (Feb. 2008).
Andor, CCD, EMCCD and ICCDComparions & Minimizing Clock Induced Charge.
Howell, S. B., Everett, M. E., Sherry, W., Horch, E., and Ciardi, D. R., AJ 142, 19 (July 2011).

Furlan, E., Ciardi, D. R., Everett, M. E., Saylors, M., Teske, J. K., Horch, E. P., Howell, S. B., van Belle, G. T., Hirsch, L. A., Gautier, III, T. N., Adams, E. R., Barrado, D., Cartier, K. M. S., Dressing C. D., Dupree, A. K., Gilliland, R. L., Lillo-Box, J., Lucas, P. W., and Wang, J.,  AJ 153, 71 (Feb. 2017).
Furlan, E. and Howell, S. B., AJ 154, 66 (Aug. 2017).
Teske, J. K., Ciardi, D. R., Howell, S. B., Hirsch, L. A., and Johnson, R. A.,  ArXiv e-prints (Apr. 2018).
Fulton, B. J., Petigura, E. A., Howard, A. W., Isaacson, H., Marcy, G. W., Cargile, P. A., Hebb, L., Weiss, L. M., Johnson, J. A., Morton, T. D.,Sinuko , E., Cross eld, I. J. M., and Hirsch, L. A., AJ 154, 109 (Sept. 2017).
Hirsch, L. A., Ciardi, D. R., Howard, A. W., Everett, M. E., Furlan, E., Saylors, M., Horch, E. P., Howell, S. B., Teske, J., and Marcy, G. W., AJ 153, 117 (Mar. 2017).

Matson, R. A., Howell, S. B., Horch, E. P., Everett, M. E., ArXiv e-prints (May. 2018)