A conclusive model for the formation of dwarf spheroidal (dSph) galaxies still remains elusive. Owing to their proximity to the massive spirals Milky Way (MW) and M31, various environmental processes have been invoked to explain their origin. In this context, the tidal stirring model postulates that interactions with MW-sized hosts can transform rotationally supported dwarfs, resembling present-day dwarf irregular (dIrr) galaxies, into systems with the kinematic and structural properties of dSphs. Using N-body+SPH simulations, we investigate the dependence of this transformation mechanism on the gas fraction, fgas, in the disk of the progenitor dwarf. Our numerical experiments incorporate for the first time the combined effects of radiative cooling, ram-pressure stripping, star formation, supernova (SN) winds, and a cosmic UV background. For a given orbit inside the primary galaxy, rotationally supported dwarfs with gas fractions akin to those of observed dIrrs (fgas >= 0.5), demonstrate a substantially enhanced likelihood and efficiency of transformation into dSphs relative to their collisionless (fgas = 0) counterparts. We argue that the combination of ram-pressure stripping and SN winds causes the gas-rich dwarfs to respond more impulsively to tides, augmenting their transformation. When fgas >= 0.5, disky dwarfs on previously unfavorable low-eccentricity or large-pericenter orbits are still able to transform. On the widest orbits, the transformation is incomplete; the dwarfs retain significant rotational support, a relatively flat shape, and some gas, naturally resembling transition-type systems. We conclude that tidal stirring constitutes a prevalent evolutionary mechanism for shaping the structure of dwarf galaxies within the currently favored CDM cosmological paradigm.

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 < 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 "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.

We study the statistics and cosmic evolution of massive black hole seeds formed during major mergers of gas-rich late-type galaxies. Generalizing the results of the hydro-simulations from Mayer et al. 2010, we envision a scenario in which a supermassive star can form at the center of galaxies that just experienced a major merger owing to a multi-scale powerful gas inflow, provided that such galaxies live in haloes with masses above 10^{11} Msun, are gas-rich and disc-dominated, and do not already host a massive black hole. We assume that the ultimate collapse of the supermassive star leads to the rapid formation of a black hole of 10^5 Msun following a quasi-star stage. Using a model for galaxy formation applied to the outputs of the Millennium Simulation, we show that the conditions required for this massive black hole formation route to take place in the concordance LambdaCDM model are actually common at high redshift, and can be realized even at low redshift. Most major mergers above z~4 in haloes with mass > 10^{11} Msun can lead to the formation of a massive seed and, at z~2, the fraction of favourable mergers decreases to about half. Interestingly, we find that even in the local universe a fraction (~20%) of major mergers in massive haloes still satisfy the conditions for our massive black hole formation route. Those late events take place in galaxies with a markedly low clustering amplitude, that have lived in isolation for most of their life, and that are experiencing a major merger for the first time. We predict that massive black hole seeds from galaxy mergers can dominate the massive end of the mass function at high (z>4) and intermediate (z~2) redshifts relative to lighter seeds formed at higher redshift, for example, by the collapse of Pop III stars. Finally, a fraction of these massive seeds could lie, soon after formation, above the MBH-MBulge relation.

