Search

As Large Language Models (LLMs) are integrated into electronic health record (EHR) workflows, validated instruments are essential to evaluate their performance before implementation. Existing instruments for provider documentation quality are often unsuitable for the complexities of LLM-generated text and lack validation on real-world data. The Provider Documentation Summarization Quality Instrument (PDSQI-9) was developed to evaluate LLM-generated clinical summaries. Multi-document summaries were generated from real-world EHR data across multiple specialties using several LLMs (GPT-4o, Mixtral 8x7b, and Llama 3-8b). Validation included Pearson correlation for substantive validity, factor analysis and Cronbach's alpha for structural validity, inter-rater reliability (ICC and Krippendorff's alpha) for generalizability, a semi-Delphi process for content validity, and comparisons of high-versus low-quality summaries for discriminant validity. Seven physician raters evaluated 779 summaries and answered 8,329 questions, achieving over 80% power for inter-rater reliability. The PDSQI-9 demonstrated strong internal consistency (Cronbach's alpha = 0.879; 95% CI: 0.867-0.891) and high inter-rater reliability (ICC = 0.867; 95% CI: 0.867-0.868), supporting structural validity and generalizability. Factor analysis identified a 4-factor model explaining 58% of the variance, representing organization, clarity, accuracy, and utility. Substantive validity was supported by correlations between note length and scores for Succinct (rho = -0.200, p = 0.029) and Organized ($\rho = -0.190$, $p = 0.037$). Discriminant validity distinguished high- from low-quality summaries ($p < 0.001$). The PDSQI-9 demonstrates robust construct validity, supporting its use in clinical practice to evaluate LLM-generated summaries and facilitate safer integration of LLMs into healthcare workflows.

For decades, it has been theorized that a tenuous but detectable intracluster medium should be present in globular clusters, which is continuously replenished by the gas and dust ejected by bright giants and periodically cleared by interactions with the Galactic disk. However, dedicated searches, especially in infrared and radio wavelengths, have returned mostly upper limits, which are lower than theoretical expectations by several orders of magnitude. We profited from recent wide-field photometry for 48 Galactic globular clusters to compute high-resolution maps of differential reddening, which can be used to correct any photometric catalog in these areas for reddening variations. Using 3D reddening maps from the literature, we evaluated the amount of foreground extinction. This allowed us to estimate the masses of the intracluster medium in our sample clusters, with an accuracy of one order of magnitude. Our estimates agree with the few available literature detections and with theoretical expectations. Because the discrepancy between observations and expectations only concerns literature upper limits, we explored possible reasons why they could be underestimated and we show that two recent discoveries can explain the discrepancy. The first is the recent discovery that the intracluster medium in 47 Tuc is not centrally concentrated. This is also supported by our maps, which in the majority of cases do not show a central reddening concentration. The second is the discovery that the dust in metal-poor ([Fe/H] less than about -1 dex) globular clusters is dominated by iron grains rather than silicates, which undermines previous dust mass estimates from observed upper limits. We conclude that current evidence, including our maps, does not contradict theoretical expectations and the problem of the missing intracluster medium is no longer an issue., Comment: Accepted by A&A, 20 pages, 14 figures

Circulating tumour DNA (ctDNA) detection of molecular residual disease (MRD) in solid tumours correlates strongly with patient outcomes and is being adopted as a new clinical standard. ctDNA levels are known to correlate with tumor volume, and although the absolute levels vary across indication and histology, its analysis is driving the adoption of MRD. MRD assays must detect tumor when imaging cannot and, as such, require very high sensitivity to detect the low levels of ctDNA found after curative intent therapy. The minimum threshold is 0.01% Tumour Fraction but current methods like Archer and Signatera are limited by detection sensitivity resulting in some patients receiving a false negative call thereby missing out on earlier therapeutic intervention. Multiple vendors are increasing the number of somatic variants tracked in tumour-informed and personalized NGS assays, from tens to thousands of variants. Most recently, assays using other biological features of ctDNA, e.g methylation or fragmentome, have been developed at the LOD required for clinical utility. These uniformed, or tumour-naive and non-personalised assays may be more easily, and therefore more rapidly, adopted in the clinic. However, this rapid development in MRD assay technology results in significant challenges in benchmarking these new technologies for use in clinical trials. This is further complicated by the fact that previous reference materials have focused on somatic variants, and do not retain all of the epigenomic features assessed by newer technologies. In this Comments and Controversy paper, we detail what is known and what remains to be determined for optimal reference materials of MRD methods and provide opinions generated during three-years of MRD technology benchmarking in AstraZeneca Translational Medicine to help guide the community conversation.

We present a photometric analysis of globular cluster 47 Tuc (NGC\,104), using near-IR imaging data from the GeMS/GSAOI Galactic Globular Cluster Survey (G4CS) which is in operation at Gemini-South telescope.~Our survey is designed to obtain AO-assisted deep imaging with near diffraction-limited spatial resolution of the central fields of Milky Way globular clusters.~The G4CS near-IR photometry was combined with an optical photometry catalog obtained from Hubble Space Telescope survey data to produce a high-quality color-magnitude diagram that reaches down to K$_s\approx$ 21 Vega mag.~We used the software suite BASE-9, which uses an adaptive Metropolis sampling algorithm to perform a Markov chain Monte Carlo (MCMC) Bayesian analysis, and obtained probability distributions and precise estimates for the age, distance and extinction cluster parameters.~Our best estimate for the age of 47 Tuc is 12.42$^{+0.05}_{-0.05}$ $\pm$ 0.08 Gyr, and our true distance modulus estimate is (m$-$M)$_0$=13.250$^{+0.003}_{-0.003}$ $\pm$ 0.028 mag, in tight agreement with previous studies using Gaia DR2 parallax and detached eclipsing binaries., Comment: 25 pages, 15 figures. Accepted for publication in The Astrophysical Journal

We investigate the properties of the mixed-mode (RRd) RR Lyrae (RRL) variables in the Fornax dwarf spheroidal (dSph) galaxy by using $B$- and $V$-band time series collected over twenty-four years. We compare the properties of the RRds in Fornax with those in the Magellanic Clouds and in nearby dSphs, with special focus on Sculptor. We found that the ratio of RRds over the total number of RRLs decreases with metallicity. Typically, dSphs have very few RRds with 0.49$\ltsim P_0 \ltsim $0.53 days, but Fornax fills this period gap in the Petersen diagram (ratio between first overtone over fundamental period versus fundamental period). We also found that the distribution in the Petersen diagram of Fornax RRds is similar to SMC RRds, thus suggesting that their old stars have a similar metallicity distribution. We introduce the Period-Amplitude RatioS (PARS) diagram, a new pulsation diagnostics independent of distance and reddening. We found that LMC RRds in this plane are distributed along a short- and a long-period sequence that we identified as the metal-rich and the metal-poor component. These two groups are also clearly separated in the Petersen and Bailey (luminosity amplitude versus logarithmic period) diagrams. These circumstantial evidence indicates that the two groups have different evolutionary properties. All the pulsation diagnostics adopted in this investigation suggest that old stellar populations in Fornax and Sculptor dSphs underwent different chemical enrichment histories. Fornax RRds are similar to SMC RRds, while Sculptor RRds are more similar to the metal-rich component of the LMC RRds., Comment: Accepted for publication by MNRAS

Internet analysis is a major challenge due to the volume and rate of network traffic. In lieu of analyzing traffic as raw packets, network analysts often rely on compressed network flows (netflows) that contain the start time, stop time, source, destination, and number of packets in each direction. However, many traffic analyses benefit from temporal aggregation of multiple simultaneous netflows, which can be computationally challenging. To alleviate this concern, a novel netflow compression and resampling method has been developed leveraging GraphBLAS hyperspace traffic matrices that preserve anonymization while enabling subrange analysis. Standard multitemporal spatial analyses are then performed on each subrange to generate detailed statistical aggregates of the source packets, source fan-out, unique links, destination fan-in, and destination packets of each subrange which can then be used for background modeling and anomaly detection. A simple file format based on GraphBLAS sparse matrices is developed for storing these statistical aggregates. This method is scale tested on the MIT SuperCloud using a 50 trillion packet netflow corpus from several hundred sites collected over several months. The resulting compression achieved is significant (<0.1 bit per packet) enabling extremely large netflow analyses to be stored and transported. The single node parallel performance is analyzed in terms of both processors and threads showing that a single node can perform hundreds of simultaneous analyses at over a million packets/sec (roughly equivalent to a 10 Gigabit link)., Comment: arXiv admin note: text overlap with arXiv:2203.13934, arXiv:2108.06653, arXiv:2008.00307

We present accurate and deep multi-band ($g,r,i$) photometry of the Local Group dwarf irregular galaxy NGC 6822. The images were collected with wide field cameras at 2m/4m- (INT,CTIO,CFHT) and 8m-class telescopes (SUBARU) covering a 2 square degrees FoV across the center of the galaxy. We performed PSF photometry of $\approx$7,000 CCD images and the final catalog includes more than 1 million objects. We developed a new approach to identify candidate field and galaxy stars, and performed a new estimate of the galaxy center by using old stellar tracers finding that it differs by 1.15 (RA) and 1.53 (DEC) arcmin from previous estimates. We also found that young (Main Sequence, Red Supergiants), intermediate (Red Clump, Asymptotic Giant Branch [AGB]) and old (Red Giant Branch [RGB]) stars display different radial distributions. Old stellar population is spherically distributed and extends to radial distances larger than previously estimated ($\sim$1 degree). The young population shows a well defined bar and a disk-like distribution, as suggested by radio measurements, that is off-center compared with old population. We discuss pros and cons of the different diagnostics adopted to identify AGB stars and develop new ones based on optical-NIR-MIR color-color diagrams (CCDs) to characterize Oxygen and Carbon (C) rich stars. We found a mean population ratio between Carbon and M-type (C/M) stars of 0.67$\pm$0.08 (optical/NIR/MIR) and we used the observed C/M ratio with empirical C/M-metallicity relations to estimate a mean iron abundance of [Fe/H]$\sim$-1.25 ($\sigma$=0.04 dex) that agrees quite well with literature estimates., Comment: Accepted for publication in ApJ, 34 pages, 22 figures, 6 tables

We present the first results of eight Globular Clusters (GCs) from the AstroSat/UVIT Legacy Survey program GlobULeS based on the observations carried out in two FUV filters (F148W and F169M). The FUV-optical and FUV-FUV color-magnitude diagrams (CMDs) of GCs with the proper motion membership were constructed by combining the UVIT data with HST UV Globular Cluster Survey (HUGS) data for inner regions and Gaia Early Data Release (EDR3) for regions outside the HST's field. We detect sources as faint as F148W $\sim$ 23.5~mag which are classified based on their locations in CMDs by overlaying stellar evolutionary models. The CMDs of 8 GCs are combined with the previous UVIT studies of 3 GCs to create stacked FUV-optical CMDs to highlight the features/peculiarities found in the different evolutionary sequences. The FUV (F148W) detected stellar populations of 11 GCs comprises 2,816 Horizontal Branch (HB) stars (190 Extreme HB candidates), 46 post-HB (pHB), 221 Blue Straggler Stars (BSS), and 107 White Dwarf (WD) candidates. We note that the blue HB color extension obtained from F148W$-$G color and the number of FUV detected EHB candidates are strongly correlated with the maximum internal Helium (He) variation within each GC, suggesting that the FUV-optical plane is the most sensitive to He abundance variations in the HB. We discuss the potential science cases that will be addressed using these catalogues including HB morphologies, BSSs, pHB, and, WD stars., Comment: Accepted for publication in MNRAS; 24 pages, 20 Figures, 6 Tables

We present new period-$\phi_{31}$-[Fe/H] relations for first overtone RRL stars (RRc), calibrated over a broad range of metallicities ($-2.5 < \textrm{[Fe/H]}< 0.0$) utilizing the largest currently available set of Galactic halo field RRL with homogeneous spectroscopic metallicities. Our relations are defined in the optical (ASAS-SN $V$-band) and, inaugurally, in the infrared (WISE $W1$ and $W2$ bands). Our $V$-band relation can reproduce individual RRc spectroscopic metallicities with a dispersion of 0.30 dex over the entire metallicity range of our calibrator sample (an RMS smaller than what we found for other relations in literature including non-linear terms). Our infrared relation has a similar dispersion in the low and intermediate metallicity range ($\textrm{[Fe/H]} < -0.5$) but tends to underestimate the [Fe/H] abundance around solar metallicity. We tested our relations by measuring both the metallicity of the Sculptor dSph and a sample of Galactic globular clusters, rich in both RRc and RRab stars. The average metallicity we obtain for the combined RRL sample in each cluster is within $\pm 0.08$ dex of their spectroscopic metallicities. The infrared and optical relations presented in this work will enable deriving reliable photometric RRL metallicities in conditions where spectroscopic measurements are not feasible; e.g., in distant galaxies or reddened regions (observed with upcoming Extremely Large Telescopes and the James Webb Space Telescope), or in the large sample of new RRL that will be discovered in large-area time-domain photometric surveys (such as LSST and the Roman space telescope)., Comment: Accepted by ApJ, 22 pages, 9 Figures, 3 Tables

Long range detection is a cornerstone of defense in many operating domains (land, sea, undersea, air, space, ..,). In the cyber domain, long range detection requires the analysis of significant network traffic from a variety of observatories and outposts. Construction of anonymized hypersparse traffic matrices on edge network devices can be a key enabler by providing significant data compression in a rapidly analyzable format that protects privacy. GraphBLAS is ideally suited for both constructing and analyzing anonymized hypersparse traffic matrices. The performance of GraphBLAS on an Accolade Technologies edge network device is demonstrated on a near worse case traffic scenario using a continuous stream of CAIDA Telescope darknet packets. The performance for varying numbers of traffic buffers, threads, and processor cores is explored. Anonymized hypersparse traffic matrices can be constructed at a rate of over 50,000,000 packets per second; exceeding a typical 400 Gigabit network link. This performance demonstrates that anonymized hypersparse traffic matrices are readily computable on edge network devices with minimal compute resources and can be a viable data product for such devices., Comment: Accepted to IEEE HPEC, Outstanding Paper Award, 8 pages, 8 figures, 1 table, 70 references. arXiv admin note: text overlap with arXiv:2108.06653, arXiv:2008.00307, arXiv:2203.10230

The Extremely Large Telescope and the Thirty Meter Telescope will use state of the art multiconjugate adaptive optics (MCAO) systems to obtain the full D4 advantage that their apertures can provide. However, to reach the full astrometric potential of these facilities for on-sky science requires understanding any residual astrometric distortions from these systems and find ways to measure and eliminate them. In this work, we use Gemini multiconjugate adaptive optic system (GeMS) observations of the core of NGC 6723 to better understand the on-sky astrometric performance of MCAO. We develop new methods to measure the astrometric distortion fields of the observing system, which probe the distortion at the highest possible spatial resolution. We also describe methods for examining the time-variable and static components of the astrometric distortion. When applied to the GeMS Gemini South Adaptive Optics Imager (GSAOI) data, we are able to see the effect of the field rotator at the subpixel level, and we are able to empirically derive the distortion due to the optical design of GeMS-GSAOI. We argue that the resulting distortion maps are a valuable tool to measure and monitor the on-sky astrometric performance of future instrumentation. Our overall astrometry pipeline produces high-quality proper motions with an uncertainty floor of 45 uas per year. We measure the proper motion dispersion profile of NGC 6723 from a radius of 10 arcsec out to 1 arcmin based on 12000 stars. We also produce a high-quality optical-near infrared color-magnitude diagram, which clearly shows the extreme horizontal branch and main-sequence knee of this cluster.

The Internet has become a critical component of modern civilization requiring scientific exploration akin to endeavors to understand the land, sea, air, and space environments. Understanding the baseline statistical distributions of traffic are essential to the scientific understanding of the Internet. Correlating data from different Internet observatories and outposts can be a useful tool for gaining insights into these distributions. This work compares observed sources from the largest Internet telescope (the CAIDA darknet telescope) with those from a commercial outpost (the GreyNoise honeyfarm). Neither of these locations actively emit Internet traffic and provide distinct observations of unsolicited Internet traffic (primarily botnets and scanners). Newly developed GraphBLAS hyperspace matrices and D4M associative array technologies enable the efficient analysis of these data on significant scales. The CAIDA sources are well approximated by a Zipf-Mandelbrot distribution. Over a 6-month period 70\% of the brightest (highest frequency) sources in the CAIDA telescope are consistently detected by coeval observations in the GreyNoise honeyfarm. This overlap drops as the sources dim (reduce frequency) and as the time difference between the observations grows. The probability of seeing a CAIDA source is proportional to the logarithm of the brightness. The temporal correlations are well described by a modified Cauchy distribution. These observations are consistent with a correlated high frequency beam of sources that drifts on a time scale of a month., Comment: 8 pages, 8 figures, 2 tables, 59 references; accepted to GrAPL 2022. arXiv admin note: substantial text overlap with arXiv:2108.06653

Adversarial Internet robots (botnets) represent a growing threat to the safe use and stability of the Internet. Botnets can play a role in launching adversary reconnaissance (scanning and phishing), influence operations (upvoting), and financing operations (ransomware, market manipulation, denial of service, spamming, and ad click fraud) while obfuscating tailored tactical operations. Reducing the presence of botnets on the Internet, with the aspirational target of zero, is a powerful vision for galvanizing policy action. Setting a global goal, encouraging international cooperation, creating incentives for improving networks, and supporting entities for botnet takedowns are among several policies that could advance this goal. These policies raise significant questions regarding proper authorities/access that cannot be answered in the abstract. Systems analysis has been widely used in other domains to achieve sufficient detail to enable these questions to be dealt with in concrete terms. Defeating botnets using an observe-pursue-counter architecture is analyzed, the technical feasibility is affirmed, and the authorities/access questions are significantly narrowed. Recommended next steps include: supporting the international botnet takedown community, expanding network observatories, enhancing the underlying network science at scale, conducting detailed systems analysis, and developing appropriate policy frameworks., Comment: 26 pages, 13 figures, 2 tables, 72 references, submitted to PlosOne

With the recognition of cyberspace as an operating domain, concerted effort is now being placed on addressing it in the whole-of-domain manner found in land, sea, undersea, air, and space domains. Among the first steps in this effort is applying the standard supporting concepts of security, defense, and deterrence to the cyber domain. This paper presents an architecture that helps realize forward defense in cyberspace, wherein adversarial actions are repulsed as close to the origin as possible. However, substantial work remains in making the architecture an operational reality including furthering fundamental research cyber science, conducting design trade-off analysis, and developing appropriate public policy frameworks.

We extend the approach Abbott, Kedlaya and Roe to computation of the zeta function of a projective hypersurface with $\tau$ isolated ordinary double points over a finite field $\mathbb{F}_q$ given by the reduction of a homogeneous polynomial $f \in \mathbb{Z}[x_0, \ldots, x_n]$, under the assumption of equisingularity over $\mathbb{Z}_q$. The algorithm is based on the results of Dimca and Saito (over the field $\mathbb{C}$ of complex numbers) on the pole order spectral sequence in the case of ordinary double points. We give some examples of explicit computations for surfaces in $\mathbb{P}^3$., Comment: Versions 2: some typos fixed and a reference [DiSt2] added

We present a detailed study of the variable star population of Eridanus II (Eri II), an ultra-faint dwarf galaxy that lies close to the Milky Way virial radius. We analyze multi-epoch $g,r,i$ ground-based data from Goodman and the Dark Energy Camera, plus $F475W, F606W, F814W$ space data from the Advanced Camera for Surveys. We report the detection of 67 RR Lyrae (RRL) stars and 2 Anomalous Cepheids, most of them new discoveries. With the RRL stars, we measure the distance modulus of Eri II, $\mu_0=22.84\pm 0.05$ mag (D$_{\odot}=370\pm9$ kpc) and derive a metallicity spread of 0.3 dex (0.2 dex intrinsic). The colour distribution of the horizontal branch (HB) and the period distribution of the RRL stars can be nicely reproduced by a combination of two stellar models of [Fe/H]=($-2.62$, $-2.14$). The overall low metallicity is consistent with the red giant branch bump location, 0.65 mag brighter than the HB. These results are in agreement with previous spectroscopic studies. The more metal-rich RRL and the RRab stars have greater central concentration than the more metal-poor RRL and the RRc stars that are mainly located outside $\sim 1$ r$_{\rm h}$. This is similar to what is found in larger dwarf galaxies such as Sculptor, and in agreement with an outside-in galaxy formation scenario. This is remarkable in such a faint dwarf galaxy with an apparently single and extremely short ($<1$ Gyr) star formation burst. Finally, we have derived new and independent structural parameters for Eri II and its star cluster using our new data that are in very good agreement with previous estimates., Comment: 20 pages, 14 figures, 8 tables. Accepted for publication in MNRAS

The Internet has never been more important to our society, and understanding the behavior of the Internet is essential. The Center for Applied Internet Data Analysis (CAIDA) Telescope observes a continuous stream of packets from an unsolicited darkspace representing 1/256 of the Internet. During 2019 and 2020 over 40,000,000,000,000 unique packets were collected representing the largest ever assembled public corpus of Internet traffic. Using the combined resources of the Supercomputing Centers at UC San Diego, Lawrence Berkeley National Laboratory, and MIT, the spatial temporal structure of anonymized source-destination pairs from the CAIDA Telescope data has been analyzed with GraphBLAS hierarchical hypersparse matrices. These analyses provide unique insight on this unsolicited Internet darkspace traffic with the discovery of many previously unseen scaling relations. The data show a significant sustained increase in unsolicited traffic corresponding to the start of the COVID19 pandemic, but relatively little change in the underlying scaling relations associated with unique sources, source fan-outs, unique links, destination fan-ins, and unique destinations. This work provides a demonstration of the practical feasibility and benefit of the safe collection and analysis of significant quantities of anonymized Internet traffic., Comment: 8 pages, 9 figures, 2 tables, 43 references, accepted to IEEE HPEC 2021. arXiv admin note: substantial text overlap with arXiv:2008.00307

Accurate astrometry is a key deliverable for the next generation of multi-conjugate adaptive optics (MCAO) systems. The MCAO Visible Imager and Spectrograph (MAVIS) is being designed for the Very Large Telescope Adaptive Optics Facility and must achieve 150 $\mu$as astrometric precision (50 $\mu$as goal). To test this before going on-sky, we have created MAVISIM, a tool to simulate MAVIS images. MAVISIM accounts for three major sources of astrometric error, high- and low-order point spread function (PSF) spatial variability, tip-tilt residual error and static field distortion. When exploring the impact of these three error terms alone, we recover an astrometric accuracy of 50 $\mu$as for all stars brighter than $m=19$ in a 30s integration using PSF-fitting photometry. We also assess the feasibility of MAVIS detecting an intermediate mass black hole (IMBH) in a Milky Way globular cluster. We use an N-body simulation of an NGC 3201-like cluster with a central 1500 M$_{\odot}$ IMBH as input to MAVISIM and recover the velocity dispersion profile from proper motion measurements. Under favourable astrometric conditions, the dynamical signature of the IMBH is detected with a precision of ~0.20 km/s in the inner ~4" of the cluster where HST is confusion-limited. This precision is comparable to measurements made by Gaia, HST and MUSE in the outer ~60" of the cluster. This study is the first step towards building a science-driven astrometric error budget for an MCAO system and a prediction of what MAVIS could do once on sky., Comment: 17 pages, 13 figures, accepted for publication in MNRAS, main result is shown in Figure 12. Fixed typo in introduction

We discuss the largest and most homogeneous spectroscopic dataset of field RR Lyrae variables (RRLs) available to date. We estimated abundances using both high-resolution and low-resolution ({\Delta S} method) spectra for fundamental (RRab) and first overtone (RRc) RRLs. The iron abundances for 7,941 RRLs were supplemented with similar literature estimates available, ending up with 9,015 RRLs (6,150 RRab, 2,865 RRc). The metallicity distribution shows a mean value of <[Fe/H]> = -1.51\pm0.01, and {\sigma}(standard deviation)= 0.41 dex with a long metal-poor tail approaching [Fe/H] = -3 and a sharp metal-rich tail approaching solar iron abundance. The RRab variables are more metal-rich (<[Fe/H]>ab = -1.48\pm0.01, {\sigma} = 0.41 dex) than RRc variables (<[Fe/H]>c = -1.58\pm0.01, {\sigma} = 0.40 dex). The relative fraction of RRab variables in the Bailey diagram (visual amplitude vs period) located along the short-period (more metal-rich) and the long-period (more metal-poor) sequences are 80% and 20\%, while RRc variables display an opposite trend, namely 30\% and 70\%. We found that the pulsation period of both RRab and RRc variables steadily decreases when moving from the metal-poor to the metal-rich regime. The visual amplitude shows the same trend, but RRc amplitudes are almost two times more sensitive than RRab amplitudes to metallicity. We also investigated the dependence of the population ratio (Nc/Ntot) of field RRLs on the metallicity and we found that the distribution is more complex than in globular clusters. The population ratio steadily increases from ~0.25 to ~0.36 in the metal-poor regime, it decreases from ~0.36 to ~0.18 for -1.8 < [Fe/H] < -0.9 and it increases to a value of ~0.3 approaching solar iron abundance., Comment: 22 pages, 13 figures, 3 tables. Accepted for publication in ApJ

We collected the largest spectroscopic catalog of RR Lyrae (RRLs) including $\approx$20,000 high-, medium- and low-resolution spectra for $\approx$10,000 RRLs. We provide the analytical forms of radial velocity curve (RVC) templates. These were built using 36 RRLs (31 fundamental -- split into three period bins -- and 5 first overtone pulsators) with well-sampled RVCs based on three groups of metallic lines (Fe, Mg, Na) and four Balmer lines (H$_\alpha$, H$_\beta$, H$_\gamma$, H$_\delta$). We tackled the long-standing problem of the reference epoch to anchor light curve and RVC templates. For the $V$-band, we found that the residuals of the templates anchored to the phase of the mean magnitude along the rising branch are $\sim$35\% to $\sim$45\% smaller than those anchored to the phase of maximum light. For the RVC, we used two independent reference epochs for metallic and Balmer lines and we verified that the residuals of the RVC templates anchored to the phase of mean RV are from 30\% (metallic lines) up to 45\% (Balmer lines) smaller than those anchored to the phase of minimum RV. We validated our RVC templates by using both the single- and the three-phase points approach. We found that barycentric velocities based on our RVC templates are two-three times more accurate than those available in the literature. We applied the current RVC templates to Balmer lines RVs of RRLs in the globular NGC~3201 collected with MUSE at VLT. We found the cluster barycentric RV of $V_{\gamma}$=496.89$\pm$8.37(error)$\pm$3.43 (standard deviation) km/s, which agrees well with literature estimates.

We provide the largest and most homogeneous sample of $\alpha$-element (Mg, Ca, Ti) and iron abundances for field RR Lyrae (RRLs, 162 variables) by using high-resolution spectra. The current measurements were complemented with similar abundances available in the literature for 46 field RRLs brought to our metallicity scale. We ended up with a sample of old (t$\ge$ 10 Gyr), low-mass stellar tracers (208 RRLs: 169 fundamental, 38 first overtone, 1 mixed mode) covering three dex in iron abundance (-3.00$\le$[Fe/H]$\le$0.24). We found that field RRLs are $\sim$0.3 dex more $\alpha$-poor than typical Halo tracers in the metal-rich regime, ([Fe/H]$\ge$-1.2) while in the metal-poor regime ([Fe/H]$\le$-2.2) they seem to be on average $\sim$0.1 dex more $\alpha$-enhanced. This is the first time that the depletion in $\alpha$-elements for solar iron abundances is detected on the basis of a large, homogeneous and coeval sample of old stellar tracers. Interestingly, we also detected a close similarity in the [$\alpha$/Fe] trend between $\alpha$-poor, metal-rich RRLs and red giants (RGs) in the Sagittarius dwarf galaxy as well as between $\alpha$-enhanced, metal-poor RRLs and RGs in ultra faint dwarf galaxies. These results are supported by similar elemental abundances for 46 field Horizontal Branch (HB) stars. These stars share with RRLs the same evolutionary phase and the same progenitors. This evidence further supports the key role that old stellar tracers play in constraining the early chemical enrichment of the Halo and, in particular, in investigating the impact that dwarf galaxies have had in the mass assembly of the Galaxy.

We present newly-calibrated period-$\phi_{31}$-[Fe/H] relations for fundamental mode RR Lyrae stars in the optical and, for the first time, mid-infrared. This work's calibration dataset provides the largest and most comprehensive span of parameter space to date with homogeneous metallicities from $-3<\textrm{[Fe/H]}<0.4$ and accurate Fourier parameters derived from 1980 ASAS-SN ($V$-band) and 1083 WISE (NEOWISE extension, $W1$ and $W2$ bands) RR Lyrae stars with well-sampled light curves. We compare our optical period-$\phi_{31}$-[Fe/H] with those available in the literature and demonstrate that our relation minimizes systematic trends in the lower and higher metallicity range. Moreover, a direct comparison shows that our optical photometric metallicities are consistent with both those from high-resolution spectroscopy and globular clusters, supporting the good performance of our relation. We found an intrinsic scatter in the photometric metallicities (0.41 dex in the $V$-band and 0.50 dex in the infrared) by utilizing large calibration datasets covering a broad metallicity range. This scatter becomes smaller when optical and infrared bands are used together (0.37 dex). Overall, the relations derived in this work have many potential applications, including large-area photometric surveys with JWST in the infrared and LSST in the optical., Comment: Accepted by ApJ, 29 pages, 14 Figures, 4 Tables

The separability of clusters is one of the most desired properties in clustering. There is a wide range of settings in which different clusterings of the same data set appear. We are interested in applications where there is a need for an explicit, gradual transition of one separable clustering into another one. This transition should be a sequence of simple, natural steps that upholds separability of the clusters throughout. We design an algorithm for such a transition. We exploit the intimate connection of separability and linear programming over bounded-shape partition and transportation polytopes: separable clusterings lie on the boundary of partition polytopes, form a subset of the vertices of the corresponding transportation polytopes, and circuits of both polytopes are readily interpreted as sequential or cyclical exchanges of items between clusters. This allows for a natural approach to achieve the desired transition through a combination of two walks: an edge walk between two so-called radial clusterings in a transportation polytope, computed through an adaptation of classical tools of sensitivity analysis and parametric programming; and a walk from a separable clustering to a corresponding radial clustering, computed through a tailored, iterative routine updating cluster sizes and re-optimizing the cluster assignment of items.

We performed the largest and most homogeneous spectroscopic survey of field RR Lyraes (RRLs). We secured $\approx$6,300 high resolution (HR, R$\sim$35,000) spectra for 143 RRLs (111 fundamental, RRab; 32 first overtone, RRc). The atmospheric parameters were estimated by using the traditional approach and the iron abundances were measured by using an LTE line analysis. The resulting iron distribution shows a well defined metal-rich tail approaching solar iron abundance. This suggests that field RRLs experienced a complex chemical enrichment in the early halo formation. We used these data to develop a new calibration of the $\Delta$S method. This diagnostic, based on the equivalent widths of CaII K and three Balmer (H$_{\delta,\gamma,\beta}$) lines, traces the metallicity of RRLs. For the first time the new empirical calibration: i) includes spectra collected over the entire pulsation cycle; ii) includes RRc variables; iii) relies on spectroscopic calibrators covering more than three dex in iron abundance; iv) provides independent calibrations based on one/two/three Balmer lines. The new calibrations were applied to both SEGUE-SDSS and degraded HR spectra totalling 6,451 low resolution (LR, R$\sim$2,000) spectra for 5,001 RRLs (3,439 RRab, 1,562 RRc). This resulted in an iron distribution with a median of -1.55$\pm$0.01 and $\sigma$=0.51 dex, in good agreement with literature values. We also found that RRc are 0.10 dex more metal-poor than RRab variables, and have a distribution with a smoother metal-poor tail. This finding supports theoretical prescriptions suggesting a steady decrease in the RRc number when moving from metal-poor to metal-rich stellar environments., Comment: Accepted by ApJ

