Bar Major Axis Research Articles

Bar-driven Gas Dynamics of M31

The large-scale gaseous shocks in the bulge of M31 can be naturally explained by a rotating stellar bar. We use gas dynamical models to provide an independent measurement of the bar pattern speed in M31. The gravitational potentials of our simulations are from a set of made-to-measure models constrained by stellar photometry and kinematics. If the inclination of the gas disk is fixed at i = 77°, we find that a low pattern speed of 16–20 km s−1 kpc−1 is needed to match the observed position and amplitude of the shock features, as shock positions are too close to the bar major axis in high Ω b models. The pattern speed can increase to 20–30 km s−1 kpc−1 if the inner gas disk has a slightly smaller inclination angle compared with the outer one. Including subgrid physics such as star formation and stellar feedback has minor effects on the shock amplitude, and does not change the shock position significantly. If the inner gas disk is allowed to follow a varying inclination similar to the H i and ionized gas observations, the gas models with a pattern speed of 38 km s−1 kpc−1, which is consistent with stellar-dynamical models, can match both the shock features and the central gas features.

Unveiling a Hidden Bar-like Structure in NGC 1087: Kinematic and Photometric Evidence Using MUSE/VLT, ALMA, and JWST

We report a faint nonaxisymmetric structure in NGC 1087 through the use of James Webb Space Telescope Near Infrared Camera, with an associated kinematic counterpart observed as an oval distortion in the stellar velocity map, Hα, and CO J = 2 → 1 velocity fields. This structure is not evident in the MUSE optical continuum images but only revealed in the near-IR with the F200W and F300M band filters at 2 μm and 3 μm, respectively. Due to its elongation, this structure resembles a stellar bar although with remarkable differences with respect to conventional stellar bars. Most of the near-IR emission is concentrated within 6″∼500 pc with a maximum extension up to 1.2 kpc. The spatial extension of the large-scale noncircular motions is coincident with the bar, which undoubtedly confirms the presence of a nonaxisymmetric perturbation in the potential of NGC 1087. The oval distortion is enhanced in CO due to its dynamically cold nature rather than in Hα. We found that the kinematics in all phases, including stellar, ionized, and molecular, can be described simultaneously by a model containing a bisymmetric perturbation; however, we find that an inflow model of gas along the bar major axis is also likely. Furthermore, the molecular mass inflow rate associated can explain the observed star-formation rate in the bar. This reinforces the idea that bars are mechanisms for transporting gas and triggering star formation. This work contributes to our understanding of nonaxisymmetry in galaxies using the most sophisticated data so far.

Quantifying the influence of bars on action-based dynamical modelling of disc galaxies

ABSTRACT Action-based dynamical modelling, using stars as dynamical tracers, is an excellent diagnostic to estimate the underlying axisymmetric matter distribution of the Milky Way. However, the Milky Way’s bar causes non-axisymmetric resonance features in the stellar disc. Using RoadMapping (an action-based dynamical modelling framework to estimate the gravitational potential and the stellar distribution function), we systematically quantify the robustness of action-based modelling in the presence of a bar. We construct a set of test-particle simulations of barred galaxies (with varying bar properties), and apply RoadMapping to different survey volumes (with varying azimuthal position, size) drawn from these barred models. For realistic bar parameters, the global potential parameters are still recovered to within ∼1–17 per cent. However, with increasing bar strength, the best-fitting values of the parameters progressively deviate from their true values. This happens due to a combination of radial heating, radial migration, and resonance overlap phenomena in our bar models. Furthermore, the azimuthal location and the size of the survey volumes play important roles in the successful recovery of the parameters. Survey volumes along the bar major axis produce larger (relative) errors in the best-fitting parameter values. In addition, the potential parameters are better recovered for survey volumes with larger spatial coverage. As the Sun is located just ∼28°–33° behind the bar’s major axis, an estimate for the bar-induced systematic bias – as provided by this study – is therefore crucial for future modelling attempts of the Milky Way.

Photometric Signature of Ultraharmonic Resonances in Barred Galaxies

