Showing 1–50 of 101 results for author: Sijacki, D

Weak lensing combined with the kinetic Sunyaev Zel'dovich effect: A study of baryonic feedback

Authors: L. Bigwood, A. Amon, A. Schneider, J. Salcido, I. G. McCarthy, C. Preston, D. Sanchez, D. Sijacki, E. Schaan, S. Ferraro, N. Battaglia, A. Chen, S. Dodelson, A. Roodman, A. Pieres, A. Ferte, A. Alarcon, A. Drlica-Wagner, A. Choi, A. Navarro-Alsina, A. Campos, A. J. Ross, A. Carnero Rosell, B. Yin, B. Yanny , et al. (100 additional authors not shown)

Abstract: Extracting precise cosmology from weak lensing surveys requires modelling the non-linear matter power spectrum, which is suppressed at small scales due to baryonic feedback processes. However, hydrodynamical galaxy formation simulations make widely varying predictions for the amplitude and extent of this effect. We use measurements of Dark Energy Survey Year 3 weak lensing (WL) and Atacama Cosmolo… ▽ More Extracting precise cosmology from weak lensing surveys requires modelling the non-linear matter power spectrum, which is suppressed at small scales due to baryonic feedback processes. However, hydrodynamical galaxy formation simulations make widely varying predictions for the amplitude and extent of this effect. We use measurements of Dark Energy Survey Year 3 weak lensing (WL) and Atacama Cosmology Telescope DR5 kinematic Sunyaev-Zel'dovich (kSZ) to jointly constrain cosmological and astrophysical baryonic feedback parameters using a flexible analytical model, `baryonification'. First, using WL only, we compare the $S_8$ constraints using baryonification to a simulation-calibrated halo model, a simulation-based emulator model and the approach of discarding WL measurements on small angular scales. We find that model flexibility can shift the value of $S_8$ and degrade the uncertainty. The kSZ provides additional constraints on the astrophysical parameters and shifts $S_8$ to $S_8=0.823^{+0.019}_{-0.020}$, a higher value than attained using the WL-only analysis. We measure the suppression of the non-linear matter power spectrum using WL + kSZ and constrain a mean feedback scenario that is more extreme than the predictions from most hydrodynamical simulations. We constrain the baryon fractions and the gas mass fractions and find them to be generally lower than inferred from X-ray observations and simulation predictions. We conclude that the WL + kSZ measurements provide a new and complementary benchmark for building a coherent picture of the impact of gas around galaxies across observations. △ Less

Submitted 9 April, 2024; originally announced April 2024.

A dormant, overmassive black hole in the early Universe

Authors: Ignas Juodžbalis, Roberto Maiolino, William M. Baker, Sandro Tacchella, Jan Scholtz, Francesco D'Eugenio, Raffaella Schneider, Alessandro Trinca, Rosa Valiante, Christa DeCoursey, Mirko Curti, Stefano Carniani, Jacopo Chevallard, Anna de Graaff, Santiago Arribas, Jake S. Bennett, Martin A. Bourne, Andrew J. Bunker, Stéphane Charlot, Brian Jiang, Sophie Koudmani, Michele Perna, Brant Robertson, Debora Sijacki, Hannah Übler , et al. (3 additional authors not shown)

Abstract: Recent observations have found a large number of supermassive black holes already in place in the first few hundred million years after Big Bang. The channels of formation and growth of these early, massive black holes are not clear, with scenarios ranging from heavy seeds to light seeds experiencing bursts of high accretion rate. Here we present the detection, from the JADES survey, of broad Halp… ▽ More Recent observations have found a large number of supermassive black holes already in place in the first few hundred million years after Big Bang. The channels of formation and growth of these early, massive black holes are not clear, with scenarios ranging from heavy seeds to light seeds experiencing bursts of high accretion rate. Here we present the detection, from the JADES survey, of broad Halpha emission in a galaxy at z=6.68, which traces a black hole with mass of ~ 4 * 10^8 Msun and accreting at a rate of only 0.02 times the Eddington limit. The host galaxy has low star formation rate (~ 1 Msun/yr, a factor of 3 below the star forming main sequence). The black hole to stellar mass ratio is ~ 0.4, i.e. about 1,000 times above the local relation, while the system is closer to the local relations in terms of dynamical mass and velocity dispersion of the host galaxy. This object is most likely the tip of the iceberg of a much larger population of dormant black holes around the epoch of reionisation. Its properties are consistent with scenarios in which short bursts of super-Eddington accretion have resulted in black hole overgrowth and massive gas expulsion from the accretion disk; in between bursts, black holes spend most of their life in a dormant state. △ Less

Submitted 6 March, 2024; originally announced March 2024.

Comments: 38 pages, 18 figures. Submitted

Ly$α$ emission as a sensitive probe of feedback-regulated LyC escape from dwarf galaxies

Authors: Yuxuan Yuan, Sergio Martin-Alvarez, Martin G. Haehnelt, Thibault Garel, Debora Sijacki

Abstract: Ly$α$ emission is an exceptionally informative tracer of the life cycle of evolving galaxies and the escape of ionising photons. However, theoretical studies of Ly$α$ emission are often limited by insufficient numerical resolution, incomplete sets of physical models, and poor line-of-sight (LOS) statistics. To overcome such limitations, we utilize here the novel PANDORA suite of high-resolution dw… ▽ More Ly$α$ emission is an exceptionally informative tracer of the life cycle of evolving galaxies and the escape of ionising photons. However, theoretical studies of Ly$α$ emission are often limited by insufficient numerical resolution, incomplete sets of physical models, and poor line-of-sight (LOS) statistics. To overcome such limitations, we utilize here the novel PANDORA suite of high-resolution dwarf galaxy simulations that include a comprehensive set of state-of-the-art physical models for ionizing radiation, magnetic fields, supernova feedback and cosmic rays. We post-process the simulations with the radiative transfer code \textsc{RASCAS} to generate synthetic observations and compare to observed properties of Ly$α$ emitters. Our simulated Ly$α$ haloes are more extended than the spatial region from which the intrinsic emission emanates and our spatially resolved maps of spectral parameters of the Ly$α$ emission are very sensitive to the underlying spatial distribution and kinematics of neutral hydrogen. Ly$α$ and LyC emission display strongly varying signatures along different LOS depending on how each LOS intersects low-density channels generated by stellar feedback. Comparing galaxies simulated with different physics, we find the Ly$α$ signatures to exhibit systematic offsets determined by the different levels of feedback strength and the clumpiness of the neutral gas. Despite this variance, and regardless of the different physics included in each model, we find universal correlations between Ly$α$ observables and LyC escape fraction, demonstrating a robust connection between Ly$α$ and LyC emission. Ly$α$ observations from a large sample of dwarf galaxies should thus give strong constraints on their stellar feedback-regulated LyC escape and confirm their important role for the reionization of the Universe. △ Less

Submitted 26 June, 2024; v1 submitted 4 January, 2024; originally announced January 2024.

A unified accretion disc model for supermassive black holes in galaxy formation simulations: method and implementation

Authors: Sophie Koudmani, Rachel S. Somerville, Debora Sijacki, Martin A. Bourne, Yan-Fei Jiang, Kasar Profit

Abstract: It is well established that supermassive black hole (SMBH) feedback is crucial for regulating the evolution of massive, if not all, galaxies. However, modelling the interplay between SMBHs and their host galaxies is challenging due to the vast dynamic range. Previous simulations have utilized simple subgrid models for SMBH accretion, while recent advancements track the properties of the unresolved… ▽ More It is well established that supermassive black hole (SMBH) feedback is crucial for regulating the evolution of massive, if not all, galaxies. However, modelling the interplay between SMBHs and their host galaxies is challenging due to the vast dynamic range. Previous simulations have utilized simple subgrid models for SMBH accretion, while recent advancements track the properties of the unresolved accretion disc, usually based on the thin $α$-disc model. However, this neglects accretion in the radiatively inefficient regime, expected to occur through a thick disc for a significant portion of an SMBH's lifetime. To address this, we present a novel 'unified' accretion disc model for SMBHs, harnessing results from the analytical advection-dominated inflow-outflow solution (ADIOS) model and state-of-the-art GR(R)MHD simulations. Going from low to high Eddington ratios, our model transitions from an ADIOS flow to a thin $α$-disc via a truncated disc, incorporating self-consistently SMBH spin evolution due to Lense-Thirring precession. Utilizing the moving mesh code AREPO, we perform simulations of single and binary SMBHs within gaseous discs to validate our model and assess its impact. The disc state significantly affects observable luminosities, and we predict markedly different electromagnetic counterparts in SMBH binaries. Crucially, the assumed disc model shapes SMBH spin magnitudes and orientations, parameters that gravitational wave observatories like LISA and IPTA are poised to constrain. Our simulations emphasize the importance of accurately modelling SMBH accretion discs and spin evolution, as they modulate the available accretion power, profoundly shaping the interaction between SMBHs and their host galaxies. △ Less

Submitted 3 July, 2024; v1 submitted 13 December, 2023; originally announced December 2023.

Comments: Accepted for publication in MNRAS; 29 pages, 17 figures, 1 appendix; doi:10.1093/mnras/stae1422

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 532, Issue 1, July 2024, Pages 60-88

Dynamics and spin alignment in massive, gravito-turbulent circumbinary discs around supermassive black hole binaries

Authors: Martin A. Bourne, Davide Fiacconi, Debora Sijacki, Joanna M. Piotrowska, Sophie Koudmani

Abstract: Parsec-scale separation supermassive black hole binaries in the centre of gas-rich galaxy merger remnants could be surrounded by massive circumbinary discs (CBDs). Black hole mass and spin evolution during the gas-rich binary inspiral are crucial in determining the direction and power of relativistic jets that radio observations with LOFAR and SKAO will probe, and for predicting gravitational wave… ▽ More Parsec-scale separation supermassive black hole binaries in the centre of gas-rich galaxy merger remnants could be surrounded by massive circumbinary discs (CBDs). Black hole mass and spin evolution during the gas-rich binary inspiral are crucial in determining the direction and power of relativistic jets that radio observations with LOFAR and SKAO will probe, and for predicting gravitational wave (GW) emission that IPTA and LISA will measure. We present 3D hydrodynamic simulations capturing gas-rich, self-gravitating CBDs around a $2\times 10^6$M$_{\odot}$ supermassive black hole binary, that probe different mass ratios, eccentricities and inclinations. We employ a sub-grid Shakura-Sunyaev accretion disc to self-consistently model black hole mass and spin evolution together with super-Lagrangian refinement techniques to resolve gas flows, streams and mini-discs within the cavity, which play a fundamental role in torquing and feeding the binary. We find that higher mass ratio and eccentric binaries result in larger cavities, while retrograde binaries result in smaller cavities. All of the simulated binaries are expected to shrink with net gravitational torques being negative. Unlike previous simulations, we do not find preferential accretion onto the secondary black hole. This implies smaller chirp masses at coalescence and hence a weaker GW background. Critically this means that spin-alignment is faster than the binary inspiral timescale even for low mass ratios. However, we find that mini-disc and hence spin alignment is not guaranteed in initially misaligned systems, potentially leading to a significant fraction of recoiled remnants displaced from their host galaxies if chaotic accretion is the dominant feeding channel. △ Less

Submitted 28 November, 2023; originally announced November 2023.

Comments: 27 pages, 14 figures, submitted to MNRAS. Comments and feedback welcome

The impact of cosmic rays on the interstellar medium and galactic outflows of Milky Way analogues

Authors: Francisco Rodríguez Montero, Sergio Martin-Alvarez, Adrianne Slyz, Julien Devriendt, Yohan Dubois, Debora Sijacki

Abstract: During the last decade, cosmological simulations have managed to reproduce realistic and morphologically diverse galaxies, spanning the Hubble sequence. Central to this success was a phenomenological calibration of the few included feedback processes, whilst glossing over higher complexity baryonic physics. This approach diminishes the predictive power of such simulations, preventing to further ou… ▽ More During the last decade, cosmological simulations have managed to reproduce realistic and morphologically diverse galaxies, spanning the Hubble sequence. Central to this success was a phenomenological calibration of the few included feedback processes, whilst glossing over higher complexity baryonic physics. This approach diminishes the predictive power of such simulations, preventing to further our understanding of galaxy formation. To tackle this fundamental issue, we investigate the impact of cosmic rays (CRs) and magnetic fields on the interstellar medium (ISM) and the launching of outflows in a cosmological zoom-in simulation of a Milky Way-like galaxy. We find that including CRs decreases the stellar mass of the galaxy by a factor of 10 at high redshift and $\sim 4$ at cosmic noon, leading to a stellar mass to halo mass ratio in good agreement with abundance matching models. Such decrease is caused by two effects: i) a reduction of cold, high-density, star-forming gas, and ii) a larger fraction of SN events exploding at lower densities, where they have a higher impact. SN-injected CRs produce enhanced, multi-phase galactic outflows, which are accelerated by CR pressure gradients in the circumgalactic medium of the galaxy. While the mass budget of these outflows is dominated by the warm ionised gas, warm neutral and cold gas phases contribute significantly at high redshifts. Importantly, our work shows that future JWST observations of galaxies and their multi-phase outflows across cosmic time have the ability to constrain the role of CRs in regulating star formation. △ Less

