The Widefield ASKAP L-band Legacy All-sky Blind surveY (WALLABY) is a next-generation survey of neutral hydrogen (HI) in the Local Universe. It uses the widefield, high-resolution capability of the Australian Square Kilometer Array Pathfinder (ASKAP), a radio interferometer consisting of 36 x 12-m dishes equipped with Phased-Array Feeds (PAFs), located in an extremely radio-quiet zone in Western Australia. WALLABY aims to survey three-quarters of the sky (-90 degr < Dec < +30 degr) to a redshift of z < 0.26, and generate spectral line image cubes at ~30 arcsec resolution and ~1.6 mJy/beam per 4 km/s channel sensitivity. ASKAP's instantaneous field of view at 1.4 GHz, delivered by the PAF's 36 beams, is about 30 sq deg. At an integrated signal-to-noise ratio of five, WALLABY is expected to detect over half a million galaxies with a mean redshift of z ~ 0.05 (~200 Mpc). The scientific goals of WALLABY include: (a) a census of gas-rich galaxies in the vicinity of the Local Group; (b) a study of the HI properties of galaxies, groups and clusters, in particular the influence of the environment on galaxy evolution; and (c) the refinement of cosmological parameters using the spatial and redshift distribution of low-bias gas-rich galaxies. For context we provide an overview of previous large-scale HI surveys. Combined with existing and new multi-wavelength sky surveys, WALLABY will enable an exciting new generation of panchromatic studies of the Local Universe. - First results from the WALLABY pilot survey are revealed, with initial data products publicly available in the CSIRO ASKAP Science Data Archive (CASDA).

We report the direct observation of neutrinos produced in the carbon-nitrogen-oxygen (CNO) fusion cycle in the Sun with the Borexino detector at the Laboratori Nazionali del Gran Sasso in Italy. This is the first experimental evidence of the existence of such reaction sequence in a star. The CNO solar neutrino interaction rate is $7.2^{+3.0}_{-1.7}$ counts per day per $100$ tonnes of target at 68% C.L., corresponding to a flux of neutrinos on Earth of $7.0^{+3.0}_{-2.0}\times10^8$ cm$^{-2}$ s$^{-1}$. The absence of CNO signal is disfavoured at 5$\sigma$.

We use SHARK, a semi-analytic galaxy formation model, to investigate the physical processes involved in dictating the shape, scatter and evolution of the HI-halo mass relation at $0\leq z \leq 2$. We compare SHARK with HI clustering and spectral stacking of the HI-halo mass relation derived from observations finding excellent agreement with the former and a deficiency of HI in SHARK at $M_{\rm vir}\approx 10^{12-13} M_{\odot}$ in the latter, but otherwise great agreement below and above that mass threshold. In SHARK, we find that the HI mass increases with the halo mass up to a critical mass of $\approx 10^{11.8} M_{\odot}$; between $\sim 10^{11.8}-10^{13}M_{\odot}$, the scatter in the relation increases by 0.7 dex and the HI mass decreases with the halo mass on average; at $M_{\rm vir} \geq 10^{13} M_{\odot}$, the HI content continues to increase with halo mass. We find that the critical halo mass of $\approx 10^{12} M_{\odot}$ is largely set by feedback from Active Galactic Nuclei (AGN), and the exact shape and scatter of the HI-halo mass relation around that mass is extremely sensitive to how AGN feedback is modelled, with other physical processes playing a less significant role. We determine the main secondary parameters responsible for the scatter of the HI-halo mass relation, namely the halo spin parameter at $M_{\rm vir}\leq 10^{11.8} M_{\odot}$, and the fractional contribution from substructure to the total halo mass for $M_{\rm vir}\geq 10^{13} M_{\odot}$. The scatter at $10^{11.8}<M_{\rm vir}<10^{13} M_{\odot}$ is best described by the black-hole-to-stellar mass ratio of the central galaxy, reflecting the AGN feedback relevance. We present a numerical model to populate dark matter-only simulations with HI at $0\leq z \leq 2$ based solely on halo parameters that are measurable in such simulations.

During a core-collapse supernova (SN), axionlike particles (ALPs) could be produced through the Primakoff process and subsequently convert into $\gamma$ rays in the magnetic field of the Milky Way. We do not find evidence for such a $\gamma$-ray burst in observations of extragalactic SNe with the Fermi Large Area Telescope (LAT). The SN explosion times are estimated from optical light curves and we find a probability of about $\sim$90% that the LAT observed at least one SN at the time of the core collapse. Under the assumption that at least one SN was contained within the LAT field of view, we exclude photon-ALP couplings $\gtrsim2.6\times10^{-12}$ GeV$^{-1}$ for ALP masses $m_a\lesssim3\times10^{-10}$ eV, improving previous limits from SN1987A by a factor of 2.

Context. Detecting and characterizing substellar companions for which the luminosity, mass, and age can be determined independently is of utter importance to test and calibrate the evolutionary models due to uncertainties in their formation mechanisms. HD 19467 is a bright and nearby star hosting a cool brown dwarf companion detected with RV and imaging, making it a valuable object for such studies. Aims. We aim to further characterize the orbital, spectral, and physical properties of the HD 19467 system. Methods. We present new high-contrast imaging data with the SPHERE and NaCo instruments. We also analyze archival data from HARPS, NaCo, HIRES, UVES, and ASAS. We also use proper motion data of the star from Hipparcos and Gaia. Results. We refine the properties of the host star and derive an age of 8.0$^{+2.0}_{-1.0}$ Gyr based on isochrones, gyrochronology, and chemical and kinematic arguments. This estimate is slightly younger than previous estimates of ~9-11 Gyr. No orbital curvature is seen in the current imaging, RV, and astrometric data. From a joint fit of the data, we refine the orbital parameters for HD 19467B: period 398$^{+95}_{-93}$ yr, inclination 129.8$^{+8.1}_{-5.1}$ deg, eccentricity 0.56$\pm$0.09, longitude of the ascending node 134.8$\pm$4.5 deg, and argument of the periastron 64.2$^{+5.5}_{-6.3}$ deg. We assess a dynamical mass of 74$^{+12}_{-9}$ MJ. The fit with atmospheric models of the spectrophotometric data of HD 19467B indicates an atmosphere without clouds or with very thin clouds, an effective temperature of 1042$^{+77}_{-71}$ K, and a large surface gravity of 5.34$^{+0.08}_{-0.09}$ dex. The comparison to model predictions of the bolometric luminosity and dynamical mass of HD 19467B, assuming our system age estimate, indicates a better agreement with the Burrows et al. models; whereas the other evolutionary models used tend to underestimate its cooling rate.

Neutrinos emitted in the carbon, nitrogen, oxygen (CNO) fusion cycle in the Sun are a sub-dominant, yet crucial component of solar neutrinos whose flux has not been measured yet. The Borexino experiment at the Laboratori Nazionali del Gran Sasso (Italy) has a unique opportunity to detect them directly thanks to the detector's radiopurity and the precise understanding of the detector backgrounds. We discuss the sensitivity of Borexino to CNO neutrinos, which is based on the strategies we adopted to constrain the rates of the two most relevant background sources, $pep$ neutrinos and $^{210}$Bi beta decays. Finally, the relevance of a detection of CNO neutrinos is extensively discussed in the context of the solar metallicity problem.

Protoplanetary disks in dense, massive star-forming regions are strongly affected by their environment. How this environmental impact changes over time is an important constraint on disk evolution and external photoevaporation models. We characterize the dust emission from 179 disks in the core of the young (0.5 Myr) NGC 2024 cluster. By studying how the disk mass varies within the cluster, and comparing these disks to those in other regions, we aim to determine how external photoevaporation influences disk properties over time. Using the Atacama Large Millimeter/submillimeter Array, a 2.9' x 2.9' mosaic centered on NGC 2024 FIR 3 was observed at 225 GHz with a resolution of 0.25'', or ~100 AU. The imaged region contains 179 disks identified at IR wavelengths, seven new disk candidates, and several protostars. The overall detection rate of disks is $32 \pm 4\%$. Few of the disks are resolved, with the exception of a giant (R = 300 AU) transition disk. Serendipitously, we observe a millimeter flare from an X-ray bright young stellar object (YSO), and resolve continuum emission from a Class 0 YSO in the FIR 3 core. Two distinct disk populations are present: a more massive one in the east, along the dense molecular ridge hosting the FIR 1-5 YSOs, with a detection rate of $45 \pm 7\%$. In the western population, towards IRS 1, only $15 \pm 4\%$ of disks are detected. NGC 2024 hosts two distinct disk populations. Disks along the dense molecular ridge are young (0.2 - 0.5 Myr) and partly shielded from the far ultraviolet radiation of IRS 2b; their masses are similar to isolated 1 - 3 Myr old SFRs. The western population is older and at lower extinctions, and may be affected by external photoevaporation from both IRS 1 and IRS 2b. However, it is possible these disks had lower masses to begin with.

Recent high-contrast imaging surveys, looking for planets in young, nearby systems showed evidence of a small number of giant planets at relatively large separation beyond typically 20 au where those surveys are the most sensitive. Access to smaller physical separations between 5 and 20 au is the next step for future planet imagers on 10 m telescopes and ELTs in order to bridge the gap with indirect techniques (radial velocity, transit, astrometry with Gaia). In that context, we recently proposed a new algorithm, Keplerian-Stacker, combining multiple observations acquired at different epochs and taking into account the orbital motion of a potential planet present in the images to boost the ultimate detection limit. We showed that this algorithm is able to find planets in time series of simulated images of SPHERE even when a planet remains undetected at one epoch. Here, we validate the K-Stacker algorithm performances on real SPHERE datasets, to demonstrate its resilience to instrumental speckles and the gain offered in terms of true detection. This will motivate future dedicated multi-epoch observation campaigns in high-contrast imaging to search for planets in emitted and reflected light. Results. We show that K-Stacker achieves high success rate when the SNR of the planet in the stacked image reaches 7. The improvement of the SNR ratio goes as the square root of the total exposure time. During the blind test and the redetection of HD 95086 b, and betaPic b, we highlight the ability of K-Stacker to find orbital solutions consistent with the ones derived by the state of the art MCMC orbital fitting techniques, confirming that in addition to the detection gain, K-Stacker offers the opportunity to characterize the most probable orbital solutions of the exoplanets recovered at low signal to noise.

Axion-like particles (ALPs), predicted in theories beyond the Standard Model, can have observational effects on the transparency of the Universe to $\gamma$ rays in the presence of magnetic fields. In this work, we search for effects compatible with the existence of ALPs with 80 months of data from the ${\it Fermi}$ Large Area Telescope, by comparing the distributions of observed highest energy photons from sources beyond redshifts of z $\geq$ 0.1 with theoretical predictions in the presence of ALPs. We find no evidence for an increased $\gamma$-ray transparency due to ALPs and therefore we set limits on the ALPs parameters assuming a value of the intergalactic magnetic field strength of 1 nG. Photon-ALP couplings above $10^{-11}$ GeV$^{-1}$ are excluded for ALP masses $m_{a}$ $\lesssim3.0$ neV. These limits exclude a region of the parameter space not covered by other $\gamma$-ray telescopes and are compatible with constraints imposed by other experiments.

The most widely-studied mechanism of mass loss from extrasolar planets is photoevaporation via XUV ionization, primarily in the context of highly irradiated planets. However, the EUV dissociation of hydrogen molecules can also theoretically drive atmospheric evaporation on low-mass planets. For temperate planets such as the early Earth, impact erosion is expected to dominate in the traditional planetesimal accretion model, but it would be greatly reduced in pebble accretion scenarios, allowing other mass loss processes to be major contributors. We apply the same prescription for photoionization to this photodissociation mechanism and compare it to an analysis of other possible sources of mass loss in pebble accretion scenarios. We find that there is not a clear path to evaporating the primordial atmosphere accreted by an early Earth analog in a pebble accretion scenario. Impact erosion could remove ~2,300 bars of hydrogen if 1% of the planet's mass is accreted as planetesimals, while the combined photoevaporation processes could evaporate ~750 bars of hydrogen. Photodissociation is likely a subdominant, but significant component of mass loss. Similar results apply to super-Earths and mini-Neptunes. This mechanism could also preferentially remove hydrogen from a planet's primordial atmosphere, thereby leaving a larger abundance of primordial water compared to standard dry formation models. We discuss the implications of these results for models of rocky planet formation including Earth's formation and the possible application of this analysis to mass loss from observed exoplanets.

Proxima Centauri is known to host an earth-like planet in its habitable zone; very recently a second candidate planet was proposed based on radial velocities. At quadrature, the expected projected separation of this new candidate is larger than 1 arcsec, making it a potentially interesting target for direct imaging. While difficult, identification of the optical counterpart of this planet would allow detailed characterization of the closest planetary system. We searched for a counterpart in SPHERE images acquired during four years through the SHINE survey. In order to account for the large orbital motion of the planet, we used a method that assumes the circular orbit obtained from radial velocities and exploits the sequence of observations acquired close to quadrature in the orbit. We checked this with a more general approach that considers keplerian motion, K-stacker. We did not obtain a clear detection. The best candidate has S/N=6.1 in the combined image. A statistical test suggests that the probability that this detection is due to random fluctuation of noise is < 1% but this result depends on the assumption that distribution of noise is uniform over the image. The position of this candidate and the orientation of its orbital plane fit well with observations in the ALMA 12m array image. However, the astrometric signal expected from the orbit of the candidate we detected is 3-sigma away from the astrometric motion of Proxima as measured from early Gaia data. This, together with the unexpectedly high flux associated with our direct imaging detection, means we cannot confirm that our candidate is indeed Proxima c. On the other hand, if confirmed, this would be the first observation in imaging of a planet discovered from radial velocities and the second one (after Fomalhaut b) of reflecting circumplanetary material. Further confirmation observations should be done as soon as possible.

