Source author record

Alexander Hobbs

Alexander Hobbs appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

12works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2020arXiv

Positive feedback at the disc-halo interface

The flat star formation (SF) history of the Milky Way (MW) requires gas in the Galactic disc to be replenished, most likely from a reservoir outside the Galaxy. Such a replenishment may be achieved by a form of `positive' feedback, whereby SF feedback creates a Galactic fountain cycle that collects and cools additional gas from the hot halo surrounding the Galaxy. In this paper we present a model of this process for the MW. A section of the Galactic disc is allowed to form stars which subsequently explode as supernovae (SNe) and send gas out into the hot halo. The gas that is sent out is colder than the hot halo gas and, as it mixes, the halo gas is cooled, providing fuel for further SF as the mixture falls back onto the Galactic disc. We find that this process can be sufficient to maintain a roughly-constant cold gas mass in the MW over at least 3 Gyr. Our results further suggest that there is a positive feedback trend whereby increasing SF leads to an increase in the cold gas budget at average SF rates below 0.5 solar masses per year but above which becomes negative, where further increasing the SFR causes the cold gas budget to decrease. We have constructed an analytical model for this that reproduces the data well and could have profound implications for galaxy evolution in feedback-dominated regimes.

preprint2016arXiv

Spurious Small-Scale Structure & Discreteness-Driven Relaxation in Cosmological Simulations

There is strong evidence that cosmological N-body simulations dominated by Warm Dark Matter (WDM) contain spurious or unphysical haloes, most readily apparent as regularly spaced low-mass haloes strung along filaments. We show that spurious haloes are a feature of traditional N-body simulations of cosmological structure formation models, including WDM and Cold Dark Matter (CDM) models, in which gravitational collapse proceeds in an initially anisotropic fashion, and arises naturally as a consequence of discreteness-driven relaxation. We demonstrate this using controlled N-body simulations of plane-symmetric collapse and show that spurious haloes are seeded at shell crossing by localised velocity perturbations induced by the discrete nature of the density field, and that their characteristic separation should be approximately the mean inter-particle separation of the N-body simulation, which is fixed by the mass resolution within the volume. Using cosmological N-body simulations in which particles are split into two collisionless components of fixed mass ratio, we find that the spatial distribution of the two components show signatures of discreteness-driven relaxation in their spatial distribution on both large and small scales. Adopting a spline kernel gravitational softening that is of order the comoving mean inter-particle separation helps to suppress the effect of discreteness-driven relaxation, but cannot eliminate it completely. These results provide further motivation for recent developments of new algorithms, which include, for example, revisions of the traditional N-body approach by means of spatially adaptive anistropric gravitational softenings or explicit solutions for the evolution of dark matter in phase space.

preprint2015arXiv

Novel Adaptive softening for collisionless N-body simulations: Eliminating spurious halos

We describe a NOVel form of Adaptive softening (NovA) for collisionless $N$-body simulations, implemented in the Ramses adaptive mesh refinement code. We introduce a refinement criterion that the particle distribution within each cell be sufficiently isotropic, as measured by its moment of inertia tensor. In this way, collapse is only refined if it occurs along all three axes, ensuring that the softening $ε$ is always of order twice the largest inter-particle spacing in a cell. This more conservative force softening criterion is designed to minimise spurious two-body effects, while maintaining high force resolution in collapsed regions of the flow. We test NovA using an antisymmetric perturbed plane wave collapse (`Valinia' test) before applying it to warm dark matter (WDM) simulations. For the Valinia test, we show that -- unlike the standard $N$-body method -- NovA produces no numerical fragmentation while still being able to correctly capture fine caustics and shells around the collapsing regions. For the WDM simulations, we find that NovA converges significantly more rapidly than standard $N$-body, producing little or no spurious halos on small scales. We show, however, that determining whether or not halos exist below the free streaming mass $M_{\rm fs}$ is complicated by the fact that our halo finder (AHF) likely incorrectly labels some caustics and criss-crossing filaments as halos, while one or two particularly massive filaments appear to fragment in any version of NovA where refinement is allowed. Such massive filaments may be physically unstable to collapse, as is the case for infinite, static, self-gravitating cylinders. We will use NovA in forthcoming papers to study the issue of halo formation below $M_{\rm fs}$; filament stability; and to obtain new constraints on the temperature of dark matter.

preprint2014arXiv

AGN Feedback models: Correlations with star formation and observational implications of time evolution

We examine the correlation between the star formation rate (SFR) and black hole accretion rate (BHAR) across a suite of different AGN feedback models, using the time evolution of a merger simulation. By considering three different stages of evolution, and a distinction between the nuclear and outer regions of star formation, we consider 63 different cases. Despite many of the feedback models fitting the M-σ relationship well, there are often distinct differences in the SFR-BHAR correlations, with close to linear trends only being present after the merger. Some of the models also show evolution in the SFR-BHAR parameter space that is at times directly across the long-term averaged SFR-BHAR correlation. This suggests that the observational SFR-BHAR correlation found for ensembles of galaxies is an approximate statistical trend, as suggested by Hickox et al. Decomposing the SFR into nuclear and outer components also highlights notable differences between models and there is only modest agreement with observational studies examining this in Seyfert galaxies. For the fraction of the black hole mass growth from the merger event relative to the final black hole mass, we find as much as a factor of three variation among models. This also translates into a similar variation in the post-starburst black hole mass growth. Overall, we find that while qualitative features are often similar amongst models, precise quantitative analysis shows there can be quite distinct differences.

