We present the first nonparametric morphological analysis of a set of spiral galaxies from UV to submm wavelengths. Our study is based on high-quality multi-wavelength imaging for nine well-resolved spiral galaxies from the DustPedia database, combined with nonparametric morphology indicators calculated in a consistent way using the {\tt{StatMorph}} package. We measure the half-light radius, the concentration index, the asymmetry index, the smoothness index, the Gini coefficient and the $M_{20}$ indicator in various wavebands from UV to submm wavelengths, as well as in stellar mass, dust mass and star formation rate maps. We find that the interstellar dust in galaxies is distributed in a more extended, less centrally concentrated, more asymmetric, and more clumpy way than the stars. This is particularly evident when comparing morphological indicators based on the stellar mass and dust mass maps. This should serve as a warning sign against treating the dust in galaxies as a simple smooth component. We argue that the nonparametric galaxy morphology of galaxies from UV to submm wavelengths is an interesting test for cosmological hydrodynamics simulations.

Tilted and warped discs inside tilted dark matter haloes are predicted from numerical and semi-analytical studies. In this paper, we use deep imaging to demonstrate the likely existence of tilted outer structures in real galaxies. We consider two SB0 edge-on galaxies, NGC4469 and NGC4452, which exhibit apparent tilted outer discs with respect to the inner structure. In NGC4469, this structure has a boxy shape, inclined by $\Delta$PA$\approx$3$^{\circ}$ with respect to the inner disc, whereas NGC4452 harbours a discy outer structure with $\Delta$PA$\approx$6$^{\circ}$. In spite of the different shapes, both structures have surface brightness profiles close to exponential and make a large contribution ($\sim30$%) to the total galaxy luminosity. In the case of NGC4452, we propose that its tilted disc likely originates from a former fast tidal encounter (probably with IC3381). For NGC4469, a plausible explanation may also be galaxy harassment, which resulted in a tilted or even a tumbling dark matter halo. A less likely possibility is accretion of gas-rich satellites several Gyr ago. New deep observations may potentially reveal more such galaxies with tilted outer structures, especially in clusters. We also consider galaxies, mentioned in the literature, where a central component (a bar or a bulge) is tilted with respect to the stellar disc. According to our numerical simulations, one of the plausible explanations of such observed "tilts" of the bulge/bar is a projection effect due to a not exactly edge-on orientation of the galaxy coupled with a skew angle of the triaxial bulge/bar.

Dust emission, an important diagnostic of star formation and ISM mass throughout the Universe, can be powered by sources unrelated to ongoing star formation. In the framework of the DustPedia project we have set out to disentangle the radiation of the ongoing star formation from that of the older stellar populations. This is done through detailed, 3D radiative transfer simulations of face-on spiral galaxies. In this particular study, we focus on NGC 1068, which contains an active galactic nucleus (AGN). The effect of diffuse dust heating by AGN (beyond the torus) was so far only investigated for quasars. This additional dust heating source further contaminates the broadband fluxes on which classic galaxy modelling tools rely to derive physical properties. We aim to fit a realistic model to the observations of NGC 1068 and quantify the contribution of the several dust heating sources. Our model is able to reproduce the global spectral energy distribution of the galaxy. It matches the resolved optical and infrared images fairly well, but deviates in the UV and the submm. We find a strong wavelength dependency of AGN contamination to the broadband fluxes. It peaks in the MIR, drops in the FIR, but rises again at submm wavelengths. We quantify the contribution of the dust heating sources in each 3D dust cell and find a median value of 83% for the star formation component. The AGN contribution is measurable at the percentage level in the disc, but quickly increases in the inner few 100 pc, peaking above 90%. This is the first time the phenomenon of an AGN heating the diffuse dust beyond its torus is quantified in a nearby star-forming galaxy. NGC 1068 only contains a weak AGN, meaning this effect can be stronger in galaxies with a more luminous AGN. This could significantly impact the derived star formation rates and ISM masses for such systems.

Interstellar dust absorbs stellar light very efficiently and thus shapes the energetic output of galaxies. Studying the impact of different stellar populations on the dust heating remains hard because it requires decoupling the relative geometry of stars and dust, and involves complex processes as scattering and non-local dust heating. We aim to constrain the relative distribution of dust and stellar populations in the spiral galaxy M81 and create a realistic model of the radiation field that describes the observations. Investigating the dust-starlight interaction on local scales, we want to quantify the contribution of young and old stellar populations to the dust heating. We aim to standardise the setup and model selection of such inverse radiative transfer simulations so this can be used for comparable modelling of other nearby galaxies. We present a semi-automated radiative transfer modelling pipeline that implements the necessary steps such as the geometric model construction and the normalisation of the components through an optimisation routine. We use the Monte Carlo radiative transfer code SKIRT to calculate a self-consistent, panchromatic model of the interstellar radiation field. By looking at different stellar populations independently, we can quantify to what extent different stellar age populations contribute to the dust heating. Our method takes into account the effects of non-local heating. We obtain a realistic 3D radiative transfer model of the face-on galaxy M81. We find that only 50.2\% of the dust heating can be attributed to young stellar populations. We confirm a tight correlation between the specific star formation rate and the heating fraction by young stellar populations, both in sky projection and in 3D, also found for radiative transfer models of M31 and M51. We conclude that... (abridged)

Context: Dust in late-type galaxies in the local Universe is responsible for absorbing approximately one third of the energy emitted by stars. It is often assumed that dust heating is mainly attributable to the absorption of UV and optical photons emitted by the youngest (<= 100 Myr) stars. Consequently, thermal re-emission by dust at FIR wavelengths is often linked to the star-formation activity of a galaxy. However, several studies argue that the contribution to dust heating by much older stars might be more significant. Advances in radiation transfer (RT) simulations finally allow us to actually quantify the heating mechanisms of diffuse dust by the stellar radiation field. Aims: As one of the main goals in the DustPedia project, we have constructed detailed 3D stellar and dust RT models for nearby galaxies. We analyse the contribution of the different stellar populations to the dust heating in four face-on barred galaxies: NGC1365, M83, M95, and M100. We aim to quantify the fraction directly related to young stars, both globally and on local scales, and to assess the influence of the bar on the heating fraction. Results: We derive global attenuation laws for each galaxy and confirm that galaxies of high sSFR have shallower attenuation curves and weaker UV bumps. On average, 36.5% of the bolometric luminosity is absorbed by dust. We report a clear effect of the bar structure on the radial profiles of the dust-heating fraction by the young stars, and the dust temperature. We find that the young stars are the main contributors to the dust heating, donating, on average ~59% of their luminosity to this purpose throughout the galaxy. This dust-heating fraction drops to ~53% in the bar region and ~38% in the bulge region where the old stars are the dominant contributors to the dust heating. We also find a strong link between the heating fraction by the young stars and the sSFR.

