We present a multi-wavelength analysis of the galaxies in nine clusters selected from the WINGS dataset, examining how galaxy structure varies as a function of wavelength and environment using the state of the art software GALAPAGOSII. We simultaneously fit single S\'ersic functions on three optical (u, B and V) and two near-infrared (J and K) bands thus creating a wavelength-dependent model of each galaxy. We measure the magnitudes, effective radius ($R_{e}$) the S\'ersic index ($n$), axis ratio and position angle in each band. The sample contains 790 cluster members (located close to the cluster center < 0.64 R$_{200}$ and 254 non-member galaxies that we further separate based on their morphology into ellipticals, lenticulars and spirals. We find that the S\'ersic index of all galaxies inside clusters remains nearly constant with wavelength while $R_{e}$ decreases as wavelength increases for all morphological types. We do not observe a significant variation on n and $R_{e}$ as a function of projected local density and distance from the clusters center. Comparing the n and $R_{e}$ of bright cluster galaxies with a subsample of non-member galaxies we find that bright cluster galaxies are more concentrated (display high $n$ values) and are more compact (low $R_{e}$). Moreover, the light profile ($\mathcal{N}$) and size ($\mathcal{R}$) of bright cluster galaxies does not change as a function of wavelength in the same manner as non-member galaxies.

The Southern Photometric Local Universe Survey (S-PLUS) is imaging ~9300 deg^2 of the celestial sphere in twelve optical bands using a dedicated 0.8 m robotic telescope, the T80-South, at the Cerro Tololo Inter-American Observatory, Chile. The telescope is equipped with a 9.2k by 9.2k e2v detector with 10 um pixels, resulting in a field-of-view of 2 deg^2 with a plate scale of 0.55"/pixel. The survey consists of four main subfields, which include two non-contiguous fields at high Galactic latitudes (8000 deg^2 at |b| > 30 deg) and two areas of the Galactic plane and bulge (for an additional 1300 deg^2). S-PLUS uses the Javalambre 12-band magnitude system, which includes the 5 u, g, r, i, z broad-band filters and 7 narrow-band filters centered on prominent stellar spectral features: the Balmer jump/[OII], Ca H+K, H-delta, G-band, Mg b triplet, H-alpha, and the Ca triplet. S-PLUS delivers accurate photometric redshifts (delta_z/(1+z) = 0.02 or better) for galaxies with r < 20 AB mag and redshift < 0.5, thus producing a 3D map of the local Universe over a volume of more than 1 (Gpc/h)^3. The final S-PLUS catalogue will also enable the study of star formation and stellar populations in and around the Milky Way and nearby galaxies, as well as searches for quasars, variable sources, and low-metallicity stars. In this paper we introduce the main characteristics of the survey, illustrated with science verification data highlighting the unique capabilities of S-PLUS. We also present the first public data release of ~336 deg^2 of the Stripe-82 area, which is available at http://datalab.noao.edu/splus.

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.

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 present the results of two-component (disc+bar) and three-component (disc+bar+bulge) multiwavelength 2D photometric decompositions of barred galaxies in five SDSS bands ($ugriz$). This sample of $\sim$3,500 nearby ($z<0.06$) galaxies with strong bars selected from the Galaxy Zoo citizen science project is the largest sample of barred galaxies to be studied using photometric decompositions which include a bar component. With detailed structural analysis we obtain physical quantities such as the bar- and bulge-to-total luminosity ratios, effective radii, S\'ersic indices and colours of the individual components. We observe a clear difference in the colours of the components, the discs being bluer than the bars and bulges. An overwhelming fraction of bulge components have S\'ersic indices consistent with being pseudobulges. By comparing the barred galaxies with a mass-matched and volume-limited sample of unbarred galaxies, we examine the connection between the presence of a large-scale galactic bar and the properties of discs and bulges. We find that the discs of unbarred galaxies are significantly bluer compared to the discs of barred galaxies, while there is no significant difference in the colours of the bulges. We find possible evidence of secular evolution via bars that leads to the build-up of pseudobulges and to the quenching of star formation in the discs. We identify a subsample of unbarred galaxies with an inner lens/oval and find that their properties are similar to barred galaxies, consistent with an evolutionary scenario in which bars dissolve into lenses. This scenario deserves further investigation through both theoretical and observational work.

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.

