Using SPITZER 3.6$\mu$m imaging, we investigate the physical and data-driven origins of up-bending (Type III) disk breaks. We apply a robust new break-finding algorithm to 175 low-inclination disk galaxies previously identified as containing Type III breaks, classify each galaxy by its outermost re-classified (via our new algorithm) break type, and compare the local environments of each resulting subgroup. Using three different measures of the local density of galaxies, we find that galaxies with extended outer spheroids (Type IIIs) occupy the highest density environments in our sample, while those with extended down-bending (Type II) disks and symmetric outskirts occupy the lowest density environments. Among outermost breaks, the most common origin of Type III breaks in our sample is methodological; the use of elliptical apertures to measure the radial profiles of asymmetric galaxies usually results in features akin to Type III breaks.

Aims: Using a sample of 175 low-inclination galaxies from the S$^{4}$G, we investigate the origins of up-bending (Type III) breaks in the 3.6 $\mu$m surface brightness profiles of disk galaxies. Methods: We re-analyze a sample of previously identified Type III disk break-hosting galaxies using a new, unbiased break-finding algorithm, which uncovered many new, sometimes subtle disk breaks across the whole sample. We classify each break by its likely origin through close examination of the galaxy images across wavelengths, and compare samples of galaxies separated by their outermost identified break types in terms of their stellar populations and local environments. Results: We find that more than half of the confirmed Type III breaks in our sample can be attributed to morphological asymmetry in the host galaxies. As these breaks are mostly an artifact of the azimuthal averaging process, their status as physical "breaks" is questionable. Such galaxies occupy some of the highest density environments in our sample, implying that much of this asymmetry is the result of tidal disturbance. Additionally, we find that Type III breaks related to extended spiral arms or star formation often host down-bending (Type II) breaks at larger radius which were previously unidentified. Such galaxies reside in the lowest density environments in our sample, in line with previous studies that found a lack of Type II breaks in clusters. Galaxies occupying the highest density environments most often show Type III breaks associated with outer spheroidal components. Conclusions: We find that Type III breaks in the outer disks of galaxies arise most often through environmental influence: either tidal disturbance (resulting in disk asymmetry) or heating through, e.g., galaxy harrassment (leading to spheroidal components). Galaxies hosting the latter break types also show... (abstract continues)

We investigate barlenses in the Calar Alto Legacy Integral Field Area (CALIFA) survey galaxies, studying their morphologies, stellar populations and metallicities. Multi-component decompositions are made using the Sloan Digital Sky Survey (SDSS) images and making GALFIT models, fitting besides bulges, disks and bars, also barlenses, which are the face-on counterparts of Boxy/Peanut bulges. Similar GALFIT models as were made for the galaxies, were made also for the simulation snap-shots. For the stellar populations and metallicities in the various structure components CALIFA IFU-observations are used. We show that, when present, barlenses account for a significant portion of light of photometric bulges (i.e. the excess light on top of the disks), which highlights the importance of bars in accumulating the central galaxy mass concentrations in the cosmic timescale. Bars and barlenses are found to have similar cumulative stellar age and metallicity distributions. This is the first time that a combined morphological and stellar population analysis is used to study barlenses. We show that their stars are accumulated in a prolonged time period, concurrently with the evolution of the narrow bar.

Using a sample of 472 local Universe (z<0.06) galaxies in the stellar mass range 10.25 < log M*/M_sun < 10.75, we explore the variation in galaxy structure as a function of morphology and galaxy colour. Our sample of galaxies is sub-divided into red, green and blue colour groups and into elliptical and non-elliptical (disk-type) morphologies. Using KiDS and VIKING derived postage stamp images, a group of eight volunteers visually classified bars, rings, morphological lenses, tidal streams, shells and signs of merger activity for all systems. We find a significant surplus of rings ($2.3\sigma$) and lenses ($2.9\sigma$) in disk-type galaxies as they transition across the green valley. Combined, this implies a joint ring/lens green valley surplus significance of $3.3\sigma$ relative to equivalent disk-types within either the blue cloud or the red sequence. We recover a bar fraction of ~44% which remains flat with colour, however, we find that the presence of a bar acts to modulate the incidence of rings and (to a lesser extent) lenses, with rings in barred disk-type galaxies more common by ~20-30 percentage points relative to their unbarred counterparts, regardless of colour. Additionally, green valley disk-type galaxies with a bar exhibit a significant $3.0\sigma$ surplus of lenses relative to their blue/red analogues. The existence of such structures rules out violent transformative events as the primary end-of-life evolutionary mechanism, with a more passive scenario the favoured candidate for the majority of galaxies rapidly transitioning across the green valley.

