Update on the correlation of the highest energy cosmic rays with nearby extragalactic matter

Data collected by the Pierre Auger Observatory through 31 August 2007 showed evidence for anisotropy in the arrival directions of cosmic rays above the Greisen-Zatsepin-Kuz'min energy threshold, 6 x 10(19) eV. The anisotropy was measured by the fraction of arrival directions that are less than 3.1 degrees from the position of an active galactic nucleus within 75 Mpc (using the Veron-Cetty and Veron 12th catalog). An updated measurement of this fraction is reported here using the arrival directions of cosmic rays recorded above the same energy threshold through 31 December 2009. The number of arrival directions has increased from 27 to 69, allowing a more precise measurement. The correlating fraction is (38(-6)(+7))%, compared with 21% expected for isotropic cosmic rays. This is down from the early estimate of (69-(+11)(13))%. The enlarged set of arrival directions is examined also in relation to other populations of nearby extragalactic objects: galaxies in the 2 Microns All Sky Survey and active galactic nuclei detected in hard X-rays by the Swift Burst Alert Telescope. A celestial region around the position of the radiogalaxy Cen A has the largest excess of arrival directions relative to isotropic expectations. The 2-point autocorrelation function is shown for the enlarged set of arrival directions and compared to the isotropic expectation. (C) 2010 Elsevier B.V. All rights reserved.

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.

X-ray analysis of local mid-infrared selected Compton-thick AGN candidates

Cosmic X-ray background synthesis models (Gilli 2007) require a significant fraction of obscured AGN, some of which are expected to be heavily obscured (Compton-thick), but the number density of observationally found obscured sources is still an open issue (Vignali 2010, 2014). This thesis work takes advantage of recent NuSTAR data and is based on a multiwavelength research approach. Gruppioni et al. 2016 compared the AGN bolometric luminosity, for a sample of local 12 micron Seyfert galaxies, derived from the SED decomposition to the same quantity obtained by the 2-10 keV luminosity (IPAC-NED). A difference up to two orders of magnitude resulted between these quantities for some sources. Thus, the intrinsic X-ray luminosity obtained correcting for the obscuration may be underestimated. In this thesis we have tested this hypothesis by re-analysing the X-ray spectra of three of the sources (UGC05101, NGC1194 and NGC3079), for which observations from NuSTAR and Chandra and/or XMM-Newton were available. This is meant to extend our analysis to energies above 10 keV and thus estimate the AGN column density as reliable as possible. For spectral fitting we made use of both the commonly used XSPEC package and the two very recent MYtorus and BNtorus physical models. The available wide bandpass allowed us to achieve new and more solid insights into the X-ray spectral properties of these sources. The measured absorption column densities are highly suggestive of heavy obscuration. Once corrected the X-ray AGN luminosity for the obscuration estimated through our spectral analysis, we compared the L(X) values in the 2-10 keV band with those derived from the MIR band, by means of the relation by Gandhi, 2009. As expected, the values derived from this relation are in good agreement with those we measured, indicating that the column densities were underestimated in the previous literature works.

The Califa and Hipass circular velocity function for all morphological galaxy types

The velocity function (VF) is a fundamental observable statistic of the galaxy population that is similar to the luminosity function in importance, but much more difficult to measure. In this work we present the first directly measured circular VF that is representative between 60 < v_circ < 320 km s^-1 for galaxies of all morphological types at a given rotation velocity. For the low-mass galaxy population (60 < v_circ < 170 km s^-1), we use the HI Parkes All Sky Survey VF. For the massive galaxy population (170 < v_circ < 320 km s^-1), we use stellar circular velocities from the Calar Alto Legacy Integral Field Area Survey (CALIFA). In earlier work we obtained the measurements of circular velocity at the 80% light radius for 226 galaxies and demonstrated that the CALIFA sample can produce volume-corrected galaxy distribution functions. The CALIFA VF includes homogeneous velocity measurements of both late and early-type rotation-supported galaxies and has the crucial advantage of not missing gas-poor massive ellipticals that HI surveys are blind to. We show that both VFs can be combined in a seamless manner, as their ranges of validity overlap. The resulting observed VF is compared to VFs derived from cosmological simulations of the z = 0 galaxy population. We find that dark-matter-only simulations show a strong mismatch with the observed VF. Hydrodynamic simulations fare better, but still do not fully reproduce observations.

