GALAXIES BEHIND THE GALACTIC PLANE: FIRST RESULTS AND PERSPECTIVES FROM THE VVV SURVEY

VISTA Variables in the Via Lactea (VVV) is an ESO variability survey that is performing observations in near-infrared bands (ZY JHK(s)) toward the Galactic bulge and part of the disk with the completeness limits at least 3 mag deeper than Two Micron All Sky Survey. In the present work, we searched in the VVV survey data for background galaxies near the Galactic plane using ZY JHK(s) photometry that covers 1.636 deg(2). We identified 204 new galaxy candidates by analyzing colors, sizes, and visual inspection of multi-band (ZY JHK(s)) images. The galaxy candidate colors were also compared with the predicted ones by star count models considering a more realistic extinction model at the same completeness limits observed by VVV. A comparison of the galaxy candidates with the expected one by Millennium simulations is also presented. Our results increase the number density of known galaxies behind the Milky Way by more than one order of magnitude. A catalog with galaxy properties including ellipticity, Petrosian radii, and ZY JHK(s) magnitudes is provided, as well as comparisons of the results with other surveys of galaxies toward the Galactic plane.

The Universe 10 Billion years ago: Morphologies, Star-formation rates and Galactic-scale winds in z ∼ 2 Galaxies

This thesis is devoted to the study of the properties of high-redsfhit galaxies in the epoch 1 < z < 3, when a substantial fraction of galaxy mass was assembled, and when the evolution of the star-formation rate density peaked. Following a multi-perspective approach and using the most recent and high-quality data available (spectra, photometry and imaging), the morphologies and the star-formation properties of high-redsfhit galaxies were investigated. Through an accurate morphological analyses, the built up of the Hubble sequence was placed around z ~ 2.5. High-redshift galaxies appear, in general, much more irregular and asymmetric than local ones. Moreover, the occurrence of morphological k-­correction is less pronounced than in the local Universe. Different star-formation rate indicators were also studied. The comparison of ultra-violet and optical based estimates, with the values derived from infra-red luminosity showed that the traditional way of addressing the dust obscuration is problematic, at high-redshifts, and new models of dust geometry and composition are required. Finally, by means of stacking techniques applied to rest-frame ultra-violet spectra of star-forming galaxies at z~2, the warm phase of galactic-scale outflows was studied. Evidence was found of escaping gas at velocities of ~ 100 km/s. Studying the correlation of inter-­stellar absorption lines equivalent widths with galaxy physical properties, the intensity of the outflow-related spectral features was proven to depend strongly on a combination of the velocity dispersion of the gas and its geometry.

Tracing AGN accretion and star formation in Herschel galaxies

In this Thesis, we study the physical properties and the cosmic evolution of AGN and their host galaxies since z∼3. Our analysis exploits samples of star forming galaxies detected with Herschel at far-IR wavelengths (from 70 up to 500 micron) in different extragalactic surveys, such as COSMOS and the deep GOODS (South and North) fields. The broad-band ancillary data available in COSMOS and the GOODS fields, allows us to implement Herschel and Spitzer photometry with multi-wavelength ancillary data. We perform a multicomponent SED-fitting decomposition to decouple the emission due to star formation from that due to AGN accretion, and to estimate both host-galaxy parameters (such as stellar mass, M* and star formation rate, SFR), and nuclear intrinsic bolometric luminosities. We use the individual estimates of AGN bolometric luminosity obtained through SED-fitting decomposition to reconstruct the redshit evolution of the AGN bolometric luminosity function since z∼3. The resulting trends are used to estimate the overall AGN accretion rate density at different cosmic epochs and to trace the first ever estimate of the AGN accretion history from an IR survey. Later on, we focus our study on the connection between AGN accretion and integrated galaxy properties. We analyse the relationships of AGN accretion with galaxy properties in the SFR-M* plane and at different cosmic epochs. Finally, we infer what is the parameter that best correlates with AGN accretion, comparing our results with previous studies and discussing their physical implications in the context of current scenarios of AGN/galaxy evolution.

CANDELS: The progenitors of compact quiescent galaxies at z ~ 2

