Measurements of the cosmic far-infrared background (CIB) indicate that emission from many extragalactic phenomena, including star formation and black hole accretion, in the Universe can be obscured by dust. Resolving the CIB to study the population of galaxies in which this activity takes place is a major goal of submillimetre astronomy. Here, we present interferometric 650$\mu$m submillimetre number counts. Using the Band 8 data from the ALMACAL survey, we have analysed 81 ALMA calibrator fields together covering a total area of 5.5~arcmin$^2$. The typical central rms in these fields is $\sim 100 \mu$Jy~beam$^{-1}$ with the deepest maps reaching $\sigma = 47 \mu$Jy~beam$^{-1}$ at sub-arcsec resolution. Multi-wavelength coverage from ALMACAL allows us to exclude contamination from jets associated with the calibrators. However, residual contamination by jets and lensing remain a possibility. Using a signal-to-noise threshold of $4.5\sigma$, we find 21 dusty, star-forming galaxies with 650$\mu$m flux densities of $\geq 0.7 $mJy. At the detection limit we resolve $\simeq 100$ per cent of the CIB at 650$\mu$m, a significant improvement compared to low resolution studies at similar wavelength. We have therefore identified all the sources contributing to the EBL at 650 microns and predict that the contribution from objects with flux 0.7<mJy will be small.

We present results of the MUSE-ALMA Halos, an ongoing study of the Circum-Galactic Medium (CGM) of low redshift galaxies (z < 1.4), currently comprising 14 strong HI absorbers in five quasar fields. We detect 43 galaxies associated with absorbers down to star formation rate (SFR) limits of 0.01-0.1 solar masses/yr, found within impact parameters (b) of 250 kpc from the quasar sightline. Excluding the targeted absorbers, we report a high detection rate of 89 per cent and find that most absorption systems are associated with pairs or groups of galaxies (three to eleven members). We note that galaxies with the smallest impact parameters are not necessarily the closest to the absorbing gas in velocity space. Using a multi-wavelength dataset (UVES/HIRES, HST, MUSE), we combine metal and HI column densities, allowing for derivation of the lower limits of neutral gas metallicity as well as emission line diagnostics (SFR, metallicities) of the ionised gas in the galaxies. We find that groups of associated galaxies follow the canonical relations of N(HI) -- b and W_r(2796) -- b, defining a region in parameter space below which no absorbers are detected. The metallicity of the ISM of associated galaxies, when measured, is higher than the metallicity limits of the absorber. In summary, our findings suggest that the physical properties of the CGM of complex group environments would benefit from associating the kinematics of individual absorbing components with each galaxy member.

Gaia16aye was a binary microlensing event discovered in the direction towards the northern Galactic disc and was one of the first microlensing events detected and alerted to by the Gaia space mission. Its light curve exhibited five distinct brightening episodes, reaching up to I=12 mag, and it was covered in great detail with almost 25,000 data points gathered by a network of telescopes. We present the photometric and spectroscopic follow-up covering 500 days of the event evolution. We employed a full Keplerian binary orbit microlensing model combined with the motion of Earth and Gaia around the Sun to reproduce the complex light curve. The photometric data allowed us to solve the microlensing event entirely and to derive the complete and unique set of orbital parameters of the binary lensing system. We also report on the detection of the first-ever microlensing space-parallax between the Earth and Gaia located at L2. The properties of the binary system were derived from microlensing parameters, and we found that the system is composed of two main-sequence stars with masses 0.57$\pm$0.05 $M_\odot$ and 0.36$\pm$0.03 $M_\odot$ at 780 pc, with an orbital period of 2.88 years and an eccentricity of 0.30. We also predict the astrometric microlensing signal for this binary lens as it will be seen by Gaia as well as the radial velocity curve for the binary system. Events such as Gaia16aye indicate the potential for the microlensing method of probing the mass function of dark objects, including black holes, in directions other than that of the Galactic bulge. This case also emphasises the importance of long-term time-domain coordinated observations that can be made with a network of heterogeneous telescopes.

