A complete map of the 3D distribution of molecular (CO) gas was constructed using a realistic dynamical model of the gas flow in the barred potential of the Milky Way. The map shows two prominent spiral arms starting at the bar ends connecting smoothly to the 4-armed spiral pattern observed in the atomic hydrogen gas in the outer Galaxy. Unlike previous attempts, our new map uncovers the gas distribution in the bar region of the Galaxy and the far side of the disk. For the first time, we can follow spiral arms in gas as they pass behind the galactic centre.

We present a new model of the three-dimensional distribution of molecular gas in the Milky Way Galaxy, based on CO line data. Our analysis is based on a gas-flow simulation of the inner Galaxy using smoothed-particle hydrodynamics (SPH) using a realistic barred gravitional potential derived from the observed COBE/DIRBE near-IR light distribution. The gas model prescribes the gas orbits much better than a simple circular rotation model and is highly constrained by observations, but it cannot predict local details. In this study, we provide a 3D map of the observed molecular gas distribution using the velocity field from the SPH model. A comparison with studies of the Galactic Center region suggests that the main structures are reproduced but somewhat stretched along the line-of-sight, probably on account of limited resolution of the underlying SPH simulation. The gas model will be publicly available and may prove useful in a number of applications, among them the analysis of diffuse gamma-ray emission as measured with GLAST.

We present a model of the global magnetic field in the barred galaxy NGC 1365 based jointly on the large-scale velocity field of interstellar gas fitted to HI and CO observations of this galaxy and on mean-field dynamo theory. The aim of the paper is to present a detailed quantitative comparison of a galactic dynamo model with independent radio observations. We consider several gas dynamical and nonlinear dynamo models that include plausible variations of parameters that are poorly known. Models of cosmic ray distribution in the galaxy are introduced to produce synthetic radio polarization maps allowing direct comparison with those observed at 3.5cm and 6.2cm. We show that the dynamo model is robust in that the most important magnetic features are controlled by the relatively well established properties of the density distribution and gas velocity field. The optimal agreement between the synthetic polarization maps and observations is obtained when a uniform cosmic ray distribution is adopted. We find some indirect evidence for enhanced turbulence in the regions of strong velocity shear and within 1-2kpc of the galactic centre. We confirm that magnetic stresses can drive an inflow of gas into the inner 1kpc of the galaxy at a rate of a few solar masses per year. The dynamo models are successful to some extent in modelling the large scale regular magnetic field in this galaxy. Our results demonstrate that dynamo models and synthetic polarization maps can provide information about both the gas dynamical models and conditions in the interstellar medium. We demonstrate that the dynamical effects of magnetic fields cannot be everywhere ignored in galaxy modelling.

The Milky Way is made up of a central bar, a disk with embedded spiral arms, and a dark matter halo. Observational and theoretical constraints for the characteristic parameters of these components will be presented, with emphasis on the constraints from the dynamics of the Milky Way gas. In particular, the fraction of dark matter inside the solar radius, the location of the main resonances, and the evidence for multiple pattern speeds will be discussed.

We present the first interferometric observations of CO(1-0) and CO(2-1) line emission from the warped LINER NGC3718, obtained with the IRAM PdBI. This L1.9 galaxy has a prominent dust lane and on kpc scales, a strongly warped atomic and molecular gas disk. The molecular gas is closely associated with the dust lane across the nucleus and its kinematic center is consistent with the mm continuum AGN. A comparison of our interferometric mosaic data, which fully cover the ~9kpc warped disk, with a previously obtained IRAM 30m single dish CO(1-0) map shows that the molecular gas distribution in the disk is heavily resolved by the PdBI map. After applying a short-spacing correction with the IRAM 30m data, we find in total six main source components within the dust lane: one associated with the nucleus, four symmetrically positioned on either side at galactocentric distances of about 1.3kpc and 4.0kpc from the center, and a sixth on the western side at ~3kpc with only a very weak eastern counterpart. In the framework of a kinematic model using tilted rings, we interpret the five symmetric source components as locations of strong orbital crowding. We further find indications that the warp appears not only on kpc scales, but continues down to 250pc. Besides the sixth feature on the western side, the lower flux of the eastern components compared to the western ones indicates an intrinsic large scale asymmetry in NGC3718 that cannot be explained by the warp. Indications for a small scale asymmetry are also seen in the central 600pc. These asymmetries might be evidence for a tidal interaction with a companion galaxy (large scales) and gas accretion onto the nucleus (small scales). Our study of NGC3718 is part of the NUGA project that aims at investigating the different processes of gas accretion onto AGN.

