The centers of most galaxies contain massive black holes surrounded by dense star clusters. The structure of these clusters determines the rate and properties of observable transient events, such as flares from tidally disrupted stars and gravitational-wave signals from stars spiraling into the black hole. Most estimates of these rates enforce spherical symmetry on the cluster. Here we show that, in the course of generic evolutionary processes, a star cluster surrounding a black hole can undergo a robust phase transition from a spherical thermal equilibrium to a lopsided equilibrium, in which most stars are on high-eccentricity orbits with aligned orientations. The rate of transient events is expected to be much higher in the ordered phase. Better models of cluster formation and evolution are needed to determine whether clusters should be found in the ordered or disordered phase.

The energy transfer among the components in a gas determines its fate. Especially at low temperatures, inelastic collisions drive the cooling and the heating mechanisms. In the early Universe as well as in zero- or low- metallicity environments the major contribution comes from the collisions among atomic and molecular hydrogen, also in its deuterated version. The present work shows some updated calculations of the H$_2$ cooling function based on novel collisional data which explicitely take into account the reactive pathway at low temperatures. Deviations from previous calculations are discussed and a multivariate data analysis is performed to provide a fit depending on both the gas temperature and the density of the gas.

We present idealized models of a razor-thin, axisymmetric, Keplerian stellar disc around a massive black hole, and study non-axisymmetric secular instabilities in the absence of either counter-rotation or loss cones. These discs are prograde mono-energetic waterbags, whose phase space distribution functions are constant for orbits within a range of eccentricities (e) and zero outside this range. The linear normal modes of waterbags are composed of sinusoidal disturbances of the edges of distribution function in phase space. Waterbags which include circular orbits (polarcaps) have one stable linear normal mode for each azimuthal wavenumber m. The m = 1 mode always has positive pattern speed and, for polarcaps consisting of orbits with e < 0.9428, only the m = 1 mode has positive pattern speed. Waterbags excluding circular orbits (bands) have two linear normal modes for each m, which can be stable or unstable. We derive analytical expressions for the instability condition, pattern speeds, growth rates and normal mode structure. Narrow bands are unstable to modes with a wide range in m. Numerical simulations confirm linear theory and follow the non-linear evolution of instabilities. Long-time integration suggests that instabilities of different m grow, interact non-linearly and relax collisionlessly to a coarse-grained equilibrium with a wide range of eccentricities.

Non-thermal photons deriving from radiative transitions among the internal ladder of atoms and molecules are an important source of photons in addition to thermal and stellar sources in many astrophysical environments. In the present work the calculation of reaction rates for the direct photodissociation of some molecules relevant in early Universe chemistry is presented; in particular, the calculations include non-thermal photons deriving from the recombination of primordial hydrogen and helium atoms for the cases of H2, HD and HeH+. New effects on the fractional abundances of chemical species are investigated and the fits for the HeH+ photodissociation rates by thermal photons are provided.

We apply photo-dissociation region (PDR) molecular line emission models, that have varying degrees of enhanced mechanical heating rates, to the gaseous component of simulations of star-forming galaxies taken from the literature. Snapshots of these simulations are used to produce line emission maps for the rotational transitions of the CO molecule and its 13CO isotope up to J = 4-3. We consider two galaxy models: a small disk galaxy of solar metallicity and a lighter dwarf galaxy with 0.2 \zsun metallicity. Elevated excitation temperatures for CO(1 - 0) correlate positively with mechanical feedback, that is enhanced towards the central region of both model galaxies. The emission maps of these model galaxies are used to compute line ratios of CO and 13CO transitions. These line ratios are used as diagnostics where we attempt to match them These line ratios are used as diagnostics where we attempt to match them to mechanically heated single component (i.e. uniform density, Far-UV flux, visual extinction and velocity gradient) equilibrium PDR models. We find that PDRs ignoring mechanical feedback in the heating budget over-estimate the gas density by a factor of 100 and the far-UV flux by factors of ~10 - 1000. In contrast, PDRs that take mechanical feedback into account are able to fit all the line ratios for the central < 2 kpc of the fiducial disk galaxy quite well. The mean mechanical heating rate per H atom that we recover from the line ratio fits of this region varies between $10^{-27}$ -- $10^{-26}$~erg s$^{-1}$.

