Probing the 2D kinematic structure of early-type galaxies out to three effective radii

We detail an innovative new technique for measuring the two-dimensional (2D) velocity moments (rotation velocity, velocity dispersion and Gauss-Hermite coefficients h(3) and h(4)) of the stellar populations of galaxy haloes using spectra from Keck DEIMOS (Deep Imaging Multi-Object Spectrograph) multi-object spectroscopic observations. The data are used to reconstruct 2D rotation velocity maps. Here we present data for five nearby early-type galaxies to similar to three effective radii. We provide significant insights into the global kinematic structure of these galaxies, and challenge the accepted morphological classification in several cases. We show that between one and three effective radii the velocity dispersion declines very slowly, if at all, in all five galaxies. For the two galaxies with velocity dispersion profiles available from planetary nebulae data we find very good agreement with our stellar profiles. We find a variety of rotation profiles beyond one effective radius, i.e. rotation speed remaining constant, decreasing and increasing with radius. These results are of particular importance to studies which attempt to classify galaxies by their kinematic structure within one effective radius, such as the recent definition of fast- and slow-rotator classes by the Spectrographic Areal Unit for Research on Optical Nebulae project. Our data suggest that the rotator class may change when larger galactocentric radii are probed. This has important implications for dynamical modelling of early-type galaxies. The data from this study are available on-line.

The 6dF Galaxy Survey: stellar population trends across and through the Fundamental Plane

We present results from an analysis of stellar population parameters for 7132 galaxies in the 6dF Galaxy Survey Fundamental Plane (FP) sample. We bin the galaxies along the axes, v1, v2 and v3, of the tri-variate Gaussian to which we have fitted the galaxy distribution in effective radius, surface brightness and central velocity dispersion (FP space), and compute median values of stellar age, [Fe/H], [Z/H] and [a/Fe]. We determine the directions of the vectors in FP space along which each of the binned stellar population parameters vary most strongly. In contrast to previous work, we find stellar population trends not just with velocity dispersion and FP residual, but with radius and surface brightness as well. The most remarkable finding is that the stellar population parameters vary through the plane (v1 direction) and across the plane (v3 direction), but show no variation at all along the plane (v2 direction). The v2 direction in FP space roughly corresponds to luminosity density. We interpret a galaxys position along this vector as being closely tied to its merger history, such that early-type galaxies with lower luminosity density are more likely to have undergone major mergers. This conclusion is reinforced by an examination of the simulations of Kobayashi, which show clear trends of merger history with v2.

THE SLOAN DIGITAL SKY SURVEY CO-ADD: A GALAXY PHOTOMETRIC REDSHIFT CATALOG

We present and describe a catalog of galaxy photometric redshifts (photo-z) for the Sloan Digital Sky Survey (SDSS) Co-add Data. We use the artificial neural network (ANN) technique to calculate the photo-z and the nearest neighbor error method to estimate photo-z errors for similar to 13 million objects classified as galaxies in the co-add with r < 24.5. The photo-z and photo-z error estimators are trained and validated on a sample of similar to 83,000 galaxies that have SDSS photometry and spectroscopic redshifts measured by the SDSS Data Release 7 (DR7), the Canadian Network for Observational Cosmology Field Galaxy Survey, the Deep Extragalactic Evolutionary Probe Data Release 3, the VIsible imaging Multi-Object Spectrograph-Very Large Telescope Deep Survey, and the WiggleZ Dark Energy Survey. For the best ANN methods we have tried, we find that 68% of the galaxies in the validation set have a photo-z error smaller than sigma(68) = 0.031. After presenting our results and quality tests, we provide a short guide for users accessing the public data.

Development of a new galaxy classification technique and its application to the zCOSMOS survey

