Showing 1–50 of 59 results for author: Sparre, M

Comparing the interstellar and circumgalactic origin of gas in the tails of jellyfish galaxies

Authors: Martin Sparre, Christoph Pfrommer, Ewald Puchwein

Abstract: Simulations and observations find long tails in jellyfish galaxies, which are commonly thought to originate from ram-pressure stripped gas of the interstellar medium (ISM) in the immediate galactic wake. While at larger distances from the galaxy, they have been claimed to form in situ owing to thermal instability and fast radiative cooling of mixed ISM and intracluster medium (ICM). In this paper,… ▽ More Simulations and observations find long tails in jellyfish galaxies, which are commonly thought to originate from ram-pressure stripped gas of the interstellar medium (ISM) in the immediate galactic wake. While at larger distances from the galaxy, they have been claimed to form in situ owing to thermal instability and fast radiative cooling of mixed ISM and intracluster medium (ICM). In this paper, we use magneto-hydrodynamical windtunnel simulations of a galaxy with the arepo code to study the origin of gas in the tails of jellyfish galaxies. To this end, we model the galaxy orbit in a cluster by accounting for a time-varying galaxy velocity, ICM density and turbulent magnetic field. Tracking gas flows between the ISM, the circumgalactic medium (CGM) and the ICM, we find contrary to popular opinion that the majority of gas in the tail originated in the CGM. Prior to the central passage of the jellyfish galaxy in the cluster, the CGM is directly transported to the clumpy jellyfish tail that has been shattered into small cloudlets. After the central cluster passage, gas in the tail originates both from the initial ISM and the CGM, but that from the latter was accreted to the galactic ISM before being ram-pressure stripped to form filamentary tentacles in the tail. Our simulation shows a declining gas metallicity in the tail as a function of downstream distance from the galaxy. We conclude that the CGM plays an important role in shaping the tails of jellyfish galaxies. △ Less

Submitted 1 May, 2024; originally announced May 2024.

Comments: 9 pages, 5 figures, submitted to A&A (Main Journal)

The magnetised and thermally unstable tails of jellyfish galaxies

Authors: Martin Sparre, Christoph Pfrommer, Ewald Puchwein

Abstract: Jellyfish galaxies are promising laboratories for studying radiative cooling and magnetic fields in multiphase gas flows. Their long, dense tails are observed to be magnetised, and they extend up to 100 kpc into the intracluster medium (ICM), suggesting that their gas is thermally unstable so that the cold gas mass grows with time rather than being fully dissolved in the hot wind as a result of hy… ▽ More Jellyfish galaxies are promising laboratories for studying radiative cooling and magnetic fields in multiphase gas flows. Their long, dense tails are observed to be magnetised, and they extend up to 100 kpc into the intracluster medium (ICM), suggesting that their gas is thermally unstable so that the cold gas mass grows with time rather than being fully dissolved in the hot wind as a result of hydrodynamical interface instabilities. In this paper we use the AREPO code to perform magnetohydrodynamical windtunnel simulations of a jellyfish galaxy experiencing ram-pressure stripping by interacting with an ICM wind. The ICM density, temperature and velocity that the galaxy encounters are time-dependent and comparable to what a real jellyfish galaxy experiences while orbiting the ICM. In simulations with a turbulent magnetised wind we reproduce observations, which show that the magnetic field is aligned with the jellyfish tails. During the galaxy infall into the cluster with a near edge-on geometry, the gas flow in the tail is fountain-like, implying preferential stripping of gas where the rotational velocity vectors add up with the ram pressure while fall-back occurs in the opposite case. Hence, the tail velocity shows a memory of the rotation pattern of the disc. At the time of the nearest cluster passage, ram-pressure stripping is so strong that the fountain flow is destroyed and instead the tail is dominated by removal of gas. We show that gas in the tail is very fragmentative, which is a prediction of shattering due to radiative cooling. △ Less

Submitted 20 November, 2023; v1 submitted 9 November, 2023; originally announced November 2023.

Comments: 14 pages, 7 figures, Accepted for publication in MNRAS

Unveiling the distant Universe: Characterizing $z\ge9$ Galaxies in the first epoch of COSMOS-Web

Authors: Maximilien Franco, Hollis B. Akins, Caitlin M. Casey, Steven L. Finkelstein, Marko Shuntov, Katherine Chworowsky, Andreas L. Faisst, Seiji Fujimoto, Olivier Ilbert, Anton M. Koekemoer, Daizhong Liu, Christopher C. Lovell, Claudia Maraston, Henry Joy McCracken, Jed McKinney, Brant E. Robertson, Micaela B. Bagley, Jaclyn B. Champagne, Olivia R. Cooper, Xuheng Ding, Nicole E. Drakos, Andrea Enia, Steven Gillman, Christopher C. Hayward, Michaela Hirschmann , et al. (25 additional authors not shown)

