Researcher profile

Stelios Kazantzidis

Stelios Kazantzidis contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
33works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

33 published item(s)

preprint2015arXiv

Black Hole Starvation and Bulge Evolution in a Milky Way-like Galaxy

We present a new zoom-in hydrodynamical simulation, "Erisbh", which follows the cosmological evolution and feedback effects of a supermassive black hole at the center of a Milky Way-type galaxy. ErisBH shares the same initial conditions, resolution, recipes of gas cooling, star formation and feedback, as the close Milky Way-analog "Eris", but it also includes prescriptions for the formation, growth and feedback of supermassive black holes. We find that the galaxy's central black hole grows mainly through mergers with other black holes coming from infalling satellite galaxies. The growth by gas accretion is minimal because very little gas reaches the sub-kiloparsec scales. The final black hole is, at z=0, about 2.6 million solar masses and it sits closely to the position of SgrA* on the MBH-MBulge and MBH-sigma planes, in a location consistent with what observed for pseudobulges. Given the limited growth due to gas accretion, we argue that the mass of the central black hole should be above 10^5 solar masses already at z~8. The effect of AGN feedback on the host galaxy is limited to the very central few hundreds of parsecs. Despite being weak, AGN feedback seems to be responsible for the limited growth of the central bulge with respect to the original Eris, which results in a significantly flatter rotation curve in the inner few kiloparsecs. Moreover, the disk of ErisBH is more prone to instabilities, as its bulge is smaller and its disk larger then Eris. As a result, the disk of ErisBH undergoes a stronger dynamical evolution relative to Eris and around z=0.3 a weak bar grows into a strong bar of a few disk scale lengths in size. The bar triggers a burst of star formation in the inner few hundred parsecs, provides a modest amount of new fuel to the central black hole, and causes the bulge of ErisBH to have, by z=0, a box/peanut morphology.(Abridged)

preprint2013arXiv

Can we measure the slopes of density profiles in dwarf spheroidal galaxies?

Using collisionless N-body simulations of dwarf galaxies orbiting the Milky Way (MW) we construct realistic models of dwarf spheroidal (dSph) galaxies of the Local Group. The dwarfs are initially composed of stellar disks embedded in dark matter haloes with different inner density slopes and are placed on an eccentric orbit typical for MW subhaloes. After a few Gyr of evolution the stellar component is triaxial as a result of bar instability induced by tidal forces. Observing the simulated dwarfs along three principal axes of the stellar component we create mock data sets and determine the their half-light radii and line-of-sight velocity dispersions. Using the estimator proposed by Wolf et al. we calculate masses within half-light radii. The masses obtained this way are over(under)estimated by up to a factor of two when the line of sight is along the longest (shortest) axis of the stellar component. We then divide the initial stellar distribution into an inner and outer population and trace their evolution in time. The two populations, although affected by tidal forces, retain different density profiles even after a few Gyr. We measure the half-light radii and velocity dispersions of the stars in the two populations along different lines of sight and use them to estimate the slope of the mass distribution in the dwarfs following the method proposed by Walker & Penarrubia. The inferred slopes are systematically over- or underestimated, depending on the line of sight. In particular, when the dwarf is seen along the longest axis of the stellar component, a significantly shallower density profile is inferred than the real one measured from the simulations. Since most dSphs are non-spherical and their orientation with respect to our line of sight is unknown, the method can be reliably applied only to a large sample of dwarfs when these systematic errors are expected to be diminished.

preprint2013arXiv

Disk Assembly and the M_BH-sigma Relation of Supermassive Black Holes

Recent Hubble Space Telescope (HST) observations have revealed that a majority of active galactic nuclei (AGN) at z ~ 1-3 are resident in isolated disk galaxies, contrary to the usual expectation that AGN are triggered by mergers. Here we develop a new test of the cosmic evolution of supermassive black holes (SMBHs) in disk galaxies by considering the local population of SMBHs. We show that substantial SMBH growth in spiral galaxies is required as disks assemble. SMBHs exhibit a tight relation between their mass and the velocity dispersion of the spheroid within which they reside, the M_BH-sigma relation. In disk galaxies the bulge is the spheroid of interest. We explore the evolution of the M_BH-sigma relation when bulges form together with SMBHs on the M_BH-sigma relation and then slowly reform a disk around them. The formation of the disk compresses the bulge raising its sigma. We present evidence for such compression in the form of larger velocity dispersion of classical bulges compared with elliptical galaxies at the same mass. This compression leads to an offset in the M_BH-sigma relation if it is not accompanied by an increased M_BH. We quantify the expected offset based on photometric data and show that, on average, SMBHs must grow by ~ 50-65% just to remain on the M_BH-sigma relation. We find no significant offset in the M_BH-sigma relations of classical bulges and of ellipticals, implying that SMBHs have been growing along with disks. Our simulations demonstrate that SMBH growth is necessary for the local population of disk galaxies to have remained on the M_BH-sigma relation.