Submitted 25 July, 2023; originally announced July 2023.

Comments: 20 pages, 11 figures, submitted MNRAS

Simulations of spin-driven AGN jets in gas-rich galaxy mergers

Authors: Rosie Y. Talbot, Debora Sijacki, Martin A. Bourne

Abstract: In this work, we use hydrodynamical simulations to explore the effects of kinetic AGN jet feedback on the progression and outcome of the major merger of two isolated, gas-rich galaxies. We present simulations that use the moving-mesh code AREPO to follow the progression of the merger through first passage and up to the final coalescence, modelling the black holes at the centres of both of the merg… ▽ More In this work, we use hydrodynamical simulations to explore the effects of kinetic AGN jet feedback on the progression and outcome of the major merger of two isolated, gas-rich galaxies. We present simulations that use the moving-mesh code AREPO to follow the progression of the merger through first passage and up to the final coalescence, modelling the black holes at the centres of both of the merging galaxies using our prescription for black hole accretion via an $α$-disc and feedback in the form of a spin-driven jet. We find that the jets drive large-scale, multiphase outflows which launch large quantities of cold gas out to distances greater than 100 kpc and with velocities that reach $\sim 2500 \, {\rm km \, s^{-1}}$. Gas in the outflows that decelerates, cools and falls back on the galaxies can provide a rich source of fuel for the black hole, leading to intense episodes of jet activity in which the jet can become significantly misaligned. The presence of AGN jets affects the growth of the stellar component: star formation is moderately suppressed at all times during the merger and the peak of the star formation rate, attained during the final coalescence of the galaxies, is reduced by a factor of $\sim 2$. Analysis of simulations such as these will play a central role in making precise predictions for multimessenger investigations of dual radio-AGN, which next-generation observational facilities such as LISA, Athena and SKA will make possible. △ Less

Submitted 16 February, 2024; v1 submitted 12 June, 2023; originally announced June 2023.

Comments: 20 pages, 10 figures, accepted for publication by MNRAS

The growth of the gargantuan black holes powering high-redshift quasars and their impact on the formation of early galaxies and protoclusters

Authors: Jake S. Bennett, Debora Sijacki, Tiago Costa, Nicolas Laporte, Callum Witten