Mass assembly and morphological transformations since z ∼ 3 from CANDELS

We quantify the evolution of the stellar mass functions (SMFs) of star-forming and quiescent galaxies as a function of morphology from z ∼ 3 to the present. Our sample consists of ∼50 000 galaxies in the CANDELS fields (∼880 arcmin^2), which we divide into four main morphological types, i.e. pure bulge-dominated systems, pure spiral disc-dominated, intermediate two-component bulge+disc systems and irregular disturbed galaxies. At z ∼ 2, 80 per cent of the stellar mass density of star-forming galaxies is in irregular systems. However, by z ∼ 0.5, irregular objects only dominate at stellar masses below 10^9 M_⊙. A majority of the star-forming irregulars present at z ∼ 2 undergo a gradual transformation from disturbed to normal spiral disc morphologies by z ∼ 1 without significant interruption to their star formation. Rejuvenation after a quenching event does not seem to be common except perhaps for the most massive objects, because the fraction of bulge-dominated star-forming galaxies with M^*/M_⊙ > 10^10.7 reaches 40 per cent at z < 1. Quenching implies the presence of a bulge: the abundance of massive red discs is negligible at all redshifts over 2 dex in stellar mass. However, the dominant quenching mechanism evolves. At z > 2, the SMF of quiescent galaxies above M^* is dominated by compact spheroids. Quenching at this early epoch destroys the disc and produces a compact remnant unless the star-forming progenitors at even higher redshifts are significantly more dense. At 1 < z < 2, the majority of newly quenched galaxies are discs with a significant central bulge. This suggests that mass quenching at this epoch starts from the inner parts and preserves the disc. At z < 1, the high-mass end of the passive SMF is globally in place and the evolution mostly happens at stellar masses below 10^10 M_⊙. These low-mass galaxies are compact, bulge-dominated systems, which were environmentally quenched: destruction of the disc through ram-pressure stripping is the likely process.

SAFARI optical system architecture and design concept

SpicA FAR infrared Instrument, SAFARI, is one of the instruments planned for the SPICA mission. The SPICA mission is the next great leap forward in space-based far-infrared astronomy and will study the evolution of galaxies, stars and planetary systems. SPICA will utilize a deeply cooled 2.5m-class telescope, provided by European industry, to realize zodiacal background limited performance, and high spatial resolution. The instrument SAFARI is a cryogenic grating-based point source spectrometer working in the wavelength domain 34 to 230 μm, providing spectral resolving power from 300 to at least 2000. The instrument shall provide low and high resolution spectroscopy in four spectral bands. Low Resolution mode is the native instrument mode, while the high Resolution mode is achieved by means of a Martin-Pupplet interferometer. The optical system is all-reflective and consists of three main modules; an input optics module, followed by the Band and Mode Distributing Optics and the grating Modules. The instrument utilizes Nyquist sampled filled linear arrays of very sensitive TES detectors. The work presented in this paper describes the optical design architecture and design concept compatible with the current instrument performance and volume design drivers.

Modeling chemistry in massive star forming regions with internal PDRs

