A spectroscopic and photometric study of MSP companions in Galactic Globular Clusters

This Thesis is devoted to the study of the optical companions of Millisecond Pulsars in Galactic Globular Clusters (GCs) as a part of a large project started at the Department of Astronomy of the Bologna University, in collaboration with other institutions (Astronomical Observatory of Cagliari and Bologna, University of Virginia), specifically dedicated to the study of the environmental effects on passive stellar evolution in galactic GCs. Globular Clusters are very efficient “Kilns” for generating exotic object, such as Millisecond Pulsars (MSP), low mass X-ray binaries(LMXB) or Blue Straggler Stars (BSS). In particular MSPs are formed in binary systems containing a Neutron Star which is spun up through mass accretion from the evolving companion (e.g. Bhattacharia & van den Heuvel 1991). The final stage of this recycling process is either the core of a peeled star (generally an Helium white dwarf) or a very light almos exhausted star, orbiting a very fast rotating Neutron Star (a MSP). Despite the large difference in total mass between the disk of the Galaxy and the Galactic GC system (up a factor 103), the percentage of fast rotating pulsar in binary systems found in the latter is very higher. MSPs in GCs show spin periods in the range 1.3 ÷ 30ms, slowdown rates ˙P 1019 s/s and a lower magnetic field, respect to ”normal” radio pulsars, B 108 gauss . The high probability of disruption of a binary systems after a supernova explosion, explain why we expect only a low percentage of recycled millisecond pulsars respect to the whole pulsar population. In fact only the 10% of the known 1800 radio pulsars are radio MSPs. Is not surprising, that MSP are overabundant in GCs respect to Galactic field, since in the Galactic Disk, MSPs can only form through the evolution of primordial binaries, and only if the binary survives to the supernova explosion which lead to the neutron star formation. On the other hand, the extremely high stellar density in the core of GCs, relative to most of the rest of the Galaxy, favors the formation of several different binary systems, suitable for the recycling of NSs (Davies at al. 1998). In this thesis we will present the properties two millisecond pulsars companions discovered in two globular clusters, the Helium white dwarf orbiting the MSP PSR 1911-5958A in NGC 6752 and the second case of a tidally deformed star orbiting an eclipsing millisecond pulsar, PSR J1701-3006B in NGC6266

A panchromatic study of stellar populations in Galactic Globular Clusters

The aim of this Thesis is to investigate (i) how common the bimodal Blue Straggler Stars (BSS) radial distribution is in stellar clusters and (ii) which are the physical processes that can produce this bimodality. We discuss possible relations between the properties of the BSS radial distribution and the dynamical state of the hosting clusters by making use of dynamical models and simulations. When relevant, we also discuss the possible links with some cluster "anomalies" and the effects of a massive object (like Imtermediate Mass Black Hole) in the cluster center. To this purpose we present the observational multiwavelength studies of the BSS populations and their radial distributions in 5 GGCs: M5, M55, M2, NGC 2419 and NGC 6388.

Mass Loss in Population II Red Giants: an IRAC@Spitzer Survey of Galactic Globular Clusters

With the goal of studying ML along the RGB, mid-IR observations of a carefully selected sample of 17 Galactic globular clusters (GGCs) with different metallicity and horizontal branch (HB) morphology have been secured with IRAC on board Spitzer: a global sample counting about 8000 giant has been obtained. Suitable complementary photometry in the optical and near-IR has been also secured in order to properly characterize the stellar counterparts to the Spitzer sources and their photospheric parameters. Stars with color (i.e. dust) excess have been identified, their likely circumstellar emission quantified and modelled, and empirical estimates of mass loss rates and timescales obtained. We find that mass loss rates increases with increasing stellar luminosity and decreasing metallicity. For a given luminosity, we find that ML rates are systematically higher than the prediction by extrapolating the Reimers law. CMDs constructed from ground based near-IR and IRAC bands show that at a given luminosity some stars have dusty envelopes and others do not. From this, we deduce that the mass loss is episodic and is ``on'' for some fraction of the time. The total mass lost on the RGB can be easily computed by multiplying ML rates by the ML timescales and integrating over the evolutionary timescale. The average total mass lost moderately increases with increasing metallicity, and for a given metallicity is systematically higher in clusters with extended blue HB.

