The Milky Way's central molecular zone (CMZ) has emerged in recent years as a unique laboratory for the study of star formation. Here we use the simulations presented in Tress et al. 2020 to investigate star formation in the CMZ. These simulations resolve the structure of the interstellar medium at sub-parsec resolution while also including the large-scale flow in which the CMZ is embedded. Our main findings are as follows. (1) While most of the star formation happens in the CMZ ring at $R\gtrsim100 \pc$, a significant amount also occurs closer to SgrA* at $R \lesssim 10\pc$. (2) Most of the star formation in the CMZ happens downstream of the apocentres, consistent with the "pearls-on-a-string" scenario, and in contrast to the notion that an absolute evolutionary timeline of star formation is triggered by pericentre passage. (3) Within the timescale of our simulations ($\sim100$~Myr), the depletion time of the CMZ is constant within a factor of $\sim2$. This suggests that variations in the star formation rate are primarily driven by variations in the mass of the CMZ, caused for example by AGN feedback or externally-induced changes in the bar-driven inflow rate, and not by variations in the depletion time. (4) We study the trajectories of newly born stars in our simulations. We find several examples that have age and 3D velocity compatible with those of the Arches and Quintuplet clusters. Our simulations suggest that these prominent clusters originated near the collision sites where the bar-driven inflow accretes onto the CMZ, at symmetrical locations with respect to the Galactic centre, and that they have already decoupled from the gas in which they were born.

It has been recently suggested that (i) nuclear rings in barred galaxies (including our own Milky Way) form at the radius where the shear parameter of the rotation curve reaches a minimum; (ii) the acoustic instability of Montenegro et al. is responsible for driving the turbulence and angular momentum transport in the central regions of barred galaxies. Here we test these suggestions by running simple hydrodynamical simulations in a logarithmic barred potential. Since the rotation curve of this potential is scale-free, the shear minimum theory predicts that no ring should form. We find that in contrast to this prediction, a ring does form in the simulation, with morphology consistent with that of nuclear rings in real barred galaxies. This proves that the presence of a shear-minimum is not a necessary condition for the formation of a ring. We also find that perturbations that are predicted to be acoustically unstable wind up and eventually propagate off to infinity, so that the system is actually stable. We conclude that (i) the shear-minimum theory is an unlikely mechanism for the formation of nuclear rings in barred galaxies; (ii) the acoustic instability is a spurious result and may not be able to drive turbulence in the interstellar medium, at least for the case without self-gravity. The question of the role of turbulent viscosity remains open.

We use hydrodynamical simulations to study the Milky Way's central molecular zone (CMZ), i.e. the star-forming nuclear ring at Galactocentric radii $R\lesssim200$ pc. The simulations comprise the gas flow in a Milky Way barred potential out to $R=5$ kpc, which is necessary in order to capture the large-scale environment in which the CMZ is embedded and with which it is strongly interacting through the bar-driven inflow. The simulations also include a non-equilibrium time-dependent chemical network, gas self-gravity, and a sub-grid model for star formation and supernova feedback, all while reaching sub-parsec resolution in the densest regions. Our main findings are as follows: (1) The distinction between inner ($R\lesssim120$ pc) and outer ($120\lesssim R\lesssim450$ pc) CMZ that is sometimes proposed in the literature is unnecessary. Instead, the CMZ is best described as single structure, namely a star-forming ring with outer radius $R\simeq 200$ pc which is interacting directly with the dust lanes that mediate the bar-driven inflow. (2) This accretion can induce a significant tilt of the CMZ out of the plane. A tilted CMZ might provide an alternative explanation to the $\infty$-shaped structure identified in Herschel data by Molinari et al. 2011. (3) The bar in our simulation efficiently drives an inflow from the Galactic disc ($R\simeq 3$ kpc) down to the CMZ ($R\simeq200$ pc) of the order of $1\rm\,M_\odot\,yr^{-1}$, consistent with observational determinations. (4) Self-gravity and supernovae feedback can drive an inflow from the CMZ inwards towards the circumnuclear disc of the order of $\sim0.03\,\rm M_\odot\,yr^{-1}$. (5) We give a new interpretation for the 3D placement of the 20 and 50 km s$^{-1}$ clouds, according to which they are close ($R\lesssim30$ pc) to the Galactic centre, but are also connected to the larger-scale streams at $R\gtrsim100$ pc.

