Context. Island effect (IE) aberrations are induced by differential pistons, tips, and tilts between neighboring pupil segments on ground-based telescopes, which severely limit the observations of circumstellar environments on the recently deployed exoplanet imagers (e.g., VLT/SPHERE, Gemini/GPI, Subaru/SCExAO) during the best observing conditions. Caused by air temperature gradients at the level of the telescope spiders, these aberrations were recently diagnosed with success on VLT/SPHERE, but so far no complete calibration has been performed to overcome this issue. Aims. We propose closed-loop focal plane wavefront control based on the asymmetric Fourier pupil wavefront sensor (APF-WFS) to calibrate these aberrations and improve the image quality of exoplanet high-contrast instruments in the presence of the IE. Methods. Assuming the archetypal four-quadrant aperture geometry in 8 m class telescopes, we describe these aberrations as a sum of the independent modes of piston, tip, and tilt that are distributed in each quadrant of the telescope pupil. We calibrate these modes with the APF-WFS before introducing our wavefront control for closed-loop operation. We perform numerical simulations and then experimental tests on a real system using Subaru/SCExAO to validate our control loop in the laboratory and on-sky. Results. Closed-loop operation with the APF-WFS enables the compensation for the IE in simulations and in the laboratory for the small aberration regime. Based on a calibration in the near infrared, we observe an improvement of the image quality in the visible range on the SCExAO/VAMPIRES module with a relative increase in the image Strehl ratio of 37%. Conclusions. Our first IE calibration paves the way for maximizing the science operations of the current exoplanet imagers. Such an approach and its results prove also very promising in light of the Extremely Large Telescopes (ELTs) and the presence of similar artifacts with their complex aperture geometry.

We present the Vortex Image Processing (VIP) library, a python package dedicated to astronomical high-contrast imaging. Our package relies on the extensive python stack of scientific libraries and aims to provide a flexible framework for high-contrast data and image processing. In this paper, we describe the capabilities of VIP related to processing image sequences acquired using the angular differential imaging (ADI) observing technique. VIP implements functionalities for building high-contrast data processing pipelines, encompassing pre- and post-processing algorithms, potential source. position and flux estimation, and sensitivity curve. generation. Among the reference point-spread. function subtraction techniques for ADI post-processing, VIP includes several flavors of principal component analysis (PCA) based algorithms, such as annular PCA and incremental PCA algorithms capable of processing big datacubes (of several gigabytes) on a computer with limited memory. Also, we present a novel ADI algorithm based on non-negative matrix factorization, which comes from the same family of low-rank matrix approximations as PCA and provides fairly similar results. We showcase the ADI capabilities of the VIP library using a deep sequence on HR 8799 taken with the LBTI/LMIRCam and its recently commissioned L-band vortex coronagraph. Using VIP, we investigated the presence of additional companions around HR 8799 and did not find any significant additional point source beyond the four known planets. VIP is available at http://github. com/vortex-exoplanet/VIP and is accompanied with Jupyter notebook tutorials illustrating the main functionalities of the library.

An investigation is presented of the dust scattering in the coma of 67P/Churyumov-Gerasimenko for the dates of 2015 February 28, March 15 and April 27. A comparison of the morphology of dust continuum maps at 1.1 mu m and gas emission shows that for the above dates the spatial distribution of the dust is strongly correlated with H2O but not with CO2 emission. For April 27, the radial profile on the illuminated side of the nucleus in the inner coma agrees well with the direct simulation Monte Carlo (DSMC) calculations as the dust is accelerating and flows outwards distribution of the dust is narrower than the broader emission of the gas. Af. values are 1.13 m for 2015 February 28, 2.02 m for April 27, while local values for March 15 are 2.3-5.3 m, depending on the nucleus illumination. In the inner coma, the spectral reflectivity from 0.35 to 3.5 mu m displays a red slope with a change at around 1 mu m. From 0.35 to 0.8 mu m, the values range from 9 to 12 +/- 1 per cent per 100 nm both on the sunlit side and on the dark side. From 1 to 2.5 mu m, the values are 1.7 +/- 0.2 per cent per 100 nm on the sunlit side and 3 +/- 1 per cent per 100 nm on the dark side. For the August 26 jet, no significant increase of the colour gradient with distance from the nucleus could be observed, nor any significant difference detected between the jet and the background coma.