Abstract: We report the identification of 15 galaxy candidates at $z\ge9$ using the initial COSMOS-Web JWST observations over 77 arcmin$^2$ through four NIRCam filters (F115W, F150W, F277W, F444W) with an overlap with MIRI (F770W) of 8.7 arcmin$^2$. We fit the sample using several publicly-available SED fitting and photometric redshift codes and determine their redshifts between $z=9.3$ and $z=10.9$ (… ▽ More We report the identification of 15 galaxy candidates at $z\ge9$ using the initial COSMOS-Web JWST observations over 77 arcmin$^2$ through four NIRCam filters (F115W, F150W, F277W, F444W) with an overlap with MIRI (F770W) of 8.7 arcmin$^2$. We fit the sample using several publicly-available SED fitting and photometric redshift codes and determine their redshifts between $z=9.3$ and $z=10.9$ ($\langle z\rangle=10.0$), UV-magnitudes between M$_{\rm UV}$ = $-$21.2 and $-$19.5 (with $\langle $M$_{\rm UV}\rangle=-20.2$) and rest-frame UV slopes ($\langle β\rangle=-2.4$). These galaxies are, on average, more luminous than most $z\ge9$ candidates discovered by JWST so far in the literature, while exhibiting similar blue colors in their rest-frame UV. The rest-frame UV slopes derived from SED-fitting are blue ($β\sim$[$-$2.0, $-$2.7]) without reaching extremely blue values as reported in other recent studies at these redshifts. The blue color is consistent with models that suggest the underlying stellar population is not yet fully enriched in metals like similarly luminous galaxies in the lower redshift Universe. The derived stellar masses with $\langle \log_{\rm 10} ($M$_\star/$M$_\odot)\rangle\approx8-9$ are not in tension with the standard $Λ$CDM model and our measurement of the volume density of such UV luminous galaxies aligns well with previously measured values presented in the literature at $z\sim9-10$. Our sample of galaxies, although compact, are significantly resolved. △ Less

Submitted 1 August, 2023; originally announced August 2023.

Comments: Submitted to ApJ

Structure, Kinematics, and Observability of the Large Magellanic Cloud's Dynamical Friction Wake in Cold vs. Fuzzy Dark Matter

Authors: Hayden R. Foote, Gurtina Besla, Philip Mocz, Nicolás Garavito-Camargo, Lachlan Lancaster, Martin Sparre, Emily C. Cunningham, Mark Vogelsberger, Facundo A. Gómez, Chervin F. P. Laporte

Abstract: The Large Magellanic Cloud (LMC) will induce a dynamical friction (DF) wake on infall to the Milky Way (MW). The MW's stellar halo will respond to the gravity of the LMC and the dark matter (DM) wake, forming a stellar counterpart to the DM wake. This provides a novel opportunity to constrain the properties of the DM particle. We present a suite of high-resolution, windtunnel-style simulations of… ▽ More The Large Magellanic Cloud (LMC) will induce a dynamical friction (DF) wake on infall to the Milky Way (MW). The MW's stellar halo will respond to the gravity of the LMC and the dark matter (DM) wake, forming a stellar counterpart to the DM wake. This provides a novel opportunity to constrain the properties of the DM particle. We present a suite of high-resolution, windtunnel-style simulations of the LMC's DF wake that compare the structure, kinematics, and stellar tracer response of the DM wake in cold DM (CDM), with and without self-gravity, vs. fuzzy DM (FDM) with $m_a = 10^{-23}$ eV. We conclude that the self-gravity of the DM wake cannot be ignored. Its inclusion raises the wake's density by $\sim 10\%$, and holds the wake together over larger distances ($\sim$ 50 kpc) than if self-gravity is ignored. The DM wake's mass is comparable to the LMC's infall mass, meaning the DM wake is a significant perturber to the dynamics of MW halo tracers. An FDM wake is more granular in structure and is $\sim 20\%$ dynamically colder than a CDM wake, but with comparable density. The granularity of an FDM wake increases the stars' kinematic response at the percent level compared to CDM, providing a possible avenue of distinguishing a CDM vs. FDM wake. This underscores the need for kinematic measurements of stars in the stellar halo at distances of 70-100 kpc. △ Less

Submitted 8 September, 2023; v1 submitted 30 June, 2023; originally announced July 2023.

Comments: 38 pages, 31 figures. Updated to published version

Journal ref: ApJ 954 163 (2023)

Resolving galactic-scale obscuration of X-ray AGN at $z\gtrsim1$ with COSMOS-Web

Authors: John D. Silverman, Vincenzo Mainieri, Xuheng Ding, Daizhong Liu, Knud Jahnke, Michaela Hirschmann, Jeyhan Kartaltepe, Erini Lambrides, Masafusa Onoue, Benny Trakhtenbrot, Eleni Vardoulaki, Angela Bongiorno, Caitlin Casey, Francesca Civano, Andreas Faisst, Maximilien Franco, Steven Gillman, Ghassem Gozaliasl, Christopher C. Hayward, Anton M. Koekemoer, Vasily Kokorev, Georgios Magdis, Stefano Marchesi, Robert Michael Rich, Martin Sparre , et al. (3 additional authors not shown)

Abstract: A large fraction of the accreting supermassive black hole population is shrouded by copious amounts of gas and dust, particularly in the distant ($z\gtrsim1$) Universe. While much of the obscuration is attributed to a parsec-scale torus, there is a known contribution from the larger-scale host galaxy. Using JWST/NIRCam imaging from the COSMOS-Web survey, we probe the galaxy-wide dust distribution… ▽ More A large fraction of the accreting supermassive black hole population is shrouded by copious amounts of gas and dust, particularly in the distant ($z\gtrsim1$) Universe. While much of the obscuration is attributed to a parsec-scale torus, there is a known contribution from the larger-scale host galaxy. Using JWST/NIRCam imaging from the COSMOS-Web survey, we probe the galaxy-wide dust distribution in X-ray selected AGN up to $z\sim2$. Here, we focus on a sample of three AGNs with their host galaxies exhibiting prominent dust lanes, potentially due to their edge-on alignment. These represent 27% (3 out of 11 with early NIRCam data) of the heavily obscured ($N_H>10^{23}$ cm$^{-2}$) AGN population. With limited signs of a central AGN in the optical and near-infrared, the NIRCam images are used to produce reddening maps $E(B-V)$ of the host galaxies. We compare the mean central value of $E(B-V)$ to the X-ray obscuring column density along the line-of-sight to the AGN ($N_H\sim10^{23-23.5}$ cm$^{-2}$). We find that the extinction due to the host galaxy is present ($0.6\lesssim E(B-V) \lesssim 0.9$; $1.9 \lesssim A_V \lesssim 2.8$) and significantly contributes to the X-ray obscuration at a level of $N_H\sim10^{22.5}$ cm$^{-2}$ assuming an SMC gas-to-dust ratio which amounts to $\lesssim$30% of the total obscuring column density. These early results, including three additional cases from CEERS, demonstrate the ability to resolve such dust structures with JWST and separate the different circumnuclear and galaxy-scale obscuring structures. △ Less

Submitted 5 June, 2023; originally announced June 2023.

Comments: 12 pages, 8 figures, Accepted for publication in ApJL

The impact of magnetic fields on cosmological galaxy mergers -- II. Modified angular momentum transport and feedback

Authors: Joseph Whittingham, Martin Sparre, Christoph Pfrommer, Rüdiger Pakmor

Abstract: The role of magnetic fields in galaxy evolution is still an unsolved question in astrophysics. We have previously shown that magnetic fields play a crucial role in major mergers between disc galaxies; in hydrodynamic simulations of such mergers, the Auriga model produces compact remnants with a distinctive bar and ring morphology. In contrast, in magnetohydrodynamic (MHD) simulations, remnants for… ▽ More The role of magnetic fields in galaxy evolution is still an unsolved question in astrophysics. We have previously shown that magnetic fields play a crucial role in major mergers between disc galaxies; in hydrodynamic simulations of such mergers, the Auriga model produces compact remnants with a distinctive bar and ring morphology. In contrast, in magnetohydrodynamic (MHD) simulations, remnants form radially-extended discs with prominent spiral arm structure. In this paper, we analyse a series of cosmological "zoom-in" simulations of major mergers and identify exactly \textit{how} magnetic fields are able to alter the outcome of the merger. We find that magnetic fields modify the transport of angular momentum, systematically hastening the merger progress. The impact of this altered transport depends on the orientation of the field, with a predominantly non-azimuthal (azimuthal) orientation increasing the central baryonic concentration (providing support against collapse). Both effects act to suppress an otherwise existent bar-instability, which in turn leads to a fundamentally different morphology and manifestation of feedback. We note, in particular, that stellar feedback is substantially less influential in MHD simulations, which allows for the later accretion of higher angular momentum gas and the subsequent rapid radial growth of the remnant disc. A corollary of the increased baryonic concentration in MHD simulations is that black holes are able to grow twice as large, although this turns out to have little impact on the remnant's development. Our results show that galaxy evolution cannot be modelled correctly without including magnetic fields. △ Less

Submitted 26 September, 2023; v1 submitted 30 January, 2023; originally announced January 2023.

Comments: 22 pages, 12 figures. Published in MNRAS on 07 September 2023

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 526, Issue 1, November 2023, Pages 224-245

COSMOS-Web: An Overview of the JWST Cosmic Origins Survey

Authors: Caitlin M. Casey, Jeyhan S. Kartaltepe, Nicole E. Drakos, Maximilien Franco, Santosh Harish, Louise Paquereau, Olivier Ilbert, Caitlin Rose, Isabella G. Cox, James W. Nightingale, Brant E. Robertson, John D. Silverman, Anton M. Koekemoer, Richard Massey, Henry Joy McCracken, Jason Rhodes, Hollis B. Akins, Aristeidis Amvrosiadis, Rafael C. Arango-Toro, Micaela B. Bagley, Angela Bongiorno, Peter L. Capak, Jaclyn B. Champagne, Nima Chartab, Oscar A. Chavez Ortiz , et al. (60 additional authors not shown)

Abstract: We present the survey design, implementation, and outlook for COSMOS-Web, a 255 hour treasury program conducted by the James Webb Space Telescope in its first cycle of observations. COSMOS-Web is a contiguous 0.54 deg$^2$ NIRCam imaging survey in four filters (F115W, F150W, F277W, and F444W) that will reach 5$σ$ point source depths ranging $\sim$27.5-28.2 magnitudes. In parallel, we will obtain 0.… ▽ More We present the survey design, implementation, and outlook for COSMOS-Web, a 255 hour treasury program conducted by the James Webb Space Telescope in its first cycle of observations. COSMOS-Web is a contiguous 0.54 deg$^2$ NIRCam imaging survey in four filters (F115W, F150W, F277W, and F444W) that will reach 5$σ$ point source depths ranging $\sim$27.5-28.2 magnitudes. In parallel, we will obtain 0.19 deg$^2$ of MIRI imaging in one filter (F770W) reaching 5$σ$ point source depths of $\sim$25.3-26.0 magnitudes. COSMOS-Web will build on the rich heritage of multiwavelength observations and data products available in the COSMOS field. The design of COSMOS-Web is motivated by three primary science goals: (1) to discover thousands of galaxies in the Epoch of Reionization ($6<z<11$) and map reionization's spatial distribution, environments, and drivers on scales sufficiently large to mitigate cosmic variance, (2) to identify hundreds of rare quiescent galaxies at $z>4$ and place constraints on the formation of the Universe's most massive galaxies ($M_\star>10^{10}$\,M$_\odot$), and (3) directly measure the evolution of the stellar mass to halo mass relation using weak gravitational lensing out to $z\sim2.5$ and measure its variance with galaxies' star formation histories and morphologies. In addition, we anticipate COSMOS-Web's legacy value to reach far beyond these scientific goals, touching many other areas of astrophysics, such as the identification of the first direct collapse black hole candidates, ultracool sub-dwarf stars in the Galactic halo, and possibly the identification of $z>10$ pair-instability supernovae. In this paper we provide an overview of the survey's key measurements, specifications, goals, and prospects for new discovery. △ Less

Submitted 8 March, 2023; v1 submitted 14 November, 2022; originally announced November 2022.

Comments: 46 pages, 16 figures, ApJ accepted

Kinematics of the Local Group gas and galaxies in the Hestia simulations

Authors: Luis Biaus, Sebastián E. Nuza, Philipp Richter, Martin Sparre, Cecilia Scannapieco, Mitali Damle, Jenny G. Sorce, Robert J. J. Grand, Elmo Tempel, Noam I. Libeskind, Maan H. Hani

Abstract: We investigate the kinematic properties of gas and galaxies in the Local Group (LG) using high-resolution simulations performed by the {\sc Hestia} (High-resolution Environmental Simulations of The Immediate Area) collaboration. Our simulations include the correct cosmography surrounding LG-like regions consisting of two main spiral galaxies of $\sim 10^{12}$~M$_\odot$, their satellites and minor… ▽ More We investigate the kinematic properties of gas and galaxies in the Local Group (LG) using high-resolution simulations performed by the {\sc Hestia} (High-resolution Environmental Simulations of The Immediate Area) collaboration. Our simulations include the correct cosmography surrounding LG-like regions consisting of two main spiral galaxies of $\sim 10^{12}$~M$_\odot$, their satellites and minor isolated galaxies, all sharing the same large-scale motion within a volume of a few Mpc. We characterise the gas and galaxy kinematics within the simulated LGs, from the perspective of the Sun, to compare with observed trends from recent HST/COS absorption-line observations and LG galaxy data. To analyse the velocity pattern of LG gas and galaxies seen in the observational data, we build sky maps from the local standard of rest, and the galactic and local group barycentre frames. Our findings show that the establishment of a radial velocity dipole at low/high latitudes, near the preferred barycentre direction, is a natural outcome of simulation kinematics for material {\it outside} the Milky Way virial radius after removing galaxy rotation when the two main LG galaxies are approaching. Our results favour a scenario where gas and galaxies stream towards the LG barycentre producing a velocity dipole resembling observations. While our study shows in a qualitative way the global matter kinematics in the LG as part of its on-going assembly, quantitative estimates of gas-flow rates and physical conditions of the LG gas have to await a more detailed modeling of the ionization conditions, which will be presented in a follow-up paper. △ Less

Submitted 27 October, 2022; originally announced October 2022.

Comments: 13 pages, 9 figures, accepted for publication in MNRAS

The stellar halo in Local Group Hestia simulations III. Chemical abundance relations for accreted and in-situ stars

Authors: Sergey Khoperskov, Ivan Minchev, Noam Libeskind, Vasily Belokurov, Matthias Steinmetz, Facundo A. Gomez, Robert J. J. Grand, Yehuda Hoffman, Alexander Knebe, Jenny G. Sorce, Martin Sparre, Elmo Tempel, Mark Vogelsberger

Abstract: Since the chemical abundances of stars are the fossil records of the physical conditions in galaxies, they provide the key information for recovering the assembly history of galaxies. In this work, we explore the chemo-chrono-kinematics of accreted and in-situ stars, by analyzing six M31/MW analogues from the HESTIA suite of cosmological hydrodynamics zoom-in simulations of the Local Group. We fou… ▽ More Since the chemical abundances of stars are the fossil records of the physical conditions in galaxies, they provide the key information for recovering the assembly history of galaxies. In this work, we explore the chemo-chrono-kinematics of accreted and in-situ stars, by analyzing six M31/MW analogues from the HESTIA suite of cosmological hydrodynamics zoom-in simulations of the Local Group. We found that the merger debris are chemically distinct from the survived dwarf galaxies. The mergers debris have abundances expected for stars originating from dwarfs that had their star formation activity quenched at early times. Accreted stellar haloes, including individual debris, reveal abundance gradients in the ELz, where the most metal-rich stars have formed in the inner parts of the disrupted systems before the merger and mainly contribute to the central regions of the hosts. Therefore, we suggest that abundance measurements in the inner MW will allow constraining better the parameters of building blocks of the MW stellar halo. The MDFs of the individual debris show several peaks and the majority of debris have lower metallicity than the in-situ stars for Lz>0, while non-rotating and retrograde accreted stars are similar to the in-situ. Prograde accreted stars show a prominent knee in the [Fe/H]-[Mg/Fe] plane while the retrograde stars typically deposit to a high-[Mg/Fe] sequence. We found that the metal-poor stars ([Fe/H]<-1) of the HESTIA galaxies exhibit between zero to 80 km/s net rotation which is consistent with the Aurora population. At higher metallicities, we detect a sharp transition (spin-up) from the turbulent phase to a disk-like rotation. Mergers debris are similar in the [Fe/H]-[Mg/Fe] plane. However, combining a set of abundances allows to capture chemical patterns corresponding to different debris, which are the most prominent as a function of stellar age. △ Less

Submitted 12 September, 2023; v1 submitted 11 June, 2022; originally announced June 2022.

Comments: Paper III in series of III, 21 pages, 17 figures, A&A in press

Journal ref: A&A 677, A91 (2023)

The stellar halo in Local Group Hestia simulations II. The accreted component

Authors: Sergey Khoperskov, Ivan Minchev, Noam Libeskind, Misha Haywood, Paola Di Matteo, Vasily Belokurov, Matthias Steinmetz, Facundo A. Gomez, Robert J. J. Grand, Yehuda Hoffman, Alexander Knebe, Jenny G. Sorce, Martin Sparre, Elmo Tempel, Mark Vogelsberger

Abstract: In the Milky Way, recent progress in the exploration of its assembly history is driven by the tremendous amount of high-quality data delivered by Gaia, which has revealed a number of substructures potentially linked to several ancient accretion events. In this work, aiming to explore the phase-space structure of accreted stars, we analyze six M31/MW analogues from the HESTIA suite of cosmological… ▽ More In the Milky Way, recent progress in the exploration of its assembly history is driven by the tremendous amount of high-quality data delivered by Gaia, which has revealed a number of substructures potentially linked to several ancient accretion events. In this work, aiming to explore the phase-space structure of accreted stars, we analyze six M31/MW analogues from the HESTIA suite of cosmological hydrodynamics zoom-in simulations of the Local Group. We found that all the HESTIA galaxies experience a few dozen mergers but only 1-4 mergers have the stellar mass ratio >0.2 where, depending on the halo definition, the most massive merger contributes from 20% to 70% of the total stellar halo. Individual merger remnants show diverse density distributions at z=0, significantly overlapping with each other and with the in-situ stars in the ELz, UV and RVphi coordinates. The mergers debris often change their position in the ELz with time due to the galactic mass growth and the non-axisymmetry of the potential. In agreement with previous works, we show that even individual merger debris exhibit a number of distinct ELz features. In the UV plane, all HESTIA galaxies reveal radially hot, non-rotating or weakly counter-rotating, Gaia-Sausage-like features. We found an age gradient in Elz space for the individual debris, where the youngest stars, formed in the inner regions of accreting systems, deposit to the innermost regions of the host. The bulk of these stars is being formed during the last stages of accretion, making it possible to date the merger. In actions space (Jr, Jz, Jφ), the mergers debris do not appear as isolated substructures but are instead scattered over a large parameters area and overlapping with the in-situ stars. We also introduce a purely kinematic space (Jz/Jr-eccentricity), where different merger debris can be disentangled better from each other and from the in-situ stars. △ Less

Submitted 12 September, 2023; v1 submitted 9 June, 2022; originally announced June 2022.

Comments: Paper II in a series of III, 22 pages, 16 figures, A&A in press

Journal ref: A&A 677, A90 (2023)

The stellar halo in Local Group Hestia simulations I. The in-situ component and the effect of mergers

Authors: Sergey Khoperskov, Ivan Minchev, Noam Libeskind, Misha Haywood, Paola Di Matteo, Vasily Belokurov, Matthias Steinmetz, Facundo A. Gomez, Robert J. J. Grand, Yehuda Hoffman, Alexander Knebe, Jenny G. Sorce, Martin Sparre, Elmo Tempel, Mark Vogelsberger

Abstract: Theory suggests that mergers play an important role in shaping galactic discs and stellar haloes, which was observationally confirmed in the MW thanks to Gaia data. In this work, aiming to probe the contribution of mergers to the in situ stellar halo formation, we analyse six M31/MW analogues from the HESTIA suite of cosmological hydrodynamical zoom-in simulations of the LG. We found that all the… ▽ More Theory suggests that mergers play an important role in shaping galactic discs and stellar haloes, which was observationally confirmed in the MW thanks to Gaia data. In this work, aiming to probe the contribution of mergers to the in situ stellar halo formation, we analyse six M31/MW analogues from the HESTIA suite of cosmological hydrodynamical zoom-in simulations of the LG. We found that all the HESTIA galaxies experience between one to four mergers with stellar mass ratios between 0.2 and 1 relative to the host at the time of the merger. These significant mergers, with a single exception, happened 7-11Gyr ago. The overall impact of the most massive mergers in HESTIA is clearly seen as a sharp increase in the orbital eccentricity (and a corresponding decrease in the rotational velocity Vphi of pre-existing disc stars of the main progenitor, thus nicely reproducing the Splash-, Plume-like feature that was discovered in the MW. We do find a correlation between mergers and close pericentric passages of massive satellites and bursts of star formation in the in situ component. Massive mergers sharply increase the disc velocity dispersion of the in situ stars; however, the latest significant merger often heats up the disc up to the numbers when the contribution of the previous ones is less prominent in the age-velocity dispersion relation. In HESTIA galaxies, the in situ halo is an important component of the inner stellar halo where its fraction is about 30-40%, while in the outer parts it typically does not exceed ~5% beyond 15 kpc. The simulations suggest that this component of the stellar haloes continues to grow well after mergers conclude; however, the most significant contribution comes from stars that formed recently before the merger. The orbital analysis of the HESTIA galaxies suggests that wedges in Rmax-Zmax space are mainly populated by the stars born between significant mergers. △ Less

Submitted 12 September, 2023; v1 submitted 9 June, 2022; originally announced June 2022.

Comments: Paper I in a series of III, 17 pages, 15 figures, A&A in press

Journal ref: A&A 677, A89 (2023)

Constraining SIDM with halo shapes: revisited predictions from realistic simulations of early-type galaxies

Authors: Giulia Despali, Levi G. Walls, Simona Vegetti, Martin Sparre, Mark Vogelsberger, Jesús Zavala

Abstract: We study the effect of self-interacting dark matter (SIDM) and baryons on the shape of early-type galaxies (ETGs) and their dark matter haloes, comparing them to the predictions of the cold dark matter (CDM) scenario. We use five hydrodynamical zoom-in simulations of haloes hosting ETGs ($M_{\rm vir}\sim 10^{13}M_{\odot}$ and $M_{*}\sim10^{11}M_{\odot}$), simulated in CDM and a SIDM model with con… ▽ More We study the effect of self-interacting dark matter (SIDM) and baryons on the shape of early-type galaxies (ETGs) and their dark matter haloes, comparing them to the predictions of the cold dark matter (CDM) scenario. We use five hydrodynamical zoom-in simulations of haloes hosting ETGs ($M_{\rm vir}\sim 10^{13}M_{\odot}$ and $M_{*}\sim10^{11}M_{\odot}$), simulated in CDM and a SIDM model with constant cross-section of $σ_T/m_χ= 1\ \mathrm{cm}^2 \mathrm{g}^{-1}$. We measure the three-dimensional and projected shapes of the dark matter haloes and their baryonic content using the inertia tensor and compare our measurements to the results of three $HST$ samples of gravitational lenses and $Chandra$ and $XMM-Newton$ X-ray observations. We find that the inclusion of baryons greatly reduces the differences between CDM and a SIDM, together with the ability to draw constraints based on shapes. Lensing measurements reject the predictions of CDM dark-matter-only simulations and prefer one of the hydro scenarios. When we consider the total sample of lenses, observational data prefer the CDM hydro scenario. The shapes of the X-ray emitting gas are compatible with observational results in both hydro runs, with CDM predicting higher elongations only in the very centre. Contrary to previous claims at the scale of elliptical galaxies, we conclude that both CDM and our SIDM model can still explain observed shapes once we include baryons in the simulations. Our results demonstrate that this is essential to derive realistic constraints and that new simulations are needed to confirm and extend our findings. △ Less

Submitted 5 September, 2022; v1 submitted 26 April, 2022; originally announced April 2022.

Comments: 17 pages, 13 figures, revised version accepted for publication in MNRAS

Cold and Hot gas distribution around the Milky-Way-M31 system in the HESTIA simulations

Authors: Mitali Damle, Martin Sparre, Philipp Richter, Maan H. Hani, Sebastián E. Nuza, Christoph Pfrommer, Robert J. J. Grand, Yehuda Hoffman, Noam Libeskind, Jenny G. Sorce, Matthias Steinmetz, Elmo Tempel, Mark Vogelsberger, Peng Wang

Abstract: Recent observations have revealed remarkable insights into the gas reservoir in the circumgalactic medium (CGM) of galaxy haloes. In this paper, we characterise the gas in the vicinity of Milky Way and Andromeda analogues in the HESTIA (High resolution Environmental Simulations of The Immediate Area) suite of constrained Local Group (LG) simulations. The HESTIA suite comprise of a set of three hig… ▽ More Recent observations have revealed remarkable insights into the gas reservoir in the circumgalactic medium (CGM) of galaxy haloes. In this paper, we characterise the gas in the vicinity of Milky Way and Andromeda analogues in the HESTIA (High resolution Environmental Simulations of The Immediate Area) suite of constrained Local Group (LG) simulations. The HESTIA suite comprise of a set of three high-resolution {\sc arepo}-based simulations of the LG, run using the Auriga galaxy formation model. For this paper, we focus only on the $z = 0$ simulation datasets and generate mock skymaps along with a power spectrum analysis to show that the distributions of ions tracing low-temperature gas (HI and SiIII) are more clumpy in comparison to warmer gas tracers (OVI, OVII and OVIII). We compare to the spectroscopic CGM observations of M31 and low-redshift galaxies. HESTIA under-produces the column densities of the M31 observations, but the simulations are consistent with the observations of low-redshift galaxies. A possible explanation for these findings is that the spectroscopic observations of M31 are contaminated by gas residing in the CGM of the Milky Way. △ Less

Submitted 16 March, 2022; v1 submitted 26 January, 2022; originally announced January 2022.

Comments: 22 pages, 9 figures, revised version, Accepted in MNRAS

Gas flows in galaxy mergers: supersonic turbulence in bridges, accretion from the circumgalactic medium, and metallicity dilution

Authors: Martin Sparre, Joseph Whittingham, Mitali Damle, Maan H. Hani, Philipp Richter, Sara L. Ellison, Christoph Pfrommer, Mark Vogelsberger

Abstract: In major galaxy mergers, the orbits of stars are violently perturbed, and gas is torqued to the centre, diluting the gas metallicity and igniting a starburst. In this paper, we study the gas dynamics in and around merging galaxies using a series of cosmological magneto-hydrodynamical (MHD) zoom-in simulations. We find that the gas bridge connecting the merging galaxies pre-coalescence is dominated… ▽ More In major galaxy mergers, the orbits of stars are violently perturbed, and gas is torqued to the centre, diluting the gas metallicity and igniting a starburst. In this paper, we study the gas dynamics in and around merging galaxies using a series of cosmological magneto-hydrodynamical (MHD) zoom-in simulations. We find that the gas bridge connecting the merging galaxies pre-coalescence is dominated by turbulent pressure, with turbulent Mach numbers peaking at values of 1.6-3.3. This implies that bridges are dominated by supersonic turbulence, and are thus ideal candidates for studying the impact of extreme environments on star formation. We also find that gas accreted from the circumgalactic medium (CGM) during the merger significantly contributes (27-51 per cent) to the star formation rate (SFR) at the time of coalescence and drives the subsequent reignition of star formation in the merger remnant. Indeed, 19-53 per cent of the SFR at z=0 originates from gas belonging to the CGM prior the merger. Finally, we investigate the origin of the metallicity-diluted gas at the centre of merging galaxies. We show that this gas is rapidly accreted onto the galactic centre with a time-scale much shorter than that of normal star-forming galaxies. This explains why coalescing galaxies are not well-captured by the fundamental metallicity relation. △ Less

Submitted 7 October, 2021; originally announced October 2021.

Comments: 15 pages, 13 figures, submitted to MNRAS

The impact of magnetic fields on cosmological galaxy mergers. I: Reshaping gas and stellar discs

Authors: Joseph Whittingham, Martin Sparre, Christoph Pfrommer, Rüdiger Pakmor

Abstract: Mergers play an important role in galaxy evolution. In particular, major mergers are able to have a transformative effect on galaxy morphology. In this paper, we investigate the role of magnetic fields in gas-rich major mergers. To this end, we run a series of high-resolution magnetohydrodynamic (MHD) zoom-in simulations with the moving-mesh code Arepo and compare the outcome with hydrodynamic sim… ▽ More Mergers play an important role in galaxy evolution. In particular, major mergers are able to have a transformative effect on galaxy morphology. In this paper, we investigate the role of magnetic fields in gas-rich major mergers. To this end, we run a series of high-resolution magnetohydrodynamic (MHD) zoom-in simulations with the moving-mesh code Arepo and compare the outcome with hydrodynamic simulations run from the same initial conditions. This is the first time that the effect of magnetic fields in major mergers has been investigated in a cosmologically-consistent manner. In contrast to previous non-cosmological simulations, we find that the inclusion of magnetic fields has a substantial impact on the production of the merger remnant. Whilst magnetic fields do not strongly affect global properties, such as the star formation history, they are able to significantly influence structural properties. Indeed, MHD simulations consistently form remnants with extended discs and well-developed spiral structure, whilst hydrodynamic simulations form more compact remnants that display distinctive ring morphology. We support this work with a resolution study and show that whilst global properties are broadly converged across resolution and physics models, morphological differences only develop given sufficient resolution. We argue that this is due to the more efficient excitement of a small-scale dynamo in higher resolution simulations, resulting in a more strongly amplified field that is better able to influence gas dynamics. △ Less

Submitted 21 March, 2022; v1 submitted 27 November, 2020; originally announced November 2020.

Comments: 27 pages, 18 figures. Published in MNRAS on 19 May 2021

Highly ordered magnetic fields in the tail of the jellyfish galaxy JO206

Authors: Ancla Müller, Bianca Poggianti, Christoph Pfrommer, Björn Adebahr, Paolo Serra, Alessandro Ignesti, Martin Sparre, Myriam Gitti, Ralf-Jürgen Dettmar, Benedetta Vulcani, Alessia Moretti

Abstract: Jellyfish galaxies have long tails of gas that is stripped from the disc by ram pressure due to the motion of galaxies in the intracluster medium in galaxy clusters. We present the first measurement of the magnetic field strength and orientation within the disc and the (90$\,$kpc-long) $\rm Hα$-emitting tail of the jellyfish galaxy JO206. The tail has a large-scale magnetic field ($>4.1\,μ$G), a s… ▽ More Jellyfish galaxies have long tails of gas that is stripped from the disc by ram pressure due to the motion of galaxies in the intracluster medium in galaxy clusters. We present the first measurement of the magnetic field strength and orientation within the disc and the (90$\,$kpc-long) $\rm Hα$-emitting tail of the jellyfish galaxy JO206. The tail has a large-scale magnetic field ($>4.1\,μ$G), a steep radio spectral index ($α\sim -2.0$), indicating an aging of the electrons propagating away from the star-forming regions, and extremely high fractional polarisation ($>50\,$%), indicating low turbulent motions. The magnetic field vectors are aligned with (parallel to) the direction of the ionised-gas tail and stripping direction. High-resolution simulations of a large, cold gas cloud that is exposed to a hot, magnetised turbulent wind show that the high fractional polarisation and the ordered magnetic field can be explained by accretion of draped magnetised plasma from the hot wind that condenses onto the external layers of the tail, where it is adiabatically compressed and sheared. The ordered magnetic field, preventing heat and momentum exchange, may be a key factor in allowing in-situ star formation in the tail. △ Less

Submitted 26 October, 2020; v1 submitted 28 September, 2020; originally announced September 2020.

Comments: accepted for publication in Nature Astronomy, Sharedlt: https://rdcu.be/b9c0L (free read only)

Interaction of a cold cloud with a hot wind: the regimes of cloud growth and destruction and the impact of magnetic fields

Authors: Martin Sparre, Christoph Pfrommer, Kristian Ehlert

Abstract: Multiphase galaxy winds, the accretion of cold gas through galaxy haloes, and gas stripping from jellyfish galaxies are examples of interactions between cold and hot gaseous phases. There are two important regimes in such systems. A sufficiently small cold cloud is destroyed by the hot wind as a result of Kelvin-Helmholtz instabilities, which shatter the cloud into small pieces that eventually mix… ▽ More Multiphase galaxy winds, the accretion of cold gas through galaxy haloes, and gas stripping from jellyfish galaxies are examples of interactions between cold and hot gaseous phases. There are two important regimes in such systems. A sufficiently small cold cloud is destroyed by the hot wind as a result of Kelvin-Helmholtz instabilities, which shatter the cloud into small pieces that eventually mix and dissolve in the hot wind. On the contrary, stripped cold gas from a large cloud mixes with the hot wind to intermediate temperatures, and then becomes thermally unstable and cools, causing a net accretion of hot gas to the cold tail. Using the magneto-hydrodynamical code AREPO, we perform cloud crushing simulations and test analytical criteria for the transition between the growth and destruction regimes to clarify a current debate in the literature. We find that the hot-wind cooling time sets the transition radius and not the cooling time of the mixed phase. Magnetic fields modify the wind-cloud interaction. Draping of wind magnetic field enhances the field upstream of the cloud and fluid instabilities are suppressed by a turbulently magnetised wind beyond what is seen for a wind with a uniform magnetic field. We furthermore predict jellyfish galaxies to have ordered magnetic fields aligned with their tails. We finally discuss how the results of idealised simulations can be used to provide input to subgrid models in cosmological (magneto-)hydrodynamical simulations, which cannot resolve the detailed small-scale structure of cold gas clouds in the circum-galactic medium. △ Less

Submitted 12 October, 2020; v1 submitted 20 August, 2020; originally announced August 2020.

Comments: Accepted for publication in MNRAS

The time-scales probed by star formation rate indicators for realistic, bursty star formation histories from the FIRE simulations

Authors: José A. Flores Velázquez, Alexander B. Gurvich, Claude-André Faucher-Giguère, James S. Bullock, Tjitske K. Starkenburg, Jorge Moreno, Alexandres Lazar, Francisco J. Mercado, Jonathan Stern, Martin Sparre, Christopher C. Hayward, Andrew Wetzel, Kareem El-Badry

Abstract: Understanding the rate at which stars form is central to studies of galaxy formation. Observationally, the star formation rates (SFRs) of galaxies are measured using the luminosity in different frequency bands, often under the assumption of a time-steady SFR in the recent past. We use star formation histories (SFHs) extracted from cosmological simulations of star-forming galaxies from the FIRE pro… ▽ More Understanding the rate at which stars form is central to studies of galaxy formation. Observationally, the star formation rates (SFRs) of galaxies are measured using the luminosity in different frequency bands, often under the assumption of a time-steady SFR in the recent past. We use star formation histories (SFHs) extracted from cosmological simulations of star-forming galaxies from the FIRE project to analyze the time-scales to which the H$α$ and far-ultraviolet (FUV) continuum SFR indicators are sensitive. In these simulations, the SFRs are highly time variable for all galaxies at high redshift, and continue to be bursty to z=0 in dwarf galaxies. When FIRE SFHs are partitioned into their bursty and time-steady phases, the best-fitting FUV time-scale fluctuates from its ~10 Myr value when the SFR is time-steady to >~100 Myr immediately following particularly extreme bursts of star formation during the bursty phase. On the other hand, the best-fitting averaging time-scale for H$α$ is generally insensitive to the SFR variability in the FIRE simulations and remains ~5 Myr at all times. These time-scales are shorter than the 100 Myr and 10 Myr time-scales sometimes assumed in the literature for FUV and H$α$, respectively, because while the FUV emission persists for stellar populations older than 100 Myr, the time-dependent luminosities are strongly dominated by younger stars. Our results confirm that the ratio of SFRs inferred using H$α$ vs. FUV can be used to probe the burstiness of star formation in galaxies. △ Less

Submitted 15 December, 2020; v1 submitted 19 August, 2020; originally announced August 2020.

Comments: 14 pages, 10 figures, accepted to MNRAS

The Hestia project: simulations of the Local Group

Authors: Noam I. Libeskind, Edoardo Carlesi, Rob J. J. Grand, Arman Khalatyan, Alexander Knebe, Ruediger Pakmor, Sergey Pilipenko, Marcel S. Pawlowski, Martin Sparre, Elmo Tempel, Peng Wang, Helene M. Courtois, Stefan Gottloeber, Yehuda Hoffman, Ivan Minchev, Christoph Pfrommer, Jenny G. Sorce, Volker Springel, Matthias Steinmetz, R. Brent Tully, Mark Vogelsberger, Gustavo Yepes

Abstract: We present the Hestia simulation suite: High-resolutions Environmental Simulations of The Immediate Area, a set of cosmological simulations of the Local Group. Initial conditions constrained by the observed peculiar velocity of nearby galaxies are employed to accurately simulate the local cosmography. Halo pairs that resemble the Local Group are found in low resolutions constrained, dark matter on… ▽ More We present the Hestia simulation suite: High-resolutions Environmental Simulations of The Immediate Area, a set of cosmological simulations of the Local Group. Initial conditions constrained by the observed peculiar velocity of nearby galaxies are employed to accurately simulate the local cosmography. Halo pairs that resemble the Local Group are found in low resolutions constrained, dark matter only simulations, and selected for higher resolution magneto hydrodynamic simulation using the Arepo code. Baryonic physics follows the Auriga model of galaxy formation. The simulations contain a high resolution region of 3-5 Mpc in radius from the Local Group midpoint embedded in the correct cosmographic landscape. Within this region a simulated Local Group consisting of a Milky Way and Andromeda like galaxy forms, whose description is in excellent agreement with observations. The simulated Local Group galaxies resemble the Milky Way and Andromeda in terms of their halo mass, mass ratio, stellar disc mass, morphology separation, relative velocity, rotation curves, bulge-disc morphology, satellite galaxy stellar mass function, satellite radial distribution and in some cases, the presence of a Magellanic cloud like object. Because these simulations properly model the Local Group in their cosmographic context, they provide a testing ground for questions where environment is thought to play an important role. △ Less

Submitted 11 August, 2020; originally announced August 2020.

Comments: 20 pages, 13 figures, submitted to MNRAS

Submillimetre galaxies in cosmological hydrodynamical simulations -- an opportunity for constraining feedback models

Authors: Christopher C. Hayward, Martin Sparre, Scott C. Chapman, Lars Hernquist, Dylan Nelson, Rüdiger Pakmor, Annalisa Pillepich, Volker Springel, Paul Torrey, Mark Vogelsberger, Rainer Weinberger

Abstract: 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. We use a scaling relation derived from previous simulations plus radiative transfer calculations to predict the submm flux densities of simulated SMGs drawn from cosmologic… ▽ More 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. We use a scaling relation derived from previous simulations plus radiative transfer calculations to predict the submm flux densities of simulated SMGs drawn from cosmological simulations from the Illustris and IllustrisTNG projects based on the simulated galaxies' star formation rates (SFRs) and dust masses 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. The redshift distribution of SMGs with $S_{850} > 3$ mJy in IllustrisTNG is consistent with the observed distribution, whereas the Illustris redshift distribution peaks at significantly lower redshift (1.5 vs. 2.8). We demonstrate that IllustrisTNG hosts fewer SMGs than Illustris because in the former, high-mass ($M_{\star} \sim 10^{11} \, \text{M}_{\odot}$) $z \sim 2-3$ galaxies have lower dust masses and 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. △ Less

Submitted 10 February, 2021; v1 submitted 3 July, 2020; originally announced July 2020.

Comments: MNRAS, in press; typo corrected

Galaxy pairs in the Sloan Digital Sky Survey -- XIV. Galaxy mergers do not lie on the Fundamental Metallicity Relation

Authors: Sebastián Bustamante, Sara L. Ellison, David R. Patton, Martin Sparre

Abstract: 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 qualit… ▽ More 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$. △ Less

Submitted 13 April, 2020; originally announced April 2020.

Comments: 13 pages, 10 figures. Accepted for publication in MNRAS

Quiescent galaxies 1.5 billion years after the Big Bang and their progenitors

Authors: Francesco Valentino, Masayuki Tanaka, Iary Davidzon, Sune Toft, Carlos Gomez-Guijarro, Mikkel Stockmann, Masato Onodera, Gabriel Brammer, Daniel Ceverino, Andreas L. Faisst, Anna Gallazzi, Christopher C. Hayward, Olivier Ilbert, Mariko Kubo, Georgios E. Magdis, Jonatan Selsing, Rhythm Shimakawa, Martin Sparre, Charles Steinhardt, Kiyoto Yabe, Johannes Zabl

Abstract: 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 alter… ▽ More 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. △ Less

Submitted 20 December, 2019; v1 submitted 23 September, 2019; originally announced September 2019.

Comments: 20 pages (+appendix), 11 figures, accepted for publication in the Astrophysical Journal (ApJ). This work comes with a companion letter by Tanaka, Valentino, Toft et al. 2019, ApJL, 885, L34

The diversity of the circumgalactic medium around z = 0 Milky Way-mass galaxies from the Auriga simulations

Authors: Maan H. Hani, Sara L. Ellison, Martin Sparre, Robert J. J. Grand, Rüediger Pakmor, Facundo A. Gomez, Volker Springel

Abstract: 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 f… ▽ More 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. △ Less

Submitted 9 July, 2019; originally announced July 2019.

Journal ref: Monthly Notices of the Royal Astronomical Society 488 (2019) 135-152

The interplay of Self-Interacting Dark Matter and baryons in shaping the halo evolution

Authors: Giulia Despali, Martin Sparre, Simona Vegetti, Mark Vogelsberger, Jesús Zavala, Federico Marinacci

Abstract: 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 $σ^\text{T}/m_χ=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… ▽ More 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 $σ^\text{T}/m_χ=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. △ Less

Submitted 28 January, 2019; v1 submitted 6 November, 2018; originally announced November 2018.

Comments: 11 pages, 11 figures, accepted for publication in MNRAS

The physics of multiphase gas flows: fragmentation of a radiatively cooling gas cloud in a hot wind

Authors: Martin Sparre, Christoph Pfrommer, Mark Vogelsberger

Abstract: 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… ▽ More 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-$α$ haloes, high-velocity clouds, and broad-line regions of active galactic nuclei. △ Less

Submitted 7 November, 2018; v1 submitted 20 July, 2018; originally announced July 2018.

Comments: Accepted for publication in MNRAS

Highly-ionized metals as probes of the circumburst gas in the natal regions of gamma-ray bursts

Authors: K. E. Heintz, D. Watson, P. Jakobsson, J. P. U. Fynbo, J. Bolmer, M. Arabsalmani, Z. Cano, S. Covino, V. D'Elia, A. Gomboc, J. Japelj, L. Kaper, J. -K. Krogager, G. Pugliese, R. Sánchez-Ramírez, J. Selsing, M. Sparre, N. R. Tanvir, C. C. Thöne, A. de Ugarte Postigo, S. D. Vergani

Abstract: 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… ▽ More 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. △ Less

Submitted 4 June, 2018; originally announced June 2018.

Comments: 13 pages, 8 figures + Appendix. Accepted in MNRAS

The fraction of ionizing radiation from massive stars that escapes to the intergalactic medium

Authors: N. R. Tanvir, J. P. U. Fynbo, A. de Ugarte Postigo, J. Japelj, K. Wiersema, D. Malesani, D. A. Perley, A. J. Levan, J. Selsing, S. B. Cenko, D. A. Kann, B. Milvang-Jensen, E. Berger, Z. Cano, R. Chornock, S. Covino, A. Cucchiara, V. D'Elia, P. Goldoni, A. Gomboc, K. E. Heintz, J. Hjorth, L. Izzo, P. Jakobsson, L. Kaper , et al. (16 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 18 May, 2018; originally announced May 2018.

Comments: 31 pages

The X-shooter GRB afterglow legacy sample (XS-GRB)

Authors: J. Selsing, D. Malesani, P. Goldoni, J. P. U. Fynbo, T. Krühler, L. A. Antonelli, M. Arabsalmani, J. Bolmer, Z. Cano, L. Christensen, S. Covino, P. D'Avanzo, V. D'Elia, A. De Cia, A. de Ugarte Postigo, H. Flores, M. Friis, A. Gomboc, J. Greiner, P. Groot, F. Hammer, O. E. Hartoog, K. E. Heintz, J. Hjorth, P. Jakobsson , et al. (31 additional authors not shown)

Abstract: In this work we present spectra of all $γ$-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 observation… ▽ More In this work we present spectra of all $γ$-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. △ Less

Submitted 21 February, 2018; originally announced February 2018.

Comments: 41 pages, 10 Figures, 4 Tables. Submitted to A&A. Paper and code also available at https://github.com/jselsing/XSGRB-sample-paper

Journal ref: A&A 623, A92 (2019)

Galaxy mergers moulding the circum-galactic medium - I. The impact of a major merger

Authors: Maan H. Hani, Martin Sparre, Sara L. Ellison, Paul Torrey, Mark Vogelsberger

Abstract: 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.… ▽ More 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. △ Less

Submitted 18 January, 2018; originally announced January 2018.

Journal ref: Monthly Notices of the Royal Astronomical Society 475 (2018) 1160 - 1176

An alternate approach to measure specific star formation rates at 2<z<7

Authors: Iary Davidzon, Olivier Ilbert, Andreas L. Faisst, Martin Sparre, Peter L. Capak

Abstract: 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 se… ▽ More 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. △ Less

Submitted 11 December, 2017; originally announced December 2017.

Comments: ApJ accepted

Merger-Induced Metallicity Dilution in Cosmological Galaxy Formation Simulations

Authors: Sebastian Bustamante, Martin Sparre, Volker Springel, Robert J. J. Grand

Abstract: 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 stat… ▽ More 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. △ Less

Submitted 16 July, 2018; v1 submitted 8 December, 2017; originally announced December 2017.

Comments: Submitted to MNRAS

Mass and metallicity scaling relations of high redshift star-forming galaxies selected by GRBs

Authors: M. Arabsalmani, P. Møller, D. A. Perley, W. Freudling, J. P. U. Fynbo, E. Le Floc'h, M. A. Zwaan, S. Schulze, N. R. Tanvir, L. Christensen, A. J. Levan, P. Jakobsson, D. Malesani, Z. Cano, S. Covino, V. D'Elia, P. Goldoni, A. Gomboc, K. E. Heintz, M. Sparre, A. de Ugarte Postigo, S. D. Vergani

Abstract: 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 der… ▽ More 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. △ Less

Submitted 19 September, 2017; originally announced September 2017.

Comments: Accepted for publication in MNRAS

Origin of chemically distinct discs in the Auriga cosmological simulations

Authors: Robert J. J. Grand, Sebastián Bustamante, Facundo A. Gómez, Daisuke Kawata, Federico Marinacci, Rüdiger Pakmor, Hans-Walter Rix, Christine M. Simpson, Martin Sparre, Volker Springel

Abstract: 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 cle… ▽ More 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 α/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[α/Fe]$ star formation phase in the inner region ($R\lesssim 5$ kpc) induced by gas-rich mergers, followed by more quiescent low-$\rm[α/Fe]$ star formation; and b) an early phase of high-$\rm[α/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[α/Fe]$ sequence stars. We discuss the implication that both processes occurred in the Milky Way. △ Less

Submitted 20 November, 2017; v1 submitted 25 August, 2017; originally announced August 2017.

Comments: 12 pages, MNRAS accepted

The host galaxy of the short GRB111117A at $z = 2.211$: impact on the short GRB redshift distribution and progenitor channels

Authors: J. Selsing, T. Krühler, D. Malesani, P. D'Avanzo, S. Schulze, J. Palmerio, S. D. Vergani, J. Japelj, B. Milvang-Jensen, D. Watson, P. Jakobsson, J. Bolmer, Z. Cano, S. Covino, L. B. Christensen, V. D'Elia, A. de Ugarte Postigo, J. P. U. Fynbo, A. Gomboc, K. E. Heintz, L. Kaper, A. J. Levan, S. Piranomonte, G. Pugliese, R. Sánchez-Ramírez , et al. (4 additional authors not shown)

Abstract: It is notoriously difficult to localize short $γ$-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$. T… ▽ More It is notoriously difficult to localize short $γ$-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. △ Less

Submitted 3 July, 2017; originally announced July 2017.

Comments: 8 pages, 5 figures, 2 tables. Submitted for publication in A&A. All data, code and calculations related to the paper along with the paper itself are available at https://github.com/jselsing/GRB111117A

Journal ref: A&A 616, A48 (2018)

Log-normal star formation histories in simulated and observed galaxies

Authors: Benedikt Diemer, Martin Sparre, Louis E. Abramson, Paul Torrey

Abstract: 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 feature… ▽ More 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. △ Less

Submitted 10 April, 2017; v1 submitted 9 January, 2017; originally announced January 2017.

Comments: 20 pages, 18 figures. Minor changes to match the version published in ApJ

Journal ref: 2017 ApJ 839, 26

Cosmic evolution and metal aversion in super-luminous supernova host galaxies

Authors: S. Schulze, T. Krühler, G. Leloudas, J. Gorosabel, A. Mehner, J. Buchner, S. Kim, E. Ibar, R. Amorín, R. Herrero-Illana, J. P. Anderson, F. E. Bauer, L. Christensen, M. de Pasquale, A. de Ugarte Postigo, A. Gallazzi, J. Hjorth, N. Morrell, D. Malesani, M. Sparre, B. Stalder, A. A. Stark, C. C. Thöne, J. C. Wheeler

Abstract: 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 com… ▽ More 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σ$ 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. △ Less

Submitted 3 November, 2017; v1 submitted 18 December, 2016; originally announced December 2016.

Comments: Accepted on 8 September 2017. Currently in press at MNRAS. Substantially revised version. 61 pages, 31 pages main text, 30 pages appendices, 22 figures and 9 tables

The unorthodox evolution of major merger remnants into star-forming spiral galaxies

Authors: Martin Sparre, Volker Springel

Abstract: 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.… ▽ More 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. △ Less

Submitted 15 June, 2017; v1 submitted 12 October, 2016; originally announced October 2016.

Comments: 15 pages, 9 figures, Accepted for publication in MNRAS

GRB 110715A: The peculiar multiwavelength evolution of the first afterglow detected by ALMA

Authors: R. Sánchez-Ramírez, P. J. Hancock, G. Jóhannesson, Tara Murphy, A. de Ugarte Postigo, J. Gorosabel, D. A. Kann, T. Krühler, S. R. Oates, J. Japelj, C. C. Thöne, A. Lundgren, D. A. Perley, D. Malesani, I. de Gregorio Monsalvo, A. J. Castro-Tirado, V. D'Elia, J. P. U. Fynbo, D. Garcia-Appadoo, P. Goldoni, J. Greiner, Y. -D. Hu, M. Jelínek, S. Jeong, A. Kamble , et al. (13 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 11 October, 2016; v1 submitted 6 October, 2016; originally announced October 2016.

Comments: Accepted for publication in MNRAS

Journal ref: MNRAS, 464, 4624 (2017)

Zooming in on major mergers: dense, starbursting gas in cosmological simulations

Authors: Martin Sparre, Volker Springel

Abstract: We introduce the `Illustris zoom simulation project', which allows the study of selected galaxies forming in the $Λ$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 bur… ▽ More We introduce the `Illustris zoom simulation project', which allows the study of selected galaxies forming in the $Λ$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. △ Less

Submitted 20 July, 2016; v1 submitted 27 April, 2016; originally announced April 2016.

Comments: 13 pages, 7 figures, Accepted for publication in MNRAS

The optical identifcation of events with poorly defined locations: The case of the Fermi GBM GRB140801A

Authors: V. M. Lipunov, J. Gorosabel, M. V. Pruzhinskaya, A. de Ugarte Postigo, V. Pelassa, A. E. Tsvetkova, I. V. Sokolov, D. A. Kann, Dong Xu, E. S. Gorbovskoy, V. V. Krushinski, V. G. Kornilov, P. V. Balanutsa, S. V. Boronina, N. M. Budnev, Z. Cano, A. J. Castro-Tirado, V. V. Chazov, V. Connaughton, C. Delvaux, D. D. Frederiks, J. F. U. Fynbo, A. V. Gabovich, A. Goldstein, J. Greiner , et al. (28 additional authors not shown)

Abstract: We report the early discovery of the optical afterglow of gamma-ray burst (GRB) 140801A in the 137 deg$^2$ 3-$σ$ 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… ▽ More We report the early discovery of the optical afterglow of gamma-ray burst (GRB) 140801A in the 137 deg$^2$ 3-$σ$ 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 $θ=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. △ Less

Submitted 27 October, 2015; originally announced October 2015.

Comments: in press MNRAS, 2015

Journal ref: MNRAS, 455, 712 (2016)

(Star)bursts of FIRE: observational signatures of bursty star formation in galaxies

Authors: Martin Sparre, Christopher C. Hayward, Robert Feldmann, Claude-André Faucher-Giguère, Alexander L. Muratov, Dušan Kereš, Philip F. Hopkins

Abstract: 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 a… ▽ More 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$α$-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$α$- and FUV-derived SFR--$M_*$ relations at $z=2$ is consistent with observational constraints. We also find that the H$α$/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$α$/FUV ratios than observed. We suggest that future cosmological simulations should compare the H$α$/FUV ratios of their galaxies with observations to constrain the feedback models employed. △ Less

Submitted 18 November, 2016; v1 submitted 13 October, 2015; originally announced October 2015.

Comments: 18 pages, 12 figures. Accepted for publication in MNRAS. Figure 3,4,10 and A1 are new. Figure 5 now includes an estimation of stochastic IMF sampling

Journal ref: MNRAS 466 (2017) 88-104

The mysterious optical afterglow spectrum of GRB140506A at z=0.889

Authors: J. P. U. Fynbo, T. Krühler, K. Leighly, C. Ledoux, P. M. Vreeswijk, S. Schulze, P. Noterdaeme, D. Watson, R. A. M. J. Wijers, J. Bolmer, Z. Cano, L. Christensen, S. Covino, V. D'Elia, H. Flores, M. Friis, P. Goldoni, J. Greiner, F. Hammer, J. Hjorth, P. Jakobsson, J. Japelj, L. Kaper, S. Klose, F. Knust , et al. (18 additional authors not shown)

Abstract: 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-shoo… ▽ More 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. △ Less

Submitted 18 September, 2014; v1 submitted 17 September, 2014; originally announced September 2014.

Comments: 8 pages, 4 figures. Accepted for publications in A&A

Journal ref: A&A 572, A12 (2014)

VLT/X-shooter spectroscopy of the afterglow of the Swift GRB 130606A: Chemical abundances and reionisation at $z\sim6$

Authors: O. E. Hartoog, D. Malesani, J. P. U. Fynbo, T. Goto, T. Krühler, P. M. Vreeswijk, A. De Cia, D. Xu, P. Møller, S. Covino, V. D'Elia, H. Flores, P. Goldoni, J. Hjorth, P. Jakobsson, J. -K. Krogager, L. Kaper, C. Ledoux, A. J. Levan, B. Milvang-Jensen, J. Sollerman, M. Sparre, G. Tagliaferri, N. R. Tanvir, A. de Ugarte Postigo , et al. (6 additional authors not shown)

Abstract: 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… ▽ More 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$α$, $β$ and $γ$ 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$α$-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$σ$). 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.) △ Less