Bars may induce morphological features, such as rings, through their resonances. Previous studies suggested that the presence of “dark gaps,” or regions of a galaxy where the difference between the surface brightness along the bar major axis and that along the bar minor axis is maximal, can be attributed to the location of bar corotation. Here, using GALAKOS, a high-resolution N-body simulation of a barred galaxy, we test this photometric method’s ability to identify the bar corotation resonance. Contrary to previous work, our results indicate that “dark gaps” are a clear sign of the location of the 4:1 ultraharmonic resonance instead of bar corotation. Measurements of the bar corotation can indirectly be inferred using kinematic information, e.g., by measuring the shape of the rotation curve. We demonstrate our concept on a sample of 578 face-on barred galaxies with both imaging and integral field observations and find that the sample likely consists primarily of fast bars.

Local variations of the stellar velocity ellipsoid – II. The effect of the bar in the inner regions of Auriga galaxies

ABSTRACT Theoretical works have shown that off-plane motions of bars can heat stars in the vertical direction during buckling but is not clear how do they affect the rest of components of the stellar velocity ellipsoid (SVE). We study the 2D spatial distribution of the vertical, σz, azimuthal, σϕ, and radial, σr velocity dispersions in the inner regions of Auriga galaxies, a set of high-resolution magneto-hydrodynamical cosmological zoom-in simulations, to unveil the influence of the bar on the stellar kinematics. σz and σϕ maps exhibit non-axisymmetric features that closely match the bar light distribution with low-σ regions along the bar major axis and high values in the perpendicular direction. On the other hand, σr velocity dispersion maps present more axisymmetric distributions. We show that isophotal profile differences best capture the impact of the bar on the three SVE components providing strong correlations with bar morphology proxies although there is no relation with individual σ. Time evolution analysis shows that these differences are a consequence of the bar formation and that they tightly coevolve with the strength of the bar. We discuss the presence of different behaviours of σz and its connection with observations. This work helps us understand the intrinsic σ distribution and motivates the use of isophotal profiles as a mean to quantify the effect of bars.

Molecular gas and star formation within 12 strong galactic bars observed with IRAM-30 m

Context.While some galactic bars show recent massive star formation (SF) along them, some others do not. Whether bars with low level of SF are a consequence of low star formation efficiency, low gas inflow rate, or dynamical effects remains a matter of debate.Aims.In order to study the physical conditions that enable or prevent SF, we perform a multi-wavelength analysis of 12 strongly barred galaxies with total stellar masses log10(M⋆/M⊙)∈[10.2, 11], chosen to host different degrees of SF along the bar major axis without any prior condition on gas content. We observe the CO(1–0) and CO(2–1) emission within bars with the IRAM-30 m telescope (beam sizes of 1.7–3.9 kpc and 0.9–2.0 kpc, respectively; 7–8 pointings per galaxy on average).Methods.We estimated molecular gas masses (Mmol) from the CO(1–0) and CO(2–1) emissions. SF rates (SFRs) were calculated from GALEX near-ultraviolet (UV) and WISE 12 μm images within the beam-pointings, covering the full bar extent (SFRs were also derived from far-UV and 22 μm).Results.We detect molecular gas along the bars of all probed galaxies. Molecular gas and SFR surface densities span the ranges log10(Σmol/[M⊙pc−2]) ∈ [0.4,2.4] and log10(ΣSFR/[M⊙pc−1kpc−2]]) ∈ [−3.25, −0.75], respectively. The star formation efficiency (SFE; i.e., SFR/Mmol) in bars varies between galaxies by up to an order of magnitude (SFE ∈[0.1, 1.8] Gyr−1). On average, SFEs are roughly constant along bars. SFEs are not significantly different from the mean value in spiral galaxies reported in the literature (∼0.43 Gyr−1), regardless of whether we estimateMmolfrom CO(1–0) or CO(2–1). Interestingly, the higher the total stellar mass of the host galaxy, the lower the SFE within their bars. In particular, the two galaxies in our sample with the lowest SFE and ΣSFR(NGC 4548 and NGC 5850, SFE ≲ 0.25 Gyr−1, ΣSFR ≲ 10−2.25 M⊙yr−1kpc−2,M⋆ ≳ 1010.7 M⊙) are also those hosting massive bulges and signs of past interactions with nearby companions.Conclusions.We present a statistical analysis of the SFE in bars for a sample of 12 galaxies. The SFE in strong bars is not systematically inhibited (either in the central, middle, or end parts of the bar). Both environmental and internal quenching are likely responsible for the lowest SFEs reported in this work.