preprint2013arXiv

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 &#34;Eris&#34; simulation suite, dissecting it into cohorts of stars formed at different epochs of cosmic history. At z=0, stars with t_form < 2 Gyr mainly occupy the stellar spheroid, with the oldest (earliest forming) stars having more centrally concentrated profiles. The younger age cohorts populate disks of progressively longer radial scale length and shorter vertical scale height. At a given radius, the vertical density profiles and velocity dispersions of stars vary smoothly as a function of age, and the superposition of old, vertically-extended and young, vertically-compact cohorts gives rise to a double-exponential profile like that observed in the MW. Turning to formation history, we find that the trends of spatial structure and kinematics with stellar age are largely imprinted at birth, or immediately thereafter. Stars that form during the active merger phase at z>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 &#34;inside-out&#34; in a radial sense and &#34;upside-down&#34; 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.

preprint2013arXiv

Tidal Stirring of Disky Dwarfs with Shallow Dark Matter Density Profiles: Enhanced Transformation into Dwarf Spheroidals

(Abridged) The origin of dSphs in the Local Group (LG) remains an enigma. The tidal stirring model posits that late-type, rotationally-supported dwarfs resembling present-day dwarf irregular (dIrr) galaxies can transform into dSphs via interactions with Milky Way-sized hosts. Using collisionless N-body simulations, we investigate for the first time how tidal stirring depends on the dark matter (DM) density distribution in the central stellar region of the progenitor disky dwarf. Specifically, we explore various asymptotic inner slopes gamma of the dwarf DM density profiles (rho \propto r^{-gamma} as r -> 0). For a given orbit inside the primary, rotationally-supported dwarfs embedded in DM halos with core-like density distributions (gamma = 0.2) and mild density cusps (gamma = 0.6) demonstrate a substantially enhanced likelihood and efficiency of transformation into dSphs compared to their counterparts with steeper DM density profiles (gamma = 1). Such shallow DM distributions are akin to those of observed dIrrs, highlighting tidal stirring as a plausible model for the LG morphology-density relation. When gamma <1, a single pericentric passage can induce dSph formation and disky dwarfs on low-eccentricity or large-pericenter orbits are able to transform into dSphs; these new results allow the tidal stirring model to explain the existence of virtually all known dSphs across a wide range of distances from their hosts. A subset of rotationally-supported dwarfs with gamma <1 are eventually disrupted by the primary; those that survive as dSphs are generally on orbits that are biased towards lower eccentricities and/or larger pericenters relative to those of typical CDM satellites. The latter could explain the rather peculiar orbits of several classic LG dSphs such as Fornax, Leo I, Tucana, and Cetus.

preprint2013arXiv

Tidal tails of dwarf galaxies on different orbits around the Milky Way

We present a phenomenological description of the properties of tidal tails forming around dwarf galaxies orbiting the Milky Way. For this purpose we use collisionless N-body simulations of dwarfs initially composed of a disk embedded in an NFW dark matter halo. The dwarfs are placed on seven orbits around the Milky Way-like host, differing in size and eccentricity, and their evolution is followed for 10 Gyr. In addition to the well-studied morphological and dynamical transformation of the dwarf&#39;s main body, the tidal stripping causes them to lose a substantial fraction of mass both in dark matter and stars which form pronounced tidal tails. We focus on the properties of the stellar component of the tidal tails thus formed. We first discuss the break radii in the stellar density profile defining the transition to tidal tails as the radii where the profile becomes shallower and relate them to the classically defined tidal radii. We then calculate the relative density and velocity of the tails at a few break radii as a function of the orbital phase. Next, we measure the orientation of the tails with respect to an observer placed at the centre of the Milky Way. The tails are perpendicular to this line of sight only for a short period of time near the pericentre. For most of the time the angles between the tails and this line of sight are low, with orbit-averaged medians below 42 degrees for all, even the almost circular orbit. The median angle is typically lower while the maximum relative density higher for more eccentric orbits. The combined effects of relative density and orientation of the tails suggest that they should be easiest to detect for dwarf galaxies soon after their pericentre passage.

preprint2012arXiv

How to make an ultra-faint dwarf spheroidal galaxy: tidal stirring of disky dwarfs with shallow dark matter density profiles

In recent years the Sloan Digital Sky Survey has unraveled a new population of ultra-faint dwarf galaxies (UFDs) in the vicinity of the Milky Way (MW) whose origin remains a puzzle. Using a suite of collisionless N-body simulations, we investigate the formation of UFDs in the context of the tidal stirring model for the formation of dwarf spheroidal galaxies in the Local Group (LG). Our simulations are designed to reproduce the tidal interactions between MW-sized host galaxies and rotationally supported dwarfs embedded in 10^9 M_sun dark matter (DM) halos. We explore a variety of inner density slopes ρ\propto r^{-α} for the dwarf DM halos, ranging from core-like (α= 0.2) to cuspy (α= 1), and different dwarf orbital configurations. Our experiments demonstrate that UFDs can be produced via tidal stirring of disky dwarfs on relatively tight orbits, consistent with a redshift of accretion by the host galaxy of z \sim 1, and with intermediate values for the halo inner density slopes (ρ\propto r^{-0.6}). The inferred slopes are in excellent agreement with those resulting from both the modeling of the rotation curves of dwarf galaxies and recent cosmological simulations of dwarf galaxy formation. Comparing the properties of observed UFDs with those of their simulated counterparts, we find remarkable similarities in terms of basic observational parameters. We conclude that tidal stirring of rotationally supported dwarfs represents a viable mechanism for the formation of UFDs in the LG environment.

