Showing 1–42 of 42 results for author: Escala, A

Efficiency of black hole formation via collisions in stellar systems: An analysis of data from simulations and observations

Authors: M. C. Vergara, D. R. G. Schleicher, A. Escala, B. Reinoso, F. Flammini Dotti, A. W. H. Kamlah, M. Liempi, N. Hoyer, N. Neumayer, R. Spurzem

Abstract: This paper explores the theoretical relation between star clusters and black holes within, focusing on the potential role of nuclear star clusters (NSCs), globular clusters (GCs), and ultra compact dwarf galaxies (UCDs) as environments that allow for black hole formation via stellar collisions. This study aims to identify the optimal conditions for stellar collisions across a range of stellar syst… ▽ More This paper explores the theoretical relation between star clusters and black holes within, focusing on the potential role of nuclear star clusters (NSCs), globular clusters (GCs), and ultra compact dwarf galaxies (UCDs) as environments that allow for black hole formation via stellar collisions. This study aims to identify the optimal conditions for stellar collisions across a range of stellar systems leading to the formation of very massive stars that subsequently collapse into black holes. We analyze data from numerical simulations and observations of diverse stellar systems, encompassing various initial conditions, initial mass functions, and evolution scenarios. We computed a critical mass, determined by the interplay of collision time, system age, and initial properties of the star cluster. The efficiency of black hole formation ($ε_{\mathrm{BH}}$) is defined as the ratio of initial stellar mass divided by critical mass. We find that stellar systems with a ratio of initial stellar mass over critical mass above 1 exhibit a high efficiencies in terms of black hole formation, ranging from $30-100\%$. While there is some scatter, potentially attributed to complex system histories and the presence of gas, the results highlight the potential for achieving high efficiencies via a purely collisional channel in black hole formation. In conclusion, this theoretical exploration elucidates the connection between star clusters and black hole formation. The study underscores the significance of UCDs, GCs, and NSCs as environments conducive to the black hole formation scenario via stellar collisions. The defined black hole formation efficiency ($ε_{\mathrm{BH}}$) is shown to be influenced by the ratio of the initial stellar mass to the critical mass. △ Less

Submitted 11 June, 2024; v1 submitted 20 May, 2024; originally announced May 2024.

Comments: 16 pages, 5 figures

Origin of supermassive black holes: predictions for the black hole population

Authors: M. Liempi, L. Almonacid, D. R. G. Schleicher, A. Escala

Abstract: The presence of supermassive black holes at redshift z > 6 raises some questions about their formation and growth in the early universe. Due to the construction of new telescopes like the ELT to observe and detect SMBHs, it will be useful to derive theoretical estimates for the population and to compare observations and model predictions in the future. In consequence our main goal is to estimate t… ▽ More The presence of supermassive black holes at redshift z > 6 raises some questions about their formation and growth in the early universe. Due to the construction of new telescopes like the ELT to observe and detect SMBHs, it will be useful to derive theoretical estimates for the population and to compare observations and model predictions in the future. In consequence our main goal is to estimate the population of SMBHs using a semi-analytic code known as Galacticus which is a code for the formation and evolution of galaxies where we are about to include different scenarios for SMBHs formation indicating the initial mass of the black hole seed, its formation conditions and recipes for the evolution of the components of the galaxies. We found that the principal mechanism of growing SMBHs is is via galaxy mergers and accretion of matter. For the comparison of our results with observations, we calculate the radius of influence of the black hole to estimate which part of the population could be detected, leading to relations similar to the observed ones. △ Less

Submitted 29 May, 2023; originally announced May 2023.

Comments: 3 pages, 3 figures

Global instability by runaway collisions in nuclear stellar clusters: Numerical tests of a route for massive black hole formation

Authors: M. C. Vergara, A. Escala, D. R. G. Schleicher, B. Reinoso

Abstract: The centres of galaxies host nuclear stellar clusters, supermassive black holes, or both. The origin of this dichotomy is still a mystery. Nuclear stellar clusters are the densest stellar system in the Universe, so they are ideal places for runaway collisions to occur. Previous studies have proposed the possible existence of a critical mass scale in such clusters, for which the occurrence of colli… ▽ More The centres of galaxies host nuclear stellar clusters, supermassive black holes, or both. The origin of this dichotomy is still a mystery. Nuclear stellar clusters are the densest stellar system in the Universe, so they are ideal places for runaway collisions to occur. Previous studies have proposed the possible existence of a critical mass scale in such clusters, for which the occurrence of collisions becomes very frequent and leads to the formation of a very massive object. While it is difficult to directly probe this scenario with simulations, we here aim for a proof of concept using toy models where the occurrence of such a transition is shown based on simplified compact systems, where the typical evolution time-scales will be faster compared to the real Universe. Indeed our simulations confirm that such a transition takes place and that up to 50 per cent of the cluster mass can go into the formation of a central massive object for clusters that are above the critical mass scale. Our results thus support the proposed new scenario on the basis of idealized simulations. A preliminary analysis of observed nuclear star clusters shows similar trends related to the critical mass as in our simulations. We further discuss the caveats for the application of the proposed scenario in real nuclear star clusters. △ Less

Submitted 22 May, 2023; v1 submitted 29 September, 2022; originally announced September 2022.

Comments: 15 pages, 9 figures

Spiral Arms are Metal Freeways: Azimuthal Gas-Phase Metallicity Variations in Simulated Cosmological Zoom-in Flocculent Disks

Authors: Matthew E. Orr, Blakesley Burkhart, Andrew Wetzel, Philip F. Hopkins, Ivanna A. Escala, Allison L. Strom, Paul F. Goldsmith, Jorge L. Pineda, Christopher C. Hayward, Sarah R. Loebman

Abstract: We examine the azimuthal variations in gas-phase metallicity profiles in simulated Milky Way mass disk galaxies from the Feedback in Realistic Environments (FIRE-2) cosmological zoom-in simulation suite, which includes a sub-grid turbulent metal mixing model. We produce spatially resolved maps of the disks at $z \approx 0$ with pixel sizes ranging from 250 to 750~pc, analogous to modern integral f… ▽ More We examine the azimuthal variations in gas-phase metallicity profiles in simulated Milky Way mass disk galaxies from the Feedback in Realistic Environments (FIRE-2) cosmological zoom-in simulation suite, which includes a sub-grid turbulent metal mixing model. We produce spatially resolved maps of the disks at $z \approx 0$ with pixel sizes ranging from 250 to 750~pc, analogous to modern integral field unit (IFU) galaxy surveys, mapping the gas-phase metallicities in both the cold & dense gas and the ionized gas correlated with HII regions. We report that the spiral arms alternate in a pattern of metal rich and metal poor relative to the median metallicity on the order of $\lesssim 0.1$~dex, appearing generally in this sample of flocculent spirals. The pattern persists even in a simulation with different strengths of metal mixing, indicating that the pattern emerges from physics above the sub-grid scale. Local enrichment does not appear to be the dominant source of the azimuthal metallicity variations at $z \approx 0$: there is no correlation with local star formation on these spatial scales. Rather, the arms are moving inwards and outwards relative to each other, carrying their local metallicity gradients with them radially before mixing into the larger-scale interstellar medium. We propose that the arms act as freeways channeling relatively metal poor gas radially inwards, and relatively enriched gas radially outwards. △ Less

Submitted 28 September, 2022; originally announced September 2022.

Comments: 16 Pages, 14 Figures, submitted to MNRAS

Black hole fueling in galaxy mergers: A high-resolution analysis

Authors: Joaquín Prieto, Andrés Escala, George Privon, Juan d'Etigny

Abstract: Using parsec scale resolution hydrodynamical adaptive mesh refinement simulations we have studied the mass transport process throughout a galactic merger. The aim of such study is to connect both the peaks of mass accretion rate onto the BHs and star formation bursts with both gravitational and hydrodynamic torques acting on the galactic gaseous component. Our merger initial conditions were chosen… ▽ More Using parsec scale resolution hydrodynamical adaptive mesh refinement simulations we have studied the mass transport process throughout a galactic merger. The aim of such study is to connect both the peaks of mass accretion rate onto the BHs and star formation bursts with both gravitational and hydrodynamic torques acting on the galactic gaseous component. Our merger initial conditions were chosen to mimic a realistic system. The simulations include gas cooling, star formation, supernovae feedback, and AGN feedback. Gravitational and hydrodynamic torques near pericenter passes trigger gas funneling to the nuclei which is associated with bursts of star formation and black hole growth. Such episodes are intimately related with both kinds of torques acting on the galactic gas. Pericenters trigger both star formation and mass accretion rates of $\sim$ few $(1-10)$ $M_\odot$/yr. Such episodes last $\sim$ $(50-75)$ Myrs. Close passes also can produce black hole accretion that approaches and reaches the Eddington rate, lasting $\sim$ few Myrs. Our simulation shows that both gravitational and hydrodynamic torques are enhanced at pericenter passes with gravitational torques tending to have higher values than the hydrodynamic torques throughout the merger. We also find that in the closest encounters, hydrodynamic and gravitational torques can be comparable in their effect on the gas, the two helping in the redistribution of both angular momentum and mass in the galactic disc. Such phenomena allow inward mass transport onto the BH influence radius, fueling the compact object and lighting up the galactic nuclei. △ Less

Submitted 22 January, 2021; originally announced January 2021.

Comments: 13 pages, 9 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 508, Issue 3, December 2021, Pages 3672-3683

A kpc-scale resolved study of unobscured and obscured star-formation activity in normal galaxies at z = 1.5 and 2.2 from ALMA and HiZELS

