We present a revised version of Peters' (1964) time-scale for the gravitational-wave (GW) induced decay of two point masses. The new formula includes the effects of the first-order post-Newtonian perturbation and additionally provides a simple fit to account for the Newtonian self-consistent evolution of the eccentricity. The revised time-scale is found by multiplying Peters' estimate by two factors, $R(e_0)= 8^{1-\sqrt{1-e_0}}$ and $Q_{\rm f}(p_0) = \exp \left(2.5 (r_{\rm S}/p_0) \right)$, where $e_0$ and $p_0$ are the initial eccentricity and periapsis, respectively, and $r_{\rm S}$ the Schwarzschild radius of the system. Their use can correct errors of a factor of 1-10 that arise from using the original Peters' formula. We apply the revised time-scales to a set of typical sources for existing ground-based laser interferometers and for the future Laser Interferometer Space Antenna (LISA), at the onset of their GW driven decay. We argue that our more accurate model for the orbital evolution will affect current event- and detection-rate estimates for mergers of compact object binaries, with stronger deviations for eccentric LISA sources, such as extreme and intermediate mass-ratio inspirals. We propose the correction factors $R$ and $Q_{\rm f}$ as a simple prescription to quantify decay time-scales more accurately in future population synthesis models. We also suggest that the corrected time-scale may be used as a computationally efficient alternative to numerical integration in other applications that include the modelling of radiation reaction for eccentric sources.

The forthcoming Laser Interferometer Space Antenna (LISA) will probe the population of coalescing massive black hole (MBH) binaries up to the onset of structure formation. Here we simulate the galactic-scale pairing of $\sim10^6 M_\odot$ MBHs in a typical, non-clumpy main-sequence galaxy embedded in a cosmological environment at $z = 7-6$. In order to increase our statistical sample, we adopt a strategy that allows us to follow the evolution of six secondary MBHs concomitantly. We find that the magnitude of the dynamical-friction induced torques is significantly smaller than that of the large-scale, stochastic gravitational torques arising from the perturbed and morphologically evolving galactic disc, suggesting that the standard dynamical friction treatment is inadequate for realistic galaxies at high redshift. The dynamical evolution of MBHs is very stochastic, and a variation in the initial orbital phase can lead to a drastically different time-scale for the inspiral. Most remarkably, the development of a galactic bar in the host system either significantly accelerates the inspiral by dragging a secondary MBH into the centre, or ultimately hinders the orbital decay by scattering the MBH in the galaxy outskirts. The latter occurs more rarely, suggesting that galactic bars overall promote MBH inspiral and binary coalescence. The orbital decay time can be an order of magnitude shorter than what would be predicted relying on dynamical friction alone. The stochasticity, and the important role of global torques, have crucial implications for the rates of MBH coalescences in the early Universe: both have to be accounted for when making predictions for the upcoming LISA observatory.

We report on high-resolution simulations that explore the orbital decay of massive black hole (MBH) pairs with masses between $10^5$ and $10^7 M_{\odot}$ embedded in a circumnuclear gas disk (CND). An adiabatic equation of state is adopted, with a range of adiabatic indices, which maintains a smooth flow. Mergers between MBHs in this mass range would be detectable by the upcoming Laser Inteferometer Space Antenna (LISA). The orbital evolution is followed from the CND scale ($100$~pc) down to separations of $0.1$--$0.01$~pc at which a circumbinary disk (CBD) could form. The decay is erratic and strongly dependent on the gas flow within the disk, that ultimately determines the net torques experienced by the sinking MBH. Overall, we can identify three different evolutionary stages: (i) an initially slow decay that leads to no significant change in the orbital angular momentum, resulting in some circularization; (ii) a fast migration phase in which the orbital angular momentum decreases rapidly; and (iii) a final, very slow decay phase, in which orbital angular momentum can even increase, and a CBD can form. The fast migration phase owes to disk-driven torques originating primarily from the co-orbital region of the secondary MBH, at a distance of 1--3 Hill radii. We find strong analogies with fast Type III migration for massive planets in protoplanetary disks. The CBD forms only when the decay rate becomes small enough to allow it enough time to carve a cavity around the primary MBH, at scales $\lesssim 1$~pc; when this happens, the MBH separation nearly stalls in our higher-resolution run. We suggest an empirically modified gap-opening criterion that takes into account such timescale effects as well as other deviations from standard assumptions made in the literature. [Abriged]

We review the current status of knowledge concerning the early phases of star formation during cosmic dawn. This includes the first generations of stars forming in the lowest mass dark matter halos in which cooling and condensation of gas with primordial composition is possible at very high redshift ($z > 20$), namely metal-free Population III stars, and the first generation of massive black holes forming at such early epochs, the so-called black hole seeds. The formation of black hole seeds as end states of the collapse of Population III stars, or via direct collapse scenarios, is discussed. In particular, special emphasis is given to the physics of supermassive stars as potential precursors of direct collapse black holes, in light of recent results of stellar evolution models, and of numerical simulations of the early stages of galaxy formation. Furthermore, we discuss the role of the cosmic radiation produced by the early generation of stars and black holes at high redshift in the process of reionization.

In the early stages of a protoplanetary disk, when its mass is a significant fraction of its star's, turbulence generated by gravitational instability (GI) should feature significantly in the disk's evolution. At the same time, the disk may be sufficiently ionised for magnetic fields to play some role in the dynamics. Though usually neglected, the impact of magnetism on the GI may be critical, with consequences for several processes: the efficiency of accretion, spiral structure formation, fragmentation, and the dynamics of solids. In this paper, we report on global three-dimensional magnetohydrodynamical simulations of a self-gravitating protoplanetary disk using the meshless finite mass (MFM) Lagrangian technique. We confirm that GI spiral waves trigger a dynamo that amplifies an initial magnetic field to nearly thermal amplitudes (plasma beta < 10), an order of magnitude greater than that generated by the magneto-rotational instability alone. We also determine the dynamo's nonlinear back reaction on the gravitoturbulent flow: the saturated state is substantially hotter, with an associated larger Toomre parameter and weaker, more 'flocculent' spirals. But perhaps of greater import is the dynamo's boosting of accretion via a significant Maxwell stress; mass accretion is enhanced by factors of several relative to either pure GI or pure MRI. Our simulations use ideal MHD, an admittedly poor approximation in protoplanetary disks, and thus future studies should explore the full gamut of non-ideal MHD. In preparation for that, we exhibit a small number of Ohmic runs that reveal that the dynamo, if anything, is stronger in a non-ideal environment. This work confirms that magnetic fields are a potentially critical ingredient in gravitoturbulent young disks, possibly controlling their evolution, especially via their enhancement of (potentially episodic) accretion.

Soon after the recent first ever detection of gravitational waves from merging black holes it has been suggested that their origin is primordial. Appealingly, a sufficient number of primordial black holes (PBHs) could also partially or entirely constitute the dark matter (DM) in the Universe. However, recent studies on PBHs in ultra-faint dwarf galaxies (UFDGs) suggest that they would dynamically heat up the stellar component due to two-body relaxation processes. From the comparison with the observed stellar velocity dispersions and the stellar half-light radii it was claimed that only PBHs with masses $\lesssim10\,M_\odot$ can significantly contribute to the DM. In this work, we improve the latter constraints by considering the largest observational sample of UFDGs and by allowing the PBH masses to follow an extended (log-normal) distribution. By means of collisional Fokker-Planck simulations, we explore a wide parameter space of UFDGs containing PBHs. The analysis of the half-light radii and velocity dispersions resulting from the simulations leads to three general findings that exclude PBHs with masses $\sim\mathcal{O}(1$-$100)\,M_\odot$ from constituting all of the DM: (i) We identify a critical sub-sample of UFDGs that only allows for $\sim\mathcal{O}(1)\,M_\odot$ PBH masses; (ii) for any PBH mass, there is an UFDG in our sample that disfavours it; (iii) the spatial extensions of a majority of simulated UFDGs containing PBHs are too large to match the observed.

The quest for binary and dual supermassive black holes (SMBHs) at the dawn of the multi-messenger era is compelling. Detecting dual active galactic nuclei (AGN) -- active SMBHs at projected separations larger than several parsecs -- and binary AGN -- probing the scale where SMBHs are bound in a Keplerian binary -- is an observational challenge. The study of AGN pairs (either dual or binary) also represents an overarching theoretical problem in cosmology and astrophysics. The AGN triggering calls for detailed knowledge of the hydrodynamical conditions of gas in the imminent surroundings of the SMBHs and, at the same time, their duality calls for detailed knowledge on how galaxies assemble through major and minor mergers and grow fed by matter along the filaments of the cosmic web. This review describes the techniques used across the electromagnetic spectrum to detect dual and binary AGN candidates and proposes new avenues for their search. The current observational status is compared with the state-of-the-art numerical simulations and models for formation of dual and binary AGN. Binary SMBHs are among the loudest sources of gravitational waves (GWs) in the Universe. The search for a background of GWs at nHz frequencies from inspiralling SMBHs at low redshifts, and the direct detection of signals from their coalescence by the Laser Interferometer Space Antenna in the next decade, make this a theme of major interest for multi-messenger astrophysics. This review discusses the future facilities and observational strategies that are likely to significantly advance this fascinating field.

We investigated the formation and evolution of satellite systems in a cold, extended circumplanetary disc around a 10 $M_{\rm{Jupiter}}$ gas giant which was generated by gravitational instability. We use a population synthesis approach, in which we seed the disc with satellite embryos and let them migrate, accrete mass, collide, with subsequent generations being created until the disc is dissipated. In each run we choose randomly the dust-to-gas ratio, dispersion- and refilling time-scales, the number of embryos and their starting locations in a plausible range. The disc structure is the result of a 3D global disc SPH simulation \citep{Szulagyi16b}. We also investigate the effect of the planet's semi-major axis on the resulting satellite systems, taking 50 AU as the nominal case. In the nominal case we find that most satellites are close in mass to the Galilean ones, with a maximum at around 3 $M_{\rm{Earth}}$ and form typically on time-scales comparable to the dispersion time-scale. We also find that about 10 $M_{\rm{Earth}}$ worth of satellites migrate into the planet, polluting it with metals. The influence of a different planet semi-major axis is mainly felt through the change in disc size. We find that for the discs closer to the star, the satellites are lighter, formation time-scales are longer, and more moons are lost into the planet, as a smaller disc means easier migration for the massive satellites. The probability of detecting satellites like the ones in this work is very low ($\leq$ 3%) even with upcoming powerful telescopes like E-ELT.

Chondrules are crystallised droplets of silicates formed by rapid heating to high temperatures (> 1800 K) of solid precursors followed by hours or days of cooling. Dating of chondrules is consistent with the formation timescale of Jupiter in the core-accretion model (1-4 Myrs). Here we investigate if the shocks generated by a massive planet could generate flash heating episodes necessary to form chondrules using high resolution 2D simulations with the multi-fluid code ROSSBI. We use different radiative cooling prescriptions, i.e. different cooling rates and models, and vary planet mass, orbit and disk models. For long disk cooling rates (> 1000 orbits) and a massive protoplanet (> 0.75 jupiter masses), we obtain hot enough temperatures for chondrule formation, while using more realistic thermodynamics is not successful in the Minimum Mass Solar Nebula (MMSN) model. However, sufficient flash heating is achieved by a Jupiter mass planet in a 5 times more massive disk, which is a conceivable scenario for the young solar nebula and exoplanetary systems. Interestingly, a gap-forming massive planet triggers vortices, which can trap dust, i.e. chondrule precursors, and generates a high pressure environment that is consistent with cosmochemical evidence from chondrules. A massive gas giant can thus, in principle, both stimulate the concentration of chondrule precursors in the vicinity of the shocking regions, and flash-heat them via the shocks.

Dark matter haloes play a fundamental role in cosmological structure formation. The most common approach to model their assembly mechanisms is through N-body simulations. In this work we present an innovative pathway to predict dark matter halo formation from the initial density field using a Deep Learning algorithm. We implement and train a Deep Convolutional Neural Network (DCNN) to solve the task of retrieving Lagrangian patches from which dark matter halos will condense. The volumetric multi-label classification task is turned into a regression problem by means of the euclidean distance transformation. The network is complemented by an adaptive version of the watershed algorithm to form the entire protohalo identification pipeline. We show that splitting the segmentation problem into two distinct sub-tasks allows for training smaller and faster networks, while the predictive power of the pipeline remains the same. The model is trained on synthetic data derived from a single full N-body simulation and achieves deviations of ~10% when reconstructing the dark matter halo mass function at z=0. This approach represents a promising framework for learning highly non-linear relations in the primordial density field. As a practical application, our method can be used to produce mock dark matter halo catalogues directly from the initial conditions of N-body simulations.

Major mergers of gas-rich galaxies provide promising conditions for the formation of supermassive black holes (SMBHs; $\gtrsim10^5$ M$_\odot$) by direct collapse because they can trigger mass inflows as high as $10^4-10^5$ M$_\odot$ yr$^{-1}$ on sub-parsec scales. However, the channel of SMBH formation in this case, either dark collapse (direct collapse without prior stellar phase) or supermassive star (SMS; $\gtrsim10^4$ M$_\odot$), remains unknown. Here, we investigate the limit in accretion rate up to which stars can maintain hydrostatic equilibrium. We compute hydrostatic models of SMSs accreting at $1-1000$ M$_\odot$ yr$^{-1}$, and estimate the departures from equilibrium a posteriori by taking into account the finite speed of sound. We find that stars accreting above the atomic cooling limit ($\gtrsim10$ M$_\odot$ yr$^{-1}$) can only maintain hydrostatic equilibrium once they are supermassive. In this case, they evolve adiabatically with a hylotropic structure, that is, entropy is locally conserved and scales with the square root of the mass coordinate. Our results imply that stars can only become supermassive by accretion at the rates of atomically cooled haloes ($\sim0.1-10$ M$_\odot$ yr$^{-1}$). Once they are supermassive, larger rates are possible.

The cold molecular gas in contemporary galaxies is structured in discrete cloud complexes. These giant molecular clouds (GMCs), with $10^4$-$10^7$ solar masses and radii of 5-100 parsecs, are the seeds of star formation. Highlighting the molecular gas structure at such small scales in distant galaxies is observationally challenging. Only a handful of molecular clouds were reported in two extreme submillimetre galaxies at high redshift. Here we search for GMCs in a typical Milky Way progenitor at z = 1.036. Using the Atacama Large Millimeter/submillimeter Array (ALMA), we mapped the CO(4-3) emission of this gravitationally lensed galaxy at high resolution, reading down to 30 parsecs, which is comparable to the resolution of CO observations of nearby galaxies. We identify 17 molecular clouds, characterized by masses, surface densities and supersonic turbulence all of which are 10-100 times higher than present-day analogues. These properties question the universality of GMCs and suggest that GMCs inherit their properties from ambient interstellar medium. The measured cloud gas masses are similar to the masses of stellar clumps seen in the galaxy in comparable numbers. This corroborates the formation of molecular clouds by fragmentation of distant turbulent galactic gas disks, which then turn into stellar clumps ubiquitously observed in galaxies at cosmic noon.

We present a novel study of dust-vortex evolution in global two-fluid disk simulations to find out if evolution toward high dust-to-gas ratios can occur in a regime of well-coupled grains with low Stokes numbers ($St=10^{-3}-{4\times 10^{-2}}$). We design a new implicit scheme in the code RoSSBi, to overcome the short timesteps occurring for small grain sizes. We discover that the linear capture phase occurs self-similarly for all grain sizes, with an intrinsic timescale (characterizing the vortex lifetime) scaling as $1/St$. After vortex dissipation, the formation of a {{global active dust ring}} is a generic outcome confirming our previous results obtained for larger grains. We propose a scenario in which, irrespective of grain size, multiple pathways can lead to local dust-to-gas ratios of order unity and above on relatively short timescales, $< 10^5$ yr, in the presence of a vortex, even with $St=10^{-3}$. When $St>10^{-2}$, the vortex is quickly dissipated by two-fluid instabilities, and large dust density enhancements form in the global dust ring. When $St<10^{-2}$, the vortex is resistant to destabilization. As a result, dust concentrations occur locally due to turbulence developing inside the vortex. Whatever the Stokes number, dust-to-gas ratios in the range $1-10$, a necessary condition to trigger a subsequent streaming instability, or even a direct gravitational instability of the dust clumps, appears to be an inevitable outcome. Although quantitative connections with other instabilities still need to be made, we argue that our results support a new scenario of vortex-driven planetesimal formation.

We propose a space-based interferometer surveying the gravitational wave (GW) sky in the milli-Hz to $\mu$-Hz frequency range. By the 2040s', the $\mu$-Hz frequency band, bracketed in between the Laser Interferometer Space Antenna (LISA) and pulsar timing arrays, will constitute the largest gap in the coverage of the astrophysically relevant GW spectrum. Yet many outstanding questions related to astrophysics and cosmology are best answered by GW observations in this band. We show that a $\mu$-Hz GW detector will be a truly overarching observatory for the scientific community at large, greatly extending the potential of LISA. Conceived to detect massive black hole binaries from their early inspiral with high signal-to-noise ratio, and low-frequency stellar binaries in the Galaxy, this instrument will be a cornerstone for multimessenger astronomy from the solar neighbourhood to the high-redshift Universe.

Bars are a key factor in the long-term evolution of spiral galaxies, in their unique role in redistributing angular momentum and transporting gas and stars on large scales. The Eris-suite simulations are cosmological zoom-in, N-body, smoothed-particle hydrodynamic simulations built to follow the formation and evolution of a Milky Way-sized galaxy across the build-up of the large scale structure. Here we analyse and describe the outcome of two particular simulations taken from the Eris suite - ErisBH and Eris2k - which mainly differ in the prescriptions employed for gas cooling, star formation, and feedback from supernovae and black holes. Our study shows that the enhanced effective feedback in Eris2k, due to the collective effect of the different micro-physics implementations, results in a galaxy which is less massive than its ErisBH counterpart till z~1. However, when the stellar content is large enough so that global dynamical instabilities can be triggered, the galaxy in Eris2k develops a stronger and more extended bar with respect to ErisBH. We demonstrate that he structural properties and time evolution of the two bars are very different. Our results highlight the importance of accurate sub-grid prescriptions in cosmological zoom-in simulations of the process of galaxy formation and evolution, and the possible use of a statistical sample of barred galaxies to assess the strength of the stellar feedback.

In the local universe, black holes of $10^{5-6}$ M$_{\odot}$ are hosted in galaxies displaying a variety of stellar profiles and morphologies. These black holes are the anticipated targets of LISA, the Laser Interferometer Space Antenna that will detect the low-frequency gravitational-wave signal emitted by binary black holes in this mass interval. In this paper, we infer upper limits on the lifetime of binary black holes of $10^{5-6}$ M$_{\odot}$ and up to $10^8$ M$_{\odot}$, forming in galaxy mergers, exploring two underlying stellar density profiles, by Dehnen and by Prugniel & Simien, and by exploiting local scaling relations between the mass of the black holes and several quantities of their hosts. We focus on the phase of the dynamical evolution when the binary is transitioning from the hardening phase ruled by the interaction with single stars to the phase driven by the emission of gravitational waves. We find that different stellar profiles predict very distinct trends with binary mass, with lifetimes ranging between fractions of a Gyr to more than 10 Gyr, and with a spread of about one order of magnitude, given by the uncertainties in the observed correlations, which are larger in the low-mass tail of the observed black hole population.

In recent years, the discovery of massive quasars at z~7 has provided a striking challenge to our understanding of the origin and growth of supermassive black holes in the early Universe. Mounting observational and theoretical evidence indicates the viability of massive seeds, formed by the collapse of supermassive stars, as a progenitor model for such early, massive accreting black holes. Although considerable progress has been made in our theoretical understanding, many questions remain regarding how (and how often) such objects may form, how they live and die, and how next generation observatories may yield new insight into the origin of these primordial titans. This review focusses on our present understanding of this remarkable formation scenario, based on discussions held at the Monash Prato Centre from November 20--24, 2017, during the workshop "Titans of the Early Universe: The Origin of the First Supermassive Black Holes."

The advent of new deep+wide photometric lensing surveys will open up the possibility of direct measurements of the dark matter halos of dwarf galaxies. The HSC wide survey will be the first with the statistical capability of measuring the lensing signal with high signal-to-noise at log(M*)=8. At this same mass scale, LSST will have the most overall constraining power with a predicted signal-to-noise for the galaxy-galaxy lensing signal around dwarfs of SN=200. WFIRST and LSST will have the greatest potential to push below the log(M*) = 7 mass scale thanks to the depth of their imaging data. Studies of the dark matter halos of dwarf galaxies at z=0.1 with gravitational lensing are soon within reach. However, further work will be required to develop optimized strategies for extracting dwarfs samples from these surveys, determining redshifts, and accurately measuring lensing on small radial scales. Dwarf lensing will be a new and powerful tool to constrain the halo masses and inner density slopes of dwarf galaxies and to distinguish between baryonic feedback and modified dark matter scenarios.

The giant impact hypothesis for Moon formation successfully explains the dynamic properties of the Earth-Moon system but remains challenged by the similarity of isotopic fingerprints of the terrestrial and lunar mantles. Moreover, recent geochemical evidence suggests that the Earth's mantle preserves ancient (or "primordial") heterogeneity that predates the Moon-forming giant impact. Using a new hydrodynamical method, we here show that Moon-forming giant impacts lead to a stratified starting condition for the evolution of the terrestrial mantle. The upper layer of the Earth is compositionally similar to the disk, out of which the Moon evolves, whereas the lower layer preserves proto-Earth characteristics. As long as this predicted compositional stratification can at least partially be preserved over the subsequent billions of years of Earth mantle convection, the compositional similarity between the Moon and the accessible Earth's mantle is a natural outcome of realistic and high-probability Moon-forming impact scenarios. The preservation of primordial heterogeneity in the modern Earth not only reconciles geochemical constraints but is also consistent with recent geophysical observations. Furthermore, for significant preservation of a proto-Earth reservoir, the bulk composition of the Earth-Moon system may be systematically shifted towards chondritic values.

Observational evidence has been mounting for the existence of intermediate mass black holes (IMBHs, 10^2-10^5 Msun), but observing them at all, much less constraining their masses, is very challenging. In one theorized formation channel, IMBHs are the seeds for supermassive black holes in the early universe. As a result, IMBHs are predicted to exist in the local universe in dwarf galaxies, as well as wandering in more massive galaxy halos. However, these environments are not conducive to the accretion events or dynamical signatures that allow us to detect IMBHs. The Laser Interferometer Space Antenna (LISA) will demystify IMBHs by detecting the mergers of these objects out to extremely high redshifts, while measuring their masses with extremely high precision. These observations of merging IMBHs will allow us to constrain the formation mechanism and subsequent evolution of massive black holes, from the 'dark ages' to the present day, and reveal the role that IMBHs play in hierarchical galaxy evolution.

The magneto-rotational instability (MRI) is one of the most important processes in sufficiently ionized astrophysical disks. Grid-based simulations, especially those using the local shearing box approximation, provide a powerful tool to study the ensuing nonlinear turbulence. On the other hand, while meshless methods have been widely used in both cosmology, galactic dynamics, and planet formation they have not been fully deployed on the MRI problem. We present local unstratified and vertically stratified MRI simulations with two meshless MHD schemes: a recent implementation of SPH MHD (Price2012), and a MFM MHD scheme with a constrained gradient divergence cleaning scheme, as implemented in the GIZMO code \citep{Hopkins2017}. Concerning variants of the SPH hydro force formulation we consider both the "vanilla" SPH and the PSPH variant included in GIZMO. We find, as expected, that the numerical noise inherent in these schemes affects turbulence significantly. A high order kernel, free of the pairing instability, is necessary. Both schemes can adequately simulate MRI turbulence in unstratified shearing boxes with net vertical flux. The turbulence, however, dies out in zero-net-flux unstratified boxes, probably due to excessive and numerical dissipation. In zero-net-flux vertically stratified simulations, MFM can reproduce the MRI dynamo and its characteristic butterfly diagram for several tens of orbits before ultimately decaying. In contrast, extremely strong toroidal fields, as opposed to sustained turbulence, develop in equivalent simulations using SPH MHD. This unphysical state in SPH MHD is likely caused by a combination of excessive artificial viscosity, numerical resistivity, and the relatively large residual errors in the divergence of the magnetic field remaining even after cleaning procedures are applied.

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

Gas clumps formed within massive gravitationally unstable circumstellar discs are potential seeds of gas giant planets, brown dwarfs and companion stars. Simulations show that competition between three processes -- migration, gas accretion and tidal disruption -- establishes what grows from a given seed. Here we investigate the robustness of numerical modelling of clump migration and accretion with the codes PHANTOM, GADGET, SPHINX, SEREN, GIZMO-MFM, SPHNG and FARGO. The test problem comprises a clump embedded in a massive disc at an initial separation of 120 AU. There is a general qualitative agreement between the codes, but the quantitative agreement in the planet migration rate ranges from $\sim 10$% to $\sim 50$%, depending on the numerical setup. We find that the artificial viscosity treatment and the sink particle prescription may account for much of the differences between the codes. In order to understand the wider implications of our work, we also attempt to reproduce the planet evolution tracks from our hydrodynamical simulations with prescriptions from three previous population synthesis studies. We find that the disagreement amongst the population synthesis models is far greater than that between our hydrodynamical simulations. The results of our code comparison project are therefore encouraging in that uncertainties in the given problem are probably dominated by the physics not yet included in the codes rather than by how hydrodynamics is modelled in them.

Giant impacts (GIs) are common in the late stage of planet formation. The Smoothed Particle Hydrodynamics (SPH) method is widely used for simulating the outcome of such violent collisions, one prominent example being the formation of the Moon. However, a decade of numerical studies in various areas of computational astrophysics has shown that the standard formulation of SPH suffers from several shortcomings such as artificial surface tension and its tendency to promptly damp turbulent motions on scales much larger than the physical dissipation scale, both resulting in the suppression of mixing. In order to quantify how severe these limitations are when modeling GIs we carried out a comparison of simulations with identical initial conditions performed with the standard SPH as well as with the novel Lagrangian Meshless Finite Mass (MFM) method in the GIZMO code. We confirm the lack of mixing between the impactor and target when SPH is employed, while MFM is capable of driving vigorous sub-sonic turbulence and leads to significant mixing between the two bodies. Modern SPH variants with artificial conductivity, a different formulation of the hydro force or reduced artificial viscosity, do not improve mixing as significantly. Angular momentum is conserved similarly well in both methods, but MFM does not suffer from spurious transport induced by artificial viscosity, resulting in a slightly higher angular momentum of the proto-lunar disk. Furthermore, SPH initial conditions exhibit an unphysical density discontinuity at the core-mantle boundary which is easily removed in MFM.