We report on an improved measurement of the $^8$B solar neutrino interaction rate with the Borexino experiment at the Laboratori Nazionali del Gran Sasso. Neutrinos are detected via their elastic scattering on electrons in a large volume of liquid scintillator. The measured rate of scattered electrons above 3 MeV of energy is $0.223\substack{+0.015 \\ -0.016}\,(stat)\,\substack{+0.006 \\ -0.006}\,(syst)$ cpd/100 t, which corresponds to an observed solar neutrino flux assuming no neutrino flavor conversion of $\Phi\substack{\rm ES \\ ^8\rm B}=2.57\substack{+0.17 \\ -0.18}(stat)\substack{+0.07\\ -0.07}(syst)\times$10$^6$ cm$^{-2}\,$s$^{-1}$. This measurement exploits the active volume of the detector in almost its entirety for the first time, and takes advantage of a reduced radioactive background following the 2011 scintillator purification campaign and of novel analysis tools providing a more precise modeling of the background. Additionally, we set a new limit on the interaction rate of solar $hep$ neutrinos, searched via their elastic scattering on electrons as well as their neutral current-mediated inelastic scattering on carbon, $^{12}$C($\nu,\nu'$)$^{12}$C* ($E_{\gamma}$= 15.1 MeV).

Directly imaged planets and substellar companions are key targets for the characterization of self-luminous atmospheres. Their photometric appearance at 4-5 $\mu$m is sensitive to the chemical composition and cloud content of their atmosphere. We aim at systematically characterizing the atmospheres of directly imaged low-mass companions at 4-5 $\mu$m. We want to homogeneously process the data and compile a photometric library at thermal wavelengths of these mostly young, low-gravity objects. In this way, we want to find trends related to their spectral type and surface gravity by comparing with isolated brown dwarfs and predictions from atmospheric models. We have used the high-resolution, high-contrast capabilities of NACO at the VLT to directly image the companions of HIP 65426, PZ Tel, and HD 206893 in the NB4.05 and/or $M'$ filters. For the same targets, and additionally $\beta$ Pic, we have also analyzed six archival VLT/NACO datasets which were taken with the NB3.74, $L'$, NB4.05, and $M'$ filters. The $L'$-NB4.05 and $L'$-$M'$ colors of the studied sample are all red while the NB4.05-$M'$ color is blue for $\beta$ Pic b, gray for PZ Tel B, and red for HIP 65426 b and HD 206893 B (although typically with low significance). The absolute NB4.05 and $M'$ fluxes of our sample are all larger than those of field dwarfs with similar spectral types. Finally, the surface gravity of $\beta$ Pic b has been constrained to $\log{g} = 4.17_{-0.13}^{+0.10}$ dex from its photometry and dynamical mass. A red color at 3-4 $\mu$m and a blue color at 4-5 $\mu$m might be (partially) caused by H$_2$O and CO absorption, respectively. The red characteristics of $\beta$ Pic b, HIP 65426 b, and HD 206893 B at 3-5$\mu$m, as well as their higher fluxes in NB4.05 and $M'$ compared to field dwarfs, indicate that cloud densities are enhanced close to the photosphere as a result of their low surface gravity.

Almost 10 yr of $\gamma$-ray observations with the Fermi Large Area Telescope (LAT) have revealed extreme $\gamma$-ray outbursts from flat spectrum radio quasars (FSRQs), temporarily making these objects the brightest $\gamma$-ray emitters in the sky. Yet, the location and mechanisms of the $\gamma$-ray emission remain elusive. We characterize long-term $\gamma$-ray variability and the brightest $\gamma$-ray flares of six FSRQs. Consecutively zooming in on the brightest flares, which we identify in an objective way through Bayesian blocks and a hill-climbing algorithm, we find variability on subhour time scales and as short as minutes for two sources in our sample (3C279, CTA102) and weak evidence for variability at time scales less than the Fermi satellite's orbit of 95 minutes for PKS1510-089 and 3C454.3. This suggests extremely compact emission regions in the jet. We do not find any signs for $\gamma$-ray absorption in the broad-line region (BLR), which indicates that $\gamma$-rays are produced at distances greater than hundreds of gravitational radii from the central black hole. This is further supported by a cross-correlation analysis between $\gamma$-ray and radio/millimeter light curves, which is consistent with $\gamma$-ray production at the same location as the millimeter core for 3C273, CTA102, and 3C454.3. The inferred locations of the $\gamma$-ray production zones are still consistent with the observed decay times of the brightest flares if the decay is caused by external Compton scattering with BLR photons. However, the minute-scale variability is challenging to explain in such scenarios.

We report on the observations of gamma-ray burst (GRB) 190114C by the Fermi Gamma-ray Space Telescope and the Neil Gehrels Swift Observatory. The early-time observations reveal multiple emission components that evolve independently, with a delayed power-law component that exhibits significant spectral attenuation above 40 MeV in the first few seconds of the burst. This power-law component transitions to a harder spectrum that is consistent with the afterglow emission observed at later times. This afterglow component is clearly identifiable in the GBM and BAT light curves as a slowly fading emission component on which the rest of the prompt emission is superimposed. As a result, we are able to constrain the transition from internal shock to external shock dominated emission. We find that the temporal and spectral evolution of the broadband afterglow emission can be well modeled as synchrotron emission from a forward shock propagating into a wind-like circumstellar environment and find that high-energy photons observed by Fermi LAT are in tension with the theoretical maximum energy that can be achieved through synchrotron emission from a shock. These violations of the maximum synchrotron energy are further compounded by the detection of very high energy (VHE) emission above 300 GeV by MAGIC concurrent with our observations. We conclude that the observations of VHE photons from GRB 190114C necessitates either an additional emission mechanism at very high energies that is hidden in the synchrotron component in the LAT energy range, an acceleration mechanism that imparts energy to the particles at a rate that is faster than the electron synchrotron energy loss rate, or revisions of the fundamental assumptions used in estimating the maximum photon energy attainable through the synchrotron process.

The Borexino detector measures solar neutrino fluxes via neutrino-electron elastic scattering. Observed spectra are determined by the solar-$\nu_{e}$ survival probability $P_{ee}(E)$, and the chiral couplings of the neutrino and electron. Some theories of physics beyond the Standard Model postulate the existence of Non-Standard Interactions (NSI's) which modify the chiral couplings and $P_{ee}(E)$. In this paper, we search for such NSI's, in particular, flavor-diagonal neutral current interactions that modify the $\nu_e e$ and $\nu_\tau e$ couplings using Borexino Phase II data. Standard Solar Model predictions of the solar neutrino fluxes for both high- and low-metallicity assumptions are considered. No indication of new physics is found at the level of sensitivity of the detector and constraints on the parameters of the NSI's are placed. In addition, with the same dataset the value of $\sin^2\theta_W$ is obtained with a precision comparable to that achieved in reactor antineutrino experiments.

We performed observations of the Sco-Cen F star HD 117214 aiming at a search for planetary companions and the characterization of the debris disk structure. HD 117214 was observed with the SPHERE subsystems IRDIS, IFS and ZIMPOL at optical and near-IR wavelengths using angular and polarimetric differential imaging techniques. This provided the first images of scattered light from the debris disk with a spatial resolution reaching 25 mas and an inner working angle $< 0.1''$. With the observations with IRDIS and IFS we derive detection limits for substellar companions. The geometrical parameters of the detected disk are constrained by fitting 3D models for the scattering of an optically thin dust disk. Investigating the possible origin of the disk gap, we introduced putative planets therein and modeled the planet-disk and planet-planet dynamical interactions. The obtained planetary architectures are compared with the detection limit curves. The debris disk has an axisymmetric ring structure with a radius of $0.42(\pm 0.01)''$ or $\sim45$ au and an inclination of $71(\pm 2.5)^\circ$ and exhibits a $0.4''$ ($\sim40$ au) wide inner cavity. From the polarimetric data, we derive a polarized flux contrast for the disk of $(F_{\rm pol})_{\rm disk}/F_{\rm \ast}> (3.1 \pm 1.2)\cdot 10^{-4}$ in the RI band. The fractional scattered polarized flux of the disk is eight times smaller than the fractional infrared flux excess. This ratio is similar to the one obtained for the debris disk HIP 79977 indicating that dust radiation properties are not very different between these two disks. Inside the disk cavity we achieve the high sensitivity limits on planetary companions with a mass down to $\sim 4 M_{\rm J}$ at projected radial separations between $0.2''$ and $0.4''$. We can exclude the stellar companions at a radial separation larger than 75 mas from the star.

We present the first simultaneous measurement of the interaction rates of $pp$, $^7$Be, and $pep$ solar neutrinos performed with a global fit to the Borexino data in an extended energy range (0.19-2.93)$\,$MeV. This result was obtained by analyzing 1291.51$\,$days of Borexino Phase-II data, collected between December 2011 and May 2016 after an extensive scintillator purification campaign. We find: rate($pp$)$\,$=$\,$$134$$\,$$\pm$$\,$$10$$\,$($stat$)$\,$$^{\rm +6}_{\rm -10}$$\,$($sys$)$\,$cpd/100$\,$t, rate($^7$Be)$\,$=$\,$$48.3$$\,$$\pm$$\,$$1.1$$\,$($stat$)$\,$$^{\rm +0.4}_{\rm -0.7}$$\,$($sys$)$\,$cpd/100$\,$t, and rate($pep$)$\,$=$\,$$2.43$$\pm$$\,$$0.36$$\,$($stat$)$^{+0.15}_{-0.22}$$\,$($sys$)$\,$cpd/100$\,$t. These numbers are in agreement with and improve the precision of our previous measurements. In particular, the interaction rate of $^7$Be $\nu$'s is measured with an unprecedented precision of 2.7%, showing that discriminating between the high and low metallicity solar models is now largely dominated by theoretical uncertainties. The absence of $pep$ neutrinos is rejected for the first time at more than 5$\,$$\sigma$. An upper limit of $8.1$$\,$cpd/100$\,$t (95%$\,$C.L.) on the CNO neutrino rate is obtained by setting an additional constraint on the ratio between the $pp$ and $pep$ neutrino rates in the fit. This limit has the same significance as that obtained by the Borexino Phase-I (currently providing the tightest bound on this component), but is obtained by applying a less stringent constraint on the $pep$ $\nu$ flux.

RefPlanets is a guaranteed time observation (GTO) programme that uses the Zurich IMaging POLarimeter (ZIMPOL) of SPHERE/VLT for a blind search for exoplanets in wavelengths from 600-900 nm. The goals of this study are the characterization of the unprecedented high polarimetic contrast and polarimetric precision capabilities of ZIMPOL for bright targets, the search for polarized reflected light around some of the closest bright stars to the Sun and potentially the direct detection of an evolved cold exoplanet for the first time. For our observations of Alpha Cen A and B, Sirius A, Altair, Eps Eri and Tau Ceti we used the polarimetric differential imaging (PDI) mode of ZIMPOL which removes the speckle noise down to the photon noise limit for angular separations >0.6". We describe some of the instrumental effects that dominate the noise for smaller separations and explain how to remove these additional noise effects in post-processing. We then combine PDI with angular differential imaging (ADI) as a final layer of post-processing to further improve the contrast limits of our data at these separations. For good observing conditions we achieve polarimetric contrast limits of 15.0-16.3 mag at the effective inner working angle of about 0.13", 16.3-18.3 mag at 0.5" and 18.8-20.4 mag at 1.5". The contrast limits closer in (<0.6") depend significantly on the observing conditions, while in the photon noise dominated regime (>0.6"), the limits mainly depend on the brightness of the star and the total integration time. We compare our results with contrast limits from other surveys and review the exoplanet detection limits obtained with different detection methods. For all our targets we achieve unprecedented contrast limits. Despite the high polarimetric contrasts we are not able to find any additional companions or extended polarized light sources in the data that has been taken so far.

The All-sky Medium Energy Gamma-ray Observatory (AMEGO) is a probe class mission concept that will provide essential contributions to multimessenger astrophysics in the late 2020s and beyond. AMEGO combines high sensitivity in the 200 keV to 10 GeV energy range with a wide field of view, good spectral resolution, and polarization sensitivity. Therefore, AMEGO is key in the study of multimessenger astrophysical objects that have unique signatures in the gamma-ray regime, such as neutron star mergers, supernovae, and flaring active galactic nuclei. The order-of-magnitude improvement compared to previous MeV missions also enables discoveries of a wide range of phenomena whose energy output peaks in the relatively unexplored medium-energy gamma-ray band.

Context. PDS 70 is a young (5.4 Myr), nearby (~113 pc) star hosting a known transition disk with a large gap. Recent observations with SPHERE and NACO in the near-infrared (NIR) allowed us to detect a planetary mass companion, PDS70b, within the disk cavity. Moreover, observations in H_alpha with MagAO and MUSE revealed emission associated to PDS70b and to another new companion candidate, PDS70c, at a larger separation from the star. Aims. Our aim is to confirm the discovery of the second planet PDS70c using SPHERE at VLT, to further characterize its physical properties, and search for additional point sources in this young planetary system. Methods. We re-analyzed archival SPHERE NIR observations and obtained new data in Y, J, H and K spectral bands for a total of four different epochs. The data were reduced using the data reduction and handling pipeline and the SPHERE data center. We then applied custom routines (e.g. ANDROMEDA and PACO) to subtract the starlight. Results. We re-detect both PDS 70 b and c and confirm that PDS70c is gravitationally bound to the star. We estimate this second planet to be less massive than 5 M Jup and with a T_eff around 900 K. Also, it has a low gravity with log g between 3.0 and 3.5 dex. In addition, a third object has been identified at short separation (~0.12") from the star and gravitationally bound to the star. Its spectrum is however very blue, so that we are probably seeing stellar light reflected by dust and our analysis seems to demonstrate that it is a feature of the inner disk. We, however, cannot completely exclude the possibility that it is a planetary mass object enshrouded by a dust envelope. In this latter case, its mass should be of the order of few tens of M_Earth. Moreover, we propose a possible structure for the planetary system based on our data that, however, cannot be stable on a long timescale.

We present the results of a binary population study in the Orion Nebula Cluster (ONC) using archival Hubble Space Telescope (HST) data obtained with the Advanced Camera for Surveys (ACS) in Johnson V filter (HST Proposal 10246, PI M. Robberto). Young clusters and associations hold clues to the origin and properties of multiple star systems. Binaries with separations $< 100 $ AU are useful as tracers of the initial binary population since they are not as likely to be destroyed through dynamical interactions. Low mass, low stellar density star-forming regions such as Taurus-Auriga, reveal an excess of multiples compared to the Galactic Field. Studying the binary population of higher mass, higher stellar density star-forming regions like the ONC provides useful information concerning the origin of the Galactic Field star population. In this survey, we characterize the previously unexplored (and incomplete) separation parameter space of binaries in the ONC (15 - 160 AU) by fitting a double-PSF model built from empirical PSFs. We identified 14 candidate binaries (11 new detections) and find that 8$_{-2\%}^{+4\%}$ of our observed sample are in binary systems, complete over mass ratios and separations of 0.6 $< $ q $< $ 1.0 and 30 $< $ a $< $ 160 AU. This is consistent with the Galactic Field M-dwarf population over the same parameter ranges, 6.5\% $\pm$ 3\%. Therefore, high mass star forming regions like the ONC would not require further dynamical evolution for their binary population to resemble the Galactic Field, as some models have hypothesized for young clusters.

Blazars and in particular the subclass of high synchrotron peaked Active Galactic Nuclei are among the main targets for the present generation of Imaging Atmospheric Cherenkov Telescopes (IACTs) and will remain of great importance for very high-energy $\gamma$-ray science in the era of the Cherenkov Telescope Array (CTA). Observations by IACTs, which have relatively small fields of view ($\sim$ few degrees), are limited by viewing conditions; therefore, it is important to select the most promising targets in order to increase the number of detections. The aim of this paper is to search for unclassified blazars among known $\gamma$-ray sources from the Fermi Large Area Telescope (LAT) third source catalog that are likely detectable with IACTs or CTA. We use an artificial neural network algorithm and updated analysis of Fermi-LAT data. We found 80 $\gamma$-ray source candidates, and for the highest-confidence candidates, we calculate their potential detectability with IACTs and CTA based on an extrapolation of their energy spectra. Follow-up observations of our source candidates could significantly increase the current TeV source population sample and could ultimately confirm the efficiency of our algorithm to select TeV sources.

PKS 2023-07 is a flat spectrum radio quasar located at a redshift $z=1.388$, farther than any source currently detected at very high energies ($E>100$ GeV). At such energies, absorption by the extragalactic background light (EBL) renders the detection of distant sources particularly challenging. The High Energy Stereoscopic System (H.E.S.S.) observed the source following reports from AGILE (April 2016) and Fermi-LAT (April 2016, September and October 2017) on high-flux states in gamma rays. During each of the three flaring periods, near-simultaneous observations were obtained with H.E.S.S., Fermi-LAT and multiple telescopes at other wavelengths. Though the source was not significantly detected by H.E.S.S., upper limits were derived for each observation period. Through constraints given by Fermi-LAT in the MeV--GeV domain and differential upper limits by H.E.S.S., we searched for an intrinsic cutoff in the EBL-corrected gamma ray spectrum of PKS 2023-07.

Exoplanet science is one of the most thriving fields of modern astrophysics. A major goal is the atmospheric characterization of dozens of small, terrestrial exoplanets in order to search for signatures in their atmospheres that indicate biological activity, assess their ability to provide conditions for life as we know it, and investigate their expected atmospheric diversity. None of the currently adopted projects or missions, from ground or in space, can address these goals. In this White Paper we argue that a large space-based mission designed to detect and investigate thermal emission spectra of terrestrial exoplanets in the MIR wavelength range provides unique scientific potential to address these goals and surpasses the capabilities of other approaches. While NASA might be focusing on large missions that aim to detect terrestrial planets in reflected light, ESA has the opportunity to take leadership and spearhead the development of a large MIR exoplanet mission within the scope of the "Voyage 2050" long-term plan establishing Europe at the forefront of exoplanet science for decades to come. Given the ambitious science goals of such a mission, additional international partners might be interested in participating and contributing to a roadmap that, in the long run, leads to a successful implementation. A new, dedicated development program funded by ESA to help reduce development and implementation cost and further push some of the required key technologies would be a first important step in this direction. Ultimately, a large MIR exoplanet imaging mission will be needed to help answer one of mankind's most fundamental questions: "How unique is our Earth?"

Using a spectral stacking technique, we measure the neutral hydrogen (HI) properties of a sample of galaxies at $z < 0.11$ across 35 pointings of the Westerbork Synthesis Radio Telescope (WSRT). The radio data contains 1,895 galaxies with redshifts and positions known from the Sloan Digital Sky Survey (SDSS). We carefully quantified the effects of sample bias, aperture used to extract spectra, sidelobes and weighting technique and use our data to provide a new estimate for the cosmic HI mass density. We find a cosmic HI mass density of $\Omega_{\rm HI} = (4.02 \pm 0.26)\times 10^{-4} h_{70}^{-1}$ at $\langle z\rangle = 0.066$, consistent with measurements from blind HI surveys and other HI stacking experiments at low redshifts. The combination of the small interferometer beam size and the large survey volume makes our result highly robust against systematic effects due to confusion at small scales and cosmic variance at large scales. Splitting into three sub-samples with $\langle z\rangle$ = 0.038, 0.067 and 0.093 shows no significant evolution of the HI gas content at low redshift.

The Cherenkov Telescope Array (CTA) is the next generation ground-based observatory for $\gamma$-ray astronomy at energies above 30 GeV. Thanks to its unique capabilities, CTA observations will address a plethora of open questions in astrophysics ranging from the origin of cosmic messengers to the exploration of the frontiers of physics. In this note, we present a comprehensive sensitivity study to assess the potential of CTA to measure the $\gamma$-ray absorption on the extragalactic background light (EBL), to constrain or detect intergalactic magnetic fields (IGMFs), and probe physics beyond the Standard Model such as axion-like particles (ALPs) and Lorentz invariance violation (LIV), which could modify the $\gamma$-ray spectra features expected from EBL absorption. Our results suggest that CTA will have unprecedented sensitivity to detect IGMF signatures and will probe so-far unexplored regions of the LIV and ALP parameter space. Furthermore, an indirect measurement of the EBL and of its evolution will be performed with unrivaled precision.

Open-access telescopes of all apertures are needed to operate a competitive and efficient national science program. While larger facilities contribute light-gathering power and angular resolution, smaller ones dominate for field of view, time-resolution, and especially, total available observing time, thereby enabling our entire, diversely-expert community. Smaller aperture telescopes therefore play a critical and indispensable role in advancing science. Thus, the divestment of NSF support for modest-aperture (1 - 4 m) public telescopes poses a serious threat to U.S. scientific leadership, which is compounded by the unknown consequences of the shift from observations driven by individual investigators to survey-driven science. Given the much higher cost efficiency and dramatic science returns for investments in modest aperture telescopes, it is hard to justify funding only the most expensive facilities. We therefore urge the Astro2020 panel to explicitly make the case for modest aperture facilities, and to recommend enhancing this funding stream to support and grow this critical component of the OIR System. Further study is urgently needed to prioritize the numerous exciting potential capabilities of smaller facilities,and to establish sustainable, long-term planning for the System.

Jets are rarely associated with pre-main-sequence intermediate-mass stars. Optical and near-IR observations of jet-driving sources are often hindered by the presence of a natal envelope. Jets around partly embedded sources are a useful diagnostic to constrain the geometry of the concealed protoplanetary disk. In fact, the jet-driving mechanisms are affected by both spatial anisotropies and episodic variations at the (sub-)au scale from the star. We obtained a rich set of high-contrast VLT/SPHERE observations from 0.6 micron to 2.2 micron of the young intermediate-mass star RY Tau. Given the proximity to the Sun of this source, our images have the highest spatial resolution ever obtained for an atomic jet. Optical observations in polarized light show no sign of the protoplanetary disk detected by ALMA. Instead, we observed a diffuse signal resembling a remnant envelope with an outflow cavity. The jet is detected in four spectral lines. The jet appears to be wiggling and its radial width increasing with the distance is complementary to the shape of the outflow cavity suggesting a strong jet/envelope interaction. Through the estimated tangential velocity, we revealed a possible connection between the launching time of the jet sub-structures and the stellar activity of RY Tau. RY Tau is at an intermediate stage toward the dispersal of the natal envelope. This source shows episodic increases of mass accretion/ejection similarly to other known intermediate-mass stars. The amount of observed jet wiggle is consistent with the presence of a precessing disk warp or misaligned inner disk that would be induced by an unseen planetary/sub-stellar companion at sub-/few-au scales. The high disk mass of RY Tau and of two other jet-driving intermediate-mass stars, HD163296 and MWC480, suggests that massive, full disks are more efficient at launching prominent jets.

HD163296 is a Herbig Ae/Be star known to host a protoplanetary disk with a ringed structure. To explain the disk features, previous works proposed the presence of planets embedded into the disk. We have observed HD163296 with the near-infrared (NIR) branch of SPHERE composed by IRDIS and IFS with the aim to put tight constraints on the presence of substellar companions around this star. Despite the low rotation of the field of view during our observation we were able to put upper mass limits of few M_Jup around this object. These limits do not allow to give any definitive conclusion about the planets proposed through the disk characteristics. On the other hand, our results seem to exclude the presence of the only candidate proposed until now using direct imaging in the NIR even if some caution has to be taken considered the different wavelength bands of the two observations.

The flat spectrum radio quasar 3C 279 is known to exhibit pronounced variability in the high-energy ($100\,$MeV$<E<100\,$GeV) $\gamma$-ray band, which is continuously monitored with Fermi-LAT. During two periods of high activity in April 2014 and June 2015 Target-of-Opportunity observations were undertaken with H.E.S.S. in the very-high-energy (VHE, $E>100\,$GeV) $\gamma$-ray domain. While the observation in 2014 provides an upper limit, the observation in 2015 results in a signal with $8.7\,\sigma$ significance above an energy threshold of $66\,$GeV. No VHE variability has been detected during the 2015 observations. The VHE photon spectrum is soft and described by a power-law index of $4.2\pm 0.3$. The H.E.S.S. data along with a detailed and contemporaneous multiwavelength data set provide constraints on the physical parameters of the emission region. The minimum distance of the emission region from the central black hole is estimated using two plausible geometries of the broad-line region and three potential intrinsic spectra. The emission region is confidently placed at $r\gtrsim 1.7\times10^{17}\,$cm from the black hole, i.e., beyond the assumed distance of the broad-line region. Time-dependent leptonic and lepto-hadronic one-zone models are used to describe the evolution of the 2015 flare. Neither model can fully reproduce the observations, despite testing various parameter sets. Furthermore, the H.E.S.S. data are used to derive constraints on Lorentz invariance violation given the large redshift of 3C 279.

This paper reports on the atomic hydrogen gas (HI) observations of the spiral galaxy NGC 1566 using the newly commissioned Australian Square Kilometre Array Pathfinder (ASKAP) radio telescope. We measure an integrated HI flux density of $180.2$ Jy km s$^{-1}$ emanating from this galaxy, which translates to an HI mass of $1.94\times10^{10}$M$_\circ$ at an assumed distance of $21.3$ Mpc. Our observations show that NGC 1566 has an asymmetric and mildly warped HI disc. The HI-to-stellar mass fraction of NGC 1566 is $0.29$, which is high in comparison with galaxies that have the same stellar mass ($10^{10.8}$M$_\circ$). We also derive the rotation curve of this galaxy to a radius of $50$ kpc and fit different mass models to it. The NFW, Burkert and pseudo-isothermal dark matter halo profiles fit the observed rotation curve reasonably well and recover dark matter fractions of $0.62$, $0.58$ and $0.66$, respectively. Down to the column density sensitivity of our observations ($N_{HI} = 3.7\times10^{19}$ cm$^{-2}$), we detect no HI clouds connected to, or in the nearby vicinity of, the HI disc of NGC 1566 nor nearby interacting systems. We conclude that, based on a simple analytic model, ram pressure interactions with the IGM can affect the HI disc of NGC 1566 and is possibly the reason for the asymmetries seen in the HI morphology of NGC 1566.

Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We discuss how LSST will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational couplings to the Standard Model, and compact object abundances. Additionally, we discuss the ways that LSST will complement other experiments to strengthen our understanding of the fundamental characteristics of dark matter. More information on the LSST dark matter effort can be found at https://lsstdarkmatter.github.io/ .

Planetary systems hold the imprint of the formation and of the evolution of planets especially at young ages, and in particular at the stage when the gas has dissipated leaving mostly secondary dust grains. The dynamical perturbation of planets in the dust distribution can be revealed with high-contrast imaging in a variety of structures. SPHERE, the high-contrast imaging device installed at the VLT, was designed to search for young giant planets in long period, but is also able to resolve fine details of planetary systems at the scale of astronomical units in the scattered-light regime. As a young and nearby star, NZ Lup was observed in the course of the SPHERE survey. A debris disk had been formerly identified with HST/NICMOS. We observed this system in the near-infrared with the camera in narrow and broad band filters and with the integral field spectrograph. High contrasts are achieved by the mean of pupil tracking combined with angular differential imaging algorithms. The high angular resolution provided by SPHERE allows us to reveal a new feature in the disk which is interpreted as a superimposition of two belts of planetesimals located at stellocentric distances of $\sim$85 and $\sim$115\,au, and with a mutual inclination of about 5$\degb$. Despite the very high inclination of the disk with respect to the line of sight, we conclude that the presence of a gap, that is, a void in the dust distribution between the belts, is likely. We discuss the implication of the existence of two belts and their relative inclination with respect to the presence of planets.

The Cherenkov Telescope Array (CTA) is the major next-generation observatory for ground-based very-high-energy gamma-ray astronomy. It will improve the sensitivity of current ground-based instruments by a factor of five to twenty, depending on the energy, greatly improving both their angular and energy resolutions over four decades in energy (from 20 GeV to 300 TeV). This achievement will be possible by using tens of imaging Cherenkov telescopes of three successive sizes. They will be arranged into two arrays, one per hemisphere, located on the La Palma island (Spain) and in Paranal (Chile). We present here the optimised and final telescope arrays for both CTA sites, as well as their foreseen performance, resulting from the analysis of three different large-scale Monte Carlo productions.

WAVES is designed to study the growth of structure, mass and energy on scales of ~1 kpc to ~10 Mpc over a 7 Gyr timeline. On the largest length scales (1-10 Mpc) WAVES will measure the structures defined by groups, filaments and voids, and their emergence over recent times. Comparisons with bespoke numerical simulations will be used to confirm, refine or refute the Cold Dark Matter paradigm. At intermediate length scales (10 kpc-1 Mpc) WAVES will probe the size and mass distribution of galaxy groups, as well as the galaxy merger rates, in order to directly measure the assembly of dark matter halos and stellar mass. On the smallest length scales (1-10 kpc) WAVES will provide accurate distance and environmental measurements to complement high-resolution space-based imaging to study the mass and size evolution of galaxy bulges, discs and bars. In total, WAVES will provide a panchromatic legacy dataset of ~1.6 million galaxies, firmly linking the very low ($z < 0.1$) and intermediate ($z \sim 0.8$) redshift Universe.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

Context. The 51 Eridani system harbors a complex architecture with its primary star forming a hierarchical system with the binary GJ 3305AB at a projected separation of 2000 au, a giant planet orbiting the primary star at 13 au, and a low-mass debris disk around the primary star with possibly a cold component and a warm component inferred from the spectral energy distribution. Aims. We aim to better constrain the orbital parameters of the known giant planet. Methods. We monitored the system over three years from 2015 to 2018 with the VLT/SPHERE exoplanet imaging instrument. Results. We measure an orbital motion for the planet of ~130 mas with a slightly decreasing separation (~10 mas) and find a hint of curvature. This potential curvature is further supported at 3$\sigma$ significance when including literature GPI astrometry corrected for calibration systematics. Fits of the SPHERE and GPI data using three complementary approaches provide broadly similar results. The data suggest an orbital period of 32$^{+17}_{-9}$ yr (i.e. 12$^{+4}_{-2}$ au in semi-major axis), an inclination of 133$^{+14}_{-7}$ deg, an eccentricity of 0.45$^{+0.10}_{-0.15}$, and an argument of periastron passage of 87$^{+34}_{-30}$ deg [mod 180 deg]. The time at periastron passage and the longitude of node exhibit bimodal distributions because we do not detect yet if the planet is accelerating or decelerating along its orbit. Given the inclinations of the planet's orbit and of the stellar rotation axis (134-144 deg), we infer alignment or misalignment within 18 deg for the star-planet spin-orbit. Further astrometric monitoring in the next 3-4 years is required to confirm at a higher significance the curvature in the planet's motion, determine if the planet is accelerating or decelerating on its orbit, and further constrain its orbital parameters and the star-planet spin-orbit.

Future instruments like NIRCam and MIRI on JWST or METIS at the ELT will be able to image exoplanets that are too faint for current direct imaging instruments. Evolutionary models predicting the planetary intrinsic luminosity as a function of time have traditionally concentrated on gas-dominated giant planets. We extend these cooling curves to Saturnian and Neptunian planets. We simulate the cooling of isolated core-dominated and gas giant planets with masses of 5 Earthmasses to 2 Jupitermasses. The luminosity includes the contribution from the cooling and contraction of the core and of the H/He envelope, as well as radiogenic decay. For the atmosphere we use grey, AMES-Cond, petitCODE, and HELIOS models. We consider solar and non-solar metallicities as well as cloud-free and cloudy atmospheres. The most important initial conditions, namely the core-to-envelope ratio and the initial luminosity are taken from planet formation simulations based on the core accretion paradigm. We first compare our cooling curves for Uranus, Neptune, Jupiter, Saturn, GJ 436b, and a 5 Earthmass-planet with a 1% H/He envelope with other evolutionary models. We then present the temporal evolution of planets with masses between 5 Earthmasses and 2 Jupitermasses in terms of their luminosity, effective temperature, radius, and entropy. We discuss the impact of different post formation entropies. For the different atmosphere types and initial conditions magnitudes in various filter bands between 0.9 and 30 micrometer wavelength are provided. Using black body fluxes and non-grey spectra, we estimate the detectability of such planets with JWST. It is found that a 20 (100) Earthmass-planet can be detected with JWST in the background limit up to an age of about 10 (100) Myr with NIRCam and MIRI, respectively.

The flat spectrum radio quasar PKS 1510-089 is a monitored target in many wavelength bands due to its high variability. It was detected as a very-high-energy (VHE) $\gamma$-ray emitter with H.E.S.S. in 2009, and has since been a regular target of VHE observations by the imaging Cherenkov observatories H.E.S.S. and MAGIC. In this paper, we summarize the current state of results focusing on the monitoring effort with H.E.S.S. and the discovery of a particularly strong VHE flare in 2016 with H.E.S.S. and MAGIC. While the source has now been established as a weak, but regular emitter at VHE, no correlation with other energy bands has been established. This is underlined by the 2016 VHE flare, where the detected optical and high-energy $\gamma$-ray counterparts evolve differently than the VHE flux.

The Kepler, K2 and TESS transit surveys are revolutionizing our understanding of planets orbiting close to their host stars and our understanding of exoplanet systems in general, but there remains a gap in our understanding of wide-orbit planets. This gap in our understanding must be filled if we are to understand planet formation and how it affects exoplanet habitability. We summarize current and planned exoplanet detection programs using a variety of methods: microlensing (including WFIRST), radial velocities, Gaia astrometry, and direct imaging. Finally, we discuss the prospects for joint analyses using results from multiple methods and obstacles that could hinder such analyses. We endorse the findings and recommendations published in the 2018 National Academy report on Exoplanet Science Strategy. This white paper extends and complements the material presented therein.

Surveys reveal that terrestrial- to Neptune-sized planets (1 $< R <$ 4 R$_{\rm{Earth}}$) are the most common type of planets in our galaxy. Detecting and characterizing such small planets around nearby stars holds the key to understanding the diversity of exoplanets and will ultimately address the ubiquitousness of life in the universe. The following fundamental questions will drive research in the next decade and beyond: (1) how common are terrestrial to Neptune-sized planets within a few AU of their host star, as a function of stellar mass? (2) How does planet composition depend on planet mass, orbital radius, and host star properties? (3) What are the energy budgets, atmospheric dynamics, and climates of the nearest worlds? Addressing these questions requires: a) diffraction-limited spatial resolution; b) stability and achievable contrast delivered by adaptive optics; and c) the light-gathering power of extremely large telescopes (ELTs), as well as multi-wavelength observations and all-sky coverage enabled by a comprehensive US ELT Program. Here we provide an overview of the challenge, and promise of success, in detecting and comprehensively characterizing small worlds around the very nearest stars to the Sun with ELTs. This white paper extends and complements the material presented in the findings and recommendations published in the National Academy reports on Exoplanet Science Strategy and Astrobiology Strategy for the Search for Life in the Universe.