We are just starting to understand the physical processes driving the dramatic change in cosmic star-formation rate between $z\sim 2$ and the present day. A quantity directly linked to star formation is the molecular gas density, which should be measured through independent methods to explore variations due to cosmic variance and systematic uncertainties. We use intervening CO absorption lines in the spectra of mm-bright background sources to provide a census of the molecular gas mass density of the Universe. The data used in this work are taken from ALMACAL, a wide and deep survey utilizing the ALMA calibrator archive. While we report multiple Galactic absorption lines and one intrinsic absorber, no extragalactic intervening molecular absorbers are detected. However, thanks to the large redshift path surveyed ($\Delta z=182$), we provide constraints on the molecular column density distribution function beyond $z\sim 0$. In addition, we probe column densities of N(H$_2$) > 10$^{16}$ atoms~cm$^{-2}$, five orders of magnitude lower than in previous studies. We use the cosmological hydrodynamical simulation IllustrisTNG to show that our upper limits of $\rho ({\rm H}_2)\lesssim 10^{8.3} \text{M}_{\odot} \text{Mpc}^{-3}$ at $0 < z \leq 1.7$ already provide new constraints on current theoretical predictions of the cold molecular phase of the gas. These results are in agreement with recent CO emission-line surveys and are complementary to those studies. The combined constraints indicate that the present decrease of the cosmic star-formation rate history is consistent with an increasing depletion of molecular gas in galaxies compared to $z\sim 2$.

Galaxy halos appear to be missing a large fraction of their baryons, most probably hiding in the circumgalactic medium (CGM), a diffuse component within the dark matter halo that extends far from the inner regions of the galaxies. A powerful tool to study the CGM gas is offered by absorption lines in the spectra of background quasars. Here, we present optical (MUSE) and mm (ALMA) observations of the field of the quasar Q1130-1449 which includes a log [N(H I)/cm^-2]=21.71+/-0.07 absorber at z=0.313. Ground-based VLT/MUSE 3D spectroscopy shows 11 galaxies at the redshift of the absorber down to a limiting SFR>0.01 M_sun yr^-1 (covering emission lines of [OII], Hbeta, [OIII], [NII] and Halpha), 7 of which are new discoveries. In particular, we report a new emitter with a smaller impact parameter to the quasar line-of-sight (b=10.6 kpc) than the galaxies detected so far. Three of the objects are also detected in CO(1-0) in our ALMA observations indicating long depletion timescales for the molecular gas and kinematics consistent with the ionised gas. We infer from dedicated numerical cosmological RAMSES zoom-in simulations that the physical properties of these objects qualitatively resemble a small group environment, possibly part of a filamentary structure. Based on metallicity and velocity arguments, we conclude that the neutral gas traced in absorption is only partly related to these emitting galaxies while a larger fraction is likely the signature of gas with surface brightness almost four orders of magnitude fainter that current detection limits. Together, these findings challenge a picture where strong-HI quasar absorbers are associated with a single bright galaxy and favour a scenario where the HI gas probed in absorption is related to far more complex galaxy structures.

Accreting supermassive black holes (SMBHs) can exhibit variable emission across the electromagnetic spectrum and over a broad range of timescales. The variability of active galactic nuclei (AGNs) in the ultraviolet and optical is usually at the few tens of per cent level over timescales of hours to weeks. Recently, rare, more dramatic changes to the emission from accreting SMBHs have been observed, including tidal disruption events, 'changing look' AGNs and other extreme variability objects. The physics behind the 're-ignition', enhancement and 'shut-down' of accretion onto SMBHs is not entirely understood. Here we present a rapid increase in ultraviolet-optical emission in the centre of a nearby galaxy, marking the onset of sudden increased accretion onto a SMBH. The optical spectrum of this flare, dubbed AT 2017bgt, exhibits a mix of emission features. Some are typical of luminous, unobscured AGNs, but others are likely driven by Bowen fluorescence - robustly linked here with high-velocity gas in the vicinity of the accreting SMBH. The spectral features and increased ultraviolet flux show little evolution over a period of at least 14 months. This disfavours the tidal disruption of a star as their origin, and instead suggests a longer-term event of intensified accretion. Together with two other recently reported events with similar properties, we define a new class of SMBH-related flares. This has important implications for the classification of different types of enhanced accretion onto SMBHs.