Chemical and kinematical properties of Blue Straggler stars in Galactic globular clusters

Blue straggler stars (BSSs) are brighter and bluer (hotter) than the main-sequence (MS) turnoff and they are known to be more massive than MS stars.Two main scenarios for their formation have been proposed:collision-induced stellar mergers (COL-BSSs),or mass-transfer in binary systems (MT-BSSs).Depleted surface abundances of C and O are expected for MT-BSSs,whereas no chemical anomalies are predicted for COL-BSSs.Both MT- and COL-BSSs should rotate fast, but braking mechanisms may intervene with efficiencies and time-scales not well known yet,thus preventing a clear prediction of the expected rotational velocities.Within this context,an extensive survey is ongoing by using the multi-object spectrograph FLAMES@VLT,with the aim to obtain abundance patterns and rotational velocities for representative samples of BSSs in several Galactic GCs.A sub-population of CO-depleted BSSs has been identified in 47 Tuc,with only one fast rotating star detected.For this PhD Thesis work I analyzed FLAMES spectra of more than 130 BSSs in four GCs:M4,NGC 6397,M30 and ω Centauri.This is the largest sample of BSSs spectroscopically investigated so far.Hints of CO depletion have been observed in only 4-5 cases (in M30 and ω Centauri),suggesting either that the majority of BSSs have a collisional origin,or that the CO-depletion is a transient phenomenon.Unfortunately,no conclusions in terms of formation mechanism could be drawn in a large number of cases,because of the effects of radiative levitation. Remarkably,however,this is the first time that evidence of radiative levitation is found in BSSs hotter than 8200 K.Finally, we also discovered the largest fractions of fast rotating BSSs ever observed in any GCs:40% in M4 and 30% in ω Centauri.While not solving the problem of BSS formation,these results provide invaluable information about the BSS physical properties,which is crucial to build realistic models of their evolution.

Multiple stellar populations in globular clusters with photometry and low resolution spectroscopy

Our view of Globular Clusters has deeply changed in the last decade. Modern spectroscopic and photometric data have conclusively established that globulars are neither coeval nor monometallic, reopening the issue of the formation of such systems. Their formation is now schematized as a two-step process, during which the polluted matter from the more massive stars of a first generation gives birth, in the cluster innermost regions, to a second generation of stars with the characteristic signature of fully CNO-processed matter. To date, star-to-star variations in abundances of the light elements (C, N, O, Na) have been observed in stars of all evolutionary phases in all properly studied Galactic globular clusters. Multiple or broad evolutionary sequences have also been observed in nearly all the clusters that have been observed with good signal-to-noise in the appropriate photometric bands. The body of evidence suggests that spreads in light-element abundances can be fairly well traced by photometric indices including near ultraviolet passbands, as CNO abundance variations affect mainly wavelengths shorter than ~400 nm owing to the rise of some NH and CN molecular absorption bands. Here, we exploit this property of near ultraviolet photometry to trace internal chemical variations and combined it with low resolution spectroscopy aimed to derive carbon and nitrogen abundances in order to maximize the information on the multiple populations. This approach has been proven to be very effective in (i) detecting multiple population, (ii) characterizing their global properties (i.e., relative fraction of stars, location in the color-magnitude diagram, spatial distribution, and trends with cluster parameters) and (iii) precisely tagging their chemical properties (i.e., extension of the C-N anticorrelation, bimodalities in the N content).

COSMIC-LAB: Unexpected Results from High-resolution Spectra of AGB Stars in Globular Clusters

