(abridged) A fundamental subject in Extragalactic Astronomy concerns the formation and evolution of late-type galaxies (LTGs). The standard scenario comprises the early assembly of the bulge followed by disk accretion. However, recent observational evidence points to a joint formation and perpetual co-evolution of these structural components. Our current knowledge on the properties of bulge and disk is mostly founded on photometric decomposition studies, which sensitively depend on the adopted methodology and enclosed assumptions on the structure of LTGs. A critical assumption whose validity was never questioned is that galactic disks conserve their exponential nature up to the galactic center. This implies that bulge and disk co-exist without significant dynamical interaction and mass exchange over nearly the entire Hubble time. Our goal is to examine the validity of the standard assumption that galactic disks preserve their exponential intensity profile inside the bulge radius all the way to the galactic center. We developed a spectrophotometric bulge-disk decomposition technique that provides an estimation for the net spectrum of the bulge. A systematic application of our spectrophotometric bulge-disk decomposition tool to a representative sample of 135 local LTGs from the CALIFA Survey yields a significant fraction (up to ~30%) of unphysical net-bulge spectra when a purely exponential intensity profile is assumed for the disk. The obtained results suggest that, for a significant fraction of LTGs, the disk component shows a down-bending beneath the bulge. If proven to be true, such result will call for a substantial revision of structural decomposition studies for LTGs and have far-reaching implications in our understanding of the photometric properties of their bulges.

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.

Taking advantage of HST CANDELS data, we analyze the lowest redshift (z<0.5) massive galaxies in order to disentangle their structural constituents and study possible faint non-axis-symmetric features. Due to the excellent HST spatial resolution for intermediate-z objects, they are hard to model by purely automatic parametric fitting algorithms. We performed careful single and double S\'ersic fits to their galaxy surface brightness profiles. We also compare the model color profiles with the observed ones and also derive multi-component global effective radii attempting to obtain a better interpretation of the mass-size relation. Additionally, we test the robustness of our measured structural parameters via simulations. We find that the S\'ersic index does not offer a good proxy for the visual morphological type for our sample of massive galaxies. Our derived multi-component effective radii give a better description of the size of our sample galaxies than those inferred from single S\'ersic models with GALFIT. Our galaxy population lays on the scatter of the local mass-size relation, indicating that these massive galaxies do not experience a significant growth in size since z~0.5. Interestingly the few outliers are late-type galaxies, indicating that spheroids must reach the local mass-size relation earlier. For most of our sample galaxies, both single and multi-component S\'ersic models with GALFIT show substantial systematic deviations from the observed SBPs in the outskirts. These residuals may be partly due to several factors, namely a non-optimal data reduction for low surface brightness features, the existence of prominent stellar haloes for massive galaxies and could also arise from conceptual shortcomings of parametric 2D image decomposition tools. They consequently propagate into galaxy color profiles.

The S\'ersic law (SL) offers a versatile functional form for the structural characterization of galaxies near and far. Whereas applying it to galaxies with a genuine SL luminosity distribution yields a robust determination of the S\'ersic exponent eta and effective surface brightness $\mu_{\rm eff}$, this is not necessarily the case for galaxies whose surface brightness profiles (SBPs) appreciably deviate from the SL (eg, early-type galaxies with a depleted core and nucleated dwarf ellipticals, or most late-type galaxies-LTGs). In this general case of "imperfect" SL profiles, the best-fitting solution may significantly depend on the radius (or surface brightness) interval fit and corrections for point spread function (PSF) convolution effects. Such uncertainties may then affect, in a non-easily predictable manner, automated structural studies of galaxies. We present a fitting concept (iFIT) that permits a robust determination of the equivalent SL model for the general case of galaxies with imperfect SL profiles. iFIT has been extensively tested on synthetic data with a S\'ersic index 0.3<${\eta}$<4.2 and an effective radius 1<$\rm{R}_{eff}$ (arcs)<20. Applied to non PSF-convolved data, iFIT can infer the S\'ersic exponent eta with an absolute error of <0.2 even for shallow SBPs. As for PSF-degraded data, iFIT can recover the input SL model parameters with a satisfactorily accuracy almost over the entire considered parameter space as long as FWHM(PSF)<$\rm{R}_{eff}$. Tests indicate that iFIT shows little sensitivity on PSF corrections and the SBP limiting surface brightness, and that subtraction of the best-fitting SL model in two different bands yields a good match to the observed radial color profile. The publicly available iFIT offers an efficient tool for the non-supervised structural characterization of large galaxy samples, as those expected to become available with Euclid and LSST.

