One of the most important properties of a galaxy is the total stellar mass, or equivalently the stellar mass-to-light ratio (M/L). It is not directly observable, but can be estimated from stellar population synthesis. Currently, a galaxy's M/L is typically estimated from global fluxes. For example, a single global g - i colour correlates well with the stellar M/L. Spectral energy distribution (SED) fitting can make use of all available fluxes and their errors to make a Bayesian estimate of the M/L. We want to investigate the possibility of using morphology information to assist predictions of M/L. Our first goal is to develop and train a method that only requires a g-band image and redshift as input. This will allows us to study the correlation between M/L and morphology. Next, we can also include the i-band flux, and determine if morphology provides additional constraints compared to a method that only uses g- and i-band fluxes. We used a machine learning pipeline that can be split in two steps. First, we detected morphology features with a convolutional neural network. These are then combined with redshift, pixel size and g-band luminosity features in a gradient boosting machine. Our training target was the M/L acquired from the GALEX-SDSS-WISE Legacy Catalog, which uses global SED fitting and contains galaxies with z ~ 0.1. Morphology is a useful attribute when no colour information is available, but can not outperform colour methods on its own. When we combine the morphology features with global g- and i-band luminosities, we find an improved estimate compared to a model which does not make use of morphology. While our method was trained to reproduce global SED fitted M/L, galaxy morphology gives us an important additional constraint when using one or two bands. Our framework can be extended to other problems to make use of morphological information.

We present results of the detailed dust energy balance study for the seven large edge-on galaxies in the HEROES sample using 3D radiative transfer (RT) modelling. Based on available optical and near-infrared observations of the HEROES galaxies, we derive the 3D distribution of stars and dust in these galaxies. For the sake of uniformity, we apply the same technique to retrieve galaxy properties for the entire sample: we use a stellar model consisting of a S\'ersic bulge and three double-exponential discs (a superthin disc for a young stellar population and thin and thick discs for old populations). For the dust component, we adopt a double-exponential disc with the new THEMIS dust-grain model. We fit oligochromatic radiative transfer (RT) models to the optical and near-infrared images with the fitting algorithm FitSKIRT and do panchromatic simulations with the SKIRT code at wavelengths ranging from ultraviolet to submillimeter. We confirm the previously stated dust energy balance problem in galaxies: for the HEROES galaxies, the dust emission derived from our RT calculations underestimates the real observations by a factor 1.5-4 for all galaxies except NGC 973 and NGC 5907 (apparently, the latter galaxy has a more complex geometry than we used). The comparison between our RT simulations and the observations at mid-infrared-submillimeter wavelengths shows that most of our galaxies exhibit complex dust morphologies (possible spiral arms, star-forming regions, more extended dust structure in the radial and vertical directions). We suggest that, in agreement with the results from Saftly et al. (2015), the large- and small-scale structure is the most probable explanation for the dust energy balance problem.

We present optical and near-infrared archival observations of eight edge-on galaxies. These observations are used to model the stellar content of each galaxy using the FitSKIRT software package. Using FitSKIRT, we can self-consistently model a galaxy in each band simultaneously while treating for dust. This allows us to accurately measure both the scale length and scale height of the stellar disc, plus the shape parameters of the bulge. By combining this data with the previously reported integrated magnitudes of each galaxy, we can infer their true luminosities. We have successfully modelled seven out of the eight galaxies in our sample. We find that stellar discs can be modelled correctly, but have not been able to model the stellar bulge reliably. Our sample consists for the most part of slow rotating galaxies, and we find that the average dust layer is much thicker than what is reported for faster rotating galaxies.

We investigate the dust energy balance for the edge-on galaxy IC 2531, one of the seven galaxies in the HEROES sample. We perform a state-of-the-art radiative transfer modelling based, for the first time, on a set of optical and near-infrared galaxy images. We show that taking into account near-infrared imaging in the modelling significantly improves the constraints on the retrieved parameters of the dust content. We confirm the result from previous studies that including a young stellar population in the modelling is important for explaining the observed stellar energy distribution. However, the discrepancy between the observed and modelled thermal emission at far-infrared wavelengths, the so-called dust energy balance problem, is still present: the model underestimates the observed fluxes by a factor of about two. We compare two different dust models, and find that dust parameters and thus the spectral energy distribution in the infrared domain are sensitive to the adopted dust model. In general, the THEMIS model reproduces the observed emission in the infrared wavelength domain better than the popular Zubko et al. BARE-GR-S model. Our study of IC 2531 is a pilot case for detailed and uniform radiative transfer modelling of the entire HEROES sample, which will shed more light on the strength and origins of the dust energy balance problem.

