Local luminous infrared galaxies: Spatially resolved mid-infrared observations with Spitzer/IRS

Luminous Infrared (IR) Galaxies (LIRGs, L_IR=10^11-10 L_⨀) are an important cosmological class of galaxies as they are the main contributors to the co-moving star formation rate density of the universe at z=1. In this paper we present a guaranteed time observation (GTO) Spitzer InfraRed Spectrograph (IRS) program aimed to obtain spectral mapping of a sample of 14 local d<76Mpc LIRGs. The data cubes map, at least, the central 20arcsec X 20arcsec to 30 arcsec X 30 arcsec regions of the galaxies, and use all four IRS modules covering the full 5-38 μ m spectral range. The final goal of this project is to characterize fully the mid-IR properties of local LIRGs as a first step to understanding their more distant counterparts. In this paper we present the first results of this GTO program. The IRS spectral mapping data allow us to build spectral maps of the bright mid-IR emission lines (e.g., [Ne II] 12.81 μ m, [Ne III]15.56 μ m, [S III] 18.71 μ m, H_2 at 17 μ m), continuum, the 6.2 and 11.3 μ m polycyclic aromatic hydrocarbon (PAH) features, and the 9.7 μ m silicate feature, as well as to extract 1D spectra for regions of interest in each galaxy. The IRS data are used to obtain spatially resolved measurements of the extinction using the 9.7 μ m silicate feature, and to trace star forming regions using the neon lines and the PAH features. We also investigate a number of active galactic nuclei (AGN) indicators, including the presence of high excitation emission lines and a strong dust continuum emission at around 6 9.7 μ m . We finally use the integrated Spitzer/IRS spectra as templates of local LIRGs. We discuss several possible uses for these templates, including the calibration of the star formation rate of IR-bright galaxies at high redshift. We also predict the intensities of the brightest mid-IR emission lines for LIRGs as a function of redshift, and compare them with the expected sensitivities of future space IR missions.

The host galaxies and black holes of typical z~0.5-1.4 AGNs

We study the stellar and star formation properties of the host galaxies of 58 X-ray-selected AGNs in the GOODS portion of the Chandra Deep Field South (CDF-S) region at z ~ 0.5-1.4. The AGNs are selected such that their rest-frame UV to near-infrared spectral energy distributions (SEDs) are dominated by stellar emission; i.e., they show a prominent 1.6 μm bump, thus minimizing the AGN emission "contamination." This AGN population comprises approximately 50% of the X-ray-selected AGNs at these redshifts. We find that AGNs reside in the most massive galaxies at the redshifts probed here. Their characteristic stellar masses (M_* ~ 7.8 × 10^10 and M_* ~ 1.2 × 10^11 M_☉ at median redshifts of 0.67 and 1.07, respectively) appear to be representative of the X-ray-selected AGN population at these redshifts and are intermediate between those of local type 2 AGNs and high-redshift (z ~ 2) AGNs. The inferred black hole masses (M_BH ~ 2 × 10^8 M_☉) of typical AGNs are similar to those of optically identified quasars at similar redshifts. Since the AGNs in our sample are much less luminous (L_2–10 keV < 10^44 erg s^−1) than quasars, typical AGNs have low Eddington ratios (η ~ 0.01-0.001). This suggests that, at least at intermediate redshifts, the cosmic AGN "downsizing" is due to both a decrease in the characteristic stellar mass of typical host galaxies and less efficient accretion. Finally, there is no strong evidence in AGN host galaxies for either highly suppressed star formation (expected if AGNs played a role in quenching star formation) or elevated star formation when compared to mass-selected (i.e., IRAC-selected) galaxies of similar stellar masses and redshifts.

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.

Near-infrared Instrumentation For Rapid-response Astronomy

