Submillimetre galaxies (SMGs) have long posed a challenge for theorists, and self-consistently reproducing the properties of the SMG population in a large-volume cosmological hydrodynamical simulation has not yet been achieved. In this work, we employ a relatively simple method based on previous work to predict the submm flux densities of simulated SMGs drawn from cosmological simulations from the Illustris and IllustrisTNG projects and compare the predicted number counts with observations. We find that the predicted SMG number counts based on IllustrisTNG are significantly less than observed (more than 1 dex at $S_{850} \gtrsim 4$ mJy). The simulation from the original Illustris project yields more SMGs than IllustrisTNG: the predicted counts are consistent with those observed at both $S_{850} \lesssim 5$ mJy and $S_{850} \gtrsim 9$ mJy and only a factor of $\sim 2$ lower than observed at intermediate flux densities. We demonstrate that IllustrisTNG hosts fewer SMGs than Illustris because in the former, high-mass ($M_{\star} \sim 10^{11} \, {\rm M}_{\odot}$) $z \sim 2-3$ galaxies have lower dust masses and star formation rates (SFRs) than in Illustris owing to differences in the sub-grid models for stellar or/and active galactic nucleus (AGN) feedback between the two simulations (we unfortunately cannot isolate the specific cause(s) post hoc). Our results demonstrate that because our method enables predicting SMG number counts in post-processing with a negligible computational expense, SMGs can provide useful constraints for tuning sub-grid models in future large-volume cosmological simulations. Higher resolution, which would lead to stronger SFR enhancements in starbursts, could at least partially reconcile the discrepancy between the IllustrisTNG SMG number counts and those observed.

In recent observational studies, star-forming galaxies have been shown to follow a relation often dubbed the fundamental metallicity relation (FMR). This relation links the stellar mass of a galaxy with its star formation rate (SFR) and its gas-phase metallicity. Specifically, the FMR predicts that galaxies, at a given stellar mass, exhibit lower metallicities for higher SFRs. This trend is qualitatively consistent with observations of galaxy pairs, which have been robustly shown to experience increasing gas-phase metallicity dilution and enhanced star formation activity with decreasing projected separation. In this work, we show that, despite the qualitative consistency with FMR expectations, the observed O/H dilution in galaxy pairs of the Sloan Digital Sky Survey is stronger than what is predicted by the FMR. We conclude that the evolutionary phase of galaxies interacting with companions is not encoded in the FMR, and thus, mergers constitute a clearly defined population of outliers. We find that galaxies in pairs are consistent with the FMR only when their separation is larger than 110 kpc. Finally, we also quantify the local environment of the pairs using the number of galaxy neighbours within $2\, \mathrm{Mpc}$, $N_2$, and the projected separation to the second closest galaxy, $r_2$. We find that pairs are more sensitive to a second companion than to the local galaxy density, displaying less elevated SFRs with smaller values of $r_2$.

We report two secure ($z=3.775, 4.012$) and one tentative ($z\approx3.767$) spectroscopic confirmations of massive and quiescent galaxies through $K$-band observations with Keck/MOSFIRE and VLT/X-Shooter. The stellar continuum emission, the absence of strong nebular emission lines and the lack of significant far-infrared detections confirm the passive nature of these objects, disfavoring the alternative solution of low-redshift dusty star-forming interlopers. We derive stellar masses of $\mathrm{log}(M_{\star}/M_\odot)\sim11$ and ongoing star formation rates placing these galaxies $\gtrsim 1-2$ dex below the main sequence at their redshifts. The adopted parametrization of the star formation history suggests that these sources experienced a strong ($\langle \rm SFR \rangle \sim 1200-3500\,M_\odot\,\mathrm{yr}^{-1}$) and short ($\sim 50$ Myr) burst of star formation, peaking $\sim 150-500$ Myr before the time of observation, all properties reminiscent of the characteristics of sub-millimeter galaxies (SMGs) at $z>4$. We investigate this connection by comparing the comoving number densities and the properties of these two populations. We find a fair agreement only with the deepest sub-mm surveys detecting not only the most extreme starbursts, but also more normal galaxies. We support these findings by further exploring the Illustris-TNG cosmological simulation, retrieving populations of both fully quenched massive galaxies at $z\sim3-4$ and SMGs at $z\sim4-5$, with number densities and properties in agreement with the observations at $z\sim3$, but in increasing tension at higher redshift. Nevertheless, as suggested by the observations, not all the progenitors of quiescent galaxies at these redshifts shine as bright SMGs in their past and, similarly, not all bright SMGs quench by $z\sim3$, both fractions depending on the threshold assumed to define the SMGs themselves.

Galaxies are surrounded by massive gas reservoirs (i.e. the circumgalactic medium; CGM) which play a key role in their evolution. The properties of the CGM, which are dependent on a variety of internal and environmental factors, are often inferred from absorption line surveys which rely on a limited number of single lines-of-sight. In this work we present an analysis of 28 galaxy haloes selected from the Auriga project, a cosmological magneto-hydrodynamical zoom-in simulation suite of isolated Milky Way-mass galaxies, to understand the impact of CGM diversity on observational studies. Although the Auriga haloes are selected to populate a narrow range in halo mass, our work demonstrates that the CGM of L* galaxies is extremely diverse: column densities of commonly observed species span ~3-4 dex and their covering fractions range from ~5 to 90 per cent. Despite this diversity, we identify the following correlations: 1) the covering fractions (CF) of hydrogen and metals of the Auriga haloes positively correlate with stellar mass, 2) the CF of H I, C IV, and Si II anticorrelate with active galactic nucleus luminosity due to ionization effects, and 3) the CF of H I, C IV, and Si II positively correlate with galaxy disc fraction due to outflows populating the CGM with cool and dense gas. The Auriga sample demonstrates striking diversity within the CGM of L* galaxies, which poses a challenge for observations reconstructing CGM characteristics from limited samples, and also indicates that long-term merger assembly history and recent star formation are not the dominant sculptors of the CGM.