We model for the first time the complete orbital evolution of a pair of Supermassive Black Holes (SMBHs) in a 1:10 galaxy merger of two disk dominated gas-rich galaxies, from the stage prior to the formation of the binary up to the onset of gravitational wave emission when the binary separation has shrunk to 1 milli parsec. The high-resolution smoothed particle hydrodynamics (SPH) simulations used for the first phase of the evolution include star formation, accretion onto the SMBHs as well as feedback from supernovae explosions and radiative heating from the SMBHs themselves. Using the direct N-body code \phi-GPU we evolve the system further without including the effect of gas, which has been mostly consumed by star formation in the meantime. We start at the time when the separation between two SMBHs is ~ 700 pc and the two black holes are still embedded in their galaxy cusps. We use 3 million particles to study the formation and evolution of the SMBH binary till it becomes hard. After a hard binary is formed, we reduce (reselect) the particles to 1.15 million and follow the subsequent shrinking of the SMBH binary due to 3-body encounters with the stars. We find approximately constant hardening rates and that the SMBH binary rapidly develops a high eccentricity. Similar hardening rates and eccentricity values are reported in earlier studies of SMBH binary evolution in the merging of dissipation-less spherical galaxy models. The estimated coalescence time is ~ 2.9 Gyr, significantly smaller than a Hubble time. We discuss why this timescale should be regarded as an upper limit. Since 1:10 mergers are among the most common interaction events for galaxies at all cosmic epochs, we argue that several SMBH binaries should be detected with currently planned space-borne gravitational wave interferometers, whose sensitivity will be especially high for SMBHs in the mass range considered here.

Supermassive black holes (SMBHs) have been detected in the centers of most nearby massive galaxies. Galaxies today are the products of billions of years of galaxy mergers, but also billions of years of SMBH activity as active galactic nuclei (AGNs) that is connected to galaxy mergers. In this context, detection of AGN pairs should be relatively common. Observationally, however, dual AGN are scant, being just a few percent of all AGN. In this Letter we investigate the triggering of AGN activity in merging galaxies via a suite of high resolution hydrodynamical simulations. We follow the dynamics and accretion onto the SMBHs as they move from separations of tens of kiloparsecs to tens of parsecs. Our resolution, cooling and star formation implementation produce an inhomogeneous, multi-phase interstellar medium, allowing us to accurately trace star formation and accretion onto the SMBHs. We study the impact of gas content, morphology, and mass ratio, allowing us to study AGN activity and dynamics across a wide range of relevant conditions. We test when the two AGN are simultaneously detectable, for how long and at which separations. We find that strong dual AGN activity occurs during the late phases of the mergers, at small separations (<1-10 kpc) below the resolution limit of most surveys. Much of the SMBH accretion is not simultaneous, limiting the dual AGN fraction detectable through imaging and spectroscopy to a few percent, in agreement with observational samples.

Simulations of the formation of late-type spiral galaxies in a cold dark matter LCDM universe have traditionally failed to yield realistic candidates. Here we report a new cosmological N-body/SPH simulation of extreme dynamic range in which a close analog of a Milky Way disk galaxy arises naturally. Termed Eris, the simulation follows the assembly of a galaxy halo of mass Mvir=7.9x10^11 Msun with a total of N=18.6 million particles (gas + dark matter + stars) within the final virial radius, and a force resolution of 120 pc. It includes radiative cooling, heating from a cosmic UV field and supernova explosions, a star formation recipe based on a high gas density threshold (nSF=5 atoms cm^-3 rather than the canonical nSF=0.1 atoms cm^-3), and neglects AGN feedback. At the present epoch, the simulated galaxy has an extended rotationally-supported disk with a radial scale length Rd=2.5 kpc, a gently falling rotation curve with circular velocity at 2.2 disk scale lenghts of V2.2=214 km/s, a bulge-to-disk ratio B/D=0.35, and a baryonic mass fraction that is 30% below the cosmic value. The disk is thin, is forming stars in the region of the Sigma_SFR - Sigma_HI plane occupied by spiral galaxies, and falls on the photometric Tully-Fisher and the stellar mass-halo virial mass relations. Hot (T>3x10^5 K), X-ray luminous halo gas makes only 26% of the universal baryon fraction and follows a flattened density profile proportional to r^-1.13 out to r=100 kpc. Eris appears then to be the first cosmological hydrodynamic simulation in which the galaxy structural properties, the mass budget in the various components, and the scaling relations between mass and luminosity are all consistent with a host of observational constraints. (Abridged)

