A new fitting concept for the robust determination of Sérsic model parameters

• Published in: Astronomy and astrophysics (Berlin) 2019-12, Vol.632, p.A128
• Main Authors: Breda, Iris, Papaderos, Polychronis, Gomes, Jean Michel, Amarantidis, Stergios
• Format: Article
• Language: English
• Subjects: Convolution; Data points; Galactic evolution; galaxies: elliptical and lenticular; galaxies: fundamental parameters; galaxies: photometry; galaxies: starburst; galaxies: structure; Galaxy distribution; Luminosity; Mathematical models; Model accuracy; Parameters; Photometry; Point spread functions; Red shift; Robustness; Starburst galaxies; Stars & galaxies; Structural analysis; Subtraction; Surface brightness; techniques: photometric; Uncertainty
• ISSN: 0004-6361; 1432-0746
• DOI: 10.1051/0004-6361/201935144

Context. The Sérsic law (SL) offers a versatile, widely used functional form for the structural characterization of galaxies near and far. Whereas fitting this three-parameter function to galaxies with a genuine SL luminosity distribution (e.g., several local early-type galaxies–ETGs) yields a robust determination of the Sérsic exponent η and effective surface brightness μeff, this is not necessarily the case for galaxies whose surface brightness profiles (SBPs) appreciably deviate, either in their centers or over an extended radius interval, from the SL (e.g., ETGs 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, the photometric passbands considered and the specifics of the fitting procedure (photometric uncertainties of SBP data points or image pixels, and corrections for point spread function (PSF) convolution effects). Such uncertainties may then affect, in a non-easily predictable manner, automated structural studies of large heterogeneous galaxy samples and introduce a scatter, if not a bias, in galaxy scaling relations and their evolution across redshift (z). Aims. Our goal is to devise a fitting concept that permits a robust determination of the equivalent SL model for the general case of galaxies with imperfect SL profiles. Methods. The distinctive feature of the concept proposed here (iFIT) is that the fit is not constrained through standard χ2 minimization between an observed SBP and the SL model of it, but instead through the search for the best match between the observationally determined and theoretically expected radial variation of the mean surface brightness and light growth curve. This approach ensures quick convergence to a unique solution for both perfect and imperfect Sérsic profiles, even shallow and resolution-degraded SBPs. iFIT allows for correction of PSF convolution effects, offering the user the option of choosing between a Moffat, Gaussian, or user-supplied PSF. iFIT, which is a standalone FORTRAN code, can be applied to any SBP that is provided in ASCII format and it has the capability of convenient graphical storage of its output. The iFIT distribution package is supplemented with an auxiliary SBP derivation tool in python. Results. iFIT has been extensively tested on synthetic data with a Sérsic index 0.3 ≤ η ≤ 4.2 and an effective rad