Context. Edge-on galaxies can offer important insights in galaxy evolution as they are the only systems where the distribution of the different components can be studied both radially and vertically. The HEROES project was designed to investigate the interplay between the gas, dust, stars and dark matter (DM) in a sample of 7 massive edge-on spiral galaxies. Aims. In this second HEROES paper we present an analysis of the atomic gas content of 6 out of 7 galaxies in our sample. The remaining galaxy was recently analysed according to the same strategy. The primary aim of this work is to constrain the surface density distribution, the rotation curve and the geometry of the gas disks in a homogeneous way. In addition we identify peculiar features and signs of recent interactions. Methods. We construct detailed tilted-ring models of the atomic gas disks based on new GMRT 21-cm observations of NGC 973 and UGC 4277 and re-reduced archival HI data of NGC 5907, NGC 5529, IC 2531 and NGC 4217. Potential degeneracies between different models are resolved by requiring a good agreement with the data in various representations of the data cubes. Results. From our modelling we find that all but one galaxy are warped along the major axis. In addition, we identify warps along the line of sight in three galaxies. A flaring gas layer is required to reproduce the data only for one galaxy, but (moderate) flares cannot be ruled for the other galaxies either. A coplanar ring-like structure is detected outside the main disk of NGC 4217, which we suggest could be the remnant of a recent minor merger event. We also find evidence for a radial inflow of 15 +- 5 km/s in the disk of NGC 5529, which might be related to the ongoing interaction with two nearby companions. (Abridged)

A fraction of the early-type galaxy population hosts a prominent dust lane. Methods to quantify the dust content of these systems based on optical imaging data usually yield dust masses which are an order of magnitude lower than dust masses derived from the observed FIR emission. High-quality optical data from the Next Generation Virgo cluster Survey (NGVS) and FIR/submm observations from the Herschel Virgo Cluster Survey (HeViCS) allow us to revisit previous methods to determine the dust content in galaxies and explore new ones. We aim to derive the dust mass in NGC 4370 from both optical and FIR data, and investigate the need to invoke a putative diffuse dust component. We create color and attenuation maps, which are converted to approximate dust mass maps based on simple dust geometries. Dust masses are also derived from SED fits to FIR/submm observations. Finally, inverse radiative transfer fitting is performed to investigate more complex dust geometries. The empirical methods applied to the optical data yield lower limits of 3.4e5 solar masses, an order of magnitude below the total dust masses derived from SED fitting. In contrast, radiative transfer models yield dust masses which are slightly lower, but fully consistent with the FIR-derived mass. Dust is more likely to be distributed in a ring around the centre of NGC 4370 as opposed to an exponential disc or a simple foreground screen. Moreover, using inverse radiative transfer fitting, we are able to constrain most of the parameters describing these geometries. The resulting dust masses are high enough to account for the dust observed at FIR/submm wavelengths, so that no diffuse dust component needs to be invoked. We furthermore caution for the interpretation of dust masses and optical depths based on optical data alone, using overly simplistic star-dust geometries and the neglect of scattering effects. [ABRIDGED]

Interstellar dust in galaxies can be traced either through its extinction effects on the star light, or through its thermal emission at infrared wavelengths. Recent radiative transfer studies of several nearby edge-on galaxies have found an apparent inconsistency in the dust energy balance: the radiative transfer models that successfully explain the optical extinction underestimate the observed fluxes by an average factor of three. We investigate the dust energy balance for IC4225 and NGC5166, two edge-on spiral galaxies observed by the Herschel Space Observatory in the frame of the H-ATLAS survey. We start from models which were constrained from optical data and extend them to construct the entire spectral energy distribution of our galaxies. These predicted values are subsequently compared to the observed far-infrared fluxes. We find that including a young stellar population in the modelling is necessary as it plays a non-negligible part in the heating of the dust grains. While the modelling approach for both galaxies is nearly identical, we find two very different results. As is often seen in other edge-on spiral galaxies, the far-infrared emission of our radiative transfer model of IC4225 underestimates the observed fluxes by a factor of about three. For NGC5166 on the other hand, we find that both the predicted spectral energy distribution as well as the simulated images match the observations particularly well. We explore possible reasons for this difference and conclude that it is unlikely that one single mechanism is the cause of the dust energy balance problem in spiral galaxies. We discuss the different approaches that can be considered in order to get a conclusive answer on the origin this discrepancy.