Substellar multiplicity is a key outcome of the formation process. The biggest challenge for the next decade will be to distinguish between the formation history, environmental conditions, and dynamical evolution leading to the least massive brown dwarfs and the most massive planets at the tail ends of their mass functions. In this white paper, we advocate for a comprehensive characterization of both the statistical distributions of the population of ultracool dwarf multiple systems and the fundamental properties of their individual components as a function of age. A space-based precision astrometry mission in near-infrared wavelengths would provide the necessary measurements to identify and characterize age-calibrated populations of multiple systems.

The era of precision cosmology has revealed that about 85% of the matter in the universe is dark matter. Two well-motivated candidates are weakly interacting massive particles (WIMPs) and weakly interacting sub-eV particles (WISPs) (e.g. axions). Both WIMPs and WISPs possess distinct {\gamma}-ray signatures. Over the last decade, data taken between 50 MeV to >300 GeV by the Fermi Large Area Telescope (Fermi-LAT) have provided stringent constraints on both classes of dark matter models. Thus far, there are no conclusive detections. However, there is an intriguing {\gamma}-ray excess associated with the Galactic center that could be explained by WIMP annihilation. At lower energies, the poor angular resolution of the Fermi-LAT makes source identification challenging, inhibiting our ability to more sensitively probe both the Galactic center excess, as well as lower-mass WIMP and WISP models. Additionally, targeted WISP searches (e.g., those probing supernovae and blazars) would greatly benefit from enhanced energy resolution and polarization measurements in the MeV range. To address these issues, a new telescope that is optimized for MeV observations is needed. Such an instrument would allow us to explore new areas of dark matter parameter space and provide unprecedented access to its particle nature.

