SPHERE (Beuzit et al,. 2019) has now been in operation at the VLT for more than 5 years, demonstrating a high level of performance. SPHERE has produced outstanding results using a variety of operating modes, primarily in the field of direct imaging of exoplanetary systems, focusing on exoplanets as point sources and circumstellar disks as extended objects. The achievements obtained thus far with SPHERE (~200 refereed publications) in different areas (exoplanets, disks, solar system, stellar physics...) have motivated a large consortium to propose an even more ambitious set of science cases, and its corresponding technical implementation in the form of an upgrade. The SPHERE+ project capitalizes on the expertise and lessons learned from SPHERE to push high contrast imaging performance to its limits on the VLT 8m-telescope. The scientific program of SPHERE+ described in this document will open a new and compelling scientific window for the upcoming decade in strong synergy with ground-based facilities (VLT/I, ELT, ALMA, and SKA) and space missions (Gaia, JWST, PLATO and WFIRST). While SPHERE has sampled the outer parts of planetary systems beyond a few tens of AU, SPHERE+ will dig into the inner regions around stars to reveal and characterize by mean of spectroscopy the giant planet population down to the snow line. Building on SPHERE's scientific heritage and resounding success, SPHERE+ will be a dedicated survey instrument which will strengthen the leadership of ESO and the European community in the very competitive field of direct imaging of exoplanetary systems. With enhanced capabilities, it will enable an even broader diversity of science cases including the study of the solar system, the birth and death of stars and the exploration of the inner regions of active galactic nuclei.

In this work we explore the possibility of using Recurrence Quantification Analysis (RQA) in astronomical high-contrast imaging to statistically discriminate the signal of faint objects from speckle noise. To this end, we tested RQA on a sequence of high frame rate (1 kHz) images acquired with the SHARK-VIS forerunner at the Large Binocular Telescope. Our tests show promising results in terms of detection contrasts at angular separations as small as $50$ mas, especially when RQA is applied to a very short sequence of data ($2$ s). These results are discussed in light of possible science applications and with respect to other techniques like, for example, Angular Differential Imaging and Speckle-Free Imaging.

Residual speckles in adaptive optics (AO) images represent a well-known limitation on the achievement of the contrast needed for faint source detection. Speckles in AO imagery can be the result of either residual atmospheric aberrations, not corrected by the AO, or slowly evolving aberrations induced by the optical system. We take advantage of the high temporal cadence (1 ms) of the data acquired by the System for Coronagraphy with High-order Adaptive Optics from R to K bands-VIS forerunner experiment at the Large Binocular Telescope to characterize the AO residual speckles at visible wavelengths. An accurate knowledge of the speckle pattern and its dynamics is of paramount importance for the application of methods aimed at their mitigation. By means of both an automatic identification software and information theory, we study the main statistical properties of AO residuals and their dynamics. We therefore provide a speckle characterization that can be incorporated into numerical simulations to increase their realism and to optimize the performances of both real-time and postprocessing techniques aimed at the reduction of the speckle noise.

A robust post processing technique is mandatory to analyse the coronagraphic high contrast imaging data. Angular Differential Imaging (ADI) and Principal Component Analysis (PCA) are the most used approaches to suppress the quasi-static structure in the Point Spread Function (PSF) in order to revealing planets at different separations from the host star. The focus of this work is to apply these two data reduction techniques to obtain the best limit detection for each coronagraphic setting that has been simulated for the SHARK-NIR, a coronagraphic camera that will be implemented at the Large Binocular Telescope (LBT). We investigated different seeing conditions ($0.4"-1"$) for stellar magnitude ranging from R=6 to R=14, with particular care in finding the best compromise between quasi-static speckle subtraction and planet detection.