The purpose of this Thesis is to develop a robust and powerful method to classify galaxies from large surveys, in order to establish and confirm the connections between the principal observational parameters of the galaxies (spectral features, colours, morphological indices), and help unveil the evolution of these parameters from $z \sim 1$ to the local Universe. Within the framework of zCOSMOS-bright survey, and making use of its large database of objects ($\sim 10\,000$ galaxies in the redshift range $0 < z \lesssim 1.2$) and its great reliability in redshift and spectral properties determinations, first we adopt and extend the \emph{classification cube method}, as developed by Mignoli et al. (2009), to exploit the bimodal properties of galaxies (spectral, photometric and morphologic) separately, and then combining together these three subclassifications. We use this classification method as a test for a newly devised statistical classification, based on Principal Component Analysis and Unsupervised Fuzzy Partition clustering method (PCA+UFP), which is able to define the galaxy population exploiting their natural global bimodality, considering simultaneously up to 8 different properties. The PCA+UFP analysis is a very powerful and robust tool to probe the nature and the evolution of galaxies in a survey. It allows to define with less uncertainties the classification of galaxies, adding the flexibility to be adapted to different parameters: being a fuzzy classification it avoids the problems due to a hard classification, such as the classification cube presented in the first part of the article. The PCA+UFP method can be easily applied to different datasets: it does not rely on the nature of the data and for this reason it can be successfully employed with others observables (magnitudes, colours) or derived properties (masses, luminosities, SFRs, etc.). The agreement between the two classification cluster definitions is very high. ``Early'' and ``late'' type galaxies are well defined by the spectral, photometric and morphological properties, both considering them in a separate way and then combining the classifications (classification cube) and treating them as a whole (PCA+UFP cluster analysis). Differences arise in the definition of outliers: the classification cube is much more sensitive to single measurement errors or misclassifications in one property than the PCA+UFP cluster analysis, in which errors are ``averaged out'' during the process. This method allowed us to behold the \emph{downsizing} effect taking place in the PC spaces: the migration between the blue cloud towards the red clump happens at higher redshifts for galaxies of larger mass. The determination of $M_{\mathrm{cross}}$ the transition mass is in significant agreement with others values in literature.

The dynamics of early-type galaxies as a tool to understand their hot coronae and their IMF

Dynamical models of galaxies are a powerful tool to study and understand several astrophysical problems related to galaxy formation and evolution. This thesis is focussed on a particular type of dynamical models, that are widely used in literature, and are based on the solution of the Jeans equations. By means of a numerical Jeans solver code, developed on purpose and able to build state-of-the-art advanced axisymmetric galaxy models, two of the main currently investigated issues in the field of research of early-type galaxies (ETGs) are addressed. The first topic concerns the hot and X-ray emitting gaseous coronae that surround ETGs. The main goal is to explain why flat and rotating galaxies generally exhibit haloes with lower gas temperatures and luminosities with respect to rounder and velocity dispersion supported systems. The second astrophysical problem addressed concerns instead the stellar initial mass function (IMF) of ETGs. Nowadays, this is a very controversial issue due to a growing number of works on ETGs, based on different and independent techniques, that show evidences of a systematic variation of the IMF normalization as a function of galaxy velocity dispersion or mass. These studies are changing the previous opinion that the IMF of ETGs was the same as that of spiral galaxies, and hence universal throughout the whole large family of galaxies.

Extracting evolutionary information from the spectral decomposition of early-type galaxies.