We provide the basic integrated physical properties of all the galaxies contained in the full CASSIS with available broad-band photometry from UV to 22 $\mu$m. We have collected photometric measurements in 14 wavelengths from available public surveys in order to study the SED of each galaxy in CASSIS, thus constructing a final sample of 1,146 galaxies in the redshift range 0<z<2.5. The SEDs are modelled with the CIGALE code which relies on the energy balance between the absorbed stellar and the dust emission while taking into account the possible contribution due to the presence of an AGN. We split the galaxies in three groups, a low-redshift, a mid-redshift and a high-redshift and find that the vast majority of the Spitzer/IRS galaxies are star-forming and lie on or above the star-forming MS of the corresponding redshift. Moreover, the emission of Spitzer/IRS galaxies with z<0.1 is mostly dominated by star-formation, galaxies in the mid-redshift bin are a mixture of star forming and AGN galaxies, while half of the galaxies with z<0.5 show moderate or high AGN activity. Additionally, using rest-frame NUV-r colour, S\'ersic indices, optical OIII and NII emission lines we explore the nature of these galaxies by investigating their structure as well as their star-formation and AGN activity. Using a colour magnitude diagram we confirm that 97% of the galaxies with redshift smaller than 0.5 have experienced a recent star-formation episode. For a sub-sample of galaxies with available structural information and redshift smaller than 0.3 we find that early-type galaxies are placed below the MS, while late-type galaxies are found on the MS as expected. Finally, for all the galaxies with z<0.5 and available optical spectral line measurements we find that galaxies with high AGN luminosity, as calculated by CIGALE, are most likely to be classified as composite or AGNs by optical spectral lines.

Using single-component fits to SDSS/UKIDSS images of galaxies in the G09 region of the GAMA survey we study radial colour gradients across the galaxy population. We use the multiwavelength information provided by MegaMorph analysis of galaxy light profiles to calculate intrinsic colour gradients, and divide into six subsamples split by overall S\'{e}rsic index ($n$) and galaxy colour. We find a bimodality in the colour gradients of high- and low-$n$ galaxies in all wavebands, which varies with overall galaxy luminosity. Global trends in colour gradients therefore result from combining the contrasting behaviour of a number of different galaxy populations. The ubiquity of strong negative colour gradients supports the picture of inside-out growth through gas accretion for blue, low-$n$ galaxies, and through dry minor mergers for red, high-$n$ galaxies. An exception is the blue high-n population, with properties indicative of dissipative major mergers.

With a large sample of bright, low-redshift galaxies with optical$-$near-IR imaging from the GAMA survey we use bulge-disc decompositions to understand the wavelength-dependent behavior of single-S\'ersic structural measurements. We denote the variation in single-S\'ersic index with wavelength as $\mathcal{N}$, likewise for effective radius we use $\mathcal{R}$. We find that most galaxies with a substantial disc, even those with no discernable bulge, display a high value of $\mathcal{N}$. The increase in S\'ersic index to longer wavelengths is therefore intrinsic to discs, apparently resulting from radial variations in stellar population and/or dust reddening. Similarly, low values of $\mathcal{R}$ ($<$ 1) are found to be ubiquitous, implying an element of universality in galaxy colour gradients. We also study how bulge and disc colour distributions vary with galaxy type. We find that, rather than all bulges being red and all discs being blue in absolute terms, both components become redder for galaxies with redder total colours. We even observe that bulges in bluer galaxies are typically bluer than discs in red galaxies, and that bulges and discs are closer in colour for fainter galaxies. Trends in total colour are therefore not solely due to the colour or flux dominance of the bulge or disc.