The separation between RR Lyrae (RRLs) and Type II Cepheid (T2Cs) variables based on their period is debated. Both types of variable stars are distance indicators and we aim to promote the use of T2Cs as distance indicators in synergy with RRLs. We adopted new and existing optical and Near-Infrared (NIR) photometry of \wcen~to investigate several diagnostics (colour-magnitude diagram, Bailey diagram, Fourier decomposition of the light curve, amplitude ratios) for their empirical separation. We found that the classical period threshold at 1 day is not universal and does not dictate the evolutionary stage: V92 has a period of 1.3 days but is likely to be still in its core Helium-burning phase, typical of RRLs. We also derived NIR Period-Luminosity relations and found a distance modulus of 13.65$\pm$0.07 (err.)$\pm$0.01 ($\sigma$) mag, in agreement with the recent literature. We also found that RRLs and T2Cs obey the same PL relations in the NIR. This equivalence gives the opportunity to adopt RRLs+T2Cs as an alternative to classical Cepheids to calibrate the extragalactic distance scale.

Our society has never been more dependent on computer networks. Effective utilization of networks requires a detailed understanding of the normal background behaviors of network traffic. Large-scale measurements of networks are computationally challenging. Building on prior work in interactive supercomputing and GraphBLAS hypersparse hierarchical traffic matrices, we have developed an efficient method for computing a wide variety of streaming network quantities on diverse time scales. Applying these methods to 100,000,000,000 anonymized source-destination pairs collected at a network gateway reveals many previously unobserved scaling relationships. These observations provide new insights into normal network background traffic that could be used for anomaly detection, AI feature engineering, and testing theoretical models of streaming networks., Comment: 6 pages, 6 figures,3 tables, 49 references, accepted to IEEE HPEC 2020