A substantial fraction of barred galaxies host additional nuclear bars which tumble with pattern speeds exceeding those of the large-scale (primary) stellar bars. We have investigated the mechanism of formation and dynamical decoupling in such nested bars which include gaseous (secondary) nuclear bars within the full size galactic disks, hosting a double inner Lindblad resonance. Becoming increasingly massive and self-gravitating, the nuclear bars lose internal (circulation) angular momentum to the primary bars and increase their strength. Developing chaos within these bars triggers a rapid gas collapse -- bar contraction. During this time period, the secondary bar pattern speed Omega_s~a^{-1}, where "a" stands for the bar size. As a result, Omega_s increases dramatically until a new equilibrium is reached (if at all), while the gas specific angular momentum decreases -- demonstrating the dynamical decoupling of nested bars. Viscosity, and therefore the gas presence, appears to be a necessary condition for the prograde decoupling of nested bars. This process maintains an inflow rate of ~1 M_o/yr over ~10^8 yrs across the central 200 pc and has important implications for fueling the nuclear starbursts and AGN.

The NUGA project is a high-resolution (0.5''-1'') CO survey of low luminosity AGN including the full sequence of activity types (Seyferts, LINERs and transition objects). NUGA aims to systematically study the different mechanisms for gas fueling of AGNs in the Local Universe. In this paper we discuss the latest results of this recently completed survey, which now includes newly acquired subarcsec resolution observations for all targets of the sample. The large variety of circumnuclear disk morphologies found in NUGA galaxies (m=1, m=2 and stochastic instabilities) is a challenging result that urges the refinement of current dynamical models. In this paper we report on new results obtained in 4 study cases for NUGA: NGC4826, NGC7217, NGC4579 and NGC6951

The NUGA project is the first high-resolution (0.5''-1'') CO survey of low luminosity AGN including the full sequence of activity types (Seyfert 1--2 and LINERs). NUGA aims to study systematically the different mechanisms for gas fueling of the AGN. In order to reach the critical spatial scales for AGN feeding (<100 pc), we are undertaking high-resolution observations, made at the Plateau de Bure Interferometer, to achieve a sharp view of the distribution and kinematics of molecular gas. In this talk we discuss the latest results of this survey which now includes newly acquired subarcsec resolution observations for 5 targets of our sample. The large variety of circumnuclear disk morphologies found in NUGA galaxies (m=1, m=2 and m=0-type gravitational instabilities) is a challenging result that urges the refinement of current dynamical models.