We have used high-resolution spectra, acquired with UVES@ESO-VLT, to determine the chemical abundances of different samples of AGB and RGB stars in 4 Galactic globular clusters, namely 47Tuc, NGC3201, M22 and M62. For almost all the analyzed AGB stars we found a clear discrepancy between the iron abundance measured from neutral lines and that obtained from single ionized lines, while this discrepancy is not obtained for the RGB samples observed in the same clusters and analyzed with the same procedure. Such a behavior exactly corresponds to what expected in the case of Non-Local Thermodynamical Equilibrium (NLTE) in the star atmosphere. These results have a huge impact on the proper determination of GC chemistry. In fact, one of the most intriguing consequences is that, at odds with previous claims, no iron spread is found in NGC3201 and M22 if the iron abundance is obtained from ionized lines only.

Exploring the internal kinematics of Galactic Globular Cluster cores

Globular clusters (GCs) are traditionally described as simple quasi-relaxed non-rotating stellar systems, characterized by spherical symmetry and isotropy in velocity space. However, recent studies have shown deviations from isotropic velocity distributions and significant internal rotation in many GCs, suggesting that their internal structure and kinematics are more complex than previously thought. The aim of this thesis is to investigate the internal kinematics of Galactic Globular Clusters (GGCs) as part of the Multi-Instrument Kinematic Survey (MIKiS), which exploits the capabilities of different ESO-VLT spectrographs to obtain comprehensive velocity dispersion (VD) and rotation profiles of GGCs. Moreover, this thesis has the particular goal of unraveling the kinematics of GC cores, which are still largely unexplored, by taking advantage of the exceptional spatial resolution of the adaptive-optics assisted integral-field spectrograph MUSE/NFM. The thesis presents a thorough kinematic study of three GGCs NGC 1904, NGC 6440, and NGC 6569. By combining the data sets acquired with four different spectrographs, we obtained the radial velocity (RV) of more than 1000 individual stars in each cluster, sampling from the innermost to the outermost regions. This allowed us to obtain the entire VD profile of each cluster and exclude the presence of an intermediate-mass black hole in the core of NGC 1904, at odds with previous findings obtained from integrated-light spectra. The studies also revealed signatures of internal rotation in each of the GCs studied. These results, supported by those of N-body simulations, prove that GCs were born with a significant initial rotation that they gradually lost through internal two-body relaxation and angular momentum loss carried away by escaping stars. Furthermore, we derived the structural parameters of NGC 6440 and NGC 6569, obtaining a comprehensive overview of the internal kinematics and structure of these GCs, which is necessary to properly reconstruct the evolutionary history of these systems.

Exploring the white dwarf cooling sequence in Globular Clusters

White dwarfs (WDs) are electron-degenerate structures that are commonly assumed to evolve via a pure cooling process, with no stable thermonuclear activity at work. Their cooling rate is adopted as a cosmic chronometer to constrain the age of several Galactic populations, including the disk, Globular Clusters (GCs) and open clusters. This thesis work is aimed at the study of the WD populations in globular clusters and is articulated in two branches. The first was focused on the study of the bright portion of the WD cooling sequence. By analyzing high resolution UV data acquired with the Hubble Space Telescope (HST), we compared the WD luminosity functions (LFs) in four Galactic GCs (namely M13, M3, NGC6752, and M5) finding an unexpected over-abundance of WDs in M13 and NGC6752 with respect to M3 and M5. Theoretical models suggest that, consistently with the blue-tail horizontal branch (HB) morphology of M13 and NGC6752, this overabundance is due to a population of slowly cooling WDs, i.e., WDs fading more slowly than in a pure cooling process thanks to an extra-energy source provided by stable thermonuclear burning in their residual hydrogen-rich envelope. This is the first empirical evidence of WDs fading at a slower rate than usually assumed, and has a crucial impact on the use of the cooling sequence as a cosmic chronometer. The second branch was focused on the search for the companion star to binary millisecond Pulsars (MSP) in the globular clusters M13 and NGC 6652: the identified companions turned out to be helium-core WDs, and provided a invaluable constraints on the mass of the neutron star and the epoch of the MSP formation.

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.

A spectroscopic and photometric study of stellar populations in a sample of clusters in the Magellanic Clouds