preprint2012arXiv

Stellar populations in tidally stirred dwarf galaxies

Using N-body simulations we study the evolution of separate stellar populations in dwarf galaxies in the context of the tidal stirring scenario for the formation of dwarf spheroidal (dSph) galaxies in the Local Group. The dwarf galaxies, initially composed of a stellar disk and a dark matter halo, are placed on seven different orbits around the Milky Way. The stars are divided into two populations, within and outside the half-light radius, and their positions are followed for 10 Gyr. We find that the populations retain different density distributions even over such long timescales. Some of the stars of the outer population migrate to the central part of the dwarf forming an extended core while the stars of the inner population develop a tail in the outer parts. In addition, the outer population is more heavily stripped by tidal forces from the Milky Way and may become subdominant at all radii on tight enough orbits. We conclude that the tidal stirring model is fully compatible with the presence of multiple stellar populations in dSph galaxies.

preprint2012arXiv

The shapes of Milky Way satellites: looking for signatures of tidal stirring

We study the shapes of Milky Way satellites in the context of the tidal stirring scenario for the formation of dwarf spheroidal galaxies. The standard procedures used to measure shapes involve smoothing and binning of data and thus may not be sufficient to detect structural properties like bars, which are usually subtle in low surface brightness systems. Taking advantage of the fact that in nearby dwarfs photometry of individual stars is available we introduce discrete measures of shape based on the two-dimensional inertia tensor and the Fourier bar mode. We apply these measures of shape first to a variety of simulated dwarf galaxies formed via tidal stirring of disks embedded in dark matter halos and orbiting the Milky Way. In addition to strong mass loss and randomization of stellar orbits, the disks undergo morphological transformation that typically involves the formation of a triaxial bar after the first pericenter passage. These tidally induced bars persist for a few Gyr before being shortened towards a more spherical shape if the tidal force is strong enough. We test this prediction by measuring in a similar way the shape of nearby dwarf galaxies, satellites of the Milky Way. We detect inner bars in Ursa Minor, Sagittarius, LMC and possibly Carina. In addition, six out of eleven studied dwarfs show elongated stellar distributions in the outer parts that may signify transition to tidal tails. We thus find the shapes of Milky Way satellites to be consistent with the predictions of the tidal stirring model.

preprint2011arXiv

Evolutionary tracks of tidally stirred disky dwarf galaxies

Using collisionless N-body simulations, we investigate the tidal evolution of late-type, rotationally supported dwarfs inside Milky Way-sized host galaxies. Our study focuses on a wide variety of dwarf orbital configurations and initial structures. During the evolution, the disky dwarfs undergo strong mass loss, the stellar disks are transformed into spheroids, and rotation is replaced by random motions of the stars. Thus, the late-type progenitors are transformed into early-type dwarfs as envisioned by the tidal stirring model for the formation of dwarf spheroidal (dSph) galaxies in the Local Group. We determine the photometric properties of the dwarfs, including the total visual magnitude, the half-light radius and the central surface brightness as they would be measured by an observer near the galactic center. Special emphasis is also placed on studying their kinematics and shapes. We demonstrate that the measured values are biased by a number of observational effects including the increasing angle of the observation cone near the orbital pericenter, the fact that away from the pericenter the tidal tails are typically oriented along the line of sight, and the fact that for most of the evolution the stellar components of the dwarfs are triaxial ellipsoids whose major axis tumbles with respect to the line of sight. Finally, we compare the measured properties of the simulated dwarfs to those of dwarf galaxies in the Local Group. The evolutionary tracks of the dwarfs in different parameter planes and the correlations between their different properties, especially the total magnitude and the surface brightness, strongly suggest that present-day dSph galaxies may have indeed formed from late-type progenitors as proposed by the tidal stirring scenario.

preprint2011arXiv

Formation of Dwarf Spheroidal Galaxies Via Mergers of Disky Dwarfs

We perform collisionless N-body simulations to investigate whether binary mergers between rotationally-supported dwarfs can lead to the formation of dwarf spheroidal galaxies (dSphs). Our simulation campaign is based on a hybrid approach combining cosmological simulations and controlled numerical experiments. We select merger events from a Constrained Local UniversE (CLUES) simulation of the Local Group (LG) and record the properties of the interacting dwarf-sized halos. This information is subsequently used to seed controlled experiments of binary encounters between dwarf galaxies consisting of exponential stellar disks embedded in cosmologically-motivated dark matter halos. These simulations are designed to reproduce eight cosmological merger events, with initial masses of the interacting systems in the range ~ (5-60) x 10^7 Mo, occurring quite early in the history of the LG, more than 10 Gyr ago. We compute the properties of the merger remnants as a distant observer would and demonstrate that at least three of the simulated encounters produce systems with kinematic and structural properties akin to those of the classic dSphs in the LG. Tracing the history of the remnants in the cosmological simulation to z=0, we find that two dSph-like objects remain isolated at distances larger than 800 kpc from either the Milky Way or M31. These systems constitute plausible counterparts of the remote dSphs Cetus and Tucana which reside in the LG outskirts, far from the tidal influence of the primary galaxies. We conclude that merging of rotationally-supported dwarfs represents a viable mechanism for the formation of dSphs in the LG and similar environments.