Ɣ-ray bursts (GRBs) are the Universe's most luminous transient events. Since the discovery of GRBs was announced in 1973, efforts have been ongoing to obtain data over a broader range of the electromagnetic spectrum at the earliest possible times following the initial detection. The discovery of the theorized ``afterglow'' emission in radio through X-ray bands in the late 1990s confirmed the cosmological nature of these events. At present, GRB afterglows are among the best probes of the early Universe (z ≳ 9). In addition to informing theories about GRBs themselves, observations of afterglows probe the circum-burst medium (CBM), properties of the host galaxies and the progress of cosmic reionization. To explore the early-time variability of afterglows, I have developed a generalized analysis framework which models near-infrared (NIR), optical, ultra-violet (UV) and X-ray light curves without assuming an underlying model. These fits are then used to construct the spectral energy distribution (SED) of afterglows at arbitrary times within the observed window. Physical models are then used to explore the evolution of the SED parameter space with time. I demonstrate that this framework produces evidence of the photodestruction of dust in the CBM of GRB 120119A, similar to the findings from a previous study of this afterglow. The framework is additionally applied to the afterglows of GRB 140419A and GRB 080607. In these cases the evolution of the SEDs appears consistent with the standard fireball model. Having introduced the scientific motivations for early-time observations, I introduce the Rapid Infrared Imager-Spectrometer (RIMAS). Once commissioned on the 4.3 meter Discovery Channel Telescope (DCT), RIMAS will be used to study the afterglows of GRBs through photometric and spectroscopic observations beginning within minutes of the initial burst. The instrument will operate in the NIR, from 0.97 μm to 2.37 μm, permitting the detection of very high redshift (z ≳ 7) afterglows which are attenuated at shorter wavelengths by Lyman-ɑ absorption in the intergalactic medium (IGM). A majority of my graduate work has been spent designing and aligning RIMAS's cryogenic (~80 K) optical systems. Design efforts have included an original camera used to image the field surrounding spectroscopic slits, tolerancing and optimizing all of the instrument's optics, thermal modeling of optomechanical systems, and modeling the diffraction efficiencies for some of the dispersive elements. To align the cryogenic optics, I developed a procedure that was successfully used for a majority of the instrument's sub-assemblies. My work on this cryogenic instrument has necessitated experimental and computational projects to design and validate designs of several subsystems. Two of these projects describe simple and effective measurements of optomechanical components in vacuum and at cryogenic temperatures using an 8-bit CCD camera. Models of heat transfer via electrical harnesses used to provide current to motors located within the cryostat are also presented.

Faint Galaxies in the Megaparsec-Scale Environments of Hyperluminous QSOs at Redshifts 2 < z < 3

This thesis presents detailed observational studies of the extended distributions of gas, galaxies, and dark matter around hyperluminous quasars (HLQSOs) at high redshift. Taken together, these works aim to coherently describe the relationships between these massive, accreting black holes and their environments: the nature of the regions that give rise to such massive black holes, the effect of HLQSO radiation on their surrounding galaxies and gas, and the ability of both galaxies and black holes to shed new light on the formation and evolution of the other.

Chapter 2 focuses on the continuum-color-selected galaxies drawn from the Keck Baryonic Structure Survey (KBSS). The KBSS is a uniquely deep spectroscopic survey of star-forming galaxies in the same volumes of space as 15 HLQSOs at 2.5 < z < 2.9. The three-dimensional distribution of these galaxies among themselves and the nearby HLQSOs is used to infer the extent to which these black holes are associated with overdense peaks in the dark matter and galaxy distribution as quantified by clustering statistics. In conjunction with recent dark-matter simulations, these data provide the first estimates of the host dark-matter halo masses for HLQSOs, providing new insight into the formation and evolution of the most massive black holes at high redshift.

Chapter 3 describes the first results from a new survey (KBSS-Lyα) conducted for this thesis. The KBSS-Lyα survey uses narrowband imaging to identify Lyα-emitters (LAEs) in the ~Mpc regions around eight of the KBSS HLQSOs. Many of these LAEs show the effect of reprocessed HLQSO radiation in their emission through the process known as Lyα fluorescence. In this chapter, these fluorescent LAEs are used to generate a coarse map of the average HLQSO ionizing emission on Mpc scales, thereby setting the first direct constraints of the lifetime and angular distribution of activity for a population of these uniquely luminous black holes.

