We constructed a Hubble Space Telescope (HST) astro-photometric catalog of the central region of the Galactic globular cluster NGC 1261. This catalog, complemented with Gaia DR2 data sampling the external regions, has been used to estimate the structural parameters of the system (i.e., core, half-mass, tidal radii and concentration) from its resolved star density profile. We computed high-precision proper motions thanks to multi-epoch HST data and derived the cluster velocity dispersion profile in the plane of the sky for the innermost region, finding that the system is isotropic. The combination with line-of-sight information collected from spectroscopy in the external regions provided us with the cluster velocity dispersion profile along the entire radial extension. We also measured the absolute proper motion of NGC 1261 using a few background galaxies as a reference. The radial distribution of the Blue Straggler Star population shows that the cluster is in a low/intermediate phase of dynamical evolution.

We used a set of moderately-deep and high-resolution optical observations obtained with the Hubble Space Telescope to investigate the properties of the stellar population in the heavily obscured bulge globular cluster NGC 6256. The analysis of the color-magnitude diagram revealed a stellar population with an extended blue horizontal branch and severely affected by differential reddening, which was corrected taking into account color excess variations up to \delta E(B-V) ~ 0.51. We implemented a Monte Carlo Markov Chain technique to perform the isochrone fitting of the observed color-magnitude diagram in order to derive the stellar age, the cluster distance and the average color excess in the cluster direction. Using different set of isochrones we found that NGC 6256 is characterized by a very old stellar age around 13.0 Gyr, with a typical uncertainty of ~ 0.5 Gyr. We also found an average color excess E(B-V) = 1.19 and a distance from the Sun of 6.8 kpc. We then derived the cluster gravitational center and measured its absolute proper motion using the Gaia-DR2 catalog. All this was used to back-integrate the cluster orbit in a Galaxy-like potential and measure its integrals of motion. It turned out that NGC 6256 is currently in a low-eccentricity orbit entirely confined within the bulge and its integrals of motion are fully compatible with a cluster purely belonging to the Galaxy native globular cluster population. All these pieces of evidence suggest that NGC 6256 is an extremely old relic of the past history of the Galaxy, formed during the very first stages of its assembly.

We used resolved star counts from Hubble Space Telescope images to determine the center of gravity and the projected density profiles of 6 old globular clusters in the Large Magellanic Cloud (LMC), namely NGC 1466, NGC 1841, NGC 1898, NGC 2210, NGC 2257 and Hodge 11. For each system, the LMC field contribution was properly taken into account by making use, when needed, of parallel HST observations. The derived values of the center of gravity may differ by several arcseconds (corresponding to more than 1 pc at the distance of the LMC) from previous determinations. The cluster density profiles are all well fit by King models, with structural parameters that may differ from the literature ones by even factors of two. Similarly to what observed for Galactic globular clusters, the ratio between the effective and the core radii has been found to anti-correlate with the cluster dynamical age.

The distribution of size as a function of age observed for star clusters in the Large Magellanic Cloud (LMC) is very puzzling: young clusters are all compact, while the oldest systems show both small and large sizes. It is commonly interpreted as due to a population of binary black holes driving a progressive expansion of cluster cores. Here we propose, instead, that it is the natural consequence of the fact that only relatively low-mass clusters have formed in the last ~3 Gyr in the LMC and only the most compact systems survived and are observable. The spread in size displayed by the oldest (and most massive) clusters, instead, can be explained in terms of initial conditions and internal dynamical evolution. To quantitatively explore the role of the latter, we selected a sample of five coeval and old LMC clusters with different sizes, and we estimated their dynamical age from the level of central segregation of blue straggler stars (the so-called dynamical clock). Similarly to what found in the Milky Way, we indeed measure different levels of dynamical evolution among the selected coeval clusters, with large-core systems being dynamically younger than those with small size. This behaviour is fully consistent with what expected from internal dynamical evolution processes over timescales mainly set by the structure of each system at formation.

We used optical images acquired with the UVIS channel of the Wide Field Camera 3 on board of the Hubble Space Telescope to construct the first high-resolution extinction map in the direction of NGC 6440, a globular cluster located in the bulge of our Galaxy. The map has a spatial resolution of 0.5" over a rectangular region of about 160" X 240" around the cluster center, with the long side in the North-West/South-East direction. We found that the absorption clouds show patchy and filamentary sub-structures with extinction variations as large as $\delta {\rm E}(B-V)\sim0.5$ mag. We also performed a first-order proper motion analysis to distinguish cluster members from field interlopers. After the field decontamination and the differential reddening correction, the cluster sequences in the color-magnitude diagram appear much better defined, providing the best optical color-magnitude diagram so far available for this cluster.