Satellites of giant planets thought to form in gaseous circumplanetary disks (CPDs) during the late planet-formation phase, but it was unknown so far whether smaller mass planets, such as the ice giants could form such disks, thus moons there. We combined radiative hydrodynamical simulations with satellite population synthesis to investigate the question in the case of Uranus and Neptune. For both ice giants we found that a gaseous CPD is created at the end of their formation. The population synthesis confirmed that Uranian-like, icy, prograde satellite-system could form in these CPDs within a couple of $10^5$ years. This means that Neptune could have a Uranian-like moon-system originally that was wiped away by the capture of Triton. Furthermore, the current moons of Uranus can be reproduced by our model without the need for planet-planet impact to create a debris disk for the moons to grow. These results highlight that even ice giants -- that among the most common mass-category of exoplanets -- can also form satellites, opening a way to a potentially much larger population of exomoons than previously thought.

We report the results of high-resolution subpixel imaging of the hard continuum and Fe K{\alpha} line of the Compton Thick (CT) Active Galactic Nucleus (AGN) ESO 428-G014, observed with Chandra ACIS. While the 3-4 keV emission is dominated by an extended component, a single nuclear point source is prominent in the 4-6 keV range. Instead, two peaks of similar intensity, separated by ~36 pc in projection on the plane of the sky are detected in the Fe K{\alpha} emission. The SE knot could be marginally associated with the heavily obscured hard continuum source. We discuss four possible interpretations of the nuclear morphology. (1) Given the bolometric luminosity and likely black hole (BH) mass of ESO 428-G014, we may be imaging two clumps of the CT obscuring torus in the Fe K{\alpha} line. (2) The Fe K{\alpha} knots may be connected with the fluorescent emission from the dusty bicone, or (3) with the light echo of a nuclear outburst. (4) We also explore the less likely possibility that we may be detecting the rare signature of merging nuclei. Considering the large-scale kpc-size extent of the hard continuum and Fe K{\alpha} emission (Papers I and II), we conclude that the AGN in ESO 428-G014 has been active for at least 104 yrs. Comparison with the models of Czerny et al (2009) suggests high accretion rates during this activity.

Theoretical models for the expected merger rates of intermediate-mass black holes (IMBHs) are vital for planned gravitational-wave detection experiments such as the Laser Interferometer Space Antenna (LISA). Using collisionless $N$-body simulations of dwarf galaxy (DG) mergers, we examine how the orbital decay of IMBHs and the efficiency of IMBH binary formation depend on the central dark matter (DM) density profile of the merging DGs. Specifically, we explore various asymptotic inner slopes $\gamma$ of the DG's DM density distribution, ranging from steep cusps ($\gamma=1$) to shallower density profiles ($\gamma<1$), motivated by well-known baryonic-feedback effects as well as by DM models that differ from cold DM at the scales of DGs. We find that the inner DM slope is crucial for the formation (or lack thereof) of an IMBH binary; only mergers between DGs with cuspy DM profiles ($\gamma=1$) are favourable to forming a hard IMBH binary, whereas when $\gamma<1$ the IMBHs stall at a separation of 50-100 pc. Consequently, the rate of LISA signals from IMBH coalescence will be determined by the fraction of DGs with a cuspy DM profile. Conversely, the LISA event rates at IMBH mass scales offer in principle a novel way to place constraints on the inner structure of DM halos in DGs and address the core-cusp controversy. We also show that, with spatial resolutions of $\sim$0.1 kpc, as often adopted in cosmological simulations, all IMBHs stall, independent of $\gamma$. This suggests caution in employing cosmological simulations of galaxy formation to study BH dynamics in DGs.

We use high-resolution cosmological hydrodynamical simulations of Milky Way-sized galaxies with varying supernovae feedback strengths and merger histories to investigate the formation of their gaseous halos and especially their hot ($>10^6$~K) X-ray luminous coronae. Our simulations predict the presence of significant hot gas in the halos as early as $z=3-4$, well before the halos ought to be able to sustain hot mode accretion in the conventional picture. The nascent coronae grow inside-out and initially do so primarily as a result of outflows from the central galaxies powered by merger-induced shock heating and strong supernovae feedback, both of which are elemental features of today's successful galaxy formation models. Furthermore, the outflows and the forming coronae also accelerate the transition from cold to hot mode accretion by contributing to the conditions for sustaining stable accretion shocks. They also disrupt the filamentary streams funneling cold gas onto the central galaxies by causing their mouths to fray into a broad delta, detach from the galaxies, and be pushed away to larger radii. And even though at early times the filaments repeatedly re-form, the hot gas and the outflows act to weaken the filaments and accelerate their ultimate disruption. Although galactic outflows are generally thought of as ejective feedback, we find that their action on the filaments suggests a preventive role as well.

The standard model for giant planet formation is based on the accretion of solids by a growing planetary embryo, followed by rapid gas accretion once the planet exceeds a so-called critical mass. The dominant size of the accreted solids (cm-size particles named pebbles or km to hundred km-size bodies named planetesimals) is, however, unknown. Recently, high-precision measurements of isotopes in meteorites provided evidence for the existence of two reservoirs in the early Solar System. These reservoirs remained separated from ~1 until ~ 3 Myr after the beginning of the Solar System's formation. This separation is interpreted as resulting from Jupiter growing and becoming a barrier for material transport. In this framework, Jupiter reached ~20 Earth masses within ~1 Myr and slowly grew to ~50 Earth masses in the subsequent 2 Myr before reaching its present-day mass. The evidence that Jupiter slowed down its growth after reaching 20 Earth masses for at least 2 Myr is puzzling because a planet of this mass is expected to trigger fast runaway gas accretion. Here, we use theoretical models to describe the conditions allowing for such a slow accretion and show that Jupiter grew in three distinct phases. First, rapid pebble accretion brought the major part of Jupiter's core mass. Second, slow planetesimal accretion provided the energy required to hinder runaway gas accretion during 2 Myr. Third, runaway gas accretion proceeded. Both pebbles and planetesimals therefore have an important role in Jupiter's formation.

We present the first 2D hydrodynamical finite volume simulations in which dust is fully coupled with the gas, including its back-reaction onto it, and at the same time the dust size is evolving according to coagulation and fragmentation based on a sub-grid model. The aim of this analysis is to present the differences occurring when dust evolution is included relative to simulations with fixed dust size, with and without an embedded Jupiter-mass planet that triggers gap formation. We use the two-fluid polar Godunov-type code RoSSBi developed by Surville et al. 2016 combined with a new local sub-grid method for dust evolution based on the model by Birnstiel et al. 2012. We find striking differences between simulations with variable and fixed dust sizes. The timescales for dust depletion differ significantly and yield a completely different evolution of the dust surface density. In general sharp features such as pile-ups of dust in the inner disk and near gap edges, when a massive planet is present, become much weaker. This has important implications on the interpretation of observed substructure in disks, suggesting that the presence of a massive planet does not necessarily cause sharp gaps and rings in the dust component. Also, particles with different dust sizes show a different distribution, pointing to the importance of multi-wavelength synthetic observations in order to compare with observations by ALMA and other instruments. We also find that simulations adopting fixed intermediate particle sizes, in the range $10^{-2} - 10^{-1}$ cm, best approximate the surface density evolution seen in simulations with dust evolution.

The satellites of Jupiter are thought to form in a circumplanetary disc. Here we address their formation and orbital evolution with a population synthesis approach, by varying the dust-to-gas ratio, the disc dispersal timescale and the dust refilling timescale. The circumplanetary disc initial conditions (density and temperature) are directly drawn from the results of 3D radiative hydrodynamical simulations. The disc evolution is taken into account within the population synthesis. The satellitesimals were assumed to grow via streaming instability. We find that the moons form fast, often within $10^4$ years, due to the short orbital timescales in the circumplanetary disc. They form in sequence, and many are lost into the planet due to fast type I migration, polluting Jupiter's envelope with typically 15 Earth-masses of metals. The last generation of moons can form very late in the evolution of the giant planet, when the disc has already lost more than the 99% of its mass. The late circumplanetary disc is cold enough to sustain water ice, hence not surprisingly the 85% of the moon population has icy composition. The distribution of the satellite-masses is peaking slightly above Galilean masses, up until a few Earth-masses, in a regime which is observable with the current instrumentation around Jupiter-analog exoplanets orbiting sufficiently close to their host stars. We also find that systems with Galilean-like masses occur in 20% of the cases and they are more likely when discs have long dispersion timescales and high dust-to-gas ratios.

$\Lambda$-Warm Dark Matter (WDM) has been proposed as alternative scenario to $\Lambda$ cold dark matter (CDM), motivated by discrepancies at the scale of dwarf galaxies, with less small-scale power and realized by collisionless particles with energies in the range $1-3$ keV. We present a new approach to constrain the viability of such WDM models using star formation histories of the dwarf spheroidal galaxies (dSphs) in the Local Group. We compare their high time-resolution star formation histories (SFHs) obtained with HST-based color magnitude diagrams with the range of possible collapse redshifts of their dark matter halos expected in CDM and in different WDM scenarios. The collapse redshift is inferred after determining a plausible infall mass of the subhalo. This is based on the current mass of individual dwarf inferred from stellar kinematics combined with results of cosmological simulations providing information on the subhalo evolution. Since WDM subhalos close to the filtering mass scale form significantly later than CDM, we show that they are in the first place difficult to reconcile with a truncation of star formation occurring as early as $z\geq 3$. The Ultra-Faint Dwarfs (UFDs) provide the most stringent constraints. Using 6 UFDs with the best determination of the SFHs, we show that we can exclude a 1 keV warm particle to a 2-$\sigma$ confidence interval consistently with other methods reported in the literature. For some objects the $2$ keV model is also excluded. We discuss the various caveats of the method, most notably the low number of dwarfs with accurately determined star formation histories and the uncertainties in the determination of the infall mass of the subhalos. Our preliminary analysis serves as a pathfinder for future investigations that will combine upcoming accurate SFHs for more local dSphs with direct analysis of WDM cosmological simulations with baryons.

The Laser Interferometer Space Antenna (LISA) will detect gravitational-wave (GW) signals from merging supermassive black holes (BHs) with masses below $10^7$~M$_{\odot}$. It is thus of paramount importance to understand the orbital dynamics of these relatively light central BHs, which typically reside in disc-dominated galaxies, in order to produce reliable forecasts of merger rates. To this aim, realistic simulations probing BH dynamics in unequal-mass disc galaxy mergers, into and beyond the binary hardening stage, are performed by combining smooth particle hydrodynamics and direct $N$-body codes. The structural properties and orbits of the galaxies are chosen to be consistent with the results of galaxy formation simulations. Stellar and dark matter distributions are triaxial down to the central 100 pc of merger remnant. In all cases, a BH binary forms and hardens on time-scales of at most 100~Myr, coalescing on another few hundred Myr time-scale, depending on the characteristic density and orbital eccentricity. Overall, the sinking of the BH binary takes no more than $\sim$0.5~Gyr after the merger of the two galaxies is completed, but can be much faster for very plunging orbits. Comparing with previous numerical simulations following the decay of BHs in massive early-type galaxies at $z \sim 3$, we confirm that the characteristic density is the most crucial parameter determining the overall BH merging time-scale, despite the structural diversity of the host galaxies. Our results lay down the basis for robust forecasts of LISA event rates in the case of merging BHs.

The emergence of a large-scale stellar bar is one of the most striking features in disc galaxies. By means of state-of-the-art cosmological zoom-in simulations, we study the formation and evolution of bars in Milky Way-like galaxies in a fully cosmological context, including the physics of gas dissipation, star formation, and supernova feedback. Our goal is to characterise the actual trigger of the non-axisymmetric perturbation that leads to the strong bar observable in the simulations at z=0, discriminating between an internal/secular versus an external/tidal origin. To this aim, we run a suite of cosmological zoom-in simulations altering the original history of galaxy-satellite interactions at a time when the main galaxy, though already bar-unstable, does not feature any non-axisymmetric structure yet. We find that the main effect of a late minor merger and of a close fly-by is to delay the time of bar formation and those two dynamical events are not directly responsible for the development of the bar and do not alter significantly its global properties (e.g. its final extension). We conclude that, once the disc has grown to a mass large enough to sustain global non-axisymmetric modes, then bar formation is inevitable.

It has been proposed that close interactions with satellite galaxies can significantly perturb the morphology of the main galaxy. However, the dynamics of an already formed bar following the interaction with the external environment has not been studied in detail in a fully cosmological context. In this work, analysing the cosmological zoom-in simulation Eris2k, we study the effects that a very unequal-mass flyby crossing the stellar disc has on the stability of the pre-existing bar. We characterize the evolution of the bar strength and length showing that the perturbation exerted by the flyby shuffles the orbits of stars for less than one Gyr. After this time, the bar shows a remarkable resilience, reforming with properties comparable to those it had before the interaction. Our work shows that close unequal-mass encounters, the most frequent interactions occurring during the evolution of cosmic structures, have (i) an overall minor impact on the global evolution of the bar in the long term, still (ii) the effect is destructive and (iii) a very weak interaction is sufficient to dismantle a strong bar leading to its "apparent death". As a consequence, due to the non-negligible duration of the bar-less period, a fraction of observed spiral galaxies classified as non-barred could be prone to bar formation.

We review and discuss theoretical studies addressing the possibility of gas accretion onto black holes occurring at rates exceeding the Eddington limit. Our focus is on the applications to the growth of black hole seeds at high redshift. We first present the general notion of Super-Eddington accretion, and then summarize the different models and numerical simulations developed to study such regime. We consider optically thick flows in accretion disks as well as in spherically symmetric envelopes, and devote particular attention to the widely adopted model based on the SLIM disk solution. While attractive for its simplicity, the SLIM disk solution is challenged by the latest generation of three-dimensional radiation (magneto)-hydrodynamical simulations, in which radiative losses can be an order of magnitude higher, and the mechanisms of radiation transport is more complex than straight advection as it takes place in a complex turbulent regime. We then discuss the gas supply rate to the sub-pc scale accretion disk or envelope from larger scales, revisiting gas inflow rates in protogalaxies under various conditions. We conclude that in the dense gaseous nuclei of high-z galaxies the conditions necessary for the onset of Super Eddington accretion regimes, such as a high optical depth and high gas supply rates from large scales, should be naturally met. Feedback from the growing BH seed should not alter significantly such conditions according to the results of radiation magneto-hydrodynamical simulations of super-critical flows in accretion disks. Furthermore, based on the required nuclear gas inflow rates and the tendency of stellar feedback to remove efficiently gas in low mass halos, we argue that super-critical accretion will be more easily achieved in relatively sizable halos, with virial masses $M_{\rm vir} > 10^{10}$ M$_{\odot}$, which become more common at $z < 15$.

According to star formation histories (SFHs), Local Group dwarf galaxies can be broadly classified in two types: those forming most of their stars before $z=2$ (${\it fast}$) and those with more extended SFHs (${\it slow}$). The most precise SFHs are usually derived from deep but not very spatially extended photometric data; this might alter the ratio of old to young stars when age gradients are present. Here we correct for this effect and derive the mass formed in stars by $z=2$ for a sample of 16 Local Group dwarf galaxies. We explore early differences between ${\it fast}$ and ${\it slow}$ dwarfs, and evaluate the impact of internal feedback by supernovae (SN) on the baryonic and dark matter (DM) component of the dwarfs. ${\it Fast}$ dwarfs assembled more stellar mass at early times and have larger amounts of DM within the half-light radius than ${\it slow}$ dwarfs. By imposing that ${\it slow}$ dwarfs cannot have lost their gas by $z=2$, we constrain the maximum coupling efficiency of SN feedback to the gas and to the DM to be $\sim$10%. We find that internal feedback alone appears insufficient to quench the SFH of ${\it fast}$ dwarfs by gas deprivation, in particular for the fainter systems. Nonetheless, SN feedback can core the DM halo density profiles relatively easily, producing cores of the sizes of the half-light radius in ${\it fast}$ dwarfs by $z=2$ with very low efficiencies. Amongst the "classical" Milky Way satellites, we predict that the smallest cores should be found in Draco and Ursa Minor, while Sculptor and Fornax should host the largest ones.

