We report on the determination of electron densities, and their impact on the outflow masses and rates, measured in the central few hundred parsecs of 11 local luminous active galaxies. We show that the peak of the integrated line emission in the AGN is significantly offset from the systemic velocity as traced by the stellar absorption features, indicating that the profiles are dominated by outflow. In contrast, matched inactive galaxies are characterised by a systemic peak and weaker outflow wing. We present three independent estimates of the electron density in these AGN, discussing the merits of the different methods. The electron density derived from the [SII] doublet is significantly lower than than that found with a method developed in the last decade using auroral and transauroral lines, as well as a recently introduced method based on the ionization parameter. The reason is that, for gas photoionized by an AGN, much of the [SII] emission arises in an extended partially ionized zone where the implicit assumption that the electron density traces the hydrogen density is invalid. We propose ways to deal with this situation and we derive the associated outflow rates for ionized gas, which are in the range 0.001--0.5 M$_{\odot}$yr$^{-1}$ for our AGN sample. We compare these outflow rates to the relation between $\dot{M}_{out}$ and $L_{AGN}$ in the literature, and argue that it may need to be modified and rescaled towards lower mass outflow rates.

The $M_{BH}$ - $\sigma_{\star}$ relation is considered a result of co-evolution between the host galaxies and their super-massive black holes. For elliptical-bulge hosting inactive galaxies, this relation is well established, but there is still a debate whether active galaxies follow the same relation. In this paper, we estimate black hole masses for a sample of 19 local luminous AGNs (LLAMA) in order to test their location on the $M_{BH}$ - $\sigma_{\star}$ relation. Super-massive black hole masses ($M_{BH}$) were derived from the broad-line based relations for H$\alpha$, H$\beta$ and Pa$\beta$ emission line profiles for the Type 1 AGNs. We compare the bulge stellar velocity dispersion ($\sigma_{\star}$) as determined from the Ca II triplet (CaT) with the dispersion measured from the near-infrared CO (2-0) absorption features for each AGN and find them to be consistent with each other. We apply an extinction correction to the observed broad line fluxes and we correct the stellar velocity dispersion by an average rotation contribution as determined from spatially resolved stellar kinematic maps. The H$\alpha$-based black hole masses of our sample of AGNs were estimated in the range 6.34 $\leq$ $\log{M_{BH}}$ $\leq$ 7.75 M$_\odot$ and the $\sigma_{\star CaT}$ estimates range between 73 $\leq$ $\sigma_{\star CaT}$ $\leq$ 227 km s$^{-1}$. From the so-constructed $M_{BH}$ - $\sigma_{\star}$ relation for our Type 1 AGNs, we estimate the black hole masses for the Type 2 AGNs and the inactive galaxies in our sample. In conclusion, we find that our sample of local luminous AGNs is consistent with the $M_{BH}$ - $\sigma_{\star}$ relation of lower luminosity AGNs and inactive galaxies, after correcting for dust extinction and the rotational contribution to the stellar velocity dispersion.

We report on our combined analysis of HST, VLT/MUSE, VLT/SINFONI, and ALMA observations of the local Seyfert 2 galaxy, NGC 5728 to investigate in detail the feeding and feedback of the AGN. The datasets simultaneously probe the morphology, excitation, and kinematics of the stars, ionized gas, and molecular gas over a large range of spatial scales (10 pc--10 kpc). NGC 5728 contains a large stellar bar which is driving gas along prominent dust lanes to the inner 1 kpc where the gas settles into a circumnuclear ring. The ring is strongly star forming and contains a substantial population of young stars as indicated by the lowered stellar velocity dispersion and gas excitation consistent with HII regions. We model the kinematics of the ring using the velocity field of the CO (2--1) emission and stars and find it is consistent with a rotating disk. The outer regions of the disk, where the dust lanes meet the ring, show signatures of inflow at a rate of 1 M$_{\sun}$ yr$^{-1}$. Inside the ring, we observe three molecular gas components corresponding to the circular rotation of the outer ring, a warped disk, and the nuclear stellar bar. The AGN is driving an ionized gas outflow that reaches a radius of 250 pc with a mass outflow rate of 0.08 M$_{\sun}$ yr$^{-1}$ consistent with its luminosity and scaling relations from previous studies. While we observe distinct holes in CO emission which could be signs of molecular gas removal, we find that largely the AGN is not disrupting the structure of the circumnuclear region.

We present ALMA observations of the CO(2-1) and CO(3-2) molecular gas transitions and associated (sub)-mm continua of the nearby Seyfert 1.5 galaxy NGC3227 with angular resolutions 0.085-0.21" (7-15pc). On large scales the cold molecular gas shows circular motions as well as streaming motions on scales of a few hundred parsecs associated with a large scale bar. We fitted the nuclear ALMA 1.3mm emission with an unresolved component and an extended component. The 850$\mu$m emission shows at least two extended components, one along the major axis of the nuclear disk and the other along the axis of the ionization cone. The molecular gas in the central region (1" ~73pc) shows several CO clumps with complex kinematics which appears to be dominated by non-circular motions. While we cannot demonstrate conclusively the presence of a warped nuclear disk, we also detected non-circular motions along the kinematic minor axis. They reach line-of-sight velocities of v-vsys =150-200km/s. Assuming that the radial motions are in the plane of the galaxy, then we interpret them as a nuclear molecular outflow due to molecular gas in the host galaxy being entrained by the AGN wind. We derive molecular outflow rates of $5\,M_\odot\,{\rm yr}^{-1}$ and $0.6\,M_\odot\,{\rm yr}^{-1}$ at projected distances of up to 30pc to the northeast and southwest of the AGN, respectively. At the AGN location we estimate a mass in molecular gas of $5\times 10^{5}\,M_\odot$ and an average column density $N({\rm H}_2) = 2-3\times 10^{23}\,{\rm cm}^{-2}$ in the inner 15pc. The nuclear molecular gas and sub-mm continuum emission of NGC3227 do not resemble the classical compact torus. Rather, these emissions extend for several tens of parsecs and appear connected with the circumnuclear ring in the host galaxy disk, as found in other local AGN. (Abridged)