In this work we determine the parameters of spiral structure for a sample of face-on spiral galaxies. In practice, the solution of this problem is a hard task because of the diversity of the observed characteristics of spiral structure, such as the arm number, their shape, arm contrast etc. In this work we study spiral structure in galaxies based on an analysis of photometric cuts perpendicular to the arm direction. The method is based on an approximation of these slices with an analytical function and derivation of the parameters of spiral structure (arm width, asymmetry, pitch angle) using the fitted parameters of this approximation. The algorithm has been applied to a sample of 155 galaxies selected from the Sloan Digital Sky Survey in different passbands. In this paper we only consider the results on the arm width: most spirals show an increase of their width with galactocentric distance. Only 14 per cent of galaxies in our sample show an opposite trend or have an almost constant arm width at all radii.

We compare the spectral energy distributions (SEDs) and inferred physical properties for simulated and observed galaxies at low redshift. We exploit UV-submillimetre mock fluxes of ~7000 z=0 galaxies from the EAGLE suite of cosmological simulations, derived using the radiative transfer code SKIRT. We compare these to ~800 observed galaxies in the UV-submillimetre range, from the DustPedia sample of nearby galaxies. To derive global properties, we apply the SED fitting code CIGALE consistently to both data sets, using the same set of ~80 million models. The results of this comparison reveal overall agreement between the simulations and observations, both in the SEDs and in the derived physical properties, with a number of discrepancies. The optical and far-infrared regimes, and the scaling relations based upon the global emission, diffuse dust and stellar mass, show high levels of agreement. However, the mid-infrared fluxes of the EAGLE galaxies are overestimated while the far-UV domain is not attenuated enough, compared to the observations. We attribute these discrepancies to a combination of galaxy population differences between the samples, and limitations in the subgrid treatment of star-forming regions in the EAGLE-SKIRT post-processing recipe. Our findings show the importance of detailed radiative transfer calculations and consistent comparison, and provide suggestions for improved numerical models.

The HERON project is aimed at studying halos and low surface brightness details near galaxies. In this second HERON paper we consider in detail deep imaging (down to surface brightness of ~28 mag/arcsec$^2$ in the r band) for 35 galaxies, viewed edge-on. We confirm a range of low surface brightness features previously described in the literature but also report new ones. We classify the observed outer shapes of the galaxies into three main types (and their prototypes): disc/diamond-like (NGC891), oval (NGC4302), and boxy (NGC3628). We show that the shape of the outer disc in galaxies does not often follow the general 3D model of an exponential disc: 17 galaxies in our sample exhibit oval or even boxy isophotes at the periphery. Also, we show that the less flattened the outer disc is, the more oval or boxy its structure. Many galaxies in our sample have an asymmetric outer structure. We propose that the observed diversity of the galaxy outer shapes is defined by the merger history and its intensity: if no recent multiple minor or single major merging took place, the outer shape is diamond-like or discy. On the contrary, interacting galaxies show oval outer shapes, whereas recent merging appears to transform the outer shape to boxy.

Different spiral generation mechanisms are expected to produce different morphological and kinematic features. In this first paper in a series we carefully study the parameters of spiral structure in 155 face-on spiral galaxies, selected from the Sloan Digital Sky Survey, in the three $gri$ bands. We use a method for deriving a set of parameters of spiral structure, such as the width of the spiral arms, their fraction to the total galaxy luminosity and their colour, which have not been properly studied before. Our method is based on an analysis of a set of photometric cuts perpendicular to the direction of a spiral arm. Based on the results of our study, we compare the main three classes of spirals: grand design, multi-armed and flocculent. We conclude that: i) for the vast majority of galaxies (86\%) we observe an increase of their arm width with galactocentric distance; ii) more luminous spirals in grand design galaxies exhibit smaller variations of the pitch angle with radius than less luminous grand design spirals; iii) grand design galaxies show less difference between the pitch angles of individual arms than multi-armed galaxies. Apart from these distinctive features, all three spiral classes do not differ significantly by their pitch angle, arm width, width asymmetry, and environment. Wavelength dependence is found only for the arm fraction. Therefore, observationally we find no strong difference (except for the view and number of arms) between grand design, multi-armed and flocculent spirals in the sample galaxies.

The purpose of this work is the characterization of the main scaling relations between all the ISM components (dust, atomic/molecular/total gas), gas-phase metallicity, and other galaxy properties, such as Mstar and galaxy morphology, for late-type galaxies in the Local Universe. This study is performed by extracting late-type galaxies from the entire DustPedia sample and by exploiting the large and homogeneous dataset available thanks to the DustPedia project. The sample consists of 436 galaxies with morphological stage from T = 1 to 10, Mstar from 6 x 10^7 to 3 x 10^11 Msun, SFR from 6 x 10^(-4) to 60 Msun/yr, and 12 + log(O/H) from 8 to 9.5. The scaling relations involving the molecular gas are studied by assuming both a constant and a metallicity-dependent XCO. The analysis has been performed by means of the survival analysis technique. We confirm that the dust mass correlates very well with the total gas mass, and find -- for the first time -- that the dust mass correlates better with the atomic gas mass than the molecular one. We characterize important mass ratios such as gas fraction, molecular-to-atomic gas mass ratio, dust-to-total gas mass ratio (DGR), and dust-to-stellar mass ratio. Only the assumption of a metallicity-dependent XCO reproduces the expected decrease of the DGR with increasing morphological stage and decreasing gas-phase metallicity, with a slope of about 1. DGR, gas-phase metallicity, and the dust-to-stellar mass ratio are, for our galaxy sample, directly linked to the galaxy morphology. The molecular-to-atomic gas mass ratio and the DGR show a positive correlation for low molecular gas fractions, but for molecular gas rich galaxies this trend breaks down. This trend has never been found previously, to our knowledge. It provides new constraints for theoretical models of galaxy evolution and a reference for high-redshift studies.