In this paper we review a new scenario for the formation of massive black hole seeds that relies on multi-scale gas inflows initiated by the merger of massive gas-rich galaxies at $z > 6$, where gas has already achieved solar composition. Hydrodynamical simulations undertaken to explore our scenario show that supermassive, gravitationally bound gaseous disks, weighing a billion solar masses and of a few pc in size, form in the nuclei of merger remnants in less than $10^5$ yr. These could later produce a supermassive protostar or supermassive star at their center via various mechanisms. Moreover, we present a new analytical model, based on angular momentum transport in mass-loaded gravitoturbulent disks. This naturally predicts that a nuclear disk accreting at rates exceeding $1000 M_{\odot}$/yr, as seen in the simulations, is stable against fragmentation irrespective of its metallicity. This is at variance with conventional direct collapse scenarios, which require the suppression of gas cooling in metal-free protogalaxies for gas collapse to take place. Such high accretion rates reflect the high free-fall velocities in massive halos appearing at $z < 10$, and occur naturally as a result of the efficient angular momentum loss provided by mergers. We discuss the implications of our scenario on the observed population of high-z quasars and on its evolution to lower redshifts using a semi-analytical galaxy formation model. Finally, we consider the intriguing possibility that the secondary gas inflows in the unstable disks might drive gas to collapse into a supermassive black hole directly via the General Relativistic radial instability. Such {\it dark collapse} route could generate gravitational wave emission detectable via the future Laser Interferometer Space Antenna (LISA). [Abridged]

Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt-Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly-formed stellar clump mass functions show more significant variation (difference by up to a factor of ~3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low density region, and between more diffusive and less diffusive schemes in the high density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes.

We present a semi-analytical population synthesis model of protoplanetary clumps formed by disk instability at radial distances of 80 - 120 AU. Various clump density profiles, initial mass functions, protoplanetary disk models, stellar masses, and gap opening criteria are considered. When we use more realistic gap opening criteria, we find that gaps open only rarely which strongly affects clump survival rates and their physical properties (mass, radius and radial distance). The inferred surviving population is then shifted towards less massive clumps at smaller radial distances. We also find that populations of surviving clumps are very sensitive to the model assumptions and used parameters. Depending on the chosen parameters, the protoplanets occupy a mass range between 0.01 and 16 $M_{\text{J}}$ and may either orbit close to the central star or as far out as 75 AU, with a sweet spot at 10-30 AU for the massive ones. However, in all the cases we consider we find that massive giant planets at very large radial distances are rare in qualitative agreement with current direct imaging surveys. We conclude that caution should be taken in deriving population synthesis models as well as when comparing the models' results with observations.

One of the most computationally demanding aspects of the hydrodynamical modelling of Astrophysical phenomena is the transport of energy by radiation or relativistic particles. Physical processes involving energy transport are ubiquitous and of capital importance in many scenarios ranging from planet formation to cosmic structure evolution, including explosive events like core collapse supernova or gamma-ray bursts. Moreover, the ability to model and hence understand these processes has often been limited by the approximations and incompleteness in the treatment of radiation and relativistic particles. The DIAPHANE project has focused in developing a portable and scalable library that handles the transport of radiation and particles (in particular neutrinos) independently of the underlying hydrodynamic code. In this work, we present the computational framework and the functionalities of the first version of the DIAPHANE library, which has been successfully ported to three different smoothed-particle hydrodynamic codes, GADGET2, GASOLINE and SPHYNX. We also present validation of different modules solving the equations of radiation and neutrino transport using different numerical schemes.

Giant stellar clumps are ubiquitous in high-redshift galaxies. They are thought to play an important role in the build-up of galactic bulges and as diagnostics of star formation feedback in galactic discs. Hubble Space Telescope (HST) blank field imaging surveys have estimated that these clumps have masses up to 10$^{9.5}$ M$_{\odot}$ and linear sizes larger than ~1 kpc. Recently, gravitational lensing has also been used to get higher spatial resolution. However, both recent lensed observations and models suggest that the clumps properties may be overestimated by the limited resolution of standard imaging techniques. A definitive proof of this observational bias is nevertheless still missing. Here we investigate directly the effect of resolution on clump properties by analysing multiple gravitationally-lensed images of the same galaxy at different spatial resolutions, down to 30 pc. We show that the typical mass and size of giant clumps, generally observed at $\sim$1 kpc resolution in high-redshift galaxies, are systematically overestimated. The high spatial resolution data, only enabled by strong gravitational lensing using currently available facilities, support smaller scales of clump formation by fragmentation of the galactic gas disk via gravitational instabilities.

We carry out simulations of gravitationally unstable disks using smoothed particle hydrodynamics(SPH) and the novel Lagrangian meshless finite mass (MFM) scheme in the GIZMO code (Hopkins 2015). Our aim is to understand the cause of the non-convergence of the cooling boundary for fragmentation reported in the literature. We run SPH simulations with two different artificial viscosity implementations, and compare them with MFM, which does not employ any artificial viscosity. With MFM we demonstrate convergence of the critical cooling time scale for fragmentation at \beta_{crit} =3.. Non-convergence persists in SPH codes, although it is significantly mitigated with schemes having reduced artificial viscosity such as inviscid SPH (ISPH) (Cullen & Dehnen 2010). We show how the non-convergence problem is caused by artificial fragmentation triggered by excessive dissipation of angular momentum in domains with large velocity derivatives. With increased resolution such domains become more prominent. Vorticity lags behind density due to numerical viscous dissipation in these regions, promoting collapse with longer cooling times. Such effect is shown to be dominant over the competing tendency of artificial viscosity to diminish with increasing resolution. When the initial conditions are first relaxed for several orbits, the flow is more regular, with lower shear and vorticity in non-axisymmetric regions, aiding convergence. Yet MFM is the only method that converges exactly. Our findings are of general interest as numerical dissipation via artificial viscosity or advection errors can also occur in grid-based codes. Indeed for the FARGO code values of \beta_{crit} significantly higher than our converged estimate have been reported in the literature. Finally, we discuss implications for giant planet formation via disk instability.

We present mock observations of forming planets with ALMA. The possible detections of circumplanetary disks (CPDs) were investigated around planets of Saturn, 1, 3, 5, and 10 Jupiter-masses that are placed at 5.2 AU from their star. The radiative, three dimensional hydrodynamic simulations were then post-processed with RADMC3D and the ALMA Observation Simulator. We found that even though the CPDs are too small to be resolved, they are hot due to the accreting planet in the optically thick limit, therefore the best chance to detect them with continuum observations in this case is at the shortest ALMA wavelengths, such as Band 9 (440 microns). Similar fluxes were found in the case of Saturn and Jupiter-mass planets, as for the 10 $\mathrm{M_{Jup}}$ gas-giant, due to temperature weighted optical depth effects: when no deep gap is carved, the planet region is blanketed by the optically thick circumstellar disk leading to a less efficient cooling there. A test was made for a 52 AU orbital separation, showed that optically thin CPDs are also detectable in band 7 but they need longer integration times ($>$5hrs). Comparing the gap profiles of the same simulation at various ALMA bands and the hydro simulation confirmed that they change significantly, first because the gap is wider at longer wavelengths due to decreasing optical depth; second, the beam convolution makes the gap shallower and at least 25% narrower. Therefore, caution has to be made when estimating planet masses based on ALMA continuum observations of gaps.

Using the Eris zoom-in cosmological simulation of assembly of a Milky Way analog, we study chemical enrichment of stars due to accretion of metal-enriched gas from the interstellar medium during the Galaxy's development. We consider metal-poor and old stars in both Galactic halo and bulge and make use of stellar orbits, gas density and metallicity distributions in Eris. Assuming spherically symmetric Bondi-Hoyle accretion, we find that halo and bulge stars accrete metals at the rate of about 10^-24 solar mass per year and 10^-22 solar mass per year, respectively, at redshifts z < 3, but this accretion rate increases hundred-fold to about 10^-20 solar mass per year at higher redshifts due to increased gas density. Bulge and halo stars accrete similar amounts of metals at high redshifts as kinematically distinct bulge and halo are not yet developed at these redshifts and both sets of stars encounter similar metal distribution in the ISM on average. Accretion alone can enrich main-sequence stars up to [Fe/H] -2 in extreme cases.Median enrichment level due to accretion in these stars is about [Fe/H]~-6 to -5.Because accretion mostly takes place at high redshifts, it is alpha-enriched to [alpha/Fe]~0.5. We find that accretive metal enrichment is significant enough to affect the predicted metallicity distribution function of halo stars at [Fe/H] < -5.This suggests that attempts to infer the natal chemical environment of the most metal-poor stars from their observed enrichment today can be hindered due to metal accretion. Peculiar enrichment patterns such as those predicted to arise from pair-instability supernovae could help in disentangling natal and accreted metal content of stars.

We investigate the simultaneous triggering of active galactic nuclei (AGN) in merging galaxies, using a large suite of high-resolution hydrodynamical simulations. We compute dual-AGN observability time-scales using bolometric, X-ray, and Eddington-ratio thresholds, confirming that dual activity from supermassive black holes (BHs) is generally higher at late pericentric passages, before a merger remnant has formed, especially at high luminosities. For typical minor and major mergers, dual activity lasts ~20-70 and ~100-160 Myr, respectively. We also explore the effects of X-ray obscuration from gas, finding that the dual-AGN time decreases at most by a factor of ~2, and of contamination from star formation. Using projected separations and velocity differences rather than three-dimensional quantities can decrease the dual-AGN time-scales by up to ~4, and we apply filters which mimic current observational-resolution limitations. In agreement with observations, we find that, for a sample of major and minor mergers hosting at least one AGN, the fraction harbouring dual AGN is ~20-30 and ~1-10 per cent, respectively. We quantify the effects of merger mass ratio (0.1 to 1), geometry (coplanar, prograde, retrograde, and inclined), disc gas fraction, and BH properties, finding that the mass ratio is the most important factor, with the difference between minor and major mergers varying between factors of a few to orders of magnitude, depending on the luminosity and filter used. We also find that a shallow imaging survey will require very high angular resolution, whereas a deep imaging survey will be less resolution-dependent.

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 present results from the "Ponos" simulation suite on the early evolution of a massive, $M_{\rm vir}(z=0)=1.2\times 10^{13}$ M$_{\odot}$ galaxy. At $z\gtrsim6$, before feedback from a central supermassive black hole becomes dominant, the main galaxy has a stellar mass $\sim 2\times 10^{9}$ M$_{\odot}$ and a star formation rate $\sim 20$ M$_{\odot}$ yr$^{-1}$. The galaxy sits near the expected main sequence of star-forming galaxies at those redshifts, and resembles moderately star-forming systems observed at $z>5$. The high specific star formation rate results in vigorous heating and stirring of the gas by supernovae feedback, and the galaxy develops a thick and turbulent disc, with gas velocity dispersion $\sim 40$ km s$^{-1}$, rotation to dispersion ratio $\sim 2$, and with a significant amount of gas at $\sim 10^5$ K. The Toomre parameter always exceeds the critical value for gravito-turbulence, $Q\sim 1.5-2$, mainly due to the contribution of warm/hot gas inside the disc. Without feedback, a nearly gravito-turbulent regime establishes with similar gas velocity dispersion and lower $Q$. We propose that the "hot and turbulent" disc regime seen in our simulations, unlike the "cold and turbulent" gravito-turbulent regime of massive clumpy disc galaxies at $z\sim 1-2$, is a fundamental characterisation of main sequence galaxies at $z\gtrsim 6$, as they can sustain star formation rates comparable to those of low-mass starbursts at $z=0$. This results in no sustained coherent gas inflows through the disc, and in fluctuating and anisotropic mass transport, possibly postponing the assembly of the bulge and causing the initial feeding of the central black hole to be highly intermittent.

We discuss the latest results of numerical simulations following the orbital decay of massive black hole pairs in galaxy mergers. We highlight important differences between gas-poor and gas-rich hosts, and between orbital evolution taking place at high redshift as opposed to low redshift. Two effects have a huge impact and are rather novel in the context of massive black hole binaries. The first is the increase in characteristic density of galactic nuclei of merger remnants as galaxies are more compact at high redshift due to the way dark halo collapse depends on redshift. This leads naturally to hardening timescales due to 3-body encounters that should decrease by two orders of magnitude up to $z=4$. It explains naturally the short binary coalescence timescale, $\sim 10$ Myr, found in novel cosmological simulations that follow binary evolution from galactic to milliparsec scales. The second one is the inhomogeneity of the interstellar medium in massive gas-rich disks at high redshift. In the latter star forming clumps 1-2 orders of magnitude more massive than local Giant Molecular Clouds (GMCs) can scatter massive black holes out of the disk plane via gravitational perturbations and direct encounters. This renders the character of orbital decay inherently stochastic, often increasing orbital decay timescales by as much as a Gyr. At low redshift a similar regime is present at scales of $1-10$ pc inside Circumnuclear Gas Disks (CNDs). In CNDs only massive black holes with masses below $10^7 M_{\odot}$ can be significantly perturbed. They decay to sub-pc separations in up to $\sim 10^8$ yr rather than the in just a few million years as in a smooth CND. Finally implications for building robust forecasts of LISA event rates are discussed

We study the orbital decay of a pair of massive black holes (BHs) with masses 5 * 10^5 and 10^7 M_sun, using hydrodynamical simulations of circumnuclear disks (CNDs) with the alternating presence of sub-grid physics such as radiative cooling, star formation, supernova feedback, BH accretion and feedback. In the absence of such processes, the orbit of the secondary BH decays over timescales of ~10 Myr to the center of the CND, where the primary BH resides. When strong dissipation operates in CNDs, fragmentation into massive objects the size of giant molecular clouds and with densities in the range 10^4 - 10^7 amu / cm^3 occurs, causing stochastic torques and hits that can eject the secondary BH from the midplane. Outside the plane, the low-density medium provides only weak drag, and the BH return is governed by inefficient dynamical friction. In rare cases, clump-BH interactions can lead to a faster decay. Feedback processes lead to outflows, but do not change significantly the overall density of the CND midplane. However, with a spherically distributed BH feedback a hot bubble is generated behind the secondary, which almost shuts off dynamical friction, a phenomenon we dub "wake evacuation", leading to delays in the decay of possibly ~0.3 Gyr. We discuss the non-trivial implications on the discovery space of the eLISA telescope. Our results suggest the largest uncertainty in predicting BH merger rates lies in the potentially wide variety of galaxy host systems, with different degrees of gas dissipation and heating, yielding decay timescales from ~10 to ~300 Myr.

We analyse stellar masses of clumps drawn from a compilation of star-forming galaxies at 1.1<z<3.6. Comparing clumps selected in different ways, and in lensed or blank field galaxies, we examine the effects of spatial resolution and sensitivity on the inferred stellar masses. Large differences are found, with median stellar masses ranging from ~10^9 Msun for clumps in the often-referenced field galaxies to ~10^7 Msun for fainter clumps selected in deep-field or lensed galaxies. We argue that the clump masses, observed in non-lensed galaxies with a limited spatial resolution of ~1 kpc, are artificially increased due to the clustering of clumps of smaller mass. Furthermore, we show that the sensitivity threshold used for the clump selection affects the inferred masses even more strongly than resolution, biasing clumps at the low mass end. Both improved spatial resolution and sensitivity appear to shift the clump stellar mass distribution to lower masses, qualitatively in agreement with clump masses found in recent high-resolution simulations of disk fragmentation. We discuss the nature of the most massive clumps, and we conclude that it is currently not possible to properly establish a meaningful clump stellar mass distribution from observations and to infer the existence and value of a characteristic clump mass scale.

We investigate the angular momentum evolution of four disk galaxies residing in Milky Way-sized halos formed in cosmological zoom-in simulations with various sub-grid physics and merging histories. We decompose these galaxies kinematically and photometrically, into their disk and bulge components. The simulated galaxies and their components lie on the observed sequences in the $j_*$--$M_*$ diagram relating the specific angular momentum and mass of the stellar component. We find that galaxies in low-density environments follow the relation $j_* \propto M_*^{\alpha}$ past major mergers, with $\alpha \sim 0.6$ in the case of strong feedback, when bulge-to-disk ratios are relatively constant, and $\alpha \sim 1.4$ in the other cases, when secular processes operate on shorter timescales. We compute the retention factors (i.e. the ratio of the specific angular momenta of stars and dark matter) for both disks and bulges and show that they vary relatively slowly after averaging over numerous but brief fluctuations. For disks, the retention factors are usually close to unity, while for bulges, they are a few times smaller. Our simulations therefore indicate that galaxies and their halos grow in a quasi-homologous way.

High-resolution simulations of star-forming massive galactic discs have shown that clumps form with a characteristic baryonic mass in the range $10^7-10^8~M_{\odot}$, with a small tail exceeding $10^9~M_{\odot}$ produced by clump-clump mergers. This is in contrast with the observed kpc-size clumps with masses up to $10^{10}~M_{\odot}$ in high-redshift star-forming galaxies. In this paper we show that the comparison between simulated and observed star-forming clumps is hindered by limited observational spatial resolution and sensitivity. We post-process high-resolution hydrodynamical simulations of clumpy discs using accurate radiative transfer to model the effect of ionizing radiation from young stars and to compute H$\alpha$ emission maps. By comparing the intrinsic clump size and mass distributions with those inferred from convolving the H$\alpha$ maps with different gaussian apertures, we mimick the typical resolution used in observations. We found that with 100 pc resolution, mock observations can recover the intrinsic clump radii and stellar masses, in agreement with those found by lensing observations. Instead, using a 1 kpc resolution smears out individual clumps, resulting in their apparent merging. This causes significant overestimations of the clump radii and, therefore, masses derived using methods that use their observed sizes. We show that limited sensitivity can also force observations to significantly overestimate the clump masses. We conclude that a significant fraction of giant clumps detected in the observations may result from artificially inflated radii and masses, and that $\approx 100$ pc spatial resolution is required to capture correctly the physical characteristics of star-forming clumps if they are coherent structures produced by disc fragmentation.

Circumplanetary disks can be found around forming giant planets, regardless of whether core accretion or gravitational instability built the planet. We carried out state-of-the-art hydrodynamical simulations of the circumplanetary disks for both formation scenarios, using as similar initial conditions as possible to unveil possible intrinsic differences in the circumplanetary disk mass and temperature between the two formation mechanisms. We found that the circumplanetary disks mass linearly scales with the circumstellar disk mass. Therefore, in an equally massive protoplanetary disk, the circumplanetary disks formed in the disk instability model can be only a factor of eight more massive than their core-accretion counterparts. On the other hand, the bulk circumplanetary disk temperature differs by more than an order of magnitude between the two cases. The subdisks around planets formed by gravitational instability have a characteristic temperature below 100 K, while the core accretion circumplanetary disks are hot, with temperatures even greater than 1000 K when embedded in massive, optically thick protoplanetary disks. We explain how this difference can be understood as the natural result of the different formation mechanisms. We argue that the different temperatures should persist up to the point when a full-fledged gas giant forms via disk instability, hence our result provides a convenient criteria for observations to distinguish between the two main formation scenarios by measuring the bulk temperature in the planet vicinity.

Massive gas-rich galaxy discs at $z \sim 1-3$ host massive star-forming clumps with typical baryonic masses in the range $10^7-10^8$~M$_{\odot}$ which can affect the orbital decay and concurrent growth of supermassive black hole (BH) pairs. Using a set of high-resolution simulations of isolated clumpy galaxies hosting a pair of unequal-mass BHs, we study the interaction between massive clumps and a BH pair at kpc scales, during the early phase of the orbital decay. We find that both the interaction with massive clumps and the heating of the cold gas layer of the disc by BH feedback tend to delay significantly the orbital decay of the secondary, which in many cases is ejected and then hovers for a whole Gyr around a separation of 1--2 kpc. In the envelope, dynamical friction is weak and there is no contribution of disc torques: these lead to the fastest decay once the orbit of the secondary BH has circularised in the disc midplane. In runs with larger eccentricities the delay is stronger, although there are some exceptions. We also show that, even in discs with very sporadic transient clump formation, a strong spiral pattern affects the decay time-scale for BHs on eccentric orbits. We conclude that, contrary to previous belief, a gas-rich background is not necessarily conducive to a fast BH decay and binary formation, which prompts more extensive investigations aimed at calibrating event-rate forecasts for ongoing and future gravitational-wave searches, such as with Pulsar Timing Arrays and the future evolved Laser Interferometer Space Antenna.

Supermassive black holes (SMBHs) are ubiquitous in galaxies with a sizable mass. It is expected that a pair of SMBHs originally in the nuclei of two merging galaxies would form a binary and eventually coalesce via a burst of gravitational waves. So far theoretical models and simulations have been unable to predict directly the SMBH merger timescale from ab-initio galaxy formation theory, focusing only on limited phases of the orbital decay of SMBHs under idealized conditions of the galaxy hosts. The predicted SMBH merger timescales are long, of order Gyrs, which could be problematic for future gravitational wave searches. Here we present the first multi-scale $\Lambda$CDM cosmological simulation that follows the orbital decay of a pair of SMBHs in a merger of two typical massive galaxies at $z\sim3$, all the way to the final coalescence driven by gravitational wave emission. The two SMBHs, with masses $\sim10^{8}$ M$_{\odot}$, settle quickly in the nucleus of the merger remnant. The remnant is triaxial and extremely dense due to the dissipative nature of the merger and the intrinsic compactness of galaxies at high redshift. Such properties naturally allow a very efficient hardening of the SMBH binary. The SMBH merger occurs in only $\sim10$ Myr after the galactic cores have merged, which is two orders of magnitude smaller than the Hubble time.

We analyse the output of the hi-res cosmological zoom-in simulation ErisBH to study self-consistently the formation of a strong stellar bar in a Milky Way-type galaxy and its effect on the galactic structure, on the central gas distribution and on star formation. The simulation includes radiative cooling, star formation, SN feedback and a central massive black hole which is undergoing gas accretion and is heating the surroundings via thermal AGN feedback. A large central region in the ErisBH disk becomes bar-unstable after z~1.4, but a clear bar-like structure starts to grow significantly only after z~0.4, possibly triggered by the interaction with a massive satellite. At z~0.1 the bar reaches its maximum radial extent of l~2.2 kpc. As the bar grows, it becomes prone to buckling instability, which we quantify based on the anisotropy of the stellar velocity dispersion. The actual buckling event is observable at z~0.1, resulting in the formation of a boxy-peanut bulge clearly discernible in the edge-on view of the galaxy at z=0. The bar in ErisBH does not dissolve during the formation of the bulge but remains strongly non-axisymmetric down to the resolution limit of ~100 pc at z=0. During its early growth, the bar exerts a strong torque on the gas within its extent and drives gas inflows that enhance the nuclear star formation on sub-kpc scales. Later on the infalling gas is nearly all consumed into stars and, to a lesser extent, accreted onto the central black hole, leaving behind a gas-depleted region within the central ~2 kpc. Observations would more likely identify a prominent, large-scale bar at the stage when the galactic central region has already been quenched. Bar-driven quenching may play an important role in disk-dominated galaxies at all redshift. [Abridged]

We study the dependence of fragmentation in massive gas-rich galaxy disks at $z > 1$ on feedback model and hydrodynamical method, employing the GASOLINE2 SPH code and the lagrangian mesh-less code GIZMO in finite mass mode. We compare non-cosmological galaxy disk runs with standard blastwave supernovae (SN)feedback, which introduces delayed cooling in order to drive winds, and runs with the new superbubble SN feedback, which produces winds naturally by modelling the detailed physics of SN-driven bubbles and leads to efficient self-regulation of star formation. We find that, with blastwave feedback, massive star forming clumps form in comparable number and with very similar masses in GASOLINE2 and GIZMO. The typical masses are in the range $10^7-10^8 M_{\odot}$, lower than in most previous works, while giant clumps with masses above $10^9 M_{\odot}$ are exceedingly rare. With superbubble feedback, instead, massive bound star forming clumps do not form because galaxies never undergo a phase of violent disk instability. Only sporadic, unbound star forming overdensities lasting only a few tens of Myr can arise that are triggered by perturbations of massive satellite companions. We conclude that there is a severe tension between explaining massive star forming clumps observed at $z > 1$ primarily as the result of disk fragmentation driven by gravitational instability and the prevailing view of feedback-regulated galaxy formation. The link between disk stability and star formation efficiency should thus be regarded as a key testing ground for galaxy formation theory.

The Initial Star formation and Lifetimes of Andromeda Satellites (ISLAndS) project uses Hubble Space Telescope imaging to study a representative sample of six Andromeda dSph satellite companion galaxies. The main goal of the program is to determine whether the star formation histories (SFHs) of the Andromeda dSph satellites demonstrate significant statistical differences from those of the Milky Way, which may be attributable to the different properties of their local environments. Our observations reach the oldest main sequence turn-offs, allowing a time resolution at the oldest ages of ~ 1 Gyr, which is comparable to the best achievable resolution in the MW satellites. We find that the six dSphs present a variety of SFHs that are not strictly correlated with luminosity or present distance from M31. Specifically, we find a significant range in quenching times (lookback times from 9 to 6 Gyr), but with all quenching times more than ~ 6 Gyr ago. In agreement with observations of Milky Way companions of similar mass, there is no evidence of complete quenching of star formation by the cosmic UV background responsible for reionization, but the possibility of a degree of quenching at reionization cannot be ruled out. We do not find significant differences between the SFHs of the three members of the vast, thin plane of satellites and the three off-plane dSphs. The primary difference between the SFHs of the ISLAndS dSphs and Milky Way dSph companions of similar luminosities and host distances is the absence of very late quenching (< 5 Gyr ago) dSphs in the ISLAndS sample. Thus, models that can reproduce satellite populations with and without late quenching satellites will be of extreme interest.

The "too big to fail" problem is revisited by studying the tidal evolution of populations of dwarf satellites with different density profiles. The high resolution cosmological $\rm \Lambda CDM$ "ErisMod" set of simulations is used. These simulations can model both the stellar and dark matter components of the satellites, and their evolution under the action of the tides of a MW-sized host halo at a force resolution better than 10 pc. The stronger tidal mass loss and re-shaping of the mass distribution induced in satellites with $\gamma=0.6$ dark matter density distributions, as those resulting from the effect of feedback in hydrodynamical simulations of dwarf galaxy formation, is sufficient to bring the circular velocity profiles in agreement with the kinematics of MW's dSphs. In contrast, in simulations in which the satellites retain cusps at $z=0$ there are several "massive failures" with circular velocities in excess of the observational constraints. Various sources of deviations in the conventionally adopted relation between the circular velocity at the half light radius and the one dimensional line-of-sight velocity dispersions are found. Such deviations are caused by the response of circular velocity profiles to tidal effects, which also varies depending on the initially assumed inner density profile, and by the complexity of the stellar kinematics, which include residual rotation and anisotropy. In addition tidal effects naturally induce large deviations in the stellar mass-halo mass relation for halo masses below $\rm 10^9 ~ M_{\odot}$, preventing any reliable application of the abundance matching technique to dwarf galaxy satellites.

The analysis of the early star formation history (SFH) of nearby galaxies, obtained from their resolved stellar populations is relevant as a test for cosmological models. However, the early time resolution of observationally derived SFHs is limited by several factors. Thus, direct comparison of observationally derived SFHs with those derived from theoretical models of galaxy formation is potentially biased. Here we investigate and quantify this effect. For this purpose, we analyze the duration of the early star formation activity in a sample of four Local Group dwarf galaxies and test whether they are consistent with being true fossils of the pre-reionization era; i.e., if the quenching of their star formation occurred before cosmic reionization by UV photons was completed. Two classical dSph (Cetus and Tucana) and two dTrans (LGS-3 and Phoenix) isolated galaxies with total stellar masses between $1.3\times 10^6$ to $7.2\times 10^6$ M$_\odot$ have been studied. Accounting for time resolution effects, the SFHs peak as much as 1.25 Gyr earlier than the optimal solutions. Thus, this effect is important for a proper comparison of model and observed SFHs. It is also shown that none of the analyzed galaxies can be considered a true-fossil of the pre-reionization era, although it is possible that the {\it outer regions} of Cetus and Tucana are consistent with quenching by reionization.

Based on data aquired in 13 orbits of HST time, we present a detailed evolutionary history of the M31 dSph satellite Andromeda XVI, including its life-time star formation history, the spatial distribution of its stellar populations, and the properties of its variable stars. And XVI is characterized by prolonged star formation activity from the oldest epochs until star formation was quenched ~6 Gyr ago, and, notably, only half of the mass in stars of And XVI was in place 10 Gyr ago. And XVI appears to be a low mass galaxy for which the early quenching by either reionization or starburst feedback seems highly unlikely, and thus, is most likely due to an environmental effect (e.g., an interaction), possibly connected to a late infall in the densest regions of the Local Group. Studying the star formation history as a function of galactocentric radius, we detect a mild gradient in the star formation history: the star formation activity between 6 and 8 Gyr ago is significantly stronger in the central regions than in the external regions, although the quenching age appears to be the same, within 1 Gyr. We also report the discovery of 9 RR Lyrae stars, 8 of which belong to And XVI. The RR Lyrae stars allow a new estimate of the distance, (m-M)0= 23.72+/-0.09 mag, which is marginally larger than previous estimates based on the tip of the red giant branch.

Phases of gravitational instability are expected in the early phases of disk evolution, when the disk mass is still a substantial fraction of the mass of the star. Disk fragmentation into sub-stellar objects could occur in the cold exterior part of the disk. Direct detection of massive gaseous clumps on their way to collapse into gas giant planets would offer an unprecedented test of the disk instability model. Here we use state-of-the-art 3D radiation-hydro simulations of disks undergoing fragmentation into massive gas giants, post-processed with the RADMC-3D ray-tracing code to produce dust continuum emission maps. These are then fed into the Common Astronomy Software Applications (CASA) ALMA simulator. The synthetic maps show that both overdense spiral arms and actual clumps at different stages of collapse can be detected with the Atacama Large Millimetre/sub-millimetre Array (ALMA) in the full configuration at the distance of the Ophiuchus star forming region (125 pc). The detection of clumps is particularly effective at shorter wavelengths (690 GHz) combining two resolutions with multi-scale clean. Furthermore, we show that a flux-based estimate of the mass of a protoplanetary clump can be from comparable to a factor of 3 higher than the gravitationally bound clump mass. The estimated mass depends on the assumed opacity, and on the gas temperature, which should be set using the input of radiation-hydro simulations. We conclude that ALMA has the capability to detect "smoking gun" systems that are a signpost of the disk instability model for gas giant planet formation.

We study dust capture by vortices and its long-term consequences in global two-fluid inviscid disk simulations using a new polar grid code RoSSBi. We perform the longest integrations so far, several hundred disk orbits, at the highest resolution attainable in global simulations of disks with dust, namely 2048x4096 grid points. This allows to study the dust evolution well beyond vortex dissipation. We vary a wide range of parameters, most notably the dust-to-gas ratio in the initial setup varies in the range $10^{-3}$ to $0.1$. Irrespective of the initial dust-to-gas ratio we find rapid concentration of the dust inside vortices, reaching dust-to-gas ratios of order unity inside the vortex. We present an analytical model that describes very well the dust capture process inside vortices, finding consistent results for all dust-to-gas ratios. A vortex streaming instability develops which causes invariably vortex destruction. After vortex dissipation large-scale dust-rings encompassing a disk annulus form in most cases, which sustain very high dust concentration, approaching ratios of order unity. The rings are long lived lasting as long as the duration of the simulations. They also develop a streaming instability, which manifests itself in eddies at various scales within which the dust forms compact high density clumps. Such clumps would be unstable to gravitational collapse in absence of strong dissipation by viscous forces. When vortices are particularly long lived, rings do not form but dust clumps inside vortices become then long lived features and would likely undergo collapse by gravitational instability. Rings encompass almost an Earth mass of solid material, while even larger masses of dust do accumulate inside vortices in the earlier stage. We argue that rapid planetesimal formation would occur in the dust clumps inside the vortices as well as in the post-vortex ring.

We investigate the properties of halo gas using three cosmological `zoom-in' simulations of realistic Milky Way-galaxy analogs with varying sub-grid physics. In all three cases, the mass of hot ($T > 10^6$ K) halo gas is $\sim 1\%$ of the host's virial mass. The X-ray luminosity of two of the runs is consistent with observations of the Milky Way, while the third simulation is X-ray bright and resembles more closely a very massive, star-forming spiral. Hot halos extend to 140 kpc from the galactic center and are surrounded by a bubble of warm-hot ($T = 10^5 - 10^6$K) gas that extends to the virial radius. Simulated halos agree well outside 20-30 kpc with the $\beta$-model of Miller & Bregman (2013) based on OVII absorption and OVIII emission measurements. Warm--hot and hot gas contribute up to $80\%$ of the total gas reservoir, and contain nearly the same amount of baryons as the stellar component. The mass of warm-hot and hot components falls into the range estimated for $L^*$ galaxies. With key observational constraints on the density of the Milky Way corona being satisfied, we show that concealing of the ubiquitous warm-hot baryons, along with the ejection of just $20-30 \%$ of the diffuse gas out of the potential wells by supernova-driven outflows, can solve the "missing baryon problem". The recovered baryon fraction within 3 virial radii is $90\%$ of the universal value. With a characteristic density of $\sim 10^{-4}$ cm$^{-3}$ at $50-80$ kpc, diffuse coronae meet the requirement for fast and complete ram--pressure stripping of the gas reservoirs in dwarf galaxy satellites.

The rapid assembly of the massive black holes that power the luminous quasars observed at $z \sim 6-7$ remains a puzzle. Various direct collapse models have been proposed to head-start black hole growth from initial seeds with masses $\sim 10^5\,\rm M_\odot$, which can then reach a billion solar mass while accreting at the Eddington limit. Here we propose an alternative scenario based on radiatively inefficient super-critical accretion of stellar-mass holes embedded in the gaseous circum-nuclear discs (CNDs) expected to exist in the cores of high redshift galaxies. Our sub-pc resolution hydrodynamical simulations show that stellar-mass holes orbiting within the central 100 pc of the CND bind to very high density gas clumps that arise from the fragmentation of the surrounding gas. Owing to the large reservoir of dense cold gas available, a stellar-mass black hole allowed to grow at super-Eddington rates according to the "slim disc" solution can increase its mass by 3 orders of magnitudes within a few million years. These findings are supported by simulations run with two different hydro codes, RAMSES based on the Adaptive Mesh Refinement technique and GIZMO based on a new Lagrangian Godunov-type method, and with similar, but not identical, sub-grid recipes for star formation, supernova feedback, black hole accretion and feedback. The low radiative efficiency of super-critical accretion flows are instrumental to the rapid mass growth of our black holes, as they imply modest radiative heating of the surrounding nuclear environment.

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)

We study the effect of sub-grid physics, galaxy mass, structural parameters and resolution on the fragmentation of gas-rich galaxy discs into massive star forming clumps. The initial conditions are set up with the aid of the ARGO cosmological hydrodynamical simulation. Blast-wave feedback does not suppress fragmentation, but reduces both the number of clumps and the duration of the unstable phase. Once formed, bound clumps cannot be destroyed by our feedback model. Widespread fragmentation is promoted by high gas fractions and low halo concentrations. Yet giant clumps $M > 10^8 M_{\odot}$ lasting several hundred Myr are rare and mainly produced by clump-clump mergers. They occur in massive discs with maximum rotational velocities $V_{max} > 250$ km/s at $z \sim 2$, at the high mass end of the observed galaxy population at those redshifts. The typical gaseous and stellar masses of clumps in all runs are in the range $\sim 10^7-10^8 M_{\odot}$ for galaxies with disc mass in the range $10^{10}-8\times 10^{10} M_{\odot}$. Clumps sizes are usually in the range $100-400$ pc, in agreement with recent clump observations in lensed high-z galaxies. \\ We argue that many of the giant clumps identified in observations are not due to in-situ fragmetation, or are the result of blending of smaller structures owing to insufficient resolution. Using an analytical model describing local collapse inside spiral arms, we can predict the characteristic gaseous masses of clumps at the onset of fragmentation ($\sim 3-5 \times 10^7 M_{\odot}$) quite accurately, while the conventional Toomre mass overestimates them. Due to their moderate masses, clumps which migrate to the centre have marginal effect on bulge growth.

We discuss how knowledge of the whole evolutionary history of dwarf galaxies, including details on the early star formation events, can provide insight on the origin of the different dwarf galaxy types. We suggest that these types may be imprinted by the early conditions of formation rather than being only the result of a recent morphological transformation driven by environmental effects. We present precise star formation histories of a sample of Local Group dwarf galaxies, derived from colour-magnitude diagrams reaching the oldest main-sequence turnoffs. We argue that these galaxies can be assigned to two basic types: fast dwarfs that started their evolution with a dominant and short star formation event, and slow dwarfs that formed a small fraction of their stars early and have continued forming stars until the present time (or almost). These two different evolutionary paths do not map directly onto the present-day morphology (dwarf spheroidal vs dwarf irregular). Slow and fast dwarfs also differ in their inferred past location relative to the Milky Way and/or M31, which hints that slow dwarfs were generally assembled in lower density environments than fast dwarfs. We propose that the distinction between a fast and slow dwarf galaxy reflects primarily the characteristic density of the environment where they form. At a later stage, interaction with a large host galaxy may play a role in the final gas removal and ultimate termination of star formation.

We present the first simulations of tidal stirring of dwarf galaxies in the Local Group carried out in a cosmological context. We use the ErisDARK simulation of a MW-sized galaxy to identify some of the most massive subhalos ($M_{vir} > 10^8 M_{\odot}$) that fall into the main host before $z=2$. Subhalos are replaced before infall with high-resolution models of dwarf galaxies comprising a faint stellar disk embedded in a dark matter halo. The set of models contains cuspy halos as well as halos with "cored" profiles (with asymptotic inner slope $\gamma = 0.6$). The simulations are then run to $z=0$ with as many as 54 million particles and resolution as small as $\sim 4$ pc using the N-Body code ChaNGa. The stellar components of all satellites are significantly affected by tidal stirring, losing stellar mass and undergoing a morphological transformation towards a pressure supported spheroidal system. However, while some remnants with cuspy halos maintain significant rotational flattening and disk-like features, all the shallow halo models achieve $v/\sigma < 0.5$ and round shapes typical of dSph satellites of the MW and M31. Mass loss is also enhanced in the latter, and remnants can reach luminosities and velocity dispersions as low as those of Ultra Faint Dwarfs (UFDs). We argue that cuspy progenitors must be the exception rather than the rule among satellites of the MW since all the MW and M31 satellites in the luminosity range of our remnants are dSphs, a result matched only in the simulation with "cored" models.

We investigate the production sites and the enrichment history of $r$-process elements in the Galaxy, as traced by the [Eu/Fe] ratio, using the high resolution, cosmological zoom-in simulation `Eris'. At $z=0$, Eris represents a close analog to the Milky Way, making it the ideal laboratory to understand the chemical evolution of our Galaxy. Eris formally traces the production of oxygen and iron due to Type-Ia and Type-II supernovae. We include in post-processing the production of $r$-process elements from compact binary mergers. Unlike previous studies, we find that the nucleosynthetic products from compact binary mergers can be incorporated into stars of very low metallicity and at early times, even with a minimum delay time of 100 Myr. This conclusion is relatively insensitive to modest variations in the merger rate, minimum delay time, and the delay time distribution. By implementing a first-order prescription for metal-mixing, we can further improve the agreement between our model and the data for the chemical evolution of both [$\alpha$/Fe] and [Eu/Fe].We argue that compact binary mergers could be the dominant source of $r$-process nucleosynthesis in the Galaxy.

We perform two-dimensional hydrodynamical simulations to quantitatively explore the torque balance criterion for gap-opening (as formulated by Crida et al. 2006) in a variety of disks when considering a migrating planet. We find that even when the criterion is satisfied, there are instances when planets still do not open gaps. We stress that gap-opening is not only dependent on whether a planet has the ability to open a gap, but whether it can do so quickly enough. This can be expressed as an additional condition on the gap-opening timescale versus the crossing time, i.e. the time it takes the planet to cross the region which it is carving out. While this point has been briefly made in the previous literature, our results quantify it for a range of protoplanetary disk properties and planetary masses, demonstrating how crucial it is for gap-opening. This additional condition has important implications for the survival of planets formed by core accretion in low mass disks as well as giant planets or brown dwarfs formed by gravitational instability in massive disks. It is particularly important for planets with intermediate masses susceptible to Type III-like migration. For some observed transition disks or disks with gaps, we expect that estimates on the potential planet masses based on the torque balance gap-opening criterion alone may not be sufficient. With consideration of this additional timescale criterion theoretical studies may find a reduced planet survivability or that planets may migrate further inwards before opening a gap.

We study supermassive black holes (BHs) in merging galaxies, using a suite of hydrodynamical simulations with very high spatial (~10 pc) and temporal (~1 Myr) resolution, where we vary the initial mass ratio, the orbital configuration, and the gas fraction. (i) We address the question of when and why, during a merger, increased BH accretion occurs, quantifying gas inflows and BH accretion rates. (ii) We also quantify the relative effectiveness in inducing AGN activity of merger-related versus secular-related causes, by studying different stages of the encounter: the stochastic (or early) stage, the (proper) merger stage, and the remnant (or late) stage. (iii) We assess which galaxy mergers preferentially enhance BH accretion, finding that the initial mass ratio is the most important factor. (iv) We study the evolution of the BH masses, finding that the BH mass contrast tends to decrease in minor mergers and to increase in major mergers. This effect hints at the existence of a preferential range of mass ratios for BHs in the final pairing stages. (v) In both merging and dynamically quiescent galaxies, the gas accreted by the BH is not necessarily the gas with $low$ angular momentum, but the gas that $loses$ angular momentum.

We use high resolution Hydro$+$N-Body cosmological simulations to compare the assembly and evolution of a small field dwarf (stellar mass ~ 10$^{6-7}$ M$\odot$, total mass 10$^{10}$ M$\odot$ in $\Lambda$ dominated CDM and 2keV WDM cosmologies. We find that star formation (SF) in the WDM model is reduced and delayed by 1-2 Gyr relative to the CDM model, independently of the details of SF and feedback. Independent of the DM model, but proportionally to the SF efficiency, gas outflows lower the central mass density through `dynamical heating', such that all realizations have circular velocities $<$ 20kms at 500$~$pc, in agreement with local kinematic constraints. As a result of dynamical heating, older stars are less centrally concentrated than younger stars, similar to stellar population gradients observed in nearby dwarf galaxies. Introducing an important diagnostic of SF and feedback models, we translate our simulations into artificial color-magnitude diagrams and star formation histories in order to directly compare to available observations. The simulated galaxies formed most of their stars in many $\sim$10 Myr long bursts. The CDM galaxy has a global SFH, HI abundance and Fe/H and alpha-elements distribution well matched to current observations of dwarf galaxies. These results highlight the importance of directly including `baryon physics' in simulations when 1) comparing predictions of galaxy formation models with the kinematics and number density of local dwarf galaxies and 2) differentiating between CDM and non-standard models with different DM or power spectra.

We present novel 3D multi-scale SPH simulations of gas-rich galaxy mergers between the most massive galaxies at $z \sim 8 - 10$, designed to scrutinize the direct collapse formation scenario for massive black hole seeds proposed in \citet{mayer+10}. The simulations achieve a resolution of 0.1 pc, and include both metallicity-dependent optically-thin cooling and a model for thermal balance at high optical depth. We consider different formulations of the SPH hydrodynamical equations, including thermal and metal diffusion. When the two merging galaxy cores collide, gas infall produces a compact, optically thick nuclear disk with densities exceeding $10^{-10}$ g cm$^3$. The disk rapidly accretes higher angular momentum gas from its surroundings reaching $\sim 5$ pc and a mass of $\gtrsim 10^9$ $M_{\odot}$ in only a few $10^4$ yr. Outside $\gtrsim 2$ pc it fragments into massive clumps. Instead, supersonic turbulence prevents fragmentation in the inner parsec region, which remains warm ($\sim 3000-6000$ K) and develops strong non-axisymmetric modes that cause prominent radial gas inflows ($> 10^4$ $M_{\odot}$ yr$^{-1}$), forming an ultra-dense massive disky core. Angular momentum transport by non-axisymmetric modes should continue below our spatial resolution limit, quickly turning the disky core into a supermassive protostar which can collapse directly into a massive black hole of mass $10^8-10^9$ $M_{\odot}$ via the relativistic radial instability. Such a "cold direct collapse"' explains naturally the early emergence of high-z QSOs. Its telltale signature would be a burst of gravitational waves in the frequency range $10^{-4} - 10^{-1}$ Hz, possibly detectable by the planned eLISA interferometer.

The Hubble sequence is a common classification scheme for the structure of galaxies. Despite the tremendous usefulness of this diagnostic, we still do not fully understand when, where, and how this morphological ordering was put in place. Here, we investigate the morphological evolution of a sample of 22 high redshift ($z\geq3$) galaxies extracted from the Argo simulation. Argo is a cosmological zoom-in simulation of a group-sized halo and its environment. It adopts the same high resolution ($\sim10^4$ M$_\odot$, $\sim100$ pc) and sub-grid physical model that was used in the Eris simulation but probes a sub-volume almost ten times bigger with as many as 45 million gas and star particles in the zoom-in region. Argo follows the early assembly of galaxies with a broad range of stellar masses ($\log M_{\star}/{\rm M}_{\odot}\sim8-11$ at $z\simeq3$), while resolving properly their structural properties. We recover a diversity of morphologies, including late-type/irregular disc galaxies with flat rotation curves, spheroid dominated early-type discs, and a massive elliptical galaxy, already established at $z\sim3$. We identify major mergers as the main trigger for the formation of bulges and the steepening of the circular velocity curves. Minor mergers and non-axisymmetric perturbations (stellar bars) drive the bulge growth in some cases. The specific angular momenta of the simulated disc components fairly match the values inferred from nearby galaxies of similar $M_{\star}$ once the expected redshift evolution of disc sizes is accounted for. We conclude that morphological transformations of high redshift galaxies of intermediate mass are likely triggered by processes similar to those at low redshift and result in an early build-up of the Hubble sequence.

Observations show a prevalence of high redshift galaxies with large stellar masses and predominantly passive stellar populations. A variety of processes have been suggested that could reduce the star formation in such galaxies to observed levels, including quasar mode feedback, virial shock heating, or galactic winds driven by stellar feedback. However, the main quenching mechanisms have yet to be identified. Here we study the origin of star formation quenching using Argo, a cosmological, hydrodynamical zoom-in simulation that follows the evolution of a massive galaxy at $z\geq{}2$. This simulation adopts the same sub-grid recipes of the Eris simulations, which have been shown to form realistic disk galaxies, and, in one version, adopts also a mass and spatial resolution identical to Eris. The resulting galaxy has properties consistent with those of observed, massive (M_* ~ 1e11 M_sun) galaxies at z~2 and with abundance matching predictions. Our models do not include AGN feedback indicating that supermassive black holes likely play a subordinate role in determining masses and sizes of massive galaxies at high z. The specific star formation rate (sSFR) of the simulated galaxy matches the observed M_* - sSFR relation at early times. This period of smooth stellar mass growth comes to a sudden halt at z=3.5 when the sSFR drops by almost an order of magnitude within a few hundred Myr. The suppression is initiated by a leveling off and a subsequent reduction of the cool gas accretion rate onto the galaxy, and not by feedback processes. This "cosmological starvation" occurs as the parent dark matter halo switches from a fast collapsing mode to a slow accretion mode. Additional mechanisms, such as perhaps radio mode feedback from an AGN, are needed to quench any residual star formation of the galaxy and to maintain a low sSFR until the present time.

We model the dynamics of dwarf early-type galaxies in the Virgo cluster when subject to a variety of environmental processes. We focus on how these processes imprint trends in dynamical state (rotational vs. pressure support as measured by the $\lambda^*_{\rm Re/2}$ statistic) with projected distance from the cluster center, and compare these results to observational estimates. We find a large scatter in the gradient of $\lambda^*_{\rm Re/2}$ with projected radius. A statistical analysis shows that models with no environmental effects produce gradients as steep as those observed in none of the 100 cluster realizations we consider, while in a model incorporating tidal stirring by the cluster potential 34% of realizations produce gradients as steep as that observed. Our results suggest that tidal stirring may be the cause of the observed radial dependence of dwarf early-type dynamics in galaxy clusters.

We analyze the formation and evolution of the stellar components in "Eris", a 120 pc-resolution cosmological hydrodynamic simulation of a late-type spiral galaxy. The simulation includes the effects of a uniform UV background, a delayed-radiative-cooling scheme for supernova feedback, and a star formation recipe based on a high gas density threshold. It allows a detailed study of the relative contributions of "in-situ" (within the main host) and "ex-situ" (within satellite galaxies) star formation to each major Galactic component in a close Milky Way analog. We investigate these two star-formation channels as a function of galactocentric distance, along different lines of sight above and along the disk plane, and as a function of cosmic time. We find that: 1) approximately 70 percent of today's stars formed in-situ; 2) more than two thirds of the ex-situ stars formed within satellites after infall; 3) the majority of ex-situ stars are found today in the disk and in the bulge; 4) the stellar halo is dominated by ex-situ stars, whereas in-situ stars dominate the mass profile at distances < 5 kpc from the center at high latitudes; and 5) approximately 25% of the inner, r < 20 kpc, halo is composed of in-situ stars that have been displaced from their original birth sites during Eris' early assembly history.

We present results from a fully cosmological, very high-resolution, LCDM "zoom-in" simulation of a group of seven field dwarf galaxies with present-day virial masses in the range M_vir=4.4e8-3.6e10 Msun. The simulation includes a blastwave scheme for supernova feedback, a star formation recipe based on a high gas density threshold, metal-dependent radiative cooling, a scheme for the turbulent diffusion of metals and thermal energy, and a uniform UV background that modifies the ionization and excitation state of the gas. The properties of the simulated dwarfs are strongly modulated by the depth of the gravitational potential well. All three halos with M_vir < 1e9 Msun are devoid of stars, as they never reach the density threshold for star formation of 100 atoms/cc. The other four, M_vir > 1e9 Msun dwarfs have blue colors, low star formation efficiencies, high cold gas to stellar mass ratios, and low stellar metallicities. Their bursty star formation histories are characterized by peak specific star formation rates in excess of 50-100 1/Gyr, far outside the realm of normal, more massive galaxies, and in agreement with observations of extreme emission-line starbursting dwarfs by the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey. Metal-enriched galactic outflows produce sub-solar effective yields and pollute with heavy elements a Mpc-size region of the intergalactic medium, but are not sufficient to completely quench star formation activity and are not ubiquitous in our dwarfs. Within the limited size of the sample, our simulations appear to simultaneously reproduce the observed stellar mass and cold gas content, resolved star formation histories, stellar kinematics, and metallicities of field dwarfs in the Local Volume.

We simulate an equal-mass merger of two Milky Way-size galaxy discs with moderate gas fractions at parsec-scale resolution including a new model for radiative cooling and heating in a multi-phase medium, as well as star formation and feedback from supernovae. The two discs initially have a $2.6\times10^6\mathrm{~M_{\odot}}$ supermassive black hole (SMBH) embedded in their centers. As the merger completes and the two galactic cores merge, the SMBHs form a a pair with a separation of a few hundred pc that gradually decays. Due to the stochastic nature of the system immediately following the merger, the orbital plane of the binary is significantly perturbed. Furthermore, owing to the strong starburst the gas from the central region is completely evacuated, requiring $\sim10$~Myr for a nuclear disc to rebuild. Most importantly, the clumpy nature of the interstellar medium has a major impact on the the dynamical evolution of the SMBH pair, which undergo gravitational encounters with massive gas clouds and stochastic torquing by both clouds and spiral modes in the disk. These effects combine to greatly delay the decay of the two SMBHs to separations of a few parsecs by nearly two orders of magnitude, $\sim 10^8$ yr, compared to previous work. In mergers of more gas-rich, clumpier galaxies at high redshift stochastic torques will be even more pronounced and potentially lead to stronger modulation of the orbital decay. This suggests that SMBH pairs at separations of several tens of parsecs should be relatively common at any redshift.

We present the first comparison between the lifetime star formation histories (SFHs) of M31 and Milky Way (MW) satellites. Using the Advanced Camera for Surveys aboard the Hubble Space Telescope, we obtained deep optical imaging of Andromeda II (M$_{V} = -$12.0; log(M$_{\star}$/M$_{\odot}$) $\sim$ 6.7) and Andromeda XVI (M$_{V} = -$7.5; log(M$_{\star}$/M$_{\odot}$) $\sim$ 4.9) yielding color-magnitude diagrams that extend at least 1 magnitude below the oldest main sequence turnoff, and are similar in quality to those available for the MW companions. And II and And XVI show strikingly similar SFHs: both formed 50-70% of their total stellar mass between 12.5 and 5 Gyr ago (z$\sim$5-0.5) and both were abruptly quenched $\sim$ 5 Gyr ago (z$\sim$0.5). The predominance of intermediate age populations in And XVI makes it qualitatively different from faint companions of the MW and clearly not a pre-reionization fossil. Neither And II nor And XVI appears to have a clear analog among MW companions, and the degree of similarity in the SFHs of And II and And XVI is not seen among comparably faint-luminous pairs of MW satellites. These findings provide hints that satellite galaxy evolution may vary substantially among hosts of similar stellar mass. Although comparably deep observations of more M31 satellites are needed to further explore this hypothesis, our results underline the need for caution when interpreting satellite galaxies of an individual system in a broader cosmological context.

We present an analysis of the star formation history (SFH) of a field near the half light radius in the Local Group dwarf irregular galaxy IC 1613 based on deep Hubble Space Telescope Advanced Camera for Surveys imaging. Our observations reach the oldest main sequence turn-off, allowing a time resolution at the oldest ages of ~1 Gyr. Our analysis shows that the SFH of the observed field in IC 1613 is consistent with being constant over the entire lifetime of the galaxy. These observations rule out an early dominant episode of star formation in IC 1613. We compare the SFH of IC 1613 with expectations from cosmological models. Since most of the mass is in place at early times for low mass halos, a naive expectation is that most of the star formation should have taken place at early times. Models in which star formation follows mass accretion result in too many stars formed early and gas mass fractions which are too low today (the "over-cooling problem"). The depth of the present photometry of IC 1613 shows that, at a resolution of ~1 Gyr, the star formation rate is consistent with being constant, at even the earliest times, which is difficult to achieve in models where star formation follows mass assembly.

We present near-infrared integral field spectroscopy of the central kiloparsec of 17 nearby luminous and ultra-luminous infrared galaxies undergoing major mergers. These observations were taken with OSIRIS assisted by the Keck I and II Adaptive Optics systems, providing spatial resolutions of a few tens of parsecs. The resulting kinematic maps reveal gas disks in at least 16 out of 19 nuclei and stellar disks in 11 out of 11 nuclei observed in these galaxy merger systems. In our late-stages mergers, these disks are young (stellar ages $<30$ Myr) and likely formed as gas disks which became unstable to star formation during the merger. On average, these disks have effective radii of a few hundred parsecs, masses between $10^{8}$ and $10^{10} M_{Sun}$, and $v/\sigma$ between 1 and 5. These disks are similar to those created in high-resolution hydrodynamical simulations of gas-rich galaxy mergers, and favor short coalescence times for binary black holes. The few galaxies in our sample in earlier stages of mergers have disks which are larger ($r_{eff}\sim200-1800$ pc) and likely are remnants of the galactic disks that have not yet been completely disrupted by the merger.

We study the dynamical evolution of supermassive black holes (BHs) in merging galaxies on scales of hundreds of kpc to 10 pc, to identify the physical processes that aid or hinder the orbital decay of BHs. We present hydrodynamical simulations of galaxy mergers with a resolution of $\leq$20 pc, chosen to accurately track the motion of the nuclei and provide a realistic environment for the evolution of the BHs. We find that, during the late stages of the merger, tidal shocks inject energy in the nuclei, causing one or both nuclei to be disrupted and leaving their BH `naked', without any bound gas or stars. In many cases, the nucleus that is ultimately disrupted is that of the larger galaxy (`nuclear coup'), as star formation grows a denser nuclear cusp in the smaller galaxy. We supplement our simulations with an analytical estimate of the orbital-decay time required for the BHs to form a binary at unresolved scales, due to dynamical friction. We find that, when a nuclear coup occurs, the time-scale is much shorter than when the secondary's nucleus is disrupted, as the infalling BH is more massive, and it also finds itself in a denser stellar environment.

We introduce the AGORA project, a comprehensive numerical study of well-resolved galaxies within the LCDM cosmology. Cosmological hydrodynamic simulations with force resolutions of ~100 proper pc or better will be run with a variety of code platforms to follow the hierarchical growth, star formation history, morphological transformation, and the cycle of baryons in and out of 8 galaxies with halo masses M_vir ~= 1e10, 1e11, 1e12, and 1e13 Msun at z=0 and two different ("violent" and "quiescent") assembly histories. The numerical techniques and implementations used in this project include the smoothed particle hydrodynamics codes GADGET and GASOLINE, and the adaptive mesh refinement codes ART, ENZO, and RAMSES. The codes will share common initial conditions and common astrophysics packages including UV background, metal-dependent radiative cooling, metal and energy yields of supernovae, and stellar initial mass function. These are described in detail in the present paper. Subgrid star formation and feedback prescriptions will be tuned to provide a realistic interstellar and circumgalactic medium using a non-cosmological disk galaxy simulation. Cosmological runs will be systematically compared with each other using a common analysis toolkit, and validated against observations to verify that the solutions are robust - i.e., that the astrophysical assumptions are responsible for any success, rather than artifacts of particular implementations. The goals of the AGORA project are, broadly speaking, to raise the realism and predictive power of galaxy simulations and the understanding of the feedback processes that regulate galaxy "metabolism." The proof-of-concept dark matter-only test of the formation of a galactic halo with a z=0 mass of M_vir ~= 1.7e11 Msun by 9 different versions of the participating codes is also presented to validate the infrastructure of the project.

The large number of detected giant exoplanets offers the opportunity to improve our understanding of the formation mechanism, evolution, and interior structure of gas giant planets. The two main models for giant planet formation are core accretion and disk instability. There are substantial differences between these formation models, including formation timescale, favorable formation location, ideal disk properties for planetary formation, early evolution, planetary composition, etc. First, we summarize the two models including their substantial differences, advantages, and disadvantages, and suggest how theoretical models should be connected to available (and future) data. We next summarize current knowledge of the internal structures of solar- and extrasolar- giant planets. Finally, we suggest the next steps to be taken in giant planet exploration.

We study the dynamics of massive black hole pairs in clumpy gaseous circumnuclear disks. We track the orbital decay of the light, secondary black hole $M_{\bullet2}$ orbiting around the more massive primary at the center of the disk, using $N$-body/smoothed particle hydrodynamic simulations. We find that the gravitational interaction of $M_{\bullet2}$ with massive clumps $M_{\rm cl}$ erratically perturbs the otherwise smooth orbital decay. In close encounters with massive clumps, gravitational slingshots can kick the secondary black hole out of the disk plane. The black hole moving on an inclined orbit then experiences the weaker dynamical friction of the stellar background, resulting in a longer orbital decay timescale. Interactions between clumps can also favor orbital decay when the black hole is captured by a massive clump which is segregating toward the center of the disk. The stochastic behavior of the black hole orbit emerges mainly when the ratio $M_{\bullet2}/M_{\rm cl}$ falls below unity, with decay timescales ranging from $\sim1$ to $\sim50$ Myr. This suggests that describing the cold clumpy phase of the inter-stellar medium in self-consistent simulations of galaxy mergers, albeit so far neglected, is important to predict the black hole dynamics in galaxy merger remnants.

Although it is fairly established that Gravitational Instability (GI) should occur in the early phases of the evolution of a protoplanetary disk, the fate of the clumps resulting from disk fragmentation and their role in planet formation is still unclear. In the present study we investigate semi-analytically their evolution following the contraction of a synthetic population of clumps with varied initial structure and orbits coupled with the surrounding disk and the central star. Our model is based on recently published state-of-the-art 3D collapse simulations of clumps with varied thermodynamics. Various evolutionary mechanisms are taken into account, and their effect is explored both individually and in combination with others: migration and tidal disruption, mass accretion, gap opening and disk viscosity. It is found that, in general, at least 50% of the initial clumps survive tides, leaving behind potential gas giant progenitors after ~10^5 yr of evolution in the disk. The rest might be either disrupted or produce super-Earths and other low mass planets provided that a solid core can be assembled on a sufficiently short timescale, a possibility that we do not address in this paper. Extrapolating to million year timescales, all our surviving protoplanets would lead to close-in gas giants. This outcome might in part reflect the limitations of the migration model adopted, and is reminiscent of the analogous result found in core-accretion models in absence of fine-tuning of the migration rate. Yet it suggests that a significant fraction of the clumps formed by gravitational instability could be the precursors of Hot Jupiters.

Based on Hubble Space Telescope observations from the Local Cosmology from Isolated Dwarfs project, we present the star formation histories, as a function of galactocentric radius, of four isolated Local Group dwarf galaxies: two dSph galaxies, Cetus and Tucana, and two transition galaxies (dTrs), LGS-3 and Phoenix. The oldest stellar populations of the dSphs and dTrs are, within the uncertainties, coeval ($\sim 13 Gyr$) at all galactocentric radii. We find that there are no significative differences between the four galaxies in the fundamental properties (such as the normalized star formation rate or age-metallicity relation) of their outer regions (radii greater than four exponential scale lengths); at large radii, these galaxies consist exclusively of old ($\geq 10.5 Gyr$) metal-poor stars. The duration of star formation in the inner regions vary from galaxy to galaxy, and the extended central star formation in the dTrs produces the dichotomy between dSph and dTr galaxy types. The dTr galaxies show prominent radial stellar population gradients: the centers of these galaxies host young ($\leq 1 Gyr$) populations while the age of the last formation event increases smoothly with increasing radius. This contrasts with the two dSph galaxies. Tucana shows a similar, but milder, gradient, but no gradient in age is detected Cetus. For the three galaxies with significant stellar population gradients, the exponential scale length decreases with time. These results are in agreement with outside-in scenarios of dwarf galaxy evolution, in which a quenching of the star formation toward the center occurs as the galaxy runs out of gas in the outskirts.

In some scenarios for the formation of the Milky Way bulge the stellar population at the edges of the boxy bulge may be younger than those on the minor axis, or close to the Galactic center. So far the only bulge region where deep color-magnitude diagrams have been obtained is indeed along the minor axis. To overcome this limitation, we aim at age-dating the bulge stellar populations far away from the bulge minor axis. Color-magnitude diagrams and luminosity functions have been obtained from deep near-IR VLT/HAWK-I images taken at the two Southern corners of the boxy bulge, i.e., near the opposite edges of the Galactic bar. The foreground disk contamination has been statistically removed using a pure disk field observed with the same instrument and located approximately at similar Galactic latitudes of the two bulge fields, and 30deg in longitude away from the Galactic center. For each bulge field, mean reddening and distance are determined using the position of red clump stars, and the metallicity distribution is derived photometrically using the color distribution of stars in the upper red giant branch. The resulting metallicity distribution function of both fields peaks around [Fe/H] -0.1 dex, with the bulk of the stellar population having a metallicity within the range: -1 dex < [Fe/H] <+0.4 dex, quite similar to that of other inner bulge fields. Like for the inner fields previously explored, the color-magnitude diagrams of the two bar fields are consistent with their stellar population being older than ~10 Gyr, with no obvious evidence for the presence of a younger population.

We revisit the phases of the pairing and sinking of BHs in galaxy mergers and circunmunclear disks in light of the results of recent simulations with massive BHs embedded in predominantly gaseous backgrounds. After a general overview we discuss the importance of a fast orbital decay regime dominated by global disk torques rather than by the local dynamical friction wake. This regime can dominate at BH binary separations of a few tens of parsecs and below, following a phase of orbital circularization dominated by local dynamical friction. It is similar to Type-I migration in planetary evolution. It can bring the black holes to separations small enough for gravitational waves to take over on a timescale ranging from less than $\sim 10^7$ yr to up to $10^8$ yr, depending on whether the interstellar medium is smooth or clumpy. Eventual gap opening at sub-pc scale separations slows down but does not interrupt the orbital decay.Subsequently, we discuss a new intriguing connection between the conditions required for rapid orbital decay of massive BH binaries and those required for prominent gas inflows in gas-rich galaxies. We derive a condition for the maximum inflow rate that a circumnuclear disk can host while still maintaining a sufficiently high gas density at large radii to sustain the decay of a BH binary. We find that gas inflows rates exceeding 10 $M_{\odot}$/yr, postulated to form massive BH seeds in some direct collapse models, would stifle the sinking of massive BH binaries in gas-dominated galactic nuclei. Vice-versa, lower inflow rates, below a solar mass per year, as required to feed typical AGNs, are compatible with a fast orbital decay of BH binaries across a wide range of masses.

We discuss a novel approach to "weighing" the Milky Way dark matter halo, one that combines the latest samples of halo stars selected from the Sloan Digital Sky Survey (SDSS) with state-of-the-art numerical simulations of Milky Way analogs. The fully cosmological runs employed in the present study include "Eris", one of the highest-resolution hydrodynamical simulations of the formation of a M_vir=8e11 M_sun late-type spiral, and the dark-matter only M_vir=1.7e12 M_sun "Via Lactea II" simulation. Eris provides an excellent laboratory for creating mock SDSS samples of tracer halo stars, and we successfully compare their density, velocity anisotropy, and radial velocity dispersion profiles with the observational data. Most mock SDSS realizations show the same "cold veil" recently observed in the distant stellar halo of the Milky Way, with tracers as cold as sigma_los ~ 50 km/s between 100-150 kpc. Controlled experiments based on the integration of the spherical Jeans equation as well as a particle tagging technique applied to Via Lactea II show that a "heavy" M_vir 2e12 M_sun realistic host produces a poor fit to the kinematic SDSS data. We argue that these results offer added evidence for a "light," centrally-concentrated Milky Way halo.

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.

The last century has seen enormous progress in our understanding of the Universe. We know the life cycles of stars, the structure of galaxies, the remnants of the big bang, and have a general understanding of how the Universe evolved. We have come remarkably far using electromagnetic radiation as our tool for observing the Universe. However, gravity is the engine behind many of the processes in the Universe, and much of its action is dark. Opening a gravitational window on the Universe will let us go further than any alternative. Gravity has its own messenger: Gravitational waves, ripples in the fabric of spacetime. They travel essentially undisturbed and let us peer deep into the formation of the first seed black holes, exploring redshifts as large as z ~ 20, prior to the epoch of cosmic re-ionisation. Exquisite and unprecedented measurements of black hole masses and spins will make it possible to trace the history of black holes across all stages of galaxy evolution, and at the same time constrain any deviation from the Kerr metric of General Relativity. eLISA will be the first ever mission to study the entire Universe with gravitational waves. eLISA is an all-sky monitor and will offer a wide view of a dynamic cosmos using gravitational waves as new and unique messengers to unveil The Gravitational Universe. It provides the closest ever view of the early processes at TeV energies, has guaranteed sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales around singularities and black holes, all the way to cosmological dimensions.

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.

We present new results on the kinematics, thermal and ionization state, and spatial distribution of metal-enriched gas in the circumgalactic medium (CGM) of massive galaxies at redshift 3, using the "Eris" suite of cosmological "zoom-in" simulations. The reference run adopts a blastwave scheme for supernova feedback that produces galactic outflows, a star formation recipe based on a high gas density threshold, metal-dependent radiative cooling, and a model for the diffusion of metals and thermal energy. Synthetic spectra through the multiphase CGM produce interstellar absorption line strengths of Lya, CII, CIV, SiII, and SiIV as a function of galactocentric impact parameter (scaled to the virial radius) that are in broad agreement with those observed at high-redshift by Steidel et al. (2010). Only about one third of all the gas within R_vir is outflowing. The fraction of sightlines within one virial radius that intercept optically thick material is 27%, in agreement with recent observations by Rudie et al. (2012). Such optically thick absorption is shown to trace inflowing "cold" streams that penetrate deep inside the virial radius. The streams, enriched to metallicities above 0.01 solar, give origin to strong (log N > 13) CII absorption with a covering factor of 22% within R_vir and 10% within 2 R_vir. Galactic outflows do not cause any substantial suppression of the cold accretion mode. The central galaxy is surrounded by a large OVI halo, with a typical column density log N>14 and a near unity covering factor maintained all the way out to 150 kpc. This matches the trends recently observed in star-forming galaxies at low redshift by Tumlinson et al. (2011). Our zoom-in simulations of this single system appear to reproduce quantitatively the complex baryonic processes that determine the exchange of matter, energy, and metals between galaxies and their surroundings. (Abridged)

We show that the position of the central dark matter density peak may be expected to differ from the dynamical center of the Galaxy by several hundred parsec. In Eris, a high resolution cosmological hydrodynamics simulation of a realistic Milky-Way-analog disk galaxy, this offset is 300 - 400 pc (~3 gravitational softening lengths) after z=1. In its dissipationless dark-matter-only twin simulation ErisDark, as well as in the Via Lactea II and GHalo simulations, the offset remains below one softening length for most of its evolution. The growth of the DM offset coincides with a flattening of the central DM density profile in Eris inwards of ~1 kpc, and the direction from the dynamical center to the point of maximum DM density is correlated with the orientation of the stellar bar, suggesting a bar-halo interaction as a possible explanation. A dark matter density offset of several hundred parsec greatly affects expectations of the dark matter annihilation signals from the Galactic Center. It may also support a dark matter annihilation interpretation of recent reports by Weniger (2012) and Su & Finkbeiner (2012) of highly significant 130 GeV gamma-ray line emission from a region 1.5 degrees (~200 parsec projected) away from Sgr A* in the Galactic plane.

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

We study the formation of a supermassive black hole (SMBH) binary and the shrinking of the separation of the two holes to sub-pc scales starting from a realistic major merger between two gas-rich spiral galaxies with mass comparable to our Milky Way. The simulations, carried out with the Adaptive Mesh Refinement (AMR) code RAMSES, are capable of resolving separations as small as 0.1 pc. The collision of the two galaxies produces a gravo-turbulent rotating nuclear disk with mass (10^9 Msun) and size (60 pc) in excellent agreement with previous SPH simulations with particle splitting that used a similar setup (Mayer et al. 2007) but were limited to separations of a few parsecs. The AMR results confirm that the two black holes sink rapidly as a result of dynamical friction onto the gaseous background, reaching a separation of 1 pc in less than 10^7 yr. We show that the dynamical friction wake is well resolved by our model and we find good agreement with analytical predictions of the drag force as a function of the Mach number. Below 1 pc, black hole pairing slows down significantly, as the relative velocity between the sinking SMBH becomes highly subsonic and the mass contained within their orbit falls below the mass of the binary itself, rendering dynamical friction ineffective. In this final stage, the black holes have not opened a gap as the gaseous background is highly pressurized in the center. Non-axisymmetric gas torques do not arise to restart sinking in absence of efficient dynamical friction, at variance with previous calculations using idealized equilibrium nuclear disk models. (abridged)

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 investigate the formation and evolution of the pseudobulge in "Eris", a high-resolution N-body + smoothed particle hydrodynamics (SPH) cosmological simulation that successfully reproduces a Milky Way-like massive late-type spiral in a cold dark matter (LCDM) universe. At the present epoch, Eris has a virial mass Mvir=8x10^11 Msun, a photometric stellar mass M*=3.2x10^10 Msun, a bulge-to-total ratio B/T = 0.26, and a weak nuclear bar. We find that the bulk of the pseudobulge forms quickly at high redshift via a combination of non-axisymmetric disk instabilities and tidal interactions or mergers both occurring on dynamical timescales, not through slow secular processes at lower redshift. Its subsequent evolution is not strictly secular either, and is closely intertwined with the evolution of the stellar bar. In fact, the structure that we recognize as a pseudobulge today evolves from a stellar bar that formed at high redshift, was destroyed by minor mergers at z~3, reformed shortly after, and weakened again following a steady gas inflow at z<1. The gradual dissolution of the bar ensues at z=1 and continues until the present without increasing the stellar velocity dispersion in the inner regions. In this scenario the pseudobulge is not a separate component from the inner disk in terms of formation path, rather it is the first step in the inside-out formation of the baryonic disk, in agreement with the fact that pseudobulges of massive spiral galaxies have typically a dominant old stellar population. If our simulations do indeed reproduce the formation mechanisms of massive spirals, then the progenitors of late-type galaxies should have strong bars and small photometric pseudobulges at high redshift.

We present a detailed study of the metal-enriched circum-galactic medium of a massive galaxy at z=3 using "ErisMC", a new cosmological hydrodynamic "zoom-in" simulations of a disk galaxy with mass comparable to the Milky Way. The run adopts a blastwave scheme for supernova feedback that generates galactic outflows without explicit wind particles, a star formation recipe based on a high gas density threshold, and high temperature metal cooling. ErisMC's main progenitor at z=3 resembles a "Lyman break" galaxy of mass M_vir=2.4e11 M_sun, virial radius R_vir=48 kpc, and star formation rate 18 M_sun/yr, and its metal-enriched CGM extends to 200 (physical) kpc from its center. Approximately 41, 9, and 50 percent of all gas-phase metals at z=3 are locked in a hot (T> 3e5 K), warm (3e5 >T> 3e4 K), and cold (T< 3e4 K) medium, respectively. We identify three sources of heavy elements: 1) the main host, responsible for 60% of all the metals found within 3R_vir; 2) its satellite progenitors, responsible for 28% of all the metals within 3R_vir, and for only 5% of those beyond 3R_vir; and nearby dwarfs, which give origin to 12% of all the metals within 3R_vir and 95% of those beyond 3R_vir. Late (z<5) galactic "superwinds" account for only 9% of all the metals observed beyond 2R_vir, the bulk having been released at redshifts 5< z < 8 by early star formation and outflows. In the CGM, lower overdensities are typically enriched by `older', colder metals. Heavy elements are accreted onto Eris along filaments via low-metallicity cold inflows, and are ejected hot via galactic outflows at a few hundred km/s. The outflow mass-loading factor is of order unity for the main halo, but can exceed 10 for nearby dwarfs. We stress that our "zoom-in" simulation focuses on the CGM of a single massive system and cannot describe the enrichment history of the IGM as a whole. (abridged)

We present 3D smoothed particle hydrodynamics simulations of the collapse of clumps formed through gravitational instability in the outer part of a protoplanetary disc. The initial conditions are taken directly from a global disc simulation, and a realistic equation of state is used to follow the clumps as they contract over several orders of magnitude in density, approaching the molecular hydrogen dissociation stage. The effects of clump rotation, asymmetries, and radiative cooling are studied. Rotation provides support against fast collapse, but non-axisymmetric modes develop and efficiently transport angular momentum outward, forming a circumplanetary disc. This transport helps the clump reach the dynamical collapse phase, resulting from molecular hydrogen dissociation, on a thousand-year timescale, which is smaller than timescales predicted by some previous spherical 1D collapse models. Extrapolation to the threshold of the runaway hydrogen dissociation indicates that the collapse timescales can be shorter than inward migration timescales, suggesting that clumps could survive tidal disruption and deliver a proto-gas giant to distances of even a few AU from the central star.

The point-like X-ray source HLX-1 is the brightest known ultraluminous X-ray source and likely the strongest intermediate-mass black hole candidate. HLX-1 is hosted in the S0 galaxy ESO 243-49, but offset with respect to the nucleus, and its optical counterpart was identified with a massive star cluster. In this paper, we study, through N-body/smoothed particle hydrodynamics simulations, the scenario where ESO 243-49 is undergoing (or just underwent) a minor merger with a gas-rich low-mass late-type galaxy. The simulations suggest that the observed star formation rate (SFR) in ESO 243-49 is a consequence of the interaction and that the companion galaxy already underwent the second pericentre passage. We propose that the counterpart of HLX-1 coincides with the nucleus (and possibly with the nuclear star cluster) of the secondary galaxy. We estimate that, if the minor merger scenario is correct, the number density of X-ray sources similar to HLX-1 is ~10^-6 Mpc^-3.

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 \rho \propto r^{-\alpha} for the dwarf DM halos, ranging from core-like (\alpha = 0.2) to cuspy (\alpha = 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 (\rho \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.

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.

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.

The formation of the massive young stars surrounding SgrA* is still an open question. In this paper, we simulate the infall of a turbulent molecular cloud towards the Galactic Center (GC). We adopt two different cloud masses (4.3x10^4 and 1.3x10^5 solar masses). We run five simulations: the gas is assumed to be isothermal in four runs, whereas radiative cooling is included in the fifth run. In all the simulations, the molecular cloud is tidally disrupted, spirals towards the GC, and forms a small, dense and eccentric disk around SgrA*. With high resolution simulations, we follow the fragmentation of the gaseous disk. Star candidates form in a ring at ~0.1-0.4 pc from the super-massive black hole (SMBH) and have moderately eccentric orbits (~0.2-0.4), in good agreement with the observations. The mass function of star candidates is top-heavy only if the local gas temperature is high (>~100 K) during the star formation and if the parent cloud is sufficiently massive (>~10^5 solar masses). Thus, this study indicates that the infall of a massive molecular cloud is a viable scenario for the formation of massive stars around SgrA*, provided that the gas temperature is kept sufficiently high (>~100 K).

Massive compact galaxies seem to be more common at high redshift than in the local universe, especially in denser environments. To investigate the fate of such massive galaxies identified at z~2 we analyse the evolution of their properties in three cosmological hydrodynamical simulations that form virialised galaxy groups of mass ~10^13 Msun hosting a central massive elliptical/S0 galaxy by redshift zero. We find that at redshift ~2 the population of galaxies with M_*> 2 10^10 Msun is diverse in terms of mass, velocity dispersion, star formation and effective radius, containing both very compact and relatively extended objects. In each simulation all the compact satellite galaxies have merged into the central galaxy by redshift 0 (with the exception of one simulation where one of such satellite galaxy survives). Satellites of similar mass at z = 0 are all less compact than their high redshift counterparts. They form later than the galaxies in the z = 2 sample and enter the group potential at z < 1, when dynamical friction times are longer than the Hubble time. Also, by z = 0 the central galaxies have increased substantially their characteristic radius via a combination of in situ star formation and mergers. Hence in a group environment descendants of compact galaxies either evolve towards larger sizes or they disappear before the present time as a result of the environment in which they evolve. Since the group-sized halos that we consider are representative of dense environments in the LambdaCDM cosmology, we conclude that the majority of high redshift compact massive galaxies do not survive until today as a result of the environment.

We investigate the formation of RE galaxies (i.e. of collisional ring galaxies with an empty ring), with N-body/SPH simulations. The simulations employ a recipe for star formation (SF) and feedback that has been shown to be crucial to produce realistic galaxies in a cosmological context. We show that RE galaxies can form via off-centre collisions (i.e. with a non-zero impact parameter), even for small inclination angles. The ring can be either a complete ring or an arc, depending on the initial conditions (especially on the impact parameter). In our simulations, the nucleus of the target galaxy is displaced from the dynamical centre of the galaxy and is buried within the ring, as a consequence of the off-centre collision. We find that the nucleus is not vertically displaced from the plane of the ring. We study the kinematics of the ring, finding agreement with the predictions by the analytic theory. The SF history of the simulated galaxies indicates that the interaction enhances the SF rate. We compare the results of our simulations with the observations of Arp 147, that is the prototype of RE galaxies.

We present the first investigation of the Blue Straggler star (BSS) population in two isolated dwarf spheroidal galaxies of the Local Group, Cetus and Tucana. Deep HST/ACS photometry allowed us to identify samples of 940 and 1214 candidates, respectively. The analysis of the star formation histories of the two galaxies suggests that both host a population of BSSs. Specifically, if the BSS candidates are interpreted as young main sequence stars, they do not conform to their galaxy's age-metallicity relationship. The analysis of the luminosity function and the radial distributions support this conclusion, and suggest a non-collisional mechanism for the BSS formation, from the evolution of primordial binaries. This scenario is also supported by the results of new dynamical simulations presented here. Both galaxies coincide with the relationship between the BSS frequency and the absolute visual magnitude Mv found by Momany et al (2007). If this relationship is confirmed by larger sample, then it could be a valuable tool to discriminate between the presence of BSSs and galaxies hosting truly young populations.

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.

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 present an analysis of the neutral hydrogen (HI) properties of a fully cosmological hydrodynamical dwarf galaxy, run with varying simulation parameters. As reported by Governato et al. (2010), the high resolution, high star formation density threshold version of this galaxy is the first simulation to result in the successful reproduction of a (dwarf) spiral galaxy without any associated stellar bulge. We have set out to compare in detail the HI distribution and kinematics of this simulated bulgeless disk with what is observed in a sample of nearby dwarfs. To do so, we extracted the radial gas density profiles, velocity dispersion (e.g., velocity ellipsoid, turbulence), and the power spectrum of structure within the cold interstellar medium from the simulations. The highest resolution dwarf, when using a high density star formation threshold comparable to densities of giant molecular clouds, possesses bulk characteristics consistent with those observed in nature, though the cold gas is not as radially extended as that observed in nearby dwarfs, resulting in somewhat excessive surface densities. The lines-of-sight velocity dispersion radial profiles have values that are in good agreement with observed dwarf galaxies, but due to the fact that only the streaming velocities of particles are tracked, a correction to include the thermal velocities can lead to profiles that are quite flat. The ISM power spectra of the simulations appear to possess more power on smaller spatial scales than that of the SMC. We conclude that unavoidable limitations remain due to the unresolved physics of star formation and feedback within pc-scale molecular clouds.

We present high-resolution 3D smoothed particle hydrodynamics simulations of the formation and evolution of protostellar discs in a turbulent molecular cloud. Using a piecewise polytropic equation of state, we perform two sets of simulations. In both cases we find that isolated systems undergo a fundamentally different evolution than members of binary or multiple systems. When formed, isolated systems must accrete mass and increase their specific angular momentum, leading to the formation of massive, extended discs, which undergo strong gravitational instabilities and are susceptible to disc fragmentation. Fragments with initial masses of 5.5 M_jup, 7.4 M_jup and 12 M_jup are produced in our simulations. In binaries and small clusters, we observe that due to competition for material from the parent core, members do not accrete significant amounts of high specific angular momentum gas relative to isolated systems. We find that discs in multiple systems are strongly self-gravitating but that they are stable against fragmentation due to disc truncation and mass profile steeping by tides, accretion of high specific angular momentum gas by other members, and angular momentum being redirected into members' orbits. In general, we expect disc fragmentation to be less likely in clusters and to be more a feature of isolated systems.

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.

The physical mechanisms and timescales that determine the morphological signatures and the quenching of star formation of typical (~L*) elliptical galaxies are not well understood. To address this issue, we have simulated the formation of a group of galaxies with sufficient resolution to track the evolution of gas and stars inside about a dozen galaxy group members over cosmic history. Galaxy groups, which harbor many elliptical galaxies in the universe, are a particularly promising environment to investigate morphological transformation and star formation quenching, due to their high galaxy density, their relatively low velocity dispersion, and the presence of a hot intragroup medium. Our simulation reproduces galaxies with different Hubble morphologies and, consequently, enables us to study when and where the morphological transformation of galaxies takes place. The simulation does not include feedback from active galactic nuclei showing that it is not an essential ingredient for producing quiescent, red elliptical galaxies in galaxy groups. Ellipticals form, as suspected, through galaxy mergers. In contrast with what has often been speculated, however, these mergers occur at z>1, before the merging progenitors enter the virial radius of the group and before the group is fully assembled. The simulation also shows that quenching of star formation in the still star-forming elliptical galaxies lags behind their morphological transformation, but, once started, is taking less than a billion years to complete. As long envisaged the star formation quenching happens as the galaxies approach and enter the finally assembled group, due to quenching of gas accretion and (to a lesser degree) stripping. A similar sort is followed by unmerged, disk galaxies, which, as they join the group, are turned into the red-and-dead disks that abound in these environments.

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.

We first review the results of the tidal stirring model for the transformation of gas-rich dwarf irregulars into dwarf spheroidals, which turns rotationally supported stellar systems into pressure supported ones. We emphasize the importance of the combined effect of ram pressure stripping and heating from the cosmic ultraviolet background in removing the gas and converting the object into a gas poor system as dSphs. We discuss how the timing of infall of dwarfs into the primary halo determines the final mass-to-light ratio and star formation history. Secondly we review the results of recent cosmological simulations of the formation of gas-rich dwarfs. These simulations are finally capable to produce a realistic object with no bulge, an exponential profile and a slowly rising rotation curve. The result owes to the inclusion of an inhomogeneous ISM and a star formation scheme based on regions having the typical density of molecular cloud complexes. Supernovae-driven winds become more effective in such mode, driving low angular momentum baryons outside the virial radius at high redshift and turning the dark matter cusp into a core. Finally we show the first tidal stirring experiments adopting dwarfs formed in cosmological simulations as initial conditions. The latter are gas dominated and have have turbulent thick gaseous and stellar disks disks that cannot develop strong bars, yet they are efficiently heated into spheroids by tidal shocks.

We make a direct comparison of the derived dark matter (DM) distributions between hydrodynamical simulations of dwarf galaxies assuming a LCDM cosmology and the observed dwarf galaxies sample from the THINGS survey in terms of (1) the rotation curve shape and (2) the logarithmic inner density slope alpha of mass density profiles. The simulations, which include the effect of baryonic feedback processes, such as gas cooling, star formation, cosmic UV background heating and most importantly physically motivated gas outflows driven by supernovae (SNe), form bulgeless galaxies with DM cores. We show that the stellar and baryonic mass is similar to that inferred from photometric and kinematic methods for galaxies of similar circular velocity. Analyzing the simulations in exactly the same way as the observational sample allows us to address directly the so-called "cusp/core" problem in the LCDM model. We show that the rotation curves of the simulated dwarf galaxies rise less steeply than CDM rotation curves and are consistent with those of the THINGS dwarf galaxies. The mean value of the logarithmic inner density slopes alpha of the simulated galaxies' dark matter density profiles is ~ -0.4 +- 0.1, which shows good agreement with \alpha = -0.29 +- 0.07 of the THINGS dwarf galaxies. The effect of non-circular motions is not significant enough to affect the results. This confirms that the baryonic feedback processes included in the simulations are efficiently able to make the initial cusps with \alpha ~ -1.0 to -1.5 predicted by dark-matter-only simulations shallower, and induce DM halos with a central mass distribution similar to that observed in nearby dwarf galaxies.

We investigate the evolution of grains composed of an ice shell surrounding an olivine core as they pass through a spiral shock in a protoplanetary disk. We use published three-dimensional radiation-hydrodynamics simulations of massive self-gravitating protoplanetary disks to extract the thermodynamics of spiral shocks in the region between 10 and 20 AU from the central star. As the density wave passes, it heats the grains, causing them to lose their ice shell and resulting in a lowering of the grain opacity. In addition, since grains of different sizes will have slightly different temperatures, there will be a migration of ice from the hotter grains to the cooler ones. The rate of migration depends on the temperature of the background gas, so a grain distribution that is effectively stable for low temperatures, can undergo an irreversible change in opacity if the gas is temporarily heated to above $\sim 150$\,K. We find that the opacity can drop more, and at a significantly faster rate throughout the spiral shocks relative to the prediction of standard dust grains models adopted in hydrodynamical calculations of protoplanetary disks. This would lead to faster gas cooling within spiral arms. We discuss the implications of our results on the susceptibility of disk fragmentation into sub-stellar objects at distances of a few tens of astronomical units.

We present an analysis of the star formation history (SFH) of the transition-type (dIrr/dSph) Local Group galaxy \objectname[]{LGS-3} (Pisces) based on deep photometry obtained with the {\it Advanced Camera for Surveys} onboard the {\it Hubble Space Telescope}. Our analysis shows that the SFH of \objectname[]{LGS-3} is dominated by a main episode $\sim 11.7$ Gyr ago with a duration of $\sim$ 1.4 Gyr which formed $\sim 90%$ of the stars. Subsequently, \objectname[]{LGS-3} continued forming stars until the present, although at a much lower rate. The lack of early chemical enrichment is in contrast to that observed in the isolated dSph galaxies of comparable luminosity, implying that the dSphs were more massive and subjected to more tidal stripping. We compare the SFH of \objectname[]{LGS-3} with expectations from cosmological models. Most or all the star formation was produced in \objectname[]{LGS-3} after the reionization epoch, assumed to be completed at $z\sim6$ or $\sim 12.7$ Gyr ago. The total mass of the galaxy is estimated to be between 2 and $4\times 10^8$ M$_\odot$, corresponding to circular velocities between 28 km\ s$^{-1}$ to 36 km\ s$^{-1}$. These values are close to but somewhat above the limit of 30 km\ s$^{-1}$ below which the UV background is expected to prevent any star formation after reionization. Feedback from SNe associated with the initial episode of star formation (mechanical luminosity from SNe $L_w=5.3\times 10^{38}$ erg s$^{-1}$) is probably inadequate to completely blow away the gas. However, the combined effects of SN feedback and UV background heating might be expected to completely halt star formation at the reionization epoch for the low mass of \objectname[]{LGS-3}; this suggests that self-shielding is important to the early evolution of galaxies in this mass range.

We present a detailed study of the dwarf galaxy VV124, recently recognized as a isolated member of the Local Group. We have obtained deep (r=26.5) wide-field g,r photometry of individual stars with the LBT under sub-arcsec seeing conditions. The Color-Magnitude Diagram suggests that the stellar content of the galaxy is dominated by an old, metal-poor population, with a significant metallicity spread. A very clean detection of the RGB tip allows us to derive an accurate distance of D=1.3 +/- 0.1 Mpc. Combining surface photometry with star counts, we are able to trace the surface brightness profile of VV124 out to ~ 5' = 1.9 kpc radius (where mu_r=30 mag/arcsec^2), showing that it is much more extended than previously believed. Moreover, the surface density map reveals the presence of two symmetric flattened wings emanating from the central elongated spheroid and aligned with its major axis, resembling a stellar disk seen nearly edge-on. We also present HI observations obtained with the WSRT, the first ever of this object. A total amount of 10^6 M_sun of HI gas is detected in VV124. Compared to the total luminosity, this gives a value of M_HI/L_V=0.11, which is particularly low for isolated Local Group dwarfs. The spatial distribution of the gas does not correlate with the observed stellar wings. The systemic velocity of the HI in the region superposed to the stellar main body of the galaxy is V_h=-25 km/s. The velocity field shows substructures typical of galaxies of this size but no sign of rotation. The HI spectra indicates the presence of a two-phase interstellar medium, again typical of many dwarf galaxies.

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.

We use 1 kpc resolution cosmological AMR simulations of a Virgo-like galaxy cluster to investigate the effect of feedback from supermassive black holes (SMBH) on the mass distribution of dark matter, gas and stars. We compared three different models: (i) a standard galaxy formation model featuring gas cooling, star formation and supernovae feedback, (ii) a "quenching" model for which star formation is artificially suppressed in massive halos and finally (iii) the recently proposed AGN feedback model of Booth & Schaye (2009). Without AGN feedback (even in the quenching case), our simulated cluster suffers from a strong overcooling problem, with a stellar mass fraction significantly above observed values in M87. The baryon distribution is highly concentrated, resulting in a strong adiabatic contraction (AC) of dark matter. With AGN feedback, on the contrary, the stellar mass in the bright central galaxy (BCG) lies below observational estimates and the overcooling problem disappears. The stellar mass of the BCG is seen to increase with increasing mass resolution, suggesting that our stellar masses converges to the correct value from below. The gas and total mass distributions are in striking agreement with observations. We also find a slight deficit (~10%) of baryons at the virial radius, due to the effect of AGN-driven shock waves pushing gas to Mpc scales and beyond. This baryon deficit results in a slight adiabatic expansion of the dark matter distribution, that can be explained quantitatively by AC theory.

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.

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.

We present a detailed study of the star formation history (SFH) of the Tucana dwarf spheroidal galaxy. High quality, deep HST/ACS data, allowed us to obtain the deepest color-magnitude diagram to date, reaching the old main sequence turnoff (F814 ~ 29) with good photometric accuracy. Our analysis, based on three different SFH codes, shows that Tucana is an old and metal-poor stellar system, which experienced a strong initial burst of star formation at a very early epoch (~ 13 Gyr ago) which lasted a maximum of 1 Gyr (sigma value). We are not able to unambiguously answer the question of whether most star formation in Tucana occurred before or after the end of the reionization era, and we analyze alternative scenarios that may explain the transformation of Tucana from a gas-rich galaxy into a dSph. Current measurements of its radial velocity do not preclude that Tucana may have crossed the inner regions of the Local Group once, and so gas stripping by ram pressure and tides due to a close interaction cannot be ruled out. On the other hand, the high star formation rate measured at early times may have injected enough energy into the interstellar medium to blow out a significant fraction of the initial gas content. Gas that is heated but not blown out would also be more easily stripped via ram pressure. We compare the SFH inferred for Tucana with that of Cetus, the other isolated LG dSph galaxy in the LCID sample. We show that the formation time of the bulk of star formation in Cetus is clearly delayed with respect to that of Tucana. This reinforces the conclusion of Monelli et al. (2010) that Cetus formed the vast majority of its stars after the end of the reionization era implying, therefore, that small dwarf galaxies are not necessarily strongly affected by reionization, in agreement with many state-of-the-art cosmological models. [abridged]

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

Massive young clusters (YCs) are expected to host intermediate-mass black holes (IMBHs) born via runaway collapse. These IMBHs are likely in binaries and can undergo mergers with other compact objects, such as stellar mass black holes (BHs) and neutron stars (NSs). We derive the frequency of such mergers starting from information available in the Local Universe. Mergers of IMBH-NS and IMBH-BH binaries are sources of gravitational waves (GWs), which might allow us to reveal the presence of IMBHs. We thus examine their detectability by current and future GW observatories, both ground- and space-based. In particular, as representative of different classes of instruments we consider Initial and Advanced LIGO, the Einstein gravitational-wave Telescope (ET) and the Laser Interferometer Space Antenna (LISA). We find that IMBH mergers are unlikely to be detected with instruments operating at the current sensitivity (Initial LIGO). LISA detections are disfavored by the mass range of IMBH-NS and IMBH-BH binaries: less than one event per year is expected to be observed by such instrument. Advanced LIGO is expected to observe a few merger events involving IMBH binaries in a 1-year long observation. Advanced LIGO is particularly suited for mergers of relatively light IMBHs (~100 Msun) with stellar mass BHs. The number of mergers detectable with ET is much larger: tens (hundreds) of IMBH-NS (IMBH-BH) mergers might be observed per year, according to the runaway collapse scenario for the formation of IMBHs. We note that our results are affected by large uncertainties, produced by poor observational constraints on many of the physical processes involved in this study, such as the evolution of the YC density with redshift.[abridged]

(Abridged) We use N-body simulations to study the effects that a divergent (i.e. "cuspy") dark matter (DM) profile introduces on the tidal evolution of dwarf spheroidal galaxies (dSphs). Our models assume cosmologically-motivated initial conditions where dSphs are DM-dominated systems on eccentric orbits about a host galaxy composed of a dark halo and a baryonic disc. We find that the resilience of dSphs to tidal stripping is extremely sensitive to the halo cuspiness; whereas dwarfs with a cored profile can be easily destroyed by the host disc, those with cusps always retain a bound remnant. For a given halo profile the evolution of the structural parameters as driven by tides is controlled solely by the total amount of mass lost. This information is used to construct a semi-analytic code that simulates the hierarchical build-up of spiral galaxies assuming different halo profiles and disc masses. We find that tidal encounters with discs tend to decrease the average mass of satellites at all galactocentric radii. Interestingly, satellites accreted before re-ionization (z>6), which may be singled out by anomalous metallicity patterns, survive only if haloes are cuspy. We show that the size-mass relation established from Milky Way (MW) dwarfs strongly supports the presence of cusps in the majority of these systems, as cored models systematically underestimate the masses of the known Ultra-Faint dSphs. Our models also indicate that a massive M31 disc may explain why many of its dSphs fall below the size-mass relationship derived from MW dSphs. We use our models to constrain the mass threshold below which star formation is suppressed in DM haloes, finding that luminous satellites must be accreted with masses above 10^8--10^9 M_sol in order to explain the size-mass relation observed in MW dwarfs.

Making robust predictions for the phase space distribution of dark matter at the solar neighbourhood is vital for dark matter direct detection experiments. To date, almost all such predictions have been based on simulations that model the dark matter alone. Here, we use three cosmological hydrodynamics simulations of bright, disc dominated galaxies to include the effects of baryonic matter self-consistently for the first time. We find that the addition of baryonic physics drastically alters the dark matter profile in the vicinity of the Solar neighbourhood. A stellar/gas disc, already in place at high redshift, causes merging satellites to be dragged preferentially towards the disc plane where they are torn apart by tides. This results in an accreted dark matter disc that contributes ~0.25 - 1.5 times the non-rotating halo density at the solar position. The dark disc, unlike dark matter streams, is an equilibrium structure that must exist in disc galaxies that form in a hierarchical cosmology. Its low rotation lag with respect to the Earth significantly boosts WIMP capture in the Earth and Sun, boosts the annual modulation signal, and leads to distinct variations in the flux as a function of recoil energy that allow the WIMP mass to be determined.

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.

We use deep HST/ACS observations to calculate the star formation history (SFH) of the Cetus dwarf spheroidal (dSph) galaxy. Our photometry reaches below the oldest main sequence turn-offs, which allows us to estimate the age and duration of the main episode of star formation in Cetus. This is well approximated by a single episode that peaked roughly 12+/-0.5 Gyr ago and lasted no longer than about 1.9+/-0.5 Gyr (FWHM). Our solution also suggests that essentially no stars formed in Cetus during the past 8 Gyrs. This makes Cetus' SFH comparable to that of the oldest Milky Way dSphs. Given the current isolation of Cetus in the outer fringes of the Local Group, this implies that Cetus is a clear outlier in the morphology-Galactocentric distance relation that holds for the majority of Milky Way dwarf satellites. Our results also show that Cetus continued forming stars through z ~ 1, long after the Universe was reionized, and that there is no clear signature of the epoch of reionization in Cetus' SFH. We discuss briefly the implications of these results for dwarf galaxy evolution models. Finally, we present a comprehensive account of the data reduction and analysis strategy adopted for all galaxies targeted by the LCID (Local Cosmology from Isolated Dwarfs) project. We employ two different photometry codes (DAOPHOT/ALLFRAME and DOLPHOT), three different SFH reconstruction codes (IAC-pop/MinnIAC, MATCH, COLE), and two stellar evolution libraries (BaSTI and Padova/Girardi), allowing for a detailed assessment of the modeling and observational uncertainties.

We explore the initial conditions for fragments in the extended regions ($r\gtrsim 50$ AU) of gravitationally unstable disks. We combine analytic estimates for the fragmentation of spiral arms with 3D SPH simulations to show that initial fragment masses are in the gas giant regime. These initial fragments will have substantial angular momentum, and should form disks with radii of a few AU. We show that clumps will survive for multiple orbits before they undergo a second, rapid collapse due to H$_2$ dissociation and that it is possible to destroy bound clumps by transporting them into the inner disk. The consequences of disrupted clumps for planet formation, dust processing, and disk evolution are discussed. We argue that it is possible to produce Earth-mass cores in the outer disk during the earliest phases of disk 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.

We trace the evolution of central galaxies in three ~10^13 M_sun galaxy groups simulated at high resolution in cosmological hydrodynamical simulations. The evolution in the group potential leads, at z=0, to central galaxies that are massive, gas-poor early-type systems supported by stellar velocity dispersion resembling either elliptical or S0 galaxies. Their z~2-2.5 main progenitors are massive M* ~ 3-10 x 10^10 M_sun, star forming (20-60 M_sun/yr) galaxies which host substantial reservoirs of cold gas (~5 x 10^9 M_sun) in extended gas disks. Our simulations thus show that star forming galaxies observed at z~2 are likely the main progenitors of central galaxies in galaxy groups at z=0. Their central stellar densities stay approximately constant from z~1.5 down to z=0. Instead, the galaxies grow inside-out, by acquiring a stellar envelope outside the innermost ~2 kpc. Consequently the density within the effective radius decreases by up to two orders of magnitude. Both major and minor mergers contribute to most of the mass accreted outside the effective radius and thus drive the evolution of the half-mass radii. In one of the three simulated groups the short central cooling time leads to a dramatic rejuvenation of the central group galaxy at z<1, affecting its morphology, kinematics and colors. This episode is eventually terminated by a group-group merger. Our simulations demonstrate that, in galaxy groups, the interplay between halo mass assembly, galaxy merging and gas accretion has a substantial influence on the star formation histories and z=0 morphologies of central galaxies.[Abridged]

For almost two decades the properties of "dwarf" galaxies have challenged the Cold Dark Matter (CDM) paradigm of galaxy formation. Most observed dwarf galaxies consists of a rotating stellar disc embedded in a massive DM halo with a near constant-density core. Yet, models based on the CDM scenario invariably form galaxies with dense spheroidal stellar "bulges" and steep central DM profiles, as low angular momentum baryons and DM sink to the center of galaxies through accretion and repeated mergers. Processes that decrease the central density of CDM halos have been identified, but have not yet reconciled theory with observations of present day dwarfs. This failure is potentially catastrophic for the CDM model, possibly requiring a different DM particle candidate. This Letter presents new hydrodynamical simulations in a Lambda$CDM framework where analogues of dwarf galaxies, bulgeless and with a shallow central DM profile, are formed. This is achieved by resolving the inhomogeneous interstellar medium, resulting in strong outflows from supernovae explosions which remove low angular momentum gas. This inhibits the formation of bulges and decreases the dark-matter density to less than half within the central kiloparsec. Realistic dwarf galaxies are thus shown to be a natural outcome of galaxy formation in the CDM scenario.

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 review numerical work carried out over the last decade on the role of environmental mechanisms in shaping the nature of the faintest galaxies known, dwarf spheroidals (dSphs). In particular we discuss a model in which dSphs originated from gas dominated, disky dwarfs that were accreted by massive galaxies several billions of years ago. Dwarfs accreting at $z > 1$, when the cosmic ultraviolet ionizing flux was much higher than today and was thus able to keep the gas in the dwarfs warm and diffuse, were rapidly stripped of their baryons via ram pressure and tidal forces, producing very dark matter dominated objects with truncated star formation histories, such as the Draco dSph. The low star formation efficiency expected in such low-metallicity objects prior to their infall was crucial for keeping their disks gas dominated until stripping took over. {\it Therefore gas stripping along with inefficient star formation provides a new feedback mechanism, alternative to photoevaporation or supernovae feedback, playing a crucial role in dwarf galaxy formation and evolution}. Tidally induced instabilities, termed "tidal stirring", turned the surviving stellar disk into a spheroid. Moreover, we examine the properties of the newly discovered ultra-faint dSphs in light of this scenario and argue that they likely belong to a different population of lower mass dwarf satellites. These were mostly shaped by reionization rather than by environmental mechanisms and are thus good candidates for being "reionization fossils". We discuss implications of the morphological transformation scenario on the substructure problem. (Abridged)

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': 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.

Using high resolution SPH simulations in a fully cosmological Lambda CDM context we study the formation of a bright disk dominated galaxy that originates from a "wet" major merger at z=0.8. The progenitors of the disk galaxy are themselves disk galaxies that formed from early major mergers between galaxies with blue colors. A substantial thin stellar disk grows rapidly following the last major merger and the present day properties of the final remnant are typical of early type spiral galaxies, with an i band B/D ~0.65, a disk scale length of 7.2 kpc, g-r = 0.5 mag, an HI line width (W_{20}/2) of 238 km/sec and total magnitude i = -22.4. The key ingredients for the formation of a dominant stellar disk component after a major merger are: i) substantial and rapid accretion of gas through cold flows followed at late times by cooling of gas from the hot phase, ii) supernova feedback that is able to partially suppress star formation during mergers and iii) relative fading of the spheroidal component. The gas fraction of the progenitors' disks does not exceed 25% at z<3, emphasizing that the continuous supply of gas from the local environment plays a major role in the regrowth of disks and in keeping the galaxies blue. The results of this simulation alleviate the problem posed for the existence of disk galaxies by the high likelihood of interactions and mergers for galaxy sized halos at relatively low z.