Far-infrared molecular emission is an important tool used to understand the excitation mechanisms of the gas in the inter-stellar medium of star-forming galaxies. In the present work, we model the emission from rotational transitions with critical densities n >~ 10^4 cm-3. We include 4-3 < J <= 15-14 transitions of CO and 13CO, in addition to J <= 7-6 transitions of HCN, HNC, and HCO+ on galactic scales. We do this by re-sampling high density gas in a hydrodynamic model of a gas-rich disk galaxy, assuming that the density field of the interstellar medium of the model galaxy follows the probability density function (PDF) inferred from the resolved low density scales. We find that in a narrow gas density PDF, with a mean density of ~10 cm-3 and a dispersion \sigma = 2.1 in the log of the density, most of the emission of molecular lines, emanates from the 10-1000 cm-3 part of the PDF. We construct synthetic emission maps for the central 2 kpc of the galaxy and fit the line ratios of CO and 13CO up to J = 15-14, as well as HCN, HNC, and HCO+ up to J = 7-6, using one photo-dissociation region (PDR) model. We attribute the goodness of the one component fits for our model galaxy to the fact that the distribution of the luminosity, as a function of density, is peaked at gas densities between 10 and 1000 cm-3. We explore the impact of different log-normal density PDFs on the distribution of the line-luminosity as a function of density, and we show that it is necessary to have a broad dispersion, corresponding to Mach numbers >~ 30 in order to obtain significant emission from n > 10^4 cm-3 gas. Such Mach numbers are expected in star-forming galaxies, LIRGS, and ULIRGS. By fitting line ratios of HCN(1-0), HNC(1-0), and HCO+(1-0) for a sample of LIRGS and ULIRGS using mechanically heated PDRs, we constrain the Mach number of these galaxies to 29 < M < 77.

CO observations in active galactic nuclei and star-bursts reveal high kinetic temperatures. Those environments are thought to be very turbulent due to dynamic phenomena such as outflows and high supernova rates. We investigate the effect of mechanical heating (MH) on atomic fine-structure and molecular lines, and their ratios. We use those ratios as a diagnostic to constrain the amount of MH in an object and also study its significance on estimating the H2 mass. Equilibrium PDRs models were used to compute the thermal and chemical balance for the clouds. The equilibria were solved for numerically using the optimized version of the Leiden PDR-XDR code. Large velocity gradient calculations were done as post-processing on the output of the PDR models using RADEX. High-J CO line ratios are very sensitive to MH. Emission becomes at least one order of magnitude brighter in clouds with n~10^5~cm^-3 and a star formation rate of 1 Solar Mass per year (corresponding to a MH rate of 2 * 10^-19 erg cm^-3 s^-1). Emission of low-J CO lines is not as sensitive to MH, but they do become brighter in response to MH. Generally, for all of the lines we considered, MH increases excitation temperatures and decreases the optical depth at the line centre. Hence line ratios are also affected, strongly in some cases. Ratios involving HCN are a good diagnostic for MH, such as HCN(1-0)/CO(1-0) and HCN(1-0)/HCO^+(1-0). Both ratios increase by a factor 3 or more for a MH equivalent to > 5 percent of the surface heating, as opposed to pure PDRs. The first major conclusion is that low-J to high-J intensity ratios will yield a good estimate of the MH rate (as opposed to only low-J ratios). The second one is that the MH rate should be taken into account when determining A_V or equivalently N_H, and consequently the cloud mass. Ignoring MH will also lead to large errors in density and radiation field estimates.