preprint2011arXiv

Radial Mixing in Galactic Disks: The Effects of Disk Structure and Satellite Bombardment

We use a suite of numerical simulations to investigate the mechanisms and effects of radial migration of stars in disk galaxies like the Milky Way (MW). An isolated, collisionless stellar disk with a MW-like scale-height shows only the radial &#34;blurring&#34; expected from epicyclic orbits. Reducing the disk thickness or adding gas to the disk substantially increases the level of radial migration, induced by interaction with transient spiral arms and/or a central bar. We also examine collisionless disks subjected to gravitational perturbations from a cosmologically motivated satellite accretion history. In the perturbed disk that best reproduces the observed properties of the MW, 20% of stars that end up in the solar annulus 7 kpc < R < 9 kpc started at R < 6 kpc, and 7% started at R > 10 kpc. This level of migration would add considerable dispersion to the age-metallicity relation of solar neighborhood stars. In the isolated disk models, the probability of migration traces the disk&#39;s radial mass profile, but in perturbed disks migration occurs preferentially at large radii, where the disk is more weakly bound. The orbital dynamics of migrating particles are also different in isolated and perturbed disks: satellite perturbations drive particles to lower angular momentum for a given change in radius. Thus, satellite perturbations appear to be a distinct mechanism for inducing radial migration, which can operate in concert with migration induced by bars and spiral structure. We investigate correlations between changes in radius and changes in orbital circularity or vertical energy, identifying signatures that might be used to test models and distinguish radial migration mechanisms in chemo-dynamical surveys of the MW disk.

preprint2010arXiv

Biases in mass estimates of dSph galaxies

Using a high resolution N-body simulation of a two-component dwarf galaxy orbiting in the potential of the Milky Way, we study two effects that lead to significant biases in mass estimates of dwarf spheroidal galaxies. Both are due to the strong tidal interaction of initially disky dwarfs with their host. The tidal stripping of dwarf stars leads to the formation of strong tidal tails that are typically aligned with the line of sight of an observer positioned close to the host center. The stars from the tails contaminate the kinematic samples leading to a velocity dispersion profile increasing with the projected radius and resulting in an overestimate of mass. The tidal stirring of the dwarf leads to the morphological transformation of the initial stellar disk into a bar and then a spheroid. The distribution of stars in the dwarf remains non-spherical for a long time leading to an overestimate of its mass if it is observed along the major axis and an underestimate if it seen in the perpendicular direction.

preprint2010arXiv

Heated Disc Stars in the Stellar Halo

Minor accretion events with mass ratio M_sat : M_host ~ 1:10 are common in the context of LCDM cosmology. We use high-resolution simulations of Galaxy-analogue systems to show that these mergers can dynamically eject disk stars into a diffuse light component that resembles a stellar halo both spatially and kinematically. For a variety of orbital configurations, we find that ~3-5e8 M_sun of primary stellar disk material is ejected to a distance larger than 5 kpc above the galactic plane. This ejected contribution is similar to the mass contributed by the tidal disruption of the satellite galaxy itself, though it is less extended. If we restrict our analysis to the approximate solar neighborhood in the disk plane, we find that ~1% of the initial disk stars in that region would be classified kinematically as halo stars. Our results suggest that the inner parts of galactic stellar halos contain ancient disk stars and that these stars may have been liberated in the very same events that delivered material to the outer stellar halo.

preprint2010arXiv

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

(Abridged) The tidal stirring model posits the formation of dSph galaxies via the tidal interactions between rotationally-supported dwarfs and MW-sized host galaxies. Using a set of collisionless N-body simulations, we investigate the efficiency of the tidal stirring mechanism. We explore a wide variety of dwarf orbital configurations and initial structures and demonstrate that in most cases the disky dwarfs experience significant mass loss and their stellar components undergo a dramatic morphological and dynamical transformation: 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). 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 demonstrate that the combination of short orbital times and small pericenters, characteristic of dwarfs being accreted at high redshift, induces the strongest transformations. Our models also indicate that the transformation efficiency is affected significantly by the structure of the progenitor disky dwarfs. Lastly, we find that the dwarf remnants satisfy the relation Vmax = \sqrt{3} * sigma, where sigma is the 1D, central stellar velocity dispersion and Vmax is the maximum halo circular velocity, with intriguing implications for the missing satellites problem. 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 processes of this type should thus be included as ingredients in models of dwarf galaxy formation and evolution.

preprint2010arXiv

The inner structure and kinematics of the Sagittarius dwarf galaxy as a product of tidal stirring

