Dimitri A. Gadotti

Structural properties of pseudo-bulges, classical bulges and elliptical galaxies: a sloan digital sky survey perspective

We have performed 2D bulge/bar/disc decompositions using g, r and i-band images of a representative sample of nearly 1000 galaxies from the Sloan Digital Sky Survey. We show that the Petrosian concentration index is a better proxy for bulge-to-total ratio than the global Sersic index. We show that pseudo-bulges can be distinguished from classical bulges as outliers in the Kormendy relation. We provide the structural parameters and distributions of stellar masses of ellipticals, classical bulges, pseudo- bulges, discs and bars, and find that 32 per cent of the total stellar mass in massive galaxies in the local universe is contained in ellipticals, 36 per cent in discs, 25 per cent in classical bulges, 3 per cent in pseudo-bulges and 4 per cent in bars. Pseudo-bulges are currently undergoing intense star formation activity and populate the blue cloud of the colour-magnitude diagram. Most (though not all) classical bulges are quiescent and populate the red sequence of the diagram. Classical bulges follow a correlation between bulge Sersic index and bulge-to-total ratio, while pseudo-bulges do not. In addition, for a fixed bulge-to-total ratio, pseudo-bulges are less concentrated than classical bulges. Pseudo-bulges follow a mass-size relation similar to that followed by bars, and different from that followed by classical bulges. In the fundamental plane, pseudo-bulges occupy the same locus as discs. While these results point out different formation processes for classical and pseudo-bulges, we also find a significant overlap in their properties, indicating that the different processes might happen concomitantly. Finally, classical bulges and ellipticals follow offset mass-size relations, suggesting that high-mass bulges might not be simply high-mass ellipticals surrounded by discs.

RECONSTRUCTING THE STELLAR MASS DISTRIBUTIONS OF GALAXIES USING S4G IRAC 3.6 AND 4.5 μm IMAGES. I. CORRECTING FOR CONTAMINATION BY POLYCYCLIC AROMATIC HYDROCARBONS, HOT DUST, AND INTERMEDIATE-AGE STARS

With the aim of constructing accurate two-dimensional maps of the stellar mass distribution in nearby galaxies from Spitzer Survey of Stellar Structure in Galaxies 3.6 and 4.5 μm images, we report on the separation of the light from old stars from the emission contributed by contaminants. Results for a small sample of six disk galaxies (NGC 1566, NGC 2976, NGC 3031, NGC 3184, NGC 4321, and NGC 5194) with a range of morphological properties, dust content, and star formation histories are presented to demonstrate our approach. To isolate the old stellar light from contaminant emission (e.g., hot dust and the 3.3 μm polycyclic aromatic hydrocarbon (PAH) feature) in the IRAC 3.6 and 4.5 μm bands we use an independent component analysis (ICA) technique designed to separate statistically independent source distributions, maximizing the distinction in the [3.6]-[4.5] colors of the sources. The technique also removes emission from evolved red objects with a low mass-to-light ratio, such as asymptotic giant branch (AGB) and red supergiant (RSG) stars, revealing maps of the underlying old distribution of light with [3.6]-[4.5] colors consistent with the colors of K and M giants. The contaminants are studied by comparison with the non-stellar emission imaged at 8 μm, which is dominated by the broad PAH feature. Using the measured 3.6 μm/8 μm ratio to select individual contaminants, we find that hot dust and PAHs together contribute between ~5% and 15% to the integrated light at 3.6 μm, while light from regions dominated by intermediate-age (AGB and RSG) stars accounts for only 1%-5%. Locally, however, the contribution from either contaminant can reach much higher levels; dust contributes on average 22% to the emission in star-forming regions throughout the sample, while intermediate-age stars contribute upward of 50% in localized knots. The removal of these contaminants with ICA leaves maps of the old stellar disk that retain a high degree of structural information and are ideally suited for tracing stellar mass, as will be the focus in a companion paper.

Image decomposition of barred galaxies and AGN hosts