Breaks in the surface brightness profiles in the outer regions of galactic discs are thought to have formed by various internal and external processes, and by studying the breaks we aim to better understand what processes are responsible for the evolution of the outer discs. We use a large well-defined sample to study how common the breaks are, and whether their properties depend on galaxy stellar mass or observed wavelength. We study radial surface brightness profiles of 753 galaxies, obtained from the $3.6 \mu m$ images of the Spitzer Survey of Stellar Structure in Galaxies (S$^4$G), and the $K_s$-band data from the Near InfraRed S0-Sa galaxy Survey (NIRS0S), covering a wide range of galaxy morphologies and stellar masses. Optical SDSS or Liverpool telescope data was used for 480 of these galaxies. We find that in low-mass galaxies the single exponential discs (Type I) are most common, and that their fraction decreases with increasing galaxy stellar mass. The fraction of down-bending (Type II) discs increases with stellar mass, possibly due to more common occurrence of bar resonance structures. The up-bending (Type III) discs are also more common in massive galaxies. The observed wavelength affects the scalelength of the disc of every profile type. Especially the scalelength of the inner disc of Type II profiles increases from infrared to u-band on average by a factor of $\sim 2.2$. Consistent with the previous studies, we find that Type II outer disc scalelengths ($h_o$) in late-type and low-mass galaxies are shorter in bluer wavelengths, possibly due to stellar radial migration populating the outer discs with old stars. In Type III discs $h_o$ are larger in the u-band, hinting to the presence of young stellar population in the outer disc. While the observed wavelength affects the disc parameters, it does not significantly affect the profile type in our sample. (Abridged)

We study the surface brightness profiles of disc galaxies in the 3.6 micron images from the Spitzer Survey of Stellar Structure in Galaxies (S4G) and K_s-band images from the Near Infrared S0-Sa galaxy Survey (NIRS0S). We particularly connect properties of single exponential (type I), downbending double exponential (type II), and upbending double exponential (type III) disc profile types, to structural components of galaxies by using detailed morphological classifications, and size measurements of rings and lenses. We also study how the local environment of the galaxies affects the profile types by calculating parameters describing the environmental density and the tidal interaction strength. We find that in majority of type II profiles the break radius is connected with structural components such as rings, lenses, and spirals. The exponential disc sections of all three profile types, when considered separately, follow the disc scaling relations. However, the outer discs of type II, and the inner discs of type III, are similar in scalelength to the single exponential discs. Although the different profile types have similar mean environmental parameters, the scalelengths of the type III profiles show a positive correlation with the tidal interaction strength.

We analyse structural decompositions of 500 late-type galaxies (Hubble $T$-type $\ge 6$) from the Spitzer Survey of Stellar Structure in Galaxies (S$^4$G), spanning a stellar mass range of about $10^7$ to a few times $10^{10}$ M$_\odot$. Their decomposition parameters are compared with those of the early-type dwarfs in the Virgo cluster from Janz et al. They have morphological similarities, including the fact that the fraction of simple one-component galaxies in both samples increases towards lower galaxy masses. We find that in the late-type two-component galaxies both the inner and outer structures are by a factor of two larger than those in the early-type dwarfs, for the same stellar mass of the component. While dividing the late-type galaxies to low and high density environmental bins, it is noticeable that both the inner and outer components of late types in the high local galaxy density bin are smaller, and lie closer in size to those of the early-type dwarfs. This suggests that, although structural differences between the late and early-type dwarfs are observed, environmental processes can plausibly transform their sizes sufficiently, thus linking them evolutionarily.

Using 3.6 and 4.5$\mu$m images of 73 late-type, edge-on galaxies from the S$^4$G survey, we compare the richness of the globular cluster populations of these galaxies to those of early type galaxies that we measured previously. In general, the galaxies presented here fill in the distribution for galaxies with lower stellar mass, M$_*$, specifically $\log({\rm M}_*/{\rm M}_\odot) < 10$, overlap the results for early-type galaxies of similar masses, and, by doing so, strengthen the case for a dependence of the number of globular clusters per $10^9\ {\rm M}_\odot$ of galaxy stellar mass, T$_{\rm N}$, on M$_*$. For $8.5 < \log ({\rm M}_*/{\rm M}_\odot) < 10.5$ we find the relationship can be satisfactorily described as T$_{\rm N} = ({\rm M}_*/10^{6.7})^{-0.56}$ when M$_*$ is expressed in solar masses. The functional form of the relationship is only weakly constrained and extrapolation outside this range is not advised. Our late-type galaxies, in contrast to our early-types, do not show the tendency for low mass galaxies to split into two T$_{\rm N}$ families. Using these results and a galaxy stellar mass function from the literature, we calculate that in a volume limited, local Universe sample, clusters are most likely to be found around fairly massive galaxies (M$_* \sim 10^{10.8}$ M$_\odot$) and present a fitting function for the volume number density of clusters as a function of parent galaxy stellar mass. We find no correlation between T$_{\rm N}$ and large-scale environment, but do find a tendency for galaxies of fixed M$_*$ to have larger T$_{\rm N}$ if they have converted a larger proportion of their baryons into stars.

Resonance rings are the consequence of secular evolution processes that redistribute material and angular momentum in discs. We produced a Catalogue and an Atlas of the rings detected in the Spitzer Survey of Stellar Structure in Galaxies (S4G) and to conduct a statistical study of the data in the Catalogue. We traced the contours of rings previously identified by Buta et al. (in preparation) and fitted them with ellipses. We found the orientation of bars by studying the galaxy ellipse fits from S4G's Pipeline 4. We used the galaxy orientation data obtained by S4G's Pipeline 4 to obtain intrinsic ellipticities and orientations of rings and the bars. ARRAKIS contains data on 724 ringed galaxies in the S4G. The frequency of resonance rings in the S4G is of 16+-1% and 35+-1% for outer and inner features, respectively. Outer rings are mostly found in Hubble stages -1<=T<=4. Inner rings are found in a distribution that covers the range -1<=T<=7. We confirm that outer rings have two preferred orientations, parallel and perpendicular to the bar. We confirm a tendency for inner rings to be oriented parallel to the bar, but we find that a significant fraction (~50%) of them have random orientations with respect to the bar. These misaligned inner rings are mostly found in late-type galaxies (T>=4). This may be due to spiral modes decoupled from the bar dominating the Fourier amplitude spectrum at the radius of the inner ring. We find that the fraction of barred galaxies hosting outer (inner) rings is ~1.7 times (~1.3 times) that in unbarred galaxies. The fact that rings are only mildly favoured by bars suggests that those in unbarred galaxies either formed due to weak departures from the axisymmetry of the galactic potential or that they are born because of bars that have been destroyed after the ring formation.