We study how the sizes and radial profiles of galaxies vary with wavelength, by fitting Sersic functions simultaneously to imaging in nine optical and near-infrared bands. To quantify the wavelength dependence of effective radius we use the ratio, $\mathcal{R}$, of measurements in two restframe bands. The dependence of Sersic index on wavelength, $\mathcal{N}$, is computed correspondingly. Vulcani et al. (2014) have demonstrated that different galaxy populations present sharply contrasting behaviour in terms of $\mathcal{R}$ and $\mathcal{N}$. Here we study the luminosity dependence of this result. We find that at higher luminosities, early-type galaxies display a more substantial decrease in effective radius with wavelength, whereas late-types present a more pronounced increase in Sersic index. The structural contrast between types thus increases with luminosity. By considering samples at different redshifts, we demonstrate that lower data quality reduces the apparent difference between the main galaxy populations. However, our conclusions remain robust to this effect. We show that accounting for different redshift and luminosity selections partly reconciles the size variation measured by Vulcani et al. with the weaker trends found by other recent studies. Dividing galaxies by visual morphology confirms the behaviour inferred using morphological proxies, although the sample size is greatly reduced. Finally, we demonstrate that varying dust opacity and disc inclination can account for features of the joint distribution of $\mathcal{R}$ and $\mathcal{N}$ for late-type galaxies. However, dust does not appear to explain the highest values of $\mathcal{R}$ and $\mathcal{N}$. The bulge-disc nature of galaxies must also contribute to the wavelength-dependence of their structure.