Submitted 9 July, 2015; v1 submitted 16 September, 2014; originally announced September 2014.

Comments: 15 pages, 12 figures

Journal ref: A&A 580, A139 (2015)

The star formation main sequence and stellar mass assembly of galaxies in the Illustris simulation

Authors: Martin Sparre, Christopher C. Hayward, Volker Springel, Mark Vogelsberger, Shy Genel, Paul Torrey, Dylan Nelson, Debora Sijacki, Lars Hernquist

Abstract: 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… ▽ More 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. △ Less

Submitted 18 December, 2014; v1 submitted 29 August, 2014; originally announced September 2014.

Comments: 16 pages, 16 figures, Accepted for MNRAS, Section 6 and Figure 16 are new. Data is available at http://martinsparre.com/dataproducts

Journal ref: MNRAS 2015, 447, 3548

Particle ejection during mergers of dark matter halos

Authors: Isabella P. Carucci, Martin Sparre, Steen H. Hansen, Michael Joyce

Abstract: 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 signifi… ▽ More 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. △ Less

Submitted 26 May, 2014; originally announced May 2014.

Comments: 18 pages, 12 figures. Accepted for publication in JCAP

Journal ref: JCAP06(2014)057

GRB 120422A/SN 2012bz: Bridging the Gap between Low- And High-Luminosity GRBs

