Alma Mater Studiorum – Università di Bologna DOTTORATO DI RICERCA IN ASTROFISICA Ciclo XXX Settore Concorsuale: 02/C1 Settore Scientifico Disciplinare: FIS/05 “The Exotic Zoo of Millisecond Pulsars in Globular Clusters: a multi-wavelength study” Presentata da: Mario Cadelano Coordinatore Dottorato Supervisore Chiar.mo Prof. F. R. Ferraro Chiar.mo Prof. F. R. Ferraro Co-Relatrici arXiv:1901.04212v1 [astro-ph.HE] 14 Jan 2019 Chiar.ma Prof.ssa B. Lanzoni Dr.ssa C. Pallanca Esame finale anno 2018 Abstract This thesis work is focused on the study of millisecond pulsars in globular cluster by using multi-wavelength observations obtained with radio and optical telescopes. Radio observations have been used to search for and timing the pulsars. While classical search routines are based on the analysis of single and long time sequence of data, we present here an alternative method. This method exploits the large amount of available archival radio data to search for very faint pulsars by stacking all the daily power spectra. This method led to the discovery of three new pulsars in the stellar system Terzan 5. Optical observations have been exploited to search for millisecond pulsar optical counterparts, whose emission is totally dominated by the companion stars. Six new companion stars have been discovered, thus increasing by ~40% the total number of companions identified in globular clusters. In particular, four companions turned out to be He white dwarfs, as expected from the canonical formation scenario. One companion turned out to be a very faint and non-degenerate obJect, showing a strong variability which is likely the result of a strong heating of the stellar side exposed to the pulsar inJected flux. Its position in the color-magnitude diagram, not compatible with that of any evolutionary sequence of the cluster, also confirms that this system underwent a non-canonical evolution. Finally, one companion is a main-sequence star which shows a strong Hα emission likely due to a low-level accretion probably from a residual disk. This system can be located in an early evolutive stage, immediately preceding the reactivation of an already re-accelerated pulsar. Furthermore, we identified the companion star to a transient low- mass X-ray binary. This companion turned out to be a sub-giant branch star and therefore this system is likely in a very early phase of the mass accretion stage, when an expanding secondary star is filling its Roche-Lobe and transferring material that eventually spins-up the NS. Conclusion are drawn in the final chapter, where the evolution of millisecond pulsars is discussed and possible future developments are suggested. To my Walden... When we become the enfolders of those orbs, and the pleasure and knowledge of every thing in them, shall we be fill’d and satisfied then? No, we but level that lift to pass and continue beyond. - Walt Whitman, Leaves of Grass - TABLE OF CONTENTS Page List of Tables xiii List of Figures xv Introduction 1 1 From Neutron Stars to Millisecond Pulsars5 1.1 Neutron Stars . .6 1.1.1 The Internal Structure of a Neutron Star . .8 1.2 Phenomenology of Pulsars . .9 1.2.1 Pulse Properties . .9 1.2.2 Derivation of the Pulsar Physical Parameters . 10 1.3 Pulsar Magnetosphere and Emission Mechanism . 12 1.4 Pulsar Evolution . 15 1.5 Millisecond Pulsars . 17 1.5.1 Canonical Formation Scenarios . 18 1.5.2 The Spiders: Black-Widows and Redbacks . 21 1.5.3 Accreting and Transitional Millisecond Pulsars . 22 1.6 Globular Clusters: the Ideal Millisecond Pulsar Factories . 24 1.6.1 Dynamical Interactions . 26 1.7 Millisecond Pulsars in Globular Clusters: State of the Art . 28 2 Searching, Timing and Optical Identification of Millisecond Pulsars in Globu- lar Clusters 33 2.1 Searching for Pulsars . 34 2.1.1 Dispersion and Dedispersion of the Signal . 34 2.1.2 Standard Search Procedures . 36 2.1.3 Flux Density Estimates . 39 2.2 Pulsar Timing . 40 2.2.1 Isolated Pulsars . 41 ix TABLE OF CONTENTS 2.2.2 Binary Pulsars . 41 2.3 Identification of Optical Counterparts in Globular Clusters . 43 2.3.1 Photometry of the Companion Stars . 43 2.3.2 Spectroscopy of the Companion Stars . 47 3 Discovery of Three New Millisecond Pulsars in Terzan 5 49 3.1 Introduction . 50 3.2 Observations and Data Analysis . 51 3.2.1 Dataset and initial data reduction . 51 3.2.2 Stacking search procedures . 51 3.3 Results . 53 3.3.1 Ter5aj . 54 3.3.2 Ter5ak . 56 3.3.3 Ter5al . 57 3.4 Accelerations and physical parameters . 58 4 The Black-Widow PSR J1953+1846A in the Globular Cluster M71 61 4.1 Introduction . 62 4.2 Radio Timing . 62 4.3 Optical Photometry of the Companion Star . 65 4.3.1 Observations and data analysis . 65 4.3.2 The companion to M71A . 66 4.4 Discussion . 68 4.4.1 Reprocessing efficiency and Roche Lobe filling factor . 72 4.4.2 A comparison between M71A and M5C . 74 4.4.3 Comparing X-ray and optical light curves . 75 5 The He White Dwarfs Orbiting the Millisecond Pulsars in the Globular Clus- ter 47 Tucanae 77 5.1 Introduction . 78 5.2 Optical Photometry of the Star Cluster . 78 5.2.1 Observations and data analysis . 78 5.2.2 Astrometry . 79 5.2.3 Identification of the MSP companions . 80 5.3 Discussion . 83 5.3.1 The physical properties of the He WD companions . 83 5.3.2 Possible candidate companion stars . 87 5.3.3 Non detections . 88 5.3.4 The companion to the RB 47TucW . 89 x TABLE OF CONTENTS 6 The Optical Counterpart to the Accreting Millisecond X-ray Pulsar SAX J1748.9- 2021 in the Globular Cluster NGC 6440 91 6.1 Introduction . 92 6.2 Observations and Data Reduction . 93 6.3 Results . 94 6.4 Discussion . 97 7 The Optical Identification of the X-ray Burster EXO 1745-248 in Terzan 5 103 7.1 Introduction . 104 7.2 Observations and Data Reduction . 105 7.3 Results . 106 Summary and Conclusions 111 A Proper Motions and Structural Parameters of the Galactic Globular Cluster M71 119 A.1 Introduction . 120 A.2 Observations and Data Reduction . 121 A.3 Relative Proper Motions . 124 A.4 Absolute Proper Motions . 126 A.4.1 The cluster orbit . 130 A.5 Gravitational Center, Structural Parameters and Initial Mass . 131 A.5.1 Gravitational center . 131 A.5.2 Stellar density profile . 132 A.5.3 Cluster initial mass . 135 A.6 Summary . 136 Bibliography 139 Acknowledgements 183 xi LIST OF TABLES TABLE Page 3.1 Timing parameters for the new Terzan 5 MSPs. 56 4.1 Timing parameters for PSR J1953+1847 (M71A). Numbers in parentheses are uncer- tainties in the last digits quoted. 63 4.2 Optical observations of COM-M71A . 69 4.3 Optical properties of COM-M71A . 70 5.1 Radio timing ephemeris of the analyzed MSPs. From left to right: MSP name, posi- tion, offset from the GC center, orbital period, mass function and characteristic age. Numbers in parentheses are uncertainties in the last digits quoted. Reference: Ridolfi et al.(2016); Freire et al.(2017). ................................ 79 5.2 Optical properties of the companion stars. From left to right: MSP companion name, position, distance from the radio MSP nominal position, F300X and F390W magnitudes and the relatives uncertainties. 83 5.3 Derived properties of the five MSPs with He WD companions. From top to bottom: companion mass, age, temperature, luminosity, PSR mass and inclination angle. 85 xiii LIST OF FIGURES FIGURE Page 1.1 Schematic representation of the NS internal structure. .9 1.2 Integrated pulse profiles for a sample of PSRs. 10 1.3 Schematic representation of a NS and its magnetosphere . 13 1.4 PSR magnetosphere in the Goldreich-Julian model. 13 1.5 The P P˙ for a large sample of radio PSRs. 16 ¡ 1.6 Surfaces of equal gravitational potential in a binary system. 18 1.7 MSPs formation scenarios. ..