Six decades and counting, the formation of hot ~20,000-30,000 K Extreme Horizontal Branch (EHB) stars in Galactic Globular Clusters remains one of the most elusive quests in stellar evolutionary theory. Here we report on two discoveries shattering their currently alleged stable luminosity. The first EHB variability is periodic and cannot be ascribed to binary evolution nor pulsation. Instead, we here attribute it to the presence of magnetic spots: superficial chemical inhomogeneities whose projected rotation induces the variability. The second EHB variability is aperiodic and manifests itself on time-scales of years. In two cases, the six-year light curves display superflare events a mammoth several million times more energetic than solar analogs. We advocate a scenario where the two spectacular EHB variability phenomena are different manifestations of diffuse, dynamo-generated, weak magnetic fields. Ubiquitous magnetic fields, therefore, force an admittance into the intricate matrix governing the formation of all EHBs, and traverse to their Galactic field counterparts. The bigger picture is one where our conclusions bridge similar variability/magnetism phenomena in all radiative-enveloped stars: young main-sequence stars, old EHBs and defunct white dwarfs., Comment: Author's version of the main article (23 pages) and Supplementary Information (22 pages) combined into a single pdf (45 pages). Readers invited to read the Nature Astronomy Published version available at this url: https://www.nature.com/articles/s41550-020-1113-4