Authors: S. Schulze, D. Malesani, A. Cucchiara, N. R. Tanvir, T. Krühler, A. de Ugarte Postigo, G. Leloudas, J. Lyman, D. Bersier, K. Wiersema, D. A. Perley, P. Schady, J. Gorosabel, J. P. Anderson, A. J. Castro-Tirado, S. B. Cenko, A. De Cia, L. E. Ellerbroek, J. P. U. Fynbo, J. Greiner, J. Hjorth, D. A. Kann, L. Kaper, S. Klose, A. J. Levan , et al. (40 additional authors not shown)

Abstract: 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 ultrarelati… ▽ More 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 $γ$-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$α$ 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 $Γ_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] △ Less

Submitted 15 January, 2014; originally announced January 2014.

Comments: 30 pages, 17 figures, 9 tables; abstract is abridged; images are shown at reduced resolution; comments are welcome

The Metallicity of High-$z$ GRBs: the Case of GRB 111008A

Authors: Martin Sparre

Abstract: 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-$α$ absorbers have a rapid decrease in metallicity at this redshift (Rafelski et al. 2013). In th… ▽ More 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-$α$ 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. △ Less

Submitted 12 December, 2013; originally announced December 2013.

Comments: 3 pages. Proceeding for the workshop "Metal Production and Distribution in a Hierarchical Universe", Paris, 21 - 25 Oct. 2013