Gas giants are the most readily detectable exoplanets but fundamental questions about their system architectures, formation, migration, and atmospheres have been unanswerable with the current generation of ground- and space-based facilities. The dominant techniques to detect and characterize giant planets $-$ radial velocities, transits, direct imaging, microlensing, and astrometry $-$ are each isolated to a limited range of planet masses, separations, ages, and temperatures. These windows into the arrangement and physical properties of giant planets have spawned new questions about the timescale and location of their assembly; the distributions of planet mass and orbital separation at young and old ages; the composition and structure of their atmospheres; and their orbital and rotational angular momentum architectures. The ELTs will address these questions by building bridges between these islands of mass, orbital distance, and age. The angular resolution, collecting area, all-sky coverage, and novel instrumentation suite of these facilities are needed to provide a complete map of the orbits and atmospheric evolution of gas giant planets (0.3$-$10 $M_\mathrm{Jup}$) across space (0.1$-$100 AU) and time (1 Myr to 10 Gyr). This white paper highlights the scientific potential of the GMT and TMT to address these outstanding questions, with a particular focus on the role of direct imaging and spectroscopy of large samples of giant planets that will soon be made available with $Gaia$.

Having discovered that Earth-sized planets are common, we are now embarking on a journey to determine if Earth-like planets are also common. Finding Earth-like planets is one of the most compelling endeavors of the 21st century - leading us toward finally answering the question: Are we alone? To achieve this forward-looking goal, we must determine the masses of the planets; the sizes of the planets, by themselves, are not sufficient for the determination of the bulk and atmospheric compositions. Masses, coupled with the radii, are crucial constraints on the bulk composition and interior structure of the planets and the composition of their atmospheres, including the search for biosignatures. Precision radial velocity is the most viable technique for providing essential mass and orbit information for spectroscopy of other Earths. The development of high quality precision radial velocity instruments coupled to the building of the large telescope facilities like TMT and GMT or space-based platforms like EarthFinder can enable very high spectral resolution observations with extremely precise radial velocities on minute timescales to allow for the modeling and removal of radial velocity jitter. Over the next decade, the legacy of exoplanet astrophysics can be cemented firmly as part of humankind's quest in finding the next Earth - but only if we can measure the masses and orbits of Earth-sized planets in habitable zone orbits around Sun-like stars.

