Tiziana Di Matteo

Bio: Tiziana Di Matteo is an academic researcher from Carnegie Mellon University. The author has contributed to research in topics: Galaxy & Supermassive black hole. The author has an hindex of 58, co-authored 233 publications receiving 25089 citations. Previous affiliations of Tiziana Di Matteo include Australian National University & University of Oxford.

Energy input from quasars regulates the growth and activity of black holes and their host galaxies

Abstract: In the early Universe, while galaxies were still forming, black holes as massive as a billion solar masses powered quasars. Supermassive black holes are found at the centres of most galaxies today, where their masses are related to the velocity dispersions of stars in their host galaxies and hence to the mass of the central bulge of the galaxy. This suggests a link between the growth of the black holes and their host galaxies, which has indeed been assumed for a number of years. But the origin of the observed relation between black hole mass and stellar velocity dispersion, and its connection with the evolution of galaxies, have remained unclear. Here we report simulations that simultaneously follow star formation and the growth of black holes during galaxy-galaxy collisions. We find that, in addition to generating a burst of star formation, a merger leads to strong inflows that feed gas to the supermassive black hole and thereby power the quasar. The energy released by the quasar expels enough gas to quench both star formation and further black hole growth. This determines the lifetime of the quasar phase (approaching 100 million years) and explains the relationship between the black hole mass and the stellar velocity dispersion.

Modelling feedback from stars and black holes in galaxy mergers

Abstract: We describe techniques for incorporating feedback from star formation and black hole (BH) accretion into simulations of isolated and merging galaxies. At present, the details of these processes cannot be resolved in simulations on galactic scales. Our basic approach therefore involves forming coarse-grained representations of the properties of the interstellar medium (ISM) and BH accretion starting from basic physical assumptions, so that the impact of these effects can be included on resolved scales. We illustrate our method using a multiphase description of star-forming gas. Feedback from star formation pressurizes highly overdense gas, altering its effective equation of state (EOS). We show that this allows the construction of stable galaxy models with much larger gas fractions than possible in earlier numerical work. We extend the model by including a treatment of gas accretion onto central supermassive BHs in galaxies. Assuming thermal coupling of a small fraction of the bolometric luminosity of accreting BHs to the surrounding gas, we show how this feedback regulates the growth of BHs. In gas-rich mergers of galaxies, we observe a complex interplay between starbursts and central active galactic nuclei (AGN) activity when the tidal interaction triggers intense nuclear inflows of gas. Once an accreting supermassive BH has grown to a critical size, feedback terminates its further growth and expels gas from the central region in a powerful quasar-driven wind. Our simulation methodology is therefore able to address the coupled processes of gas dynamics, star formation and BH accretion during the formation of galaxies.

Modeling feedback from stars and black holes in galaxy mergers

Abstract: We describe techniques for incorporating feedback from star formation and black hole accretion into simulations of isolated and merging galaxies. At present, the details of these processes cannot be resolved in simulations on galactic scales. Our basic approach therefore involves forming coarse-grained representations of the properties of the interstellar medium and black hole accretion starting from basic physical assumptions, so that the impact of these effects can be included on resolved scales. We illustrate our method using a multiphase description of star-forming gas. Feedback from star formation pressurises highly overdense gas, altering its effective equation of state. We show that this allows the construction of stable galaxy models with much larger gas fractions than possible in earlier numerical work. We extend the model by including a treatment of gas accretion onto central supermassive black holes in galaxies. Assuming thermal coupling of a small fraction of the bolometric luminosity of accreting black holes to the surrounding gas, we show how this feedback regulates the growth of black holes. In gas-rich mergers of galaxies, we observe a complex interplay between starbursts and central AGN activity when the tidal interaction triggers intense nuclear inflows of gas. Once an accreting supermassive black hole has grown to a critical size, feedback terminates its further growth, and expels gas from the central region in a powerful quasar-driven wind. Our simulation methodology is therefore able to address the coupled processes of gas dynamics, star formation, and black hole accretion during the formation of galaxies.

A Unified, Merger-driven Model of the Origin of Starbursts, Quasars, the Cosmic X-Ray Background, Supermassive Black Holes, and Galaxy Spheroids

Abstract: We present an evolutionary model for starbursts, quasars, and spheroidal galaxies in which mergers between gas-rich galaxies drive nuclear inflows of gas, producing starbursts and feeding the buried growth of supermassive black holes (BHs) until feedback expels gas and renders a briefly visible optical quasar. The quasar lifetime and obscuring column density depend on both the instantaneous and peak quasar luminosity, and we determine this dependence using a large set of galaxy merger simulations varying galaxy properties, orbital geometry, and gas physics. We use these fits to deconvolve observed quasar luminosity functions and obtain the evolution of the formation rate of quasars with peak luminosity, (Lpeak, z). Quasars spend extended periods at luminosities well below peak, so (Lpeak) has a maximum corresponding to the break in the observed luminosity function. From (Lpeak) and our simulations, we obtain self-consistent hard and soft X-ray and optical luminosity functions and predict many observables at multiple redshifts, including column density distributions of optical and X-ray samples, the luminosity function of broad-line quasars in X-ray samples and broad-line fraction versus luminosity, active BH mass functions, the distribution of Eddington ratios, the mass function of relic BHs and total BH mass density, and the cosmic X-ray background. In every case, our predictions agree well with observed estimates, without invoking ad hoc assumptions about source properties or distributions. We provide a library of Monte Carlo realizations of our models for comparison with observations.