Abstract: High-redshift quasars ($z\gtrsim6$), powered by black holes (BHs) with large inferred masses, imply rapid BH growth in the early Universe. The most extreme examples have inferred masses of $\sim \! 10^9\,$M$_\odot$ at $z = 7.5$ and $\sim \! 10^{10}\,$M$_\odot$ at $z = 6.3$. Such dramatic growth via gas accretion likely leads to significant energy input into the quasar host galaxy and its surroundi… ▽ More High-redshift quasars ($z\gtrsim6$), powered by black holes (BHs) with large inferred masses, imply rapid BH growth in the early Universe. The most extreme examples have inferred masses of $\sim \! 10^9\,$M$_\odot$ at $z = 7.5$ and $\sim \! 10^{10}\,$M$_\odot$ at $z = 6.3$. Such dramatic growth via gas accretion likely leads to significant energy input into the quasar host galaxy and its surroundings, however few theoretical predictions of the impact of such objects currently exist. We present zoom-in simulations of a massive high-redshift protocluster, with our fiducial FABLE model incapable of reproducing the brightest quasars. With modifications to this model to promote early BH growth, such as earlier seeding and mildly super-Eddington accretion, such `gargantuan' BHs can be formed. With this new model, simulated host dust masses and star formation rates are in good agreement with existing JWST and ALMA data from ultraluminous quasars. We find the quasar is often obscured as it grows, and that strong, ejective feedback is required to have a high probability of detecting the quasar in the rest-frame UV. Fast and energetic quasar-driven winds expel metal-enriched gas, leading to significant metal pollution of the circumgalactic medium (CGM) out to twice the virial radius. As central gas densities and pressures are reduced, we find weaker signals from the CGM in mock X-ray and Sunyaev-Zeldovich maps, whose detection - with proposed instruments such as Lynx, and even potentially presently with ALMA - can constrain quasar feedback. △ Less

Submitted 8 November, 2023; v1 submitted 19 May, 2023; originally announced May 2023.

Comments: 22 pages, 11 figures, published in MNRAS

Deciphering Lyman-$α$ Emission Deep into the Epoch of Reionisation

Authors: Callum Witten, Nicolas Laporte, Sergio Martin-Alvarez, Debora Sijacki, Yuxuan Yuan, Martin G. Haehnelt, William M. Baker, James S. Dunlop, Richard S. Ellis, Norman A. Grogin, Garth Illingworth, Harley Katz, Anton M. Koekemoer, Daniel Magee, Roberto Maiolino, William McClymont, Pablo G. Pérez-González, David Puskás, Guido Roberts-Borsani, Paola Santini, Charlotte Simmonds

Abstract: During the epoch of reionisation the first galaxies were enshrouded in pristine neutral gas, with one of the brightest emission lines in star-forming galaxies, Lyman-$α$ (Ly$α$), expected to remain undetected until the Universe became ionised. Providing an explanation for the surprising detection of Ly$α$ in these early galaxies is a major challenge for extra-galactic studies. Recent JWST observat… ▽ More During the epoch of reionisation the first galaxies were enshrouded in pristine neutral gas, with one of the brightest emission lines in star-forming galaxies, Lyman-$α$ (Ly$α$), expected to remain undetected until the Universe became ionised. Providing an explanation for the surprising detection of Ly$α$ in these early galaxies is a major challenge for extra-galactic studies. Recent JWST observations have reignited the debate on whether residence in an overdensity of galaxies is a it sufficient and necessary condition for Ly$α$ to escape. Here, we take unique advantage of both high-resolution and high-sensitivity images from the JWST instrument NIRCam to reveal that all galaxies in a sample of z>7 Ly$α$ emitters have close companions. We exploit novel on-the-fly radiative transfer magnetohydrodynamical simulations with cosmic ray feedback to show that galaxies with frequent mergers have very bursty star formation which drives episodes of high intrinsic Ly$α$ emission and facilitates the escape of Ly$α$ photons along channels cleared of neutral gas. We conclude that the rapid build up of stellar mass through mergers presents a compelling solution to the long-standing puzzle of the detection of Ly$α$ emission deep into the epoch of reionisation. △ Less

Submitted 20 January, 2024; v1 submitted 28 March, 2023; originally announced March 2023.

Comments: Published in Nature Astronomy. 45 pages, 9 figures, 2 tables

Shattering and growth of cold clouds in galaxy clusters: the role of radiative cooling, magnetic fields and thermal conduction

Authors: Fred Jennings, Ricarda Beckmann, Debora Sijacki, Yohan Dubois

Abstract: In galaxy clusters, the hot intracluster medium (ICM) can develop a striking multi-phase structure around the brightest cluster galaxy. Much work has been done on understanding the origin of this central nebula, but less work has studied its eventual fate after the originally filamentary structure is broken into individual cold clumps. In this paper we perform a suite of 30 (magneto-)hydrodynamica… ▽ More In galaxy clusters, the hot intracluster medium (ICM) can develop a striking multi-phase structure around the brightest cluster galaxy. Much work has been done on understanding the origin of this central nebula, but less work has studied its eventual fate after the originally filamentary structure is broken into individual cold clumps. In this paper we perform a suite of 30 (magneto-)hydrodynamical simulations of kpc-scale cold clouds with typical parameters as found by galaxy cluster simulations, to understand whether clouds are mixed back into the hot ICM or can persist. We investigate the effects of radiative cooling, small-scale heating, magnetic fields, and (anisotropic) thermal conduction on the long-term evolution of clouds. We find that filament fragments cool on timescales shorter than the crushing timescale, fall out of pressure equilibrium with the hot medium, and shatter, forming smaller clumplets. These act as nucleation sites for further condensation, and mixing via Kelvin-Helmholtz instability, causing cold gas mass to double within 75 Myr. Cloud growth depends on density, as well as on local heating processes, which determine whether clouds undergo ablation- or shattering-driven evolution. Magnetic fields slow down but don't prevent cloud growth, with the evolution of both cold and warm phase sensitive to the field topology. Counter-intuitively, anisotropic thermal conduction increases the cold gas growth rate compared to non-conductive clouds, leading to larger amounts of warm phase as well. We conclude that dense clumps on scales of $500$ pc or more cannot be ignored when studying the long-term cooling flow evolution of galaxy clusters. △ Less

Submitted 30 November, 2022; v1 submitted 16 November, 2022; originally announced November 2022.

Comments: 21 pages, 26 figures, Accepted for Publication in MNRAS

The Pandora project. I: the impact of radiation and cosmic rays on baryonic and dark matter properties of dwarf galaxies

Authors: Sergio Martin-Alvarez, Debora Sijacki, Martin G. Haehnelt, Marion Farcy, Yohan Dubois, Vasily Belokurov, Joakim Rosdahl, Enrique Lopez-Rodriguez

Abstract: Enshrouded in several well-known controversies, dwarf galaxies have been extensively studied to learn about the underlying cosmology, notwithstanding that physical processes regulating their properties are poorly understood. To shed light on these processes, we introduce the Pandora suite of 17 high-resolution (3.5 parsec half-cell side) dwarf galaxy formation cosmological simulations. Commencing… ▽ More Enshrouded in several well-known controversies, dwarf galaxies have been extensively studied to learn about the underlying cosmology, notwithstanding that physical processes regulating their properties are poorly understood. To shed light on these processes, we introduce the Pandora suite of 17 high-resolution (3.5 parsec half-cell side) dwarf galaxy formation cosmological simulations. Commencing with thermo-turbulent star formation and mechanical supernova feedback, we gradually increase the complexity of physics incorporated leading to full-physics models combining magnetism, on-the-fly radiative transfer and the corresponding stellar photoheating, and SN-accelerated cosmic rays. We investigate combinations of these processes, comparing them with observations to constrain what are the main mechanisms determining dwarf galaxy properties. We find hydrodynamical `SN feedback-only' simulations struggle to produce realistic dwarf galaxies, leading either to overquenched or too centrally concentrated, dispersion dominated systems when compared to observed field dwarfs. Accounting for radiation with cosmic rays results in extended and rotationally-supported systems. Spatially `distributed' feedback leads to realistic stellar and HI masses as well as kinematics. Furthermore, resolved kinematic maps of our full-physics models predict kinematically distinct clumps and kinematic misalignments of stars, HI and HII after star formation events. Episodic star formation combined with its associated feedback induces more core-like dark matter central profiles, which our `SN feedback-only' models struggle to achieve. Our results demonstrate the complexity of physical processes required to capture realistic dwarf galaxy properties, making tangible predictions for integral field unit surveys, radio synchrotron emission, and for galaxy and multi-phase interstellar medium properties that JWST will probe. △ Less

Submitted 25 September, 2023; v1 submitted 16 November, 2022; originally announced November 2022.

Comments: Replaced to match the MNRAS published version

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 525, Issue 3, November 2023, Pages 3806-3830

Two can play at that game: constraining the role of supernova and AGN feedback in dwarf galaxies with cosmological zoom-in simulations

Authors: Sophie Koudmani, Debora Sijacki, Matthew C. Smith

Abstract: There is growing observational evidence for dwarf galaxies hosting active galactic nuclei (AGN), including hints of AGN-driven outflows in dwarfs. However, in the common theoretical model of galaxy formation, efficient supernova (SN) feedback is the tool of choice for regulating star formation in the low-mass regime. In this paper, we present a suite of high-resolution cosmological dwarf zoom-in s… ▽ More There is growing observational evidence for dwarf galaxies hosting active galactic nuclei (AGN), including hints of AGN-driven outflows in dwarfs. However, in the common theoretical model of galaxy formation, efficient supernova (SN) feedback is the tool of choice for regulating star formation in the low-mass regime. In this paper, we present a suite of high-resolution cosmological dwarf zoom-in simulations relaxing the assumption of strong SN feedback, with the goal to determine whether more moderate SN feedback in combination with an efficient AGN could be a suitable alternative. Importantly, we find that there are sufficient amounts of gas to power brief Eddington-limited accretion episodes in dwarfs. This leads to a variety of outcomes depending on the AGN accretion model: from no additional suppression to moderate regulation of star formation to catastrophic quenching. Efficient AGN can drive powerful outflows, depleting the gas reservoir of their hosts via ejective feedback and then maintaining a quiescent state through heating the circumgalactic medium. Moderate AGN outflows can be as efficient as the strong SN feedback commonly employed, leading to star formation regulation and HI gas masses in agreement with observations of field dwarfs. All efficient AGN set-ups are associated with overmassive black holes (BHs) compared to the (heavily extrapolated) observed BH mass - stellar mass scaling relations, with future direct observational constraints in this mass regime being crucially needed. Efficient AGN activity is mostly restricted to high redshifts, with hot, accelerated outflows and high X-ray luminosities being the clearest tell-tale signs for future observational campaigns. △ Less

Submitted 1 September, 2022; v1 submitted 22 June, 2022; originally announced June 2022.

Comments: Accepted for publication in MNRAS; 30 pages, 14 figures (including 3 appendices)

Towards convergence of turbulent dynamo amplification in cosmological simulations of galaxies

Authors: Sergio Martin-Alvarez, Julien Devriendt, Adrianne Slyz, Debora Sijacki, Mark L. A. Richardson, Harley Katz

Abstract: Our understanding of the process through which magnetic fields reached their observed strengths in present-day galaxies remains incomplete. One of the advocated solutions is a turbulent dynamo mechanism that rapidly amplifies weak magnetic field seeds to the order of ${\sim}μ$G. However, simulating the turbulent dynamo is a very challenging computational task due to the demanding span of spatial s… ▽ More Our understanding of the process through which magnetic fields reached their observed strengths in present-day galaxies remains incomplete. One of the advocated solutions is a turbulent dynamo mechanism that rapidly amplifies weak magnetic field seeds to the order of ${\sim}μ$G. However, simulating the turbulent dynamo is a very challenging computational task due to the demanding span of spatial scales and the complexity of the required numerical methods. In particular, turbulent velocity and magnetic fields are extremely sensitive to the spatial discretisation of simulated domains. To explore how refinement schemes affect galactic turbulence and amplification of magnetic fields in cosmological simulations, we compare two refinement strategies. A traditional quasi-Lagrangian adaptive mesh refinement approach focusing spatial resolution on dense regions, and a new refinement method that resolves the entire galaxy with a high resolution quasi-uniform grid. Our new refinement strategy yields much faster magnetic energy amplification than the quasi-Lagrangian method, which is also significantly greater than the adiabatic compressional estimate indicating that the extra amplification is produced through stretching of magnetic field lines. Furthermore, with our new refinement the magnetic energy growth factor scales with resolution following $\propto \Dres^{-1/2}$, in much better agreement with small-scale turbulent box simulations. Finally, we find evidence suggesting most magnetic amplification in our simulated galaxies occurs in the warm phase of their interstellar medium, which has a better developed turbulent field with our new refinement strategy. △ Less

Submitted 21 April, 2022; v1 submitted 12 November, 2021; originally announced November 2021.

Comments: Replaced to match version accepted for publication in MNRAS. 20 pages, 17 figures and 2 appendices

Blandford-Znajek jets in galaxy formation simulations: exploring the diversity of outflows produced by spin-driven AGN jets in Seyfert galaxies

Authors: Rosie Y. Talbot, Debora Sijacki, Martin A. Bourne

Abstract: Recent observations of Seyfert galaxies indicate that low power, misaligned jets can undergo significant interaction with the gas in the galactic disc and may be able to drive large-scale, multiphase outflows. We apply our novel sub-grid model for Blandford-Znajek jets to simulations of the central regions of Seyferts, in which a black hole is embedded in a dense, sub-kpc circumnuclear disc (CND)… ▽ More Recent observations of Seyfert galaxies indicate that low power, misaligned jets can undergo significant interaction with the gas in the galactic disc and may be able to drive large-scale, multiphase outflows. We apply our novel sub-grid model for Blandford-Znajek jets to simulations of the central regions of Seyferts, in which a black hole is embedded in a dense, sub-kpc circumnuclear disc (CND) and surrounded by a dilute circumgalactic medium (CGM). We find that the variability of the accretion flow is highly sensitive both to the jet power and to the CND thermodynamics and, ultimately, is determined by the complex interplay between jet-driven outflows and backflows. Even at moderate Eddington ratios, AGN jets are able to significantly alter the thermodynamics and kinematics of CNDs and entrain up to 10% of their mass in the outflow. Mass outflow rates and kinetic powers of the warm outflowing component are in agreement with recent observations for black holes with similar bolometric luminosities, with outflow velocities that are able to reach 500 km/s. Depending on their power and direction, jets are able to drive a wide variety of large-scale outflows, ranging from light, hot and collimated structures to highly mass-loaded, multiphase, bipolar winds. This diversity of jet-driven outflows highlights the importance of applying physically motivated models of AGN feedback to realistic galaxy formation contexts. Such simulations will play a crucial role in accurately interpreting the wealth of data that next generation facilities such as JWST, SKA and Athena will provide. △ Less

Submitted 7 June, 2022; v1 submitted 2 November, 2021; originally announced November 2021.

Comments: 25 pages, 12 figures, 2 appendices, MNRAS accepted

Momentum deposition of supernovae with cosmic rays

Authors: Francisco Rodríguez Montero, Sergio Martin-Alvarez, Debora Sijacki, Adrianne Slyz, Julien Devriendt, Yohan Dubois

Abstract: The cataclysmic explosions of massive stars as supernovae are one of the key ingredients of galaxy formation. However, their evolution is not well understood in the presence of magnetic fields or cosmic rays (CRs). We study the expansion of individual supernova remnants (SNRs) using our suite of 3D hydrodynamical (HD), magnetohydrodynamical (MHD) and CRMHD simulations generated using RAMSES. We ex… ▽ More The cataclysmic explosions of massive stars as supernovae are one of the key ingredients of galaxy formation. However, their evolution is not well understood in the presence of magnetic fields or cosmic rays (CRs). We study the expansion of individual supernova remnants (SNRs) using our suite of 3D hydrodynamical (HD), magnetohydrodynamical (MHD) and CRMHD simulations generated using RAMSES. We explore multiple ambient densities, magnetic fields and fractions of supernova energy deposited as CRs ($χ_{\rm CR}$), accounting for cosmic ray anisotropic diffusion and streaming. All our runs have comparable evolutions until the end of the Sedov-Taylor phase. However, our CRMHD simulations experience an additional CR pressure-driven snowplough phase once the CR energy dominates inside the SNR. We present a model for the final momentum deposited by supernovae that captures this new phase: $p_{\rm SNR} = 2.87\times 10^{5} (χ_{\text{CR}} + 1)^{4.82}\left(\frac{n}{\text{cm}^{-3}}\right)^{-0.196} M_{\odot}$ km s$^{-1}$. Assuming a 10% fraction of SN energy in CRs leads to a 50% boost of the final momentum, with our model predicting even higher impacts at lower ambient densities. The anisotropic diffusion of CRs assuming an initially uniform magnetic field leads to extended gas and cosmic ray outflows escaping from the supernova poles. We also study a tangled initial configuration of the magnetic field, resulting instead in a quasi-isotropic diffusion of CRs and earlier momentum deposition. Finally, synthetic synchrotron observations of our simulations using the POLARIS code show that the local magnetic field configuration in the interstellar medium modifies the overall radio emission morphology and polarisation. △ Less

Submitted 19 October, 2021; originally announced October 2021.

Comments: 17 pages, 11 figures, submitted to MNRAS

A disturbing FABLE of mergers, feedback, turbulence, and mass biases in simulated galaxy clusters

Authors: Jake S. Bennett, Debora Sijacki

Abstract: The use of galaxy clusters as cosmological probes often relies on understanding the properties and evolution of the intracluster medium (ICM). However, the ICM is a complex plasma, regularly stirred by mergers and feedback, with non-negligible bulk and turbulent motions and a non-thermal pressure component, making it difficult to construct a coherent and comprehensive picture. To this end, we use… ▽ More The use of galaxy clusters as cosmological probes often relies on understanding the properties and evolution of the intracluster medium (ICM). However, the ICM is a complex plasma, regularly stirred by mergers and feedback, with non-negligible bulk and turbulent motions and a non-thermal pressure component, making it difficult to construct a coherent and comprehensive picture. To this end, we use the FABLE simulations to investigate how the hydrostatic mass bias is affected by mergers, turbulence, and feedback. Following in detail a single, massive cluster we find the bias varies significantly over cosmic time, rarely staying at the average value found at a particular epoch. Variations of the bias at a given radius are contemporaneous with periods where outflows dominate the mass flux, either due to mergers or interestingly, at high redshift, AGN feedback. The $z=0$ ensemble median mass bias in FABLE is $\sim\!13$ per cent at $R_\mathrm{500}$ and $\sim\!15$ per cent at $R_\mathrm{200}$, but with a large scatter in individual values. In halo central regions, we see an increase in temperature and a decrease in non-thermal pressure support with cosmic time as turbulence thermalises, leading to a reduction in the mass bias within $\sim\!0.2 \, R_\mathrm{200}$. When using a fitted pressure profile, instead of the simulation data, to estimate the bias, we find there can be significant differences, particularly at larger radii and higher redshift. We therefore caution over the use of such fits in future work when comparing with the next generation of X-ray and SZ observations. △ Less

Submitted 28 April, 2022; v1 submitted 14 October, 2021; originally announced October 2021.

Comments: 16 pages, 7 figures, MNRAS accepted

Morphological evolution of supermassive black hole merger hosts and multimessenger signatures

Authors: Colin DeGraf, Debora Sijacki, Tiziana Di Matteo, Kelly Holley-Bockelmann, Greg Snyder, Volker Springel

Abstract: With projects such as Laser Interferometer Space Antenna (LISA) and Pulsar Timing Arrays expected to detect gravitational waves from supermassive black hole mergers in the near future, it is key that we understand what we expect those detections to be, and maximize what we can learn from them. To address this, we study the mergers of supermassive black holes in the Illustris simulation, the overal… ▽ More With projects such as Laser Interferometer Space Antenna (LISA) and Pulsar Timing Arrays expected to detect gravitational waves from supermassive black hole mergers in the near future, it is key that we understand what we expect those detections to be, and maximize what we can learn from them. To address this, we study the mergers of supermassive black holes in the Illustris simulation, the overall rate of mergers, and the correlation between merging black holes and their host galaxies. We find that these mergers occur in typical galaxies along the $M_{\rm{BH}}-M_*$ relation, and that between LISA and PTAs we expect to probe the full range of galaxy masses. As galaxy mergers can trigger increased star formation, we find that galaxies hosting low-mass black hole mergers tend to show a slight increase in star formation rates compared to a mass-matched sample. However, high-mass merger hosts have typical star formation rates, due to a combination of low gas fractions and powerful AGN feedback. Although minor black hole mergers do not correlate with disturbed morphologies, major mergers (especially at high-masses) tend to show morphological evidence of recent galaxy mergers which survives for ~500 Myr. This is on the same scale as the infall/hardening time of the merging black holes, suggesting that electromagnetic followups to gravitational wave signals may not be able to observe this correlation. We further find that incorporating a realistic timescale delay for the black hole mergers could shift the distribution of merger masses toward higher-masses, decreasing the rate of LISA detections while increasing the rate of PTA detections. △ Less

Submitted 28 July, 2021; v1 submitted 1 December, 2020; originally announced December 2020.

Comments: 14 pages, 17 figures. Published in MNRAS

Journal ref: MNRAS (2021), 503, 3629D

Unraveling the origin of magnetic fields in galaxies

Authors: Sergio Martin-Alvarez, Harley Katz, Debora Sijacki, Julien Devriendt, Adrianne Slyz

Abstract: Despite their ubiquity, there are many open questions regarding galactic and cosmic magnetic fields. Specifically, current observational constraints cannot rule out if magnetic fields observed in galaxies were generated in the Early Universe or are of astrophysical nature. Motivated by this we use our magnetic tracers algorithm to investigate whether the signatures of primordial magnetic fields pe… ▽ More Despite their ubiquity, there are many open questions regarding galactic and cosmic magnetic fields. Specifically, current observational constraints cannot rule out if magnetic fields observed in galaxies were generated in the Early Universe or are of astrophysical nature. Motivated by this we use our magnetic tracers algorithm to investigate whether the signatures of primordial magnetic fields persist in galaxies throughout cosmic time. We simulate a Milky Way-like galaxy in four scenarios: magnetised solely by primordial magnetic fields, magnetised exclusively by SN-injected magnetic fields, and two combined primordial + SN magnetisation cases. We find that once primordial magnetic fields with a comoving strength $B_0 >10^{-12}$ G are considered, they remain the primary source of galaxy magnetisation. Our magnetic tracers show that, even combined with galactic sources of magnetisation, when primordial magnetic fields are strong, they source the large-scale fields in the warm metal-poor phase of the simulated galaxy. In this case, the circumgalactic and intergalactic medium can be used to probe $B_0$ without risk of pollution by magnetic fields originated in the galaxy. Furthermore, whether magnetic fields are primordial or astrophysically-sourced can be inferred by studying local gas metallicity. As a result, we predict that future state-of-the-art observational facilities of magnetic fields in galaxies will have the potential to unravel astrophysical and primordial magnetic components of our Universe. △ Less

Submitted 28 July, 2021; v1 submitted 23 November, 2020; originally announced November 2020.

Comments: Accepted for publication in MNRAS. 19 pages, 13 figures and 1 appendix

Journal ref: MNRAS,504, 2517 (2021)

Blandford-Znajek jets in galaxy formation simulations: method and implementation

Authors: Rosie Y. Talbot, Martin A. Bourne, Debora Sijacki

Abstract: Jets launched by active galactic nuclei (AGN) are believed to play a significant role in shaping the properties of galaxies and provide an energetically viable mechanism through which galaxies can become quenched. Here we present a novel AGN feedback model, which we have incorporated into the AREPO code, that evolves the black hole mass and spin as the accretion flow proceeds through a thin $α$-di… ▽ More Jets launched by active galactic nuclei (AGN) are believed to play a significant role in shaping the properties of galaxies and provide an energetically viable mechanism through which galaxies can become quenched. Here we present a novel AGN feedback model, which we have incorporated into the AREPO code, that evolves the black hole mass and spin as the accretion flow proceeds through a thin $α$-disc which we self-consistently couple to a Blandford-Znajek jet. We apply our model to the central region of a typical radio-loud Seyfert galaxy embedded in a hot circumgalactic medium (CGM). We find that jets launched into high pressure environments thermalise efficiently due to the formation of recollimation shocks and the vigorous instabilities that these shocks excite increase the efficiency of the mixing of CGM and jet material. The beams of more overpressured jets, however, are not as readily disrupted by instabilities so the majority of the momentum flux at the jet base is retained out to the head, where the jet terminates in a reverse shock. All jets entrain a significant amount of cold circumnuclear disc material which, while energetically insignificant, dominates the lobe mass together with the hot, entrained CGM material. The jet power evolves significantly due to effective self-regulation by the black hole, fed by secularly-driven, intermittent mass flows. The direction of jets launched directly into the circumnuclear disc changes considerably due to effective Bardeen-Petterson torquing. Interestingly, these jets obliterate the innermost regions of the disc and drive large-scale, multi-phase, turbulent, bipolar outflows. △ Less

Submitted 29 July, 2021; v1 submitted 20 November, 2020; originally announced November 2020.

Comments: 33 pages, 12 figures, 2 appendices, published in MNRAS

Journal ref: MNRAS, Volume 504, Issue 3, July 2021

AGN jet feedback on a moving mesh: gentle cluster heating by weak shocks and lobe disruption

Authors: Martin A. Bourne, Debora Sijacki

Abstract: While there is overwhelming observational evidence of AGN-driven jets in galaxy clusters and groups, if and how the jet energy is delivered to the ambient medium remains unanswered. Here we perform very high resolution AGN jet simulations within a live, cosmologically evolved cluster with the moving mesh code AREPO. We find that mock X-ray and radio lobe properties are in good agreement with obser… ▽ More While there is overwhelming observational evidence of AGN-driven jets in galaxy clusters and groups, if and how the jet energy is delivered to the ambient medium remains unanswered. Here we perform very high resolution AGN jet simulations within a live, cosmologically evolved cluster with the moving mesh code AREPO. We find that mock X-ray and radio lobe properties are in good agreement with observations with different power jets transitioning from FR-I to FR-II-like morphologies. During the lobe inflation phase, heating by both internal and bow shocks contributes to lobe energetics, and $\sim 40$ per cent of the feedback energy goes into the $PdV$ work done by the expanding lobes. Low power jets are more likely to simply displace gas during lobe inflation, but higher power jets become more effective at driving shocks and heating the intracluster medium (ICM), although shocks rarely exceed $\mathcal{M}\sim 2-3$. Once the lobe inflation phase ceases, cluster weather significantly impacts the lobe evolution. Lower power jet lobes are more readily disrupted and mixed with the ICM, depositing up to $\sim 70$ per cent of the injected energy, however, ultimately the equivalent of $\gtrsim 50$ per cent of the feedback energy ends up as potential energy of the system. Even though the mean ICM entropy is increased up to $80$~Myr after the jets switch off, AGN heating is gentle, inducing no large variations in cluster radial profiles in accord with observations. △ Less

Submitted 14 June, 2021; v1 submitted 28 August, 2020; originally announced August 2020.

Comments: 29 pages, 14 figures, accepted for publication in MNRAS

A little FABLE: exploring AGN feedback in dwarf galaxies with cosmological simulations

Authors: Sophie Koudmani, Nicholas A. Henden, Debora Sijacki

Abstract: Contrary to the standard lore, there is mounting observational evidence that feedback from active galactic nuclei (AGN) may also play a role at the low-mass end of the galaxy population. We investigate this using the cosmological simulation suite FABLE, with a particular focus on the dwarf regime ($M_\mathrm{stellar} < 10^{9.5} \ \mathrm{M_{\odot}}$). We find that overmassive black holes (BHs), wi… ▽ More Contrary to the standard lore, there is mounting observational evidence that feedback from active galactic nuclei (AGN) may also play a role at the low-mass end of the galaxy population. We investigate this using the cosmological simulation suite FABLE, with a particular focus on the dwarf regime ($M_\mathrm{stellar} < 10^{9.5} \ \mathrm{M_{\odot}}$). We find that overmassive black holes (BHs), with respect to the mean scaling relations with their host galaxies, drive hotter and faster outflows and lead to significantly reduced gas mass fractions. They are also more likely to display a kinematically misaligned ionized gas component in our mock MaNGA velocity maps, although we caution that cosmic inflows and mergers contribute to misalignments as well. While in the local Universe the majority of AGN in dwarfs are much dimmer than the stellar component, for $z \geq 2$ there is a significant population that outshines their hosts. These high-redshift overmassive BHs contribute to the quenching of dwarfs, whereas at late cosmic times supernova (SN) feedback is more efficient. While our results are overall in good agreement with X-ray observations of AGN in dwarfs, the lack of high-luminosity X-ray AGN in FABLE at low redshifts highlights an interesting possibility that SN feedback could be too strong in FABLE's dwarfs, curtailing AGN growth and feedback. We predict that future observations may uncover many more AGN in dwarfs with lower luminosities and at higher redshifts. △ Less

Submitted 18 March, 2021; v1 submitted 20 July, 2020; originally announced July 2020.

Comments: Accepted for publication in MNRAS; 24 pages, 12 figures, 1 appendix

Resolving shocks and filaments in galaxy formation simulations: effects on gas properties and star formation in the circumgalactic medium

Authors: Jake S. Bennett, Debora Sijacki

Abstract: There is an emerging consensus that large amounts of gas do not shock heat in the circumgalactic medium (CGM) of massive galaxies, but instead pierce deep into haloes from the cosmic web via filaments. To better resolve this process numerically, we have developed a novel `shock refinement' scheme within the moving mesh code AREPO that adaptively improves resolution around shocks on-the-fly in gala… ▽ More There is an emerging consensus that large amounts of gas do not shock heat in the circumgalactic medium (CGM) of massive galaxies, but instead pierce deep into haloes from the cosmic web via filaments. To better resolve this process numerically, we have developed a novel `shock refinement' scheme within the moving mesh code AREPO that adaptively improves resolution around shocks on-the-fly in galaxy formation simulations. We apply this to a massive $\sim10^{12}$ M$_\odot$ halo at $z=6$ using the successful FABLE model, increasing the mass resolution by a factor of 512. With better refinement there are significantly more dense, metal-poor and fast-moving filaments and clumps flowing into the halo, leading to a more multiphase CGM. We find a $\sim50$ per cent boost in cool-dense gas mass and a 25 per cent increase in inflowing mass flux. Better resolved accretion shocks cause turbulence to increase dramatically, leading to a doubling in the halo's non-thermal pressure support. Despite much higher thermalisation at shocks with higher resolution, increased cooling rates suppress the thermal energy of the halo. In contrast, the faster and denser filaments cause a significant jump in the bulk kinetic energy of cool-dense gas, while in the hot phase turbulent energy increases by up to $\sim150$ per cent. Moreover, HI covering fractions within the CGM increase by up to 60 per cent. Consequently star formation is spread more widely and we predict a population of metal-poor stars forming within primordial filaments that deep JWST observations may be able to probe. △ Less