Over the past decades star formation has been a very attractive field because knowledge of star formation leads to a better understanding of the formation of planets and thus of our solar system but also of the evolution of galaxies. Conditions leading to the formation of high-mass stars are still under investigation but an evolutionary scenario has been proposed: As a cold pre-stellar core collapses under gravitational force, the medium warms up until it reaches a temperature of 100 K and enters the hot molecular core (HMC) phase. The forming central proto-star accretes materials, increasing its mass and luminosity and eventually it becomes sufficiently evolved to emit UV photons which irradiate the surrounding environment forming a hyper compact (HC) and then a ultracompact (UC) HII region. At this stage, a very dense and very thin internal photon-dominated region (PDR) forms between the HII region and the molecular core. Information on the chemistry allows to trace the physical processes occurring in these different phases of star formation. Formation and destruction routes of molecules are influenced by the environment as reaction rates depend on the temperature and radiation field. Therefore, chemistry also allows the determination of the evolutionary stage of astrophysical objects through the use of chemical models including the time evolution of the temperature and radiation field. Because HMCs host a very rich chemistry with high abundances of complex organic molecules (COMs), several astrochemical models have been developed to study the gas phase chemistry as well as grain chemistry in these regions. In addition to HMCs models, models of PDRs have also been developed to study in particular photo-chemistry. So far, few studies have investigated internal PDRs and only in the presence of outflows cavities. Thus, these unique regions around HC/UCHII regions remain to be examined thoroughly. My PhD thesis focuses on the spatio-temporal chemical evolution in HC/UC HII regions with internal PDRs as well as in HMCs. The purpose of this study is first to understand the impact and effects of the radiation field, usually very strong in these regions, on the chemistry. Secondly, the goal is to study the emission of various tracers of HC/UCHII regions and compare it with HMCs models, where the UV radiation field does not impact the region as it is immediately attenuated by the medium. Ultimately we want to determine the age of a given region using chemistry in combination with radiative transfer.

PRISM (Polarized Radiation Imaging and Spectroscopy Mission): an extended white paper

PRISM (Polarized Radiation Imaging and Spectroscopy Mission) was proposed to ESA in May 2013 as a large-class mission for investigating within the framework of the ESA Cosmic Vision program a set of important scientific questions that require high res- olution, high sensitivity, full-sky observations of the sky emission at wavelengths ranging from millimeter-wave to the far-infrared. PRISM’s main objective is to explore the distant universe, probing cosmic history from very early times until now as well as the structures, distribution of matter, and velocity flows throughout our Hubble volume. PRISM will survey the full sky in a large number of frequency bands in both intensity and polarization and will measure the absolute spectrum of sky emission more than three orders of magnitude bet- ter than COBE FIRAS. The data obtained will allow us to precisely measure the absolute sky brightness and polarization of all the components of the sky emission in the observed frequency range, separating the primordial and extragalactic components cleanly from the galactic and zodiacal light emissions. The aim of this Extended White Paper is to provide a more detailed overview of the highlights of the new science that will be made possible by PRISM, which include: (1) the ultimate galaxy cluster survey using the Sunyaev-Zeldovich (SZ) e↵ect, detecting approximately 106 clusters extending to large redshift, including a char- acterization of the gas temperature of the brightest ones (through the relativistic corrections to the classic SZ template) as well as a peculiar velocity survey using the kinetic SZ e↵ect that comprises our entire Hubble volume; (2) a detailed characterization of the properties and evolution of dusty galaxies, where the most of the star formation in the universe took place, the faintest population of which constitute the di↵use CIB (Cosmic Infrared Background); (3) a characterization of the B modes from primordial gravity waves generated during inflation and from gravitational lensing, as well as the ultimate search for primordial non-Gaussianity using CMB polarization, which is less contaminated by foregrounds on small scales than thetemperature anisotropies; (4) a search for distortions from a perfect blackbody spectrum, which include some nearly certain signals and others that are more speculative but more informative; and (5) a study of the role of the magnetic field in star formation and its inter- action with other components of the interstellar medium of our Galaxy. These are but a few of the highlights presented here along with a description of the proposed instrument.

Analysis of the high-resolution ro-vibrational spectrum of DC3N in the far and mid infrared regions