Chapter 4 contains a more detailed description of the KBSS-Lyα survey itself and the detailed properties of the star-forming and fluorescent objects selected therein. Using imaging and spectroscopic data covering rest-frame UV and optical wavelengths, including spectra from the new near-infrared spectrometer MOSFIRE, we characterize this population of nascent galaxies in terms of their kinematics, enrichment, gas properties, and luminosity distribution while comparing and contrasting them with previously-studied populations of continuum-selected galaxies and LAEs far from the effects of HLQSO emission.

At the conclusion of this thesis, I briefly present future directions for the continuation of this research. In Appendix A, I provide background information on the instrumentation used in this thesis, including my own contributions to MOSFIRE.

The J1030 field: a new window on early large scale structures and faint radio-galaxy populations

The correlations between the evolution of the Super Massive Black Holes (SMBHs) and their host galaxies suggests that the SMBH accretion on sub-pc scales (active galactice nuclei, AGN) is linked to the building of the galaxy over kpc scales, through the so called AGN feedback. Most of the galaxy assembly occurs in overdense large scale structures (LSSs). AGN residing in powerful sources in LSSs, such as the proto-brightest cluster galaxies (BCGs), can affect the evolution of the surrounding intra-cluster medium (ICM) and nearby galaxies. Among distant AGN, high-redshift radio-galaxies (HzRGs) are found to be excellent BCG progenitor candidates. In this Thesis we analyze novel interferometric observations of the so-called "J1030" field centered around the z = 6.3 SDSS Quasar J1030+0524, carried out with the Atacama large (sub-)millimetre array (ALMA) and the Jansky very large array (JVLA). This field host a LSS assembling around a powerful HzRG at z = 1.7 that shows evidence of positive AGN feedback in heating the surrounding ICM and promoting star-formation in multiple galaxies at hundreds kpc distances. We report the detection of gas-rich members of the LSS, including the HzRG. We showed that the LSS is going to evolve into a local massive cluster and the HzRG is the proto-BCG. we unveiled signatures of the proto-BCG's interaction with the surrounding ICM, strengthening the positive AGN feedback scenario. From the JVLA observations of the "J1030" we extracted one of the deepest extra-galactic radio surveys to date (~12.5 uJy at 5 sigma). Exploiting the synergy with the X-ray deep survey (~500 ks) we investigated the relation of the X-ray/radio emission of a X-ray-selected sample, unveiling that the radio emission is powered by different processes (star-formation and AGN), and that AGN-driven sample is mostly composed by radio-quiet objects that display a significant X-ray/radio correlation.

Hydrodynamic simulations of multiple stellar populations in globular clusters

The presence of multiple stellar populations in globular clusters (GCs) is now well accepted, however, very little is known regarding their origin. In this Thesis, I study how multiple populations formed and evolved by means of customized 3D numerical simulations, in light of the most recent data from spectroscopic and photometric observations of Local and high-redshift Universe. Numerical simulations are the perfect tool to interpret these data: hydrodynamic simulations are suited to study the early phases of GCs formation, to follow in great detail the gas behavior, while N-body codes permit tracing the stellar component. First, we study the formation of second-generation stars in a rotating massive GC. We assume that second-generation stars are formed out of asymptotic giant branch stars (AGBs) ejecta, diluted by external pristine gas. We find that, for low pristine gas density, stars mainly formed out of AGBs ejecta rotate faster than stars formed out of more diluted gas, in qualitative agreement with current observations. Then, assuming a similar setup, we explored whether Type Ia supernovae affect the second- generation star formation and their chemical composition. We show that the evolution depends on the density of the infalling gas, but, in general, an iron spread is developed, which may explain the spread observed in some massive GCs. Finally, we focused on the long-term evolution of a GC, composed of two populations and orbiting the Milky Way disk. We have derived that, for an extended first population and a low-mass second one, the cluster loses almost 98 percent of its initial first population mass and the GC mass can be as much as 20 times less after a Hubble time. Under these conditions, the derived fraction of second-population stars reproduces the observed value, which is one of the strongest constraints of GC mass loss.

