ABSTRACT A new class of high-contrast image analysis algorithms that empirically fit and subtract systematic noise has lead to recent discoveries of faint exoplanet/substellar companions and ...scattered light images of circumstellar disks. These methods are extremely efficient at enhancing the detectability of a faint astrophysical signal, but they generally create systematic biases in their observed properties. This paper provides a general solution for this outstanding problem. We present an analytical derivation of a linear expansion that captures the impact of astrophysical over-subtraction or self-subtraction in current image analysis techniques. We examine the general case for which the reference images of the astrophysical scene move azimuthally and/or radially across the field of view as a result of the observation strategy. Our new method is based on perturbing the covariance matrix underlying any least-squares speckles problem, and propagating this perturbation through the data analysis algorithm. Most of the work in this paper is presented in the Principal Component Analysis framework, but it can be easily generalized to methods relying on the linear combination of images (instead of eigenmodes). Based on this linear expansion, which is obtained in the most general case, we then demonstrate practical applications of this new algorithm. We first consider the spectral extraction of faint point sources in IFS data and illustrate, using public Gemini Planet Imager commissioning data, that our novel perturbation-based Forward Modeling, which we named Karhunen Loeve Image Processing (KLIP-FM), can indeed alleviate algorithmic biases. We then apply KLIP-FM to the detection of point sources and show how it decreases the rate of false negatives while keeping the rate of false positives unchanged when compared to classical KLIP. This can potentially have important consequences on the design of follow-up strategies of ongoing direct imaging surveys.
We apply the vectorized non-negative matrix factorization (NMF) method to the post-processing of the direct imaging data of exoplanetary systems such as circumstellar disks. NMF is an iterative ...approach, which first creates a nonorthogonal and non-negative basis of components using the given reference images and then models a target with the components. The constructed model is then rescaled with a factor to compensate for the contribution from the disks. We compare NMF with existing methods (classical reference differential imaging method, and the Karhunen-Loève image projection algorithm) using synthetic circumstellar disks and demonstrate the superiority of NMF: with no need of prior selection of references, NMF not only can detect fainter circumstellar disks but also better preserves their morphology and does not require forward modeling. As an application to a well-known disk example, we process the archival Hubble Space Telescope STIS coronagraphic observations of HD 181327 with different methods and compare them, and NMF is able to extract some circumstellar materials inside the primary ring for the first time. In an appendix, we mathematically investigate the stability of NMF components during the iteration and the linearity of NMF modeling.
The Apodized Pupil Lyot Coronagraph (APLC) is a diffraction suppression system installed in the recently deployed instruments Palomar/P1640, Gemini/GPI, and VLT/SPHERE to allow direct imaging and ...spectroscopy of circumstellar environments. Using a prolate apodization, the current implementations offer raw contrasts down to 10 super(-7) at 0.2 arcsec from a star over a wide bandpass (20%), in the presence of central obstruction and struts, enabling the study of young or massive gaseous planets. Observations of older or lighter companions at smaller separations would require improvements in terms of the inner working angle (IWA) and contrast, but the methods originally used for these designs were not able to fully explore the parameter space. We propose a novel approach to improve the APLC performance. Our method relies on the linear properties of the coronagraphic electric field with the apodization at any wavelength to develop numerical solutions producing coronagraphic star images with high-contrast region in broadband light. We explore the parameter space by considering different aperture geometries, contrast levels, dark-zone sizes, bandpasses, and focal plane mask sizes. We present an application of these solutions to the case of Gemini/GPI with a design delivering a 10 super(-8) raw contrast at 0.19 arcsec and offering a significantly reduced sensitivity to low-order aberrations compared to the current implementation. Optimal solutions have also been found to reach 10 super(-10) contrast in broadband light regardless of the aperture shape, with effective IWA in the 2-3.5 lambda/D range, therefore making the APLC a suitable option for the future exoplanet direct imagers on the ground or in space.
We present a statistical analysis of the first 300 stars observed by the Gemini Planet Imager Exoplanet Survey. This subsample includes six detected planets and three brown dwarfs; from these ...detections and our contrast curves we infer the underlying distributions of substellar companions with respect to their mass, semimajor axis, and host stellar mass. We uncover a strong correlation between planet occurrence rate and host star mass, with stars M* > 1.5 M more likely to host planets with masses between 2 and 13MJup and semimajor axes of 3-100 au at 99.92% confidence. We fit a double power-law model in planet mass (m) and semimajor axis (a) for planet populations around high-mass stars (M* > 1.5 M ) of the form , finding = −2.4 0.8 and β = −2.0 0.5, and an integrated occurrence rate of % between 5-13MJup and 10-100 au. A significantly lower occurrence rate is obtained for brown dwarfs around all stars, with % of stars hosting a brown dwarf companion between 13-80MJup and 10-100 au. Brown dwarfs also appear to be distributed differently in mass and semimajor axis compared to giant planets; whereas giant planets follow a bottom-heavy mass distribution and favor smaller semimajor axes, brown dwarfs exhibit just the opposite behaviors. Comparing to studies of short-period giant planets from the radial velocity method, our results are consistent with a peak in occurrence of giant planets between ∼1 and 10 au. We discuss how these trends, including the preference of giant planets for high-mass host stars, point to formation of giant planets by core/pebble accretion, and formation of brown dwarfs by gravitational instability.
Abstract We have obtained new images of the protoplanetary disk orbiting TW Hya in visible, total intensity light with the Space Telescope Imaging Spectrograph (STIS) on the Hubble Space Telescope ( ...HST ), using the newly commissioned BAR5 occulter. These HST /STIS observations achieved an inner working angle of ∼0.″2, or 11.7 au, probing the system at angular radii coincident with recent images of the disk obtained by ALMA and in polarized intensity near-infrared light. By comparing our new STIS images to those taken with STIS in 2000 and with NICMOS in 1998, 2004, and 2005, we demonstrate that TW Hya’s azimuthal surface brightness asymmetry moves coherently in position angle. Between 50 au and 141 au we measure a constant angular velocity in the azimuthal brightness asymmetry of 22.°7 yr −1 in a counterclockwise direction, equivalent to a period of 15.9 yr assuming circular motion. Both the (short) inferred period and lack of radial dependence of the moving shadow pattern are inconsistent with Keplerian rotation at these disk radii. We hypothesize that the asymmetry arises from the fact that the disk interior to 1 au is inclined and precessing owing to a planetary companion, thus partially shadowing the outer disk. Further monitoring of this and other shadows on protoplanetary disks potentially opens a new avenue for indirectly observing the sites of planet formation.
First light of the Gemini Planet Imager Macintosh, Bruce; Graham, James R.; Ingraham, Patrick ...
Proceedings of the National Academy of Sciences - PNAS,
09/2014, Letnik:
111, Številka:
35
Journal Article
Recenzirano
Odprti dostop
The Gemini Planet Imager is a dedicated facility for directly imaging and spectroscopically characterizing extrasolar planets. It combines a very high-order adaptive optics system, a ...diffraction-suppressing coronagraph, and an integral field spectrograph with low spectral resolution but high spatial resolution. Every aspect of the Gemini Planet Imager has been tuned for maximum sensitivity to faint planets near bright stars. During first-light observations, we achieved an estimated H band Strehl ratio of 0.89 and a 5-σ contrast of 10 ⁶ at 0.75 arcseconds and 10 ⁵ at 0.35 arcseconds. Observations of Beta Pictoris clearly detect the planet, Beta Pictoris b, in a single 60-s exposure with minimal postprocessing. Beta Pictoris b is observed at a separation of 434 ± 6 milliarcseconds (mas) and position angle 211.8 ± 0.5°. Fitting the Keplerian orbit of Beta Pic b using the new position together with previous astrometry gives a factor of 3 improvement in most parameters over previous solutions. The planet orbits at a semimajor axis of Formula near the 3:2 resonance with the previously known 6-AU asteroidal belt and is aligned with the inner warped disk. The observations give a 4% probability of a transit of the planet in late 2017.
We report the results of a ∼4 yr direct imaging survey of 104 stars to resolve and characterize circumstellar debris disks in scattered light as part of the Gemini Planet Imager (GPI) Exoplanet ...Survey. We targeted nearby ( 150 pc), young ( 500 Myr) stars with high infrared (IR) excesses (LIR/L > 10−5), including 38 with previously resolved disks. Observations were made using the GPI high-contrast integral field spectrograph in H-band (1.6 m) coronagraphic polarimetry mode to measure both polarized and total intensities. We resolved 26 debris disks and 3 protoplanetary/transitional disks. Seven debris disks were resolved in scattered light for the first time, including newly presented HD 117214 and HD 156623, and we quantified basic morphologies of five of them using radiative transfer models. All of our detected debris disks except HD 156623 have dust-poor inner holes, and their scattered-light radii are generally larger than corresponding radii measured from resolved thermal emission and those inferred from spectral energy distributions. To assess sensitivity, we report contrasts and consider causes of nondetections. Detections were strongly correlated with high IR excess and high inclination, although polarimetry outperformed total intensity angular differential imaging for detecting low-inclination disks ( 70°). Based on postsurvey statistics, we improved upon our presurvey target prioritization metric predicting polarimetric disk detectability. We also examined scattered-light disks in the contexts of gas, far-IR, and millimeter detections. Comparing H-band and ALMA fluxes for two disks revealed tentative evidence for differing grain properties. Finally, we found no preference for debris disks to be detected in scattered light if wide-separation substellar companions were present.
We describe a Bayesian rejection-sampling algorithm designed to efficiently compute posterior distributions of orbital elements for data covering short fractions of long-period exoplanet orbits. Our ...implementation of this method, Orbits for the Impatient (OFTI), converges up to several orders of magnitude faster than two implementations of Markov Chain Monte Carlo (MCMC) in this regime. We illustrate the efficiency of our approach by showing that OFTI calculates accurate posteriors for all existing astrometry of the exoplanet 51 Eri b up to 100 times faster than a Metropolis-Hastings MCMC. We demonstrate the accuracy of OFTI by comparing our results for several orbiting systems with those of various MCMC implementations, finding the output posteriors to be identical within shot noise. We also describe how our algorithm was used to successfully predict the location of 51 Eri b six months in the future based on less than three months of astrometry. Finally, we apply OFTI to 10 long-period exoplanets and brown dwarfs, all but one of which have been monitored over less than 3% of their orbits, producing fits to their orbits from astrometric records in the literature.
Gravitational deflection of starlight around the Sun during the 1919 total solar eclipse provided measurements that confirmed Einstein’s general theory of relativity. We have used the Hubble Space ...Telescope to measure the analogous process of astrometric microlensing caused by a nearby star, the white dwarf Stein 2051 B. As Stein 2051 B passed closely in front of a background star, the background star’s position was deflected. Measurement of this deflection at multiple epochs allowed us to determine the mass of Stein 2051 B—the sixth-nearest white dwarf to the Sun—as 0.675 ± 0.051 solar masses. This mass determination provides confirmation of the physics of degenerate matter and lends support to white dwarf evolutionary theory.
Abstract We have performed a comprehensive study of the Orion Nebula Cluster (ONC) combining the photometric data obtained by the two Hubble Space Telescope Treasury programs that targeted this ...region. To consistently analyze the rich data set obtained in a wide variety of filters, we adopted a Bayesian approach to fit the spectral energy distribution of the sources, deriving mass, age, extinction, distance, and accretion for each source in the region. The three-dimensional study of mass distribution for bona fide cluster members shows that mass segregation in the ONC extends to subsolar masses, while the age distribution strongly supports the idea that star formation in the ONC is best described by a major episode of star formation that happened ∼1 Myr ago. For masses ≳0.1 M ⊙ , our derived empirical initial mass function (IMF) is in good agreement with a Chabrier system IMF. Both the accretion luminosity ( L acc ) and mass accretion rates ( M ̇ acc ) are best described by broken power-law relations. This suggests that for the majority of young circumstellar disks in this cluster the excess emission may be dominated by X-ray-driven photoevaporation by the central star rather than external photoevaporation. If this is the case, the slopes of the power-law relations may be largely determined by the initial conditions set at the onset of the star formation process, which may be quite similar between regions that eventually form clusters of different sizes.