In this article, we present a simulator conceived for the conceptual study of an AO-fed high-contrast coronagraphic imager. The simulator implements physical optics: a complex disturbance (the electric field) is Fresnel-propagated through any user-defined optical train, in an end-to-end fashion. The effect of atmospheric residual aberrations and their evolution with time can be reproduced by introducing in input a temporal sequence of phase screens: synthetic images are then generated by co-adding instantaneous PSFs. This allows studying with high accuracy the impact of AO correction on image quality for different integration times and observing conditions. In addition, by conveniently detailing the optical model, the user can easily implement any coronagraphic set-up and introduce optical aberrations at any position. Furthermore, generating multiple images can allow exploring detection limits after a differential post-processing algorithm is applied (e.g. Angular Differential Imaging). The simulator has been developed in the framework of the design of SHARK-NIR, the second-generation high contrast imager selected for the Large Binocular Telescope.

SHARK-NIR is one of the two coronagraphic instruments proposed for the Large Binocular Telescope. Together with SHARK-VIS (performing coronagraphic imaging in the visible domain), it will offer the possibility to do binocular observations combining direct imaging, coronagraphic imaging and coronagraphic low resolution spectroscopy in a wide wavelength domain, going from 0.5{\mu}m to 1.7{\mu}m. Additionally, the contemporary usage of LMIRCam, the coronagraphic LBTI NIR camera, working from K to L band, will extend even more the covered wavelength range. In January 2017 SHARK-NIR underwent a successful final design review, which endorsed the instrument for construction and future implementation at LBT. We report here the final design of the instrument, which foresees two intermediate pupil planes and three focal planes to accomodate a certain number of coronagraphic techniques, selected to maximize the instrument contrast at various distances from the star. Exo-Planets search and characterization has been the science case driving the instrument design, but the SOUL upgrade of the LBT AO will increase the instrument performance in the faint end regime, allowing to do galactic (jets and disks) and extra-galactic (AGN and QSO) science on a relatively wide sample of targets, normally not reachable in other similar facilities.

We present a new processing technique aimed at significantly improving the angular differential imaging method (ADI) in the context of high-contrast imaging of faint objects nearby bright stars in observations obtained with extreme adaptive optics (EXAO) systems. This technique, named "SFADI" for "Speckle-Free ADI", allows to improve the achievable contrast by means of speckles identification and suppression. This is possible in very high cadence data, which freeze the atmospheric evolution. Here we present simulations in which synthetic planets are injected into a real millisecond frame rate sequence, acquired at the LBT telescope at visible wavelength, and show that this technique can deliver low and uniform background, allowing unambiguous detection of $10^{-5}$ contrast planets, from $100$ to $300$ mas separations, under poor and highly variable seeing conditions ($0.8$ to $1.5$ arcsec FWHM) and in only $20$ min of acquisition. A comparison with a standard ADI approach shows that the contrast limit is improved by a factor of $5$. We extensively discuss the SFADI dependence on the various parameters like speckle identification threshold, frame integration time, and number of frames, as well as its ability to provide high-contrast imaging for extended sources, and also to work with fast acquisitions.

The reionization of the Universe is one of the most important topics of present day astrophysical research. The most plausible candidates for the reionization process are star-forming galaxies, which according to the predictions of the majority of the theoretical and semi-analytical models should dominate the HI ionizing background at z~3. We aim at measuring the Lyman continuum escape fraction, which is one of the key parameters to compute the contribution of star-forming galaxies to the UV background. We have used ultra-deep U-band imaging (U=30.2mag at 1sigma) by LBC/LBT in the CANDELS/GOODS-North field, as well as deep imaging in COSMOS and EGS fields, in order to estimate the Lyman continuum escape fraction of 69 star-forming galaxies with secure spectroscopic redshifts at 3.27<z<3.40 to faint magnitude limits (L=0.2L*, or equivalently M1500~-19). We have measured through stacks a stringent upper limit (<1.7% at 1sigma) for the relative escape fraction of HI ionizing photons from bright galaxies (L>L*), while for the faint population (L=0.2L*) the limit to the escape fraction is ~10%. We have computed the contribution of star-forming galaxies to the observed UV background at z~3 and we have found that it is not enough to keep the Universe ionized at these redshifts, unless their escape fraction increases significantly (>10%) at low luminosities (M1500>-19). We compare our results on the Lyman continuum escape fraction of high-z galaxies with recent estimates in the literature and discuss future prospects to shed light on the end of the Dark Ages. In the future, strong gravitational lensing will be fundamental to measure the Lyman continuum escape fraction down to faint magnitudes (M1500~-16) which are inaccessible with the present instrumentation on blank fields.