We present the stellar and gaseous kinematics of an Sb galaxy, NGC 3223, with the aim of determining the vertical and radial stellar velocity dispersion as a function of radius, which can help to constrain disk heating theories. Together with the observed NIR photometry, the vertical velocity dispersion is also used to determine the stellar mass-to-light (M/L) ratio, typically one of the largest uncertainties when deriving the dark matter distribution from the observed rotation curve. We find a vertical-to-radial velocity dispersion ratio of sigma_z/sigma_R=1.21+-0.14, significantly higher than expectations from known correlations, and a weakly-constrained Ks-band stellar M/L ratio in the range 0.5-1.7, at the high end of (but consistent with) the predictions of stellar population synthesis models. Such a weak constraint on the stellar M/L ratio, however, does not allow us to securely determine the dark matter density distribution. To achieve this, either a statistical approach or additional data (e.g. integral-field unit) are needed.

The interstellar dust content in galaxies can be traced in extinction at optical wavelengths, or in emission in the far-infrared. Several studies have found that radiative transfer models that successfully explain the optical extinction in edge-on spiral galaxies generally underestimate the observed FIR/submm fluxes by a factor of about three. In order to investigate this so-called dust energy balance problem, we use two Milky Way-like galaxies produced by high-resolution hydrodynamical simulations. We create mock optical edge-on views of these simulated galaxies (using the radiative transfer code SKIRT), and we then fit the parameters of a basic spiral galaxy model to these images (using the fitting code FitSKIRT). The basic model includes smooth axisymmetric distributions along a S\'ersic bulge and exponential disc for the stars, and a second exponential disc for the dust. We find that the dust mass recovered by the fitted models is about three times smaller than the known dust mass of the hydrodynamical input models. This factor is in agreement with previous energy balance studies of real edge-on spiral galaxies. On the other hand, fitting the same basic model to less complex input models (e.g. a smooth exponential disc with a spiral perturbation or with random clumps), does recover the dust mass of the input model almost perfectly. Thus it seems that the complex asymmetries and the inhomogeneous structure of real and hydrodynamically simulated galaxies are a lot more efficient at hiding dust than the rather contrived geometries in typical quasi-analytical models. This effect may help explain the discrepancy between the dust emission predicted by radiative transfer models and the observed emission in energy balance studies for edge-on spiral galaxies.

Context: Dust reprocesses about half of the stellar radiation in galaxies. The thermal re-emission by dust of absorbed energy is considered driven merely by young stars and, consequently, often applied to trace the star formation rate in galaxies. Recent studies have argued that the old stellar population might anticipate a non-negligible fraction of the radiative dust heating. Aims: In this work, we aim to analyze the contribution of young (< 100 Myr) and old (~ 10 Gyr) stellar populations to radiative dust heating processes in the nearby grand-design spiral galaxy M51 using radiative transfer modeling. High-resolution 3D radiative transfer (RT) models are required to describe the complex morphologies of asymmetric spiral arms and clumpy star-forming regions and model the propagation of light through a dusty medium. Methods: In this paper, we present a new technique developed to model the radiative transfer effects in nearby face-on galaxies. We construct a high-resolution 3D radiative transfer model with the Monte-Carlo code SKIRT accounting for the absorption, scattering and non-local thermal equilibrium (NLTE) emission of dust in M51. The 3D distribution of stars is derived from the 2D morphology observed in the IRAC 3.6 {\mu}m, GALEX FUV, H{\alpha} and MIPS 24 {\mu}m wavebands, assuming an exponential vertical distribution with an appropriate scale height. The dust geometry is constrained through the far-ultraviolet (FUV) attenuation, which is derived from the observed total-infrared-to-far-ultraviolet luminosity ratio. The stellar luminosity, star formation rate and dust mass have been scaled to reproduce the observed stellar spectral energy distribution (SED), FUV attenuation and infrared SED. (abridged)