We combine high-resolution Hubble Space Telescope/WFC3 images with multi-wavelength photometry to track the evolution of structure and activity of massive (M_*> 10^10 M_☉) galaxies at redshifts z = 1.4-3 in two fields of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey. We detect compact, star-forming galaxies (cSFGs) whose number densities, masses, sizes, and star formation rates (SFRs) qualify them as likely progenitors of compact, quiescent, massive galaxies (cQGs) at z = 1.5-3. At z≲2, cSFGs present SFR = 100-200 M_☉ yr^–1, yet their specific star formation rates (sSFR ~ 10^–9 yr^–1) are typically half that of other massive SFGs at the same epoch, and host X-ray luminous active galactic nuclei (AGNs) 30 times (~30%) more frequently. These properties suggest that cSFGs are formed by gas-rich processes (mergers or disk-instabilities) that induce a compact starburst and feed an AGN, which, in turn, quench the star formation on dynamical timescales (few 10^8 yr). The cSFGs are continuously being formed at z = 2-3 and fade to cQGs down to z ~ 1.5. After this epoch, cSFGs are rare, thereby truncating the formation of new cQGs. Meanwhile, down to z = 1, existing cQGs continue to enlarge to match local QGs in size, while less-gas-rich mergers and other secular mechanisms shepherd (larger) SFGs as later arrivals to the red sequence. In summary, we propose two evolutionary tracks of QG formation: an early (z≲2), formation path of rapidly quenched cSFGs fading into cQGs that later enlarge within the quiescent phase, and a late-arrival (z≳2) path in which larger SFGs form extended QGs without passing through a compact state.

Testing diagnostics of nuclear activity and star formation in galaxies at z > 1

We present some of the first science data with the new Keck/MOSFIRE instrument to test the effectiveness of different AGN/SF diagnostics at z ~ 1.5. MOSFIRE spectra were obtained in three H-band multi-slit masks in the GOODS-S field, resulting in 2 hr exposures of 36 emission-line galaxies. We compare X-ray data with the traditional emission-line ratio diagnostics and the alternative mass-excitation and color-excitation diagrams, combining new MOSFIRE infrared data with previous HST/WFC3 infrared spectra (from the 3D-HST survey) and multiwavelength photometry. We demonstrate that a high [O III]/Hβ ratio is insufficient as an active galactic nucleus (AGN) indicator at z > 1. For the four X-ray-detected galaxies, the classic diagnostics ([O III]/Hβ versus [N II]/Hα and [S II]/Hα) remain consistent with X-ray AGN/SF classification. The X-ray data also suggest that "composite" galaxies (with intermediate AGN/SF classification) host bona fide AGNs. Nearly ~2/3 of the z ~ 1.5 emission-line galaxies have nuclear activity detected by either X-rays or the classic diagnostics. Compared to the X-ray and line ratio classifications, the mass-excitation method remains effective at z > 1, but we show that the color-excitation method requires a new calibration to successfully identify AGNs at these redshifts.

Calibrating photometric redshifts of luminous red galaxies

We discuss the construction of a photometric redshift catalogue of luminous red galaxies (LRGs) from the Sloan Digital Sky Survey (SDSS), emphasizing the principal steps necessary for constructing such a catalogue: (i) photometrically selecting the sample, (ii) measuring photometric redshifts and their error distributions, and (iii) estimating the true redshift distribution. We compare two photometric redshift algorithms for these data and find that they give comparable results. Calibrating against the SDSS and SDSS-2dF (Two Degree Field) spectroscopic surveys, we find that the photometric redshift accuracy is sigma similar to 0.03 for redshifts less than 0.55 and worsens at higher redshift (similar to 0.06 for z < 0.7). These errors are caused by photometric scatter, as well as systematic errors in the templates, filter curves and photometric zero-points. We also parametrize the photometric redshift error distribution with a sum of Gaussians and use this model to deconvolve the errors from the measured photometric redshift distribution to estimate the true redshift distribution. We pay special attention to the stability of this deconvolution, regularizing the method with a prior on the smoothness of the true redshift distribution. The methods that we develop are applicable to general photometric redshift surveys.

Heavily obscured AGN: an ideal laboratory to study the early co-evolution of galaxies and black holes