We report on the discovery and follow-up of a peculiar transient, OGLE17aaj, which occurred in the nucleus of a weakly active galaxy. We investigate whether it can be interpreted as a new candidate for a tidal disruption event (TDE). We present the OGLE-IV light curve that covers the slow 60-day-long rise to maximum along with photometric, spectroscopic, and X-ray follow-up during the first year. OGLE17aaj is a nuclear transient exhibiting some properties similar to previously found TDEs, including a long rise time, lack of colour-temperature evolution, and high black-body temperature. On the other hand, its narrow emission lines and slow post-peak evolution are different from previously observed TDEs. Its spectrum and light-curve evolution is similar to F01004-2237 and AT 2017bgt. Signatures of historical low-level nuclear variability suggest that OGLE17aaj may instead be related to a new type of accretion event in active super-massive black holes.

SN 2017dio shows both spectral characteristics of a type-Ic supernova (SN) and signs of a hydrogen-rich circumstellar medium (CSM). Prominent, narrow emission lines of H and He are superposed on the continuum. Subsequent evolution revealed that the SN ejecta are interacting with the CSM. The initial SN Ic identification was confirmed by removing the CSM interaction component from the spectrum and comparing with known SNe Ic, and reversely, adding a CSM interaction component to the spectra of known SNe Ic and comparing them to SN 2017dio. Excellent agreement was obtained with both procedures, reinforcing the SN Ic classification. The light curve constrains the pre-interaction SN Ic peak absolute magnitude to be around $M_g = -17.6$ mag. No evidence of significant extinction is found, ruling out a brighter luminosity required by a SN Ia classification. These pieces of evidence support the view that SN 2017dio is a SN Ic, and therefore the first firm case of a SN Ic with signatures of hydrogen-rich CSM in the early spectrum. The CSM is unlikely to have been shaped by steady-state stellar winds. The mass loss of the progenitor star must have been intense, $\dot{M} \sim 0.02$ $(\epsilon_{H\alpha}/0.01)^{-1}$ $(v_\textrm{wind}/500$ km s$^{-1}$) $(v_\textrm{shock}/10 000$ km s$^{-1})^{-3}$ $M_\odot$~yr$^{-1}$, peaking at a few decades before the SN. Such a high mass loss rate might have been experienced by the progenitor through eruptions or binary stripping.

Gravitational waves were discovered with the detection of binary black hole mergers and they should also be detectable from lower mass neutron star mergers. These are predicted to eject material rich in heavy radioactive isotopes that can power an electromagnetic signal called a kilonova. The gravitational wave source GW170817 arose from a binary neutron star merger in the nearby Universe with a relatively well confined sky position and distance estimate. Here we report observations and physical modelling of a rapidly fading electromagnetic transient in the galaxy NGC4993, which is spatially coincident with GW170817 and a weak short gamma-ray burst. The transient has physical parameters broadly matching the theoretical predictions of blue kilonovae from neutron star mergers. The emitted electromagnetic radiation can be explained with an ejected mass of 0.04 +/- 0.01 Msol, with an opacity of kappa <= 0.5 cm2/gm at a velocity of 0.2 +/- 0.1c. The power source is constrained to have a power law slope of beta = -1.2 +/- 0.3, consistent with radioactive powering from r-process nuclides. We identify line features in the spectra that are consistent with light r-process elements (90 < A < 140). As it fades, the transient rapidly becomes red, and emission may have contribution by a higher opacity, lanthanide-rich ejecta component. This indicates that neutron star mergers produce gravitational waves, radioactively powered kilonovae, and are a nucleosynthetic source of the r-process elements.

We present a collection of 450 598 eclipsing and ellipsoidal binary systems detected in the OGLE fields toward the Galactic bulge. The collection consists of binary systems of all types: detached, semi-detached, and contact eclipsing binaries, RS CVn stars, cataclysmic variables, HW Vir binaries, double periodic variables, and even planetary transits. For all stars we provide the I- and V-band time-series photometry obtained during the OGLE-II, OGLE-III, and OGLE-IV surveys. We discuss methods used to identify binary systems in the OGLE data and present several objects of particular interest.