We present CO(1-0) and CO(2-1) maps of the LINER galaxy NGC7217, obtained with the IRAM interferometer, at 2.4"x1.9" and 1.2"x0.8" resolution respectively. The nuclear ring (at r=12"=0.8kpc) dominates the CO maps, and has a remarkable sharp surface density gradient at its inner edge. The latter is the site of the stellar/Halpha ring, while the CO emission ring extends farther or is broader (500-600pc). This means that the star formation has been more intense toward the inner edge of the CO ring, in a thin layer, just at the location of the high gas density gradient. In the CO(2-1) map, a central concentration inside 3" coincides with a spiral dust lane, clearly seen in the HST V-I color image. N-body simulations including gas dissipation and star formation are performed to better understand the nature of the nuclear ring observed. In the best model, the gas content was higher in the recent past (having been consumed via star formation), and the oval/spiral formed were more self-gravitating (abridged).

We present new gas flow models for the Milky Way inside the solar circle. Using a mass model derived from the COBE/DIRBE maps and clump giant star counts, and using a parametric model for the spiral arm pattern in the disk, we calculate the gas flow and compare with 12CO observations. We find that models with 4 spiral arms fit the observations better then 2-armed models. We also find that models with separate pattern speeds for the bar and spiral arms can explain the gas flow in the bar corotation region better than single-pattern speed models.

We present new high-resolution observations of the nucleus of the counter-rotating LINER NGC4826, made in the J=1-0 and J=2-1 lines of 12CO with the IRAM Plateau de Bure mm-interferometer(PdBI).The CO maps, which achieve 0.8''(16pc) resolution in the 2-1 line, fully resolve an inner molecular gas disk which is truncated at an outer radius of 700pc. The total molecular gas mass is distributed in a lopsided nuclear disk of 40pc radius and two one-arm spirals, which develop at different radii in the disk. The distribution and kinematics of molecular gas in the inner 1kpc of NGC4826 show the prevalence of different types of m=1 perturbations in the gas. Although dominated by rotation, the gas kinematics are perturbed by streaming motions related to the m=1 instabilities. The non-circular motions associated with the inner m=1 perturbations agree qualitatively with the pattern expected for a trailing wave developed outside corotation ('fast' wave). In contrast, the streaming motions in the outer m=1 spiral are better explained by a 'slow' wave. A paradoxical consequence is that the inner m=1 perturbations would not favour AGN feeding. An independent confirmation that the AGN is not being generously fueled at present is found in the low values of the gravitational torques exerted by the stellar potential for R<530pc. The distribution of star formation in the disk of NGC4826 is also strongly asymmetrical. Massive star formation is still vigorous, fed by the significant molecular gas reservoir at R<700pc. There is supporting evidence for a recent large mass inflow episode in NGC4826. These observations have been made in the context of the NUclei of GAlaxies (NUGA) project, aimed at the study of the different mechanisms for gas fueling of AGN.

We present new gas flow models for the Milky Way inside the solar circle. To this end we use SPH simulations in gravitational potentials determined from the NIR luminosity distribution (including spiral arms) which are based on the COBE/DIRBE maps. Gas flows in models which include massive spiral arms clearly match the observed 12CO lvplot better than if the potential does not include spiral structure. Besides single pattern speed models we investigate models with separate pattern speeds for the bar and spiral arms. The most important difference is that in the latter case the gas spiral arms go through the bar corotation region, keeping the gas aligned with the arms there. In the (l,v) plot this results in characteristic regions which appear to be nearly void of gas.In single pattern speed models these regions are filled with gas because the spiral arms dissolve in the bar corotation region. Comparing with the 12CO data we find evidence for separate pattern speeds in the Milky Way.From a series of models the preferred range for the bar pattern speed is Om_p=60\pm5 /Gyr, corresponding to corotation at 3.4\pm0.3kpc. The spiral pattern speed is less well constrained, but our preferred value is Om_sp\approx 20 /Gyr. A further series of gas models is computed for different bar angles, using separately determined luminosity models and gravitational potentials in each case. We find acceptable gas models for 20<=\phibar<=25. The model with (\phibar=20, Om_p=60 /Gyr, Om_sp=20 /Gyr) gives an excellent fit to the spiral arm ridges in the observed (l,v) plot.

The NUclei of GAlaxies (NUGA) project is a combined effort to carry out a high-resolution (<1'') interferometer CO survey of a sample of 12 nearby AGN spiral hosts, using the IRAM array. We map the distribution and dynamics of molecular gas in the inner 1 kpc of the nuclei with resolutions of 10-50 pc, and study the mechanisms for gas fueling of the different low-luminosity AGN. First results show evidence for the occurrence of strong m=1 gas instabilities in Seyferts. NUGA maps allow us to address the origin/nature of m=1 modes and their link with m=2 modes and acoustic instabilities, present in other targets.

We present CO(1-0) and CO(2-1) maps of the Seyfert galaxy NGC 7217, obtained with the IRAM interferometer, at 3" and 1.5" resolution respectively. The nuclear ring (at r=12"=0.8kpc) is predominant in the CO maps, with a remarkable surface density gradient between the depleted region inside the ring and the inner border of the ring. The CO nuclear ring is significantly broader (500-600pc) than the dust lane ring. The CO(2-1)/CO(1-0) ratio is around 1, typical of optically thick gas with high density. The overall morphology of the ring is quite circular, with no evidence of non-circular velocities. In the CO(2-1) map, a central concentration might be associated with the circumnuclear ring of ionised gas detected inside r=3" and interpreted as a polar ring by Sil'chenko and Afanasiev (2000). Our interpretation is more in terms of a bar/spiral structure, in the same plane as the global galaxy but affected by non-circular motions, which results in a characteristic S-shape of the isovels. This nuclear bar/spiral structure, clearly seen in a V-I HST colour image, is essentially gaseous and might be explained with acoustic waves.

Through observations and modeling, we demonstrate how the recently discovered large-scale bar in NGC 5248 generates spiral structure which extends from 10 kpc down to 100 pc, fuels star formation on progressively smaller scales, and drives disk evolution. Deep inside the bar, two massive molecular spirals cover nearly 180 degrees in azimuth, show streaming motions of 20--40 km/s and feed a starburst ring of super star clusters at 375 pc. They also connect to two narrow K-band spirals which delineate the UV-bright star clusters in the ring. The data suggest that the K-band spirals are young, and the starburst has been triggered by a bar-driven spiral density wave (SDW). The latter may even have propagated closer to the center where a second Halpha ring and a dust spiral are found. The molecular and HST data support a scenario where stellar winds and supernovae efficiently clear out gas from dense star-forming regions on timescales less than a few Myrs. We have investigated the properties of massive CO spirals within the framework of bar-driven SDWs, incorporating the effect of gas self-gravity. We find good agreement between the model predictions and the observed morphology, kinematics, and pitch angle of the spirals. This combination of observations and modeling provides the best evidence to date for a strong dynamical coupling between the nuclear region and the surrounding disk. It also confirms that a low central mass concentration, which may be common in late-type galaxies, is particularly favorable to the propagation of a bar-driven gaseous SDW deep into the central region of the galaxy, whereas a large central mass concentration favors other processes, such as the formation and decoupling of nuclear bars.

For decades, the grand-design SAB spiral galaxy NGC 5248 has been postulated to host a short bar of semi-major axis 22" (1.6 kpc). From dynamical and morphological arguments, we argue, however, that its spiral structure is being driven by a large-scale bar whose corotation radius lies at around 115" (8.6 kpc).Our estimate is based partially on a deep R-band image which reveals that the feature previously thought to be an inclined disk is in fact an extended stellar bar. The bar is embedded within a fainter outer disk visible out to a radius of 230" (17.2 kpc). The bar has a deprojected ellipticity of 0.44 and a semi-major axis of 95" (7.1 kpc). The classical grand-design spirals of NGC 5248, prominent in B, R, and K light, lie on the leading edge of the large-scale stellar bar and are accompanied by concave dust lanes out to at least 70". The offset between the dust and young stars is consistent with our understanding of gas flows in barred galaxies, where shocks along the leading edges of a moderately strong bar compress the gas to form massive young stars. While in many strongly barred galaxies, optical spiral arms are prominent outside the bar but not within it, NGC 5248 illustrates how intense star formation along a moderately strong bar can lead to conspicuous open spiral arms within the bar itself. NGC 5248 also provides a clear example of how a large-scale stellar bar embedded within a faint outer optical disk can be misidentified as an inclined disk when imaging studies lack the sensitivity to detect the actual outer disk. We discuss the implications for the estimated bar fraction at higher redshifts.

We find that nuclear rings in barred galaxies can be subject to a new type of non-self-gravitational dynamical instability. The instability leads to the formation of gaseous molecular bars with pattern speeds which are substantially slower than speeds of the primary stellar bars. This spectacular decoupling of nuclear bars from the underlying gravitational potential is triggered but is not driven by the gas viscosity. We find that low-viscosity systems can spend a substantial period of time in a fully decoupled state, with the nuclear bar slowly tumbling in the gravitational field of the primary bar. Higher viscosity systems form nuclear bars which librate about the primary bar. The shape of a nuclear bar, i.e., its eccentricity, correlates strongly with the angle between the bars. We also find that such decoupling, partial or full, most probably will be associated with bursts of star formation and with gas inflow across the inner (Lindblad) resonance zone towards smaller radii.

We present for the first time a two-dimensional velocity field of the central region of the grand-design spiral galaxy NGC 5248, at 0.9 arcsec spatial resolution. The H-alpha velocity field is dominated by circular rotation. While no systematic streaming motions are seen in the area of the nuclear grand-design spiral or the circumnuclear ring, the amplitude of residual velocities, after subtracting a model circular velocity field, reaches 20 km/s in projection. The rotation curve levels out at around 140 km/s, after a well-resolved and rather shallow rise. We have generated an analytical model for the nuclear spiral and fitted it to our observations to obtain estimates of the pattern speed of the spiral and the sound speed in the central region of NGC 5248. Our results are consistent with a low pattern speed, suggesting that the nuclear spiral rotates with the same rate as the main spiral structure in NGC 5248, and thus that the spiral structure is coupled from scales of a few hundred parsecs to several kiloparsecs. We have also compared the observed structure and kinematics between the nuclear regions of NGC 5248 and M100. Several similarities and differences are discussed, including the location of the peak emission regions on major and minor axes, and the spiral arm streaming motions. We find no kinematic evidence for a presence of a nuclear bar in NGC 5248.

The prominent interacting galaxy pair NGC 4038/9 contains many active star-forming regions and is continuously forming new star clusters. We present a self-consistent n-body model for this system which includes an SPH gas component. The model qualitatively explains the apparent concentration of gas in the so-called overlap region between the two nuclei as a bridge of gas connecting the two galaxies. Projected on the sky, the bridge appears as a dense spot of gas. We discuss some implications for the evolution of Ultra-luminous infrared galaxies.

The observed gas dynamics in the Milky Way can only be explained by a bar in the galactic center. Such a bar is directly visible in the near-IR maps of the bulge, where it causes a distinctive asymmetric light distribution pattern. Another large-scale structure is the grand-design 4-arm spiral pattern, most clearly observed in the spatial distribution of molecular gas and HII-regions. In order to model the observed gas flow structure, we constructed a model for the stellar mass distribution. For the inner 5 $\kpc$ we used the 3D deprojected near-IR light distribution, as observed by the COBE/DIRBE experiment, and added an analytical disk model outside the box as well as a halo model. With this frozen mass distribution, we computed the stationary gas flow for various deprojection parameters and pattern speeds. For all reasonable parameter choices, we obtain a 4-armed spiral pattern, which can be matched to the observed spiral arms. In the bar region, our model can explain the non-circular motion visible in the terminal velocity curve as well as part of the forbidden velocities. Inside the corotation, we also find 4 spiral arms, the nearest arm corresponds to the 3-kpc-arm, although only qualitatively. The missing southern 3-kpc-arm at the far end of the galaxy is explained by running parallel to another arm. Close to the center, we find gas on circular orbits forming a disk. Such a disk has been observed in emission of the CS molecule, however only part of the disk appears to be occupied by dense enough gas to be traced by CS. Further we compare our model to the distribution of OH/IR stars in the inner galaxy.

We analyze formation of grand-design two-arm spiral structure in the nuclear regions of disk galaxies. Such morphology has been recently detected in a number of objects using high-resolution near-infrared observations. Motivated by the observed (1) continuity between the nuclear and kpc-scale spiral structures, and by (2) low arm-interarm contrast, we apply the density wave theory to explain the basic properties of the spiral nuclear morphology. In particular, we address the mechanism for the formation, maintenance and the detailed shape of nuclear spirals. We find, that the latter depends mostly on the shape of the underlying gravitational potential and the sound speed in the gas. Detection of nuclear spiral arms provides diagnostics of mass distribution within the central kpc of disk galaxies. Our results are supported by 2D numerical simulations of gas response to the background gravitational potential of a barred stellar disk. We investigate the parameter space allowed for the formation of nuclear spirals using a new method for constructing a gravitational potential in a barred galaxy, where positions of resonances are prescribed.

We report simulations of the gas flow in the gravitational potential of the COBE NIR bulge and disk. These models lead to four--armed spiral structure between corotation of the bar and the Sun, in agreement with the observed spiral arm tangents. The 3-kpc-arm is identified with one of the arms emanating from the ends of the bar.

We present a new model for the gas dynamics in the galactic disk inside the Sun's orbit. Quasi-equilibrium flow solutions are determined in the gravitational potential of the deprojected COBE NIR bar and disk, complemented by a central cusp and, in some models, an outer halo. These models generically lead to four-armed spiral structure between corotation of the bar and the solar circle; their large-scale morphology is not sensitive to the precise value of the bar's pattern speed, to the orientation of the bar with respect to the observer, and to whether or not the spiral arms carry mass. Our best model provides a coherent interpretation of many observed gas dynamical features. Its four-armed spiral structure outside corotation reproduces quantitatively the directions to the five main spiral arm tangents at |l|<=60deg observed in a variety of tracers. The 3-kpc-arm is identified with one of the model arms emanating from the ends of the bar, extending into the corotation region. The model features an inner gas disk with a cusped orbit shock transition to an x_2 orbit disk of radius R~150pc. The bar's corotation radius is fairly well--constrained at R_c=3.5 +/- 0.5 kpc. The best value for the orientation angle of the bar is probably 20-25deg, but the uncertainty is large since no detailed quantitative fit to all features in the observed lv-diagrams is yet possible. The Galactic terminal velocity curve from HI and CO observations out to l=+/-45deg (=5 kpc) is approximately described by a maximal disk model with constant mass-to-light ratio for the NIR bulge and disk.

We examine the left-right asymmetry in the cleaned COBE/DIRBE near-infrared data of the inner Galaxy and show (i) that the Galactic bar is probably not seen very nearly end-on, and (ii) that even if it is, it is not highly elongated. The assumption of constant mass-to-light ratio is used to derive simulated terminal-velocity plots for the ISM from our model luminosity distributions. By comparing these plots with observed terminal velocities we determine the mass-to-light ratio of the near-IR bulge and disk. Assuming that all this mass contributes to gravitational microlensing we compute optical depths $\tau$ for microlensing in Galactic-centre fields. For three models with bar major axis between $10\deg-25\deg$ from the Sun-Galactic Center line, the resulting optical depths in Baade's window lie in the range $0.83\times10^{-6} \lta \tau \lta 0.89\times10^{-6}$ for main-sequence stars and $1.2\times10^{-6} \lta \tau \lta 1.3\times10^{-6}$ for red-clump giants. We discuss a number of uncertainties including possible variations of the near-infrared mass-to-light ratio. We conclude that, although the values predicted from analyzing the COBE and gas velocity data are inconsistent at the $2-2.5\sigma$ level with recent observational determinations of $\tau$, we believe they should be taken seriously.

We investigate stationary gas flows in a fixed, rotating barred potential. The gas is assumed to be isothermal with an effective sound speed c_s, and the equations of motion are solved with smoothed particle hydrodynamics (SPH). Since the thermal energy in cloud random motions is negligible compared to the orbital kinetic energy, no dependence of the flow on c_s is expected. However, this is not the case when shocks are involved. For low values of c_s an open, off-axis shock flow forms that is characteristic for potentials with an inner Lindblad resonance (ILR). Through this shock the gas streams inwards from x_1 to x_2-orbits. At high sound speeds the gas arranges itself in a different, on-axis shock flow pattern. In this case, there is no gas on x_2-orbits, demonstrating that the gas can behave as if there were no ILR. The critical effective sound speed dividing the two regimes is in the range of values observed in the Milky Way. We give a heuristic explanation for this effect. A possible consequence is that star formation may change the structure of the flow by which it was initiated. Low-mass galaxies should predominantly be in the on-axis regime. A brief comparison of our SPH results with those from a grid-based hydrodynamic code is also given.