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.

Recent cosmological hydrodynamical simulations suggest that baryonic processes, and in particular supernova feedback after bursts of star formation, can alter the structure of dark matter haloes and transform primordial cusps into shallower cores. To assess whether this mechanism offers a solution to the cusp-core controversy, simulated haloes must be compared to real dark matter haloes inferred from galaxy rotation curves. For this purpose, two new dark matter density profiles were recently derived from simulations of galaxies in complementary mass ranges: the DC14 halo ($10^{10} < M_{\text{halo}}/M_{\odot} < 8 \times 10^{11}$) and the coreNFW halo ($10^{7} < M_{\text{halo}}/M_{\odot} < 10^{9}$). Both models have individually been found to give good fits to observed rotation curves. For the DC14 model, however, the agreement of the predicted halo properties with cosmological scaling relations was confirmed by one study, but strongly refuted by another. A next question is whether the two models converge to the same solution in the mass range where both should be appropriate. To investigate this, we tested the DC14 and cNFW halo models on the rotation curves of a selection of galaxies with halo masses in the range $4 \times 10^{9}$ - $7 \times 10^{10}$ $M_{\odot}$. We further applied the DC14 model to a set of rotation curves at higher halo masses, up to $9 \times 10^{11}$ $M_{\odot}$, to verify the agreement with the cosmological scaling relations. We find that both models are generally able to reproduce the observed rotation curves, in line with earlier results, and the predicted dark matter haloes are consistent with the cosmological $c-M_{\text{halo}}$ and $M_{*}-M_{\text{halo}}$ relations. The DC14 and cNFW models are also in fairly good agreement with each other, even though DC14 tends to predict slightly less extended cores and somewhat more concentrated haloes than cNFW.

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)

Evidence is mounting that a significant fraction of the early-type galaxy population contains substantial reservoirs of cold interstellar gas and dust. We investigate the gas and dust in NGC 5485, an early-type galaxy with a prominent minor-axis dust lane. Using new Herschel PACS and SPIRE imaging data, we detect 3.8 x 10^6 Msun of cool interstellar dust in NGC 5485, which is in stark contrast with the non-detection of the galaxy in sensitive HI and CO observations from the ATLAS3D consortium. The resulting gas-to-dust ratio upper limit is Mgas/Md < 14.5, almost an order of magnitude lower than the canonical value for the Milky Way. We scrutinize the reliability of the dust, atomic gas and molecular gas mass estimates, but these do not show systematic uncertainties that can explain the extreme gas-to-dust ratio. Also a warm or hot ionized gas medium does not offer an explanation. A possible scenario could be that NGC 5485 merged with an SMC-type metal-poor galaxy with a substantial CO-dark molecular gas component and that the bulk of atomic gas was lost during the interaction, but it remains to be investigated whether such a scenario is possible.

We investigate the connection between dust and gas in the nearby edge-on spiral galaxy NGC 891. High resolution Herschel PACS and SPIRE 70, 100, 160, 250, 350, and 500 $\mu$m images are combined with JCMT SCUBA 850 $\mu$m observations to trace the far-infrared/submillimetre spectral energy distribution (SED). Maps of the HI 21 cm line and CO(J=3-2) emission trace the atomic and molecular hydrogen gas, respectively. We fit one-component modified blackbody models to the integrated SED, finding a global dust mass of 8.5$\times$10$^{7}$ M$_{\odot}$ and an average temperature of 23$\pm$2 K. We also fit the pixel-by-pixel SEDs to produce maps of the dust mass and temperature. The dust mass distribution correlates with the total stellar population as traced by the 3.6 $\mu$m emission. The derived dust temperature, which ranges from approximately 17 to 24 K, is found to correlate with the 24 $\mu$m emission. Allowing the dust emissivity index to vary, we find an average value of $\beta$ = 1.9$\pm$0.3. We confirm an inverse relation between the dust emissivity spectral index and dust temperature, but do not observe any variation of this relationship with vertical height from the mid-plane of the disk. A comparison of the dust properties with the gaseous components of the ISM reveals strong spatial correlations between the surface mass densities of dust and the molecular hydrogen and total gas surface densities. Observed asymmetries in the dust temperature, and the H$_{2}$-to-dust and total gas-to-dust ratios hint that an enhancement in the star formation rate may be the result of larger quantities of molecular gas available to fuel star formation in the NE compared to the SW. Whilst the asymmetry likely arises from dust obscuration due to the geometry of the line-of-sight projection of the spiral arms, we cannot exclude an enhancement in the star formation rate in the NE side of the disk.

We cross-correlate the Planck Catalogue of Compact Sources (PCCS) with the fully sampled 84 deg2 Herschel Virgo Cluster Survey (HeViCS) fields. We search for and identify the 857 and 545 GHz PCCS sources in the HeViCS fields by studying their FIR/submm and optical counterparts. We find 84 and 48 compact Planck sources in the HeViCS fields at 857 and 545 GHz, respectively. Almost all sources correspond to individual bright Virgo Cluster galaxies. The vast majority of the Planck detected galaxies are late-type spirals, with the Sc class dominating the numbers, while early-type galaxies are virtually absent from the sample, especially at 545 GHz. We compare the HeViCS SPIRE flux densities for the detected galaxies with the four different PCCS flux density estimators and find an excellent correlation with the aperture photometry flux densities, even at the highest flux density levels. We find only seven PCCS sources in the HeViCS fields without a nearby galaxy as obvious counterpart, and conclude that all of these are dominated by Galactic cirrus features or are spurious detections. No Planck sources in the HeViCS fields seem to be associated to high-redshift proto-clusters of dusty galaxies or strongly lensed submm sources. Finally, our study is the first empirical confirmation of the simulation-based estimated completeness of the PCCS, and provides a strong support of the internal PCCS validation procedure.

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.