Obscured AGN in deep fields: tracing accretion and star-formation up to z ∼ 5

The discovery of scaling relations between the mass of the SMBH and some key physical properties of the host galaxy suggests that the growth of the SMBH and that of the galaxy are coupled, with the AGN activity and the star-formation (SF) processes influencing each other. Although the mechanism of this co-evolution are still a matter of debate, all scenarios agree that a key phase of the co-evolution is represented by the obscured accretion phase. This phase is of the co-evolution is the least studied, mostly due to the challenge in detecting and recognizing such obscured AGN. My thesis aims at investigating the AGN-galaxy co-evolution paradigm by identifying and studying AGN in the obscured accretion phase. The study of obscured AGN is key for our understanding of the feedback processes and of the mutual influence of the SF and the AGN activity. Moreover, these obscured and elusive AGN are needed to explain the X-ray background spectrum and to reconcile the measurements and the theoretical prediction of the BH accretion rate density. In this thesis, we firstly investigate the synergies between IR and X-ray missions in detecting and characterizing AGN, with a particular focus on the most obscured ones. We exploited UV/optical emission lines to select high-redshift obscured AGN at the cosmic noon, where the highest SFR density and BH accretion rate density are expected. We provide X-ray spectral analysis and UV-to-far-IR SED-fitting. We show that our samples host a significant fraction of very obscured sources; many of these are highly accreting. Finally, we performe a thoughtful investigation of a galaxy at z~5 with unusual and peculiar features, that lead us to identify a second extremely young population of stars and hidden AGN activity.

The molecular gas content of star-forming galaxies in the ALPINE fields.

In this thesis, I aim to study the evolution with redshift of the gas mass fraction of a sample of 53 sources (from z ∼ 0.5 to z > 5) serendipitously detected in ALMA band 7 as part of the ALMA Large Program to INvestigate C II at Early Times (ALPINE). First, I used SED-fitting software CIGALE, which is able to implement energy balancing between the optical and the far infrared part, to produce a best-fit template of my sources and to have an estimate of some physical properties, such as the star formation rate (SFR), the total infrared luminosity and the total stellar mass. Then, using the tight correlation found by Scoville et al. (2014) between the ISM molecular gas mass and the rest-frame 850 μm luminosity, I used the latter, extrapolating it from the best-fit template using a code that I wrote in Python, as a tracer for the molecular gas. For my sample, I then derived the most important physical properties, such as molecular gas mass, gas mass fractions, specific star formation rate and depletion timescales, which allowed me to better categorize them and find them a place within the evolutionary history of the Universe. I also fitted our sources, via another code I wrote again in Python, with a general modified blackbody (MBB) model taken from the literature (Gilli et al. (2014), D’Amato et al. (2020)) to have a direct method of comparison with similar galaxies. What is evident at the end of the paper is that the methods used to derive the physical quantities of the sources are consistent with each other, and these in turn are in good agreement with what is found in the literature.

Non-thermal emission in High Frequency Peaked blazars towards the Square Kilometer Array era

In this thesis work we will explore and discuss the properties of the gamma-ray sources included in the first Fermi-LAT catalog of sources above 10 GeV (1FHL), by considering both blazars and the non negligible fraction of still unassociated gamma-ray sources (UGS, 13%). We perform a statistical analysis of a complete sample of hard gamma-ray sources, included in the 1FHL catalog, mostly composed of HSP blazars, and we present new VLBI observations of the faintest members of the sample. The new VLBI data, complemented by an extensive search of the archives for brighter sources, are essential to gather a sample as large as possible for the assessment of the significance of the correlation between radio and very high energy (E>100 GeV) emission bands. After the characterization of the statistical properties of HSP blazars and UGS, we use a complementary approach, by focusing on an intensive multi-frequency observing VLBI and gamma-ray campaign carried out for one of the most remarkable and closest HSP blazar Markarian 421.