We collected over 6000 high-resolution spectra of four dozen field RR Lyrae (RRL) variables pulsating either in the fundamental (39 RRab) or in the first overtone (9 RRc) mode. We measured radial velocities (RVs) of four strong metallic and four Balmer lines along the entire pulsational cycle and derived RV amplitudes with accuracies better than 1$-$2~\kmsec. The new amplitudes were combined with literature data for 23~RRab and 3~RRc stars (total sample 74 RRLs) which allowed us to investigate the variation of the Bailey diagram (photometric amplitude versus period) when moving from optical to mid-infrared bands and to re-cast the Bailey diagram in terms of RV amplitudes. We found that RV amplitudes for RRab are minimally affected by nonlinear phenomena (shocks) and multi-periodicity (Blazhko effect). The RV slope ($\log P$--A(V$_r$)) when compared with the visual slope ($\log P$--A($V$)) is shallower and the dispersion, at fixed period, decreases by a factor of two. We constructed homogeneous sets of Horizontal Branch evolutionary models and nonlinear, convective pulsation models of RRLs to constrain the impact of evolutionary effects on their pulsation properties. Evolution causes, on the Bailey diagram based on RV amplitudes, a modest variation in pulsation period and a large dispersion in amplitude. The broad dispersion in period of the Bailey diagram is mainly caused by variation in RRL intrinsic parameters (stellar mass, chemical composition). Empirical evidence indicates that RV amplitudes are an optimal diagnostic for tracing the mean effective temperature across the RRab instability strip., Comment: 14 pages, 4 figures, 1 table, Accepted on ApJ Letter

