SDSS-IV MaNGA: the indispensable role of bars in enhancing the central star formation of low- z galaxies

We analyse two-dimensional maps and radial proﬁles of EW(H α ), EW(H δ A ), and D n (4000) of low-redshift galaxies using integral ﬁeld spectroscopy from the MaNGA survey. Out of ≈ 1400 nearly face-on late-type galaxies with a redshift z < 0.05, we identify 121 “turnover” galaxies that each have a central upturn in EW(H α ), EW(H δ A ), and/or a central drop in D n (4000), indicative of ongoing/recent star formation. The turnover features are found mostly in galaxies with a stellar mass above ∼ 10 10 M (cid:4) and NUV – r colour less than ≈ 5. The majority of the turnover galaxies are barred, with a bar fraction of 89 ± 3 per cent. Furthermore, for barred galaxies, the radius of the central turnover region is found to tightly correlate with one-third of the bar length. Comparing the observed and the inward extrapolated star formation rate surface density, we estimate that the central SFR have been enhanced by an order of magnitude. Conversely, only half of the barred galaxies in our sample have a central turnover feature, implying that the presence of a bar is not sufﬁcient to lead to a central SF enhancement. We further examined the SF enhancement in paired galaxies, as well as the local environment, ﬁnding no relation. This implies that the environment is not a driving factor for central SF enhancement in our sample. Our results reinforce both previous ﬁndings and theoretical expectation that galactic bars play a crucial role in the secular evolution of galaxies by driving gas inﬂow and enhancing the star formation and bulge growth in the centre.


I N T RO D U C T I O N
The bulges of disc-dominated galaxies are more complicated than previously thought. Initially believed to be similar to elliptical galaxies, and thus formed via violent-merging events and dominated by old stellar populations (e.g. Whitford 1978;Renzini 1999), a large fraction of pseudobulges were later identified, built up from discs via bar-or spiral arm-driven processes (e.g. Peletier & Balcells 1996;Helfer et al. 2003;Kormendy & Kennicutt 2004). Several studies then suggested that bulges can be composite systems in which a classical bulge and a pseudobulge coexist, following a complex formation and evolutionary scenario (e.g. Athanassoula 2005;Gadotti 2009; Kormendy & Barentine 2010;Nowak et al. 2010;Méndez-Abreu et al. 2014;Erwin et al. 2015). Alternatively, E-mail: linlin@shao.ac.cn (LL); cli2015@tsinghua.edu.cn (CL) direct observations of high-redshift galaxies suggest that bulges are formed through the inward migration and coalescence of clumps. These massive gas-rich clumps can sustain the rejuvenation of bulges, and act together with stellar migrations, minor mergers, and other dynamical effects to finally build up bulges in the following several billion years (e.g. Noguchi 1999;Carollo et al. 2007;Roškar et al. 2012;Bird et al. 2013;Mandelker et al. 2014Mandelker et al. , 2017. Tracing the central star formation of galaxies is key to understand how bulges build up their stellar masses. Any mechanism that causes disc instabilities and gas inflows may trigger central star formation, thus contributing to bulge growth. In the local Universe, the merger rate decreases as the Universe expands. Galaxy evolution is thus governed by long-term and secular processes. Internal mechanisms such as bars and spiral arms (and associated resonances) have been proposed to transfer angular momentum to the outer parts of discs, driving gas inflows towards the central regions (e.g. Athanassoula 2003). External processes like minor mergers can also cause disc instabilities that lead to gas flows towards the centres of the galaxies, eventually building up bulge components (e.g. Kormendy & Kennicutt 2004).
Many observations in the local Universe suggest that the existence of a bar-like structure in a galaxy contributes to cold gas inflows and the enhancement of central star formation (e.g. Sheth et al. 2005;Ellison et al. 2011;Wang et al. 2012). Barred galaxies are indeed found to have higher molecular gas concentrations and systematically higher central star formation rates (SFRs) than unbarred galaxies (e.g. Sakamoto et al. 1999;Sheth et al. 2005). Ellison et al. (2011) used Sloan Digital Sky Survey (SDSS) spectra to estimate the central SFRs of galaxies and compared them to global SFRs at given stellar mass, finding that barred galaxies have systematically higher central star formation than unbarred galaxies. Coelho & Gadotti (2011) found that barred galaxies have a bimodal bulge age distribution, indicating an excess of young stars.
Bars are, however, thought to suppress the global or galaxyscale star formation as well. Bar-induced star formation exhausts the infalling gas quickly, helping to build-up a central bulge, that can, in turn, stabilize the disc and stop further gas falling in (e,g. Bournaud & Combes 2002;Athanassoula 2003). Observations show that barlike structures are preferentially found in massive, gas-poor galaxies rather than in gas-rich galaxies (e.g. Masters et al. 2012;Cheung et al. 2013;Cervantes Sodi 2017). Wang et al. (2012) thus found that galaxies with strong bars can either enhance central star formation or suppress it, and the stellar bars have been proposed to play an important role in star formation quenching of local isolated galaxies (e.g. Masters et al. 2011;Gavazzi et al. 2015;George, Subramanian & Paul 2019;Newnham et al. 2020).
Despite the above, the connection between bars and central nuclear activity is still unclear. Numerous studies suggest that bars are not related to the active galactic nucleus (AGN) (e.g. Cheung et al. 2015;Cisternas et al. 2015), while others claim there is a weak connection (e.g. Knapen, Shlosman & Peletier 2000;Laine et al. 2002;Galloway et al. 2015). Nuclear, or inner bars would be needed to take material from the galactic scale bar to the supermassive black hole scale, which could explain this mismatch.
Apart from bar-induced central star formation, galaxy interactions/disruptions are another well-known mechanism for enhancing star formation, as strong starbursts are predominately found in merging systems. Nevertheless, several studies suggest that enhanced SFRs may not be ubiquitous in interacting galaxies. The average SFR enhancement in pair galaxies is below a factor of 2-3 (e.g. Lin et al. 2007Lin et al. , 2008Ellison et al. 2008;Knapen, Cisternas & Querejeta 2015). Many factors, such as the separation of the companions, gas fraction, mass ratio of the systems or even the merging geometry affect the star formation enhancement (e.g. Ellison et al. 2008;Davies et al. 2015;Pan et al. 2018). Projected separation of less than ∼100 kpc is required for central enhancement (Li et al. 2008).
The recent development in the field of integral-field spectrographs (IFSs) provides us with a unique opportunity to revisit the detailed morphological substructures in the inner regions of galaxies, including bulges and bars. In a previous work (Lin et al. 2017), we used the 4000Å break (D n (4000)), the equivalent width of Hδ absorption (EW(Hδ A )) and Hα emission (EW(Hδ A )) to reveal the recent star formation history of 57 galaxies selected from Calar Alto Legacy Integral Field Area Survey (CALIFA; Sánchez et al. 2012). By analysing their two-dimensional maps and radial profiles, we identified a class of "turnover" galaxies that each have a significant inner drop in D n (4000) and/or a corresponding upturn in EW(Hα) and EW(Hδ A ). This indicates that the central regions of these galaxies have experienced star formation in the past 1-2 Gyr. We found most turnover galaxies to host a bar structure, but only half of the barred galaxies to have a central turnover. In the following work of Chown et al. (2019), we studied the link between star formation enhancements and molecular gas surface densities. Turnover galaxies are associated with higher molecular gas concentrations, supporting bar-or interaction-driven gas inflow scenarios. However, due to the limited sample size in that study, turnover features in recent star formation indicators are still not fully understood. It is therefore necessary to extend previous work with a larger sample covering a wider parameter space and enabling us to build a control sample.
In this paper, we use the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA; Bundy et al. 2015) data set to accomplish this, and confirm that there is a strong link between turnover features and the bars. This paper is organized as follows. Section 2 describes our data and sample selection, as well as our measurements of bar properties, the environment catalogue we adopted, and how we identify turnover features. We present the main result in Section 3. The discussion and conclusions are addressed in Section 4 and 5, respectively. Throughout this paper, we adopt a flat Lambda cold dark matter cosmology, with parameters = 0.3, = 0.7, and H 0 = 70 km s −1 Mpc −1 .