The tidal stirring model envisions the formation of dwarf spheroidal (dSph) galaxies in the Local Group via the tidal interaction of disky dwarf systems with a larger host galaxy like the Milky Way. These progenitor disks are embedded in extended dark halos and during the evolution both components suffer strong mass loss. In addition, the disks undergo the morphological transformation into spheroids and the transition from ordered to random motion of their stars. Using collisionless N-body simulations we construct a model for the nearby and highly elongated Sagittarius (Sgr) dSph galaxy within the framework of the tidal stirring scenario. Constrained by the present known orbit of the dwarf, the model suggests that in order to produce the majority of tidal debris observed as the Sgr stream, but not yet transform the core of the dwarf into a spherical shape, Sgr must have just passed the second pericenter of its current orbit around the Milky Way. In the model, the stellar component of Sgr is still very elongated after the second pericenter and morphologically intermediate between the strong bar formed at the first pericenter and the almost spherical shape existing after the third pericenter. This is thus the first model of the evolution of the Sgr dwarf that accounts for its observed very elliptical shape. At the present time there is very little intrinsic rotation left and the velocity gradient detected along the major axis is almost entirely of tidal origin. We model the recently measured velocity dispersion profile for Sgr assuming that mass traces light and estimate its current total mass within 5 kpc to be 5.2 x 10^8 M_sun. To have this mass at present, the model requires that the initial virial mass of Sgr must have been as high as 1.6 x 10^10 M_sun, comparable to that of the Large Magellanic Cloud, which may serve as a suitable analog for the pre-interaction, Sgr progenitor.

preprint2010arXiv

The Sphericalization of Dark Matter Halos by Galaxy Disks

Cosmological simulations indicate that cold dark matter (CDM) halos should be triaxial. Verifying observationally this theoretical prediction is, however, less than straightforward because the assembly of galaxies is expected to modify the halo shapes and to render them more axisymmetric. We use a suite of N-body simulations to investigate quantitatively the effect of the growth of a central disk galaxy on the shape of triaxial dark matter halos. As expected, the halo responds to the presence of the disk by becoming more spherical. The net effect depends only weakly on the orientation of the disk relative to the halo principal axes or the timescale of disk assembly, but strongly on the overall gravitational importance of the disk. Our results show that exponential disks whose contribution peaks at less than ~50% of their circular velocity are unable to modify noticeably the shape of the gravitational potential of their surrounding halos. Many dwarf and low surface brightness galaxies are expected to be in this regime, and therefore their detailed kinematics could be used to probe halo triaxiality, one of the basic predictions of the CDM paradigm. We argue that the complex disk kinematics of the dwarf galaxy NGC 2976 might be the reflection of a triaxial halo. Such signatures of halo triaxiality should be common in galaxies where the luminous component is subdominant.

preprint2010arXiv

Thick-disk evolution induced by the growth of an embedded thin disk

We perform collisionless N-body simulations to investigate the evolution of the structural and kinematical properties of simulated thick disks induced by the growth of an embedded thin disk. The thick disks used in the present study originate from cosmologically-common 5:1 encounters between initially-thin primary disk galaxies and infalling satellites. The growing thin disks are modeled as static gravitational potentials and we explore a variety of growing-disk parameters that are likely to influence the response of thick disks. We find that the final thick-disk properties depend strongly on the total mass and radial scale-length of the growing thin disk, and much less sensitively on its growth timescale and vertical scale-height as well as the initial sense of thick-disk rotation. Overall, the growth of an embedded thin disk can cause a substantial contraction in both the radial and vertical direction, resulting in a significant decrease in the scale-lengths and scale-heights of thick disks. Kinematically, a growing thin disk can induce a notable increase in the mean rotation and velocity dispersions of thick-disk stars. We conclude that the reformation of a thin disk via gas accretion may play a significant role in setting the structure and kinematics of thick disks, and thus it is an important ingredient in models of thick-disk formation.

preprint2010arXiv

Tidal stirring of Milky Way satellites: a simple picture with the integrated tidal force

Most of dwarf spheroidal galaxies in the Local Group were probably formed via environmental processes like the tidal interaction with the Milky Way. We study this process via N-body simulations of dwarf galaxies evolving on seven different orbits around the Galaxy. The dwarf galaxy is initially composed of a rotating stellar disk and a dark matter halo. Due to the action of tidal forces it loses mass and the disk gradually transforms into a spheroid while stellar motions become increasingly random. We measure the characteristic scale-length of the dwarf, its maximum circular velocity, mass, shape and kinematics as a function of the integrated tidal force along the orbit. The final properties of the evolved dwarfs are remarkably similar if the total tidal force they experienced was the same, independently of the actual size and eccentricity of the orbit.

preprint2009arXiv

Direct Formation of Supermassive Black Holes via Multi-Scale Gas Inflows in Galaxy Mergers

Observations of distant bright quasars suggest that billion solar mass supermassive black holes (SMBHs) were already in place less than a billion years after the Big Bang. Models in which light black hole seeds form by the collapse of primordial metal-free stars cannot explain their rapid appearance due to inefficient gas accretion. Alternatively, these black holes may form by direct collapse of gas at the center of protogalaxies. However, this requires metal-free gas that does not cool efficiently and thus is not turned into stars, in contrast with the rapid metal enrichment of protogalaxies. Here we use a numerical simulation to show that mergers between massive protogalaxies naturally produce the required central gas accumulation with no need to suppress star formation. Merger-driven gas inflows produce an unstable, massive nuclear gas disk. Within the disk a second gas inflow accumulates more than 100 million solar masses of gas in a sub-parsec scale cloud in one hundred thousand years. The cloud undergoes gravitational collapse, which eventually leads to the formation of a massive black hole. The black hole can grow to a billion solar masses in less than a billion years by accreting gas from the surrounding disk.