Submitted 24 September, 2020; v1 submitted 17 June, 2020; originally announced June 2020.

Comments: 21 pages, 14 figures. Accepted for publication in MNRAS after minor revisions

Metal enrichment in the circumgalactic medium and Lyα haloes around quasars at z $\sim$3

Authors: Yucheng Guo, Roberto Maiolino, Linhua Jiang, Kenta Matsuoka, Tohru Nagao, Oli Luiz Dors, Michele Ginolfi, Nick Henden, Jake Bennett, Debora Sijacki, Ewald Puchwein

Abstract: Deep observations have detected extended Ly$α$ emission nebulae surrounding tens of quasars at redshift 2 to 6. However, the metallicity of such extended haloes is still poorly understood. We perform a detailed analysis on a large sample of 80 quasars at $z\sim3$ based on MUSE-VLT data. We find clear evidence of extended emission of the UV nebular lines such as CIV or HeII for about 20$\%$ of the… ▽ More Deep observations have detected extended Ly$α$ emission nebulae surrounding tens of quasars at redshift 2 to 6. However, the metallicity of such extended haloes is still poorly understood. We perform a detailed analysis on a large sample of 80 quasars at $z\sim3$ based on MUSE-VLT data. We find clear evidence of extended emission of the UV nebular lines such as CIV or HeII for about 20$\%$ of the sample, while CIII] is only marginally detected in a few objects. By stacking the cubes we detect emission of CIV, HeII and CIII] out to a radius of about 45 kpc. CIV and HeII show a radial decline much steeper than Ly$α$, while CIII] shows a shallower profile similar to Ly$α$ in the inner 45 kpc. We infer that the average metallicity of the circumgalactic gas within the central 30-50~kpc is $\sim$0.5 solar, or even higher. However, we also find evidence of a component of the Ly$α$ haloes, which has much weaker metal emission lines relative to Ly$α$. We suggest that the high metallicity of the circumgalactic medium within the central 30-50 kpc is associated with chemical pre-enrichment by past quasar-driven outflows and that there is a more extended component of the CGM that has much lower metallicity and likely associated with near-pristine gas accreted from the intergalactic medium. We show that our observational results are in good agreement with the expectations of the FABLE zoom-in cosmological simulations. △ Less

Submitted 14 June, 2020; v1 submitted 15 January, 2020; originally announced January 2020.

Comments: 24 pages, 17 figures. Accepted for publication in ApJ

The baryon content of groups and clusters of galaxies in the FABLE simulations

Authors: Nicholas A. Henden, Ewald Puchwein, Debora Sijacki

Abstract: We study the gas and stellar mass content of galaxy groups and clusters in the FABLE suite of cosmological hydrodynamical simulations, including the evolution of their central brightest cluster galaxies (BCGs), satellite galaxies and intracluster light (ICL). The total gas and stellar mass of FABLE clusters are in very good agreement with observations and show negligible redshift evolution at fixe… ▽ More We study the gas and stellar mass content of galaxy groups and clusters in the FABLE suite of cosmological hydrodynamical simulations, including the evolution of their central brightest cluster galaxies (BCGs), satellite galaxies and intracluster light (ICL). The total gas and stellar mass of FABLE clusters are in very good agreement with observations and show negligible redshift evolution at fixed halo mass for $M_{500} \gtrsim 3 \times 10^{14} M_{\odot}$ at $z \lesssim 1$, in line with recent findings from Sunyaev-Zel'dovich (SZ)-selected cluster samples. Importantly, the simulations predict significant redshift evolution in these quantities in the low mass ($M_{500} \sim 10^{14} M_{\odot}$) regime, which will be testable with upcoming SZ surveys such as SPT-3G. While the stellar masses of FABLE BCGs are in reasonable agreement with observations, the total stellar mass in satellite galaxies is lower than observed and the total mass in ICL is somewhat higher. This may be caused by enhanced tidal stripping of satellite galaxies due to their large sizes. BCGs are characterised by moderate stellar mass growth at $z < 1$ coincident with a late-time development of the ICL. The level of BCG mass growth is in good agreement with recent observations, however, we caution that the inferred growth depends sensitively on the mass definition. We further show that in-situ star formation contributes more than half the mass of a BCG over its lifetime, the bulk of which is gained at $z > 1$ where star formation rates are highest. The stellar mass profiles of the BCG+ICL component are similar to observed profiles out to $\sim 100$ kpc at $z \approx 0$ and follow a close to power law shape out to several hundred kpc. We further demonstrate that the inferred size growth of BCGs can be severely biased by the choice of parametric model and the outer radius of the fit. △ Less

Submitted 25 July, 2020; v1 submitted 27 November, 2019; originally announced November 2019.

Comments: 27 pages, 13 figures, 1 appendix; MNRAS accepted

Cosmological simulations of massive black hole seeds: predictions for next generation electromagnetic and gravitational wave observations

Authors: Colin DeGraf, Debora Sijacki