Computational astrophysics has undergone unprecedented development over the last decade, becoming a field of its own. The challenge ahead of us will involve increasingly complex multi-scale simulations. These will bridge the gap between areas of astrophysics such as star and planet formation, or star formation and galaxy formation, that have evolved separately until today. A global knowledge of the physics and modeling techniques of astrophysical simulations is thus an important asset for the next generation of modelers. With the aim at fostering such a global approach, we present the Special Issue on Computational Astrophysics for the Advanced Science Letters (http://www.aspbs.com/science.htm). The Advanced Science Letters (ASL) is a new multi-disciplinary scientific journal which will cover extensively computational astrophysics and cosmology, and will act as a forum for the presentation and discussion of novel work attempting to connect different research areas. This Special Issue collects 9 reviews on 9 key topics of modern astrophysics and cosmology in which simulations have had a pivotal contribution. From asteroid collisions to the formation of stars, from the physics and origin of galaxy clusters to large-scale phenomena such as the reionization of the Universe, it is a showcase of state-of-the-art numerical work in a diverse range of topics. The 9 reviews are available on today's listing of the astro-ph archive (arXiv id numbers are provided in this foreword).

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, \sigma, 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} \sigma 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.

(abridged) We present the first study of the variable star populations in the isolated dwarf spheroidal galaxies (dSph) Cetus and Tucana. Based on Hubble Space Telescope images obtained with the Advanced Camera for Surveys in the F475W and F814W bands, we identified 180 and 371 variables in Cetus and Tucana, respectively. The vast majority are RR Lyrae stars. In Cetus we also found three anomalous Cepheids, four candidate binaries and one candidate long-period variable (LPV), while six anomalous Cepheids and seven LPV candidates were found in Tucana. Of the RR Lyrae stars, 147 were identified as fundamental mode (RRab) and only eight as first-overtone mode (RRc) in Cetus, with mean periods of 0.614 and 0.363 day, respectively. In Tucana we found 216 RRab and 82 RRc giving mean periods of 0.604 and 0.353 day. These values place both galaxies in the so-called Oosterhoff Gap, as is generally the case for dSph. We calculated the distance modulus to both galaxies using different approaches based on the properties of RRab and RRc, namely the luminosity-metallicity and period-luminosity-metallicity relations, and found values in excellent agreement with previous estimates using independent methods: (m-M)_{0,Cet}=24.46+-0.12 and (m-M)_{0,Tuc}=24.74+-0.12, corresponding to 780+-40 kpc and 890+-50 kpc. We also found numerous RR Lyrae variables pulsating in both modes simultaneously (RRd): 17 in Cetus and 60 in Tucana. Tucana is, after Fornax, the second dSph in which such a large fraction of RRd (~17%) has been observed. We provide the photometry and pulsation parameters for all the variables, and compare the latter with values from the literature for well-studied dSph of the Local Group and Galactic globular clusters.

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.

