• We present HST photometry of a selected sample of 50 long-period, low-extinction Milky Way Cepheids measured on the same WFC3 F555W, F814W, and F160W-band photometric system as extragalactic Cepheids in SN Ia hosts. These bright Cepheids were observed with the WFC3 spatial scanning mode in the optical and near-infrared to mitigate saturation and reduce pixel-to-pixel calibration errors to reach a mean photometric error of 5 millimags per observation. We use the new Gaia DR2 parallaxes and HST photometry to simultaneously constrain the cosmic distance scale and to measure the DR2 parallax zeropoint offset appropriate for Cepheids. We find a value for the zeropoint offset of -46 +/- 13 muas or +/- 6 muas for a fixed distance scale, higher than found from quasars, as expected, for these brighter and redder sources. The precision of the distance scale from DR2 has been reduced by a factor of 2.5 due to the need to independently determine the parallax offset. The best fit distance scale is 1.006 +/- 0.033, relative to the scale from Riess et al 2016 with H0=73.24 km/s/Mpc used to predict the parallaxes photometrically, and is inconsistent with the scale needed to match the Planck 2016 CMB data combined with LCDM at the 2.9 sigma confidence level (99.6%). At 96.5% confidence we find that the formal DR2 errors may be underestimated as indicated. We identify additional error associated with the use of augmented Cepheid samples utilizing ground-based photometry and discuss their likely origins. Including the DR2 parallaxes with all prior distance ladder data raises the current tension between the late and early Universe route to the Hubble constant to 3.8 sigma (99.99 %). With the final expected precision from Gaia, the sample of 50 Cepheids with HST photometry will limit to 0.5% the contribution of the first rung of the distance ladder to the uncertainty in the Hubble constant.
  • We present year-long, near-infrared Hubble Space Telescope WFC3 observations of Mira variables in the water megamaser host galaxy NGC 4258. Miras are AGB variables that can be divided into oxygen- (O-) and carbon- (C-) rich subclasses. Oxygen-rich Miras follow a tight (scatter $\sim 0.14$ mag) Period-Luminosity Relation (PLR) in the near-infrared and can be used to measure extragalactic distances. The water megamaser in NGC 4258 gives a geometric distance to the galaxy accurate to 2.6% that can serve to calibrate the Mira PLR. We develop criteria for detecting and classifying O-rich Miras with optical and NIR data as well as NIR data alone. In total, we discover 438 Mira candidates that we classify with high confidence as O-rich. Our most stringent criteria produce a sample of 139 Mira candidates that we use to measure a PLR. We use the OGLE-III sample of O-rich Miras in the LMC to obtain a relative distance modulus, $\mu_{4258} - \mu_{LMC} = 10.95 \pm 0.01 $ (statistical) $\pm 0.06 $ (systematic) mag in good agreement with the relative distance determined using Cepheids. These results demonstrate the feasibility of discovering and characterizing Miras using the near-infrared with the Hubble Space Telescope and the upcoming James Webb Space Telescope and using them to measure extragalactic distances and determine the Hubble constant.
  • We present new parallax measurements of 7 long-period (> 10 days) Milky Way Cepheids (SS CMa, XY Car, VY Car, VX Per, WZ Sgr, X Pup and S Vul) using astrometry from spatial scanning of WFC3 on HST. Observations were obtained at 6 month intervals over 4 years. The distances are 1.7--3.6 kpc with a mean precision of 45 microarcseconds and a best of 29 microarcseconds (SNR = 14). The accuracy of the parallaxes is demonstrated through independent analyses of >100 reference stars. This raises to 10 the number of long-period Cepheids with significant parallax measurements, 8 obtained from this program. We also present high-precision F555W, F814W, and F160W magnitudes of these Cepheids, allowing a direct, zeropoint-independent comparison to >1800 extragalactic Cepheids in the hosts of 19 SNeIa. This sample addresses two outstanding systematic uncertainties affecting prior comparisons of Milky Way and extragalactic Cepheids used to calibrate H_0: their dissimilarity of periods and photometric systems. Comparing the new parallaxes to their predicted values derived from reversing the distance ladder gives a ratio (or independent scale for H_0) of 1.037+/-0.036, consistent with no change and inconsistent at the 3.5 sigma level with a ratio of 0.91 needed to match the value predicted by Planck+LCDM. Using these data instead to augment the Riess et al. (2016) measurement of H_0 improves the precision to 2.3%, yielding 73.48+/-1.66 km/s/Mpc, and tension with Planck+LCDM increases to 3.7 sigma. The future combination of Gaia parallaxes and HST spatial scanning photometry of 50 Milky Way Cepheids can support a < 1% calibration of H_0.
  • We study the near-infrared properties of 690 Mira candidates in the central region of the Large Magellanic Cloud, based on time-series observations at JHKs. We use densely-sampled I-band observations from the OGLE project to generate template light curves in the near infrared and derive robust mean magnitudes at those wavelengths. We obtain near-infrared Period-Luminosity relations for Oxygen-rich Miras with a scatter as low as 0.12 mag at Ks. We study the Period-Luminosity-Color relations and the color excesses of Carbon-rich Miras, which show evidence for a substantially different reddening law.
  • We present the discovery of 1847 Mira candidates in the Local Group galaxy M33 using a novel semi-parametric periodogram technique coupled with a Random Forest classifier. The algorithms were applied to ~2.4x10^5 I-band light curves previously obtained by the M33 Synoptic Stellar Survey. We derive preliminary Period-Luminosity relations at optical, near- & mid-infrared wavelengths and compare them to the corresponding relations in the Large Magellanic Cloud.
  • We develop a non-linear semi-parametric Gaussian process model to estimate periods of Miras with sparsely-sampled light curves. The model uses a sinusoidal basis for the periodic variation and a Gaussian process for the stochastic changes. We use maximum likelihood to estimate the period and the parameters of the Gaussian process, while integrating out the effects of other nuisance parameters in the model with respect to a suitable prior distribution obtained from earlier studies. Since the likelihood is highly multimodal for period, we implement a hybrid method that applies the quasi-Newton algorithm for Gaussian process parameters and search the period/frequency parameter over a dense grid. A large-scale, high-fidelity simulation is conducted to mimic the sampling quality of Mira light curves obtained by the M33 Synoptic Stellar Survey. The simulated data set is publicly available and can serve as a testbed for future evaluation of different period estimation methods. The semi-parametric model outperforms an existing algorithm on this simulated test data set as measured by period recovery rate and quality of the resulting Period-Luminosity relations.
  • We present results of an optical search for Cepheid variable stars using the Hubble Space Telescope (HST) in 19 hosts of Type Ia supernovae (SNe Ia) and the maser-host galaxy NGC 4258, conducted as part of the SH0ES project (Supernovae and H0 for the Equation of State of dark energy). The targets include 9 newly imaged SN Ia hosts using a novel strategy based on a long-pass filter that minimizes the number of HST orbits required to detect and accurately determine Cepheid properties. We carried out a homogeneous reduction and analysis of all observations, including new universal variability searches in all SN Ia hosts, that yielded a total of 2200 variables with well-defined selection criteria -- the largest such sample identified outside the Local Group. These objects are used in a companion paper to determine the local value of H0 with a total uncertainty of 2.4%.
  • We present the results of the optical follow-up conducted by the TOROS collaboration of the first gravitational-wave event GW150914. We conducted unfiltered CCD observations (0.35-1 micron) with the 1.5-m telescope at Bosque Alegre starting ~2.5 days after the alarm. Given our limited field of view (~100 square arcmin), we targeted 14 nearby galaxies that were observable from the site and were located within the area of higher localization probability. We analyzed the observations using two independent implementations of difference-imaging algorithms, followed by a Random-Forest-based algorithm to discriminate between real and bogus transients. We did not find any bona fide transient event in the surveyed area down to a 5-sigma limiting magnitude of r=21.7 mag (AB). Our result is consistent with the LIGO detection of a binary black hole merger, for which no electromagnetic counterparts are expected, and with the expected rates of other astrophysical transients.
  • We use the Wide Field Camera 3 (WFC3) on the Hubble Space Telescope (HST) to reduce the uncertainty in the local value of the Hubble constant (H_0) from 3.3% to 2.4%. Improvements come from new, near-infrared observations of Cepheid variables in 11 new hosts of recent SNe~Ia, more than doubling the sample of SNe~Ia having a Cepheid-calibrated distance for a total of 19; these leverage the magnitude-z relation based on 300 SNe~Ia at z<0.15. All 19 hosts and the megamaser system NGC4258 were observed with WFC3, thus nullifying cross-instrument zeropoint errors. Other improvements include a 33% reduction in the systematic uncertainty in the maser distance to NGC4258, more Cepheids and a more robust distance to the LMC from late-type DEBs, HST observations of Cepheids in M31, and new HST-based trigonometric parallaxes for Milky Way (MW) Cepheids. We consider four geometric distance calibrations of Cepheids: (i) megamasers in NGC4258, (ii) 8 DEBs in the LMC, (iii) 15 MW Cepheids with parallaxes, and (iv) 2 DEBs in M31. H_0 from each is 72.25+/-2.51, 72.04+/-2.67, 76.18+/-2.37, and 74.50+/-3.27 km/sec/Mpc, respectively. Our best estimate of 73.24+/-1.74 km/sec/Mpc combines the anchors NGC4258, MW, and LMC, and includes systematic errors for a final uncertainty of 2.4%. This value is 3.4 sigma higher than 66.93+/-0.62 km/sec/Mpc predicted by LambdaCDM with 3 neutrinos with mass 0.06 eV and the Planck data, but reduces to 2.1 sigma relative to the prediction of 69.3+/-0.7 km/sec/Mpc with the combination of WMAP+ACT+SPT+BAO, suggesting systematic uncertainties in CMB measurements may play a role in the tension. If we take the conflict between Planck and H_0 at face value, one plausible explanation could involve an additional source of dark radiation in the early Universe in the range of Delta N_eff=0.4-1. We anticipate significant improvements in H_0 from upcoming parallax measurements.
  • Based on $J/\psi$ events collected with the BESIII detector, with corresponding Monte Carlo samples, the tracking efficiency and its systematic uncertainty are studied using a control sample of $J/\psi \to p \overline p \pi^+ \pi^-$. Validation methods and different factors influencing the tracking efficiency are presented in detail. The tracking efficiency and its systematic uncertainty for protons and pions with the transverse momentum and polar angle dependence are also discussed.
  • In this paper, we report the detections of stellar variabilities from the first 2-year observations of sky area of about 1300 square degrees from the Tsinghua University-NAOC Transient Survey (TNTS). A total of 1237 variable stars (including 299 new ones) were detected with brightness < 18.0 mag and magnitude variation >= 0.1 mag on a timescale from a few hours to few hundred days. Among such detections, we tentatively identified 661 RR Lyrae stars, 431 binaries, 72 Semiregular pulsators, 29 Mira stars, 11 slow irregular variables, 11 RS Canum Venaticorum stars, 7 Gamma Doradus stars, 5 long period variables, 3 W Virginis stars, 3 Delta Scuti stars, 2 Anomalous Cepheids, 1 Cepheid, and 1 nove-like star based on their time-series variability index Js and their phased diagrams. Moreover, we found that 14 RR Lyrae stars show the Blazhko effect and 67 contact eclipsing binaries exhibit the O'Connell effect. Since the period and amplitude of light variations of RR Lyrae variables depend on their chemical compositions, their photometric observations can be used to investigate distribution of metallicity along the direction perpendicular to the Galactic disk. We find that the metallicity of RR Lyrae stars shows large scatter at regions closer to the Galactic plane (e.g., -3.0 < [Fe/H] < 0) but tends to converge at [Fe/H]~ -1.7 at larger Galactic latitudes. This variation may be related to that the RRAB Lyrae stars in the Galactic halo come from globular clusters with different metallicity and vertical distances, i.e. OoI and OoII populations, favoring for the dual-halo model.