We investigate the amount and spatial distribution of interstellar dust in edge-on spiral galaxies, using detailed radiative transfer modeling of a homogeneous sample of 12 galaxies selected from the CALIFA survey. Our automated fitting routine, FitSKIRT, was first validated against artificial data. This is done by simultaneously reproducing the SDSS $g$-, $r$-, $i$- and $z$-band observations of a toy model in order to combine the information present in the different bands. We show that this combined, oligochromatic fitting, has clear advantages over standard monochromatic fitting especially regarding constraints on the dust properties. We model all galaxies in our sample using a three-component model, consisting of a double exponential disc to describe the stellar and dust discs and using a S\'ersic profile to describe the central bulge. The full model contains 19 free parameters, and we are able to constrain all these parameters to a satisfactory level of accuracy without human intervention or strong boundary conditions. Apart from two galaxies, the entire sample can be accurately reproduced by our model. We find that the dust disc is about 75% more extended but only half as high as the stellar disc. The average face-on optical depth in the V-band is $0.76$ and the spread of $0.60$ within our sample is quite substantial, which indicates that some spiral galaxies are relatively opaque even when seen face-on.

Context. Edge-on spiral galaxies with prominent dust lanes provide us with an excellent opportunity to study the distribution and properties of the dust within them. The HEROES project was set up to observe a sample of seven large edge-on galaxies across various wavelengths for this investigation. Aims. Within this first paper, we present the Herschel observations and perform a qualitative and quantitative analysis on them, and we derive some global properties of the far infrared and submillimetre emission. Methods. We determine horizontal and vertical profiles from the Herschel observations of the galaxies in the sample and describe the morphology. Modified black-body fits to the global fluxes, measured using aperture photometry, result in dust temperatures and dust masses. The latter values are compared to those that are derived from radiative transfer models taken from the literature. Results. On the whole, our Herschel flux measurements agree well with archival values. We find that the exponential horizontal dust distribution model often used in the literature generally provides a good description of the observed horizontal profiles. Three out of the seven galaxies show signatures of extended vertical emission at 100 and 160 {\mu}m at the 5{\sigma} level, but in two of these it is probably due to deviations from an exactly edge-on orientation. Only for NGC 4013, a galaxy in which vertically extended dust has already been detected in optical images, we can detect vertically extended dust, and the derived scaleheight agrees with the value estimated through radiative transfer modelling. Our analysis hints at a correlation between the dust scaleheight and its degree of clumpiness, which we infer from the difference between the dust masses as calculated from modelling of optical data and from fitting the spectral energy distribution of Herschel datapoints.

We present FitSKIRT, a method to efficiently fit radiative transfer models to UV/optical images of dusty galaxies. These images have the advantage that they have better spatial resolution compared to FIR/submm data. FitSKIRT uses the GAlib genetic algorithm library to optimize the output of the SKIRT Monte Carlo radiative transfer code. Genetic algorithms prove to be a valuable tool in handling the multi- dimensional search space as well as the noise induced by the random nature of the Monte Carlo radiative transfer code. FitSKIRT is tested on artificial images of a simulated edge-on spiral galaxy, where we gradually increase the number of fitted parameters. We find that we can recover all model parameters, even if all 11 model parameters are left unconstrained. Finally, we apply the FitSKIRT code to a V-band image of the edge-on spiral galaxy NGC4013. This galaxy has been modeled previously by other authors using different combinations of radiative transfer codes and optimization methods. Given the different models and techniques and the complexity and degeneracies in the parameter space, we find reasonable agreement between the different models. We conclude that the FitSKIRT method allows comparison between different models and geometries in a quantitative manner and minimizes the need of human intervention and biasing. The high level of automation makes it an ideal tool to use on larger sets of observed data.

We combine new dust continuum observations of the edge-on spiral galaxy NGC 4565 in all Herschel/SPIRE (250, 350, 500 micron) wavebands, obtained as part of the Herschel Reference Survey, and a large set of ancillary data (Spitzer, SDSS, GALEX) to analyze its dust energy balance. We fit a radiative transfer model for the stars and dust to the optical maps with the fitting algorithm FitSKIRT. To account for the observed UV and mid-infrared emission, this initial model was supplemented with both obscured and unobscured star-forming regions. Even though these star-forming complexes provide an additional heating source for the dust, the far-infrared/submillimeter emission long wards of 100 micron is underestimated by a factor of 3-4. This inconsistency in the dust energy budget of NGC 4565 suggests that a sizable fraction (two-thirds) of the total dust reservoir (Mdust ~ 2.9e+8 Msun) consists of a clumpy distribution with no associated young stellar sources. The distribution of those dense dust clouds would be in such a way that they remain unresolved in current far-infrared/submillimeter observations and hardly comtribute to the attenuation at optical wavelengths. More than two-thirds of the dust heating in NGC 4565 is powered by the old stellar population, with localized embedded sources supplying the remaining dust heating in NGC 4565. The results from this detailed dust energy balance study in NGC 4565 is consistent with that of similar analyses of other edge-on spirals.