We announce the discovery of a new Milky Way satellite Segue 2 found in the data of the Sloan Extension for Galactic Understanding and Exploration (SEGUE). We followed this up with deeper imaging and spectroscopy on the Multiple Mirror Telescope. From this, we derive a luminosity of M_v = -2.5, a half-light radius of 34 pc and a systemic velocity of -40$ km/s. Our MMT data also provides evidence for a stream around Segue 2 at a similar heliocentric velocity, and the SEGUE data show that it is also present in neighboring fields. We resolve the velocity dispersion of Segue 2 as 3.4 km/s and the possible stream as about 7 km/s. This object shows points of comparison with other recent discoveries, Segue 1, Boo II and Coma. We speculate that all four objects may be representatives of a population of satellites of satellites -- survivors of accretion events that destroyed their larger but less dense parents. They are likely to have formed at redshifts z > 10 and are good candidates for fossils of the reionization epoch.

Predicting the local flux of dark matter particles is vital for dark matter direct detection experiments. To date, such predictions have been based on simulations that model the dark matter alone. Here we include the influence of the baryonic matter for the first time. We use two different approaches. Firstly, we use dark matter only simulations to estimate the expected merger history for a Milky Way mass galaxy, and then add a thin stellar disc to measure its effect. Secondly, we use three cosmological hydrodynamic simulations of Milky Way mass galaxies. In both cases, we find that a stellar/gas disc at high redshift (z~1) causes merging satellites to be preferentially dragged towards the disc plane. This results in an accreted dark matter disc that contributes ~0.25 - 1 times the non-rotating halo density at the solar position. An associated thick stellar disc forms with the dark disc and shares a similar velocity distribution. If these accreted stars can be separated from those that formed in situ, future astronomical surveys will be able to infer the properties of the dark disc from these stars. The dark disc, unlike dark matter streams, is an equilibrium structure that must exist in disc galaxies that form in a hierarchical cosmology. Its low rotation lag with respect to the Earth significantly boosts WIMP capture in the Earth and Sun, increases the likelihood of direct detection at low recoil energy, boosts the annual modulation signal, and leads to distinct variations in the flux as a function of recoil energy that allow the WIMP mass to be determined (see contribution from T. Bruch this volume).

