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 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.
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, Volume:
495, Issue:
1
Journal Article
Peer reviewed
Open access
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 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.
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 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 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.
We provide and describe a library of diffuse stellar radiation fields in spiral galaxies derived using calculations of the transfer of stellar radiation from the main morphological components - ...discs, thin discs and bulges - through the dusty interstellar medium. These radiation fields are self-consistent with the solutions for the integrated panchromatic spectral energy distributions (SEDs) previously presented using the same model. Because of this, observables calculated from the radiation fields, such as gamma-ray or radio emission, can be self-consistently combined with the solutions for the ultraviolet/optical/submillimetre SEDs, thus expanding the range of applicability of the radiation transfer model to a broader range of wavelengths and physical quantities. We also give analytic solutions for radiation fields in optically thin stellar discs, in stellar discs with one dust disc and in stellar discs with two dust discs. The analytic solutions for the direct light are exact and can be used as benchmarks. The analytic solutions with scattering are only approximate, becoming exact only in the extreme optically thick limit. We find strongly contrasting solutions for the spatial distribution of the radiation fields for discs, thin discs and bulges. For bulges, we find a strong dependence of the radiation fields on the Sérsic index.
We look for correlated changes in stellar mass and star formation rate (SFR) along filaments in the cosmic web by examining the stellar masses and UV-derived SFRs of 1799 ungrouped and unpaired ...spiral galaxies that reside in filaments. We devise multiple distance metrics to characterize the complex geometry of filaments, and find that galaxies closer to the cylindrical centre of a filament have higher stellar masses than their counterparts near the periphery of filaments, on the edges of voids. In addition, these peripheral spiral galaxies have higher SFRs at a given mass. Complementing our sample of filament spiral galaxies with spiral galaxies in tendrils and voids, we find that the average SFR of these objects in different large-scale environments are similar to each other with the primary discriminant in SFR being stellar mass, in line with previous works. However, the distributions of SFRs are found to vary with large-scale environment. Our results thus suggest a model in which in addition to stellar mass as the primary discriminant, the large-scale environment is imprinted in the SFR as a second-order effect. Furthermore, our detailed results for filament galaxies suggest a model in which gas accretion from voids on to filaments is primarily in an orthogonal direction. Overall, we find our results to be in line with theoretical expectations of the thermodynamic properties of the intergalactic medium in different large-scale environments.