Galaxies within galaxies in the TIMER survey: stellar populations of inner bars are scaled replicas of main bars

Inner bars are frequent structures in the local Universe and thought to substantially influence the nuclear regions of disc galaxies. In this study we explore the structure and dynamics of inner bars by deriving maps and radial profiles of their mean stellar population content and comparing them to previous findings in the context of main bars. To this end, we exploit observations obtained with the integral-field spectrograph MUSE of three double-barred galaxies in the TIMER sample. The results indicate that inner bars can be clearly distinguished based on their stellar population properties alone. More precisely, inner bars show significantly elevated metallicities and depleted [α/Fe] abundances. Although they exhibit slightly younger stellar ages compared to the nuclear disc, the typical age differences are small, except at their outer ends. These ends of the inner bars are clearly younger compared to their inner parts, an effect known from main bars as orbital age separation. In particular, the youngest stars (i.e. those with the lowest radial velocity dispersion) seem to occupy the most elongated orbits along the (inner) bar major axis. We speculate that these distinct ends of bars could be connected to the morphological feature of ansae. Radial profiles of metallicity and [α/Fe] enhancements are flat along the inner bar major axis, but show significantly steeper slopes along the minor axis. This radial mixing in the inner bar is also known from main bars and indicates that inner bars significantly affect the radial distribution of stars. In summary, based on maps and radial profiles of the mean stellar population content and in line with previous TIMER results, inner bars appear to be scaled down versions of the main bars seen in galaxies. This suggests the picture of a ‘galaxy within a galaxy’, with inner bars in nuclear discs being dynamically equivalent to main bars in main galaxy discs.

Stellar populations across galaxy bars in the MUSE TIMER project

Stellar populations in barred galaxies save an imprint of the influence of the bar on the host galaxy’s evolution. We present a detailed analysis of star formation histories (SFHs) and chemical enrichment of stellar populations in nine nearby barred galaxies from the TIMER project. We used integral field observations with the MUSE instrument to derive unprecedented spatially resolved maps of stellar ages, metallicities, [Mg/Fe] abundances, and SFHs, as well as Hαas a tracer of ongoing star formation. We find a characteristic V-shaped signature in the SFH that is perpendicular to the bar major axis, which supports the scenario where intermediate-age stars (∼2 − 6 Gyr) are trapped on more elongated orbits shaping a thinner part of the bar, while older stars (> 8 Gyr) are trapped on less elongated orbits shaping a rounder and thicker part of the bar. We compare our data to state-of-the-art cosmological magneto-hydrodynamical simulations of barred galaxies and show that such V-shaped SFHs arise naturally due to the dynamical influence of the bar on stellar populations with different ages and kinematic properties. Additionally, we find an excess of very young stars (< 2 Gyr) on the edges of the bars, predominantly on the leading side, thus confirming typical star formation patterns in bars. Furthermore, mass-weighted age and metallicity gradients are slightly shallower along the bar than in the disc, which is likely due to orbital mixing in the bar. Finally, we find that bars are mostly more metal-rich and less [Mg/Fe]-enhanced than the surrounding discs. We interpret this as a signature that the bar quenches star formation in the inner region of discs, usually referred to as star formation deserts. We discuss these results and their implications on two different scenarios of bar formation and evolution.

High Performance Computing of Magnetized Galactic Disks