Holding the major share of stellar mass in galaxies and being also old and passively evolving, early-type galaxies (ETGs) are the primary probes in investigating these various evolution scenarios, as well as being useful means to provide insights on cosmological parameters. In this thesis work I focused specifically on ETGs and on their capability in constraining galaxy formation and evolution; in particular, the principal aims were to derive some of the ETGs evolutionary parameters, such as age, metallicity and star formation history (SFH) and to study their age-redshift and mass-age relations. In order to infer galaxy physical parameters, I used the public code STARLIGHT: this program provides a best fit to the observed spectrum from a combination of many theoretical models defined in user-made libraries. the comparison between the output and input light-weighted ages shows a good agreement starting from SNRs of ∼ 10, with a bias of ∼ 2.2% and a dispersion 3%. Furthermore, also metallicities and SFHs are well reproduced. In the second part of the thesis I performed an analysis on real data, starting from Sloan Digital Sky Survey (SDSS) spectra. I found that galaxies get older with cosmic time and with increasing mass (for a fixed redshift bin); absolute light-weighted ages, instead, result independent from the fitting parameters or the synthetic models used. Metallicities, instead, are very similar from each other and clearly consistent with the ones derived from the Lick indices. The predicted SFH indicates the presence of a double burst of star formation. Velocity dispersions and extinctiona are also well constrained, following the expected behaviours. As a further step, I also fitted single SDSS spectra (with SNR∼ 20), to verify that stacked spectra gave the same results without introducing any bias: this is an important check, if one wants to apply the method at higher z, where stacked spectra are necessary to increase the SNR. Our upcoming aim is to adopt this approach also on galaxy spectra obtained from higher redshift Surveys, such as BOSS (z ∼ 0.5), zCOSMOS (z 1), K20 (z ∼ 1), GMASS (z ∼ 1.5) and, eventually, Euclid (z 2). Indeed, I am currently carrying on a preliminary study to estabilish the applicability of the method to lower resolution, as well as higher redshift (z 2) spectra, just like the Euclid ones.

On the Luminous and Dark Matter Distribution in Early-Type Galaxies

The way mass is distributed in galaxies plays a major role in shaping their evolution across cosmic time. The galaxy's total mass is usually determined by tracing the motion of stars in its potential, which can be probed observationally by measuring stellar spectra at different distances from the galactic centre, whose kinematics is used to constrain dynamical models. A class of such models, commonly used to accurately determine the distribution of luminous and dark matter in galaxies, is that of equilibrium models. In this Thesis, a novel approach to the design of equilibrium dynamical models, in which the distribution function is an analytic function of the action integrals, is presented. Axisymmetric and rotating models are used to explain observations of a sample of nearby early-type galaxies in the Calar Alto Legacy Integral Field Area survey. Photometric and spectroscopic data for round and flattened galaxies are well fitted by the models, which are then used to get the galaxies' total mass distribution and orbital anisotropy. The time evolution of massive early-type galaxies is also investigated with numerical models. Their structural properties (mass, size, velocity dispersion) are observed to evolve, on average, with redshift. In particular, they appear to be significantly more compact at higher redshift, at fixed stellar mass, so it is interesting to investigate what drives such evolution. This Thesis focuses on the role played by dark-matter haloes: their mass-size and mass-velocity dispersion correlations evolve similarly to the analogous correlations of ellipticals; at fixed halo mass, the haloes are more compact at higher redshift, similarly to massive galaxies; a simple model, in which all the galaxy's size and velocity-dispersion evolution is due to the cosmological evolution of the underlying halo population, reproduces the observed size and velocity-dispersion of massive compact early-type galaxies up to redshift of about 2.

Observing the epoch of galaxy formation

Significant observational progress in addressing the question of the origin and early evolution of galaxies has been made in the past few years, allowing for direct comparison of the epoch when most of the stars in the universe were forming to prevailing theoretical models. There is currently broad consistency between theoretical expectations and the observations, but rapid improvement in the data will provide much more critical tests of theory in the coming years.

Galaxy formation

It is argued that within the standard Big Bang cosmological model the bulk of the mass of the luminous parts of the large galaxies likely had been assembled by redshift z ∼ 10. Galaxy assembly this early would be difficult to fit in the widely discussed adiabatic cold dark matter model for structure formation, but it could agree with an isocurvature version in which the cold dark matter is the remnant of a massive scalar field frozen (or squeezed) from quantum fluctuations during inflation. The squeezed field fluctuations would be Gaussian with zero mean, and the distribution of the field mass therefore would be the square of a random Gaussian process. This offers a possibly interesting new direction for the numerical exploration of models for cosmic structure formation.

