Millisecond pulsars are subject to accelerations in globular clusters GCs that manifest themselves in both the first and second spin period time derivatives and can be used to explore the mass distribution of the potentials they inhabit. Here we report on over 20 yr of pulsar timing observations of five millisecond radio pulsars in the core of the corecollapse GC NGC 6752 with the Parkes Murriyang and MeerKAT radio telescopes which have allowed us to measure the proper motions positions and first and second time derivatives of the pulsars. The pulsar timing parameters indicate that all the pulsars in the core experience accelerations and jerks that can be explained only if an amount of nonluminous mass of at least 2.56  10SUP3SUP M SUBSUB is present in the core of NGC 6752. On the other hand our studies highly disfavor the presence of an intermediatemass black hole at the center of the cluster with a mass equal to or greater than 3000M SUBSUB.
['Millisecond pulsars', 'Pulsars', 'Globular star clusters', '1062', '1306', '656', 'Astrophysics - High Energy Astrophysical Phenomena']
Timing of Millisecond Pulsars in NGC 6752. III. On the Presence of Nonluminous Matter in the Clusters Core
{'No Header': '<!-- image -->', "Timing of Millisecond Pulsars in NGC 6752. III. On the Presence of Nonluminous Matter in the Cluster ' s Core": 'A. Corongiu 1 , A. Ridol /uniFB01 1,2 , F. Abbate 1,2 , M. Bailes 3,4 , A. Possenti 1 , M. Geyer 5,6 , R. N. Manchester 7 , M. Kramer 2,8 ,  \n1  \nP. C. C. Freire 2 , M. Burgay 1 , S. Buchner 5 , and F. Camilo 5  \nINAF  \n-  \nOsservatorio Astronomico di Cagliari, Via Della Scienza 5, I-09047 Selargius, Italy; alessandro.corongiu@inaf.it  \n2  \nMax-Planck-Institut für Radioastronomie, Auf dem Hügel 69, D-53121 Bonn, Germany  \n3 Centre for Astrophysics and Supercomputing, Swinburne University of Technology, PO Box 218, Hawthorn, Vic, 3122, Australia  \n4 The ARC Centre of Excellence for Gravitational Wave Discovery ( OzGrav ) , Mail H29, Swinburne University Of Technology, PO Box 218, Hawthorn, Vic 3122,  \nAustralia  \n5 South African Radio Astronomy Observatory, 2 Fir Street, Black River Park, Observatory 7925, South Africa  \n6 Department of Astronomy, University of Cape Town, Rondebosch, Cape Town, 7700, South Africa 7 Australia Telescope National Facility, CSIRO, Space and Astronomy, Epping NSW 1710, Australia 8 Jodrell Bank Centre for Astrophysics, School of Physics and Astronomy, Univ. of Manchester, Manchester, M13 9PL, UK  \nReceived 2023 August 1; revised 2024 June 25; accepted 2024 July 1; published 2024 September 10', 'Abstract': 'Millisecond pulsars are subject to accelerations in globular clusters ( GCs ) that manifest themselves in both the /uniFB01 rst and second spin period time derivatives, and can be used to explore the mass distribution of the potentials they inhabit. Here we report on over 20 yr of pulsar timing observations of /uniFB01 ve millisecond radio pulsars in the core of the core-collapse GC NGC 6752 with the Parkes ( Murriyang ) and MeerKAT radio telescopes, which have allowed us to measure the proper motions, positions, and /uniFB01 rst and second time derivatives of the pulsars. The pulsar timing parameters indicate that all the pulsars in the core experience accelerations and jerks that can be explained only if an amount of nonluminous mass of at least 2.56 × 10 3 M e is present in the core of NGC 6752. On the other hand, our studies highly disfavor the presence of an intermediate-mass black hole at the center of the cluster, with a mass equal to or greater than ∼ 3000 M e .  \nUni /uniFB01 ed Astronomy Thesaurus concepts: Millisecond pulsars ( 1062 ) ; Pulsars ( 1306 ) ; Globular star clusters ( 656 )', '1. Introduction': "Globular clusters ( GCs ) are among the oldest and densest structures in a galaxy and often contain several × 10 5 stars in spherical structures with radii of just a few tens of parsecs. In recent decades, more and more powerful optical telescopes have been investigating these systems, providing new insights into their mass distribution and internal structure ( see, e.g., Baumgardt & Hilker 2018; Leitinger et al. 2023, and references therein ) . Today accurate censuses of their stellar populations can be obtained via optical data, and their structure, dynamical status, and evolution can be explored with N -body simulations. While their macroscopic properties are increasingly well understood, their dense cores remain dif /uniFB01 cult to probe. Even high spatial resolution instruments like GAIA ( Gaia Collaboration et al. 2016 ) struggle to resolve stars in the central regions due to the high number densities and distances. Another reason is the presence in the central regions of high-mass nonluminous stellar remnants like neutron stars ( NS ) and possibly stellar mass black holes ( Bahcall & Wolf 1976, Bahcall & Wolf 1977 ) , and also massive underluminous white dwarfs ( e.g., Vitral et al. 2023 ) which have migrated to the clusters ' cores, thus leaving the lightest objects at larger distances from the center due to mass segregation. This scenario /uniFB01 nds con /uniFB01 rmation via the observed spatial distribution of millisecond pulsars ( MSP ) in GCs, which are on average among the heaviest objects in a cluster and usually are located near the cores. The pulsars located outside the core have been postulated to be the  \n<!-- image -->  \nOriginal content from this work may be used under the terms of the Creative Commons Attribution 4.0 licence. Any further distribution of this work must maintain attribution to the author ( s ) and the title of the work, journal citation and DOI.  \nresult of multibody ejection mechanisms, some involving the presence of an intermediate-mass black hole ( IMBH ) of several tens of solar masses ( e.g., Colpi et al. 2002 ) , located near the cluster ' s center of gravity. In recent works, instead, it has been argued that three-body interactions involving stellar mass objects are a viable mechanism for expelling binary MSPs from the cluster core ( see, e.g., Leigh et al. 2024 ) .  \nThe MSPs in the GC NGC 6752 show characteristics that are typical of the populations of core-collapsed GCs: a large predominance of isolated pulsars ( for other examples, see the pulsar population of M15, Anderson 1993, NGC 6624, Abbate et al. 2022 and NGC 6517, Pan et al. 2021 ) and a much larger fraction of pulsars at large distances from the core. As discussed in detail by Verbunt & Freire ( 2014 ) , these facts are likely causally related: the large number of stellar interactions per binary in the core region -to which NSs migrate via dynamical friction -implies that many MSP binaries are disrupted in three-body encounters, with many being ejected from the cluster, and with several others barely staying bound to it. The ejections of NSs are more likely if they are involved in interactions with more massive stellar remnants.  \nHowever, even in comparison with the pulsar populations of other core-collapsed GCs, the population of NGC 6752 is extreme in these respects. Its peculiarities are hardly explainable without invoking the presence of a nonnegligible amount of under- or nonluminous ( hereafter, we name all of that as nonluminous ) mass in the core, under the form of massive remnants. PSR J1910 -5959A ( hereafter PSRA, D ' Amico et al. 2001 ) , the /uniFB01 rst pulsar discovered in this cluster and its only binary pulsar known to date, has the largest offset from the center of its host GC in core radius units, namely ∼ 58 r c ( this fact has cast doubts on the  \nassociation, which has, however, been demonstrated recently by Corongiu et al. 2023 ) and angular units ( ~¢ 6.3 ) ;whilePSR J1910 -5959C ( hereafter PSRC, D ' Amico et al. 2002,hereafter Paper I ) has the third angular offset in core radius units, namely ∼ 23 r c, and fourth in angular units ( ∼ 2 6 9 ) . For these two pulsars, recent episodes of ejection from the core are a very likely possibility, since dynamical friction should have taken them back to the core in less than ∼ 1 Gyr ( Colpi et al. 2002 ) . The other known pulsars, namely PSR J1910 -5959B, PSR J1910 -5959D, and PSR J1910 -5959E ( hereafter PSRB, PSRD, and PSRE respectively, Paper I ) , and the recently discovered PSR J1910 -5959F ( hereafter PSRF, Ridol /uniFB01 et al. 2021 ) are located in the core of NGC 6752. We will refer to these four pulsars as the core pulsars .  \nIn Paper I it was shown that the measured spin period derivatives of PSRB, PSRD, and PSRE are dominated by the acceleration imparted to them by the cluster potential well, and at least 1.3 × 10 4 M e should be present in the form of nonluminous objects in the core. The measured spin period derivative of PSRF ( Ridol /uniFB01 et al. 2021 ) is of the same order of magnitude as the other three core pulsars, and also its value is likely due to the gravitational effect of nonluminous matter in the core of NGC 6752.  \nIn this work, we report updated timing solutions for PSRB, PSRC, PSRD, PSRE, and PSRF 10 based on observations taken over ∼ 20 yr with the Parkes Murriyang and MeerKAT radio telescopes, and investigate the total mass and distribution of the nonluminous matter in the cluster core. This paper is organized as follows: In Section 2 we describe the observations and the data analysis, and in Section 3 we comment on the proper motions of the core pulsars. In Section 4 we place constraints on the nonluminous matter in the cluster core, and in Section 5 we discuss our results in terms of the central mass-to-light ratio. In Section 6 we comment on the main sources of uncertainty in our analysis, and we review our numerical results. Finally, in Section 7 we summarize our /uniFB01 ndings.", '2.1. Parkes Observations': "The GC NGC 6752 was observed with the Parkes Murriyang 64 m radio telescope, located in Australia, from 1999 September to 2016 March. Observations were carried out in the L -band with the Multibeam and the H -OH receivers, depending on availability, with a bandwidth of 256 MHz and using two orthogonal polarizations that were detected and summed before digitization. The Stokes I signal was processed with the 1 bit Analog Filterbank ( AFB; see, e.g., D ' Amico et al. 2001 ) digital signal processor ( hereafter ' backend ' ) up to its decommissioning in 2012 September. The AFB acquired the detected signal at the central frequency of 1390 MHz, split the observed bandwidth into 512 frequency channels, and recorded the data with a sampling time of 80 /uni03BC s ( 125 /uni03BC s in some earlier observations ) . In all observations taken after 2012 September, and in some before, the data were recorded in search mode with the Pulsar Digital Filterbank 3 ( PDFB3 ) and 4 ( PDFB4; see, e.g., Manchester et al. 2013 and references therein ) . These backends acquired the detected signal at the central frequency  \nof 1369 MHz, split the observed bandwidth in 512 ( DFB3 ) or 1024 ( DFB4 ) frequency channels, and recorded the data with a sampling time of 80 /uni03BC s ( 125 /uni03BC s in a few observations ) .  \nWe reduced and analyzed the raw data with the same method used and explained in detail in Paper I and Corongiu et al. ( 2006, hereafter Paper II ) . Brie /uniFB02 y, for each pulsar, we folded the data at the approximate topocentric period with the software package dspsr 11 ( van Straten & Bailes 2011 ) , as predicted by the ephemeris published in Paper II. We folded the data with a typical subintegration length of 1 minute and kept the same number of frequency channels as the raw data to aid in accurate dedispersion. The resulting reduced data, whose /uniFB01 les are usually described as ' timing archives, ' have the form of a matrix whose elements are the binned pro /uniFB01 le amplitudes observed in each frequency channel, one for each subintegration time, and some metadata ( header ) information.", '2.2. MeerKAT Observations': "The MeerKAT radio telescope ( Jonas & MeerKAT Team 2016; Camilo et al. 2018 ) , located in South Africa, is an array composed of 64 offset-Gregorian dishes, each equivalent to a circular aperture of approximately 13.5 m across. We used MeerKAT and its UHF-band ( 544 -1088 MHz ) and L -band ( 856 -1712 MHz ) receivers to observe NGC 6752 on several occasions between 2019 July and 2022 January, as part of the MeerTime 12 ( Bailes et al. 2020 ) and the TRansients And PUlsars with MeerKAT ( TRAPUM 13 ; Stappers & Kramer 2016 ) Large Survey Projects. Each of these projects has its own backend, tailored to its main scienti /uniFB01 c goals. MeerTime, which is mostly devoted to precision pulsar timing and polarimetry, makes use of the Pulsar Timing User Supplied Equipment ( PTUSE ) backend. This is capable of acquiring up to four Nyquist-sampled tied-array beams with up to 4096 coherently dedispersed frequency channels with four Stokes parameters. The data were recorded as 8 bit search mode PSRFITS /uniFB01 les with a time resolution of 9.57 /uni03BC sat L -Band and 7.53 /uni03BC sat UHF. Given that only pulsars A and C are located far from the core of the cluster, while all the others are less than 0 2 away from the nominal center of the cluster, three PTUSE beams ( centered on pulsar A, on pulsar C, and on the nominal center of NGC 6752, respectively ) were suf /uniFB01 cient to record the signals from all the known pulsars. Also, being the DMs of the pulsars all in the range of 33.20 -33.70 pc cm -3 , a single search mode /uniFB01 le coherently dedispersed at a DM = 33.3 pc cm -3 was suf /uniFB01 cient to allow the redetection of all the pulsars without any loss of sensitivity. TRAPUM, on the other hand, is focused on the search for pulsars and radio transients, and uses two computing clusters. First, the Filterbanking BeamFormer User Supplied Equipment ( FBFUSE ) cluster applies the beamforming technique to combine the raw signals from all the antennas and synthesize hundreds ( up to 288 for our NGC 6752 observations ) tied-array beams on the sky ( Barr 2018; Chen et al. 2021 ) . These are then recorded as search mode ' /uniFB01 lterbank ' /uniFB01 les by the Accelerated Pulsar Search User Supplied Equipment ( APSUSE ) cluster: the observing band is recorded with a typical time resolution of ∼ 60 -80 /uni03BC s and is split into 4096 frequency channels. The /uniFB01 ne channelization is particularly important in this case to remove the effects of  \nTable 1 List of the MeerKAT Observations of NGC 6752 Used for This Work  \nNote. The symbols stand for the following: f c, central frequency; /uni0394 f , observing bandwidth; CD DM: dispersion measure used for coherent dedispersion; N chan, number of frequency channels; N pol, number of Stokes parameters; t samp, sampling time; N ant, number of antennas.  \ninterstellar dispersion, as APSUSE does not bene /uniFB01 t from coherent dedispersion, due to computational constraints. The majority of the MeerKAT observations of NGC 6752 were carried out using PTUSE and APSUSE simultaneously, so as to take advantage of their complementary characteristics. The exact used setups depended on the main scienti /uniFB01 c purpose of each observing session, as well as on the availability of new observing modes ( such as the possibility of recording more than a single beam with PTUSE, or to record a different number of channels between PTUSE and APSUSE during a simultaneous MeerTime + TRAPUM session ) which were gradually implemented over the course of these projects. We excluded the APSUSE data for PSRB, because the time stamp was not correctly recorded in the corresponding data /uniFB01 les. We report all the MeerKAT observations of NGC 6752 used in this paper, along with their con /uniFB01 gurations, in Table 1.", '2.3. ToA Extraction and Timing': "We extracted the pulse times of arrival ( ToAs ) with the routines of the software suite psrchive ( Hotan et al. 2004 ) , by coherently adding in phase the pro /uniFB01 les in each archive with respect to subintegrations and channels, and convolving them with a high signal-to-noise ratio ( S / N ) template obtained by summing in phase the observed pro /uniFB01 les with the best S / N. For each pulsar, we produced a speci /uniFB01 c template for each telescope-backend-observing mode combination. The pulsars were often subject to strong scintillation, due to the ionized interstellar medium, which produces random variations of the pulse S / N with respect to both time and frequency, correlated on some timescale. For this cluster the timescale is such that amplitude variations occur not only from one observation to another, but also within single observations, and they are so dramatic that pulses can be very bright at some moments and / or in some frequency subbands, and completely undetectable in others. For these reasons, we could not determine a single optimal integration length and subband frequency width along which the pro /uniFB01 les in the archive should always be summed; instead, we visually inspected every single archive to determine the best time interval ( s ) and frequency band ( s ) where pulses were evaluated as detections. Whenever possible, we extracted more than one ToA from a single observation, either by decimating in time, frequency, or both.  \nWe determined the instrumental time jump between the different backends used for Parkes observations, and between the Parkes and MeerKAT data sets as a whole. We determined these jumps for each pulsar separately, since they are not due to technical reasons only, but mainly because templates differ, from one pulsar to another, in shape, displacement, and time extent. We used the ToAs obtained from the Parkes AFB backend as reference ( i.e., we determined the time jumps with respect to ) , since for all pulsars these sets span the largest time interval and contain the largest number of ToAs. At /uniFB01 rst we considered the ToAs for each telescope separately, and /uniFB01 tted them against the best timing model for each set. We took the resulting χ 2 and multiplied the ToA uncertainties by its square root. Once this operation is performed, a new /uniFB01 t on the same ToAs against the same timing model should return a reduced χ 2 = 1. We then combined the two telescope ' s ToAs and /uniFB01 tted them for a single time jump alone, using Parkes ToAs as the reference. Finally, we proceeded with the /uniFB01 t against the pulsar timing model, still allowing this time jump to vary for correctly estimating the uncertainties in the timing model ' s parameters.  \nWe /uniFB01 tted the obtained ToAs against the rotational and astrometric parameters with the pulsar timing software TEMPO2 ( Hobbs et al. 2006 ) , using the DE430 solar system ephemeris from the NASA Jet Propulsion Laboratory 14 for barycentering. Table 2 reports our timing results. The values for each parameter in the timing model are reported with the nominal TEMPO2 1 σ uncertainty on the last signi /uniFB01 cant digit ( in parentheses ) . The main improvement in the timing solutions with respect to Paper II is the determination of the second derivative of the spin period for all pulsars, the proper motions for PSRB, D, and E, and the re /uniFB01 nement of the latter for PSRF. Parameters whose measurement has already been reported in previous works are now determined with a similar or higher precision. Figure 1 displays the corresponding timing residuals. The amplitude of the plotted error bars is equal to the measured one for each ToA, multiplied by the amount necessary to obtain χ 2 = 1 in the preliminary /uniFB01 t of each telescope data set separately ( see above ) . No clear trend is evident in the  \nFigure 1. Timing residuals vs. MJD. ToAs obtained with Parkes and MeerKAT are plotted with blue triangles and red squares, respectively. If not visible, error bars are smaller than the marking symbol. Units and scales are the same for all pulsars for an immediate comparison of their timing rms and noise level. Pulsar names and rms values are indicated near the bottom right corner of the dedicated panel.  \n<!-- image -->  \nTable 2 Measured and Derived Parameters for the Pulsars in NGC 6752", 'Notes.': "- a Units of R.A. are hours, minutes, and seconds, and units of decl. are degrees, arcminutes, and arcseconds.\n- b Reference epoch for both the spin period and the position.\n- c The offset of the pulsars is calculated with respect to the position of the cluster ' s center of gravity reported by Ferraro et al. ( 2003; see Table 3 ) .  \nTable 3 Positional, Kinematical, and Structural Parameters for NGC 6752 Published in Literature and Used in This Work  \nReferences. ( a ) Ferraro et al. ( 2003 ) ; ( b ) Vasiliev & Baumgardt ( 2021 ) ; ( c ) Baumgardt & Vasiliev ( 2021 ) ; ( d ) Globular clusters online catalog https: // people.smp.uq.edu.au / HolgerBaumgardt / globular / , v.3, see Hilker et al. ( 2020 ) ; ( e ) Noyola & Gebhardt ( 2006 ) ; ( f ) Baumgardt & Hilker ( 2018 ) ; ( g ) Libralato et al. ( 2022 ) .  \nresiduals of any pulsar, thus indicating that the adopted timing models adequately describe the available data sets.", '3. Proper Motions of the Core Pulsars and the Central Velocity Dispersion': "The pulsars in the core of NGC6752 have excellent proper motion determinations, namely with a signi /uniFB01 cance of 11 σ ,25 σ ,31 σ ,and11 σ for PSRB, PSRD, PSRE, and PSRF in alphabetical order, respectively. The average motion of the pulsars results in m d á ñ=- a cos 3.15 0.51 mas yr -1 ( quoted uncertainties are everywhere at the 1 σ level unless explicitly indicated ) and 〈 /uni03BCδ 〉 = -4.13 ± 0.31 mas yr -1 , in excellent agreement with the cluster motion m d =- a cos 3.162 0.023 mas yr -1 and /uni03BCδ = -4.028 ± 0.022 mas yr -1 , obtained from GAIA optical observations ( all cluster ' sparameterspicked up from literature and used in this work are presented in Table 3 with all necessary references, unless explicitly indicated ) .Atthe cluster distance of 4.125 ± 0.040 kpc the relative projected transverse velocity of each pulsar with respect to the cluster center, for which the GAIA proper motion is assumed, are 15.5 ± 4.7, 7.7 ± 4.1, 4.8 ± 3.6, and 15.4 ± 11.1 km s -1 respectively, thus implying that ( barring a large radial velocity component ) all pulsars are compatible with being bound to the cluster core, since the escape velocity from the core is V esc = 34.5 km s -1 .  \nWe estimated the 1D velocity dispersion of the core pulsars, by calculating the standard deviation of the proper motion components in R.A. and decl. separately, then averaging the values in the two directions. We obtain σ V,PSRs = 8.0 ± 1.9 km s -1 . This value is compatible with the value σ 0 = 8.5 ± 0.2 km s -1 , resulting from the measured proper motion dispersion in the core of the cluster, σ/uni03BC = 0.436 ± 0.009 mas yr -1 , using optical observations.", '4. Presence of Nonluminous Matter in the Core of NGC6752': "Two pulsars in the core of the cluster, PSRB and PSRE, have a negative measured time derivative of their spin period  P ; meas that is the signature of their accelerated motion in the gravitational potential well of the cluster ( Paper I ) , since  P is always intrinsically positive because of pulsar spindown. In the case of PSRD ( Paper I ) and PSRF ( Ridol /uniFB01 et al. 2021, this work )  P meas is positive, with an absolute value at least as large as for PSRB. An insight into the Galactic population of MSPs 15 ( here de /uniFB01 ned as pulsars with a spin period not larger than 10 ms ) shows that it is highly unlikely that  P meas for PSRD and PSRF is dominated by their intrinsic spindown  P intr .  \nIn fact, the Galactic population of MSPs consists of 228 objects for which it has been measured  > P 0 . Among them, 223 MSPs have  < -P 10 19 ss -1 , with an average spindown rate  á ñ= ' -P 1.40 10 20 ss -1 , and a standard deviation  s = ' -1.22 10 P 20 ss -1 . For the four pulsars in the core,  P meas /uni2223 /uni2223 is larger than  á ñ P by factors that range from 31 ( in the case of PSRE ) to 69 ( in the case of PSRD ) . The maximum  P in this sample is 7.74 × 10 -20 ss -1 ( PSR J0218 + 4232, Desvignes et al. 2016 ) , a value that is an order of magnitude smaller than  P meas of PSRF, namely the smallest value among the core pulsars for which  > P 0 meas . For four of the remaining /uniFB01 ve Galactic MSPs,  P is comprised between 1.05 × 10 -19 ss -1 ( PSRB1937 + 21, Reardon et al. 2021 ) and 1.63 × 10 -19 ss -1 ( PSRJ1850 + 0242, Scholz et al. 2015 ) . But even in the case of PSRJ1850 + 0242, the  P value is still a factor of ∼ 4.5 smaller than  P meas for PSRF. Moreover, if one includes these four MSPs in the aforementioned calculation of the average spindown, one obtains  á ñ= ' -P 1.60 10 20 ss -1 and  s = P ' -2.00 10 20 ss -1 . It is immediately seen that the  P of these four MSPs is more than 4 σ away from the average. The /uniFB01 fth one is PSR J1402 + 13 ( P = 5.9 ms,  = ' -P 4.8 10 17 ss -1 , Abdollahi et al. 2022 ) , about which no discussion seems to be present in the literature. Because these /uniFB01 ve high  P objects represent the 2.2% only of the sample of Galactic MSPs, and their  P is signi /uniFB01 cantly higher than the average MSPs spindown rate, we evaluated these /uniFB01 ve objects as outliers in the  P distribution of the Galactic MSPs, and assumed that the distribution of the intrinsic spindown of the galactic MSPs is well represented by the 223 objects for which  < < -P 0.0 10 19 ss -1 . It could be objected that an MSP in a GC might have a  P intr that is much higher than the average of the Galactic MSPs, since the cluster ' s crowded environment might have induced a different formation path for the hosted MSPs, which results in a considerably large  P intr . Nevertheless, NGC6752 also hosts PSRA, for which Corongiu et al. ( 2023 ) obtained  = ' -P 5.02 10 intr 21 ss -1 , and PSRC, whose  P meas and projected angular separation from the cluster ' s center of gravity imply   ' -P 1.71 10 intr 20 ss -1 . Since for both PSRA and PSRC  P intr is fully consistent with the average of the Galactic /uniFB01 eld MSPs, we can safely assume that this is also true for the pulsars in the core of NGC 6752.  \nThe effects of the cluster environment on the core pulsars are also evident in the second time derivative of their spin periods. In fact, the upper limits, both observational and theoretical, on the intrinsic P /uni0308 for MSPs are so low, that it can safely be  \nassumed that the intrinsic P /uni0308 for the core pulsars are completely negligible in comparison to the measured ones, whose origin can thus be attributed to the cluster ' s dynamics.  \nIn Paper I it was investigated whether the acceleration imparted to the core pulsars ( known at that time ) could be entirely due to the observed luminous mass, or if a further amount of matter was necessary to explain their  P meas . The authors of Paper I used the argument of the mass-to-light ratio   in the central regions ( see Equation ( 1 ) in Paper I ) , and found that a minimum   of an order of 10, namely >   9 for PSRB and E, and >   13 for PSRD, is required to impart to these pulsars the accelerations, whose component along the line of sight produces the measured  P . Hereafter in this paper, in order to improve the reading of the text, the term acceleration will indicate the component along the line of sight of an acceleration, unless explicitly speci /uniFB01 ed otherwise. These   values were much larger than the one available at that time and obtained with optical observations, =   1.1 ( Pryor & Meylan 1993 ) . The authors of Paper I calculated the amount of mass that should be present in the core of the cluster in the form of nonluminous objects, and found a lower limit of 1.3 × 10 4 M e .", '4.1. Our Approach': 'Given the high precision of our determination of both the /uniFB01 rst and second derivatives of the spin period for the core pulsars, we revisited the estimate of the nonluminous mass in the core of NGC 6752 with a different approach. Phinney ( 1993 ) detailed how  P meas is related to the intrinsic spindown  P intr and all accelerations acting on a pulsar in a GC with the following relation 16 :  \n  = + + + P P P P a c a c a c 1 l l l meas intr ,MW ,SHK ,GC ( )  \nwhere c is the speed of light, al ,MW is the acceleration imparted by the Milky Way, 17 a l ,SHK is the apparent acceleration due to the pulsar proper motion, also known as the Shklovskii ( 1970 ) effect, and al ,GC is the acceleration imparted by the GC. When the second derivative of the spin period is also measured, a further equation is available, that links the measured P /uni0308 to the time derivative of the imparted acceleration, often referred to as the jerk. In fact, taking the time derivative of Equation ( 1 ) and considering relevant terms only, one obtains:  \n  = + + P P P P a c a c 2 l l meas intr ,GC ,NN /uni0308 /uni0308 ( )  \nwhere we have introduced the effect of the nearest neighbor star,  al ,NN , which cannot be in general neglected in crowded environments like GCs ( see, e.g., Abbate et al. 2019 ) , while we dropped the time derivatives of the Milky Way acceleration and that of the Shklowskii effect. These two terms can be neglected because, across the epoch range spanned by the observations, the 3D motion of the cluster in the Milky Way does not produce signi /uniFB01 cant changes on its position, hence on the acceleration imparted by the Galaxy. The pulsar proper motion  \nchanges even less, hence the apparent acceleration due to its transverse motion remains substantially constant.  \nThe acceleration and jerk imparted by the cluster can thus be determined by solving Equation ( 1 ) with respect to al ,GC and Equation ( 2 ) with respect to  al ,GC :  \n/uni23A1 /uni23A3 /uni23A2 /uni23A4 /uni23A6 /uni23A5   = ---a c P P P P a a 3 l l l ,GC meas intr ,MW ,SHK ( )  \n/uni239C /uni239F /uni239C /uni239F /uni23A1 /uni23A3 /uni23A2 /uni239B /uni239D /uni239E /uni23A0 /uni239B /uni239D /uni239E /uni23A0 /uni23A4 /uni23A6 /uni23A5   = --a c P P P P a 4 l l ,GC meas intr ,NN /uni0308 /uni0308 ( )  \nand the link with the cluster structural parameters is obtained by replacing a l ,GC in Equation ( 3 ) and  al ,GC in Equation ( 4 ) with the expression that results by assuming a given mass distribution for the cluster. In the simplest case of a stationary and spherically symmetric distribution for the mass of the cluster, al ,GC takes the well-known simple form:  \n=-a GM r r r 5 l l ,GC 3 ( ) ( )  \nwhere M ( r ) is the cluster mass enclosed in a sphere of radius r , centered at the center of gravity of the cluster, and G is the Gravitational constant ( 6.67 × 10 -11 m 3 kg -1 s -2 ) . We express a generic position, with respect to the cluster center, in terms of two coordinates perpendicular to the line of sight, a a d º -a r cos PSR GC GC ( ) and r δ ≡ δ PSR -δ GC, and rl , namely the pulsar distance from the center along the line of sight; the relation = + + a d r r r r l 2 2 2 obviously holds. In these de /uniFB01 nitions, α PSR and δ PSR are the celestial coordinates of a given pulsar, while α GC and δ GC are the celestial coordinates of the cluster gravity center, and rl is de /uniFB01 ned so that its value is positive for positions farther than the cluster center with respect to the observer, and negative otherwise. At this stage, we considered r α and r δ of each pulsar as known exactly, despite the uncertainties on the assumed coordinates of the center of gravity of the cluster ( 0 5 in both coordinates, Ferraro et al. 2003 ) . In Section 6 we will discuss the impact of these uncertainties on our results.  \nThe analytic expression for  al ,GC is obtained by taking the time derivative of Equation ( 5 ) :  \n/uni239B /uni239D /uni239E /uni23A0 /uni239B /uni239D /uni239E /uni23A0  = -= ---r v r v a d dt GM r r r G dM r dr r r GM r r v r r 3 6 l l l l ,GC 3 4 3 l 2 ( ) ( ) · ( ) · ( )  \nwhere v is the pulsar velocity with respect to the cluster center, whose components are m d º = -a a a v dr dt cos ,PSR PSR ( m d a D cos ,GC GC ) and v δ ≡ dr δ / dt = ( /uni03BCδ ,PSR -/uni03BCδ ,GC ) D , where in turn D is the cluster distance, m d a cos ,PSR PSR and /uni03BCδ ,PSR are the proper motion components of a pulsar, while m d a cos ,GC GC and /uni03BC δ ,GC are the proper motion components of the cluster, and v l ≡ dr l / dt . Given the above de /uniFB01 nition of r l , v l is positive when the pulsar is moving away from the observer, and negative otherwise. We did not consider the quantities v α and v δ to be exactly known, since their uncertainty is comparable to their amplitude and to the central velocity dispersion.  \nThe problem is so far underdetermined: for a given pulsar, only two equations are available, but the unknown terms are its  \ndepth inside the cluster rl , its 3D velocity v , and all the parameters that enter in the analytic expression of a given mass distribution for the cluster. We indicate the set of all of these parameters with the formal vector s ≡ ( s 1, s 2, K , sN ) . Moreover, the contribution to  al ,GC due to the nearest neighbor calls into play some further parameters that describe the cluster structure and dynamic properties at the pulsar position ( see below ) .', '4.2. The Bayesian Analysis': "The availability of more than one pulsar does not allow us to fully constrain the problem, since each object has its own position and velocity, and also, the cluster ' s local properties differ from a given position to another. For this reason we adopted a Bayesian approach. All terms in the right-hand side of Equation ( 3 ) but the intrinsic spindown ( see below ) are known. We can then assume that they obey a Gaussian probability distribution with a mean value equal to their value, either measured or derived, and a standard deviation equal to their 1 σ uncertainty. This in turn implies that also al ,GC ( hereafter al for notation simplicity ) follows a Gaussian probability distribution:  \n/uni23A7 /uni23A8 /uni23A9 /uni23AB /uni23AC /uni23AD s p s = --Pa a a 1 2 exp 1 2 7 l a l l a ,m 2 2 l l ,m ,m ( ) ( ) ( )  \nwhose mean value al ,m is given by Equation ( 3 ) , and whose standard deviation s al ,m is given by the usual rules for the propagation of the uncertainties. Once al is substituted with the explicit expression given by Equation ( 5 ) , where in turn M ( r ) is substituted with the analytical expression that corresponds to a chosen model for the mass distribution, Equation ( 7 ) returns the probability distribution s P r , a l ,X l ( ) , for the structural parameters in s and the still unknown pulsar depth in the cluster rl ,asit results by considering a single pulsar alone, identi /uniFB01 ed by the formal index X ( the subscript al in the symbol s P r , a l ,X l ( ) indicates that this probability distribution comes from the analysis of the acceleration, in order to distinguish it from the distribution that results from the analysis of the jerk ) .  \nWe calculated the Milky Way acceleration by applying Equation ( 16 ) in Lazaridis et al. ( 2009 ) , but using for the vertical component of the Galactic acceleration F z the analytic formula provided by Li & Widrow ( 2021 )( Equation ( 14 )) . The adopted values for the Solar motion in the Galaxy are Θ e = 240.5 ± 4.1 km s -1 and R e = 8.275 ± 0.034 kpc ( Gravity Collaboration et al. 2021 ) . At the cluster ' s distance and position ( l g = 336.4929 deg, b g = -25.628 deg in Galactic coordinates ) we obtain a l ,MW = ( + 2.7 ± 0.6 ) × 10 -11 ms -2 . The Shklovskii effect gives a contribution = al ,SHK m  ' --7.52 0.07 10 5 mas yr 11 1 2 ( ) ( ) ms -2 , where /uni03BC is the proper motion amplitude, for a pulsar at the distance of NGC6752. As already seen above,  P intr can be as high as several 10 -20 ss -1 , a value that would give a ∼ 5% contribution on the observed accelerations. Therefore, we calculated a l ,m as  = --a cP P a a l l l ,m meas ,SHK ,MW , and we took into account our poor knowledge of the intrinsic spindown by adding in quadrature the quantity   s 'áñ+ c P P P ( ) to the uncertainty on al ,m that results from the standard propagation of the uncertainties on the three known terms.  \nWe applied the method illustrated above to Equation ( 4 ) , for deriving the probability distribution obeyed by  al ,GC ( hereafter  al for notation simplicity ) , which in turn translates into the  \nprobability distribution  Pa ,X l for the structural parameters in s and the unknown quantities for a generic pulsar X. In this case it must be taken into account that the nearest neighbor contribution on the measured jerk cannot be calculated but only statistically treated and that, most important, it does not follow a Gaussian distribution, but a Lorentzian one ( see, e.g., Abbate et al. 2019 and references therein ) :  \n    p = + Pa a a a 1 . 8 NN 0 NN 2 0 2 ( ) ( )  \nThe quantity  a 0 is dubbed as characteristic jerk , and it is related to the cluster ' s structure ( Prager et al. 2017 ) by the following relation:  \n px s = á ñ a G n m 2 3 9 0 v ( )  \nwhere ξ ; 3.04 is a numerical constant, while σ v, n , and á ñ m are the cluster ' s 1D velocity dispersion, star number density and mean stellar mass at the pulsar position. The product áñ n m can be set equal to the mass density at the pulsar position, while σ v can be only constrained to be smaller than the central 1D velocity dispersion. Once a model is assumed for the cluster ' s mass distribution, the probability distribution for  a NN depends on the pulsar position and the local velocity dispersion σ v. Therefore,  Pa ,X l depends on the parameters in s , the pulsar depth in the cluster, the pulsar 3D velocity in the cluster frame, and the cluster ' s 1D velocity dispersion at the pulsar position. Because of the peculiar distribution for  a NN ,  s s v P r , , , a l ,X v l ( ) must be calculated from more general principles of statistics. Let α and β be two generic quantities, which follow probability distributions P α ( α ) and P β ( β ) respectively, and γ = α -β their difference. The probability distribution for γ , P ( γ ) , is given by the following integral 18 :  \nò g a a g a = -a b -¥ +¥ P P P d . 10 ( ) ( ) ( ) ( )  \nIn our case a = cP P /uni0308 , β is indeed  a NN , and γ is the value for  al as predicted by Equation ( 6 ) once a given mass distribution is assumed. Because of the measurement of the second-order time derivative of the spin period, P /uni0308 follows a Gaussian distribution with mean value P meas /uni0308 , and standard deviation s P meas /uni0308 . As a direct consequence, once P is kept /uniFB01 xed, also cP P /uni0308 follows a Gaussian distribution with mean value cP P meas /uni0308 and standard deviation s c P P meas /uni0308 . We completely neglected the contribution due to the intrinsic evolution of the spin period, namely P intr /uni0308 since, as we already commented on, it is negligible in comparison with P meas /uni0308 .  \nThe joint probability of simultaneously having an acceleration al and a jerk  al acting on a given pulsar X is the product of the two above-mentioned distributions. Its marginalization with respect to rl , v and σ v returns the probability distribution P X ( s ) for the parameters that describe the assumed mass model:  \n s s s = ' s s s v v v P P r P r P r P P dr d d , , , , 11 a l a l l l X ,X ,X v v 3 v l l /uni222D ( ) ( ) ( ) ˜ ( ) ˜ ( ) ˜ ( ) ( )  \nwhere the P ˜ functions are the prior probability distributions for their arguments. We assumed that Prl ˜ ( ) is /uniFB02 at and nonzero only for those values for r l that place a given pulsar within one core radius from the cluster center, i.e.,   -^ r r r 0 l c 2 2 /uni2223 /uni2223 , where we de /uniFB01 ned º + a d ^ r r r 2 2 , and whose sign is opposite to  P meas , namely r l > 0if  < P 0 meas and vice versa. In fact, both a l ,MW and a l ,SKH are positive, and their sum is so smaller than  c P P meas /uni2223 /uni2223 for all core pulsars, that Equation ( 1 ) is satis /uniFB01 ed only if al ,GC has the same sign of  P meas , and the explicit expression for a l ,GC, as given by Equation ( 5 ) , states that the sign of r l must be opposite to a l ,GC, hence to  P meas .  \nWe assumed v P ˜ ( ) to be the product of the priors a Pv ˜ ( ) , d Pv ˜ ( ) , and Pvl ˜ ( ) for the three components of the velocity. Thanks to our measurement of the proper motion for all core pulsars, we assumed that a Pv ˜ ( ) and d Pv ˜ ( ) are Gaussian distributions, whose mean and standard deviation can be obtained from their de /uniFB01 nition and the usual uncertainty propagation rules. Because no information is available about the radial velocity of the core pulsars, we assumed Pvl ˜ ( ) to be /uniFB02 at and nonzero in the range -50 km s -1 /uni0084 v l /uni0084 50 km s -1 , where the limits are a conservative rounding of the escape velocity from the core, and we allowed both positive and negative values since it is not a priori known whether a pulsar is moving away or toward the observer in the cluster frame. Finally, we assumed s P v ˜ ( ) to be /uniFB02 at in the range 0 /uni0084 σ v /uni0084 σ 0, since it cannot be larger than the central velocity dispersion. Its upper limit will be discussed later, since it is model dependent.  \nBecause the same mass distribution must be responsible for the accelerations and the jerks observed in all the core pulsars, the /uniFB01 nal probability distribution P ( s ) is given by the product of the distributions P X ( s ) , later multiplied by the prior distribution s P ˜ ( ) :  \n/uni23A7 /uni23A8 /uni23A9 /uni23AB /uni23AC /uni23AD  = s s s P P P . 12 X X ( ) ( ) ˜ ( ) ( )", '4.3. The Mass Models and the Distribution of the Luminous Mass': 'The models that we considered for the mass distribution in the core of the cluster are the sum of two components. The /uniFB01 rst one is due to the luminous mass, whose spatial distribution is described by the model that best /uniFB01 ts to the data taken with optical observations. We assumed it was relatively well known and thus we kept it /uniFB01 xed. The second component is due to the nonluminous mass, and we investigated two scenarios for its distribution. In the /uniFB01 rst one we simply assumed that all the nonluminous mass is generically contained in a sphere, whose radius is smaller than the true distance of the closest pulsar to the cluster center ( see Section 4.4 ) . This scenario allowed us to obtain a /uniFB01 rst estimate of the amount of nonluminous matter in the core of the cluster. In the second scenario, instead, we assumed that an IMBH is located at the center of gravity of the cluster. This scenario differs from the /uniFB01 rst one ( see Section 4.5 ) , since the presence of an IMBH has a substantial impact on the mass distribution in the core, and it allowed us to test whether the resulting picture represents a physically meaningful description of the mass distribution in the core of the cluster.  \nThe distribution of the luminous mass obviously plays a key role in our work. In fact, if one adopts a model that, e.g., underestimates its content, one obtains an overestimation of the  \nnonluminous mass, and vice versa. For this reason, we paid great attention to the previous works, available in the literature, on the mass distribution of NGC 6752, and we identi /uniFB01 ed a model that allows one to obtain a consistent picture of the optically observed structure of this cluster.  \nThe most recent study of the mass surface and radial pro /uniFB01 le of NGC 6752 was reported in Baumgardt & Hilker ( 2018 ) . The authors performed N -body simulations of the formation and evolution of the Galactic GCs. The results of their simulations for NGC 6752 ( see their Figure E13 ) are consistent with the data at angular radii θ ⊥ /uniF089 10 " , but they clearly underestimate the observed surface density pro /uniFB01 le at smaller radii. This evidence can be considered consistent with the results obtained by Ferraro et al. ( 2003 ) , who studied the surface star density pro /uniFB01 le of NGC 6752 using different models. In particular, they obtained a very good description of their data at all angular radii by using a combination of two King ( 1966 ) density pro /uniFB01 les 19 ( see text and Figure 5 in Ferraro et al. 2003 ) , for separately describing the cluster star density in its inner and outer regions, in projection on the sky. We recall that a King pro /uniFB01 le describes the structure of a GC in terms of its mass density, whose analytic expression is given by the following equation:  \n/uni239C /uni239F /uni23A1 /uni23A3 /uni23A2 /uni239B /uni239D /uni239E /uni23A0 /uni23A4 /uni23A6 /uni23A5 r r = + -r r r 1 13 K 0 c 2 3 2 ( ) ( )  \nwhere ρ 0 is the central mass density and r c is the core radius. It is interesting to note that the transition between the two pro /uniFB01 les occurs at θ ⊥ ∼ 10 " , i.e., at about the same angular radius at which the simulations by Baumgardt & Hilker ( 2018 ) begin to agree with the data. Since we have assumed that all core pulsars are indeed located in the core of the cluster, i.e., their true 3D distance from the center of gravity of the cluster is smaller than the core radius, which in turn is smaller than the angular radius at which the outer King pro /uniFB01 le begins to describe the data, we can safely assume that the luminous mass distribution that contributes to the acceleration imparted on the pulsars is well described by the inner King pro /uniFB01 le only. We tested whether the the double King pro /uniFB01 le by Ferraro et al. ( 2003 ) adequately describes the cluster pro /uniFB01 le at all radii, by calculating the cluster total mass that results from all the above-mentioned assumptions, and the already known cluster \' s parameters ( see Table 3 ) . We set the value for the core radius of the outer King pro /uniFB01 le to r c = 28 " , according to the results by Ferraro et al. ( 2003 ) , and we normalized it by imposing that at r = 10 " it must return the same density given by the inner one. Finally, we integrated this two-component density pro /uniFB01 le up to the radius r NGC 6752 = 12.26 pc, which corresponds to 4.27 projected half-light radii, thus following the prescription given by Leitinger et al. ( 2023; see their Section 4.1 for a discussion about this choice ) . We obtained a total mass of the cluster M TOT = 2.865 × 10 5 M e , a value just a few percent higher than the value ( 2.76 ± 0.04 ) × 10 5 M e ( see Table 3 ) , thus ensuring that the double King pro /uniFB01 le adequately describes the optically observed mass pro /uniFB01 le of the cluster at all radii.', '4.4. A First Estimate of the Nonluminous Mass in the Core': "Our /uniFB01 rst estimate of the nonluminous mass M NL in the cluster core is based on assuming that it is entirely contained in a sphere whose radius is smaller than the distance of the closest pulsar to the cluster center. We simpli /uniFB01 ed the problem by assuming that the nonluminous mass does not perturb the distribution of the luminous matter. With these assumptions, this distribution for the nonluminous mass can be mathematically treated as a point-mass M NL placed at the center of gravity of the cluster. For this reason, this model will later be referred to as the point-mass model. The formal vector s is thus the 1D vector s ≡ ( M NL ) . We assumed the prior distribution PM NL ˜ ( ) to be /uniFB02 at and nonzero in the range 0 /uni0084 M NL /uni0084 10 4 M e . In this model we set to 10 km s -1 the upper limit for the local velocity dispersion that enters in the characteristic jerk expression. This value is a conservative upper limit for the measured central velocity dispersion, namely σ 0 = 8.5 ± 0.2 km s -1 . We obtained  = ' -+ M M 3.05 10 NL 0.12 0.10 3 at the 1 σ level (  = ' -+ M M 3.05 10 NL 0.23 0.21 3 and = ' -+ M 3.05 NL 0.37 0.35  M 10 3 at the 2 σ and 3 σ level, respectively ) . This means that about 3 thousand solar masses of nonluminous mass, with a 3 σ lower limit M NL /uni0085 2.68 × 10 3 M e , are required to explain the measured accelerations and jerks acting on the core pulsars.", '4.5. Is There an IMBH in the Core of NGC 6752?': 'We explored a scenario where an IMBH resides at the center of gravity of NGC 6752. The presence of an IMBH modi /uniFB01 es the density pro /uniFB01 le of a GC, with respect to the predictions of a pure King pro /uniFB01 le, up to a radius r i dubbed as the IMBH in /uniFB02 uence radius ( Baumgardt et al. 2004a and references therein ) . Such a density distribution is called cusp , since for r /uni0084 r i it predicts higher values for the density than the ones predicted at the same radii by the King model that describes the density pro /uniFB01 le at r > r i. In this scenario the nonluminous mass is not only due to the IMBH alone, but also to the extra mass  \nOutside r i the mass density pro /uniFB01 le follows the standard King ( 1966 ) law ( see Equation ( 13 )) . In the simulations by Baumgardt et al. ( 2004b ) the radial density pro /uniFB01 le is a continuous function at r i. Therefore the condition r r r = = + r r r r 1 C i K i 0 i c 232 ( ) ( ) [ ( ) ] must hold. Baumgardt et al. ( 2004a ) also cited the further condition that the total mass of the cusp should not be larger than the mass of the IMBH, if the latter is heavier than a few percent of the cluster mass. From the point-mass model we obtain M NL /uniF088 3500 M e at the 3 σ level, i.e., about 1.4% of the cluster mass. Under the reasonable assumption that the true amount of nonluminous matter cannot be larger than twice the aforementioned upper limit, this constraint is not necessary in this case. A visual inspection of Figure 5 in Ferraro et al. ( 2003 ) led us to place an upper limit on r i. In fact, if the cusp is present, it should leave a signature in the optically observed radial pro /uniFB01 le of the cluster. In our inspection of the aforementioned /uniFB01 gure, we could not identify any deviation from the King pro /uniFB01 le that describes the data at angular radii up to 10 " . We thus deduce that, if indeed a cusp is present in the core, its radius of in /uniFB02 uence cannot be larger than the x -coordinate x 1 of the /uniFB01 rst data point in the plot, which we extrapolated to be = x log arcseconds 0.1 10 1 ( ) , i.e., 0.025 pc at the cluster distance. We conservatively assumed that the cusp can extend up to the x -coordinate x 2 of the second data point, which we extrapolated to be = x log arcseconds 0.6 10 2 ( ) , i.e., 4 0, which corresponds to 0.08 pc at the cluster distance. This value is roughly equal to the projected distance of PSRB from the cluster \' s center of gravity, the core pulsar with the second highest projected distance. One can easily verify that in the cases of PSRB, PSRD, and PSRE, the parameter space that results from all the considerations above also contains a subspace where these pulsars are located inside the cusp, thus ensuring that these three objects can be used as probes of its structure.  \nThe resulting cluster \' s mass pro /uniFB01 le M ( r ) , including the IMBH, is then:  \n/uni23A7 /uni23A8 /uni23AA /uni23A9 /uni23AA  ò ò ò p r p r p r = = + < = + + Mr Mr M r rdr r r Mr M r rdr r rdr r r 4 if 4 4 if 16 r r r r BH 0 2 C i BH 0 2 C 2 K i i i ( ) ( ) ( ) ( ) ( ) ( ) ( )  \nM CUSP -M King ( r i ) , namely the difference between the amounts of the cluster mass contained up to r = ri , as predicted by the cusp density law ( M CUSP ) , and the King model that describes the luminous mass at radii r /uni0085 r i ( M King ( r i )) . Prager et al. ( 2017 ) explicitly obtained the analytic expression for the density pro /uniFB01 le of the cusp ρ C ( r ) , as it results from the works by Baumgardt et al. ( 2004a, 2004b ) :  \n/uni239C /uni239F /uni239B /uni239D /uni239E /uni23A0 r r = -r r r . 14 C i i 1.55 ( ) ( )  \nThe in /uniFB02 uence radius r i is related to the IMBH mass M BH and the 1D velocity dispersion of the stars in the cluster core σ 0 by the formula ( Baumgardt et al. 2004a, Equation ( 1 )) :  \ns = r GM 2 . 15 i BH 0 2 ( )  \nwhere the King pro /uniFB01 le that describes the mass distribution at r /uni0085 r i is the inner King pro /uniFB01 le from the studies by Ferraro et al. ( 2003; see Section 4.3 ) . In this model the still unconstrained structural parameter is M BH, hence s = ( M BH ) , and we assumed its prior PM BH ˜ ( ) to be /uniFB02 at and nonzero in the conservative range 0 /uni0084 M BH /uni0084 10 4 M e . We also included σ 0 among the parameters of our probability space, with a Gaussian prior whose mean and standard deviation are its measured value and 1 σ uncertainty, respectively. We obtained = M BH  \' -+ M 2.67 10 0.08 0.05 3 at the 1 σ level (  \' -+ M 2.67 10 0.16 0.12 3 and  \' -+ M 2.67 10 0.24 0.18 3 at the 2 σ and 3 σ levels, respectively ) , while the overall amount of the nonluminous mass, M NL = M BH + M CUSP -M King ( ri ) , now results in = \' -+ M 3.08 NL 0.06 0.08  M 10 3 at the 1 σ level (  \' -+ M 3.08 10 0.14 0.14 3 and \' -+ 3.08 0.25 0.20  M 10 3 at the 2 σ and 3 σ levels, respectively ) . Figure 2 displays the triangular plot for the posterior probability of the M BH -σ 0  \nFigure 2. Posterior probability for the M BH -σ 0 space resulting from the analysis aimed to constrain the possible presence of an IMBH in the core of NGC 6752. The horizontal and vertical scales have been zoomed in the ranges 2 × 10 3 M e -4 × 10 3 M e and 8 km s -1 -9kms -1 , respectively, for a better inspection of the relevant portion of each plot. The diagonal black line in the colored plot delimits the observational constraint on the size of the IMBH in /uniFB02 uence radius. The scale for the 2D color map, displayed in the vertical left panel, is normalized to the maximum value of the 2D probability distribution. Units for all probability distributions are arbitrary.  \n<!-- image -->  \nspace. The black line in the color plot marks the boundary between the two regions where the observational constraint on the IMBH radius of in /uniFB02 uence is ( left side ) and is not ( right side ) satis /uniFB01 ed. We immediately see that the points with the highest probability density lie at the boundary of the allowed region, and not inside it, and this leads us to put into question the physical validity of our result, i.e., if there is indeed a ∼ 2.7 × 10 3 M e IMBH at the center of the cluster, or whether it is simply a mathematical solution for our problem with no real physical meaning.  \nFigure 3 also displays the posterior probability of the M BH -σ 0 space, but in the case where the constraint on the IMBH radius of in /uniFB02 uence is not applied. The points with the highest probability are located in the not allowed region, and those points, that show in Figure 2 the highest probability, have now a probability not negligibly lower than the maximum. Quantitatively speaking, the probability that the true solution of the problem lies in the allowed region, i.e., that the IMBH scenario has a real physical meaning, is ∼ 0.25 only. This probability is not low enough for an immediate ruling out of this model, nor high enough to consider physically certain the presence of an IMBH in the core of NGC 6752. We could in principle relax the upper limit on the radius of in /uniFB02 uence, allowing r i to be larger than the 4 0 limit assumed above. This would shift, in Figure 3, the diagonal line toward larger values for M BH, thus increasing the probability that the IMBH scenario has a real physical meaning. But a cusp with such a  \nlarge radius would already have been detected in the analysis by Ferraro et al. ( 2003 ) , at odds with their results.  \nIf we, instead, put a more stringent upper limit on r i, assuming that its angular size at the cluster distance is not larger than the /uniFB01 rst data point of Figure 5 in Ferraro et al. ( 2003 ) , i.e., r i /uni0084 x 1 = 1 25, further considerations arise against the physical validity of the IMBH scenario. It must be noted that, with such a tight constraint on r i, the cusp radius would be smaller than the projected distance of PSRE from the cluster \' s gravity center. We recall here that PSRE is the core pulsar with the smallest angular distance from the cluster center. This means that none of the core pulsars would be located inside the cusp, hence none of them could probe the cusp structure. As a consequence, the IMBH model would be totally indistinguishable from the point-mass model, and it should return exactly the same results for the overall amount of the nonluminous mass. It is easy to check whether this was the case, namely if we could obtain an overall amount for the nonluminous of at least 2.66 × 10 3 M e ( 3 σ lower limit, see Section 4.4 ) , by applying our methods with this new tighter constraint on r i.If we conservatively impose that the central velocity dispersion σ 0 lies in the range 6.5 -10.5 km s -1 , i.e., within 10 σ from its measured value of 8.5 km s -1 ( see Table 3 ) , the overall a priori predicted amount of the nonluminous mass can be at most 1.3 × 10 3 M e , in clear contradiction with the results given by the point-mass model. The IMBH model can, again a priori, predict M NL /uni0085 2.66 × 10 3 M e , only if the cluster \' s central  \nFigure 3. Same as Figure 2, but also considering those points in the M BH -σ 0 space, whose corresponding value for the IMBH in /uniFB02 uence radius is larger than the upper limit we imposed on r i ( see Section 4.5 ) . Labels have been added for marking and easily identifying in which region our condition on r i is obeyed ( \' ALLOWED \' ) or violated ( \' NOT ALLOWED \' ) . In both the uppermost and rightmost panels, the black solid and red dashed lines are the posterior probability distributions that are obtained by considering, in the marginalization, the allowed region only and the entire M BH -σ 0 space, respectively. All 1D posteriors have been plotted before being normalized after the marginalization, in order to highlight their relative height.  \n<!-- image -->  \nvelocity dispersion is at least 15.1 km s -1 , which is too large with respect to the one inferred from the optical observations.  \nAll of these considerations point against the presence in the core of the cluster of an IMBH, either with a mass equal to or higher than few thousands solar masses, or anyway massive enough to be responsible for the presence of the extra mass that explains the measured derivatives of the spin period for the core pulsars. The results from the point-mass model /uniFB01 rmly state that an amount of nonluminous mass of at least ∼ 3 × 10 3 M e is necessary to justify the observed accelerations and jerks acting on the core pulsars, but it is very unlikely that the IMBH model is able to provide this amount of mass, consistently with all other optically observed features of NGC 6752. Most likely, the only possible scenario remaining invokes the presence in the cluster core of a population of much lighter nonluminous objects, whose nature and spatial distribution cannot be investigated with the present data, and which are so far undetected in optical observations and not yet predicted by models and simulations published by other authors.', '5. The Central Mass-to-light Ratio of NGC 6752': "The mass-to-light ratio (   ) in GCs provides a useful link between the directly measurable emitted light and the indirectly determinable distribution of mass, thus adding elements in the investigation of the structure of this kind of stellar associations, and probing models obtained by theoretical calculations and numerical simulations.  \nFollowing the method reported in Paper I, we revisited the mass-to-light ratio argument for PSRB, D and E, also including PSRF. Equation ( 1 ) in Paper I reads:  \n/uni23A1 /uni23A3 /uni23A2 /uni23A4 /uni23A6 /uni23A5    q q q p q q < < = ' S< ^ ^ ^ ^ -^ --  P P a c G c M D L 1.1 5.1 10 10 pc s 17 l V V V ,max cyl 2 2 18 4 , 2 1 ( ) ( ) ( ) ( ) ( )  \nwhere θ ⊥ is the projected angular offset of the pulsar with respect to the cluster center, q < ^ al ,max ( ) is the maximum acceleration, due to the cluster potential that an object can experience at the angular offset θ ⊥ , M cyl ( < θ ⊥ ) is the cluster mass enclosed in a cylinder of radius θ ⊥ whose main axis is parallel to the line of sight and passes through the cluster center, Σ V ( < θ ⊥ ) is the cluster ' s average surface brightness within the angular radius θ ⊥ from its center, and LV , e is the luminosity of the Sun in the V -band. Using the value Σ V = 4.025 × 10 4 L e pc -2 , assumed constant according to the work reported in Noyola & Gebhardt ( 2006 ) , we obtained >   4.6, 5.2, 4.6, and 4.3 ( here and hereafter in this section in solar units for PSRB, PSRD, PSRE, and PSRF, in the given order ) . The discrepancies between these values and the ones reported in Paper I for the /uniFB01 rst three pulsars, namely    9 for PSRB and PSRE, and    13 for PSRD, are due to our use of more recent determinations of the surface  \nbrightness pro /uniFB01 le and the coordinates for the cluster center ( see Table 3 ) . The position of the cluster center of gravity, assumed by the authors of Paper I, also led them to obtain the much larger lower limit of 1.3 × 10 4 M e for the mass responsible for the accelerations acting on PSRB, PSRD, and PSRE, than the one we obtained in our analysis ( see Section 4 ) .  \nWe also derived the   values by considering the luminous mass distribution for the cluster reported and commented in Section 4. We integrated the double King density pro /uniFB01 le ( see Section 4.3 ) along cylinders of radii equal to the angular separation of each core pulsar, and parallel to the line of sight. At /uniFB01 rst, we considered the luminous mass only, and we found =   2.0 for PSRB, PSRD, and PSRE, and 1.9 for PSRF. These values are substantially smaller than those we obtained by applying Equation ( 17 ) . We then added a mass of 2680 M e in the core, namely the 3 σ lower limit on the nonluminous mass that results by applying the point-mass model, and we obtained   = 5.65, 8.15, 9.31, and 4.35, thus con /uniFB01 rming that Equation ( 17 ) may give a lower limit that is too conservative for   , hence for the amount of mass responsible for the acceleration acting on a pulsar in a GC.  \nIt is worth mentioning a plot of the radial pro /uniFB01 le of   for NGC 6752, recently published in the online catalog of the Galactic GCs ( see reference ' d ' in Table 3 ) . This pro /uniFB01 le was obtained by Baumgardt ( 2017 ) by /uniFB01 tting the observed surface brightness pro /uniFB01 le against N -body simulations of the evolution of GCs, including the formation of nonluminous objects. We extrapolated from this plot the   values at the projected distance r ⊥ of each core pulsar, obtaining   = 3.1, 3.2, 3.3, and 2.6. These values are larger than the ones we obtained by considering the luminous mass only, yet still smaller than the ones implied by the  P P ratios and, consequently, than the ones we obtained considering our conservative lower limit on the amount of nonluminous matter. These discrepancies clearly are nonnegligible, and seem to mean that the methods used in Baumgardt ( 2017 ) still lead to an underestimation of the mass in the core of this cluster, in the form of objects that are not detectable in optical observations.", '6. The Sources of Uncertainty in Our Results': 'In our Bayesian analysis we encountered three sources of uncertainty for our results. The /uniFB01 rst one is our poor knowledge of the intrinsic  P , for the core pulsars, which we addressed with a correction of the uncertainty in the measured accelerations. The second one resides in the treatment of the contribution due to the nearest neighbor to the time derivatives of the accelerations, a contribution that, as we already commented on, can be only statistically addressed. The third one is due to the uncertainty on the coordinates of the center of gravity of NGC 6752, which we have not taken into account so far. In this section we discuss the impact of these three sources of uncertainty on our results.', '6.1. The Intrinsic Spindown': "We treated our poor knowledge of the intrinsic spindown of the core pulsars in terms of an uncertainty to be added in quadrature to the one on all other terms on the right-hand side of Equation ( 3 ) . As we already discussed in Section 4.2,webased our quanti /uniFB01 cation of this additional term on the observed distribution of the measured  P for the MSPs in the Galactic /uniFB01 eld. Table 4 details all contributions to the uncertainty in the acceleration due to the GC, i.e., the uncertainty on the value one  \nobtains by calculating al from Equation ( 3 ) , in terms of the 1 σ fractional uncertainty. The term due to the intrinsic spindown dominates the uncertainty on a l by 2 orders of magnitude, with respect to the Galactic acceleration and the Shklovskii effect, and by 3 orders of magnitude with respect to the uncertainty on the measured  P . Nevertheless, the resulting total fractional uncertainty on a l is just a few percent, with an average of 3.9%. The resulting fractional uncertainty on the nonluminous mass can be estimated as follows. Equation ( 5 ) implies that the fractional uncertainty on a l is the same as on M TOT ( r ) , being the latter the total amount of mass responsible for the observed acceleration, i.e., the sum of the amount of luminous mass M L with the amount of nonluminous mass M NL. A reasonable reference value for the amount of the luminous mass can be estimated by calculating the mass predicted by the King pro /uniFB01 le we used to describe the luminous mass in the core of the cluster ( see Section 4.3 ) , within a radius 〈 r 3D 〉 equal to the average 3D distance of the core pulsars from the center of gravity. We estimated it as á ñº á ñ= ^ r r r 3 2 0.663 3D ,PSR c , where 〈 r ⊥ , PSR 〉 is the average angular displacement of the core pulsars from the coordinates of the cluster ' s gravity center. We thus obtained M L ( 〈 r 3D 〉 ) = 460 M e . By using the value M NL = 3050 M e , namely the most probable value for the amount of nonluminous mass we obtained by applying the point-mass model, it results in M TOT ( 〈 r 3D 〉 ) = 3510 M e , whose 3.9% error is 137 M e .Because we considered known the distribution of the luminous mass, this value represents the contribution  s M P , NL to the uncertainty on the nonluminous mass that is induced from the poor knowledge of the intrinsic spindown of the core pulsars.", '6.2. The Nearest Neighbor Contribution to the Jerks': "Table 5 details all contributions to the uncertainty in the jerk due to the GC potential, namely the uncertainty on the value one obtains by calculating  al from Equation ( 4 ) , again in terms of the 1 σ fractional uncertainty. All columns, but the one labeled   s a a l NN , present the same quantities as the corresponding ones in Table 4, but for the case of the jerk instead of the acceleration, and P /uni0308 instead of  P . The contribution to the uncertainty due to the intrinsic P /uni0308 is not present since, as we already showed in Section 4, this quantity is several orders of magnitude lower than the uncertainty on P meas /uni0308 , hence it can be completely neglected. The parameter  s a NN quanti /uniFB01 es the uncertainty on the contribution to the jerk due to the nearest neighbor star. We de /uniFB01 ned it so that the probability to have    s a a NN NN /uni2223 /uni2223 is equal to the probability that a generic Gaussian distributed quantity x is within ± 1 σ x from its mean x m, namely:  \n/uni23A7 /uni23A8 /uni23A9 /uni23AB /uni23AC /uni23AD       ò ò p ps s + = ' --s s s s -+ -+ a a a da x x dx 1 1 2 exp 1 2 18 x x x x 0 0 2 NN 2 NN m 2 2 a a x x NN NN m m ( ) ( )  \nfrom which one obtains   s º a 1.837 a 0 NN . We used the reference value  = ' --a 3.2 10 ms 0 20 3 , which results by using in Equation ( 9 ) the distance á ñ= r r 0.663 3D c from the cluster center, as in Section 6.1, σ V = σ V,PSR = 8.0 km s -1 , namely the 1D velocity dispersion of the pulsars for the local velocity dispersion at =á ñ r r 3D ( see Section 3 ) , and the approximation r áñ» á ñ n m r K 3D ( ) , where ρ K ( r ) is given by Equation ( 13 ) .  \nIn this case the contribution due to the nearest neighbor dominates the overall uncertainty in  al , while the uncertainty on  \nTable 4 Contributions to the Uncertainties on the Pulsar Accelerations due to the GC PotentialTable 5  \nContributions to the Uncertainties on the Pulsar Jerks due to the GC Potential  \nP meas /uni0308 plays a nearly negligible role. The resulting contribution  s M a , NL NN on the uncertainty on the nonluminous mass cannot be obtained as simply as in the case of the /uniFB01 rst derivatives, but requires some additional considerations given the complexity of the expression for  al ( see Equation ( 6 )) . The most generic form of the relation between  s a a , l NN , namely the contribution on the uncertainty on  al due to the poor knowledge of  a NN , and  s M a , NL NN can be written as 20 :  \n   s s = ¶ ¶ a M 19 a a l M a , NL , l NN NL NN ( )  \nafter keeping /uniFB01 xed to a reference value all quantities in Equation ( 6 ) but M NL. Assuming   s h = a a a l , l NN , where η = 0.5404 is the average of the values in the fourth column of Table 5, the contribution on the uncertainty on M NL, due to the poor knowledge of  a NN only, results:  \n   s h = ¶ ¶ -a M a . 20 M a l l , NL 1 NL NN /uni2223 /uni2223 ( )  \nIt must be noted that the explicit expressions for  al and  ¶ ¶ a M l NL also depend on the orientation, with respect to the observer, of both the pulsar position r and velocity v in the cluster ' s center reference frame. Because these dependencies cannot be neglected, we averaged Equation ( 20 ) over all possible orientations of the two vectors r and v , thus obtaining:  \n s h á ñ = -r dM r dr Mr 3 3 . 21 M a , NL NN ( ) ( ) ( )  \nWe used again the results from the point-mass model ( see Section 4.4 ) , and the reference position =á ñ r r 3D , as above, thus obtaining  s á ñ = M a , NL NN 2193 M e .  \nAny comparison between the values we obtained for  s M P , NL and  s á ñ M a , l NL must consider that the overall probability distribution for M NL is the product of the two distributions that result by separately considering the acceleration and jerk. One can easily see that the product of two peaked distributions, that are peaked at values that are close to each other but with different widths of the peak, results in another function that is  \npeaked at the average value of the peaks of the two factors, with a characteristic width which is smaller than the one of both factors. Nevertheless, the value obtained for  s á ñ M a , NL NN seems so large that one may deduce that the probability distribution obtained by considering the jerks only might not be able to tighten the uncertainties on the nonluminous mass that is obtained by considering the accelerations only. For this reason we repeated the analysis reported in Section 4.4, but considering only the accelerations. We obtained M NL = ( 3.105 ± 0.137 ) × 10 3 M e at the 1 σ level ( = -+ M 3.105 NL 0.313 0.332  ' M 10 3 and  = ' -+ M M 3.105 10 NL 0.508 1.035 3 at the 2 σ and 3 σ level, respectively ) . The best value is consistent, within the uncertainties, with the one reported in Section 4.4, but in this case, the uncertainties are notably larger than the ones obtained by including the jerks. In practice, the amount of nonluminous mass is mainly determined by the accelerations, while the jerks allow a nonnegligible tightening of the uncertainties.", "6.3. The Cluster ' s Gravity Center and the Final Estimate of the Amount of Nonluminous Mass": 'All quantitative results, reported and commented so far, rely on the assumed position for the cluster \' s gravity center, which we held /uniFB01 xed at the values indicated in Table 3 and hereafter referred to as the measured coordinates, whose 1 σ uncertainties are 0 5 in both directions. Given their nonnegligible uncertainties, we investigated whether our estimation of the nonluminous mass might change if the center of gravity were not placed at the currently bona /uniFB01 de coordinates. We considered positions for the gravity center with an angular displacement r off /uni0084 1 " . For each position, we recomputed the amount of nonluminous matter by applying the point-mass model. In several cases, the resulting probability distribution showed more than one single peak for M NL. For this reason we consider in this discussion the mean mass M NL,mean ( α , δ ) so de /uniFB01 ned:  \nò ò ad ad ad = M M P M dM PM dM , , , , , 22 NL,mean NL NL NL NL NL ( ) ( ) ( ) ( )  \nwhere P ( M NL, α , δ ) is the probability density distribution obtained for the gravity center located at the celestial coordinates ( α , δ ) . We obtained values in the range 2.87 × 10 3 M e /uniF088 M NL,mean /uniF088 7.54 × 10 3 M e , with mean 〈 M NL,mean 〉 = 3.9 × 10 3 M e , standard deviation  s = \' á ñ M 1.0 10 M 3 NL,mean , and median  = \' M M 3.6 10 NL,mean 3 ˜ .  \nFigure 4 ( panel ( a )) displays a histogram of the obtained values, with a bin width of 200 M e . There is a clear peak around 3000 M e , in agreement with our previous result, but also a nonnegligible tail extendingupto ∼ 7600 M e .Figure4 also contains information about the r off values, in steps of 0 25, from which one deduces that the tail above 5000 M e receives contributions from positions at r off /uni0085 0 5, and that the peak at 3000 M e is mainly due to positions with an offset in the range 0 25 /uni0084 r off /uni0084 0 75, with also  \nFigure 4. Histograms of the amount of nonluminous matter in the core of NGC 6752 for the explored positions of the cluster center of gravity. Panel ( a ) : cumulative histogram: for each bin different colors represent the contribution at distances r off, from the coordinates of the cluster center of gravity by Ferraro et al. ( 2003 ) , r off < 0 25 ( red ) ,0 25 /uni0084 r off < 0 5 ( dark blue ) ,0 5 /uni0084 r off < 0 75 ( orange ) ,0 75 /uni0084 r off < 1 0 ( light blue ) , respectively. Panels ( b ) to ( e ) : separate histograms for each annulus de /uniFB01 ned as above.  \n<!-- image -->  \nFigure 5. Amount of the nonluminous mass in the core of NGC 6752 as a function of the position of the cluster center of gravity. Panel ( a ) : region containing the explored positions for the cluster center of gravity ( small colored circle ) , compared to the position of the core pulsars ( black bullets ) , and the core radius ( blue dashed circle ) . Panel ( b ) : a color map for the amount of the nonluminous mass in the core of the cluster as a function of the position of the cluster center of gravity. For each exploited position, the plotted value is the mean value of the posterior probability distribution ( see Section 6.3 ) . Dotted -dashed circles delimit the regions whose points have offsets r off /uni0084 0.5 σ , 0.5 σ < r off /uni0084 1.0 σ , 1.0 σ < r off /uni0084 1.5 σ , 1.5 σ < r off /uni0084 2.0 σ from the innermost circle to the outermost annulus, respectively. In both panels, coordinates are centered at the measured position for the cluster center of gravity ( RAJ = 19:10:52.04, DECJ = -59:59:04.64, Ferraro et al. 2003 ) . The reference color scale is displayed in the top horizontal band.  \n<!-- image -->  \nrelevant contributions from nearer and farther points. Such a distribution may be a consequence of our division of the exploited area in annuli, whose areas are in the ratio 1:3:5:7. In order to  \nhave a better insight into the dependence of the amount of nonluminous mass on the position of the cluster center of gravity, we plotted in Figure 5 a color map of the obtained values  \nFigure 6. Combined unweighted ( red line, top panel ) and weighted ( green line, middle panel ) probability distributions for the nonluminous mass in the core of NGC 6752, determined as in Section 6.3. The probability distribution obtained using the coordinates of Ferraro et al. ( 2003 ) for the cluster center is also presented for immediate comparison ( blue line, bottom panel, label measured center ) . Units for the probability distributions are arbitrary.  \n<!-- image -->  \nfor M NL,mean as a function of the position of the center of gravity. In panel ( a ) the obtained color map is placed in the context of the cluster core, while panel ( b ) shows an insight into the region under consideration. The position of the measured coordinates lies close to the northern boundary of a clear large area where M NL,mean ( α , δ ) is up to ∼ 3500 M e . Its displacement in decl. and elongation in R.A. can be seen as related to the fact that this area covers decl. around the mean of the core pulsars. The highest values are mostly required for positions toward the northwest, where one /uniFB01 nds the farthest points from PSRD and PSRF, which are the pulsars with the highest  P /uni2223 /uni2223 ( PSRD ) and distance ( PSRF ) from the measured coordinates.  \nA /uniFB01 nal insight can be obtained after calculating the overall probability distribution P ( M NL ) , by marginalizing P ( M NL, α GC, δ GC ) over the position of the cluster center. We performed this sum by taking into account the probability that the true center of gravity of the cluster is located in the considered position ( weighted case ) :  \n a d a d d a d = \'  23 P M P M P P d d , , cos r W NL 1 NL GC GC GC GC GC GC GC off /uni222C ( ) ( ) ( ) ( ) ( )  \nwhere P ( α GC ) and P ( δ GC ) are Gaussian distributions with mean value α GC and δ GC, respectively, and standard deviation σ = 0 5 in both cases, and also in the case of a /uniFB02 at distribution for the probability density of both coordinates ( unweighted case ) :  \n a d d a d =  P M P M d d , , cos . r U NL NL GC GC GC GC GC /uni222C ( ) ( )  \n24 1 off ( )  \nThe top and middle panels of Figure 6 display the obtained unweighted and weighted distributions, jointly with the distribution for the nonluminous mass we obtained in Section 4, i.e., in the case where the cluster center of gravity is located at the measured coordinates. Both P W and P U peak at the same value of 2.97 × 10 3 M e , and have similar lower uncertainties at 1 σ ( -0.23 × 10 3 M e and -0.21 × 10 3 M e ) , at 2 σ ( -0.33 × 10 3 M e and -0.29 × 10 3 M e , respectively ) and 3 σ ( -0.41 × 10 3 M e and -0.39 × 10 3 M e , respectively ) . The higher mass tails are clearly not negligible up to 8 × 10 3 M e but, as shown by the histograms in Figure 4 and panel ( b ) in Figure 5,an  \namount of nonluminous mass of at least 5 × 10 3 M e is required only if the true gravity center of the cluster has an offset r off /uni0085 0 5 from the measured coordinates, i.e., /uni0085 1 σ away from the position measured by Ferraro et al. ( 2003 ) . If one extracts from the probability distribution P W ( M NL ) a measure for M NL with its ( asymmetric ) uncertainties, one obtains = M NL  \' -+ M 2.97 10 0.23 0.88 3 at the 1 σ level, with a 3 σ lower limit of 2.56 × 10 3 M e . This is the most constrained value for M NL we can obtain from our data after considering all the sources of uncertainty, yet its errors are much larger than the one introduced by our poor knowledge of the intrinsic /uniFB01 rst and second-order time derivative of the spin period of the core pulsars. A decrease by a factor of 10 on the uncertainties on the gravity center coordinates is still not enough for the uncertainty on M NL to be equal to the one induced by the the poor knowledge of  P . Adopting a very conservative approach, we do not conclude with a de /uniFB01 nite value for M NL, since from both P W and P U we deduce a 3 σ upper limit on M NL at least of /uni0085 10 4 M e . At the same time, we can place the solid lower limit of M NL /uni0085 2.56 × 10 3 M e on the amount of nonluminous mass in the core of NGC 6752.  \nThe position of the cluster \' s center of gravity is thus the dominant source of uncertainty on the estimation of the nonluminous mass in the core of NGC 6752. If its coordinates were known with a precision comparable to the one of the position of the pulsars in the core, the poor knowledge of both  P and P /uni0308 for the core pulsars would become the main sources of uncertainty. But in such a situation, the resulting uncertainty on M NL would be of only a few percent in a case like NGC 6752, thus providing a good test bed for modeling and simulations.', '7. Summary': "We presented a timing analysis, spanning about 21 yr of observations, of the isolated pulsars PSR J1910 -5959B, PSR J1910-5959C, PSR J1910 -5959D, PSR J1910 -5959E, and PSR J1910 -5959F in the GC NGC 6752. The measured spin period derivative for the pulsars in the core is dominated by the gravitational pull of the mass in the central regions, and we estimated the amount of mass that is required to fully justify these measurements. We found that at least 2.56 × 10 3 M e , not yet taken into account in models and / or simulations, are present under the form of nonluminous matter, within a sphere of radius not larger than the distance of the closest pulsar to the cluster center of gravity. We also explored the scenario where an IMBH resides at the center of gravity of this cluster. We found that the presence of an IMBH with a mass of at least a few thousand solar masses can explain the observed accelerations and jerks experienced by the core pulsars, but is highly, although not completely, incompatible with the observational constraints deductible by the works of other authors, and that all these constraints hardly allow this scenario to provide all the necessary extra mass. We thus inferred that this extra amount of matter is most likely in the form of a population of much lighter sub- or nonluminous objects, whose exact nature and distribution cannot be investigated with the present data. Once this extra mass is taken into account, the inferred mass-to-light ratio in the central regions is higher than that jointly deduced from optical observations and from simulations of the structure and evolution of this GC. We have also shown that the main source of uncertainty in our results is the error in the coordinates of the cluster ' s center of gravity, and that our poor knowledge of the intrinsic spindown  \nof the core pulsars has a smaller impact, unless the precision on the coordinates of the cluster center were comparable to that of the core pulsars, as one can obtain from their timing.", 'Acknowledgments': 'The MeerKAT telescope is operated by the South African Radio Astronomy Observatory, which is a facility of the National Research Foundation, an agency of the Department of Science and Innovation. SARAO acknowledges the ongoing advice and calibration of GPS systems by the National Metrology Institute of South Africa ( NMISA ) and the time -space reference systems department of the Paris Observatory. MeerTime data is housed on the OzSTAR supercomputer at Swinburne University of Technology supported by ADACS and the Gravitational Wave Data Centre via Astronomy Australia Ltd. MeerKAT observations used the PTUSE and APSUSE computing clusters for data acquisition, storage, and analysis. These clusters were funded and installed by the MaxPlanck-Institut für Radioastronomie ( MPIfR ) and the MaxPlanck-Gesellschaft. The Parkes radio telescope is funded by the Commonwealth of Australia for operation as a National Facility managed by CSIRO. We acknowledge the Wiradjuri people as the traditional owners of the Observatory site. M.Ba. acknowledges support from the ARC Centre of Excellence for Gravitational Wave Discovery ( OzGrav ) under grant CE170100004. M.Ba. acknowledges support from the Australian Research Council via grant CE170100004, the ARC Centre of Excellence for Gravitational Wave Discovery ( OzGrav ) . Data used in this project made use of the OzSTAR supercomputer supported by ADACS and the GWDC. Part of this work has been funded using resources from the INAF Large Grant 2022 GCjewels ( P.I. Andrea Possenti ) approved with the Presidential Decree 30 / 2022. This work was also in part supported by the Italian Ministry of Foreign Affairs and International Cooperation, grant No. ZA23GR03.', 'ORCID iDs': 'A. Ridol /uniFB01 https: // orcid.org / 0000-0001-6762-2638 F. Abbate https: // orcid.org / 0000-0002-9791-7661 M. Bailes https: // orcid.org / 0000-0003-3294-3081 A. Possenti https: // orcid.org / 0000-0001-5902-3731 R. N. Manchester https: // orcid.org / 0000-0001-9445-5732 M. Kramer https: // orcid.org / 0000-0002-4175-2271 P. C. C. Freire https: // orcid.org / 0000-0003-1307-9435 M. Burgay https: // orcid.org / 0000-0002-8265-4344 F. Camilo https: // orcid.org / 0000-0002-1873-3718', 'References': "Abbate, F., Ridol /uniFB01 , A., Barr, E. D., et al. 2022, MNRAS, 513, 2292  \n- Abbate, F., Spera, M., & Colpi, M. 2019, MNRAS, 487, 769\n- Abdollahi, S., Acero, F., Baldini, L., et al. 2022, ApJS, 260, 53\n- Anderson, S. B. 1993, PhD thesis, California Institute of Technology\n- Bahcall, J. N., & Wolf, R. A. 1976, ApJ, 209, 214\n- Bahcall, J. N., & Wolf, R. A. 1977, ApJ, 216, 883\n- Bailes, M., Jameson, A., Abbate, F., et al. 2020, PASA, 37, e028\n- Barr, E. D. 2018, in IAU Symp. 337, Pulsar Astrophysics the Next Fifty\n- Years, ed. P. Weltevrede et al., 175\n- Baumgardt, H. 2017, MNRAS, 464, 2174\n- Baumgardt, H., & Hilker, M. 2018, MNRAS, 478, 1520\n- Baumgardt, H., Makino, J., & Ebisuzaki, T. 2004a, ApJ, 613, 1133\n- Baumgardt, H., Makino, J., & Ebisuzaki, T. 2004b, ApJ, 613, 1143\n- Baumgardt, H., & Vasiliev, E. 2021, MNRAS, 505, 5957\n- Camilo, F., Scholz, P., Serylak, M., et al. 2018, ApJ, 856, 180  \n- Chen, W., Barr, E., Karuppusamy, R., Kramer, M., & Stappers, B. 2021, JAI, 10, 2150013\n- Colpi, M., Possenti, A., & Gualandris, A. 2002, ApJL, 570, L85\n- Corongiu, A., Possenti, A., Lyne, A. G., et al. 2006, ApJ, 653, 1417\n- Corongiu, A., Venkatraman Krishnan, V., Freire, P. C. C., et al. 2023, A&A, 671, A72\n- D ' Amico, N., Lyne, A. G., Manchester, R. N., Possenti, A., & Camilo, F. 2001, ApJL, 548, L171\n- D ' Amico, N., Possenti, A., Fici, L., et al. 2002, ApJL, 570, L89\n- Desvignes, G., Caballero, R. N., Lentati, L., et al. 2016, MNRAS, 458, 3341 Ferraro, F. R., Possenti, A., Sabbi, E., et al. 2003, ApJ, 595, 179\n- Gaia Collaboration, Prusti, T., de Bruijne, J. H. J., et al. 2016, A&A, 595, A1 Gravity Collaboration, Abuter, R., Amorim, A., et al. 2021, A&A, 647, A59\n- Hilker, M., Baumgardt, H., Sollima, A., & Bellini, A. 2020, in IAU Symp. 351, Star Clusters: From the Milky Way to the Early Universe, ed. A. Bragaglia et al., 451\n- Hobbs, G. B., Edwards, R. T., & Manchester, R. N. 2006, MNRAS, 369, 655 Hotan, A. W., van Straten, W., & Manchester, R. N. 2004, PASA, 21, 302\n- Jonas, J. & MeerKAT Team 2016, in Proc. of MeerKAT Science: On the Pathway to the SKA ( Trieste: Proceedings of Science ) , 1\n- King, I. 1966, AJ, 71, 64\n- Lazaridis, K., Wex, N., Jessner, A., et al. 2009, MNRAS, 400, 805\n- Leigh, N. W. C., Ye, C. S., Grondin, S. M., et al. 2024, MNRAS, 527, 6913  \nLeitinger, E., Baumgardt, H., Cabrera-Ziri, I., Hilker, M., & Pancino, E. 2023,  \n- MNRAS, 520, 1456  \n- Li, H., & Widrow, L. M. 2021, MNRAS, 503, 1586 Libralato, M., Bellini, A., Vesperini, E., et al. 2022, ApJ, 934, 150 Manchester, R. N., Hobbs, G., Bailes, M., et al. 2013, PASA, 30, e017 Manchester, R. N., Hobbs, G. B., Teoh, A., & Hobbs, M. 2005, AJ, 129, 1993\n- Noyola, E., & Gebhardt, K. 2006, AJ, 132, 447\n- Pan, Z., Ma, X.-Y., Qian, L., et al. 2021, RAA, 21, 143\n- Phinney, E. S. 1993, in ASP Conf. Ser. 50, Structure and Dynamics of Globular Clusters, ed. S. G. Djorgovski & G. Meylan ( San Francisco, CA: ASP ) , 141\n- Prager, B. J., Ransom, S. M., Freire, P. C. C., et al. 2017, ApJ, 845, 148\n- Pryor, C., & Meylan, G. 1993, in ASP Conf. Ser. 50, Structure and Dynamics\n- of Globular Clusters, ed. S. G. Djorgovski & G. Meylan ( San Francisco, CA: ASP ) , 357\n- Reardon, D. J., Shannon, R. M., Cameron, A. D., et al. 2021, MNRAS, 507, 2137\n- Ridol /uniFB01 , A., Gautam, T., Freire, P. C. C., et al. 2021, MNRAS, 504, 1407\n- Scholz, P., Kaspi, V. M., Lyne, A. G., et al. 2015, ApJ, 800, 123\n- Shklovskii, I. S. 1970, SvA, 13, 562\n- Stappers, B., & Kramer, M. 2016, in Proc. of MeerKAT Science: On the Pathway to the SKA ( Trieste: Proceedings of Science ) , 9\n- van Straten, W., & Bailes, M. 2011, PASA, 28, 1\n- Vasiliev, E., & Baumgardt, H. 2021, MNRAS, 505, 5978\n- Verbunt, F., & Freire, P. C. C. 2014, A&A, 561, A11\n- Vitral, E., Libralato, M., Kremer, K., et al. 2023, MNRAS, 522, 5740"}