Abstract: We study how statistical properties of supermassive black holes depend on the frequency and conditions for massive seed formation in cosmological simulations of structure formation. We develop a novel method to recalculate detailed growth histories and merger trees of black holes within the framework of the Illustris simulation for several seed formation models, including a physically motivated mo… ▽ More We study how statistical properties of supermassive black holes depend on the frequency and conditions for massive seed formation in cosmological simulations of structure formation. We develop a novel method to recalculate detailed growth histories and merger trees of black holes within the framework of the Illustris simulation for several seed formation models, including a physically motivated model where black hole seeds only form in progenitor galaxies that conform to the conditions for direct collapse black hole formation. While all seed models considered here are in a broad agreement with present observational constraints on black hole populations from optical, UV and X-ray studies, we find that they lead to widely different black hole number densities and halo occupation fractions which are currently observationally unconstrained. In terms of future electromagnetic spectrum observations, the faint-end quasar luminosity function and the low-mass-end black hole-host galaxy scaling relations are very sensitive to the specific massive seed prescription. Specifically, the direct collapse model exhibits a seeding efficiency which decreases rapidly with cosmic time and produces much fewer black holes in low mass galaxies, in contrast to the original Illustris simulation. We further find that the total black hole merger rate varies by more than one order of magnitude for different seed models, with the redshift evolution of the chirp mass changing as well. Supermassive black hole merger detections with LISA and International Pulsar Timing Array may hence provide the most direct means of constraining massive black hole seed formation in the early Universe. △ Less

Submitted 29 July, 2021; v1 submitted 26 June, 2019; originally announced June 2019.

Comments: 20 pages, 17 figures; Published in MNRAS

Journal ref: MNRAS (2020), 491,4973

Constraints on Chameleon f(R)-Gravity from Galaxy Rotation Curves of the SPARC Sample

Authors: Aneesh P. Naik, Ewald Puchwein, Anne-Christine Davis, Debora Sijacki, Harry Desmond

Abstract: In chameleon $f(R)$-gravity, the fifth force will lead to `upturns' in galaxy rotation curves near the screening radius. The location of the upturn depends on the cosmic background value of the scalar field $f_{R0}$, as well as the mass, size and environment of the galaxy. We search for this signature of modified gravity in the SPARC sample of measured rotation curves, using an MCMC technique to d… ▽ More In chameleon $f(R)$-gravity, the fifth force will lead to `upturns' in galaxy rotation curves near the screening radius. The location of the upturn depends on the cosmic background value of the scalar field $f_{R0}$, as well as the mass, size and environment of the galaxy. We search for this signature of modified gravity in the SPARC sample of measured rotation curves, using an MCMC technique to derive constraints on $f_{R0}$. Assuming NFW dark matter haloes and with $f_{R0}$ freely varying for each galaxy, most galaxies prefer $f(R)$ gravity to $Λ$CDM, but there is a large spread of inferred $f_{R0}$ values, inconsistent with a single global value. Requiring instead a consistent $f_{R0}$ value for the whole sample, models with $\log_{10}|f_{R0}| > -6.1$ are excluded. On the other hand, models in the range $-7.5<\log_{10}|f_{R0}|<-6.5$ seem to be favoured with respect to $Λ$CDM, with a significant peak at -7. However, this signal is largely a result of galaxies for which the $f(R)$ signal is degenerate with the core/cusp problem, and when the NFW profile is replaced with a cored halo profile, $Λ$CDM gives better fits than any given $f(R)$ model. Thus, we find no convincing evidence of $f(R)$ gravity down to the level of $|f_{R0}| \sim 6 \times 10^{-8}$, with the caveat that if cored halo density profiles cannot ultimately be explained within $Λ$CDM, a screened modified gravity theory could possibly provide an alternative solution for the core/cusp problem. However, the $f(R)$ models studied here fall short of achieving this. △ Less

Submitted 31 July, 2019; v1 submitted 30 May, 2019; originally announced May 2019.

Comments: Updated to match version accepted by MNRAS. 19 pages, 11 figures (inc. appendix). All code and plotting scripts are publicly available

The redshift evolution of X-ray and Sunyaev-Zel'dovich scaling relations in the FABLE simulations

Authors: Nicholas A. Henden, Ewald Puchwein, Debora Sijacki

Abstract: We study the redshift evolution of the X-ray and Sunyaev-Zel'dovich (SZ) scaling relations for galaxy groups and clusters in the FABLE suite of cosmological hydrodynamical simulations. Using an expanded sample of $27$ high-resolution zoom-in simulations, together with a uniformly-sampled cosmological volume to sample low-mass systems, we find very good agreement with the majority of observational… ▽ More We study the redshift evolution of the X-ray and Sunyaev-Zel'dovich (SZ) scaling relations for galaxy groups and clusters in the FABLE suite of cosmological hydrodynamical simulations. Using an expanded sample of $27$ high-resolution zoom-in simulations, together with a uniformly-sampled cosmological volume to sample low-mass systems, we find very good agreement with the majority of observational constraints up to $z \sim 1$. We predict significant deviations of all examined scaling relations from the simple self-similar expectations. While the slopes are approximately independent of redshift, the normalisations evolve positively with respect to self-similarity, even for commonly-used mass proxies such as the $Y_{\mathrm{X}}$ parameter. These deviations are due to a combination of factors, including more effective AGN feedback in lower mass haloes, larger binding energy of gas at a given halo mass at higher redshifts and larger non-thermal pressure support from kinetic motions at higher redshifts. Our results have important implications for cluster cosmology from upcoming SZ surveys such as SPT-3G, ACTpol and CMB-S4, as relatively small changes in the observable--mass scaling relations (within theoretical uncertainties) have a large impact on the predicted number of high-redshift clusters and hence on our ability to constrain cosmology using cluster abundances. In addition, we find that the intrinsic scatter of the relations, which agrees well with most observational constraints, increases at lower redshifts and for lower mass systems. This calls for a more complex parametrization than adopted in current observational studies to be able to accurately account for selection biases. △ Less

Submitted 16 August, 2019; v1 submitted 30 April, 2019; originally announced May 2019.

Comments: 35 pages, 13 figures, 2 appendices; MNRAS accepted

AGN jet feedback on a moving mesh: lobe energetics and X-ray properties in a realistic cluster environment

Authors: Martin A. Bourne, Debora Sijacki, Ewald Puchwein

Abstract: Jet feedback from active galactic nuclei (AGN) harboured by brightest cluster galaxies is expected to play a fundamental role in regulating cooling in the intracluster medium (ICM). While observations and theory suggest energy within jet lobes balances ICM radiative losses, the modus operandi of energy communication with the ICM remains unclear. We present simulations of very high resolution AGN-d… ▽ More Jet feedback from active galactic nuclei (AGN) harboured by brightest cluster galaxies is expected to play a fundamental role in regulating cooling in the intracluster medium (ICM). While observations and theory suggest energy within jet lobes balances ICM radiative losses, the modus operandi of energy communication with the ICM remains unclear. We present simulations of very high resolution AGN-driven jets launching in a live, cosmological galaxy cluster, within the moving mesh code AREPO. As the jet propagates through the ICM the majority of its energy, which is initially in the kinetic form, thermalizes quickly through internal shocks and inflates lobes of very hot gas. The jets effectively heat the cluster core, with $PdV$ work and weather-aided mixing being the main channels of energy transfer from the lobes to the ICM, while strong shocks and turbulence are subdominant. We additionally present detailed mock X-ray maps at different stages of evolution, revealing clear cavities surrounded by X-ray bright rims, with lobes being detectable for up to $\sim 10^8$ yr even when magnetic draping is ineffective. We find bulk motions in the cluster can significantly affect lobe propagation, offsetting them from the jet direction and imparting bulk velocities that can dominate over the buoyantly rising motion. △ Less

Submitted 16 October, 2019; v1 submitted 30 January, 2019; originally announced January 2019.

Comments: 8 pages, 3 figures, accepted for publication in MNRAS

Journal ref: MNRAS, Volume 490, Issue 1, November 2019

Fast and energetic AGN-driven outflows in simulated dwarf galaxies

Authors: Sophie Koudmani, Debora Sijacki, Martin A. Bourne, Matthew C. Smith

Abstract: The systematic analysis of optical large-scale surveys has revealed a population of dwarf galaxies hosting AGN, which have been confirmed by X-ray follow-up observations. Recently, the MaNGA survey identified six dwarf galaxies that appear to have an AGN that is preventing on-going star formation. It is therefore timely to study the physical properties of dwarf galaxies, in particular whether the… ▽ More The systematic analysis of optical large-scale surveys has revealed a population of dwarf galaxies hosting AGN, which have been confirmed by X-ray follow-up observations. Recently, the MaNGA survey identified six dwarf galaxies that appear to have an AGN that is preventing on-going star formation. It is therefore timely to study the physical properties of dwarf galaxies, in particular whether the presence of an AGN can affect their evolution. Using the moving mesh code AREPO, we have investigated different models of AGN activity, ranging from simple energy-driven spherical winds to collimated, mass-loaded, bipolar outflows in high resolution simulations of isolated dwarf galaxies hosting an active black hole. Our simulations also include a novel implementation of star formation and mechanical supernova (SN) feedback. We find that AGN outflows have a small but systematic effect on the central star formation rates (SFRs) for all set-ups explored, while substantial effects on the global SFR are only obtained with strong SNe and a sustained high-luminosity AGN with an isotropic wind. This suggests that AGN feedback in dwarf galaxies is unlikely to directly regulate their global SFRs. There is, however, a significant effect on outflow properties, which are notably enhanced by the AGN to much higher outflow temperatures and velocities, in agreement with kinematic signatures from the MaNGA survey. This indicates that AGN may play an indirect role in regulating the baryon cycle in dwarf galaxies by hindering cosmic gas inflows. △ Less

Submitted 8 January, 2019; v1 submitted 11 December, 2018; originally announced December 2018.

Comments: Accepted for publication in MNRAS, 22 pages, 11 figures

Cosmological simulations of dwarfs: the need for ISM physics beyond SN feedback alone

Authors: Matthew C. Smith, Debora Sijacki, Sijing Shen

Abstract: The dominant feedback mechanism in low mass haloes is usually assumed to take the form of massive stars exploding as supernovae (SNe). We perform very high resolution cosmological zoom-in simulations of five dwarf galaxies to z = 4 with our mechanical SN feedback model. This delivers the correct amount of momentum corresponding to the stage of the SN remnant evolution resolved, and has been shown… ▽ More The dominant feedback mechanism in low mass haloes is usually assumed to take the form of massive stars exploding as supernovae (SNe). We perform very high resolution cosmological zoom-in simulations of five dwarf galaxies to z = 4 with our mechanical SN feedback model. This delivers the correct amount of momentum corresponding to the stage of the SN remnant evolution resolved, and has been shown to lead to realistic dwarf properties in isolated simulations. We find that in 4 out of our 5 simulated cosmological dwarfs, SN feedback has insufficient impact resulting in excessive stellar masses, extremely compact sizes and central super-solar stellar metallicities. The failure of the SN feedback in our dwarfs is physical in nature within our model and is the result of the build up of very dense gas in the early universe due to mergers and cosmic inflows prior to the first SN occurring. We demonstrate that our results are insensitive to resolution (provided that it is high enough), details of the (spatially uniform) UV background and reasonable alterations within our star formation prescription. We therefore conclude that the ability of SNe to regulate dwarf galaxy properties is dependent on other physical processes, such as turbulent pressure support, clustering and runaway of SN progenitors and other sources of stellar feedback. △ Less

Submitted 2 April, 2019; v1 submitted 11 July, 2018; originally announced July 2018.

Comments: 18 pages, 13 figures, accepted MNRAS

Journal ref: 2019 MNRAS 485, 3317

Tracing the Sources of Reionization in Cosmological Radiation Hydrodynamics Simulations

Authors: Harley Katz, Taysun Kimm, Martin G. Haehnelt, Debora Sijacki, Joakim Rosdahl, Jeremy Blaizot

Abstract: We use the photon flux and absorption tracer algorithm presented in Katz et al.~2018, to characterise the contribution of haloes of different mass and stars of different age and metallicity to the reionization of the Universe. We employ a suite of cosmological multifrequency radiation hydrodynamics AMR simulations that are carefully calibrated to reproduce a realistic reionization history and gala… ▽ More We use the photon flux and absorption tracer algorithm presented in Katz et al.~2018, to characterise the contribution of haloes of different mass and stars of different age and metallicity to the reionization of the Universe. We employ a suite of cosmological multifrequency radiation hydrodynamics AMR simulations that are carefully calibrated to reproduce a realistic reionization history and galaxy properties at $z \geq 6$. In our simulations, haloes with mass $10^9{\rm M_{\odot}}h^{-1}<M<10^{10}{\rm M_{\odot}}h^{-1}$, stars with metallicity $10^{-3}Z_{\odot}<Z<10^{-1.5}Z_{\odot}$, and stars with age $3\,{\rm Myr} < t < 10 \, {\rm Myr}$ dominate reionization by both mass and volume. We show that the sources that reionize most of the volume of the Universe by $z=6$ are not necessarily the same sources that dominate the meta-galactic UV background at the same redshift. We further show that in our simulations, the contribution of each type of source to reionization is not uniform across different gas phases. The IGM, CGM, filaments, ISM, and rarefied supernova heated gas have all been photoionized by different classes of sources. Collisional ionisation contributes at both the lowest and highest densities. In the early stages of the formation of individual HII bubbles, reionization proceeds with the formation of concentric shells of gas ionised by different classes of sources, leading to large temperature variations as a function of galacto-centric radius. The temperature structure of individual HII bubbles may thus give insight into the star formation history of the galaxies acting as the first ionising sources. Our explorative simulations highlight how the complex nature of reionization can be better understood by using our photon tracer algorithm. △ Less