The possible association of the blazar TXS 0506+056 with a high-energy neutrino detected by IceCube holds the tantalizing potential to answer three astrophysical questions: 1. Where do high-energy neutrinos originate? 2. Where are cosmic rays produced and accelerated? 3. What radiation mechanisms produce the high-energy {\gamma}-rays in blazars? The MeV gamma-ray band holds the key to these questions, because it is an excellent proxy for photo-hadronic processes in blazar jets, which also produce neutrino counterparts. Variability in MeV gamma-rays sheds light on the physical conditions and mechanisms that take place in the particle acceleration sites in blazar jets. In addition, hadronic blazar models also predict a high level of polarization fraction in the MeV band, which can unambiguously distinguish the radiation mechanism. Future MeV missions with a large field of view, high sensitivity, and polarization capabilities will play a central role in multi-messenger astronomy, since pointed, high-resolution telescopes will follow neutrino alerts only when triggered by an all-sky instrument.

The processes that transform gas and dust in circumstellar disks into diverse exoplanets remain poorly understood. One key pathway is to study exoplanets as they form in their young ($\sim$few~Myr) natal disks. Extremely Large Telescopes (ELTs) such as GMT, TMT, or ELT, can be used to establish the initial chemical conditions, locations, and timescales of planet formation, via (1)~measuring the physical and chemical conditions in protoplanetary disks using infrared spectroscopy and (2)~studying planet-disk interactions using imaging and spectro-astrometry. Our current knowledge is based on a limited sample of targets, representing the brightest, most extreme cases, and thus almost certainly represents an incomplete understanding. ELTs will play a transformational role in this arena, thanks to the high spatial and spectral resolution data they will deliver. We recommend a key science program to conduct a volume-limited survey of high-resolution spectroscopy and high-contrast imaging of the nearest protoplanetary disks that would result in an unbiased, holistic picture of planet formation as it occurs.

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

Ground-based telescopes coupled with adaptive optics (AO) have been playing a leading role in exoplanet direct imaging science and technological development for the past two decades and will continue to have an indispensable role for the next decade and beyond. Over the next decade, extreme AO systems on 8-10m telescopes will 1) mitigate risk for WFIRST-CGI by identifying numerous planets the mission can spectrally characterize, 2) validate performance requirements and motivate improvements to atmosphere models needed to unambiguously characterize solar system-analogues from space, and 3) mature novel technological innovations useful for space. Extremely Large Telescopes can deliver the first thermal infrared (10 $\mu m$) images of rocky planets around Sun-like stars and identify biomarkers. These data provide a future NASA direct imaging flagship mission (i.e. HabEx, LUVOIR) with numerous exo-Earth candidates and critical ancillary information to help clarify whether these planets are habitable.

We present the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept selected for study by NASA. It combines huge collecting area, high throughput, broad energy coverage, and excellent spectral and temporal resolution in a single facility. STROBE-X offers an enormous increase in sensitivity for X-ray spectral timing, extending these techniques to extragalactic targets for the first time. It is also an agile mission capable of rapid response to transient events, making it an essential X-ray partner facility in the era of time-domain, multi-wavelength, and multi-messenger astronomy. Optimized for study of the most extreme conditions found in the Universe, its key science objectives include: (1) Robustly measuring mass and spin and mapping inner accretion flows across the black hole mass spectrum, from compact stars to intermediate-mass objects to active galactic nuclei. (2) Mapping out the full mass-radius relation of neutron stars using an ensemble of nearly two dozen rotation-powered pulsars and accreting neutron stars, and hence measuring the equation of state for ultradense matter over a much wider range of densities than explored by NICER. (3) Identifying and studying X-ray counterparts (in the post-Swift era) for multiwavelength and multi-messenger transients in the dynamic sky through cross-correlation with gravitational wave interferometers, neutrino observatories, and high-cadence time-domain surveys in other electromagnetic bands. (4) Continuously surveying the dynamic X-ray sky with a large duty cycle and high time resolution to characterize the behavior of X-ray sources over an unprecedentedly vast range of time scales. STROBE-X's formidable capabilities will also enable a broad portfolio of additional science.

Aims: We want to detect and quantify observables related to accretion processes occurring locally in circumstellar disks, which could be attributed to young forming planets. We focus on objects known to host protoplanet candidates and/or disk structures thought to be the result of interactions with planets. Methods: We analyzed observations of 6 young stars (age $3.5-10$ Myr) and their surrounding environments with the SPHERE/ZIMPOL instrument on the VLT in the H$\alpha$ filter (656 nm) and a nearby continuum filter (644.9 nm). Results: We re-detect the known accreting M-star companion HD142527 B with the highest published signal to noise to date in both H$\alpha$ and the continuum. We derive new astrometry ($r = 62.8^{+2.1}_{-2.7}$ mas and $\text{PA} = (98.7\,\pm1.8)^\circ$) and photometry ($\Delta$N_Ha=$6.3^{+0.2}_{-0.3}$ mag, $\Delta$B_Ha=$6.7\pm0.2$ mag and $\Delta$Cnt_Ha=$7.3^{+0.3}_{-0.2}$ mag) for the companion in agreement with previous studies, and estimate its mass accretion rate ($\dot{M}\approx1-2\,\times10^{-10}\,M_\odot\text{ yr}^{-1}$). A faint point-like source around HD135344 B (SAO206462) is also investigated, but a second deeper observation is required to reveal its nature. No other companions are detected. In the framework of our assumptions we estimate detection limits at the locations of companion candidates around HD100546, HD169142 and MWC758 and calculate that processes involving H$\alpha$ fluxes larger than $\sim8\times10^{-14}-10^{-15}\,\text{erg/s/cm}^2$ ($\dot{M}>10^{-10}-10^{-12}\,M_\odot\text{ yr}^{-1}$) can be excluded. Furthermore, flux upper limits of $\sim10^{-14}-10^{-15}\,\text{erg/s/cm}^2$ ($\dot{M}<10^{-11}-10^{-12}\,M_\odot \text{ yr}^{-1}$) are estimated within the gaps identified in the disks surrounding HD135344B and TW Hya.

In contrast to the water-poor inner solar system planets, stochasticity during planetary formation and order of magnitude deviations in exoplanet volatile contents suggest that rocky worlds engulfed in thick volatile ice layers are the dominant family of terrestrial analogues among the extrasolar planet population. However, the distribution of compositionally Earth-like planets remains insufficiently constrained, and it is not clear whether the solar system is a statistical outlier or can be explained by more general planetary formation processes. Here we employ numerical models of planet formation, evolution, and interior structure, to show that a planet's bulk water fraction and radius are anti-correlated with initial $^{26}$Al levels in the planetesimal-based accretion framework. The heat generated by this short-lived radionuclide rapidly dehydrates planetesimals prior to accretion onto larger protoplanets and yields a system-wide correlation of planet bulk abundances, which, for instance, can explain the lack of a clear orbital trend in the water budgets of the TRAPPIST-1 planets. Qualitatively, our models suggest two main scenarios of planetary systems' formation: high-$^{26}$Al systems, like our solar system, form small, water-depleted planets, whereas those devoid of $^{26}$Al predominantly form ocean worlds, where the mean planet radii between both scenarios deviate by up to about 10%.

Aims. R Coronae Australis (R CrA) is the brightest star of the Coronet nebula of the Corona Australis (CrA) star forming region. It has very red colors, probably due to dust absorption and it is strongly variable. High contrast instruments allow for an unprecedented direct exploration of the immediate circumstellar environment of this star. Methods. We observed R CrA with the near-IR channels (IFS and IRDIS) of SPHERE at VLT. In this paper, we used four different epochs, three of them from open time observations while one is from the SPHERE guaranteed time. The data were reduced using the DRH pipeline and the SPHERE Data Center. On the reduced data we implemented custom IDL routines with the aim to subtract the speckle halo.We have also obtained pupil-tracking H-band (1.45-1.85 micron) observations with the VLT/SINFONI near-infrared medium-resolution (R~3000) spectrograph. Results. A companion was found at a separation of 0.156" from the star in the first epoch and increasing to 0.18400 in the final one. Furthermore, several extended structures were found around the star, the most noteworthy of which is a very bright jet-like structure North-East from the star. The astrometric measurements of the companion in the four epochs confirm that it is gravitationally bound to the star. The SPHERE photometry and the SINFONI spectrum, once corrected for extinction, point toward an early M spectral type object with a mass between 0.3 and 0.55 M?. The astrometric analyis provides constraints on the orbit paramenters: e~0.4, semi-major axis at 27-28 au, inclination of ~ 70{\deg} and a period larger than 30 years. We were also able to put constraints of few MJup on the mass of possible other companions down to separations of few tens of au.

Young planets are expected to cause perturbations in protostellar disks that may be used to infer their presence. Clear detection of still-forming planets embedded within gas-rich disks is rare. HD 169142 is a very young Herbig Ae-Be star surrounded by a pre-transitional disk, composed of at least three rings. While claims of sub-stellar objects around this star have been made previously, follow-up studies remain inconclusive. We used SPHERE at ESO VLT to obtain a sequence of high-contrast images of the immediate surroundings of this star over about three years. This enables a photometric and astrometric analysis of the structures in the disk. While we were unable to definitively confirm the previous claims of a massive sub-stellar object at 0.1-0.15 arcsec from the star, we found both spirals and blobs within the disk. The spiral pattern may be explained as due to the presence of a primary, a secondary, and a tertiary arm excited by a planet of a few Jupiter masses lying along the primary arm, likely in the cavities between the rings. The blobs orbit the star consistently with Keplerian motion, allowing a dynamical determination of the mass of the star. While most of these blobs are located within the rings, we found that one of them lies in the cavity between the rings, along the primary arm of the spiral design. This blob might be due to a planet that might also be responsible for the spiral pattern observed within the rings and for the cavity between the two rings. The planet itself is not detected at short wavelengths, where we only see a dust cloud illuminated by stellar light, but the planetary photosphere might be responsible for the emission observed in the K band. The mass of this putative planet may be constrained using photometric and dynamical arguments; it should be between 1 and 4 Jupiter masses. The brightest blobs are found at the 1:2 resonance with this putative planet

We detect the accreting low-mass companion HD142527B at a separation of 73 mas (11.4 au) from the star. No other companions with mass greater than 10 MJ are visible in the field of view of IFS (\sim 100 au centered on the star) or in the IRDIS field of view (\sim 400 au centered on the star). Measurements from IFS, SAM IFS, and IRDIS suggest an M6 spectral type for HD142527B, with an uncertainty of one spectral subtype, compatible with an object of M=0.11 \pm 0.06 MSun and R=0.15 \pm 0.07 RSun. The determination of the mass remains a challenge using contemporary evolutionary models, as they do not account for the energy input due to accretion from infalling material. We consider that the spectral type of the secondary may also be earlier than the type we derived from IFS spectra. From dynamical considerations, we further constrain the mass to 0.26^{+0.16}_{-0.14} MSun , which is consistent with both our spectroscopic analysis and the values reported in the literature. Following previous methods, the lower and upper dynamical mass values correspond to a spectral type between M2.5 and M5.5 for the companion. By fitting the astrometric points, we find the following orbital parameters: a period of P=35-137 yr; an inclination of i=121-130 deg.; , a value of Omega=124-135 deg for the longitude of node, and an 68% confidence interval of \sim 18 - 57 au for the separation at periapsis. Eccentricity and time at periapsis passage exhibit two groups of values: \sim0.2-0.45 and \sim0.45-0.7 for e, and \sim 2015-2020 and \sim2020-2022 for T_0. While these orbital parameters might at first suggest that HD142527B is not the companion responsible for the outer disk truncation, a previous hydrodynamical analysis of this system showed that they are compatible with a companion that is able to produce the large cavity and other observed features.

We present an analysis of the integrated neutral hydrogen (HI) properties for 27 galaxies within nine low mass, gas-rich, late-type dominated groups which we denote "Choirs". We find that majority of the central Choir galaxies have average HI content: they have a normal gas-mass fraction with respect to isolated galaxies of the same stellar mass. In contrast, we find more satellite galaxies with a lower gas-mass fraction than isolated galaxies of the same stellar mass. A likely reason for the lower gas content in these galaxies is tidal stripping. Both the specific star formation rate and the star formation efficiency of the central group galaxies are similar to galaxies in isolation. The Choir satellite galaxies have similar specific star formation rate as galaxies in isolation, therefore satellites that exhibit a higher star formation efficiency simply owe it to their lower gas-mass fractions. We find that the most HI massive galaxies have the largest HI discs and fall neatly onto the HI size-mass relation, while outliers are galaxies that are experiencing interactions. We find that high specific angular momentum could be a reason for galaxies to retain the large fraction of HI gas in their discs. This shows that for the Choir groups with no evidence of interactions, as well as those with traces of minor mergers, the internal galaxy properties dominate over the effects of residing in a group. The probed galaxy properties strengthen evidence that the Choir groups represent the early stages of group assembly.