We have constructed the Spectral Energy Distributions (SEDs) of a sample of Blue Straggler Stars (BSSs) in the core of the globular cluster 47 Tucanae, taking advantage of the large set of high resolution images, ranging from the ultraviolet to the near infrared, obtained with the ACS/HRC camera of the Hubble Space Telescope. Our final BSS sample consists of 22 objects, spanning the whole color and magnitude extension of the BSS sequence in 47 Tucanae. We fitted the BSS broadband SEDs with models to derive temperature, surface gravity, radius, luminosity and mass. We show that BSSs indeed define a mass sequence, where the mass increases for increasing luminosity. Interestingly, the BSS masses estimates from the SED fitting turn out to be comparable to those derived from the projection of the stellar position in the color-magnitude diagram onto standard star evolutionary tracks. We compare our results with previous, direct mass estimates of a few BSSs in 47 Tucanae. We also find a couple of supermassive BSS candidates, i.e., BSSs with a mass larger than twice the turn-off mass, the formation of which must have involved more than two progenitors.

We report on the discovery of the companion star to the millisecond pulsar J1342+2822B in the globular cluster M3. We exploited a combination of near-ultraviolet and optical observations acquired with the Hubble Space Telescope in order to search for the optical counterparts to the known millisecond pulsars in this cluster. At a position in excellent agreement with that of the radio pulsar J1342+2822B (M3B), we have identified a blue and faint object (mF275W \approx 22.45) that, in the color-magnitude diagram of the cluster, is located in the region of He core white dwarfs. From the comparison of the observed magnitudes with theoretical cooling tracks we have estimated the physical properties of the companion star: it has a mass of only 0.19 \pm 0.02 Msun, a surface temperature of 12 \pm 1 x 10^3 K and a cooling age of 1.0\pm0.2 Gyr. Its progenitor was likely a ~ 0.84 M star and the bulk of the mass-transfer activity occurred during the sub-giant branch phase. The companion mass, combined with the pulsar mass function, implies that this system is observed almost edge-on and that the neutron star has a mass of 1.1 \pm 0.3 Msun, in agreement with the typical values measured for recycled neutron stars in these compact binary systems. We have also identified a candidate counterpart to the wide and eccentric binary millisecond pulsar J1342+2822D. It is another white dwarf with a He core and a mass of 0.22 \pm 0.2 Msun, implying that the system is observed at a high inclination angle and hosts a typical NS with a mass of 1.3 \pm 0.3 Msun. At the moment, the large uncertainty on the radio position of this millisecond pulsar prevents us from robustly concluding that the detected star is its optical counterpart.

As part of the ESO-VLT Multi-Instrument Kinematic Survey (MIKiS) of Galactic globular clusters, we present a detailed investigation of the internal kinematics of NGC 5986. The analysis is based on about 300 individual radial velocities of stars located at various distances from the cluster center, up to 300 arcseconds (about 4 half-mass radii). Our analysis reveals the presence of a solid-body rotation extending from the cluster center to the outermost regions probed by the data, and a velocity dispersion profile initially declining with the distance from the cluster's center, but flattening and staying constant at ~5 km/s for distances larger than about one half-mass radius. This is the first globular cluster for which evidence of the joint presence of solid-body rotation and flattening in the outer velocity dispersion profile is found. The combination of these two kinematical features provides a unique opportunity to shed light on fundamental aspects of globular cluster dynamics and probe the extent to which internal relaxation, star escape, angular momentum transport and loss, and the interaction with the Galaxy tidal field can affect a cluster's dynamical evolution and determine its current kinematical properties. We present the results of a series of N-body simulations illustrating the possible dynamical paths leading to kinematic features like those observed in this cluster and the fundamental dynamical processes that underpin them.

The parameter A+, defined as the area enclosed between the cumulative radial distribution of blue straggler stars (BSSs) and that of a reference population, is a powerful indicator of the level of BSS central segregation. As part of the Hubble Space Telescope UV Legacy Survey of Galactic globular clusters (GCs), here we present the BSS population and the determination of A+ in 27 GCs observed out to about one half-mass radius. In combination with 21 additional clusters discussed in a previous paper this provides us with a global sample of 48 systems (corresponding to \sim 32\% of the Milky Way GC population), for which we find a strong correlation between A+ and the ratio of cluster age to the current central relaxation time. Tight relations have been found also with the core radius and the central luminosity density, which are expected to change with the long-term cluster dynamical evolution. An interesting relation is emerging between A+ and the ratio of the BSS velocity dispersion relative to that of main sequence turn-off stars, which measures the degree of energy equipartition experienced by BSSs in the cluster. These results provide further confirmation that BSSs are invaluable probes of GC internal dynamics and A+ is a powerful dynamical clock.