The ionized gas in the CALIFA early-type galaxies I. Mapping two representative cases: NGC 6762 and NGC 5966

As part of the ongoing CALIFA survey, we have conducted a thorough bidimensional analysis of the ionized gas in two E/S0 galaxies, NGC 6762 and NGC 5966, aiming to shed light on the nature of their warm ionized ISM. Specifically, we present optical (3745–7300 Å) integral field spectroscopy obtained with the PMAS/PPAK integral field spectrophotometer. Its wide field-of-view (1′ × 1′) covers the entire optical extent of each galaxy down to faint continuum surface brightnesses. To recover the nebular lines, we modeled and subtracted the underlying stellar continuum from the observed spectra using the STARLIGHT spectral synthesis code. The pure emission-line spectra were used to investigate the gas properties and determine the possible sources of ionization. We show the advantages of IFU data in interpreting the complex nature of the ionized gas in NGC 6762 and NGC 5966. In NGC 6762, the ionized gas and stellar emission display similar morphologies, while the emission line morphology is elongated in NGC 5966, spanning ~6 kpc, and is oriented roughly orthogonal to the major axis of the stellar continuum ellipsoid. Whereas gas and stars are kinematically aligned in NGC 6762, the gas is kinematically decoupled from the stars in NGC 5966. A decoupled rotating disk or an “ionization cone” are two possible interpretations of the elongated ionized gas structure in NGC 5966. The latter would be the first “ionization cone” of such a dimension detected within a weak emission-line galaxy. Both galaxies have weak emission-lines relative to the continuum[EW(Hα) ≲ 3 Å] and have very low excitation, log([OIII]λ5007/Hβ) ≲ 0.5. Based on optical diagnostic ratios ([OIII]λ5007/Hβ, [NII]λ6584/Hα, [SII]λ6717, 6731/Hα, [OI]λ6300/Hα), both objects contain a LINER nucleus and an extended LINER-like gas emission. The emission line ratios do not vary significantly with radius or aperture, which indicates that the nebular properties are spatially homogeneous. The gas emission in NGC 6762 can be best explained by photoionization by pAGB stars without the need of invoking any other excitation mechanism. In the case of NGC 5966, the presence of a nuclear ionizing source seems to be required to shape the elongated gas emission feature in the “ionization cone” scenario, although ionization by pAGB stars cannot be ruled out. Further study of this object is needed to clarify the nature of its elongated gas structure.

Caught in the act: gas and stellar velocity dispersions in a fast Quenching compact star-forming galaxy at z ~ 1.7

We present Keck I MOSFIRE spectroscopy in the Y and H bands of GDN-8231, a massive, compact, star-forming galaxy at a redshift of z ~ 1.7. Its spectrum reveals both Hα and [Nii] emission lines and strong Balmer absorption lines. The Hα and Spitzer MIPS 24 μm fluxes are both weak, thus indicating a low star-formation rate of SFR≲5-10 M_⨀ yr−1. This, added to a relatively young age of ~700 Myr measured from the absorption lines, provides the first direct evidence for a distant galaxy being caught in the act of rapidly shutting down its star formation. Such quenching allows GDN-8231 to become a compact, quiescent galaxy, similar to three other galaxies in our sample, by z ~ 1.5. Moreover, the color profile of GDN-8231 shows a bluer center, consistent with the predictions of recent simulations for an early phase of inside-out quenching. Its line-of-sight velocity dispersion for the gas, σ_LOG^gas = 127 ± 32 km s^−1, is nearly 40% smaller than that of its stars, σ_LOG^* = 215 ± 35 km s^−1. High-resolution hydro-simulations of galaxies explain such apparently colder gas kinematics of up to a factor of ~1.5 with rotating disks being viewed at different inclinations and/or centrally concentrated star-forming regions. A clear prediction is that their compact, quiescent descendants preserve some remnant rotation from their star-forming progenitors.

A new algorithm for the detection of intercluster galaxy filaments using galaxy orientation alignments