We use high-resolution hydrodynamical simulation to test the difference of halo properties in cold dark matter (CDM) and a self-interacting dark matter (SIDM) scenario with a constant cross-section of $\sigma^\text{T}/m_{\chi}=1\;\text{cm}^{2}\text{g}^{-1}$. We find that the interplay between dark matter self-interaction and baryonic physics induces a complex evolution of the halo properties, which depends on the halo mass and morphological type, as well as on the halo mass accretion history. While high mass haloes, selected as analogues of early-type galaxies, show cored profiles in the SIDM run, systems of intermediate mass and with a significant disk component can develop a profile that is similar or cuspier than in CDM. The final properties of SIDM haloes - measured at z=0.2 - correlate with the halo concentration and formation time, suggesting that the differences between different systems are due to the fact that we are observing the impact self-interaction. We also search for signatures of self-interacting dark matter in the lensing signal of the main haloes and hints of potential differences in the distribution of Einstein radii, which suggests that future wide-field survey might be able to distinguish between CDM and SIDM models on this basis. Finally, we find that the subhalo abundances are not altered in the adopted SIDM model with respect to CDM.

Galactic winds exhibit a multiphase structure that consists of hot-diffuse and cold-dense phases. Here we present high-resolution idealised simulations of the interaction of a hot supersonic wind with a cold cloud with the moving-mesh code arepo in setups with and without radiative cooling. We demonstrate that cooling causes clouds with sizes larger than the cooling length to fragment in two- and three-dimensional simulations (2D and 3D). We confirm earlier 2D simulations by McCourt et al. 2018 and highlight differences of the shattering processes of 3D clouds that are exposed to a hot wind. The fragmentation process is quantified with a friends-of-friends analysis of shattered cloudlets and density power spectra. Those show that radiative cooling causes the power spectral index to gradually increase when the initial cloud radius is larger than the cooling length and with increasing time until the cloud is fully dissolved in the hot wind. A resolution of around 1 pc is required to reveal the effect of cooling-induced fragmentation of a 100 pc outflowing cloud. Thus, state-of-the-art cosmological zoom simulations of the circumgalactic medium (CGM) fall short by orders of magnitudes from resolving this fragmentation process. This physics is, however, necessary to reliably model observed column densities and covering fractions of Lyman-$\alpha$ haloes, high-velocity clouds, and broad-line regions of active galactic nuclei.

Observational studies have revealed that galaxy pairs tend to have lower gas-phase metallicity than isolated galaxies. This metallicity deficiency can be caused by inflows of low-metallicity gas due to the tidal forces and gravitational torques associated with galaxy mergers, diluting the metal content of the central region. In this work we demonstrate that such metallicity dilution occurs in state-of-the-art cosmological simulations of galaxy formation. We find that the dilution is typically 0.1 dex for major mergers, and is noticeable at projected separations smaller than $40$ kpc. For minor mergers the metallicity dilution is still present, even though the amplitude is significantly smaller. Consistent with previous analysis of observed galaxies we find that mergers are outliers from the \emph{fundamental metallicity relation}, with deviations being larger than expected for a Gaussian distribution of residuals. Our large sample of mergers within full cosmological simulations also makes it possible to estimate how the star formation rate enhancement and gas consumption timescale behave as a function of the merger mass ratio. We confirm that strong starbursts are likely to occur in major mergers, but they can also arise in minor mergers if more than two galaxies are participating in the interaction, a scenario that has largely been ignored in previous work based on idealised isolated merger simulations.

We present here a survey of high-ionization absorption lines in the afterglow spectra of long-duration gamma-ray bursts (GRBs) obtained with the VLT/X-shooter spectrograph. Our main goal is to investigate the circumburst medium in the natal regions of GRBs. Our primary focus is on the NV 1238,1242 line transitions, but we also discuss other high-ionization lines such as OVI, CIV and SiIV. We find no correlation between the column density of NV and the neutral gas properties such as metallicity, HI column density and dust depletion, however the relative velocity of NV, typically a blueshift with respect to the neutral gas, is found to be correlated with the column density of HI. This may be explained if the NV gas is part of an HII region hosting the GRB, where the region's expansion is confined by dense, neutral gas in the GRB's host galaxy. We find tentative evidence (at 2-sigma significance) that the X-ray derived column density, N_H,X, may be correlated with the column density of NV, which would indicate that both measurements are sensitive to the column density of the gas located in the vicinity of the GRB. We investigate the scenario where NV (and also OVI) is produced by recombination after the corresponding atoms have been stripped entirely of their electrons by the initial prompt emission, in contrast to previous models where highly-ionized gas is produced by photoionization from the GRB afterglow.

The part played by stars in the ionization of the intergalactic medium remains an open question. A key issue is the proportion of the stellar ionizing radiation that escapes the galaxies in which it is produced. Spectroscopy of gamma-ray burst afterglows can be used to determine the neutral hydrogen column-density in their host galaxies and hence the opacity to extreme ultra-violet radiation along the lines-of-sight to the bursts. Thus, making the reasonable assumption that long-duration GRB locations are representative of the sites of massive stars that dominate EUV production, one can calculate an average escape fraction of ionizing radiation in a way that is independent of galaxy size, luminosity or underlying spectrum. Here we present a sample of NH measures for 138 GRBs in the range 1.6<z<6.7 and use it to establish an average escape fraction at the Lyman limit of <fesc>~0.005, with a 98% confidence upper limit of ~0.015. This analysis suggests that stars provide a small contribution to the ionizing radiation budget of the IGM at z<5, where the bulk of the bursts lie. At higher redshifts, z>5, firm conclusions are limited by the small size of the GRB sample, but any decline in average HI column-density seems to be modest. We also find no indication of a significant correlation of NH with galaxy UV luminosity or host stellar mass, for the subset of events for which these are available. We discuss in some detail a number of selection effects and potential biases. Drawing on a range of evidence we argue that such effects, while not negligible, are unlikely to produce systematic errors of more than a factor ~2, and so would not affect the primary conclusions. Given that many GRB hosts are low metallicity, high specific star-formation rate, dwarf galaxies, these results present a particular problem for the hypothesis that such galaxies dominated the reionization of the universe.