Using a sample of 472 local Universe (z<0.06) galaxies in the stellar mass range 10.25 < log M*/M_sun < 10.75, we explore the variation in galaxy structure as a function of morphology and galaxy colour. Our sample of galaxies is sub-divided into red, green and blue colour groups and into elliptical and non-elliptical (disk-type) morphologies. Using KiDS and VIKING derived postage stamp images, a group of eight volunteers visually classified bars, rings, morphological lenses, tidal streams, shells and signs of merger activity for all systems. We find a significant surplus of rings ($2.3\sigma$) and lenses ($2.9\sigma$) in disk-type galaxies as they transition across the green valley. Combined, this implies a joint ring/lens green valley surplus significance of $3.3\sigma$ relative to equivalent disk-types within either the blue cloud or the red sequence. We recover a bar fraction of ~44% which remains flat with colour, however, we find that the presence of a bar acts to modulate the incidence of rings and (to a lesser extent) lenses, with rings in barred disk-type galaxies more common by ~20-30 percentage points relative to their unbarred counterparts, regardless of colour. Additionally, green valley disk-type galaxies with a bar exhibit a significant $3.0\sigma$ surplus of lenses relative to their blue/red analogues. The existence of such structures rules out violent transformative events as the primary end-of-life evolutionary mechanism, with a more passive scenario the favoured candidate for the majority of galaxies rapidly transitioning across the green valley.

Using new APEX and JCMT spectroscopy of the CO 2-1 line, we undertake a controlled study of cold molecular gas in moderately luminous Active Galactic Nuclei (AGN) and inactive galaxies from the Luminous Local AGN with Matched Analogs (LLAMA) survey. We use spatially resolved infrared photometry of the LLAMA galaxies from 2MASS, WISE, IRAS & Herschel, corrected for nuclear emission using multi-component spectral energy distribution (SED) fits, to examine the dust-reprocessed star-formation rates (SFRs), molecular gas fractions and star formation efficiencies (SFEs) over their central 1 - 3 kpc. We find that the gas fractions and central SFEs of both active and inactive galaxies are similar when controlling for host stellar mass and morphology (Hubble type). The equivalent central molecular gas depletion times are consistent with the discs of normal spiral galaxies in the local Universe. Despite energetic arguments that the AGN in LLAMA should be capable of disrupting the observable cold molecular gas in their central environments, our results indicate that nuclear radiation only couples weakly with this phase. We find a mild preference for obscured AGN to contain higher amounts of central molecular gas, which suggests a connection between AGN obscuration and the gaseous environment of the nucleus. Systems with depressed SFEs are not found among the LLAMA AGN. We speculate that the processes that sustain the collapse of molecular gas into dense pre-stellar cores may also be a prerequisite for the inflow of material on to AGN accretion disks.

In a complete sample of local 14-195 keV selected AGNs and inactive galaxies, matched by their host galaxy properties, we study the spatially resolved stellar kinematics and luminosity distributions at near-infrared wavelengths on scales of 10-150 pc, using SINFONI on the VLT. In this paper, we present the first half of the sample, which comprises 13 galaxies, 8 AGNs and 5 inactive galaxies. The stellar velocity fields show a disk-like rotating pattern, for which the kinematic position angle is in agreement with the photometric position angle obtained from large scale images. For this set of galaxies, the stellar surface brightness of the inactive galaxy sample is generally comparable to the matched sample of AGN but extends to lower surface brightness. After removal of the bulge contribution, we find a nuclear stellar light excess with an extended nuclear disk structure, and which exhibits a size-luminosity relation. While we expect the excess luminosity to be associated with a dynamically cooler young stellar population, we do not typically see a matching drop in dispersion. This may be because these galaxies have pseudo-bulges in which the intrinsic dispersion increases towards the centre. And although the young stars may have an impact in the observed kinematics, their fraction is too small to dominate over the bulge and compensate the increase in dispersion at small radii, so no dispersion drop is seen. Finally, we find no evidence for a difference in the stellar kinematics and nuclear stellar luminosity excess between these active and inactive galaxies.

We discuss the environment of local hard X-ray selected active galaxies, with reference to two independent group catalogues. We find that the fraction of these AGN in S0 host galaxies decreases strongly as a function of galaxy group size (halo mass) - which contrasts with the increasing fraction of galaxies of S0 type in denser environments. However, there is no evidence for an environmental dependence of AGN in spiral galaxies. Because most AGN are found in spiral galaxies, this dilutes the signature of environmental dependence for the population as a whole. We argue that the differing results for AGN in disk-dominated and bulge-dominated galaxies is related to the source of the gas fuelling the AGN, and so may also impact the luminosity function, duty cycle, and obscuration. We find that there is a significant difference in the luminosity function for AGN in spiral and S0 galaxies, and tentative evidence for some difference in the fraction of obscured AGN.

We use high spectral resolution (R > 8000) data covering 3800-13000\r{A} to study the physical conditions of the broad line region (BLR) of nine nearby Seyfert 1 galaxies. Up to six broad HI lines are present in each spectrum. A comparison - for the first time using simultaneous optical to near-infrared observations - to photoionisation calculations with our devised simple scheme yields the extinction to the BLR at the same time as determining the density and photon flux, and hence distance from the nucleus, of the emitting gas. This points to a typical density for the HI emitting gas of 10$^{11}$cm$^{-3}$ and shows that a significant amount of this gas lies at regions near the dust sublimation radius, consistent with theoretical predictions. We also confirm that in many objects the line ratios are far from case B, the best-fit intrinsic broad-line H$\alpha$/H$\beta$ ratios being in the range 2.5-6.6 as derived with our photoionization modeling scheme. The extinction to the BLR, based on independent estimates from HI and HeII lines, is A$_V$ $\le$ 3 for Seyfert 1-1.5s, while Seyfert 1.8-1.9s have A$_V$ in the range 4-8. A comparison of the extinction towards the BLR and narrow line region (NLR) indicates that the structure obscuring the BLR exists on scales smaller than the NLR. This could be the dusty torus, but dusty nuclear spirals or filaments could also be responsible. The ratios between the X-ray absorbing column N$_H$ and the extinction to the BLR are consistent with the Galactic gas-to-dust ratio if N$_H$ variations are considered.

We describe a complete volume limited sample of nearby active galaxies selected by their 14-195keV luminosity, and outline its rationale for studying the mechanisms regulating gas inflow and outflow. We describe also a complementary sample of inactive galaxies, selected to match the AGN host galaxy properties. The active sample appears to have no bias in terms of AGN type, the only difference being the neutral absorbing column which is two orders of magnitude greater for the Seyfert 2s. In the luminosity range spanned by the sample, log L_{14-195keV} [erg/s] = 42.4-43.7, the optically obscured and X-ray absorbed fractions are 50-65%. The similarity of these fractions to more distant spectroscopic AGN samples, although over a limited luminosity range, suggests that the torus does not strongly evolve with redshift. Our sample confirms that X-ray unabsorbed Seyfert 2s are rare, comprising not more than a few percent of the Seyfert 2 population. At higher luminosities, the optically obscured fraction decreases (as expected for the increasing dust sublimation radius), but the X-ray absorbed fraction changes little. We argue that the cold X-ray absorption in these Seyfert 1s can be accounted for by neutral gas in clouds that also contribute to the broad line region (BLR) emission; and suggest that a geometrically thick neutral gas torus co-exists with the BLR and bridges the gap to the dusty torus.