A parallel implementation of the magneto-hydrodynamical code for global modeling of the galactic evolution is reported. The code is parallelized by using MPI interface, and it shows ideal scaling up to 200–300 cores on Lomonosov supercomputer with fast interconnect. In the benchmarking of this code, we study the dynamics of a magnetized gaseous disk of a galaxy with a bar. We run a high-resolution 3D magnetohydrodynamic simulation taking into account the Milky Way-like gravitational potential, gas self-gravity and a network of cooling and heating processes in the interstellar medium. By using this simulation the evolution of morphology and enhancement of the magnetic field are explored. In agreement to hydrodynamical models, when the bar is strong enough, the gas develops sharp shocks at the leading side of the bar. In such a picture we found that when typically the magnetic field strength traces the location of the largescale shocks along the bar major axis, the magnetic field pressure weakens the shocks and reduces the inflow of gas towards the galactic center.

Models of bars - I. Flattish profiles for early-type galaxies

We introduce a simple family of barred galaxy models with flat rotation curves. They are built by convolving the axisymmetric logarithmic model with a needle density. The density contours in the bar region are highly triaxial and elongated, but become spheroidal in the outer parts. The mass distribution differs markedly from the elliptical shape assumed in other analytical models, like Ferrers or Freeman bars. The surface density profile along the bar major axis is flattish, as befits models for bars in early-type galaxies (SB0, SBa). The two-dimensional orbital structure of the models is analysed with surfaces of section and characteristic diagrams and it reveals qualitatively new features. For some pattern speeds, additional Lagrange points can exist along the major axis, and give rise to off-centered, trapped orbits. When the bar is weak, the orbital structure is very simple, comprising just prograde, aligned $x_1$ orbits and retrograde anti-aligned $x_4$ orbits. As the bar strength increases, the $x_1$ family becomes unstable and vanishes, with propeller orbits dominating the characteristic diagram.

The BaLROG project – II. Quantifying the influence of bars on the stellar populations of nearby galaxies

We continue the exploration of the BaLROG (Bars in Low Redshift Optical Galaxies) sample: 16 large mosaics of barred galaxies observed with the integral field unit SAURON. We quantify the influence of bars on the composition of the stellar component. We derive linestrength indices of H${\beta}$, Fe5015 and Mgb. Based on single stellar population (SSP) models, we calculate ages, metallicities and [Mg/Fe] abundances and their gradients along the bar major and minor axes. The high spatial resolution of our data allows us to identify breaks among index and SSP profiles, commonly at 0.13$\pm$0.06 bar length, consistent with kinematic features. Inner gradients are about ten times steeper than outer gradients and become larger when there is a central rotating component, implying that the gradients are not independent of dynamics and orbits. Central ages appear to be younger for stronger bars. Yet, the bar regions are usually old. We find a flattening of the iron (Fe5015) and magnesium (Mgb) outer gradients along the bar major axis, translating into a flattening of the metallicity gradient. This gradient is found to be 0.03$\pm$0.07 dex/kpc along the bar major axis while the mean value of the bar minor axis compares well with that of an unbarred control sample and is significantly steeper, namely -0.20$\pm$0.04 dex/kpc. These results confirm recent simulations and discern the important localized influence of bars. The elevated [Mg/Fe] abundances of bars and bulges compared to the lower values of discs suggest an early formation, in particular for early type galaxies.

Shepherding tidal debris with the Galactic bar: the Ophiuchus stream

The dynamics of stellar streams in rotating barred potentials is explained for the first time. Naturally, neighbouring stream stars reach pericentre at slightly different times. In the presence of a rotating bar, these neighbouring stream stars experience different bar orientations during pericentric passage and hence each star receives a different torque from the bar. These differing torques reshape the angular momentum and energy distribution of stars in the stream, which in turn changes the growth rate of the stream. For a progenitor orbiting in the same sense as the bar's rotation and satisfying a resonance condition, the resultant stream can be substantially shorter or longer than expected, depending on whether the pericentric passages of the progenitor occur along the bar's minor or major axis respectively. We present a full discussion of this phenomenon focusing mainly on streams confined to the Galactic plane. In stark contrast with the evolution in static potentials, which give rise to streams that grow steadily in time, rotating barred potentials can produce dynamically old, short streams. This challenges the traditional viewpoint that the inner halo consists of well phase-mixed material whilst the tidally-disrupted structures in the outer halo are more spatially coherent. We argue that this mechanism plays an important role in explaining the mysteriously short Ophiuchus stream that was recently discovered near the bulge region of the Milky Way.