The metallicity and dust content of a redshift 5 gamma-ray burst host galaxy

Authors: M. Sparre, O. E. Hartoog, T. Krühler, J. P. U. Fynbo, D. J. Watson, K. Wiersema, V. D'Elia, T. Zafar, P. M. J. Afonso, S. Covino, A. de Ugarte Postigo, H. Flores, P. Goldoni, J. Greiner, J. Hjorth, P. Jakobsson, L. Kaper, S. Klose, A. J. Levan, D. Malesani, B. Milvang-Jensen, M. Nardini, S. Piranomonte, J. Sollerman, R. Sánchez-Ramírez , et al. (4 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 12 March, 2014; v1 submitted 11 September, 2013; originally announced September 2013.

Comments: Accepted for publication in ApJ

Journal ref: 2014 ApJ 785 150

Spectroscopy of the short-hard GRB 130603B: The host galaxy and environment of a compact object merger

Authors: A. de Ugarte Postigo, C. C. Thoene, A. Rowlinson, R. Garcia-Benito, A. J. Levan, J. Gorosabel, P. Goldoni, S. Schulze, T. Zafar, K. Wiersema, R. Sanchez-Ramirez, A. Melandri, P. D'Avanzo, S. Oates, V. D'Elia, M. De Pasquale, T. Kruehler, A. J. van der Horst, D. Xu, D. Watson, S. Piranomonte, S. Vergani, B. Milvang-Jensen, L. Kaper, D. Malesani , et al. (18 additional authors not shown)

Abstract: 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, u… ▽ More 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. △ Less

Submitted 23 December, 2013; v1 submitted 13 August, 2013; originally announced August 2013.

Comments: Accepted for publication in Astronomy and Astrophysics

Journal ref: Astron. Astrophys. 563 (2014) A62

Discovery of the broad-lined Type Ic SN 2013cq associated with the very energetic GRB 130427A

Authors: D. Xu, A. de Ugarte Postigo, G. Leloudas, T. Kruhler, Z. Cano, J. Hjorth, D. Malesani, J. P. U. Fynbo, C. C. Thoene, R. Sanchez-Ramirez, S. Schulze, P. Jakobsson, L. Kaper, J. Sollerman, D. J. Watson, A. Cabrera-Lavers, C. Cao, S. Covino, H. Flores, S. Geier, J. Gorosabel, S. M. Hu, B. Milvang-Jensen, M. Sparre, L. P. Xin , et al. (3 additional authors not shown)

Abstract: 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 reds… ▽ More 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. △ Less

Submitted 20 August, 2013; v1 submitted 29 May, 2013; originally announced May 2013.

Comments: 12 pages, 4 figures, 1 table, accepted for publication in ApJ

Journal ref: ApJ 776 (2013) 98