Starburst galaxies are undergoing intense periods of star formation. Understanding the heating and cooling mechanisms in these galaxies can give us insight to the driving mechanisms that fuel the starburst. Molecular emission lines play a crucial role in the cooling of the excited gas. With SPIRE on the Herschel Space Observatory we have observed the rich molecular spectrum towards the central region of NGC 253. CO transitions from J=4-3 to 13-12 are observed and together with low-J line fluxes from ground based observations, these lines trace the excitation of CO. By studying the CO excitation ladder and comparing the intensities to models, we investigate whether the gas is excited by UV radiation, X-rays, cosmic rays, or turbulent heating. Comparing the $^{12}$CO and $^{13}$CO observations to large velocity gradient models and PDR models we find three main ISM phases. We estimate the density, temperature,and masses of these ISM phases. By adding $^{13}$CO, HCN, and HNC line intensities, we are able to constrain these degeneracies and determine the heating sources. The first ISM phase responsible for the low-J CO lines is excited by PDRs, but the second and third phases, responsible for the mid to high-J CO transitions, require an additional heating source. We find three possible combinations of models that can reproduce our observed molecular emission. Although we cannot determine which of these are preferable, we can conclude that mechanical heating is necessary to reproduce the observed molecular emission and cosmic ray heating is a negligible heating source. We then estimate the mass of each ISM phase; $6\times 10^7$ M$_\odot$ for phase 1 (low-J CO lines), $3\times 10^7$ M$_\odot$ for phase 2 (mid-J CO lines), and $9\times 10^6$ M$_\odot$ for phase 3 (high-J CO lines) for a total system mass of $1\times10^{8}$ M$_\odot$.

It is strongly believed that Andromeda's double nucleus signals a disk of stars revolving around its central super-massive black hole on eccentric Keplerian orbits with nearly aligned apsides. A self-consistent stellar dynamical origin for such apparently long-lived alignment has so far been lacking, with indications that cluster self-gravity is capable of sustaining such lopsided configurations if and when stimulated by external perturbations. Here, we present results of N-body simulations which show unstable counter-rotating stellar clusters around super-massive black holes saturating into uniformly precessing lopsided nuclei. The double nucleus in our featured experiment decomposes naturally into a thick eccentric disk of apo-apse aligned stars which is embedded in a lighter triaxial cluster. The eccentric disk reproduces key features of Keplerian disk models of Andromeda's double nucleus; the triaxial cluster has a distinctive kinematic signature which is evident in HST observations of Andromeda's double nucleus, and has been difficult to reproduce with Keplerian disks alone. Our simulations demonstrate how the combination of eccentric disk and triaxial cluster arises naturally when a star cluster accreted over a pre-existing and counter-rotating disk of stars, drives disk and cluster into a mutually destabilizing dance. Such accretion events are inherent to standard galaxy formation scenarios. They are here shown to double stellar black hole nuclei as they feed them.

Context. Multitransition CO observations of galaxy centers have revealed that significant fractions of the dense circumnuclear gas have high kinetic temperatures, which are hard to explain by pure photon excitation, but may be caused by dissipation of turbulent energy. Aims. We aim to determine to what extent mechanical heating should be taken into account while modelling PDRs. To this end, the effect of dissipated turbulence on the thermal and chemical properties of PDRs is explored. Methods. Clouds are modelled as 1D semi-infinite slabs whose thermal and chemical equilibrium is solved for using the Leiden PDR-XDR code. Results. In a steady-state treatment, mechanical heating seems to play an important role in determining the kinetic temperature of the gas in molecular clouds. Particularly in high-energy environments such as starburst galaxies and galaxy centers, model gas temperatures are underestimated by at least a factor of two if mechanical heating is ignored. The models also show that CO, HCN and H2 O column densities increase as a function of mechanical heating. The HNC/HCN integrated column density ratio shows a decrease by a factor of at least two in high density regions with n \sim 105 cm-3, whereas that of HCN/HCO+ shows a strong dependence on mechanical heating for this same density range, with boosts of up to three orders of magnitude. Conclusions. The effects of mechanical heating cannot be ignored in studies of the molecular gas excitation whenever the ratio of the star formation rate to the gas density is close to, or exceeds, 7 \times 10-6 M yr-1 cm4.5 . If mechanical heating is not included, predicted column densities are underestimated, sometimes even by a few orders of magnitude. As a lower bound to its importance, we determined that it has non-negligible effects already when mechanical heating is as little as 1% of the UV heating in a PDR.

We show that the algorithm proposed by Gauss to compute the secular evolution of gravitationally interacting Keplerian rings extends naturally to softened gravitational interactions. The resulting tool is ideal for the study of the secular dynamical evolution of nearly Keplerian systems such as stellar clusters surrounding black holes in galactic nuclei, cometary clouds, or planetesimal discs. We illustrate its accuracy, efficiency and versatility on a variety of configurations. In particular, we examine a secularly unstable unstable system of counter-rotating disks, and follow the unfolding and saturation of the instability into a global, uniformly precessing, lopsided (m=1) mode.