Obscured AGN are a crucial ingredient to understand the full growth history of super massive black holes and the coevolution with their host galaxies, since they constitute the bulk of the BH accretion. In the distant Universe, many of them are hosted by submillimeter galaxies (SMGs), characterized by a high production of stars and a very fast consumption of gas. Therefore, the analysis of this class of objects is fundamental to investigate the role of the ISM in the early coevolution of galaxies and black holes. We present a multiwavelength study of a sample of six obscured X-ray selected AGN at z>2.5 in the CDF-S, detected in the far-IR/submm bands. We performed the X-ray spectral analysis based on the 7Ms Chandra dataset, which provides the best X-ray spectral information currently available for distant AGN. We were able to place constraints on the obscuring column densities and the intrinsic luminosities of our targets. Moreover, we built up the UV to FIR spectral energy distributions (SEDs) by combining the broad-band photometry from CANDELS and the Herschel catalogs, and analyzed them by means of an SED decomposition technique. Therefore, we derived important physical parameters of both the host galaxy and the AGN. In addition, we obtained, through an empirical calibration, the gas mass in the host galaxy and assessed the galaxy sizes in order to estimate the column density associated with the host ISM. The comparison of the ISM column densities with the values measured from the X-ray spectral analysis pointed out that the contribution of the host ISM to the obscuration of the AGN emission can be substantial, ranging from ~10% up to ~100% of the value derived from the X-ray spectra. The absorption may occur at different physical scales in these sources and, in particular, the medium in the host galaxy is an ingredient that should be taken into account, since it may have a relevant role in driving the early co-evolution of galaxies with their black holes.

The relation between cool cluster cores and Herschel-detected star formation in brightest cluster galaxies

We present far-infrared (FIR) analysis of 68 brightest cluster galaxies (BCGs) at 0.08 < z < 1.0. Deriving total infrared luminosities directly from Spitzer and Herschel photometry spanning the peak of the dust component (24-500 μm), we calculate the obscured star formation rate (SFR). 22^+6.2 _–5.3% of the BCGs are detected in the far-infrared, with SFR = 1-150 M ☉ yr^–1. The infrared luminosity is highly correlated with cluster X-ray gas cooling times for cool-core clusters (gas cooling time <1 Gyr), strongly suggesting that the star formation in these BCGs is influenced by the cluster-scale cooling process. The occurrence of the molecular gas tracing Hα emission is also correlated with obscured star formation. For all but the most luminous BCGs (L_TIR > 2 × 10^11 L_☉), only a small (≤0.4 mag) reddening correction is required for SFR(Hα) to agree with SFR_FIR. The relatively low Hα extinction (dust obscuration), compared to values reported for the general star-forming population, lends further weight to an alternate (external) origin for the cold gas. Finally, we use a stacking analysis of non-cool-core clusters to show that the majority of the fuel for star formation in the FIR-bright BCGs is unlikely to originate from normal stellar mass loss.

Far-infrared constraints on the contamination by dust-obscured galaxies of high-z dropout searches

The spectral energy distributions (SED) of dusty galaxies at intermediate redshift may look similar to very high-redshift galaxies in the optical/near infrared (NIR) domain. This can lead to the contamination of high-redshift galaxy searches based on broad-band optical/NIR photometry by lower redshift dusty galaxies because both kind of galaxies cannot be distinguished. The contamination rate could be as high as 50%. This work shows how the far-infrared (FIR) domain can help to recognize likely low-z interlopers in an optical/NIR search for high-z galaxies. We analyze the FIR SEDs of two galaxies that are proposed to be very high-redshift (z > 7) dropout candidates based on deep Hawk-I/VLT observations. The FIR SEDs are sampled with PACS/Herschel at 100 and 160 μm, with SPIRE/Herschel at 250, 350 and 500 μm and with LABOCA/APEX at 870 μm. We find that redshifts > 7 would imply extreme FIR SEDs (with dust temperatures >100 K and FIR luminosities >10^13 L_⊙). At z ~ 2, instead, the SEDs of both sources would be compatible with those of typical ultra luminous infrared galaxies or submillimeter galaxies. Considering all available data for these sources from visible to FIR we re-estimate the redshifts and find z ~ 1.6–2.5. Owing to the strong spectral breaks observed in these galaxies, standard templates from the literature fail to reproduce the visible-to-near-IR part of the SEDs even when additional extinction is included. These sources strongly resemble dust-obscured galaxies selected in Spitzer observations with extreme visible-to-FIR colors, and the galaxy GN10 at z = 4. Galaxies with similar SEDs could contaminate other high-redshift surveys.

Determining star formation rates for infrared galaxies