In the context of the ESO-VLT Multi-Instrument Kinematic Survey (MIKiS) of Galactic globular clusters, we present the line-of-sight rotation curve and velocity dispersion profile of M5 (NGC 5904), as determined from the radial velocity of more than 800 individual stars observed out to 700" (~ 5 half-mass radii) from the center. We find one of the cleanest and most coherent rotation patterns ever observed for globular clusters, with a very stable rotation axis (having constant position angle of 145^o at all surveyed radii) and a well-defined rotation curve. The density distribution turns out to be flattened in the direction perpendicular to the rotation axis, with a maximum ellipticity of 0.15. The rotation velocity peak (~3 km/s in projection) is observed at ~0.6 half-mass radii, and its ratio with respect to the central velocity dispersion (~0.3-0.4 at 4 projected half-mass radii) indicates that ordered motions play a significant dynamical role. This result strengthens the growing empirical evidence of the kinematic complexity of Galactic globular clusters and motivates the need of fundamental investigations of the role of angular momentum in collisional stellar dynamics.

We present the first results of the Multi-Instrument Kinematic Survey of Galactic Globular Clusters, a project aimed at exploring the internal kinematics of a representative sample of Galactic globular clusters from the radial velocity of individual stars, covering the entire radial extension of each system. This is achieved by exploiting the formidable combination of multi-object and integral field unit spectroscopic facilities of the ESO Very Large Telescope. As a first step, here we discuss the results obtained for 11 clusters from high and medium resolution spectra acquired through a combination of FLAMES and KMOS observations. We provide the first kinematical characterization of NGC 1261 and NGC 6496. In all the surveyed systems, the velocity dispersion profile declines at increasing radii, in agreement with the expectation from the King model that best fits the density/luminosity profile. In the majority of the surveyed systems we find evidence of rotation within a few half-mass radii from the center. These results are in general overall agreement with the predictions of recent theoretical studies, suggesting that the detected signals could be the relic of significant internal rotation set at the epoch of the cluster's formation.

We report on the discovery of three new millisecond pulsars (namely J1748-2446aj, J1748-2446ak and J1748-2446al) in the inner regions of the dense stellar system Terzan 5. These pulsars have been discovered thanks to a method, alternative to the classical search routines, that exploited the large set of archival observations of Terzan 5 acquired with the Green Bank Telescope over 5 years (from 2010 to 2015). This technique allowed the analysis of stacked power spectra obtained by combining ~206 hours of observation. J1748-2446aj has a spin period of ~2.96 ms, J1748-2446ak of ~1.89 ms (thus it is the fourth fastest pulsar in the cluster) and J1748-2446al of ~5.95 ms. All the three millisecond pulsars are isolated and currently we have timing solutions only for J1748-2446aj and J1748-2446ak. For these two systems, we evaluated the contribution to the measured spin-down rate of the acceleration due to the cluster potential field, thus estimating the intrinsic spin-down rates, which are in agreement with those typically measured for millisecond pulsars in globular clusters. Our results increase to 37 the number of pulsars known in Terzan 5, which now hosts 25% of the entire pulsar population identified, so far, in globular clusters.

We used a dataset of archival Hubble Space Telescope images obtained through the F555W, F814W and F656N filters, to perform a complete search for objects showing $H\alpha$ emission in the globular cluster NGC 6397. As photometric diagnostic, we used the $(V-H\alpha)_0$ color excess in the $(V-H\alpha)_0$-$(V-I)_0$ color-color diagram. In the analysed field of view, we identified 53 $H\alpha$ emitters. In particular, we confirmed the optical counterpart to 20 X-ray sources (7 cataclysmic variables, 2 millisecond pulsars and 11 active binaries) and identified 33 previously unknown sources, thus significantly enlarging the population of known active binaries in this cluster. We report the main characteristics for each class of objects. Photometric estimates of the equivalent width of the $H\alpha$ emission line, were derived from the $(V-H\alpha)_0$-excess and, for the first time, compared to the spectroscopic measurements obtained from the analysis of MUSE spectra. The very good agreement between the spectroscopic and photometric measures fully confirmed the reliability of the proposed approach to measure the $H\alpha$ emission. The search demonstrated the efficiency of this novel approach to pinpoint and measure $H\alpha$-emitters, thus offering a powerful tool to conduct complete census of objects whose formation and evolution can be strongly affected by dynamical interactions in star clusters.