Observations of evolution in the dust-to-metal ratio allow us to constrain the dominant dust processing mechanisms. In this work, we present a study of the dust-to-metal and dust-to-gas ratios in a sub-sample of ~500 DustPedia galaxies. Using literature and MUSE emission line fluxes, we derived gas-phase metallicities (oxygen abundances) for over 10000 individual regions and determine characteristic metallicities for each galaxy. We study how the relative dust, gas, and metal contents of galaxies evolve by using metallicity and gas fraction as proxies for evolutionary state. The global oxygen abundance and nitrogen-to-oxygen ratio are found to increase monotonically as galaxies evolve. Additionally, unevolved galaxies (gas fraction > 60%, metallicity 12 + log(O/H) < 8.2) have dust-to-metal ratios that are about a factor of 2.1 lower (a factor of six lower for galaxies with gas fraction > 80%) than the typical dust-to-metal ratio (Md/MZ ~ 0.214) for more evolved sources. However, for high gas fractions, the scatter is larger due to larger observational uncertainties as well as a potential dependence of the dust grain growth timescale and supernova dust yield on local conditions and star formation histories. We find chemical evolution models with a strong contribution from dust grain growth describe these observations reasonably well. The dust-to-metal ratio is also found to be lower for low stellar masses and high specific star formation rates (with the exception of some sources undergoing a starburst). Finally, the metallicity gradient correlates weakly with the HI-to-stellar mass ratio, the effective radius and the dust-to-stellar mass ratio, but not with stellar mass.

The dust mass absorption coefficient, $\kappa_{d}$, is the conversion function used to infer physical dust masses from observations of dust emission. However, it is notoriously poorly constrained, and it is highly uncertain how it varies, either between or within galaxies. Here we present the results of a proof-of concept study, using the DustPedia data for two nearby face-on spiral galaxies M74 (NGC 628) and M83 (NGC 5236), to create the first ever maps of $\kappa_{d}$ in galaxies. We determine $\kappa_{d}$ using an empirical method that exploits the fact that the dust-to-metals ratio of the interstellar medium is constrained by direct measurements of the depletion of gas-phase metals. We apply this method pixel-by-pixel within M74 and M83, to create maps of $\kappa_{d}$. We also demonstrate a novel method of producing metallicity maps for galaxies with irregularly-sampled measurements, using the machine learning technique of Gaussian process regression. We find strong evidence for significant variation in $\kappa_{d}$. We find values of $\kappa_{d}$ at 500 $\mu$m spanning the range 0.11-0.25 ${\rm m^{2}\,kg^{-1}}$ in M74, and 0.15-0.80 ${\rm m^{2}\,kg^{-1}}$ in M83. Surprisingly, we find that $\kappa_{d}$ shows a distinct inverse correlation with the local density of the interstellar medium. This inverse correlation is the opposite of what is predicted by standard dust models. However, we find this relationship to be robust against a large range of changes to our method - only the adoption of unphysical or highly unusual assumptions would be able to suppress it.

Dust plays an important role in shaping a galaxy's spectral energy distribution (SED). It absorbs ultraviolet (UV) to near-infrared (NIR) radiation and re-emits this energy in the far-infrared (FIR). The FIR is essential to understand dust in galaxies. However, deep FIR observations require a space mission, none of which are still active today. We aim to infer the FIR emission across six Herschel bands, along with dust luminosity, mass, and effective temperature, based on the available UV to mid-infrared (MIR) observations. We also want to estimate the uncertainties of these predictions, compare our method to energy balance SED fitting, and determine possible limitations of the model. We propose a machine learning framework to predict the FIR fluxes from 14 UV-MIR broadband fluxes. We used a low redshift sample by combining DustPedia and H-ATLAS, and extracted Bayesian flux posteriors through SED fitting. We trained shallow neural networks to predict the far-infrared fluxes, uncertainties, and dust properties. We evaluated them on a test set using a root mean square error (RMSE) in log-space. Our results (RMSE = 0.19 dex) significantly outperform UV-MIR energy balance SED fitting (RMSE = 0.38 dex), and are inherently unbiased. We can identify when the predictions are off, for example when the input has large uncertainties on WISE 22, or when the input does not resemble the training set. The galaxies for which we have UV-FIR observations can be used as a blueprint for galaxies that lack FIR data. This results in a 'virtual FIR telescope', which can be applied to large optical-MIR galaxy samples. This helps bridge the gap until the next FIR mission.

[ABRIDGED] We derive the dust properties for 753 local galaxies and examine how these relate to some of their physical properties. We model their global dust-SEDs, treated statistically as an ensemble within a hierarchical Bayesian dust-SED modeling approach. The model-derived properties are the dust masses (Mdust), the average interstellar radiation field intensities (Uav), the mass fraction of very small dust grains ('QPAH' fraction), as well as their standard deviations. In addition, we use mid-IR observations to derive SFR and Mstar, quantities independent of the modeling. We derive distribution functions of the properties for the galaxy ensemble and per galaxy type. The mean value of Mdust for the ETGs is lower than that for the LTGs and IRs, despite ETGs and LTGs having Mstar spanning across the whole range observed. The Uav and 'QPAH' fraction show no difference among different galaxy types. When fixing Uav to the Galactic value, the derived 'QPAH' fraction varies across the Galactic value (0.071). The sSFR increases with galaxy type, while this is not the case for the dust-sSFR (=SFR/Mdust), showing an almost constant SFE per galaxy type. The galaxy sample is characterised by a tight relation between Mdust and Mstar for the LTGs and Irs, while ETGs scatter around this relation and tend towards smaller Mdust. While the relation indicates that Mdust may fundamentally be linked to Mstar, metallicity and Uav are the second parameter driving the scatter, which we investigate in a forthcoming work. We use the extended KS law to estimate Mgas and the GDR. The Mgas derived from the extended KS law is on average ~20% higher than that derived from the KS law, and a large standard deviation indicates the importance of the average SF present to regulate star formation and gas supply. The average GDR for the LTGs and IRs is 370, while including the ETGs gives an average of 550. [ABRIDGED]