The mid-infrared is an optimal window to trace stellar mass in nearby galaxies and the 3.6$\mu 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$\mu m$ and obtain stellar mass maps for more than 1600 galaxies available from the Spitzer Survey of Stellar Structure in Galaxies (S$^{4}$G). This survey consists of images in two infrared bands (3.6 and 4.5$\mu m$), and we use the Independent Component Analysis (ICA) method presented in Meidt et al. (2012) to separate the dominant light from old stars and the dust emission that can significantly contribute to the observed 3.6$\mu 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$\mu 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: $\log(M/L)=-0.339 (\pm 0.057) \times ([3.6]-[4.5]) -0.336 (\pm 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.

We use mid-infrared 3.6 and 4.5microns imaging of NGC 3906 from the Spitzer Survey of Stellar Structure in Galaxies (S4G) to understand the nature of an unusual offset between its stellar bar and the photometric center of an otherwise regular, circular outer stellar disk. We measure an offset of ~720 pc between the center of the stellar bar and photometric center of the stellar disk; the bar center coincides with the kinematic center of the disk determined from previous HI observations. Although the undisturbed shape of the disk suggests that NGC 3906 has not undergone a significant merger event in its recent history, the most plausible explanation for the observed offset is an interaction. Given the relatively isolated nature of NGC 3906 this interaction could be with dark matter sub structure in the galaxy's halo or from a recent interaction with a fast moving neighbor which remains to be identified. Simulations aimed at reproducing the observed offset between the stellar bar / kinematic center of the system and the photometric center of the disk are necessary to confirm this hypothesis and constrain the interaction history of the galaxy.

The Spitzer Survey of Stellar Structure in Galaxies (S4G) is a volume, magnitude, and size-limited survey of 2352 nearby galaxies with deep imaging at 3.6 and 4.5um. In this paper we describe our surface photometry pipeline and showcase the associated data products that we have released to the community. We also identify the physical mechanisms leading to different levels of central stellar mass concentration for galaxies with the same total stellar mass. Finally, we derive the local stellar mass-size relation at 3.6um for galaxies of different morphologies. Our radial profiles reach stellar mass surface densities below 1 Msun pc-2. Given the negligible impact of dust and the almost constant mass-to-light ratio at these wavelengths, these profiles constitute an accurate inventory of the radial distribution of stellar mass in nearby galaxies. From these profiles we have also derived global properties such as asymptotic magnitudes (and the corresponding stellar masses), isophotal sizes and shapes, and concentration indices. These and other data products from our various pipelines (science-ready mosaics, object masks, 2D image decompositions, and stellar mass maps), can be publicly accessed at IRSA (http://irsa.ipac.caltech.edu/data/SPITZER/S4G/).

The Spitzer Survey of Stellar Structure in Galaxies (S$^4$G, Sheth et. al. 2010) is a deep 3.6 and 4.5 $\mu$m imaging survey of 2352 nearby ($< 40$ Mpc) galaxies. We describe the S$^4$G data analysis pipeline 4, which is dedicated to 2-dimensional structural surface brightness decompositions of 3.6 $\mu$m images, using GALFIT3.0 \citep{peng2010}. Besides automatic 1-component S\'ersic fits, and 2-component S\'ersic 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 S\'ersic index and bulge-to-total light ratio ($B/T$), confirming earlier results \citep{laurikainen2007, gadotti2008, weinzirl2009}. In this first paper, 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 ({\bf \tt 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. In the second paper we will present a detailed analysis of the bulge, disk, and bar parameter derived from multi-component decompositions.

We obtained GALEX FUV, NUV, and Spitzer/IRAC 3.6$\mu$m photometry for > 2000 galaxies, available for 90% of the S4G sample. We find a very tight "GALEX Blue Sequence (GBS)" in the (FUV-NUV) versus (NUV-[3.6]) color-color diagram which is populated by irregular and spiral galaxies, and is mainly driven by changes in the formation timescale ($\tau$) and a degeneracy between $\tau$ and dust reddening. The tightness of the GBS provides an unprecedented way of identifying star-forming galaxies and objects that are just evolving to (or from) what we call the "GALEX Green Valley (GGV)". At the red end of the GBS, at (NUV-[3.6]) > 5, we find a wider "GALEX Red Sequence (GRS)" mostly populated by E/S0 galaxies that has a perpendicular slope to that of the GBS and of the optical red sequence. We find no such dichotomy in terms of stellar mass (measured by $\rm{M}_{[3.6]}$), since both massive ($M_{\star} > 10^{11} M_{\odot}$) blue and red sequence galaxies are identified. The type that is proportionally more often found in the GGV are the S0-Sa's and most of these are located in high-density environments. We discuss evolutionary models of galaxies that show a rapid transition from the blue to the red sequence on timescale of $10^{8}$years.

The Spitzer Survey of Stellar Structure in Galaxies (S4G) is the largest available database of deep, homogeneous middle-infrared (mid-IR) images of galaxies of all types. The survey, which includes 2352 nearby galaxies, reveals galaxy morphology only minimally affected by interstellar extinction. This paper presents an atlas and classifications of S4G galaxies in the Comprehensive de Vaucouleurs revised Hubble-Sandage (CVRHS) system. The CVRHS system follows the precepts of classical de Vaucouleurs (1959) morphology, modified to include recognition of other features such as inner, outer, and nuclear lenses, nuclear rings, bars, and disks, spheroidal galaxies, X patterns and box/peanut structures, OLR subclass outer rings and pseudorings, bar ansae and barlenses, parallel sequence late-types, thick disks, and embedded disks in 3D early-type systems. We show that our CVRHS classifications are internally consistent, and that nearly half of the S4G sample consists of extreme late-type systems (mostly bulgeless, pure disk galaxies) in the range Scd-Im. The most common family classification for mid-IR types S0/a to Sc is SA while that for types Scd to Sm is SB. The bars in these two type domains are very different in mid-IR structure and morphology. This paper examines the bar, ring, and type classification fractions in the sample, and also includes several montages of images highlighting the various kinds of "stellar structures" seen in mid-IR galaxy morphology.

Using 3.6$\mu$m images of 97 early-type galaxies, we develop and verify methodology to measure globular cluster populations from the S$^4$G survey images. We find that 1) the ratio, T$_{\rm N}$, of the number of clusters, N$_{\rm CL}$, to parent galaxy stellar mass, M$_*$, rises weakly with M$_*$ for early-type galaxies with M$_* > 10^{10}$ M$_\odot$ when we calculate galaxy masses using a universal stellar initial mass function (IMF), but that the dependence of T$_{\rm N}$ on M$_*$ is removed entirely once we correct for the recently uncovered systematic variation of IMF with M$_*$, and 2) for M$_* < 10^{10}$ M$_\odot$ there is no trend between N$_{\rm CL}$ and M$_*$, the scatter in T$_{\rm N}$ is significantly larger (approaching 2 orders of magnitude), and there is evidence to support a previous, independent suggestion of two families of galaxies. The behavior of N$_{\rm CL}$ in the lower mass systems is more difficult to measure because these systems are inherently cluster poor, but our results may add to previous evidence that large variations in cluster formation and destruction efficiencies are to be found among low mass galaxies. The average fraction of stellar mass in clusters is $\sim$ 0.0014 for M$_* > 10^{10}$ M$_\odot$ and can be as large as $\sim 0.02$ for less massive galaxies. These are the first results from the S$^4$G sample of galaxies, and will be enhanced by the sample of early-type galaxies now being added to S$^4$G and complemented by the study of later type galaxies within S$^4$G.