In February 2014, the SHARK-VIS (System for High contrast And coronography from R to K at VISual bands) Forerunner, a high contrast experimental imager operating at visible wavelengths, was installed at LBT (Large Binocular Telescope). Here we report on the first results obtained by recent on-sky tests. These results show the extremely good performance of the LBT ExAO (Extreme Adaptive Optics) system at visible wavelengths, both in terms of spatial resolution and contrast achieved. Similarly to what was done by (Amara et al. 2012), we used the SHARK-VIS Forerunner data to quantitatively assess the contrast enhancement. This is done by injecting several different synthetic faint objects in the acquired data and applying the ADI (angular differential imaging) technique. A contrast of the order of $5 \times 10^{-5}$ is obtained at 630 nm for angular separations from the star larger than 100 mas. These results are discussed in light of the future development of SHARK-VIS and compared to those obtained by other high contrast imagers operating at similar wavelengths.

The demonstration of efficient single-mode fiber (SMF) coupling is a key requirement for the development of a compact, ultra-precise radial velocity (RV) spectrograph. iLocater is a next generation instrument for the Large Binocular Telescope (LBT) that uses adaptive optics (AO) to inject starlight into a SMF. In preparation for commissioning iLocater, a prototype SMF injection system was installed and tested at the LBT in the Y-band (0.970-1.065 $\mu$m). This system was designed to verify the capability of the LBT AO system as well as characterize on-sky SMF coupling efficiencies. SMF coupling was measured on stars with variable airmasses, apparent magnitudes, and seeing conditions for six half-nights using the Large Binocular Telescope Interferometer. We present the overall optical and mechanical performance of the SMF injection system, including details of the installation and alignment procedure. A particular emphasis is placed on analyzing the instrument's performance as a function of telescope elevation to inform the final design of the fiber injection system for iLocater.

We present the results of a study along with a first prototype of a high precision system (? 1 arcsec) for pointing and tracking light (near-infrared) telescopes on board stratospheric balloons. Such a system is essentially composed by a star sensor and by a star tracker, able to recognize the field and to adequately track the telescope, respectively. We present the software aimed at processing the star sensor image and the predictive algorithm that allows the fine tracking of the source at a sub-pixel level. The laboratory tests of the system are described and its performance is analyzed. We demonstrate how such a device, when used at the focal plane of enough large telescopes (2-4m, F/10), is capable to provide (sub-)arcsec diffraction limited images in the near infrared bands.

Aims. The Large Binocular Cameras (LBC) are two twin wide field cameras (FOV ~ 23'x 25') mounted at the prime foci of the 8.4m Large Binocular Telescope (LBT). We performed a weak lensing analysis of the z=0.288 cluster Abell 611 on g-band data obtained by the blue-optimized Large Binocular Camera in order to estimate the cluster mass. Methods. Due to the complexity of the PSF of LBC, we decided to use two different approaches, KSB and Shapelets, to measure the shape of background galaxies and to derive the shear signal produced by the cluster. Then we estimated the cluster mass with both aperture densitometry and parametric model fits. Results. The combination of the large aperture of the telescope and the wide field of view allowed us to map a region well beyond the expected virial radius of the cluster and to get a high surface density of background galaxies (23 galaxies/arcmin^2). This made possible to estimate an accurate mass for Abell 611. We find that the mass within 1.5 Mpc is: $(8\pm3)\times 10^{14} M_\odot$ from the aperture mass technique and $(5\pm1)\times 10^{14} M_\odot$ using the model fitting by a NFW mass density profile, for both Shapelets and KSB methods. This analysis demonstrates that LBC is a powerful instrument for weak gravitational lensing studies.