In this work we present spectra of all $\gamma$-ray burst (GRB) afterglows that have been promptly observed with the X-shooter spectrograph until 31-03-2017. In total, we obtained spectroscopic observations of 103 individual GRBs observed within 48 hours of the GRB trigger. Redshifts have been measured for 97 per cent of these, covering a redshift range from 0.059 to 7.84. Based on a set of observational selection criteria that minimize biases with regards to intrinsic properties of the GRBs, the follow-up effort has been focused on producing a homogeneous sample of 93 afterglow spectra for GRBs discovered by the Swift satellite. We here provide a public release of all the reduced spectra, including continuum estimates and telluric absorption corrections. For completeness, we also provide reductions for the 18 late-time observations of the underlying host galaxies. We provide an assessment of the degree of completeness with respect to the parent GRB population, in terms of the X-ray properties of the bursts in the sample and find that the sample presented here is representative of the full Swift sample. We constrain the fraction of dark bursts to be < 28 per cent and we confirm previous results that higher optical darkness is correlated with increased X-ray absorption. For the 42 bursts for which it is possible, we provide a measurement of the neutral hydrogen column density, increasing the total number of published HI column density measurements by $\sim$ 33 per cent. This dataset provides a unique resource to study the ISM across cosmic time, from the local progenitor surroundings to the intervening universe.

Galaxies are surrounded by sizeable gas reservoirs which host a significant amount of metals: the circum-galactic medium (CGM). The CGM acts as a mediator between the galaxy and the extra-galactic medium. However, our understanding of how galaxy mergers, a major evolutionary transformation, impact the CGM remains deficient. We present a theoretical study of the effect of galaxy mergers on the CGM. We use hydrodynamical cosmological zoom-in simulations of a major merger selected from the Illustris project such that the z=0 descendant has a halo mass and stellar mass comparable to the Milky Way. To study the CGM we then re-simulated this system at a 40 times better mass resolution, and included detailed post-processing ionization modelling. Our work demonstrates the effect the merger has on the characteristic size of the CGM, its metallicity, and the predicted covering fraction of various commonly observed gas-phase species, such as H I, C IV, and O VI. We show that merger-induced outflows can increase the CGM metallicity by 0.2-0.3 dex within 0.5 Gyr post-merger. These effects last up to 6 Gyr post-merger. While the merger increases the total metal covering fractions by factors of 2-3, the covering fractions of commonly observed UV ions decrease due to the hard ionizing radiation from the active galactic nucleus, which we model explicitly. Our study of the single simulated major merger presented in this work demonstrates the significant impact that a galaxy interaction can have on the size, metallicity, and observed column densities of the CGM.

We trace the specific star formation rate (sSFR) of massive star-forming galaxies ($\gtrsim\!10^{10}\,\mathcal{M}_\odot$) from $z\sim2$ to 7. Our method is substantially different from previous analyses, as it does not rely on direct estimates of star formation rate, but on the differential evolution of the galaxy stellar mass function (SMF). We show the reliability of this approach by means of semi-analytical and hydrodynamical cosmological simulations. We then apply it to real data, using the SMFs derived in the COSMOS and CANDELS fields. We find that the sSFR is proportional to $(1+z)^{1.1\pm0.2}$ at $z>2$, in agreement with other observations but in tension with the steeper evolution predicted by simulations from $z\sim4$ to 2. We investigate the impact of several sources of observational bias, which however cannot account for this discrepancy. Although the SMF of high-redshift galaxies is still affected by significant errors, we show that future large-area surveys will substantially reduce them, making our method an effective tool to probe the massive end of the main sequence of star-forming galaxies.

The stellar disk of the Milky Way shows complex spatial and abundance structure that is central to understanding the key physical mechanisms responsible for shaping our Galaxy. In this study, we use six very high resolution cosmological zoom simulations of Milky Way-sized haloes to study the prevalence and formation of chemically distinct disc components. We find that our simulations develop a clearly bimodal distribution in the $[\rm \alpha/Fe]$ -- $[\rm Fe/H]$ plane. We find two main pathways to creating this dichotomy which operate in different regions of the galaxies: a) an early ($z>1$) and intense high-$\rm[\alpha/Fe]$ star formation phase in the inner region ($R\lesssim 5$ kpc) induced by gas-rich mergers, followed by more quiescent low-$\rm[\alpha/Fe]$ star formation; and b) an early phase of high-$\rm[\alpha/Fe]$ star formation in the outer disc followed by a shrinking of the gas disc owing to a temporarily lowered gas accretion rate, after which disc growth resumes. In process b), a double-peaked star formation history around the time and radius of disc shrinking accentuates the dichotomy. If the early star formation phase is prolonged (rather than short and intense), chemical evolution proceeds as per process a) in the inner region, but the dichotomy is less clear. In the outer region, the dichotomy is only evident if the first intense phase of star formation covers a large enough radial range before disc shrinking occurs; otherwise, the outer disc consists of only low-$\rm[\alpha/Fe]$ sequence stars. We discuss the implication that both processes occurred in the Milky Way.

The SUperluminous Supernova Host galaxIES (SUSHIES) survey aims to provide strong new constraints on the progenitors of superluminous supernovae (SLSNe) by understanding the relationship to their host galaxies. We present the photometric properties of 53 H-poor and 16 H-rich SLSN host galaxies out to $z\sim4$. We model their spectral energy distributions to derive physical properties, which we compare with other galaxy populations. At low redshift, H-poor SLSNe are preferentially found in very blue, low-mass galaxies with high average specific star-formation rates. As redshift increases, the host population follows the general evolution of star-forming galaxies towards more luminous galaxies. After accounting for secular evolution, we find evidence for differential evolution in galaxy mass, but not in the $B$-band and the far UV luminosity ($3\sigma$ confidence). Most remarkable is the scarcity of hosts with stellar masses above $10^{10}~M_\odot$ for both classes of SLSNe. In the case of H-poor SLSNe, we attribute this to a stifled production efficiency above $\sim0.4$ solar metallicity. However, we argue that, in addition to low metallicity, a short-lived stellar population is also required to regulate the SLSN production. H-rich SLSNe are found in a very diverse population of star-forming galaxies. Still, the scarcity of massive hosts suggests a stifled production efficiency above $\sim0.8$ solar metallicity. The large dispersion of the H-rich SLSNe host properties is in stark contrast to those of gamma-ray burst, regular core-collapse SN, and H-poor SLSNe host galaxies. We propose that multiple progenitor channels give rise to this sub-class.