We have measured the radial light profiles and global shapes of bars using two-dimensional 3.6 $\mu m $ image decompositions for 144 face-on barred galaxies from the Spitzer Survey of Stellar Structure in Galaxies (S4G). The bar surface brightness profile is correlated with the stellar mass and bulge-to-total (B/T) ratio of their host galaxies. Bars in massive and bulge-dominated galaxies (B/T$>$0.2) show a flat profile, while bars in less massive, disk-dominated galaxies (B/T$\sim$0) show an exponential, disk-like profile with a wider spread in the radial profile than in the bulge-dominated galaxies. The global two-dimensional shapes of bars, however, are rectangular/boxy, independent of the bulge or disk properties. We speculate that because bars are formed out of disk, bars initially have an exponential (disk-like) profile which evolves over time, trapping more stars into the boxy bar orbits. This leads bars to become stronger and have flatter profiles. The narrow spread of bar radial profiles in more massive disks suggests that these bars formed earlier (z$>$1), while the disk-like profiles and a larger spread in the radial profile in less massive systems imply a later and more gradual evolution, consistent with the cosmological evolution of bars inferred from observational studies. Therefore, we expect that the flatness of the bar profile can be used as a dynamical age indicator of the bar to measure the time elapsed since the bar formation. We argue that cosmic gas accretion is required to explain our results on bar profile and the presence of gas within the bar region.

We present a catalogue and images of visually detected features, such as asymmetries, extensions, warps, shells, tidal tails, polar rings, and obvious signs of mergers or interactions, in the faint outer regions (at and outside of R_25) of nearby galaxies. This catalogue can be used in future quantitative studies that examine galaxy evolution due to internal and external factors. We are able to reliably detect outer region features down to a brightness level of 0.03 MJy/sr per pixel at 3.6 microns in the Spitzer Survey of Stellar Structure in Galaxies (S4G). We also tabulate companion galaxies. We find asymmetries in the outer isophotes in 22+/-1 per cent of the sample. The asymmetry fraction does not correlate with galaxy classification as an interacting galaxy or merger remnant, or with the presence of companions. We also compare the detected features to similar features in galaxies taken from cosmological zoom re-simulations. The simulated images have a higher fraction (33 per cent) of outer disc asymmetries, which may be due to selection effects and an uncertain star formation threshold in the models. The asymmetries may have either an internal (e.g., lopsidedness due to dark halo asymmetry) or external origin.