We introduce a new suite of simulations, "The Cloud Factory", which self-consistently forms molecular cloud complexes at high enough resolution to resolve internal substructure (up to 0.25 Msol in mass) all while including galactic-scale forces. We use a version of the Arepo code modified to include a detailed treatment of the physics of the cold molecular ISM, and an analytical galactic gravitational potential for computational efficiency. The simulations have nested levels of resolution, with the lowest layer tied to tracer particles injected into individual cloud complexes. These tracer refinement regions are embedded in the larger simulation so continue to experience forces from outside the cloud. This allows the simulations to act as a laboratory for testing the effect of galactic environment on star formation. Here we introduce our method and investigate the effect of galactic environment on filamentary clouds. We find that cloud complexes formed after a clustered burst of feedback, have shorter lengths and are less likely to fragment compared to quiescent clouds (e.g. the Musca filament) or those dominated by the galactic potential (e.g. Nessie). Spiral arms and differential rotation preferentially align filaments, but strong feedback randomises them. Long filaments formed within the cloud complexes are necessarily coherent with low internal velocity gradients, which has implications for the formation of filamentary star-clusters. Cloud complexes formed in regions dominated by supernova feedback have fewer star-forming cores, and these are more widely distributed. These differences show galactic-scale forces can have a significant impact on star formation within molecular clouds.

We present here the first of a series of papers aimed at better understanding the evolution and properties of giant molecular clouds (GMCs) in a galactic context. We perform high resolution, three-dimensional {\sc arepo} simulations of an interacting galaxy inspired by the well-observed M51 galaxy. Our fiducial simulations include a non-equilibrium, time-dependent, chemical network that follows the evolution of atomic and molecular hydrogen as well as carbon and oxygen self-consistently. Our calculations also treat gas self-gravity and subsequent star formation (described by sink particles), and coupled supernova feedback. In the densest parts of the simulated interstellar medium (ISM) we reach sub-parsec resolution, granting us the ability to resolve individual GMCs and their formation and destruction self-consistently throughout the galaxy. In this initial work we focus on the general properties of the ISM with a particular focus on the cold star-forming gas. We discuss the role of the interaction with the companion galaxy in generating cold molecular gas and controlling stellar birth. We find that while the interaction drives large-scale gas flows and induces spiral arms in the galaxy, it is of secondary importance in determining gas fractions in the different ISM phases and the overall star-formation rate. The behaviour of the gas on small GMC scales instead is mostly controlled by the self-regulating property of the ISM driven by coupled feedback.

We explore the origin of phase-space substructures revealed by the second Gaia data release in the disc of the Milky Way, such as the ridges in the $V_{\phi}$-$r$ plane, the undulations in the $V_{\phi}$-$r$-$V_r$ space and the streams in the $V_{\phi}$-$V_r$ plane. We use a collisionless N-body simulation with co-spatial thin and thick discs, along with orbit integration, to study the orbital structure close to the Outer Lindblad Resonance (OLR) of the bar. We find that a prominent, long-lived ridge is formed in the $V_{\phi}$-$r$ plane due to the OLR which translates to streams in the $V_{\phi}$-$V_r$ plane and examine which closed periodic and trapped librating orbits are responsible for these features. We find that orbits which carry out small librations around the $x_1(1)$ family are preferentially found at negative $V_r$, giving rise to a `horn'-like feature, while orbits with larger libration amplitudes, trapped around the $x_1(2)$ and $x_1(1)$ families, constitute the positive $V_r$ substructure, i.e. the Hercules-like feature. This changing libration amplitude of orbits will translate to a changing ratio of thin/thick disc stars, which could have implications on the metallicity distribution in this plane. We find that a scenario in which the Sun is placed close to the OLR gives rise to a strong asymmetry in $V_r$ in the $V_{\phi}$-$V_r$ plane (i.e. Hercules vs. `the horn') and subsequently to undulations in the $V_{\phi}$-$r$-$V_r$ space. We also explore a scenario in which the Sun is placed closer to the bar corotation and find that the bar perturbation $alone$ cannot give rise to the these features.