Observations of turbulent velocity dispersions in the HI component of galactic disks show a characteristic floor in galaxies with low star formation rates and within individual galaxies the dispersion profiles decline with radius. We carry out several high resolution adaptive mesh simulations of gaseous disks embedded within dark matter haloes to explore the roles of cooling, star-formation, feedback, shearing motions and baryon fraction in driving turbulent motions. In all simulations the disk slowly cools until gravitational and thermal instabilities give rise to a multi-phase medium in which a large population of dense self-gravitating cold clouds are embedded within a warm gaseous phase that forms through shock heating. The diffuse gas is highly turbulent and is an outcome of large scale driving of global non-axisymmetric modes as well as cloud-cloud tidal interactions and merging. At low star-formation rates these processes alone can explain the observed HI velocity dispersion profiles and the characteristic value of ~10 km/s observed within a wide range of disk galaxies. Supernovae feedback creates a significant hot gaseous phase and is an important driver of turbulence in galaxies with a star-formation rate per unit area >10^-3 M_sun/yr/kpc^2.

The final inspiral and coalescence of a black hole binary can produce highly beamed gravitational wave radiation. To conserve linear momentum, the black hole remnant can recoil with "kick" velocity as high as 4000 km/s. We present two sets of full N-body simulations of recoiling massive black holes (MBH) in high-resolution, non-axisymmetric potentials. The host to the first set of simulations is the main halo of the Via Lactea I simulation (Diemand et al. 2007). The nature of the resulting orbits is investigated through a numerical model where orbits are integrated assuming an evolving, triaxial NFW potential, and dynamical friction is calculated directly from the velocity dispersion along the major axes of the main halo of Via Lactea I. By comparing the triaxial case to a spherical model, we find that the wandering time spent by the MBH is significantly increased due to the asphericity of the halo. For kicks larger than 200 km/s, the remnant MBH does not return to the inner 200 pc within 1 Gyr, a timescale an order of magnitude larger than the upper limit of the estimated QSO lifetime. The second set of simulations is run using the outcome of a high-resolution gas-rich merger (Mayer et al. 2007) as host potential. In this case, a recoil velocity of 500 km/s cannot remove the MBH from the nuclear region.

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.

We use a set of high-resolution N-body simulations of binary galaxy mergers to show that the morphologies of the tidal features that are seen around a large fraction of nearby, massive ellipticals in the field, cannot be reproduced by equal-mass dissipationless mergers; rather, they are well explained by the accretion of disk-dominated galaxies. In particular, the arm- and looplike morphologies of the observed tidal debris can only be produced by the kinematically cold material of the disk components of the accreted galaxies. The tidal features that arise from such "cold-accretion" events onto a massive elliptical are visible for significantly longer timescales than the features produced by elliptical-elliptical mergers (about 1-2 Gyr vs. a few hundred million years). Mass ratios of the order of 1:10 between the accreting elliptical and the accreted disk are sufficient to match the brightness of the observed debris. Furthermore, stellar population synthesis models and simple order-of-magnitude calculations indicate that the colors of the tidal features generated in such minor cold-accretion events are relatively red, in agreement with the observations. The minor cold-accretion events that explain the presence, brightness, and structural and color properties of the tidal debris cause only a modest mass and luminosity increase in the accreting massive elliptical. These results, coupled with the relative statistical frequencies of disk- and bulge-dominated galaxies in the field, suggest that massive ellipticals assemble most of their mass well before their tidal debris forms through the accretion of relatively little, kinematically cold material rather than in very recent, dissipationless major mergers.

We study the evolution of a massive black hole pair in a rotationally supported nuclear disc. The distributions of stars and gas mimic the nuclear region of a gas-rich galaxy merger remnant. Using high-resolution SPH simulations, we follow the black hole dynamics and trace the evolution of the underlying background, until the black holes form a binary. We find that the gravitational perturbation of the pair creates a core in the disc density profile, hence decreasing the gas-dynamical drag. This leads the newly formed binary to stall at a separation of ~5 pc. In the early phases of the sinking, black holes lose memory of their initial orbital eccentricity if they co-rotate with the disc, as rotation of the gaseous background promotes circularization of the black hole orbits. Circularization is efficient until the black holes bind in a binary, though in the latest stages of the simulations a residual eccentricity > 0.1 is still present. Black holes are treated as sink particles, allowing for gas accretion. We find that accretion strongly depends on the dynamical properties of the black holes, and occurs preferentially after circularization.

(Abridged) We review the results of the first multi-scale, hydrodynamical simulations of mergers between galaxies with central supermassive black holes (SMBHs) to investigate the formation of SMBH binaries in galactic nuclei. We demonstrate that strong gas inflows produce nuclear disks at the centers of merger remnants whose properties depend sensitively on the details of gas thermodynamics. In numerical simulations with parsec-scale spatial resolution in the gas component and an effective equation of state appropriate for a starburst galaxy, we show that a SMBH binary forms very rapidly, less than a million years after the merger of the two galaxies. Binary formation is significantly suppressed in the presence of a strong heating source such as radiative feedback by the accreting SMBHs. We also present preliminary results of numerical simulations with ultra-high spatial resolution of 0.1 pc in the gas component. These simulations resolve the internal structure of the resulting nuclear disk down to parsec scales and demonstrate the formation of a central massive object (~ 10^8 Mo) by efficient angular momentum transport. This is the first time that a radial gas inflow is shown to extend to parsec scales as a result of the dynamics and hydrodynamics involved in a galaxy merger, and has important implications for the fueling of SMBHs. Due to the rapid formation of the central clump, the density of the nuclear disk decreases significantly in its outer region, reducing dramatically the effect of dynamical friction and leading to the stalling of the two SMBHs at a separation of ~1 pc. We discuss how the orbital decay of the black holes might continue in a more realistic model which incorporates star formation and the multi-phase nature of the ISM.