We present a comprehensive study of the relations between gas kinematics, metallicity, and stellar mass in a sample of 82 GRB-selected galaxies using absorption and emission methods. We find the velocity widths of both emission and absorption profiles to be a proxy of stellar mass. We also investigate the velocity-metallicity correlation and its evolution with redshift and find the correlation derived from emission lines to have a significantly smaller scatter compared to that found using absorption lines. Using 33 GRB hosts with measured stellar mass and metallicitiy, we study the mass-metallicity relation for GRB host galaxies in a stellar mass range of $10^{8.2} M_{\odot}$ to $10^{11.1} M_{\odot}$ and a redshift range of $ z\sim 0.3-3.4$. The GRB-selected galaxies appear to track the mass-metallicity relation of star forming galaxies but with an offset of 0.15 towards lower metallicities. This offset is comparable with the average error-bar on the metallicity measurements of the GRB sample and also the scatter on the MZ relation of the general population. It is hard to decide whether this relatively small offset is due to systematic effects or the intrinsic nature of GRB hosts. We also investigate the possibility of using absorption-line metallicity measurements of GRB hosts to study the mass-metallicity relation at high redshifts. Our analysis shows that the metallicity measurements from absorption methods can significantly differ from emission metallicities and assuming identical measurements from the two methods may result in erroneous conclusions.

It is notoriously difficult to localize short $\gamma$-ray bursts (sGRBs) and their hosts to measure their redshifts. These measurements, however, are critical to constrain the nature of sGRB progenitors, their redshift distribution and the $r$-process element enrichment history of the universe. Here, we present spectroscopy of the host galaxy of GRB111117A and measure its redshift to be $z = 2.211$. This makes GRB111117A the most distant high-confidence short duration GRB detected to date. Our spectroscopic redshift supersedes a lower, previously estimated photometric redshift value for this burst. We use the spectroscopic redshift, as well as new imaging data to constrain the nature of the host galaxy and the physical parameters of the GRB. The rest-frame X-ray derived hydrogen column density, for example, is the highest compared to a complete sample of sGRBs and seems to follow the evolution with redshift as traced by the hosts of long GRBs (lGRBs). The host lies in the brighter end of the expected sGRB host brightness distribution at $z = 2.211$, and is actively forming stars. Using the host as a benchmark for redshift determination, we find that between 43 and 71 per cent of all sGRB redshifts should be missed due to host faintness for hosts at $z\sim2$. The high redshift of GRB111117A is evidence against a lognormal delay-time model for sGRBs through the predicted redshift distribution of sGRBs, which is very sensitive to high-$z$ sGRBs. From the age of the universe at the time of GRB explosion, an initial neutron star (NS) separation of $a_0 < 3.2~R_\odot$ is required in the case where the progenitor system is a circular pair of inspiralling NSs. This constraint excludes some of the longest sGRB formation channels for this burst.

Galaxy mergers are believed to play a key role in transforming star-forming disk galaxies into quenched ellipticals. Most of our theoretical knowledge about such morphological transformations does, however, rely on idealised simulations where processes such as cooling of hot halo gas into the disk and gas accretion in the post-merger phase are not treated in a self-consistent cosmological fashion. In this paper we study the morphological evolution of the stellar components of four major mergers occurring at z=0.5 in cosmological hydrodynamical zoom-simulations. In all simulations the merger reduces the disk mass-fraction, but all galaxies simulated at our highest resolution regrow a significant disk by z=0 (with a disk fraction larger than 24%). For runs with our default physics model, which includes galactic winds from star formation and black hole feedback, none of the merger remnants are quenched, but in a set of simulations with stronger black hole feedback we find that major mergers can indeed quench galaxies. We conclude that major merger remnants commonly evolve into star-forming disk galaxies, unless sufficiently strong AGN feedback assists in the quenching of the remnant.

Gladders et al. have recently suggested that the star formation histories (SFHs) of individual galaxies are characterized by a log-normal function in time, implying a slow decline rather than rapid quenching. We test their conjecture on theoretical SFHs from the cosmological simulation Illustris and on observationally inferred SFHs. While the log-normal form necessarily ignores short-lived features such as starbursts, it fits the overall shape of the majority of SFHs very well. In particular, 85% of the cumulative SFHs are fitted to within a maximum error of 5% of the total stellar mass formed, and 99% to within 10%. The log-normal performs systematically better than the commonly used delayed-tau model, and is superseded only by functions with more than three free parameters. Poor fits are mostly found in galaxies that were rapidly quenched after becoming satellites. We explore the log-normal parameter space of normalization, peak time, and full width at half maximum, and find that the simulated and observed samples occupy similar regions, though Illustris predicts wider, later-forming SFHs on average. The ensemble of log-normal fits correctly reproduces complex metrics such as the evolution of Illustris galaxies across the star formation main sequence, but overpredicts their quenching timescales. SFHs in Illustris are a diverse population not determined by any one physical property of galaxies, but follow a tight relation, where $\mathrm{width}\propto\mathrm{(peak\ time)}^{3/2}$. We show that such a relation can be explained qualitatively (though not quantitatively) by a close connection between the growth of dark matter halos and their galaxies.

Galaxy formation models are now able to reproduce observed relations such as the relation between galaxies' star formation rates (SFRs) and stellar masses ($M_*$) and the stellar mass--halo mass relation. We demonstrate that comparisons of the short-timescale variability in galaxy SFRs with observational data provide an additional useful constraint on the physics of galaxy formation feedback. We apply SFR indicators with different sensitivity timescales to galaxies from the Feedback in Realistic Environments (FIRE) simulations. We find that the SFR--$M_*$ relation has a significantly greater scatter when the H$\alpha$-derived SFR is considered compared with when the far-ultraviolet (FUV)-based SFR is used. This difference is a direct consequence of bursty star formation because the FIRE galaxies exhibit order-of-magnitude SFR variations over timescales of a few Myr. We show that the difference in the scatter between the simulated H$\alpha$- and FUV-derived SFR--$M_*$ relations at $z=2$ is consistent with observational constraints. We also find that the H$\alpha$/FUV ratios predicted by the simulations at $z=0$ are similar to those observed for local galaxies except for a population of low-mass ($M_* \lesssim 10^{9.5} {\rm M}_\odot$) simulated galaxies with lower H$\alpha$/FUV ratios than observed. We suggest that future cosmological simulations should compare the H$\alpha$/FUV ratios of their galaxies with observations to constrain the feedback models employed.