We used a combination of deep optical and Halpha images of the Galactic globular cluster NGC 6440, acquired with the Hubble Space Telescope, to identify the optical counterpart to the accreting millisecond X-ray pulsar SAX J1748.9-2021during quiescence. A strong Halpha emission has been detected from a main sequence star (hereafter COM-SAX J1748.9-2021) located at only 0.15" from the nominal position of the X-ray source. The position of the star also agrees with the optical counterpart found by Verbunt et al. (2000) during an outburst. We propose this star as the most likely optical counterpart to the binary system. By direct comparison with isochrones, we estimated that COM-SAX J1748.9-2021 has a mass of 0.70 Msun - 0.83 Msun, a radius of 0.88 pm 0.02 Rsun and a superficial temperature of 5250pm80 K. These parameters combined with the orbital characteristics of the binary suggest that the system is observed at a very low inclination angle (~8 deg -14 deg) and that the star is filling or even overflowing its Roche Lobe. This, together with the equivalent width of the Halpha emission (~20 Ang), suggest possible on-going mass transfer. The possibile presence of such a on-going mass transfer during a quiescence state also suggests that the radio pulsar is not active yet and thus this system, despite its similarity with the class of redback millisecond pulsars, is not a transitional millisecond pulsar.

About one third of the 3033 $\gamma$-ray sources in the Third Fermi-LAT Gamma-ray Source Catalogue (3FGL) are unidentified and do not have even a tentative association with a known object, hence they are defined as unassociated. Among Galactic $\gamma$-ray sources, pulsars represent the largest class, with over 200 identifications to date. About one third of them are milli-second pulsars (MSPs) in binary systems. Therefore, it is plausible that a sizeable fraction of the unassociated Galactic $\gamma$-ray sources belong to this class. We collected X-ray and optical observations of the fields of twelve unassociated Fermi sources that have been classified as likely MSPs according to statistical classification techniques. To find observational support for the proposed classification, we looked for periodic modulations of the X-ray and optical flux of these sources, which could be associated with the orbital period of a MSP in a tight binary system. Four of the observed sources were identified as binary MSPs, or proposed as high-confidence candidates, while this work was in progress. For these sources, we present the results of our follow-up investigations, whereas for the others we present possible evidence of new MSP identifications. In particular, we discuss the case of 3FGL J0744.1-2523 that we proposed as a possible binary MSP based upon the preliminary detection of a 0.115 d periodicity in the flux of its candidate optical counterpart. We also found very marginal evidence of periodicity in the candidate optical counterpart to 3FGL J0802.3-5610, at a period of 0.4159 d, which needs to be confirmed by further observations.

By exploiting two ACS/HST datasets separated by a temporal baseline of ~7 years, we have determined the relative stellar proper motions (providing membership) and the absolute proper motion of the Galactic globular cluster M71. The absolute proper motion has been used to reconstruct the cluster orbit within a Galactic, three-component, axisymmetric potential. M71 turns out to be in a low latitude disk-like orbit inside the Galactic disk, further supporting the scenario in which it lost a significant fraction of its initial mass. Since large differential reddening is known to affect this system, we took advantage of near-infrared, ground-based observations to re-determine the cluster center and density profile from direct star counts. The new structural parameters turn out to be significantly different from the ones quoted in the literature. In particular, M71 has a core and a half-mass radii almost 50% larger than previously thought. Finally we estimate that the initial mass of M71 was likely one order of magnitude larger than its current value, thus helping to solve the discrepancy with the observed number of X-ray sources.

We present the results of the analysis of deep photometric data for a sample of three Galactic globular clusters (NGC5466, NGC6218 and NGC6981) with the aim of estimating their degree of mass segregation and testing the predictions of analytic dynamical models. The adopted dataset, composed by both Hubble Space Telescope and ground based data, reaches the low-mass end of the mass functions of these clusters from the center up to their tidal radii allowing to derive their radial distribution of stars with different masses. All the analysed clusters show evidence of mass segregation with the most massive stars more concentrated than low-mass ones. The structures of NGC5466 and NGC6981 are well reproduced by multimass dynamical models adopting a lowered-Maxwellian distribution function and the prescription for mass segregation given by Gunn & Griffin (1979). Instead, NGC6218 appears to be more mass segregated than model predictions. By applying the same technique to mock observations derived from snapshots selected from suitable N-body simulations we show that the deviation from the behaviour predicted by these models depends on the particular stage of dynamical evolution regardless of initial conditions.

