Source author record

Andres Escala

Andres Escala 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

18works
8topics
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

18 published item(s)

preprint2025arXiv

Rapid formation of a very massive star >50000 $M_\odot$ and subsequently an IMBH from runaway collisions. Direct N-body and Monte Carlo simulations of dense star clusters

Context. We present simulations of a massive young star cluster using \textsc{Nbody6++GPU} and \textsc{MOCCA}. The cluster is initially more compact than previously published models, with one million stars, a total mass of $5.86 \times 10^5~\mathrm{M}_{\odot}$, and a half-mass radius of $0.1~\mathrm{pc}$. Aims. We analyse the formation and growth of a very massive star (VMS) through successive stellar collisions and investigate the subsequent formation of an intermediate-mass black hole (IMBH) in the core of a dense star cluster. Methods. We use both direct \textit{N}-body and Monte Carlo simulations, incorporating updated stellar evolution prescriptions (SSE/BSE) tailored to massive stars and VMSs. These include revised treatments of stellar radii, rejuvenation, and mass loss during collisions. While the prescriptions represent reasonable extrapolations into the VMS regime, the internal structure and thermal state of VMSs formed through stellar collisions remain uncertain, and future work may require further refinement. Results. We find that runaway stellar collisions in the cluster core produce a VMS exceeding $5 \times 10^4~\mathrm{M}_{\odot}$ within 5 Myr, which subsequently collapses into an IMBH. Conclusions. Our model suggests that dense stellar environments may enable the formation of very massive stars and massive black hole seeds through runaway stellar collisions. These results provide a potential pathway for early black hole growth in star clusters and offer theoretical context for interpreting recent JWST observations of young, compact clusters at high redshift.

preprint2020arXiv

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

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.

preprint2020arXiv

Universal Relation for Life-span Energy Consumption in Living Organisms: Insights for the origin of ageing

Metabolic energy consumption has long been thought to play a major role in the aging process ({\it 1}). Across species, a gram of tissue on average expends about the same amount of energy during life-span ({\it 2}). Energy restriction has also been shown that increases maximum life-span ({\it 3}) and retards age-associated changes ({\it 4}). However, there are significant exceptions to a universal energy consumption during life-span, mainly coming from the inter-class comparison ({\it 5, 6}). Here we present a unique relation for life-span energy consumption, valid for $\sim$300 species representing all classes of living organisms, from unicellular ones to the largest mammals. The relation has an average scatter of only 0.3 dex, with 95\% ($\rm 2-σ$) of the organisms having departures less than a factor of $π$ from the relation, despite the $\sim$20 orders of magnitude difference in body mass, reducing any possible inter-class variation in the relation to only a geometrical factor. This result can be interpreted as supporting evidence for the existence of an approximately constant total number $\rm N_r \sim 10^8$ of respiration cycles per lifetime for all organisms, effectively predetermining the extension of life by the basic energetics of respiration, being an incentive for future studies that investigate the relation of such constant $\rm N_r$ cycles per lifetime with the production rates of free radicals and oxidants, which may give definite constraints on the origin of ageing.

preprint2020arXiv

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

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.

preprint2016arXiv

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

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.

preprint2016arXiv

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

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.

preprint2015arXiv

On the Functional Form of the Universal Star Formation Law

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.

preprint2015arXiv

Super massive black holes in star forming gaseous circumnuclear discs

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.

preprint2014arXiv

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

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.

preprint2013arXiv

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

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.

preprint2013arXiv

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

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.

preprint2013arXiv

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

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.

preprint2012arXiv

Binary-disk interaction: Gap-Opening criteria

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

preprint2012arXiv

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

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.

preprint2012arXiv

Gravitational Fragmentation in Galaxy Mergers: A Stability Criteria

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.

preprint2011arXiv

A law for star formation in galaxies

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.

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.

preprint2006arXiv

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

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.