On the nature of the barlens component in barred galaxies: what do boxy/peanut bulges look like when viewed face-on?

Barred galaxies have interesting morphological features whose presence and properties set constraints on galactic evolution. Here we examine barlenses, i.e. lens-like components whose extent along the bar major axis is shorter than that of the bar and whose outline is oval or circular. We identify and analyse barlenses in $N$-body plus SPH simulations, compare them extensively with those from the NIRS0S (Near-IR S0 galaxy survey) and the S$^4$G samples (Spitzer Survey of Stellar Structure in Galaxies) and find very good agreement. We observe barlenses in our simulations from different viewing angles. This reveals that barlenses are the vertically thick part of the bar seen face-on, i.e. a barlens seen edge-on is a boxy/peanut/X bulge. In morphological studies, and in the absence of kinematics or photometry, a barlens, or part of it, may be mistaken for a classical bulge. Thus the true importance of classical bulges, both in numbers and mass, is smaller than currently assumed, which has implications for galaxy formation studies. Finally, using the shape of the isodensity curves, we propose a rule of thumb for measuring the barlens extent along the bar major axis of moderately inclined galaxies, thus providing an estimate of which part of the bar is thicker.

The X-shaped Milky Way bulge in OGLE-III★ photometry and in N-body models

We model the split red clump of the Galactic bulge in OGLE-III photometry, and compare the results to predictions from two N-body models. Our analysis yields precise maps of the brightness of the two red clumps, the fraction of stars in the more distant peak, and their combined surface density. We compare the observations to predictions from two N-body models previously used in the literature. Both models correctly predict several features as long as one assumes an angle $\alpha_{\rm{Bar}} \approx 30^{\circ}$ between the Galactic bar's major axis and the line of sight to the Galactic centre. In particular that the fraction of stars in the faint red clump should decrease with increasing longitude. The biggest discrepancies between models and data are in the rate of decline of the combined surface density of red clump stars toward negative longitudes and of the brightness difference between the two red clumps toward positive longitudes, with neither discrepancy exceeding $\sim$25% in amplitude. Our analysis of the red giant luminosity function also yields an estimate of the red giant branch bump parameters toward these high-latitude fields, and evidence for a high rate ($\sim$25%) of disk contamination in the bulge at the colour and magnitude of the red clump, with the disk contamination rate increasing toward sightlines further distant from the plane.

Metallicity and kinematics of the bar in situ

Constraints on the Galactic bulge/bar structure and formation history from stellar kinematics and metallicities mainly come from relatively high-latitude fields (|b|>4) where a complex mix of stellar population is seen. We aim here to constrain the formation history of the Galactic bar by studying the radial velocity and metallicity distributions of stars in-situ (|b|<1). We observed red clump stars in four fields along the bar's major axis (l=10,6,-6 and b=0 plus a field at l=0,b=1) with low-resolution spectroscopy from VLT/FLAMES, observing around the CaII triplet. We developed robust methods for extracting radial velocity and metallicity estimates from these low signal-to-noise spectra. We derived distance probability distributions using Bayesian methods rigorously handling the extinction law. We present radial velocities and metallicity distributions, as well as radial velocity trends with distance. We observe an increase in the radial velocity dispersion near the Galactic plane. We detect the streaming motion of the stars induced by the bar in fields at l=+/-6, the highest velocity components of this bar stream being metal-rich ([Fe/H]~0.2 dex). Our data is consistent with a bar inclined at 26+/-3 from the Sun-Galactic centre line. We observe a significant fraction of metal-poor stars, in particular in the field at l=0,b=1. We confirm the flattening of the metallicity gradient along the minor axis when getting closer to the plane, with a hint that it could actually be inverted. Our stellar kinematics corresponds to the expected behaviour of a bar issued from the secular evolution of the Galactic disc. The mix of several populations, seen further away from the plane, is also seen in the bar in-situ since our metallicity distributions highlight a different spatial distribution between metal-poor and metal-rich stars, the more metal-poor stars being more centrally concentrated.

A new photometric model of the Galactic bar using red clump giants