We analyse the 2-dimensional distribution and kinematics of the stars as well as molecular and ionised gas in the central few hundred parsecs of 5 active and 5 matched inactive galaxies. The equivalent widths of the Br-gamma line indicate there is no on-going star formation in their nuclei, although recent (terminated) starbursts are possible in the active galaxies. The stellar velocity fields show no signs of non-circular motions, while the 1-0S(1) H_2 kinematics exhibit significant deviations from simple circular rotation. In the active galaxies the H_2 kinematics reveal inflow and outflow superimposed on disk rotation. Steady-state circumnuclear inflow is seen in three AGN, and hydrodynamical models indicate it can be driven by a large scale bar. In three of the five AGN, molecular outflows are spatially resolved. The outflows are oriented such that they intersect, or have an edge close to, the disk - which may be the source of molecular gas in the outflow. The relatively low speeds imply the gas will fall back onto the disk; and with moderate outflow rates, they will have only a local impact on the host galaxy. H_2 was detected in two inactive galaxies. These exhibit chaotic circumnuclear dust morphologies and have molecular structures that are counter-rotating with respect to the main gas component, which could lead to gas inflow in the near future. In our sample, all four galaxies with chaotic dust morphology in the circumnuclear region exist in moderately dense groups with 10-15 members where accretion of stripped gas can easily occur.

We present an analysis of the H2 emission-line gas kinematics in the inner < 4 arcsec radius of six nearby spiral galaxies, based on AO-assisted integral-field observations obtained in the K-band with SINFONI/VLT. Four of the six galaxies in our sample display ordered H2 velocity fields, consistent with gas moving in the plane of the galaxy and rotating in the same direction as the stars. However, the gas kinematics is typically far from simple circular motion. We can classify the observed velocity fields into four different types of flows, ordered by increasing complexity: (1) circular motion in a disc (NGC3351); (2) oval motion in the galaxy plane (NGC3627 and NGC4536); (3) streaming motion superimposed on circular rotation (NGC4501); and (4) disordered streaming motions (NGC4569 and NGC4579). The H2 velocity dispersion in the galaxies is usually higher than 50 km/s in the inner 1-2 arcsec radii. The four galaxies with ordered kinematics have v/sigma < 1 at radii less than 40-80 pc. The radius at which v/sigma = 1 is independent of the type of nuclear activity. While the low values of v/sigma could be taken as an indication of a thick disc in the innermost regions of the galaxies, other lines of evidence (e.g. H2 morphologies and velocity fields) argue for a thin disc interpretation in the case of NGC3351 and NGC4536. We discuss the implications of the high values of velocity dispersion for the dynamics of the gaseous disc and suggest caution when interpreting the velocity dispersion of ionized and warm tracers as being entirely dynamical. Understanding the nature and role of the velocity dispersion in the gas dynamics, together with the full 2D information of the gas, is essential for obtaining accurate black hole masses from gas kinematics.

Although nearly one-third of barred galaxies host an inner, secondary bar, the formation and evolution of double barred galaxies remain unclear. We show here an example model of a galaxy, dominated by a live dark matter halo, in which double bars form naturally, without requiring gas, and we follow its evolution for a Hubble time. The inner bar in our model galaxy rotates almost as slowly as the outer bar, and it can reach up to half of its length. The route to the formation of a double bar may be different from that of a single strong bar. Massive dark matter halo or dynamically hot stellar disc may play an important role in the formation of double bars and their subsequent evolution.

In a matched sample of Seyfert and quiescent galaxies we simultaneously probe the stellar and molecular gas kinematics from 1 kpc down to 50 pc with the aim of identifying the dynamical processes dictating black hole accretion rates. This first paper compares the global characteristics of a sample of ten galaxies. We find several differences within a radius of 500 pc that are correlated with Active Galactic Nucleus (AGN) activity. The Seyferts have: (1) a more centrally concentrated nuclear stellar surface brightness with lower stellar luminosities beyond a radius of 100 pc, (2) a lower stellar velocity dispersion within a radius of 200 pc, (3) elevated H2 luminosity out to a radius of at least 250 pc, and (4) more centrally concentrated H2 surface brightness profiles. These observed differences can be interpreted as evidence for Seyfert galaxies having a dynamically cold (in comparison to the bulge) nuclear structure composed of a significant gas reservoir and a relatively young stellar population. This structure is undetected (and possibly does not exist) in quiescent galaxies. The presence of such a nuclear structure in Seyfert galaxies provides evidence for inflow of the surrounding interstellar medium since the nuclear stellar population requires a supply of gas from which to form. The fueling of a Seyfert AGN is thus associated with the formation of a dynamically cold component of gas and stars on scales of hundreds of parsecs.

We present a new method that derives both velocity components in the equatorial plane of a barred stellar disc from the observed line-of-sight velocity, assuming geometry of a thin disc. The method can be applied to large departures from circular motion, and does not require multipole decomposition. It is based on assumptions that the bar is close to steady-state (i.e. does not evolve fast), and that both morphology and kinematics are symmetrical with respect to the major axis of the bar. We derive the equations used in the method, and analyze the effect of observational errors on the inferred velocity fields. We show that this method produces meaningful results via a simple toy model. We also apply the method on integral-field data of NGC 936, for which we recover both velocity components in the disc. Knowing both velocity components in the disc, i.e. the non-observable transverse velocity in addition to the line-of-sight velocity, puts additional constraints on dynamical models and allows for new ways of determining parameters that are crucial in characterizing galaxies.