(abridged) This study revolves around dmB, a new distance- and extinction-independent measure of the contribution by stellar populations older than 9 Gyr to the mean r-band surface brightness of the bulge component in 135 late-type galaxies (LTGs) from the CALIFA survey, spanning a range of 2.6 dex and 3 dex in total and bulge stellar mass (M*T~10^(8.9-11.5) M_solar and M*B~10^(8.3-11.3) M_solar, respectively). The main insight from this study is that LTG bulges form a continuous sequence of increasing dmB with increasing M*T, M*B, stellar mass surface density S* and mass-weighted age and metallicity: high-dmB bulges are the oldest, densest and most massive ones, and vice versa. Furthermore, we find that the bulge-to-disk age and metallicity contrast, as well as the bulge-to-disk mass ratio increase with M*T, raising from, respectively, ~0 Gyr, 0 dex and 0.25 to ~3 Gyr, ~0.3 dex and 0.67 across the mass range covered by our sample. Whereas gas excitation in lower-mass bulges is invariably dominated by star formation (SF), LINER- and Seyfert-specific emission-line ratios were exclusively documented in high-mass, high-S* bulges. The continuity both in the properties of LTG bulges themselves and in their age and metallicity contrast to their parent disks suggests that these components evolve alongside in a concurrent process that leads to a continuum of physical and evolutionary characteristics. Our results are consistent with a picture where bulge growth in LTGs is driven by a superposition of quick-early and slow-secular processes, the relative importance of which increases with M*T. These processes, which presumably combine in situ SF in the bulge and inward migration of material from the disk, are expected to lead to a non-homologous radial growth of S* and a trend for an increasing Sersic index with increasing galaxy mass.

The morphological, spectroscopic and kinematical properties of the warm interstellar medium (wim) in early-type galaxies (ETGs) hold key observational constraints to nuclear activity and the buildup history of these massive, quiescent systems. High-quality integral field spectroscopy (IFS) data with a wide spectral and spatial coverage, such as those from the CALIFA survey, offer an unprecedented opportunity for advancing our understanding of the wim in ETGs. This article centers on a 2D investigation of the wim component in 32 nearby (<~150Mpc) ETGs from CALIFA, complementing a previous 1D analysis of the same sample (Papaderos et al. 2013; P13). We include here H\alpha\ intensity and equivalent width (EW) maps and radial profiles, diagnostic emission-line ratios, besides ionized-gas and stellar kinematics. This study is supplemented by \tau-ratio maps as an efficient means to quantify the role of photoionization by pAGB stars, as compared to other mechanisms (e.g., AGN, low-level star formation). Additionally, we extend the tentative classification proposed in P13 by the type i+, which is assigned to a subset of type i ETGs exhibiting ongoing low-level star-formation (SF) in their periphery. This finding along with faint traces of localized SF in the extranuclear component of several of our sample ETGs points to a non-negligible contribution by OB stars to the total ionizing budget. We also demonstrate that, at the typical emission-line detection threshold of ~2\AA in previous studies, most of the extranuclear wim emission in an ETG may evade detection, which could in turn prompt its classification as an entirely gas-devoid system. This study adds further observational evidence for a considerable heterogeneity among ETGs with regard to the physical properties and 2D kinematics of the wim component, and underscores the importance of IFS studies over their entire optical extent.

Integral field spectroscopy studies based on CALIFA data have recently revealed the presence of ongoing low-level star formation (SF) in the periphery of ~10% of local early-type galaxies (ETGs), witnessing a still ongoing inside-out galaxy growth process. A distinctive property of the nebular component in these ETGs, classified i+, is a two-radial-zone structure, with the inner zone displaying LINER emission with a H\alpha equivalent width EW~1{\AA}, and the outer one (3{\AA}<EW<~20{\AA}) showing HII-region characteristics. Using CALIFA IFS data, we empirically demonstrate that the confinement of nebular emission to the galaxy periphery leads to a strong aperture (or, redshift) bias in spectroscopic single-fiber studies of type i+ ETGs: At low redshift (<~0.45), SDSS spectroscopy is restricted to the inner (SF-devoid LINER) zone, thereby leading to their erroneous classification as "retired" galaxies (systems lacking SF and whose faint emission is powered by pAGB stars). Only at higher z's the SDSS aperture can encompass the outer SF zone, permitting their unbiased classification as "composite SF/LINER". We also demonstrate that the principal effect of a decreasing aperture on the classification of i+ ETGs via standard BPT emission-line ratios consists in a monotonic up-right shift precisely along the upper-right wing of the "seagull" distribution. Motivated by these insights, we also investigate theoretically these biases in aperture-limited studies of inside-out growing galaxies as a function of z. To this end, we devise a simple model, which involves an outwardly propagating SF process, that reproduces the radial extent and two-zone EW distribution of i+ ETGs. By simulating on this model the spectroscopic SDSS aperture, we find that SDSS studies at z<~1 are progressively restricted to the inner LINER-zone, and miss an increasingly large portion of the H\alpha-emitting periphery.

Based on a combined analysis of SDSS imaging and CALIFA integral field spectroscopy data, we report on the detection of faint (24 < {\mu}$_r$ mag/arcsec$^2$ < 26) star-forming spiral-arm-like features in the periphery of three nearby early-type galaxies (ETGs). These features are of considerable interest because they document the still ongoing inside-out growth of some local ETGs and may add valuable observational insight into the origin and evolution of spiral structure in triaxial stellar systems. A characteristic property of the nebular component in the studied ETGs, classified i+, is a two-radial-zone structure, with the inner zone that displays faint (EW(H\alpha)$\simeq$1{\AA}) low-ionization nuclear emission-line region (LINER) properties, and the outer one (3{\AA}<EW(H\alpha)<~20{\AA}) HII-region characteristics. This spatial segregation of nebular emission in two physically distinct concentric zones calls for an examination of aperture effects in studies of type i+ ETGs with single-fiber spectroscopic data.