I present the results of multicomponent decomposition of V and R broad-band images of a sample of 17 nearby galaxies, most of them hosting bars and active galactic nuclei (AGN). I use BUDDA v2.1 to produce the fits, allowing the inclusion of bars and AGN in the models. A comparison with previous results from the literature shows a fairly good agreement. It is found that the axial ratio of bars, as measured from ellipse fits, can be severely underestimated if the galaxy axisymmetric component is relatively luminous. Thus, reliable bar axial ratios can only be determined by taking into account the contributions of bulge and disc to the light distribution in the galaxy image. Through a number of tests, I show that neglecting bars when modelling barred galaxies can result in an overestimation of the bulge-to-total luminosity ratio of a factor of 2. Similar effects result when bright, type 1 AGN are not considered in the models. By artificially redshifting the images, I show that the structural parameters of more distant galaxies can in general be reliably retrieved through image fitting, at least up to the point where the physical spatial resolution is ≈1.5 kpc. This corresponds, for instance, to images of galaxies at z = 0.05 with a seeing full width at half-maximum (FWHM) of 1.5 arcsec, typical of the Sloan Digital Sky Survey (SDSS). In addition, such a resolution is also similar to what can be achieved with the Hubble Space Telescope (HST), and ground-based telescopes with adaptive optics, at z ∼ 1‐2. Thus, these results also concern deeper studies such as COSMOS and SINS. This exercise shows that disc parameters are particularly robust, but bulge parameters are prone to errors if its effective radius is small compared to the seeing radius, and might suffer from systematic effects. For instance, the bulge-to-total luminosity ratio is systematically overestimated, on average, by 0.05 (i.e. 5 per cent of the galaxy total luminosity). In this low-resolution regime, the effects of ignoring bars are still present, but AGN light is smeared out. I briefly discuss the consequences of these results to studies of the structural properties of galaxies, in particular on the stellar mass budget in the local Universe. With reasonable assumptions, it is possible to show that the stellar content in bars can be similar to that in classical bulges and elliptical galaxies. Finally, I revisit the cases of NGC 4608 and 5701 and show that the lack of stars in the disc region inside the bar radius is significant. Accordingly, the best-fitting model for the former uses a Freeman type II disc.

THICK DISKS OF EDGE-ON GALAXIES SEEN THROUGH THE SPITZER SURVEY OF STELLAR STRUCTURE IN GALAXIES (S4G): LAIR OF MISSING BARYONS?

Most, if not all, disk galaxies have a thin (classical) disk and a thick disk. In most models thick disks are thought to be a necessary consequence of the disk formation and/or evolution of the galaxy. We present the results of a study of the thick disk properties in a sample of carefully selected edge-on galaxies with types ranging from T = 3 to T = 8. We fitted one-dimensional luminosity profiles with physically motivated functions—the solutions of two stellar and one gaseous isothermal coupled disks in equilibrium—which are likely to yield more accurate results than other functions used in previous studies. The images used for the fits come from the Spitzer Survey of Stellar Structure in Galaxies (S^4G). We found that thick disks are on average more massive than previously reported, mostly due to the selected fitting function. Typically, the thin and thick disks have similar masses. We also found that thick disks do not flare significantly within the observed range in galactocentric radii and that the ratio of thick-to-thin disk scale heights is higher for galaxies of earlier types. Our results tend to favor an in situ origin for most of the stars in the thick disk. In addition, the thick disk may contain a significant amount of stars coming from satellites accreted after the initial buildup of the galaxy and an extra fraction of stars coming from the secular heating of the thin disk by its own overdensities. Assigning thick disk light to the thin disk component may lead to an underestimate of the overall stellar mass in galaxies because of different mass-to-light ratios in the two disk components. On the basis of our new results, we estimate that disk stellar masses are between 10% and 50% higher than previously thought and we suggest that thick disks are a reservoir of “local missing baryons.”

The Spitzer Survey of Stellar Structure in Galaxies (S4G): Precise Stellar Mass Distributions from Automated Dust Correction at 3.6 μm