We present the collection of eclipsing binaries in the Large and Small Magellanic Clouds, based on the OGLE survey. It contains 48605 systems, 40204 belonging to the LMC and 8401 to the SMC. Out of the total number of the presented binaries, 16374 are the new ones. We present the time-series photometry obtained for the selected objects during the fourth phase of the OGLE project. The catalog has been created using a two step machine learning procedure based on the Random Forest algorithm.

We present the discovery and first three months of follow-up observations of a currently on-going unusual transient detected by the OGLE-IV survey, located in the centre of a galaxy at redshift z=0.1655. The long rise to absolute magnitude of -20.5 mag, slow decline, very broad He and H spectral features make OGLE16aaa similar to other optical/UV Tidal Disruption Events (TDEs). Weak narrow emission lines in the spectrum and archival photometric observations suggest the host galaxy is a weak-line Active Galactic Nucleus (AGN), which has been accreting at higher rate in the past. OGLE16aaa, along with SDSS J0748, seems to form a sub-class of TDEs by weakly or recently active super-massive black holes (SMBHs). This class might bridge the TDEs by quiescent SMBHs and flares observed as "changing-look QSOs", if we interpret the latter as TDEs. If this picture is true, the previously applied requirement for identifying a flare as a TDE that it had to come from an inactive nucleus, could be leading to observational bias in TDE selection, thus affecting TDE-rate estimations.

We use the fundamental-mode RR Lyr-type variable stars (RRab) from OGLE-IV survey to draw a 3D picture of the central part of the tidally disrupted Sagittarius Dwarf Spheroidal (Sgr dSph) galaxy. We estimate the line-of-sight thickness of the Sgr dSph stream to be FWHM_cen=2.42 kpc. Based on OGLE-IV observations collected in seasons 2011-2014 we conduct a comprehensive study of stellar variability in the field of the globular cluster M54 (NGC 6715) residing in the core of this dwarf galaxy. Among the total number of 268 detected variable stars we report the identification of 174 RR Lyr stars, four Type II Cepheids, 51 semi-regular variable red giants, three SX Phe-type stars, 18 eclipsing binary systems. Eighty-three variable stars are new discoveries. The distance to the cluster determined from RRab stars is d_M54 = 26.7 +/-0.03(stat) +/-1.3(sys) kpc. From the location of RRab stars in the period-amplitude (Bailey) diagram we confirm the presence of two old populations, both in the cluster and the Sgr dSph stream.

We report the discovery and characterisation of a deeply eclipsing AM CVn-system, Gaia14aae (= ASSASN-14cn). Gaia14aae was identified independently by the All-Sky Automated Survey for Supernovae (ASAS-SN; Shappee et al. 2014) and by the Gaia Science Alerts project, during two separate outbursts. A third outburst is seen in archival Pan-STARRS-1 (PS1; Schlafly et al. 2012; Tonry et al. 2012; Magnier et al. 2013) and ASAS-SN data. Spectroscopy reveals a hot, hydrogen-deficient spectrum with clear double-peaked emission lines, consistent with an accreting double degenerate classification. We use follow-up photometry to constrain the orbital parameters of the system. We find an orbital period of 49.71 min, which places Gaia14aae at the long period extremum of the outbursting AM CVn period distribution. Gaia14aae is dominated by the light from its accreting white dwarf. Assuming an orbital inclination of 90 degrees for the binary system, the contact phases of the white dwarf lead to lower limits of 0.78 M solar and 0.015 M solar on the masses of the accretor and donor respectively and a lower limit on the mass ratio of 0.019. Gaia14aae is only the third eclipsing AM CVn star known, and the first in which the WD is totally eclipsed. Using a helium WD model, we estimate the accretor's effective temperature to be 12900+-200 K. The three out-burst events occurred within 4 months of each other, while no other outburst activity is seen in the previous 8 years of Catalina Real-time Transient Survey (CRTS; Drake et al. 2009), Pan-STARRS-1 and ASAS-SN data. This suggests that these events might be rebrightenings of the first outburst rather than individual events.