Source author record

Lucio Mayer

Lucio Mayer appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

95works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

95 published item(s)

preprint2026arXiv

Disc fragmentation. II. Ejection of low mass Free Floating Planets from growing binary systems

Observations indicate that disc fragmentation due to Gravitational Instability (GI) is the likely origin of massive companions to stars, such as giant planets orbiting M-dwarf stars, Brown Dwarf (BD) companions to FGK stars, and binary stars with separations smaller than 100 au. Additionally, we have recently showed that disc fragmentation in young rapidly evolving binary systems ejects an abundant population of massive Jupiter-mass Free-Floating Planets (FFPs). In this model, a massive disc around an initially single protostar undergoes GI and hatches a number of fragments; the most massive oligarch grows by runaway accretion into the secondary star. As the system rearranges itself from a single to a binary star configuration, a dramatic "pincer movement" by the binary ejects planets through dynamical interactions with the stars. Here we propose that the same scenario applies to an even more abundant population of smaller FFPs discovered by the microlensing surveys. Although disc fragmentation is usually believed to form only massive objects, several pathways for forming small core-dominated planets at separations of tens of au exist. We present results from three complementary simulation approaches, all of which confirm planet ejection efficiency as high as 0.5 for secondaries more massive than $\sim 10$\% of the primary star mass. On the other hand, Jovian mass planets migrate through the region of tens of au too rapidly to eject planets from that region. We discuss implications of this scenario, concluding that microlensing FFPs may be the most convincing evidence yet that disc fragmentation forms planets much less massive than Jupiter.

preprint2025arXiv

Semi-empirical constraints on the HI mass function of star-forming galaxies and $Ω_{\rm HI}$ at $z\sim 0.37$ from interferometric surveys

The HI mass function is a crucial tool to understand the evolution of the HI content in galaxies over cosmic times. We aim to derive semi-empirical constraints at $z\sim 0.37$ by combining literature results on the $M_\star$ function from optical surveys with recent findings on the $M_{\rm HI}-M_\star$ scaling relation derived via spectral stacking analysis applied to 21-cm line interferometric data from the MIGHTEE and CHILES surveys, conducted with the MeerKAT and VLA radio telescopes, respectively. We draw synthetic $M_\star$ samples directly from the publicly-available results underlying the analysis of the COSMOS2020 galaxy photometric sample. Afterwards, we convert $M_\star$ into $M_{\rm HI}$ using analytical fitting functions to the data points from HI stacking. We then fit a Schechter function to the median HIMF from all the samples via MCMC. We finally derive the posterior distribution for $Ω_{\rm HI}$ by integrating the models for the HIMF built from the posteriors samples of the Schechter parameters. We find a deviation of the HIMF at $z\sim 0.37$ from the results at $z\sim 0$ from the ALFALFA survey and at $z\sim 1$ from uGMRT data. Our results for $Ω_{\rm HI}$ are in broad agreement with other literature results, and follow the overall trend on $Ω_{\rm HI}$ as a function of redshift. The derived value $Ω_{\rm HI}=\left(7.02^{+0.59}_{-0.52}\right)\times10^{-4}$ at $z\sim 0.37$ from the combined analysis deviates at $\sim 2.9σ$ from the ALFALFA result at $z\sim 0$. Our findings about the HIMF and $Ω_{\rm HI}$ differ from previous literature results at $z\sim0$ and $z\sim1$, although we are unable to confirm at this stage whether these differences are due to cosmic evolution consistent with a smooth transition of the HI content of galaxies over the last 8 Gyr or due to selection biases and systematics.

preprint2023arXiv

Direct collapse of exceptionally heavy black holes in the merger-driven scenario