Overview of MaNGA
As one of the three key programs of the fourth-generation of SDSS (SDSS-IV; Blanton et al. 2017), MaNGA is an integral-field unit (IFU) spectroscopic survey (Bundy et al. 2015), aiming to achieve kpc-scale spatially resolved IFU datacubes for about 10 000 nearby galaxies. The MaNGA instrument has 29 IFUs including 12 sevenfibre mini-bundles for flux calibration and 17 science bundles with multiple fields of view ranging from 12 to 32 arcsec on sky (Drory et al. 2015). The spectra are obtained using the two dual-channel BOSS spectrographs at the 2.5-m Sloan Telescope (Gunn et al. 2006;Smee et al. 2013), covering wavelengths from 3622 to 10 354 Å with a spectral resolution of R ≈ 2000, and reaching a target r-band signalto-noise ratio SNR =4-8 (per Å per 2-fibre arcsec) at 1-2 effective radius (R e ) with a typical exposure time of three hours.
MaNGA target galaxies are selected from NASA-Sloan Atlas catalogue 1 (NSA; Blanton et al. 2011), which is an SDSS-based galaxy catalogue including physical parameters of 640 000 galaxies from GALEX, SDSS, and 2MASS. The target sample is well defined to have a statistically significant size, a uniform spatial coverage in units of R e , and an approximately flat total stellar mass distribution between 10 9 and 10 11 M . Details of the sample design and optimization can be found in Wake et al. (2017). The MaNGA sample includes three subsamples. The primary and secondary samples are selected to have a flat distribution in i-band absolute magnitude, and a spatial coverage out to 1.5 and 2.5R e , respectively. The third colour enhanced sample increases the fraction of rare populations in the galaxy colour-magnitude diagram, such as high-mass blue galaxies and low-mass red galaxies. The median redshift of the MaNGA galaxies is z ≈ 0.03.
The MaNGA observing strategy, survey execution and flux calibration are described in Yan et al. (2016aYan et al. ( , 2016b. The absolute flux calibration is better than 5 per cent for more than 80 per cent of the wavelength range. Raw data are reduced with the Data Reduction Figure 1. Distribution of MPL-7 galaxies on the diagram of redshift versus the r-band absolute magnitude. The spatial resolution, corresponding to the physical scale of a typical PSF FWHM of 2.5 arcsec, is also shown along the ordinate. The galaxies populate two branches, corresponding to the MaNGA primary and secondary samples. To ensure sufficient spatial resolution, we limit our sample galaxies at redshifts z < 0.05 with axial ratios b/a > 0.5. Pipeline (DRP; Law et al. 2016), which produces a datacube for each galaxy with a spaxel size of 0.5 arcsec and an effective spatial resolution of ≈2.5 arcsec. The MaNGA team has developed a Data Analysis Pipeline (DAP) that performs full spectral fitting to the spectra in DRP datacubes and calculates a variety of physical parameters including stellar kinematics, emission-line properties and spectral line indices (Belfiore et al. 2019;Westfall et al. 2019) In this work, we use the internal data release of MaNGA, the MaNGA Product Launch 7 (MPL-7), that includes DRP and DAP products for 4672 galaxies and is identical to the MaNGA data released with the SDSS Data Release 15 (DR15; Aguado et al. 2019).

Sample selection
In this work, we will focus on the central regions of galaxies, thus requiring a statistical sample of data with sufficient spatial resolution. Fig. 1 shows the distribution of the MPL-7 galaxies on the diagram of redshift versus the r-band absolute magnitude. The spatial resolution in kpc is plotted on the ordinate, corresponding to the 2.5 arcsec full width at half-maximum (FWHM) of the MaNGA effective point spread function (PSF). The galaxies are located in two separate sequences in the figure, corresponding to MaNGA Primary and Secondary samples described above.
For our study, we select a sample of galaxies from the MPL-7 as follows. First, we select nearly face-on galaxies, with minor-tomajor axial ratio b/a > 0.5, where the semimajor (a) and semiminor (b) axes are taken from the NSA and are measured at the 25 mag arcsec −2 isophote in the r-band. Secondly, we restrict ourselves to relatively low redshifts z < 0.05. As can be seen from Fig. 1, this redshift cut corresponds to a worst spatial resolution of ≈2.3 kpc, comparable to the typical size of galactic bulges and bars (Gadotti 2011). This choice of the upper limit on redshift is a trade-off between having a substantially large sample for good statistics and ensuring a relatively high spatial resolution for measuring central structures. More discussion on the spatial resolution can be found in Appendix A. Third, we select late-type galaxies by requiring T-type > 0. The morphological classification was taken from Domínguez . Finally, we make use of the MPL-7 version of the close pair galaxy catalogue (see Section 2.4 for details) to exclude 78 galaxies that are merging systems at the final coalescence stage. These galaxies do not have regular morphologies/structures/centres, and it is not straightforward to define radial profiles and make comparisons to other types of galaxies.
These restrictions give rise to a sample of 1452 late-type galaxies, which forms the main sample to be studied below. The numbers of galaxies in different categories are listed in Table 1. In this work, we use data produces provided by the MaNGA DAP. Specifically, we use the "hybrid" binning (HYB10) products, for which the stellar spectrum (continuum plus absorption lines) and stellar kinematics are obtained by applying the Penalized PiXel-Fitting (pPXF) code (Cappellari 2017) to spectra binned to S/N ∼ 10. For each spaxel in the datacube, the best-fitting stellar spectrum is subtracted from the observed spectrum and the starlight-subtracted spectrum is then used to measure emission lines by fitting each line with a Gaussian profile (Belfiore et al. 2019). We remove spaxels that might be contaminated by foreground stars or have critical failures during data reduction, using the same bitmask as adopted in the DAP.
Based on the EW(Hα), EW(Hδ A ), and D n (4000) maps from the MaNGA DAP, we estimate the one-dimensional radial profiles of these quantities by adopting a constant spatial sampling of 0.5 arcsec along the major axis, with fixed ellipticity and position angle (PA) and thus perfectly concentric elliptical annuli for all radial bins. The global ellipticity, PA of the major axis, and the major-to-minor axial ratio are taken from the NSA catalogue. Uncertainties on the derived quantities are given by the 1σ scatter of the spaxels within each radial annulus. We note that the one-dimensional profile is mixing non-axisymmetric structures, so it is important to examine the twodimensional map simultaneously.

Bar classification and measurement
We identify galactic bars in our sample and measure their lengths and ellipticities by applying the commonly used Image Reduction and Analysis Facility (IRAF) task ELLIPSE to the backgroundsubtracted r-band SDSS images. The ELLIPSE task fits twodimensional ellipses to the isophotes for each galaxy elliptical annuli, at logarithmically spaced radii along a galaxy major axis. For each galaxy in our sample, we adopt the sampling radius with a step of 1.1, and the centre of each ellipse is allow to vary. We apply the ELLIPSE task twice, adopting R 50 and R 90 , which taken from the NSA catalogue, as the starting radius, respectively. R 50 and R 90 are the elliptical Petrosian radii enclosing 50 and 90 per cent of the total light in the r-band. Then the two profiles are merged in order to fully cover the radial profile. Following Lin et al. (2017), we identify the bar structure according to the presence of an abrupt decrease of the ellipticity profile and an associated change of PA -a galaxy is classified as having a bar if the ellipticity increases above 0.25 as the radius increase from the galactic centre outwards, before decreasing by at least 0.1 at some radius. Meanwhile, the PA is required to change more than 10 • from the end of the bar to the outer disc region.
We visually examined the SDSS gri images of all 1425 galaxies in our sample. As noted by Barazza, Jogee & Marinova (2008) and Wang et al. (2012), some barred galaxies (though few in number) may be missed by the ellipse-fitting procedure, probably due to the existence of a bulge or the alignment of the bar with the minor axis at a moderate inclination angle (making the bar appear rounder and smoother). In these cases, the ellipticity would be underestimated and so would not fulfil the requirement of ellipticity > 0.25. Galaxies in  which the bar smoothly connects with the spiral arms or ring would also be missed. After visually examining all the images, we include 15 additional galaxies with a weak bar, 6 per cent of our final barred galaxy sample.
Our final barred galaxy sample consists of 252 galaxies. The halflength (r bar in arcseconds), ellipticity (e bar ) and PA (PA bar in degrees) of each bar is determined from the best-fitting ellipse at the semimajor axis radius at which the ellipticity reaches a local maximum. 2 We further calculate the deprojected bar length and ellipticity following Gadotti et al. (2007). We adopt uncertainties of 10 per cent in both the bar length and ellipticity measurements, taking into account the uncertainties due to different bar measurement methods (Zou, Shen & Li 2014). Our bar measurements are listed in Table 2.
The bar fraction in the main sample is 18 per cent (252/1425), which seems much lower than the known fraction of one-third for strong bars or the fraction of two-thirds for all bars. We note that the bar fraction is a strong function of stellar mass and sSFR. After the selection criteria, our main sample is restricted to late-type galaxies 2 The bar length determined by the location of the maximum ellipticity is thought to be a lower limit (Erwin 2005). However, it is well correlated with the true value. The correlation we find in the following should not be affected. and have a median stellar mass of ∼10 10 M . The bar fraction in this type of galaxies is ∼20 per cent in Cheung et al. (2013), close to our bar fraction. In addition, our bar fraction in the control sample is 34 per cent (41/121). The increase in the bar fraction can be understood from the larger stellar mass in the control sample. Overall, our bar fraction is consistent with the results from the more general population.

Environment catalogues
To investigate potential environmental effects on our galaxies at different physical scales, we make use of the MaNGA galaxy pair catalogue of Pan et al. (2019), the SDSS galaxy group catalogue of Yang et al. (2007), and the local environment density measure inferred from the density field reconstructed by Wang et al. (2016). We briefly describe these catalogues below.
The close galaxy pair catalogue was constructed using the MPL-7 (or SDSS DR15) version of the MaNGA sample as described in detail in Pan et al. (2019). A galaxy pair is defined as a system of two galaxies with a projected separation less than 50 kpc h −1 and a line-of-sight velocity difference less than 500 km s −1 . The pairs are further classified into four different merger stages based on a visual inspection of their SDSS images. Galaxies in pairs at Stage 1 are those without any morphology distortion. Galaxies in pairs at Stage 2 have obvious tidal tails or bridges. Galaxies in pairs at Stage 3 have only weak morphology distortions. Stage 4 encompasses galaxy mergers close to the final coalescence stage, for which the two galaxies are largely overlapping, meaning it is hard to measure radial profiles. Therefore, we exclude Stage 4, in this study and focus on the well-separated pairs at Stages 1-3.
The galaxy group catalogue of Yang et al. (2007) has been widely used in previous studies. Galaxy groups are identified by applying a modified version of the halo-based group-finding algorithm of Yang et al. (2005) to the galaxy sample of SDSS Data Release 7 with 0.01 < z < 0.2, with redshift completeness greater than 70 per cent. Following common practice, for each group we take the most massive galaxy as the central galaxy and consider the other member galaxies as satellites. Stellar masses are taken from the NSA catalogue (Blanton et al. 2011).
We also use the three-dimensional reconstructed local mass density (δ = ρ/ρ, where ρ is the local matter density andρ is the mean matter density) to quantify the density of the local environment for each galaxy in our sample. The mass density field is taken from the Exploring the Local Universe with the reConstructed Initial Density Field project (ELUCID; Wang et al. 2016), that by construction can very well reproduce the observed distributions of both galaxies and groups of galaxies. ELUCID provides local densities measured at different scales. In this work, we use the density estimated by smoothing the density field with a Gaussian kernel at a scale of 1 Mpc h −1 .

Identification of 'turnover' galaxies
We firstly identify "turnover" galaxies in our sample that present a significant star formation enhancement in their centres, applying a selection procedure similar to that in our previous work (Lin et al. 2017, hereafter Paper I). In short, the identification procedure comprises two steps. In the first step, we select potential turnover galaxies that each show a central upturn or drop in the onedimensional radial profiles of three parameters EW(Hα), EW(Hδ A ), and D n (4000). Next, we visually examine the SDSS gri-image, twodimensional maps, and the radial profiles of each galaxy, and decide whether or not the turnover feature indeed occurs in the inner regions. Here, 'inner regions' mean the inner parts of the disc, including both the bulge and bar-dominated regions (if present). In the rest of this subsection, we describe the identification process in more detail.
For a given radial profile, we start from the outermost region by fitting a linear function to the three data points at the largest radii. We then determine whether the next data point at the neighbouring smaller radius follows the same trend, by evaluating the deviation of that data point from the best-fitting line. If this deviation is less than the 1σ scatter of all the spaxels at the same radius as the data point, we perform the linear fit again but now include the new data point. In this way, the data points at smaller radii are gradually added to the fit, one by one. If there is a newly added data point deviated significantly from the best-fitting line by more than 1σ , indicating a sudden break in the slope of the profile. From this radius, we repeat the above process, start fitting a new line with a different slope to the next three data points and gradually adding more data points at smaller radii. We thus consider the whole profile from the outside in by adding data points one-by-one with a slightly adjusted slope or an entirely new slope until the radius reaches half of the PSF FWHM (1.25 arcsec).
This procedure was applied to the profiles of EW(Hα), EW(Hδ A ), and D n (4000) for every galaxy in our sample. We then select turnover galaxy candidates by requiring at least one of the three parameters to have a slope break in its profile [increment for EW(Hα) and EW(Hδ A ), and turn-down for D n (4000) as the radius decreases] within R 50 or R bar if R bar >R 50 . This restriction to the inner regions is based on our empirical finding from Paper I that nearly all turnover radii are smaller than R 50 or R bar . We thereby select 223 turnover galaxy candidates out of the 1425 galaxies in our sample.
We then visually inspect the SDSS gri-image of each turnover galaxy candidate, overplotting the isophotal ellipse at the break radius on to the image. Much attention is paid to the location of the break relative to the bulge, bar and disc components. A candidate is selected as a real turnover galaxy if the innermost break happens within the bulge radius if the galaxy is unbarred, or within the bar length if it is a barred galaxy. This visual inspection was carried out independently by two of the authors (LL and CD). In this step, 102 candidates are excluded for various reasons, as illustrated in Fig. 2 showing for three example galaxies the optical image, the two-dimensional maps of EW(Hα), EW(Hδ A ), D n (4000), log( SFR ), and Baldwin, Phillips and Telervich diagrams (BPT; Baldwin, Phillips & Terlevich 1981) , and the corresponding radial profiles. The SFRs are calculated from dust-corrected Hα fluxes (see Section 3.3 for details). In the majority of these galaxies (72/102), the slope break of the radial profile reflects the transition between bulge and disc, rather than a central turnover feature within the bulge. This is the case for the first galaxy (MaNGA ID 7962-12704) in Fig. 2. This galaxy shows a break at ≈4 arscsec in EW(Hα), EW(Hδ A ), and D n (4000) at a radius that is much larger than the central bright core, which is likely a pseudo-bulge with the size of ∼1 arcsec but shows no obvious turnover features in the profiles. There are 24 candidates are even bulgeless, thus they are also excluded. The breaks are due to fluctuations among spiral arms. For example, the second galaxy (MaNGA ID 7962-12701) in Fig. 2 is a bulgeless barred galaxy. It has a break within R 50 , indicated by the vertical black dashed line, but the profiles within the bar (indicated by the green dotted line) are roughly flat, with no obvious turnover feature. Finally, a few galaxies (6/102) are excluded because this slope changes are too gentle and smooth, and are thus unlikely caused by significant changes of stellar populations. As shown in the bottom panels, the selected galaxy (MaNGA ID 8313-12704) does have a variation of the slope of its EW(Hα) profile, but the signal is no obvious.
Our final sample includes 121 galaxies selected following the above two-step scheme based on EW(Hα), EW(Hδ A ), and D n (4000) independently, and each presenting significant central turnover feature in at least one of the three parameters. However, we find that in most cases the turnover feature is simultaneously present in all the three parameters. The number of turnover galaxies identified from each parameter is listed in Table 1. More than 80 per cent are identified from at least two parameters. In addition, the turnover radii measured from these different parameters are highly consistent with each other. This can be seen from Fig. 3, comparing the turnover radii determined from the three different parameters. The typical scatter is 0.11 dex. In what follows, if one galaxy has different turnover radii measurements, we thus adopt the median of the different turnover radii as our turnover radius and the standard deviation as its uncertainty. If one galaxy has only one turnover radii measurement, we adopt the typical scatter of 0.11 dex as its uncertainty. The measurements of the turnover features are listed in Table 3 for every turnover galaxy. Fig. 4 shows four examples of the turnover galaxies, analogously to Fig. 2. The first two are barred galaxies with regular morphologies, while the last two have slightly disturbed morphologies and are included in the pair catalogue. The first galaxy, MaNGA 8602-12701, has obvious breaks at ≈5 arcsec, and its central region is likely to be low-ionization nuclear emission-line region (LINER) in the BPT diagram. The second galaxy, MaNGA 9487-12701, has similar profiles to the first one, but its central region is classified as star-forming according to the BPT diagram. The third galaxy, MaNGA 9507-12705 is an unbarred galaxy in a pair of Stage 3. It shows a break at 3.5 arcsec in both EW(Hα) and D n (4000). The break radius is smaller than the bulge radius ≈5 arcsec. The last galaxy, MaNGA 8484-12703, is a barred galaxy in a pair of Stage of 2. It has a break in all three profiles, with a star-forming centre.

Global properties of turnover galaxies
We start by examining the global properties of our turnover galaxies. For comparison, we also select a control sample of galaxies from the MaNGA/MPL-7 by requiring the control sample galaxies and turnover sample galaxy to have similar stellar masses and NUVr colours. For this purpose, we make use of the stellar masses and  (4000), and SFR surface density profile (measured from dust-corrected Hα emission), and resolved BPT diagram. The red hexagons correspond to the field of view of the MaNGA bundle. In each map, the blue and green ellipses indicate the potential turnover radii and (if present) the bar, respectively. In each radial profile, the solid and dotted black lines are the median value and 1σ scatter in each annulus, the segment fitting is plotted as solid red lines, the green (if present) and black vertical dashed lines indicate, respectively, the end of the bar (R bar ) and the effective radius (R 50 ). The first galaxy has a break radius in all the three indicators, but it happens at the bulge-disc transition (by visual inspection). We thus conclude it is not a genuine turnover within bulge. The second galaxy has break radius within R 50 , but the profiles are straight within R bar . The last galaxy is excluded because it lacks of a robust feature. More details are provided in the text.
NUVr colours from the NSA catalogue. Typical uncertainties of the stellar mass is about 0.1 dex, 3 NUVr uncertainties are 0.06 dex (Blanton & Roweis 2007). In practice, for each turnover galaxy, we search for a control galaxy from the MPL-7 with a stellar mass difference log M * /M ≤ 0.2 and a colour difference (NUVr) ≤ 0.3. We randomly select one of the many galaxies that meet these criteria, regardless of the presence or absence of a turnover feature. We note that the control sample includes 41 barred galaxies, i.e. it has a bar fraction of about one-third, consistent with the typical strong bar fraction of galaxies in similar mass and colour ranges (Barazza et al. 2008;Cheung et al. 2013). In addition, the control sample . Comparison between the different turnover radii of our final turnover sample galaxies. The black dotted line is the one-to-one line, the blue and red solid lines are the linear fittings between R t, D4000 and R t, EW(Hα) , R t, D4000 and R t, EW(Hδ) , respectively. Different measurements are highly consistent with each other, with a typical 1σ scatter of 0.11 dex.
includes 37 turnover galaxies, indicating a typical turnover galaxy fraction (30 ± 4 per cent) for galaxies with M * > 10 10 M and NUV r <5. Table 1 lists the number of barred and unbarred galaxies, the number of galaxies at different pair stages, and the number of galaxies in different classes in the BPT diagram, for both the turnover galaxy sample and the control sample. The most striking difference between the turnover and control samples is that between their barred galaxy fraction, which is 89 ± 3 (108/121) and 34 ± 4 per cent (41/121), respectively. The high fraction of bars in turnover galaxies is expected, given the similarly high fraction of bars in the turnover galaxies from CALIFA (88 ± 8 per cent, 15/17; see Paper I). This, again, strongly suggests a close relationship between the presence of a bar and a central star formation enhancement (as indicated by a turnover).
In Fig. 5, we examine the global properties of our turnover galaxies by plotting them on four diagrams: NUVr versus log M * /M , R 90 /R 50 versus log M * /M , log 1 /M kpc −2 versus log M * /M , and the BPT diagram. The stellar surface mass density within a radius of 1 kpc, 1 , is calculated using the total light within the radius of 1 kpc in i-band and stellar mass-to-light ratio M/L i from g − i colour following Fang et al. (2013). We adopt the typical uncertainties of colour-based stellar M/L ratio estimation are 0.1 dex (Bell et al. 2003), the typical uncertainties of R 50 and R 90 measurements are 4 per cent (Blanton et al. 2011). Emission-line measurements and the corresponding uncertainties for the BPT diagram are taken from the MPA-JHU database 4 (Brinchmann et al. 2004), based on the central 3-arcsec-fibre SDSS spectra. In Fig. 5, the barred and unbarred turnover galaxies are shown as blue and red circles, separately. The control galaxies are shown as green crosses. For comparison with the general galaxy population, we also select a volume-limited galaxy sample from the NSA, consisting of 35 070 galaxies with the r-band absolute magnitude M r < −17.2 and redshift 0.01 < z < 0.03. The distribution of this sample is illustrated as the greyscale background in the figure.
In the colour-mass diagram (top left-hand panel of Fig. 5), galaxies can be divided into three regions: blue cloud, green valley, and red sequence, as demonstrated by the volume-limited sample in the background. The turnover galaxies are found mainly in the blue cloud and green valley (NUVr 5), with stellar masses M * ranging from ∼10 10 to ∼10 11.5 M . After visually inspecting their SDSS images, we find most of the turnover galaxies to have morphological types from Sab to Sbc, with very few later than Sc. This result suggests that the turnover features indicative of enhanced star formation in galactic centres occur preferentially in massive spiral galaxies. The ranges of stellar masses and colours found here for the turnover galaxies are consistent with those found from the CALIFA sample in Paper I. By construction, the control sample covers similar ranges of stellar masses and NUVr colours as the turnover galaxy sample.
The upper right-hand panel of Fig. 5 shows the distribution of the different samples in the plane of concentration versus stellar mass, where concentration is quantified by the ratio of R 90 to R 50 in the r-band (R 90 /R 50 ). After being matched with similar stellar mass, NUVr colour and morphological type, the turnover sample is found to have similar concentrations with control sample. Majority of them have intermediate-to-low concentrations (R 90 /R 50 3). We further examine the bulge-to-total ratios (B/T) in turnover and control galaxies using the B/T measurements from the bulge-disc Sérsic+exponential profiles decomposition catalogue of Fischer, Domínguez Sánchez & Bernardi (2019). Similar distributions are seen in turnover and control galaxies, with a median B/T of 0.2. We should point out that a three-component decomposition including bulge, bar, and disc components would be preferable, given that the majority of turnover galaxies are barred. If bars are neglected in the decomposition, this usually results in an overestimation of the B/T ratio (Gadotti 2008). Thus, one may expect a smaller B/T in turnover galaxies.
The central stellar mass surface density within the radius of central 1kpc, 1 , has been suggested to be a good indicator of the prominence of a central bulge (e.g. Fang et al. 2013;Luo et al. 2020). In the lower left-hand panel of Fig. 5, the black dashed line indicates the dividing line between quenched and unquenched galaxies proposed by Fang et al. (2013): at fixed mass, it has been found that most of the red (quenched) galaxies lie above this line. As can be seen in Fig. 5, a few (11 ± 3 per cent, 14/121) turnover galaxies are above this line, consistent with their relatively blue NUVr colours (NUVr 5) and relatively small concentration parameters (R 90 /R 50 3). In contrast, more (17 ± 3 per cent, 21/121) control galaxies are above the dividing line, suggesting the presence of a high-density central bulge. In addition, we have compared the central r-band surface brightnesses of the turnover and control sample galaxies, finding the former to be brighter on average than the latter. This implies a lower average mass-to-light ratio (M/L), thus younger stellar populations, in the centres of turnover galaxies compared to the centres of control galaxies.
The lower right-hand panel of Fig. 5 shows the sample galaxies in the BPT diagram. The turnover and control sample galaxies have very similar distributions in this diagram, with galaxies distributed over all the different ionization classes. However, when compared to the general galaxy population, the turnover galaxies (and the corresponding control galaxies) are limited to particular regions. In the star-forming region, the turnover and control sample galaxies are located in the region of low [O III]/Hβ and high [N II]/Hα. This can be understood from the relatively high stellar masses of the turnover galaxies. In the composite and AGN regions, the turnover, and control Notes. Columns are as follows: Column (1) -MaNGA ID; Column (2) -Turnover radius and uncertainty; Column (3) -log EW(Hα) enhancement and uncertainty; Column (4) -EW(Hδ) enhancement and uncertainty; Column (5) -D n (4000) drop and uncertainty; Column (6) -log SFR enhancement within turnover radius and uncertainty. Full table is available online in the machinereadable format.
galaxies are similarly found in the Seyfert, LINER, and composite regions. The number of galaxies of each type is listed in Table 1. We find that 25 ± 4 per cent (30/121) of turnover galaxies are classified as star-forming galaxies, 16 ± 3 per cent (19/121) are composite, 32 ± 4 per cent (39/121) are central LINERs and 21 ± 4 per cent (25/121) are Seyferts. Turnover galaxies are likely to have slightly higher fractions for LINERs and Seyferts. We will return to these populations in Section 4.6.
To summarize, central turnover features are found in massive spiral galaxies with M * 10 10 M and the majority of them are barred, falling in the blue cloud, and green valley with NUVr 5. Turnover galaxies have intermediate-to-low concentration parameters (R 90 /R 50 3) and relatively low central surface mass densities, suggesting relatively low B/T. Their centres harbour all the different ionization mechanisms in the BPT diagram, with slightly higher fractions of LINERs and Seyferts to those of the control sample.

Enhanced star formation in turnover galaxies
The central turnovers suggest young stellar populations are enhanced in the centres of turnover galaxies, due to recent and/or ongoing star formation. To quantify these enhancements, we follow Paper I and obtain an extrapolated value of EW(Hα), EW(Hδ A ), and D n (4000) for the central region of each galaxy. This is derived by extrapolating the linear fit to the radial profile from radii beyond the turnover radius down to the centre. The extrapolated parameters can be considered as the values expected if the central star formation were not enhanced. In the top panels of Fig. 6, we plot both the observed (purple solid circles) and extrapolated (green open circles) EW(Hα) as a function of stellar mass (lefthand panel) and NUVr (middle panel). The similar but larger symbols indicate the EW(Hα) running averages for three intervals of stellar mass or colour. The observed and extrapolated values at a given mass or colour are connected by a black arrow. As can be seen in the figure, the observed Hα equivalent widths are enhanced by about 1 dex compared to the extrapolated ones, and this appears to be independent of stellar mass. The reddest colour bin is less enhanced, which may be reflecting the dependence of the central star formation enhancement on cold gas content at the centre. However, given the relatively small number of galaxies in the reddest colour bin, this slightly weak enhancement should not be overemphasized. Fig. 6(c) shows the observed and extrapolated EW(Hα) as a function of d SF , the effective distance from the star-forming sequence in the BPT diagram. Here, d SF is calculated by using a set of curves parallel to but offset from the diagnostic line of Kauffmann et al. (2003) as a metric, to characterize the position of galaxies along the mixing sequence. As can be seen, the centres of turnover galaxies can be broadly divided into two groups, i.e. those star-forming centres at low d SF and those with AGN-like centres at high d SF (with just a few galaxies in between, with d SF ≈ 0.5). The Hα equivalent widths are enhanced by about one dex for star-forming galaxies, more than for galaxies with an AGN-like centre, that on average are enhanced by ≈0.5 dex.
Hα emission is a proxy for the SFR for star-forming regions. The data in the top panels of Fig. 6 thus provide clear evidence for an enhancement of central star formation in turnover galaxies. To quantify the central SF enhancements more directly, we use dust-corrected Hα flux maps to infer SFR surface densities ( SFR ), and then integrate the profile within the turnover region. Following common practice, the dust extinctions are calculated from the Hα/Hβ ratios by assuming an intrinsic ratio of 2.86 and adopting a Milky-Way-like extinction law from Cardelli, Clayton & Mathis (1989). The SFR calibration is given by Calzetti (2013) (2001) stellar initial mass function (IMF) and a constant star formation time-scale from 0.1 to 100 M . After deriving the SFR surface density profiles, we apply the same analysis as done previously for EW(Hα), EW(Hδ A ), and D n (4000), fitting the SFR profiles and deriving extrapolated SFR for the central regions. Then, we integrate the profiles from the centres to the turnover radii to get the observed and extrapolated SFRs, respectively, denoted as SFR <Rt .
The bottom panels of Fig. 6 show the observed and extrapolated SFRs of the central regions of the turnover galaxies as a function of log M * /M , NUVr and d SF . The symbols/colours are the same as top panels, except that galaxies with LINER or Seyfert centres are shown in grey. We caution that the Hα emission in these galaxies could be contaminated by nuclear activity. The SFRs running averages indicated by the large symbols are calculated using SF galaxies only. Overall, turnover galaxies have SFR enhancements of about 0.5-1 dex, similar to the trends seen in EW(Hα). . Four examples of turnover galaxies. From the left-to right-hand panel: Optical gri image, maps, and profiles of EW(Hα), EW(Hδ A ), D n (4000), and SFR surface density (measured from dust-corrected Hα emission), and resolved BPT diagram. Labels and symbols are the same as Fig. 2. The first two are barred galaxies, which are selected to demonstrate central SF and LINER-emission. The last two are disturbed galaxies in pair systems. The central SFRs could be contaminated by nuclear activity if they have LINER or Syfert centres. (red open circles) galaxies. The control sample galaxies (green crosses) are matched in stellar mass and NUV-r colour. The volume-limited SDSS sample (greyscale) is shown as background. In the top left-hand panel, the black dashed lines are for NUV−r = 4 and 5, commonly used to divide galaxies into blue-cloud, green-valley, and red-sequence galaxies. In the bottom left-hand panel, the black dashed line indicates the 1 threshold for quenching (Fang et al. 2013). In the bottom right-hand panel, the three black lines separate galaxies into the star-forming, composite, LINER, and Seyfert categories (Kauffmann et al. 2003;Kewley et al. 2006;Cid Fernandes et al. 2010). Typical uncertainties are presented at the lower right-hand corner in each panel.

Relationship between turnovers and bars
Now we examine the relationship between the turnover features (including the turnover radius R t and the SFR enhancement within the turnover radius [log (SFR <Rt /M yr −1 )]) and bar properties (including the bar length R bar and bar ellipticity e bar ). The results are shown in Fig. 7. In each panel, the turnover barred galaxies of star-forming and composite types and those of LINER and Seyfert types are shown as solid blue and empty blue circles, respectively. In some panels, we normalize the bar length and turnover radius by R 25 , the radius at which the r-band surface brightness is 25 mag arcsec −2 .
The top left-hand panel of Fig. 7 presents the tightest correlation -an increasing of the turnover radius with increasing bar length, a trend holding for both SF/composite galaxies and AGN. This correlation can be well described by a linear function, R t = (0.35 ± 0.02)R bar +(0.61 ± 0.13), plotted as a dashed line in the panel. A similar correlation holds, albeit with more scatter, when R t and R bar are scaled by R 25 (Fig. 7d). The star formation enhancement within the turnover region, as quantified by [log (SFR <Rt /M yr −1 )] is also positively correlated with R bar and R bar /R 25 (Fig. 7g and h), but with with larger scatter than the correlations between turnover radius and bar length.
Unlike bar length, bar ellipticity shows no correlation with turnover radius or SFR enhancement. Both bar length and bar ellipticity are commonly used as proxies of bar strength, because of their correlations with gravitational bar torques (e.g. Combes & Sanders 1981;Block et al. 2001). Bar length and ellipticity are indeed correlated, as shown by Wang et al. (2012) with a large sample of barred galaxies from SDSS. However, the turnover galaxies in our sample show no correlation between R bar and e bar , as can be seen in Fig. 7(c). Therefore, the lack of correlation between the turnover features and bar ellipticity is likely due to the limited size of our sample, or it may be suggesting that bar length is really more fundamental to the turnover feature. To clarity this, one would need a larger sample of turnover galaxies as well as theoretical studies (such as numerically simulating the co-evolution of galactic centres and bars). Besides, the measurement of bar ellipticity using the ellipse fitting technique would be affected by B/T, bar strength and other features like spiral arms or rings. We matched our sample with the catalogue in Kruk et al. (2018) in which they considered bar in photometric decomposition. In total, 44 galaxies are included in their measurements. We still found no obvious correlations between turnover radius versus ellipticity or B/T in their catalogue. Hopefully this can be examined again with larger samples in the future.

Central versus global star formation enhancements
Central turnovers were identified by comparing the observed and inward extrapolated values of the three parameters closely related to recent star formation that we considered. These turnovers reflect relative SF enhancements with respect to the preexisting underlying old populations at the galaxy centres. Therefore, the central turnover galaxies selected in this way are expected to have recently experienced bulge growth and/or rejuvenation (see also Thomas & Davies 2006;Coelho & Gadotti 2011;Méndez-Abreu et al. 2014;Chown et al. 2019) In the literature, there have been various attempts to estimate SFR enhancements in the central regions of galaxies, driven by instabilities induced by bars or tidal interactions. For instance, Li et al. (2008) quantified the enhancements of SFR as a function of projected distances to neighbouring galaxies, finding the SFRs can be enhanced by a factor of 2-3 at scales comparable to the sizes of individual galaxies. Ellison et al. (2011) estimated the average central log SFR of barred galaxies relative to unbarred galaxies at fixed stellar mass, finding the former to be 0.2-0.3 dex higher. Wang et al. (2012) calculated the central-to-total SFR ratios of a sample of barred galaxies and a carefully selected sample of control galaxies, finding the average SFR ratio of galaxies with strong bars is about 0.2-0.3 dex higher than that of the control galaxies (see table 1 of their paper).
The typical SF enhancement of our turnover galaxies is ≈1.0 dex (see Fig. 6), higher than those previously found (0.2-0.4 dex as described above). We argue that this difference can be mainly attributed to the different definitions of star formation enhancement. In our case, the enhancement is quantified for the central region of the each galaxy, comparing the observed and expected (extrapolated) SFRs. In fact, in a recent paper, Chown et al. (2019) estimated the enhancements of the central SFRs of barred galaxies in the same way as in our work, also finding a high enhancement of ≈1 dex (see their Fig. 5). In previous studies, the enhancement of each galaxy was calculated by comparing the central SFR with either the total SFR of the whole galaxy (e.g. Wang et al. 2012) or the central SFR of control galaxies Figure 7. Turnover quantities as a function of bar properties. Bar properties are described by the bar length (R bar ), normalized bar length (R bar /R 25 ), and bar ellipticity (e bar ) along the abscissa. Turnover features are characterized by the turnover radius R t , the normalized turnover radius R t /R 25 and integral SFR enhancement within R t ( [log SFR <Rt /M yr −1 ]) along the ordinate. Only barred turnover galaxies are shown. Data points are separated into SF (solid blue circles) and AGN (open blue circles) according to their locations in the BPT diagram. We caution that the SFR estimate in LINER/Syfert galaxies could be contaminated by nuclear activity. The dashed line in the top left-hand panel shows a linear fit between turnover radius and bar length. The Spearman correlation coefficients r are listed at the top left-hand corner. Typical uncertainties are presented at the lower right-hand corner.
(e.g. Li et al. 2008;Ellison et al. 2011). In fact, we re-calculated log SFR by comparing the central SFR of each turnover galaxy with the central SFR of its control galaxy, and found an average SFR enhancement of 0.28 dex, thus consistent with previous studies.
Nevertheless, we would like to emphasize that the turnover galaxies studied here are not necessarily the same objects as those studied in previous studies, that usually present globally higher SFRs, with peak SFRs in the centre and weaker SFRs at larger radii. By investigating the spatially resolved SFR profiles, galaxies above (or below) the star formation main sequence (SFMS) are found to have enhanced (or suppressed) star formation at all galactocentric radii (e.g. Ellison et al. 2018;Guo et al. 2019), though with stronger effects at the centre. This can be interpreted as variations of the overall gas inflow rates (e.g. Lilly et al. 2013;Wang et al. 2019). However, our turnover galaxies are expected to have star formation enhanced at their centres compared to the outskirts of the bulge or bar regions. According to our sample selection criteria, turnover galaxies do not necessarily have globally higher SFRs compared to the SFMS. Therefore, the turnover features are related to radial gas inflows to the galactic centres, rather than enhancements of the overall gas inflow rates.
Besides, instead of considering central star formation enhancement, one could think about turnover behaviour in another way: The central region is surrounded by a "desert" ring, where the star formation has been suppressed (James & Percival 2018). According to the bar evolution scenario proposed in James & Percival (2018), the star formation enhancement would be an initial response after the radial gas inflow, after that, the bar would suppress SF and form a desert ring. In this study, we define turnover galaxies according to their radial profiles. Thus our sample could include SF enhancement at the beginning phase, as well as objects at the later phase that the surrounding area is being suppressed, while the central SF still exist. The presence of a star formation desert potentially increases the contrast between central and surrounding region in the radial profiles, which could bias the ability to detect turnover features towards strong-barred galaxies.
4.2 Lack of "turnover" galaxies below M * ∼ 10 10 M It is interesting to note that the turnover features identified here are mostly found in relatively high-mass galaxies, with stellar masses above 10 10 M . Although adopting a different definition of star formation enhancement, the earlier study of barred galaxies in SDSS by Ellison et al. (2011) also found a lack of star formation enhancement at stellar masses below 10 10 M . The lack of central turnover features in low-mass galaxies indicates a lack of central star formation enhancements, that probably is not unexpected for low-mass galaxies. First, most low-mass spirals do not have a clear bulge structure or are even bulgeless. Furthermore, in most cases, the low-mass galaxies in our sample are actively forming stars (and are thus entirely dominated by young stellar populations) according to their EW(Hα), EW(Hδ A ), and/or D n (4000) profiles, that are typically flat or slightly negative, and could thus be well described by a single slope. Therefore, low-mass galaxies may be globally enhanced, i.e. with star formation enhancement at all radii (e.g. Fraser-McKelvie et al. 2019). In this work, we only considered the central regions of our galaxies. In a following paper, we will examine the characteristic radial profile of star formation enhancement (or suppression) as a function of galaxy mass and other properties.

Role of bars
We find 89 ± 3 per cent of the turnover galaxies are barred, a very high fraction compared to the bar fraction of 34 ± 4 per cent in the control sample. A similarly high fraction of 88 ± 8 per cent was found in our previous work (Lin et al. 2017) from the CALIFA survey. Combined with the finding of Chown et al. (2019) on centrally concentrated molecular gas in turnover galaxies, our results thus reinforce the conclusion that bars must play a dominated role in enhancing central star formation in (barred) galaxies, by driving cold gas inflows.
In this work, we have further examined the correlation of turnover features with bar properties, as quantified by bar length and ellipticity. It is striking that turnover radius is tightly correlated with bar length, the turnover occurring at about one-third of the bar length. This is, again, strongly suggesting that bars are indeed the main driving factor of turnover features. In a few previous studies, similar turnover positions have been found, e.g. Seidel et al. (2015Seidel et al. ( , 2016 identified breaks in the radial profiles of kinematic features and stellar population indices of barred galaxies. Those authors found the breaks usually happen at 0.13-0.2R bar . The turnover radii in our galaxies are slightly larger, likely due to the fact that the turnover features in our work are traced by young stellar populations. As the central bulges grow, one would expect to find younger populations at larger radii, consistent with the inside-out formation of bars as predicted by numerical simulations (e.g. Athanassoula et al. 2016). A recent study that looks at the stellar populations in the barred galaxies also find a similar break at ≈0.3R bar (Neumann et al. 2020).
We can also compare our bar classifications with the barred galaxies of the Galaxy Zoo (GZ) catalogue (Willett et al. 2013). In the MPL-7 sample, ≈90 per cent galaxies can be matched with the current GZ catalogue. Using the criteria bar probability p bar > 0.5 and non-edge-on probability p edge < 0.5, 269 objects are selected as barred galaxies. The bar fraction (19 ± 1 per cent) is consistent with our own (18 ± 1 per cent). Comparing the two barred samples, 192 objects (≈70 per cent) overlap. The differences are mainly due to weak bars, whose identification is somewhat ambiguous and subjective, both GZ and our classification miss or mis-identify some weak bars, and we identified more lenticular objects as barred galaxies. Nevertheless, regarding our turnover sample, if we use the GZ bar classifications, there are 89 barred turnover galaxies, 26 unbarred turnover galaxies, and 6 galaxies lacking a classification. These numbers yield a similarly high bar fraction of 77 ± 4 per cent among turnover galaxies. So our results would not change much if we used GZ bar classifications.

Role of interactions/mergers and local density
In addition to bar-driven instabilities, tidal interactions with close companions are also known to be able to drive cold gas from the outer parts of discs to galactic centres, thus leading to similar star formation enhancements in the central regions of galaxies (e.g. Ellison et al. 2008;Li et al. 2008). For instance, a recent MaNGA-based study by Pan et al. (2019) revealed significant star formation enhancements in paired galaxies with a separation less than 20 kpc, an effect that is strongest at the galaxy centres. We have therefore also examined the correlation of our turnover galaxies with galaxy pairs/mergers, but find no obvious link. Of the 121 turnover galaxies, only 20 are in a pair/merger according to the classification of Pan et al. (2019). This fraction (17 ± 3 per cent) is consistent with the fraction of pair/merger galaxies in the control sample (24/121). On the other hand, of the whole 228 pairs/mergers galaxies in our parent sample, there are only 20 galaxies are identified as turnover galaxies, the fraction (9 ± 2 per cent) is even lower.
In our turnover galaxies, except the 108 barred galaxies, there are 16 turnover galaxies with no bar. One may expect the 16 unbarred turnover galaxies to be dominated by galaxy pairs/mergers, so that bars and tidal interactions combined can completely explain the central star formation enhancements of our sample galaxies. Indeed, in a recent study, Chown et al. (2019) examined the central star formation enhancement in a sample of 58 nearby galaxies from CALIFA, finding the enhanced star formation are either barred (13/19, 68 ± 11 per cent) or in pairs/mergers (6/19, 32 ± 11 per cent). However, this is not the case in our sample, where we find a similarly high bar fraction in the turnover galaxies (89 ± 3 per cent), but only 4 of the 13 unbarred turnover galaxies are associated with pairs/mergers (31 ± 13 per cent in our sample versus 100 per cent in theirs). Most (16/21) of our turnover galaxies with galaxy pairs/mergers are also barred.
In conclusion, in contrast to our initial expectation, galaxy interactions/mergers are not obviously related to the central star formation enhancement as indicated by the central turnover features. This is true for the majority of the turnover galaxies, regardless of whether they are barred or unbarred. This result may be caused by multiple Figure 8. Cumulative fraction of local density in 1 Mpc for turnover galaxies (black dashed line) and control sample galaxies (red solid line). The shaded areas represent the uncertainties computed by assuming the Poisson distribution. The K-S test probability is 0.15, suggesting no significant difference between the two samples.
reasons. First, several paired galaxies are dry mergers, thus lacking cold gas to support the central star formation. We found that about 13 per cent (30/228) of the paired galaxies have NUVr > 5, with little cold gas to drive and/or support central star formation. Second, as pointed out in Section 4.2, low-mass galaxies usually have weak or no bulge component, showing no turnover features in their centres, although their star formation may be globally enhanced. The rest paired galaxies (198/228) with NUVr < 5, more than half (54 ± 3 per cent, 124/228) have stellar masses less than 10 10 M . For these galaxies, interactions/mergers probably have a more important role in enhancing star formation over large areas of the galaxies, rather than only driving central star formation enhancements. We indeed see several paired galaxies in which enhanced star formation extends to large radii, far beyond the central regions. In our work, they are not selected as turnover galaxies as the turnover radii are defined to be smaller than the size of the bulge or bar. Third, recent studies suggest that star formation enhancements are significant only when the paired galaxies are close enough. For instance, Pan et al. (2019) found significant star formation enhancement when pair separation is less than 20 kpc. Finally, other parameters like the mass ratios of the paired galaxies, their merging geometries and the gas fractions should also affect the resulting star formation enhancements.
We have also examined the environments of our galaxies over larger physical scales. First, using the group catalogue of Yang et al. (2007), we find that the ratio of central and satellite galaxies in the turnover galaxy sample is very similar to that of control sample, with a central galaxy fraction of ≈70 per cent in both cases. Next, based on the density field of the local Universe reconstructed by Wang et al. (2016), Fig. 8 shows the cumulative distribution of the local mass density (δ = ρ/ρ) measured over a scale of 1 Mpc for both the turnover galaxies and the control sample galaxies. The two samples show very similar cumulative density distributions, with a Kolmogorov-Smirnov (K-S) test probability larger than 10 per cent indicating no significant difference between the two distributions. Therefore, we conclude that the large-scale environment is not obviously linked to the turnover features of our galaxies.

Role of cold gas
Central star formation enhancements are expected to be associated with higher cold gas contents in galactic centres. In fact, Chown et al. (2019) jointly analysed the integral-field spectroscopy from CALIFA and CO line imaging from the Extragalatic Database for Galaxy Evolution (EDGE) survey for a sample of nearby galaxies. They found the central enhancements of SFR to be positively correlated with the concentrations of molecular gas, an effect only present in barred galaxies and those in pairs/mergers. This provides direct evidence that cold gas inflows due to baror merger-driven instabilities enhance star formation in the central regions of galaxies. On the other hand, Ellison et al. (2020) studied the CO content of MaNGA galaxies with a central starburst, finding that the elevated SFRs in their centres are instead driven by enhancements of the star formation efficiency, defined as the SFR divided by the H 2 mass. Nevertheless, they also found lower gas-phase metallicities in the central regions of starburst galaxies compared to their outer parts, consistent with a metal-poor gas inflow scenario.
In our previous work (Paper I), a weak correlation was also found between the turnover features and the H I gas mass fraction of their host galaxies. Here we attempted to examine the H I gas mass of our MaNGA galaxies using both the ALFALFA survey (Haynes et al. 2018) and the H I follow-up observations being carried out by the MaNGA team , which is aiming to observe all MaNGA galaxies with z < 0.05. We thereby obtained H I gas masses for 90 turnover galaxies, but there is no significant correlation of neither turnover radius or [log (SFR <Rt /M yr −1 )] with the H I masses. We also compared the H I fraction in turnover barred and non-turnover barred galaxies, and we find no obvious difference either. It is thus very likely that the weak correlation found in our previous work is not real, but was due to the small sample size. Having said that, the sample from MaNGA is larger but still rather limited. Further analysis of the H I gas content of our galaxies is beyond the scope of this paper, and we leave it to future works when larger samples of H I observations are available.
Nevertheless, even with a more complete sample, the absence of a correlation between central star formation enhancement and H I gas mass may not be surprising, for the following two considerations. First, it is widely believed that atomic gas is less closely related to star formation than molecular gas (Bigiel et al. 2008;Schruba et al. 2011). Secondly, the cold gas is expected to be driven towards the centres of the galaxies, but the global gas contents are not necessarily increased. This is indeed what was found by Chown et al. (2019) considering the spatially resolved CO maps. Unfortunately, spatially resolved H I maps are not yet available for our galaxies. In a sample with spatially resolved H I, selected to have both strong bars and significant H I content, Newnham et al. (2020) demonstrated that central H I holes are common, but not ubiquitous in barred galaxies. In that sample, the presence or absence of a H I hole correlated with total H I gas fraction, as well as the global star formation properties, with barred galaxies without central H I holes more likely to have the highest gas fractions and ongoing global SF.

Turnover galaxies with central LINER/Seyfert emission
It is notable that nearly half of our turnover galaxies are located in the AGN region of the BPT diagram, including 32 ± 4 per cent (39/121) LINERs and 21 ± 4 per cent (25/121) Seyferts. This appears to be in conflict with our expectation that the central regions of these galaxies are dominated by recent/ongoing star formation, as indicated by the turnover of EW(Hα), EW(Hδ A ), and D n (4000).
The origin of LINER emission has been widely discussed in recent years. Rather than connecting it with central low ionization nuclear activity as traditionally assumed, recent studies suggest LINER-like emission may be produced by evolved stars, in particular post-AGB stars (Stasińska et al. 2008;Cid Fernandes et al. 2011;Belfiore et al. 2016). This again appears to be in conflict with our expectations. Is it possible that the Hα emission produced by star formation is mixed with diffuse LINER emission? To answer this question, we have examined the locations of our turnover regions on the spatially resolved SFMS (Rosales-Ortega et al. 2012;Sánchez et al. 2013). As suggested by Hsieh et al. (2017), the location on the spatially resolved main-sequence diagram may indicate whether the Hα emission is powered by young or old stellar populations. We find that most of the turnover regions lie between the star-forming sequence and the quiescent sequence. This supports our conjecture that the turnover regions are a mixture of underlying old stellar populations and recently formed young populations.
In addition, Belfiore et al. (2016) separated LINER-like galaxies into central-LIER (cLIER) and extended-LIER (eLIER) according to the location of the LINER emission within each host galaxy. Our LINER-like turnover galaxies are likely to be cLIERs based on their morphologies. The central D n (4000) of these cLIERs range from 1.6 to 1.9, systematically lower than the D n (4000) of eLIERs (see fig. 12 in Belfiore et al. (2016)), consistent with our suggestion that their central regions contain significant fractions of young stellar populations.
For the turnover galaxies with Seyfert-type emission, the Hα enhancement could actually be contributed by AGN activity. Communicating privately with colleagues working on AGN sources, we note that there are no any broad emission line in these galaxies. If this is true, these Seyfert-type turnover galaxies should be lowluminosity AGN. We find that most of them have turnovers in both EW(Hα) and D n (4000), while only five of them have only one EW(Hα) turnover. Since AGN can contribute to Hα emission but not to D n (4000), the associated D n (4000) turnovers indicate that these could be systems with AGN surrounded by circumnuclear star formation. The Seyfert-type fraction in the turnover sample (21 ± 4 per cent, 25/121) is slightly higher than that in the control sample (13 ± 3 per cent, 16/121), implying that the physical processes (e.g. bar-driven instability) driving central star formation also tend to boost nuclear activity. Given that the turnover galaxies are predominately barred, the association of the turnover feature with Seyfert-type emission also suggests a physical link between the presence of bars and nuclear activity.
Although many studies agree that bars indeed cause an inflow of gas toward the central regions of galaxies (e.g. Regan & Teuben 2004;Sheth et al. 2005;Hunt et al. 2008) as summarized in the Introduction, the connection between bars and nuclear activity is still under debate. A few high-resolution observations see fueling of nuclei through bars (e.g. Fathi et al. 2013), while many statistical studies do not find a different incidence of AGN activity associated with a bar (e.g. Ho, Filippenko & Sargent 1997;Hao et al. 2009). Numerical simulations suggest that additional effects, such as a secondary bar or a nuclear ring, are required to drive the gas further inwards (e.g. Regan & Teuben 2004). Recent high-resolution observations also suggest that there are many different ways to fuel galactic nuclei (e.g. Hunt et al. 2008;Gadotti et al. 2019). To figure out the connection between central turnover features and nuclear activity, high-resolution observations of gas and kinematics will be needed.

Future perspectives
In this work, we focus on the turnover feature happened within the inner region where a bulge is usually present. However, the star formation enhancement can be more extended and beyond the entire bulge. For example, the first galaxy  in Fig. 2 indeed shows a break in EW(Hα), EW(Hδ A ), and D n (4000) profiles. The break happens at ∼4 arcsec, while the central bright core (likely a pseudobulge) is much smaller, ∼1 arcsec, according to the SDSS image. Therefore, the "turnover" seen in the profiles cannot be a feature of the central core. Rather, the break happens at the transition between the inner region and the spirals. The whole inner region is young and without break in the profile. Besides, we also see some interacting galaxies that have extreme star formation across the entire galaxy. There is no turnover in the profile either because the entire galaxy is young and the profile can be fitted by single slope.
These phenomena are also interesting by themselves. They probably result from similar physics as turnover galaxies, such as non-axisymmetric-induced gas inflow through spiral arms or tidal interactions with close companions. Different selection criteria are needed to search for these galaxies, which can be an additional topic for future studies.
Many nuclear structures (such as nuclear ring, secondary bar and nuclear spiral) are likely to be built by the infalling gas and intense star formation. They help build up the central bulge and also play an important role in the secular evolution. A number of studies had been done from both theoretical and observational (Shlosman, Frank & Begelman 1989;Knapen et al. 1995;Piner, Stone & Teuben 1995;Comerón et al. 2010;Cole et al. 2014;Seo et al. 2019) to understand their physical properties. High spatial resolution observations from MUSE or ALMA are essential for fully understanding the formation and evolution of these substructures.

C O N C L U S I O N S
In this work, we analysed spatially resolved maps and radial profiles of EW(Hα), EW(Hδ A ), and D n (4000) for a sample of ≈1400 nearly face-on galaxies at z < 0.05, using integral field spectroscopy from the current sample (MPL-7) of the MaNGA survey. We identified 121 galaxies that presenting an upturns of EW(Hα), EW(Hδ A ), and/or a drop of D n (4000) with decreasing radius in their central regions. Such "turnover" features indicate recent/ongoing star formation in the central regions of these galaxies. We also examined global properties of the turnover galaxies, such as stellar mass and NUVr colour, and structural properties, particularly the presence of a bar as well as environment. We quantified the central turnover features and studied their correlation with global galaxy properties and bar properties.
Our conclusions can be summarized as follows: (i) Central turnover features are found primarily in galaxies with stellar masses above ∼10 10 M and NUVr colour less than five.
(ii) The majority of the turnover galaxies are barred, with a bar fraction of 89 ± 3 per cent (108/121), much higher than the bar fraction of 34 ± 4 per cent of a control sample of galaxies that are closely matched in stellar mass and colour. This reinforces previous findings that bar-driven instabilities lead to star formation enhancements in the central regions of galaxies.
(iii) For turnover galaxies with a bar, the radius of the central turnover region is found to positively correlate with the radius of the bar: R t = (0.35 ± 0.02)R bar +(0.61 ± 0.13). This provides further evidence supporting a close link of central star formation enhancement and bars.
(iv) The central star formation enhancements can be quantified by comparing the observed and inward-extrapolated SFRs within turnover radii, and can be as high as ≈1 dex in our turnover galaxies. These SFR enhancements show no obvious correlation with global properties such as stellar mass, colour, or location in the BPT diagram.
(v) There is no significant correlation of the turnover features with large-scale environment or the galaxy-galaxy interactions/mergers.

S U P P O RT I N G I N F O R M AT I O N
Supplementary data are available at MNRAS online. Table 2. Measurements of barred galaxies. Table 3. Turnover property measurements.
Please note: Oxford University Press is not responsible for the content or functionality of any supporting materials supplied by the authors. Any queries (other than missing material) should be directed to the corresponding author for the article.

A P P E N D I X A : T H E I M PAC T O F S PAT I A L R E S O L U T I O N O N M E A S U R E M E N T O F T U R N OV E R R A D I U S
One may worry about the limited spatial resolution of the MaNGA data which may not be high enough for all the turnover radii to be measured reliably. In Fig. A1 (upper panel), we show the turnover radius R t in unit of kpc as a function of redshift for the turnover galaxies in our sample (blue circles). The black dotted line corresponds to 2.5 arcsec, i.e. the spatial resolution of MaNGA. As can be seen, all the turnover galaxies are well above the dotted line, thus ensuring that the turnover features are well resolved in our galaxies. On the other hand, however, we see that the turnover radius is clearly correlated with the spatial resolution, or equivalent the redshift. This is not a real correlation, but a result of the correlation Figure A1. Turnover radius R t (upper panel) and the ratio of R t /R e (lower panel) as a function of redshift. Blue circles are turnover galaxies from this work based on MaNGA, and red triangles are turnover galaxies from our previous work based on CALIFA. The dotted line indicates the physical scale at given redshift corresponding to 2.5 arcsec, the MaNGA spatial resolution. of galaxy size (or mass) with redshift, which is a combined effect of the fixed field of view of MaNGA IFUs and the sample selection strategy. The IFUs of MaNGA have five different sizes, and the galaxy targets are selected so as to have their angular size matched to the IFU, so that the IFU covers out to 1.5 or 2.5R e of each galaxy. Due to this selection, the average luminosity, stellar mass, and size of the MaNGA sample is an increasing function of redshift (see Fig. 1). This leads to the lack of turnover galaxies with large R t at lower redshifts as seen in Fig. A1, where the physical size of the galaxies and their central features (bulges, bars, turnovers, etc.) are smaller. In the figure, we also show the turnover galaxies from the CALIFA survey (Paper I), which has better resolution than MaNGA and is limited to even lower redshifts. The same correlation is seen between R t and redshift. The CALIFA also uses IFU with a fixed FoV of ∼1 arcmin, so it is also expected to see correlations of the (physical) sizes of galaxies and their internal features/structures with the redshift.
In the lower panel of Fig. A1, we show the ratio of R t relative to the effective radius R e . We find no obvious correlation with the redshift, demonstrating that the correlation of R t with redshift in the upper panel is not real. In addition, we have done a analysis by re-constructing the control sample by closely matching the control sample with the turnover galaxy sample in redshift, in addition to stellar mass and colour. Our results and conclusions remain unchanged. Therefore, we conclude that our results are not affected by the spatial resolution of the data. This paper has been typeset from a T E X/L A T E X file prepared by the author.