Aims: We compare the far-infrared to sub-millimetre dust emission properties measured in high Galactic latitude cirrus with those determined in a sample of 204 late-type DustPedia galaxies. The aim is to verify if it is appropriate to use Milky Way dust properties to derive dust masses in external galaxies. Methods: We used Herschel observations and atomic and molecular gas masses to estimate the disc-averaged dust emissivity at 250 micrometres, and from this, the absorption cross section per H atom and per dust mass. The emissivity requires one assumption, which is the CO-to-H_2 conversion factor, and the dust temperature is additionally required for the absorption cross section per H atom; yet another constraint on the dust-to-hydrogen ratio D/H, depending on metallicity, is required for the absorption cross section dust mass. Results: We find epsilon(250) = 0.82 +/- 0.07 MJy sr^-1 (1E20 H cm^-2)^-1 for galaxies with 4 < F(250)/F(500) < 5. This depends only weakly on the adopted CO-to-H_2 conversion factor. The value is almost the same as that for the Milky Way at the same colour ratio. Instead, for F(250)/F(500) > 6, epsilon(250) is lower than predicted by its dependence on the heating conditions. The reduction suggests a variation in dust emission properties for spirals of earlier type, higher metallicity, and with a higher fraction of molecular gas. When the standard emission properties of Galactic cirrus are used for these galaxies, their dust masses might be underestimated by up to a factor of two. Values for the absorption cross sections at the Milky Way metallicity are also close to those of the cirrus. Mild trends of the absorption cross sections with metallicity are found, although the results depend on the assumptions made.

We use a dedicated 0.7-m telescope to image the halos of 119 galaxies in the Local Volume to $\mu_r \sim 28-30$ mag/arcsec$^2$. The sample is primarily from the 2MASS Large Galaxy Atlas and extended to include nearby dwarf galaxies and more distant giant ellipticals, and spans fully the galaxy colour-magnitude diagram including the blue cloud and red sequence. We present an initial overview, including deep images of our galaxies. Our observations reproduce previously reported low surface brightness structures, including extended plumes in M51, and a newly discovered tidally extended dwarf galaxy in NGC7331. Low surface brightness structures, or "envelopes", exceeding 50 kpc in diameter are found mostly in galaxies with $M_V<-20.5$, and classic interaction signatures are infrequent. Defining a halo diameter at the surface brightness 28 mag/arcsec$^2$, we find that halo diameter is correlated with total galaxy luminosity. Extended signatures of interaction are found throughout the galaxy colour-magnitude diagram without preference for the red or blue sequences, or the green valley. Large envelopes may be found throughout the colour-magnitude diagram with some preference for the bright end of the red sequence. Spiral and S0 galaxies have broadly similar sizes, but ellipticals extend to notably greater diameters, reaching 150 kpc. We propose that the extended envelopes of disk galaxies are dominated by an extension of the disk population rather than by a classical population II halo.

We use a sub-set of the DustPedia galaxy sample (461 galaxies) to investigate the effect the environment has had on galaxies. We consider Virgo cluster and field samples and also assign a density contrast parameter to each galaxy, as defined by the local density of SDSS galaxies. We consider their chemical evolution (using M_{Dust}/M_{Baryon} and M_{Gas}/M_{Baryon}), their specific star formation rate (SFR/M_{Stars}), star formation efficiency (SFR/M_{Gas}), stars-to-dust mass ratio (M_{Stars}/M_{Dust}), gas-to-dust mass ratio (M_{Gas}/M_{Dust}) and the relationship between star formation rate per unit mass of dust and dust temperature (SFR/M_{Dust} and T_{Dust}). Late type galaxies (later than Sc) in all of the environments can be modelled using simple closed box chemical evolution and a simple star formation history (SFR(t) \propto t\exp{-t/\tau}). For earlier type galaxies the physical mechanisms that give rise to their properties are clearly much more varied and require a more complicated model (mergers, gas in or outflow). However, we find little or no difference in the properties of galaxies of the same morphological type within the cluster, field or with different density contrasts. It appears that it is morphology, how and whenever this is laid down, and consistent internal physical processes that primarily determine the derived properties of galaxies in the DustPedia sample and not processes related to differences in the local environment.