HIP 65426 b is a recently discovered exoplanet imaged during the course of the SPHERE-SHINE survey. Here we present new $L'$ and $M'$ observations of the planet from the NACO instrument at the VLT from the NACO-ISPY survey, as well as a new $Y-H$ spectrum and $K$-band photometry from SPHERE-SHINE. Using these data, we confirm the nature of the companion as a warm, dusty planet with a mid-L spectral type. From comparison of its SED with the BT-Settl atmospheric models, we derive a best-fit effective temperature of $T_{\text{eff}}=1618\pm7$ K, surface gravity $\log g=3.78^{+0.04}_{-0.03}$ and radius $R=1.17\pm0.04$ $R_{\text{J}}$ (statistical uncertainties only). Using the DUSTY and COND isochrones we estimate a mass of $8\pm1$ $M_{\text{J}}$. Combining the astrometric measurements from our new datasets and from the literature, we show the first indications of orbital motion of the companion (2.6$\sigma$ significance) and derive preliminary orbital constraints. We find a highly inclined orbit ($i=107^{+13}_{-10}$ deg) with an orbital period of $800^{+1200}_{-400}$ yr. We also report SPHERE sparse aperture masking observations that investigate the possibility that HIP 65426 b was scattered onto its current orbit by an additional companion at a smaller orbital separation. From this data we rule out the presence of brown dwarf companions with masses greater than 16 $M_{\text{J}}$ at separations larger than 3 AU, significantly narrowing the parameter space for such a companion.

With an orbital distance comparable to that of Saturn in the solar system, \bpic b is the closest (semi-major axis $\simeq$\,9\,au) exoplanet that has been imaged to orbit a star. Thus it offers unique opportunities for detailed studies of its orbital, physical, and atmospheric properties, and of disk-planet interactions. With the exception of the discovery observations in 2003 with NaCo at the Very Large Telescope (VLT), all following astrometric measurements relative to \bpic have been obtained in the southwestern part of the orbit, which severely limits the determination of the planet's orbital parameters. We aimed at further constraining \bpic b orbital properties using more data, and, in particular, data taken in the northeastern part of the orbit. We used SPHERE at the VLT to precisely monitor the orbital motion of beta \bpic b since first light of the instrument in 2014. We were able to monitor the planet until November 2016, when its angular separation became too small (125 mas, i.e., 1.6\,au) and prevented further detection. We redetected \bpic b on the northeast side of the disk at a separation of 139\,mas and a PA of 30$^{\circ}$ in September 2018. The planetary orbit is now well constrained. With a semi-major axis (sma) of $a = 9.0 \pm 0.5$ au (1 $\sigma $), it definitely excludes previously reported possible long orbital periods, and excludes \bpic b as the origin of photometric variations that took place in 1981. We also refine the eccentricity and inclination of the planet. From an instrumental point of view, these data demonstrate that it is possible to detect, if they exist, young massive Jupiters that orbit at less than 2 au from a star that is 20 pc away.

We present a study of 16 HI-detected galaxies found in 178 hours of observations from Epoch 1 of the COSMOS HI Large Extragalactic Survey (CHILES). We focus on two redshift ranges between 0.108 <= z <= 0.127 and 0.162 <= z <= 0.183 which are among the worst affected by radio frequency interference (RFI). While this represents only 10% of the total frequency coverage and 18% of the total expected time on source compared to what will be the full CHILES survey, we demonstrate that our data reduction pipeline recovers high quality data even in regions severely impacted by RFI. We report on our in-depth testing of an automated spectral line source finder to produce HI total intensity maps which we present side-by-side with significance maps to evaluate the reliability of the morphology recovered by the source finder. We recommend that this become a common place manner of presenting data from upcoming HI surveys of resolved objects. We use the COSMOS 20k group catalogue, and we extract filamentary structure using the topological DisPerSE algorithm to evaluate the \hi\ morphology in the context of both local and large-scale environments and we discuss the shortcomings of both methods. Many of the detections show disturbed HI morphologies suggesting they have undergone a recent interaction which is not evident from deep optical imaging alone. Overall, the sample showcases the broad range of ways in which galaxies interact with their environment. This is a first look at the population of galaxies and their local and large-scale environments observed in HI by CHILES at redshifts beyond the z=0.1 Universe.

We present long baseline Atacama Large Millimeter/submillimeter Array (ALMA) observations of the 870$\,\mu$m dust continuum emission and CO (3-2) from the protoplanetary disk around the Herbig Ae/Be star HD 100546, which is one of the few systems claimed to have two young embedded planets. These observations achieve a resolution of 4 au (3.8 mas), an rms noise of 66$\mu$Jy/beam, and reveal an asymmetric ring between $\sim$20-40 au with largely optically thin dust continuum emission. This ring is well fit by two concentric and overlapping Gaussian rings of different widths and a Vortex. In addition, an unresolved component is detected at a position consistent with the central star, which may trace the central inner disk ($<$2au in radius). We report a lack of compact continuum emission at the positions of both claimed protoplanets. We use this result to constrain the circumplanetary disk (CPD) mass and size of 1.44M$_{\rm Earth}$ and 0.44au in the optically thin and thick regime, respectively, for the case of the previously directly imaged protoplanet candidate at $\sim$55 au (HD100546 b). We compare these empirical CPD constraints to previous numerical simulations. This suggests that HD100546 b is inconsistent with several planet accretion models, while gas-starved models are also still compatible. We estimate the planetary mass as 1.65 M$_J$ by using the relation between planet, circumstellar, and circumplanetary masses derived from numerical simulations. Finally, the CO integrated intensity map shows a possible spiral arm feature that would match the spiral features identified in Near-Infrared scattered light polarized emission, which suggests a real spiral feature in the disk surface that needs to be confirmed with further observations.

We present Widefield ASKAP L-band Legacy All-sky Blind Survey (WALLABY) early science results from the Australian Square Kilometre Array Pathfinder (ASKAP) observations of the NGC 7162 galaxy group. We use archival HIPASS and Australia Telescope Compact Array (ATCA) observations of this group to validate the new ASKAP data and the data reduction pipeline ASKAPsoft. We detect six galaxies in the neutral hydrogen (HI) 21-cm line, expanding the NGC 7162 group membership from four to seven galaxies. Two of the new detections are also the first HI detections of the dwarf galaxies, AM 2159-434 and GALEXASC J220338.65-431128.7, for which we have measured velocities of $cz=2558$ and $cz=2727$ km s$^{-1}$, respectively. We confirm that there is extended HI emission around NGC 7162 possibly due to past interactions in the group as indicated by the $40^{\circ}$ offset between the kinematic and morphological major axes for NGC 7162A, and its HI richness. Taking advantage of the increased resolution (factor of $\sim1.5$) of the ASKAP data over archival ATCA observations, we fit a tilted ring model and use envelope tracing to determine the galaxies' rotation curves. Using these we estimate the dynamical masses and find, as expected, high dark matter fractions of $f_{\mathrm{DM}}\sim0.81-0.95$ for all group members. The ASKAP data are publicly available.

In this paper, we address two issues related to primordial disk evolution in three clusters (NGC 1333, IC 348, and Orion A) observed by the INfrared Spectra of Young Nebulous Clusters (IN-SYNC) project. First, in each cluster, averaged over the spread of age, we investigate how disk lifetime is dependent on stellar mass. The general relation in IC 348 and Orion A is that primordial disks around intermediate mass stars (2--5$M_{\odot}$) evolve faster than those around loss mass stars (0.1--1$M_{\odot}$), which is consistent with previous results. However, considering only low mass stars, we do not find a significant dependence of disk frequency on stellar mass. These results can help to better constrain theories on gas giant planet formation timescales. Secondly, in the Orion A molecular cloud, in the mass range of 0.35--0.7$M_{\odot}$, we provide the most robust evidence to date for disk evolution within a single cluster exhibiting modest age spread. By using surface gravity as an age indicator and employing 4.5 $\mu m$ excess as a primordial disk diagnostic, we observe a trend of decreasing disk frequency for older stars. The detection of intra-cluster disk evolution in NGC 1333 and IC 348 is tentative, since the slight decrease of disk frequency for older stars is a less than 1-$\sigma$ effect.

The nearby Herbig Be star HD100546 is known to be a laboratory for the study of protoplanets and their relation with the circumstellar disk that is carved by at least 2 gaps. We observed the HD100546 environment with high contrast imaging exploiting several different observing modes of SPHERE, including datasets with/without coronagraphs, dual band imaging, integral field spectroscopy and polarimetry. The picture emerging from these different data sets is complex. Flux-conservative algorithms images clearly show the disk up to 200au. More aggressive algorithms reveal several rings and warped arms overlapping the main disk. The bright parts of this ring lie at considerable height over the disk mid-plane at about 30au. Our images demonstrate that the brightest wings close to the star in the near side of the disk are a unique structure, corresponding to the outer edge of the intermediate disk at ~40au. Modeling of the scattered light from the disk with a geometrical algorithm reveals that a moderately thin structure can well reproduce the light distribution in the flux-conservative images. We suggest that the gap between 44 and 113 au span between the 1:2 and 3:2 resonance orbits of a massive body located at ~70au that might coincide with the candidate planet HD100546b detected with previous thermal IR observations. In this picture, the two wings can be the near side of a ring formed by disk material brought out of the disk at the 1:2 resonance with the same massive object. While we find no clear evidence confirming detection of the planet candidate HD100546c in our data, we find a diffuse emission close to the expected position of HD100546b. This source can be described as an extremely reddened substellar object surrounded by a dust cloud or its circumplanetary disk. Its astrometry is broadly consistent with a circular orbital motion on the disk plane.

BL Lac Objects (BL Lacs) and Flat Spectrum Radio Quasars (FSRQs) are radio-loud active galaxies (AGNs) whose jets are seen at a small viewing angle (blazars), while Misaligned Active Galactic Nuclei (MAGNs) are mainly radiogalaxies of type FRI or FRII and Steep Spectrum Radio Quasars (SSRQs), which show jets of radiation oriented away from the observer's line of sight. MAGNs are very numerous and well studied in the lower energies of the electromagnetic spectrum but are not commonly observed in the gamma-ray energy range, because their inclination leads to the loss of relativistic boosting of the jet emission. The Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope in the 100 MeV -300 GeV energy range detected only 18 MAGNs (15 radio galaxies and 3 SSRQs) compared to 1144 blazars. Studying MAGNs and their environment in the gamma-ray sky is extremely interesting, because FRI and FRII radio galaxies are respectively considered the parent populations of BL Lacs and FSRQs, and these account for more than 50% of the known gamma-ray sources. The aim of this study is to hunt new gamma-ray MAGN candidates among the remaining blazars of uncertain type and unassociated AGNs, using machine learning techniques and other physical constraints when strict classifications are not available. We found 10 new MAGN candidates associated with gamma-ray sources. Their features are consistent with a source with a misaligned jet of radiation. This study reinforces the need for more systematic investigation of MAGNs in order to improve understanding of the radiation emission mechanisms and and the disparity of detection between more powerful and weaker gamma-ray AGNs.

Context. A low-mass brown dwarf has been recently imaged around HR 2562 (HD 50571), a star hosting a debris disk resolved in the far infrared. Interestingly, the companion location is compatible with an orbit coplanar with the disk and interior to the debris belt. This feature makes the system a valuable laboratory to analyze the formation of substellar companions in a circumstellar disk and potential disk-companion dynamical interactions. Aims. We aim to further characterize the orbital motion of HR 2562 B and its interactions with the host star debris disk. Methods. We performed a monitoring of the system over ~10 months in 2016 and 2017 with the VLT/SPHERE exoplanet imager. Results. We confirm that the companion is comoving with the star and detect for the first time an orbital motion at high significance, with a current orbital motion projected in the plane of the sky of 25 mas (~0.85 au) per year. No orbital curvature is seen in the measurements. An orbital fit of the SPHERE and literature astrometry of the companion without priors on the orbital plane clearly indicates that its orbit is (quasi-)coplanar with the disk. To further constrain the other orbital parameters, we used empirical laws for a companion chaotic zone validated by N-body simulations to test the orbital solutions that are compatible with the estimated disk cavity size. Non-zero eccentricities (>0.15) are allowed for orbital periods shorter than 100 yr, while only moderate eccentricities up to ~0.3 for orbital periods longer than 200 yr are compatible with the disk observations. A comparison of synthetic Herschel images to the real data does not allow us to constrain the upper eccentricity of the companion.

Binary stars constitute a large percentage of the stellar population, yet relatively little is known about the planetary systems orbiting them. Most constraints on circumbinary planets (CBPs) so far come from transit observations with the Kepler telescope, which is sensitive to close-in exoplanets but does not constrain planets on wider orbits. However, with continuous developments in high-contrast imaging techniques, this population can now be addressed through direct imaging. We present the full survey results of the Search for Planets Orbiting Two Stars (SPOTS) survey, which is the first direct imaging survey targeting CBPs. The SPOTS observational program comprises 62 tight binaries that are young and nearby, and thus suitable for direct imaging studies, with VLT/NaCo and VLT/SPHERE. Results from SPOTS include the resolved circumbinary disk around AK Sco, the discovery of a low-mass stellar companion in a triple packed system, the relative astrometry of up to 9 resolved binaries, and possible indications of non-background planetary-mass candidates around HIP 77911. We did not find any CBP within 300 AU, which implies a frequency upper limit on CBPs (1--15 $M_{\rm Jup}$) of 6--10 % between 30-300 AU. Coupling these observations with an archival dataset for a total of 163 stellar pairs, we find a best-fit CBP frequency of 1.9 % (2--15 $M_{\rm Jup}$) between 1--300 AU with a 10.5 % upper limit at a 95 % confidence level. This result is consistent with the distribution of companions around single stars.