Stacked searches for high-energy neutrinos from blazars with IceCube

In the year 2013, the detection of a diffuse astrophysical neutrino flux with the IceCube neutrino telescope – constructed at the geographic South Pole – was announced by the IceCube collaboration. However, the origin of these neutrinos is still unknown as no sources have been identified to this day. Promising neutrino source candidates are blazars, which are a subclass of active galactic nuclei with radio jets pointing towards the Earth. In this thesis, the neutrino flux from blazars is tested with a maximum likelihood stacking approach, analyzing the combined emission from uniform groups of objects. The stacking enhances the sensitivity w.r.t. the still unsuccessful single source searches. The analysis utilizes four years of IceCube data including one year from the completed detector. As all results presented in this work are compatible with background, upper limits on the neutrino flux are given. It is shown that, under certain conditions, some hadronic blazar models can be challenged or even rejected. Moreover, the sensitivity of this analysis – and any other future IceCube point source search – was enhanced by the development of a new angular reconstruction method. It is based on a detailed simulation of the photon propagation in the Antarctic ice. The median resolution for muon tracks, induced by high-energy neutrinos, is improved for all neutrino energies above IceCube’s lower threshold at 0.1TeV. By reprocessing the detector data and simulation from the year 2010, it is shown that the new method improves IceCube’s discovery potential by 20% to 30% depending on the declination.

The Fornax spectroscopic survey I. Survey strategy and preliminary results on the redshift distribution of a complete sample of stars and galaxies

The Fornax Spectroscopic Survey will use the Two degree Field spectrograph (2dF) of the Angle-Australian Telescope to obtain spectra for a complete sample of all 14000 objects with 16.5 less than or equal to b(j) less than or equal to 19.7 in a 12 square degree area centred on the Fornax Cluster. The aims of this project include the study of dwarf galaxies in the cluster (both known low surface brightness objects and putative normal surface brightness dwarfs) and a comparison sample of background field galaxies. We will also measure quasars and other active galaxies, any previously unrecognised compact galaxies and a large sample of Galactic stars. By selecting all objects-both stars and galaxies-independent of morphology, we cover a much larger range of surface brightness and scale size than previous surveys. In this paper we first describe the design of the survey. Our targets are selected from UK Schmidt Telescope sky survey plates digitised by the Automated Plate Measuring (APM) facility. We then describe the photometric and astrometric calibration of these data and show that the APM astrometry is accurate enough for use with the 2dF. We also describe a general approach to object identification using cross-correlations which allows us to identify and classify both stellar and galaxy spectra. We present results from the first 2dF field. Redshift distributions and velocity structures are shown for all observed objects in the direction of Fornax, including Galactic stars? galaxies in and around the Fornax Cluster, and for the background galaxy population. The velocity data for the stars show the contributions from the different Galactic components, plus a small tail to high velocities. We find no galaxies in the foreground to the cluster in our 2dF field. The Fornax Cluster is clearly defined kinematically. The mean velocity from the 26 cluster members having reliable redshifts is 1560 +/- 80 km s(-1). They show a velocity dispersion of 380 +/- 50 km s(-1). Large-scale structure can be traced behind the cluster to a redshift beyond z = 0.3. Background compact galaxies and low surface brightness galaxies are found to follow the general galaxy distribution.

Very-High-Energy gamma rays from a Distant Quasar: How Transparent Is the Universe?