We present a study of the luminosity density distribution of the Galactic bar using number counts of red clump giants (RCGs) from the OGLE-III survey. The data were recently published by Nataf et al. (2013) for 9019 fields towards the bulge and have $2.94\times 10^6$ RC stars over a viewing area of $90.25 \,\textrm{deg}^2$. The data include the number counts, mean distance modulus ($\mu$), dispersion in $\mu$ and full error matrix, from which we fit the data with several tri-axial parametric models. We use the Markov Chain Monte Carlo (MCMC) method to explore the parameter space and find that the best-fit model is the $E_3$ model, with the distance to the GC is 8.13 kpc, the ratio of semi-major and semi-minor bar axis scale lengths in the Galactic plane $x_{0},y_{0}$, and vertical bar scale length $z_0$, is $x_0:y_0:z_0 \approx 1.00:0.43:0.40$ (close to being prolate). The scale length of the stellar density profile along the bar's major axis is $\sim$ 0.67 kpc and has an angle of $29.4^\circ$, slightly larger than the value obtained from a similar study based on OGLE-II data. The number of estimated RC stars within the field of view is $2.78 \times 10^6$, which is systematically lower than the observed value. We subtract the smooth parametric model from the observed counts and find that the residuals are consistent with the presence of an X-shaped structure in the Galactic centre, the excess to the estimated mass content is $\sim 5.8%$. We estimate the total mass of the bar is $\sim 1.8 \times 10^{10} M_\odot$. Our results can be used as a key ingredient to construct new density models of the Milky Way and will have implications on the predictions of the optical depth to gravitational microlensing and the patterns of hydrodynamical gas flow in the Milky Way.

REDDENING AND EXTINCTION TOWARD THE GALACTIC BULGE FROM OGLE-III: THE INNER MILKY WAY'SRV∼ 2.5 EXTINCTION CURVE

We combine VI photometry from OGLE-III with VVV and 2MASS measurements of E(J-K_{s}) to resolve the longstanding problem of the non-standard optical extinction toward the Galactic bulge. We show that the extinction is well-fit by the relation A_{I} = 0.7465*E(V-I) + 1.3700*E(J-K_{s}), or, equivalently, A_{I} = 1.217*E(V-I)(1+1.126*(E(J-K_{s})/E(V-I)-0.3433)). The optical and near-IR reddening law toward the inner Galaxy approximately follows an R_{V} \approx 2.5 extinction curve with a dispersion {\sigma}_{R_{V}} \approx 0.2, consistent with extragalactic investigations of the hosts of type Ia SNe. Differential reddening is shown to be significant on scales as small as as our mean field size of 6', with the 1{\sigma} dispersion in reddening averaging 9% of total reddening for our fields. The intrinsic luminosity parameters of the Galactic bulge red clump (RC) are derived to be (M_{I,RC}, \sigma_{I,RC,0}, (V-I)_{RC,0}, \sigma_{(V-I)_{RC}}, (J-K_{s})_{RC,0}) = (-0.12, 0.09, 1.06, 0.121, 0.66). Our measurements of the RC brightness, brightness dispersion and number counts allow us to estimate several Galactic bulge structural parameters. We estimate a distance to the Galactic center of 8.20 kpc, resolving previous discrepancies in distance determinations to the bulge based on I-band observations. We measure an upper bound on the tilt {\alpha} \approx 40{\deg}. between the bar's major axis and the Sun-Galactic center line of sight, though our brightness peaks are consistent with predictions of an N-body model oriented at {\alpha} \approx 25{\deg}. The number of RC stars suggests a total stellar mass for the Galactic bulge of 2.0*10^{10} M_{\odot}, if one assumes a Salpeter IMF.

Peanuts at an angle: detecting and measuring the three-dimensional structure of bars in moderately inclined galaxies