We revisit the conditions present in supermassive discs (SMDs) formed by the merger of gas-rich, metal-enriched galaxies at red-shift $z\sim 10$. We find that SMDs naturally form hydrostatic cores which go through a rapidly accreting supermassive star phase, before directly collapsing into massive black holes via the general relativistic instability. The growth and collapse of the cores occurs within $\sim 5\times 10^5$ yr from the formation of the SMD, producing bright electromagnetic, neutrino and gravitational wave transients with a typical duration of a few minutes and, respectively, a typical flux and a typical strain amplitude at Earth of $\sim 10^{-8}$ erg s$^{-1}$ cm$^{-2}$ and $\sim4\times 10^{-21}$. We provide a simple fitting formula for the the resulting black hole masses, which range from a few $10^6$ M$_{\odot}$ to $10^8$ M$_{\odot}$ depending on the initial SMD configuration. Crucially, our analysis does not require any specific assumption on the thermal properties of the gas, nor on the angular momentum loss mechanisms within the SMD. Led by these findings, we argue that the merger-driven scenario provides a robust pathway for the rapid formation of supermassive black holes at $z > 6$. It provides an explanation for the origin of the brightest and oldest quasars without the need of a sustained growth phase from a much smaller seed. Its smoking gun signatures can be tested directly via multi-messenger observations.

preprint2022arXiv

Dirty waveforms: multiband harmonic content of gas-embedded gravitational wave sources

We analyse the effect of stochastic torque fluctuations on the orbital evolution and the gravitational wave (GW)emission of gas-embedded sources with intermediate and extreme mass ratios. We show that gas-driven fluctuations imprint additional harmonic content in the GWs of the binary system, which we dub dirty waveforms(DWs). We find three interesting observational prospects for DWs, provided that torque fluctuations do indeed persist beyond the resolution limit of current hydrodynamical simulations. Firstly, DWs can produce a significant stochastic GW background, comparable to other GW noise sources. Secondly, the energy flux implied by the additional harmonics can cause a detectable secular phase shift in Laser Interferometer Space Antenna (LISA) sources, even if the net torque fluctuations vanish when averaged over orbital time-scales. Lastly, the DWs of moderate-redshift nHz supermassive binaries detectable by pulsar timing arrays (PTAs) could be detectable in the mHz range, producing a new type of PTA-LISA multiband gravitational source. Our results suggest that searching for DWs and their effects can potentially be a novel way to probe the heaviest of black holes and the physics of the accretion discs surrounding them. We find these results to be a further confirmation of the many exciting prospects of actively searching for environmental effects within the data stream of future GW detectors.

preprint2022arXiv

From EMBER to FIRE: predicting high resolution baryon fields from dark matter simulations with Deep Learning

Hydrodynamic simulations provide a powerful, but computationally expensive, approach to study the interplay of dark matter and baryons in cosmological structure formation. Here we introduce the EMulating Baryonic EnRichment (EMBER) Deep Learning framework to predict baryon fields based on dark-matter-only simulations thereby reducing computational cost. EMBER comprises two network architectures, U-Net and Wasserstein Generative Adversarial Networks (WGANs), to predict two-dimensional gas and HI densities from dark matter fields. We design the conditional WGANs as stochastic emulators, such that multiple target fields can be sampled from the same dark matter input. For training we combine cosmological volume and zoom-in hydrodynamical simulations from the Feedback in Realistic Environments (FIRE) project to represent a large range of scales. Our fiducial WGAN model reproduces the gas and HI power spectra within 10% accuracy down to ~10 kpc scales. Furthermore, we investigate the capability of EMBER to predict high resolution baryon fields from low resolution dark matter inputs through upsampling techniques. As a practical application, we use this methodology to emulate high-resolution HI maps for a dark matter simulation of a L=100 Mpc/h comoving cosmological box. The gas content of dark matter haloes and the HI column density distributions predicted by EMBER agree well with results of large volume cosmological simulations and abundance matching models. Our method provides a computationally efficient, stochastic emulator for augmenting dark matter only simulations with physically consistent maps of baryon fields.

preprint2022arXiv

Multiply lensed star forming clumps in the A521-sys1 galaxy at redshift 1