We present an observational estimate of the fraction and distribution of dark mass in the innermost region of the two Galactic globular clusters NGC 6218 (M12) and NGC 288. Such an assessment has been made by comparing the dynamical and luminous mass profiles derived from an accurate analysis of the most extensive spectroscopic and photometric surveys performed on these stellar systems. We find that non-luminous matter constitutes more than 60% of the total mass in the region probed by our data (R<1.6 arcmin~r_h) in both clusters. We have carefully analyzed the effects of binaries and tidal heating on our estimate and ruled out the possibility that our result is a spurious consequence of these effects. The dark component appears to be more concentrated than the most massive stars suggesting that it is likely composed of dark remnants segregated in the cluster core.

Globular clusters are known to host peculiar objects, named Blue Straggler Stars (BSSs), significantly heavier than the normal stellar population. While these stars can be easily identified during their core hydrogen-burning phase, they are photometrically indistinguishable from their low-mass sisters in advanced stages of the subsequent evolution. A clear-cut identification of these objects would require the direct measurement of the stellar mass. We used the detailed comparison between chemical abundances derived from neutral and from ionized spectral lines as a powerful stellar "weighing device" to measure stellar mass and to identify an evolved BSS in 47 Tucanae. In particular, high-resolution spectra of three bright stars located slightly above the level of the "canonical" horizontal branch sequence in the color-magnitude diagram of 47 Tucanae, have been obtained with UVES spectrograph. The measurements of iron and titanium abundances performed separately from neutral and ionized lines reveal that two targets have stellar parameters fully consistent with those expected for low-mass post-horizontal branch objects, while for the other target the elemental ionization balance is obtained only by assuming a mass of ~1.4Msol, which is significantly larger than the main sequence turn-off mass of the cluster (~0.85 Msol). The comparison with theoretical stellar tracks suggests that this is a BSS descendant possibly experiencing its core helium-burning phase. The large applicability of the proposed method to most of the globular clusters in our Galaxy opens the possibility to initiate systematic searches for evolved BSSs, thus giving access to still unexplored phases of their evolution.

We present multi-wavelength observations of the unassociated gamma-ray source 3FGL J2039.6-5618 detected by the Fermi Large Area Telescope. The source gamma-ray properties suggest that it is a pulsar, most likely a millisecond pulsar, for which neither radio nor $\gamma$-ray pulsations have been detected yet. We observed 3FGL J2039.6-5618 with XMM-Newton and discovered several candidate X-ray counterparts within/close to the gamma-ray error box. The brightest of these X-ray sources is variable with a period of 0.2245$\pm$0.0081 d. Its X-ray spectrum can be described by a power law with photon index $\Gamma_X =1.36\pm0.09$, and hydrogen column density $N_{\rm H} < 4 \times 10^{20}$ cm$^{-2}$, which gives an unabsorbed 0.3--10 keV X-ray flux of $1.02 \times 10^{-13}$ erg cm$^{-2}$ s$^{-1}$. Observations with the Gamma-Ray Burst Optical/Near-Infrared Detector (GROND) discovered an optical counterpart to this X-ray source, with a time-average magnitude $g'\sim 19.5$. The counterpart features a flux modulation with a period of 0.22748$\pm$0.00043 d that coincides, within the errors, with that of the X-ray source, confirming the association based on the positional coincidence. We interpret the observed X-ray/optical periodicity as the orbital period of a close binary system where one of the two members is a neutron star. The light curve profile of the companion star, with two asymmetric peaks, suggests that the optical emission comes from two regions at different temperatures on its tidally-distorted surface. Based upon its X-ray and optical properties, we consider this source as the most likely X-ray counterpart to 3FGL J2039.6-5618, which we propose to be a new redback system.

We used ultra-deep UV observations obtained with the Hubble Space Telescope to search for optical companions to binary millisecond pulsars (MSPs) in the globular cluster 47 Tucanae. We identified four new counterparts (to MSPs 47TucQ, 47TucS, 47TucT and 47TucY) and confirmed those already known (to MSPs 47TucU and 47TucW). In the color magnitude diagram, the detected companions are located in a region between the main sequence and the CO white dwarf cooling sequences, consistent with the cooling tracks of He white dwarfs of mass between 0.15 Msun and 0.20 Msun. For each identified companion, mass, cooling age, temperature and pulsar mass (as a function of the inclination angle) have been derived and discussed. For 47TucU we also found that the past accretion history likely proceeded in a sub-Eddington rate. The companion to the redback 47TucW is confirmed to be a non degenerate star, with properties particularly similar to those observed for black widow systems. Two stars have been identified within the 2-sigma astrometric uncertainty from the radio positions of 47TucH and 47TucI, but the available data prevent us from firmly assessing whether they are the true companions of these two MSPs.

