Barred Spiral Research Articles

Enhanced Star Formation in Both Disks and Ram-pressure-stripped Tails of GASP Jellyfish Galaxies

Abstract Exploiting the data from the GAs Stripping Phenomena in galaxies with MUSE (GASP) program, we compare the integrated star formation rate–mass relation (SFR–M *) relation of 42 cluster galaxies undergoing ram-pressure stripping (RPS; “stripping galaxies”) to that of 32 field and cluster undisturbed galaxies. Theoretical predictions have so far led to contradictory conclusions about whether or not ram pressure can enhance the star formation (SF) in the gas disks and tails, and until now a statistically significant observed sample of stripping galaxies was lacking. We find that stripping galaxies occupy the upper envelope of the control sample SFR–M * relation, showing a systematic enhancement of the SFR at any given mass. The star formation enhancement occurs in the disk (0.2 dex), and additional SF takes place in the tails. Our results suggest that strong RPS events can moderately enhance the SF also in the disk prior to gas removal.

Time Inference with MUSE in Extragalactic Rings (TIMER): properties of the survey and high-level data products

The Time Inference with MUSE in Extragalactic Rings (TIMER) project is a survey with the VLT-MUSE integral-field spectrograph of 24 nearby barred galaxies with prominent central structures (e.g., nuclear rings or inner discs). The main goals of the project are: (i) estimating the cosmic epoch when discs of galaxies settle, leading to the formation of bars; (ii) testing the hypothesis whereby discs in more massive galaxies are assembled first; and (iii) characterising the history of external gas accretion in disc galaxies. We present details on the sample selection, observations, data reduction, and derivation of high-level data products, including stellar kinematics, ages and metallicities. We also derive star formation histories and physical properties and kinematics of ionised gas. We illustrate how this dataset can be used for a plethora of scientific applications, e.g., stellar feedback, outflows, nuclear and primary bars, stellar migration and chemical enrichment, and the gaseous and stellar dynamics of nuclear spiral arms, barlenses, box/peanuts and bulges. Amongst our first results - based on a few selected galaxies -, we show that the dynamics of nuclear rings and inner discs is consistent with the picture in which they are formed by bars, that the central few hundred parsecs in massive disc galaxies tend to show a pronounced peak in stellar metallicity, and that nuclear rings can efficiently prevent star formation in this region. Finally, we present evidence that star-bursting nuclear rings can be fed with low-metallicity gas from low-mass companions.

Dust in the eye of Andromeda

We present new Herschel-derived images of warm dust in the Andromeda Galaxy, M31, with unprecedented spatial resolution (~ 30 pc), column density accuracy, and constraints on the three-dimensional distributions of dust temperature and dust opacity index (hence grain size and composition), based on the new PPMAP Bayesian analysis procedure. We confirm the overall radial variation of dust opacity index reported by other recent studies, including the central decrease within ~ 3 kpc of the nucleus. We also investigate the detailed distribution of dust in the nuclear region, a prominent feature of which is a ~ 500 pc bar-like structure seen previously in H{\alpha}. The nature of this feature has been the subject of some debate. Our maps show it to be the site of the warmest dust, with a mean line-of-sight temperature ~ 30 K. A comparison with the stellar distribution, based on 2MASS data, provides strong evidence that it is a gravitationally induced bar. A comparison with radial velocity maps suggests the presence of an inflow towards the nucleus from opposite directions along this bar, fed presumably by the nuclear spiral with which it appears to connect. Such behaviour is common in large-scale bars in spiral galaxies, as is the phenomenon of nested bars whereby a subkiloparsec nuclear bar exists within a large-scale primary bar. We suggest that M31 represents an example of such nesting.

ALMA 12CO (J = 1–0) imaging of the nearby galaxy M 83: Variations in the efficiency of star formation in giant molecular clouds