We study the population of star-forming clumps in A521-sys1, a $\rm z=1.04$ system gravitationally lensed by the foreground ($\rm z=0.25$) cluster Abell 0521. The galaxy presents one complete counter--image with a mean magnification of $\rm μ\sim4$ and a wide arc containing two partial images of A521-sys1 with magnifications reaching $\rm μ>20$, allowing the investigations of clumps down to scales of $\rm R_{eff}<50$ pc. We identify 18 unique clumps with a total of 45 multiple images. Intrinsic sizes and UV magnitudes reveal clumps with elevated surface brightnesses, comparable to similar systems at redshifts $\rm z\gtrsim1.0$. Such clumps account for $\sim40\%$ of the galaxy UV luminosity, implying that a significant fraction of the recent star-formation activity is taking place there. Clump masses range from $\rm 10^6\ M_\odot$ to $\rm 10^9\ M_\odot$ and sizes from tens to hundreds of parsec, resulting in mass surface densities from $10$ to $\rm 10^3\ M_\odot\ pc^{-2}$, with a median of $\rm \sim10^2\ M_\odot\ pc^{-2}$. These properties suggest that we detect star formation taking place across a wide range of scale, from cluster aggregates to giant star-forming complexes. We find ages of less than $100$ Myr, consistent with clumps being observed close to their natal region. The lack of galactocentric trends with mass, mass density, or age and the lack of old migrated clumps can be explained either by dissolution of clumps after few $\sim100$ Myr or by stellar evolution making them fall below the detectability limits of our data.

preprint2022arXiv

Presolar grain dynamics: creating nucleosynthetic variations through a combination of drag and viscous evolution

Meteoritic studies of solar system objects show evidence of nucleosynthetic heterogeneities that are inherited from small presolar grains (< 10 $μ$m) formed in stellar environments external to our own. The initial distribution and subsequent evolution of these grains are currently unconstrained. Using 3D, gas-dust simulations, we find that isotopic variations on the order of those observed in the solar system can be generated and maintained by drag and viscosity. Small grains are dragged radially outwards without size/density sorting by viscous expansion and backreaction, enriching the outer disc with presolar grains. Meanwhile large aggregates composed primarily of silicates drift radially inwards due to drag, further enriching the relative portion of presolar grains in the outer disc and diluting the inner disc. The late accumulation of enriched aggregates outside Jupiter could explain some of the isotopic variations observed in solar system bodies, such as the enrichment of supernovae derived material in carbonaceous chondrites. We also see evidence for isotopic variations in the inner disc that may hold implications for enstatite and ordinary chondrites that formed closer to the Sun. Initial heterogeneities in the presolar grain distribution that are not continuously reinforced are dispersed by diffusion, radial surface flows, and/or planetary interactions over the entire lifetime of the disc. For younger, more massive discs we expect turbulent diffusion to be even more homogenising, suggesting that dust evolution played a more central role in forming the isotopic anomalies in the solar system than originally thought.

preprint2022arXiv

The Dawn of Disk Formation in a Milky Way-sized Galaxy Halo: Thin Stellar Disks at $z > 4$

We present results from \textsc{GigaEris}, a cosmological, $N$-body hydrodynamical "zoom-in" simulation of the formation of a Milky Way-sized galaxy halo with unprecedented resolution, encompassing of order a billion particles within the refined region. The simulation employs a modern implementation of smoothed-particle hydrodynamics, including metal-line cooling and metal and thermal diffusion. We focus on the early assembly of the galaxy, down to redshift $z=4.4$. The simulated galaxy has properties consistent with extrapolations of the main sequence of star-forming galaxies to higher redshifts and levels off to a star formation rate of $\sim$60$\, M_{\odot}$~yr$^{-1}$ at $z=4.4$. A compact, thin rotating stellar disk with properties analogous to those of low-redshift systems arises already at $z \sim 8$. The galaxy rapidly develops a multi-component structure, and the disk, at least at these early stages, does not grow "upside-down" as often reported in the literature. Rather, at any given time, newly born stars contribute to sustain a thin disk. The kinematics reflect the early, ubiquitous presence of a thin disk, as a stellar disk component with $v_ϕ/σ_R$ larger than unity is already present at $z \sim 9$--10. Our results suggest that high-resolution spectro-photometric observations of very high-redshift galaxies should find thin rotating disks, consistent with the recent discovery of cold rotating gas disks by ALMA. Finally, we present synthetic images for the JWST NIRCam camera, showing how the early disk would be easily detectable already at those early times.

