Stelios Kazantzidis

The Effect of Gas Cooling on the Shapes of Dark Matter Halos

We analyze the effect of dissipation on the shapes of dark matter (DM) halos using high-resolution cosmological gasdynamics simulations of clusters and galaxies in the ΛCDM cosmology. We find that halos formed in simulations with gas cooling are significantly more spherical than corresponding halos formed in adiabatic simulations. Gas cooling results in an average increase of the principle axis ratios of halos by ~0.2-0.4 in the inner regions. The systematic difference decreases slowly with radius but persists almost to the virial radius. We argue that the differences in simulations with and without cooling arise both during periods of quiescent evolution, when gas cools and condenses toward the center, and during major mergers. We perform a series of high-resolution N-body simulations to study the shapes of remnants in major mergers of DM halos and halos with embedded stellar disks. In the DM halo-only mergers, the shape of the remnants depends only on the orbital angular momentum of the encounter and not on the internal structure of the halos. However, significant shape changes in the DM distribution may result if stellar disks are included. In this case the shape of the DM halos is correlated with the morphology of the stellar remnants.

Generating Equilibrium Dark Matter Halos: Inadequacies of the Local Maxwellian Approximation

We describe an algorithm for constructing N-body realizations of equilibrium spherical systems. A general form for the mass density ? is used, making it possible to represent most of the popular density profiles found in the literature, including the cuspy density profiles found in high-resolution cosmological simulations. We demonstrate explicitly that our models are in equilibrium. In contrast, many existing N-body realizations of isolated systems have been constructed under the assumption that the local velocity distribution is Maxwellian. We show that a Maxwellian halo with an initial ?(r) r-1 central density cusp immediately develops a constant-density core. Moreover, after just one crossing time the orbital anisotropy has changed over the entire system, and the initially isotropic model becomes radially anisotropic. These effects have important implications for many studies, including the survival of substructure in cold dark matter (CDM) models. Comparing the evolution and mass-loss rate of isotropic Maxwellian and self-consistent Navarro, Frenk, & White (NFW), satellites orbiting inside a static host CDM potential, we find that the former are unrealistically susceptible to tidal disruption. Thus, recent studies of the mass-loss rate and disruption timescales of substructure in CDM models may be compromised by using the Maxwellian approximation. We also demonstrate that a radially anisotropic, self-consistent NFW satellite loses mass at a rate several times higher than that of its isotropic counterpart on the same external tidal field and orbit.

Rapid Formation of Supermassive Black Hole Binaries in Galaxy Mergers with Gas

Supermassive black holes (SMBHs) are a ubiquitous component of the nuclei of galaxies. It is normally assumed that after the merger of two massive galaxies, a SMBH binary will form, shrink because of stellar or gas dynamical processes, and ultimately coalesce by emitting a burst of gravitational waves. However, so far it has not been possible to show how two SMBHs bind during a galaxy merger with gas because of the difficulty of modeling a wide range of spatial scales. Here we report hydrodynamical simulations that track the formation of a SMBH binary down to scales of a few light years after the collision between two spiral galaxies. A massive, turbulent, nuclear gaseous disk arises as a result of the galaxy merger. The black holes form an eccentric binary in the disk in less than 1 million years as a result of the gravitational drag from the gas rather than from the stars.

THE CAUSES OF HALO SHAPE CHANGES INDUCED BY COOLING BARYONS : DISKS VERSUS SUBSTRUCTURES

Cold dark matter cosmogony predicts triaxial dark matter halos, whereas observations find quite round halos. This is most likely due to the condensation of baryons leading to rounder halos. We examine the halo phase space distribution basis for such shape changes. Triaxial halos are supported by box orbits, which pass arbitrarily close to the density center. The decrease in triaxiality caused by baryons is thought to be due to the scattering of these orbits. We test this hypothesis with simulations of disks grown inside triaxial halos. After the disks are grown we check whether the phase space structure has changed by evaporating the disks and comparing the initial and final states. While the halos are substantially rounder when the disk is at full mass, their final shape after the disk is evaporated is not much different from the initial. Likewise, the halo becomes (more) radially anisotropic when the disk is grown, but the final anisotropy is consistent with the initial. Only if the baryons are unreasonably compact or massive does the halo change irreversibly. We show that the character of individual orbits is not generally changed by the growing mass. Thus, the central condensation of baryons does not destroy enough box orbits to cause the shape change. Rather, box orbits merely become rounder along with the global potential. However, if angular momentum is transferred to the halo, either via satellites or via bars, a large irreversible change in the halo distribution occurs. The ability of satellites to alter the phase space distribution of the halo is of particular concern to galaxy formation simulations since halo triaxiality can profoundly influence the evolution of disks.

INSIDE OUT AND UPSIDE DOWN: TRACING THE ASSEMBLY OF A SIMULATED DISK GALAXY USING MONO-AGE STELLAR POPULATIONS