We combine data from the Spitzer Survey for Stellar Structure in Galaxies (S$^4$G), a recently calibrated empirical stellar mass estimator from Eskew et al., and an extensive database of HI spectral line profiles to examine the baryonic Tully-Fisher (BTF) relation. We find 1) that the BTF has lower scatter than the classic Tully-Fisher (TF) relation and is better described as a linear relationship, confirming similar previous results, 2) that the inclusion of a radial scale in the BTF decreases the scatter but only modestly, as seen previously for the TF relation, and 3) that the slope of the BTF, which we find to be $3.5\pm 0.2$ ($\Delta$ log $M_{baryon}/\Delta$ log $v_c$), implies that on average a nearly constant fraction ($\sim 0.4$) of all baryons expected to be in a halo are "condensed" onto the central region of rotationally supported galaxies. The condensed baryon fraction, $M_{baryon}/M_{total}$, is, to our measurement precision, nearly independent of galaxy circular velocity (our sample spans circular velocities, $v_c$, between 60 and 250 km s$^{-1}$, but is extended to $v_c\sim 10$ km s$^{-1}$ using data from the literature). The observed galaxy-to-galaxy scatter in this fraction is generally $\le$ a factor of 2 despite fairly liberal selection criteria. These results imply that cooling and heating processes, such as cold vs. hot accretion, mass loss due to stellar winds, and AGN driven feedback, to the degree that they affect the global galactic properties involved in the BTF, are independent of halo mass for galaxies with $10 < v_c < 250$ km/s and typically introduce no more than a factor of two range in the resulting $M_{baryon}/M_{total}$. Recent simulations by Aumer et al. of a small sample of disk galaxies are in excellent agreement with our data, suggesting that current simulations are capable of reproducing the global properties of individual disk galaxies.

We present a new approach for estimating the 3.6 micron stellar mass-to-light ratio in terms of the [3.6]-[4.5] colors of old stellar populations. Our approach avoids several of the largest sources of uncertainty in existing techniques. By focusing on mid-IR wavelengths, we gain a virtually dust extinction-free tracer of the old stars, avoiding the need to adopt a dust model to correctly interpret optical or optical/NIR colors normally leveraged to assign M/L. By calibrating a new relation between NIR and mid-IR colors of GLIMPSE giant stars we also avoid discrepancies in model predictions for the [3.6]-[4.5] colors of old stellar populations due to uncertainties in molecular line opacities. We find that the [3.6]-[4.5] color, which is driven primarily by metallicity, provides a tight constraint on M/L_3.6, which varies intrinsically less than at optical wavelengths. The uncertainty on M/L_3.6 of ~0.07 dex due to unconstrained age variations marks a significant improvement on existing techniques for estimating the stellar M/L with shorter wavelength data. A single M/L_3.6=0.6 (assuming a Chabrier IMF), independent of [3.6]-[4.5] color, is also feasible as it can be applied simultaneously to old, metal-rich and young, metal-poor populations, and still with comparable (or better) accuracy (~0.1 dex) as alternatives. We expect our M/L_3.6 to be optimal for mapping the stellar mass distributions in S4G galaxies, for which we have developed an Independent Component Analysis technique to first isolate the old stellar light at 3.6 micron from non-stellar emission (e.g. hot dust and the 3.3 PAH feature). Our estimate can also be used to determine the fractional contribution of non-stellar emission to global (rest-frame) 3.6 micron fluxes, e.g. in WISE imaging, and establishes a reliable basis for exploring variations in the stellar IMF.

We have performed two-dimensional multicomponent decomposition of 144 local barred spiral galaxies using 3.6 $\mu {\rm m}$ images from the Spitzer Survey of Stellar Structure in Galaxies. Our model fit includes up to four components (bulge, disk, bar, and a point source) and, most importantly, takes into account disk breaks. We find that ignoring the disk break and using a single disk scale length in the model fit for Type II (down-bending) disk galaxies can lead to differences of 40% in the disk scale length, 10% in bulge-to-total luminosity ratio (B/T), and 25% in bar-to-total luminosity ratios. We find that for galaxies with B/T $\geq$ 0.1, the break radius to bar radius, $r_{\rm br}/R_{\rm bar}$, varies between 1 and 3, but as a function of B/T the ratio remains roughly constant. This suggests that in bulge-dominated galaxies the disk break is likely related to the outer Lindblad Resonance (OLR) of the bar, and thus moves outwards as the bar grows. For galaxies with small bulges, B/T $<$ 0.1, $r_{\rm br}/R_{\rm bar}$ spans a wide range from 1 to 6. This suggests that the mechanism that produces the break in these galaxies may be different from that in galaxies with more massive bulges. Consistent with previous studies, we conclude that disk breaks in galaxies with small bulges may originate from bar resonances that may be also coupled with the spiral arms, or be related to star formation thresholds.