Submitted 20 November, 2018; v1 submitted 19 June, 2018; originally announced June 2018.

Comments: 13 pages, 8 figures, Accepted to MNRAS

The FABLE simulations: A feedback model for galaxies, groups and clusters

Authors: Nicholas A. Henden, Ewald Puchwein, Sijing Shen, Debora Sijacki

Abstract: We present the Feedback Acting on Baryons in Large-scale Environments (FABLE) suite of cosmological hydrodynamical simulations of galaxies, groups and clusters. The simulations use the AREPO moving-mesh code with a set of physical models for galaxy formation based on the successful Illustris simulation, but with updated AGN and supernovae feedback models. This allows us to simultaneously reproduce… ▽ More We present the Feedback Acting on Baryons in Large-scale Environments (FABLE) suite of cosmological hydrodynamical simulations of galaxies, groups and clusters. The simulations use the AREPO moving-mesh code with a set of physical models for galaxy formation based on the successful Illustris simulation, but with updated AGN and supernovae feedback models. This allows us to simultaneously reproduce the observed redshift evolution of the galaxy stellar mass function together with the stellar and gas mass fractions of local groups and clusters across a wide range of halo masses. Focusing on the properties of groups and clusters, we find very good agreement with a range of observed scaling relations, including the X-ray luminosity--total mass and gas mass relations as well as the total mass--temperature and Sunyaev-Zel'dovich flux--mass relations. Careful comparison of our results with scaling relations based on X-ray hydrostatic masses as opposed to weak lensing-derived masses reveals some discrepancies, which hint towards a non-negligible X-ray mass bias in observed samples. We further show that radial profiles of density, pressure and temperature of the simulated intracluster medium are in very good agreement with observations, in particular for $r > 0.3\ r_{500}$. In the innermost regions however we find too large entropy cores, which indicates that a more sophisticated modelling of the physics of AGN feedback may be required to accurately reproduce the observed populations of cool-core and non-cool-core clusters. △ Less

Submitted 6 July, 2018; v1 submitted 13 April, 2018; originally announced April 2018.

Comments: 31 pages, 20 figures, 2 appendices; MNRAS accepted

A Census of the LyC Photons that Form the UV Background During Reionization

Authors: Harley Katz, Taysun Kimm, Martin Haehnelt, Debora Sijacki, Joakim Rosdahl, Jeremy Blaizot

Abstract: We present a new, on-the-fly photon flux and absorption tracer algorithm designed to directly measure the contribution of different source populations to the metagalactic UV background and to the ionisation fraction of gas in the Universe. We use a suite of multifrequency radiation hydrodynamics simulations that are carefully calibrated to reproduce a realistic reionization history and galaxy prop… ▽ More We present a new, on-the-fly photon flux and absorption tracer algorithm designed to directly measure the contribution of different source populations to the metagalactic UV background and to the ionisation fraction of gas in the Universe. We use a suite of multifrequency radiation hydrodynamics simulations that are carefully calibrated to reproduce a realistic reionization history and galaxy properties at $z \ge 6$, to disentangle the contribution of photons emitted by different mass haloes and by stars with different metallicities and ages to the UV background during reionization. While at very early cosmic times low mass, metal poor haloes provide most of the LyC photons, their contribution decreases steadily with time. At $z = 6$ it is the photons emitted by massive systems (${\rm M_{halo}}/{\rm M_\odot} > 10^{10} \, {\rm h ^{-1}}$) and by the metal enriched stars ($10^{-3} < Z/Z_{\rm \odot} < 10^{-1.5}$) that provide the largest contribution to the ionising UV background. We demonstrate that there are large variations in the escape fraction depending on the source, with the escape fraction being highest ($\sim 45-60\%$) for photons emitted by the oldest stars that penetrate into the IGM via low opacity channels carved by the ionising photons and supernova from younger stars. Before HII regions begin to overlap, the photoionisation rate strongly fluctuates between different, isolated HII bubbles, depending on the embedded ionising source, which we suggest may result in spatial variations in the properties of dwarf galaxies. △ Less

Submitted 31 May, 2018; v1 submitted 5 February, 2018; originally announced February 2018.

Comments: 19 pages, 14 figures, Accepted MNRAS

Galactic nuclei evolution with spinning black holes: method and implementation

Authors: Davide Fiacconi, Debora Sijacki, J. E. Pringle

Abstract: Supermassive black holes at the centre of galactic nuclei mostly grow in mass through gas accretion over cosmic time. This process also modifies the angular momentum (or spin) of black holes, both in magnitude and in orientation. Despite being often neglected in galaxy formation simulations, spin plays a crucial role in modulating accretion power, driving jet feedback, and determining recoil veloc… ▽ More Supermassive black holes at the centre of galactic nuclei mostly grow in mass through gas accretion over cosmic time. This process also modifies the angular momentum (or spin) of black holes, both in magnitude and in orientation. Despite being often neglected in galaxy formation simulations, spin plays a crucial role in modulating accretion power, driving jet feedback, and determining recoil velocity of coalescing black hole binaries. We present a new accretion model for the moving-mesh code {\sc arepo} that incorporates (i) mass accretion through a thin $α$-disc, and (ii) spin evolution through the Bardeen-Petterson effect. We use a diverse suite of idealised simulations to explore the physical connection between spin evolution and larger scale environment. We find that black holes with mass $\lesssim 10^{7}$ M$_{\odot}$ experience quick alignment with the accretion disc. This favours prolonged phases of spin-up, and the spin direction evolves according to the gas inflow on timescales as short as $\lesssim 100$ Myr, which might explain the observed jet direction distribution in Seyfert galaxies. Heavier black holes ($\gtrsim 10^{8}$ M$_{\odot}$) are instead more sensitive to the local gas kinematic. Here we find a wider distribution in spin magnitudes: spin-ups are favoured if gas inflow maintains a preferential direction, and spin-downs occur for nearly isotropic infall, while the spin direction does not change much over short timescales $\sim 100$ Myr. We therefore conclude that supermassive black holes with masses $\gtrsim 5 \times 10^{8}$ M$_{\odot}$ may be the ideal testbed to determine the main mode of black hole fuelling over cosmic time. △ Less

Submitted 5 April, 2018; v1 submitted 30 November, 2017; originally announced December 2017.

Comments: accepted for publication in MNRAS (new version after acceptance, some changes to the text, results unchanged)

Quenching star formation with quasar outflows launched by trapped IR radiation

Authors: Tiago Costa, Joakim Rosdahl, Debora Sijacki, Martin G. Haehnelt

Abstract: We present cosmological radiation-hydrodynamic simulations, performed with the code Ramses-RT, of radiatively-driven outflows in a massive quasar host halo at $z = 6$. Our simulations include both single- and multi-scattered radiation pressure on dust from a quasar and are compared against simulations performed with thermal feedback. For radiation pressure-driving, we show that there is a critical… ▽ More We present cosmological radiation-hydrodynamic simulations, performed with the code Ramses-RT, of radiatively-driven outflows in a massive quasar host halo at $z = 6$. Our simulations include both single- and multi-scattered radiation pressure on dust from a quasar and are compared against simulations performed with thermal feedback. For radiation pressure-driving, we show that there is a critical quasar luminosity above which a galactic outflow is launched, set by the equilibrium of gravitational and radiation forces. While this critical luminosity is unrealistically high in the single-scattering limit for plausible black hole masses, it is in line with a $\approx 3 \times 10^9 \, \rm M_\odot$ black hole accreting at its Eddington limit, if infrared (IR) multi-scattering radiation pressure is included. The outflows are fast ($v \, \gtrsim \, 1000 \, \rm km \, s^{-1}$) and strongly mass-loaded with peak mass outflow rates $\approx 10^3 - 10^4 \, \rm M_\odot \, yr^{-1}$, but short-lived ($< 10 \, \rm Myr$). Outflowing material is multi-phase, though predominantly composed of cool gas, forming via a thermal instability in the shocked swept-up component. Radiation pressure- and thermally-driven outflows both affect their host galaxies significantly, but in different, complementary ways. Thermally-driven outflows couple more efficiently to diffuse halo gas, generating more powerful, hotter and more volume-filling outflows. IR radiation, through its ability to penetrate dense gas via diffusion, is more efficient at ejecting gas from the bulge. The combination of gas ejection through outflows with internal pressurisation by trapped IR radiation leads to a complete shut down of star formation in the bulge. We hence argue that radiation pressure-driven feedback may be an important ingredient in regulating star formation in compact starbursts, especially during the quasar's `obscured' phase. △ Less

Submitted 29 May, 2018; v1 submitted 25 September, 2017; originally announced September 2017.

Comments: 34 pages (26 main paper + 8 pages of references and appendices). Accepted by MNRAS

Impact of Cosmological Satellites on Stellar Discs: Dissecting One Satellite at a Time

Authors: Shaoran Hu, Debora Sijacki

Abstract: Within the standard hierarchical structure formation scenario, Milky Way-mass dark matter haloes have hundreds of dark matter subhaloes with mass $\gtrsim 10^8 \, {\rm M_{\odot}}$. Over the lifetime of a galactic disc a fraction of these may pass close to the central region and interact with the disc. We extract the properties of subhaloes, such as their mass and trajectories, from a realistic cos… ▽ More Within the standard hierarchical structure formation scenario, Milky Way-mass dark matter haloes have hundreds of dark matter subhaloes with mass $\gtrsim 10^8 \, {\rm M_{\odot}}$. Over the lifetime of a galactic disc a fraction of these may pass close to the central region and interact with the disc. We extract the properties of subhaloes, such as their mass and trajectories, from a realistic cosmological simulation to study their potential effect on stellar discs. We find that massive subhalo impacts can generate disc heating, rings, bars, warps, lopsidedness as wells as spiral structures in the disc. Specifically, strong counter-rotating single-armed spiral structures form each time a massive subhalo passes through the disc. Such single-armed spirals wind up relatively quickly (over $1-2$ Gyrs) and are generally followed by co-rotating two-armed spiral structures that both develop and wind up more slowly. In our simulations self-gravity in the disc is not very strong and these spiral structures are found to be kinematic density waves. We demonstrate that there is a clear link between each spiral mode in the disc and a given subhalo that caused it, and by changing the mass of the subhalo we can modulate the strength of the spirals. Furthermore, we find that the majority of subhaloes interact with the disc impulsively, such that the strength of spirals generated by subhaloes is proportional to the total torque they exert. We conclude that only a handful of encounters with massive subhaloes is sufficient for re-generating and sustaining spiral structures in discs over their entire lifetime. △ Less

Submitted 10 June, 2018; v1 submitted 18 September, 2017; originally announced September 2017.

Comments: 21 pages, 15 figures. MNRAS accepted

Journal ref: MNRAS, Volume 478, Issue 2, 1 August 2018, Page 1576

Supernova feedback in numerical simulations of galaxy formation: separating physics from numerics

Authors: Matthew C. Smith, Debora Sijacki, Sijing Shen