Milli-second pulsars (MSPs) are rapidly spinning neutron stars, with spin periods P_s <= 10 ms, which have been most likely spun up after a phase of matter accretion from a companion star. In this work we present the results of the search for the companion stars of four binary milli-second pulsars, carried out with archival data from the Gemini South telescope. Based upon a very good positional coincidence with the pulsar radio coordinates, we likely identified the companion stars to three MSPs, namely PSRJ0614-3329 (g=21.95 +- 0.05), J1231-1411 (g=25.40 +-0.23), and J2017+0603 (g=24.72 +- 0.28). For the last pulsar (PSRJ0613-0200) the identification was hampered by the presence of a bright star (g=16 +- 0.03) at \sim 2" from the pulsar radio coordinates and we could only set 3-sigma upper limits of g=25.0, r= 24.3, and i= 24.2 on the magnitudes of its companion star. The candidate companion stars to PSRJ0614-3329, J1231-1411, and J2017+0603 can be tentatively identified as He white dwarfs (WDs) on the basis of their optical colours and brightness and the comparison with stellar model tracks. From the comparison of our multi-band photometry with stellar model tracks we also obtained possible ranges on the mass, temperature, and gravity of the candidate WD companions to these three MSPs. Optical spectroscopy observations are needed to confirm their possible classification as He WDs and accurately measure their stellar parameters.

We report on the optical identification of the neutron star burster EXO 1745-248 in Terzan 5. The identification was performed by exploiting HST/ACS images acquired in Director's Discretionary Time shortly after (approximately 1 month) the Swift detection of the X-ray burst. The comparison between these images and previous archival data revealed the presence of a star that currently brightened by ~3 magnitudes, consistent with expectations during an X-ray outburst. The centroid of this object well agrees with the position, in the archival images, of a star located in the Turn-Off/Sub Giant Branch region of Terzan 5. This supports the scenario that the companion should has recently filled its Roche Lobe. Such a system represents the pre-natal stage of a millisecond pulsar, an evolutionary phase during which heavy mass accretion on the compact object occurs, thus producing X-ray outbursts and re-accelerating the neutron star.

We report on the determination of the astrometric, spin and orbital parameters for PSR J1953+1846A, a "black widow" binary millisecond pulsar in the globular cluster M71. By using the accurate position and orbital parameters obtained from radio timing, we identified the optical companion in ACS/Hubble Space Telescope images. It turns out to be a faint (m_F606W>=24, m_F814W>=23) and variable star located at only ~0.06" from the pulsar timing position. The light curve shows a maximum at the pulsar inferior conjunction and a minimum at the pulsar superior conjunction, thus confirming the association with the system. The shape of the optical modulation suggests that the companion star is heated, likely by the pulsar wind. The comparison with the X-ray light curve possibly suggests the presence of an intra-binary shock due to the interaction between the pulsar wind and the material released by the companion. This is the second identification (after COM-M5C) of an optical companion to a black widow pulsar in a globular cluster. Interestingly, the two companions show a similar light curve and share the same position in the color magnitude diagram.

We report on the determination of astrometric, spin and orbital parameters for PSR J1518+0204C, a "black widow" binary millisecond pulsar in the globular cluster M5. The accurate position and orbital parameters obtained from radio timing allowed us to search for the optical companion. By using WFC3/HST images we identified a very faint variable star (m_F390W > 24.8, m_F606W > 24.3, m_F814W > 23.1) located at only 0.25" from the pulsar's timing position. Due to its strong variability, this star is visible only in a sub-sample of images. However, the light curve obtained folding the available data with the orbital parameters of the pulsar shows a maximum at the pulsar inferior conjunction and a possible minimum at the pulsar superior conjunction. Furthermore, the shape of the optical modulation indicates a heating process possibly due to the pulsar wind. This is the first identification of an optical companion to a black widow pulsar in the dense stellar environment of a globular cluster.

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 PSR J1740-5340A in the GC NGC 6397, confirming that it is a deeply peeled star descending from a ~0.8$M_{\odot}$ progenitor. This nicely confirms the theoretical expectations about the formation and evolution of MSPs.

We report on the identification of the optical counterpart to Star1, the exotic object serendipitously discovered by Deutsch et al. in the core of the Galactic globular cluster NGC 6624. Star1 has been classified by Deutsch et al. as either a quiescent Cataclysmic Variable or a low-mass X-ray binary. Deutsch et al. proposed StarA as possible optical counterpart to this object. We used high-resolution images obtained with the Hubble Space Telescope to perform a variability analysis of the stars close to the nominal position of Star1. While no variability was detected for StarA, we found another star, here named COM_Star1, showing a clear sinusoidal light modulation with amplitude \Delta m_F435W~0.7 mag and orbital period of P_orb~98 min. The shape of the light curve is likely caused by strong irradiation by the primary heating one hemisphere of the companion, thus suggesting a quite hot primary.