Within the framework of the DustPedia project we investigate the properties of cosmic dust and its interaction with the stellar radiation (originating from different stellar populations) for 814 galaxies in the nearby Universe, all observed by the Herschel Space Observatory. We take advantage of the widely used galaxy SED fitting code CIGALE, properly adapted to include the state-of-the-art dust model THEMIS. Using the DustPedia photometry we determine the physical properties of the galaxies, such as, the dust and stellar mass, the star-formation rate, the bolometric luminosity as well as the unattenuated and the absorbed by dust stellar light, for both the old (> 200 Myr) and young (<= 200 Myr) stellar populations. We show how the mass of stars, dust, and atomic gas, as well as the star-formation rate and the dust temperature vary between galaxies of different morphologies and provide recipes to estimate these parameters given their Hubble stage (T). We find a mild correlation between the mass fraction of the small a-C(:H) grains with the specific star-formation rate. On average, young stars are very efficient in heating the dust, with absorption fractions reaching as high as ~77% of the total, unattenuated luminosity of this population. On the other hand, the maximum absorption fraction of old stars is ~24%. Dust heating in early-type galaxies is mainly due to old stars, up to a level of ~90%. Young stars progressively contribute more for `typical' spiral galaxies and they become the dominant source of dust heating for Sm type and irregular galaxies, donating up to ~60% of their luminosity to this purpose. Finally, we find a strong correlation of the dust heating fraction by young stars with morphology and the specific star-formation rate.

We present new isochrone fits to colour-magnitude diagrams of the Galactic globular cluster NGC 5904 (M5). We utilise 29 photometric bands from the ultraviolet to mid-infrared by use of the data from the {\it Hubble Space Telescope}, {\it Gaia} DR2, {\it Wide-field Infrared Survey Explorer}, Sloan Digital Sky Survey (SDSS), and other photometric data. In our isochrone fitting we use the PAdova and TRieste Stellar Evolution Code, the MESA Isochrones and Stellar Tracks, the Dartmouth Stellar Evolution Program, and a Bag of Stellar Tracks and Isochrones both for the solar-scaled and enhanced He and $\alpha$ abundances with a metallicity about [Fe/H]$=-1.33$ adopted from the literature. All tools provide us with estimates of the distance, age, and extinction law to the cluster. The best-fit distance, true distance modulus, and age are $7.4\pm0.3$ kpc, $14.34\pm0.09$ mag, and $12.15\pm1.00$ Gyr, respectively. The derived distance agrees with the literature, including the {\it Gaia} DR2 parallax with its known global zero-point correction. All the data and models, except some UV and SDSS data, agree with the extinction law of Cardelli-Clayton-Mathis with $R_\mathrm{V}=3.60\pm0.05$ and $A_\mathrm{V}=0.20\pm0.02$ mag. This extinction is twice as high as generally accepted due to a rather high extinction between 625 and 2000 nm. An offset of the model colours instead of the high extinction in this range is a less likely, yet possible explanation of the discovered large deviations of the isochrones from the data.

Polar-ring galaxies (PRGs) are fascinating systems in which the central object (typically, an early-type galaxy) is encircled by a large-scale ring of stars, gas, and dust with almost perpendicular spin. PRGs are rare objects and their formation mechanism is not entirely clear. In this paper, we present a new sample of 31 candidates to PRGs identified in the Sloan Digital Sky Survey (SDSS). Using their stacked $gri$ images, we determined geometrical parameters of these galaxies (the position angle and the size of the host galaxy and the ring). We compare our sample objects to the previously known PRGs and discuss their general characteristics (the frequency of faint outer structures, the luminosity--size relation, and the distribution by the apparent angle between the ring and the host galaxy). Our main results are: (i) Central galaxies of PRGs follow the luminosity--size relation for ordinary galaxies. The ring components are located along the similar relation but with a larger scatter. This increasing scatter may be the result of a secondary origin of polar structures. (ii) At least half of PRGs show a ring component within 20$^\circ$ from the perpendicular orientation.

Most radiative transfer models assume that dust in spiral galaxies is distributed exponentially. In this paper our goal is to verify this assumption by analysing the two-dimensional large-scale distribution of dust in galaxies from the DustPedia sample. For this purpose, we make use of Herschel imaging in five bands, from 100 to 500{\mu}m, in which the cold dust constituent is primarily traced and makes up the bulk of the dust mass in spiral galaxies. For a subsample of 320 disc galaxies, we successfully perform a simultaneous fitting with a single S\'ersic model of the Herschel images in all five bands using the multiband modelling code GALFITM. We report that the S\'ersic index $n$, which characterises the shape of the S\'ersic profile, lies systematically below 1 in all Herschel bands and is almost constant with wavelength. The average value at 250{\mu}m is $0.67\pm0.37$ (187 galaxies are fitted with $n_{250}\leq0.75$, 87 galaxies have $0.75<n_{250}\leq1.25$, and 46 - with $n_{250}>1.25$). Most observed profiles exhibit a depletion in the inner region (at $r<0.3-0.4$ of the optical radius $r_{25}$ ) and are more or less exponential in the outer part. We also find breaks in the dust emission profiles at longer distances $(0.5-0.6)r_{25}$ which are associated with the breaks in the optical and near-infrared. We assume that the observed deficit of dust emission in the inner galaxy region is related to the depression in the radial profile of the HI surface density in the same region because the atomic gas reaches high enough surface densities there to be transformed into molecular gas. If a galaxy has a triggered star formation in the inner region (for example, because of a strong bar instability, which transfers the gas inwards to the centre, or a pseudobulge formation), no depletion or even an excess of dust emission in the centre is observed.

We analyse an all-sky compilation of thirteen data sources with optical interstellar linear polarisation for 3871 {\it Gaia} DR2 and {\it Hipparcos} stars within 500 pc, which do not exhibit a considerable intrinsic polarisation. These data are analysed together with five 3D maps and models of the reddening $E(B-V)$. We consider variations of the polarisation degree $P$, position angle $\theta$, and polarisation efficiency $P/E(B-V)$ with Galactic coordinates, distance $R$, and dereddened colour. $P$ shows a maximum at the Gould Belt mid-plane. $P/R$ drops in the Local Bubble by several times defining a boundary of the Bubble at $P=0.1$ per cent. All the data sources of the reddening, except Lallement et al., show a drop of $P/E(B-V)$ in the Bubble. A significant number of stars with a too high $P/E(B-V)$, where the reddening is taken from Arenou et al., Drimmel et al. or Lallement et al., either reject the Serkowski's limit $P/E(B-V)<9.3$ per cent per mag or manifest a systematic underestimation of the reddening by these data sources. The spatial and colour dependent variations of $P$ and $E(B-V)$ outside the Bubble compensate each other, resulting in a nearly constant $P/E(B-V)$. A giant envelope of aligned dust dominates at middle and high latitudes outside the Bubble. The Markkanen's cloud, the North Polar Spur, and some other filaments are parts of this envelope.

We study the fraction of stellar radiation absorbed by dust, f_abs, in 814 galaxies of different morphological types. The targets constitute the vast majority (93%) of the DustPedia sample, including almost all large (optical diameter larger than 1'), nearby (v <= 3000 km/s) galaxies observed with the Herschel Space Observatory. For each object, we model the spectral energy distribution from the ultraviolet to the sub-millimetre using the dedicated, aperture-matched DustPedia photometry and the fitting code CIGALE. The value of f_abs is obtained from the total luminosity emitted by dust and from the bolometric luminosity, which are estimated by the fit. On average, 19% of the stellar radiation is absorbed by dust in DustPedia galaxies. The fraction rises to 25% if only late-type galaxies are considered. The dependence of f_abs on morphology, showing a peak for Sb-Sc galaxies, is weak; it reflects a stronger, yet broad, positive correlation with the bolometric luminosity, which is identified for late-type, disk-dominated, high-specific-star-formation rate, gas-rich objects. We find no variation of f_abs with inclination, at odds with radiative transfer models of edge-on galaxies. These results call for a self-consistent modelling of the evolution of the dust mass and geometry along the build-up of the stellar content. We also provide template spectral energy distributions in bins of morphology and luminosity and study the variation of f_abs with stellar mass and specific star formation rate. We confirm that the local Universe is missing the high f_abs}, luminous and actively star-forming objects necessary to explain the energy budget in observations of the extragalactic background light.

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.

Gaia DR1 Tycho-Gaia Astrometric Solution parallaxes, Tycho-2 photometry and reddening/extinction estimates from nine data sources for 38074 giants within 415 pc from the Sun are used to compare their position in the Hertzsprung-Russell diagram with theoretical estimates, which are based on the PARSEC and MIST isochrones and the TRILEGAL model of the Galaxy with its parameters being widely varied. We conclude that (1) some systematic errors of the reddening/extinction estimates are the main uncertainty in this study; (2) any emission-based 2D reddening map cannot give reliable estimates of reddening within 415 pc due to a complex distribution of dust; (3) if a TRILEGAL's set of the parameters of the Galaxy is reliable and if the solar metallicity is Z<0.021, then the reddening at high Galactic latitudes behind the dust layer is underestimated by all 2D reddening maps based on the dust emission observations of IRAS, COBE, and Planck and by their 3D followers (we also discuss some explanations of this underestimation); 4) the reddening/extinction estimates from recent 3D reddening map by Gontcharov, including the median reddening E(B-V)=0.06 mag at |b|>50 deg, give the best fit of the empirical and theoretical data with each other.

We compare eight sources of reddening and extinction estimates for approximately 60,000 {\it Gaia} DR1 Tycho--Gaia Astrometric Solution (TGAS) main sequence stars younger than 3~Gyr with a relative error of the {\it Gaia} parallax less than 0.1. For the majority of the stars, the best 2D dust emission-based reddening maps show considerable differences between the reddening to infinity and the one calculated to the stellar distance using the barometric law of the dust distribution. This proves that the majority of the TGAS stars are embedded in to the Galactic dust layer and a proper 3D treatment of the reddening/extinction is required to reliably calculate their de-reddened colors and absolute magnitudes. The sources with the 3D estimates of reddening are tested in their ability to put the stars among the PARSEC and MIST theoretical isochrones in the Hertzsprung--Russell diagram based on the precise {\it Gaia}, {\it Tycho-2}, 2MASS and {\it WISE} photometry. Only the reddening/extinction estimates by Arenou et al. (1992), Gontcharov (2012), and Gontcharov (2017), being appropriate for nearby stars within 280~pc, provide both the minimal number of outliers bluer than any reasonable isochrone and the correct number of stars younger than 3~Gyr in agreement with the Besan\c{c}on Galaxy model.

The DustPedia project is capitalising on the legacy of the Herschel Space Observatory, using cutting-edge modelling techniques to study dust in the 875 DustPedia galaxies - representing the vast majority of extended galaxies within 3000 km s$^{-1}$ that were observed by Herschel. This work requires a database of multiwavelength imagery and photometry that greatly exceeds the scope (in terms of wavelength coverage and number of galaxies) of any previous local-Universe survey. We constructed a database containing our own custom Herschel reductions, along with standardised archival observations from GALEX, SDSS, DSS, 2MASS, WISE, Spitzer, and Planck. Using these data, we performed consistent aperture-matched photometry, which we combined with external supplementary photometry from IRAS and Planck. We present our multiwavelength imagery and photometry across 42 UV-microwave bands for the 875 DustPedia galaxies. Our aperture-matched photometry, combined with the external supplementary photometry, represents a total of 21,857 photometric measurements. A typical DustPedia galaxy has multiwavelength photometry spanning 25 bands. We also present the Comprehensive & Adaptable Aperture Photometry Routine (CAAPR), the pipeline we developed to carry out our aperture-matched photometry. CAAPR is designed to produce consistent photometry for the enormous range of galaxy and observation types in our data. In particular, CAAPR is able to determine robust cross-compatible uncertainties, thanks to a novel method for reliably extrapolating the aperture noise for observations that cover a very limited amount of background. Our rich database of imagery and photometry is being made available to the community

The vast majority of Gaia DR1 Tycho-Gaia astrometric solution (TGAS) stars are located within a distance of 1200 pc from the Sun and up to 600 pc within the Galactic plane. Complete usage of their parallaxes is possible only if their reddening and interstellar extinction are known. We obtained those from a comparison of the eight most accurate and complete 3D maps of reddening and extinction in the considered space. We analyzed the statistics of the maps and their differences in the cells of $20\times20\times20$ pc in the whole considered space. The comparison of the maps in the cells outside the dust layer allows us to obtain a formula of the correction for both overestimated and minimal reddenings of the SFD map through the dust half-layer. We point out the need for a further study of the minimal reddening through the dust half-layer and the overestimation of higher reddenings by use of the most precise data. Based on the data of all the maps under consideration we compile two catalogues of the reddenings: for 630,109 cells of $20\times20\times20$ pc in the considered space and for 1,758,723 TGAS stars with the most accurate distances.

We present a detailed photometric study of a sample of 22 edge-on galaxies with clearly visible X-shaped structures. We propose a novel method to derive geometrical parameters of these features, along with the parameters of their host galaxies based on the multi-component photometric decomposition of galactic images. To include the X-shaped structure into our photometric model, we use the IMFIT package, in which we implement a new component describing the X-shaped structure. This method is applied for a sample of galaxies with available SDSS and Spitzer IRAC 3.6 $\mu$m observations. In order to explain our results, we perform realistic $N$-body simulations of a Milky Way-type galaxy and compare the observed and the model X-shaped structures. Our main conclusions are as follows: (1) galaxies with strong X-shaped structures reside in approximately the same local environments as field galaxies; (2) the characteristic size of the X-shaped structures is about 2/3 of the bar size; (3) there is a correlation between the X-shaped structure size and its observed flatness: the larger structures are more flattened; (4) our $N$-body simulations qualitatively confirm the observational results and support the bar-driven scenario for the X-shaped structure formation.

The purpose of this work is the characterization of the radial distribution of dust, stars, gas, and star-formation rate (SFR) in a sub-sample of 18 face-on spiral galaxies extracted from the DustPedia sample. This study is performed by exploiting the multi-wavelength, from UV to sub-mm bands, DustPedia database, in addition to molecular (12CO) and atomic (HI) gas maps and metallicity abundance information available in the literature. We fitted the surface brightness profiles of the tracers of dust and stars, the mass surface density profiles of dust, stars, molecular gas, and total gas, and the SFR surface density profiles with an exponential curve and derived their scale-lengths. We also developed a method to solve for the CO-to-H2 conversion factor (alpha_CO) per galaxy by using dust and gas mass profiles. Although each galaxy has its own peculiar behaviour, we identified a common trend of the exponential scale-lengths vs. wavelength. On average, the scale-lengths normalized to the B-band 25 mag/arcsec^2 radius decrease from UV to 70 micron, from 0.4 to 0.2, and then increase back up to 0.3 at 500 microns. The main result is that, on average, the dust mass surface density scale-length is about 1.8 times the stellar one derived from IRAC data and the 3.6 micron surface brightness, and close to that in the UV. We found a mild dependence of the scale-lengths on the Hubble stage T: the scale-lengths of the Herschel bands and the 3.6 micron scale-length tend to increase from earlier to later types, the scale-length at 70 micron tends to be smaller than that at longer sub-mm wavelength with ratios between longer sub-mm wavelengths and 70 micron that decrease with increasing T. The scale-length ratio of SFR and stars shows a weak increasing trend towards later types.

Recently, a deviation of the Gaia TGAS parallaxes from the asteroseismic ones for giants was found. We show that for parallaxes $\varpi<1.5$ mas it can be explained by a selection effect in favour of bright and luminous giants in the Tycho-2 and TGAS catalogues. Another explanation of this deviation seems to be valid for $\varpi>1.5$ mas based on the best extinction estimates: the deviation may be caused not by a bias of parallax, but by an underestimation of the extinction (and, consequently, an overestimation of the calculated absolute magnitude) in the asteroseismic results. We demonstrate that the reliable estimates of the reddening and extinction (about 0.22 mag of the visual extinction for the Kepler field) better fit both the giants and main sequence stars to the PARSEC, MIST and YaPSI isochrones.

We study the properties of galaxies with very thin discs using a sample of 85 objects whose stellar disc radial-to-vertical scale ratio determined from photometric decomposition, exceeds nine. We present evidences of similarities between the very thin disc galaxies (VTD galaxies) and low surface brightness (LSB) disc galaxies, and conclude that both small and giant LSB galaxies may reveal themselves as VTD, edge-on galaxies. Our VTD galaxies are mostly bulgeless, and those with large radial scale length tend to have redder colors. We performed spectral observations of 22 VTD galaxies with the Dual Imaging Spectrograph on the 3.5m telescope at the Apache Point Observatory. The spectra with good resolution (R ~ 5000) allow us to determine the distance and the ionized gas rotation curve maximum for the galaxies. Our VTD galaxies have low dust content, in contrast to regular disc galaxies. Apparently, VTD galaxies reside in specific cosmological low-density environments and tend to have less connection with filaments. Comparing a toy model that assumes marginally low star formation in galactic discs with obtained gas kinematics data, we conclude that there is a threshold central surface density of about 88 Mo/pc**2, which we observe in the case of very thin, rotationally supported galactic discs.

The European Space Agency has invested heavily in two cornerstones missions; Herschel and Planck. The legacy data from these missions provides us with an unprecedented opportunity to study cosmic dust in galaxies so that we can answer fundamental questions about, for example: the origin of the chemical elements, physical processes in the interstellar medium (ISM), its effect on stellar radiation, its relation to star formation and how this relates to the cosmic far infrared background. In this paper we describe the DustPedia project, which is enabling us to develop tools and computer models that will help us relate observed cosmic dust emission to its physical properties (chemical composition, size distribution, temperature), to its origins (evolved stars, super novae, growth in the ISM) and the processes that destroy it (high energy collisions and shock heated gas). To carry out this research we will combine the Herschel/Planck data with that from other sources of data, providing observations at numerous wavelengths (< 41) across the spectral energy distribution, thus creating the DustPedia database. To maximise our spatial resolution and sensitivity to cosmic dust we limit our analysis to 4231 local galaxies (v < 3000 km/s) selected via their near infrared luminosity (stellar mass). To help us interpret the data we have developed a new physical model for dust (THEMIS), a new Bayesian method of fitting and interpreting spectral energy distributions (HerBIE) and a state-of-the-art Monte Carlo photon tracing radiative transfer model (SKIRT). In this the first of the DustPedia papers we describe the project objectives, data sets used and provide an insight into the new scientific methods we plan to implement.

In this paper, we present results of a photometric and kinematic study for a sample of 13 edge-on spiral galaxies with pronounced integral-shape warps of their stellar discs. The global structure of the galaxies is analyzed on the basis of the Sloan Digital Sky Survey (SDSS) imaging, in the g, r and i passbands. Spectroscopic observations are obtained with the 6-m Special Astrophysical Observatory telescope. In general, galaxies of the sample are typical bright spiral galaxies satisfying the Tully-Fisher relation. Most of the galaxies reside in dense spatial environments and, therefore, tidal encounters are the most probable mechanism for generating their stellar warps. We carried out a detailed analysis of the galaxies and their warps and obtained the following main results: (i) maximum angles of stellar warps in our sample are about 20{\deg}; (ii) warps start, on average, between 2 and 3 exponential scale lengths of a disc; (iii) stronger warps start closer to the center, weak warps start farther; (iv) warps are asymmetric, with the typical degree of asymmetry of about several degrees (warp angle); (v) massive dark halo is likely to preclude the formation of strong and asymmetric warps.

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.

Based on SDSS data, we have selected a sample of nine edge-on spiral galaxies with bulges whose major axes show a high inclination to the disk plane. Such objects are called polar-bulge galaxies. They are similar in their morphology to polar-ring galaxies, but the central objects in them have small size and low luminosity. We have performed a photometric analysis of the galaxies in the g and r bands and determined the main characteristics of their bulges and disks. We show that the disks of such galaxies are typical for the disks of spiral galaxies of late morphological types. The integrated characteristics of their bulges are similar to the parameters of normal bulges. The stellar disks of polar-bulge galaxies often show large-scale warps, which can be explained by their interaction with neighboring galaxies or external accretion from outside.

We analyze the dependence of the stellar disc flatness on the galaxy morphological type using 2D decomposition of galaxies from the reliable subsample of the Edge-on Galaxies in SDSS (EGIS) catalogue. Combining these data with the retrieved models of the edge-on galaxies from the Two Micron All Sky Survey (2MASS) and the Spitzer Survey of Stellar Structure in Galaxies (S$^4$G) catalogue, we make the following conclusions: (1) The disc relative thickness $z_0/h$ in the near- and mid-infrared passbands correlates weakly with morphological type and does not correlate with the bulge-to-total luminosity ratio $B/T$ in all studied bands. (2) Applying an 1D photometric profile analysis overestimates the disc thickness in galaxies with large bulges making an illusion of the relationship between the disc flattening and the ratio $B/T$. (3) In our sample the early-type disc galaxies (S0/a) have both flat and "puffed" discs. The early spirals and intermediate-type galaxies have a large scatter of the disc flatness, which can be caused by the presence of a bar: barred galaxies have thicker stellar discs, on average. On the other hand, the late-type spirals are mostly thin galaxies, whereas irregular galaxies have puffed stellar discs.

We present a catalog of true edge-on disk galaxies automatically selected from the Seventh Data Release (DR7) of the Sloan Digital Sky Survey. A visual inspection of the $g$, $r$ and $i$ images of about 15000 galaxies allowed us to split the initial sample of edge-on galaxy candidates into 4768 (31.8% of the initial sample) genuine edge-on galaxies, 8350 (55.7%) non-edge-ons, and 1865 (12.5%) edge-on galaxies not suitable for simple automatic analysis because these objects show signs of interaction, warps, or nearby bright stars project on it. We added more candidate galaxies from RFGC, EFIGI, RC3, and Galaxy Zoo catalogs found in the SDSS footprints. Our final sample consists of 5747 genuine edge-on galaxies. We estimate the structural parameters of the stellar disks (the stellar disk thickness, radial scale length, and central surface brightness) in the galaxies by analyzing photometric profiles in each of the g, r, and i images. We also perform simplified 3-D modeling of the light distribution in the stellar disks of edge-on galaxies from our sample. Our large sample is intended to be used for studying scaling relations in the stellar disks and bulges and for estimating parameters of the thick disks in different types of galaxies via the image stacking. In this paper we present the sample selection procedure and general description of the sample.

We analyzed several basic correlations between structural parameters of galaxies. The data were taken from various samples in different passbands which are available in the literature. We discuss disc scaling relations as well as some debatable issues concerning the so-called Photometric Plane for bulges and elliptical galaxies in different forms and various versions of the famous Kormendy relation. We show that some of the correlations under discussion are artificial (self-correlations), while others truly reveal some new essential details of the structural properties of galaxies. Our main results are as follows: (1) At present, we can not conclude that faint stellar discs are, on average, more thin than discs in high surface brightness galaxies. The ``central surface brightness -- thickness'' correlation appears only as a consequence of the transparent exponential disc model to describe real galaxy discs. (2) The Photometric Plane appears to have no independent physical sense. Various forms of this plane are merely sophisticated versions of the Kormendy relation or of the self-relation involving the central surface brightness of a bulge/elliptical galaxy and the Sersic index n. (3) The Kormendy relation is a physical correlation presumably reflecting the difference in the origin of bright and faint ellipticals and bulges. We present arguments that involve creating artificial samples to prove our main idea.

A new method for selecting stars in the Galactic bar based on 2MASS infrared photometry in combination with stellar proper motions from the Kharkiv XPM catalogue has been implemented. In accordance with this method, red clump and red giant branch stars are preselected on the color -- magnitude diagram and their photometric distances are calculated. Since the stellar proper motions are indicators of a larger velocity dispersion toward the bar and the spiral arms compared to the stars with circular orbits, applying the constraints on the proper motions of the preselected stars that take into account the Galactic rotation has allowed the background stars to be eliminated. Based on a joint analysis of the velocities of the selected stars and their distribution on the Galactic plane, we have confidently identified the segment of the Galactic bar nearest to the Sun with an orientation of 20$^\circ$--25$^\circ$ with respect to the Galactic center -- Sun direction and a semimajor axis of no more than 3 kpc.

We present a sample of edge-on spiral galaxies both of early and late types.The sample consists of 175 galaxies in the Ks-filter, 169 galaxies in the H-filter and 165 galaxies in the J-filter. Bulge and disc decompositions of each galaxy image, taken from the Two Micron All Sky Survey (2MASS), were performed. We discuss several scaling relations for bulges and discs which indicate a tight link between their formation and evolution. We show that galaxies with bulges fitted by the Sersic index n<2 (pseudobulges) have quite different distributions of their structural parameters than galaxies with n>=2 bulges (classical bulges). First of all, the distribution of the apparent bulge axis ratio q_b for the subsample with n<2 can be attributed to triaxial, nearly prolate bulges, while n>=2 bulges seem to be oblate spheroids with moderate flattening. Secondly, the Photometric Plane of the sample bulges is not flat and has a prominent curvature towards small values of n. Thirdly, despite of the existence of a clear relation between the flattening of stellar discs h/z_0 and the relative mass of a spherical component, the distributions over both parameters are quite different for galaxies possesing bulges and pseudobulges.

A sample of edge-on spiral galaxies aimed at a study of the main structural and photometric parameters of edge-on galaxies both of early and late types is presented. The data were taken from the 2MASS in the J, H and K_s filters. The sample consists of 175 galaxies in the K_s-filter, 169 galaxies in the H-filter and 165 galaxies in the J-filter. We present bulge and disc decompositions of each galaxy image. All galaxies have been modelled with a Sersic bulge and exponential disc with the BUDDA v2.1 package. The main conclusions of our general statistical analysis of the sample are: (1) The distribution of the apparent bulge axis ratio q_b for the subsample with n < 2 can be attributed to triaxial, nearly prolate bulges that are seen from different projections, while n > 2 bulges seem to be oblate spheroids with moderate flattening. (2) For the sample galaxies, the effective radius of the bulge r_{e,b}, the disc scalelength h and the disc scaleheight z_0 are well correlated. However, there is a clear trend for the ratio r_{e,b}/h to increase with n. (3) There is a hint that the fundamental planes of discs, which links only disc parameters and the maximum rotational velocity of gas, are different for galaxies with different bulges. (4) The investigation of the Photometric Plane of sample bulges shows that the plane is not flat and has a prominent curvature towards small values of n. For bulges this fact was not noticed earlier. (5) The clear relation between the flattening of stellar discs h/z_0 and the relative mass of a spherical component, including a dark halo, is confirmed not for bulgeless galaxies but for galaxies with massive bulges. (Abridged)