Cyanoacetylene HC3N is a molecule of great astronomical importance and it has been observed in many interstellar environments. Its deuterated form DC3N has been detected in number of sources from external galaxies to Galactic interstellar clouds, star-forming regions and planetary atmospheres. All these detections relied on previous laboratory investigations, which however still lack some essential information concerning its infrared spectrum. In this project, high-resolution ro-vibrational spectra of DC3N have been recorded in two energy regions: 150 – 450 cm-1 and 1800 – 2800 cm-1. In the first window the ν7← GS, 2ν7 ← ν7, ν5 ← ν7, ν5+ν7 ← 2ν7, ν6+ν7 → 2v7, 4ν7 ← 2ν7 bands have been assigned, while in the second region the three stretching fundamental bands ν1, ν2, ν3 have been observed and analysed. The 150 – 450 cm-1 region spectra have been recorded at the AILES beamline at the SOLEIL synchrotron (France), the 1800 – 2800 cm-1 spectra at the Department of Industrial Chemistry “Toso Montanari” in Bologna. In total, 2299 transitions have been assigned. Such experimental transition, together with data previously recorded for DC3N, were included in a least-squares fitting procedure from which several spectroscopic parameters have been determined with high precision and accuracy. They include rotational, vibrational and resonance constants. The spectroscopic data of DC3N have been included in a line catalog for this molecule in order to assist future astronomical observations and data interpretation. A paper which includes this research work has been published (M. Melosso, L. Bizzocchi, A. Adamczyk, E. Cane, P. Caselli, L. Colzid, L. Dorea, B. M. Giulianob, J.-C. Guillemine, M-A. Martin-Drumel, O. Piralif, A. Pietropolli Charmet , D. Prudenzano, V. M. Rivillad, F. Tamassia, Extensive ro-vibrational analysis of deuterated-cyanoacetylene (DC3N) from millimeter wavelengths to the infrared domain, Jour. of Quant. Spectr. and Rad. Tran. 254, 107221, 2020).

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.

Chemical composition of Milky Way satellites: Magellanic Clouds and Sagittarius dwarf galaxy

This PhD project is aimed at investigating the chemical composition of the stellar populations in the closest satellites of the Milky Way (MW), namely the Large and Small Magellanic Cloud (LMC and SMC, respectively) and the remnant of the Sagittarius (Sgr) dwarf spheroidal galaxy. Their proximity allows us to resolve their individual stars both with spectroscopy and photometry, studying in detail the characteristics of their stellar populations. All these objects are interacting galaxies: LMC and SMC are in an early stage of a minor merger event, and Sgr is being disrupted by the tidal field of the MW. There is a plenty of literature regarding the chemical composition of these systems, however, the extension of these galaxies prevents a complete and homogeneous analysis. Therefore, we homogeneously analysed stellar spectra belonging to MW and its satellites galaxies and we derived their chemical compositions. We highlighted the importance of a homogeneous analysis in the comparison among different galaxies or different samples, to avoid systematics due to different methods or physical assumptions.

Probing the innermost regions of AGN via time-resolved spectral analysis

The dynamics and geometry of the material inflowing and outflowing close to the supermassive black hole in active galactic nuclei are still uncertain. X-rays are the most suitable way to study the AGN innermost regions because of the Fe Kα emission line, a proxy of accretion, and Fe absorption lines produced by outflows. Winds are typically classified as Warm Absorbers (slow and mildly ionized) and Ultra Fast Outflows (fast and highly ionized). Transient Obscurers -optically thick winds that produce strong spectral hardening in X-rays, lasting from days to months- have been observed recently. Emission and absorption features vary on time-scales from hours to years, probing phenomena at different distances from the SMBH. In this work, we use time-resolved spectral analysis to investigate the accretion and ejection flows, to characterize them individually and search for correlations. We analyzed XMM-Newtomn data of a set of the brightest Seyfert 1 galaxies that went through an obscuration event: NGC 3783, NGC 3227, NGC 5548, and NGC 985. Our aim is to search for emission/absorption lines in short-duration spectra (∼ 10ks), to explore regions as close as the SMBH as the statistics allows for, and possibly catch transient phenomena. First we run a blind search to detect emission/absorption features, then we analyze their evolution with Residual Maps: we visualize simultaneously positive and negative residuals from the continuum in the time-energy plane, looking for patterns and relative time-scales. In NGC 3783 we were able to ascribe variations of the Fe Kα emission line to absorptions at the same energy due to clumps in the obscurer, whose presence is detected at >3σ, and to determine the size of the clumps. In NGC 3227 we detected a wind at ∼ 0.2c at ∼ 2σ, briefly appearing during an obscuration event.

Towards a full cosmological exploitation of cosmic voids