Authors: Cheng Cheng, Edo Ibar, Ian Smail, Juan Molina, David Sobral, Andres Escala, Philip Best, Rachel Cochrane, Steven Gillman, Mark Swinbank, R. J. Ivison, Jia-Sheng Huang, Thomas M. Hughes, Eric Villard, Michele Cirasuolo

Abstract: We present Atacama Large Millimeter/Submillimeter Array (ALMA) continuum observations of a sample of nine star-forming galaxies at redshifts 1.47 and 2.23 selected from the High-$z$ Emission Line Survey (HiZELS). Four galaxies in our sample are detected at high significance by ALMA at a resolution of 0.25'' at rest-frame 355 $μ$m. Together with the previously observed H$α$ emission, from adaptive… ▽ More We present Atacama Large Millimeter/Submillimeter Array (ALMA) continuum observations of a sample of nine star-forming galaxies at redshifts 1.47 and 2.23 selected from the High-$z$ Emission Line Survey (HiZELS). Four galaxies in our sample are detected at high significance by ALMA at a resolution of 0.25'' at rest-frame 355 $μ$m. Together with the previously observed H$α$ emission, from adaptive optics-assisted integral-field-unit spectroscopy (0.15'' resolution), and F606W and F140W imaging from the Hubble Space Telescope (0.2'' resolution), we study the star-formation activity, stellar and dust mass in these high-redshift galaxies at $\sim$kpc-scale resolution. We find that ALMA detection rates are higher for more massive galaxies ($M_*>10^{10.5}$ M$_\odot$) and higher [N {\sc ii}]/H$α$ ratios ($>0.25$, a proxy for gas-phase metallicity). The dust extends out to a radius of 8 kpc, with a smooth structure, even for those galaxies presenting clumpy H$α$ morphologies. The half-light radii ($R_{\rm dust}$) derived for the detected galaxies are of the order $\sim$4.5 kpc, more than twice the size of submillimetre-selected galaxies at a similar redshift. Our global star-formation rate estimates -- from far-IR and extinction-corrected H$α$ luminosities -- are in good agreement. However, the different morphologies of the different phases of the interstellar medium suggest complex extinction properties of the high-redshift normal galaxies. △ Less

Submitted 1 November, 2020; originally announced November 2020.

Comments: 13 pages, 7 figures, MNRAS accepted

VALES VII: Molecular and ionized gas properties in pressure balanced interstellar medium of starburst galaxies at z ~ 0.15

Authors: Juan Molina, Edo Ibar, Nicolas Godoy, Andres Escala, Tomonari Michiyama, Cheng Cheng, Thomas M. Hughes, Maarten Baes, Yongquan Xue, Michał J. Michałowski, Paul van der Werf, Xue-Jian Jiang

Abstract: Context. Spatially resolved observations of the ionized and molecular gas are critical for understanding the physical processes that govern the interstellar medium (ISM) in galaxies. Aims. To study the morpho-kinematic properties of the ionized and molecular gas in three dusty starburst galaxies at $z = 0.12-0.17$ to explore the relation between molecular ISM gas phase dynamics and the star-format… ▽ More Context. Spatially resolved observations of the ionized and molecular gas are critical for understanding the physical processes that govern the interstellar medium (ISM) in galaxies. Aims. To study the morpho-kinematic properties of the ionized and molecular gas in three dusty starburst galaxies at $z = 0.12-0.17$ to explore the relation between molecular ISM gas phase dynamics and the star-formation activity. Methods. We analyse $\sim$kpc-scale ALMA CO(1--0) and seeing limited SINFONI Paschen-$α$ observations. We use a dynamical mass model, which accounts for beam-smearing effects, to constrain the CO-to-H$_2$ conversion factor. Results. One starburst galaxy shows irregular morphology which may indicate a major merger, while the other two systems show disc-like morpho-kinematics. The two disc-like starbursts show molecular gas velocity dispersion values comparable with that seen in local LIRG/ULIRGs, but in an ISM with molecular gas fraction and surface density values consistent to that reported for local star-forming galaxies. These molecular gas velocity dispersion values can be explained by assuming vertical pressure equilibrium. The star-formation activity is correlated with the molecular gas content suggesting depletion times of the order of $\sim 0.1-1$ Gyr. The star formation rate surface density ($Σ_{\rm SFR}$) correlates with the ISM pressure set by self-gravity ($P_{\rm grav}$) following a power law with an exponent close to 0.8. Conclusions. In dusty disc-like starburst galaxies, our data support the scenario in which the molecular gas velocity dispersion values are driven by the ISM pressure set by self-gravity, responsible to maintain the vertical pressure balance. The correlation between $Σ_{\rm SFR}$ and $P_{\rm grav}$ suggests that, in these dusty starbursts galaxies, the star formation activity arises as a consequence of the ISM pressure balance. △ Less

Submitted 17 September, 2020; originally announced September 2020.

Comments: Accepted for publication in A&A; 18 pages, 9 figures, 5 tables plus appendix. Figure 3 may not be correctly visualized in web browser

Journal ref: A&A 643, A78 (2020)

Elemental Abundances in the Local Group: Tracing the Formation History of the Great Andromeda Galaxy

Authors: Ivanna A. Escala

Abstract: The Local Group (LG) is an environment accessible to detailed studies of galaxy formation, providing a complement to the early universe. In particular, spectroscopy of resolved stellar populations in the LG provides kinematical and chemical information for individual stars that can be used to infer the history of $L_\star$ galaxies like the Milky Way (MW) and Andromeda (M31). The Gaia revolution i… ▽ More The Local Group (LG) is an environment accessible to detailed studies of galaxy formation, providing a complement to the early universe. In particular, spectroscopy of resolved stellar populations in the LG provides kinematical and chemical information for individual stars that can be used to infer the history of $L_\star$ galaxies like the Milky Way (MW) and Andromeda (M31). The Gaia revolution in the MW has enabled comprehensive observational studies of the MW's formation history. The exquisite detail in which the MW has been studied is currently not achievable in any other $L_\star$ galaxy. For this reason, the MW is a template for our understanding of galaxy formation. M31 is the only external galaxy that we can currently hope study in a level of detail approaching the MW. Studies of M31 have recently taken on greater significance, given the growing body of evidence that its formation history differs substantially from that of the MW. In an era of limited information about elemental abundances in M31, I have developed a technique to apply spectral synthesis to low-resolution stellar spectroscopy in order to measure abundances for individual giant stars in distant LG galaxies. Through undertaking the largest deep, spectroscopic survey of M31 to date with my collaborators, this has resulted in the first measurements of the elemental abundances in the inner stellar halo and stellar disk of M31, and the largest homogeneous catalog of elemental abundances in M31. With this foundational work, we have opened the doors to detailed studies of the chemical composition of M31. Now, we can begin to ask--and answer--what differences in the elemental abundances of the M31 and the MW imply for our knowledge of galaxy formation in the broader universe. △ Less

Submitted 4 August, 2020; originally announced August 2020.

Comments: Caltech PhD thesis (Jun 2020). See Table 1.1 for a summary of relevant M31 vs. MW halo and disk properties

Observational Support for Massive Black Hole Formation Driven by Runaway Stellar Collisions in Galactic Nuclei

Authors: Andres Escala

Abstract: We explore here an scenario for massive black hole formation driven by stellar collisions in galactic nuclei, proposing a new formation regime of global instability in nuclear stellar clusters triggered by runaway stellar collisions. Using order of magnitude estimations, we show that observed nuclear stellar clusters avoid the regime where stellar collisions are dynamically relevant over the whole… ▽ More We explore here an scenario for massive black hole formation driven by stellar collisions in galactic nuclei, proposing a new formation regime of global instability in nuclear stellar clusters triggered by runaway stellar collisions. Using order of magnitude estimations, we show that observed nuclear stellar clusters avoid the regime where stellar collisions are dynamically relevant over the whole system, while resolved detections of massive black holes are well into such collision-dominated regime. We interpret this result in terms of massive black holes and nuclear stellar clusters being different evolutionary paths of a common formation mechanism, unified under the standard terminology of being both central massive objects. We propose a formation scenario where central massive objects more massive than $\rm \sim 10^8 \, Msun$, which also have relaxation times longer that their collision times, will be too dense (in virial equilibrium) to be globally stable against stellar collisions and most of its mass will collapse towards the formation of a massive black hole. Contrarily, this will only be the case at the core of less dense central massive objects leading to the formation of black holes with much lower black hole efficiencies $\rm ε_{BH} = \frac{M_{BH}}{M_{CMO}}$, with these efficiencies $\rm ε_{BH}$ drastically growing for central massive objects more massive than $\rm \sim 10^7 \, Msun$, approaching unity around $\rm M_{CMO} \sim 10^8 \, Msun$. We show that the proposed scenario successfully explains the relative trends observed in the masses, efficiencies, and scaling relations between massive black holes and nuclear stellar clusters. △ Less

Submitted 9 September, 2020; v1 submitted 2 June, 2020; originally announced June 2020.

Comments: comments welcome (aescala@das.uchile.cl)

When Gas Dynamics Decouples from Galactic Rotation: Characterizing ISM Circulation in Disk Galaxies

Authors: Jose Utreras, Guillermo A. Blanc, Andres Escala, Sharon Meidt, Eric Emsellem, Frank Bigiel, Simon C. O. Glover, Jonathan Henshaw, Alex Hygate, J. M. Diederik Kruijssen, Erik Rosolowsky, Eva Schinnerer, Andreas Schruba