Abstract We present results of the 12CO (1–0) mosaic observations of the nearby barred-spiral galaxy M 83 obtained with the Atacama Large Millimeter/submillimeter Array (ALMA). The total flux is recovered by combining the ALMA data with single-dish data obtained using the Nobeyama 45 m telescope. The combined map covers a ∼13 kpc2 field that includes the galactic center, eastern bar, and spiral arm with a resolution of 2${^{\prime\prime}_{.}}$03 × 1${^{\prime\prime}_{.}}$15 (∼45 pc × ∼25 pc). With a resolution comparable to typical sizes of giant molecular clouds (GMCs), the CO distribution in the bar and arm is resolved into many clumpy peaks that form ridge-like structures. Remarkably, in the eastern arm, the CO peaks form two arc-shaped ridges that run along the arm and exhibit a distinct difference in the activity of star formation: the one on the leading side has numerous H ii regions associated with it, whereas the other one on the trailing side has only a few. To see whether GMCs form stars with uniform star formation efficiency (SFE) per free-fall time (SFEff), GMCs are identified from the data cube and then cross-matched with the catalog of ${H}\,\, \small {II}$ regions to estimate the star formation rate for each of them. 179 GMCs with a median mass of 1.6 × 106 M⊙ are identified. The mass-weighted average SFEff of the GMCs is ∼9.4 × 10−3, which is in agreement with models of turbulence regulated star formation. Meanwhile, we find that SFEff is not universal within the mapped region. In particular, one of the arm ridges shows a high SFEff with a mass-weighted value of ∼2.7 × 10−2, which is higher by more than a factor of 5 compared to the inter-arm regions. This large regional variation in SFEff favors the recent interpretation that GMCs do not form stars at a constant rate within their lifetime.

Ram pressure stripping made easy: an analytical approach

The removal of gas by ram pressure stripping of galaxies is treated by a purely kinematic description. The solution has two asymptotic limits: if the duration of the ram pressure pulse exceeds the period of vertical oscillations perpendicular to the galactic plane, the commonly used quasi-static criterion of Gunn & Gott is obtained which uses the maximum ram pressure that the galaxy has experienced along its orbit. For shorter pulses the outcome depends on the time-integrated ram pressure. This parameter pair fully describes the gas mass fraction that is stripped from a given galaxy. This approach closely reproduces results from SPH simulations. We show that typical galaxies follow a very tight relation in this parameter space corresponding to a pressure pulse length of about 300 Myr. Thus, the Gunn & Gott criterion provides a good description for galaxies in larger clusters. Applying the analytic description to a sample of 232 Virgo galaxies from the GoldMine database, we show that the ICM provides indeed the ram pressures needed to explain the deficiencies. We also can distinguish current and past strippers, including objects whose stripping state was unknown.

Multi-scale simulations of black hole accretion in barred galaxies

Due to the non-axisymmetric potential of the central bar, in addition to their characteristic arms and bar, barred spiral galaxies form a variety of structures within the thin gas disk, such as nuclear rings, inner spirals, and dust lanes. These structures in the inner kiloparsec are extremely important in order to explain and understand the rate of black hole feeding. The aim of this work is to investigate the influence of stellar bars in spiral galaxies on the thin self-gravitating gas disk. We focus on the accretion of gas onto the central supermassive black hole and its time-dependent evolution. We conducted multi-scale simulations simultaneously resolving the galactic disk and the accretion disk around the central black hole. In all the simulations we varied the initial gas disk mass. As an additional parameter we chose either the gas temperature for isothermal simulations or the cooling timescale for non-isothermal simulations. Accretion was either driven by a gravitationally unstable or clumpy accretion disk or by energy dissipation in strong shocks. Most of the simulations show a strong dependence of the accretion rate at the outer boundary of the central accretion disk (r< 300 pc) on the gas flow at kiloparsec scales. The final black hole masses reach up to ~109 M⊙after 1.6 Gyr. Our models show the expected influence of the Eddington limit and a decline in growth rate at the corresponding sub-Eddington limit.

Chemical enrichment in isolated barred spiral galaxies

