Radial age gradients hold the cumulative record of the multitude of physical processes driving the build-up of stellar populations and the ensuing star formation (SF) quenching process in galaxy bulges, therefore potentially sensitive discriminators between competing theoretical concepts on bulge formation and evolution. Based on spectral modeling of integral field spectroscopy data from the CALIFA survey, we derive mass- and light-weighted stellar age gradients ($\nabla$(t,B)L,M) within the photometrically determined bulge radius (RB) of a representative sample of local face-on late-type galaxies that span 2.6 dex in stellar mass. Our analysis documents a trend for decreasing $\nabla$(t,B)L,M with increasing M,T, with high-mass bulges predominantly showing negative age gradients and vice versa. The inversion from positive to negative $\nabla$(t,B)L,M occurs at logM,T ~ 10, which roughly coincides with the transition from lower-mass bulges whose gas excitation is powered by SF to bulges classified as Composite, LINER or Seyfert. We discuss two limiting cases for the origin of radial age gradients in massive LTG bulges. The first assumes that the stellar age in the bulge is initially spatially uniform, thus the observed age gradients arise from an inside-out SF quenching (ioSFQ) front that is radially expanding with a mean velocity vq. In this case, the age gradients translate into a slow ioSFQ that lasts until z~2, suggesting mild negative feedback by SF or an AGN. If negative age gradients in massive bulges are not due to ioSFQ but primarily due to their inside-out formation process, then the standard hypothesis of quasi-monolithic bulge formation has to be discarded in favor of a scenario that involves gradual buildup of stellar mass over 2-3 Gyr through, e.g., inside-out SF and inward migration of SF clumps from the disk. In this case, rapid AGN-driven ioSFQ cannot be ruled out.

Spectral population synthesis (PS) is a fundamental tool in extragalactic research that aims to decipher the assembly history of galaxies from their SED. However, until recently all PS codes were restricted to purely stellar fits, neglecting the essential contribution of nebular emission (NE). With the advent of FADO, the now possible self-consistent modelling of stellar and NE opens new routes to the exploration of galaxy SFHs. The main goal of this study is to quantitatively explore the accuracy to which FADO can recover physical and evolutionary properties of galaxies and compare its output with that from purely stellar PS codes. With this in mind, FADO and STARLIGHT were applied to synthetic SEDs that track the spectral evolution of stars and gas in extinction-free mock galaxies that form their stellar mass ($M_\star$) according to different parametric SFHs. Spectral fits were computed for two different set-ups that approximate the spectral range of SDSS and CALIFA data. Our analysis indicates that FADO can recover the key physical and evolutionary properties of galaxies, such as $M_\star$ and mass- and light-weighted mean age and metallicity, with an accuracy better than 0.2 dex. This is the case even in phases of strongly elevated sSFR and thus with considerable NE contamination. As for STARLIGHT, our analysis documents a moderately good agreement with theoretical values only for evolutionary phases for which NE drops to low levels. Indeed, fits with STARLIGHT during phases of high sSFR severely overestimate both $M_\star$ and the mass-weighted stellar age, whereas strongly underestimate the light-weighted age and metallicity. The insights from this study suggest that the neglect of nebular continuum emission in STARLIGHT and similar purely stellar PS codes could systematically impact $M_\star$ and SFH estimates for star-forming galaxies.

The effect of the featureless power-law (PL) continuum of an active galactic nucleus (AGN) on the estimation of physical properties of galaxies with optical population spectral synthesis (PSS) remains largely unknown. With this in mind, we fit synthetic galaxy spectra representing a wide range of galaxy star formation histories (SFHs) and including distinct PL contributions of the form $F_{\nu} \propto \nu^{-\alpha}$ with the PSS code STARLIGHT to study to which extent various inferred quantities (e.g. stellar mass, mean age, and mean metallicity) match the input. The synthetic spectral energy distributions (SEDs) computed with our evolutionary spectral synthesis code include an AGN PL component with $0.5 \leq \alpha \leq 2$ and a fractional contribution $0.2 \leq x_{\mathrm{AGN}} \leq 0.8$ to the monochromatic flux at 4020 \AA. At the empirical AGN detection threshold $x_{\mathrm{AGN}}\simeq 0.26$ that we previously inferred in a pilot study on this subject, our results show that the neglect of a PL component in spectral fitting can lead to an overestimation by $\sim$2 dex in stellar mass and by up to $\sim$1 and $\sim$4 dex in the light- and mass-weighted mean stellar age, respectively, whereas the light- and mass-weighted mean stellar metallicity are underestimated by up to $\sim$0.3 and $\sim$0.6 dex, respectively. Other fitting set-ups including either a single PL or multiple PLs in the base reveal, on average, much lower unsystematic uncertainties of the order of those typically found when fitting purely stellar SEDs with stellar templates, however, reaching locally up to $\sim$1, 3 and 0.4 dex in mass, age and metallicity, respectively. Our results underscore the importance of an accurate modelling of the AGN spectral contribution in PSS fits as a minimum requirement for the recovery of the physical and evolutionary properties of stellar populations in active galaxies.