preprint2014arXiv

Black hole feedback in a multiphase interstellar medium

Ultrafast outflows (UFOs) from supermassive black holes (SMBHs) are thought to regulate the growth of SMBHs and host galaxies, resulting in a number of observational correlations. We present high-resolution numerical simulations of the impact of a thermalized UFO on the ambient gas in the inner part of the host galaxy. Our results depend strongly on whether the gas is homogeneous or clumpy. In the former case all of the ambient gas is driven outward rapidly as expected based on commonly used energy budget arguments, while in the latter the flows of mass and energy decouple. Carrying most of the energy, the shocked UFO escapes from the bulge via paths of least resistance, taking with it only the low-density phase of the host. Most of the mass is however in the high-density phase, and is affected by the UFO much less strongly, and may even continue to flow inwards. We suggest that the UFO energy leakage through the pores in the multiphase interstellar medium (ISM) may explain why observed SMBHs are so massive despite their overwhelmingly large energy production rates. The multiphase ISM effects reported here are probably under-resolved in cosmological simulations but may be included in prescriptions for active galactic nuclei feedback in future simulations and in semi-analytical models.

preprint2014arXiv

Growing galaxies via superbubble-driven accretion flows

We use a suite a cooling halo simulations to study a new mechanism for rapid accretion of hot halo gas onto star-forming galaxies. Correlated supernovae events create converging 'superbubbles' in the halo gas. Where these collide, the density increases, driving cooling filaments of low metallicity gas that feed the disc. At our current numerical resolution (20 pc) we are only able to resolve the most dramatic events; these could be responsible for the build-up of galaxy discs after the most massive gas-rich mergers have completed (z < 1). As we increase the numerical resolution, we find that the filaments persist for longer, driving continued late-time star formation. This suggests that SNe-driven accretion could act as an efficient mechanism for extracting cold gas from the hot halo, driving late-time star formation in disc galaxies. We show that such filament feeding leads to a peak star formation rate (SFR) of $\sim 3$ M$_{\rm sun}$ yr$^{-1}$, consistent with estimates for the Milky Way. By contrast, direct cooling from the hot halo ('hot-mode' accretion, not present in the simulations that show filament feeding) falls short of the SNe-driven SFR by a factor of 3-4, and is sustained over a shorter time period. The filaments we resolve extend to $\sim$ 50 kpc, reaching column densities of $\sim 10^{18}$ cm$^{-2}$. We show that such structures can plausibly explain the broad dispersion in Mg II absorption seen along sight lines to quasars. Our results suggest a dual role for stellar feedback in galaxy formation, suppressing hot-mode accretion while promoting cold-mode accretion along filaments. This ultimately leads to more star formation, suggesting that the positive feedback effect outweighs the negative. Finally, since the filamentary gas has higher angular momentum than that coming from hot-mode accretion, we show that this leads to the formation of substantially larger gas discs.

preprint2012arXiv

Growing supermassive black holes: sub-grid modelling and intermediate-scale processes

The sheer range of scales in the Universe makes it impossible to model all at once. It is necessary, therefore, when conducting numerical experiments, that we employ sub-resolution prescriptions that can represent the scales we are unable to model directly. In this article we present a prescription for black hole growth that incorporates a different accretion regime from the standard approach used in the literature, and discuss the results of dedicated simulations of intermediate processes between small-scale accretion flows and large-scale cosmological volumes that can strongly enhance the accretion rate onto the black hole at the centre of a galaxy.

preprint2012arXiv

Modelling supermassive black hole growth: towards an improved sub-grid prescription

Accretion onto supermassive black holes (SMBHs) in galaxy formation simulations is frequently modelled by the Bondi-Hoyle formalism. Here we examine the validity of this approach analytically and numerically. We argue that the character of the flow where one evaluates the gas properties is unlikely to satisfy the simple Bondi-Hoyle model. Only in the specific case of hot virialised gas with zero angular momentum and negligible radiative cooling is the Bondi-Hoyle solution relevant. In the opposite extreme, where the gas is in a state of free-fall at the evaluation radius due to efficient cooling and the dominant gravity of the surrounding halo, the Bondi-Hoyle formalism can be erroneous by orders of magnitude in either direction. This may impose artificial trends with halo mass in cosmological simulations by being wrong by different factors for different halo masses. We propose an expression for the sub-grid accretion rate which interpolates between the free-fall regime and the Bondi-Hoyle regime, therefore taking account of the contribution of the halo to the gas dynamics.

preprint2012arXiv