To investigate the role of bars in the chemical evolution of isolated disc galaxies, we performed a series of 39 gas dynamical simulations of isolated barred and unbarred galaxies with various masses, initial gas fractions, and AGN feedback models. The presence of a bar drives a substantial amount of gas toward the central region of the galaxy. In the most massive galaxies, this results in a violent starburst, followed by a drop in star formation resulting from gas exhaustion. The time delay between Type Ia and Type II supernovae explosions means that barred galaxies experience a rapid increase in [O/H] in the central region, and a much more gradual increase in [Fe/H]. In unbarred galaxies, star formation proceeds at a slow and steady rate, and oxygen and iron are produced at steady rates which are similar except for a time offset. Comparing the abundance ratios in barred and unbarred galaxies with the same central stellar mass $M_*$, we find in barred galaxies an enhancement of 0.07 dex in [O/H], 0.05 dex in [Fe/H], and 0.05 dex in [O/Fe]. The [O/H] enhancement is in excellent agreement with observations from the SDSS. The initial gas fraction has very little effect on the abundance ratios in barred and unbarred galaxies, unless the galaxies experience a starburst. We considered AGN-host galaxies located near the bottom of the AGN regime, $M_*>3\times10^{10}M_\odot$, where AGN feedback dominates over supernovae feedback. We found that the impact of AGN feedback on the central abundances is marginal.

Collisionless Boltzmann equation approach for the study of stellar discs within barred galaxies

We have studied the kinematics of stellar disc populations within the solar neighbourhood in order to find the imprints of the Galactic bar. We carried out the analysis by developing a numerical resolution of the 2D2V (two-dimensional in the physical space, 2D, and two-dimensional in the velocity motion, 2V) collisionless Boltzmann equation and modelling the stellar motions within the plane of the Galaxy within the solar neighbourhood. We recover similar results to those obtained by other authors using N-body simulations, but we are also able to numerically identify faint structures thanks to the cancelling of the Poisson noise. We find that the ratio of the bar pattern speed to the local circular frequency is in the range ΩB∕Ω = 1.77 to 1.91. If the Galactic bar angle orientation is within the range from 24 to 45 degrees, the bar pattern speed is between 46 and 49 km s−1 kpc−1.

Detection of bars in galaxies using a deep convolutional neural network

We present an automated method for the detection of bar structure in optical images of galaxies using a deep convolutional neural network which is easy to use and provides good accuracy. In our study we use a sample of 9346 galaxies in the redshift range 0.009-0.2 from the Sloan Digital Sky Survey, which has 3864 barred galaxies, the rest being unbarred. We reach a top precision of ~94 per cent in identifying bars in galaxies using the trained network. This accuracy matches the accuracy reached by human experts on the same data without additional information about the images. Since Deep Convolutional Neural Networks can be scaled to handle large volumes of data, the method is expected to have great relevance in an era where astronomy data is rapidly increasing in terms of volume, variety, volatility and velocity along with other V's that characterize big data. With the trained model we have constructed a catalogue of barred galaxies from SDSS and made it available online.

Orbits in elementary, power-law galaxy bars – 1. Occurrence and role of single loops

Orbits in galaxy bars are generally complex, but simple closed loop orbits play an important role in our conceptual understanding of bars. Such orbits are found in some well-studied potentials, provide a simple model of the bar in themselves, and may generate complex orbit families. The precessing, power ellipse (p-ellipse) orbit approximation provides accurate analytic orbit orbits in symmetric galaxy potentials. It remains useful for finding and fitting simple loop orbits in the frame of a rotating bar with bar-like and symmetric power-law potentials. Second order perturbation theory yields two or fewer simple loop solutions in these potentials. Numerical integrations in the parameter space neighborhood of perturbation solutions reveal zero or one ac- tual loops in a range of such potentials with rising rotation curves. These loops are embedded in a small parameter region of similar, but librating orbits, which have a subharmonic frequency superimposed on the basic loop. These loops and their librat- ing companions support annular bars. Solid bars can be produced in more complex potentials, as shown by an example with power-law indices varying with radius. The power-law potentials can be viewed as the elementary constituents of more complex potentials. Numerical integrations also reveal interesting classes of orbits with multiple loops. In two-dimensional, self-gravitating bars, with power-law potentials, single loop orbits are very rare. This result suggests that gas bars or oval distortions are unlikely to be long-lived, and that complex orbits or three-dimensional structure must support self-gravitating stellar bars.

Molecular Gas Properties in M83 from CO PDFs