Centaurus A (Cen A) is the nearest radio galaxy discovered as a very-high-energy (VHE; 100 GeV-100 TeV) $\gamma$-ray source by the High Energy Stereoscopic System (H.E.S.S.). It is a faint VHE $\gamma$-ray emitter, though its VHE flux exceeds both the extrapolation from early Fermi-LAT observations as well as expectations from a (misaligned) single-zone synchrotron-self Compton (SSC) description. The latter satisfactorily reproduces the emission from Cen A at lower energies up to a few GeV. New observations with H.E.S.S., comparable in exposure time to those previously reported, were performed and eight years of Fermi-LAT data were accumulated to clarify the spectral characteristics of the $\gamma$-ray emission from the core of Cen A. The results allow us for the first time to achieve the goal of constructing a representative, contemporaneous $\gamma$-ray core spectrum of Cen A over almost five orders of magnitude in energy. Advanced analysis methods, including the template fitting method, allow detection in the VHE range of the core with a statistical significance of 12$\sigma$ on the basis of 213 hours of total exposure time. The spectrum in the energy range of 250 GeV-6 TeV is compatible with a power-law function with a photon index $\Gamma=2.52\pm0.13_{\mathrm{stat}}\pm0.20_{\mathrm{sys}}$. An updated Fermi-LAT analysis provides evidence for spectral hardening by $\Delta\Gamma\simeq0.4\pm0.1$ at $\gamma$-ray energies above $2.8^{+1.0}_{-0.6}$ GeV at a level of $4.0\sigma$. The fact that the spectrum hardens at GeV energies and extends into the VHE regime disfavour a single-zone SSC interpretation for the overall spectral energy distribution (SED) of the core and is suggestive of a new $\gamma$-ray emitting component connecting the high-energy emission above the break energy to the one observed at VHE energies.

Young circumstellar disks are of prime interest to understand the physical and chemical conditions under which planet formation takes place. Only very few detections of planet candidates within these disks exist, and most of them are currently suspected to be disk features. In this context, the transition disk around the young star PDS 70 is of particular interest, due to its large gap identified in previous observations, indicative of ongoing planet formation. We aim to search for the presence of planets and search for disk structures indicative for disk-planet interactions and other evolutionary processes. We analyse new and archival near-infrared (NIR) images of the transition disk PDS 70 obtained with the VLT/SPHERE, VLT/NaCo and Gemini/NICI instruments in polarimetric differential imaging (PDI) and angular differential imaging (ADI) modes. We detect a point source within the gap of the disk at about 195 mas (about 22 au) projected separation. The detection is confirmed at five different epochs, in three filter bands and using different instruments. The astrometry results in an object of bound nature, with high significance. The comparison of the measured magnitudes and colours to evolutionary tracks suggests that the detection is a companion of planetary mass. We confirm the detection of a large gap of about 54 au in size within the disk in our scattered light images, and detect a signal from an inner disk component. We find that its spatial extent is very likely smaller than about 17 au in radius. The images of the outer disk show evidence of a complex azimuthal brightness distribution which may in part be explained by Rayleigh scattering from very small grains. Future observations of this system at different wavelengths and continuing astrometry will allow us to test theoretical predictions regarding planet-disk interactions, planetary atmospheres and evolutionary models.

The G-type star GJ504A is known to host a 3 to 35 MJup companion whose temperature, mass, and projected separation all contribute to make it a test case for the planet formation theories and for atmospheric models of giant planets and light brown dwarfs. We collected data from the CHARA interferometer, SOPHIE spectrograph, and VLT/SPHERE high contrast imager to revisit the properties of the system. We measure a radius of 1.35+/- 0.04Rsun for GJ504A which yields isochronal ages of 21+/-2Myr or 4.0+/-1.8Gyr for the system and line-of-sight stellar rotation axis inclination of $162.4_{-4.3}^{+3.8}$ degrees or $18.6_{-3.8}^{+4.3}$ degrees. We re-detect the companion in the Y2, Y3, J3, H2, and K1 dual band SPHERE images. The complete 1-4 $\mu$m SED shape of GJ504b is best reproduced by T8-T9.5 objects with intermediate ages ($\leq1.5$Gyr), and/or unusual dusty atmospheres and/or super-solar metallicities. All six atmospheric models used yield $\mathrm{T_{eff}=550 \pm 50}$K for GJ504b and point toward a low surface gravity (3.5-4.0 dex). The accuracy on the metallicity value is limited by model-to-model systematics. It is not degenerate with the C/O ratio. We derive $\mathrm{log\:L/L_{\odot}=-6.15\pm0.15}$ dex for the companion compatible with masses of $\mathrm{M=1.3^{+0.6}_{-0.3}M_{Jup}}$ and $\mathrm{M=23^{+10}_{-9} M_{Jup}}$ for the young and old age ranges, respectively. The semi-major axis (sma) is above 27.8 au and the eccentricity lower than 0.55. The posterior on GJ~504b's orbital inclination suggests a misalignment with GJ~504A rotation axis. We combine the radial velocity and multi-epoch imaging data to exclude additional objects (90\% prob.) more massive than 2.5 and 30 $\mathrm{M_{Jup}}$ with sma in the range 0.01-80 au for the young and old system ages, respectively. The companion is in the envelope of the population of planets synthetized with our core-accretion model.

Aims: We aim to characterize the orbital and atmospheric properties of PDS 70 b, which was first identified on May 2015 in the course of the SHINE survey with SPHERE, the extreme adaptive-optics instrument at the VLT. Methods: We obtained new deep SPHERE/IRDIS imaging and SPHERE/IFS spectroscopic observations of PDS 70 b. The astrometric baseline now covers 6 years which allows us to perform an orbital analysis. For the first time, we present spectrophotometry of the young planet which covers almost the entire near-infrared range (0.96 to 3.8 micrometer). We use different atmospheric models covering a large parameter space in temperature, log(g), chemical composition, and cloud properties to characterize the properties of the atmosphere of PDS 70 b. Results: PDS 70 b is most likely orbiting the star on a circular and disk coplanar orbit at ~22 au inside the gap of the disk. We find a range of models that can describe the spectrophotometric data reasonably well in the temperature range between 1000-1600 K and log(g) no larger than 3.5 dex. The planet radius covers a relatively large range between 1.4 and 3.7 R_jupiter with the larger radii being higher than expected from planet evolution models for the age of the planet of 5.4 Myr. Conclusions: This study provides a comprehensive dataset on the orbital motion of PDS 70 b, indicating a circular orbit and a motion coplanar with the disk. The first detailed spectral energy distribution of PDS 70 b indicates a temperature typical for young giant planets. The detailed atmospheric analysis indicates that a circumplanetary disk may contribute to the total planet flux.

We present observations with the planet finder SPHERE of a selected sample of the most promising radial velocity (RV) companions for high-contrast imaging. Using a Monte Carlo simulation to explore all the possible inclinations of the orbit of wide RV companions, we identified the systems with companions that could potentially be detected with SPHERE. We found the most favorable RV systems to observe are : HD\,142, GJ\,676, HD\,39091, HIP\,70849, and HD\,30177 and carried out observations of these systems during SPHERE Guaranteed Time Observing (GTO). To reduce the intensity of the starlight and reveal faint companions, we used Principle Component Analysis (PCA) algorithms alongside angular and spectral differential imaging. We injected synthetic planets with known flux to evaluate the self-subtraction caused by our data reduction and to determine the 5$\sigma$ contrast in the J band $vs$ separation for our reduced images. We estimated the upper limit on detectable companion mass around the selected stars from the contrast plot obtained from our data reduction. Although our observations enabled contrasts larger than 15 mag at a few tenths of arcsec from the host stars, we detected no planets. However, we were able to set upper mass limits around the stars using AMES-COND evolutionary models. We can exclude the presence of companions more massive than 25-28 \MJup around these stars, confirming the substellar nature of these RV companions.

For the first time both H$\alpha$ and far-ultraviolet (FUV) observations from an HI-selected sample are used to determine the dust-corrected star formation rate density (SFRD: $\dot{\rho}$) in the local Universe. Applying the two star formation rate indicators on 294 local galaxies we determine log($\dot{\rho}$$ _{H\alpha}) = -1.68~^{+0.13}_{-0.05}$ [M$_{\odot} $ yr$^{-1} $ Mpc$^{-3}]$ and log($\dot{\rho}_{FUV}$) $ = -1.71~^{+0.12}_{-0.13}$ [M$_\odot $ yr$^{-1} $ Mpc$^{-3}]$. These values are derived from scaling H$\alpha$ and FUV observations to the HI mass function. Galaxies were selected to uniformly sample the full HI mass (M$_{HI}$) range of the HI Parkes All-Sky Survey (M$_{HI} \sim10^{7}$ to $\sim10^{10.7}$ M$_{\odot}$). The approach leads to relatively larger sampling of dwarf galaxies compared to optically-selected surveys. The low HI mass, low luminosity and low surface brightness galaxy populations have, on average, lower H$\alpha$/FUV flux ratios than the remaining galaxy populations, consistent with the earlier results of Meurer. The near-identical H$\alpha$- and FUV-derived SFRD values arise with the low H$\alpha$/FUV flux ratios of some galaxies being offset by enhanced H$\alpha$ from the brightest and high mass galaxy populations. Our findings confirm the necessity to fully sample the HI mass range for a complete census of local star formation to include lower stellar mass galaxies which dominate the local Universe.

The Deep Extragalactic VIsible Legacy Survey (DEVILS) is a large spectroscopic campaign at the Anglo-Australian Telescope (AAT) aimed at bridging the near and distant Universe by producing the highest completeness survey of galaxies and groups at intermediate redshifts ($0.3<z<1.0$). Our sample consists of $\sim$60,000 galaxies to Y$<$21.2mag, over $\sim$6deg$^{2}$ in three well-studied deep extragalactic fields (Cosmic Origins Survey field, COSMOS, Extended Chandra Deep Field South, ECDFS and the X-ray Multi-Mirror Mission Large-Scale Structure region, XMM-LSS - all Large Synoptic Survey Telescope deep-drill fields). This paper presents the broad experimental design of DEVILS. Our target sample has been selected from deep Visible and Infrared Survey Telescope for Astronomy (VISTA) Y-band imaging (VISTA Deep Extragalactic Observations, VIDEO and UltraVISTA), with photometry measured by ProFound. Photometric star/galaxy separation is done on the basis of NIR colours, and has been validated by visual inspection. To maximise our observing efficiency for faint targets we employ a redshift feedback strategy, which continually updates our target lists, feeding back the results from the previous night's observations. We also present an overview of the initial spectroscopic observations undertaken in late 2017 and early 2018.

Evolutionary models are widely used to infer the mass of stars, brown dwarfs, and giant planets. Their predictions are thought to be less reliable at young ages ($<$ 200 Myr) and in the low-mass regime ($\mathrm{<1~M_{\odot}}$). GJ 2060 AB and TWA 22 AB are two rare astrometric M-dwarf binaries respectively members of the AB Doradus and Beta Pictoris moving groups. As their dynamical mass can be measured within a few years, they can be used to calibrate the evolutionary tracks and set new constraints on the age of young moving groups. We find a total mass of $\mathrm{0.18\pm 0.02~M_\odot}$ for TWA 22. That mass is in good agreement with model predictions at the age of the Beta Pic moving group. We obtain a total mass of $\mathrm{1.09 \pm 0.10~M_{\odot}}$ for GJ 2060. We estimate a spectral type of M$1\pm0.5$, $\mathrm{L/L_{\odot}=-1.20\pm0.05}$ dex, and $\mathrm{T_{eff}=3700\pm100}$ K for GJ 2060 A. The B component is a M$3\pm0.5$ dwarf with $\mathrm{L/L_{\odot}=-1.63\pm0.05}$ dex and $\mathrm{T_{eff}=3400\pm100}$ K. The dynamical mass of GJ 2060 AB is inconsistent with the most recent models predictions (BCAH15, PARSEC) for an ABDor age in the range 50-150 Myr. It is 10 to 20\% (1-2 sigma, depending on the assumed age) above the models predictions, corresponding to an underestimation of $0.10$ to $0.20~\mathrm{M_\odot}$. Coevality suggests a young age for the system ($\sim$ 50 Myr) according to most evolutionary models. TWA 22 validates the predictions of recent evolutionary tracks at $\sim$20 Myr. On the other hand, we evidence a 1-2 sigma mismatch between the predicted and observed mass of GJ 2060 AB. This slight departure may indicate that one of the star hosts a tight companion. Alternatively, this would confirm the models tendency to underestimate the mass of young low-mass stars.