We perform a suite of high-resolution smoothed particle hydrodynamics simulations to investigate the orbital decay and mass evolution of massive black hole (MBH) pairs down to scales of ~30 pc during minor mergers of disk galaxies. Our simulation set includes star formation and accretion onto the MBHs, as well as feedback from both processes. We consider 1:10 merger events starting at z~3, with MBH masses in the sensitivity window of the Laser Interferometer Space Antenna, and we follow the coupling between the merger dynamics and the evolution of the MBH mass ratio until the satellite galaxy is tidally disrupted. While the more massive MBH accretes in most cases as if the galaxy were in isolation, the satellite MBH may undergo distinct episodes of enhanced accretion, owing to strong tidal torques acting on its host galaxy and to orbital circularization inside the disk of the primary galaxy. As a consequence, the initial 1:10 mass ratio of the MBHs changes by the time the satellite is disrupted. Depending on the initial fraction of cold gas in the galactic disks and the geometry of the encounter, the mass ratios of the MBH pairs at the time of satellite disruption can stay unchanged or become as large as 1:2. Remarkably, the efficiency of MBH orbital decay correlates with the final mass ratio of the pair itself: MBH pairs that increase significantly their mass ratio are also expected to inspiral more promptly down to nuclear-scale separations. These findings indicate that the mass ratios of MBH pairs in galactic nuclei do not necessarily trace the mass ratios of their merging host galaxies, but are determined by the complex interplay between gas accretion and merger dynamics.

The asymmetric emission of gravitational waves produced during the coalescence of a massive black hole (MBH) binary imparts a velocity "kick" to the system that can displace the hole from the center of its host. Here we study the trajectories and observability of MBHs recoiling in three (one major, two minor) gas-rich galaxy merger remnants that were previously simulated at high resolution, and in which the pairing of the MBHs had been shown to be successful. We run new simulations of MBHs recoiling in the major merger remnant with Mach numbers in the range 1<M<6, and use simulation data to construct a semi-analytical model for the orbital evolution of MBHs in gas-rich systems. We show that: 1) in major merger remnants the energy deposited by the moving hole into the rotationally supported, turbulent medium makes a negligible contribution to the thermodynamics of the gas. This contribution becomes significant in minor merger remnants, potentially allowing for an electromagnetic signature of MBH recoil; 2) in major merger remnants, the combination of both deeper central potential well and drag from high-density gas confines even MBHs with kick velocities as high as 1200 km/s within 1 kpc from the host's center; 3) kinematically offset nuclei may be observable for timescales of a few Myr in major merger remnants in the case of recoil velocities in the range 700-1,000 km/s; 4) in minor mergers remnants the effect of gas drag is weaker, and MBHs with recoil speeds in the range 300-600 km/s will wander through the host halo for longer timescales. When accounting for the probability distribution of kick velocities, however, we find that the likelihood of observing recoiling MBHs in gas-rich galaxy mergers is very low, typically below 10^-5 - 10^-6.

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.