Abstract We have obtained 12CO(1–0) data of the nearby barred spiral galaxy M83 from Atacama Large Millimeter/submillimeter Array and Nobeyama 45 m observations. By combining these two data sets, the total CO flux has been recovered, and a high angular resolution ( corresponding to ∼40 pc at the distance of M83) has been achieved. The field of view is corresponding to ∼3.4 kpc and covers the galactic center, bar, and spiral arm regions. In order to investigate how these galactic structures affect gas properties, we have created a probability distribution function (PDF) of the CO integrated intensity ( ), peak temperature, and velocity dispersion for a region with each structure. We find that the PDF for the bar shows a bright-end tail while that for the arm does not. As the star formation efficiency is lower in the bar, this difference in PDF shape is contrary to the trend in Milky Way studies where the bright-end tail is found for star-forming molecular clouds. While the peak temperature PDFs are similar for the bar and arm regions, velocity dispersion in the bar is systematically larger than in the arm. This large velocity dispersion is likely a major cause of the bright-end tail and of suppressed star formation. We also investigate an effect of stellar feedback to PDF profiles and find that the different PDFs between bar and arm regions cannot be explained by the feedback effect, at least at the current spatial scale.

The dependence of bar frequency on galaxy mass, colour, and gas content – and angular resolution – in the local universe

I use distance- and mass-limited subsamples of the Spitzer Survey of Stellar Structure in Galaxies (S4G) to investigate how the presence of bars in spiral galaxies depends on mass, colour, and gas content and whether large, SDSS-based investigations of bar frequencies agree with local data. Bar frequency reaches a maximum of $\approx 0.70$ at $M_{\star} \sim 10^{9.7} M_{\odot}$, declining to both lower and higher masses. It is roughly constant over a wide range of colours ($g - r \approx 0.1$-0.8) and atomic gas fractions ($\log (M_{HI} / M_{\star}) \approx -2.5$ to 1). Bars are thus as common in blue, gas-rich galaxies are they are in red, gas-poor galaxies. This is in sharp contrast to many SDSS-based studies of $z \sim 0.01$-0.1 galaxies, which report f_bar increasing strongly to higher masses (from $M_{\star} \sim 10^{10}$ to $10^{11} M_{\odot}$), redder colours, and lower gas fractions. The contradiction can be explained if SDSS-based studies preferentially miss bars in, and underestimate the bar fraction for, lower-mass (bluer, gas-rich) galaxies due to poor spatial resolution and the correlation between bar size and stellar mass. Simulations of SDSS-style observations using the S4G galaxies as a parent sample, and assuming that bars below a threshold angular size of twice the PSF FWHM cannot be identified, successfully reproduce typical SDSS f_bar trends for stellar mass and gas mass ratio. Similar considerations may affect high-redshift studies, especially if bars grow in length over cosmic time; simulations suggest that high-redshift bar fractions may thus be systematically underestimated.

Star formation suppression and bar ages in nearby barred galaxies

We present new spectroscopic data for 21 barred spiral galaxies, which we use to explore the effect of bars on disk star formation, and to place constraints on the characteristic lifetimes of bar episodes. The analysis centres on regions of heavily suppressed star formation activity, which we term 'star formation deserts'. Long-slit optical spectroscopy is used to determine H beta absorption strengths in these desert regions, and comparisons with theoretical stellar population models are used to determine the time since the last significant star formation activity, and hence the ages of the bars. We find typical ages of approx. 1 Gyr, but with a broad range, much larger than would be expected from measurement errors alone, extending from about 0.25 Gyr to more than 4 Gyr. Low-level residual star formation, or mixing of stars from outside the 'desert' regions, could result in a doubling of these age estimates. The relatively young ages of the underlying populations coupled with the strong limits on the current star formation rule out a gradual exponential decline in activity, and hence support our assumption of an abrupt truncation event.

Detection of Enhanced Central Mass-to-light Ratios in Low-mass Early-type Galaxies: Evidence for Black Holes?

