Establishing the role of active galactic nuclei (AGN) during the formation of galaxies remains one of the greatest challenges of galaxy formation theory. Towards addressing this, we summarise our recent work investigating: (1) the physical drivers of ionised outflows and (2) observational signatures of the impact by jets/outflows on star formation and molecular gas content in AGN host galaxies. We confirm a connection between radio emission and extreme ionised gas kinematics in AGN hosts. Emission-line selected AGN are significantly more likely to exhibit ionised outflows (as traced by the [O III] emission line) if the projected linear extent of the radio emission is confined within the spectroscopic aperture. Follow-up high resolution radio observations and integral field spectroscopy of 10 luminous Type 2 AGN reveal moderate power, young (or frustrated) jets interacting with the interstellar medium. We find that these sources live in highly star forming and gas rich galaxies. Additionally, by combining ALMA-derived dust maps with integral field spectroscopy for eight host galaxies of z~2 X-ray AGN, we show that H-alpha emission is an unreliable tracer of star formation. For the five targets with ionised outflows we find no dramatic in-situ shut down of the star formation. Across both of these studies we find that if these AGN do have a negative impact upon their host galaxies, it must be happening on small (unresolved) spatial scales and/or an observable galaxy-wide impact has yet to occur.

As part of our KMOS AGN Survey at High-redshift (KASHz), we present spatially-resolved VLT/KMOS and VLT/SINFONI spectroscopic data and ALMA 870$\mu$m continuum imaging of eight $z$=1.4--2.6 moderate AGN ($L_{\rm 2-10 \rm kev}$ = $10^{42} - 10^{45}$ ergs s$^{-1}$). We map [OIII], H$\alpha$ and rest-frame FIR emission to search for any spatial anti-correlation between ionised outflows (traced by the [OIII] line) and star formation (SF; traced by H$\alpha$ and FIR), that has previously been claimed for some high-z AGN and used as evidence for negative and/or positive AGN feedback. Firstly, we conclude that H$\alpha$ is unreliable to map SF inside our AGN host galaxies based on: (i) SF rates inferred from attenuation-corrected H$\alpha$ can lie below those inferred from FIR; (ii) the FIR continuum is more compact than the H$\alpha$ emission by a factor of $\sim 2$ on average; (iii) in half of our sample, we observe significant spatial offsets between the FIR and H$\alpha$ emission, with an average offset of $1.4\pm0.6$ kpc. Secondly, for the five targets with outflows we find no evidence for a spatial anti-correlation between outflows and SF using either H$\alpha$ or FIR as a tracer. This holds for our re-analysis of a famous $z$=1.6 X-ray AGN (`XID 2028') where positive and negative feedback has been previously claimed. Based on our results, any impact on SF by ionised outflows must be subtle, either occurring on scales below our resolution, or on long timescales.

We present sensitive 850$\mu$m imaging of the COSMOS field using 640hr of new and archival observations taken with SCUBA-2 at the East Asian Observatory's James Clerk Maxwell Telescope. The SCUBA-2 COSMOS survey (S2COSMOS) achieves a median noise level of $\sigma_{850\mu{\mathrm{m}}}$=1.2mJy/beam over an area of 1.6 sq. degree (MAIN; HST/ACS footprint), and $\sigma_{850\mu{\mathrm{m}}}$=1.7mJy/beam over an additional 1 sq. degree of supplementary (SUPP) coverage. We present a catalogue of 1020 and 127 sources detected at a significance level of >4$\sigma$ and >4.3$\sigma$ in the MAIN and SUPP regions, respectively, corresponding to a uniform 2% false-detection rate. We construct the single-dish 850$\mu$m number counts at $S_{850}$>2mJy and show that these S2COSMOS counts are in agreement with previous single-dish surveys, demonstrating that degree-scale fields are sufficient to overcome the effects of cosmic variance in the $S_{850}$=2-10mJy population. To investigate the properties of the galaxies identified by S2COSMOS sources we measure the surface density of near-infrared-selected galaxies around their positions and identify an average excess of 2.0$\pm$0.2 galaxies within a 13$''$ radius (~100kpc at $z$~2). The bulk of these galaxies represent near-infrared-selected SMGs and/or spatially-correlated sources and lie at a median photometric redshift of $z$=2.0$\pm$0.1. Finally, we perform a stacking analysis at sub-millimeter and far-infrared wavelengths of stellar-mass-selected galaxies ($M_{\star}$=10$^{10}$-10$^{12}{\rm M_{\odot}}$) from $z$=0-4, obtaining high-significance detections at 850um in all subsets (SNR=4-30), and investigate the relation between far-infrared luminosity, stellar mass, and the peak wavelength of the dust SED. The publication of this survey adds a new deep, uniform sub-millimeter layer to the wavelength coverage of this well-studied COSMOS field.

Active galactic nuclei (AGN) are believed to regulate star formation inside their host galaxies through "AGN feedback". We summarise our on-going study of luminous AGN (z~0.2-3; L_(AGN,bol)>~10^43 erg/s), which is designed to search for observational signatures of feedback by combining observed star-formation rate (SFR) measurements from statistical samples with cosmological model predictions. Using the EAGLE hydrodynamical cosmological simulations, in combination with our Herschel+ALMA surveys, we show that - even in the presence of AGN feedback - we do not necessarily expect to see any relationships between average galaxy-wide SFRs and instantaneous AGN luminosities. We caution that the correlation with stellar mass for both SFR and AGN luminosity can contribute to apparent observed positive trends between these two quantities. On the other hand, the EAGLE simulations, which reproduce our observations, predict that a signature of AGN feedback can be seen in the wide specific SFR distributions of $all$ massive galaxies (not just AGN hosts). Overall, whilst we can not rule out that AGN have an immediate small-scale impact on in-situ star-formation, all of our results are consistent with a feedback model where galaxy-wide in-situ star formation is not rapidly suppressed by AGN, but where the feedback likely acts over a longer timescale than a single AGN episode.

We identify multi-wavelength counterparts to 1,147 submillimeter sources from the S2COSMOS SCUBA-2 survey of the COSMOS field by employing a recently developed radio$+$machine-learning method trained on a large sample of ALMA-identified submillimeter galaxies (SMGs), including 260 SMGs identified in the AS2COSMOS pilot survey. In total, we identify 1,222 optical/near-infrared(NIR)/radio counterparts to the 897 S2COSMOS submillimeter sources with S$_{850}$>1.6mJy, yielding an overall identification rate of ($78\pm9$)%. We find that ($22\pm5$)% of S2COSMOS sources have multiple identified counterparts. We estimate that roughly 27% of these multiple counterparts within the same SCUBA-2 error circles very likely arise from physically associated galaxies rather than line-of-sight projections by chance. The photometric redshift of our radio$+$machine-learning identified SMGs ranges from z=0.2 to 5.7 and peaks at $z=2.3\pm0.1$. The AGN fraction of our sample is ($19\pm4$)%, which is consistent with that of ALMA SMGs in the literature. Comparing with radio/NIR-detected field galaxy population in the COSMOS field, our radio+machine-learning identified counterparts of SMGs have the highest star-formation rates and stellar masses. These characteristics suggest that our identified counterparts of S2COSMOS sources are a representative sample of SMGs at z<3. We employ our machine-learning technique to the whole COSMOS field and identified 6,877 potential SMGs, most of which are expected to have submillimeter emission fainter than the confusion limit of our S2COSMOS surveys (S$_{850}$<1.5mJy). We study the clustering properties of SMGs based on this statistically large sample, finding that they reside in high-mass dark matter halos ($(1.2\pm0.3)\times10^{13}\,h^{-1}\,\rm M_{\odot}$), which suggests that SMGs may be the progenitors of massive ellipticals we see in the local Universe.

We highlight that the anomalous orbits of Trans-Neptunian Objects (TNOs) and an excess in microlensing events in the 5-year OGLE dataset can be simultaneously explained by a new population of astrophysical bodies with mass several times that of Earth ($M_\oplus$). We take these objects to be primordial black holes (PBHs) and point out the orbits of TNOs would be altered if one of these PBHs was captured by the Solar System, inline with the Planet 9 hypothesis. Capture of a free floating planet is a leading explanation for the origin of Planet 9 and we show that the probability of capturing a PBH instead is comparable. The observational constraints on a PBH in the outer Solar System significantly differ from the case of a new ninth planet. This scenario could be confirmed through annihilation signals from the dark matter microhalo around the PBH.

We investigate the contribution of clumps and satellites to the galaxy mass assembly. We analyzed spatially-resolved Hubble Space Telescope observations (imaging and slitless spectroscopy) of 53 star-forming galaxies at z ~ 1 - 3. We created continuum and emission line maps and pinpointed residual "blobs" detected after subtracting the galaxy disk. Those were separated into compact (unresolved) and extended (resolved) components. Extended components have sizes ~ 2 kpc and comparable stellar mass and age as the galaxy disks, whereas the compact components are 1.5 dex less massive and 0.4 dex younger than the disks. Furthermore the extended blobs are typically found at larger distances from the galaxy barycenter than the compact ones. Prompted by these observations and by the comparison with simulations, we suggest that compact blobs are in-situ formed clumps, whereas the extended ones are accreting satellites. Clumps and satellites enclose respectively ~ 20% and ~< 80% of the galaxy stellar mass, ~ 30% and ~ 20% of its star formation rate. Considering the compact blobs, we statistically estimated that massive clumps (Mstar >~ 10^9 Msun) have lifetimes of ~ 650 Myr, and the less massive ones (10^8 < Mstar < 10^9 Msun) of ~ 145 Myr. This supports simulations predicting long-lived clumps (lifetime > 100 Myr). Finally, ~< 30% (13%) of our sample galaxies are undergoing single (multiple) merger(s), they have a projected separation ~< 10 kpc, and the typical mass ratio of our satellites is 1:5 (but ranges between 1:10 and 1:1), in agreement with literature results for close pair galaxies.

We predict the range of proper motions of 19 satellite galaxies of M31 that would rotationally stabilise the M31 plane of satellites consisting of 15-20 members as identified by Ibata et al. (2013). Our prediction is based purely on the current positions and line-of-sight velocities of these satellites and the assumption that the plane is not a transient feature. These predictions are therefore independent of the current debate about the formation history of this plane. We further comment on the feasibility of measuring these proper motions with future observations by the THEIA satellite mission as well as the currently planned observations by HST and JWST.

We investigate the production of circularly polarised X and gamma-ray signals in cosmic accelerators such as supernova remnants and AGN jets. Proton-proton and proton-photon collisions within these sites produce a charge asymmetry in the distribution of mesons and muons that eventually leads to a net circular polarisation signal as these particles decay radiatively. We find that the fraction of circular polarisation thus produced is at the level of $ 5 \times 10^{-4}$, regardless of the exact beam spectrum, as long as it is made predominantly of protons. While this fraction is very small, the detection of circular polarisation signals in conjunction with high-energy neutrinos would provide an unambiguous signature of the presence of high-energy protons in cosmic accelerators. In Supernovae shocks in particular, this would indicate the presence of relativistic protons hitting stationary protons and/or low-energy photons in the intergalactic or interstellar medium.

We recently discovered that the luminous radio-quiet QSO LBQS0302-0019 at z=3.286 is likely accompanied by an obscured AGN at 20kpc projected distance, which we dubbed Jil. It represents the tightest candidate obscured/unobscured dual AGN system at z>3. To verify the dual AGN scenario we obtained deep $K_s$ band (rest-frame $V$ band) imaging with the VLT/HAWK-I+GRAAL instrument at 0.4" resolution during science verification in January 2018. Indeed, we detect the individual host galaxies of the QSO and Jil with estimated stellar masses of $\log(M_\star/M_{\odot})=11.4\pm0.5$ and $\log(M_\star/M_{\odot})=10.9\pm0.5$, respectively. Near-IR spectra obtained with VLT-KMOS reveal a clear [OIII] $\lambda$5007 line detection at the location of Jil which does not contribute significantly to the $K_s$ band flux. Both observations therefore corroborate the dual AGN scenario. A comparison to Illustris cosmological simulations suggests a parent halo mass of $\log(M_\mathrm{halo}/M_{\odot})=13.2\pm0.5$ for this interacting galaxy system, corresponding to a very massive dark matter halo at that epoch.

We present the star formation rates (SFRs) of a sample of 109 galaxies with X-ray selected active galactic nuclei (AGN) with moderate to high X-ray luminosities (L(2-8keV)= 10^42-10^45 erg/s), at redshifts 1 < z < 4.7, that were selected to be faint or undetected in the Herschel bands. We combine our deep ALMA continuum observations with deblended 8-500{\mu}m photometry from Spitzer and Herschel, and use infrared (IR) SED fitting and AGN - star formation decomposition methods. The addition of the ALMA photometry results in an order of magnitude more X-ray AGN in our sample with a measured SFR (now 37%). The remaining 63% of the sources have SFR upper limits that are typically a factor of 2-10 times lower than the pre-ALMA constraints. With the improved constraints on the IR SEDs, we can now identify a mid-IR (MIR) AGN component in 50% of our sample, compared to only ~1% previously. We further explore the F870{\mu}m/F24{\mu}m-redshift plane as a tool for the identification of MIR emitting AGN, for three different samples representing AGN dominated, star formation dominated, and composite sources. We demonstrate that the F870{\mu}m/F24{\mu}m-redshift plane can successfully split between AGN and star formation dominated sources, and can be used as an AGN identification method.

We show that the Higgsplosion mechanism makes a prediction for the mass and coupling of a WIMP-like minimal scalar dark matter model. In particular the currently favoured minimal value for the Higgsplosion scale, $E_\mathrm{H}\sim 25$ TeV, implies a dark matter mass $m_\mathrm{DM} \sim 1.25$ TeV and a moderate quartic coupling with the Standard Model Higgs field $\lambda_\mathrm{H,DM} \sim 0.4$. This point in the parameter space is still allowed by all current experimental bounds, including direct detection (XENON), indirect detection (HESS, Fermi, Planck) and collider searches. We have updated the scalar dark matter bounds to reflect the latest results from XENON and HESS experiments. We also comment on vacuum stability and dark matter self-interactions in this model.

We present sensitive 870$\mu$m continuum measurements from our ALMA programmes of 114 X-ray selected AGN in the CDF-S and COSMOS fields. We use these observations in combination with data from Spitzer and Herschel to construct a sample of 86 X-ray selected AGN, 63 with ALMA constraints at $z=1.5-3.2$ with stellar mass $>2\times10^{10}M_{\odot}$. We constructed broad-band spectral energy distributions in the infrared band (8 - 1000$\mu$m) and constrain star-formation rates (SFRs) uncontaminated by the AGN. Using a hierarchical Bayesian method that takes into account the information from upper limits, we fit SFR and specific SFR (sSFR) distributions. We explore these distributions as a function of both X-ray luminosity and stellar mass. We compare our measurements to two versions of the EAGLE hydrodynamical simulations: the reference model with AGN feedback and the model without AGN. We find good agreement between the observations and that predicted by the EAGLE reference model for the modes and widths of the sSFR distributions as a function of both X-ray luminosity and stellar mass; however, we found that the EAGLE model without AGN feedback predicts a significantly narrower width when compared to the data. Overall, from the combination of the observations with the model predictions, we conclude that (1) even with AGN feedback, we expect no strong relationship between the sSFR distribution parameters and instantaneous AGN luminosity and (2) a signature of AGN feedback is a broad distribution of sSFRs for all galaxies (not just those hosting an AGN) with stellar masses above $\approx 10^{10}$M$_{\odot}$.

Recently it was shown that dark matter with mass of order the weak scale can be charged under a new long-range force, decoupled from the Standard Model, with only weak constraints from early Universe cosmology. Here we consider the implications of an additional charged particle $C$ that is light enough to lead to significant dissipative dynamics on galactic times scales. We highlight several novel features of this model, which can be relevant even when the $C$ particle constitutes only a small fraction of the number density (and energy density). We assume a small asymmetric abundance of the $C$ particle whose charge is compensated by a heavy $X$ particle so that the relic abundance of dark matter consists mostly of symmetric $X$ and $\bar{X}$, with a small asymmetric component made up of $X$ and $C$. As the universe cools, it undergoes asymmetric recombination binding the free $C$s into $(XC)$ dark atoms efficiently. Even with a tiny asymmetric component, the presence of $C$ particles catalyzes tight coupling between the heavy dark matter $X$ and the dark photon plasma that can lead to a significant suppression of the matter power spectrum on small scales and lead to some of the strongest bounds on such dark matter theories. We find a viable parameter space where structure formation constraints are satisfied and significant dissipative dynamics can occur in galactic haloes but show a large region is excluded. Our model shows that subdominant components in the dark sector can dramatically affect structure formation.

We revisit constraints on dark matter that is charged under a $U(1)$ gauge group in the dark sector, decoupled from Standard Model forces. We find that the strongest constraints in the literature are subject to a number of mitigating factors. For instance, the naive dark matter thermalization timescale in halos is corrected by saturation effects that slow down isotropization for modest ellipticities. The weakened bounds uncover interesting parameter space, making models with weak-scale charged dark matter viable, even with electromagnetic strength interaction. This also leads to the intriguing possibility that dark matter self-interactions within small dwarf galaxies are extremely large, a relatively unexplored regime in current simulations. Such strong interactions suppress heat transfer over scales larger than the dark matter mean free path, inducing a dynamical cutoff length scale above which the system appears to have only feeble interactions. These effects must be taken into account to assess the viability of darkly-charged dark matter. Future analyses and measurements should probe a promising region of parameter space for this model.

In the context of the ESA M5 (medium mission) call we proposed a new satellite mission, Theia, based on relative astrometry and extreme precision to study the motion of very faint objects in the Universe. Theia is primarily designed to study the local dark matter properties, the existence of Earth-like exoplanets in our nearest star systems and the physics of compact objects. Furthermore, about 15 $\%$ of the mission time was dedicated to an open observatory for the wider community to propose complementary science cases. With its unique metrology system and "point and stare" strategy, Theia's precision would have reached the sub micro-arcsecond level. This is about 1000 times better than ESA/Gaia's accuracy for the brightest objects and represents a factor 10-30 improvement for the faintest stars (depending on the exact observational program). In the version submitted to ESA, we proposed an optical (350-1000nm) on-axis TMA telescope. Due to ESA Technology readiness level, the camera's focal plane would have been made of CCD detectors but we anticipated an upgrade with CMOS detectors. Photometric measurements would have been performed during slew time and stabilisation phases needed for reaching the required astrometric precision.

We show that Fermi repulsion can lead to cored density profiles in dwarf galaxies for sub-keV fermionic dark matter. We treat the dark matter as a quasi-degenerate self-gravitating Fermi gas and calculate its density profile assuming hydrostatic equilibrium. We find that suitable dwarf galaxy cores of larger than 130 pc can be achieved for fermion dark matter with mass in the range 70 eV - 400 eV. While in conventional dark matter scenarios, such sub-keV thermal dark matter would be excluded by free streaming bounds, the constraints are ameliorated in models with dark matter at lower temperature than conventional thermal scenarios, such as the Flooded Dark Matter model that we have previously considered. Modifying the arguments of Tremaine and Gunn we derive a conservative lower bound on the mass of fermionic dark matter of 70 eV and a stronger lower bound from Lyman-$\alpha$ clouds of about 470 eV, leading to slightly smaller cores than have been observed. We comment on this result and how the tension is relaxed in dark matter scenarios with non-thermal momentum distributions.

Most dark matter (DM) models set the DM relic density by some interaction with Standard Model particles. Such models generally assume the existence of Standard Model particles early on, with the DM relic density a later consequence of those interactions. Perhaps a more compelling assumption is that DM is not part of the Standard Model sector and a population of DM too is generated at the end of inflation. This democratic assumption does not necessarily provide a natural value for the DM relic density, and superficially leads to too much entropy in the dark sector. We address the latter issue by the late decay of heavy particles produced at early times, associating the DM relic density with the lifetime of a long-lived state. We ask what it would take for this scenario to be compatible with observations in what we call Flooded Dark Matter (FDM) and discuss several interesting consequences. One is that DM can be very light and furthermore, light DM is in some sense the most natural scenario in FDM as it is compatible with larger couplings of the decaying particle. Moreover, the decay of the field with the smallest coupling and hence the longest lifetime dominates the entropy and possibly the matter content of the Universe, a principle we refer to as 'Maximum Baroqueness'. We also show that the dark sector should be colder than the ordinary sector, relaxing the free-streaming constraints on light DM. We will discuss the implications for the core-cusp problem in a follow-up paper. FDM also has interesting baryogenesis implications. One possibility is that both DM and baryon asymmetries are simultaneously diluted by a late entropy dump. Alternatively, FDM is compatible with an elegant non-thermal leptogenesis implementation in which decays of a heavy RH neutrino lead to late time reheating of the Standard Model and provide suitable conditions for creation of a lepton asymmetry.

We show that dissipative dark matter can potentially explain the large observed mass to light ratio of the dwarf satellite galaxies that have been observed in the recently identified planar structure around Andromeda, which are thought to result from tidal forces during a galaxy merger. Whereas dwarf galaxies created from ordinary disks would be dark matter poor, dark matter inside the galactic plane not only provides a source of dark matter, but one that is more readily bound due to the dark matter's lower velocity. This initial N-body study shows that with a thin disk of dark matter inside the baryonic disk, mass-to-light ratios as high as O(30) can be generated when tidal forces pull out patches of sizes similar to the scales of Toomre instabilities of the dark disk. A full simulation will be needed to confirm this result.

We explore the possibility that the dark matter is a condensate of a very light vector boson. Such a condensate could be produced during inflation, provided the vector mass arises via the Steuckelberg mechanism. We derive bounds on the kinetic mixing of the dark matter boson with the photon, and point out several potential signatures of this model.