Using the Inamori Magellan Areal Camera and Spectrograph (IMACS) integral-field unit (IFU) on the 6.5m Magellan telescope, we have designed the first statistically significant investigation of the two-dimensional distribution and kinematics of ionized gas and stars in the central kiloparsec regions of a well-matched sample of Seyfert and inactive control galaxies selected from the Sloan Digital Sky Survey. The goals of the project are to use the fine spatial sampling (0.2 arcsec/pixel) and large wavelength coverage (4000-7000A) of the IMACS-IFU to search for dynamical triggers of nuclear activity in the central region where active galactic nucleus (AGN) activity and dynamical timescales become comparable, to identify and assess the impact of AGN-driven outflows on the host galaxy and to provide a definitive sample of local galaxy kinematics for comparison with future three-dimensional kinematic studies of high-redshift systems. In this paper, we provide the first detailed description of the procedure to reduce and calibrate data from the IMACS-IFU in `long mode' to obtain two-dimensional maps of the distribution and kinematics of ionized gas and stars. The sample selection criteria are presented, observing strategy described and resulting maps of the sample galaxies presented along with a description of the observed properties of each galaxy and the overall observed properties of the sample.

We report the discovery of a two-armed mini-spiral structure within the inner kiloparsec of the E0 LINER/Seyfert 1 galaxy Arp102B. The arms are observed in H-alpha emission and located East and West of the nucleus, extending up to about 1 kpc from it. We use narrow-band imaging from the Hubble Space Telescope Advanced Camera for Surveys, in combination with archival VLA radio images at 3.6 and 6 cm to investigate the origin of the nuclear spiral. From the H-alpha luminosity of the spiral, we obtain an ionized gas mass of the order of one million solar masses. One possibility is that the nuclear spiral represents a gas inflow triggered by a recent accretion event which has replenished the accretion disk, giving rise to the double-peaked emission-line profiles characteristic of Arp102B. However, the radio images show a one-sided curved jet which correlates with the eastern spiral arm observed in the H-alpha image. A published milliarcsecond radio image also shows one-sided structure at position angle about 40 degrees, approximately aligned with the inner part of the eastern spiral arm. The absence of a radio counter-part to the western spiral arm is tentatively interpreted as indicating that the jet is relativistic, with an estimated speed of 0.45c. Estimates of the jet kinetic energy and the ionizing luminosity of the active nucleus indicate that both are capable of ionizing the gas along the spiral arms. We conclude that, although the gas in the nuclear region may have originated in an accretion event, the mini-spiral is most likely the result of a jet-cloud interaction rather than an inflowing stream.

We analyze how the orbital support of the inner bar in a double-barred galaxy (nested bars) depends on the angular velocity (i.e. pattern speed) of this bar. We study orbits in seven models of double bars using the method of invariant loops. The range of pattern speed is covered exhaustively. We find that not all pattern speeds are allowed when the inner bar rotates in the same direction as the outer bar. Below a certain minimum pattern speed orbital support for the inner bar abruptly disappears, while at high values of this speed the orbits indicate an increasingly round bar that looks more like a twist in the nuclear isophotes than a dynamically independent component. For values between these two extremes, orbits supporting the inner bar extend further out as the bar's pattern speed decreases, their corresponding loops become more eccentric, pulsate more, and their rotation becomes increasingly non-uniform, as they speed up and slow down in their motion. Lower pattern speeds also lead to a less coherent bar, as the pulsation and acceleration increasingly varies among the loops supporting the inner bar. The morphologies of fast and slow inner bars expected from the orbital structure studied here are recently recovered observationally by decomposition of double barred galaxies. Our findings allow us to link the observed morphology to the dynamics of the inner bar.

Our adaptive optics observations of nearby AGN at spatial resolutions as small as 0.085arcsec show strong evidence for recent, but no longer active, nuclear star formation. We begin by describing observations that highlight two contrasting methods by which gas can flow into the central tens of parsecs. Gas accumulation in this region will inevitably lead to a starburst, and we discuss the evidence for such events. We then turn to the impact of stellar evolution on the further inflow of gas by combining a phenomenological approach with analytical modelling and hydrodynamic simulations. These complementary perspectives paint a picture in which all the processes are ultimately regulated by the mass accretion rate into the central hundred parsecs, and the ensuing starburst that occurs there. The resulting supernovae delay accretion by generating a starburst wind, which leaves behind a clumpy interstellar medium. This provides an ideal environment for slower stellar outflows to accrete inwards and form a dense turbulent disk on scales of a few parsecs. Such a scenario may resolve the discrepancy between the larger scale structure seen with adaptive optics and the small scale structure seen with VLTI.

We present spatially resolved distributions and kinematics of the stars and molecular gas in the central 320pc of NGC1097. The stellar continuum confirms the previously reported 3-arm spiral pattern extending into the central 100pc. The stellar kinematics and the gas distribution imply this is a shadowing effect due to extinction by gas and dust in the molecular spiral arms. The molecular gas kinematics show a strong residual (i.e. non-circular) velocity, which is manifested as a 2-arm kinematic spiral. Linear models indicate that this is the line-of-sight velocity pattern expected for a density wave in gas that generates a 3-arm spiral morphology. We estimate the inflow rate along the arms. Using hydrodynamical models of nuclear spirals, we show that when deriving the accretion rate into the central region, outflow in the disk plane between the arms has to be taken into account. For NGC1097, despite the inflow rate along the arms being ~1.2Msun/yr, the net gas accretion rate to the central few tens of parsecs is much smaller. The numerical models indicate that the inflow rate could be as little as ~0.06Msun/yr. This is sufficient to generate recurring starbursts, similar in scale to that observed, every 20-150Myr. The nuclear spiral represents a mechanism that can feed gas into the central parsecs of the galaxy, with the gas flow sustainable for timescales of a Gigayear.

The method to study oscillating potentials of double bars, based on invariant loops, is introduced here in a new way, intended to be more intelligible. Using this method, I show how the orbital structure of a double-barred galaxy (nested bars) changes with the variation of nuclear bar's pattern speed. Not all pattern speeds are allowed when the inner bar rotates in the same direction as the outer bar. Below certain minimum pattern speed orbital support for the inner bar abruptly disappears, while high values of this speed lead to loops that are increasingly round. For values between these two extremes, loops supporting the inner bar extend further out as its pattern speed decreases, and they become more eccentric and pulsate more. These findings do not apply to counter-rotating inner bars.

We show that stable double-frequency orbits form the backbone of double bars, because they trap around themselves regular orbits, as stable closed periodic orbits do in single bars, and in both cases the trapped orbits occupy similar volume of phase-space. We perform a global search for such stable double-frequency orbits in a model of double bars by constructing maps of trajectories with initial conditions well sampled over the available phase-space. We use the width of a ring sufficient to enclose a given map as the indicator of how tightly the trajectory is trapped around a double-frequency orbit. We construct histograms of these ring widths in order to determine the fraction of phase-space occupied by ordered motions. We build 22 further models of double bars, and we construct histograms showing the fraction of the phase-space occupied by regular orbits in each model. Our models indicate that resonant coupling between the bars may not be the dominant factor reducing chaos in the system.

We report on our attempts to achieve a nearly steady-state gas flow in hydrodynamical simulations of doubly barred galaxies. After exploring the parameter space, we construct two models, for which we evaluate the photometric and the kinematic integrals, present in the Tremaine-Weinberg method, in search of observational signatures of two rotating patterns. We show that such signatures are often present, but a direct fit to data points is likely to return incorrect pattern speeds. However, for a particular distribution of the tracer, presented here, the values of the pattern speeds can be retrieved reliably even with the direct fit.

The backbone of double bars is made out of double-frequency orbits, and loops, their maps, indicate the bars' extent, morphology and dynamics.

In central regions of non-axisymmetric galaxies high-resolution hydrodynamical simulations indicate spiral shocks, which are capable of transporting gas inwards. The efficiency of transport is lower at smaller radii, therefore instead of all gas dropping onto the galactic centre, a roughly uniform distribution of high-density gas develops in the gaseous nuclear spiral downstream from the shock, and the shear in gas is very low there. These are excellent conditions for star formation. This mechanism is likely to contribute to the process of (pseudo-) bulge formation.

Bars in galaxies are mainly supported by particles trapped around stable periodic orbits. These orbits represent oscillatory motion with only one frequency, which is the bar driving frequency, and miss free oscillations. We show that a similar situation takes place in double bars: particles get trapped around parent orbits, which in this case represent oscillatory motion with two frequencies of driving by the two bars, and which also lack free oscillations. Thus the parent orbits, which constitute the backbone of an oscillating potential of two independently rotating bars, are the double-frequency orbits. These orbits do not close in any reference frame, but they map onto loops, first introduced by Maciejewski & Sparke (1997). Trajectories trapped around the parent double-frequency orbit map onto a set of points confined within a ring surrounding the loop.

We present new Space Telescope Imaging Spectrograph (STIS) observations of three spiral galaxies, NGC 4303, NGC 3310 and NGC 4258. The bright optical emission lines H$\alpha$ $\lambda$ $6564 \AA$, [NII] $\lambda$$\lambda$ $6549,6585 \AA$ and [SII] $\lambda$$\lambda$ $ 6718,6732 \AA$ were used to study the kinematics of the ionized gas in the nuclear region of each galaxy with a $\sim 0.07\arcsec$ spatial resolution. In NGC 3310, the observed gas kinematics is well matched by a circularly rotating disk model but we are only able to set an upper limit to the BH mass which, taking into account the allowed disk inclinations, varies in the range $5.0 \times 10^{6} - 4.2 \times 10^{7} M_{\odot}$ at the 95% confidence level. In NGC 4303 the kinematical data require the presence of a BH with mass $M_{BH}=(5.0)^{+0.87}_{-2.26}\times 10^{6}M_{\odot}$ (for a disk inclination $i=70$ deg).In NGC 4258, the observed kinematics require the presence of a black hole with $M_{BH}= (7.9)^{+6.2}_{-3.5} \times 10^{7}M_{\odot}$ ($i=60$ deg). This result is in good agreement with the published value $(3.9 \pm 0.1) \times 10^{7} M_{\odot}$, derived from $H_{2}O$-maser observations. Our attempt at measuring BH masses in these 3 late type Sbc spiral galaxies has shown that these measurements are very challenging and at the limit of the highest spatial resolution currently available. Nonetheless our estimates are in good agreement with the scaling relations between black holes and their host spheroids suggesting that (i) they are reliable and (ii) black holes in spiral galaxies follows the same scaling relations as those in more massive early-type galaxies. A crucial test for the gas kinematical method, the correct recovery of the known BH mass in NGC 4258, has been successful. [abridged]

Recent high-resolution observations indicate that nuclear spirals are often present in the innermost few hundred parsecs of disc galaxies. My models show that nuclear spirals form naturally as a gas response to non-axisymmetry in the gravitational potential. Some nuclear spirals take the form of spiral shocks, resulting in streaming motions in the gas, and in inflow comparable to the accretion rates needed to power local Active Galactic Nuclei. Recently streaming motions of amplitude expected from the models have been observed in nuclear spirals, confirming the role of nuclear spirals in feeding of the central massive black holes.

When integrals in the standard Tremaine-Weinberg method are evaluated for the case of a realistic model of a doubly barred galaxy, their modifications introduced by the second rotating pattern are in accord with what can be derived from a simple extension of that method, based on separation of tracer's density. This extension yields a qualitative argument that discriminates between prograde and retrograde inner bars. However, the estimate of the value of inner bar's pattern speed requires further assumptions. When this extension of the Tremaine-Weinberg method is applied to the recent observation of the doubly barred galaxy NGC 2950, it indicates that the inner bar there is counter-rotating, possibly with the pattern speed of -140 +/- 50 km/s/arcsec. The occurrence of counter-rotating inner bars can constrain theories of galaxy formation.

If supermassive black holes in centres of galaxies form by merging of black-hole remnants of massive Population III stars, then there should be a few black holes of mass one or two orders of magnitude smaller than that of the central ones, orbiting around the centre of a typical galaxy. These black holes constitute a weak perturbation in the gravitational potential, which can generate wave phenomena in gas within a disc close to the centre of the galaxy. Here we show that a single orbiting black hole generates a three-arm spiral pattern in the central gaseous disc. The density excess in the spiral arms in the disc reaches values of 3-12% when the orbiting black hole is about ten times less massive than the central black hole. Therefore the observed density pattern in gas can be used as a signature in detecting the most massive orbiting black holes.

Near-infrared images of the prototype LINER / Seyfert 1 galaxy NGC 1097 observed with the Very Large Telescope (VLT) using adaptive optics disclose with unprecedented detail a complex central network of filamentary structure spiralling down to the center of the galaxy. The structure, consisting of several spiral arms, some almost completing a revolution about the center, is most prominent within the radius of about 300 pc. Some filaments can be traced further out, where they seem to connect with the nuclear star-forming ring at 0.7 kpc radius. Straight principal shocks running along the primary large-scale bar of NGC 1097, seen in the optical images as prominent dust lanes, curve into this ring, but radio polarization vectors cross the nuclear ring under a rather large angle. Here we attempt to explain this morphology in terms of three-dimensional gas flow in a barred galaxy. In our scenario, parts of the principal shock, which propagate in the off-plane gas, can cross the nuclear star-forming ring, and excite waves inward from it. If the dispersion relation of the excited waves allows for their propagation, they will naturally take the shape of the observed central spiral. The nuclear region of NGC 1097 remains unresolved at sub-arcsec scales in the near-IR, with an upper size limit of <10 pc FWHM. Thus, any putative central dusty torus or gaseosus disk envisaged by the AGN unified schemes has to be smaller than 10 pc in diameter at near-IR wavelengths. The extinction in the region between the nuclear star-forming ring and the nucleus increases very moderately, reaching A_v~1 at the immediate surrounding of the nucleus. Thus, if the nuclear filaments are tracing cold dust, they contribute to a very low extinction in the line of sight and are likely to be distributed in a rather thin disk.

We investigate the central regions of 23 spiral galaxies using archival NICMOS imaging and STIS spectroscopy. The sample is taken from our program to determine the masses of central massive black holes (MBH) in 54 nearby spiral galaxies. Stars are likely to contribute significantly to any dynamical central mass concentration that we find in our MBH program and this paper is part of a series to investigate the nuclear properties of these galaxies. We use the Nuker law to fit surface brightness profiles, derived from the NICMOS images, to look for nuclear star clusters and find possible extended sources in 3 of the 23 galaxies studied (13 per cent). The fact that this fraction is lower than that inferred from optical Hubble Space Telescope studies is probably due to the greater spatial resolution of those studies. Using R-H and J-H colors and equivalent widths of H-alpha emission (from the STIS spectra) we investigate the nature of the stellar population with evolutionary models. Under the assumption of hot stars ionizing the gas, as opposed to a weak AGN, we find that there are young stellar populations (~10-20 Myr) however these data do not allow us to determine what percentage of the total nuclear stellar population they form. Also, in an attempt to find any unknown AGN we use [N II] and [S II] line flux ratios (relative to H-alpha) and find tentative evidence for weak AGN in NGC 1300 and NGC 4536.

We present Space Telescope Imaging Spectrograph emission-line spectra of the central regions of the spiral galaxies NGC 1300 and NGC 2748. From the derived kinematics of the nuclear gas we have found evidence for central supermassive black holes in both galaxies. The estimated mass of the black hole in NGC 1300 is 6.6 (+6.3, -3.2) x 10^7 solar masses and in NGC 2748 is 4.4 (+3.5, -3.6) x 10^7 solar masses (both at the 95% confidence level). These two black hole mass estimates contribute to the poorly sampled low-mass end of the nuclear black hole mass spectrum.

We present surface photometry for the central regions of a sample of 48 spiral galaxies (mostly unbarred and barred of types Sbc or Sc) observed with the Space Telescope Imaging Spectrograph on board the Hubble Space Telescope. Surface brightness profiles were derived and modeled with a Nuker law. We also analyzed archival Wide Field Planetary Camera~2 images with a larger field of view, available for 18 galaxies in our sample. We modeled the extracted bulge surface brightness profiles with an exponential, a de Vaucouleurs or a Sersic profile. In agreement with previous studies, we find that bulges of Sbc galaxies fall into two categories: bulges well described by an exponential profile and those well described by an de Vaucouleurs profile. Only one galaxy requires the use of a more general Sersic profile to properly describe the bulge. Nuclear photometrically distinct components are found in ~55% of the galaxies. For those that we classify as star clusters based on their resolved extent we find absolute magnitudes that are brighter on average than those previously identified in spiral galaxies. This might be due to a bias in our sample toward star forming galaxies, combined with a trend for star forming galaxies to host brighter central clusters.

Nuclear spirals naturally form as a gas response to non-axisymmetry in the galactic potential, even if the degree of this asymmetry is very small. Linear wave theory well describes weak nuclear spirals, but spirals induced by stronger asymmetries in the potential are clearly beyond the linear regime. Hydrodynamical models indicate spiral shocks in this latter case that, depending on how the spiral intersects the x2 orbits, either get damped, leading to the formation of the nuclear ring, or get strengthened, and propagate towards the galaxy centre. Central massive black hole of sufficient mass can allow the spiral shocks to extend all the way to its immediate vicinity, and to generate gas inflow up to 0.03 M_sun/yr, which coincides with the accretion rates needed to power luminous local Active Galactic Nuclei.

Nuclear spirals can provide a wealth of information about the nuclear potential in disc galaxies. They are unlikely to form in nuclei with solid-body rotation, yet they are present in a majority of galactic centres. Their morphology varies depending on whether a central massive black hole (MBH) is present or absent in the galaxy. In this paper I consider predictions of the linear theory for waves induced in gas by asymmetric gravitational potential, which are applicable to the nuclear spirals observed in galaxies. Generation and propagation of waves is governed by dynamical resonances, and inclusion of a MBH can move or even create resonances, greatly altering the extent and shape of the nuclear spiral. I will use predictions of the linear theory presented here as a guideline when interpreting hydrodynamical models in the second paper of this series. I also comment on modifications that self-gravity in gas imposes on the induced waves.

Nuclear spirals can provide a wealth of information about the nuclear potential in disc galaxies. They form naturally as a gas response to non-axisymmetry in the gravitational potential, even if the degree of this asymmetry is very small. Linear wave theory well describes weak nuclear spirals, but stronger asymmetries in the potential induce waves beyond the linear regime, which appear as spiral shocks. If a central massive black hole (MBH) is present, spiral shocks can extend all the way to its immediate vicinity, and generate gas inflow up to 0.03 Msun/yr. This coincides with the accretion rates needed to power local Active Galactic Nuclei.

Bars in galaxies are mainly supported by particles trapped around closed periodic orbits. These orbits respond to the bar's forcing frequency only and lack free oscillations. We show that a similar situation takes place in double bars: particles get trapped around orbits which only respond to the forcing from the two bars and lack free oscillations. We find that writing the successive positions of a particle on such an orbit every time the bars align generates a closed curve, which we call a loop. Loops allow us to verify consistency of the potential. As maps of doubly periodic orbits, loops can be used to search the phase-space in double bars in order to determine the fraction occupied by ordered motions.

Asymmetries in galactic potentials, either self-induced, or caused by a passing companion, play an important role in global gas dynamics in galaxies. In particular, they are able to trigger gas inflow, which in turn feeds nuclear activity. In the inner kiloparsec, the inflowing gas becomes subject to various resonances induced by asymmetry in the potential. This produces complicated gas morphology and dynamics seen in the nuclear rings and spirals. Formation of nuclear rings is related to barred galaxies, and explained by the orbital structure in bars. One can get basic understanding of weak nuclear spirals from the linear wave theory, but stronger spirals are likely to be out of the linear regime, and they appear as spiral shocks in hydrodynamical models. In addition, nuclear stellar bars are observed in a considerable fraction of disc galaxies, often nested inside the large bar; their orbital structure further modifies the nuclear gas dynamics. If a massive black hole (MBH) is present in the galactic centre, it governs the resonances even beyond its classically defined sphere of influence. Since resonances shape gas dynamics in the nuclear region, the central black hole should be able to regulate gas flow around itself.

We present HST/STIS spectra of the Sbc spiral galaxy NGC 4041 which were used to map the velocity field of the gas in its nuclear region. We detect the presence of a compact (r~0.4" ~40 pc), high surface brightness, rotating nuclear disk co-spatial with a nuclear star cluster. The disk is characterized by a rotation curve with a peak to peak amplitude of ~40 km/s and is systematically blueshifted by 10 - 20 km/s with respect to the galaxy systemic velocity. With the standard assumption of constant mass-to-light ratio and with the nuclear disk inclination taken from the outer disk, we find that a dark point mass of 1(-0.7;+0.6) 10^7 Msun is needed to reproduce the observed rotation curve. However the observed blueshift suggests the possibility that the nuclear disk could be dynamically decoupled. Following this line of reasoning we relax the standard assumptions and find that the kinematical data can be accounted for by the stellar mass provided that either the central mass-to-light ratio is increased by a factor of ~2 or that the inclination is allowed to vary. This model results in a 3 sigma upper limit of 6 10^6 Msun on the mass of any nuclear black hole. Overall, our analysis only allows us to set an upper limit of 2 10^7 Msun on the mass of the nuclear black hole. If this upper limit is taken in conjunction with an estimated bulge B magnitude of -17.7 and with a central stellar velocity dispersion of 95 km/s, then these results are not inconsistent with both the MBH-Lsph and the MBH-sigma correlations. Constraints on black hole masses in spiral galaxies of types as late as Sbc are still very scarce and therefore the present result adds an important new datapoint to our understanding of black hole demography.

The bulk of the molecular line emission in the double barred galaxy NGC4303 as observed in its CO(1-0) line with the OVRO mm-interferometer comes from two straight gas lanes which run north-south along the leading sides of the large-scale primary bar. Inside a radius of ~ 400 pc the molecular gas forms a spiral pattern which, for the northern arm, can be traced to the nucleus. Comparison of the OVRO and archival HST data with dynamical models of gas flow in the inner kiloparsec of single- and double-barred galaxies shows that the observed global properties of the molecular gas are in agreement with models for the gas flow in a strong, large-scale bar, and the two-arm spiral structure seen in CO in the inner kiloparsec can already be explained by a density wave initiated by the potential of that bar. Only a weak correlation between the molecular gas distribution and the extinction seen in the HST V-H map is found in the inner 400 pc of NGC4303: The innermost part of one arm of the nuclear CO spiral correlates with a weak dust filament in the color map, while the overall dust distribution follows a ring or single-arm spiral pattern well correlated with the UV continuum. This complicated nuclear geometry of the stellar and gaseous components allows for two scenarios: (A) A self-gravitating m=1 mode is present forming the spiral structure seen in the UV continuum. In this case the gas kinematics would be unaffected by the small (~ 4'') inner bar. (B) The UV continuum traces a complete ring which is heavily extincted north of the nucleus. Such a ring forms in hydrodynamic models of double bars, but the models cannot account for the UV emission observed on the leading side of the inner bar. (abridged)

A wide-spread belief that nested bars enhance gas inflow to the galactic centre has recently been contradicted by dynamical models in which inner bars seem to prohibit such inflow. Can the existing models of dynamically possible double bars be modified to enable strong inflow in the secondary bar? I present here simple dynamical arguments which imply that in general, double bars in resonant coupling do not enhance gas inflow. However, stronger inflow with straight shocks in the inner bar can occur if there is no resonant coupling of the commonly assumed form between the bars.

We report here results of high-resolution hydrodynamical simulations of gas flows in barred galaxies, with a focus on gas dynamics in the central kiloparsec. In a single bar with an Inner Lindblad Resonance, we find either near-circular motion of gas in the nuclear ring, or a spiral shock extending towards the galaxy center, depending on the sound speed in the gas. From a simple model of a dynamically-possible doubly barred galaxy with resonant coupling, we infer that the secondary bar is likely to end well inside its corotation. Such a bar cannot create shocks in the gas flow, and therefore will not reveal itself in color maps through straight dust lanes: the gas flows induced by it are different from those caused by the rapidly rotating main bars. In particular, we find that secondary stellar bars are unlikely to increase the mass inflow rate into the galactic nucleus.

With a help of hydrodynamical models, we explore fundamental modes of gas flow in central parts of a galaxy potential with non-axisymmetric perturbations of various strengths. Our grid-based algorithm allows us to get a detailed picture of gas dynamics down to about 10 pc from the galaxy center, where the central black hole becomes dominant. Gas inflow in a bar with an Inner Lindblad Resonance often stops at the nuclear ring. We find that embedded bars in resonant coupling are unlikely to increase the inflow, but for certain gas parameters in a single bar, inflow to the center develops in a spiral shock, similar to the patterns observed in nuclear discs of real galaxies. The spiral pattern persists for very weak non-axisymmetries of the potential.

We describe an on-going HST program aimed at determining the relationship between the nuclear black hole mass and bulge mass in spiral galaxies. We have selected a volume limited sample of 54 nearby spiral galaxies for which we already have ground based emission line rotation curves, CCD surface photometry and radio maps. We are now obtaining HST/STIS longslit observations of each of the galaxies in the sample in order to determine the nuclear Halpha rotation curve at high (~0.1") spatial resolution. We will use these data to measure the unresolved dark mass concentration at the nucleus of each object. Here we show the first results from observations of objects in the sample.

Motivated by STIS observations of more than 50 nearby galactic nuclei, we consider long-slit emission-line spectra when the slit is wider than the instrumental PSF, and the target has arbitrarily large velocity gradients. The finite width of the slit generates complex patterns in the spectra that can be misinterpreted as coming from various physically distinct nuclear components, but when interpreted correctly, they can have considerable diagnostic power. For a thin disk in circular motion around a central galactic black hole (BH), a characteristic artifact occurs in the spectrum at the outer edge of the BH's sphere of influence. It betrays the presence of a BH, and allows us to develop a new method for estimating its mass, which gives higher sensitivity to BH detection than traditional methods.

We consider long-slit emission-line spectra of galactic nuclei when the slit is wider than the instrumental PSF, and the target has large velocity gradients. The finite width of the slit generates complex distributions of brightness at a given spatial point in the measured spectrum, which can be misinterpreted as coming from additional physically distinct nuclear components. We illustrate this phenomenon for the case of a thin disc in circular motion around a nuclear black hole (BH). We develop a new method for estimating the mass of the BH that exploits a feature in the spectrum at the outer edge of the BH's sphere of influence, and therefore gives higher sensitivity to BH detection than traditional methods. Moreover, with this method we can determine the black hole mass and the inclination of the surrounding disc separately, whereas the traditional approach to black-hole estimation requires two long-slit spectra to be taken. We show that with a given spectrograph, the detectability of a BH depends on the sense of rotation of the nuclear disc. We apply our method to estimate the BH mass in M84 from a publicly available spectrum, and recover a value 4 times lower than that published previously from the same data.

Motivated by STIS observations of more than 50 nearby galactic nuclei, we consider long-slit emission-line spectra when the slit is wider than the instrumental PSF. The finite width of the slit generates complex patterns in the spectra that can have considerable diagnostic power if they are modeled with adequate sophistication. They allowed us to develop a new method for estimating the BH mass that gives fuller information and higher sensitivity to BH detection than traditional methods.

High-resolution observations of the inner regions of barred disk galaxies have revealed many asymmetrical, small-scale central features, some of which are best described as secondary bars. Because orbital time-scales in the galaxy center are short, secondary bars are likely to be dynamically decoupled from the main kiloparsec-scale bars. Here, we show that regular orbits exist in such doubly-barred potentials and that they can support the bars in their motion. We find orbits in which particles remain on loops: closed curves which return to their original positions after two bars have come back to the same relative orientation. Stars trapped around stable loops could form the building blocks for a long-lived, doubly-barred galaxy. Using the loop representation, we can find which orbits support the bars in their motion, and what are the constraints on the sizes and shapes of self-consistent double bars. In particular, it appears that a long-lived secondary bar may exist only when an Inner Lindblad Resonance is present in the primary bar, and that it would not extend beyond this resonance.

Active Galactic Nuclei (AGN) require mass accretion onto the central engine. On large galactic scales, torques from a stellar bar can efficiently remove angular momentum from gas, and cause it to move inwards along two hydrodynamical shocks on the leading edges of the bar. Inflowing gas settles on near-circular orbits around the Inner Lindblad Resonance (ILR), and forms a nuclear ring, about 1 kpc in size. A secondary bar inside the main one, with its own pair of shocks, has been proposed to drive further inflow, and thus feed the AGN in a manner similar to the inflow on large scales. Here, we report results of high resolution hydrodynamical simulations, where we examine the nature of the nuclear ring, and check how efficient double bars can be in fueling AGNs.

We propose simple, explicit, analytical approximations for the kinematics of an adiabatic blast wave propagating in an exponentially stratified ambient medium, and for the onset of radiative cooling, which ends the adiabatic era. Our method, based on the Kompaneets implicit solution and the Kahn approximation for the radiative cooling coefficient, gives straightforward estimates for the size, expansion velocity, and progression of cooling times over the surface, when applied to supernova remnants (SNRs). The remnant shape is remarkably close to spherical for moderate density gradients, but even a small gradient in ambient density causes the cooling time to vary substantially over the remnant's surface, so that for a considerable period there will be a cold dense expanding shell covering only a part of the remnant. Our approximation provides an effective tool for identifying the approximate parameters when planning 2-dimensional numerical models of SNRs, the example of W44 being given in a subsequent paper.

Inner regions of barred disk galaxies often include asymmetrical, small-scale central features, some of which are best described as secondary bars. Because orbital timescales in the galaxy center are short, secondary bars are likely to be dynamically decoupled from the main kiloparsec-scale bars. We found that non-chaotic multiply-periodic particle orbits can exist in potentials with two dynamically decoupled bars.Stars trapped around these orbits could form the building blocks for a long-lived, doubly-barred galaxy. A self-consistent secondary bar appears to induce formation of inner gaseous rings rather than shocks in gas flow.

(shortened version) We show that many observations of W44, a supernova remnant in the galactic plane at a distance of about 2500 pc, are remarkably consistent with the simplest realistic model. The model remnant is evolving in a smooth ambient medium of fairly high density, about 6 cm^-3 on average, with a substantial density gradient. At the observed time it has an age of about 20,000 years, consistent with the age of the associated pulsar, and a radius of 11 to 13 pc. Over most of the outer surface, radiative cooling has become important in the post shock gas; on the denser end there has been sufficient compression of the cooled gas to develop a very thin dense half shell of about 450 M_sun, supported against further compression by nonthermal pressure. The half shell has an expansion velocity of about 150 km s^-1, and is bounded on the outer surface by a radiative shock with that speed. We provide several analytic tools for the assembly of models of this type. We review the early evolution and shell formation analyses and their generalizations to evolution in a density gradient. We also calculate the density and temperature that should be present in the hot interior of a remnant with thermal conduction. We supply the van der Laan mechanism in a particularly useful form for the calculation of radio continuum from radiative remnants. Finally, we demonstrate a simple technique for estimating the optical emission expected. These tools are employed to choose parameters of models which we then explore with our 1d and 2d hydrocodes, providing, respectively, the detailed x-ray spectra and dynamical characteristics.

We show that non-chaotic multiply-periodic particle orbits can exist in a galaxy-like potential where a small fast-tumbling nuclear bar is nested inside a main bar which has a slower pattern speed. We introduce the concept of a loop: a one-dimensional curve such that particles distributed along it at some initial instant return to the same curve (as viewed in the rotating frame of one of the bars) after the bars return to the same relative position. Just as particles following regular orbits in a simply-barred potential are trapped around closed periodic orbits, so regular orbits in a doubly-barred potential oscillate about stable loops. We find both loops which remain aligned with the inner bar and loops aligned with the outer bar: particles trapped around these loops could form the building blocks for a long-lived, self-consistent, doubly-barred galaxy. In a realistic doubly-barred galaxy potential, we demonstrate the existence of stable loops which support the structure of the inner bar. We use the epicyclic approximation to preview the pattern of gas flow in a doubly-barred potential.