Abstract We present dynamical measurements of the central mass-to-light ratio (M/L) of a sample of 27 low-mass early-type galaxies. We consider all galaxies with 9.7 < log( < 10.5 in our analysis, selecting out galaxies with available high-resolution Hubble Space Telescope (HST) data, and eliminating galaxies with significant central color gradients or obvious dust features. We use the HST images to derive mass models for these galaxies and combine these with the central velocity dispersion values from data to obtain a central dynamical M/L estimate. These central dynamical are higher than dynamical derived at larger radii and stellar population estimates of the galaxy centers in ∼80% of galaxies, with a median enhancement of ∼14% and a statistical significance of 3.3σ. We show that the enhancement in the central M/L is best described either by the presence of black holes in these galaxies or by radial initial mass function variations. Assuming a black hole model, we derive black hole masses for the sample of galaxies. In two galaxies, NGC 4458 and NGC 4660, the data suggest significantly overmassive black holes, while in most others only upper limits are obtained. We also show that the level of M/L enhancements we see in these early-type galaxy nuclei are consistent with the larger enhancements seen in ultracompact dwarf galaxies (UCDs), supporting the scenario where massive UCDs are created by stripping galaxies of these masses.

Orbital tori for non-axisymmetric galaxies

Our Galaxy's bar makes the Galaxy's potential distinctly non-axisymmetric. All orbits are affected by non-axisymmetry, and significant numbers are qualitatively changed by being trapped at a resonance with the bar. Orbital tori are used to compute these effects. Thick-disc orbits are no less likely to be trapped by corotation or a Lindblad resonance than thin-disc orbits. Perturbation theory is used to create non-axisymmetric orbital tori from standard axisymmetric tori, and both trapped and untrapped orbits are recovered to surprising accuracy. Code is added to the TorusModeller library that makes it as easy to manipulate non-axisymmetric tori as axisymmetric ones. The augmented TorusModeller is used to compute the velocity structure of the solar neighbourhood for bars of different pattern speeds and a simple action-based distribution function. The technique developed here can be applied to any non-axisymmetric potential that is stationary in a rotating from - hence also to classical spiral structure.

Intracluster light at the Frontier – II. The Frontier Fields Clusters

Multiwavelength deep observations are a key tool to understand the origin of the diffuse light in clusters of galaxies: the intra-cluster light (ICL). For this reason, we take advantage of the Hubble Frontier Fields survey to investigate the properties of the stellar populations of the ICL of its 6 massive intermediate redshift (0.3<z<0.6) clusters. We carry on this analysis down to a radial distance of ~120 kpc from the brightest cluster galaxy. We found that the average metallicity of the ICL is [Fe/H] ~-0.5, compatible with the value of the outskirts of the Milky Way. The mean stellar ages of the ICL are between 2 to 6 Gyr younger than the most massive galaxies of the clusters. Those results suggest that the ICL of these massive (> 10^15 Msol) clusters is formed by the stripping of MW-like objects that have been accreted at z<1, in agreement with current simulations. We do not find any significant increase in the fraction of light of the ICL with cosmic time, although the redshift range explored is narrow to derive any strong conclusion. When exploring the slope of the stellar mass density profile, we found that the ICL of the HFF clusters follows the shape of their underlying dark matter haloes, in agreement with the idea that the ICL is the result of the stripping of galaxies at recent times.

Distribution and kinematics of atomic and molecular gas inside the solar circle

The detailed distribution and kinematics of the atomic and the CO-bright molecular hydrogen in the disc of the Milky Way inside the Solar circle are derived under the assumptions of axisymmetry and pure circular motions. We divide the Galactic disc into a series of rings, and assume that the gas in each ring is described by four parameters: its rotation velocity, velocity dispersion, midplane density and its scale height. We fit these parameters to the Galactic HI and CO (J=1-0) data by producing artificial HI and CO line-profiles and comparing them with the observations. Our approach allows us to fit all parameters to the data simultaneously without assuming a-priori a radial profile for one of the parameters. We present the distribution and kinematics of the HI and H2 in both the approaching (QIV) and the receding (QI) regions of the Galaxy. Our best-fit models reproduces remarkably well the observed HI and CO longitude-velocity diagrams up to a few degrees of distance from the midplane. With the exception of the innermost 2.5 kpc, QI and QIV show very similar kinematics. The rotation curves traced by the HI and H2 follow closely each other, flattening beyond R=6.5 kpc. Both the HI and the H2 surface densities show a) a deep depression at 0.5<R<2.5 kpc, analogous to that shown by some nearby barred galaxies, b) local overdensities that can be interpreted in terms of spiral arms or ring-like features in the disk. The HI (H2) properties are fairly constant in the region outside the depression, with typical velocity dispersion of 8.9+/-1.1 (4.4+/-1.2) km/s, density of 0.43+/-0.11 (0.42+/-0.22) cm-3 and HWHM scale height of 202+/-28 (64+/-12) pc. We also show that the HI opacity in the LAB data can be accounted for by using an `effective' spin temperature of about 150 K: assuming an optically thin regime leads to underestimate the HI mass by about 30%.

Most pseudo-bulges can be formed at later stages of major mergers

Most giant spiral galaxies have pseudo or disk-like bulges that are considered to be the result of purely secular processes. This may challenge the hierarchical scenario predicting about one major merger per massive galaxy ($>$$3\times 10^{10} M_{\odot}$) since the last $\sim$ 9 billion years. Here we verify whether or not the association between pseudo-bulges and secular processes is irrevocable. Using GADGET2 N-body/SPH simulations, we have conducted a systematic study of remnants of major mergers which progenitors have been selected (1) to follow the gas richness-look back time relationship, and (2) with a representative distribution of orbits and spins in a cosmological frame. Analyzing the surface-mass density profile of both nearby galaxies and merger remnants with two components, we find that most of them show pseudo-bulges or bar dominated centers. Even if some orbits lead to classical bulges just after the fusion, the contamination by the additional gas that gradually accumulates to the center and forming stars later on, leads to remnants apparently dominated by pseudo-bulges. We also found that simple SPH simulations should be sufficient to form realistic spiral galaxies as remnants of ancient gas-rich mergers without need for specifically tuned feedback conditions. We then conclude that pseudo-bulges and bars in spiral galaxies are natural consequences of major mergers when they are realized in a cosmological context, i.e., with gas-rich progenitors as expected when selected in the distant Universe.

On stellar mass loss from galaxies in groups and clusters

We estimate the stellar mass that satellite galaxies lose once they enter groups (and clusters) by identifying groups in a high-resolution cosmological N-body simulation, assigning entry masses to satellite galaxies with halo abundance matching at the entry time, and comparing the predicted conditional stellar mass function of satellite galaxies at $z\simeq0$ with observations. Our results depend on the mass of the stars that form in satellite galaxies after the entry time. A model in which star formation shuts down completely as soon a galaxy enters a group environment is ruled out because it underpredicts the stellar masses of satellite galaxies even in the absence of tidal stripping. The greater is the stellar mass that is allowed to form, the greater the fraction that needs to be tidally stripped. The stellar mass fraction lost by satellite galaxies after entering a group or cluster environment is consistent with any value in the range $0-40\%$. To place stronger constraints, we consider a more refined model of tidal stripping of galaxies on elongated orbits (where stripping occurs at orbit pericentres). Our model predicts less tidal stripping: satellite galaxies lose $\sim 20-25\%$ of their stellar mass since their entry into the group. This finding is consistent with a slow-starvation delayed-quenching picture, in which galaxies that enter a group or cluster environment keep forming stars until at least the first pericentric passage.

Dissecting galaxies: separating star formation, shock excitation and AGN activity in the central region of NGC 613

The most rapidly evolving regions of galaxies often display complex optical spectra with emission lines excited by massive stars, shocks and accretion onto supermassive black holes. Standard calibrations (such as for the star formation rate) cannot be applied to such mixed spectra. In this paper we isolate the contributions of star formation, shock excitation and active galactic nucleus (AGN) activity to the emission line luminosities of individual spatially resolved regions across the central 3 $\times$ 3 kpc$^2$ region of the active barred spiral galaxy NGC$\sim$613. The star formation rate and AGN luminosity calculated from the decomposed emission line maps are in close agreement with independent estimates from data at other wavelengths. The star formation component traces the B-band stellar continuum emission, and the AGN component forms an ionization cone which is aligned with the nuclear radio jet. The optical line emission associated with shock excitation is cospatial with strong $H_2$ and [Fe II] emission and with regions of high ionized gas velocity dispersion ($\sigma > 100$ km s$^{-1}$). The shock component also traces the outer boundary of the AGN ionization cone and may therefore be produced by outflowing material interacting with the surrounding interstellar medium. Our decomposition method makes it possible to determine the properties of star formation, shock excitation and AGN activity from optical spectra, without contamination from other ionization mechanisms.