Various lines of evidence suggest that the cores of a large portion of early-type galaxies (ETGs) are virtually evacuated of warm ionised gas. This implies that the Lyman-continuum (LyC) radiation produced by an assumed active galactic nucleus (AGN) can escape from the nuclei of these systems without being locally reprocessed into nebular emission, which would prevent their reliable spectroscopic classification as Seyfert galaxies with standard diagnostic emission-line ratios. The spectral energy distribution (SED) of these ETGs would then lack nebular emission and be essentially composed of an old stellar component and the featureless power-law (PL) continuum from the AGN. A question that arises in this context is whether the AGN component can be detected with current spectral population synthesis in the optical, specifically, whether these techniques effectively place an AGN detection threshold in LyC-leaking galaxies. To quantitatively address this question, we took a combined approach that involves spectral fitting with STARLIGHT of synthetic SEDs composed of stellar emission that characterises a 10 Gyr old ETG and an AGN power-law component that contributes a fraction $0\leq x_{\mathrm{AGN}} < 1$ of the monochromatic luminosity at $\lambda_0=$ 4020 \AA. In addition to a set of fits for PL distributions $F_{\nu} \propto \nu^{-\alpha}$ with the canonical $\alpha=1.5$, we used a base of multiple PLs with $0.5 \leq \alpha \leq 2$ for a grid of synthetic SEDs with a signal-to-noise ratio of 5-$10^3$. Our analysis indicates an effective AGN detection threshold at $x_{\mathrm{AGN}}\simeq 0.26$, which suggests that a considerable fraction of ETGs hosting significant accretion-powered nuclear activity may be missing in the AGN demographics.

We observe the radio galaxy PKS 1934-63 (at $z=0.1825$) using MUSE (Multi Unit Spectroscopic Explorer) on the Very Large Telescope (VLT). The radio source is GigaHertz Peaked Spectrum and compact (0.13 kpc), implying an early stage of evolution ($\leq 10^4$ yr). Our data show an interacting pair of galaxies, projected separation 9.1 kpc, velocity difference $\Delta(v)=216$ km $\rm s^{-1}$. The larger galaxy is a $\rm M_{*}\simeq 10^{11}M_{\odot}$ spheroidal with the emission-line spectrum of a high-excitation young radio AGN, e.g. strong [OI]6300 and [OIII]5007. Emission-line ratios indicate a large contribution to the line luminosity from high-velocity shocks ($\simeq 550$ km $\rm s^{-1}$) . The companion is a non-AGN disk galaxy, with extended $\rm H\alpha$ emission from which its star-formation rate is estimated as $\rm 0.61~M_{\odot}yr^{-1}$. Both galaxies show rotational velocity gradients in $\rm H\alpha$ and other lines, with the interaction being prograde-prograde. The SE-NW velocity gradient of the AGN host is misaligned from the E-W radio axis, but aligned with a previously discovered central ultraviolet source, and a factor 2 greater in amplitude in $\rm H\alpha$ than in other (forbidden) lines (e.g. [OIII]5007). This could be produced by a fast rotating (100-150 km $\rm s^{-1}$) disk with circumnuclear star-formation. We also identify a broad component of [OIII]5007 emission, blueshifted with a velocity gradient aligned with the radio jets, and associated with outflow. However, the broad component of [OI]6300 is redshifted. In spectral fits, both galaxies have old stellar populations plus $\sim 0.1\%$ of very young stars, consistent with the galaxies undergoing first perigalacticon, triggering infall and star-formation from $\sim 40$ Myr ago followed by the radio outburst.