We show that measures of star formation rates (SFRs) for infrared galaxies using either single-band 24 μm or extinction-corrected Paα luminosities are consistent in the total infrared luminosity = L(TIR) ~ 10^10 L_☉ range. MIPS 24 μm photometry can yield SFRs accurately from this luminosity upward: SFR(M_☉ yr^–1) = 7.8 × 10^–10 L(24 μm, L_☉) from L(TIR) = 5× 10^9 L_☉ to 10^11 L_☉ and SFR = 7.8 × 10^–10 L(24 μm, L_☉)(7.76 × 10^–11 L(24))^0.048 for higher L(TIR). For galaxies with L(TIR) ≥ 10^10 L_☉, these new expressions should provide SFRs to within 0.2 dex. For L(TIR) ≥ 10^11 L_☉, we find that the SFR of infrared galaxies is significantly underestimated using extinction-corrected Paα (and presumably using any other optical or near-infrared recombination lines). As a part of this work, we constructed spectral energy distribution templates for eleven luminous and ultraluminous purely star forming infrared galaxies and over the spectral range 0.4 μm to 30 cm. We use these templates and the SINGS data to construct average templates from 5 μm to 30 cm for infrared galaxies with L(TIR) = 5× 10^9 to 10^13 L_☉. All of these templates are made available online.

The star formation and extinction coevolution of UV-selected galaxies over 0.05 < z < 1.2

We use a new stacking technique to obtain mean mid-IR and far-IR to far-UV flux ratios over the rest-frame near-UV, near-IR color-magnitude diagram. We employ COMBO-17 redshifts and COMBO-17 optical, GALEX far- and near-UV, and Spitzer IRAC and MIPS mid-IR photometry. This technique permits us to probe the infrared excess (IRX), the ratio of far-IR to far-UV luminosity, and the specific star formation rate (SSFR) and their coevolution over 2 orders of magnitude of stellar mass and over redshift 0.1 < z < 1.2. We find that the SSFR and the characteristic mass (Script M_0) above which the SSFR drops increase with redshift (downsizing). At any given epoch, the IRX is an increasing function of mass up to Script M_0. Above this mass the IRX falls, suggesting gas exhaustion. In a given mass bin below Script M_0, the IRX increases with time in a fashion consistent with enrichment. We interpret these trends using a simple model with a Schmidt-Kennicutt law and extinction that tracks gas density and enrichment. We find that the average IRX and SSFR follow a galaxy age parameter ξ, which is determined mainly by the galaxy mass and time since formation. We conclude that blue-sequence galaxies have properties which show simple, systematic trends with mass and time such as the steady buildup of heavy elements in the interstellar media of evolving galaxies and the exhaustion of gas in galaxies that are evolving off the blue sequence. The IRX represents a tool for selecting galaxies at various stages of evolution.

The SCUBA HAlf degree extragalactic survey - III. Identification of radio and mid-infrared counterparts to submillimetre galaxies

Determining an accurate position for a submillimetre (submm) galaxy (SMG) is the crucial step that enables us to move from the basic properties of an SMG sample - source counts and 2D clustering - to an assessment of their detailed, multiwavelength properties, their contribution to the history of cosmic star formation and their links with present-day galaxy populations. In this paper, we identify robust radio and/or infrared (IR) counterparts, and hence accurate positions, for over two-thirds of the SCUBA HAlf-Degree Extragalactic Survey (SHADES) Source Catalogue, presenting optical, 24-μm and radio images of each SMG. Observed trends in identification rate have given no strong rationale for pruning the sample. Uncertainties in submm position are found to be consistent with theoretical expectations, with no evidence for significant additional sources of error. Employing the submm/radio redshift indicator, via a parametrization appropriate for radio-identified SMGs with spectroscopic redshifts, yields a median redshift of 2.8 for the radio-identified subset of SHADES, somewhat higher than the median spectroscopic redshift. We present a diagnostic colour-colour plot, exploiting Spitzer photometry, in which we identify regions commensurate with SMGs at very high redshift. Finally, we find that significantly more SMGs have multiple robust counterparts than would be expected by chance, indicative of physical associations. These multiple systems are most common amongst the brightest SMGs and are typically separated by 2-6 arcsec, similar to 15-20/sin i kpc at z~ 2, consistent with early bursts seen in merger simulations.

The Most Complete Photometric Analysis of 548 CALIFA Galaxies