We show that direct detection and measurement of the vertically thickened parts of bars (so-called ‘boxy’ or ‘peanut-shaped’ bulges) are possible not only for edge-on galaxies but also for galaxies with moderate inclinations (i < 70°), and that examples are relatively common in the nearby Universe. The analysis of a sample of 78 nearby, moderately inclined (i ≲ 65°) early-type (S0–Sb) barred galaxies shows that the isophotal signature of the box/peanut can usually be detected for inclinations as low as i ∼ 40° – and in exceptional cases down to i ∼ 30°. In agreement with the predictions from N-body simulations, the signature is most easily detectable when the bar's position angle is within ∼50° of the galaxy major axis; in particular, galaxies where the bar lies very close to the minor axis do not show the signature clearly or at all. For galaxies with i = 40°–65° and relative angles <45°, we find evidence for the signature ≈2/3 of the time; the true frequency of box/peanut structures in bars may be higher. Comparison with N-body models also allows us to link observed photometric morphology with 3D physical structures, and thus estimate the relative sizes of box/peanut structures and bars. For our local sample, we find that box/peanut structures range in radial size (measured along the bar major axis) from 0.4 to 3.8 kpc (mean =1.5 ± 0.9 kpc) and span 0.26–0.58 of the bar length (mean of 0.38 ± 0.08). This is a clear observational confirmation that when bars thicken, it is not the entire bar which does so, but only the inner part. This technique can also be used to identify galaxies with bars which have not vertically thickened. We suggest that NGC 3049 and IC 676 may be particularly good examples, and that the fraction of S0–Sb bars which lack box/peanut structures is at least ∼13 per cent.

CENTRAL REGIONS OF BARRED GALAXIES: TWO-DIMENSIONAL NON-SELF-GRAVITATING HYDRODYNAMIC SIMULATIONS

The inner regions of barred galaxies contain substructures such as off-axis shocks, nuclear rings, and nuclear spirals. These substructure may affect star formation, and control the activity of a central black hole (BH) by determining the mass inflow rate. We investigate the formation and properties of such substructures using high-resolution, grid-based hydrodynamic simulations. The gaseous medium is assumed to be infinitesimally-thin, isothermal, and non-self-gravitating. The stars and dark matter are represented by a static gravitational potential with four components: a stellar disk, the bulge, a central BH, and the bar. To investigate various galactic environments, we vary the gas sound speed c_s as well as the mass of the central BH M_BH. Once the flow has reached a quasi-steady state, off-axis shocks tend to move closer to the bar major axis as c_s increases. Nuclear rings shrink in size with increasing c_s, but are independent of M_BH, suggesting that ring position is not determined by the Lindblad resonances. Rings in low-c_s models are narrow since they are occupied largely by gas on x2-orbits and well decoupled from nuclear spirals, while they become broad because of large thermal perturbations in high-c_s models. Nuclear spirals persist only when either c_s is small or M_BH is large; they would otherwise be destroyed completely by the ring material on eccentric orbits. The shape and strength of nuclear spirals depend sensitively on c_s and M_BH such that they are leading if both c_s and M_BH are small, weak trailing if c_s is small and M_BH is large, and strong trailing if both c_s and M_BH are large. While the mass inflow rate toward the nucleus is quite small in low-c_s models because of the presence of a narrow nuclear ring, it becomes larger than 0.01 Msun/yr when c_s is large, providing a potential explanation of nuclear activity in Seyfert galaxies.

Rings and spirals in barred galaxies – II. Ring and spiral morphology

In this series of papers, we propose a theory to explain the formation and properties of rings and spirals in barred galaxies. The building blocks of these structures are orbits guided by the manifolds emanating from the unstable Lagrangian points located near the ends of the bar. In this paper we focus on a comparison of the morphology of observed and of theoretical spirals and rings and we also give some predictions for further comparisons. Our theory can account for spirals as well as both inner and outer rings. The model outer rings have the observed $R_1$, $R_1'$, $R_2$, $R_2'$ and $R_1R_2$ morphologies, including the dimples near the direction of the bar major axis. We explain why the vast majority of spirals in barred galaxies are two armed and trailing, and discuss what it would take for higher multiplicity arms to form. We show that the shapes of observed and theoretical spirals agree and we predict that stronger non-axisymmetric forcings at and somewhat beyond corotation will drive more open spirals. We compare the ratio of ring diameters in theory and in observations and predict that more elliptical rings will correspond to stronger forcings. We find that the model potential may influence strongly the numerical values of these ratios.