Aims. With an observed and rest-frame duration of < 2s and < 0.5s, respectively, GRB090426 could be classified as a short GRB. The prompt detection, both from space and ground-based telescopes, of a bright optical counterpart to this GRB offered a unique opportunity to complete a detailed study. Methods. Based on an extensive ground-based observational campaign, we obtained the spectrum of the optical afterglow of GRB090426, measuring its redshift and obtaining information about the medium in which the event took place. We completed follow-up observation of the afterglow optical light curve down to the brightness level of the host galaxy that we firmly identified and studied. We also retrieved and analyzed all the available high-energy data of this event, and compared the results with our findings in the optical. This represents one of the most detailed studies of a short-duration event presented so far. Results. The time properties qualify GRB090426 as a short burst. In this case, its redshift of z = 2.61 would be the highest yet found for a GRB of this class. On the other hand, the spectral and energy properties are more similar to those of long bursts. LBT late-time deep imaging identifies a star-forming galaxy at a redshift consistent with that of the GRB. The afterglow lies within the light of its host and shows evidence of local absorption.

Deep multicolour surveys are the main tool to explore the formation and evolution of the faint galaxies which are beyond the spectroscopic limit with the present technology. The photometric properties of these faint galaxies are usually compared with current renditions of semianalytical models to provide constraints on the fundamental physical processes involved in galaxy formation and evolution, namely the mass assembly and the star formation. Galaxy counts over large sky areas in the near-UV band are important because they are difficult to obtain given the low efficiency of near-UV instrumentation, even at 8m class telescopes. A large instrumental field of view helps in minimizing the biases due to the cosmic variance. We have obtained deep images in the 360nm U band provided by the blue channel of the Large Binocular Camera at the prime focus of the Large Binocular Telescope. We have derived over an area of ~0.4 sq. deg. the galaxy number counts down to U=27 in the Vega system (corresponding to U=27.86 in the AB system) at a completeness level of 30% reaching the faintest current limit for this wavelength and sky area. The shape of the galaxy counts in the U band can be described by a double power-law, the bright side being consistent with the shape of shallower surveys of comparable or greater areas. The slope bends over significantly at U>23.5 ensuring the convergence of the contribution by star forming galaxies to the EBL in the near-UV band to a value which is more than 70% of the most recent upper limits derived for this band. We have jointly compared our near-UV and K band counts collected from the literature with few selected hierarchical CDM models emphasizing critical issues in the physical description of the galaxy formation and evolution.

We present analysis of the asteroid surface density distribution of main belt asteroids (mean perihelion $\Delta \simeq 2.404$ AU) in five ecliptic latitude fields, $-17 \gtsimeq \beta(\degr) \ltsimeq +15$, derived from deep \textit{Large Binocular Telescope} (LBT) $V-$band (85% completeness limit $V = 21.3$ mag) and \textit{Spitzer Space Telescope} IRAC 8.0 \micron (80% completeness limit $\sim 103 \mu$Jy) fields enabling us to probe the 0.5--1.0 km diameter asteroid population. We discovered 58 new asteroids in the optical survey as well as 41 new bodies in the \textit{Spitzer} fields. The derived power law slopes of the number of asteroids per square degree are similar within each $\sim 5$\degr{} ecliptic latitude bin with a mean value of $ -0.111 \pm 0.077$. For the 23 known asteroids detected in all four IRAC channels mean albedos range from $0.24 \pm 0.07$ to $0.10 \pm 0.05$. No low albedo asteroids ($p_{V}$ $\ltsimeq$ 0.1) were detected in the \textit{Spitzer} FLS fields, whereas in the SWIRE fields they are frequent. The SWIRE data clearly samples asteroids in the middle and outer belts providing the first estimates of these km-sized asteroids' albedos. Our observed asteroid number densities at optical wavelengths are generally consistent with those derived from the Standard Asteroid Model within the ecliptic plane. However, we find an over density at $\beta \gtsimeq 5$\degr{} in our optical fields, while the infrared number densities are under dense by factors of 2 to 3 at all ecliptic latitudes.