We model the recently published kinematic data set for Leo I dSph galaxy by fitting the solutions of the Jeans equations to the velocity dispersion and kurtosis profiles measured from the data. We demonstrate that when the sample is cleaned of interlopers the data are consistent with the assumption that mass follows light and isotropic stellar orbits with no need for an extended dark matter halo. Our interloper removal scheme does not clean the data of contamination completely, as demonstrated by the rotation curve of Leo I. When moving away from the centre of the dwarf, the rotation appears to be reversed. We interpret this behaviour using the results of an N-body simulation of a dwarf galaxy possessing some intrinsic rotation, orbiting in the Milky Way potential and show that it can be reproduced if the galaxy is viewed almost along the tidal tails so that the leading (background) tail contaminates the western part of Leo I while the trailing (foreground) tail the eastern one. We show that this configuration leads to a symmetric and Gaussian distribution of line-of-sight velocities. The simulation is also applied to test our modelling method on mock data sets. We demonstrate that when the data are cleaned of interlopers and the fourth velocity moment is used the true parameters of the dwarf are typically within 1 \sigma errors of the best-fitting parameters. Restricting the fitting to the inner part of Leo I our best estimate for the anisotropy is \beta = -0.2^{+0.3}_{-0.4} and the total mass M = (4.5 +/- 0.7) x 10^7 M_sun. The mass-to-light ratio including the errors in mass, brightness and distance is M/L_V = 8.2 +/- 4.5 solar units.

We report on the results of novel global high-resolution three-dimensional simulations of disk-planet interaction which incorporate simultaneously realistic radiation physics and the self-gravity of the gas, as well as allowing the planet to move. We find that thermodynamics and radiative physics have a remarkable effect on both migration and accretion of Jupiter mass planets. In simulations with radiative transfer adopting flux-limited diffusion, inward migration can be decreased by about 30% relative to the isothermal case, while in adiabatic runs migration nearly shuts off after a few tens of orbits. Migration varies because the relative strength of the inner and outer spiral perturbations is affected by thermodynamics, thus changing the net torque acting on the planet. Mass accretion rates on the planet can be reduced by more than an order of magnitude going from isothermal to radiative transfer and adiabatic simulations. A circumplanetary disk always forms except in adiabatic runs. With radiative transfer the disk is sub-keplerian (Vrot/Vkep ~ 0.7) owing to significant pressure support. We discuss the effect of circumplanetary disk structure on the drift of embedded dust grains and planetesimals and thus on the formation of the rocky satellites of giant planets.

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.

The formation of the massive young stars surrounding SgrA* is still an open question. In this paper, we simulate the infall of an isothermal, turbulent molecular cloud towards the Galactic Centre (GC). As it spirals towards the GC, the molecular cloud forms a small and dense disc around SgrA*. Efficient star formation (SF) is expected to take place in such a dense disc. We model SF by means of sink particles. At ~6x10^5 yr, ~6000 solar masses of stars have formed, and are confined within a thin disc with inner and outer radius of 0.06 and 0.5 pc, respectively. Thus, this preliminary study shows that the infall of a molecular cloud is a viable scenario for the formation of massive stars around SgrA*. Further studies with more realistic radiation physics and SF will be required to better constrain this intriguing scenario.

The formation of disk galaxies is one of the most outstanding problems in modern astrophysics and cosmology. We review the progress made by numerical simulations carried out on large parallel supercomputers. Recent progress stems from a combination of increased resolution and improved treatment of the astrophysical processes modeled in the simulations, such as the phenomenological description of the interstellar medium and of the process of star formation. High mass and spatial resolution is a necessary condition in order to obtain large disks comparable with observed spiral galaxies avoiding spurious dissipation of angular momentum. A realistic model of the star formation history. gas-to-stars ratio and the morphology of the stellar and gaseous component is instead controlled by the phenomenological description of the non-gravitational energy budget in the galaxy. We show that simulations of gas collapse within cold dark matter halos including a phenomenological description of supernovae blast-waves allow to obtain stellar disks with nearly exponential surface density profiles as those observed in real disk galaxies, counteracting the tendency of gas collapsing in such halos to form cuspy baryonic profiles. However, the ab-initio formation of a realistic rotationally supported disk galaxy with a pure exponential disk in a fully cosmological simulation is still an open problem. We argue that the suppression of bulge formation is related to the physics of galaxy formation during the merger of the most massive protogalactic lumps at high redshift, where the reionization of the Universe likely plays a key role. A sufficiently high resolution during this early phase of galaxy formation is also crucial to avoid artificial angular momentum loss (Abridged).

We present a new set of multi-million particle SPH simulations of the formation of disk dominated galaxies in a cosmological context. Some of these galaxies are higher resolution versions of the models already described in Governato et al (2007). To correctly compare simulations with observations we create artificial images of our simulations and from them measure photometric Bulge to Disk (B/D) ratios and disk scale lengths. We show how feedback and high force and mass resolution are necessary ingredients to form galaxies that have flatter rotation curves, larger I band disk scale lengths and smaller B/D ratios. A new simulated disk galaxy has an I-band disk scale length of 9.2 kpc and a B/D flux ratio of 0.64 (face on, dust reddened).

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 simulate the collisional formation of a ring galaxy and we integrate its evolution up to 1.5 Gyr after the interaction. About 100-200 Myr after the collision, the simulated galaxy is very similar to observed ring galaxies (e.g. Cartwheel). After this stage, the ring keeps expanding and fades. Approximately 0.5-1 Gyr after the interaction, the disc becomes very large (~100 kpc) and flat. Such extended discs have been observed only in giant low surface brightness galaxies (GLSBs). We compare various properties of our simulated galaxies (surface brightness profile, morphology, HI spectrum and rotation curve) with the observations of four well-known GLSBs (UGC6614, Malin 1, Malin 2 and NGC7589). The simulations match quite well the observations, suggesting that ring galaxies could be the progenitors of GLSBs. This result is crucial for the cold dark matter (CDM) model, as it was very difficult, so far, to explain the formation of GLSBs within the CDM scenario.

If massive black holes (BHs) are ubiquitous in galaxies and galaxies experience multiple mergers during their cosmic assembly, then BH binaries should be common albeit temporary features of most galactic bulges. Observationally, the paucity of active BH pairs points toward binary lifetimes far shorter than the Hubble time, indicating rapid inspiral of the BHs down to the domain where gravitational waves lead to their coalescence. Here, we review a series of studies on the dynamics of massive BHs in gas-rich galaxy mergers that underscore the vital role played by a cool, gaseous component in promoting the rapid formation of the BH binary. The BH binary is found to reside at the center of a massive self-gravitating nuclear disc resulting from the collision of the two gaseous discs present in the mother galaxies. Hardening by gravitational torques against gas in this grand disc is found to continue down to sub-parsec scales. The eccentricity decreases with time to zero and when the binary is circular, accretion sets in around the two BHs. When this occurs, each BH is endowed with it 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 < 1 Myr. The double nuclear point-like sources that may appear have typical separation of < 10 pc, and are likely to be embedded in the still ongoing starburst. We note that a potential threat of binary stalling, in a gaseous environment, may come from radiation and/or mechanical energy injections by the BHs. Only short-lived or sub-Eddington accretion episodes can guarantee the persistence of a dense cool gas structure around the binary necessary for continuing BH inspiral.

We review recent results of SPH simulations of gravitational instability in gaseous protoplanetary disks,emphasizing the role of thermodynamics in both isolated and binary systems. Contradictory results appeared in the literature regarding disk fragmentation at tens of AU from the central star are likely due to the different treatment of radiation physics as well as reflecting different initial conditions. Further progress on the subject requires extensive comparisons between different codes with the requirement that the same initial conditions are adopted. It is discussed how the local conditions of the disks undergoing fragmentation at $R < 25$ AU in recent SPH simulations are in rough agreement with the prediction of analytical models, with small differences being likely related to the inability of analytical models to account for the dynamics and thermodynamics of three-dimensional spiral shocks. We report that radically different adaptive hydrodynamical codes, SPH and adaptive mesh refinement (AMR), yield very similar results on disk fragmentation at comparable resolution in the simple case of an isothermal equation of state. A high number of refinements in AMR codes is necessary but not sufficient to correctly follow fragmentation, rather an initial resolution of the grid high enough to capture the wavelength of the strongest spiral modes when they are still barely nonlinear is essential. These tests represent a useful benchmark and a starting point for a forthcoming code comparison with realistic radiation physics.

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.

The dynamics of two massive black holes in a rotationally supported nuclear disc of 10^8 solar masses is explored using N-Body/SPH simulations. Gas and star particles are co-present in the disc. Described by a Mestel profile, the disc has a vertical support provided by turbulence of the gas, and by stellar velocity dispersion. A primary black hole of 4 million solar masses is placed at the centre of the disc, while a secondary black hole is set initially on an eccentric co-rotating orbit in the disc plane. Its mass is in a 1 to 1, 1 to 4, and 1 to 10 ratio, relative to the primary. With this choice, we mimic the dynamics of black hole pairs released in the nuclear region at the end of a gas-rich galaxy merger. It is found that, under the action of dynamical friction, the two black holes form a close binary in ~10 Myrs. The inspiral process is insensitive to the mass fraction in stars and gas present in the disc and is accompanied by the circularization of the orbit. We detail the gaseous mass profile bound to each black hole that can lead to the formation of two small Keplerian discs, weighing ~2 % of the black hole mass, and of size ~0.01 pc. The mass of the tightly (loosely) bound particles increases (decreases) with time as the black holes spiral into closer and closer orbits. Double AGN activity is expected to occur on an estimated timescale of ~10 Myrs, comparable to the inspiral timescale. The double nuclear point-like sources that may appear during dynamical evolution have typical separations of ~10 pc.

Supermassive black holes (SMBHs) are a ubiquitous component of the nuclei of galaxies. It is normally assumed that, following the merger of two massive galaxies, a SMBH binary will form, shrink due to 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 due to 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 following 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 a million years as a result of the gravitational drag from the gas rather than from the stars.

We review the models and results of simulations of self-gravitating, gaseous protoplanetary disks in binary star systems. These models have been calculated by three different groups with three different computational methods, two particle-based and one grid-based. We show that interactions with the companion star can affect the temperature distribution and structural evolution of disks, and discuss the implications for giant planet formation by gravitational instability as well as by core accretion. Complete consensus has not been reached yet on whether the formation of giant planets is promoted or suppressed by tidal interactions with a companion star. While systems with binary separations of order 100 AU or larger appear to behave more or less as in isolation, systems with smaller separations exhibit an increased or decreased susceptibility to fragmentation, depending on the details of thermodynamics, in particular on the inclusion or absence of artificial viscosity, and on the initial conditions. While code comparisons on identical problems need to be carried out (some of which are already in progress) to decide which computer models are more realistic, it is already clear that relatively close binary systems, with separations of order 60 AU or less, should provide strong constraints on how giant planets form in these systems.

We study the origin and properties of the population of unbound stars in the kinematic samples of dwarf spheroidal galaxies. For this purpose we have run a high resolution N-body simulation of a two-component dwarf galaxy orbiting in a Milky Way potential. The dwarf is placed on a highly eccentric orbit, its initial stellar component is in the form of an exponential disk and it has a NFW-like dark matter halo. After 10 Gyrs of evolution the dwarf produces a spheroidal stellar component and is strongly tidally stripped so that mass follows light and the stars are on almost isotropic orbits. From this final state, we create mock kinematic data sets for 200 stars by observing the dwarf in different directions. We find that when the dwarf is observed along the tidal tails the kinematic samples are strongly contaminated by unbound stars from the tails. We also study another source of possible contamination by adding stars from the Milky Way. We demonstrate that most of the unbound stars can be removed by the method of interloper rejection proposed by den Hartog & Katgert. We model the cleaned-up kinematic samples using solutions of the Jeans equation with constant mass-to-light ratio and velocity anisotropy parameter. We show that the Jeans analysis, when applied to cleaned samples, allows us to reproduce the mass and mass-to-light ratio of the dwarf with accuracy typically better than 25 percent. The analysis was applied to the new data for the Fornax dSph galaxy. We show that after careful removal of interlopers the velocity dispersion profile of Fornax can be reproduced by a model in which mass traces light with a mass-to-light ratio of 11 solar units and isotropic orbits. We demonstrate that most of the contamination in the kinematic sample of Fornax probably originates from the Milky Way. (Abridged)

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.

We present preliminary results concerning the search for short-period variable stars in Tucana and LGS3 based on very deep HST/ACS imaging. In the fraction of the observed field we studied in each galaxy, a total of 133 and 30 variables were found, respectively. For Tucana, we identified 76 of them as RR Lyrae (RRL) stars pulsating in the fundamental mode (RRab) and 32 in the first-overtone mode (RRc), as well as 2 anomalous Cepheids (AC). The mean period of the RRab and RRc is 0.59 and 0.35 days, respectively. In the case of LGS3, we found 24 RRab and 4 RRc, with mean periods 0.61and 0.39 days, respectively, plus two candidate ACs. These values place both galaxies in the Oosterhoff gap.

CCD photometry on the intermediate-band vbyCaHbeta system is presented for the metal-rich, old open cluster, NGC 6791. Preliminary analysis led to [Fe/H] above +0.4 with an anomalously high reddening and an age below 5 Gyr. A revised calibration between (b-y)_0 and [Fe/H] at a given temperature shows that the traditional color-metallicity relations underestimate the color of the turnoff stars at high metallicity. With the revised relation, the metallicity from hk and the reddening for NGC 6791 become [Fe/H] = +0.45 +/- 0.04 and E(b-y) = 0.113 +/- 0.012 or E(B-V) = 0.155 +/- 0.016. Using the same technique, reanalysis of the photometry for NGC 6253 produces [Fe/H] = +0.58 +/-0.04 and E(b-y) = 0.120 +/- 0.018 or E(B-V) = 0.160 +/- 0.025. The errors quoted include both the internal and external errors. For NGC 6791, the metallicity from m_1 is a factor of two below that from hk, a result that may be coupled to the consistently low metal abundance from DDO photometry of the cluster and the C-deficiency found from high dispersion spectroscopy. E(B-V) is the same value predicted from Galactic reddening maps. With E(B-V) = 0.15 and [Fe/H] = +0.45, the available isochrones predict an age of 7.0 +/- 1.0 Gyr and an apparent modulus of (m-M) = 13.60 +/- 0.15, with the dominant source of the uncertainty arising from inconsistencies among the isochrones. The reanalysis of NGC 6253 with the revised lower reddening confirms that on both the hk and m_1 metallicity scales, NGC 6253, while less than half the age of NGC 6791, remains at least as metal-rich as NGC 6791, if not richer.

We perform controlled N-Body/SPH simulations of disk galaxy formation by cooling a rotating gaseous mass distribution inside equilibrium cuspy spherical and triaxial dark matter halos. We systematically study the angular momentum transport and the disk morphology as we increase the number of dark matter and gas particles from 10^4 to 10^6, and decrease the gravitational softening from 2 kpc to 50 parsecs. The angular momentum transport, disk morphology and radial profiles depend sensitively on force and mass resolution. At low resolution, similar to that used in most current cosmological simulations, the cold gas component has lost half of its initial angular momentum via different mechanisms. The angular momentum is transferred primarily to the hot halo component, by resolution-dependent hydrodynamical and gravitational torques, the latter arising from asymmetries in the mass distribution. In addition, disk-particles can lose angular momentum while they are still in the hot phase by artificial viscosity. In the central disk, particles can transfer away over 99% of their initial angular momentum due to spiral structure and/or the presence of a central bar. The strength of this transport also depends on force and mass resolution - large softening will suppress the bar instability, low mass resolution enhances the spiral structure. This complex interplay between resolution and angular momentum transfer highlights the complexity of simulations of galaxy formation even in isolated haloes. With 10^6 gas and dark matter particles, disk particles lose only 10-20% of their original angular momentum, yet we are unable to produce pure exponential profiles.

We used fully cosmological, high resolution N-body + SPH simulations to follow the formation of disk galaxies with rotational velocities between 135 and 270 km/sec in a Lambda CDM universe. The simulations include gas cooling, star formation, the effects of a uniform UV background and a physically motivated description of feedback from supernovae. The host dark matter halos have a spin and last major merger redshift typical of galaxy sized halos as measured in recent large scale N--Body simulations. The simulated galaxies form rotationally supported disks with realistic exponential scale lengths and fall on both the I-band and baryonic Tully Fisher relations. An extended stellar disk forms inside the Milky Way sized halo immediately after the last major merger. The combination of UV background and SN feedback drastically reduces the number of visible satellites orbiting inside a Milky Way sized halo, bringing it in fair agreement with observations. Our simulations predict that the average age of a primary galaxy's stellar population decreases with mass, because feedback delays star formation in less massive galaxies. Galaxies have stellar masses and current star formation rates as a function of total mass that are in good agreement with observational data. We discuss how both high mass and force resolution and a realistic description of star formation and feedback are important ingredients to match the observed properties of galaxies.

ABRIDGED: We use HSTACS and NICMOS imaging to study the structure and colors of a sample of nine late-type spirals. We find: (1) A correlation between bulge and disks scale-lengths, and a correlation between the colors of the bulges and those of the inner disks. Our data show a trend for bulges to be more metal-enriched than their surrounding disks, but otherwise no simple age-metallicity connection between these systems; (2) A large range in bulge stellar population properties, and, in particular, in stellar ages. Specifically, in about a half of the late-type bulges in our sample the bulk of the stellar mass was produced recently. Thus, in a substantial fraction of the z=0 disk-dominated bulged galaxies, bulge formation occurs after the formation/accretion of the disk; (3) In about a half of the late-type bulges in our sample, however, the bulk of the stellar mass was produced at early epochs; (4) Even these "old" late-type bulges host a significant fraction of stellar mass in a young(er) c component; (5) A correlation for bulges between stellar age and stellar mass, in the sense that more massive late-type bulges are older than less massive late-type bulges. Since the overall galaxy luminosity (mass) also correlates with the bulge luminosity (mass), it appears that the galaxy mass regulates not only what fraction of itself ends up in the bulge component, but also "when" bulge formation takes place. We show that dynamical friction of massive clumps in gas-rich disks is a plausible disk-driven mode for the formation of "old" late-type bulges. If disk evolutionary processes are responsible for the formation of the entire family of late-type bulges, CDM simulations need to produce a similar number of initially bulgeless disks in addition to the disk galaxies that are observed to be bulgeless at z=0.

We have carried out a hydrodynamical code comparison study of interacting multiphase fluids. The two commonly used techniques of grid and smoothed particle hydrodynamics (SPH) show striking differences in their ability to model processes that are fundamentally important across many areas of astrophysics. Whilst Eulerian grid based methods are able to resolve and treat important dynamical instabilities, such as Kelvin-Helmholtz or Rayleigh-Taylor, these processes are poorly or not at all resolved by existing SPH techniques. We show that the reason for this is that SPH, at least in its standard implementation, introduces spurious pressure forces on particles in regions where there are steep density gradients. This results in a boundary gap of the size of the SPH smoothing kernel over which information is not transferred.

We study the formation of disks via the cooling flow of gas within galactic haloes using smoothed particle hydrodynamics simulations. These simulations resolve mass scales of a few thousand solar masses in the gas component for the first time. Thermal instabilities result in the formation of numerous warm clouds that are pressure confined by the hot ambient halo gas. The clouds fall slowly onto the disk through non-spherical accretion from material flowing preferentially down the angular momentum axis. The rotational velocity of the infalling cold gas decreases as a function of height above the disk, closely resembling that of the extra-planar gas recently observed around the spiral galaxy NGC 891.

We perform N-Body/SPH simulations of disk galaxy formation inside equilibrium spherical and triaxial cuspy dark matter halos. We systematically study the disk properties and morphology as we increase the numbers of dark matter and gas particles from 10^4 to 10^6 and change the force resolution. The force resolution influences the morphological evolution of the disk quite dramatically. Unless the baryon fraction is significantly lower than the universal value, with high force resolution a gaseous bar always forms within a billion years after allowing cooling to begin. The bar interacts with the disk, transferring angular momentum and increasing its scale length. In none of the simulations does the final mass distribution of the baryons obey a single exponential profile. Indeed within a few hundred parsecs to a kiloparsec from the center the density rises much more steeply than in the rest of the disk, and this is true irrespective of the presence of the bar.

We present a comprehensive series of $N$-body as well as $N$-body + SPH simulations to study the secular evolution of the structure of disk galaxies. Our simulations are organized in a hierarchy of increasing complexity, ranging from rigid-halo collisionless simulations to fully live simulations with gas and star formation. Comparisons between the different types of simulations allow us to isolate the role of various physical mechanisms. We focus on the evolution of systems expected in a LCDM universe. Our goal is to examine which structural properties of disk galaxies may result from secular evolution rather than from direct hierarchical assembly. In the vertical direction, we find that various mechanisms can lead to heating. The strongest heating occurs during the vertical buckling instability of a bar. Among the consequences of this instability is the formation of peanut-shaped bulges which produce clear kinematic signatures when observed face-on. We find that bars are robust structures that are not destroyed by buckling. They can be destroyed instead by a central mass concentration but we find that this mass needs to be a large fraction of the total mass of the disk. We then study the evolution of stellar surface density profiles showing how angular momentum redistribution leads to increasing central densities and disk scale lengths and to profile breaks at large radii. The breaks in these simulations are in excellent agreement with observed breaks, even when the evolution is purely collisionless. Disk scale-lengths increase even when the total disk angular momentum is conserved; thus mapping halo angular momenta to scale-lengths is non-trivial. [Abridged]

We report on the results of the first 3D SPH simulation of massive, gravitationally unstable protoplanetary disks with radiative transfer. We adopt a flux-limited diffusion scheme justified by the high opacity of most of the disk. The optically thin surface of the disk cools as a blackbody. The disks grow slowly in mass starting from a Toomre-stable initial condition to the point at which they become marginally unstable. We find that gravitationally bound clumps with masses close to a Jupiter mass can arise. Fragmentation appears to be driven by vertical convective-like motions capable of transporting the heat from the disk midplane to its surface on a timescale of only about 40 years at 10 AU. A larger or smaller cooling efficiency of the disk at the optically thin surface can promote or stifle fragmentation by affecting the vertical temperature profile, which determines whether convection can happen or not, and by regulating the accretion flow from optically thin regions towards overdense regions. We also find that the chances of fragmentation increase for a higher mean molecular weight $\mu$ since compressional heating is reduced. Around a star with mass $1 M_{\odot}$ only disks with $\mu \ge 2.4$, as expected for gas with a metallicity comparable to solar or higher, fragment. This suggests that disk instability, like core-accretion, should be more effective in forming gas giants at higher gas metallicities, consistent with the observed correlation between metallicity of the planet-hosting stars and frequency of planets