We used a proper combination of high-resolution HST observations and wide-field ground based data to derive the radial star density profile of 26 Galactic globular clusters from resolved star counts (which can be all freely downloaded on-line). With respect to surface brightness (SB) profiles (which can be biased by the presence of sparse, bright stars), star counts are considered to be the most robust and reliable tool to derive cluster structural parameters. For each system a detailed comparison with both King and Wilson models has been performed and the most relevant best-fit parameters have been obtained. This is the largest homogeneous catalog collected so far of star count profiles and structural parameters derived therefrom. The analysis of the data of our catalog has shown that: (1) the presence of the central cusps previously detected in the SB profiles of NGC 1851, M13 and M62 is not confirmed; (2) the majority of clusters in our sample are fitted equally well by the King and the Wilson models; (3) we confirm the known relationship between cluster size (as measured by the effective radius) and galactocentric distances; (4) the ratio between the core and the effective radii shows a bimodal distribution, with a peak at ~ 0.3 for about 80% of the clusters, and a secondary peak at ~ 0.6 for the remaining 20%. Interestingly, the main peak turns out to be in agreement with what expected from simulations of cluster dynamical evolution and the ratio between these two radii well correlates with an empirical dynamical age indicator recently defined from the observed shape of blue straggler star radial distribution, thus suggesting that no exotic mechanisms of energy generation are needed in the cores of the analyzed clusters.

By using XSHOOTER spectra acquired at the ESO Very Large Telescope, we have studied the surface chemical composition of the companion star to the binary millisecond pulsar PSR J1740-5340 in the globular cluster NGC 6397. The measured abundances of Fe, Mg, Al and Na confirm that the star belongs to the cluster. On the other hand, the measured surface abundance of nitrogen ([N/Fe]=+0.53 +- 0.15 dex) combined with the carbon upper limit ([C/Fe] <-2 dex) previously obtained from UVES spectra allow us to put severe constraints on its nature, strongly suggesting that the pulsar companion is a deeply peeled star. In fact, the comparison with theoretical stellar models indicates that the matter currently observed at the surface of this star has been processed by the hydrogen-burning CN-cycle at equilibrium. In turn, this evidence suggests that the pulsar companion is a low mass (~0.2 Msun) remnant star, descending from a ~0.8 Msun progenitor which lost ~70-80 % of its original material because of mass transfer activity onto the pulsar.

We report on the identification of the optical counterpart to the recently detected INTEGRAL transient IGR J18245-2452 in the Galactic globular cluster M28. From the analysis of a multi epoch HST dataset we have identified a strongly variable star positionally coincident with the radio and Chandra X-ray sources associated to the INTEGRAL transient. The star has been detected during both a quiescent and an outburst state. In the former case it appears as a faint, unperturbed main sequence star, while in the latter state it is about two magnitudes brighter and slightly bluer than main sequence stars. We also detected Halpha excess during the outburst state, suggestive of active accretion processes by the neutron star.

We present the identification of the companion star to the intermediate mass binary pulsar J1439-5501 obtained by means of ground-based deep images in the B, V and I bands, acquired with FORS2 mounted at the ESO-VLT. The companion is a massive white dwarf (WD) with B=23.57+-0.02, V=23.21+-0.01 and I=22.96+-0.01, located at only ~0.05" from the pulsar radio position. Comparing the WD location in the (B, B-V) and (V, V-I) Color-Magnitude diagrams with theoretical cooling sequences we derived a range of plausible combinations of companion masses (1<~Mcom<~1.3 Msun), distances (d<~1200 pc), radii (<~7.8 10^3 Rsun) and temperatures (T=31350^{+21500}_{-7400}). From the PSR mass function and the estimated mass range we also constrained the inclination angle i >~ 55 degrees and the pulsar mass (Mpsr <~2.2 Msun). The comparison between the WD cooling age and the spin down age suggests that the latter is overestimated by a factor of about ten.

By combining high spatial resolution and wide-field spectroscopy performed, respectively, with SINFONI and FLAMES at the ESO/VLT we measured the radial velocities of more than 600 stars in the direction of NGC 6388, a Galactic globular cluster which is suspected to host an intermediate-mass black hole. Approximately 55% of the observed targets turned out to be cluster members. The cluster velocity dispersion has been derived from the radial velocity of individual stars: 52 measurements in the innermost 2", and 276 stars located between 18" and 600". The velocity dispersion profile shows a central value of ~13 km/s, a flat behavior out to ~60" and a decreasing trend outwards. The comparison with spherical and isotropic models shows that the observed density and velocity dispersion profiles are inconsistent with the presence of a central black hole more massive than ~2000 Msol. These findings are at odds with recent results obtained from integrated light spectra, showing a velocity dispersion profile with a steep central cusp of 23-25 km/s at r<2" and suggesting the presence of a black hole with a mass of 17,000 Msol (Lutzgendorf et al. 2011). We also found some evidence of systemic rotation with amplitude Arot ~8 km/s in the innermost 2" (0.13 pc), decreasing to Arot= 3.2 km/s at 18"<r<160".