We present the extensive follow-up campaign on the afterglow of GRB 110715A at 17 different wavelengths, from X-ray to radio bands, starting 81 seconds after the burst and extending up to 74 days later. We performed for the first time a GRB afterglow observation with the ALMA observatory. We find that the afterglow of GRB 110715A is very bright at optical and radio wavelengths. We use optical and near infrared spectroscopy to provide further information about the progenitor's environment and its host galaxy. The spectrum shows weak absorption features at a redshift z = 0.8225, which reveal a host galaxy environment with low ionization, column density and dynamical activity. Late deep imaging shows a very faint galaxy, consistent with the spectroscopic results. The broadband afterglow emission is modelled with synchrotron radiation using a numerical algorithm and we determine the best fit parameters using Bayesian inference in order to constrain the physical parameters of the jet and the medium in which the relativistic shock propagates. We fitted our data with a variety of models, including different density profiles and energy injections. Although the general behaviour can be roughly described by these models, none of them are able to fully explain all data points simultaneously. GRB 110715A shows the complexity of reproducing extensive multi-wavelength broadband afterglow observations, and the need of good sampling in wavelength and time and more complex models to accurately constrain the physics of GRB afterglows.

We introduce the `Illustris zoom simulation project', which allows the study of selected galaxies forming in the $\Lambda$CDM cosmology with a 40 times better mass resolution than in the parent large-scale hydrodynamical Illustris simulation. We here focus on the starburst properties of the gas in four cosmological simulations of major mergers. The galaxies in our high-resolution zoom runs exhibit a bursty mode of star formation with gas consumption timescales 10 times shorter than for the normal star formation mode. The strong bursts are only present in the simulations with the highest resolution, hinting that a too low resolution is the reason why the original Illustris simulation showed a dearth of starburst galaxies. Very pronounced bursts of star formation occur in two out of four major mergers we study. The high star formation rates, the short gas consumption timescales and the morphology of these systems strongly resemble observed nuclear starbursts. This is the first time that a sample of major mergers is studied through self-consistent cosmological hydrodynamical simulations instead of using isolated galaxy models setup on a collision course. We also study the orbits of the colliding galaxies and find that the starbursting gas preferentially appears in head-on mergers with very high collision velocities. Encounters with large impact parameters do typically not lead to the formation of starbursting gas.

We report the early discovery of the optical afterglow of gamma-ray burst (GRB) 140801A in the 137 deg$^2$ 3-$\sigma$ error-box of the Fermi Gamma-ray Burst Monitor (GBM). MASTER is the only observatory that automatically react to all Fermi alerts. GRB 140801A is one of the few GRBs whose optical counterpart was discovered solely from its GBM localization. The optical afterglow of GRB 140801A was found by MASTER Global Robotic Net 53 sec after receiving the alert, making it the fastest optical detection of a GRB from a GBM error-box. Spectroscopy obtained with the 10.4-m Gran Telescopio Canarias and the 6-m BTA of SAO RAS reveals a redshift of $z=1.32$. We performed optical and near-infrared photometry of GRB 140801A using different telescopes with apertures ranging from 0.4-m to 10.4-m. GRB 140801A is a typical burst in many ways. The rest-frame bolometric isotropic energy release and peak energy of the burst is $E_\mathrm{iso} = 5.54_{-0.24}^{+0.26} \times 10^{52}$ erg and $E_\mathrm{p, rest}\simeq280$ keV, respectively, which is consistent with the Amati relation. The absence of a jet break in the optical light curve provides a lower limit on the half-opening angle of the jet $\theta=6.1$ deg. The observed $E_\mathrm{peak}$ is consistent with the limit derived from the Ghirlanda relation. The joint Fermi GBM and Konus-Wind analysis shows that GRB 140801A could belong to the class of intermediate duration. The rapid detection of the optical counterpart of GRB 140801A is especially important regarding the upcoming experiments with large coordinate error-box areas.

The reionisation of the Universe is thought to have ended around z~6, as inferred from spectroscopy of distant bright background sources, such as quasars (QSO) and gamma-ray burst (GRB) afterglows. Furthermore, spectroscopy of a GRB afterglow provides insight in its host galaxy, which is often too dim and distant to study otherwise. We present the high S/N VLT/X-shooter spectrum of GRB130606A at z=5.913. We aim to measure the degree of ionisation of the IGM between 5.02<z<5.84 and to study the chemical abundance pattern and dust content of its host galaxy. We measured the flux decrement due to absorption at Ly$\alpha$, $\beta$ and $\gamma$ wavelength regions. The hydrogen and metal absorption lines formed in the host galaxy were fitted with Voigt profiles to obtain column densities. Our measurements of the Ly$\alpha$-forest optical depth are consistent with previous measurements of QSOs, but have a much smaller uncertainty. The analysis of the red damping wing yields a neutral fraction $x_{HI}<0.05$ (3$\sigma$). We obtain column density measurements of several elements. The ionisation corrections due to the GRB is estimated to be negligible (<0.03 dex), but larger corrections may apply due to the pre-existing radiation field (up to 0.4 dex based on sub-DLA studies). Our measurements confirm that the Universe is already predominantly ionised over the redshift range probed in this work, but was slightly more neutral at z>5.6. GRBs are useful probes of the ionisation state of the IGM in the early Universe, but because of internal scatter we need a larger statistical sample to draw robust conclusions. The high [Si/Fe] in the host can be due to dust depletion, alpha-element enhancement, or a combination of both. The very high value of [Al/Fe]=2.40+/-0.78 might connected to the stellar population history. We estimate the host metallicity to be -1.7<[M/H]<-0.9 (2%-13% of solar). (trunc.)

Understanding the physical processes that drive star formation is a key challenge for galaxy formation models. In this article we study the tight correlation between the star formation rate (SFR) and stellar mass of galaxies at a given redshift, how halo growth influences star formation, and star formation histories of individual galaxies. We study these topics using Illustris, a state-of-the-art cosmological hydrodynamical simulation of galaxy formation. Illustris reproduces the observed relation (the star formation main sequence; SFMS) between SFR and stellar mass at redshifts z=0 and z=4, but at intermediate redshifts of z~2, the simulated SFMS has a significantly lower normalisation than reported by observations. The scatter in the relation is consistent with the observed scatter. However, the fraction of outliers above the SFR-stellar mass relation in Illustris is less than that observed. Galaxies with halo masses of ~10^{12} solar masses dominate the SFR density of the Universe, in agreement with the results of abundance matching. Furthermore, more-massive galaxies tend to form the bulk of their stars at high redshift, which indicates that `downsizing' occurs in Illustris. We also studied the star formation histories of individual galaxies, including the use of a principal component analysis decomposition. We find that for fixed stellar mass, galaxies that form earlier have more-massive black holes at z=0, indicating that star formation and black hole growth are tightly linked processes in Illustris. While many of the properties of normal star-forming galaxies are well-reproduced in the Illustris simulation, forming a realistic population of starbursts will likely require higher resolution and probably a more sophisticated treatment of star formation and feedback from stars and black holes.