Star-forming regions that are visible at 3.6 microns and Halpha but not in the u,g,r,i,z bands of the Sloan Digital Sky survey (SDSS), are measured in five nearby spiral galaxies to find extinctions averaging ~3.8 mag and stellar masses averaging ~5x10^4 Msun. These regions are apparently young star complexes embedded in dark filamentary shock fronts connected with spiral arms. The associated cloud masses are ~10^7 Msun. The conditions required to make such complexes are explored, including gravitational instabilities in spiral shocked gas and compression of incident clouds. We find that instabilities are too slow for a complete collapse of the observed spiral filaments, but they could lead to star formation in the denser parts. Compression of incident clouds can produce a faster collapse but has difficulty explaining the semi-regular spacing of some regions along the arms. If gravitational instabilities are involved, then the condensations have the local Jeans mass. Also in this case, the near-simultaneous appearance of equally spaced complexes suggests that the dust lanes, and perhaps the arms too, are relatively young.

The morphology of galaxies can be quantified to some degree using a set of scale-invariant parameters. Concentration (C), Asymmetry (A), Smoothness (S), the Gini index (G), relative contribution of the brightest pixels to the second order moment of the flux (M20), ellipticity (E), and the Gini index of the second order moment (GM) have all been applied to morphologically classify galaxies at various wavelengths. Here we present a catalog of these parameters for the Spitzer Survey of Stellar Structure in Galaxies (S4G), a volume-limited near-infrared imaging survey of nearby galaxies using the 3.6 and 4.5 micron channels of the IRAC camera. Our goal is to provide a reference catalog of near-infrared quantified morphology for high-redshift studies and galaxy evolution models with enough detail to resolve stellar mass morphology. We explore where normal, non-interacting galaxies -those typically found on the Hubble tuning fork- lie in this parameter space and show that there is a tight relation between Concentration and M20 for normal galaxies. M20 can be used to classify galaxies into earlier and later types (e.g., to separate spirals from irregulars). Several criteria using these parameters exist to select systems with a disturbed morphology, i.e., those that appear to be undergoing a tidal interaction. We examine the applicability of these criteria to Spitzer near-infrared imaging. We find that four relations, based on the parameters A&S, G&M20, GM, and C&M20, respectively, select outliers in the morphological parameter space, but each selects different subsets of galaxies. Two criteria (GM > 0.6, G > -0.115 x M20 + 0.384) seem most appropriate to identify possible mergers and the merger fraction in near-infrared surveys. We find no strong relation between lopsidedness and most of these morphological parameters, except for a weak dependence of lopsidedness on Concentration and M20.

Stellar bars can lead to gas inflow toward the center of a galaxy and stimulate nuclear star formation. However, there is no compelling evidence on whether they also feed a central supermassive black hole: by measuring the fractions of barred active and inactive galaxies, previous studies have yielded conflicting results. In this paper, we aim to understand the lack of observational evidence for bar-driven active galactic nucleus (AGN) activity by studying a sample of 41 nearby (d < 35 Mpc) barred galaxies from the Spitzer Survey for Stellar Structure in Galaxies. We use Chandra observations to measure nuclear 2--10 keV X-ray luminosities and estimate Eddington ratios, together with Spitzer 3.6um imaging to quantify the strength of the stellar bar in two independent ways: (1) from its structure, as traced by its ellipticity and boxiness, and (2) from its gravitational torque Q_b, taken as the maximum ratio of the tangential force to the mean background radial force. In this way, rather than discretizing the presence of both stellar bars and nuclear activity, we are able to account for the continuum of bar strengths and degrees of AGN activity. We find nuclear X-ray sources in 31 out of 41 galaxies with median X-ray luminosity and Eddington ratio of L_X=4.3x10^{38} erg/s and L_bol/L_Edd=6.9x10^{-6} respectively, consistent with low-luminosity AGN activity. Including upper limits for those galaxies without nuclear detections, we find no significant correlation between any of the bar strength indicators and the degree of nuclear activity, irrespective of galaxy luminosity, stellar mass, Hubble type, or bulge size. Strong bars do not favor brighter or more efficient nuclear activity, implying that at least for the low-luminosity regime, supermassive black hole fueling is not closely connected to large scale features.

We study the m=1 distortions (lopsidedness) in the stellar components of 167 nearby galaxies that span a wide range of morphologies and luminosities. We confirm the previous findings of 1) a high incidence of lopsidedness in the stellar distributions, 2) increasing lopsidedness as a function of radius out to at least 3.5 exponential scale lengths, and 3) greater lopsidedness, over these radii, for galaxies of later type and lower surface brightness. Additionally, the magnitude of the lopsidedness 1) correlates with the character of the spiral arms (stronger arm patterns occur in galaxies with less lopsidedness), 2) is not correlated with the presence or absence of a bar, or the strength of the bar when one is present, 3) is inversely correlated to the stellar mass fraction, f_*, within one radial scale length, and 4) correlates directly with f_* measured within the radial range over which we measure lopsidedness. We interpret these findings to mean that lopsidedness is a generic feature of galaxies and does not, generally, depend on a rare event, such as a direct accretion of a satellite galaxy onto the disk of the parent galaxy. While lopsidedness may be caused by several phenomena, moderate lopsidedness (<A_1>_i + <A_1>_o)/2 < 0.3) is likely to reflect halo asymmetries to which the disk responds or a gravitationally self-generated mode . We hypothesize that the magnitude of the stellar response depends both on how centrally concentrated the stars are with respect to the dark matter and whether there are enough stars in the region of the lopsidedness that self-gravity is dynamically important.