PSR J1811-1736 (P=104 ms) is an old (~1.89 Gyrs) binary pulsar (P_orb=18.8 d) in a highly eccentric orbit (e=0.828) with an unidentified companion. Interestingly enough, the pulsar timing solution yields an estimated companion mass 0.93 M_{\odot}<M_C<1.5 M_{\odot}, compatible with that of a neutron star. As such, it is possible that PSR J1811-1736 is a double neutron star (DNS) system, one of the very few discovered so far. This scenario can be investigated through deep optical/infrared (IR) observations. We used J, H, K-band images, obtained as part of the UK Infrared Telescope (UKIRT) Infrared Deep Sky Survey (UKIDSS), and available in the recent Data Release 9 Plus, to search for its undetected companion of the PSR J1811-1736 binary pulsar. We detected a possible companion star to PSR J1811-1736 within the 3 sigma radio position uncertainty (1.32 arcsec), with magnitudes J=18.61+/-0.07, H=16.65+/-0.03, and K=15.46+/-0.02. The star colours are consistent with either a main sequence (MS) star close to the turn-off or a lower red giant branch (RGB) star, at a pulsar distance of ~5.5 kpc and with a reddening of E(B-V)~4.9. The star mass and radius would be compatible with the constraints on the masses and orbital inclination of the binary system inferred from the mass function and the lack of radio eclipses near superior conjunction. Thus, it is possible that it is the companion to PSR J1811-1736. However, based on the star density in the field, we estimated a quite large chance coincidence probability of ~0.27 between the pulsar and the star, which makes the association unlikely. No other star is detected within the 3 sigma pulsar radio position down to J~20.5, H~19.4$ and K~18.6, which would allow us to rule out a MS companion star earlier than a mid-to-late M spectral type.

We have used deep V and R images acquired at the ESO Very Large Telescope to identify the optical companion to the binary pulsar PSR J0610-2100, one of the black-widow millisecond pulsars recently detected by the Fermi Gamma-ray Telescope in the Galactic plane. We found a faint star (V~26.7) nearly coincident (\delta r ~0".28) with the pulsar nominal position. This star is visible only in half of the available images, while it disappears in the deepest ones (those acquired under the best seeing conditions), thus indicating that it is variable. Although our observations do not sample the entire orbital period (P=0.28 d) of the pulsar, we found that the optical modulation of the variable star nicely correlates with the pulsar orbital period and describes a well defined peak (R~25.6) at \Phi=0.75, suggesting a modulation due to the pulsar heating. We tentatively conclude that the companion to PSR J0610-2100 is a heavily ablated very low mass star (~ 0.02Msun) that completely filled its Roche Lobe.

We present the first optical observations of the unique system J0737-3039 (composed of two pulsars, hereafter PSR-A and PSR-B). Ultra-deep optical observations, performed with the High Resolution Camera of the Advanced Camera for Surveys on board the Hubble Space Telescope could not detect any optical emission from the system down to m_F435W=27.0 and m_F606W=28.3. The estimated optical flux limits are used to constrain the three-component (two thermal and one non-thermal) model recently proposed to reproduce the XMM-Newton X-ray spectrum. They suggest the presence of a break at low energies in the non-thermal power law component of PSR-A and are compatible with the expected black-body emission from the PSR-B surface. The corresponding efficiency of the optical emission from PSR-A's magnetosphere would be comparable to that of other Myr-old pulsars, thus suggesting that this parameter may not dramatically evolve over a time-scale of a few Myr.

We report on the optical identification of the companion star to the eclipsing millisecond pulsar PSR J1824-2452H in the galactic globular cluster M28 (NGC 6626). This star is at only 0.2" from the nominal position of the pulsar and it shows optical variability (~ 0.25 mag) that nicely correlates with the pulsar orbital period. It is located on the blue side of the cluster main sequence, ~1.5 mag fainter than the turn-off point. The observed light curve shows two distinct and asymmetric minima, suggesting that the companion star is suffering tidal distortion from the pulsar. This discovery increases the number of non-degenerate MSP companions optically identified so far in globular clusters (4 out of 7), suggesting that these systems could be a common outcome of the pulsar recycling process, at least in dense environments where they can be originated by exchange interactions.