We present single-Sérsic two-dimensional (2D) model fits to 167 600 galaxies modelled independently in the ugrizYJHK bandpasses using reprocessed Sloan Digital Sky Survey Data Release Seven (SDSS ...DR7) and UKIRT Infrared Deep Sky Survey Large Area Survey imaging data available from the Galaxy And Mass Assembly (GAMA) data base. In order to facilitate this study we developed Structural Investigation of Galaxies via Model Analysis (sigma), an r wrapper around several contemporary astronomy software packages including source extractor, psf extractor and galfit 3. sigma produces realistic 2D model fits to galaxies, employing automatic adaptive background subtraction and empirical point spread function measurements on the fly for each galaxy in GAMA. Using these results, we define a common coverage area across the three GAMA regions containing 138 269 galaxies. We provide Sérsic magnitudes truncated at 10r
e which show good agreement with SDSS Petrosian and GAMA photometry for low Sérsic index systems (n < 4), and much improved photometry for high Sérsic index systems (n > 4), recovering as much as Δm= 0.5 mag in the r band. We employ a K-band Sérsic index/u−r colour relation to delineate the massive (n > ∼2) early-type galaxies (ETGs) from the late-type galaxies (LTGs). The mean Sérsic index of these ETGs shows a smooth variation with wavelength, increasing by 30 per cent from g through K. LTGs exhibit a more extreme change in Sérsic index, increasing by 52 per cent across the same range. In addition, ETGs and LTGs exhibit a 38 and 25 per cent decrease, respectively, in half-light radius from g through K. These trends are shown to arise due to the effects of dust attenuation and stellar population/metallicity gradients within galaxy populations.
In this paper, we demonstrate a new method for fitting galaxy profiles which makes use of the full multiwavelength data provided by modern large optical-near-infrared imaging surveys. We present a ...new version of galapagos, which utilizes a recently developed multiwavelength version of galfit, and enables the automated measurement of wavelength-dependent Sérsic profile parameters for very large samples of galaxies. Our new technique is extensively tested to assess the reliability of both pieces of software, galfit and galapagos on both real ugrizY JHK imaging data from the Galaxy And Mass Assembly survey and simulated data made to the same specifications. We find that fitting galaxy light profiles with multiwavelength data increases the stability and accuracy of the measured parameters, and hence produces more complete and meaningful multiwavelength photometry than has been available previously. The improvement is particularly significant for magnitudes in low-S/N bands and for structural parameters like half-light radius r
e and Sérsic index n for which a prior is used by constraining these parameters to a polynomial as a function of wavelength. This allows the fitting routines to push the magnitude of galaxies for which sensible values can be derived to fainter limits. The technique utilizes a smooth transition of galaxy parameters with wavelength, creating more physically meaningful transitions than single-band fitting and allows accurate interpolation between passbands, perfect for derivation of rest-frame values.
We describe the calculation of the stochastically heated dust emission using the 3D ray-tracing dust radiative transfer code dart-ray, which is designed to solve the dust radiative transfer problem ...for galaxies with arbitrary geometries. In order to reduce the time required to derive the non-equilibrium dust emission spectra from each volume element within a model, we implemented an adaptive spectral energy distribution library approach, which we tested for the case of axisymmetric galaxy geometries. To show the capabilities of the code, we applied dart-ray to a high-resolution N-body+SPH galaxy simulation to predict the appearance of the simulated galaxy at a set of wavelengths from the UV to the sub-mm. We analyse the results to determine the effect of dust on the observed radial and vertical profiles of the stellar emission as well as on the attenuation and scattering of light from the constituent stellar populations. We also quantify the proportion of dust re-radiated stellar light powered by young and old stellar populations, both bolometrically and as a function of infrared wavelength.
We report the morphological classification of 3727 galaxies from the Galaxy and Mass Assembly survey with M
r
< −17.4 mag and in the redshift range 0.025 < z < 0.06 (2.1 × 105 Mpc3) into E, S0-Sa, ...SB0-SBa, Sab-Scd, SBab-SBcd, Sd-Irr and little blue spheroid classes. Approximately 70 per cent of galaxies in our sample are disc-dominated systems, with the remaining ∼30 per cent spheroid dominated. We establish the robustness of our classifications, and use them to derive morphological-type luminosity functions and luminosity densities in the ugrizYJHK passbands, improving on prior studies that split by global colour or light profile shape alone. We find that the total galaxy luminosity function is best described by a double-Schechter function while the constituent morphological-type luminosity functions are well described by a single-Schechter function. These data are also used to derive the star formation rate densities for each Hubble class, and the attenuated and unattenuated (corrected for dust) cosmic spectral energy distributions, i.e. the instantaneous energy production budget. While the observed optical/near-IR energy budget is dominated 58:42 by galaxies with a significant spheroidal component, the actual energy production rate is reversed, i.e. the combined disc-dominated populations generate ∼1.3 times as much energy as the spheroid-dominated populations. On the grandest scale, this implies that chemical evolution in the local Universe is currently largely confined to mid-type spiral classes like our Milky Way.
We present results of a study to quantify the effects of dust on the derived photometric parameters of disk and bulges obtained from bulge-disk decomposition: scale-length, effective radius, Sérsic ...index, disk axis-ratio, and bulge-to-disk ratio. The dust induced changes in these parameters were obtained by fitting simulated images of composite systems (containing a disk and a bulge) produced using radiative transfer calculations. The simulations were fitted with the GALFIT 3.0.2 data analysis algorithm. Fits were done with both a combination of an exponential plus a variable-index Sérsic function as well as with a combination of two variable-index Sérsic functions. We find that dust is biasing the derived exponential scale-length of decomposed disks towards smaller values than would be otherwise derived if the galaxy were to have no bulge. Similarly, the derived bulge-to-disk ratio is biased towards smaller values. However, the derived axis-ratio of the disk is not changed in the decomposition process. The derived effective radius of decomposed disks of systems having exponential bulges is found to be less affected by dust when fits are done with two variable-index Sérsic functions. For the same type of fits dust is found to bias the value of the derived effective radius of decomposed disks towards lower values for systems having de Vaucouleurs bulges.
A radiative transfer model for the spiral galaxy M33 Thirlwall, Jordan J; Popescu, Cristina C; Tuffs, Richard J ...
Monthly Notices of the Royal Astronomical Society,
06/2020, Letnik:
495, Številka:
1
Journal Article
Recenzirano
Odprti dostop
ABSTRACT
We present the first radiative transfer (RT) model of a non-edge-on disc galaxy in which the large-scale geometry of stars and dust is self-consistently derived through the fitting of ...multiwavelength imaging observations from the ultraviolet to the submm. To this end, we used the axisymmetric RT model of Popescu et al. and a new methodology for deriving geometrical parameters, and applied this to decode the spectral energy distribution (SED) of M33. We successfully account for both the spatial and spectral energy distribution, with residuals typically within $7{{\ \rm per\ cent}}$ in the profiles of surface brightness and within $8{{\ \rm per\ cent}}$ in the spatially integrated SED. We predict well the energy balance between absorption and re-emission by dust, with no need to invoke modified grain properties, and we find no submm emission that is in excess of our model predictions. We calculate that $80\pm 8{{\ \rm per\ cent}}$ of the dust heating is powered by the young stellar populations. We identify several morphological components in M33, a nuclear, an inner, a main and an outer disc, showing a monotonic trend in decreasing star formation surface density (ΣSFR) from the nuclear to the outer disc. In relation to surface density of stellar mass, the ΣSFR of these components defines a steeper relation than the ‘main sequence’ of star-forming galaxies, which we call a ‘structurally resolved main sequence’. Either environmental or stellar feedback mechanisms could explain the slope of the newly defined sequence. We find the star formation rate to be ${\rm SFR}=0.28^{+0.02}_{-0.01}{\rm M}_{\odot }{\rm yr}^{-1}$.
We present an estimate of the galaxy stellar mass function and its division by morphological type in the local (0.025 < z < 0.06) Universe. Adopting robust morphological classifications as previously ...presented (Kelvin et al.) for a sample of 3727 galaxies taken from the Galaxy And Mass Assembly survey, we define a local volume and stellar mass limited sub-sample of 2711 galaxies to a lower stellar mass limit of
$\mathcal {M}=10^{9.0}\,{\rm M}_{{\odot }}$
. We confirm that the galaxy stellar mass function is well described by a double-Schechter function given by
$\mathcal {M}^{*}=10^{10.64}\,{\rm M}_{{\odot }}$
, α1 = −0.43,
$\phi _{1}^{*}=4.18\;\mathrm{dex}^{-1}\,\mathrm{Mpc}^{-3}$
, α2 = −1.50 and
$\phi _{2}^{*}=0.74\;\mathrm{dex}^{-1}\,\mathrm{Mpc}^{-3}$
. The constituent morphological-type stellar mass functions are well sampled above our lower stellar mass limit, excepting the faint little blue spheroid population of galaxies. We find approximately
$71{}_{-4}^{+3}$
per cent of the stellar mass in the local Universe is found within spheroid-dominated galaxies; ellipticals and S0-Sas. The remaining
$29{}_{-3}^{+4}$
per cent falls predominantly within late-type disc-dominated systems, Sab-Scds and Sd-Irrs. Adopting reasonable bulge-to-total ratios implies that approximately half the stellar mass today resides in spheroidal structures, and half in disc structures. Within this local sample, we find approximate stellar mass proportions for E : S0-Sa : Sab-Scd : Sd-Irr of 34 : 37 : 24 : 5.
We build, as far as theory will permit, self-consistent model H II regions around central clusters of aging stars. These produce strong emission line diagnostics applicable to either individual H II ...regions in galaxies or to the integrated emission line spectra of disk or starburst galaxies. The models assume that the expansion and internal pressure of individual H II regions is driven by the net input of mechanical energy from the central cluster, be it through winds or supernova events. This eliminates the ionization parameter as a free variable, replacing it with a parameter that depends on the ratio of the cluster mass to the pressure in the surrounding interstellar medium. These models explain why H II regions with low abundances have high excitation and demonstrate that at least part of the warm ionized medium is the result of overlapping faint, old, large, and low-pressure H II regions. We present line ratios (at both optical and IR wavelengths) that provide reliable abundance diagnostics for both single H II regions or for integrated galaxy spectra, and we find a number that can be used to estimate the mean age of the cluster stars exciting individual H II regions.
We investigate the relationship between colour and structure within galaxies using a large, volume-limited sample of bright, low-redshift galaxies with optical–near-infrared imaging from the Galaxy ...And Mass Assembly survey. We fit single-component, wavelength-dependent, elliptical Sérsic models to all passbands simultaneously, using software developed by the MegaMorph project. Dividing our sample by n and colour, the recovered wavelength variations in effective radius (R
e) and Sérsic index (n) reveal the internal structure, and hence formation history, of different types of galaxies. All these trends depend on n; some have an additional dependence on galaxy colour. Late-type galaxies (n
r
< 2.5) show a dramatic increase in Sérsic index with wavelength. This might be a result of their two-component (bulge–disc) nature, though stellar population gradients within each component and dust attenuation are likely to play a role. All galaxies show a substantial decrease in R
e with wavelength. This is strongest for early types (n
r
> 2.5), even though they maintain constant n with wavelength, revealing that ellipticals are a superimposition of different stellar populations associated with multiple collapse and merging events. Processes leading to structures with larger R
e must be associated with lower metallicity or younger stellar populations. This appears to rule out the formation of young cores through dissipative gas accretion as an important mechanism in the recent lives of luminous elliptical galaxies.
Based on our sample of 10 095 galaxies with bulge–disc decompositions we derive the empiricalBMGC-band internal attenuation–inclination relation for galaxy discs and their associated central bulges. ...Our results agree well with the independently derived dust models of Tuffs et al., leading to a direct constraint on the mean opacity of spiral discs of τfB= 3.8 ± 0.7 (central face-on BMGC-band opacity). Depending on inclination, the BMGC-band attenuation correction varies from 0.2 to 1.1 mag for discs and from 0.8 to 2.6 mag for bulges. We find that, overall, 37 per cent of all BMGC-band photons produced in discs in the nearby Universe are absorbed by dust, a figure that rises to 71 per cent for bulge photons. The severity of internal dust extinction is such that one must incorporate internal dust corrections in all optical studies of large galaxy samples. This is particularly pertinent for optical Hubble Space Telescope comparative evolutionary studies as the dust properties will also be evolving. We use the new results to revise our recent estimates of the spheroid and disc luminosity functions. The implied stellar mass densities at redshift zero are somewhat higher than our earlier estimates: ρdiscs= (3.8 ± 0.6) → (4.4 ± 0.6) × 108 h M⊙ Mpc−3 and ρbulges= (1.6 ± 0.4) → (2.2 ± 0.4) × 108 h M⊙ Mpc−3. From our best-fitting dust models we derive a redshift zero cosmic dust density of ρdust≈ (5.3 ± 1.7) × 105 h M⊙ Mpc−3. This implies that (0.0083 ± 0.0027)h per cent of the baryons in the Universe are in the form of dust and (11.9 ± 1.7) h per cent (Salpeter-‘lite’ initial mass function) are in the form of stars (∼58 per cent reside in galaxy discs, ∼10 per cent in red elliptical galaxies, ∼29 per cent in classical galaxy bulges and the remainder in low-luminosity blue spheroid systems/components).