Context. Gamma-ray burst (GRBs) afterglows probe sightlines to star-forming regions in distant star-forming galaxies. Here we present a study of the peculiar afterglow spectrum of the z = 0.889 Swift GRB 140506A. Aims. Our aim is to understand the origin of the very unusual properties of the absorption along the line-of-sight. Methods. We analyse spectroscopic observations obtained with the X-shooter spectrograph mounted on the ESO/VLT at two epochs 8.8 h and 33 h after the burst as well as imaging from the GROND instrument. We also present imaging and spectroscopy of the host galaxy obtained with the Magellan telescope. Results. The underlying afterglow appears to be a typical afterglow of a long-duration GRB. However, the material along the line-of- sight has imprinted very unusual features on the spectrum. Firstly, there is a very broad and strong flux drop below 8000 AA (4000 AA in the rest frame), which seems to be variable between the two spectroscopic epochs. We can reproduce the flux-drops both as a giant 2175 AA extinction bump and as an effect of multiple scattering on dust grains in a dense environment. Secondly, we detect absorption lines from excited H i and He i. We also detect molecular absorption from CH+ . Conclusions. We interpret the unusual properties of these spectra as reflecting the presence of three distinct regions along the line-of-sight: the excited He i absorption originates from an H ii-region, whereas the Balmer absorption must originate from an associated photodissociation region. The strong metal line and molecular absorption and the dust extinction must originate from a third, cooler region along the line-of-sight. The presence of (at least) three separate regions is reflected in the fact that the different absorption components have different velocities relative to the systemic redshift of the host galaxy.

Dark matter halos are built from accretion and merging. During merging some of the dark matter particles may be ejected with velocities higher than the escape velocity. We use both N-body simulations and single-particle smooth-field simulations to demonstrate that rapid changes to the mean field potential are responsible for such ejection, and in particular that dynamical friction plays no significant role in it. Studying a range of minor mergers, we find that typically between 5-15% of the particles from the smaller of the two merging structures are ejected. We also find that the ejected particles originate essentially from the small halo, and more specifically are particles in the small halo which pass later through the region in which the merging occurs.

Observations of the afterglows of long gamma-ray bursts (GRBs) allow the study of star-forming galaxies across most of cosmic history. Here we present observations of GRB 111008A from which we can measure metallicity, chemical abundance patterns, dust-to-metals ratio and extinction of the GRB host galaxy at z=5.0. The host absorption system is a damped Lyman-alpha absorber (DLA) with a very large neutral hydrogen column density of log N(HI)/cm^(-2) = 22.30 +/- 0.06, and a metallicity of [S/H]= -1.70 +/- 0.10. It is the highest redshift GRB with such a precise metallicity measurement. The presence of fine-structure lines confirms the z=5.0 system as the GRB host galaxy, and makes this the highest redshift where Fe II fine-structure lines have been detected. The afterglow is mildly reddened with A_V = 0.11 +/- 0.04 mag, and the host galaxy has a dust-to-metals ratio which is consistent with being equal to or lower than typical values in the Local Group.

At low redshift, a handful of gamma-ray bursts (GRBs) have been discovered with peak luminosities ($L_{\rm iso} < 10^{48.5}~\rm{erg\,s}^{-1}$) substantially lower than the average of the more distant ones ($L_{\rm iso} > 10^{49.5}~\rm{erg\,s}^{-1}$). The properties of several low-luminosity (low-$L$) GRBs indicate that they can be due to shock break-out, as opposed to the emission from ultrarelativistic jets. Owing to this, it is highly debated how both populations are connected, and whether there is a continuum between them. The burst at redshift $z=0.283$ from 2012 April 22 is one of the very few examples of intermediate-$L$ GRBs with a $\gamma$-ray luminosity of $L\sim10^{48.9}~\rm{erg\,s}^{-1}$ that have been detected up to now. Together with the robust detection of its accompanying supernova SN 2012bz, it has the potential to answer important questions on the origin of low- and high-$L$ GRBs and the GRB-SN connection. We carried out a spectroscopy campaign using medium- and low-resolution spectrographs at 6--10-m class telescopes, covering the time span of 37.3 days, and a multi-wavelength imaging campaign from radio to X-ray energies over a duration of $\sim270$ days. Furthermore, we used a tuneable filter centred at H$\alpha$ to map star formation in the host galaxy and the surrounding galaxies. We used these data to extract and model the properties of different radiation components and incorporate spectral-energy-distribution fitting techniques to extract the properties of the host galaxy. Modelling the light curve and spectral energy distribution from the radio to the X-rays revealed the blast-wave to expand with an initial Lorentz factor of $\Gamma_0\sim60$, low for a high-$L$ GRB, and that the afterglow had an exceptional low peak luminosity-density of $\lesssim2\times10^{30}~\rm{erg\,s}^{-1}\,\rm{Hz}^{-1}$ in the sub-mm. [Abridged]