This Ph.D. Thesis has been carried out in the framework of a long-term and large project devoted to describe the main photometric, chemical, evolutionary and integrated properties of a representative sample of Large and Small Magellanic Cloud (LMC and SMC respectively) clusters. The globular clusters system of these two Irregular galaxies provides a rich resource for investigating stellar and chemical evolution and to obtain a detailed view of the star formation history and chemical enrichment of the Clouds. The results discussed here are based on the analysis of high-resolution photometric and spectroscopic datasets obtained by using the last generation of imagers and spectrographs. The principal aims of this project are summarized as follows: • The study of the AGB and RGB sequences in a sample of MC clusters, through the analysis of a wide near-infrared photometric database, including 33 Magellanic globulars obtained in three observing runs with the near-infrared camera SOFI@NTT (ESO, La Silla). • The study of the chemical properties of a sample of MCs clusters, by using optical and near-infrared high-resolution spectra. 3 observing runs have been secured to our group to observe 9 LMC clusters (with ages between 100 Myr and 13 Gyr) with the optical high-resolution spectrograph FLAMES@VLT (ESO, Paranal) and 4 very young (<30 Myr) clusters (3 in the LMC and 1 in the SMC) with the near-infrared high-resolution spectrograph CRIRES@VLT. • The study of the photometric properties of the main evolutive sequences in optical Color- Magnitude Diagrams (CMD) obtained by using HST archive data, with the final aim of dating several clusters via the comparison between the observed CMDs and theoretical isochrones. The determination of the age of a stellar population requires an accurate measure of the Main Sequence (MS) Turn-Off (TO) luminosity and the knowledge of the distance modulus, reddening and overall metallicity. For this purpose, we limited the study of the age just to the clusters already observed with high-resolution spectroscopy, in order to date only clusters with accurate estimates of the overall metallicity.

Cosmic-Lab: Optical companions to binary Millisecond Pulsars

Millisecond Pulsars (MSPs) are fast rotating, highly magnetized neutron stars. According to the "canonical recycling scenario", MSPs form in binary systems containing a neutron star which is spun up through mass accretion from the evolving companion. Therefore, the final stage consists of a binary made of a MSP and the core of the deeply peeled companion. In the last years, however an increasing number of systems deviating from these expectations has been discovered, thus strongly indicating that our understanding of MSPs is far to be complete. The identification of the optical companions to binary MSPs is crucial to constrain the formation and evolution of these objects. In dense environments such as Globular Clusters (GCs), it also allows us to get insights on the cluster internal dynamics. By using deep photometric data, acquired both from space and ground-based telescopes, we identified 5 new companions to MSPs. Three of them being located in GCs and two in the Galactic Field. The three new identifications in GCs increased by 50% the number of such objects known before this Thesis. They all are non-degenerate stars, at odds with the expectations of the "canonical recycling scenario". These results therefore suggest either that transitory phases should also be taken into account, or that dynamical processes, as exchange interactions, play a crucial role in the evolution of MSPs. We also performed a spectroscopic follow-up of the companion to PSRJ1740-5340A in the GC NGC 6397, confirming that it is a deeply peeled star descending from a ~0.8Msun progenitor. This nicely confirms the theoretical expectations about the formation and evolution of MSPs.

New indicators of star cluster dynamical evolution

This thesis presents a study of globular clusters (GCs), based on analysis of Monte Carlo simulations of globular clusters (GCs) with the aim to define new empirical parameters measurable from observations and able to trace the different phases of their dynamical evolution history. During their long term dynamical evolution, due to mass segregation and and dynamical friction, massive stars transfer kinetic energy to lower-mass objects, causing them to sink toward the cluster center. This continuous transfer of kinetic energy from the core to the outskirts triggers the runaway contraction of the core, known as "core collapse" (CC), followed by episodes of expansion and contraction called gravothermal oscillations. Clearly, such an internal dynamical evolution corresponds to significant variations also of the structure of the system. Determining the dynamical age of a cluster can be challenging as it depends on various internal and external properties. The traditional classification of GCs as CC or post-CC systems relies on detecting a steep power-law cusp in the central density profile, which may not always be reliable due to post-CC oscillations or other processes. In this thesis, based on the normalized cumulative radial distribution (nCRD) within a fraction of the half-mass radius is analyzed, and three diagnostics (A5, P5, and S2.5) are defined. These diagnostics show sensitivity to dynamical evolution and can distinguish pre-CC clusters from post-CC clusters.The analysis performed using multiple simulations with different initial conditions, including varying binary fractions and the presence of dark remnants showed the time variations of the diagnostics follow distinct patterns depending on the binary fraction and the retention or ejection of black holes. This analysis is extended to a larger set of simulations matching the observed properties of Galactic GCs, and the parameters show a potential to distinguish the dynamical stages of the observed clusters as well.