Thermal instabilities in cooling galactic coronae: fuelling star formation in galactic discs

We investigate the means by which cold gas can accrete onto Milky Way mass galaxies from a hot corona of gas, using a new smoothed particle hydrodynamics code, 'SPHS'. We find that the 'cold clumps' seen in many classic SPH simulations in the literature are not present in our SPHS simulations. Instead, cold gas condenses from the halo along filaments that form at the intersection of supernovae-driven bubbles from previous phases of star formation. This positive feedback feeds cold gas to the galactic disc directly, fuelling further star formation. The resulting galaxies in the SPH and SPHS simulations differ greatly in their morphology, gas phase diagrams, and stellar content. We show that the classic SPH cold clumps owe to a numerical thermal instability caused by an inability for cold gas to mix in the hot halo. The improved treatment of mixing in SPHS suppresses this instability leading to a dramatically different physical outcome. In our highest resolution SPHS simulation, we find that the cold filaments break up into bound clumps that form stars. The filaments are overdense by a factor of 10-100 compared to the surrounding gas, suggesting that the fragmentation results from a physical non-linear instability driven by the overdensity. This 'fragmenting filament' mode of disc growth has important implications for galaxy formation, in particular the role of star formation in bringing cold gas into disc galaxies.

preprint2011arXiv

Feeding SMBHs through supersonic turbulence and ballistic accretion

It has long been recognised that the main obstacle to accretion of gas onto supermassive black holes (SMBHs) is large specific angular momentum. It is feared that the gas settles in a large scale disc, and that accretion would then proceed too inefficiently to explain the masses of the observed SMBHs. Here we point out that, while the mean angular momentum in the bulge is very likely to be large, the deviations from the mean can also be significant. Indeed, cosmological simulations show that velocity and angular momentum fields of gas flows onto galaxies are very complex. Furthermore, inside bulges the gas velocity distribution can be further randomised by the velocity kicks due to feedback from star formation. We perform hydrodynamical simulations of gaseous rotating shells infalling onto an SMBH, attempting to quantify the importance of velocity dispersion in the gas at relatively large distances from the black hole. We implement this dispersion by means of a supersonic turbulent velocity spectrum. We find that, while in the purely rotating case the circularisation process leads to efficient mixing of gas with different angular momentum, resulting in a low accretion rate, the inclusion of turbulence increases this accretion rate by up to several orders of magnitude. We show that this can be understood based on the notion of "ballistic" accretion, whereby dense filaments, created by convergent turbulent flows, travel through the ambient gas largely unaffected by hydrodynamical drag. We derive a simple analytical formula that captures the numerical results to within a factor of a few. Rescaling our results to astrophysical bulges, we argue that this "ballistic" mode of accretion could provide the SMBHs with a sufficient supply of fuel without the need to channel the gas via large-scale discs or bars, and therefore that star formation in bulges can be a strong catalyst for SMBH accretion.

preprint2009arXiv

Dynamic Monte Carlo radiation transfer in SPH. Radiation pressure force implementation

We present a new framework for radiation hydrodynamics simulations. Gas dynamics is modelled by the Smoothed Particle Hydrodynamics (SPH) method, whereas radiation transfer is simulated via a time-dependent Monte-Carlo approach that traces photon packets. As a first step in the development of the method, in this paper we consider the momentum transfer between radiation field and gas, which is important for systems where radiation pressure is high. There is no fundamental limitations on the number of radiation sources, geometry or the optical depth of the problems that can be studied with the method. However, as expected for any Monte-Carlo transfer scheme, stochastic noise presents a serious limitation. We present a number of tests that show that the errors of the method can be estimated accurately by considering Poisson noise fluctuations in the number of photon packets that SPH particles interact with per dynamical time. It is found that for a reasonable accuracy the momentum carried by photon packets must be much smaller than a typical momentum of SPH particles. We discuss numerical limitations of the code, and future steps that can be taken to improve performance and applicability of the method.

preprint2008arXiv

Simulations of the formation of a gaseous disc and young stars near Sgr A* via cloud-cloud collisions

Young massive stars in the central parsec of our Galaxy are best explained by star formation within at least one, and possibly two, massive self-gravitating gaseous discs. With help of numerical simulations, we here consider whether the observed population of young stars could have originated from a large angle collision of two massive gaseous clouds at R ~ 1 pc from Sgr A*. In all the simulations performed, the post-collision gas flow forms an inner nearly circular gaseous disc and one or two eccentric outer filaments, consistent with the observations. Furthermore, the radial stellar mass distribution is always very steep, sigma proportional to R^-2, again consistent with the data. The 3D velocity structure of the stellar distribution is however sensitive to initial conditions (e.g., the impact parameter of the clouds) and gas cooling details. In all the cases the amount of gas accreted by our inner boundary condition is large, enough to allow Sgr A* to radiate near its Eddington limit during ~ 10^5 years. This suggests that a refined model would have physically larger clouds (or a cloud and a disc such as the CND) colliding at a distance of a few pc rather than 1 pc as in our simulations.