preprint2022arXiv

The role of bars on the dynamical-friction driven inspiral of massive perturbers

In this paper, we systematically explore the impact of a galactic bar on the inspiral time-scale of a massive object (MO) within a Milky Way-like galaxy. We integrate the orbit of MOs in a multi-component galaxy model via a semi-analytical approach that accounts for dynamical friction generalized to rotationally supported backgrounds. We compare the MO evolution in a galaxy featuring a Milky Way-like rotating bar to the evolution within an analogous axisymmetric galaxy without the bar. In agreement with previous studies, we find that the bar presence may significantly affect the inspiral, sometimes making it shorter by a factor of a few, sometimes hindering it for a Hubble time. The erratic behaviour is mainly impacted by the relative phase at which the MO encounters the stronger bar-induced resonances. In particular, the effect of the bar is more prominent for initially in-plane, prograde MOs, especially those crossing the bar co-rotation radius or outer Lindblad resonance. In the barred galaxy, we find the sinking of the most massive MOs (>~10^7.5 Msun) approaching the galaxy from large separations (>~8 kpc) to be most efficiently hampered. Neglecting the effect of global torques associated with the non-symmetric mass distribution is thus not advisable even within an idealized, smooth galaxy model; we further note that spiral patterns are unlikely to affect the inspiral due to their transient and fluctuating nature. We speculate that the sinking efficiency of massive black holes involved in minor galaxy mergers may be hampered in barred galaxies, making them less likely to host a gravitational wave signal accessible to low-frequency detectors.

preprint2021arXiv

A wide-orbit giant planet in the high-mass b Centauri binary system

Planet formation occurs around a wide range of stellar masses and stellar system architectures. An improved understanding of the formation process can be achieved by studying it across the full parameter space, particularly toward the extremes. Earlier studies of planets in close-in orbits around high-mass stars have revealed an increase in giant planet frequency with increasing stellar mass until a turnover point at 1.9 solar masses, above which the frequency rapidly decreases. This could potentially imply that planet formation is impeded around more massive stars, and that giant planets around stars exceeding 3 solar masses may be rare or non-existent. However, the methods used to detect planets in small orbits are insensitive to planets in wide orbits. Here we demonstrate the existence of a planet at 560 times the Sun-Earth distance from the 6-10 solar mass binary b Centauri through direct imaging. The planet-to-star mass ratio of 0.10-0.17% is similar to the Jupiter-Sun ratio, but the separation of the detected planet is ~100 times wider than that of Jupiter. Our results show that planets can reside in much more massive stellar systems than what would be expected from extrapolation of previous results. The planet is unlikely to have formed in-situ through the conventional core accretion mechanism, but might have formed elsewhere and arrived to its present location through dynamical interactions, or might have formed via gravitational instability.

preprint2021arXiv

BEAST begins: Sample characteristics and survey performance of the B-star Exoplanet Abundance Study

While the occurrence rate of wide giant planets appears to increase with stellar mass at least up through the A-type regime, B-type stars have not been systematically studied in large-scale surveys so far. It therefore remains unclear up to what stellar mass this occurrence trend continues. The B-star Exoplanet Abundance Study (BEAST) is a direct imaging survey with the extreme adaptive optics instrument SPHERE, targeting 85 B-type stars in the young Scorpius-Centaurus (Sco-Cen) region with the aim to detect giant planets at wide separations and constrain their occurrence rate and physical properties. The statistical outcome of the survey will help determine if and where an upper stellar mass limit for planet formation occurs. In this work, we describe the selection and characterization of the BEAST target sample. Particular emphasis is placed on the age of each system, which is a central parameter in interpreting direct imaging observations. We implement a novel scheme for age dating based on kinematic sub-structures within Sco-Cen, which complements and expands upon previous age determinations in the literature. We also present initial results from the first epoch observations, including the detections of ten stellar companions, of which six were previously unknown. All planetary candidates in the survey will need follow up in second epoch observations, which are part of the allocated observational programme and will be executed in the near future.

preprint2020arXiv

Concurrent effects of gas drag and planet migration on pebble accretion