We present an extensive photometric catalog for 548 CALIFA galaxies observed as of the summer of 2015. CALIFA is currently lacking photometry matching the scale and diversity of its spectroscopy; this work is intended to meet all photometric needs for CALIFA galaxies while also identifying best photometric practices for upcoming integral field spectroscopy surveys such as SAMI and MaNGA. This catalog comprises gri surface brightness profiles derived from Sloan Digital Sky Survey (SDSS) imaging, a variety of non-parametric quantities extracted from these pro files, and parametric models fitted to the i-band pro files (1D) and original galaxy images (2D). To compliment our photometric analysis, we contrast the relative performance of our 1D and 2D modelling approaches. The ability of each measurement to characterize the global properties of galaxies is quantitatively assessed, in the context of constructing the tightest scaling relations. Where possible, we compare our photometry with existing photometrically or spectroscopically obtained measurements from the literature. Close agreement is found with Walcher et al. (2014), the current source of basic photometry and classifications of CALIFA galaxies, while comparisons with spectroscopically derived quantities reveals the effect of CALIFA's limited field of view compared to broadband imaging surveys such as the SDSS. The colour-magnitude diagram, star formation main sequence, and Tully-Fisher relation of CALIFA galaxies are studied, to give a small example of the investigations possible with this rich catalog. We conclude with a discussion of points of concern for ongoing integral field spectroscopy surveys and directions for future expansion and exploitation of this work.

Sub-kiloparsec alma imaging of compact star-forming galaxies at z ~ 2.5: revealing the formation of dense galactic cores in the progenitors of compact quiescent galaxies

We present spatially resolved Atacama Large Millimeter/submillimeter Array (ALMA) 870 μm dust continuum maps of six massive, compact, dusty star-forming galaxies at z ~ 2.5. These galaxies are selected for their small rest-frame optical sizes (r_e,F160W ~ 1.6 kpc) and high stellar mass densities that suggest that they are direct progenitors of compact quiescent galaxies at z ~ 2. The deep observations yield high far-infrared (FIR) luminosities of L_IR = 10^12.3-12.8 L_⨀ and star formation rates (SFRs) of SFR = 200–700 M_⊙ yr^−1, consistent with those of typical star-forming "main sequence" galaxies. The high spatial resolution (FWHM ~ 0 12–0 18) ALMA and Hubble Space Telescope photometry are combined to construct deconvolved, mean radial profiles of their stellar mass and (UV+IR) SFR. We find that the dusty, nuclear IR–SFR overwhelmingly dominates the bolometric SFR up to r ~ 5 kpc, by a factor of over 100× from the unobscured UV–SFR. Furthermore, the effective radius of the mean SFR profile (r_e,SFR ~ 1 kpc) is ~30% smaller than that of the stellar mass profile. The implied structural evolution, if such nuclear starburst last for the estimated gas depletion time of Δt = ±100 Myr, is a 4×increase of the stellar mass density within the central 1 kpc and a 1.6× decrease of the half-mass–radius. This structural evolution fully supports dissipation-driven, formation scenarios in which strong nuclear starbursts transform larger, star-forming progenitors into compact quiescent galaxies.

Space density distribution of galaxies in the absolute magnitude - rotation velocity plane: a volume-complete Tully-Fisher relation from CALIFA stellar kinematics

We measured the distribution in absolute magnitude - circular velocity space for a well-defined sample of 199 rotating galaxies of the Calar Alto Legacy Integral Field Area Survey (CALIFA) using their stellar kinematics. Our aim in this analysis is to avoid subjective selection criteria and to take volume and large-scale structure factors into account. Using stellar velocity fields instead of gas emission line kinematics allows including rapidly rotating early-type galaxies. Our initial sample contains 277 galaxies with available stellar velocity fields and growth curve r-band photometry. After rejecting 51 velocity fields that could not be modelled because of the low number of bins, foreground contamination, or significant interaction, we performed Markov chain Monte Carlo modelling of the velocity fields, from which we obtained the rotation curve and kinematic parameters and their realistic uncertainties. We performed an extinction correction and calculated the circular velocity v_circ accounting for the pressure support of a given galaxy. The resulting galaxy distribution on the M-r - v(circ) plane was then modelled as a mixture of two distinct populations, allowing robust and reproducible rejection of outliers, a significant fraction of which are slow rotators. The selection effects are understood well enough that we were able to correct for the incompleteness of the sample. The 199 galaxies were weighted by volume and large-scale structure factors, which enabled us to fit a volume-corrected Tully-Fisher relation (TFR). More importantly, we also provide the volume-corrected distribution of galaxies in the M_r - v_circ plane, which can be compared with cosmological simulations. The joint distribution of the luminosity and circular velocity space densities, representative over the range of -20 > M_r > -22 mag, can place more stringent constraints on the galaxy formation and evolution scenarios than linear TFR fit parameters or the luminosity function alone.