Protoplanetary gas disks are likely to experience gravitational instabilites (GI's) during some phase of their evolution. Density perturbations in an unstable disk grow on a dynamic time scale into spiral arms that produce efficient outward transfer of angular momentum and inward transfer of mass through gravitational torques. In a cool disk with rapid enough cooling, the spiral arms in an unstable disk form self-gravitating clumps. Whether gas giant protoplanets can form by such a disk instability process is the primary question addressed by this review. We discuss the wide range of calculations undertaken by ourselves and others using various numerical techniques, and we report preliminary results from a large multi-code collaboration. Additional topics include -- triggering mechanisms for GI's, disk heating and cooling, orbital survival of dense clumps, interactions of solids with GI-driven waves and shocks, and hybrid scenarios where GI's facilitate core accretion. The review ends with a discussion of how well disk instability and core accretion fare in meeting observational constraints.

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 $\Lambda$CDM model. Our findings have important implications for planned gravitational wave detection experiments such as {\it LISA}.

We used fully cosmological, high resolution N-body + SPH simulations to follow the formation of disk galaxies with a rotational velocity between 140 and 280 Km/sec in a Lambda CDM universe. The simulations include gas cooling, star formation, the effects of a uniform UV background and a physically motivated description of feedback from supernovae. The host dark matter halos have a spin and last major merger redshift typical of galaxy sized halos as measured in recent large scale N-Body simulations. Galaxies formed rotationally supported disks with realistic exponential scale lengths and fall on the I-band and baryonic Tully Fisher relations. The combination of UV background and SN feedback drastically reduced the number of visible satellites orbiting inside a Milky Way sized halo, bringing it fair agreement with observations. Feedback delays SF in small galaxies and more massive ones contain older stellar populations. The current star formation rates as a function of galaxy stellar mass are in good agreement with those measured by the SDSS.

Chandra images of the elliptical galaxy NGC4261 have revealed an anisotropy in the spatial distribution of the off-nuclear X-ray sources, interpreted by Zezas et al. as evidence of an association with a young stellar population. Our independent analysis of archive X-ray (Chandra) and optical (INT and HST) observations confirms the anisotropy of the X-ray sources but conducts to a different interpretation for their origin. We find that nearly half of the X-ray sources are associated to a globular cluster (optical counterpart) suggesting that they are accreting low-mass X-ray binaries. Where color index information is available, the X-ray sources are found to reside in red (metal-rich) systems. The luminosity function of the X-ray sources is also consistent with the one drawn from a population of Low Mass X-ray Binaries. We further investigate the properties of the sample of point-like sources obtained from archival optical images that we suggest are good globular cluster candidates and for which we find that the projected spatial distribution is non-homogeneous. In addition, we show that the distributions of the optical and X-ray populations are very similar, which leads us to conclude that the spatial anisotropy of the X-ray sources in NGC4261 is mostly a reflection of the anisotropy of the globular cluster population.

We discuss the role of environmental mechanisms in the evolution of dwarf galaxy satellites using high-resolution N-Body+SPH simulations that include simultaneously tidal forces, ram pressure and heating from ionizing radiation fields. Tidally induced bar-buckling instabilities can transform a rotating disky dwarf into a pressure supported spheroidal. Efficient gas removal requires instead a combination of tidal mass loss and ram pressure stripping in a diffuse gaseous corona around the primary system. The efficiency of ram pressure depends strongly on how extended the gas remains during the evolution. Bar driven inflows that tend to drive the gas to the bottom of the potential well can be opposed by the heating from ultraviolet radiation. We show that even fairly massive dwarfs ($V_{peak} >$ 30 km/s) would be stripped of their gas over a few Gyr if they enter the Milky Way halo at $z > 2$ thanks to the effect of the cosmic UV background. Gas mass loss can be much faster, occurring in less than 1 Gyr, if dwarf satellites have their first close approach with the primary at the epoch of bulge formation. Indeed at that time the primary galaxy should have a FUV luminosity comparable to that of major present-day starbursts, resulting in a local radiation field even more intense than the cosmic background.

We present a kinematic diagnostic for peanut-shaped bulges in nearly face-on galaxies. The face-on view provides a novel perspective on peanuts which would allow study of their relation to bars and disks in greater detail than hitherto possible. The diagnostic is based on the fact that peanut shapes are associated with a flat density distribution in the vertical direction. We show that the kinematic signature corresponding to such a distribution is a minimum in the fourth-order Gauss-Hermite moment $s_4$. We demonstrate our method on $N$-body simulations of varying peanut strength, showing that strong peanuts can be recognized to inclinations $i \simeq 30\degrees$, regardless of the strength of the bar. We also consider compound systems in which a bulge is present in the initial conditions as may happen if bulges form at high redshift through mergers. We show that in this case, because the vertical structure of the bulge is not derived from that of the disk, that the signature of a peanut in $s_4$ is weakened. Thus the same kinematic signature of peanuts can be used to explore bulge formation mechanisms. The observational requirements for identifying peanuts with this method are challenging, but feasible.

We perform high-resolution N-Body+SPH simulations of gas-rich dwarf galaxy satellites orbiting within a Milky Way-sized halo and study for the first time the combined effects of tides and ram pressure. The structure of the galaxy models and the orbital configurations are chosen in accordance to those expected in a LCDM Universe.While tidal stirring of disky dwarfs produces objects whose stellar structure and kinematics resembles that of dwarf spheroidals after a few orbits, ram pressure stripping is needed to entirely remove their gas component. Gravitational tides can aid ram pressure stripping by diminishing the overall potential of the dwarf, but tides also induce bar formation which funnels gas inwards making subsequent stripping more difficult. This inflow is particularly effective when the gas can cool radiatively. Assuming a low density of the hot Galactic corona consistent with observational constraints, dwarfs with V_{peak} < 30 km/s can be completely stripped of their gas content on orbits with pericenters of 50 kpc or less. Instead, dwarfs with more massive dark haloes and V_{peak} > 30 km/s lose most or all of their gas content only if a heating source keeps the gas extended, partially counteracting the bar-driven inflow. We show that the ionizing radiation from the cosmic UV background at z > 2 can provide the required heating. In these objects most of the gas is removed or becomes ionized at the first pericenter passage,explaining the early truncation of the star formation observed in Draco and Ursa Minor. The stripped gas breaks up into individual clouds pressure confined by the outer gaseous medium that have masses, sizes and densities comparable to the HI clouds recently discovered around M31.(abridged)

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.

We use high resolution N-Body/SPH simulations to study the hydrodynamical and gravitational interaction between the Large Magellanic Cloud and the Milky Way. We model the dark and hot extended halo components as well as the stellar/gaseous disks of the two galaxies. Both galaxies are embedded in extended cuspy LCDM dark matter halos. We follow the previous four Gyrs of the LMC's orbit such that it ends up with the correct location and orientation on the sky. Tidal forces elongate the LMC's disk, forcing a bar and creating a strong warp and diffuse stellar halo, although very few stars become unbound. The stellar halo may account for some of the microlensing events. Ram-pressure from a low density ionised halo is then sufficient to remove 1.4e8 Msolar of gas from the LMC's disk forming a great circle trailing stream around the Galaxy. The column density of stripped gas falls by two orders of magnitude 100 degrees from LMC. The LMC does not induce any response in the Milky Way disk. On the contrary, the tides raised by the Milky Way determine the truncation of the satellite at about 11 kpc. After several Gyrs of interaction the gas disk of the LMC is smaller than the stellar disk due to ram pressure and its size compares well with the observational data.

The recent discovery of hidden non-axysimmetric and disk-like structures in bright Virgo dwarf elliptical and lenticular galaxies (dE/dSph/dS0) indicates that they may have late-type progenitors. Using N-body simulations we follow the evolution of disk galaxies within a LambdaCDM cluster simulated with 10e7 particles, where the hierarchical growth and galaxy harassment are modeled self-consistently. Most of the galaxies undergo significant morphological transformation, even at the outskirts of the cluster, and move through the Hubble sequence from late type disks to dwarf spheroidals. None of the disks are completely destroyed therefore they can not be the progenitors of ultra compact dwarf galaxies (UCDs). The time evolution of the simulated galaxies is compared with unsharp-masked images obtained from VLT data and the projected kinematics of our models with the latest high resolution spectroscopic studies from the Keck and Palomar telescopes.

We discuss the role that numerical effects have in cosmological N-Body/SPH simulations of disk galaxy formation. We show that the resolution of current state-of-the-art calculations, about $10^5$ SPH and dark matter particles within the virial radius of a Milky Way-sized halo at $z=0$, is just enough to avoid artificial losses of angular momentum and dramatic numerical disk heating during the phase of disk assembly. Instead both effects will still be an issue at earlier epochs, when the progenitors of the final galaxies are resolved by significantly less than $10^5$ particles.This must be at least partially responsible for why simulations show a ubiquitous low angular momentum spheroid that is too massive compared to that of typical disk galaxies

The structural evolution of substructure in cold dark matter (CDM) models is investigated combining ``low-resolution'' satellites from cosmological N-body simulations of parent halos with N=10^7 particles with high-resolution individual subhalos orbiting within a static host potential. We show that, as a result of mass loss, convergence in the central density profiles requires the initial satellites to be resolved with N=10^7 particles and parsec-scale force resolution. We find that the density profiles of substructure halos can be well fitted with a power-law central slope that is unmodified by tidal forces even after the tidal stripping of over 99% of the initial mass and an exponential cutoff in the outer parts. The solution to the missing-satellites problem advocated by Stoehr et al. in 2002 relied on the flattening of the dark matter (DM) halo central density cusps by gravitational tides, enabling the observed satellites to be embedded within DM halos with maximum circular velocities as large as 60 km/s. In contrast, our results suggest that tidal interactions do not provide the mechanism for associating the dwarf spheroidal satellites (dSphs) of the Milky Way with the most massive substructure halos expected in a CDM universe. We compare the predicted velocity dispersion profiles of Fornax and Draco to observations, assuming that they are embedded in CDM halos. Models with isotropic and tangentially anisotropic velocity distributions for the stellar component fit the data only if the surrounding DM halos have maximum circular velocities in the range 20-35 km/s. If the dSphs are embedded within halos this large then the overabundance of satellites within the concordance LCDM cosmological model is significantly alleviated, but this still does not provide the entire solution.

We use high resolution 3D SPH simulations to study the evolution of self-gravitating binary protoplanetary disks. Heating by shocks and cooling are included. We consider different orbital separations and masses of the disks and central stars. Isolated massive disks ($M \sim 0.1 M_{{\o}dot}$) fragment into protoplanets as a result of gravitational instability for cooling times comparable to the orbital time. Fragmentation does not occur in binary systems with a separation of about 60 AU. This is because efficient heating owing to strong tidally induced spiral shocks damps any overdensity. The resulting temperatures, above 200 K, would vaporize water ice in the outer disk, posing a problem even for the other model of giant planet formation, core-accretion. Light disks ($M \sim 0.01 M_{\odot}$) do not fragment but remain cold because their low self-gravity inhibits strong shocks. Core accretion would not be hampered in the latter. At separations of about 120 AU the efficiency of fragmentation by disk instability rises and approaches that in isolated systems. If disk instability is the main formation mechanism for giant planets, on going surveys targeting binary systems should find considerably fewer planets in systems with separations below 100 AU.

We use high resolution collisionless $N$-body simulations to study the secular evolution of disk galaxies and in particular the final properties of disks that suffer a bar and perhaps a bar-buckling instability. Although we find that bars are not destroyed by the buckling instability, when we decompose the radial density profiles of the secularly-evolved disks into inner S\'ersic and outer exponential components, for favorable viewing angles, the resulting structural parameters, scaling relations and global kinematics of the bar components are in good agreement with those obtained for bulges of late-type galaxies. Round bulges may require a different formation channel or dissipational processes.

We simulate the formation of a realistic disk galaxy within the hierarchical scenario of structure formation and study its internal properties to the present epoch. We compare results from a LambdaCDM simulation with a LambdaWDM (2keV) simulation that forms significantly less small scale structure. We show how high mass and force resolution in both the gas and dark matter components play an important role in solving the angular momentum catastrophe claimed from previous simulations of galaxy formation within the hierarchical framework. The stellar material in the disk component has a final specific angular momentum equal to 40% and 90% of that of the dark halo in the LambdaCDM and LambdaWDM models respectively. The LambdaWDM galaxy has a drastically reduced satellite population and a negligible stellar spheroidal component. Encounters with satellites play only a minor role in disturbing the disk. Satellites possess a variety of star formation histories linked to mergers and pericentric passages along their orbit around the primary galaxy. In both cosmologies, the galactic halo retains most of the baryons accreted and builds up a hot gas phase with a substantial X-ray emission. Therefore, while we have been successful in creating a realistic stellar disk in a massive galaxy within the LambdaCDM scenario, energy injection emerges as necessary ingredient to reduce the baryon fraction in galactic halos, independent of the cosmology adopted. (abridged)

We present galaxy circular velocity functions (GCVFs) for 34 low redshift (z < 0.15) clusters identified in the Sloan Digital Sky Survey (SDSS), for fifteen clusters drawn from dark matter simulations of hierarchical structure growth in a LambdaCDM cosmology, and for ~29,000 SDSS field galaxies. We find that the observed and simulated cluster GCVFs take the form of a power law. The cumulative GCVFs of the simulated clusters are very similar across a wide range of cluster masses, provided individual subhalo circular velocities are scaled by the circular velocities of the parent cluster. Once all sources of error are accounted for, the intrinsic scatter in the cumulative, scaled observed cluster GCVF is consistent with the simulations. The slope of the observed cluster GCVF is -2.5, independent of cluster velocity dispersion. The average slope of the simulated GCVFs is somewhat steeper, although formally consistent given the errors. Using our highest resolution cluster, we find that the effects of baryons on galaxy rotation curves is to flatten the simulated cluster GCVF into better agreement with observations. Finally, we find that the field GCVF deviates significantly from a power law, being flatter than the cluster GCVF at circular velocities less than 200 km/s, and steeper at circular velocities greater than 200 km/s.

We present the results of high resolution SPH simulations of the evolution of gravitationally unstable protoplanetary disks. We report on calculations in which the disk is evolved using a locally isothermal or adiabatic equation of state (with shock heating), and also on new simulations in which cooling and heating by radiation are explicitly modeled. We find that disks with a minimum Toomre parameter $< 1.4$ fragment into several gravitationally bound protoplanets with masses from below to a few Jupiter masses. This is confirmed also in runs where the disk is given a quiet start, growing gradually in mass over several orbital times. A cooling time comparable to the orbital time is needed to achieve fragmentation, for disk masses in the range $0.08-0.1 M_{\odot}$. After about 30 orbital times, merging between the bound condensations always leads to 2-3 protoplanets on quite eccentric orbits.

We carry out a large set of very high resolution, three dimensional smoothed particle hydrodynamics (SPH) simulations describing the evolution of gravitationally unstable gaseous protoplanetary disks. We consider a broad range of initial disk parameters. Disk masses out to 20 AU range from 0.075 to 0.125 $M_{\odot}$, roughly consistent with the high-end of the mass distribution inferred for disks around T Tauri stars. The initial disks span minimum $Q$ parameters between 0.8 and 2, with most models being around $\sim 1.4$. The disks are evolved assuming either a locally isothermal equation of state or an adiabatic equation of state with varying $\gamma$. When overdensities above a specific threshold appear as a result of gravitational instability in a locally isothermal calculation, the equation of state is switched to adiabatic to account for the increased optical depth. We show that when a disk has a minimum $Q$ parameter less than 1.4 strong trailing spiral instabilities, typically three or four armed modes, form and grow until fragmentation occurs along the arms after about 5 mean disk orbital times. The resulting clumps contract quickly to densities several orders of magnitude higher than the initial disk density, and the densest of them survive even under adiabatic conditions. These clumps are stable to tidal disruption and merge quickly, leaving 2-3 protoplanets on fairly eccentric orbits (the mean eccentricity being around 0.2) with masses between 0.7 and more than $7 M_{Jup}$, well in agreement with those of detected extrasolar planets.after $\sim 10^3$ years. Fragmentation is not strongly dependent on whether the disk starts from a marginally unstable state or gradually achieves it, as shown by a test where the disk gradually achieves the critical $Q$ by growing in mass.

We show that a global relation between baryonic mass and virial velocity can be constructed from the scale of dwarf galaxies up to that of rich galaxy clusters. The slope of this relation is close to that expected if dark matter halos form in the standard hierarchical cosmogony and capture a universal baryon fraction, once the details of halo structure and the adiabatic contraction of halos due to cooling gas are taken into account. The scatter and deficiency of baryons within low mass halos ($V_{vir} < 50$ km/s) is consistent with the expected suppression of gas accretion by photo-evaporation due to the cosmic UV background at high redshift. The data are not consistent with significant gas removal from strong supernovae winds unless the velocities of galaxies measured from their gas kinematics are significantly lower than the true halo velocities for objects with $V_{vir} < 100$ km/s. Thus models such as $\Lambda$CDM with a steep mass function of halos may find it difficult to reproduce both the baryonic mass-velocity relation presented here whilst at the same time reproducing the flat luminosity/HI function of galaxies. Galaxies hold about 10% of the baryons in the Universe, which is close to the collapsed mass fraction expected within hierarchical models on these scales, suggesting a high efficiency for galaxy formation. Most of the baryons are expected to be evenly distributed between diffuse intergalactic gas in low density environments and the intra-galactic medium within galaxy groups.

We use high resolution N-Body/SPH simulations to study the hydro-dynamical and gravitational interaction between the Large Magellanic Cloud and the Milky Way. We model the dark and hot extended halo components as well as the stellar/gaseous disks of the two galaxies. Tidal forces distort the LMC's disk, forcing a bar and creating a diffuse stellar halo and a strong warp, although very few stars are unbound from the LMC. Ram-pressure from a low density ionised halo is then sufficient to remove $1.4 \times 10^8M_\odot$ of gas from the LMC's disk forming a great circle trailing stream around the Galaxy.

We present the first high-resolution N-Body/SPH simulations that follow the evolution of low surface brightness disk satellites in a primary halo containing both dark matter and a hot gas component. Tidal shocks turn the stellar disk into a spheroid with low $v/\sigma$ and remove most of the outer dark and baryonic mass. In addition, by weakening the potential well of the dwarf, tides enhance the effect of ram pressure, and the gas is stripped down to radius three times smaller than the stellar component A very low gas/stars ratio results after several Gyr, similarly to what seen in dwarf spheroidal satellites of the Milky Way and M31.

The Ultra Compact Dwarf (UCD) galaxies recently discovered in the Fornax and Virgo clusters exhibit structural similarity to the dense nuclei of nucleated dEs indicating that the progenitor galaxy and its halo have been entirely tidally disrupted. Using high resolution $N$-body simulations with up to ten million particles we investigate the evolution and tidal stripping of substructure halos orbiting within a host potential. We find that complete disruption of satellite halos modeled following the NFW density profile occurs only for very low values of concentration in disagreement with the theoretical predictions of CDM models. This discrepancy is further exacerbated when we include the effect of baryons since disk formation increases the central density.

We perform several high resolution N-Body/SPH simulations of low surface brightness galaxies (LSBs) embedded in cold dark matter halos to study how likely is bar formation in such systems. The behavior of various collisionless galaxy models is studied both in isolation and in the presence of a large perturbing satellite. We also consider models with a dominant gaseous component in the disk. We find that in general bar formation requires disk masses at least a factor of 2 higher than those inferred for LSBs under the assumption of a normal stellar mass-to-light ratio. Low surface density stellar disks contributing less than 10% of the total virial mass are stable within NFW halos spanning a range of concentrations. However, a purely gaseous disk can form a bar even for quite low masses and for realistic temperatures provided that cooling is very efficient (we adopt an isothermal equation of state) and that the halo has a very low concentration, c < 5. The bars that form in these LSB models are shorter than the typical halo scale radius - their overall angular momentum content might be too low to affect significantly the inner dark halo structure. Once formed, all the bars evolve into bulge-like structures in a few Gyr and can excite spiral patterns in the surrounding disk component. The recently discovered red LSBs show significant non-axisymmetric structure and bulge-like components and share many of their structural properties with the final states of our LSB models with massive disks. Our results imply that a bulge-like component must be present in any low surface brightness galaxy that ever went bar unstable in the past.

We study the evolution of cold, gravitationally unstable protoplanetary gaseous disks performing 3D SPH simulations with up to a million particles on large parallel machines. We show that self-gravitating protoplanets can form in disks with initial $Q_{min} \sim 1.4$ provided that disks are evolved with a locally isothermal equation of state until localized regions reach a sufficiently high density. Lower resolution simulations carried out for as long as a thousand years show that systems of giant protoplanets with masses and orbits comparable to those of the observed extrasolar planets arise naturally.

We study the evolution of dark matter satellites orbiting inside more massive haloes using semi-analytical tools coupled with high-resolution N-Body simulations. We select initial satellite sizes, masses, orbital energies, and eccentricities as predicted by hierarchical models of structure formation. Both the satellite and the main halo are described by a Navarro, Frenk & White density profile with various concentrations.

The evolution of gravitationally unstable protoplanetary gaseous disks has been studied with the use of three-dimensional smoothed particle hydrodynamics simulations with unprecedented resolution. We have considered disks with initial masses and temperature profiles consistent with those inferred for the protosolar nebula and for other protoplanetary disks. We show that long-lasting, self-gravitating protoplanets arise after a few disk orbital periods if cooling is efficient enough to maintain the temperature close to 50 K. The resulting bodies have masses and orbital eccentricities similar to those of detected extrasolar planets.

In a 6h Halpha exposure of the N-W region of the cluster of galaxies A1367 we discovered a 75 kpc cometary emission of ionized gas trailing behind two Irregular galaxies. The Halpha trails correspond in position and length with tails of syncrotron radiation. At the galaxy side opposite to the tails the two galaxies show bright HII regions aligned along arcs, where the star formation takes place at the prodigeous rate of about 1 solar mass per yr. From the morphology of the galaxies and of the trailing material, we infer that the two galaxies are suffering from ram pressure due to their high velocity motion through the cluster IGM. We estimate that about 10^9 solar masses of gas, probably ionized in the giant HII regions, is swept out forming the tails. The tails cross each other at some 100 kpc from the present galaxy location, indicating that a major tidal event occurred some about 5x10^7 yr ago. We exclude that mutual harassment produced the observed morphology and we show with numerical simulations that it could have marginally aided ram pressure stripping by loosening the potential well of the galaxies.

Recent observations suggest that Carina and other nearby dwarf spheroidal galaxies are surrounded by unbound stars tidally stripped by the Milky Way. We run high-resolution N-Body simulations of dwarf galaxies orbiting within the Milky Way halo to determine if such observations can be explained with dark matter potentials as those implied by current structure formation models. We show that tidal forces acting on dwarfs with constant density cores or with cuspy profiles having a low concentration parameter ($c < 5$) lead to flat outer stellar density profiles like that of Carina for a variety of orbital configurations. On the contrary, it is more difficult to remove stars from cuspy dark matter halos with concentrations as high as predicted by CDM models at the mass scale of dwarf galaxies ($c \simgt 10$) and the data can only be reproduced assuming nearly radial orbits. Our simulations show that Carina is losing mass at a fractional rate $< 0.1$ Gyr$^{-1}$ and its mass-to-light ratio could be inflated by at most a factor of 2 due to unbound stars projected along the line of sight. We follow the evolution of the tidal debris within a triaxial clumpy cold dark matter Milky Way halo which causes differential precession and small scale heating of the stellar streams. This renders their use as a dynamical tracer of the Galactic potential practically useless, but does provide a novel test of the nature of the dark matter. Models with warm dark matter (WDM) or fluid dark matter (FDM) produce dwarf halos with lower central densities than CDM and would be consistent with the observed tidal tails even for orbits with eccentricities as low as indicated by current data on nearby dwarf spheroidals. Galactic halos in FDM are smooth and spherical and would be favored by the detection of coherent streams.

We investigate how the survival of dark matter satellites inside virialized halos depends on tidal stripping and dynamical friction. We use an analytic approach and then compare the results with N-Body simulations. Both the satellites and the primary halos are similar to cosmological haloes and have NFW density profiles. Satellites can either merge with the primary halo or continue to move on barely perturbed orbits, eventually being disrupted, depending on the relative strength of friction and tidal forces. We propose that their actual fate depends simply on their mass ratio relative to the primary halo.

Recent N-Body simulations have shown that dIrrs evolve into dSphs owing to the strong tidal perturbation of the Milky Way. Satellites whose dark matter halos have a core or an NFW profile with a concentration $c < 5$ undergo severe stripping even on low eccentricity orbits and their remnants have projected stellar profiles that become flat at large radii as observed in some of the dSphs. Satellites with higher central densities, as predicted by CDM models, are more robust to tides and can reproduce the observations only on nearly radial orbits. Accurate proper motions for LG dSphs will allow to test current models of structure formation.

We present results from high-resolution N-Body/SPH simulations of rotationally supported dwarf irregular galaxies moving on bound orbits in the massive dark matter halo of the Milky Way.The dwarf models span a range in disk surface density and the masses and sizes of their dark halos are consistent with the predictions of cold dark matter cosmogonies. We show that the strong tidal field of the Milky Way determines severe mass loss in their halos and disks and induces bar and bending instabilities that transform low surface brightness dwarfs (LSBs) into dwarf spheroidals (dSphs) and high surface brightness dwarfs (HSBs) into dwarf ellipticals (dEs) in less than 10 Gyr. The final central velocity dispersions of the remnants are in the range 8-30 km/s and their final $v/\sigma$ falls to values $< 0.5$, matching well the kinematics of early-type dwarfs. The transformation requires the orbital time of the dwarf to be $\simlt 3-4$ Gyr, which implies a halo as massive and extended as predicted by hierarchical models of galaxy formation to explain the origin of even the farthest dSph satellites of the Milky Way, Leo I and Leo II. Only dwarfs with central dark matter densities as high as those of Draco and Ursa Minor can survive for 10 Gyr in the proximity of the Milky Way: this is naturally achieved within hierarchical models, where the densest objects should have small orbital times due to their early formation epochs. Part of the gas is stripped and part is funneled to the center due to the bar, generating one strong burst of star formation in HSBs and smaller, multiple bursts in LSBs. Extended low-surface brightness stellar and gaseous streams originate from LSBs and, when projected along the line of sight, can lead to overestimate the mass-to-light ratio of the bound remnant by a factor $\simlt 2$,

N-Body/SPH simulations are used to study the evolution of dwarf irregular galaxies (dIrrs) entering the dark matter halo of the Milky Way or M31 on plunging orbits. We propose a new dynamical mechanism driving the evolution of gas rich, rotationally supported dIrrs, mostly found at the outskirts of the Local Group (LG), into gas free, pressure supported dwarf spheroidals (dSphs) or dwarf ellipticals (dEs), observed to cluster around the two giant spirals. The initial model galaxies are exponential disks embedded in massive dark matter halos and reproduce nearby dIrrs. Repeated tidal shocks at the pericenter of their orbit partially strip their halo and disk and trigger dynamical instabilities that dramatically reshape their stellar component. After only 2-3 orbits low surface brightness (LSB) dIrrs are transformed into dSphs, while high surface brightness (HSB) dIrrs evolve into dEs. This evolutionary mechanism naturally leads to the morphology-density relation observed for LG dwarfs. Dwarfs surrounded by very dense dark matter halos, like the archetypical dIrr GR8, are turned into Draco or Ursa Minor, the faintest and most dark matter dominated among LG dSphs. If disks include a gaseous component, this is both tidally stripped and consumed in periodic bursts of star formation. The resulting star formation histories are in good qualitative agreement with those derived using HST color-magnitude diagrams for local dSphs.

The evolution of a small satellite inside a more massive truncated isothermal spherical halo is studied using both the Theory of Linear Response for dynamical friction and N-Body simulations. The analytical approach includes the effects of the gravitational wake, of the tidal deformation and the shift of the barycenter of the primary, so unifying the local versus global interpretation of dynamical friction. Sizes, masses, orbital energies and eccentricities are chosen as expected in hierarchical clustering models. We find that in general the drag force in self-gravitating backgrounds is weaker than in uniform media and that the orbital decay is not accompanied by a significant circularization. We also show that the dynamical friction time scale is weakly dependent on the initial circularity. We provide a fitting formula for the decay time that includes the effect of mass and angular momentum loss. Live satellites with dense cores can survive disruption up to an Hubble time within the primary, notwithstanding the initial choice of orbital parameters. Dwarf spheroidal satellites of the Milky Way, like Sagittarius A and Fornax, have already suffered mass stripping and, with their present masses, the sinking times exceed 10 Gyr even if they are on very eccentric orbits.

We study the evolution of galaxy satellites with high resolution N-body simulations. Satellites are modeled as replicas of typical low and high surface brightness galaxies (LSBs and HSBs). Encounters on high eccentricity orbits (as typical in hierarchical models of galaxy formation) strip LSBs of most of their stars and tend to decrease their surface brightness. On the contrary, bar instability in HSBs leads to substantial loss of angular momentum of the stellar component and to an increase of central surface brightness. In both cases the remnant resembles a spheroidal galaxy with an exponential surface brightness profile. A simple modeling of color evolution and interactions driven star formation gives M/L ratios for the remnants that are roughly consistent with observations. These results suggest an evolutionary scenario for the dwarf galaxies in our Local Group, faint dSphs being the descendents of LSBs and brighter dSphs/dEs being the final state of HSB satellites.