We study the effect of a migrating planet ($10<M_p<20$ Earth mass) on the dynamics of pebbles in a radiative disk using 2D two-fluid simulations carried out with the RoSSBi code. The combined action of the waves induced by the migrating planet and drag back-reaction on the gas produces a complex evolution of the flux of pebbles. The waves excite zonal flows which act as dust traps, and accumulate pebbles. Owing to drag back-reaction, a Kelvin-Helmholtz instability develops, generating turbulent dust rings containing several Earth masses of solids, where planetesimal formation is likely. The more massive the planet, the faster the dust rings develop. In competition with planet migration, this process triggers several dust rings, with a radial separation scaling as $1/M_p$. We discover a transition around $13$ Earth masses. Below that, pebbles are stopped at the inner side of the planet's orbit, but pebble accretion on the planet is sustained. Above that, drag back-reaction in the dust ring region flattens the pressure profile while the planet migrates, which limits further growth by pebble accretion. This new {\it pebble isolation mass} is about a factor of 2 lower than often reported in the literature. This reduces the overall formation timescale of Super-Earth planets, while favoring their survival after the disk's dissipation because further accretion is stifled. Finally, our results support an hybrid model for the formation of Jupiter via both pebbles and planetesimals accretion.

preprint2020arXiv

Formation of the first stars and black holes

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.

preprint2020arXiv

Global simulations of self-gravitating magnetized protoplanetary disks

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.

preprint2020arXiv

Global torques and stochasticity as the drivers of massive black hole pairing in the young Universe

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.

preprint2020arXiv

Improved gravitational radiation time-scales: significance for LISA and LIGO-Virgo sources

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.

preprint2020arXiv

The erratic path to coalescence of LISA massive black hole binaries in sub-pc resolution simulations of smooth circumnuclear gas disks

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]

preprint2020arXiv

The Quest for Dual and Binary Supermassive Black Holes: A Multi-Messenger View

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.

preprint2019arXiv

Predicting dark matter halo formation in N-body simulations with deep regression networks

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.

preprint2019arXiv

Primordial Earth mantle heterogeneity caused by the Moon-forming giant impact

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.

preprint2018arXiv

Super-Eddington accretion; flow regimes and conditions in high-z galaxies

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

preprint2016arXiv

Bar-driven evolution and quenching of spiral galaxies in cosmological simulations

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]

preprint2016arXiv

Clumpy high-z galaxies as a testbed for feedback-regulated galaxy formation

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.

preprint2016arXiv

Diffuse coronae in cosmological simulations of Milky Way-sized galaxies

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 $β$-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.

preprint2016arXiv

Direct detection of precursors of gas giants formed by gravitational instability with the Atacama Large Millimetre/sub-millimetre Array

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.

preprint2016arXiv

Dust capture and long-lived density enhancements triggered by vortices in 2D protoplanetary disks

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.

preprint2016arXiv

Supermassive black hole pairs in clumpy galaxies at high redshift: delayed binary formation and concurrent mass growth

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.

preprint2016arXiv

Swift coalescence of supermassive black holes in cosmological mergers of massive galaxies

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 $Λ$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.

preprint2016arXiv

The ACS LCID Project XI. On the early time resolution of LG dwarf galaxy SFHs: Comparing the effects of reionization in models with observations

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.

preprint2016arXiv

The ISLANDS project I: Andromeda XVI, An Extremely Low Mass Galaxy not Quenched 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.

preprint2016arXiv

Tidal stirring of satellites with shallow density profiles prevents them from being too big to fail

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 Λ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 $γ=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.

preprint2015arXiv

A lower fragmentation mass scale in high redshift galaxies and its implications on giant clumps: a systematic numerical study

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.

preprint2015arXiv

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

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

preprint2015arXiv

On the gap-opening criterion of migrating planets in protoplanetary disks

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.

preprint2015arXiv

The evolution of dwarf galaxy satellites with different dark matter density profiles in the ErisMod simulations. I. The early infalls

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 $γ= 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/σ< 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.

preprint2015arXiv

