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, $\Gamma$, 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 $\lambda_c$. We define the spatial scale $\lambda_{\rm eq}$ at which galactic rotation and non-circular motions contribute equally to $\Gamma$. 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 $\lambda_c$ and this transition is necessary to reproduce the circulation distribution. We find that $\lambda_{\rm eq}$, and therefore the structure of the gas velocity field, is set by the local conditions of gravitational stability and stellar feedback.

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 collision 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 \, M_O$ 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 \epsilon_{BH} = \frac{M_{BH}}{M_{CMO}}$, with these efficiencies $\rm \epsilon_{BH}$ drastically growing for central massive objects more massive than $\rm \sim 10^7 \,M_O$, approaching unity around $\rm M_{CMO} \sim 10^8 \, M_O$. We show that the proposed scenario successfully explain the relative trends observed in the masses, efficiencies and scaling relations between massive black holes and nuclear stellar cluster

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.

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$\alpha$-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$\alpha$. 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 \sigma_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 $\tau_{\rm dep}=2.3 \pm 1.2$\,Gyr for the galaxy as a whole. This value is in agreement with the average $\tau_{\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.

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

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\'iso 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}/\sigma_v$, corrected for inclination). We find a dependence between $V_{\rm rot}/\sigma_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}/\sigma_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.

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.

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.

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) $\Sigma_{\rm SFR} \propto \sqrt{G/L}\Sigma_{\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 $\Omega$ 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\Omega 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.

We present adaptive optics (AO) assisted SINFONI integral field unit (IFU) spectroscopy of eleven H$\alpha$ 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$\alpha$ emission line ratio show a median negative gradient for the SHiZELS survey of $\Delta$log(O/H)/$\Delta$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.

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.

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 $\alpha$ parameter is dominated by the Reynolds term, $\alpha_R$, with $\alpha_R\gg 1$. The gravitational part of the $\alpha$ parameter, $\alpha_G$, has an increasing trend toward the galactic center at higher redshifts, with values $\alpha_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.

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.

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{\Sigma}_{\star} \propto \sqrt{\frac{G}{L}}\Sigma_{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{\Sigma}_{\star}$ scaling, that are widely studied in the literature: $\rm \Sigma_{gas}^{3/2}, \Sigma_{gas}/t_{orb}, \Sigma_{gas}/t_{ff} \, and \, \Sigma_{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.

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=\epsilon\rho_{\rm gas}\Delta 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.

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.

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.

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.

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.

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

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.

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.

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.

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

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

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.

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.

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.

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

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

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.

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.

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.

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.