Using the 8.4m Large Binocular Telescope, we observed six GRB afterglows from 2.8 hours to 30.8 days after the burst triggers to systematically probe the late time behaviors of afterglows including jet breaks, flares, and supernova bumps. We detected five afterglows with Sloan r' magnitudes ranging from 23.0-26.3 mag. The depth of our observations allows us to extend the temporal baseline for measuring jet breaks by another decade in time scale. We detected two jet breaks and a third candidate, all of which are not detectable without deep, late time optical observations. In the other three cases, we do not detect the jet breaks either because of contamination from the host galaxy light, the presence of a supernova bump, or the intrinsic faintness of the optical afterglow. This suggests that the basic picture that GRBs are collimated is still valid and that the apparent lack of Swift jet breaks is due to poorly sampled afterglow light curves, particularly at late times.

We present the characteristics and some early scientific results of the first instrument at the Large Binocular Telescope (LBT), the Large Binocular Camera (LBC). Each LBT telescope unit will be equipped with similar prime focus cameras. The blue channel is optimized for imaging in the UV-B bands and the red channel for imaging in the VRIz bands. The corrected field-of-view of each camera is approximately 30 arcminutes in diameter, and the chip area is equivalent to a 23x23 arcmin2 field. In this paper we also present the commissioning results of the blue channel. The scientific and technical performance of the blue channel was assessed by measurement of the astrometric distortion, flat fielding, ghosts, and photometric calibrations. These measurements were then used as input to a data reduction pipeline applied to science commissioning data. The measurements completed during commissioning show that the technical performance of the blue channel is in agreement with original expectations. Since the red camera is very similar to the blue one we expect similar performance from the commissioning that will be performed in the following months in binocular configuration. Using deep UV image, acquired during the commissioning of the blue camera, we derived faint UV galaxy-counts in a ~500 sq. arcmin. sky area to U(Vega)=26.5. These galaxy counts imply that the blue camera is the most powerful UV imager presently available and in the near future in terms of depth and extent of the field-of-view. We emphasize the potential of the blue camera to increase the robustness of the UGR multicolour selection of Lyman break galaxies at redshift z~3.