The mid-infrared is an optimal window to trace stellar mass in nearby galaxies and the 3.6 m m IRAC band has been exploited to this effect, but such mass estimates can be biased by dust emission. We present our pipeline to reveal the old stellar flux at 3.6 μm and obtain stellar mass maps for more than 1600 galaxies available from the Spitzer Survey of Stellar Structure in Galaxies (SG). This survey consists of images in two infrared bands (3.6 and 4.5 m m ), and we use the Independent Component Analysis (ICA) method presented in Meidt et al. to separate the dominant light from old stars and the dust emission that can significantly contribute to the observed 3.6 m m flux. We exclude from our ICA analysis galaxies with low signal-to-noise ratio (S N 10 < ) and those with original [3.6]–[4.5] colors compatible with an old stellar population, indicative of little dust emission (mostly early Hubble types, which can directly provide good mass maps). For the remaining 1251 galaxies to which ICA was successfully applied, we find that as much as 10%–30% of the total light at 3.6 m m typically originates from dust, and locally it can reach even higher values. This contamination fraction shows a correlation with specific star formation rates, confirming that the dust emission that we detect is related to star formation. Additionally, we have used our large sample of mass estimates to calibrate a relationship of effective mass-to-light ratio (M/L) as a function of observed [3.6]–[4.5] color: M L log( ) = 0.339( 0.057)  ́ ([3.6] [4.5]) 0.336( 0.002)  . Our final pipeline products have been made public through IRSA, providing the astronomical community with an unprecedentedly large set of stellar mass maps ready to use for scientific applications.

An observer's view of simulated galaxies: disc‐to‐total ratios, bars and (pseudo‐)bulges

We use cosmological hydrodynamical simulations of the formation of Milky Way-mass galaxies to study the relative importance of the main stellar components, i.e. discs, bulges and bars, at redshift zero. The main aim of this Letter is to understand if estimates of the structural parameters of these components determined from kinematics (as is usually done in simulations) agree well with those obtained using a photometric bulge/disc/bar decomposition (as done in observations). To perform such a comparison, we have produced synthetic observations of the simulation outputs with the Monte Carlo radiative transfer code sunrise and used the budda code to make 2D photometric decompositions of the resulting images (in the i and g bands). We find that the kinematic disc-to-total (D/T) ratio estimates are systematically and significantly lower than the photometric ones. While the maximum D/T ratios obtained with the former method are of the order of 0.2, they are typically >0.4, and can be as high as 0.7, according to the latter. The photometric decomposition shows that many of the simulated galaxies have bars, with Bar/T ratios in the range 0.2–0.4, and that bulges have in all cases low Sersic indices, resembling observed pseudo-bulges instead of classical ones. Simulated discs, bulges and bars generally have similar g−i colours, which are in the blue tail of the distribution of observed colours. This is not due to the presence of young stars, but rather due to low metallicities and poor gas content in the simulated galaxies, which makes dust extinction low. Photometric decompositions thus match the component ratios usually quoted for spiral galaxies better than kinematic decompositions, but the shift is insufficient to make the simulations consistent with observed late-type systems.

Herschel-ATLAS: The dust energy balance in the edge-on spiral galaxy UGC 4754

We use Herschel PACS and SPIRE observations of the edge-on spiral galaxy UGC4754, taken as part of the H-ATLAS SDP observations, to investigate the dust energy balance in this galaxy. We build detailed SKIRT radiative models based on SDSS and UKIDSS maps and use these models to predict the far-infrared emission. We find that our radiative transfer model underestimates the observed FIR emission by a factor of two to three. Similar discrepancies have been found for other edge-on spiral galaxies based on IRAS, ISO, and SCUBA data. Thanks to the good sampling of the SED at FIR wavelengths, we can rule out an underestimation of the FIR emissivity as the cause for this discrepancy. Instead we support highly obscured star formation that contributes little to the optical extinction as a more probable explanation.

THE IMPACT OF BARS ON DISK BREAKS AS PROBED BY S 4 G IMAGING.