Abstract: While feedback from massive stars exploding as supernovae (SNe) is thought to be one of the key ingredients regulating galaxy formation, theoretically it is still unclear how the available energy couples to the interstellar medium and how galactic scale outflows are launched. We present a novel implementation of six sub-grid SN feedback schemes in the moving-mesh code Arepo, including injections o… ▽ More While feedback from massive stars exploding as supernovae (SNe) is thought to be one of the key ingredients regulating galaxy formation, theoretically it is still unclear how the available energy couples to the interstellar medium and how galactic scale outflows are launched. We present a novel implementation of six sub-grid SN feedback schemes in the moving-mesh code Arepo, including injections of thermal and/or kinetic energy, two parametrizations of delayed cooling feedback and a `mechanical' feedback scheme that injects the correct amount of momentum depending on the relevant scale of the SN remnant resolved. All schemes make use of individually time-resolved SN events. Adopting isolated disk galaxy setups at different resolutions, with the highest resolution runs reasonably resolving the Sedov-Taylor phase of the SN, we aim to find a physically motivated scheme with as few tunable parameters as possible. As expected, simple injections of energy overcool at all but the highest resolution. Our delayed cooling schemes result in overstrong feedback, destroying the disk. The mechanical feedback scheme is efficient at suppressing star formation, agrees well with the Kennicutt-Schmidt relation and leads to converged star formation rates and galaxy morphologies with increasing resolution without fine tuning any parameters. However, we find it difficult to produce outflows with high enough mass loading factors at all but the highest resolution, indicating either that we have oversimplified the evolution of unresolved SN remnants, require other stellar feedback processes to be included, require a better star formation prescription or most likely some combination of these issues. △ Less

Submitted 26 April, 2018; v1 submitted 11 September, 2017; originally announced September 2017.

Comments: 32 pages, 22 figures. Accepted for publication in MNRAS

The effect of stellar and AGN feedback on the low redshift Lyman-$α$ forest in the Sherwood simulation suite

Authors: Fahad Nasir, James S. Bolton, Matteo Viel, Tae-Sun Kim, Martin G. Haehnelt, Ewald Puchwein, Debora Sijacki

Abstract: We study the effect of different feedback prescriptions on the properties of the low redshift ($z\leq1.6$) Ly$α$ forest using a selection of hydrodynamical simulations drawn from the Sherwood simulation suite. The simulations incorporate stellar feedback, AGN feedback and a simplified scheme for efficiently modelling the low column density Ly$α$ forest. We confirm a discrepancy remains between Cos… ▽ More We study the effect of different feedback prescriptions on the properties of the low redshift ($z\leq1.6$) Ly$α$ forest using a selection of hydrodynamical simulations drawn from the Sherwood simulation suite. The simulations incorporate stellar feedback, AGN feedback and a simplified scheme for efficiently modelling the low column density Ly$α$ forest. We confirm a discrepancy remains between Cosmic Origins Spectrograph (COS) observations of the Ly$α$ forest column density distribution function (CDDF) at $z \simeq 0.1$ for high column density systems ($N_{\rm HI}>10^{14}\rm\,cm^{-2}$), as well as Ly$α$ velocity widths that are too narrow compared to the COS data. Stellar or AGN feedback -- as currently implemented in our simulations -- have only a small effect on the CDDF and velocity width distribution. We conclude that resolving the discrepancy between the COS data and simulations requires an increase in the temperature of overdense gas with $Δ=4$--$40$, either through additional He$\,\rm \scriptstyle II\ $ photo-heating at $z>2$ or fine-tuned feedback that ejects overdense gas into the IGM at just the right temperature for it to still contribute significantly to the Ly$α$ forest. Alternatively a larger, currently unresolved turbulent component to the line width could resolve the discrepancy. △ Less

Submitted 28 June, 2017; v1 submitted 15 June, 2017; originally announced June 2017.

Comments: 16 pages, 8 figures, Accepted to MNRAS

AGN jet feedback on a moving mesh: cocoon inflation, gas flows and turbulence

Authors: Martin A. Bourne, Debora Sijacki

Abstract: In many observed galaxy clusters, jets launched by the accretion process on to supermassive black holes, inflate large-scale cavities filled with energetic, relativistic plasma. This process is thought to be responsible for regulating cooling losses, thus moderating the inflow of gas on to the central galaxy, quenching further star formation and maintaining the galaxy in a red and dead state. In t… ▽ More In many observed galaxy clusters, jets launched by the accretion process on to supermassive black holes, inflate large-scale cavities filled with energetic, relativistic plasma. This process is thought to be responsible for regulating cooling losses, thus moderating the inflow of gas on to the central galaxy, quenching further star formation and maintaining the galaxy in a red and dead state. In this paper, we implement a new jet feedback scheme into the moving mesh-code arepo, contrast different jet injection techniques and demonstrate the validity of our implementation by comparing against simple analytical models. We find that jets can significantly affect the intracluster medium (ICM), offset the overcooling through a number of heating mechanisms, as well as drive turbulence, albeit within the jet lobes only. Jet-driven turbulence is, however, a largely ineffective heating source and is unlikely to dominate the ICM heating budget even if the jet lobes efficiently fill the cooling region, as it contains at most only a few per cent of the total injected energy. We instead show that the ICM gas motions, generated by orbiting substructures, while inefficient at heating the ICM, drive large-scale turbulence and when combined with jet feedback, result in line-of-sight velocities and velocity dispersions consistent with the Hitomi observations of the Perseus cluster. △ Less

Submitted 1 November, 2017; v1 submitted 22 May, 2017; originally announced May 2017.

Comments: 31 pages, 19 figures, MNRAS accepted

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 472, Issue 4, 21 December 2017, Pages 4707-4735

Driving gas shells with radiation pressure on dust in radiation-hydrodynamic simulations

Authors: Tiago Costa, Joakim Rosdahl, Debora Sijacki, Martin Haehnelt

Abstract: We present radiation-hydrodynamic simulations of radiatively-driven gas shells launched by bright active galactic nuclei (AGN) in isolated dark matter haloes. Our goals are (1) to investigate the ability of AGN radiation pressure on dust to launch galactic outflows and (2) to constrain the efficiency of infrared (IR) multi-scattering in boosting outflow acceleration. Our simulations are performed… ▽ More We present radiation-hydrodynamic simulations of radiatively-driven gas shells launched by bright active galactic nuclei (AGN) in isolated dark matter haloes. Our goals are (1) to investigate the ability of AGN radiation pressure on dust to launch galactic outflows and (2) to constrain the efficiency of infrared (IR) multi-scattering in boosting outflow acceleration. Our simulations are performed with the radiation-hydrodynamic code RAMSES-RT and include both single- and multi-scattered radiation pressure from an AGN, radiative cooling and self-gravity. Since outflowing shells always eventually become transparent to the incident radiation field, outflows that sweep up all intervening gas are likely to remain gravitationally bound to their halo even at high AGN luminosities. The expansion of outflowing shells is well described by simple analytic models as long as the shells are mildly optically thick to IR radiation. In this case, an enhancement in the acceleration of shells through IR multi-scattering occurs as predicted, i.e. a force dP/dt = tau_IR L/c is exerted on the gas. For high optical depths tau_IR > 50, however, momentum transfer between outflowing optically thick gas and IR radiation is rapidly suppressed, even if the radiation is efficiently confined. At high tau_IR, the characteristic flow time becomes shorter than the required trapping time of IR radiation such that the momentum flux dP/dt << tau_IR L/c. We argue that while unlikely to unbind massive galactic gaseous haloes, AGN radiation pressure on dust could play an important role in regulating star formation and black hole accretion in the nuclei of massive compact galaxies at high redshift. △ Less

Submitted 24 May, 2018; v1 submitted 16 March, 2017; originally announced March 2017.

Comments: 25 pages, 12 figures, 3 tables (accepted by MNRAS)

The Origin and Evolution of Fast and Slow Rotators in the Illustris Simulation

Authors: Zephyr Penoyre, Benjamin P. Moster, Debora Sijacki, Shy Genel

Abstract: Using the Illustris simulation, we follow thousands of elliptical galaxies back in time to identify how the dichotomy between fast and slow rotating ellipticals (FRs and SRs) develops. Comparing to the $\textrm{ATLAS}^\textrm{3D}$ survey, we show that Illustris reproduces similar elliptical galaxy rotation properties, quantified by the degree of ordered rotation, $λ_\textrm{R}$. There is a clear s… ▽ More Using the Illustris simulation, we follow thousands of elliptical galaxies back in time to identify how the dichotomy between fast and slow rotating ellipticals (FRs and SRs) develops. Comparing to the $\textrm{ATLAS}^\textrm{3D}$ survey, we show that Illustris reproduces similar elliptical galaxy rotation properties, quantified by the degree of ordered rotation, $λ_\textrm{R}$. There is a clear segregation between low-mass ($M_{\rm *} < 10^{11} M_{\rm \odot}$) ellipticals, which form a smooth distribution of FRs, and high-mass galaxies ($M_{\rm *} > 10^{11.5} M_{\rm \odot}$), which are mostly SRs, in agreement with observations. We find that SRs are very gas poor, metal rich and red in colour, while FRs are generally more gas rich and still star forming. We suggest that ellipticals begin naturally as FRs and, as they grow in mass, lose their spin and become SRs. While at $z = 1$, the progenitors of SRs and FRs are nearly indistinguishable, their merger and star formation histories differ thereafter. We find that major mergers tend to disrupt galaxy spin, though in rare cases can lead to a spin-up. No major difference is found between the effects of gas-rich and gas-poor mergers and the amount of minor mergers seem to have little correlation with galaxy spin. In between major mergers, lower-mass ellipticals, which are mostly gas-rich, tend to recover their spin by accreting gas and stars. For galaxies with $M_{\rm *}$ above $\sim 10^{11} M_{\rm \odot}$, this trend reverses; galaxies only retain or steadily lose their spin. More frequent mergers, accompanied by an inability to regain spin, lead massive ellipticals to lose most of ordered rotation and transition from FRs to SRs. △ Less

Submitted 28 March, 2017; v1 submitted 1 March, 2017; originally announced March 2017.

Comments: 26 pages, 21 figures, MNRAS accepted

Constraints on AGN feedback from its Sunyaev-Zel'dovich imprint on the cosmic background radiation

Authors: Bjoern Soergel, Tommaso Giannantonio, George Efstathiou, Ewald Puchwein, Debora Sijacki

Abstract: We derive constraints on feedback by active galactic nuclei (AGN) by setting limits on their thermal Sunyaev-Zel'dovich (SZ) imprint on the cosmic microwave background (CMB). The amplitude of any SZ signature is small and degenerate with the poorly known sub-mm spectral energy distribution of the AGN host galaxy and other unresolved dusty sources along the line of sight. Here we break this degener… ▽ More We derive constraints on feedback by active galactic nuclei (AGN) by setting limits on their thermal Sunyaev-Zel'dovich (SZ) imprint on the cosmic microwave background (CMB). The amplitude of any SZ signature is small and degenerate with the poorly known sub-mm spectral energy distribution of the AGN host galaxy and other unresolved dusty sources along the line of sight. Here we break this degeneracy by combining microwave and sub-mm data from Planck with all-sky far-infrared maps from the AKARI satellite. We first test our measurement pipeline using the Sloan Digital Sky Survey (SDSS) redMaPPer catalogue of galaxy clusters, finding a highly significant detection ($>$$20σ$) of the SZ effect together with correlated dust emission. We then constrain the SZ signal associated with spectroscopically confirmed quasi-stellar objects (QSOs) from SDSS data release 7 (DR7) and the Baryon Oscillation Spectroscopic Survey (BOSS) DR12. We obtain a low-significance ($1.6σ$) hint of an SZ signal, pointing towards a mean thermal energy of $\simeq 5 \times 10^{60}$ erg, lower than reported in some previous studies. A comparison of our results with high-resolution hydrodynamical simulations including AGN feedback suggests QSO host masses of $M_{200c} \sim 4 \times 10^{12}~h^{-1}M_\odot$, but with a large uncertainty. Our analysis provides no conclusive evidence for an SZ signal specifically associated with AGN feedback. △ Less

Submitted 3 May, 2017; v1 submitted 19 December, 2016; originally announced December 2016.

Comments: 21 pages, 14 figures, matches version published by MNRAS

Journal ref: MNRAS (2017) 468 (1): 577-596

Interpreting ALMA Observations of the ISM During the Epoch of Reionisation

Authors: Harley Katz, Taysun Kimm, Debora Sijacki, Martin Haehnelt

Abstract: We present cosmological, radiation-hydrodynamics simulations of galaxy formation during the epoch of reionization in an effort towards modelling the interstellar medium (ISM) and interpreting ALMA observations. Simulations with and without stellar radiation are compared at large (Mpc), intermediate (tens of kpc), and small (sub kpc) scales. At large scales, the dense regions around galaxies reioni… ▽ More We present cosmological, radiation-hydrodynamics simulations of galaxy formation during the epoch of reionization in an effort towards modelling the interstellar medium (ISM) and interpreting ALMA observations. Simulations with and without stellar radiation are compared at large (Mpc), intermediate (tens of kpc), and small (sub kpc) scales. At large scales, the dense regions around galaxies reionize first before UV photons penetrate the voids; however, considerable amounts of neutral gas remain present within the haloes. The spatial distribution of neutral gas is highly dynamic and is anti-correlated with the presence of stars older than a few Myrs. For our specific feedback implementation, most of the metals remain inside the virial radii of haloes and they are proportionally distributed over the ionized and neutral medium by mass. For our most massive galaxy with ${\rm M_h}\sim10^{11}$M$_{\odot}$, the majority of the CII and OI mass are associated with cold neutral clumps. NII is more diffuse and arises in warmer gas while OIII arises in hotter gas with a higher ionization parameter, produced by photo-heating and supernovae. If smaller pockets of high metallicity gas exist in the ISM, the emission from these ions may be observable by ALMA while the low metallicity of the galaxy may cause these systems to fall below the local [CII]-SFR relation. The presence of dust can cause spatial offsets between UV/Ly$α$ and [CII] emission as suggested by the recent observations of Maiolino et al. [OIII] may be spatially offset from both of these components since it arises from a different part of density-temperature phase-space. △ Less

Submitted 21 March, 2017; v1 submitted 6 December, 2016; originally announced December 2016.

Comments: 34 Pages, 26 figures, MNRAS Accepted

Black hole clustering and duty cycles in the Illustris simulation

Authors: Colin DeGraf, Debora Sijacki

Abstract: We use the high-resolution cosmological simulation Illustris to investigate the clustering of supermassive black holes across cosmic time, the link between black hole clustering and host halo masses, and the implications for black hole duty cycles. Our predicted black hole correlation length and bias match the observational data very well across the full redshift range probed. Black hole clusterin… ▽ More We use the high-resolution cosmological simulation Illustris to investigate the clustering of supermassive black holes across cosmic time, the link between black hole clustering and host halo masses, and the implications for black hole duty cycles. Our predicted black hole correlation length and bias match the observational data very well across the full redshift range probed. Black hole clustering is strongly luminosity-dependent on small, 1-halo scales, with some moderate dependence on larger scales of a few Mpc at intermediate redshifts. We find black hole clustering to evolve only weakly with redshift, initially following the behaviour of their hosts. However below z ~ 2 black hole clustering increases faster than that of their hosts, which leads to a significant overestimate of the clustering-predicted host halo mass. The full distribution of host halo masses is very wide, including a low-mass tail extending up to an order of magnitude below the naive prediction for minimum host mass. Our black hole duty cycles follow a power-law dependence on black hole mass and decrease with redshift, and we provide accurate analytic fits to these. The increase in clustering amplitude at late times, however, means that duty cycle estimates based on black hole clustering can overestimate duty cycles substantially, by more than two orders of magnitude. We find the best agreement when the minimum host mass is assumed to be $10^{11.2} M_\odot$ , which provides an accurate measure across all redshifts and luminosity ranges probed by our simulation. △ Less

Submitted 20 January, 2017; v1 submitted 21 September, 2016; originally announced September 2016.

Comments: 15 pages, 14 figures; MNRAS accepted

Resolving flows around black holes: the impact of gas angular momentum

Authors: Michael Curtis, Debora Sijacki

Abstract: Cosmological simulations almost invariably estimate the accretion of gas on to supermassive black holes using a Bondi-Hoyle-like prescription. Doing so ignores the effects of the angular momentum of the gas, which may prevent or significantly delay accreting material falling directly on to the black hole. We outline a black hole accretion rate prescription using a modified Bondi-Hoyle formulation… ▽ More Cosmological simulations almost invariably estimate the accretion of gas on to supermassive black holes using a Bondi-Hoyle-like prescription. Doing so ignores the effects of the angular momentum of the gas, which may prevent or significantly delay accreting material falling directly on to the black hole. We outline a black hole accretion rate prescription using a modified Bondi-Hoyle formulation that takes into account the angular momentum of the surrounding gas. Meaningful implementation of this modified Bondi-Hoyle formulation is only possible when the inner vorticity distribution is well resolved, which we achieve through the use of a super-Lagrangian refinement technique around black holes within our simulations. We then investigate the effects on black hole growth by performing simulations of isolated as well as merging disc galaxies using the moving-mesh code AREPO. We find that the gas angular momentum barrier can play an important role in limiting the growth of black holes, leading also to a several Gyr delay between the starburst and the quasar phase in major merger remnants. We stress, however, that the magnitude of this effect is highly sensitive to the thermodynamical state of the accreting gas and to the nature of the black hole feedback present. △ Less

Submitted 10 October, 2016; v1 submitted 8 June, 2016; originally announced June 2016.

Comments: 16 pages, 8 figures, MNRAS accepted

Journal ref: MNRAS (November 21, 2016) 463 (1): 63-77

The Sherwood simulation suite: overview and data comparisons with the Lyman-alpha forest at redshifts 2 < z < 5

Authors: James S. Bolton, Ewald Puchwein, Debora Sijacki, Martin G. Haehnelt, Tae-Sun Kim, Avery Meiksin, John A. Regan, Matteo Viel

Abstract: We introduce a new set of large scale, high resolution hydrodynamical simulations of the intergalactic medium: the Sherwood simulation suite. These are performed in volumes 10^3-160^3 h^-3 comoving Mpc^{3}, span almost four orders of magnitude in mass resolution with up to 17.2 billion particles, and employ a variety of physics variations including warm dark matter and galactic outflows. We undert… ▽ More We introduce a new set of large scale, high resolution hydrodynamical simulations of the intergalactic medium: the Sherwood simulation suite. These are performed in volumes 10^3-160^3 h^-3 comoving Mpc^{3}, span almost four orders of magnitude in mass resolution with up to 17.2 billion particles, and employ a variety of physics variations including warm dark matter and galactic outflows. We undertake a detailed comparison of the simulations to high resolution, high signal-to-noise observations of the Lyman-alpha forest over the redshift range 2 < z < 5. The simulations are in very good agreement with the observational data, lending further support to the paradigm that the Lyman-alpha forest is a natural consequence of the web-like distribution of matter arising in LCDM cosmological models. Only a small number of minor discrepancies remain with respect to the observational data. Saturated Lyman-alpha absorption lines with column densities N_HI > 10^14.5 cm^-2 at 2 < z < 2.5 are underpredicted in the models. An uncertain correction for continuum placement bias is required to match the distribution and power spectrum of the transmitted flux, particularly at z > 4. Finally, the temperature of intergalactic gas in the simulations may be slightly too low at z=2.7 and a flatter temperature-density relation is required at z=2.4, consistent with the expected effects of non-equilibrium ionisation during He-II reionisation. △ Less

Submitted 19 September, 2016; v1 submitted 11 May, 2016; originally announced May 2016.

Comments: 19 pages, 21 figures, accepted to MNRAS. Images and movies may be viewed at http://www.nottingham.ac.uk/astronomy/sherwood/

Shock finding on a moving-mesh: II. Hydrodynamic shocks in the Illustris universe

Authors: Kevin Schaal, Volker Springel, Rüdiger Pakmor, Christoph Pfrommer, Dylan Nelson, Mark Vogelsberger, Shy Genel, Annalisa Pillepich, Debora Sijacki, Lars Hernquist

Abstract: Hydrodynamical shocks are a manifestation of the non-linearity of the Euler equations and play a fundamental role in cosmological gas dynamics. In this work, we identify and analyse shocks in the Illustris simulation, and contrast the results with those of non-radiative runs. We show that simulations with more comprehensive physical models of galaxy formation pose new challenges for shock finding… ▽ More Hydrodynamical shocks are a manifestation of the non-linearity of the Euler equations and play a fundamental role in cosmological gas dynamics. In this work, we identify and analyse shocks in the Illustris simulation, and contrast the results with those of non-radiative runs. We show that simulations with more comprehensive physical models of galaxy formation pose new challenges for shock finding algorithms due to radiative cooling and star-forming processes, prompting us to develop a number of methodology improvements. We find in Illustris a total shock surface area which is about 1.4 times larger at the present epoch compared to non-radiative runs, and an energy dissipation rate at shocks which is higher by a factor of around 7. Remarkably, shocks with Mach numbers above and below $\mathcal{M}\approx10$ contribute about equally to the total dissipation across cosmic time. This is in sharp contrast to non-radiative simulations, and we demonstrate that a large part of the difference arises due to strong black hole radio-mode feedback in Illustris. We also provide an overview of the large diversity of shock morphologies, which includes complex networks of halo-internal shocks, shocks on to cosmic sheets, feedback shocks due to black holes and galactic winds, as well as ubiquitous accretion shocks. In high redshift systems more massive than $10^{12}\,\mathrm{M}_\odot$ we discover the existence of a double accretion shock pattern in haloes. They are created when gas streams along filaments without being shocked at the outer accretion shock, but then forms a second, roughly spherical accretion shock further inside. △ Less

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

Comments: 26 pages, 15 figures, 1 movie, published in mnras

Journal ref: MRNAS 461, 4441-4465 (2016)

Powerful quasar outflow in a massive disc galaxy at $z \sim 5$

Authors: Michael Curtis, Debora Sijacki

Abstract: There is growing observational evidence of high-redshift quasars launching energetic, fast outflows, but the effects that these have on their host galaxies is poorly understood. We employ the moving-mesh code AREPO to study the feedback from a quasar that has grown to $\sim 10^9 M_\odot$ by $z \sim 5$ and the impact that this has on its host galaxy. Our simulations use a super-Lagrangian refinemen… ▽ More There is growing observational evidence of high-redshift quasars launching energetic, fast outflows, but the effects that these have on their host galaxies is poorly understood. We employ the moving-mesh code AREPO to study the feedback from a quasar that has grown to $\sim 10^9 M_\odot$ by $z \sim 5$ and the impact that this has on its host galaxy. Our simulations use a super-Lagrangian refinement technique to increase the accuracy with which the interface of the quasar-driven wind and the surrounding gas is resolved. We find that the feedback injected in these simulations is less efficient at removing gas from the galaxy than in an identical simulation with no super-Lagrangian refinement. This leads to the growth of a massive, rotationally supported, star-forming disc, co-existing with a powerful quasar-driven outflow. The properties of our host galaxy, including the kinematical structure of the gaseous disc and of the outflow, are in good agreement with current observations. Upcoming ALMA and JWST observations will be an excellent test of our model and will provide further clues as to the variance in properties of high-redshift quasar hosts. △ Less

Submitted 21 January, 2016; v1 submitted 21 October, 2015; originally announced October 2015.

Comments: 6 pages, 4 figures, MNRAS submitted, figures and videos are available at http://www.ast.cam.ac.uk/~mc636/

Journal ref: Monthly Notices of the Royal Astronomical Society: Letters 2016 457 (1): L34-L38

Recoiling black holes: prospects for detection and implications of spin alignment

Authors: Laura Blecha, Debora Sijacki, Luke Zoltan Kelley, Paul Torrey, Mark Vogelsberger, Dylan Nelson, Volker Springel, Gregory Snyder, Lars Hernquist

Abstract: Supermassive black hole (BH) mergers produce powerful gravitational wave (GW) emission. Asymmetry in this emission imparts a recoil kick to the merged BH, which can eject the BH from its host galaxy altogether. Recoiling BHs could be observed as offset active galactic nuclei (AGN). Several candidates have been identified, but systematic searches have been hampered by large uncertainties regarding… ▽ More Supermassive black hole (BH) mergers produce powerful gravitational wave (GW) emission. Asymmetry in this emission imparts a recoil kick to the merged BH, which can eject the BH from its host galaxy altogether. Recoiling BHs could be observed as offset active galactic nuclei (AGN). Several candidates have been identified, but systematic searches have been hampered by large uncertainties regarding their observability. By extracting merging BHs and host galaxy properties from the Illustris cosmological simulations, we have developed a comprehensive model for recoiling AGN. Here, for the first time, we model the effects of BH spin alignment and recoil dynamics based on the gas-richness of host galaxies. We predict that if BH spins are not highly aligned, seeing-limited observations could resolve offset AGN, making them promising targets for all-sky surveys. For randomly-oriented spins, less than about 10 spatially-offset AGN may be detectable in HST-COSMOS, and > 10^3 could be found with Pan-STARRS, LSST, Euclid, and WFIRST. Nearly a thousand velocity-offset AGN are predicted within the SDSS footprint; the rarity of large broad-line offsets among SDSS quasars is likely due in part to selection effects but suggests that spin alignment plays a role in suppressing recoils. Nonetheless, in our most physically motivated model where alignment occurs only in gas-rich mergers, hundreds of offset AGN should be found in all-sky surveys. Our findings strongly motivate a dedicated search for recoiling AGN. △ Less

Submitted 10 December, 2015; v1 submitted 6 August, 2015; originally announced August 2015.

Comments: 30 pages, 19 figures. Accepted to MNRAS after minor revisions

Stellar Spiral Structures in Triaxial Dark Matter Haloes

Authors: Shaoran Hu, Debora Sijacki

Abstract: We employ very high resolution simulations of isolated Milky Way-like galaxies to study the effect of triaxial dark matter haloes on exponential stellar discs. Non-adiabatic halo shape changes can trigger two-armed grand-design spiral structures which extend all the way to the edge of the disc. Their pattern speed coincides with the inner Lindblad resonance indicating that they are kinematic densi… ▽ More We employ very high resolution simulations of isolated Milky Way-like galaxies to study the effect of triaxial dark matter haloes on exponential stellar discs. Non-adiabatic halo shape changes can trigger two-armed grand-design spiral structures which extend all the way to the edge of the disc. Their pattern speed coincides with the inner Lindblad resonance indicating that they are kinematic density waves which can persist up to several Gyr. In dynamically cold discs, grand-design spirals are swing amplified and after a few Gyr can lead to the formation of (multi-armed) transient recurrent spirals. Stellar discs misaligned to the principal planes of the host triaxial halo develop characteristic integral shaped warps, but otherwise exhibit very similar spiral structures as aligned discs. For the grand-design spirals in our simulations, their strength dependence with radius is determined by the torque on the disc, suggesting that by studying grand-design spirals without bars it may be possible to set constraints on the tidal field and host dark matter halo shape. △ Less

Submitted 4 July, 2016; v1 submitted 6 July, 2015; originally announced July 2015.

Comments: 21 pages, 20 figures. MNRAS accepted