Observations of molecular gas near the Galactic centre ($| l | < 10^\circ$, $| b | < 1^\circ$) reveal the presence of a distinct population of enigmatic compact clouds which are characterised by extreme velocity dispersions ($\Delta v > 100\, \rm km/s$). These Extended Velocity Features (EVFs) are very prominent in the datacubes and dominate the kinematics of molecular gas just outside the Central Molecular Zone (CMZ). The prototypical example of such a cloud is Bania Clump 2. We show that similar features are naturally produced in simulations of gas flow in a realistic barred potential. We analyse the structure of the features obtained in the simulations and use this to interpret the observations. We find that the features arise from collisions between material that has been infalling rapidly along the dust lanes of the Milky Way bar and material that belongs to one of the following two categories: (i) material that has `overshot' after falling down the dust lanes on the opposite side; (ii) material which is part of the CMZ. Both types of collisions involve gas with large differences in the line-of-sight velocities, which is what produces the observed extreme velocity dispersions. Examples of both categories can be identified in the observations. If our interpretation is correct, we are directly witnessing (a) collisions of clouds with relative speeds of $\sim 200\, \rm km/s$ and (b) the process of accretion of fresh gas onto the CMZ.

Observations show that (i) multiple gas phases can coexist in the atmospheres of galaxies and clusters; (ii) these atmospheres may be significantly rotating in the inner parts, with typical velocities that approach or even exceed the local sound speed. The thermal instability is a natural candidate to explain the formation of cold structures via condensation of a hotter gas phase. Here we systematically study the effect of rotation on the thermal stability of stratified plane-parallel atmospheres, using both analytical arguments and numerical simulations. We find that the formation of cold structures starting from small isobaric perturbations is enhanced in the regions where the rotation of the system is dynamically important (i.e. when the rotational velocity becomes comparable to the sound speed). In particular, the threshold value of the ratio between the cooling and dynamical time $t_{\rm cool}/t_{\rm dyn}$ below which condensations can form is increased by a factor up to $\sim 10$ in the presence of significant rotation. We briefly discuss the implications of our results for galaxies and clusters.

Whether High Velocity Clouds (HVCs) can form by condensation of the hot ($T \sim 10^6 \, {\rm K}$) Galactic corona as a consequence of thermal instabilities has been controversial. Here we re-examine this problem and we suggest that rotation of the corona might be a missing key ingredient. We do this by studying the evolution of the models of rotating galactic coronae presented in Sormani et al. (2018) under the presence of cooling and thermal conduction. We combine a linear stability analysis with the results of local and global hydrodynamical simulations. We find that condensations are likely to occur in regions where the corona has substantial rotational support. Under reasonably general assumptions on the rotation profile of the corona, the locations where condensations are expected are in remarkable agreement with the observed location of the major non-magellanic HVCs complexes in our Galaxy (namely, at distances $< 15 \, \rm kpc$ from the Sun and within $30^\circ$ from the disc plane). We conclude that HVCs can form by thermal instabilities provided that (i) the corona is rotating substantially in the inner ($R < 50 \, \rm kpc$) parts, as suggested by current observational data and predicted by cosmological simulations of galaxy formation; (ii) close to the disc the corona is well-represented by a nearly-equilibrium stratified rotating structure (as opposed to a fast cooling flow). Our results also suggest that a better understanding of the disc-halo interface, including supernova feedback, is critical to understand the origin of HVCs.

It is well known that the Galactic bar drives a gas inflow into the Central Molecular Zone, which fuels star formation, accretion onto the central super-massive black hole, and large-scale outflows. This inflow happens mostly through two symmetrical dust-lanes, similar to those often seen in external barred galaxies. Here we use the fact that the Milky Way dust-lanes have been previously identified in $^{12}$CO datacubes and a simple geometrical model to derive the first observational determination of the mass inflow rate into the Central Molecular Zone. We find that the time-averaged inflow rate along the near-side dust lane is $1.2^{+0.7}_{-0.8}\, \rm M_\odot yr^{-1}$ and along the far-side dust lane is $1.5^{+0.9}_{-1.0}\, \rm M_\odot yr^{-1}$, which gives a total inflow of $2.7^{+1.5}_{-1.7}\, \rm M_\odot yr^{-1}$. We also provide the time series of the inflow rate $\dot{M}$ for the future few Myr. The latter shows that the inflow rate is variable with time, supporting a scenario of episodic accretion onto the Central Molecular Zone.

Fitting equilibrium dynamical models to observational data is an essential step in understanding the structure of the gaseous hot haloes that surround our own and other galaxies. However, the two main categories of models that are used in the literature are poorly suited for this task: (i) simple barotropic models are analytic and can therefore be adjusted to match the observations, but are clearly unrealistic because the rotational velocity $v_\phi(R,z)$ does not depend on the distance $z$ from the galactic plane, while (ii) models obtained as a result of cosmological galaxy formation simulations are more realistic, but are impractical to fit to observations due to high computational cost. Here we bridge this gap by presenting a general method to construct axisymmetric baroclinic equilibrium models of rotating galactic coronae in arbitrary external potentials. We consider in particular a family of models whose equipressure surfaces in the $(R,z)$ plane are ellipses of varying axis ratio. These models are defined by two one-dimensional functions, the axial ratio of pressure $q_{\rm axis}(z)$ and the value of the pressure $P_{\rm axis}(z)$ along the galaxy's symmetry axis. These models can have a rotation speed $v_\phi(R,z)$ that realistically decreases as one moves away from the galactic plane, and can reproduce the angular momentum distribution found in cosmological simulations. The models are computationally cheap to construct and can thus be used in fitting algorithms. We provide a python code that given $q_{\rm axis}(z)$, $P_{\rm axis}(z)$ and $\Phi(R,z)$ returns $\rho(R,z)$, $T(R,z)$, $P(R,z)$, $v_\phi(R,z)$. We show a few examples of these models using the Milky Way as a case study.

We develop a dynamical theory for the origin of nuclear rings in barred galaxies. In analogy with the standard theory of accretion discs, our theory is based on shear viscous forces among nested annuli of gas. However, the fact that gas follows non circular orbits in an external barred potential has profound consequences: it creates a region of reverse shear in which it is energetically favourable to form a stable ring which does not spread despite dissipation. Our theory allows us to approximately predict the size of the ring given the underlying gravitational potential. The size of the ring is loosely related to the location of the Inner Lindblad Resonance in the epicyclic approximation, but the predicted location is more accurate and is also valid for strongly barred potentials. By comparing analytical predictions with the results of hydrodynamical simulations, we find that our theory provides a viable mechanism for ring formation if the effective sound speed of the gas is low ($\cs\lesssim1\kms$), but that nuclear spirals/shocks created by pressure destroy the ring when the sound speed is high ($\cs\simeq10\kms$). We conclude that whether this mechanism for ring formation is relevant for real galaxies ultimately depends on the effective equation of state of the ISM. Promising confirmation comes from simulations in which the ISM is modelled using state-of-the-art cooling functions coupled to live chemical networks, but more tests are needed regarding the role of turbulence driven by stellar feedback. If the mechanism is relevant in real galaxies, it could provide a powerful tool to constrain the gravitational potential, in particular the bar pattern speed.

We report the discovery from our IRSF/SIRIUS Near-Infrared (NIR) variability survey of five new classical Cepheids located in the inner Galactic Disk, at longitude $l\simeq -40^{\circ}$. The new Cepheids are unique in probing the kinematics and metallicity of young stars at the transition between the inner Disk and the minor axis of the central Bar, where they are expected to be less affected by its dynamical influence.This is also the first time that metallicity of Cepheids is estimated on the basis of medium-resolution ($R\sim3,000$) NIR spectra, and we validated our results with data in the literature, finding a minimal dependence on the adopted spectroscopic diagnostics. This result is very promising for using Cepheids as stellar proxy of the present-time chemical content of the obscured regions in the Disk. We found that the three Cepheids within 8--10 kpc from us have metallicities consistent with the mean radial metallicity gradient, and kinematics consistent with the Galactic rotation curve. Instead, the closest ($\sim$4 kpc)/farthest ($\sim$12 kpc) Cepheids have significant negative/positive residuals, both in velocity and in iron content. We discuss the possibility that such residuals are related to large-scale dynamical instabilities, induced by the bar/spiral-arm pattern, but the current sample is too limited to reach firm conclusions.

It has been known for more than thirty years that the distribution of molecular gas in the innermost 300 parsecs of the Milky Way, the Central Molecular Zone, is strongly asymmetric. Indeed, approximately three quarters of molecular emission comes from positive longitudes, and only one quarter from negative longitudes. However, despite much theoretical effort, the origin of this asymmetry has remained a mystery. Here we show that the asymmetry can be neatly explained by unsteady flow of gas in a barred potential. We use high-resolution 3D hydrodynamical simulations coupled to a state-of-the-art chemical network. Despite the initial conditions and the bar potential being point-symmetric with respect to the Galactic Centre, asymmetries develop spontaneously due to the combination of a hydrodynamical instability known as the "wiggle instability" and the thermal instability. The observed asymmetry must be transient: observations made tens of megayears in the past or in the future would often show an asymmetry in the opposite sense. Fluctuations of amplitude comparable to the observed asymmetry occur for a large fraction of the time in our simulations, and suggest that the present is not an exceptional moment in the life of our Galaxy.

We run hydrodynamical simulations of a 2D isothermal non self-gravitating inviscid gas flowing in a rigidly rotating externally imposed potential formed by only two components: a monopole and a quadrupole. We explore systematically the effects of varying the quadrupole while keeping fixed the monopole and discuss the consequences for the interpretation of longitude-velocity diagrams in the Milky Way. We find that the gas flow can constrain the quadrupole of the potential and the characteristics of the bar that generates it. The exponential scale length of the bar must be at least $1.5\rm\, kpc$. The strength of the bar is also constrained. Our global interpretation favours a pattern speed of $\Omega=40\,\rm km s^{-1} {kpc}^{-1}$. We find that for most observational features, there exist a value of the parameters that matches each individual feature well, but is difficult to reproduce all the important features at once. Due to the intractably high number of parameters involved in the general problem, quantitative fitting methods that can run automatic searches in parameter space are necessary.

Long-duration, spectrally-soft Gamma-Ray Bursts (GRBs) are associated with Type Ic Core Collapse (CC) Supernovae (SNe), and thus arise from the death of massive stars. In the collapsar model, the jet launched by the central engine must bore its way out of the progenitor star before it can produce a GRB. Most of these jets do not break out, and are instead "choked" inside the star, as the central-engine activity time, $t_{\rm e}$, is not long enough. Modelling the long-soft GRB duration distribution assuming a power-law distribution for their central-engine activity times, $\propto t_{\rm e}^{-\alpha}$ for $t_{\rm e}>t_{\rm b}$, we find a steep distribution ($\alpha\sim4$) and a typical GRB jet breakout time of $t_{\rm b}\sim 60\text{ s}$ in the star's frame. The latter suggests the presence of a low-density, extended envelope surrounding the progenitor star, similar to that previously inferred for low-luminosity GRBs. Extrapolating the range of validity of this power law below what is directly observable, to $t_{\rm e}<t_{\rm b}$, by only a factor of $\sim$4-5 produces enough events to account for all Type Ib/c SNe. Such extrapolation is necessary to avoid fine-tuning the distribution of central engine activity times with the breakout time, which are presumably unrelated. We speculate that central engines launching relativistic jets may operate in all Type Ib/c SNe. In this case, the existence of a common central engine would imply that (i) the jet may significantly contribute to the energy of the SN; (ii) various observational signatures, like the asphericity of the explosion, could be directly related to jet's interaction with the star.

We study the dynamical stability of stationary galactic spiral shocks. The steady-state equilibrium flow contains a shock of the type derived by Roberts in the tightly wound approximation. We find that boundary conditions are critical in determining whether the solutions are stable or not. Shocks are unstable if periodic boundary conditions are imposed. For intermediate strengths of the spiral potential, the instability disappears if boundary conditions are imposed such that the upstream flow is left unperturbed as in the classic analysis of D'yakov and Kontorovich. This reconciles apparently contradictory findings of previous authors regarding the stability of spiral shocks. This also shows that the instability is distinct from the Kelvin-Helmholtz instability, confirming the findings of Kim et al. We suggest that instability is a general characteristics of periodic shocks, regardless of the presence of shear, and provide a physical picture as to why this is the case. For strong spiral potentials, high post-shock shear makes the system unstable also to parasitic Kelvin-Helmholtz instability regardless of the boundary conditions. Our analysis is performed in the context of a simplified problem that, while preserving all the important characteristics of the original problem, strips it from unnecessary complications, and assumes that the gas is isothermal, non self-gravitating, non-magnetised.

We use hydrodynamical simulations to construct a new coherent picture for the gas flow in the Central Molecular Zone (CMZ), the region of our Galaxy within $R\leq 500\, \mathrm{pc}$. We relate connected structures observed in $(l,b,v)$ data cubes of molecular tracers to nuclear spiral arms. These arise naturally in hydrodynamical simulations of barred galaxies, and are similar to those that can be seen in external galaxies such as NGC4303 or NGC1097. We discuss a face-on view of the CMZ including the position of several prominent molecular clouds, such as Sgr B2, the $20\,{\rm km\, s^{-1}}$ and $50\,{\rm km\, s^{-1}}$ clouds, the polar arc, Bania Clump 2 and Sgr C. Our model is also consistent with the larger scale gas flow, up to $R\simeq 3\,\rm kpc$, thus providing a consistent picture of the entire Galactic bar region.

In the paradigm of magnetic acceleration of relativistic jets, one of the key points is identifying a viable mechanism to convert the Poynting flux into the kinetic energy of the plasma beyond equipartition. A promising candidate is the kink instability, which deforms the body of the jet through helical perturbations. Since the detailed structure of real jets is unknown, we explore a large family of cylindrical, force-free equilibria to get robust conclusions. We find that the growth rate of the instability depends primarily on two parameters: (i) the gradient of the poloidal magnetic field; (ii) the Lorentz factor of the perturbation, which is closely related to the velocity of the plasma. We provide a simple fitting formula for the growth rate of the instability. As a tentative application, we use our results to interpret the dynamics of the jet in the nearby active galaxy M87. We show that the kink instability becomes non-linear at a distance from the central black hole comparable to where the jet stops accelerating. Hence (at least for this object), the kink instability of the jet is a good candidate to drive the transition from a Poynting-dominated to a kinetic-energy-dominated flow.

Hydrodynamical simulations of star formation often do not possess the dynamic range needed to fully resolve the build-up of individual stars and star clusters, and thus have to resort to subgrid models. A popular way to do this is by introducing Lagrangian sink particles, which replace contracting high density regions at the point where the resolution limit is reached. A common problem then is how to assign fundamental stellar properties to sink particles, such as the distribution of stellar masses. We present a new and simple statistical method to assign stellar contents to sink particles. Once the stellar content is specified, it can be used to determine a sink particle's radiative output, supernovae rate or other feedback parameters that may be required in the calculations. Advantages of our method are (i) it is simple to implement, (ii) it guarantees that the obtained stellar populations are good samples of the initial mass function, (iii) it can easily deal with infalling mass accreted at later times, and (iv) it does not put restrictions on the sink particles' masses in order to be used. The method works very well for sink particles that represent large star clusters and for which the stellar mass function is well sampled, but can also handle the transition to sink particles that represent a small number of stars.

We describe a scenario to explain blazar periodicities with timescales of $\sim$ few years. The scenario is based on a binary super-massive black hole (SMBH) system in which one of the two SMBH carries a jet. We discuss the various mechanisms that can cause the jet to precess and produce corkscrew patterns through space with a scale of $\sim$ few pc. It turns out that the dominant mechanism responsible for the precession is simply the imprint of the jet-carrying SMBH orbital speed on the jet. Gravitational deflection and Lense-Thirring precession (due to the gravitational field of the other SMBH) are second order effects. We complement the scenario with a kinematical jet model which is inspired to the spine-sheath structure observed in M87. One of the main advantages of such a structure is that it allows the peak of the synchrotron emission to scale with frequency according to $\nu F\propto \nu^{\xi}$ as the viewing angle is changed, where $\xi$ is not necessarily 3 or 4 as in the case of jets with uniform velocity, but can be $\xi \sim 1$. Finally, we apply the model to the source PG1553+113, which has been recently claimed to show a $T_{\rm obs}=(2.18\pm 0.08)\text{ yr}$ periodicity. We are able to reproduce the optical and gamma-ray light curves and multiple synchrotron spectra simultaneously. We also give estimates of the source mass and size.

Spiral arms that emerge from the ends of a galactic bar are important in interpreting observations of our and external galaxies. It is therefore important to understand the physical mechanism that causes them. We find that these spiral arms can be understood as kinematic density waves generated by librations around underlying ballistic closed orbits. This is even true in the case of a strong bar, provided the librations are around the appropriate closed orbits and not around the circular orbits that form the basis of the epicycle approximation. An important consequence is that it is a potential's orbital structure that determines whether a bar should be classified as weak or strong, and not crude estimates of the potential's deviation from axisymmetry.

We use a Cartesian grid to simulate the flow of gas in a barred Galactic potential and investigate the effects of varying the sound speed in the gas and the resolution of the grid. For all sound speeds and resolutions, streamlines closely follow closed orbits at large and small radii. At intermediate radii shocks arise and the streamlines shift between two families of closed orbits. The point at which the shocks appear and the streamlines shift between orbit families depends strongly on sound speed and resolution. For sufficiently large values of these two parameters, the transfer happens at the cusped orbit as hypothesised by Binney et al. over two decades ago. For sufficiently high resolutions the flow downstream of the shocks becomes unsteady. If this unsteadiness is physical, as appears to be the case, it provides a promising explanation for the asymmetry in the observed distribution of CO.

We present a new method for fitting simple hydrodynamical models to the (l,v) distribution of atomic and molecular gas observed in the Milky Way. The method works by matching features found in models and observations. It is based on the assumption that the large-scale features seen in (l,v) plots, such as ridgelines and the terminal velocity curve, are influenced primarily by the underlying large-scale Galactic potential and are only weakly dependent on local ISM heating and cooling processes. In our scheme one first identifies by hand the features in the observations: this only has to be done once. We describe a procedure for automatically extracting similar features from simple hydrodynamical models and quantifying the "distance" between each model's features and the observations. Application to models of the Galactic Bar region (|l|<30deg) shows that our feature-fitting method performs better than \chi^2 or envelope distances at identifying the correct underlying galaxy model.

A re-investigation of the gravothermal catastrophe is presented. By means of a linear perturbation analysis, we study the dynamical stability of a spherical self-gravitating isothermal fluid of finite volume and find that the conditions for the onset of the gravothermal catastrophe, under different external conditions, coincide with those obtained from thermodynamical arguments. This suggests that the gravothermal catastrophe may reduce to Jeans instability, rediscovered in an inhomogeneous framework. We find normal modes and frequencies for the fluid system and show that instability develops on the dynamical time scale. We then discuss several related issues. In particular: (1) For perturbations at constant total energy and constant volume, we introduce a simple heuristic term in the energy budget to mimic the role of binaries. (2) We outline the analysis of the two-component case and show how linear perturbation analysis can be carried out also in this more complex context in a relatively straightforward way. (3) We compare the behavior of the fluid model with that of the collisionless sphere. In the collisionless case the instability seems to disappear, which is at variance with the linear Jeans stability analysis in the homogeneous case; we argue that a key ingredient to understand the difference (a spherical stellar system is expected to undergo the gravothermal catastrophe only in the presence of some collisionality, which suggests that the instability is dissipative and not dynamical) lies in the role of the detailed angular momentum in a collisionless system. Finally, we briefly comment on the meaning of the Boltzmann entropy and its applicability to the study of the dynamics of self- gravitating inhomogeneous gaseous systems.