We used a proper combination of multiband high-resolution and wide field multi-wavelength observations collected at three different telescopes (HST, LBT and CFHT) to probe Blue Straggler Star (BSS) populations in the globular cluster M53. Almost 200 BSS have been identified over the entire cluster extension. The radial distribution of these stars has been found to be bimodal (similarly to that of several other clusters) with a prominent dip at ~60'' (~2 r_c) from the cluster center. This value turns out to be a factor of two smaller than the radius of avoidance (r_avoid, the radius within which all the stars of ~1.2 M_sun have sunk to the core because of dynamical friction effects in an Hubble time). While in most of the clusters with a bimodal BSS radial distribution, r_avoid has been found to be located in the region of the observed minimum, this is the second case (after NGC6388) where this discrepancy is noted. This evidence suggests that in a few clusters the dynamical friction seems to be somehow less efficient than expected. We have also used this data base to construct the radial star density profile of the cluster: this is the most extended and accurate radial profile ever published for this cluster, including detailed star counts in the very inner region. The star density profile is reproduced by a standard King Model with an extended core (~25'') and a modest value of the concentration parameter (c=1.58). A deviation from the model is noted in the most external region of the cluster (at r>6.5' from the center). This feature needs to be further investigated in order to address the possible presence of a tidal tail in this cluster.

We present the first deep color-magnitude diagram of the Canes Venatici I (CVnI) dwarf galaxy from observations with the wide field Large Binocular Camera on the Large Binocular Telescope. Reaching down to the main-sequence turnoff of the oldest stars, it reveals a dichotomy in the stellar populations of CVnI: it harbors an old (> 10 Gyr), metal-poor ([Fe/H] ~ -2.0) and spatially extended population along with a much younger (~ 1.4-2.0 Gyr), 0.5 dex more metal-rich, and spatially more concentrated population. These young stars are also offset by 64_{-20}^{+40} pc to the East of the galaxy center. The data suggest that this young population, which represent ~ 3-5 % of the stellar mass of the galaxy within its half-light radius, should be identified with the kinematically cold stellar component found by Ibata et al. (2006). CVnI therefore follows the behavior of the other remote MW dwarf spheroidals which all contain intermediate age and/or young populations: a complex star formation history is possible in extremely low-mass galaxies.

We present a 900 sec, wide-field U image of the inner region of the Andromeda galaxy obtained during the commissioning of the blue channel of the Large Binocular Camera mounted on the prime focus of the Large Binocular Telescope. Relative photometry and absolute astrometry of individual sources in the image was obtained along with morphological parameters aimed at discriminating between stars and extended sources, e.g. globular clusters. The image unveils the near-ultraviolet view of the inner ring of star formation recently discovered in the infrared by the Spitzer Space Telescope and shows in great detail the fine structure of the dust lanes associated with the galaxy inner spiral arms. The capabilities of the blue channel of the Large Binocular Camera at the Large Binocular Telescope (LBC-Blue) are probed by direct comparison with ultraviolet GALEX observations of the same region in M31. We discovered 6 new candidate stellar clusters in this high-background region of M31. We also recovered 62 bona-fide globulars and 62 previously known candidates from the Revised Bologna Catalogue of the M31 globular clusters, and firmly established the extended nature of 19 of them.

We present optical and near-infrared photometry, along with optical spectra, of the Type Ia supernova SN 2000E in the spiral galaxy NGC 6951. It was discovered by the staff of the Teramo Observatory during the monitoring of the SN 1999el. The observations span a time interval of 234 days in the optical and 134 days in the near-infrared (starting 16 days and 7 days before maximum B light, respectively). Optical spectra are available from 6 days before maximum B light to 122 days after it. SN 2000E exhibits a Dm15(B) = 0.94, thus being classifiable as a slow-declining Type Ia SN and showing the distinctive features of such a class of objects. Spectroscopically, SN 2000E appears as a normal Type Ia SN, like SN 1990N. We could constrain reddening [E(B-V) ~0.5 mag] and distance (mu0 ~32.14 mag) using a number of different methods. The bolometric luminosity curve of SN 2000E, allows a determination of the Ni56 mass amounting to 0.9Msun.

Optical and near-infrared light curves of the Type IIn supernova 1999el in NGC 6951 are presented. A period of 220 days (416 days in the near-infrared) is covered from the first observation obtained a few days before maximum light. Spectroscopic observations are also discussed. Using as a distance calibrator the Type Ia SN 2000E, which occurred some months later in the same galaxy, and fitting a blackbody law to the photometric data we obtain a maximum bolometric luminosity for SN 1999el of $\sim 10^{44}$ erg s$^{-1}$. In general, the photometric properties of SN 1999el are very similar to those of SN 1998S, a bright and well studied Type IIn SN, showing a fast decline in all observed bands similar to those of Type II-L SNe. The differences with SN 1998S are analyzed and ascribed to the differences in a pre-existing circumstellar envelope in which dust was already present at the moment of the SN outburst. We infer that light echoes may play a possibly significant role in affecting the observed properties of the light curves, although improved theoretical models are needed to account for the data. We conclude that mass loss in the progenitor RG stars is episodic and occurs in an asymmetric way. This implies that collapsing massive stars appear as normal Type II SN if this occurs far from major mass loss episodes, whereas they appear as Type IIn SNe if a large mass loss episode is in progress.

Until recently the positions of gamma ray bursts were not sufficiently well known within a short timescale to localize and identify them with known celestial sources. Following the historical detection of the X-ray afterglow of the burst GRB970228, extending from 30 s to 3 days after the main peak, by the Beppo-SAX satellite and that of an optical transient 21 hr after the burst, we report here the detection of the same optical transient, in images obtained only 16 hours after the burst.