We present a new algorithm for detecting intercluster galaxy filaments based upon the assumption that the orientations of constituent galaxies along such filaments are non-isotropic. We apply the algorithm to the 2dF Galaxy Redshift Survey catalogue and find that it readily detects many straight filaments between close cluster pairs. At large intercluster separations (> 15 h(-1) Mpc), we find that the detection efficiency falls quickly, as it also does with more complex filament morphologies. We explore the underlying assumptions and suggest that it is only in the case of close cluster pairs that we can expect galaxy orientations to be significantly correlated with filament direction.

NGC922 - a New Drop-through Ring Galaxy

We have found the peculiar galaxy NGC 922 to be a new drop-through ring galaxy using multiwavelength (ultraviolet-radio) imaging and spectroscopic observations. Its 'C'-shaped morphology and tidal plume indicate a recent strong interaction with its companion which was identified with these observations. Using numerical simulations we demonstrate that the main properties of the system can be generated by a high-speed off-axis drop-through collision of a small galaxy with a larger disc system, thus making NGC 922 one of the nearest known collisional ring galaxies. While these systems are rare in the local Universe, recent deep Hubble Space Telescope images suggest they were more common in the early Universe.

IRAS 00183-7111: ALMA and X-ray view of an ultra luminous infrared galaxy

The aim of this Thesis work is to study the multi-frequency properties of the Ultra Luminous Infrared Galaxy (ULIRG) IRAS 00183-7111 (I00183) at z = 0.327, connecting ALMA sub-mm/mm observations with those at high energies in order to place constraints on the properties of its central power source and verify whether the gas traced by the CO may be responsible for the obscuration observed in X-rays. I00183 was selected from the so-called Spoon diagnostic diagram (Spoon et al. 2007) for mid-infrared spectra of infrared galaxies based on the equivalent width of the 6.2 μm Polycyclic Aromatic Hydrocarbon (PAH) emission feature versus the 9.7 μm silicate strength. Such features are a powerful tool to investigate the contribution of star formation and AGN activity in this class of objects. I00183 was selected from the top-left region of the plot where the most obscured sources, characterized by a strong Si absorption feature, are located. To link the sub-mm/mm to the X-ray properties of I00183, ALMA archival Cycle 0 data in Band 3 (87 GHz) and Band 6 (270 GHz) have been calibrated and analyzed, using CASA software. ALMA Cycle 0 was the Early Science program for which data reprocessing is strongly suggested. The main work of this Thesis consisted in reprocessing raw data to provide an improvement with respect to the available archival products and results, which were obtained using standard procedures. The high-energy data consists of Chandra, XMM-Newton and NuSTAR observations which provide a broad coverage of the spectrum in the energy range 0.5 − 30 keV. Chandra and XMM archival data were used, with an exposure time of 22 and 22.2 ks, respectively; their reduction was carried out using CIAO and SAS software. The 100 ks NuSTAR are still private and the spectra were obtained by courtesy of the PI (K. Iwasawa). A detailed spectral analysis was done using XSPEC software; the spectral shape was reproduced starting from simple phenomenological models, and then more physical models were introduced to account for the complex mechanisms that involve this source. In Chapter 1, an overview of the scientific background is discussed, with a focus on the target, I00183, and the Spoon diagnostic diagram, from which it was originally selected. In Chapter 2, the basic principles of interferometry are briefly introduced, with a description of the calibration theory applied to interferometric observations. In Chapter 3, ALMA and its capabilities, both current and future, are shown, explaining also the complex structure of the ALMA archive. In Chapter 4, the calibration of ALMA data is presented and discussed, showing also the obtained imaging products. In Chapter 5, the analysis and discussion of the main results obtained from ALMA data is presented. In Chapter 6, the X-ray observations, data reduction and spectral analysis are reported, with a brief introduction to the basic principle of X-ray astronomy and the instruments from which the observations were carried out. Finally, the overall work is summarized, with particular emphasis on the main obtained results and the possible future perspectives.

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.