Aims. This work investigates the potential of using the wavelength-dependence of galaxy structural parameters (S\'ersic index, n, and effective radius, Re) to separate galaxies into distinct types. Methods. A sample of nearby galaxies with reliable visual morphologies is considered, for which we measure structural parameters by fitting multi-wavelength single-S\'ersic models. Additionally, we use a set of artificially redshifted galaxies to test how these classifiers behave when the signal-to-noise decreases. Results. We show that the wavelength-dependence of n may be employed to separate visually-classified early- and late-type galaxies, in a manner similar to the use of colour and n. Furthermore, we find that the wavelength variation of n can recover galaxies that are misclassified by these other morphological proxies. Roughly half of the spiral galaxies that contaminate an early-type sample selected using (u-r) versus n can be correctly identified as late-types by N, the ratio of n measured in two different bands. Using a set of artificially-redshifted images, we show that this technique remains effective up to z ~ 0.1. N can therefore be used to achieve purer samples of early-types and more complete samples of late-types than using a colour-n cut alone. We also study the suitability of R, the ratio of Re in two different bands, as a morphological classifier, but find that the average sizes of both early- and late-type galaxies do not change substantially over optical wavelengths.

Bulge-disc decomposition is a valuable tool for understanding galaxies. However, achieving robust measurements of component properties is difficult, even with high quality imaging, and it becomes even more so with the imaging typical of large surveys. In this paper we consider the advantages of a new, multi-band approach to galaxy fitting. We perform automated bulge-disc decompositions for 163 nearby galaxies, by simultaneously fitting multiple images taken in five photometric filters. We show that we are able to recover structural measurements that agree well with various other works, and confirm a number of key results. We additionally use our results to illustrate the link between total S\'ersic index and bulge-disc structure, and demonstrate that the visually classification of lenticular galaxies is strongly dependent on the inclination of their disc component. By simulating the same set of galaxies as they would appear if observed at a range of redshifts, we are able to study the behaviour of bulge-disc decompositions as data quality diminishes. We examine how our multi-band fits perform, and compare to the results of more conventional, single-band methods. Multi-band fitting improves the measurement of all parameters, but particularly the bulge-to total flux ratio and component colours. We therefore encourage the use of this approach with future surveys.

We investigate the relationship between colour and structure within galaxies using a large, volume-limited sample of bright, low-redshift galaxies with optical to near-infrared imaging from the GAMA survey. We fit single-component, wavelength-dependent, elliptical Sersic models to all passbands simultaneously, using software developed by the MegaMorph project. Dividing our sample by Sersic index and colour, the recovered wavelength variations in effective radius (R_e) and Sersic index (n) reveal the internal structure, and hence formation history, of different types of galaxies. All these trends depend on n; some have an additional dependence on galaxy colour. Late-type galaxies (n_r < 2.5) show a dramatic increase in Sersic index with wavelength. This might be a result of their two-component (bulge-disk) nature, though stellar population gradients within each component and dust attenuation are likely to play a role. All galaxies show a substantial decrease in R_e with wavelength. This is strongest for early-types (n_r > 2.5), even though they maintain constant n with wavelength, revealing that ellipticals are a superimposition of different stellar populations associated with multiple collapse and merging events. Processes leading to structures with larger R_e must be associated with lower metallicity or younger stellar populations. This appears to rule out the formation of young cores through dissipative gas accretion as an important mechanism in the recent lives of luminous elliptical galaxies.

We demonstrate a new multi-wavelength technique for two-dimensional parametric modelling of galaxy surface-brightness profiles, which we have incorporated into the widely used software GALFIT. Our new method, named GALFITM, extends GALFIT3's current single-band fitting process by simultaneously using multiple images of the same galaxy to constrain a wavelength-dependent model. Each standard profile parameter may vary as a function of wavelength, with a user-definable degree of smoothness, from constant to fully free. The performance of GALFITM is evaluated by fitting elliptical S\'ersic profiles to ugriz imaging data for 4026 galaxies, comprising the original SDSS imaging for 163 low redshift (v < 7000 kms) galaxies and 3863 artificially redshifted (0.01 < z < 0.25) images of the same galaxies. Comparing results from single-band and multi-band techniques, we show that GALFITM significantly improves the extraction of information, particularly from bands with low signal-to-noise ratio (e.g., u and z SDSS bands) when combined with higher signal-to-noise images. We also study systematic trends in the recovered parameters, particularly S\'ersic index, that appear when one performs measurements of the same galaxies at successively higher redshifts. We argue that it is vital that studies investigating the evolution of galaxy structure are careful to avoid or correct for these biases. The resulting multi-band photometric structural parameters for our sample of 163 galaxies are provided. We demonstrate the importance of considering multi-band measurements by showing that the S\'ersic indices of spiral galaxies increase to redder wavelengths, as expected for composite bulge-disk systems. Finally, for the ellipticals in our sample, which should be well-represented by single-S\'ersic models, we compare our measured parameters to those from previous studies.

Scaling relations between supermassive black hole mass, M_BH, and host galaxy properties are a powerful instrument for studying their coevolution. A complete picture involving all of the black hole scaling relations, in which each relation is consistent with the others, is necessary to fully understand the black hole-galaxy connection. The relation between M_BH and the central light concentration of the surrounding bulge, quantified by the S\'ersic index n, may be one of the simplest and strongest such relations, requiring only uncalibrated galaxy images. We have conducted a census of literature S\'ersic index measurements for a sample of 54 local galaxies with directly measured M_BH values. We find a clear M_BH - n relation, despite an appreciable level of scatter due to the heterogeneity of the data. Given the current M_BH - L_sph and the L_sph - n relations, we have additionally derived the expected M_BH - n relations, which are marginally consistent at the 2 sigma level with the observed relations. Elliptical galaxies and the bulges of disc galaxies are each expected to follow two distinct bent M_BH - n relations due to the S\'ersic/core-S\'ersic divide. For the same central light concentration, we predict that M_BH in the S\'ersic bulges of disc galaxies are an order magnitude higher than in S\'ersic elliptical galaxies if they follow the same M_BH - L_sph relation.

In this paper, we demonstrate a new method for fitting galaxy profiles which makes use of the full multi-wavelength data provided by modern large optical-near-infrared imaging surveys. We present a new version of GALAPAGOS, which utilises a recently-developed multi-wavelength version of GALFIT, and enables the automated measurement of wavelength dependent S\'ersic profile parameters for very large samples of galaxies. Our new technique is extensively tested to assess the reliability of both pieces of software, GALFIT and GALAPAGOS on both real ugrizY JHK imaging data from the GAMA survey and simulated data made to the same specifications. We find that fitting galaxy light profiles with multi-wavelength data increases the stability and accuracy of the measured parameters, and hence produces more complete and meaningful multi-wavelength photometry than has been available previously. The improvement is particularly significant for magnitudes in low S/N bands and for structural parameters like half-light radius re and S\'ersic index n for which a prior is used by constraining these parameters to a polynomial as a function of wavelength. This allows the fitting routines to push the magnitude of galaxies for which sensible values can be derived to fainter limits. The technique utilises a smooth transition of galaxy parameters with wavelength, creating more physically meaningful transitions than single-band fitting and allows accurate interpolation between passbands, perfect for derivation of rest-frame values.

We present near-IR surface photometry (2D-profiling) for a sample of 29 nearby galaxies for which super-massive black hole (SMBH) masses are constrained. The data is derived from the UKIDSS-LASS survey representing a significant improvement in image quality and depth over previous studies based on 2MASS data. We derive the spheroid luminosity and spheroid S\'ersic index for each galaxy with GALFIT3 and use these data to construct SMBH mass -bulge luminosity ($M_{\rm bh}$--$L$) and SMBH - S\'ersic index ($M_{\rm bh}$--$n$) relations. The best fit K-band relation for elliptical and disk galaxies is $\log(M_{\rm bh}/M_{\odot})= -0.36(\pm 0.03) (M_{\rm K} + 18) + 6.17(\pm 0.16)$ with an intrinsic scatter of 0.4$^{+0.09}_{-0.06}$dex whilst for elliptical galaxies we find $\log(M_{\rm bh}/M_{\odot})= -0.42(\pm 0.06) (M_{\rm K} + 22) + 7.5(\pm 0.15)$ with an intrinsic scatter of 0.31$^{+0.087}_{-0.047}$dex. Our revised $M_{\rm bh}$--$L$ relation agrees closely with the previous near-IR constraint by \citet{tex:G07}. The lack of improvement in the intrinsic scatter in moving to higher quality near-IR data suggests that the SMBH relations are not currently limited by the quality of the imaging data but is either intrinsic or a result of uncertainty in the precise number of required components required in the profiling process. Contrary to expectation (see \citealt{tex:GD07a}) a relation between SMBH mass and the S\'ersic index was not found at near-IR wavelengths. This latter outcome is believed to be explained by the generic inconsistencies between 1D and 2D galaxy profiling which are currently under further investigation.

Supermassive black hole mass estimates are derived for 1743 galaxies from the Millennium Galaxy Catalogue using the recently revised empirical relation between supermassive black hole mass and the luminosity of the host spheroid. The MGC spheroid luminosities are based on $R^{1/n}$-bulge plus exponential-disc decompositions. The majority of black hole masses reside between $10^6 M_{\odot}$ and an upper limit of $2\times10^9 M_{\odot}$. Using previously determined space density weights, we derive the SMBH mass function which we fit with a Schechter-like function. Integrating the black hole mass function over $10^6< M_{bh}/ M_{\odot} < 10^{10}$ gives a supermassive black hole mass density of ($3.8 \pm 0.6) \times 10^5 h^{3}_{70} M_{\odot}$ Mpc$^{-3}$ for early-type galaxies and ($0.96 \pm 0.2) \times10^5 h^{3}_{70} M_{\odot}$ Mpc$^{-3}$ for late-type galaxies. The errors are estimated from Monte Carlo simulations which include the uncertainties in the $M_{bh}$--$L$ relation, the luminosity of the host spheroid and the intrinsic scatter of the $M_{bh}$--$L$ relation. Assuming supermassive black holes form via baryonic accretion we find that ($0.008\pm0.002) h_{70}^{3}$ per cent of the Universe's baryons are currently locked up in supermassive black holes. This result is consistent with our previous estimate based on the $M_{bh}$--$n$ (S{\'e}rsic index) relation.