Star count density profiles and structural parameters of 26 Galactic globular clusters

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.Comment: Accepted for publication in The Astrophysical Journal; 19 pages (emulateapj style), 15 figures, 2 table

The binary fraction in the globular cluster M10 (NGC 6254): comparing core and outer regions

We study the binary fraction of the globular cluster M10 (NGC 6254) as a function of radius from the cluster core to the outskirts, by means of a quan- titative analysis of the color distribution of stars relative to the fiducial main sequence. By taking advantage of two data-sets, acquired with the Advanced Camera for Survey and the Wide Field Planetary Camera 2 on board the Hubble Space Telescope, we have studied both the core and the external regions of the cluster. The binary fraction is found to decrease from 14% within the core, to 1.5% in a region between 1 and 2 half-mass radii from the cluster centre. Such a trend and the derived values are in agreement with previous results ob- tained in clusters of comparable total magnitude. The estimated binary fraction is sufficient to account for the suppression of mass segregation observed in M10, without any need to invoke the presence of an intermediate-mass black hole in its centre.Comment: Accepted for publication in ApJ (22 pages, 7 figures, 3 tables

Fundamental parameters of 16 late-type stars derived from their angular diameter measured with VLTI/AMBER

Thanks to their large angular dimension and brightness, red giants and supergiants are privileged targets for optical long-baseline interferometers. Sixteen red giants and supergiants have been observed with the VLTI/AMBER facility over a two-years period, at medium spectral resolution (R=1500) in the K band. The limb-darkened angular diameters are derived from fits of stellar atmospheric models on the visibility and the triple product data. The angular diameters do not show any significant temporal variation, except for one target: TX Psc, which shows a variation of 4% using visibility data. For the eight targets previously measured by Long-Baseline Interferometry (LBI) in the same spectral range, the difference between our diameters and the literature values is less than 5%, except for TX Psc, which shows a difference of 11%. For the 8 other targets, the present angular diameters are the first measured from LBI. Angular diameters are then used to determine several fundamental stellar parameters, and to locate these targets in the Hertzsprung-Russell Diagram (HRD). Except for the enigmatic Tc-poor low-mass carbon star W Ori, the location of Tc-rich stars in the HRD matches remarkably well the thermally-pulsating AGB, as it is predicted by the stellar-evolution models. For pulsating stars with periods available, we compute the pulsation constant and locate the stars along the various sequences in the Period -- Luminosity diagram. We confirm the increase in mass along the pulsation sequences, as predicted by the theory, except for W Ori which, despite being less massive, appears to have a longer period than T Cet along the first-overtone sequence.Comment: 15 pages, 9 figures, 6 table

The extended halo of NGC 2682 (M 67) from Gaia DR2

Context: NGC 2682 is a nearby open cluster, approximately 3.5 Gyr old. Dynamically, most open clusters should dissolve on shorter timescales, of ~ 1 Gyr. Having survived until now, NGC 2682 was likely much more massive in the past, and is bound to have an interesting dynamical history. Aims: We investigate the spatial distribution of NGC 2682 stars to constrain its dynamical evolution, especially focusing on the marginally bound stars in the cluster outskirts. Methods: We use Gaia DR2 data to identify NGC 2682 members up to a distance of ~150 pc (10 degrees). Two methods (Clusterix and UPMASK) are applied to this end. We estimate distances to obtain three-dimensional stellar positions using a Bayesian approach to parallax inversion, with an appropriate prior for star clusters. We calculate the orbit of NGC 2682 using the GRAVPOT16 software. Results: The cluster extends up to 200 arcmin (50 pc) which implies that its size is at least twice as previously believed. This exceeds the cluster Hill sphere based on the Galactic potential at the distance of NGC 2682. Conclusions: The extra-tidal stars in NGC 2682 may originate from external perturbations such as disk shocking or dynamical evaporation from two-body relaxation. The former origin is plausible given the orbit of NGC 2682, which crossed the Galactic disk ~40 Myr ago.Comment: 9 pages, 5 figures, accepted for publication on A&

Probing the role of dynamical friction in shaping the BSS radial distribution. I - Semi-analytical models and preliminary N-body simulations

We present semi-analytical models and simplified $N$-body simulations with $10^4$ and $10^5$ particles aimed at probing the role of dynamical friction (DF) in determining the radial distribution of Blue Straggler Stars (BSSs) in globular clusters. The semi-analytical models show that DF (which is the only evolutionary mechanism at work) is responsible for the formation of a bimodal distribution with a dip progressively moving toward the external regions of the cluster. However, these models fail to reproduce the formation of the long-lived central peak observed in all dynamically evolved clusters. The results of $N$-body simulations confirm the formation of a sharp central peak, which remains as a stable feature over the time regardless of the initial concentration of the system. In spite of a noisy behavior, a bimodal distribution forms in many cases, with the size of the dip increasing as a function of time. In the most advanced stages the distribution becomes monotonic. These results are in agreement with the observations. Also the shape of the peak and the location of the minimum (which in most of the cases is within 10 core radii) turn out to be consistent with observational results. For a more detailed and close comparison with observations, including a proper calibration of the timescales of the dynamical processes driving the evolution of the BSS spatial distribution, more realistic simulations will be necessary.Comment: Accepted for publication by ApJ; 11 pages, 11 figure

The dynamical state of the Globular Cluster M10 (NGC 6254)

Studying the radial variation of the stellar mass function in globular clusters (GCs) has proved a valuable tool to explore the collisional dynamics leading to mass segregation and core collapse. In order to study the radial dependence of the luminosity and mass function of M 10, we used ACS/HST deep high resolution archival images, reaching out to approximately the cluster's half-mass radius (rhm), combined with deep WFPC2 images that extend our radial coverage to more than 2 rhm. From our photometry, we derived a radial mass segregation profile and a global mass function that we compared with those of simulated clusters containing different energy sources (namely hard binaries and/or an IMBH) able to halt core collapse and to quench mass segregation. A set of direct N-body simulations of GCs, with and without an IMBH of mass 1% of the total cluster mass, comprising different initial mass functions (IMFs) and primordial binary fractions, was used to predict the observed mass segregation profile and mass function. The mass segregation profile of M 10 is not compatible with cluster models without either an IMBH or primordial binaries, as a source of energy appears to be moderately quenching mass segregation in the cluster. Unfortunately, the present observational uncertainty on the binary fraction in M10 does not allow us to confirm the presence of an IMBH in the cluster, since an IMBH, a dynamically non-negligible binary fraction (~ 5%), or both can equally well explain the radial dependence of the cluster mass function.Comment: 15 pages, 8 figures, accepted for publication on Ap

A Molecular Sensor To Characterize Arenavirus Envelope Glycoprotein Cleavage by Subtilisin Kexin Isozyme 1/Site 1 Protease.

UNLABELLED: Arenaviruses are emerging viruses including several causative agents of severe hemorrhagic fevers in humans. The advent of next-generation sequencing technology has greatly accelerated the discovery of novel arenavirus species. However, for many of these viruses, only genetic information is available, and their zoonotic disease potential remains unknown. During the arenavirus life cycle, processing of the viral envelope glycoprotein precursor (GPC) by the cellular subtilisin kexin isozyme 1 (SKI-1)/site 1 protease (S1P) is crucial for productive infection. The ability of newly emerging arenaviruses to hijack human SKI-1/S1P appears, therefore, to be a requirement for efficient zoonotic transmission and human disease potential. Here we implement a newly developed cell-based molecular sensor for SKI-1/S1P to characterize the processing of arenavirus GPC-derived target sequences by human SKI-1/S1P in a quantitative manner. We show that only nine amino acids flanking the putative cleavage site are necessary and sufficient to accurately recapitulate the efficiency and subcellular location of arenavirus GPC processing. In a proof of concept, our sensor correctly predicts efficient processing of the GPC of the newly emergent pathogenic Lujo virus by human SKI-1/S1P and defines the exact cleavage site. Lastly, we employed our sensor to show efficient GPC processing of a panel of pathogenic and nonpathogenic New World arenaviruses, suggesting that GPC cleavage represents no barrier for zoonotic transmission of these pathogens. Our SKI-1/S1P sensor thus represents a rapid and robust test system for assessment of the processing of putative cleavage sites derived from the GPCs of newly discovered arenavirus by the SKI-1/S1P of humans or any other species, based solely on sequence information. IMPORTANCE: Arenaviruses are important emerging human pathogens that can cause severe hemorrhagic fevers with high mortality in humans. A crucial step in productive arenavirus infection of human cells is the processing of the viral envelope glycoprotein by the cellular subtilisin kexin isozyme 1 (SKI-1)/site 1 protease (S1P). In order to break the species barrier during zoonotic transmission and cause severe disease in humans, newly emerging arenaviruses must be able to hijack human SKI-1/S1P efficiently. Here we implement a newly developed cell-based molecular sensor for human SKI-1/S1P to characterize the processing of arenavirus glycoproteins in a quantitative manner. We further use our sensor to correctly predict efficient processing of the glycoprotein of the newly emergent pathogenic Lujo virus by human SKI-1/S1P. Our sensor thus represents a rapid and robust test system with which to assess whether the glycoprotein of any newly emerging arenavirus can be efficiently processed by human SKI-1/S1P, based solely on sequence information

Interpreting automatic AGN classifiers with saliency maps

Classification of the optical spectra of active galactic nuclei (AGN) into different types is currently based on features such as line widths and intensity ratios. Although well founded on AGN physics, this approach involves some degree of human oversight and cannot scale to large datasets. Machine learning (ML) tackles this classification problem in a fast and reproducible way, but is often (and not without reason) perceived as a black box. However, ML interpretability and are active research areas in computer science that are providing us with tools to mitigate this issue. We apply ML interpretability tools to a classifier trained to predict AGN types from spectra. Our goal is to demonstrate the use of such tools in this context, obtaining for the first time insight into an otherwise black box AGN classifier. In particular, we want to understand which parts of each spectrum most affect the predictions of our classifier, checking that the results make sense in the light of our theoretical expectations. We trained a support-vector machine on 3346 high-quality, low-redshift AGN spectra from SDSS DR15. We considered either two-class classification (type 1 versus 2) or multiclass (type 1 versus 2 versus intermediate-type). The spectra were previously and independently hand-labeled and divided into types 1 and 2, and intermediate-type (i.e., sources in which the Balmer line profile consists of a sharp narrow component superimposed on a broad component). We performed a train-validation-test split, tuning hyperparameters and independently measuring performance via a variety of metrics. On a selection of test-set spectra, we computed the gradient of the predicted class probability at a given spectrum. Regions of the spectrum were then color-coded based on the direction and the amount by which they influence the predicted class, effectively building a saliency map. We also visualized the high-dimensional space of AGN spectra using t-distributed stochastic neighbor embedding (t-SNE), showing where the spectra for which we computed a saliency map are located. Our best classifier reaches an F-score of 0.942 on our test set (with 0.948 precision and 0.936 recall). We computed saliency maps on all misclassified spectra in the test set and on a sample of randomly selected spectra. Regions that affect the predicted AGN type often coincide with physically relevant features, such as spectral lines. t-SNE visualization shows good separability of type 1 and type 2 spectra. Intermediate-type spectra either lie in-between, as expected, or appear mixed with type 2 spectra. Misclassified spectra are typically found among the latter. Some clustering structure is apparent among type 2 and intermediate-type spectra, though this may be an artifact. Saliency maps show why a given AGN type was predicted by our classifier resulting in a physical interpretation in terms of regions of the spectrum that affected its decision, making it no longer a black box. These regions coincide with those used by human experts, for example relevant spectral lines, and are even used in a similar way; the classifier effectively measures the width of a line by weighing its center and its tails oppositely

Dynamical age differences among coeval star clusters as revealed by blue stragglers

Globular star clusters that formed at the same cosmic time may have evolved rather differently from a dynamical point of view (because that evolution depends on the internal environment) through a variety of processes that tend progressively to segregate stars more massive than the average towards the cluster centre. Therefore clusters with the same chronological age may have reached quite different stages of their dynamical history (that is, they may have different dynamical ages). Blue straggler stars have masses greater than those at the turn-off point on the main sequence and therefore must be the result of either a collision or a mass-transfer event. Because they are among the most massive and luminous objects in old clusters, they can be used as test particles with which to probe dynamical evolution. Here we report that globular clusters can be grouped into a few distinct families on the basis of the radial distribution of blue stragglers. This grouping corresponds well to an effective ranking of the dynamical stage reached by stellar systems, thereby permitting a direct measure of the cluster dynamical age purely from observed properties.Comment: Published on the 20 December 2012 issue of Natur

Constructing Impactful Machine Learning Research for Astronomy: Best Practices for Researchers and Reviewers

Machine learning has rapidly become a tool of choice for the astronomical community. It is being applied across a wide range of wavelengths and problems, from the classification of transients to neural network emulators of cosmological simulations, and is shifting paradigms about how we generate and report scientific results. At the same time, this class of method comes with its own set of best practices, challenges, and drawbacks, which, at present, are often reported on incompletely in the astrophysical literature. With this paper, we aim to provide a primer to the astronomical community, including authors, reviewers, and editors, on how to implement machine learning models and report their results in a way that ensures the accuracy of the results, reproducibility of the findings, and usefulness of the method.Comment: 14 pages, 3 figures; submitted to the Bulletin of the American Astronomical Societ