Abstract: In galactic disks, galactic rotation sets the bulk motion of gas, and its energy and momentum can be transferred toward small scales. Additionally, in the interstellar medium, random and noncircular motions arise from stellar feedback, cloud-cloud interactions, and instabilities, among other processes. Our aim is to comprehend to which extent small-scale gas dynamics is decoupled from galactic rot… ▽ More In galactic disks, galactic rotation sets the bulk motion of gas, and its energy and momentum can be transferred toward small scales. Additionally, in the interstellar medium, random and noncircular motions arise from stellar feedback, cloud-cloud interactions, and instabilities, among other processes. Our aim is to comprehend to which extent small-scale gas dynamics is decoupled from galactic rotation. We study the relative contributions of galactic rotation and local noncircular motions to the circulation of gas, $Γ$, a macroscopic measure of local rotation, defined as the line integral of the velocity field around a closed path. We measure the circulation distribution as a function of spatial scale in a set of simulated disk galaxies and we model the velocity field as the sum of galactic rotation and a Gaussian random field. The random field is parameterized by a broken power law in Fourier space, with a break at the scale $λ_c$. We define the spatial scale $λ_{\rm eq}$ at which galactic rotation and non-circular motions contribute equally to $Γ$. For our simulated galaxies, the gas dynamics at the scale of molecular clouds is usually dominated by noncircular motions, but in the center of galactic disks galactic rotation is still relevant. Our model shows that the transfer of rotation from large scales breaks at the scale $λ_c$ and this transition is necessary to reproduce the circulation distribution. We find that $λ_{\rm eq}$, and therefore the structure of the gas velocity field, is set by the local conditions of gravitational stability and stellar feedback. △ Less

Submitted 10 July, 2020; v1 submitted 2 March, 2020; originally announced March 2020.

Comments: The Astrophysical Journal, 20 pages, 16 figures

Journal ref: The Astrophysical Journal, Volume 892, 2020, Issue 2, id.94

The kiloparsec-scale gas kinematics in two star-forming galaxies at $z \sim 1.47$ seen with ALMA and VLT-SINFONI

Authors: J. Molina, Edo Ibar, I. Smail, A. M. Swinbank, E. Villard, A. Escala, D. Sobral, T. M. Hughes

Abstract: We present Atacama Large Millimeter/submillimeter Array (ALMA) CO($J=2-1$) observations of two main-sequence star-forming galaxies at $z\sim1.47$ taken from the High-Z Emission Line Survey (HiZELS). These two systems have been previously reported to be molecular gas rich $f_{\rm H_2} \equiv M_{\rm H_2}/(M_{\rm H_2} + M_\star) \sim 0.8$. We carried out a follow-up study to resolve, at $\sim$kpc-sca… ▽ More We present Atacama Large Millimeter/submillimeter Array (ALMA) CO($J=2-1$) observations of two main-sequence star-forming galaxies at $z\sim1.47$ taken from the High-Z Emission Line Survey (HiZELS). These two systems have been previously reported to be molecular gas rich $f_{\rm H_2} \equiv M_{\rm H_2}/(M_{\rm H_2} + M_\star) \sim 0.8$. We carried out a follow-up study to resolve, at $\sim$kpc-scales, the CO emission. These new observations are combined with our earlier ALMA observations (sensitive to diffuse CO emission) and compared with our previous H$α$-based study at matched spatial resolution. One target is marginally resolved in CO(2-1), showing complex dynamics with respect to the ionised gas traced by H$α$. While the other source is spatially resolved, enabling a detailed exploration of its internal dynamical properties. In this system, both gaseous phases show similar spatial extension, rotation velocities and velocity dispersions ($V_{\rm rot} \sim σ_v \sim 100$\,km\,s$^{-1}$) suggesting a rotational velocity to velocity dispersion ratio consistent with unity. By comparing the ionized and molecular gas tracers through the use of a two-dimensional kinematic model, we estimate a median depletion time $τ_{\rm dep}=2.3 \pm 1.2$\,Gyr for the galaxy as a whole. This value is in agreement with the average $τ_{\rm dep}$ value observed in local star-forming galaxies at similar spatial scales. Using a thick-disk dynamical modelling, we derive a dynamical mass $M_{\rm dyn} = (1.59\pm0.19) \times 10^{11}$\,$M_\odot$ within $\approx 6$\,kpc. This suggests a dark matter fraction ($f_{\rm DM} \equiv M_{\rm DM}/M_{\rm dyn}$) of $0.59\pm0.10$, in agreement with the average $f_{\rm DM}$ value derived from stacked rotation curve analysis of galaxies at similar redshift range. △ Less

Submitted 14 June, 2019; originally announced June 2019.

Comments: 15 Pages, 9 Figures and 4 Tables. Accepted for publication in MNRAS

The role of gas fragmentation during the formation of supermassive black holes

Authors: Matías Suazo, Joaquín Prieto, Andrés Escala, Dominik Schleicher

Abstract: We have performed cosmological hydrodynamic simulations to study the effect of fragmentation on the SMBH seed mass in the direct collapse formation scenario. We considered different background UV intensities, host halo spin, and halo merger histories. Our simulations in low-spin halos, in the presence of a strong UV background are consistent with the Direct Collapse Black Hole model, in which a si… ▽ More We have performed cosmological hydrodynamic simulations to study the effect of fragmentation on the SMBH seed mass in the direct collapse formation scenario. We considered different background UV intensities, host halo spin, and halo merger histories. Our simulations in low-spin halos, in the presence of a strong UV background are consistent with the Direct Collapse Black Hole model, in which a single massive object $\sim10^5$ M$_{\odot}$ is formed in the center of a proto-galaxy. While in our simulations under the presence of a low UV background, we find fragmentation and the formation of various minor seeds. These fragments have masses of 10$^3$ - 10$^4$ M$_{\odot}$. These values are significant if we consider the potential mergers between them, and the fact that these minor objects are formed earlier in cosmic time compared to the massive single seeds. In one of our simulations, we observed gas fragmentation even in the presence of a strong UV intensity. Said structure arose in a dark matter halo that formed after various merger episodes, and the one with the highest spin value. The final mass obtained was $\sim 10^5$ M$_{\odot}$ in this run. From these results, we conclude that fragmentation in fact produces less massive objects, however, they are still prone to merge. In simulations that formed many fragments, they all approach the most massive one with time. We see no uniqueness in the strength of the UV intensity value required to achieve a DCBH, since it depends in other factors like the system dynamics in our cases. △ Less

Submitted 22 October, 2019; v1 submitted 8 March, 2019; originally announced March 2019.

VALES V: A kinematic analysis of the molecular gas content in $H$-ATLAS galaxies at $z\sim0.03-0.35$ using ALMA

Authors: J. Molina, Edo Ibar, V. Villanueva, A. Escala, C. Cheng, M. Baes, H. Messias, C. Yang, F. E. Bauer, P. van der Werf, R. Leiton, M. Aravena, A. M. Swinbank, M. J. Michałowski, A. M. Muñoz-Arancibia, G. Orellana, T. M. Hughes, D. Farrah, G. De Zotti, M. A. Lara-López, S. Eales, L. Dunne

Abstract: We present Atacama Large Millimeter/submillimeter Array (ALMA) resolved observations of molecular gas in galaxies up to $z=0.35$ to characterise the role of global galactic dynamics on the global interstellar medium (ISM) properties. These observations consist of a sub-sample of 39 galaxies taken from the Valparaíso ALMA Line Emission Survey (VALES). From the CO($J=1-0)$ emission line, we quantify… ▽ More We present Atacama Large Millimeter/submillimeter Array (ALMA) resolved observations of molecular gas in galaxies up to $z=0.35$ to characterise the role of global galactic dynamics on the global interstellar medium (ISM) properties. These observations consist of a sub-sample of 39 galaxies taken from the Valparaíso ALMA Line Emission Survey (VALES). From the CO($J=1-0)$ emission line, we quantify the kinematic parameters by modelling the velocity fields. We find that the IR luminosity increases with the rotational to dispersion velocity ratio ($V_{\rm rot}/σ_v$, corrected for inclination). We find a dependence between $V_{\rm rot}/σ_v$ and the [CII]/IR ratio, suggesting that the so-called `[CII] deficit' is related to the dynamical state of the galaxies. We find that global pressure support is needed to reconcile the dynamical mass estimates with the stellar masses in our systems with low $V_{\rm rot}/σ_v$ values. The star formation rate (SFR) is weakly correlated with the molecular gas fraction ($f_{\rm H_2}$) in our sample, suggesting that the release of gravitational energy from cold gas may not be the main energy source of the turbulent motions seen in the VALES galaxies. By defining a proxy of the `star formation efficiency' parameter as the SFR divided by the CO luminosity (SFE$'\equiv$ SFR/L$'_{\rm CO}$), we find a constant SFE$'$ per crossing time ($t_{\rm cross}$). We suggest that $t_{\rm cross}$ may be the controlling timescale in which the star formation occurs in dusty $z\sim0.03-0.35$ galaxies. △ Less

Submitted 27 September, 2018; originally announced September 2018.

Comments: Accepted for publication in MNRAS. 28 pages, 15 figures, 2 Tables