preprint2009arXiv

Evolution of Dark Matter Phase-Space Density Distributions in Equal-Mass Halo Mergers

We use dissipationless N-body simulations to investigate the evolution of the true coarse-grained phase-space density distribution f(x,v) in equal-mass mergers between dark matter (DM) halos. The halo models are constructed with various asymptotic power-law indices ranging from steep cusps to core-like profiles and we employ the phase-space density estimator ``Enbid&#39;&#39; developed by Sharma & Steinmetz to compute f(x,v). The adopted force resolution allows robust phase-space density profile estimates in the inner ~1% of the virial radii of the simulated systems. We confirm that mergers result in a decrease of the coarse-grained phase-space density in accordance with expectations from Mixing Theorems for collisionless systems. We demonstrate that binary mergers between identical DM halos produce remnants that retain excellent memories of the inner slopes and overall shapes of the phase-space density distribution of their progenitors. The robustness of the phase-space density profiles holds for a range of orbital energies, and a variety of encounter configurations including sequences of several consecutive merger events, designed to mimic hierarchical merging, and collisions occurring at different cosmological epochs. If the progenitor halos are constructed with appreciably different asymptotic power-law indices, we find that the inner slope and overall shape of the phase-space density distribution of the remnant are substantially closer to that of the initial system with the steepest central density cusp. These results explicitly demonstrate that mixing is incomplete in equal-mass mergers between DM halos, as it does not erase memory of the progenitor properties. Our results also confirm the recent analytical predictions of Dehnen (2005) regarding the preservation of merging self-gravitating central density cusps.

preprint2009arXiv

Evolution of the Dark Matter Phase-Space Density Distributions of LCDM Halos

We study the evolution of phase-space density during the hierarchical structure formation of LCDM halos. We compute both a spherically-averaged surrogate for phase-space density (Q) and the coarse-grained distribution function f(x,v) for dark matter particles that lie within~2 virial radii of four Milky-Way-sized dark matter halos. The estimated f(x,v) spans over four decades at any radius. Dark matter particles that end up within two virial radii of a Milky-Way-sized DM halo at $z=0$ have an approximately Gaussian distribution in log(f) at early redshifts, but the distribution becomes increasingly skewed at lower redshifts. The value corresponding to the peak of the Gaussian decreases as the evolution progresses and is well described by a power-law in (1+z). The highest values of f are found at the centers of dark matter halos and subhalos, where f can be an order of magnitude higher than in the center of the main halo. The power-law Q(r) profile likely reflects the distribution of entropy (K = sigma^2/rho^{2/3} \propto r^{1.2}), which dark matter acquires as it is accreted onto a growing halo. The estimated f(x, v), on the other hand, exhibits a more complicated behavior. Although the median coarse-grained phase-space density profile F(r) can be approximated by a power-law in the inner regions of halos and at larger radii the profile flattens significantly. This is because phase-space density averaged on small scales is sensitive to the high-f material associated with surviving subhalos, as well as relatively unmixed material (probably in streams) resulting from disrupted subhalos, which contribute a sizable fraction of matter at large radii. (ABRIDGED)

preprint2009arXiv

Pairing of Supermassive Black Holes in unequal-mass galaxy mergers

We examine the pairing process of supermassive black holes (SMBHs) down to scales of 20-100 pc using a set of N-body/SPH simulations of binary mergers of disk galaxies with mass ratios of 1:4 and 1:10. Our numerical experiments are designed to represent merger events occurring at various cosmic epochs. The initial conditions of the encounters are consistent with the LambdaCDM paradigm of structure formation, and the simulations include the effects of radiative cooling, star formation, and supernovae feedback. We find that the pairing of SMBHs depends sensitively on the amount of baryonic mass preserved in the center of the companion galaxies during the last phases of the merger. In particular, due to the combination of gasdynamics and star formation, we find that a pair of SMBHs can form in 1:10 minor mergers provided that galaxies are relatively gas-rich (gas fractions of 30% of the disk mass) and that the mergers occur at relatively high redshift (z~3), when dynamical friction timescales are shorter. Since 1:10 mergers are most common events during the assembly of galaxies, and mergers are more frequent at high redshift when galaxies are also more gas-rich, our results have positive implications for future gravitational wave experiments such as the Laser Interferometer Space Antenna.

preprint2009arXiv

The orientation and kinematics of inner tidal tails around dwarf galaxies orbiting the Milky Way