The technique of non-redundant masking (NRM) transforms a conventional telescope into an interferometric array. In practice, this provides a much better constrained point-spread function than a filled aperture and thus higher resolution than traditional imaging methods. Here, we describe an NRM data reduction pipeline. We discuss strategies for NRM observations regarding dithering patterns and calibrator selection. We describe relevant image calibrations and use example Large Binocular Telescope data sets to show their effects on the scatter in the Fourier measurements. We also describe the various ways to calculate Fourier quantities, and discuss different calibration strategies. We present the results of image reconstructions from simulated observations where we adjust prior images, weighting schemes, and error bar estimation. We compare two imaging algorithms and discuss implications for reconstructing images from real observations. Finally, we explore how the current state of the art compares to next-generation Extremely Large Telescopes.

We present an empirical assessment of the use of broad-band optical colours as age indicators for unresolved extragalactic clusters and investigate stochastic sampling effects on integrated colours. We use the integrated properties of Galactic open clusters (OCs) as models for unresolved extragalactic clusters. The population synthesis code Starburst99 (SB99) and four optical colours were used to estimate how well we can recover the ages of 62 well-studied Galactic OCs with published ages. We provide a method for estimating the ages of unresolved clusters and for reliably determining the uncertainties in the age estimates. Our results support earlier conclusions based on comparisons to synthetic clusters, namely the (U - B) colour is critical to the estimation of the ages of star-forming regions. We compare the observed optical colours with those obtained from SB99 using the published ages and get good agreement. The scatter in the (B - V)observed - (B - V)model is larger for lower luminosity clusters, perhaps due to stochastic effects.

Context. Since its rendezvous with comet 67P/Churyumov-Gerasimenko (67P), the Rosetta spacecraft has provided invaluable information contributing to our understanding of the cometary environment. On board, the VIRTIS and ROSINA instruments can both measure gas parameters in the rarefied cometary atmosphere, the so-called coma, and provide complementary results with remote sensing and in situ measurement techniques, respectively. The data from both ROSINA and VIRTIS instruments suggest that the source regions of H2O and CO2 are not uniformly distributed over the surface of the nucleus even after accounting for the changing solar illumination of the irregularly shaped rotating nucleus. The source regions of H2O and CO2 are also relatively different from one another. Aims. The use of a combination of a formal numerical data inversion method with a fully kinetic coma model is a way to correlate and interpret the information provided by these two instruments to fully understand the volatile environment and activity of comet 67P. Methods. In this work, the nonuniformity of the outgassing activity at the surface of the nucleus is described by spherical harmonics and constrained by ROSINA-DFMS data. This activity distribution is coupled with the local illumination to describe the inner boundary conditions of a 3D direct simulation Monte-Carlo (DSMC) approach using the Adaptive Mesh Particle Simulator (AMPS) code applied to the H2O and CO2 coma of comet 67P. Results. We obtain activity distribution of H2O and CO2 showing a dominant source of H2O in the Hapi region, while more CO2 is produced in the southern hemisphere. The resulting model outputs are analyzed and compared with VIRTIS-M/-H and ROSINADFMS measurements, showing much better agreement between model and data than a simpler model assuming a uniform surface activity. The evolution of the H2O and CO2 production rates with heliocentric distance are derived accurately from the coma model showing agreement between the observations from the different instruments and ground-based observations. Conclusions. We derive the activity distributions for H2O and CO2 at the surface of the nucleus described in spherical harmonics, which we couple to the local solar illumination to constitute the boundary conditions of our coma model. The model presented reproduces the coma observations made by the ROSINA and VIRTIS instruments on board the Rosetta spacecraft showing our understanding of the physics of 67P’s coma. This model can be used for further data analyses, such as dust modeling, in a future work.

No description available.