We have analyzed the radial distribution of old stars in a sample of 218 nearby face-on disks, using deep 3.6 mu m images from the Spitzer Survey of Stellar Structure in Galaxies. In particular, we have studied the structural properties of those disks with a broken or down-bending profile. We find that, on average, disks with a genuine single-exponential profile have a scale length and a central surface brightness which are intermediate to those of the inner and outer components of a down-bending disk with the same total stellar mass. In the particular case of barred galaxies, the ratio between the break and the bar radii (R-br/R-bar) depends strongly on the total stellar mass of the galaxy. For galaxies more massive than 10(10) M-circle dot, the distribution is bimodal, peaking at R-br/R-bar similar to 2 and similar to 3.5. The first peak, which is the most populated one, is linked to the outer Lindblad resonance of the bar, whereas the second one is consistent with a dynamical coupling between the bar and the spiral pattern. For galaxies below 10(10) M-circle dot, breaks are found up to similar to 10 R-bar, but we show that they could still be caused by resonances given the rising nature of rotation curves in these low-mass disks. While not ruling out star formation thresholds, our results imply that radial stellar migration induced by non-axisymmetric features can be responsible not only for those breaks at similar to 2 R-bar, but also for many of those found at larger radii.

The spitzer survey of stellar structure in galaxies (S^4G): multi-component decomposition strategies and data release

The Spitzer Survey of Stellar Structure in Galaxies (S^4G) is a deep 3.6 and 4.5 μm imaging survey of 2352 nearby (<40 Mpc) galaxies. We describe the S(^4)G data analysis pipeline 4, which is dedicated to two-dimensional structural surface brightness decompositions of 3.6 μm images, using GALFIT3.0. Besides automatic 1-component Sersic fits, and 2-component Sersic bulge + exponential disk fits, we present human-supervised multi-component decompositions, which include, when judged appropriate, a central point source, bulge, disk, and bar components. Comparison of the fitted parameters indicates that multi-component models are needed to obtain reliable estimates for the bulge Sersic index and bulge-to-total light ratio (B/T), confirming earlier results. Here, we describe the preparations of input data done for decompositions, give examples of our decomposition strategy, and describe the data products released via IRSA and via our web page (www.oulu.fi/astronomy/S4G_PIPELINE4/MAIN). These products include all the input data and decomposition files in electronic form, making it easy to extend the decompositions to suit specific science purposes. We also provide our IDL-based visualization tools (GALFIDL) developed for displaying/running GALFIT-decompositions, as well as our mask editing procedure (MASK_EDIT) used in data preparation. A detailed analysis of the bulge, disk, and bar parameters derived from multi-component decompositions will be published separately.

Secular evolution and structural properties of stellar bars in galaxies

I present results from the modeling of stellar bars in nearly 300 barred galaxies in the local universe through parametric multi-component multi-band image fitting. The surface brightness radial profile of bars is described using a Sersic function, and param- eters such as bar effective radius, ellipticity, boxiness, length and mass, and bar-to-total luminosity and mass ratios, are determined, which is unprecedented for a sample of this size. The properties of bars in galaxies with classical bulges and pseudo-bulges are compared. For a fixed bar-to-total mass ratio, pseudo-bulges are on average signif- icantly less massive than classical bulges, indicating that, if pseudo-bulges are formed through bars, further processes are necessary to build a classical bulge. I find a corre- lation between bar ellipticity and boxiness, and define bar strength as the product of these two quantities. I also find correlations between bar strength and normalised bar size, between the sizes of bars and bulges, and between normalised bar size and bulge- to-total ratio. Bars with different ellipticities follow parallel lines in the latter two correlations. These correlations can arise if, starting off with different normalised sizes and ellipticities, bars grow longer and stronger with dynamical age, as a result of an- gular momentum exchange from the inner to the outer parts of galaxies, in agreement with previous theoretical predictions. As a consequence, bar pattern speeds should become lower with bar dynamical age, and towards galaxies with more prominent bulges.