Using high-resolution collisionless N-body simulations we study the properties of tidal tails formed in the immediate vicinity of a two-component dwarf galaxy evolving in a static potential of the Milky Way (MW). The stellar component of the dwarf is initially in the form of a disk and the galaxy is placed on an eccentric orbit motivated by CDM-based cosmological simulations. We measure the orientation, density and velocity distribution of the stars in the tails. Due to the geometry of the orbit, in the vicinity of the dwarf, where the tails are densest and therefore most likely to be detectable, they are typically oriented towards the MW and not along the orbit. We report on an interesting phenomenon of `tidal tail flipping&#39;: on the way from the pericentre to the apocentre the old tails following the orbit are dissolved and new ones pointing towards the MW are formed over a short timescale. We also find a tight linear relation between the velocity of stars in the tidal tails and their distance from the dwarf. Using mock data sets we demonstrate that if dwarf spheroidal (dSph) galaxies in the vicinity of the MW are tidally affected their kinematic samples are very likely contaminated by tidally stripped stars which tend to artificially inflate the measured velocity dispersion. The effect is stronger for dwarfs on their way from the peri- to the apocentre due to the formation of new tidal tails after pericentre. Realistic mass estimates of dSph galaxies thus require removal of these stars from kinematic samples.

preprint2009arXiv

The stellar structure and kinematics of dwarf spheroidal galaxies formed by tidal stirring

Using high-resolution N-body simulations we study the stellar properties of dwarf spheroidal galaxies resulting from the tidally induced morphological transformation of disky dwarfs on a cosmologically motivated eccentric orbit around the Milky Way. Dwarf galaxy models initially consist of an exponential stellar disk embedded in an extended spherical dark matter halo. Depending on the initial orientation of the disk with respect to the orbital plane, different final configurations are obtained. The least evolved dwarf is triaxial and retains a significant amount of rotation. The more evolved dwarfs are prolate spheroids with little rotation. We show that the final density distribution of stars can be approximated by a simple modification of the Plummer law. The kinematics of the dwarfs is significantly different depending on the line of sight which has important implications for mapping the observed stellar velocity dispersions of dwarfs to subhalo circular velocities. When the dwarfs are observed along the long axis, the measured velocity dispersion is higher and decreases faster with radius. In the case where rotation is significant, when viewed perpendicular to the long axis, the effect of minor axis rotation is detected, as expected for triaxial systems. We model the velocity dispersion profiles and rotation curves of the dwarfs by solving the Jeans equations for spherical and axisymmetric systems and adjusting different sets of free parameters. We find that the mass is typically overestimated when the dwarf is seen along the long axis and underestimated when the observation is along the short or intermediate axis. The effect of non-sphericity cannot however bias the inferred mass by more than 60 percent in either direction, even for the most strongly stripped dwarf which is close to disruption.

preprint2009arXiv

Tidal evolution of disky dwarf galaxies in the Milky Way potential: the formation of dwarf spheroidals

We conduct high-resolution collisionless N-body simulations to investigate the tidal evolution of dwarf galaxies on an eccentric orbit in the Milky Way (MW) potential. The dwarfs originally consist of a low surface brightness stellar disk embedded in a cosmologically motivated dark matter halo. During 10 Gyr of dynamical evolution and after 5 pericentre passages the dwarfs suffer substantial mass loss and their stellar component undergoes a major morphological transformation from a disk to a bar and finally to a spheroid. The bar is preserved for most of the time as the angular momentum is transferred outside the galaxy. A dwarf spheroidal (dSph) galaxy is formed via gradual shortening of the bar. This work thus provides a comprehensive quantitative explanation of a potentially crucial morphological transformation mechanism for dwarf galaxies that operates in groups as well as in clusters. We compare three cases with different initial inclinations of the disk and find that the evolution is fastest when the disk is coplanar with the orbit. Despite the strong tidal perturbations and mass loss the dwarfs remain dark matter dominated. For most of the time the 1D stellar velocity dispersion, σ, follows the maximum circular velocity, V_{\rm max}, and they are both good tracers of the bound mass. Specifically, we find that M_{\rm bound} \propto V_{\rm max}^{3.5} and V_{\rm max} \sim \sqrt{3} σin agreement with earlier studies based on pure dark matter simulations. The latter relation is based on directly measuring the stellar kinematics of the simulated dwarf and may thus be reliably used to map the observed stellar velocity dispersions of dSphs to halo circular velocities when addressing the missing satellites problem.

preprint2008arXiv

Cold Dark Matter Substructure and Galactic Disks

We perform a set of high-resolution, dissipationless N-body simulations to investigate the influence of cold dark matter (CDM) substructure on the dynamical evolution of thin galactic disks. Our method combines cosmological simulations of galaxy-sized CDM halos to derive the properties of substructure populations and controlled numerical experiments of consecutive subhalo impacts onto initially-thin, fully-formed disk galaxies. We demonstrate that close encounters between massive subhalos and galactic disks since z~1 should be common occurrences in LCDM models. In contrast, extremely few satellites in present-day CDM halos are likely to have a significant impact on the disk structure. One typical host halo merger history is used to seed controlled N-body experiments of subhalo-disk encounters. As a result of these accretion events, the disk thickens considerably at all radii with the disk scale height increasing in excess of a factor of 2 in the solar neighborhood. We show that interactions with the subhalo population produce a wealth of distinctive morphological signatures in the disk stars including: conspicuous flares; bars; low-lived, ring-like features in the outskirts; and low-density, filamentary structures above the disk plane. We compare a resulting dynamically-cold, ring-like feature in our simulations to the Monoceros ring stellar structure in the MW. The comparison shows quantitative agreement in both spatial distribution and kinematics, suggesting that such observed complex stellar components may arise naturally as disk stars are excited by encounters with subhalos. These findings highlight the significant role of CDM substructure in setting the structure of disk galaxies and driving galaxy evolution.

preprint2008arXiv

Cold Dark Matter Substructure and Galactic Disks I: Morphological Signatures of Hierarchical Satellite Accretion

(Abridged) We conduct a series of high-resolution, dissipationless N-body simulations to investigate the cumulative effect of substructure mergers onto thin disk galaxies in the context of the LCDM paradigm of structure formation. Our simulation campaign is based on a hybrid approach. Substructure properties are culled directly from cosmological simulations of galaxy-sized cold dark matter (CDM) halos. In contrast to what can be inferred from statistics of the present-day substructure populations, accretions of massive subhalos onto the central regions of host halos, where the galactic disk resides, since z~1 should be common occurrences. One host halo merger history is subsequently used to seed controlled numerical experiments of repeated satellite impacts on an initially-thin Milky Way-type disk galaxy. We show that these accretion events produce several distinctive observational signatures in the stellar disk including: a ring-like feature in the outskirts; a significant flare; a central bar; and faint filamentary structures that (spuriously) resemble tidal streams. The final distribution of disk stars exhibits a complex vertical structure that is well-described by a standard ``thin-thick&#39;&#39; disk decomposition. We conclude that satellite-disk encounters of the kind expected in LCDM models can induce morphological features in galactic disks that are similar to those being discovered in the Milky Way, M31, and in other disk galaxies. These results highlight the significant role of CDM substructure in setting the structure of disk galaxies and driving galaxy evolution. Upcoming galactic structure surveys and astrometric satellites may be able to distinguish between competing cosmological models by testing whether the detailed structure of galactic disks is as excited as predicted by the CDM paradigm.

preprint2007arXiv

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.

preprint2006arXiv

Multi-scale simulations of merging galaxies with supermassive black holes

We present the results of the first multi-scale N-Body+SPH simulations of merging galaxies containing central supermassive black holes (SMBHs) and having a spatial resolution of only a few parsecs. Strong gas inflows associated with equal-mass mergers produce non-axisymmetric nuclear disks with masses of order $10^9 M_{\odot}$, resolved by about $10^6$ SPH particles. Such disks have sizes of several hundred parsecs but most of their mass is concentrated within less than $50$ pc. We find that a close SMBH pair forms after the merger. The separation of the two SMBHs then shrinks further owing to dynamical friction against the predominantly gaseous background. The orbits of the SMBHs decay down to the minimum resolvable scale in a few million years for an ambient gas temperature and density typical of a region undergoing a starburst. These results suggest the initial conditions necessary for the eventual coalescence of the two holes arise naturally from the merging of two equal-mass galaxies whose structure and orbits are consistent with the predictions of the $Λ$CDM model. Our findings have important implications for planned gravitational wave detection experiments such as {\it LISA}.

preprint2006arXiv

On relaxation processes in collisionless mergers

We analyze N-body simulations of halo mergers to investigate the mechanisms responsible for driving mixing in phase-space and the evolution to dynamical equilibrium. We focus on mixing in energy and angular momentum and show that mixing occurs in step-like fashion following pericenter passages of the halos. This makes mixing during a merger unlike other well known mixing processes such as phase mixing and chaotic mixing whose rates scale with local dynamical time. We conclude that the mixing process that drives the system to equilibrium is primarily a response to energy and angular momentum redistribution that occurs due to impulsive tidal shocking and dynamical friction rather than a result of chaotic mixing in a continuously changing potential. We also analyze the merger remnants to determine the degree of mixing at various radii by monitoring changes in radius, energy and angular momentum of particles. We confirm previous findings that show that the majority of particles retain strong memory of their original kinetic energies and angular momenta but do experience changes in their potential energies owing to the tidal shocks they experience during pericenter passages. Finally, we show that a significant fraction of mass (~ 40%) in the merger remnant lies outside its formal virial radius and that this matter is ejected roughly uniformly from all radii outside the inner regions. This highlights the fact that mass, in its standard virial definition, is not additive in mergers. We discuss the implications of these results for our understanding of relaxation in collisionless dynamical systems.

preprint2005arXiv

Dark Matter Halos: Shapes, The Substructure Crisis, and Indirect Detection

In this proceeding, we briefly review three recent results. First, we show that halos formed in simulations with gas cooling are significantly rounder than halos formed in dissipationless $N$-body simulations. The increase in principle axis ratios is $δ(c/a) ~ 0.2 - 0.4$ in the inner halo and remains significant at large radii. Second, we discuss the CDM substructure crisis and demonstrate the sensitivity of the crisis to the spectrum of primordial density fluctuations on small scales. Third, we assess the ability of experiments like VERITAS and GLAST to detect $γ$-rays from neutralino dark matter annihilation in dark subhalos about the MW.

preprint2005arXiv

Fate of Supermassive Black Holes and Evolution of the M(BH)-sigma 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, we show that 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 M(BH)-sigma 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.