The Principle of Similitude in Biology: From Allometry to the Formulation of Dimensionally Homogenous `Laws'

Authors: Andres Escala

Abstract: Meaningful laws of nature must be independent of the units employed to measure the variables. The principle of similitude (Rayleigh 1915) or dimensional homogeneity, states that only commensurable quantities (ones having the same dimension) may be compared, therefore, meaningful laws of nature must be homogeneous equations in their various units of measurement, a result which was formalized in the… ▽ More Meaningful laws of nature must be independent of the units employed to measure the variables. The principle of similitude (Rayleigh 1915) or dimensional homogeneity, states that only commensurable quantities (ones having the same dimension) may be compared, therefore, meaningful laws of nature must be homogeneous equations in their various units of measurement, a result which was formalized in the $\rm Π$ theorem (Vaschy 1892; Buckingham 1914). However, most relations in allometry do not satisfy this basic requirement, including the `3/4 Law' (Kleiber 1932) that relates the basal metabolic rate and body mass, which it is sometimes claimed to be the most fundamental biological rate (Brown et al. 2004) and the closest to a law in life sciences (West \& Brown 2004). Using the $\rm Π$ theorem, here we show that it is possible to construct a unique homogeneous equation for the metabolic rates, in agreement with data in the literature. We find that the variations in the dependence of the metabolic rates on body mass are secondary, coming from variations in the allometric dependence of the heart frequencies. This includes not only different classes of animals (mammals, birds, invertebrates) but also different exercise conditions (basal and maximal). Our results demonstrate that most of the differences found in the allometric exponents (White et al. 2007) are due to compare incommensurable quantities and that our dimensionally homogenous formula, unify these differences into a single formulation. We discuss the ecological implications of this new formulation in the context of the Malthusian's, Fenchel's and the total energy consumed in a lifespan relations. △ Less

Submitted 7 January, 2019; v1 submitted 7 July, 2017; originally announced July 2017.

Comments: A accepted for publication in Theoretical Ecology. Comments are welcome (andres.escala@aya.yale.edu)

How AGN and SNe feedback affect mass transport and black hole growth in high redshift galaxies

Authors: Joaquin Prieto, Andres Escala, Marta Volonteri, Yohan Dubois

Abstract: By using cosmological hydrodynamical simulations we study the effect of supernova (SN) and active galactic nuclei (AGN) feedback on the mass transport of gas on to galactic nuclei and the black hole (BH) growth down to redshift z~6. We study the BH growth in relation with the mass transport processes associated with gravity and pressure torques, and how they are modified by feedback. Cosmological… ▽ More By using cosmological hydrodynamical simulations we study the effect of supernova (SN) and active galactic nuclei (AGN) feedback on the mass transport of gas on to galactic nuclei and the black hole (BH) growth down to redshift z~6. We study the BH growth in relation with the mass transport processes associated with gravity and pressure torques, and how they are modified by feedback. Cosmological gas funelled through cold flows reaches the galactic outer region close to free-fall. Then torques associated to pressure triggered by gas turbulent motions produced in the circum-galactic medium by shocks and explosions from SNe are the main source of mass transport beyond the central ~ 100 pc. Due to high concentrations of mass in the central galactic region, gravitational torques tend to be more important at high redshift. The combined effect of almost free-falling material and both gravity and pressure torques produces a mass accretion rate of order ~ 1 M_sun/yr at ~ pc scales. In the absence of SN feedback, AGN feedback alone does not affect significantly either star formation or BH growth until the BH reaches a sufficiently high mass of $\sim 10^6$ M_sun to self-regulate. SN feedback alone, instead, decreases both stellar and BH growth. Finally, SN and AGN feedback in tandem efficiently quench the BH growth, while star formation remains at the levels set by SN feedback alone due to the small final BH mass, ~ few 10^5 M_sun. SNe create a more rarefied and hot environment where energy injection from the central AGN can accelerate the gas further. △ Less

Submitted 22 January, 2017; originally announced January 2017.

Comments: 19 pages, 11 figures

SINFONI-HiZELS: The dynamics, merger rates & metallicity gradients of `typical' star-forming galaxies at $z$ = 0.8-2.2

Authors: J. Molina, E. Ibar, A. M. Swinbank, D. Sobral, P. N. Best, I. Smail, A. Escala, M. Cirasuolo

Abstract: We present adaptive optics (AO) assisted SINFONI integral field unit (IFU) spectroscopy of eleven H$α$ emitting galaxies selected from the High-Z Emission Line Survey (HiZELS). We obtain spatially resolved dynamics on ~kpc-scales of star-forming galaxies (stellar mass M$_\star$ = 10$^{9.5-10.5}$ M$_\odot$ and star formation rate SFR = 2-30 M$_\odot$ yr$^{-1}$) near the peak of the cosmic star-form… ▽ More We present adaptive optics (AO) assisted SINFONI integral field unit (IFU) spectroscopy of eleven H$α$ emitting galaxies selected from the High-Z Emission Line Survey (HiZELS). We obtain spatially resolved dynamics on ~kpc-scales of star-forming galaxies (stellar mass M$_\star$ = 10$^{9.5-10.5}$ M$_\odot$ and star formation rate SFR = 2-30 M$_\odot$ yr$^{-1}$) near the peak of the cosmic star-formation rate history. Combining these observations with our previous SINFONI-HiZELS campaign, we construct a sample of twenty homogeneously selected galaxies with IFU AO-aided observations -- the `SHiZELS' survey, with roughly equal number of galaxies per redshift slice, at $z$ = 0.8, 1.47, and 2.23. We measure the dynamics and identify the major kinematic axis by modelling their velocity fields to extract rotational curves and infer their inclination-corrected rotational velocities. We explore the stellar mass Tully-Fisher relationship, finding that galaxies with higher velocity dispersions tend to deviate from this relation. Using kinemetry analyses we find that galaxy interactions might be the dominant mechanism controlling the star-formation activity at $z$ = 2.23 but they become gradually less important down to $z$ = 0.8. Metallicity gradients derived from the [NII]/H$α$ emission line ratio show a median negative gradient for the SHiZELS survey of $Δ$log(O/H)/$Δ$R = -0.026$\pm$0.008 dex kpc$^{-1}$. We find that metal-rich galaxies tend to show negative gradients, whereas metal-poor galaxies tend to exhibit positive metallicity gradients. This result suggests that the accretion of pristine gas in the periphery of galaxies plays an important role in replenishing the gas in `typical' star-forming galaxies. △ Less

Submitted 1 December, 2016; originally announced December 2016.

Comments: Accepted for publication in MNRAS. 16 pages, 10 figures, 2 tables

A Portrait of Cold Gas in Galaxies at 60pc Resolution and a Simple Method to Test Hypotheses That Link Small-Scale ISM Structure to Galaxy-Scale Processes

Authors: Adam K. Leroy, Annie Hughes, Andreas Schruba, Erik Rosolowsky, Guillermo Blanc, Alberto D. Bolatto, Dario Colombo, Andres Escala, Carsten Kramer, J. M. Diederik Kruijssen, Sharon Meidt, Jerome Pety, Miguel Querejeta, Karin Sandstrom, Eva Schinnerer, Kazimierz Sliwa, Antonio Usero

Abstract: The cloud-scale density, velocity dispersion, and gravitational boundedness of the interstellar medium (ISM) vary within and among galaxies. In turbulent models, these properties play key roles in the ability of gas to form stars. New high fidelity, high resolution surveys offer the prospect to measure these quantities across galaxies. We present a simple approach to make such measurements and to… ▽ More The cloud-scale density, velocity dispersion, and gravitational boundedness of the interstellar medium (ISM) vary within and among galaxies. In turbulent models, these properties play key roles in the ability of gas to form stars. New high fidelity, high resolution surveys offer the prospect to measure these quantities across galaxies. We present a simple approach to make such measurements and to test hypotheses that link small-scale gas structure to star formation and galactic environment. Our calculations capture the key physics of the Larson scaling relations, and we show good correspondence between our approach and a traditional "cloud properties" treatment. However, we argue that our method is preferable in many cases because of its simple, reproducible characterization of all emission. Using, low-J 12CO data from recent surveys, we characterize the molecular ISM at 60pc resolution in the Antennae, the Large Magellanic Cloud, M31, M33, M51, and M74. We report the distributions of surface density, velocity dispersion, and gravitational boundedness at 60pc scales and show galaxy-to-galaxy and intra-galaxy variations in each. The distribution of flux as a function of surface density appears roughly lognormal with a 1sigma width of ~0.3 dex, though the center of this distribution varies from galaxy to galaxy. The 60pc resolution line width and molecular gas surface density correlate well, which is a fundamental behavior expected for virialized or free-falling gas. Varying the measurement scale for the LMC and M31, we show that the molecular ISM has higher surface densities, lower line widths, and more self-gravity at smaller scales. △ Less

Submitted 22 June, 2016; originally announced June 2016.

Comments: Accepted for publication in the Astrophysical Journal, 36 pages (24+appendix), 21 figures (12+appendix), until publication high resolution version at http://www.astronomy.ohio-state.edu/~leroy.42/cloudscale.pdf

Unveiling the role of galactic rotation on star formation

Authors: José Utreras, Fernando Becerra, Andrés Escala