We analyze the present day structure and assembly history of a high-resolution hydrodynamic simulation of the formation of a Milky-Way-(MW)-like disk galaxy, from the Eris simulation suite, dissecting it into cohorts of stars formed at different epochs of cosmic history. At z = 0, stars with t form 3 are quickly scattered into rounded, kinematically hot configurations. The oldest disk cohorts form in structures that are radially compact and relatively thick, while subsequent cohorts form in progressively larger, thinner, colder configurations from gas with increasing levels of rotational support. The disk thus forms inside out in a radial sense and upside down in a vertical sense. Secular heating and radial migration influence the final state of each age cohort, but the changes they produce are small compared to the trends established at formation. The predicted correlations of stellar age with spatial and kinematic structure are in good qualitative agreement with the correlations observed for mono-abundance stellar populations in the MW.

ON THE EFFICIENCY OF THE TIDAL STIRRING MECHANISM FOR THE ORIGIN OF DWARF SPHEROIDALS: DEPENDENCE ON THE ORBITAL AND STRUCTURAL PARAMETERS OF THE PROGENITOR DISKY DWARFS

The tidal stirring model posits the formation of dwarf spheroidal galaxies (dSphs) via the tidal interactions between late-type, rotationally supported dwarfs and Milky-Way-sized host galaxies. Using a comprehensive set of collisionless N-body simulations, we investigate the efficiency of the tidal stirring mechanism for the origin of dSphs. In particular, we examine the degree to which the tidal field of the primary galaxy affects the sizes, masses, shapes, and kinematics of the disky dwarfs for a range of dwarf orbital and structural parameters. Our study is the first to employ self-consistent, equilibrium models for the progenitor dwarf galaxies constructed from a composite distribution function and consisting of exponential stellar disks embedded in massive, cosmologically motivated dark matter halos. Exploring a wide variety of dwarf orbital configurations and initial structures, we demonstrate that in the majority of cases the disky dwarfs experience significant mass loss and their stellar distributions undergo a dramatic morphological, as well as dynamical, transformation. Specifically, the stellar components evolve from disks to bars and finally to pressure-supported, spheroidal systems with kinematic and structural properties akin to those of the classic dSphs in the Local Group (LG) and similar environments. The self-consistency of the adopted dwarf models is crucial for confirming this complex transformation process via tidally induced dynamical instabilities and impulsive tidal heating of the stellar distribution. Our results suggest that such tidal transformations should be common occurrences within the currently favored cosmological paradigm and highlight the key factor responsible for an effective metamorphosis to be the strength of the tidal shocks at the pericenters of the orbit. We also demonstrate that the combination of short orbital times and small pericentric distances, characteristic of dwarfs being accreted by their hosts at high redshift, induces the strongest and most complete transformations. Our models also indicate that the efficiency of the transformation via tidal stirring is affected significantly by the structure of the progenitor disky dwarfs. While the mass-to-light ratios, M/L, of the dwarf galaxies typically decrease monotonically with time as the extended dark matter halos are efficiently tidally stripped, we identify a few cases where this trend is reversed later in the evolution when stellar mass loss becomes more effective. We also find that the dwarf remnants satisfy the relation , where σ* is the one-dimensional, central stellar velocity dispersion and V max is the maximum halo circular velocity, which has intriguing implications for the missing satellites problem. Assuming that the distant dSphs in the LG, such as Leo I, Tucana, and Cetus, are the products of tidal stirring, our findings suggest that these galaxies should have only been partially stirred by the tidal field of their hosts. We thus predict that these remote dwarfs should exhibit higher values of V rot/σ*, where V rot is the stellar rotational velocity, compared with those of dSphs located closer to the primary galaxies. Overall, we conclude that the action of tidal forces from the hosts constitutes a crucial evolutionary mechanism for shaping the nature of dwarf galaxies in environments such as that of the LG. Environmental mechanisms of this type should thus be included as ingredients in models of dwarf galaxy formation and evolution.

Direct formation of supermassive black holes via multi-scale gas inflows in galaxy mergers

Observations of distant quasars indicate that supermassive black holes of billions of solar masses already existed less than a billion years after the Big Bang. Models in which the ‘seeds’ of such black holes form by the collapse of primordial metal-free stars cannot explain the rapid appearance of these supermassive black holes because gas accretion is not sufficiently efficient. Alternatively, these black holes may form by direct collapse of gas within isolated protogalaxies, but current models require idealized conditions, such as metal-free gas, to prevent cooling and star formation from consuming the gas reservoir. Here we report simulations showing that mergers between massive protogalaxies naturally produce the conditions for direct collapse into a supermassive black hole with no need to suppress cooling and star formation. Merger-driven gas inflows give rise to an unstable, massive nuclear gas disk of a few billion solar masses, which funnels more than 108 solar masses of gas to a sub-parsec-scale gas cloud in only 100,000 years. The cloud undergoes gravitational collapse, which eventually leads to the formation of a massive black hole. The black hole can subsequently grow to a billion solar masses on timescales of about 108 years by accreting gas from the surrounding disk.