The galactic globular cluster Omega Centauri is the most massive of its kind, with a complex mix of multiple stellar populations and several kinematic and dynamical peculiarities. Different mean proper motions have been detected among the three main sub-populations, implying that the most metal-rich one is of accreted origin. This particular piece of evidence has been a matter of debate because the available data have either not been sufficiently precise or limited to a small region of the cluster to ultimately confirm or refute the result. Using astrometry from the second Gaia data release and recent high-quality, multi-band photometry, we are now in a position to resolve the controversy. We reproduced the original analysis using the Gaia data and found that the three populations have the same mean proper motion. Thus, there is no need to invoke an accreted origin for the most metal-rich sub-population., Comment: 6 pages, 2 tables, 5 figures. Accepted for publication on A&A

Time series observations of a single dithered field centered on the diffuse dwarf satellite galaxy Crater II were obtained with the Dark Energy Camera (DECam) at the 4m Blanco Telescope at Cerro Tololo Inter-American Observatory, Chile, uniformly covering up to two half-light radii. Analysis of the $g$ and $i$ time series results in the identification and characterization of 130 periodic variable stars, including 98 RR Lyrae stars, 7 anomalous Cepheids, and 1 SX Phoenicis star belonging to the Crater II population, and 24 foreground variables of different types. Using the large number of ab-type RR Lyrae stars present in the galaxy, we obtained a distance modulus to Crater II of $(m-M)_0=20.333\pm 0.004$ (stat) $\pm 0.07$ (sys). The distribution of the RR Lyrae stars suggests an elliptical shape for Crater II, with an ellipticity of 0.24 and a position angle of $153^\circ$. From the RR Lyrae stars we infer a small metallicity dispersion for the old population of Crater II of only 0.17 dex. There are hints that the most metal-poor stars in that narrow distribution have a wider distribution across the galaxy, while the slightly more metal rich part of the population is more centrally concentrated. Given the features in the color-magnitude diagram of Crater II, the anomalous Cepheids in this galaxy must have formed through a binary evolution channel of an old population., Comment: Accepted for publication in MNRAS