A stately fraction of the Universe volume is dominated by almost empty space. Alongside the luminous filamentary structures that make it up, there are vast and smooth regions that have remained outside the Cosmology spotlight during the past decades: cosmic voids. Although essentially devoid of matter, voids enclose fundamental information about the cosmological framework and have gradually become an effective and competitive cosmological probe. In this Thesis work we present fundamental results about the cosmological exploitation of voids. We focused on the number density of voids as a function of their radius, known as void size function, developing an effective pipeline for its cosmological usage. We proposed a new parametrisation of the most used theoretical void size function to model voids identified in the distribution of biased tracers (i.e. dark matter haloes, galaxies and galaxy clusters), a step of fundamental importance to extend the analysis to real data surveys. We then applied our built methodology to study voids in alternative cosmological scenarios. Firstly we exploited voids with the aim of breaking the degeneracies between cosmological scenarios characterised by modified gravity and the inclusion of massive neutrinos. Secondly we analysed voids in the perspective of the Euclid survey, focusing on the void abundance constraining power on dynamical dark energy models with massive neutrinos. Moreover we explored other void statistics like void profiles and clustering (i.e. the void-galaxy and the void-void correlation), providing cosmological forecasts for the Euclid mission. We finally focused on the probe combination, highlighting the incredible potential of the joint analysis of multiple void statistics and of the combination of the void size function with different cosmological probes. Our results show the fundamental role of the void analysis in constraining the fundamental parameters of the cosmological model and pave the way for future studies on this topic.

Unveiling the unknown of cool stars with high-resolution spectroscopy

Cool giant and supergiant stars are among the brightest populations in any stellar system and they are easily observable out to large distances, especially at infrared wavelengths. These stars also dominate the integrated light of star clusters in a wide range of ages, making them powerful tracers of stellar populations in more distant galaxies. High-resolution near-IR spectroscopy is a key tool for quantitatively investigating their kinematic, evolutionary and chemical properties. However, the systematic exploration and calibration of the NIR spectral diagnostics to study these cool stellar populations based on high-resolution spectroscopy is still in its pioneering stage. Any effort to make progress in the field is innovative and of impact on stellar archaeology and stellar evolution. This PhD project takes the challenge of exploring that new parameter space and characterizing the physical properties, the chemical content and the kinematics of cool giants and supergiants in selected disc fields and clusters of our Galaxy, with the ultimate goal of tracing their past and recent star formation and chemical enrichment history. By using optical HARPS-N and near-infrared GIANO-B high-resolution stellar spectra in the context of the large program SPA-Stellar Population Astrophysics: the detailed, age-resolved chemistry of the Milky Way disk” (PI L. Origlia), an extensive study of Arcturus, a standard calibrator for red giant stars, has been performed. New diagnostics of stellar parameters as well as optimal linelists for chemical analysis have been provided. Then, such diagnostics have been used to determine evolutionary properties, detailed chemical abundances of almost 30 different elements and mixing processes of a homogeneous sample of red supergiant stars in the Perseus complex.

A hydrodynamical view on the mass accretion onto massive halos

This Thesis presents the results of my work on how galaxy clusters form by the accretion of sub-clumps and diffuse materials, and how the accreted energy is distributed in the X-ray emitting plasma. Indeed, on scales larger than tens of millions of light years, the Universe is self-organised by gravity into a spiderweb, the Cosmic Web. Galaxy clusters are the knots of this Cosmic Web, but a strong definition of filaments (which link different knots) and their physical proprieties, is still uncertain. Even if this pattern was determined by studying the spatial distribution of galaxies in the optical band, recently, also in the X-rays probes of filamentary structures around galaxy clusters were obtained. Therefore, given these observational facilities, the galaxy clusters’ outskirts are the best candidate regions to detect filaments and study their physical characteristics. However, from X-rays observations, we have only a few detections of cosmic filaments to date. On the other hand, it is crucial to understand how the accreted energy is dissipated in the baryon content of galaxy clusters and groups. Indeed, it is well known that in the central region of galaxy clusters and groups, the baryon fraction increases with the halo mass. On the outer region, the lack of X-rays constraints influences our understanding of the evolution of baryons in the halos volume. The standard assumption of “closed-box” system, for which the baryon fraction should approach the cosmological ratio Omega_bar/Omega_m, for galaxy clusters and groups seems to be too strong, especially for less massive objects. Moreover, a complete redshift evolution of baryons in galaxy clusters and groups is still missing.