Homogeneous analysis of a sample of Open Clusters in the context of the BOCCE project and the Gaia-ESO Survey

The open clusters (OC) are gravitationally bound systems of a few tens or hundreds of stars. In our Galaxy, the Milky Way, we know about 3000 open clusters, of very different ages in the range of a few millions years to about 9 Gyr. OCs are mainly located in the Galactic thin disc, with distances from the Galactic centre in the range 4-22 kpc and a height scale on the disc of about 200 pc. Their chemical properties trace those of the environment in which they formed and the metallicity is in the range -0.5<[Fe/H]<+0.5 dex. Through photometry and spectroscopy it is possible to study relatively easily the properties of the OCs and estimate their age, distance, and chemistry. For these reasons they are considered primary tracers of the chemical properties and chemical evolution of the Galactic disc. The main subject of this thesis is the comprehensive study of several OCs. The research embraces two different projects: the Bologna Open Cluster Chemical Evolution project (BOCCE) and the Gaia-ESO Survey. The first is a long-term programme, aiming at studying the chemical evolution of the Milky Way disc by means of a homogeneous sample of OCs. The latter is a large public spectroscopy survey, conducted with the high-resolution spectrograph FLAMES@VLT and targeting about 10^5 stars in different part of the Galaxy and 10^4 stars in about 100 OCs. The common ground between the two projects is the study of the properties of the OCs as tracers of the disc's characteristics. The impressive scientific outcome of the Gaia-ESO Survey and the unique framework of homogeneity of the BOCCE project can propose, especially once combined together, a much more accurate description of the properties of the OCs. In turn, this will give fundamental constraints for the interpretation of the properties of the Galactic disc.

A comprehensive study of the AGN feedback cycle in galaxy clusters from high resolution X-ray and radio observations

At the center of galaxy clusters, a dramatic interplay known as feedback cycle occurs between the hot intracluster medium (ICM) and the active galactic nucleus (AGN) of the central galaxy. The footprints of this interplay are evident from X-ray observations of the ICM, where X-ray cavities and shock fronts are associated with radio lobe emission tracing energetic AGN outbursts. While such jet activity reduces the efficiency of the hot gas to cool to lower temperatures, residual cooling can generate warm and cold gas clouds around the central galaxy. The condensed gas parcels can ultimately reach the core of the galaxy and be accreted by the AGN. This picture is the result of tremendous advances over the last three decades. Yet, a deeper understanding of the details of how the heating–cooling regulation is achieved and maintained is still missing. In this Thesis, we delve into key aspects of the feedback cycle. To this end, we leverage high-resolution (sub-arcsecond), multifrequency observations (mainly X-ray and radio) of several top-level facilities (e.g., Chandra, JVLA, VLBA, LOFAR). First, we investigate which conditions trigger a feedback response to gas cooling, by studying the properties of clusters where feedback is just about to start. Then, we focus on the details of how the AGN–ICM interaction progresses by examining cavity and shock heating in the cluster RBS797, an exemplary case of the jet feedback paradigm. Furthermore, we explore the importance of shock heating and the coupling of distinct jet power regimes (i.e., FRII, FRI and FR0 radio galaxies) to the environment. Ultimately, as heating models rely on the connection between the direct evidence (the jets) and the smoking gun (the X-ray cavities) of feedback, we examine the cases in which these two are dramatically misaligned.