The History of R-Process Enrichment in the Milky Way

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 [$α$/Fe] and [Eu/Fe].We argue that compact binary mergers could be the dominant source of $r$-process nucleosynthesis in the Galaxy.

preprint2014arXiv

Building Late-Type Spiral Galaxies by In-Situ and Ex-Situ Star Formation

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.

preprint2014arXiv

Comparing M31 and Milky Way Satellites: The Extended Star Formation Histories of Andromeda II and Andromeda XVI

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.

preprint2014arXiv

Direct formation of supermassive black holes in metal-enriched gas at the heart of high-redshift galaxy mergers

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.

preprint2014arXiv

Faint dwarfs as a test of DM models: WDM vs. CDM

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 $Λ$ 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.

preprint2014arXiv

Growth and activity of black holes in galaxy mergers with varying mass ratios

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.

preprint2014arXiv

Nuclear coups: dynamics of black holes in galaxy mergers

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.

preprint2014arXiv

Stellar and Gaseous Nuclear Disks Observed in Nearby (U)LIRGs

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/σ$ 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.

preprint2014arXiv

The ACS LCID project. X. The Star Formation History of IC 1613: Revisiting the Over-Cooling Problem

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.

preprint2014arXiv

The Argo Simulation: I. Quenching of Massive Galaxies at High Redshift as a Result of Cosmological Starvation

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.

preprint2014arXiv

The Argo Simulation: II. The Early Build-up of the Hubble Sequence

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.

preprint2014arXiv

The Baryon Cycle of Dwarf Galaxies: Dark, Bursty, Gas-Rich Polluters

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.

preprint2013arXiv

A "Light," Centrally-Concentrated Milky Way Halo?

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.

preprint2013arXiv

An Off-center Density Peak in the Milky Way's Dark Matter Halo?

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.

preprint2013arXiv

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

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

preprint2013arXiv

Early evolution of clumps formed via gravitational instability in protoplanetary disks; precursors of Hot Jupiters?

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.

preprint2013arXiv

Inside Out and Upside Down: Tracing the Assembly of a Simulated Disk Galaxy Using Mono-Age Stellar Populations

We analyze the present-day structure and assembly history of a high resolution hydrodynamic simulation of the formation of a Milky Way (MW)-like disk galaxy, from the "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.

preprint2013arXiv

Massive Black Hole binaries in gas-rich galaxy mergers; multiple regimes of orbital decay and interplay with gas inflows

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.

preprint2013arXiv

Massive black hole pairs in clumpy, self-gravitating circumnuclear disks: stochastic orbital decay

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.

preprint2013arXiv

The ACS LCID project. IX. Imprints of the early Universe in the radial variation of the star formation history of dwarf galaxies

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.

preprint2013arXiv

The circumgalactic medium of massive galaxies at z~3: a test for stellar feedback, galactic outflows, and cold streams

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)

preprint2013arXiv

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

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

preprint2012arXiv

Formation and Hardening of Supermassive Black Hole Binaries in Minor Mergers of Disk Galaxies

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

preprint2012arXiv

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

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

preprint2012arXiv

Hydrodynamics of galaxy mergers with supermassive black holes: is there a last parsec problem ?

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)

preprint2012arXiv

Massive black hole seeds born via direct gas collapse in galaxy mergers: their properties, statistics and environment

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.

preprint2012arXiv

Observability of Dual Active Galactic Nuclei in Merging Galaxies

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.

preprint2012arXiv

Pseudobulge Formation as a Dynamical Rather than a Secular Process

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.

preprint2012arXiv

The fate of high redshift massive compact galaxies in dense environments

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.

preprint2012arXiv

The origin of metals in the circum-galactic medium of massive galaxies at z=3

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)

preprint2012arXiv

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

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

preprint2011arXiv

Evolutionary tracks of tidally stirred disky dwarf galaxies

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

preprint2011arXiv

Formation and evolution of dwarf galaxies in the CDM Universe

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.

preprint2011arXiv

Formation of Dwarf Spheroidal Galaxies Via Mergers of Disky Dwarfs

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

preprint2011arXiv

Forming Realistic Late-Type Spirals in a LCDM Universe: The Eris Simulation

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)