Spectral lines are among the most powerful signatures for dark matter (DM) annihilation searches in very-high-energy $\gamma$-rays. The central region of the Milky Way halo is one of the most promising targets given its large amount of DM and proximity to Earth. We report on a search for a monoenergetic spectral line from self-annihilations of DM particles in the energy range from 300 GeV to 70 TeV using a two-dimensional maximum likelihood method taking advantage of both the spectral and spatial features of signal versus background. The analysis makes use of Galactic Center (GC) observations accumulated over ten years (2004 - 2014) with the H.E.S.S. array of ground-based Cherenkov telescopes. No significant $\gamma$-ray excess above the background is found. We derive upper limits on the annihilation cross section $\langle\sigma v\rangle$ for monoenergetic DM lines at the level of $\sim4\times10^{-28}$ cm$^{3}$s$^{-1}$ at 1 TeV, assuming an Einasto DM profile for the Milky Way halo. For a DM mass of 1 TeV, they improve over the previous ones by a factor of six. The present constraints are the strongest obtained so far for DM particles in the mass range 300 GeV - 70 TeV. Ground-based $\gamma$-ray observations have reached sufficient sensitivity to explore relevant velocity-averaged cross sections for DM annihilation into two $\gamma$-ray photons at the level expected from the thermal relic density for TeV DM particles.

The study of high contrast imaged brown dwarfs and exoplanets depends strongly on evolutionary models. To estimate the mass of a directly imaged substellar object, its extracted photometry or spectrum is used and adjusted with model spectra together with the estimated age of the system. These models still need to be properly tested and constrained. HD 4747B is a brown dwarf close to the H burning mass limit, orbiting a nearby, solar-type star and has been observed with the radial velocity method over almost two decades now. Its companion was also recently detected by direct imaging, allowing a complete study of this particular object. We aim to fully characterize HD 4747B by combining a well constrained dynamical mass and a study of its observed spectral features in order to test evolutionary models for substellar objects and characterize its atmosphere. We combine the radial velocity measurements of HIRES and CORALIE taken over two decades and high contrast imaging of several epochs from NACO, NIRC2 and SPHERE to obtain a dynamical mass. From the SPHERE data we obtain a low resolution spectrum of the companion from Y to H band, as well as two narrow band-width photometric measurements in the K band. A study of the primary star allows in addition to constrain the age of the system as well as its distance. Thanks to the new SPHERE epoch and NACO archival data combined with previous imaging data and high precision radial velocity measurements, we have been able to derive a well constrained orbit. We derive a dynamical mass of mB=70.0$\pm$1.6 MJup which is higher than a previous study, but in better agreement with the models. By comparing the object with known brown dwarfs spectra, we derive a spectral type of L9 and an effective temperature of 1350$\pm$50 K. With a retrieval analysis we constrain the oxygen and carbon abundances and compare them with the ones from the HR 8799 planets.

A search for new supernova remnants (SNRs) has been conducted using TeV gamma-ray data from the H.E.S.S. Galactic plane survey. As an identification criterion, shell morphologies that are characteristic for known resolved TeV SNRs have been used. Three new SNR candidates were identified in the H.E.S.S. data set with this method. Extensive multiwavelength searches for counterparts were conducted. A radio SNR candidate has been identified to be a counterpart to HESS J1534-571. The TeV source is therefore classified as a SNR. For the other two sources, HESS J1614-518 and HESS J1912+101, no identifying counterparts have been found, thus they remain SNR candidates for the time being. TeV-emitting SNRs are key objects in the context of identifying the accelerators of Galactic cosmic rays. The TeV emission of the relativistic particles in the new sources is examined in view of possible leptonic and hadronic emission scenarios, taking the current multiwavelength knowledge into account.

Aims. We study gamma-ray emission from the shell-type supernova remnant (SNR) RX J0852.0$-$4622 to better characterize its spectral properties and its distribution over the SNR. Methods. The analysis of an extended High Energy Spectroscopic System (H.E.S.S.) data set at very high energies (E > 100 GeV) permits detailed studies, as well as spatially resolved spectroscopy, of the morphology and spectrum of the whole RX J0852.0$-$4622 region. The H.E.S.S. data are combined with archival data from other wavebands and interpreted in the framework of leptonic and hadronic models. The joint Fermi-LAT-H.E.S.S. spectrum allows the direct determination of the spectral characteristics of the parent particle population in leptonic and hadronic scenarios using only GeV-TeV data. Results. An updated analysis of the H.E.S.S. data shows that the spectrum of the entire SNR connects smoothly to the high-energy spectrum measured by Fermi-LAT. The increased data set makes it possible to demonstrate that the H.E.S.S. spectrum deviates significantly from a power law and is well described by both a curved power law and a power law with an exponential cutoff at an energy of $E_\mathrm{cut} = (6.7 \pm 1.2_\mathrm{stat} \pm 1.2_\mathrm{syst})$ TeV. The joint Fermi-LAT-H.E.S.S. spectrum allows the unambiguous identification of the spectral shape as a power law with an exponential cutoff. No significant evidence is found for a variation of the spectral parameters across the SNR, suggesting similar conditions of particle acceleration across the remnant. A simple modeling using one particle population to model the SNR emission demonstrates that both leptonic and hadronic emission scenarios remain plausible. It is also shown that at least a part of the shell emission is likely due to the presence of a pulsar wind nebula around PSR J0855$-$4644.

The next generation ground-based extremely large telescopes (ELTs) present incredible opportunities to discover and characterize diverse planetary systems, even potentially habitable worlds. Adaptive-optics assisted thermal-IR (3-14 micron) imaging is a powerful tool to study exoplanets with extant 6-12 meter telescopes. ELTs have the spatial resolution and sensitivity that offer an unparalleled expansion of the available discovery space. AO-assisted thermal-IR instruments on ELTs will be superior to JWST for high contrast imaging in the thermal-IR, and complementary to high contrast observations at shorter wavelengths, in space or with second-generation extreme AO instruments. With appropriate investments in instrumentation and pre-cursor observations, thermal-IR equipped ELTs could image the first terrestrial and super-earth planets around nearby stars, opening the door to characterization of potentially habitable planets from the ground and space.

Dozens of habitable zone, approximately earth-sized exoplanets are known today. An emerging frontier of exoplanet studies is identifying which of these habitable zone, small planets are actually habitable (have all necessary conditions for life) and, of those, which are earth-like. Many parameters and processes influence habitability, ranging from the orbit through detailed composition including volatiles and organics, to the presence of geological activity and plate tectonics. While some properties will soon be directly observable, others cannot be probed by remote sensing for the foreseeable future. Thus, statistical understanding of planetary systems' formation and evolution is a key supplement to the direct measurements of planet properties. Probabilistically assessing parameters we cannot directly measure is essential to reliably assess habitability, to prioritizing habitable-zone planets for follow-up, and for interpreting possible biosignatures.

Ground-based telescopes have been playing a leading role in exoplanet direct imaging science and technological development for the past two decades and will continue to have an indispensable role for the next decade and beyond. Extreme adaptive optics (AO) systems will advance focal-plane wavefront control and coronagraphy, augmenting the performance of and mitigating risk for WFIRST-CGI, while validating performance requirements and motivating improvements to atmosphere models needed to unambiguously characterize solar system-analogues with HabEx/LUVOIR. Specialized instruments for Extremely Large Telescopes may deliver the first thermal infrared images of rocky planets around Sun-like stars, providing HabEx/LUVOIR with numerous exo-Earth candidates and key ancillary information that can help clarify whether the planets are habitable.

We report the discovery of a bright, brown dwarf companion to the star HIP 64892, imaged with VLT/SPHERE during the SHINE exoplanet survey. The host is a B9.5V member of the Lower-Centaurus-Crux subgroup of the Scorpius Centaurus OB association. The measured angular separation of the companion ($1.2705\pm0.0023$") corresponds to a projected distance of $159\pm12$ AU. We observed the target with the dual-band imaging and long-slit spectroscopy modes of the IRDIS imager to obtain its SED and astrometry. In addition, we reprocessed archival NACO L-band data, from which we also recover the companion. Its SED is consistent with a young (<30 Myr), low surface gravity object with a spectral type of M9$_{\gamma}\pm1$. From comparison with the BT-Settl atmospheric models we estimate an effective temperature of $T_{\textrm{eff}}=2600 \pm 100$ K, and comparison of the companion photometry to the COND evolutionary models yields a mass of $\sim29-37$ M$_{\text{J}}$ at the estimated age of $16^{+15}_{-7}$ Myr for the system. HIP 64892 is a rare example of an extreme-mass ratio system ($q\sim0.01$) and will be useful for testing models relating to the formation and evolution of such low-mass objects.

PKS 0625-354 (z=0.055) was observed with the four H.E.S.S. telescopes in 2012 during 5.5 hours. The source was detected above an energy threshold of 200 GeV at a significance level of 6.1$\sigma$. No significant variability is found in these observations. The source is well described with a power-law spectrum with photon index $\Gamma =2.84 \pm 0.50_{stat} \pm 0.10_{syst}$ and normalization (at $E_0$=1.0 TeV) $N_0(E_0)=(0.58 \pm 0.22_{stat} \pm 0.12_{syst})\times10^{-12}$ TeV$^{-1}$cm$^{-2}$s$^{-1}$. Multi-wavelength data collected with Fermi-LAT, Swift-XRT, Swift-UVOT, ATOM and WISE are also analysed. Significant variability is observed only in the Fermi-LAT $\gamma$-ray and Swift-XRT X-ray energy bands. Having a good multi-wavelength coverage from radio to very high energy, we performed a broadband modelling from two types of emission scenarios. The results from a one zone lepto-hadronic, and a multi-zone leptonic models are compared and discussed. On the grounds of energetics, our analysis favours a leptonic multi-zone model. Models associated to the X-ray variability constraint supports previous results suggesting a BL Lac nature of PKS 0625-354, with, however, a large-scale jet structure typical of a radio galaxy.

Shell-type supernova remnants (SNRs) are considered prime candidates for the acceleration of Galactic cosmic rays (CRs) up to the knee of the CR spectrum at $\mathrm{E} \approx \mathrm{3}\times \mathrm{10}^\mathrm{15}$ eV. Our Milky Way galaxy hosts more than 350 SNRs discovered at radio wavelengths and at high energies, of which 220 fall into the H.E.S.S. Galactic Plane Survey (HGPS) region. Of those, only 50 SNRs are coincident with a H.E.S.S source and in 8 cases the very high-energy (VHE) emission is firmly identified as an SNR. The H.E.S.S. GPS provides us with a legacy for SNR population study in VHE $\gamma$-rays and we use this rich data set to extract VHE flux upper limits from all undetected SNRs. Overall, the derived flux upper limits are not in contradiction with the canonical CR paradigm. Assuming this paradigm holds true, we can constrain typical ambient density values around shell-type SNRs to $n\leq 7~\textrm{cm}^\textrm{-3}$ and electron-to-proton energy fractions above 10~TeV to $\epsilon_\textrm{ep} \leq 5\times 10^{-3}$. Furthermore, comparisons of VHE with radio luminosities in non-interacting SNRs reveal a behaviour that is in agreement with the theory of magnetic field amplification at shell-type SNRs.

The Laser Interferometer Gravitational Wave Observatory (LIGO) consists of two widely separated 4 km laser interferometers designed to detect gravitational waves from distant astrophysical sources in the frequency range from 10 Hz to 10 kHz. The first observation run of the Advanced LIGO detectors started in September 2015 and ended in January 2016. A strain sensitivity of better than $10^{-23}/\sqrt{\text{Hz}}$ was achieved around 100 Hz. Understanding both the fundamental and the technical noise sources was critical for increasing the observable volume in the universe. The average distance at which coalescing binary black hole systems with individual masses of 30 $M_\odot$ could be detected was 1.3 Gpc. Similarly, the range for binary neutron star inspirals was about 75 Mpc. With respect to the initial detectors, the observable volume of Universe increased respectively by a factor 69 and 43. These improvements allowed Advanced LIGO to detect the gravitational wave signal from the binary black hole coalescence, known as GW150914.

The Cherenkov Telescope Array, CTA, will be the major global observatory for very high energy gamma-ray astronomy over the next decade and beyond. The scientific potential of CTA is extremely broad: from understanding the role of relativistic cosmic particles to the search for dark matter. CTA is an explorer of the extreme universe, probing environments from the immediate neighbourhood of black holes to cosmic voids on the largest scales. Covering a huge range in photon energy from 20 GeV to 300 TeV, CTA will improve on all aspects of performance with respect to current instruments. The observatory will operate arrays on sites in both hemispheres to provide full sky coverage and will hence maximize the potential for the rarest phenomena such as very nearby supernovae, gamma-ray bursts or gravitational wave transients. With 99 telescopes on the southern site and 19 telescopes on the northern site, flexible operation will be possible, with sub-arrays available for specific tasks. CTA will have important synergies with many of the new generation of major astronomical and astroparticle observatories. Multi-wavelength and multi-messenger approaches combining CTA data with those from other instruments will lead to a deeper understanding of the broad-band non-thermal properties of target sources. The CTA Observatory will be operated as an open, proposal-driven observatory, with all data available on a public archive after a pre-defined proprietary period. Scientists from institutions worldwide have combined together to form the CTA Consortium. This Consortium has prepared a proposal for a Core Programme of highly motivated observations. The programme, encompassing approximately 40% of the available observing time over the first ten years of CTA operation, is made up of individual Key Science Projects (KSPs), which are presented in this document.