Early gas stripping as the origin of the darkest galaxies in the Universe

The known galaxies most dominated by dark matter (Draco, Ursa Minor and Andromeda IX) are satellites of the Milky Way and the Andromeda galaxies. They are members of a class of faint galaxies, devoid of gas, known as dwarf spheroidals, and have by far the highest ratio of dark to luminous matter. None of the models proposed to unravel their origin can simultaneously explain their exceptional dark matter content and their proximity to a much larger galaxy. Here we report simulations showing that the progenitors of these galaxies were probably gas-dominated dwarf galaxies that became satellites of a larger galaxy earlier than the other dwarf spheroidals. We find that a combination of tidal shocks and ram pressure swept away the entire gas content of such progenitors about ten billion years ago because heating by the cosmic ultraviolet background kept the gas loosely bound: a tiny stellar component embedded in a relatively massive dark halo survived until today. All luminous galaxies should be surrounded by a few extremely dark-matter-dominated dwarf spheroidal satellites, and these should have the shortest orbital periods among dwarf spheroidals because they were accreted early.

COLD DARK MATTER SUBSTRUCTURE AND GALACTIC DISKS. II. DYNAMICAL EFFECTS OF HIERARCHICAL SATELLITE ACCRETION

We perform a set of fully self-consistent, dissipationless N-body simulations to elucidate the dynamical response of thin galactic disks to bombardment by cold dark matter (CDM) substructure. Our method combines (1) cosmological simulations of the formation of Milky Way (MW)-sized CDM halos to derive the properties of substructure, and (2) controlled numerical experiments of consecutive subhalo impacts onto an initially thin, fully formed MW-type disk galaxy. The present study is the first to account for the evolution of satellite populations over cosmic time in such an investigation of disk structure. In contrast to what can be inferred from statistics of the z = 0 surviving substructure, we find that accretions of massive subhalos onto the central regions of host halos, where the galactic disks reside, since z ~ 1 should be common. One host halo accretion history is used to initialize the controlled simulations of satellite-disk encounters. The specific merger history involves six dark matter substructures, with initial masses in the range ~20%-60% of the disk mass and of comparable size to the disk, crossing the central regions of their host in the past ~8 Gyr. We show that these accretion events severely perturb the thin galactic disk and produce a wealth of distinctive dynamical signatures on its structure and kinematics. These include (1) considerable thickening and heating at all radii, with the disk thickness and velocity ellipsoid nearly doubling at the solar radius; (2) prominent flaring associated with an increase in disk thickness greater than a factor of 4 in the disk outskirts; (3) surface density excesses at large radii, beyond ~5 disk scale lengths, resembling those of the observed antitruncated disks; (4) long-lived, lopsidedness at levels similar to those measured in observational samples of disk galaxies; and (5) substantial tilting. The interaction with the most massive subhalo in the simulated accretion history drives the disk response while subsequent bombardment is much less efficient at disturbing the disk. We also explore a variety of disk and satellite properties that influence these responses. We conclude that substructure-disk encounters of the kind expected in the ΛCDM paradigm play a significant role in setting the structure of disk galaxies and driving galaxy evolution.

The Fate of Supermassive Black Holes and the Evolution of the MBH-σ Relation in Merging Galaxies: The Effect of Gaseous Dissipation

We analyze the effect of dissipation on the orbital evolution of supermassive black holes (SMBHs) using high- resolution self-consistent gasdynamical simulations of binary equal- and unequal-mass mergers of disk galaxies. The galaxy models are consistent with the LCDM paradigm of structure formation, and the simulations include the effects of radiative cooling and star formation. We find that equal-mass mergers always lead to the formation of a close SMBH pair at the center of the remnant, with separations limited solely by the adopted force resolution of ∼100 pc. Instead, the final SMBH separation in unequal-mass mergers depends sensitively on how the central structure of the merging galaxies is modified by dissipation. In the absence of dissipation, the satellite galaxy can be entirely disrupted before the merger is completed, leaving its SMBH wandering at a distance too far from the center of the remnant for the formation of a close pair. In contrast, gas cooling facilitates the pairing process by increasing the resilience of the companion galaxy to tidal disruption. Moreover, we demonstrate that merging disk galaxies constructed to obey the MBH-j relation move relative to it depending on whether they undergo a dissipational or collisionless merger, regardless of the mass ratio of the merging systems. Collisionless simulations reveal that remnants tend to move away from the mean relation, highlighting the role of gas-poor mergers as a possible source of scatter. In dissipational mergers, the interplay between strong gas inflows associated with the formation of massive nuclear disks and the consumption of gas by star formation provides the necessary fuel to the SMBHs and allows the merger remnants to satisfy the relation