Abstract: We study the star formation process at galactic scales and the role of rotation through numerical simulations of spiral and starburst galaxies using the adaptive mesh refinement code Enzo. We focus on the study of three integrated star formation laws found in the literature: the Kennicutt-Schmidt (KS) and Silk-Elmegreen (SE) laws, and the dimensionally homogeneous equation proposed by Escala (2015… ▽ More We study the star formation process at galactic scales and the role of rotation through numerical simulations of spiral and starburst galaxies using the adaptive mesh refinement code Enzo. We focus on the study of three integrated star formation laws found in the literature: the Kennicutt-Schmidt (KS) and Silk-Elmegreen (SE) laws, and the dimensionally homogeneous equation proposed by Escala (2015) $Σ_{\rm SFR} \propto \sqrt{G/L}Σ_{\rm gas}^{1.5}$. We show that using the last we take into account the effects of the integration along the line of sight and find a unique regime of star formation for both types of galaxies, suppressing the observed bi-modality of the KS law. We find that the efficiencies displayed by our simulations are anti-correlated with the angular velocity of the disk $Ω$ for the three laws studied in this work. Finally, we show that the dimensionless efficiency of star formation is well represented by an exponentially decreasing function of $-1.9Ωt_{\rm ff}^{\rm ini}$, where $t_{\rm ff}^{\rm ini}$ is the initial free-fall time. This leads to a unique galactic star formation relation which reduces the scatter of the bi-modal KS, SE, and Escala (2015) relations by 43%, 43%, and 35% respectively. △ Less

Submitted 30 December, 2016; v1 submitted 10 June, 2016; originally announced June 2016.

Comments: Accepted in ApJ, 22 pages, 14 figures

MSC Class: 85A04

Journal ref: The Astrophysical Journal, Volume 833, Issue 1, article id. 13, 18 pp. (2016)

Cusps in the center of galaxies: a real conflict with observations or a numerical artefact of cosmological simulations?

Authors: A. N. Baushev, L. del Valle, L. E. Campusano, A. Escala, R. R. Muñoz, G. A. Palma

Abstract: Galaxy observations and N-body cosmological simulations produce conflicting dark matter halo density profiles for galaxy central regions. While simulations suggest a cuspy and universal density profile (UDP) of this region, the majority of observations favor variable profiles with a core in the center. In this paper, we investigate the convergency of standard N-body simulations, especially in the… ▽ More Galaxy observations and N-body cosmological simulations produce conflicting dark matter halo density profiles for galaxy central regions. While simulations suggest a cuspy and universal density profile (UDP) of this region, the majority of observations favor variable profiles with a core in the center. In this paper, we investigate the convergency of standard N-body simulations, especially in the cusp region, following the approach proposed by \protect{\citep{13}}. We simulate the well known Hernquist model using the SPH code Gadget-3 and consider the full array of dynamical parameters of the particles. We find that, although the cuspy profile is stable, all integrals of motion characterizing individual particles suffer strong unphysical variations along the whole halo, revealing an effective interaction between the test bodies. This result casts doubts on the reliability of the velocity distribution function obtained in the simulations. Moreover, we find unphysical Fokker-Planck streams of particles in the cusp region. The same streams should appear in cosmological N-body simulations, being strong enough to change the shape of the cusp or even to create it. Our analysis, based on the Hernquist model and the standard SPH code, strongly suggests that the UDPs generally found by the cosmological N-body simulations may be a consequence of numerical effects. A much better understanding of the N-body simulation convergency is necessary before a 'core-cusp problem' can properly be used to question the validity of the CDM model. △ Less

Submitted 10 May, 2017; v1 submitted 9 June, 2016; originally announced June 2016.

Comments: 7 pages, 5 figures, accepted to JCAP

Journal ref: JCAP (2017) vol. 05, p. 042

Multiscale mass transport in z~6 galactic discs: fueling black holes

Authors: Joaquin Prieto, Andres Escala

Abstract: By using AMR cosmological hydrodynamic N-body zoom-in simulations, with the RAMSES code, we studied the mass transport processes onto galactic nuclei from high redshift up to $z\sim6$. Due to the large dynamical range of the simulations we were able to study the mass accretion process on scales from $\sim50[kpc]$ to $\sim$ few $1[pc]$. We studied the BH growth on to the galactic center in relation… ▽ More By using AMR cosmological hydrodynamic N-body zoom-in simulations, with the RAMSES code, we studied the mass transport processes onto galactic nuclei from high redshift up to $z\sim6$. Due to the large dynamical range of the simulations we were able to study the mass accretion process on scales from $\sim50[kpc]$ to $\sim$ few $1[pc]$. We studied the BH growth on to the galactic center in relation with the mass transport processes associated to both the Reynolds stress and the gravitational stress on the disc. Such methodology allowed us to identify the main mass transport process as a function of the scales of the problem. We found that in simulations that include radiative cooling and SNe feedback, the SMBH grows at the Eddington limit for some periods of time presenting $\langle f_{EDD}\rangle\approx 0.5$ throughout its evolution. The $α$ parameter is dominated by the Reynolds term, $α_R$, with $α_R\gg 1$. The gravitational part of the $α$ parameter, $α_G$, has an increasing trend toward the galactic center at higher redshifts, with values $α_G\sim 1$ at radii <$\sim$ few $ 10^1[pc]$ contributing to the BH fueling. In terms of torques, we also found that gravity has an increasing contribution toward the galactic center at earlier epochs with a mixed contribution above $\sim 100 [pc]$. This complementary work between pressure gradients and gravitational potential gradients allows an efficient mass transport on the disc with average mass accretion rates of the order $\sim$ few $1 [M_{\odot}/yr]$. These level of SMBH accretion rates found in our cosmological simulations are needed in all models of SMBH growth that attempt to explain the formation of redshift $6-7$ quasars. △ Less

Submitted 25 May, 2016; v1 submitted 8 December, 2015; originally announced December 2015.

Super massive black holes in star forming gaseous circumnuclear discs

Authors: Luciano del Valle, Andres Escala, Cristian Maureira-Fredes, Juan Molina, Jorge Cuadra, Pau Amaro-Seoane

Abstract: Using N-body/SPH simulations we study the evolution of the separation of a pair of SMBHs embedded in a star forming circumnuclear disk (CND). This type of disk is expected to be formed in the central kilo parsec of the remnant of gas-rich galaxy mergers. Our simulations indicate that orbital decay of the SMBHs occurs more quickly when the mean density of the CND is higher, due to increased dynamic… ▽ More Using N-body/SPH simulations we study the evolution of the separation of a pair of SMBHs embedded in a star forming circumnuclear disk (CND). This type of disk is expected to be formed in the central kilo parsec of the remnant of gas-rich galaxy mergers. Our simulations indicate that orbital decay of the SMBHs occurs more quickly when the mean density of the CND is higher, due to increased dynamical friction. However, in simulations where the CND is fragmented in high density gaseous clumps (clumpy CND), the orbits of the SMBHs are erratically perturbed by the gravitational interaction with these clumps, delaying, in some cases, the orbital decay of the SMBHs. The densities of these gaseous clumps in our simulations and in recent studies of clumpy CNDs are significantly higher than the observed density of molecular clouds in isolated galaxies or ULIRGs, thus, we expect that SMBH orbits are perturbed less in real CNDs than in the simulated CNDs of this study and other recent studies. We also find that the migration timescale has a weak dependence on the star formation rate of the CND. Furthermore, the migration timescale of a SMBH pair in a star-forming clumpy CND is at most a factor three longer than the migration timescale of a pair of SMBHs in a CND modeled with more simple gas physics. Therefore, we estimate that the migration timescale of the SMBHs in a clumpy CND is on the order of $10^7$ yrs. △ Less

Submitted 21 September, 2015; v1 submitted 5 March, 2015; originally announced March 2015.

Comments: 12 pages, 17 figures

Journal ref: 2015 Astrophys. J. 811 59

On the Functional Form of the Universal Star Formation Law

Authors: Andres Escala

Abstract: We study the functional form of the star formation law, using the Vaschy-Buckingham Pi theorem. We find that that it should have a form $\rm \dotΣ_{\star} \propto \sqrt{\frac{G}{L}}Σ_{gas}^{3/2}$, where L is a characteristic length that is related with an integration scale. With a reasonable estimation for L, we find that galaxies from different types and redshifts, including Low Surface Brightnes… ▽ More We study the functional form of the star formation law, using the Vaschy-Buckingham Pi theorem. We find that that it should have a form $\rm \dotΣ_{\star} \propto \sqrt{\frac{G}{L}}Σ_{gas}^{3/2}$, where L is a characteristic length that is related with an integration scale. With a reasonable estimation for L, we find that galaxies from different types and redshifts, including Low Surface Brightness galaxies, and individual star-forming regions in our galaxy, obey this single star formation law. We also find that depending on the assumption for L, this star formation law adopt different formulations of $\rm \dotΣ_{\star}$ scaling, that are widely studied in the literature: $\rm Σ_{gas}^{3/2}, Σ_{gas}/t_{orb}, Σ_{gas}/t_{ff} \, and \, Σ_{gas}^{2}/v_{turb}$. We also study secondary controlling parameters of the star formation law, based on the current evidence from numerical simulations and find that for galaxies, the star formation efficiency should be controlled, at least, by the turbulent Toomre parameter, the sonic and Alfvenic Mach numbers. △ Less

Submitted 19 February, 2015; v1 submitted 25 November, 2014; originally announced November 2014.

Comments: Accepted version after referee report. 20 pages, 4 figures

Super Star Clusters in Luminous Infrared Galaxies: the SUNBIRD Survey

Authors: P. Vaisanen, Z. Randriamanakoto, A. Escala, E. Kankare, A. Kniazev, J. K. Kotilainen, S. Mattila, R. Ramphul, S. Ryder, A. Tekola

Abstract: We present recent results from an adaptive optics imaging survey of 40 Luminous IR Galaxies (LIRGs) searching for obscured core collapse supernovae and studying the galaxies themselves. Here, in particular, we discuss the Super Star Clusters (SSC) populations in the LIRGs. We have constructed the first statistically significant samples of Luminosity Functions (LF) of SSCs in the near-IR, and find… ▽ More We present recent results from an adaptive optics imaging survey of 40 Luminous IR Galaxies (LIRGs) searching for obscured core collapse supernovae and studying the galaxies themselves. Here, in particular, we discuss the Super Star Clusters (SSC) populations in the LIRGs. We have constructed the first statistically significant samples of Luminosity Functions (LF) of SSCs in the near-IR, and find evidence that the LF slopes in LIRGs are shallower than in more quiescent spiral galaxies. Distance and blending effects were investigated in detail paving the way for SSC studies further out than done previously. We have also correlated the luminosities of the brightest clusters with the star formation rates (SFR) of the hosts. The relation is similar, though somewhat steeper than that found in the optical and at lower SFR levels, suggesting systematic extinction and/or age effects. We find that the characteristics of the relation suggest an underlying physical driver rather than solely a size-of-sample effect. In particular, a truncated luminosity/mass function would naturally explain the small scatter we find. Finally, we are modelling the ages and masses of our near-IR detected clusters in conjunction with HST optical data and present early results of using SSC properties to trace the histories of the target LIRG systems. △ Less

Submitted 10 February, 2014; originally announced February 2014.

Comments: Slightly expanded version of proceedings to appear in "Massive Young Star Clusters Near and Far: From the Milky Way to Reionization", 2013 Guillermo Haro Conference. Eds. Y. D. Mayya, D. Rosa-Gonzalez, & E. Terlevich, INAOE and AMC

The Interstellar Medium and Star Formation in Local Galaxies: Variations of the Star Formation Law in Simulations

Authors: Fernando Becerra, Andres Escala

Abstract: We use the Adaptive Mesh Refinement code Enzo to model the interstellar medium in isolated local disk galaxies. The simulation includes a treatment for star formation and stellar feedback. We get a highly supersonic turbulent disk, which is fragmented at multiple scales and characterized by a multi-phase interstellar medium. We show that a Kennicutt-Schmidt (KS) relation only holds when averaging… ▽ More We use the Adaptive Mesh Refinement code Enzo to model the interstellar medium in isolated local disk galaxies. The simulation includes a treatment for star formation and stellar feedback. We get a highly supersonic turbulent disk, which is fragmented at multiple scales and characterized by a multi-phase interstellar medium. We show that a Kennicutt-Schmidt (KS) relation only holds when averaging over large scales. However, values of star formation rates and gas surface densities lie close in the plot for any averaging size. This suggests an intrinsic relation between stars and gas at cell-size scales, which dominates over the global dynamical evolution. To investigate this effect, we develop a method to simulate the creation of stars based on the density field from the snapshots, without running the code again. We also investigate how the star formation law is affected by the characteristic star formation timescale, the density threshold and the efficiency considered in the recipe. We find that the slope of the law might vary from ~1.4 for a free-fall timescale, to ~1.0 for a constant depletion timescale. We further demonstrate that a power-law is recovered just by assuming that the mass of the new stars is a fraction of the mass of the cell $m_\star=ερ_{\rm gas}Δx^3$, with no other physical criteria required. We show that both efficiency and density threshold do not affect the slope, but the right combination of them can adjust the normalization of the relation, which in turn could explain a possible bi-modality in the law. △ Less

Submitted 23 April, 2014; v1 submitted 20 January, 2014; originally announced January 2014.

Comments: 30 pages, 12 figures. Replaced to match published version in ApJ

Journal ref: Astrophys.J.786:56-68,2014

Binary Disk interaction II: Gap-Opening criteria for unequal mass binaries

Authors: Luciano del Valle, Andres Escala

Abstract: We study the interaction between an unequal mass binary with an isothermal circumbinary disk, motivated by the theoretical and observational evidence that after a major merger of gas-rich galaxies, a massive gaseous disk with a SMBH binary will be formed in the nuclear region. We focus on the gravitational torques that the binary exerts onto the disk and how these torques can drive the formation o… ▽ More We study the interaction between an unequal mass binary with an isothermal circumbinary disk, motivated by the theoretical and observational evidence that after a major merger of gas-rich galaxies, a massive gaseous disk with a SMBH binary will be formed in the nuclear region. We focus on the gravitational torques that the binary exerts onto the disk and how these torques can drive the formation of a gap in the disk. This exchange of angular momentum between the binary and the disk is mainly driven by the gravitational interaction between the binary and a strong non-axisymmetric density perturbation that is produced in the disk, as response to the presence of the binary. Using SPH numerical simulations we tested two gap-opening criterion, one that assumes that the geometry of the density perturbation is an ellipsoid/thick-spirals and another that assumes a geometry of flat-spirals for the density perturbation. We find that the flat-spirals gap opening criterion successfully predicts which simulations will have a gap on the disk and which simulations will not have a gap on the disk. We also study the limiting cases predicted by the gap-opening criteria. Since the viscosity in our simulations is considerably smaller than the expected value in the nuclear regions of gas-rich merging galaxies, we conclude that in such environments the formation of a circumbinary gap is unlikely. △ Less

Submitted 16 October, 2013; originally announced October 2013.

Near-infrared adaptive optics imaging of infrared luminous galaxies: the brightest cluster magnitude - star formation rate relation

Authors: Zara Randriamanakoto, Andres Escala, Petri Vaisanen, Erkki Kankare, Jari Kotilainen, Seppo Mattila, Stuart Ryder

Abstract: We have established a relation between the brightest super star cluster magnitude in a galaxy and the host star formation rate (SFR) for the first time in the near infrared (NIR). The data come from a statistical sample of ~ 40 luminous IR galaxies (LIRGs) and starbursts utilizing K-band adaptive optics imaging. While expanding the observed relation to longer wavelengths, less affected by extincti… ▽ More We have established a relation between the brightest super star cluster magnitude in a galaxy and the host star formation rate (SFR) for the first time in the near infrared (NIR). The data come from a statistical sample of ~ 40 luminous IR galaxies (LIRGs) and starbursts utilizing K-band adaptive optics imaging. While expanding the observed relation to longer wavelengths, less affected by extinction effects, it also pushes to higher SFRs. The relation we find, M_K ~ - 2.6 log SFR, is similar to that derived previously in the optical and at lower SFRs. It does not, however, fit the optical relation with a single optical to NIR color conversion, suggesting systematic extinction and/or age effects. While the relation is broadly consistent with a size-of-sample explanation, we argue physical reasons for the relation are likely as well. In particular, the scatter in the relation is smaller than expected from pure random sampling strongly suggesting physical constraints. We also derive a quantifiable relation tying together cluster-internal effects and host SFR properties to possibly explain the observed brightest SSC magnitude vs. SFR dependency. △ Less

Submitted 28 August, 2013; originally announced August 2013.

Comments: 6 pages, 4 figures, accepted for publication in ApJ Letters

The Impact of Bound Stellar Orbits and General Relativity on the Temporal Behavior of Tidal Disruption Flares

Authors: Lixin Dai, Andres Escala, Paolo Coppi

Abstract: We have carried out general relativistic particle simulations of stars tidally disrupted by massive black holes. When a star is disrupted in a bound orbit with moderate eccentricity instead of a parabolic orbit, the temporal behavior of the resulting stellar debris changes qualitatively. The debris is initially all bound, returning to pericenter in a short time ~ the original stellar orbital times… ▽ More We have carried out general relativistic particle simulations of stars tidally disrupted by massive black holes. When a star is disrupted in a bound orbit with moderate eccentricity instead of a parabolic orbit, the temporal behavior of the resulting stellar debris changes qualitatively. The debris is initially all bound, returning to pericenter in a short time ~ the original stellar orbital timescale. The resulting fallback rate can thus be much higher than the Eddington rate. Furthermore if the star is disrupted close to the hole, in a regime where general relativity is important, the stellar and debris orbits display general relativistic precession. Apsidal precession can make the debris stream cross itself after several orbits, likely leading to fast debris energy dissipation. If the star is disrupted in an inclined orbit around a spinning hole, nodal precession reduces the probability of self-intersection, and circularization may take many dynamical timescales, delaying the onset of flare activity. An examination of the particle dynamics suggests that quasi-periodic flares with short durations, produced when the center of the tidal stream passes pericenter, may occur in the early-time light curve. The late-time light curve may still show power-law behavior which is generic to disk accretion processes. The detection triggers for future surveys should be extended to capture such "non-standard" short-term flaring activity before the event enters the asymptotic decay phase, as this activity is likely to be more sensitive to physical parameters such as the black hole spin. △ Less

Submitted 18 September, 2013; v1 submitted 20 March, 2013; originally announced March 2013.

Comments: 19 pages, 5 figures, Accepted for publication in ApJL

Binary-disk interaction: Gap-Opening criteria

Authors: Luciano del Valle, Andres Escala

Abstract: We study the interaction of an equal mass binary with an isothermal circumbinary disk motivated by the theoretical and observational evidence of the formation of massive black holes binaries surrounded by gas, after a major merger of gas-rich galaxies. We focus on the torques that the binary produces on the disk and how the exchange of angular momentum can drive the formation of a gap on it. We pr… ▽ More We study the interaction of an equal mass binary with an isothermal circumbinary disk motivated by the theoretical and observational evidence of the formation of massive black holes binaries surrounded by gas, after a major merger of gas-rich galaxies. We focus on the torques that the binary produces on the disk and how the exchange of angular momentum can drive the formation of a gap on it. We propose that the angular momentum exchange between the binary and the disk is through the gravitational interaction of the binary and a (tidally formed) global non-axisymmetric perturbation in the disk. Using this gravitational interaction we derive an analytic criterion for the formation of a gap in the disk that can be expressed on the structural parameters h/a and M(< r)/M_{bin}. Using SPH simulations we show that the simulations where the binary opens a gap in the disk and the simulations where the disk does not have a gap are distributed in two well separate regions. Our analytic gap-opening criterion predicts a shape of the threshold between this two regions that is consistent with our simulations and the other ones in the literature. We propose an analogy between the regime without (with) a gap in the disk and the Type I (Type II) migration that is observed in simulations of planet-disk interaction (binaries with extreme mass ratios), emphasizing that the interaction that drives the formation of a gap on the disk is different in the regime that we analyze (comparable mass binary). △ Less

Submitted 21 November, 2012; v1 submitted 26 September, 2012; originally announced September 2012.

Journal ref: 2012, ApJ, 761, 31

Central regions of LIRGs: rings, hidden starbursts, Supernovae and star clusters

Authors: Petri Vaisanen, Andres Escala, Erkki Kankare, Jari Kotilainen, Seppo Mattila, Vinesh Rajpaul, Zara Randriamanakoto, Juha Reunanen, Stuart Ryder, Albert Zijlstra

Abstract: We study star formation (SF) in very active environments, in luminous IR galaxies, which are often interacting. A variety of phenomena are detected, such as central starbursts, circumnuclear SF, obscured SNe tracing the history of recent SF, massive super star clusters, and sites of strong off-nuclear SF. All of these can be ultimately used to define the sequence of triggering and propagation of s… ▽ More We study star formation (SF) in very active environments, in luminous IR galaxies, which are often interacting. A variety of phenomena are detected, such as central starbursts, circumnuclear SF, obscured SNe tracing the history of recent SF, massive super star clusters, and sites of strong off-nuclear SF. All of these can be ultimately used to define the sequence of triggering and propagation of star-formation and interplay with nuclear activity in the lives of gas rich galaxy interactions and mergers. In this paper we present analysis of high-spatial resolution integral field spectroscopy of central regions of two interacting LIRGs. We detect a nuclear 3.3 um PAH ring around the core of NGC 1614 with thermal-IR IFU observations. The ring's characteristics and relation to the strong star-forming ring detected in recombination lines are presented, as well as a scenario of an outward expanding starburst likely initiated with a (minor) companion detected within a tidal feature. We then present NIR IFU observations of IRAS 19115-2124, aka the Bird, which is an intriguing triple encounter. The third component is a minor one, but, nevertheless, is the source of 3/4 of the SFR of the whole system. Gas inflows and outflows are detected at the locations of the nuclei. Finally, we briefly report on our on-going NIR adaptive optics imaging survey of several dozen LIRGs. We have detected highly obscured core-collapse SNe in the central kpc, and discuss the statistics of "missing SNe" due to dust extinction. We are also determining the characteristics of hundreds of super star clusters in and around the core regions of LIRGs, as a function of host-galaxy properties. △ Less

Submitted 28 February, 2012; originally announced February 2012.

Comments: 8 pages, 7 figures, to appear in "The Central Kiloparsec in Galactic Nuclei: Astronomy at High Angular Resolution 2011", open access Journal of Physics: Conference Series (JPCS), published by IOP Publishing

Gravitational Fragmentation in Galaxy Mergers: A Stability Criteria

Authors: Andres Escala, Fernando Becerra, Luciano del Valle, Esteban Castillo

Abstract: We study the gravitational stability of gaseous streams in the complex environment of a galaxy merger, because mergers are known to be places of ongoing massive cluster formation and bursts of star formation. We find an analytic stability parameter for case of gaseous streams orbiting around the merger remnant. We test our stability criteria using hydrodynamical simulations of galaxy mergers, obta… ▽ More We study the gravitational stability of gaseous streams in the complex environment of a galaxy merger, because mergers are known to be places of ongoing massive cluster formation and bursts of star formation. We find an analytic stability parameter for case of gaseous streams orbiting around the merger remnant. We test our stability criteria using hydrodynamical simulations of galaxy mergers, obtaining satisfactory results. We find that our criteria successfully predicts the streams that will be gravitationally unstable to fragment into clumps. △ Less

Submitted 4 October, 2012; v1 submitted 6 February, 2012; originally announced February 2012.

Comments: 10 pages, 2 figures. After 2nd referee report

A law for star formation in galaxies

Authors: Andres Escala

Abstract: We study the galactic-scale triggering of star formation. We find that the largest mass-scale not stabilized by rotation, a well defined quantity in a rotating system and with clear dynamical meaning, strongly correlates with the star formation rate in a wide range of galaxies. We find that this relation can be understood in terms of self-regulation towards marginal Toomre stability and the amount… ▽ More We study the galactic-scale triggering of star formation. We find that the largest mass-scale not stabilized by rotation, a well defined quantity in a rotating system and with clear dynamical meaning, strongly correlates with the star formation rate in a wide range of galaxies. We find that this relation can be understood in terms of self-regulation towards marginal Toomre stability and the amount of turbulence allowed to sustain the system in this self-regulated quasi-stationary state. We test such an interpretation by computing the predicted star formation rates for a galactic interstellar medium characterized by lognormal probability distribution function and find good agreement with the observed relation. △ Less

Submitted 18 April, 2011; originally announced April 2011.

Comments: 13 pages, 2 figures, Accepted in ApJ

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

Authors: Lucio Mayer, Stelios Kazantzidis, Andres Escala, Simone Callegari

Abstract: 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… ▽ More 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. △ Less

Submitted 22 December, 2009; originally announced December 2009.

Comments: 26 pages, 4 Figures, submitted to Nature (includes Supplementary Information)

Four Highly Luminous Massive Star Forming Regions in the Norma Spiral Arm II. Deep NIR imaging

Authors: L. Chavarria, D. Mardones, G. Garay, A. Escala, L. Bronfman, S. Lizano

Abstract: We present sensitive NIR (J, H and K) imaging observations toward four luminous massive star forming regions in the Norma Spiral Arm: G324.201+0.119, G328.307+0.432, G329.337+0.147 and G330.949-0.174. We identify three clusters of young stellar objects (YSO) based on surface density diagnostics. We also find that sources detected only in the H and K-bands and with colors corresponding to spectra… ▽ More We present sensitive NIR (J, H and K) imaging observations toward four luminous massive star forming regions in the Norma Spiral Arm: G324.201+0.119, G328.307+0.432, G329.337+0.147 and G330.949-0.174. We identify three clusters of young stellar objects (YSO) based on surface density diagnostics. We also find that sources detected only in the H and K-bands and with colors corresponding to spectral types earlier than B2, are likely YSOs. We analyze the spatial distribution of stars of different masses and find signatures in two clusters of primordial mass segregation which can't be explained as due to incompleteness effects. We show that dynamic interactions of cluster members with the dense gas from the parent core can explain the observed mass segregation, indicating that the gas plays an important role in the dynamics of young clusters. △ Less

Submitted 25 November, 2009; originally announced November 2009.

Comments: Submitted to ApJ in September 2009

Journal ref: Astrophys.J.710:583-596,2010

On the global triggering mechanism of star formation in galaxies

Authors: Andres Escala

Abstract: We study the large-scale triggering of star formation in galaxies. We find that the largest mass-scale not stabilized by rotation, a well defined quantity in a rotating system and with clear dynamical meaning, strongly correlates with the star formation rate in a wide range of galaxies. We find that this relation can be explained in terms of the threshold for stability and the amount of turbulen… ▽ More We study the large-scale triggering of star formation in galaxies. We find that the largest mass-scale not stabilized by rotation, a well defined quantity in a rotating system and with clear dynamical meaning, strongly correlates with the star formation rate in a wide range of galaxies. We find that this relation can be explained in terms of the threshold for stability and the amount of turbulence allowed to sustain the system in equilibrium. Using this relation we also derived the observed correlation between the star formation rate and the luminosity of the brightest young stellar cluster. △ Less

Submitted 23 September, 2009; originally announced September 2009.

Comments: 11 pages, 1 figure. Submitted

The Evolution of Swift/BAT blazars and the origin of the MeV background

Authors: M. Ajello, L. Costamante, R. M. Sambruna, N. Gehrels, J. Chiang, A. Rau, A. Escala, J. Greiner, J. Tueller, J. V. Wall, R. F. Mushotzky

Abstract: We use 3 years of data from the Swift/BAT survey to select a complete sample of X-ray blazars above 15 keV. This sample comprises 26 Flat-Spectrum Radio Quasars (FSRQs) and 12 BL Lac objects detected over a redshift range of 0.03<z<4.0. We use this sample to determine, for the first time in the 15--55 keV band, the evolution of blazars. We find that, contrary to the Seyfert-like AGNs detected by… ▽ More We use 3 years of data from the Swift/BAT survey to select a complete sample of X-ray blazars above 15 keV. This sample comprises 26 Flat-Spectrum Radio Quasars (FSRQs) and 12 BL Lac objects detected over a redshift range of 0.03<z<4.0. We use this sample to determine, for the first time in the 15--55 keV band, the evolution of blazars. We find that, contrary to the Seyfert-like AGNs detected by BAT, the population of blazars shows strong positive evolution. This evolution is comparable to the evolution of luminous optical QSOs and luminous X-ray selected AGNs. We also find evidence for an epoch-dependence of the evolution as determined previously for radio-quiet AGNs. We interpret both these findings as a strong link between accretion and jet activity. In our sample, the FSRQs evolve strongly, while our best-fit shows that BL Lacs might not evolve at all. The blazar population accounts for 10-20 % (depending on the evolution of the BL Lacs) of the Cosmic X-ray background (CXB) in the 15--55 keV band. We find that FSRQs can explain the entire CXB emission for energies above 500 keV solving the mystery of the generation of the MeV background. The evolution of luminous FSRQs shows a peak in redshift ($z_c$=4.3$\pm0.5$) which is larger than the one observed in QSOs and X-ray selected AGNs. We argue that FSRQs can be used as tracers of massive elliptical galaxies in the early Universe. △ Less

Submitted 4 May, 2009; originally announced May 2009.

Comments: Accepted for publication in The Astrophysical Journal. 59 pages, 18 figures, 5 tables

Journal ref: Astrophys.J.699:603-625,2009

Stability of Galactic Gaseous Disks and the Formation of Massive Clusters

Authors: Andres Escala, Richard B. Larson

Abstract: We study gravitational instabilities in disks, with special attention to the most massive clumps that form because they are expected to be the progenitors of globular-type clusters. The maximum unstable mass is set by rotation and depends only on the surface density and orbital frequency of the disk. We propose that the formation of massive clusters is related to this largest scale in galaxies n… ▽ More We study gravitational instabilities in disks, with special attention to the most massive clumps that form because they are expected to be the progenitors of globular-type clusters. The maximum unstable mass is set by rotation and depends only on the surface density and orbital frequency of the disk. We propose that the formation of massive clusters is related to this largest scale in galaxies not stabilized by rotation. Using data from the literature, we predict that globular-like clusters can form in nuclear starburst disks and protogalactic disks but not in typical spiral galaxies, in agreement with observations. △ Less

Submitted 8 August, 2008; originally announced August 2008.

Comments: Accepted for publication in ApJ Letters. Minor comments added due to the referee's suggestion

Formation of Nuclear Disks and Supermassive Black Hole Binaries in Multi-Scale Hydrodynamical Galaxy Mergers

Authors: Lucio Mayer, Stelios Kazantzidis, Andres Escala

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

Submitted 22 July, 2008; originally announced July 2008.

Comments: 12 pages, 5 figures. Review to appear in the proceedings of the Conference "The Central Kiloparsec: Active Galactic Nuclei and Their Hosts", Ierapetra 4-6 June 2008, Greece, Memorie della Societa Astronomica Italiana

Stability of Galactic Gas Disks and the Formation of Massive Clusters

Authors: Andres Escala, Richard B. Larson

Abstract: We study gravitational instabilities in disks, with special attention to the most massive clumps that form because they are expected to be the progenitors of globular-type clusters. The maximum unstable mass is set by rotation and depends only on the surface density and orbital frequency of the disk. We propose that the formation of massive clusters is related to this largest scale in galaxies n… ▽ More We study gravitational instabilities in disks, with special attention to the most massive clumps that form because they are expected to be the progenitors of globular-type clusters. The maximum unstable mass is set by rotation and depends only on the surface density and orbital frequency of the disk. We propose that the formation of massive clusters is related to this largest scale in galaxies not stabilized by rotation. Using data from the literature, we predict that globular-like clusters can form in nuclear starburst disks and protogalactic disks but not in typical spiral galaxies, in agreement with observations. △ Less

Submitted 4 June, 2008; originally announced June 2008.

Comments: 12 pages, 2 figures, submitted to ApjL

Towards a Comprehensive Fueling-Controlled Theory on the Growth of Massive Black Holes and Host Spheroids

Authors: Andres Escala

Abstract: We study the relation between nuclear massive black holes and their host spheroid gravitational potential. Using AMR numerical simulations, we analyze how gas is transported in the nuclear (central kpc) regions of galaxies. We study the gas fueling onto the inner accretion disk (sub-pc scale) and the star formation in a massive nuclear disk like those generally found in proto-spheroids (ULIRGs,… ▽ More We study the relation between nuclear massive black holes and their host spheroid gravitational potential. Using AMR numerical simulations, we analyze how gas is transported in the nuclear (central kpc) regions of galaxies. We study the gas fueling onto the inner accretion disk (sub-pc scale) and the star formation in a massive nuclear disk like those generally found in proto-spheroids (ULIRGs, SCUBA Galaxies). These sub-pc resolution simulation of gas fueling that is mainly depleted by star formation naturally satisfy the `M_BH - $M_virial' relation, with a scatter considerably less than the observed one. We found a generalized version of Kennicutt-Schmidt Law for starbursts is satisfied, in which the total gas depletion rate (dot{M}_gas = dot{M}_BH + dot{M}_SF) is the one that scales as M_gas/t_orbital. We also found that the `M_BH - sigma' relation is a byproduct of the `M_BH - M_virial' relation in the fueling controlled scenario. △ Less

Submitted 30 May, 2007; originally announced May 2007.

Comments: 12 pages, figures, submited to ApJ, email: andres@slac.stanford.edu

Journal ref: Astrophys.J.671:1264-1271,2007

On the Fueling of Massive Black Holes and the Properties of their Host Spheroids

Authors: Andres Escala

Abstract: We study the relation between nuclear massive black holes and their host spheroid gravitational potential. Using simple models, we analyze how gas is expected to be transported in the nuclear regions of galaxies. When we couple it with the expected gas lifetime given by the Kennicutt-Schmidt Law, this naturally leads to the `M_BH - M_virial' and `M_BH - sigma' relations. We also numerically test… ▽ More We study the relation between nuclear massive black holes and their host spheroid gravitational potential. Using simple models, we analyze how gas is expected to be transported in the nuclear regions of galaxies. When we couple it with the expected gas lifetime given by the Kennicutt-Schmidt Law, this naturally leads to the `M_BH - M_virial' and `M_BH - sigma' relations. We also numerically test, using AMR simulations, our simple models for the mass transport with satisfactory results. △ Less

Submitted 31 August, 2006; v1 submitted 24 April, 2006; originally announced April 2006.

Comments: Reduced version that fit the 4-page limit. In press ApJ Letters

Journal ref: Astrophys.J.648:L13-L16,2006

The Role of Gas in the Merging of Massive Black Holes in Galactic Nuclei. II. Black Hole Merging in a Clumpy Disk

Authors: Andres Escala, Richard B. Larson, Paolo S. Coppi, Diego Mardones

Abstract: Using SPH numerical simulations, we investigate the effects of gas on the inspiral and merger of a massive black hole binary. This study is motivated by the very massive nuclear gas disks observed in the central regions of merging galaxies. Here we present results that expand on the treatment in a previous work (Escala, Larson, Coppi & Mardones 2004; henceforth Paper I), by studying more realist… ▽ More Using SPH numerical simulations, we investigate the effects of gas on the inspiral and merger of a massive black hole binary. This study is motivated by the very massive nuclear gas disks observed in the central regions of merging galaxies. Here we present results that expand on the treatment in a previous work (Escala, Larson, Coppi & Mardones 2004; henceforth Paper I), by studying more realistic models. We run a variety of models, ranging from simulations with a relatively smooth gas disk to cases in which the gas has a more clumpy spatial distribution. We also vary the inclination angle between the plane of the binary and the plane of the disk, and the mass ratio between the MBHs and the gaseous disk. We find that as in Paper I, in the early evolution of the system the binary separation diminishes due to gravitational drag, and in the later stages the medium responds by forming an ellipsoidal density enhancement whose axis lags behind the binary axis; this offset produces a torque on the binary that causes continuing loss of angular momentum and is able to reduce the separation to distances where gravitational radiation is efficient. In the simulations that we perform, we find that the coalescence timescale for the MBH binary varies between 5 10^6 yr and 2.5 10^7 yr. For MBHs that satisfy the observed `m-sigma_c' relation, we predict that in a merger of galaxies that have at least 1% of their total mass in gas, the MBHs will coalesce soon after the galaxies merge. We also predict that if the MBHs depart from the `m-sigma_c' relation, the MBH binary can create a circumbinary gap in the disk that stalls the coalescence, but this gap formation can act as a self-regulatory mechanism on MBH growth that can help to explain the existence of the `m-sigma_c' relation. △ Less

Submitted 11 June, 2004; originally announced June 2004.

Comments: 31 pages, 14 figures, submitted to ApJ

Journal ref: Astrophys.J. 630 (2005) 152-166

The Role of Gas in the Merging of Massive Black Holes in Galactic Nuclei. I. Black Hole Merging in a Spherical Gas Cloud

Authors: Andres Escala, Richard B. Larson, Paolo S. Coppi, Diego Mardones

Abstract: Using high-resolution SPH numerical simulations, we investigate the effects of gas on the inspiral and merger of a massive black hole binary. This study is motivated by both observational and theoretical work that indicate the presence of large amounts of gas in the central regions of merging galaxies. N-body simulations have shown that the coalescence of a massive black hole binary eventually s… ▽ More Using high-resolution SPH numerical simulations, we investigate the effects of gas on the inspiral and merger of a massive black hole binary. This study is motivated by both observational and theoretical work that indicate the presence of large amounts of gas in the central regions of merging galaxies. N-body simulations have shown that the coalescence of a massive black hole binary eventually stalls in a stellar background. However, our simulations suggest that the massive black hole binary will finally merge if it is embedded in a gaseous background. Here we present results in which the gas is assumed to be initially spherical with a relatively smooth distribution. In the early evolution of the binary, the separation dimishes due to the gravitational drag exerted by the background gas. In the later stages, when the binary dominates the gravitational potential in its vicinity, the medium responds by forming an ellipsoidal density enhancement whose axis lags behind the binary axis, and this offset produces a torque on the binary that causes continuing loss of angular momentum and is able to reduce the binary separation to distances where gravitational radiation is efficient. Assuming typical parameters from observations of Ultra Luminous Infrared Galaxies, we predict that a black hole binary will merge within $10^{7}$yrs; therefore these results imply that in a merger of gas-rich galaxies, any massive central black holes will coalescence soon after the galaxies merge. Our work thus supports scenarios of massive black hole evolution and growth where hierarchical merging plays an important role. The final coalescence of the black holes leads to gravitational radiation emission that would be detectable up to high redshift by LISA. We show that similar physical effects are important for the formation of close binary stars. △ Less

Submitted 10 November, 2003; v1 submitted 29 October, 2003; originally announced October 2003.

Comments: 38 pages, 14 figures, submitted to ApJ

Journal ref: Astrophys.J. 607 (2004) 765-777