We present a deep Blanco/DECam colour-magnitude diagram (CMD) for the large but very diffuse Milky Way satellite dwarf galaxy Crater II. The CMD shows only old stars with a clearly bifurcated subgiant branch (SGB) that feeds a narrow red giant branch. The horizontal branch (HB) shows many RR Lyrae and red HB stars. Comparing the CMD with [Fe/H] = -2.0 and [$\alpha$/Fe] = +0.3 alpha-enhanced BaSTI isochrones indicates a mean age of 12.5 Gyr for the main event and a mean age of 10.5 Gyr for the brighter SGB. With such multiple star formation events Crater II shows similarity to more massive dwarfs that have intermediate age populations, however for Crater II there was early quenching of the star formation and no intermediate age or younger stars are present. The spatial distribution of Crater II stars overall is elliptical in the plane of the sky, the detailed distribution shows a lack of strong central concentration, and some inhomogeneities. The 10.5 Gyr subgiant and upper main sequence stars show a slightly higher central concentration when compared to the 12.5 Gyr population. Matching to Gaia DR2 we find the proper motion of Crater II: $\mu_{\alpha}\cos \delta$=-0.14 $\pm$ 0.07 , $\mu_{\delta}$=-0.10 $\pm$ 0.04 mas yr$^{-1}$, approximately perpendicular to the semi-major axis of Crater II. Our results provide constraints on the star formation and chemical enrichment history of Crater II, but cannot definitively determine whether or not substantial mass has been lost over its lifetime., Comment: 13 pages, 9 figures, 3 tables. Accepted for publication in MNRAS