Short duration gamma-ray bursts (SGRBs) are thought to be related to the violent merger of compact objects, such as neutron stars or black holes, which makes them promising sources of gravitational waves. The detection of a 'kilonova'-like signature associated to the Swift-detected GRB 130603B has suggested that this event is the result of a compact object merger. Our knowledge on SGRB has been, until now, mostly based on the absence of supernova signatures and the analysis of the host galaxies to which they cannot always be securely associated. Further progress has been significantly hampered by the faintness and rapid fading of their optical counterparts (afterglows), which has so far precluded spectroscopy of such events. Afterglow spectroscopy is the key tool to firmly determine the distance at which the burst was produced, crucial to understand its physics, and study its local environment. Here we present the first spectra of a prototypical SGRB afterglow in which both absorption and emission features are clearly detected. Together with multiwavelength photometry we study the host and environment of GRB 130603B. From these spectra we determine the redshift of the burst to be z = 0.3565+/-0.0002, measure rich dynamics both in absorption and emission, and a substantial line of sight extinction of A_V = 0.86+/-0.15 mag. The GRB was located at the edge of a disrupted arm of a moderately star forming galaxy with near-solar metallicity. Unlike for most long GRBs (LGRBs), N_HX / A_V is consistent with the Galactic ratio, indicating that the explosion site differs from those found in LGRBs. The merger is not associated with the most star-forming region of the galaxy; however, it did occur in a dense region, implying a rapid merger or a low natal kick velocity for the compact object binary.

With optical and near-infrared follow-up observations of gamma-ray bursts (GRBs) it is possible to study the chemical enrichment of galaxies at high redshift. Especially interesting are measurements of the metallicity of GRB host with $z\gtrsim 4.7$, since it has recently been shown that damped Lymann-$\alpha$ absorbers have a rapid decrease in metallicity at this redshift (Rafelski et al. 2013). In this paper measurements of the amount of metals and dust in the host galaxy of GRB 111008A at $z=5.0$ is presented, and it is discussed how observations of more high-$z$ GRB host galaxies can give clues about dust formation mechanisms and the distribution of metals in the early Universe.

Long-duration gamma-ray bursts (GRBs) at z < 1 are in most cases found to be accompanied by bright, broad-lined Type Ic supernovae (SNe Ic-BL). The highest-energy GRBs are mostly located at higher redshifts, where the associated SNe are hard to detect observationally. Here we present early and late observations of the optical counterpart of the very energetic GRB 130427A. Despite its moderate redshift z = 0.3399+/-0.0002, GRB 130427A is at the high end of the GRB energy distribution, with an isotropic-equivalent energy release of Eiso ~ 9.6x10^53 erg, more than an order of magnitude more energetic than other GRBs with spectroscopically confirmed SNe. In our dense photometric monitoring, we detect excess flux in the host-subtracted r-band light curve, consistent with what expected from an emerging SN, ~0.2 mag fainter than the prototypical SN 1998bw. A spectrum obtained around the time of the SN peak (16.7 days after the GRB) reveals broad undulations typical of SNe Ic-BL, confirming the presence of a SN, designated SN 2013cq. The spectral shape and early peak time are similar to those of the high expansion velocity SN 2010bh associated with GRB 100316D. Our findings demonstrate that high-energy long-duration GRBs, commonly detected at high redshift, can also be associated with SNe Ic-BL, pointing to a common progenitor mechanism.

Dark matter haloes from cosmological N-body simulations typically have triaxial shapes and anisotropic velocity distributions. Recently it has been shown that the velocity anisotropy, beta, of cosmological haloes and major merger remnants depends on direction in such a way that beta is largest along the major axis and smallest along the minor axis. In this work we use a wide range of non-cosmological N-body simulations to examine halo shapes and direction-dependence of velocity anisotropy profiles. For each of our simulated haloes we define 48 cones pointing in different directions, and from the particles inside each cone we compute velocity anisotropy profiles. We find that elongated haloes can have very distinct velocity anisotropies. We group the behaviour of haloes into three different categories, that range from spherically symmetric profiles to a much more complex behaviour, where significant differences are found for beta along the major and minor axes. We encourage future studies of velocity anisotropies in haloes from cosmological simulations to calculate beta-profiles in cones, since it reveals information, which is hidden from a spherically averaged profile. Finally, we show that spherically averaged profiles often obey a linear relation between beta and the logarithmic density slope in the inner parts of haloes, but this relation is not necessarily obeyed, when properties are calculated in cones.

In several gamma-ray bursts (GRBs) excess emission, in addition to the standard synchrotron afterglow spectrum, has been discovered in the early time X-ray observations. It has been proposed that this excess comes from black body emission, which may be related to the shock break-out of a supernova in the GRBs progenitor star. This hypothesis is supported by the discovery of excess emission in several GRBs with an associated supernova. Using mock spectra we show that it is only likely to detect such a component, similar to the one proposed in GRB 101219B, at low redshift and in low absorption environments. We also perform a systematic search for black body components in all the GRBs observed with the Swift satellite and find six bursts (GRB 061021, 061110A, 081109, 090814A, 100621A and 110715A) with possible black body components. Under the assumption that their excess emission is due to a black body component we present radii, temperatures and luminosities of the emitting components. We also show that detection of black body components only is possible in a fraction of the Swift bursts.

(Abridged) Long gamma-ray burst (GRB) host galaxies might open a short-cut to the characteristics of typical star-forming galaxies throughout the history of the Universe. Due to the absence of near-infrared (NIR) spectroscopy, however, detailed investigations, specifically a determination of the gas-phase metallicity of gamma-ray burst hosts, was largely limited to redshifts z < 1 to date. Here, we report observations of the galaxy hosting GRB 080605 at z = 1.64 using optical/NIR spectroscopy and high-resolution HST/WFC3 imaging. We avail of VLT/X-shooter spectroscopy to measure the metallicity, electron density, star-formation rate (SFR), and reddening of the host. Specifically, we use different strong-line diagnostics based on [N II] to robustly measure the gas-phase metallicity for the first time for a GRB host at this redshift. The host of the energetic (E_iso ~ 2 x 10^53 erg) GRB 080605 is a morphologically complex, vigorously star-forming galaxy with an H\alpha-derived SFR of 31 M_sun/yr. Its ISM is significantly enriched with metals with an oxygen abundance between 8.3 and 8.6 depending on the adopted strong-line calibrator. This corresponds to values in the range of 0.4-0.8 times the solar value. For its measured stellar mass of 8 x 10^9 M_sun and SFR this value is consistent with the fundamental metallicity relation defined by star-forming field galaxies. Our observations directly illustrate that GRB hosts are not necessarily metal-poor, both on absolute scales as well as relative to their stellar mass and SFR. GRB hosts could thus be fair tracers of the population of ordinary star-forming galaxies as a whole at high redshift.