We have analyzed the radial distribution of old stars in a sample of 218 nearby face-on disks, using deep 3.6um images from the Spitzer Survey of Stellar Structure in Galaxies (S4G). 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 case of barred galaxies, the ratio between the break and the bar radii (Rbr/Rbar) depends strongly on the total stellar mass of the galaxy. For galaxies more massive than 10^10 Msun, the distribution is bimodal, peaking at Rbr/Rbar~2 and ~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 Msun, breaks are found up to ~10 Rbar, 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-axysymmetric features can be responsible not only for those breaks at 2 Rbar, but also for many of those found at larger radii.

We present a study of the kinematics of the isolated spiral galaxy NGC 864, using H{\alpha} Fabry-Perot data obtained with the GH{\alpha}FaS instrument at the William Herschel Telescope in La Palma, complemented with images at 3.6 {\mu}m, in the R band and in H{\alpha} filter, and integral field spectroscopic data. The resulting data cubes and velocity maps allow the study of the kinematics of the galaxy, including in-depth investigations of the rotation curve, velocity moment maps, velocity residual maps, gradient maps and position-velocity diagrams. We find asymmetries in the velocity field in the bar zone, caused by non-circular motions, probably in response to the potential of the bar. We also find a flat-profile bar, in agreement with the strong bar, with the grand design spiral pattern, and with the gap between the ends of the bar and the start of the spiral arms. We quantify the rate of massive star formation, which is concentrated in the two spiral arms.

The mechanism causing breaks in the radial surface brightness distribution of spiral galaxies is not yet well known. Despite theoretical efforts, there is not a unique explanation for these features and the observational results are not conclusive. In an attempt to address this problem, we have selected a sample of 34 highly inclined spiral galaxies present both in the Sloan Digital Sky Survey and in the Spitzer Survey of Stellar Structure in Galaxies. We have measured the surface brightness profiles in the five Sloan optical bands and in the 3.6$\mu m$ Spitzer band. We have also calculated the color and stellar surface mass density profiles using the available photometric information, finding two differentiated features: an innermost break radius at distances of $\sim 8 \pm 1$ kpc [$0.77 \pm 0.06$ $R_{25}$] and a second characteristic radius, or truncation radius, close to the outermost optical extent ($\sim 14 \pm 2$ kpc [$1.09 \pm 0.05$ $R_{25}$]) of the galaxy. We propose in this work that the breaks might be a phenomena related to a threshold in the star formation, while truncations are more likely a real drop in the stellar mass density of the disk associated with the maximum angular momentum of the stars.

Breaks in the radial luminosity profiles of galaxies have been until now mostly studied averaged over discs. Here we study separately breaks in thin and thick discs in 70 edge-on galaxies using imaging from the Spitzer Survey of Stellar Structure in Galaxies. We built luminosity profiles of the thin and the thick discs parallel to midplanes and we found that thin discs often truncate (77%). Thick discs truncate less often (31%), but when they do, their break radius is comparable with that in the thin disc. This suggests either two different truncation mechanisms - one of dynamical origin affecting both discs simultaneously and another one only affecting the thin disc - or a single mechanism that creates a truncation in one disc or in both depending on some galaxy property. Thin discs apparently antitruncate in around 40% of galaxies. However, in many cases, these antitruncations are an artifact caused by the superposition of a thin disc and a thick disc with the latter having a longer scale length. We estimate the real thin disc antitruncation fraction to be less than 15%. We found that the ratio of the thick and thin stellar disc mass is roughly constant (0.2<M_T/M_t<0.7) for circular velocities v_c>120 km/s, but becomes much larger at smaller velocities. We hypothesize that this is due to a combination of a high efficiency of supernova feedback and a slower dynamical evolution in lower-mass galaxies causing stellar thin discs to be younger and less massive than in higher-mass galaxies.

Tidal debris around galaxies can yield important clues on their evolution. We have identified tidal debris in 11 early type galaxies (T \leq 0) from a sample of 65 early types drawn from the Spitzer Survey of Stellar Structure in Galaxies (S4G). The tidal debris includes features such as shells, ripples and tidal tails. A variety of techniques, including two-dimensional decomposition of galactic structures, was used to quantify the residual tidal features. The tidal debris contributes ~3 - 10% to the total 3.6 {\mu}m luminosity of the host galaxy. Structural parameters of the galaxies were estimated using two-dimensional profile fitting. We investigate the locations of galaxies with tidal debris in the Fundamental Plane and Kormendy relation. We find that galaxies with tidal debris lie within the scatter of early type galaxies without tidal features. Assuming that the tidal debris is indicative of recent gravitational interaction or merger, this suggests that these galaxies have either undergone minor merging events so that the overall structural properties of the galaxies are not significantly altered, or they have undergone a major merging events but already have experienced sufficient relaxation and phase-mixing so that their structural properties become similar to those of the non-interacting early type galaxies.