Context. Radial velocity and transit methods are effective for the study of short orbital period exoplanets but they hardly probe objects at large separations for which direct imaging can be used. Aims. We carried out the international deep planet survey of 292 young nearby stars to search for giant exoplanets and determine their frequency. Methods. We developed a pipeline for a uniform processing of all the data that we have recorded with NIRC2/Keck II, NIRI/Gemini North, NICI/Gemini South, and NACO/VLT for 14 yr. The pipeline first applies cosmetic corrections and then reduces the speckle intensity to enhance the contrast in the images. Results. The main result of the international deep planet survey is the discovery of the HR8799 exoplanets. We also detected 59 visual multiple systems including 16 new binary stars and 2 new triple stellar systems, as well as 2279 point-like sources. We used Monte Carlo simulations and the Bayesian theorem to determine that 1.05(-0.70)(+2.80)% of stars harbor at least one giant planet between 0.5 and 14 MJ and between 20 and 300AU. This result is obtained assuming uniform distributions of planet masses and semi-major axes. If we consider power law distributions as measured for close-in planets instead, the derived frequency is 2.30(-1.55)(+5.95)%, recalling the strong impact of assumptions on Monte Carlo output distributions. We also find no evidence that the derived frequency depends on the mass of the hosting star, whereas it does for close-in planets. Conclusions. The international deep planet survey provides a database of confirmed background sources that may be useful for other exoplanet direct imaging surveys. It also puts new constraints on the number of stars with at least one giant planet reducing by a factor of two the frequencies derived by almost all previous works.

We present new 3.6 and 4.5 mu m Spitzer phase curves for the highly irradiated hot Jupiter WASP-33b and the unusually dense Saturn-mass planet HD 149026b. As part of this analysis, we develop a new variant of pixel-level decorrelation that is effective at removing intrapixel sensitivity variations for long observations (>10 hr) where the position of the star can vary by a significant fraction of a pixel. Using this algorithm, we measure eclipse depths, phase amplitudes, and phase offsets for both planets at 3.6 and 4.5 mu m. We use a simple toy model to show that WASP-33b's phase offset, albedo, and heat recirculation efficiency are largely similar to those of other hot Jupiters despite its very high irradiation. On the other hand, our fits for HD 149026b prefer a very high albedo. We also compare our results to predictions from general circulation models, and we find that while neither planet matches the models well, the discrepancies for HD 149026b are especially large. We speculate that this may be related to its high bulk metallicity, which could lead to enhanced atmospheric opacities and the formation of reflective cloud layers in localized regions of the atmosphere. We then place these two planets in a broader context by exploring relationships between the temperatures, albedos, heat transport efficiencies, and phase offsets of all planets with published thermal phase curves. We find a striking relationship between phase offset and irradiation temperature: the former drops with increasing temperature until around 3400 K and rises thereafter. Although some aspects of this trend are mirrored in the circulation models, there are notable differences that provide important clues for future modeling efforts.

yes / A comprehensive spectral analysis of both the solar background magnetic field (SBMF) in cycles 21-23 and the sunspot magnetic field in cycle 23 reported in our recent paper showed the presence of two principal components (PCs) of SBMF having opposite polarity, e. g., originating in the northern and southern hemispheres, respectively. Over a duration of one solar cycle, both waves are found to travel with an increasing phase shift toward the northern hemisphere in odd cycles 21 and 23 and to the southern hemisphere in even cycle 22. These waves were linked to solar dynamo waves assumed to form in different layers of the solar interior. In this paper, for the first time, the PCs of SBMF in cycles 21-23 are analyzed with the symbolic regression technique using Hamiltonian principles, allowing us to uncover the underlying mathematical laws governing these complex waves in the SBMF presented by PCs and to extrapolate these PCs to cycles 24-26. The PCs predicted for cycle 24 very closely fit (with an accuracy better than 98%) the PCs derived from the SBMF observations in this cycle. This approach also predicts a strong reduction of the SBMF in cycles 25 and 26 and, thus, a reduction of the resulting solar activity. This decrease is accompanied by an increasing phase shift between the two predicted PCs (magnetic waves) in cycle 25 leading to their full separation into the opposite hemispheres in cycle 26. The variations of the modulus summary of the two PCs in SBMF reveals a remarkable resemblance to the average number of sunspots in cycles 21-24 and to predictions of reduced sunspot numbers compared to cycle 24: 80% in cycle 25 and 40% in cycle 26.