We present simultaneous optical and near-infrared (NIR) spectroscopy of 19 Swift GRB host galaxies with VLT/X-shooter with the aim of measuring their redshifts. Galaxies were selected from The Optically Unbiased GRB Host (TOUGH) survey (15 of the 19 galaxies) or because they hosted GRBs without a bright optical afterglow. Here, we provide emission-line redshifts for 13 of the observed galaxies with brightnesses between F606W > 27 mag and R=22.9 mag (median R=24.6 mag). The median redshift is z=2.1 for all, and z=2.3 for the TOUGH hosts. Our new data significantly improve the redshift completeness of the TOUGH survey, which now stands at 77% (53 out of 69 GRBs). They furthermore provide accurate redshifts for nine prototype-dark GRBs (e.g., GRBs 071021 at z=2.452 and 080207 at z=2.086), which are exemplary of GRBs where redshifts are challenging to obtain via afterglow spectroscopy. This establishes X-shooter spectroscopy as an efficient tool for redshift determination of faint, star-forming, high-redshift galaxies such as GRB hosts. It is hence a further step towards removing the bias in GRB samples that is caused by optically-dark events, and provides the basis for a better understanding of the conditions in which GRBs form. The distribution of column densities as measured from X-ray data (N_{H,X}), for example, is closely related to the darkness of the afterglow and skewed towards low N_{H, X} values in samples that are dominated by bursts with bright optical afterglows.

Dark matter haloes in cosmological N-body simulations are affected by processes such as mergers, accretion and the gravitational interaction with baryonic matter. Typically the analysis of dark matter haloes is performed in spherical or elliptical bins and the velocity distributions are often assumed to be constant within those bins. However, the velocity anisotropy, which describes differences between the radial and tangential velocity dispersion, has recently been show to have a strong dependence on direction in the triaxial halos formed in cosmological simulations. In this study we derive properties of particles in cones parallel or perpendicular to the collision axis of merger remnants. We find that the velocity anisotropy has a strong dependence on direction. The finding that the direction-dependence of the velocity anisotropy of a halo depends on the merger history, explains the existence of such trends in cosmological simulations. It also explains why a large diversity is seen in the velocity anisotropy profiles in the outer parts of high-resolution simulations of cosmological haloes.

All dark matter structures appear to follow a set of universalities, such as phase-space density or velocity anisotropy profiles, however, the origin of these universalities remains a mystery. Any equilibrated dark matter structure can be fully described by two functions, namely the radial and the tangential velocity distribution functions (VDF), and when we will understand these two then we will understand all the observed universalities. Here we demonstrate that if we know the radial VDF, then we can derive and understand the tangential VDF. This is based on simple dynamical arguments about properties of collisionless systems. We use a range of controlled numerical simulations to demonstrate the accuracy of this result. We therefore boil the question of the dark matter structural properties down to understanding the radial VDF.

Recently, the number of detected galaxy counterparts of z > 2 Damped Lyman-alpha Absorbers in QSO spectra has increased substantially so that we today have a sample of 10 detections. M{\o}ller et al. in 2004 made the prediction, based on a hint of a luminosity-metallicity relation for DLAs, that HI size should increase with increasing metallicity. In this paper we investigate the distribution of impact parameter and metallicity that would result from the correlation between galaxy size and metallicity. We compare our observations with simulated data sets given the relation of size and metallicity. The observed sample presented here supports the metallicity-size prediction: The present sample of DLA galaxies is consistent with the model distribution. Our data also show a strong relation between impact parameter and column density of HI. We furthermore compare the observations with several numerical simulations and demonstrate that the observations support a scenario where the relation between size and metallicity is driven by feedback mechanisms controlling the star-formation efficiency and outflow of enriched gas.

We report on the spectroscopic detection of supernova SN 2010ma associated with the long gamma-ray burst GRB 101219B. We observed the optical counterpart of the GRB on three nights with the X-shooter spectrograph at the VLT. From weak absorption lines, we measure a redshift of z = 0.55. The first epoch UV--near-infrared afterglow spectrum, taken 11.6 hr after the burst, is well fit by a power law consistent with the slope of the X-ray spectrum. The second and third epoch spectra (obtained 16.4 and 36.7 days after the burst), however, display clear bumps closely resembling those of the broad-lined type-Ic SN 1998bw if placed at z=0.55. Apart from demonstrating that spectroscopic SN signatures can be observed for GRBs at these large distances, our discovery makes a step forward in establishing a general connection between GRBs and SNe. In fact, unlike most previous unambiguous GRB-associated SNe, GRB 101219B has a large gamma-ray energy (E_iso = 4.2 x 10^{51} erg), a bright afterglow, and obeys the "Amati" relation, thus being fully consistent with the cosmological population of GRBs.

Cosmological simulations of dark matter structures have identified a set of universal profiles, and similar characteristics have been seen in non-cosmological simulations. It has therefore been speculated whether these profiles of collisionless systems relate to accretion and merger history, or if there is an attractor for the dark matter systems. Here we identify such a 1-dimensional attractor in the 3-dimensional space spanned by the 2 radial slopes of the density and velocity dispersion, and the velocity anisotropy. This attractor effectively removes one degree of freedom from the Jeans equation. It also allows us to speculate on a new fluid interpretation for the Jeans equation, with an effective polytropic index for the dark matter particles between 1/2 and 3/4. If this attractor solution holds for other collisionless structures, then it may hold the key to break the mass-anisotropy degeneracy, which presently prevents us from measuring the mass profiles in dwarf galaxies uniquely.