A Unified, Merger-Driven Model for the Origin of Starbursts, Quasars, the Cosmic X-Ray Background, Supermassive Black Holes and Galaxy Spheroids

Abstract: We present an evolutionary model for starbursts, quasars, and spheroidal galaxies in which mergers between gas-rich galaxies drive nuclear inflows of gas, producing intense starbursts and feeding the buried growth of supermassive black holes (BHs) until feedback expels gas and renders a briefly visible optical quasar. The quasar lifetime and obscuring column density depend on both the instantaneous and peak luminosity of the quasar, and we determine this dependence using a large set of simulations of galaxy mergers varying host galaxy properties, orbital geometry, and gas physics. We use these fits to deconvolve observed quasar luminosity functions (LFs) and obtain the evolution of the formation rate of quasars with a certain peak luminosity, n(L_peak,z). Quasars spend extended periods of time at luminosities well below peak, and so n(L_peak) has a maximum corresponding to the 'break' in the observed LF, falling off at both brighter and fainter luminosities. From n(L_peak) and our simulation results, we obtain self-consistent fits to hard and soft X-ray and optical quasar LFs and predict many observables, including: column density distributions of optical and X-ray samples, the LF of broad-line quasars in X-ray samples and the broad-line fraction as a function of luminosity, active BH mass functions, the distribution of Eddington ratios at z~0-2, the z=0 mass function of relic BHs and total mass density of BHs, and the cosmic X-ray background. In every case, our predictions agree well with observed estimates, and unlike previous modeling attempts, we are able to reproduce them without invoking any ad hoc assumptions about source properties or distributions. We provide a library of Monte Carlo realizations of our models for comparison with observations. (Abridged)

Pattern Recognition and Machine Learning

Abstract: Probability Distributions.- Linear Models for Regression.- Linear Models for Classification.- Neural Networks.- Kernel Methods.- Sparse Kernel Machines.- Graphical Models.- Mixture Models and EM.- Approximate Inference.- Sampling Methods.- Continuous Latent Variables.- Sequential Data.- Combining Models.

The Cosmological simulation code GADGET-2

Abstract: We discuss the cosmological simulation code GADGET-2, a new massively parallel TreeSPH code, capable of following a collisionless fluid with the N-body method, and an ideal gas by means of smoothed particle hydrodynamics (SPH). Our implementation of SPH manifestly conserves energy and entropy in regions free of dissipation, while allowing for fully adaptive smoothing lengths. Gravitational forces are computed with a hierarchical multipole expansion, which can optionally be applied in the form of a TreePM algorithm, where only short-range forces are computed with the ‘tree’ method while long-range forces are determined with Fourier techniques. Time integration is based on a quasi-symplectic scheme where long-range and short-range forces can be integrated with different time-steps. Individual and adaptive short-range time-steps may also be employed. The domain decomposition used in the parallelization algorithm is based on a space-filling curve, resulting in high flexibility and tree force errors that do not depend on the way the domains are cut. The code is efficient in terms of memory consumption and required communication bandwidth. It has been used to compute the first cosmological N-body simulation with more than 10 10 dark matter particles, reaching a homogeneous spatial dynamic range of 10 5 per dimension in a three-dimensional box. It has also been used to carry out very large cosmological SPH simulations that account for radiative cooling and star formation, reaching total particle numbers of more than 250 million. We present the algorithms used by the code and discuss their accuracy and performance using a number of test problems. GADGET-2 is publicly released to the research community. Ke yw ords: methods: numerical ‐ galaxies: interactions ‐ dark matter.

Simulations of the formation, evolution and clustering of galaxies and quasars

Abstract: The cold dark matter model has become the leading theoretical picture for the formation of structure in the Universe. This model, together with the theory of cosmic inflation, makes a clear prediction for the initial conditions for structure formation and predicts that structures grow hierarchically through gravitational instability. Testing this model requires that the precise measurements delivered by galaxy surveys can be compared to robust and equally precise theoretical calculations. Here we present a simulation of the growth of dark matter structure using 2,1603 particles, following them from redshift z = 127 to the present in a cube-shaped region 2.230 billion lightyears on a side. In postprocessing, we also follow the formation and evolution of the galaxies and quasars. We show that baryon-induced features in the initial conditions of the Universe are reflected in distorted form in the low-redshift galaxy distribution, an effect that can be used to constrain the nature of dark energy with future generations of observational surveys of galaxies.

The Observation of Gravitational Waves from a Binary Black Hole Merger

Abstract: On September 14, 2015 at 09:50:45 UTC the two detectors of the Laser Interferometer Gravitational-Wave Observatory simultaneously observed a transient gravitational-wave signal. The signal sweeps upwards in frequency from 35 to 250 Hz with a peak gravitational-wave strain of 1.0×10(-21). It matches the waveform predicted by general relativity for the inspiral and merger of a pair of black holes and the ringdown of the resulting single black hole. The signal was observed with a matched-filter signal-to-noise ratio of 24 and a false alarm rate estimated to be less than 1 event per 203,000 years, equivalent to a significance greater than 5.1σ. The source lies at a luminosity distance of 410(-180)(+160) Mpc corresponding to a redshift z=0.09(-0.04)(+0.03). In the source frame, the initial black hole masses are 36(-4)(+5)M⊙ and 29(-4)(+4)M⊙, and the final black hole mass is 62(-4)(+4)M⊙, with 3.0(-0.5)(+0.5)M⊙c(2) radiated in gravitational waves. All uncertainties define 90% credible intervals. These observations demonstrate the existence of binary stellar-mass black hole systems. This is the first direct detection of gravitational waves and the first observation of a binary black hole merger.