Theory and observations indicate that the distribution of stars along the horizontal branch of Galactic globular clusters mainly depends on the metal content. However, the existence of globular clusters with similar metal content and absolute age but different horizontal branch morphologies, suggests the presence of another parameter affecting the star distribution along the branch. To investigate the variation of the horizontal branch morphology in Galactic globular clusters, we define a new photometric horizontal branch morphology index, overcoming some of the limitations and degeneracies affecting similar indices available in the literature. We took advantage of a sample of 64 Galactic globular clusters, with both space-based imaging data and homogeneous ground-based photometric catalogues in five different bands ($U$,$B$,$V$, $R$, $I$). The new index, $\tau_{HB}$, is defined as the ratio between the areas subtended by the cumulative number distribution in magnitude ($I$) and in colour ($V-I$) of all stars along the horizontal branch. This new index shows a linear trend over the entire range in metallicity (-2.35 $\leq$ [Fe/H] $\leq$ -0.12) covered by our Galactic globular cluster sample. We found a linear relation between $\tau_{HB}$ and absolute cluster ages. We also found a quadratic anti-correlation with [Fe/H], becoming linear when we eliminate the age effect on $\tau_{HB}$ values. Moreover, we identified a subsample of eight clusters that are peculiar according to their $\tau_{HB}$ values. These clusters have bluer horizontal branch morphology when compared to typical ones of similar metallicity. These findings allow us to define them as the 'second parameter' clusters in the sample. A comparison with synthetic horizontal branch models suggests that they cannot be entirely explained with a spread in helium content., Comment: 18 pages, 21 figures, accepted for publication in Astronomy & Astrophysics

We directly observe the hydrodynamic linear response of a unitary Fermi gas confined in a box potential and subject to a spatially periodic optical potential that is translated into the cloud at speeds ranging from subsonic to supersonic. We show that the time-dependent change of the density profile is sensitive to the thermal conductivity, which controls the relaxation rate of the temperature gradients and hence the responses arising from adiabatic and isothermal compression., Comment: 12 pages, 7 figures