Issue |
A&A
Volume 691, November 2024
|
|
---|---|---|
Article Number | A58 | |
Number of page(s) | 18 | |
Section | Extragalactic astronomy | |
DOI | https://doi.org/10.1051/0004-6361/202450500 | |
Published online | 29 October 2024 |
Multiple-scale gas infall through gravity torques on Milky Way twins
1
Observatoire de Paris, LERMA, PSL University, CNRS, Sorbonne University, Paris, France
2
Instituto de Astronomia, Geofísica e Ciências Atmosféricas, Departamento de Astronomia, Universidade de São Paulo, São Paulo, Brazil
3
Collège de France, 11 Pl. Marcelin Berthelot, 75231 Paris, France
⋆ Corresponding authors; patricia2.silva@alumni.usp.br; francoise.combes@obspm.fr
Received:
25
April 2024
Accepted:
4
June 2024
One of the main problems raised by the feeding of super-massive black holes (SMBHs) at the centres of galaxies is the huge angular momentum of the circumnuclear gas and of the gas reservoir in the galaxy disk. Because viscous torques are not efficient at kiloparsec or 100 pc scales, the angular momentum must be exchanged through gravity torques that arise from the non-axisymmetric patterns in the disks. Our goal here is to quantify the efficiency of bars and spirals in driving the gas towards the centre at different scales in galaxies. We selected a sample of nearby galaxies considered to be analogues of the Milky Way, that is, galaxies of late morphological type Sbc. Their bar strength was variable, either SB, or SAB, or SA, so that we were able to quantify the influence of the bar. The gravitational potential was computed from deprojected red images, either from Hubble Space Telescope or Legacy survey, depending on the spatial resolution and field of view considered. The torques were computed on the gas through CO emission maps from ALMA at different resolutions. Hα maps from MUSE were used, when available. Eight out of ten galaxies are barred. The torques are found to be negative in the eight barred objects at kiloparsec scales, between corotation and the inner Lindblad resonance (ILR), with a loss of angular momentum in a few rotations. Inside the ILR, the torques are negative in only five cases, with a timescale of one to two rotations. The torques are positive for the galaxies without bars. The torques applied on the ionized gas are comparable to what is deduced from molecular gas. The bars are confirmed to be the essential pattern in the SMBH feeding at kiloparsec and 100 pc scales; higher-resolution gas maps are required to explore scales of 10 pc.
Key words: galaxies: evolution / galaxies: kinematics and dynamics / galaxies: nuclei / galaxies: spiral
© The Authors 2024
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.
1. Introduction
Super-massive black holes (SMBHs) are ubiquitous at the centres of galaxies, and they play an important role in their evolution. A possible correlation exists between the mass of the black hole and the bulge mass, suggesting a co-evolution of the two (Kormendy & Ho 2013). One of the main phenomena invoked for this co-evolution is the feedback of active galaxy nucleus (AGN) episodes on the star formation rate in galaxies. AGN winds are able to produce molecular outflows from galactic nuclei, depleting the gaseous medium and moderating or quenching star formation in disks (Fiore et al. 2017). The efficiency of AGN winds is related to the Eddington ratio, and it is only significant when this ratio exceeds 0.01 (Heckman & Best 2014). Below this, the AGN feedback is more efficient in the radio mode through radio jets (García-Burillo et al. 2014; Dasyra et al. 2015).
Driving gas towards the galaxy centre to fuel both the central black hole and central star formation is an important ingredient in the co-evolution as well. The mechanisms for these inflows widely depend on the scale: at kiloparsec (kpc) scales in galaxy disks, viscous torques are negligible (Lin & Pringle 1987), and only gravity torques due to non-axisymmetric potentials are possible (Buta & Combes 1996). Stellar bars are the most efficient and long-lived features in this domain, and they are able to remove the angular momentum of the gas in a few rotations (García-Burillo et al. 2005). By computing gravity torques from the observed images of the old stellar component, together with the observed gas maps in nearby galaxies, García-Burillo & Combes (2012) concluded that the gas is driven inwards towards the SMBH in only 30% of cases, while it is either stalled in a nuclear ring at the inner Lindblad resonance (ILR) of the bar or driven outwards, from corotation to the outer Lindblad resonance (OLR). The gas infall requires multi-scale mechanisms to reach the very centre, as was already described by Shlosman et al. (1989, 1990), Berentzen et al. (2007). For example, recent high-resolution ALMA observations have revealed that embedded nuclear bars efficiently feed the nucleus at a scale of 10 pc (Combes et al. 2014; Audibert et al. 2019). Simulations of nuclear bars that feed the nucleus were made by Englmaier & Shlosman (2004), and Laine et al. (2002) observed that more Seyfert galaxies are present in barred galaxies, proving the importance of this feature in the feeding of the nucleus.
We focus on spiral galaxies of the same type as the Milky Way (MW), that is, with an Sbc morphology, with different bar strengths (SABbc and SBbc) and non-barred counterparts. The study of morphological twins of the MW will allow us to understand the history of its secular evolution. Although our MW is the best-known galaxy because Earth lies in it, it is necessary to study its nearby twins for a grasp on the evolution it may have had and to understand its formation through gas infall, fostered by the bar. Time is very important because the environment can change drastically. For instance, the bar strength is likely not to be a long-lived characteristic (1 to 2 Gyr if there is gas in the disk) because it varies according to gas flows (Bournaud & Combes 2002) or even according to many other mechanisms during cosmological evolution (Bi et al. 2022). Regarding the AGN activity, we know that the Eddington ratio of the MW nucleus is very low ∼10−9, but it could have been quite high in the past (Sarkar et al. 2023). When we take into account that AGN episodes are short-lived in galaxies, with a duty cycle of about 40 Myr (Shankar et al. 2009), the influence must be acknowledged as variable, if it exists.
We study the molecular gas torques in order to understand the feeding processes in the MW analogues, and also compared with the non-barred galaxies to quantify the effect of the bar strength. By using data from Atacama Large Millimeter/Submillimeter Array (ALMA), Hubble Space Telescope (HST) images in the F814W filter, and the Legacy survey images in the Z and R bands, we were able to map the torque of the gas and to understand better how the molecular gas is driven towards the nuclear region of these galaxies.
Section 2 presents the sample and the selection criteria, and Section 3 describes the method we used to compute gravity torques in the galaxies. We analyse the data in Sect. 4, and the results are discussed in Sect. 5. Sect. 6 gathers out conclusions.
2. Sample
In order to estimate the gravitational potential in the disk, we needed red images of the galaxies to avoid dust extinction and to avoid biasing by the light of a small number of OB stars with low mass-to-light ratios (M/L). The red images were available with several resolutions and fields of view (FOVs): We used HST images in the F814W filter and Legacy survey DR9 images in the R (for NGC 4303 and NGC 4536 because their Z-band images are saturated at the nucleus) and Z bands. To compute the torques on the molecular gas, we used ALMA data mainly from the Phangs survey (Leroy et al. 2021) and from the ALMA archive. The data details are displayed in Table 1.
Observational information of the sample galaxies.
We searched for the data of SABbc, SBbc, and SAbc galaxies with inclinations lower than 65° because the computations cannot be performed for edge-on galaxies. We also only selected nearby galaxies (D < 18 Mpc) to have enough spatial resolution, and galaxies with a magnitude in the B-band lower than 12 for a complete sample. Based on the classification in the NASA/IPAC Extragalactic Database (NED1), these criteria led to 43 galaxies. The availability of the data limited the sample to ten galaxies, with the exceptions of the two galaxies that were already studied (NGC 1566 Combes et al. 2014 and NGC 613 Audibert et al. 2019). Another two galaxies have ALMA data, but the FOV was too small for this work : NGC 1255 and NGC 4030, and there were problems with the data of a final two galaxies: NGC 3521 and NGC 5530. When these galaxies were removed, we obtained ten galaxies. Their properties are shown in Table 2. The inclinations of NGC 2903 and NGC 4536 are higher than 65°. However, given the uncertainties on this parameter, which is about 5° at least, we chose to include these galaxies in the sample.
Properties of the sample galaxies.
3. Method
The goal of this work is to quantify the rate of gas inflow due to gravity torques exerted by a bar (or spirals) on the molecular gas. This work was done at different scales, according to the spatial resolutions and FOVs of the observations.
The gravitational potential is mainly due to stars, and because the old population is more representative of the mass, red or near-infrared images are the best tracers, while the dark matter can be neglected in the central parts of MW-like galaxies (Persic et al. 1996). The red images at high resolution are usually from the filter F814W from HST for the central region and from the Legacy survey (at bands Z or R) at large scales to benefit from a wider FOV, but the spatial resolution is lower.
The images were rotated and deprojected according to the PAs and inclination angles given in Table 2. The bulge component was separated for the deprojection, because it was assumed to be spherical; in general, the bulge does not significantly contribute to the mass, because of the late morphological types of the targets. After rotation and deprojection, the images were Fourier transformed to compute the gravitational potential and forces. The M/L was assumed to be constant throughout the disk and bulge. Its precise value was considered not to be important because we computed normalised torques, corresponding to the relative loss of angular momentum per rotation. To take the third dimension of the disk into account, we assumed a stellar exponential disk thickness of ∼1/12th of the radial scale length of the galaxies (García-Burillo et al. 2005; Audibert et al. 2019). After computing the potential and forces in each pixel, we computed the torques with an emission map of the molecular gas at the same resolution. All maps were resampled to the same number of pixels and pixel sizes.
The amount of data available for both infrared images and molecular gas varies for a given galaxy. NGC 1300 has a considerable volume of data. In this case, we opted to use four ALMA data sets: two data sets were from the program 2015.1.00925.S and were combined in order to form a larger FOV (represented by the yellow square in Fig. 1 of NGC 1300). These data were used in the analysis with a Legacy image with a corresponding FOV. We used Phangs ALMA data associated with HST and Legacy images (green and red squares) and also the ALMA data from the program 2019.1.01718.S, which has a better spatial resolution, associated with the same HST image (cyan square). For NGC 1792, we used Phangs data together with HST (green square in Fig. 1 of NGC 1792) and Legacy images (red square). NGC 2903 also had four data sets: two data sets from the program 2019.1.01730.S, which had the same goal as NGC 1300. They were combined to form a larger FOV and were used together with a Legacy image that was matched to have the same size (yellow square). This galaxy also had data from Phangs that were used together with Legacy and HST data (red and green squares, respectively). We also used high spatial resolution ALMA data (program 2018.1.00517.S) together with HST data represented by the cyan square in Fig. 1 of NGC 2903. For NGC 3059, NGC 4303, NGC 4321, NGC 4536, and NGC 4689, only Phangs ALMA data were used, together with HST and Legacy data (with the exception of NGC 3059, which lacks Legacy images). For NGC 5248, the Phangs data were used together with HST and Legacy images as well; however, ALMA data were also used from the program 2019.1.00876.S (with higher spatial resolution) together with the HST images. Finally, for NGC 5921, we used ALMA data from 2019.1.01742.S program with a high spatial resolution, together with the HST image of this galaxy. This galaxy has Legacy images; because the ALMA data have a small FOV, it would provide the same results as were obtained with HST images with a lower spatial resolution, however. For better results, we used HST data. The summary of all information is shown in Table 1 and Fig. 1.
![]() |
Fig. 1. Composite of Legacy images of the G (blue), R (green), and Z (red) bands (except for NGC 3059, which lacks Legacy images; in this case, we used the images from the Digitalized Sky Survey). The squares represent the size of the FOV that was used for each data set: green for HST and ALMA, cyan for the HST and the ALMA high spatial resolution set, red for Legacy and ALMA, and yellow for the Legacy and ALMA larger FOV (combination of two FOVs). The pink lines represent the PA of the bar. The orientation N–E is north up and east left. |
To prepare the data for the computation of gravitational torques from the infrared images, we first removed the MW bright stars in the FOV. We then rotated the images using the PA of the galaxy (Table 2) so that its major axis was horizontal. This implied the definition of a common centre. Using the infrared image, we established a centre based on the emission peak of the nuclear source. We also centred the ALMA data cube according either to the emission peak or to the centre of the ring, and we used the velocity field. The two components (infrared images and ALMA data cube) were then re-centred so that they had the same centre. The images were then deprojected by stretching the vertical axis. For the ALMA data, the same rotation and deprojection were applied to the data cube converted into velocities relative to the centre. These processes were performed with scripts in Python and Gildas environments.
4. Data analysis
4.1. Gas morphology and kinematics
Figs. 2a and b show the nuclear structure of NGC 1300. The morphology is that of a ring with a radius of about 300 pc (already detected in previous works; see Appendix A), which is decomposed into two concentric rings in the ALMA mean image. There is clear evidence of a structure that connects the nuclear source to the ring. The nuclear low-resolution image reveals that this ring is connected to the spiral structure. This is better visible in the images with a larger FOV (Figs. 2c and d). The ring and spiral obey the same regular rotation. The two spiral arms emerge outside the ring. They first correspond to the leading dust lanes in the bar, and then continue to wind out outside the bar.
![]() |
Fig. 2. NGC 1300. Left panel: Images of the F814W filter of HST (a, b), and Legacy Z band (c, d). Middle and right panels: ALMA data (mean flux and velocity map) associated with each FOV, as described in section 3. The orientation N–E is north up and east right. The cross (arbitrary size) indicates the centre of the image, which was taken from the emission peak from the infrared image and was used to match the ALMA images. |
Figs. 3a and b show NGC 1792 data in two different FOVs. They show the nuclear region, which is considerably flocculent. The nuclear source appears to be divided into two regions (as seen in the ALMA data). Because of the spatial resolution, however, we cannot conclude this with certainty. This galaxy has a short nuclear bar that might be related to this structure (see Appendix A). The gas clearly rotates in both FOVs. The spiral arms from the galaxy appear to be connected to the nuclear source seen in the nuclear ALMA image. This galaxy is classified as non-barred, but these maps are quite similar to those of the barred counterparts.
![]() |
Fig. 3. NGC 1792. Same as the left panel of Fig. 2. (a) Image of the HST filter F814W and (b) image of the Z band from the Legacy survey. |
The structure of the nuclear region of NGC 2903 (Fig. 4) is very flocculent and hard to interpret. However, a nuclear ring can be distinguished in the gas map. According to the literature, this ring has radius of about 600 pc (see Appendix A for more information). This structure is also connected to the spiral arms of the galaxy and follows the orientation of the galaxy bar (clearly visible in panel d). The gas obeys a clear rotation, and, as in NGC 1300, the ring also follows the same rotation pattern as the remaining galaxy disk.
The data from NGC 3059 (Fig. 5) do not present much detail because of the spatial resolution; the bar is very thin and strong in this almost face-on galaxy. The nuclear source is connected to the spiral arms in the same direction as the bar. The rotation is also clear in this case. This galaxy was reported to have a weak bar (see Appendix A), but we found that this galaxy has the strongest bar in the sample.
NGC 4303 (Fig. 6) presents a nuclear ring at the bar ILR in the centre that is connected to the spiral arms, leading dust lanes along the bar. These in turn are connected to the spiral arms outside the bar. The ring radius is about 310 pc. The spiral arms are very clear and rotate regularly. The spatial resolution did not allow us to see the nuclear bar that was reported in the literature (see Appendix A for more details).
The gas morphology of NGC 4321 (Fig. 7) is very similar to that of NGC 4303. The spiral arms follow the morphology typical of a barred galaxy. The nuclear ring radius is ∼500 pc. As in 4303, the spatial resolution was not high enough for us to detect the nuclear bar that was observed in previous studies (see Appendix A).
NGC 4536 has a very obscured nucleus, with a conspicuous dust lane on the near side, in the north-east. The bar orientation of this galaxy is confused with its major axis because of the high inclination. However, in Fig. 8, molecular gas images show a clear nuclear bar with a PA ∼ 136° and a diameter of about 2.1 kpc that was also detected in previous studies (see Appendix A). One gas spiral arm is also connected to the bar in the south-west. The opposite arm in the north-east is better visible in the stellar component.
NGC 4689 is another non-barred galaxy. In this case, the difference between this non-barred galaxy and the barred galaxies is clear: The nucleus, which is very flocculent, lacks any clear structure (Fig. 9). Some studies reported evidence of a weak bar, but we detected no clear structure in this galaxy (see Appendix A). The radio data show no clear nuclear source either. Moreover, the spiral structure in the molecular gas images is not as clear as in the infrared images.
NGC 5248 has a very clear spiral structure and a nuclear ring at the centre, with a radius of about 88 pc (Fig. 10). It also has an inner pseudo-ring that encircles the bar. This is expected to delineate the ultra-harmonic resonance near the corotation at a radius ∼500 pc. Then the spiral arms connect to this pseudo-ring outside the bar. We did not detect the possible nuclear bar (see Appendix A). This might be due to the low spatial resolution.
![]() |
Fig. 10. NGC 5248. Same as Fig. 2. Panels (a) and (b) show images of the HST filter F814W, and panel (c) shows an image from the Z band from the Legacy survey. |
Several structures are embedded in galaxy NGC 5921 (Fig. 11). The morphology at large scales suggests that an inner ring encircles the bar with a radius ∼5 kpc. Two spiral arms are connected to this ring/bar (see Fig. 1). At smaller scales, a nuclear ring with a radius of about 230 pc is visible. Nuclear spiral arms are also connected to this nuclear source.
![]() |
Fig. 11. NGC 5921. Left panel: F814W image from HST. Middle and right panels: ALMA data as in Fig. 2. |
In summary, all barred galaxies of the sample present a very complex nuclear region. Six galaxies present clear evidence of nuclear rings (N1300, N2903, N4303, N4321, N5248, and N5921). Two objects have a nuclear bar (N4536 and N5921). NGC 4689 alone did not show clear structures in its nucleus. This might be related to the lack of a bar. Although it was classified as a non-barred galaxy, NGC 1792 has a clear structure in its nucleus.
4.2. Gas torque analysis
The computation of the gravity potential from the deprojected red images yields the estimate of the bar and spiral strength as a function of radius. The nature of the pattern (bar or spiral) is indicated by the phase and its variation with radius of the different Fourier components into which the potential Φ(R, θ) is decomposed,
The strength of the m-Fourier component, Qm(R), is defined as Qm(R) = mΦm/R|F0(R)|, that is, by the ratio between tangential and radial forces (Combes & Sanders 1981). The strength of the total non-axisymmetric perturbation QT(R) is defined similarly, with the maximum amplitude of the tangential force FTmax(R). Their radial distributions and the radial phase variations are displayed for a typical galaxy in Fig. 12. We determine in Table 3 the radii and maximum amplitudes of QT(R) over the disk for all galaxies.
![]() |
Fig. 12. Radial distribution of the amplitudes (top) and phases (bottom) of the non-axisymmetric perturbations in each Fourier component (m = 1 to 4). QT represents the total perturbation curve in the potential for the data set of N4303 (b). This plot illustrates the typical potential component curves obtained in the sample. The maximum values of the Q2 and QT curves are shown in Table 3. |
Maximum and minimum values with their associated radii of the torque and relative maxima of non-axisymmetric potential curves (Q2 and QT).
For each galaxy, the forces per unit mass (Fx and Fy) at each pixel were obtained by deriving the potential. The torques per unit mass t(x, y) were then computed by
The sign of the torque was determined relative to the sense of rotation in the plane of the galaxy. The net effect on the gas at each radius is then the product of the torque by the gas density Σ. This quantity, t(x, y)×Σ(x, y), is shown in the left panels of Figs. 13–23.
![]() |
Fig. 13. NGC 1300 torque maps from the same set of data as in panel (a) Fig. 2a, (b) Fig. 2b, (c) Fig. 2c, and (d) Fig. 2d. The images were rotated and de-projected. The orientation N–E is indicated in the torque maps. The crosses represent the centre of the galaxy estimated from the emission peaks of the infrared images and ALMA data. The dark violet line represents the PA of the bar after rotation. Each map is associated with a graph that represents the radial distribution of the torque, averaged over the azimuth, and normalised to the angular momentum, per rotation at this radius. The colour scales in the maps are in arbitrary units (except for panels c and d), of the torques t(x,y) weighted by the gas surface density Σ(x,y). The contours of Σ(x,y) are from the mean molecular gas map (a) one to nine times 18.55 Jy/beam km/s, (b) 1, 1.1, 1.4, 2.5, 5.9, 16.5 times 0.3 Jy/beam km/s (c) 1, 5, 11, 23.8, 51.5 times 1.65 Jy/beam km/s, and (d) 1, 2.5, 12.5, 27.2, 58.9, 127.2 times 0.49 Jy/beam km/s. The velocities are relative to the centre. |
![]() |
Fig. 14. NGC 1792 torque maps from the same set of data as in panel (a) Fig. 3a and (b) Fig. 3b. Same caption as Fig. 13. The contour levels were taken from the mean molecular gas map (a) one to six times 34.13 Jy/beam km/s and (b) one to five times 32.11 Jy/beam km/s. |
![]() |
Fig. 15. NGC 2903 torque maps from the same set of data as in panel (a) Fig. 4a, (b) Fig. 4b, (c) Fig. 4c, and (d) Fig. 4d. Same caption as Fig.13. The contour levels were taken from the mean molecular gas map (a) 1, 4, 9, 16, 25, 36, and 49, times 6.8 Jy/beam km/s, (b) one to seven times 1.08 Jy/beam km/s, (c) one to five times 1.43 Jy/beam km/s, and (d) 1, 4, 9, 16, 25, and 36 times 6.3 Jy/beam km/s. |
![]() |
Fig. 16. NGC 3059 torque maps from the same set of data as in Fig. 5. Same caption as Fig. 13. The contour levels were taken from the mean molecular gas map 1, 2, 3, 4, 6, 8, 10, and 11 times 15.7 Jy/beam km/s. |
![]() |
Fig. 17. NGC 4303 torque maps from the same set of data as in panel (a) Fig. 6a and (b) Fig. 6b. Same caption as Fig. 13. The contour levels were taken from the mean molecular gas map (a) 1, 2, 4, 5, 7, 9, 10, and 12 times 14.6 Jy/beam km/s, and (b) 1, 2, 5, 6, 8, 9, and 10 times 15.9 Jy/beam km/s. |
![]() |
Fig. 18. NGC 4321 torque maps from the same set of data as in panel (a) Fig. 7a and (b) Fig. 7b. Same caption as Fig. 13. The contour levels were taken from the mean molecular gas map (a and b) 1, 2, 5, 10, 23, 34, and46 times 4.4 Jy/beam km/s. |
![]() |
Fig. 19. NGC 4536 torque maps from the same set of data as in panel (a) Fig. 8a and (b) Fig. 8b. Same caption as Fig. 13. The contour levels were taken from the mean molecular gas map (a) 1, 4.2, 14.2, 45.8, and 145.8 times 1.56 Jy/beam km/s, and (b) 1, 1.1, 1.4, 2, 3.1, 5.5, 10, 19.3, 37.6, and 74.1 times 5 Jy/beam km/s. |
![]() |
Fig. 20. NGC 4689 torque maps from the same set of data as in panel (a) Fig. 9a and (b) Fig. 9b. Same caption as Fig. 13. The contour levels were taken from the mean molecular gas map (a) 1, 1.4, 1.7, 2, 2.2, and 2.5 times 6.8 Jy/beam km/s, and (b) one to three times 4.9 Jy/beam km/s. |
![]() |
Fig. 21. NGC 5248 torque maps from the same set of data as in panel (a) Fig. 10a, (b) Fig. 10b, and (c) Fig. 10c. Same caption as Fig. 13. The contour levels were taken from the mean molecular gas map (a) 1, 1.4, 1.7, 2, 2.2, 2.5, and 2.7 times 59.3 Jy/beam km/s, and (b) one to five times 0.25 Jy/beam km/s. |
![]() |
Fig. 22. cont. of Fig. 21. The contour levels were taken from the mean molecular gas map (c) 1, 4, 9, 16, 25, and 36 times 3.48 Jy/beam km/s. |
![]() |
Fig. 23. NGC 5921 torque maps from the same set of data as in Fig. 11. Same caption as Fig. 13. The contour levels were taken from the mean molecular gas map one to seven times 0.74 Jy/beam km/s. |
The torque weighted by the gas density Σ(x, y) was then averaged over azimuth,
The quantity t(R) represents the time derivative of the specific angular momentum L = ΩR2 of the gas averaged azimuthally (García-Burillo et al. 2005). The efficiency of the gas flow was then estimated through normalising t(R) at each radius by the angular momentum L, divided by the rotation period T = 2π/Ω, as shown in the right panels of Figs. 13–23. With this normalisation, the computations are independent of the constant M/L ratios used.
These figures show some general features. Negative torques, implying gas infall between corotation and ILR, as expected, are indeed observed in several barred galaxies when the FOV is sufficient and includes enough of the bar region. This is the case for NGC 1300 in Figs. 13c and d, and also NGC 4303, NGC 4321, NGC 4536, NGC 5248, and NGC 5921. When only the nuclear ring is included in the FOV, the situation is more complex, and the torques do not show a single coherent behaviour, as for NGC 1300 in Figs. 13a and b. The torques are not significant in non-barred flocculent galaxies, as for NGC 1792 in Fig. 14 and NGC 4689 in Fig. 20.
NGC 2903 is a special case (Fig. 15). Although barred, with a coherent spiral structure, the torque values are more chaotic. This is due to several factors. The main factor is the patchiness of the CO emission, which prevents a spatial continuity in the torque computation. We note that the torques we computed from the two different fields of view are sometimes different. The reason is the different gas maps that were obtained with different configurations of the ALMA interferometer, with different spatial resolutions. Another reason is that the red image of NGC 2903 is quite dusty, with conspicuous dust lanes. This leads to an irregular potential estimation that differs from the real potential.
NGC 3059 has a very thin bar, almost on the minor axis. It is therefore even longer when deprojected, although the inclination is quite low. The torque is positive from 0.3 to 0.6 kpc, and then negative from 0.6 to 1 kpc (Fig. 16). This would mean that gas is driven to a ring at a radius of 0.6 kpc. However, in this case, the gas orbits are expected to be very elongated, and positive and negative torques might partially cancel out over the orbits. The torques are negative at radii lower than 0.3 kpc, which explains the gas concentration inside a radius of 3 arcsec, that is, 220 pc.
The quantification of the torques for each galaxy and each FOV is displayed in Table 3. The relative loss and gain of angular momentum per rotation is higher than one in absolute value only for NGC 1300, and it is not indicated in this table. It does not satisfy the main hypothesis underlying the azimuthal average of the torque, that the gas has time to remain in its radial range during one rotation. In this case, we plot only the torque map, normalized with physical quantities, computed with the rotation curve of the galaxy. The latter was derived from the potential itself, assuming a constant mass-to-light ratio, and normalised to the CO, HI, and optical observations (see Fig. 24). In this figure, the rotation curve from the potential model is compared to the observed curves from Lang et al. (2020), England (1989), and Lindblad et al. (1997).
![]() |
Fig. 24. Gas velocity rotation curves from HI (Lindblad et al. 1997) (red squares) and from CO (Lang et al. 2020) (red points) of NGC 1300 used to compute the torque map values of the data sets in panels (c) and (d). Only the highest points of each curve were kept. The blue curve represents the velocity model obtained from the computed potential (with limitations of the red image cut-off beyond 5 kpc). |
5. Discussion
We studied ten MW analogues, and only two galaxies of this sample are considered non-barred. This sample is homogeneous when we consider that galactic bars may be in a 1-2 Gyr transient phase (e.g. Bournaud & Combes 2002) and when, as previous studies argued (see Appendix A), the AGN activity is weak in these nuclei (which is also an intermittent characteristic, with a duty cycle of 40 Myr). The comparison of gas flow processes within this matched sample can help us to further understand the evolution of these types of galaxies, including our own. Although this sample is small because available data are limited, it is representative because it was selected based on objective criteria. The results are consistent for this sample.
It is important to discuss the limitations of this work. One limitation is the spatial resolution of the data. One of the most frequently encountered questions in astronomy concerns the feeding and feedback processes in galactic nuclei. In order to study this matter comprehensively, high spatial resolution data are fundamental, because these processes occur in a very small region (10–100 pc). The spatial resolution of the data we analysed here is insufficient to compute the gas dynamics at these scales in detail. As a consequence, we disregarded results below 1 arcsec (corresponding to 66 pc at an average distance of 13.7 Mpc) in general. The torque maps and values in this central region are uncertain, not only because of the spatial resolution, but sometimes also due to phenomena associated with possible AGN, which may generate ambiguous results that can only be differentiated with data with a higher spatial resolution. These are rare. Not many galaxies are observed with the necessary instrumental configuration for an analysis like this. The spatial resolution we used is considered moderate, and can only provide a general idea of what occurs in these galaxies at scales of 100 pc – 1 kpc.
The assumption that dark matter has no effect is a good approach because these galaxies are massive enough (see the stellar masses in Table 2). In massive disk-dominated galaxies like this, dark matter is not dominant in the central region and shows a cored distribution (de Blok et al. 2008). In addition, we studied features in the nuclear regions of the objects, where baryons dominate the gravitational potential.
The torque maps were computed in each pixel, weighted by the gas density. To derive the global impact on the gas inflow/outflow, we averaged on the azimuth at each radius, although the orbits of the gas, when deprojected, are not circular. We then quantified the relative angular momentum variation in one circular rotation. Although gas orbits might be elliptical in these galaxies, this average gives the correct order of magnitude over many gas streamlines at a given radius. We did not considered this average when the relative loss of angular momentum was higher than one in absolute value (e.g. NGC 1300).
In the case of NGC 1300, we have two data sets (c and d), which are related to the large FOVs that provided noisier results, but also large variations in angular momentum. As a consequence, we could not take the self-normalised torque curves into account and had to use the velocity rotation curve (Fig. 24) in order to normalise the results. Some parts of the CO maps of these data sets show no emission, and the S/N in other parts is too low, which causes the torque curves to be noisy. The same was observed in other galaxies, but with fewer fluctuations.
Another limitation of this work is the determination of the PA and the galaxy inclination. This relies on the assumption that the outer optical isophotes are circular. These determinations can be biased in barred galaxies, where the PA of the bar can easily be misinterpreted to be aligned with the galaxy major axis, causing systematic errors. The inclination angle has the same problem because these galaxies are not axisymmetric disks. The results strongly depend on these values because the deprojection is important for the computation of the torques.
We did not detect any dynamical influence of AGN feedback in these galaxies. One of the reasons, as mentioned, is the low spatial resolution and low AGN activity. We did not detect radio jets either. An optical data analysis (with a higher spatial resolution within the 500 pc scale) would be interesting in this case, in order to compare the ionized and molecular gas behaviour, and also to determine whether the AGN affects the surrounding ionized gas.
We provided the characteristics of gas flows in MW twins at 100 pc–1 kpc scales. We consider this study as preliminary results that are to be pursued in the future with data that have a better spatial resolution, and perhaps with multiwavelength data. With these data, it will be possible to connect what we observe at these scales with the smaller scales and to quantify the AGN fueling mechanisms and/or the feedback on the gas flow in these galaxies. In ringed galaxies, the infall of gas is expected to be contained within the ring, but the gas can also be driven towards the nucleus through nuclear spirals/bars.
6. Conclusions
We have quantified the gravity torques and their impact on the gas flows in a sample of eight MW analogues with two non-barred counterparts. The sample is homogeneous in mass, morphological type, and fundamental physical characteristics. The bar and spiral strengths varied, as did the AGN types, because they are considered variable during the lifetimes of typical spiral galaxies. This allowed us to compare the effect of non-axisymmetric features on the gas inflows and outflows. By using HST and Legacy images of the stellar emission and ALMA data of the molecular gas emission, we were able to compute the molecular gas torques in this sample at scales of hundreds of parsecs to kiloparsecs.
The galaxy sample presented a diversification of structures, such as rings (NGC 1300, NGC 2903, NGC 4303, NCG 4321, NGC 5248, and NGC 5921) and spiral arms connected to them. The sample also contains flocculent galaxies (the best example is NGC 4689, without any recognisable structure).
Our multi-scale study showed that torques are not significant in the inner resonant rings, even when they are not circular. They are not significant for flocculent galaxies either. For barred galaxies, torques are negative between the corotation and the inner Lindblad resonance, as expected. The torques computed on the Hα emitting gas are compatible with those on the molecular gas, but they are noisier and more patchy.
The relative loss of angular momentum can be quite high and may reach nearly 100% in one rotation for extreme cases. When the bar is very strong and thin, the gas orbits are very elongated, which prevents the formation of a ring. The gas is then effectively driven to the nucleus.
We studied intermediate scales between 100 pc and 1 kpc scales. The gas inflow revealed by the torques produces starbursts in the centre of barred galaxies. Future studies at higher spatial resolution are needed to follow-up the AGN feeding into the nucleus.
Acknowledgments
PS acknowledges Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) post-doctoral fellowships 2020/13239-5 and 2022/14382-1. The Legacy Surveys consist of three individual and complementary projects: the Dark Energy Camera Legacy Survey (DECaLS; Proposal ID #2014B-0404; PIs: David Schlegel and Arjun Dey), the Beijing-Arizona Sky Survey (BASS; NOAO Prop. ID #2015A-0801; PIs: Zhou Xu and Xiaohui Fan), and the Mayall z-band Legacy Survey (MzLS; Prop. ID # 2016A-0453; PI: Arjun Dey). DECaLS, BASS and MzLS together include data obtained, respectively, at the Blanco telescope, Cerro Tololo Inter-American Observatory, NSF’s NOIRLab; the Bok telescope, Steward Observatory, University of Arizona; and the Mayall telescope, Kitt Peak National Observatory, NOIRLab. Pipeline processing and analyses of the data were supported by NOIRLab and the Lawrence Berkeley National Laboratory (LBNL). The Legacy Surveys project is honored to be permitted to conduct astronomical research on Iolkam Du’ag (Kitt Peak), a mountain with particular significance to the Tohono O’odham Nation. The Digitized Sky Surveys were produced at the Space Telescope Science Institute under U.S. Government grant NAG W-2166. The images of these surveys are based on photographic data obtained using the Oschin Schmidt Telescope on Palomar Mountain and the UK Schmidt Telescope. The plates were processed into the present compressed digital form with the permission of these institutions. This research is based on observations made with the NASA/ESA Hubble Space Telescope obtained from the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555. These observations are associated with program(s) 15133, 15654, 9788, 9042, 15323. This paper makes use of the following ALMA data: ADS/JAO.ALMA#2015.1.00925.S, ADS/JAO.ALMA#2015.1.00956.S, ADS/JAO.ALMA#2017.1.00886.L, ADS/JAO.ALMA#2018.1.00517.S, ADS/JAO.ALMA#2018.1.01651.S, ADS/JAO.ALMA#2019.1.01718.S, ADS/JAO.ALMA#2019.1.01730.S, ADS/JAO.ALMA#2019.1.00876.S, and ADS/JAO.ALMA#2019.1.01742.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), NSTC and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. This work was supported by the Programme National Cosmology et Galaxies (PNCG) of CNRS/INSU with INP and IN2P3, co-funded by CEA and CNES.
References
- Alonso-Herrero, A., Ryder, S. D., & Knapen, J. H. 2001, MNRAS, 322, 757 [NASA ADS] [CrossRef] [Google Scholar]
- Argyle, R. W., & Eldridge, P. 1990, MNRAS, 243, 504 [NASA ADS] [Google Scholar]
- Arsenault, R. 1989, A&A, 217, 66 [NASA ADS] [Google Scholar]
- Arsenault, R., Boulesteix, J., Georgelin, Y., & Roy, J. R. 1988, A&A, 200, 29 [NASA ADS] [Google Scholar]
- Atkinson, J. W., Collett, J. L., Marconi, A., et al. 2005, MNRAS, 359, 504 [NASA ADS] [CrossRef] [Google Scholar]
- Audibert, A., Combes, F., García-Burillo, S., et al. 2019, A&A, 632, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Beck, R., Shoutenkov, V., Ehle, M., et al. 2002, A&A, 391, 83 [CrossRef] [EDP Sciences] [Google Scholar]
- Berentzen, I., Shlosman, I., Martinez-Valpuesta, I., & Heller, C. H. 2007, ApJ, 666, 189 [Google Scholar]
- Berrier, J. C., Davis, B. L., Kennefick, D., et al. 2013, ApJ, 769, 132 [NASA ADS] [CrossRef] [Google Scholar]
- Bi, D., Shlosman, I., & Romano-Díaz, E. 2022, ApJ, 934, 52 [NASA ADS] [CrossRef] [Google Scholar]
- Bonatto, C., Bica, E., & Alloin, D. 1989, A&A, 226, 23 [NASA ADS] [Google Scholar]
- Bottema, R. 1990, A&A, 233, 372 [NASA ADS] [Google Scholar]
- Bournaud, F., & Combes, F. 2002, A&A, 392, 83 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Burtscher, L., Davies, R. I., Shimizu, T. T., et al. 2021, A&A, 654, A132 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Buta, R., & Combes, F. 1996, Fund. Cosmic Phys., 17, 95 [Google Scholar]
- Buta, R., & Crocker, D. A. 1993, AJ, 105, 1344 [NASA ADS] [CrossRef] [Google Scholar]
- Canzian, B. 1993, ApJ, 414, 487 [NASA ADS] [CrossRef] [Google Scholar]
- Castillo-Morales, A., Jiménez-Vicente, J., Mediavilla, E., & Battaner, E. 2007, MNRAS, 380, 489 [NASA ADS] [CrossRef] [Google Scholar]
- Casuso, E., Beckman, J. E., & Cepa, J. 1990, Ap&SS, 170, 275 [NASA ADS] [CrossRef] [Google Scholar]
- Chiang, I.-D., Sandstrom, K. M., Chastenet, J., et al. 2024, ApJ, 964, 18 [NASA ADS] [CrossRef] [Google Scholar]
- Colina, L., & Arribas, S. 1999, ApJ, 514, 637 [NASA ADS] [CrossRef] [Google Scholar]
- Colina, L., & Wada, K. 2000, ApJ, 529, 845 [NASA ADS] [CrossRef] [Google Scholar]
- Colina, L., García Vargas, M. L., Mas-Hesse, J. M., Alberdi, A., & Krabbe, A. 1997, ApJ, 484, L41 [NASA ADS] [CrossRef] [Google Scholar]
- Colina, L., Gonzalez Delgado, R., Mas-Hesse, J. M., & Leitherer, C. 2002, ApJ, 579, 545 [NASA ADS] [CrossRef] [Google Scholar]
- Combes, F., & Sanders, R. H. 1981, A&A, 96, 164 [NASA ADS] [Google Scholar]
- Combes, F., García-Burillo, S., Casasola, V., et al. 2014, A&A, 565, A97 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Comerón, S., Knapen, J. H., Beckman, J. E., et al. 2010, MNRAS, 402, 2462 [Google Scholar]
- Comerón, S., Salo, H., Laurikainen, E., et al. 2014, A&A, 562, A121 [Google Scholar]
- Constantin, A., Shields, J. C., Ho, L. C., et al. 2015, ApJ, 814, 149 [NASA ADS] [CrossRef] [Google Scholar]
- Dahlem, M. 1992, A&A, 264, 483 [NASA ADS] [Google Scholar]
- Dahlem, M., Bomans, D. J., & Will, J.-M. 1994, ApJ, 432, 590 [NASA ADS] [CrossRef] [Google Scholar]
- Das, M., Teuben, P. J., Vogel, S. N., et al. 2003, ApJ, 582, 190 [NASA ADS] [CrossRef] [Google Scholar]
- Dasyra, K. M., Bostrom, A. C., Combes, F., & Vlahakis, N. 2015, ApJ, 815, 34 [Google Scholar]
- Davies, R. I., Sugai, H., & Ward, M. J. 1997, MNRAS, 291, 314 [NASA ADS] [CrossRef] [Google Scholar]
- Davis, B. L., Graham, A. W., & Combes, F. 2019, ApJ, 877, 64 [NASA ADS] [CrossRef] [Google Scholar]
- de Blok, W. J. G., Walter, F., Brinks, E., et al. 2008, AJ, 136, 2648 [NASA ADS] [CrossRef] [Google Scholar]
- de Vaucouleurs, G., de Vaucouleurs, A., Corwin, Jr., H. G., et al. 1991, Third Reference Catalogue of Bright Galaxies [Google Scholar]
- Díaz-García, S., Lisenfeld, U., Pérez, I., et al. 2021, A&A, 654, A135 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Dutil, Y., & Roy, J.-R. 1999, ApJ, 516, 62 [NASA ADS] [CrossRef] [Google Scholar]
- Elmegreen, B. G., Elmegreen, D. M., & Montenegro, L. 1992, ApJS, 79, 37 [NASA ADS] [CrossRef] [Google Scholar]
- Elmegreen, D. M., Elmegreen, B. G., Chromey, F. R., & Hasselbacher, D. A. 1996, ApJ, 469, 131 [NASA ADS] [CrossRef] [Google Scholar]
- Elmegreen, D. M., Chromey, F. R., Santos, M., & Marshall, D. 1997, AJ, 114, 1850 [NASA ADS] [CrossRef] [Google Scholar]
- England, M. N. 1989, ApJ, 344, 669 [NASA ADS] [CrossRef] [Google Scholar]
- Englmaier, P., & Shlosman, I. 2004, ApJ, 617, L115 [NASA ADS] [CrossRef] [Google Scholar]
- Erwin, P. 2004, A&A, 415, 941 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Eskridge, P. B., Frogel, J. A., Pogge, R. W., et al. 2002, ApJS, 143, 73 [NASA ADS] [CrossRef] [Google Scholar]
- Fathi, K., Beckman, J. E., Piñol-Ferrer, N., et al. 2009, ApJ, 704, 1657 [NASA ADS] [CrossRef] [Google Scholar]
- Fiore, F., Feruglio, C., Shankar, F., et al. 2017, A&A, 601, A143 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Font, J., Beckman, J. E., Zaragoza-Cardiel, J., et al. 2014, MNRAS, 444, L85 [NASA ADS] [CrossRef] [Google Scholar]
- Fujimoto, Y., Maeda, F., Habe, A., & Ohta, K. 2020, MNRAS, 494, 2131 [NASA ADS] [CrossRef] [Google Scholar]
- Gadotti, D. A., Sánchez-Blázquez, P., Falcón-Barroso, J., et al. 2019, MNRAS, 482, 506 [Google Scholar]
- Gadotti, D. A., Bittner, A., Falcón-Barroso, J., et al. 2020, A&A, 643, A14 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- García-Burillo, S., & Combes, F. 2012, J. Phys. Conf. Ser., 372, 012050 [CrossRef] [Google Scholar]
- Garcia-Burillo, S., Sempere, M. J., Combes, F., & Neri, R. 1998, A&A, 333, 864 [Google Scholar]
- García-Burillo, S., Combes, F., Schinnerer, E., Boone, F., & Hunt, L. K. 2005, A&A, 441, 1011 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- García-Burillo, S., Combes, F., Usero, A., et al. 2014, A&A, 567, A125 [Google Scholar]
- Garrison, R. F., & Walborn, N. R. 1974, JRASC, 68, 117 [NASA ADS] [Google Scholar]
- González-Martín, O., Masegosa, J., Márquez, I., Guainazzi, M., & Jiménez-Bailón, E. 2009, A&A, 506, 1107 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Goulding, A. D., & Alexander, D. M. 2009, MNRAS, 398, 1165 [NASA ADS] [CrossRef] [Google Scholar]
- Goulding, A. D., Alexander, D. M., Lehmer, B. D., & Mullaney, J. R. 2010, MNRAS, 406, 597 [NASA ADS] [CrossRef] [Google Scholar]
- Graham, A. W., Soria, R., & Davis, B. L. 2019, MNRAS, 484, 814 [NASA ADS] [Google Scholar]
- Haan, S., Schinnerer, E., Emsellem, E., et al. 2009, ApJ, 692, 1623 [NASA ADS] [CrossRef] [Google Scholar]
- Heckman, T. M., & Best, P. N. 2014, ARA&A, 52, 589 [Google Scholar]
- Hernandez, O., Carignan, C., Amram, P., Chemin, L., & Daigle, O. 2005, MNRAS, 360, 1201 [CrossRef] [Google Scholar]
- Ho, L. C., Filippenko, A. V., & Sargent, W. L. W. 1997, ApJS, 112, 315 [NASA ADS] [CrossRef] [Google Scholar]
- Hoyer, N., Arcodia, R., Bonoli, S., et al. 2024, A&A, 682, A36 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Hughes, M. A., Axon, D., Atkinson, J., et al. 2005, AJ, 130, 73 [NASA ADS] [CrossRef] [Google Scholar]
- Jackson, J. M., Snell, R. L., Ho, P. T. P., & Barrett, A. H. 1989, ApJ, 337, 680 [NASA ADS] [CrossRef] [Google Scholar]
- Jackson, J. M., Eckart, A., Cameron, M., et al. 1991, ApJ, 375, 105 [NASA ADS] [CrossRef] [Google Scholar]
- Jang, I., Gliozzi, M., Hughes, C., & Titarchuk, L. 2014, MNRAS, 443, 72 [CrossRef] [Google Scholar]
- Jiménez-Bailón, E., Santos-Lleó, M., Mas-Hesse, J. M., et al. 2003, ApJ, 593, 127 [CrossRef] [Google Scholar]
- Jiménez-Vicente, J., Castillo-Morales, A., Mediavilla, E., & Battaner, E. 2007, MNRAS, 382, L16 [CrossRef] [Google Scholar]
- Jogee, S., Shlosman, I., Laine, S., et al. 2002, ApJ, 575, 156 [NASA ADS] [CrossRef] [Google Scholar]
- Jogee, S., Scoville, N., & Kenney, J. D. P. 2005, ApJ, 630, 837 [NASA ADS] [CrossRef] [Google Scholar]
- Jungwiert, B., Combes, F., & Axon, D. J. 1997, A&AS, 125, 479 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Keel, W. C. 1983, ApJS, 52, 229 [NASA ADS] [CrossRef] [Google Scholar]
- Knapen, J. H., Beckman, J. E., Heller, C. H., Shlosman, I., & de Jong, R. S. 1995a, ApJ, 454, 623 [NASA ADS] [CrossRef] [Google Scholar]
- Knapen, J. H., Beckman, J. E., Shlosman, I., et al. 1995b, ApJ, 443, L73 [NASA ADS] [CrossRef] [Google Scholar]
- Knapen, J. H., Shlosman, I., Heller, C. H., et al. 2000, ApJ, 528, 219 [NASA ADS] [CrossRef] [Google Scholar]
- Knapen, J. H., Pérez-Ramírez, D., & Laine, S. 2002, MNRAS, 337, 808 [NASA ADS] [CrossRef] [Google Scholar]
- Kormendy, J., & Ho, L. C. 2013, ARA&A, 51, 511 [Google Scholar]
- Laine, S., Shlosman, I., Knapen, J. H., & Peletier, R. F. 2002, ApJ, 567, 97 [Google Scholar]
- Lang, P., Meidt, S. E., Rosolowsky, E., et al. 2020, ApJ, 897, 122 [CrossRef] [Google Scholar]
- Leitherer, C., Li, I. H., Calzetti, D., & Heckman, T. M. 2002, ApJS, 140, 303 [NASA ADS] [CrossRef] [Google Scholar]
- Leon, S., Jeyakumar, S., Pérez-Ramírez, D., et al. 2008, A&A, 491, 703 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Leroy, A. K., Schinnerer, E., Hughes, A., et al. 2021, ApJS, 257, 43 [NASA ADS] [CrossRef] [Google Scholar]
- Lin, D. N. C., & Pringle, J. E. 1987, MNRAS, 225, 607 [Google Scholar]
- Lindblad, P. A. B., & Kristen, H. 1996, A&A, 313, 733 [NASA ADS] [Google Scholar]
- Lindblad, P. A. B., Kristen, H., Joersaeter, S., & Hoegbom, J. 1997, A&A, 317, 36 [NASA ADS] [Google Scholar]
- Maeda, F., Ohta, K., Fujimoto, Y., Habe, A., & Baba, J. 2018, PASJ, 70, 37 [NASA ADS] [Google Scholar]
- Maeda, F., Ohta, K., Fujimoto, Y., & Habe, A. 2020, MNRAS, 493, 5045 [NASA ADS] [CrossRef] [Google Scholar]
- Makarov, D., Prugniel, P., Terekhova, N., Courtois, H., & Vauglin, I. 2014, A&A, 570, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Maoz, D., Barth, A. J., Ho, L. C., Sternberg, A., & Filippenko, A. V. 2001, AJ, 121, 3048 [NASA ADS] [CrossRef] [Google Scholar]
- Martin, P., & Friedli, D. 1997, A&A, 326, 449 [NASA ADS] [Google Scholar]
- Martínez-García, E. E. 2012, ApJ, 744, 92 [CrossRef] [Google Scholar]
- Martini, P., Regan, M. W., Mulchaey, J. S., & Pogge, R. W. 2003, ApJS, 146, 353 [NASA ADS] [CrossRef] [Google Scholar]
- Matsuoka, Y., Yuan, F.-T., Takeuchi, Y., & Yanagisawa, K. 2012, PASJ, 64, 44 [NASA ADS] [CrossRef] [Google Scholar]
- Mazzalay, X., Saglia, R. P., Erwin, P., et al. 2013, MNRAS, 428, 2389 [NASA ADS] [CrossRef] [Google Scholar]
- Mazzalay, X., Maciejewski, W., Erwin, P., et al. 2014, MNRAS, 438, 2036 [NASA ADS] [CrossRef] [Google Scholar]
- Mazzuca, L. M., Knapen, J. H., Veilleux, S., & Regan, M. W. 2008, ApJS, 174, 337 [NASA ADS] [CrossRef] [Google Scholar]
- Mazzuca, L. M., Swaters, R. A., Knapen, J. H., & Veilleux, S. 2011, ApJ, 739, 104 [NASA ADS] [CrossRef] [Google Scholar]
- McAlpine, W., Satyapal, S., Gliozzi, M., et al. 2011, ApJ, 728, 25 [NASA ADS] [CrossRef] [Google Scholar]
- Menezes, R. B., Steiner, J. E., Ricci, T. V., & da Silva, P. 2022, MNRAS, 513, 5935 [NASA ADS] [CrossRef] [Google Scholar]
- Murugeshan, C., Kilborn, V., Jarrett, T., et al. 2020, MNRAS, 496, 2516 [NASA ADS] [CrossRef] [Google Scholar]
- Neumann, L., Gallagher, M. J., Bigiel, F., et al. 2023, MNRAS, 521, 3348 [NASA ADS] [CrossRef] [Google Scholar]
- Pereira-Santaella, M., Diamond-Stanic, A. M., Alonso-Herrero, A., & Rieke, G. H. 2010, ApJ, 725, 2270 [Google Scholar]
- Pérez-Ramírez, D., Knapen, J. H., Peletier, R. F., et al. 2000, MNRAS, 317, 234 [CrossRef] [Google Scholar]
- Pérez-Ramírez, D., Caballero-García, M. D., Ebrero, J., & Leon, S. 2010, A&A, 522, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Persic, M., Salucci, P., & Stel, F. 1996, MNRAS, 281, 27 [Google Scholar]
- Pierce, M. J. 1986, AJ, 92, 285 [NASA ADS] [CrossRef] [Google Scholar]
- Planesas, P., Colina, L., & Perez-Olea, D. 1997, A&A, 325, 81 [NASA ADS] [Google Scholar]
- Popping, G., Pérez, I., & Zurita, A. 2010, A&A, 521, A8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Regan, M. W., Sheth, K., & Vogel, S. N. 1999, ApJ, 526, 97 [NASA ADS] [CrossRef] [Google Scholar]
- Riffel, R. A., Colina, L., Storchi-Bergmann, T., et al. 2016, MNRAS, 461, 4192 [NASA ADS] [CrossRef] [Google Scholar]
- Rosado-Belza, D., Falcón-Barroso, J., Knapen, J. H., et al. 2020, A&A, 644, A116 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Saikia, D. J., Pedlar, A., Unger, S. W., & Axon, D. J. 1994, MNRAS, 270, 46 [NASA ADS] [CrossRef] [Google Scholar]
- Sakamoto, K., Okumura, S., Minezaki, T., Kobayashi, Y., & Wada, K. 1995, AJ, 110, 2075 [NASA ADS] [CrossRef] [Google Scholar]
- Sakamoto, K., Okumura, S. K., Ishizuki, S., & Scoville, N. Z. 1999, ApJS, 124, 403 [NASA ADS] [CrossRef] [Google Scholar]
- Sánchez-Gil, M. C., Jones, D. H., Pérez, E., et al. 2011, MNRAS, 415, 753 [CrossRef] [Google Scholar]
- Sarkar, K. C., Mondal, S., Sharma, P., & Piran, T. 2023, ApJ, 951, 36 [NASA ADS] [CrossRef] [Google Scholar]
- Satyapal, S., Vega, D., Dudik, R. P., Abel, N. P., & Heckman, T. 2008, ApJ, 677, 926 [NASA ADS] [CrossRef] [Google Scholar]
- Schinnerer, E., Maciejewski, W., Scoville, N., & Moustakas, L. A. 2002, ApJ, 575, 826 [NASA ADS] [CrossRef] [Google Scholar]
- Sérsic, J. L. 1973, PASP, 85, 103 [CrossRef] [Google Scholar]
- Shankar, F., Weinberg, D. H., & Miralda-Escudé, J. 2009, ApJ, 690, 20 [NASA ADS] [CrossRef] [Google Scholar]
- Shaw, M., Axon, D., Probst, R., & Gatley, I. 1995, MNRAS, 274, 369 [NASA ADS] [CrossRef] [Google Scholar]
- She, R., Ho, L. C., & Feng, H. 2017, ApJ, 842, 131 [NASA ADS] [CrossRef] [Google Scholar]
- Sheth, K., Vogel, S. N., Regan, M. W., et al. 2002, AJ, 124, 2581 [NASA ADS] [CrossRef] [Google Scholar]
- Shlosman, I., Frank, J., & Begelman, M. C. 1989, Nature, 338, 45 [Google Scholar]
- Shlosman, I., Begelman, M. C., & Frank, J. 1990, Nature, 345, 679 [Google Scholar]
- Simons, D. A., Depoy, D. L., Becklin, E. E., et al. 1988, ApJ, 335, 126 [NASA ADS] [CrossRef] [Google Scholar]
- Sodré, L. J. & Stasińska, G., 1999, A&A, 345, 391 [Google Scholar]
- Sofue, Y., Koda, J., Nakanishi, H., & Onodera, S. 2003a, PASJ, 55, 59 [NASA ADS] [Google Scholar]
- Sofue, Y., Koda, J., Nakanishi, H., et al. 2003b, PASJ, 55, 17 [CrossRef] [Google Scholar]
- Song, Y., Linden, S. T., Evans, A. S., et al. 2021, ApJ, 916, 73 [NASA ADS] [CrossRef] [Google Scholar]
- Stuber, S. K., Saito, T., Schinnerer, E., et al. 2021, A&A, 653, A172 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Stuber, S. K., Schinnerer, E., Williams, T. G., et al. 2023, A&A, 676, A113 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Sun, J., Leroy, A. K., Schinnerer, E., et al. 2020, ApJ, 901, L8 [NASA ADS] [CrossRef] [Google Scholar]
- Turnrose, B. E. 1976, ApJ, 210, 33 [NASA ADS] [CrossRef] [Google Scholar]
- Tzanavaris, P., & Georgantopoulos, I. 2007, A&A, 468, 129 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- van den Bosch, R. C. E. 2016, ApJ, 831, 134 [NASA ADS] [CrossRef] [Google Scholar]
- Veron-Cetty, M. P., & Veron, P. 1986, A&AS, 66, 335 [Google Scholar]
- Wada, K., Sakamoto, K., & Minezaki, T. 1998, ApJ, 494, 236 [NASA ADS] [CrossRef] [Google Scholar]
- Weinzirl, T., Jogee, S., Khochfar, S., Burkert, A., & Kormendy, J. 2009, ApJ, 696, 411 [NASA ADS] [CrossRef] [Google Scholar]
- Wynn-Williams, C. G., & Becklin, E. E. 1985, ApJ, 290, 108 [NASA ADS] [CrossRef] [Google Scholar]
- Yuan, C., & Kuo, C.-L. 1998, ApJ, 497, 689 [NASA ADS] [CrossRef] [Google Scholar]
- Yukita, M., Swartz, D. A., Tennant, A. F., Soria, R., & Irwin, J. A. 2012, ApJ, 758, 105 [CrossRef] [Google Scholar]
- Zabel, N., Brown, T., Wilson, C. D., et al. 2022, ApJ, 933, 10 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Comments about individual galaxies:
NGC 1300: is a very well known MW morphological twin with a prominent strong bar and two arms rich in star forming knots (Eskridge et al. 2002; Knapen et al. 2002). Spiral arms starting at the end of the bar, extend out to the outer Lindblad resonance (England 1989), responsible for the outer pseudo-ring (major axis: 5’.71– Buta & Crocker 1993, age: 109 yrs –Elmegreen et al. 1996). The inner Lindblad resonance is the site of a circumnuclear ring. Star formation is concentrated at the ends of the bar (Elmegreen et al. 1992), and, in the central part, there is not much star formation, no molecular gas concentration, and no clear HII regions (Knapen et al. 2002; Maeda et al. 2018), possibly due to the fast collision of gas clouds (Maeda et al. 2018, 2020; Fujimoto et al. 2020). There are many classifications for this galaxy nucleus. According to some works there is no AGN (e. g. Martini et al. 2003; Stuber et al. 2021); while others see evidence of AGN activity (outflow: Saikia et al. 1994, low luminosity AGN: Jang et al. 2014; Constantin et al. 2015; Gadotti et al. 2019; Menezes et al. 2022). In any case, there is evidence of gas infall into the inner galaxy (radius < 25" Lindblad & Kristen 1996), which harbours an SMBH with (Atkinson et al. 2005). No evidence of molecular gas outflow was detected according to Stuber et al. (2021). The inner ring is more continuous in the northern part, with distinct regions of star formation and there is no evidence of a nuclear bar (Erwin 2004, Pérez-Ramírez et al. 2000, and references therein). This ring appears complete in the NIR (Knapen et al. 2002) and its semi-major axis is 400 pc (Comerón et al. 2010), with age of 4.6 ± 0.8 Gyrs (Rosado-Belza et al. 2020). The stellar velocity maps from MUSE show rotation with the blueshift towards east and redshift towards west (Gadotti et al. 2020).
NGC 1792: is classified as a non-barred SA(rs)bc galaxy (de Vaucouleurs et al. 1991; Eskridge et al. 2002; Sun et al. 2020; Neumann et al. 2023; Stuber et al. 2023), with flocculent (Eskridge et al. 2002) chaotic (Garrison & Walborn 1974) spiral arms rich in star-forming regions (Dahlem et al. 1994; Eskridge et al. 2002). The galaxy is currently interacting with NGC 1808 (Dahlem 1992). Its nuclear region, even though showing [NeV] emission (that is commonly found in AGNs –Goulding & Alexander 2009; Goulding et al. 2010, is classified as non-active in many studies (Veron-Cetty & Veron 1986; Stuber et al. 2021; Menezes et al. 2022; Neumann et al. 2023). Pereira-Santaella et al. (2010) say that the [NeV] emission detection in this galaxy might be due to supernova remnants, planetary nebulae, Wolf-Rayet stars or that its HII region classification might be due to strong obscuration of the nuclear region or to a very low luminosity AGN. The mass of the SMBH in its centre was estimated as (Goulding et al. 2010). Observations at the H-band indicate the presence of a short nuclear bar (Jungwiert et al. 1997). There is no evidence of molecular gas outflows (Stuber et al. 2021).
NGC 2903: is a very well-studied SAB(rs)bc (de Vaucouleurs et al. 1991) galaxy with a strong (in the NIR –Regan et al. 1999; Sheth et al. 2002; Hernandez et al. 2005, in short wavelengths it is similar to a patchy spiral disk –Popping et al. 2010) and young bar (Dutil & Roy 1999, with age between 200 and 600 Myr – Leon et al. 2008) immersed in a flocculent disk (Canzian 1993), somehow in a very complex structure. Bright HII regions delineate the bar and the multi-arm structure of the galaxy (Knapen et al. 2002). This galaxy is a typical example of "hot-spot" galaxy (Sérsic 1973; Turnrose 1976; Bonatto et al. 1989), rich in individual young stellar clusters (Alonso-Herrero et al. 2001), and starburst galaxy (Turnrose 1976; Wynn-Williams & Becklin 1985; Simons et al. 1988; Jackson et al. 1989; Planesas et al. 1997; Pérez-Ramírez et al. 2010). The hot-spots near the nucleus seem part of the circumnuclear starburst ring (Planesas et al. 1997; Sakamoto et al. 1999; Pérez-Ramírez et al. 2000; Alonso-Herrero et al. 2001; Mazzuca et al. 2011), which has a radius of about 600 pc (Leitherer et al. 2002). The galaxy hosts a SMBH with a mass of (van den Bosch 2016), and is commonly classified as non-active (Casuso et al. 1990; Ho et al. 1997; Sodré 1999; Pereira-Santaella et al. 2010; Neumann et al. 2023), which is supported by no evidence of molecular gas outflows (Stuber et al. 2021), no significant [NeV] emission (Goulding & Alexander 2009), and low luminosity in X-rays (0.5-8.0 keV –Yukita et al. 2012). However, NGC 2903 has NIR emission-line ratios compatible with an AGN (Matsuoka et al. 2012) and new measurements of its X-ray luminosity indicate the presence of a variable AGN (Hoyer et al. 2024). The nuclear starburst might be due to the bar that drives the molecular gas into the centre (Jackson et al. 1991; Hernandez et al. 2005; Leon et al. 2008).
NGC 3059: is not a very well-studied galaxy. It is classified as SB(rs)bc; however some studies classify it as SBcd (e. g. Eskridge et al. 2002), weakly barred (Stuber et al. 2023) or with a small bar (Beck et al. 2002). There is no evidence of an inner ring in the literature (Stuber et al. 2023), and no AGN (Stuber et al. 2021; Goulding & Alexander 2009; Neumann et al. 2023).
NGC 4303: is a grand design SAB(rs)bc (de Vaucouleurs et al. 1991) double barred galaxy (Colina et al. 1997; Colina & Arribas 1999; Colina & Wada 2000, nuclear bar diameter: 0.2 kpc, primary bar diameter: 2.64 kpc– Erwin 2004). It has a young (2.5 -25 Myrs) circumnuclear ring (200-250 pc– Riffel et al. 2016), connected to the bar (which is possibly responsible for the feeding of the nuclear region and the ring– Colina et al. 1997), by two spiral arms, and to an inner bar (Colina et al. 1997; Colina & Wada 2000; Pérez-Ramírez et al. 2000). Schinnerer et al. (2002) studied the molecular emission of the central kiloparsec of NGC 4303 and found that inflowing motions of the primary bar of the galaxy are responsible for the complex structure of the centre of this galaxy, and that this circumnuclear ring might be in an early stage of formation. The nuclear activity of this galaxy was classified as Seyfert 2 or LINER or in-between (e. g. Colina & Arribas 1999; Argyle & Eldridge 1990), or HII region (Ho et al. 1997). According to Colina et al. (2002), there is a young (∼ 4 Myr) super star cluster coexisting with the low luminosity AGN in the nucleus (within 3.1 pc), and Tzanavaris & Georgantopoulos (2007), Jiménez-Bailón et al. (2003) could not distinguish its nuclear emission from a hidden AGN, a central stellar cluster, or a ultra-luminous X-ray source. NGC 4303 harbours an SMBH with a mass of (Davis et al. 2019). An [OIII] outflow was detected north-east of the nucleus by Colina & Arribas (1999) and no molecular gas outflow was detected (Stuber et al. 2021).
NGC 4321: is a very well studied MW morphological twin (SABbc– de Vaucouleurs et al. 1991). Due to its inclination, distance, and contrasted features, it was subject to many detailed studies. For instance, this galaxy has a clumpy circumnuclear ring, whose peak radius is 458 pc (Song et al. 2021). Such starburst ring is being maintained by the gravitational influence of the bar (Arsenault et al. 1988; Knapen et al. 1995a; Pérez-Ramírez et al. 2000) and its age is about 50 Myrs (the majority of stars are younger than 6 Myr – Sánchez-Gil et al. 2011). Knapen et al. (1995a) state that there are two Inner Lindblad Resonances in the nuclear region: one at 500 pc and the other at 1.3 kpc from the nucleus. According to these authors, most of star formation happens at the inner resonance (also see Sakamoto et al. 1995; Yuan & Kuo 1998; Knapen et al. 2000). The ring is connected by two arms leading the bar to the nucleus (Knapen et al. 1995b,a, 2002 and the younger stellar populations are located near the connecting points between the ring and the larger bar (Mazzuca et al. 2008). There is also a double bar structure (Pierce 1986; Shaw et al. 1995), the nuclear bar being clearly detected in the NIR and the primary bar being weaker (Font et al. 2014) or moderate (Knapen et al. 1995a). Although Knapen et al. (1995a) do not found strong evidence for the second bar, they say both structures are part of the same bar, aligned with comparable strength; Garcia-Burillo et al. (1998) dynamical models favor that this nuclear bar is decoupled from the major bar and it is rotating faster. There are two dust arms connecting the nuclear bar to the outer bar (Knapen et al. 1995a; Sakamoto et al. 1995), and Wada et al. (1998) say that this nuclear structure is short lived (107 yrs) and recurrent. NGC 4321 nucleus was classified as HII region (Keel 1983, supported by González-Martín et al. 2009 that have not found any AGN activity); however, some studies classify it as a transition object (Ho et al. 1997) and the emission of [NeV] detected by Satyapal et al. (2008) sustains the hypothesis of a low luminosity AGN. The SMBH in NGC 4321 centre has a mass of 6.67+0.17 -6.67 (van den Bosch 2016). The high nuclear CO concentration is a result of the gas being driven by gravitational torques from the bar; however the fact that there is no strong AGN activity at the centre of this galaxy shows that this inflow is still at intermediate scale, and is not yet prolonged into the center to trigger nuclear activity (Sakamoto et al. 1995). Such remark is actually proven by the detailed torque analysis made by García-Burillo et al. (2005): there are positive torques along the spiral arms at the trailing edges of the nuclear bar and negative torques along the leading quadrants and outer disk (r > 900). They also noticed that the average torque inside 200 pc is positive, implying that there is no significant inflow of molecular gas in the nucleus. Castillo-Morales et al. (2007) detected a nuclear blueshifted region that might be related to an outflow of gas. Such outflow was also reported by Jiménez-Vicente et al. (2007) and the authors say that such wind was generated by a nuclear starburst.
NGC 4536: is a SAB(rs)bc (de Vaucouleurs et al. 1991) galaxy with two conspicuous spiral arms (Mazzalay et al. 2014). This galaxy was previously classified as starburst (Ho et al. 1997); however Hughes et al. (2005) classified it as a weak AGN, based on optical emission-line ratios. There is [NeV] emission, detected 10" north-east from the optical nucleus (Satyapal et al. 2008), and XMM-Newton data presents a power law emission typical of low luminosity AGNs (McAlpine et al. 2011), that supports the presence of a low luminosity AGN in this nucleus. A circumnuclear ring with radius of about 270 pc with age of 6.5 Myr was detected with near-infrared data, with no gradients along the ring (Mazzalay et al. 2013, 2014). There is no evidence of molecular gas outflow (Stuber et al. 2021) and, in the inner 500 pc, the gas kinematics is predominantly circular with streams of inflowing gas along bar into the centre (Jogee et al. 2005). The nucleus is connected to the circumnuclear ring and there is evidence of a nuclear bar (Mazzalay et al. 2014). The star formation is predominant in the ring and it has ceased in the nucleus (Davies et al. 1997; Jogee et al. 2005). If there is a low luminosity AGN in this galaxy, there is no influence of its emission in the circumnuclear ring (Mazzalay et al. 2013). The SMBH mass was estimated to be (Berrier et al. 2013).
NGC 4689: is a Milky-Way analogue classified non-barred SA(rs)bc (de Vaucouleurs et al. 1991). The galaxy has flocculent spiral arms (Jungwiert et al. 1997); however, Arsenault (1989) says that it has a weak bar. Its nucleus is classified as starburst (Ho et al. 1997), confirmed by the non- detection of X-ray emission, with a SMBH mass of log(MBH/M⊙) = 4.2 ± 0.8 (Graham et al. 2019). No molecular outflow was detected (Stuber et al. 2021). Molecular gas shows patchy morphology without clear structure (bar or spiral), with a regular velocity field (Sofue et al. 2003b,a). Two pseudo rings were found, one located at the bar Inner Lindblad Resonance, 0’.57 ×0’.44, and other located at the Outer Lindblad Resonance, with 1’.86 ×1’.51 (Comerón et al. 2014).
NGC5248: is an SAB(rs)bc (de Vaucouleurs et al. 1991) galaxy with strong spiral arms (Elmegreen et al. 1992). Within our field of view, there is evidence of two rings, one inner ring of 650pc diameter around the bar and a double nuclear ring (Buta & Crocker 1993; Comerón et al. 2010). The bar fuels the nuclear ring with diameter of 150 pc (Comerón et al. 2010) with two molecular trailing spirals, which connect to other two arms (Jogee et al. 2002) (Buta & Crocker 1993; Maoz et al. 2001; Comerón et al. 2010). The nuclear rings have stellar population ages between 1 and 10 Myr (Elmegreen et al. 1997). Its nucleus was classified as starburst (Ho et al. 1997) and Seyfert/LINER (Gadotti et al. 2019) from the optical emission-line ratios, and harbours a SMBH with mass of (She et al. 2017). There is no evidence of molecular gas outflows (Stuber et al. 2021). Haan et al. (2009) studied the torques of the CO molecular gas and found that NGC 5248 shows a mixed pattern of gravitational torques that corresponds to the bar and spiral arms, comparable to our results.
NGC 5921: is a SB(r)bc (de Vaucouleurs et al. 1991) galaxy with a compact nucleus and a prominent bar that ends in an outer ring (Eskridge et al. 2002; Fathi et al. 2009). Its nucleus is classified as composite (Ho et al. 1997), LINER and Seyfert-like (Veron-Cetty & Veron 1986; Arsenault 1989), and non-active (Burtscher et al. 2021). This galaxy has an inner ring-like structure (Hernandez et al. 2005; Martínez-García 2012). Its SMBH estimated mass is (van den Bosch 2016).
Appendix B: Hα gas torques with MUSE data
![]() |
Fig. B.1. Torque maps of Hα emission line from MUSE data cube of: (a)NGC 1300 with contours levels of 1 to 6 times 2.4 ×10−17 erg Å−1 cm−2 s−1, (b)NGC 4303 with contours levels of 1, 2.5, 6, and 14.3 times 1.75 ×10−17 erg Å−1 cm−2 s−1, (c) NGC 4321 with contours levels of 1, 2, 3.7, 7.1, 13.3 and 25.1 times 5.4 ×10−18 erg Å−1 cm−2 s−1, and (d) NGC 5248 with contour levels of 1, 1.8, 3.2 and 5.7 times 5.7 ×10−17 erg Å−1 cm−2 s−1. |
For the sample of 10 galaxies of this work, only four galaxies (NGC 1300, NGC 4303, NGC 4321 and NGC 5248) had Multi Unit Spectroscopic Explorer (MUSE) data cubes available. The Hα emission line can trace the star forming regions and interstellar medium in the less dense regions, giving a more complete image of the structures of the galaxy. In contrast, the molecular gas emission is more compact, giving a overall idea of the structures punctually. In order to see the distribution of torques in general in the galaxies and compare them with the obtained molecular gas torques, we made images of Hα emission line from MUSE data cubes, due to the good spatial resolution, and FOV, to compute the torques on this gas. The data were taken from the ESO public archive. The information about the observations of each galaxy is shown in table B.1.
The Hα emission-line images were taken from the data cubes by summing the spectral interval of the emission line and by subtracting the nearby continuum equivalent to the same interval. After that, the torques were computed associated to the red
Observational information of MUSE data from NGC 1300, NGC 4303, NGC 4321, and NGC 5248.
images based on the size of the FOV. In the case of NGC 1300, we used HST (cyan square in Fig. 1 or Fig. 2b), in the case of NGC 4303, the Legacy image was used (red square in Fig. 1 or Fig. 6b), for NGC 4321, we also used the Legacy image (red square in Fig. 1 or Fig. 7b), and for NCG 5248 we used the HST image (cyan square in Fig. 1 or Fig. 12b). The FOV of the MUSE data was cut to be the same size of the associated FOVs. We then computed the torques analogously to what was done with the ALMA data.
All Tables
Maximum and minimum values with their associated radii of the torque and relative maxima of non-axisymmetric potential curves (Q2 and QT).
Observational information of MUSE data from NGC 1300, NGC 4303, NGC 4321, and NGC 5248.
All Figures
![]() |
Fig. 1. Composite of Legacy images of the G (blue), R (green), and Z (red) bands (except for NGC 3059, which lacks Legacy images; in this case, we used the images from the Digitalized Sky Survey). The squares represent the size of the FOV that was used for each data set: green for HST and ALMA, cyan for the HST and the ALMA high spatial resolution set, red for Legacy and ALMA, and yellow for the Legacy and ALMA larger FOV (combination of two FOVs). The pink lines represent the PA of the bar. The orientation N–E is north up and east left. |
In the text |
![]() |
Fig. 2. NGC 1300. Left panel: Images of the F814W filter of HST (a, b), and Legacy Z band (c, d). Middle and right panels: ALMA data (mean flux and velocity map) associated with each FOV, as described in section 3. The orientation N–E is north up and east right. The cross (arbitrary size) indicates the centre of the image, which was taken from the emission peak from the infrared image and was used to match the ALMA images. |
In the text |
![]() |
Fig. 3. NGC 1792. Same as the left panel of Fig. 2. (a) Image of the HST filter F814W and (b) image of the Z band from the Legacy survey. |
In the text |
![]() |
Fig. 4. NGC 2903. Same as Fig. 2. |
In the text |
![]() |
Fig. 5. NGC 3059. Same as Fig. 2. The HST image is shown in the left panel. |
In the text |
![]() |
Fig. 6. NGC 4303. Same as Fig. 3. In this case, an R-band image from the Legacy survey was used. |
In the text |
![]() |
Fig. 7. NGC 4321. Same as Fig. 3. |
In the text |
![]() |
Fig. 8. NGC 4536. Same as Fig. 3. In this case, an R-band image from the Legacy survey was used. |
In the text |
![]() |
Fig. 9. NGC 4689. Same as Fig. 3. |
In the text |
![]() |
Fig. 10. NGC 5248. Same as Fig. 2. Panels (a) and (b) show images of the HST filter F814W, and panel (c) shows an image from the Z band from the Legacy survey. |
In the text |
![]() |
Fig. 11. NGC 5921. Left panel: F814W image from HST. Middle and right panels: ALMA data as in Fig. 2. |
In the text |
![]() |
Fig. 12. Radial distribution of the amplitudes (top) and phases (bottom) of the non-axisymmetric perturbations in each Fourier component (m = 1 to 4). QT represents the total perturbation curve in the potential for the data set of N4303 (b). This plot illustrates the typical potential component curves obtained in the sample. The maximum values of the Q2 and QT curves are shown in Table 3. |
In the text |
![]() |
Fig. 13. NGC 1300 torque maps from the same set of data as in panel (a) Fig. 2a, (b) Fig. 2b, (c) Fig. 2c, and (d) Fig. 2d. The images were rotated and de-projected. The orientation N–E is indicated in the torque maps. The crosses represent the centre of the galaxy estimated from the emission peaks of the infrared images and ALMA data. The dark violet line represents the PA of the bar after rotation. Each map is associated with a graph that represents the radial distribution of the torque, averaged over the azimuth, and normalised to the angular momentum, per rotation at this radius. The colour scales in the maps are in arbitrary units (except for panels c and d), of the torques t(x,y) weighted by the gas surface density Σ(x,y). The contours of Σ(x,y) are from the mean molecular gas map (a) one to nine times 18.55 Jy/beam km/s, (b) 1, 1.1, 1.4, 2.5, 5.9, 16.5 times 0.3 Jy/beam km/s (c) 1, 5, 11, 23.8, 51.5 times 1.65 Jy/beam km/s, and (d) 1, 2.5, 12.5, 27.2, 58.9, 127.2 times 0.49 Jy/beam km/s. The velocities are relative to the centre. |
In the text |
![]() |
Fig. 14. NGC 1792 torque maps from the same set of data as in panel (a) Fig. 3a and (b) Fig. 3b. Same caption as Fig. 13. The contour levels were taken from the mean molecular gas map (a) one to six times 34.13 Jy/beam km/s and (b) one to five times 32.11 Jy/beam km/s. |
In the text |
![]() |
Fig. 15. NGC 2903 torque maps from the same set of data as in panel (a) Fig. 4a, (b) Fig. 4b, (c) Fig. 4c, and (d) Fig. 4d. Same caption as Fig.13. The contour levels were taken from the mean molecular gas map (a) 1, 4, 9, 16, 25, 36, and 49, times 6.8 Jy/beam km/s, (b) one to seven times 1.08 Jy/beam km/s, (c) one to five times 1.43 Jy/beam km/s, and (d) 1, 4, 9, 16, 25, and 36 times 6.3 Jy/beam km/s. |
In the text |
![]() |
Fig. 16. NGC 3059 torque maps from the same set of data as in Fig. 5. Same caption as Fig. 13. The contour levels were taken from the mean molecular gas map 1, 2, 3, 4, 6, 8, 10, and 11 times 15.7 Jy/beam km/s. |
In the text |
![]() |
Fig. 17. NGC 4303 torque maps from the same set of data as in panel (a) Fig. 6a and (b) Fig. 6b. Same caption as Fig. 13. The contour levels were taken from the mean molecular gas map (a) 1, 2, 4, 5, 7, 9, 10, and 12 times 14.6 Jy/beam km/s, and (b) 1, 2, 5, 6, 8, 9, and 10 times 15.9 Jy/beam km/s. |
In the text |
![]() |
Fig. 18. NGC 4321 torque maps from the same set of data as in panel (a) Fig. 7a and (b) Fig. 7b. Same caption as Fig. 13. The contour levels were taken from the mean molecular gas map (a and b) 1, 2, 5, 10, 23, 34, and46 times 4.4 Jy/beam km/s. |
In the text |
![]() |
Fig. 19. NGC 4536 torque maps from the same set of data as in panel (a) Fig. 8a and (b) Fig. 8b. Same caption as Fig. 13. The contour levels were taken from the mean molecular gas map (a) 1, 4.2, 14.2, 45.8, and 145.8 times 1.56 Jy/beam km/s, and (b) 1, 1.1, 1.4, 2, 3.1, 5.5, 10, 19.3, 37.6, and 74.1 times 5 Jy/beam km/s. |
In the text |
![]() |
Fig. 20. NGC 4689 torque maps from the same set of data as in panel (a) Fig. 9a and (b) Fig. 9b. Same caption as Fig. 13. The contour levels were taken from the mean molecular gas map (a) 1, 1.4, 1.7, 2, 2.2, and 2.5 times 6.8 Jy/beam km/s, and (b) one to three times 4.9 Jy/beam km/s. |
In the text |
![]() |
Fig. 21. NGC 5248 torque maps from the same set of data as in panel (a) Fig. 10a, (b) Fig. 10b, and (c) Fig. 10c. Same caption as Fig. 13. The contour levels were taken from the mean molecular gas map (a) 1, 1.4, 1.7, 2, 2.2, 2.5, and 2.7 times 59.3 Jy/beam km/s, and (b) one to five times 0.25 Jy/beam km/s. |
In the text |
![]() |
Fig. 22. cont. of Fig. 21. The contour levels were taken from the mean molecular gas map (c) 1, 4, 9, 16, 25, and 36 times 3.48 Jy/beam km/s. |
In the text |
![]() |
Fig. 23. NGC 5921 torque maps from the same set of data as in Fig. 11. Same caption as Fig. 13. The contour levels were taken from the mean molecular gas map one to seven times 0.74 Jy/beam km/s. |
In the text |
![]() |
Fig. 24. Gas velocity rotation curves from HI (Lindblad et al. 1997) (red squares) and from CO (Lang et al. 2020) (red points) of NGC 1300 used to compute the torque map values of the data sets in panels (c) and (d). Only the highest points of each curve were kept. The blue curve represents the velocity model obtained from the computed potential (with limitations of the red image cut-off beyond 5 kpc). |
In the text |
![]() |
Fig. B.1. Torque maps of Hα emission line from MUSE data cube of: (a)NGC 1300 with contours levels of 1 to 6 times 2.4 ×10−17 erg Å−1 cm−2 s−1, (b)NGC 4303 with contours levels of 1, 2.5, 6, and 14.3 times 1.75 ×10−17 erg Å−1 cm−2 s−1, (c) NGC 4321 with contours levels of 1, 2, 3.7, 7.1, 13.3 and 25.1 times 5.4 ×10−18 erg Å−1 cm−2 s−1, and (d) NGC 5248 with contour levels of 1, 1.8, 3.2 and 5.7 times 5.7 ×10−17 erg Å−1 cm−2 s−1. |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.