We examine the effect of circumstellar dust extinction on the near-IR contribution of asymptotic giant branch (AGB) stars in intermediate-age clusters throughout the disk of M100. For our sample of 17 AGB-dominated clusters we extract optical-to-mid-IR SEDs and find that NIR brightness is coupled to the mid-IR dust emission in such a way that a significant reduction of AGB light, of up to 1 mag in K-band, follows from extinction by the dust shell formed during this stage. Since the dust optical depth varies with AGB chemistry (C-rich or O-rich), our results suggest that the contribution of AGB stars to the flux from their host clusters will be closely linked to the metallicity and the progenitor mass of the AGB star, to which dust chemistry and mass-loss rate are sensitive. Our sample of clusters--each the analogue of a ~1 Gyr old post-starburst galaxy--has implications within the context of mass and age estimation via SED modelling at high z: we find that the average ~0.5 mag extinction estimated here may be sufficient to reduce the AGB contribution in (rest-frame) K-band from ~70%, as predicted in the latest generation of synthesis models, to ~35%. Our technique for selecting AGB-dominated clusters in nearby galaxies promises to be effective for discriminating the uncertainties associated with AGB stars in intermediate-age populations that plague age and mass estimation in high-z galaxies.

Modern adaptive optics (AO) systems for large telescopes require tomographic techniques to reconstruct the phase aberrations induced by the turbulent atmosphere along a line of sight to a target which is angularly separated from the guide sources that are used to sample the atmosphere. Multi-object adaptive optics (MOAO) is one such technique. Here, we present a method which uses an artificial neural network (ANN) to reconstruct the target phase given off-axis references sources. We compare our ANN method with a standard least squares type matrix multiplication method and to the learn and apply method developed for the CANARY MOAO instrument. The ANN is trained with a large range of possible turbulent layer positions and therefore does not require any input of the optical turbulence profile. It is therefore less susceptible to changing conditions than some existing methods. We also exploit the non-linear response of the ANN to make it more robust to noisy centroid measurements than other linear techniques.

With the aim of constructing accurate 2D maps of the stellar mass distribution in nearby galaxies from S4G 3.6 and 4.5 micron images, we report on the separation of the light from old stars from the emission contributed by contaminants (e.g. hot dust and the 3.3 micron PAH feature). 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 contents and star formation histories are presented to demonstrate our approach. 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 intermediate-age 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. Contaminants are identified via comparison to the non-stellar emission imaged at 8 microns, which is dominated by the broad PAH feature. Using the measured 3.6/8 micron ratio to select the individual contaminants, we find that hot dust and PAH together contribute between ~5-15% to the integrated light at 3.6 microns, while light from regions dominated by intermediate-age 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 upwards 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 the stellar mass, as will be the focus in a companion paper.

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 (S4G). We found that thick disks are on average more massive than previously reported, mostly due to the selected fitting function. Typically, the thin and the thick disk 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 scaleheights 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 build-up 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".

NGC 4013 is a nearby Sb edge-on galaxy known for its "prodigious" HI warp and its "giant" tidal stream. Previous work on this unusual object shows that it cannot be fitted satisfactorily by a canonical thin+thick disk structure. We have produced a new decomposition of NGC 4013, considering three stellar flattened components (thin+thick disk plus an extra and more extended component) and one gaseous disk. All four components are considered to be gravitationally coupled and isothermal. To do so, we have used the 3.6micron images from the Spitzer Survey of Stellar Structure in Galaxies (S4G). We find evidence for NGC4013 indeed having a thin and a thick disk and an extra flattened component. This smooth and extended component (scaleheight z_EC kpc) could be interpreted as a thick disk or as a squashed ellipsoidal halo and contains ~20% of the total mass of all three stellar components. We argue it is unlikely to be related to the ongoing merger or due to the off-plane stars from a warp in the other two disk components. Instead, we favor a scenario in which the thick disk and the extended component were formed in a two-stage process, in which an initially thick disk has been dynamically heated by a merger soon enough in the galaxy history to have a new thick disk formed within it.

The Spitzer Survey of Stellar Structure in Galaxies S^4G is an Exploration Science Legacy Program approved for the Spitzer post-cryogenic mission. It is a volume-, magnitude-, and size-limited (d < 40 Mpc, |b| > 30 degrees, m_(Bcorr) < 15.5, D25>1') survey of 2,331 galaxies using IRAC at 3.6 and 4.5 microns. Each galaxy is observed for 240 s and mapped to > 1.5 x D25. The final mosaicked images have a typical 1 sigma rms noise level of 0.0072 and 0.0093 MJy / sr at 3.6 and 4.5 microns, respectively. Our azimuthally-averaged surface brightness profile typically traces isophotes at mu_3.6 (AB) (1 sigma) ~ 27 mag arcsec^-2, equivalent to a stellar mass surface density of ~ 1 Msun pc^-2. S^4G thus provides an unprecedented data set for the study of the distribution of mass and stellar structures in the local Universe. This paper introduces the survey, the data analysis pipeline and measurements for a first set of galaxies, observed in both the cryogenic and warm mission phase of Spitzer. For every galaxy we tabulate the galaxy diameter, position angle, axial ratio, inclination at mu_3.6 (AB) = 25.5 and 26.5 mag arcsec^-2 (equivalent to ~ mu_B (AB) =27.2 and 28.2 mag arcsec^-2, respectively). These measurements will form the initial S^4G catalog of galaxy properties. We also measure the total magnitude and the azimuthally-averaged radial profiles of ellipticity, position angle, surface brightness and color. Finally, we deconstruct each galaxy using GALFIT into its main constituent stellar components: the bulge/spheroid, disk, bar, and nuclear point source, where necessary. Together these data products will provide a comprehensive and definitive catalog of stellar structures, mass and properties of galaxies in the nearby Universe.