No Arabic abstract
We explore the possibility of inferring the properties of the Galactic neutron star population through machine learning. In particular, in this paper we focus on their dynamical characteristics and show that an artificial neural network is able to estimate with high accuracy the parameters which control the current positions of a mock population of pulsars. For this purpose, we implement a simplified population-synthesis framework (where selection biases are neglected at this stage) and concentrate on the natal kick-velocity distribution and the distribution of birth distances from the Galactic plane. By varying these and evolving the pulsar trajectories in time, we generate a series of simulations that are used to train and validate a suitably structured convolutional neural network. We demonstrate that our network is able to recover the parameters governing the kick-velocity and Galactic height distribution with a mean relative error of about $10^{-2}$. We discuss the limitations of our idealized approach and study a toy problem to introduce selection effects in a phenomenological way by incorporating the observed proper motions of 216 isolated pulsars. Our analysis highlights that increasing the sample of pulsars with accurate proper motion measurements by a factor of $sim$10, one of the future breakthroughs of the Square Kilometer Array, we might succeed in constraining the birth spatial and kick-velocity distribution of the neutron stars in the Milky Way with high precision through machine learning.
We study individual pulses of Vela (PSR B0833-45,/,J0835-4510) from daily observations of over three hours (around 120,000 pulses per observation), performed simultaneously with the two radio telescopes at the Argentine Institute of Radioastronomy. We select 4 days of observations in January-March 2021 and study their statistical properties with machine learning techniques. We first use density based DBSCAN clustering techniques, associating pulses mainly by amplitudes, and find a correlation between higher amplitudes and earlier arrival times. We also find a weaker (polarization dependent) correlation with the mean width of the pulses. We identify clusters of the so-called mini-giant pulses, with $sim10times$ the average pulse amplitude. We then perform an independent study, with Self-Organizing Maps (SOM) clustering techniques. We use Variational AutoEncoder (VAE) reconstruction of the pulses to separate them clearly from the noise and select one of the days of observation to train VAE and apply it to thre rest of the observations. We use SOM to determine 4 clusters of pulses per day per radio telescope and conclude that our main results are robust and self-consistent. These results support models for emitting regions at different heights (separated each by roughly a hundred km) in the pulsar magnetosphere. We also model the pulses amplitude distribution with interstellar scintillation patterns at the inter-pulses time-scale finding a characterizing exponent $n_{mathrm{ISS}}sim7-10$. In the appendices we discuss independent checks of hardware systematics with the simultaneous use of the two radio telescopes in different one-polarization / two-polarizations configurations. We also provide a detailed analysis of the processes of radio-interferences cleaning and individual pulse folding.
Through high-precision radio timing observations, we show that five recycled pulsars in the direction of the Galactic Centre (GC) have anomalous spin period time derivative ($dot P$) measurements -- PSRs J1748$-$3009, J1753$-$2819, J1757$-$2745, and J1804$-$2858 show negative values of $dot P$ and PSR J1801$-$3210 is found to have an exceptionally small value of $dot P$. We attribute these observed $dot P$ measurements to acceleration of these pulsars along their lines-of-sight (LOSs) due to the Galactic gravitational field. Using models of the Galactic mass distribution and pulsar velocities, we constrain the distances to these pulsars, placing them on the far-side of the Galaxy, providing the first accurate distance measurements to pulsars located in this region and allowing us to consider the electron density along these LOSs. We find the new electron density model YMW16 to be more consistent with these observations than the previous model NE2001. The LOS dynamics further constrain the model-dependent intrinsic $dot P$ values for these pulsars and they are consistent with measurements for other known pulsars. In the future, the independent distance measurements to these and other pulsars near the GC would allow us to constrain the Galactic gravitational potential more accurately.
We simulate the star cluster, made of stars in the main sequence and different black hole (BH) remnants, around SgrA* at the center of the Milky Way galaxy. Tracking stellar evolution, we find the BH remnant masses and construct the BH mass function. We sample 4 BH species and consider the impact of the mass-function in the dynamical evolution of system. Starting from an initial 6 dimensional family of parameters and using an MCMC approach, we find the best fits to various parameters of model by directly comparing the results of the simulations after $t = 10.5$ Gyrs with current observations of the stellar surface density, stellar mass profile and the mass of SgrA*. Using these parameters, we study the dynamical evolution of system in detail. We also explore the mass-growth of SgrA* due to tidally disrupted stars and swallowed BHs. We show that the consumed mass is dominated for the BH component with larger initial normalization as given by the BH mass-function. Assuming that about 10% of the tidally disrupted stars contribute in the growth of SgrA* mass, stars make up the second dominant effect in enhancing the mass of SgrA*. We consider the detectability of the GW signal from inspiralling stellar mass BHs around SgrA* with LISA. Computing the fraction of the lifetime of every BH species in the LISA band, with signal to noise ratio $gtrsim 8$, to their entire lifetime, and rescaling this number with the total number of BHs in the system, we find that the total expected rate of inspirals per Milky-Way sized galaxy per year is $10^{-5}$. Quite interestingly, the rate is dominated for the BH component with larger initial normalization as dictated by the BH mass-function. We interpret it as the second signature of the BH mass-function.
We use topological data analysis and machine learning to study a seminal model of collective motion in biology [DOrsogna et al., Phys. Rev. Lett. 96 (2006)]. This model describes agents interacting nonlinearly via attractive-repulsive social forces and gives rise to collective behaviors such as flocking and milling. To classify the emergent collective motion in a large library of numerical simulations and to recover model parameters from the simulation data, we apply machine learning techniques to two different types of input. First, we input time series of order parameters traditionally used in studies of collective motion. Second, we input measures based in topology that summarize the time-varying persistent homology of simulation data over multiple scales. This topological approach does not require prior knowledge of the expected patterns. For both unsupervised and supervised machine learning methods, the topological approach outperforms the one that is based on traditional order parameters.
Gamma-ray observations have shown pulsars to be efficient converters of rotational energy into GeV photons and it is of wide-ranging interest to determine their contribution to the gamma-ray background. We arrive at flux predictions from both the young (<~ Myr) and millisecond (~Gyr) Galactic pulsar populations. We find that unresolved pulsars can yield both a significant fraction of the excess GeV gamma rays near the Galactic Center and an inverse Compton flux in the inner kpc similar to that inferred by Fermi. We compare models of the young pulsar population and millisecond pulsar population to constraints from gamma-ray and radio observations. Overall, we find that the young pulsars should outnumber millisecond pulsars as unassociated gamma-ray point sources in this region. The number of young radio pulsars discovered near the Galactic Center is in agreement with our model of the young pulsar population. Deeper radio observations at higher latitudes can constrain the total gamma-ray emission from both young and millisecond pulsars from the inner galaxy. While this is a step towards better understanding of pulsars, cosmic rays in the Milky Way, and searches for dark matter, we also discuss a few interesting puzzles that arise from the underlying physics of pulsar emission and evolution.