The atmospheric Cherenkov gamma-ray telescope MAGIC, designed for a low-energy threshold, has detected very-high-energy gamma rays from a giant flare of the distant Quasi-Stellar Radio Source (in short: radio quasar) 3C 279, at a distance of more than 5 billion light-years (a redshift of 0.536). No quasar has been observed previously in very-high-energy gamma radiation, and this is also the most distant object detected emitting gamma rays above 50 gigaelectron volts. Because high-energy gamma rays may be stopped by interacting with the diffuse background light in the universe, the observations by MAGIC imply a low amount for such light, consistent with that known from galaxy counts.

Very High Energy Gamma-Rays emitting BL LAC's Population Study

Context: BL Lacs are the most numerous extragalactic objects which are detected in Very High Energy (VHE) gamma-rays band. They are a subclass of blazars. Large flux variability amplitude, sometimes happens in very short time scale, is a common characteristic of them. Significant optical polarization is another main characteristics of BL Lacs. BL Lacs' spectra have a continuous and featureless Spectral Energy Distribution (SED) which have two peaks. Among 1442 BL Lacs in the Roma-BZB catalogue, only 51 are detected in VHE gamma-rays band. BL Lacs are most numerous (more than 50% of 514 objects) objects among the sources that are detected above 10 GeV by FERMI-LAT. Therefore, many BL Lacs are expected to be discovered in VHE gamma-rays band. However, due to the limitation on current and near future technology of Imaging Air Cherenkov Telescope, astronomers are forced to predict whether an object emits VHE gamma-rays or not. Some VHE gamma-ray prediction methods are already introduced but still are not confirmed. Cross band correlations are the building blocks of introducing VHE gamma-rays prediction method. Aims: We will attempt to investigate cross band correlations between flux energy density, luminosity and spectral index of the sample. Also, we will check whether recently discovered MAGIC J2001+435 is a typical BL Lac. Methods: We select a sample of 42 TeV BL Lacs and collect 20 of their properties within five energy bands from literature and Tuorla blazar monitoring program database. All of the data are synchronized to be comparable to each other. Finally, we choose 55 pair of datasets for cross band correlations finding and investigating whether there is any correlation between each pair. For MAGIC J2001+435 we analyze the publicly available SWIFT-XRT data, and use the still unpublished VHE gamma-rays data from MAGIC collaboration. The results are compared to the other sources of the sample. Results: Low state luminosity of multiple detected VHE gamma-rays is strongly correlated luminosities in all other bands. However, the high state does not show such strong correlations. VHE gamma-rays single detected sources have similar behaviour to the low state of multiple detected ones. Finally, MAGIC J2001+435 is a typical TeV BL Lac. However, for some of the properties this source is located at the edge of the whole sample (e.g. in terms of X-rays flux). Keywords: BL Lac(s), Population study, Correlations finding, Multi wavelengths analysis, VHE gamma-rays, gamma-rays, X-rays, Optical, Radio

Transition redshift: new kinematic constraints from supernovae

The transition redshift (deceleration/acceleration) is discussed by expanding the deceleration parameter to first order around its present value. A detailed study is carried out by considering two different parametrizations, q = q(0) + q(1)z and q = q(0) + q(1)z(1 + z)(-1), and the associated free parameters (q(0), q(1)) are constrained by three different supernovae (SNe) samples. A previous analysis by Riess et al. using the first expansion is slightly improved and confirmed in light of their recent data (Gold07 sample). However, by fitting the model with the Supernova Legacy Survey (SNLS) type Ia sample, we find that the best fit to the redshift transition is z(t) = 0.61, instead of z(t) = 0.46 as derived by the High-z Supernovae Search (HZSNS) team. This result based in the SNLS sample is also in good agreement with the sample of Davis et al., z(t) = 0.60(-0.11)(+0.28) (1 sigma). Such results are in line with some independent analyses and accommodate more easily the concordance flat model (Lambda CDM). For both parametrizations, the three SNe Ia samples considered favour recent acceleration and past deceleration with a high degree of statistical confidence level. All the kinematic results presented here depend neither on the validity of general relativity nor on the matter-energy contents of the Universe.