preprint2011arXiv

The ACS LCID project. V. The Star Formation History of the Dwarf Galaxy \objectname[]{LGS-3}: Clues for Cosmic Reionization and Feedback

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.

preprint2011arXiv

The central slope of dark matter cores in dwarf galaxies: Simulations vs. THINGS

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 α= -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 α~ -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.

preprint2010arXiv

Biases in mass estimates of dSph galaxies

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

preprint2010arXiv

Clumps in the Outer Disk by Disk Instability: Why They are Initially Gas Giants and the Legacy of Disruption

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.

preprint2010arXiv

Mass Distribution in Galaxy Clusters: the Role of AGN Feedback

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.

preprint2010arXiv

On the Efficiency of the Tidal Stirring Mechanism for the Origin of Dwarf Spheroidals: Dependence on the Orbital and Structural Parameters of the Progenitor Disky Dwarfs

(Abridged) The tidal stirring model posits the formation of dSph galaxies via the tidal interactions between rotationally-supported dwarfs and MW-sized host galaxies. Using a set of collisionless N-body simulations, we investigate the efficiency of the tidal stirring mechanism. We explore a wide variety of dwarf orbital configurations and initial structures and demonstrate that in most cases the disky dwarfs experience significant mass loss and their stellar components undergo a dramatic morphological and dynamical transformation: from disks to bars and finally to pressure-supported spheroidal systems with kinematic and structural properties akin to those of the classic dSphs in the Local Group (LG). Our results suggest that such tidal transformations should be common occurrences within the currently favored cosmological paradigm and highlight the key factor responsible for an effective metamorphosis to be the strength of the tidal shocks at the pericenters of the orbit. We demonstrate that the combination of short orbital times and small pericenters, characteristic of dwarfs being accreted at high redshift, induces the strongest transformations. Our models also indicate that the transformation efficiency is affected significantly by the structure of the progenitor disky dwarfs. Lastly, we find that the dwarf remnants satisfy the relation Vmax = \sqrt{3} * sigma, where sigma is the 1D, central stellar velocity dispersion and Vmax is the maximum halo circular velocity, with intriguing implications for the missing satellites problem. Overall, we conclude that the action of tidal forces from the hosts constitutes a crucial evolutionary mechanism for shaping the nature of dwarf galaxies in environments such as that of the LG. Environmental processes of this type should thus be included as ingredients in models of dwarf galaxy formation and evolution.

preprint2010arXiv

Recoiling Massive Black Holes in Gas-Rich Galaxy Mergers

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.

preprint2010arXiv

The impact of dark matter cusps and cores on the satellite galaxy population around spiral galaxies

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

preprint2010arXiv

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

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

preprint2010arXiv

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

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

preprint2009arXiv

At the heart of the matter: the origin of bulgeless dwarf galaxies and Dark Matter cores

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.

preprint2009arXiv

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

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

preprint2009arXiv

Environmental mechanisms shaping the nature of dwarf spheroidal galaxies; the view of computer simulations

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)

preprint2009arXiv

Pairing of Supermassive Black Holes in unequal-mass galaxy mergers

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

preprint2009arXiv

The ACS LCID Project. I. Short-Period Variables in the Isolated Dwarf Spheroidal Galaxies Cetus & Tucana

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

preprint2009arXiv

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

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

preprint2009arXiv

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

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

preprint2009arXiv

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

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

preprint2008arXiv

The formation of disk galaxies in computer simulations

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

preprint2007arXiv

Early Gas Stripping as the Origin of the Darkest Galaxies in the Universe

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

preprint2006arXiv

Multi-scale simulations of merging galaxies with supermassive black holes

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

preprint2005arXiv

Fate of Supermassive Black Holes and Evolution of the M(BH)-sigma Relation in Merging Galaxies: The Effect of Gaseous Dissipation

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

preprint2003arXiv

On the Life and Death of Satellite Haloes

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.

preprint2001arXiv

The Metamorphosis of Tidally Stirred Dwarf Galaxies

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/σ$ 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$,

preprint2000arXiv

Tidal stirring and the origin of dwarf spheroidals in the Local Group

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.