(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.

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.

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.

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.

We report on key studies on the dynamics of black holes (BHs) in gas-rich galaxy mergers that underscore the vital role played by gas dissipation in promoting BH inspiral down to the smallest scales ever probed with use of high-resolution numerical simulations. In major mergers, the BHs sink rapidly under the action of gas-dynamical friction while orbiting inside the massive nuclear disc resulting from the merger. The BHs then bind and form a Keplerian binary on a scale of 5 pc. In minor mergers, BH pairing proceeds down to the minimum scale explored of 10-100 pc only when the gas fraction in the less massive galaxy is comparatively large to avoid its tidal and/or ram pressure disruption and the wandering of the light BH in the periphery of the main halo. Binary BHs enter the gravitational wave dominated inspiral only when their relative distance is typically of 0.001 pc. If the gas preserves the degree of dissipation expected in a star-burst environment, binary decay continues down to 0.1 pc, the smallest length-scale ever attained. Stalling versus hardening below 0.1 pc is still matter of deep investigations.

Giant low surface brightness galaxies (GLSBs) have flat discs extending up to ~100 kpc. Their formation is a puzzle for cosmological simulations in the cold dark matter scenario. We suggest that GLSBs might be the final product of the evolution of collisional ring galaxies. In fact, our simulations show that, approximately 0.5-1.5 Gyr after the collision which led to the formation of the ring galaxy, the ring keeps expanding and fades, while the disc becomes very large (~100 kpc) and flat. At this stage, our simulated galaxies match many properties of GLSBs (surface brightness profile, morphology, HI spectrum and rotation curve).

Disk galaxies with a spheroidal component are known to host Supermassive Black Holes (SMBHs) in their center. Unequal-mass galaxy mergers have been rarely studied despite the fact that they are the large majority of merging events by number and they are associated with the typical targets of gravitational wave experiments such as LISA. We performed N-body/SPH simulations of disk galaxy mergers with mass ratios 1:4 and 1:10 at redshifts z=0 and z=3. They have the highest resolution achieved so far for merging galaxies, and include star formation and supernova feedback. Gas dissipation is found to be necessary for the pairing of SMBHs in these minor mergers. Still, 1:10 mergers with gas allow an efficient pairing only at high z when orbital times are short enough compared to the Hubble time.

Chandra and XMM-Newton observations of the Cartwheel galaxy show ~17 bright X-ray sources (>~5x10^38 erg s^-1), all within the gas-rich outer ring. We explore the hypothesis that these X-ray sources are powered by intermediate-mass black holes (IMBHs) accreting gas or undergoing mass transfer from a stellar companion. To this purpose, we run N-body/SPH simulations of the galaxy interaction which might have led to the formation of Cartwheel, tracking the dynamical evolution of two different IMBH populations: halo and disc IMBHs. Halo IMBHs cannot account for the observed X-ray sources, as only a few of them cross the outer ring. Instead, more than half of the disc IMBHs are pulled in the outer ring as a consequence of the galaxy collision. However, also in the case of disc IMBHs, accretion from surrounding gas clouds cannot account for the high luminosities of the observed sources. Finally, more than 500 disc IMBHs are required to produce <~15 X-ray sources via mass transfer from very young stellar companions. Such number of IMBHs is very large and implies extreme assumptions. Thus, the hypothesis that all the observed X-ray sources in Cartwheel are associated with IMBHs is hardly consistent with our simulations, even if it is still possible that IMBHs account for the few (<~1-5) brightest ultraluminous X-ray sources (ULXs).

We present a study on the dynamics of massive BHs in galaxy mergers, obtained from a series of high-resolution N-Body/SPH simulations. The presence of a gaseous component is essential for the rapid formation of an eccentric (Keplerian) BH binary, that resides at the center of a massive (~10^9 Msun) turbulent nuclear disc. Using physically and/or numerically motivated recipes, we follow the accretion history of the BHs during the merger. The mass of the BHs increases as large central inflows of gas occur inside each galaxy, and their mass ratio varies with time. Given the encountered strong degeneracy between numerical resolution and physical assumptions, we suggest here three possible paths followed by the galaxies and the BHs during a merger in order to fulfill the M-sigma relation : Adjustment, Symbiosis, and BH Dominance. In an extremely high resolution run, we resolved the turbulent gas pattern down to parsec scales, and found that BH feedback is expected to be effective near the end of the merger. We then trace the BH binary orbit down to a scale of 0.1 pc modeling the nuclear disc as an equilibrium Mestel disc composed either of gas, gas and stars, or just stars. Under the action of dynamical friction against the rotating gaseous and/or stellar background the orbit circularizes. When this occurs, each BH is endowed with its own small-size (~0.01 pc) accretion disc comprising a few percent of the BH mass. Double AGN activity is expected to occur on an estimated timescale of ~10 Myrs, comparable to the inspiral time. The double nuclear point--like sources that may appear have typical separations of ~10 pc, and are likely to be embedded in the still ongoing starburst.