Search

Draco dwarf Spheroidal galaxy (dSph) is one of the nearest and the most dark matter dominated satellites of the Milky Way. We obtained multi-epoch near-infrared (NIR, $JHK_s$) observations of the central region of Draco dSph covering a sky area of $\sim 21'\times21'$ using the WIRCam instrument at the 3.6-m Canada-France-Hawaii Telescope. Homogeneous $JHK_s$ time-series photometry for 212 RR Lyrae (173 fundamental-mode, 24 first-overtone, and 15 mixed-mode variables) and 5 Anomalous Cepheids in Draco dSph is presented and used to derive their period-luminosity relations at NIR wavelengths for the first-time. The small scatter of $\sim 0.05$~mag in these empirical relations for RR Lyrae stars is consistent with those in globular clusters and suggests a very small metallicity spread, up to $\sim0.2$~dex, among these centrally located variables. Based on empirically calibrated NIR period-luminosity-metallicity relations for RR Lyrae in globular clusters, we determined a distance modulus to Draco dSph of $\mu_\textrm{RRL} = 19.557 \pm 0.026$ mag. The calibrated $K_s$-band period-luminosity relations for Anomalous Cepheids in the Draco dSph and the Large Magellanic Cloud exhibit statistically consistent slopes but systematically different zero-points, hinting at possible metallicity dependence of $\sim-0.3$ mag~dex$^{-1}$. Finally, the apparent magnitudes of the tip of the red giant branch in $I$ and $J$ bands also agree well with their absolute calibrations with the adopted RR Lyrae distance to Draco. Our recommended $\sim1.5\%$ precise RR Lyrae distance, $D_\textrm{Draco} = 81.55 \pm 0.98 \textrm{(statistical)} \pm 1.17 \textrm{(systematic)}$~kpc, is the most accurate and precise distance to Draco dSph galaxy., Comment: Accepted for publication in The Astronomical Journal

Milky Way Cepheid variables with accurate {\it Hubble Space Telescope} photometry have been established as standards for primary calibration of the cosmic distance ladder to achieve a percent-level determination of the Hubble constant ($H_0$). These 75 Cepheid standards are the fundamental sample for investigation of possible residual systematics in the local $H_0$ determination due to metallicity effects on their period-luminosity relations. We obtained new high-resolution ($R\sim81,000$), high signal-to-noise ($S/N\sim50-150$) multi-epoch spectra of 42 out of 75 Cepheid standards using ESPaDOnS instrument at the 3.6-m Canada-France-Hawaii Telescope. Our spectroscopic metallicity measurements are in good agreement with the literature values with systematic differences up to $0.1$ dex due to different metallicity scales. We homogenized and updated the spectroscopic metallicities of all 75 Milky Way Cepheid standards and derived their multiwavelength ($GVIJHK_s$) period-luminosity-metallicity and period-Wesenheit-metallicity relations using the latest {\it Gaia} parallaxes. The metallicity coefficients of these empirically calibrated relations exhibit large uncertainties due to low statistics and a narrow metallicity range ($\Delta\textrm{[Fe/H]}=0.6$~dex). These metallicity coefficients are up to three times better constrained if we include Cepheids in the Large Magellanic Cloud and range between $-0.21\pm0.07$ and $-0.43\pm0.06$ mag/dex. The updated spectroscopic metallicities of these Milky Way Cepheid standards were used in the Cepheid-Supernovae distance ladder formalism to determine $H_0=72.9~\pm 1.0$\textrm{~km~s$^{-1}$~Mpc$^{-1}$}, suggesting little variation ($\sim 0.1$ ~km~s$^{-1}$~Mpc$^{-1}$) in the local $H_0$ measurements due to different Cepheid metallicity scales., Comment: 14 pages, 4 figures, Accepted for publication in The Astrophysical Journal Letters

Based on 14 Miras located in 7 globular clusters, we derived the first gr-band period-luminosity (PL) at maximum light for the large-amplitude Mira variables using the multi-year light-curve data collected from the Zwicky Transient Facility (ZTF). Since Miras are red variables, we applied a color-term correction to subsets of ZTF light curves, and found that such corrections do not have a large impact on period determinations. We applied our derived PL relations to the known extragalactic Miras in five local galaxies (Sextans, Leo I, Leo II, NGC6822 and IC1613), and determined their Mira-based distances. We demonstrated that our PL relations can be applied to short-period (<300 days) Miras, including those in the two most distant galaxies (NGC6822 and IC1613) in our sample even when only a portion of the light-curves around maximum light have detections. We have also shown that the long-period extragalactic Miras do not follow the PL relations extrapolated to longer periods. Hence, our derived PL relations are only applicable to the short-period Miras, which will be discovered in abundance in local galaxies within the era of Vera C. Rubin Observatory's Legacy Survey of Space and Time., Comment: 12 pages, 2 Tables and 10 Figures, AJ accepted

SX Phoenicis (SXP) variables are short period pulsating stars that exhibit a period-luminosity (PL) relation. We derived the gri-band PL and extinction-free period-Wesenheit (PW) relations, as well as the period-color (PC) and reddening-free period-Q-index (PQ) relations for 47 SXP stars in located in 21 globular clusters using the optical light curves taken from Zwicky Transient Facility (ZTF). These empirically relations were derived for the first time in the gri filters except for the g-band PL relation. We used our gi band PL and PW relations to derive a distance modulus to Crater II dwarf spheroidal which hosts one SXP variable. Assuming that the fundamental and first-overtone pulsation mode for the SXP variable in Crater II, we found distance moduli of $20.03 \pm 0.23$ mag and $20.37 \pm 0.24$ mag, respectively, using the PW relation, where the latter is in excellent agreement with independent RR Lyrae based distance to Crater II dwarf galaxy., Comment: 12 pages, 1 Table and 10 Figures; AJ accepted

RR Lyrae variables are excellent population II distance indicators thanks to their well-defined period-luminosity relations (PLRs) at infrared wavelengths. We present results of near-infrared (NIR) monitoring of Galactic globular clusters to empirically quantify the metallicity dependence of NIR PLRs for RR Lyrae variables. Our sample includes homogeneous, accurate, and precise photometric data for 964 RR Lyrae variables in 11 globular clusters covering a large metallicity range ($\Delta\textrm{[Fe/H]}\sim2$~dex). We derive $JHK_s$ band period-luminosity-metallicity (PLZ) and period-Wesenheit-metallicity (PWZ) relations anchored using 346 Milky Way field RR Lyrae stars with {\it Gaia} parallaxes, and simultaneously solved for independent distances to globular clusters. We find a significant metallicity dependence of $\sim0.2$~mag/dex in $JHK_s$ band PLZ and PWZ relations for RR Lyrae stars independent of the adopted metallicity scale. The metallicity coefficients and the zero-points of the empirical PLZ and PWZ relations are in excellent agreement with the predictions from the horizontal branch evolution and pulsation models. Furthermore, RR Lyrae based distances to our sample of globular clusters are also statistically consistent with other independent measurements in the literature. Our recommended empirical $JHK_s$ band PLZ relations are also provided for RR Lyrae based distance measurements., Comment: 10 pages, 4 figures, abridged abstract, accepted in the ApJ Letters

We present a systematic analysis to determine and improve the pulsation periods of 1637 known long-period Mira variables in M33 using $gri$-band light curves spanning $\sim18$~years from several surveys, including M33 variability survey, Panoramic Survey Telescope and Rapid Response System, Palomar Transient Factory (PTF), intermediate PTF, and Zwicky Transient Facility. Based on these collections of light curves, we found that optical band light curves that are as complete as possible are crucial to determine the periods of distant Miras. We demonstrated that the machine learning techniques can be used to classify Miras into O-rich and C-rich based on the $(J-K_s)$ period--color plane. Finally, We derived the distance modulus to M33 using O-rich Miras at maximum light together with our improved periods as $24.67 \pm 0.06$~mag, which is in good agreement with the recommended value given in the literature.

Earth-sized exoplanets that transit nearby, late spectral type red dwarfs will be prime targets for atmospheric characterization in the coming decade. Such systems, however, are difficult to find via wide-field transit surveys like Kepler or TESS. Consequently, the presence of such transiting planets is unexplored and the occurrence rates of short-period Earth-sized planets around late M dwarfs remain poorly constrained. Here, we present the deepest photometric monitoring campaign of 22 nearby late M dwarf stars, using data from over 500 nights on seven 1-2 meter class telescopes. Our survey includes all known single quiescent northern late M dwarfs within 15 pc. We use transit-injection-and-recovery tests to quantify the completeness of our survey, successfully identify most ($>80\%$) transiting short-period (0.5-1 d) super-Earths ($R > 1.9 R_\oplus$), and are sensitive ($\sim50\%$) to transiting Earth-sized planets ($1.0-1.2 R_\oplus$). Our high sensitivity to transits with a near-zero false positive rate demonstrates an efficient survey strategy. Our survey does not yield a transiting planet detection, yet it provides the most sensitive upper limits on transiting planets orbiting our target stars. Finally, we explore multiple hypotheses about the occurrence rates of short-period planets (from Earth-sized planets to giant planets) around late M dwarfs. We show, for example, that giant planets at short periods ($<1$ day) are uncommon around our target stars. Our dataset provides some insight into occurrence rates of short-period planets around TRAPPIST-1-like stars, and our results can help test planetary formation and system evolution models, as well as guide future observations of nearby late M dwarfs., Comment: 27 pages, 11 figures

Multi-band light curves of two RR Lyrae variables in Segue II and Ursa Major II ultra-faint dwarf (UFD) galaxies were collected from near simultaneous observations using the Lulin One-meter Telescope in Vgri bands. Together with Gaia G-band light curves, we determined photometric metallicities using empirical relations involving pulsation period and Fourier parameter as dependent parameters. We demonstrated that the RR Lyrae photometric metallicity can be determined accurately when these empirical relations were employed at multiple wavelengths, which can potentially improve the distance determination based on RR Lyrae stars. The photometric metallicities based on our approach were found to be $-2.27\pm0.13$ dex and $-1.87\pm0.16$ dex for the RR Lyrae in Segue II and Ursa Major II UFD, respectively, with corresponding distance moduli of $17.69\pm0.15$ mag and $17.58\pm0.15$ mag, in agreement with previous literature determinations. This approach of photometric metallicity of RR Lyrae based on multi-band optical light curves will be particularly relevant for distance measurements in the era of the Vera C Rubin's Legacy Survey of Space and Time., Comment: 9 pages, 3 figures and 5 tables, AJ accepted

RR Lyrae stars are pulsating stars, many of which also show a long-term variation called the Blazhko effect which is a modulation of amplitude and phase of the lightcurve. In this work, we searched for the incidence rate of the Blazhko effect in the first-overtone pulsating RR Lyrae (RRc) stars of the Galactic halo. The focus was on the Stripe 82 region in the Galactic halo which was studied by Sesar et al using the Sloan Digital Sky Survey (SDSS) data. In their work, 104 RR Lyrae stars were classified as RRc type. We combined their SDSS light curves with Zwicky Transient Facility (ZTF) data, and use them to document the Blazhko properties of these RRc stars. Our analysis showed that among the 104 RRc stars, 8 were rather RRd stars, and were excluded from the study. Out of remaining 96, 34 were Blazhko type, 62 were non-Blazhko type, giving the incidence rate of 35.42% for Blazhko RRc stars. The shortest Blazhko period found was 12.808 +/- 0.001 d for SDSS 747380, while the longest was 3100 +/- 126 d for SDSS 3585856. Combining SDSS and ZTF data sets increased the probability of detecting the small variations due to the Blazhko effect, and thus provided a unique opportunity to find longer Blazhko periods. We found that 85% of RRc stars had the Blazhko period longer than 200 d., Comment: 9 pages, 2 tables, 8 figures, AJ accepted

In recent years, certain luminous extragalactic optical transients have been observed to last only a few days. Their short observed duration implies a different powering mechanism from the most common luminous extragalactic transients (supernovae) whose timescale is weeks. Some short-duration transients, most notably AT2018cow, display blue optical colours and bright radio and X-ray emission. Several AT2018cow-like transients have shown hints of a long-lived embedded energy source, such as X-ray variability, prolonged ultraviolet emission, a tentative X-ray quasiperiodic oscillation, and large energies coupled to fast (but subrelativistic) radio-emitting ejecta. Here we report observations of minutes-duration optical flares in the aftermath of an AT2018cow-like transient, AT2022tsd (the "Tasmanian Devil"). The flares occur over a period of months, are highly energetic, and are likely nonthermal, implying that they arise from a near-relativistic outflow or jet. Our observations confirm that in some AT2018cow-like transients the embedded energy source is a compact object, either a magnetar or an accreting black hole., Comment: 79 pages, 3 figures (main text) + 7 figures (extended data) + 2 figures (supplementary information). Published online in Nature on 15 November 2023

In recent years, certain luminous extragalactic optical transients have been observed to last only a few days 1. Their short observed duration implies a different powering mechanism from the most common luminous extragalactic transients (supernovae), whose timescale is weeks 2. Some short-duration transients, most notably AT2018cow (ref. 3), show blue optical colours and bright radio and X-ray emission 4. Several AT2018cow-like transients have shown hints of a long-lived embedded energy source 5, such as X-ray variability 6,7, prolonged ultraviolet emission 8, a tentative X-ray quasiperiodic oscillation 9,10 and large energies coupled to fast (but subrelativistic) radio-emitting ejecta 11,12. Here we report observations of minutes-duration optical flares in the aftermath of an AT2018cow-like transient, AT2022tsd (the ‘Tasmanian Devil’). The flares occur over a period of months, are highly energetic and are probably nonthermal, implying that they arise from a near-relativistic outflow or jet. Our observations confirm that, in some AT2018cow-like transients, the embedded energy source is a compact object, either a magnetar or an accreting black hole.

A black hole can launch a powerful relativistic jet after it tidally disrupts a star. If this jet fortuitously aligns with our line of sight, the overall brightness is Doppler boosted by several orders of magnitude. Consequently, such on-axis relativistic tidal disruption events have the potential to unveil cosmological (redshift z > 1) quiescent black holes and are ideal test beds for understanding the radiative mechanisms operating in super-Eddington jets. Here we present multiwavelength (X-ray, UV, optical and radio) observations of the optically discovered transient AT 2022cmc at z = 1.193. Its unusual X-ray properties, including a peak observed luminosity of ≳1048 erg s−1, systematic variability on timescales as short as 1,000 s and overall duration lasting more than 30 days in the rest frame, are traits associated with relativistic tidal disruption events. The X-ray to radio spectral energy distributions spanning 5–50 days after discovery can be explained as synchrotron emission from a relativistic jet (radio), synchrotron self-Compton (X-rays) and thermal emission similar to that seen in low-redshift tidal disruption events (UV/optical). Our modelling implies a beamed, highly relativistic jet akin to blazars but requires extreme matter domination (that is, a high ratio of electron-to-magnetic-field energy densities in the jet) and challenges our theoretical understanding of jets.

We measured the parallax of the James Webb Space Telescope based on near simultaneous observations using the Lulin One-meter Telescope and the GROWTH India Telescope, separated at a distance of ~4214 km. This serves a great demonstration for the concept of parallax commonly taught in introductory astronomy courses., Comment: 2 pages and 1 figure, RNAAS in-press

Based on time-series observations collected from Zwicky Transient Facility (ZTF), we derived period-luminosity-metallicity (PLZ) and period-Wesenheit-metallicity (PWZ) relations for RR Lyrae located in globular clusters. We have applied various selection criteria to exclude RR Lyrae with problematic or spurious light curves. These selection criteria utilized information on the number of data points per light curve, amplitudes, colors, and residuals on the period-luminosity and/or period-Wesenheit relations. Due to blending, a number of RR Lyrae in globular clusters were found to be anomalously bright and have small amplitudes of their ZTF light curves. We used our final sample of ~750 RR Lyrae in 46 globular clusters covering a wide metallicity range (-2.36 < [Fe/H] < -0.54 dex) to derive PLZ and PWZ relations in gri bands. In addition, we have also derived the period-color-metallicity (PCZ) and for the first time, the PQZ relations where the Q-index is extinction-free by construction. We have compared our various relations to empirical and theoretical relations available in literature, and found a good agreement with most studies. Finally, we applied our derived PLZ relation to a dwarf galaxy, Crater II, and found its true distance modulus should be larger than the most recent determination., Comment: 28 pages, 7 Tables and 24 Figures. AJ accepted

Mira variables are asymptotic giant branch pulsating stars with long pulsation periods and large amplitudes in optical bands. By applying the random forest algorithm to the I-band light curves for the Miras in the Magellanic Clouds, we have classified these Miras into regular Miras and non-regular Miras. Wherein non-regular Miras exhibit a long-term variation in addition to their primary pulsation periods. Our results confirm that the period-luminosity relation for maximum light has a smaller dispersion, but only occurs on the regular oxygen-rich (O-rich) Miras, which we recommend to be applied in future distance scale work. We have also collected multi-band photometry for these Miras to perform a spectral-energy-distribution (SED) fitting with stellar and dust components, showing that a significant fraction of dust is present around the non-regular Miras. According to our results, we believe that the periodic long-term variations seen in the non-regular Miras might be due to the presence of dust.

A black hole can launch a powerful relativistic jet after it tidally disrupts a star. If this jet fortuitously aligns with our line of sight, the overall brightness is Doppler boosted by several orders of magnitude. Consequently, such on-axis relativistic tidal disruption events (TDEs) have the potential to unveil cosmological (redshift $z>$1) quiescent black holes and are ideal test beds to understand the radiative mechanisms operating in super-Eddington jets. Here, we present multi-wavelength (X-ray, UV, optical, and radio) observations of the optically discovered transient \target at $z=1.193$. Its unusual X-ray properties, including a peak observed luminosity of $\gtrsim$10$^{48}$ erg s$^{-1}$, systematic variability on timescales as short as 1000 seconds, and overall duration lasting more than 30 days in the rest-frame are traits associated with relativistic TDEs. The X-ray to radio spectral energy distributions spanning 5-50 days after discovery can be explained as synchrotron emission from a relativistic jet (radio), synchrotron self-Compton (X-rays), and thermal emission similar to that seen in low-redshift TDEs (UV/optical). Our modeling implies a beamed, highly relativistic jet akin to blazars but requires extreme matter-domination, i.e, high ratio of electron-to-magnetic field energy densities in the jet, and challenges our theoretical understanding of jets., Comment: To appear in Nature Astronomy on 30th November 2022. Also see here for an animation explaining the result: https://youtu.be/MQHdSbxuznY

The transiting planet HD80606b undergoes a 1000-fold increase in insolation during its 111-day orbit due to it being highly eccentric (e=0.93). The planet's effective temperature increases from 400K to over 1400K in a few hours as it makes a rapid passage to within 0.03AU of its host star during periapsis. Spectroscopic observations during the eclipse (which is conveniently oriented a few hours before periapsis) of HD80606b with the James Webb Space Telescope (JWST) are poised to exploit this highly variable environment to study a wide variety of atmospheric properties, including composition, chemical and dynamical timescales, and large scale atmospheric motions. Critical to planning and interpreting these observations is an accurate knowledge of the planet's orbit. We report on observations of two full-transit events: 7 February 2020 as observed by the TESS spacecraft and 7--8 December 2021 as observed with a worldwide network of small telescopes. We also report new radial velocity observations which when analyzed with a coupled model to the transits greatly improve the planet's orbital ephemeris. Our new orbit solution reduces the uncertainty in the transit and eclipse timing of the JWST era from tens of minutes to a few minutes. When combined with the planned JWST observations, this new precision may be adequate to look for non-Keplerian effects in the orbit of HD80606b., Comment: Accepted for publication in AJ; in press

We present the first gri-band period-luminosity (PL) and period-Wesenheit (PW) relations for the fundamental mode anomalous Cepheids. These PL and PW relations were derived from a combined sample of five anomalous Cepheids in globular cluster M92 and the Large Magellanic Cloud, both of which have distance accurate to ~1% available from literature. Our g-band PL relation is similar to the B-band PL relation as reported in previous study. We applied our PL and PW relations to anomalous Cepheids discovered in dwarf galaxy Crater II, and found a larger but consistent distance modulus than the recent measurements based on RR Lyrae. Our calibrations of gri-band PL and PW relations, even though less precise due to small number of anomalous Cepheids, will be useful for distance measurements to dwarf galaxies., Comment: 8 pages, 1 Table and 4 Figures; AJ accepted

We present the first absolute calibration for the yellow post-asymptotic-giant-branch (PAGB) stars in the g- and r-band based on time-series observations from the Zwicky Transient Facility. These absolute magnitudes were calibrated using four yellow PAGB stars (one non-varying star and three Type II Cepheids) located in the globular clusters. We provide two calibrations of the gr-band absolute magnitudes for the yellow PAGB stars, by using an arithmetic mean and a linear regression. We demonstrate that the linear regression provides a better fit to the g-band absolute magnitudes for the yellow PAGB stars. These calibrated gr-band absolute magnitudes have a potential to be used as population II distance indicators in the era of time-domain synoptic sky surveys., Comment: 6 pages, 3 Figures and 1 Table; AJ accepted

The Vera C. Rubin Legacy Survey of Space and Time holds the potential to revolutionize time domain astrophysics, reaching completely unexplored areas of the Universe and mapping variability time scales from minutes to a decade. To prepare to maximize the potential of the Rubin LSST data for the exploration of the transient and variable Universe, one of the four pillars of Rubin LSST science, the Transient and Variable Stars Science Collaboration, one of the eight Rubin LSST Science Collaborations, has identified research areas of interest and requirements, and paths to enable them. While our roadmap is ever-evolving, this document represents a snapshot of our plans and preparatory work in the final years and months leading up to the survey's first light., Comment: 202 pages (in book format) 34 figures plus chapter heading figures (13)

We present the first gri-band period-luminosity (PL) and period-Wesenheit (PW) relations for 37 Type II Cepheids (hereafter TIIC) located in 18 globular clusters based on photometric data from the Zwicky Transient Facility. We also updated BV IJHK-band absolute magnitudes for 58 TIIC in 24 globular clusters using the latest homogeneous distances to the globular clusters. The slopes of g/r/i and B/V/I band PL relations are found to be statistically consistent when using the same sample of distance and reddening. We employed the calibration of ri-band PL/PW relations in globular clusters to estimate a distance to M31 based on a sample of ~270 TIIC from the PAndromeda project. The distance modulus to M31, obtained using calibrated ri-band PW relation, agrees well with the recent determination based on classical Cepheids. However, distance moduli derived using the calibrated r- and i-band PL relations are systematically smaller by ~0.2 mag, suggesting there are possible additional systematic error on the PL relations. Finally, we also derive the period-color (PC) relations and for the first time the period-Q-index (PQ) relations, where the Q-index is reddening-free, for our sample of TIIC. The PC relations based on (r-i) and near-infrared colors and the PQ relations are found to be relatively independent of the pulsation periods., Comment: 14 pages, 8 Figures and 5 Tables, AJ accepted

The aim of this work is to evaluate the performance of photometric metallicity [Fe/H], determined based on V-band light-curves photometrically transformed from the gr-band light-curves. We tested this by using a set of homogeneous sample of fundamental mode RR Lyrae located in the Kepler field. It was found that the color-term is necessary in such photometric transformation. We demonstrated that including the color-term the determined photometric [Fe/H] are in good agreement to the spectroscopic [Fe/H], either based on the calibrated or the transformed V-band light-curves. We also tested the impact of Blazhko RR Lyrae in determining the photometric Fe/H], and found that Blazhko RR Lyrae can give consistent photometric [Fe/H]. Finally, we derived independent gVr-band [Fe/H]-$\phi_{31}$-P relations (where $\phi_{31}$ and P are Fourier parameter and pulsation period, respectively) using our light-curves. The V-band relation is in good agreement with the most recent determination given in the literature., Comment: 12 pages, 2 Tables and 7 Figures; AJ accepted

Based on time-series observations collected from Zwicky Transient Facility (ZTF), we derived period-luminosity-metallicity (PLZ) and period-Wesenheit-metallicity (PWZ) relations for RR Lyrae located in globular clusters. We have applied various selection criteria to exclude RR Lyrae with problematic or spurious light curves. These selection criteria utilized information on the number of data points per light curve, amplitudes, colors, and residuals on the period-luminosity and/or period-Wesenheit relations. Due to blending, a number of RR Lyrae in globular clusters were found to be anomalously bright and have small amplitudes of their ZTF light curves. We used our final sample of ~750 RR Lyrae in 46 globular clusters covering a wide metallicity range (-2.36 < [Fe/H] < -0.54 dex) to derive PLZ and PWZ relations in gri bands. In addition, we have also derived the period-color-metallicity (PCZ) and for the first time, the PQZ relations where the Q-index is extinction-free by construction. We have compared our various relations to empirical and theoretical relations available in literature, and found a good agreement with most studies. Finally, we applied our derived PLZ relation to a dwarf galaxy, Crater II, and found its true distance modulus should be larger than the most recent determination., Comment: 28 pages, 7 Tables and 24 Figures. AJ accepted

We measured the parallax of the James Webb Space Telescope based on near simultaneous observations using the Lulin One-meter Telescope and the GROWTH India Telescope, separated at a distance of ~4214 km. This serves a great demonstration for the concept of parallax commonly taught in introductory astronomy courses., Comment: 2 pages and 1 figure, RNAAS in-press

Mira variables are asymptotic giant branch pulsating stars with long pulsation periods and large amplitudes in optical bands. By applying the random forest algorithm to the I-band light curves for the Miras in the Magellanic Clouds, we have classified these Miras into regular Miras and non-regular Miras. Wherein non-regular Miras exhibit a long-term variation in addition to their primary pulsation periods. Our results confirm that the period-luminosity relation for maximum light has a smaller dispersion, but only occurs on the regular oxygen-rich (O-rich) Miras, which we recommend to be applied in future distance scale work. We have also collected multi-band photometry for these Miras to perform a spectral-energy-distribution (SED) fitting with stellar and dust components, showing that a significant fraction of dust is present around the non-regular Miras. According to our results, we believe that the periodic long-term variations seen in the non-regular Miras might be due to the presence of dust.

Nova eruptions, thermonuclear explosions on the surfaces of white dwarfs (WDs), are now recognized to be among the most common shock-powered transients. We present the early discovery and rapid ultraviolet (UV), optical, and infrared (IR) temporal development of AT 2019qyl, a recent nova in NGC 300. The light curve shows a rapid rise lasting ≲1 day, reaching a peak absolute magnitude of M_V = −9.2 mag, and a very fast decline, fading by 2 mag over 3.5 days. A steep drop-off in the light curves after 71 days and the rapid decline timescale suggest a low-mass ejection from a massive WD with MWD ≳ 1.2 M_⊙. We present an unprecedented view of the early spectroscopic evolution of such an event. Three spectra prior to the peak reveal a complex, multi-component outflow giving rise to internal collisions and shocks in the ejecta of an He/N-class nova. We identify a coincident IR-variable counterpart in the extensive pre-eruption coverage of the transient location, and infer the presence of a symbiotic progenitor system with an O-rich asymptotic-giant-branch donor star, as well as evidence for an earlier UV-bright outburst in 2014. We suggest that AT 2019qyl is analogous to the subset of Galactic recurrent novae with red-giant companions such as RS Oph and other embedded nova systems like V407 Cyg. Our observations provide new evidence that internal outflow collisions likely play an important role in generating the shock-powered emission from such systems.

The TESS space telescope is collecting continuous, high-precision optical photometry of stars throughout the sky, including thousands of RR Lyrae stars. In this paper, we present results for an initial sample of 118 nearby RR Lyrae stars observed in TESS Sectors 1 and 2. We use differential-image photometry to generate light curves and analyse their mode content and modulation properties. We combine accurate light curve parameters from TESS with parallax and color information from the Gaia mission to create a comprehensive classification scheme. We build a clean sample, preserving RR Lyrae stars with unusual light curve shapes, while separating other types of pulsating stars. We find that a large fraction of RR Lyrae stars exhibit various low-amplitude modes, but the distribution of those modes is markedly different from those of the bulge stars. This suggests that differences in physical parameters have an observable effect on the excitation of extra modes, potentially offering a way to uncover the origins of these signals. However, mode identification is hindered by uncertainties when identifying the true pulsation frequencies of the extra modes. We compare mode amplitude ratios in classical double-mode stars to stars with extra modes at low amplitudes and find that they separate into two distinct groups. Finally, we find a high percentage of modulated stars among the fundamental-mode pulsators, but also find that at least 28% of them do not exhibit modulation, confirming that a significant fraction of stars lack the Blazhko effect., Comment: 35 pages, 26 figures, and a 19-page appendix. Accepted for publication in ApJS. Light curve data temporarily available: https://konkoly.hu/staff/lmolnar/tess_rrl_firstlight_all_lcs.txt

In this work, we aimed to derive the $gri$-band period-luminosity (PL) and period-luminosity-color (PLC) relations for late-type contact binaries, for the first time, located in the globular clusters, using the homogeneous light curves collected by the Zwicky Transient Factory (ZTF). We started with 79 contact binaries in 15 globular clusters, and retained 30 contact binaries in 10 globular clusters that have adequate number of data points in the ZTF light curves and unaffected by blending. Magnitudes at mean and maximum light of these contact binaries were determined using a fourth-order Fourier expansion, while extinction corrections were done using the {\tt Bayerstar2019} 3D reddening map together with adopting the homogeneous distances to their host globular clusters. After removing early-type and "anomaly" contact binaries, our derived $gri$-band PL and period-Wesenheit (PW) relations exhibit a much larger dispersion with large errors on the fitted coefficients. Nevertheless, the $gr$-band PL and PW relations based on this small sample of contact binaries in globular clusters were consistent with those based on a larger sample of nearby contact binaries. Good agreements of the PL and PW relations suggested both samples of contact binaries in the local Solar neighborhood and in the distant globular clusters can be combined and used to derive and calibrate the PL, PW and PLC relations. The final derived $gr$-band PL, PW and PLC relations were much improved than those based on the limited sample of contact binaries in the globular clusters., Comment: 26 pages, 4 Tables and 16 Figures; AJ Accepted

Nova eruptions, thermonuclear explosions on the surfaces of white dwarfs (WDs), are now recognized to be among the most common shock-powered astrophysical transients. We present the early discovery and rapid ultraviolet (UV), optical, and infrared (IR) temporal development of AT 2019qyl, a recent nova in the nearby Sculptor Group galaxy NGC 300. The light curve shows a rapid rise lasting ≲1 day, reaching a peak absolute magnitude of MV = −9.2 mag and a very fast decline, fading by 2 mag over 3.5 days. A steep dropoff in the light curves after 71 days and the rapid decline timescale suggest a low-mass ejection from a massive WD with MWD ≳ 1.2 M⊙. We present an unprecedented view of the early spectroscopic evolution of such an event. Three spectra prior to the peak reveal a complex, multicomponent outflow giving rise to internal collisions and shocks in the ejecta of an He/N-class nova. We identify a coincident IR-variable counterpart in the extensive preeruption coverage of the transient location and infer the presence of a symbiotic progenitor system with an O-rich asymptotic-giant-branch donor star, as well as evidence for an earlier UV-bright outburst in 2014. We suggest that AT 2019qyl is analogous to the subset of Galactic recurrent novae with red-giant companions such as RS Oph and other embedded nova systems like V407 Cyg. Our observations provide new evidence that internal shocks between multiple, distinct outflow components likely contribute to the generation of the shock-powered emission from such systems.

We observed the episodically active asteroid (6478) Gault in 2020 with multiple telescopes in Asia and North America and found that it is no longer active after its recent outbursts at the end of 2018 and the start of 2019. The inactivity during this apparition allowed us to measure the absolute magnitude of Gault of H ( r ) = 14.63 +/- 0.02, G ( r ) = 0.21 +/- 0.02 from our secular phase-curve observations. In addition, we were able to constrain Gault's rotation period using time-series photometric lightcurves taken over 17 hr on multiple days in 2020 August, September, and October. The photometric lightcurves have a repeating less than or similar to 0.05 mag feature suggesting that (6478) Gault has a rotation period of similar to 2.5 hr and may have a semispherical or top-like shape, much like the near-Earth asteroids Ryugu and Bennu. The rotation period of similar to 2.5 hr is near the expected critical rotation period for an asteroid with the physical properties of (6478) Gault, suggesting that its activity observed over multiple epochs is due to surface mass shedding from its fast rotation spin-up by the Yarkovsky-O'Keefe-Radzievskii-Paddack effect.

The TESS space telescope is collecting continuous, high-precision optical photometry of stars throughout the sky, including thousands of RR Lyrae stars. In this paper, we present results for an initial sample of 118 nearby RR Lyrae stars observed in TESS Sectors 1 and 2. We use differential-image photometry to generate light curves and analyse their mode content and modulation properties. We combine accurate light curve parameters from TESS with parallax and color information from the Gaia mission to create a comprehensive classification scheme. We build a clean sample, preserving RR Lyrae stars with unusual light curve shapes, while separating other types of pulsating stars. We find that a large fraction of RR Lyrae stars exhibit various low-amplitude modes, but the distribution of those modes is markedly different from those of the bulge stars. This suggests that differences in physical parameters have an observable effect on the excitation of extra modes, potentially offering a way to uncover the origins of these signals. However, mode identification is hindered by uncertainties when identifying the true pulsation frequencies of the extra modes. We compare mode amplitude ratios in classical double-mode stars to stars with extra modes at low amplitudes and find that they separate into two distinct groups. Finally, we find a high percentage of modulated stars among the fundamental-mode pulsators, but also find that at least 28% of them do not exhibit modulation, confirming that a significant fraction of stars lack the Blazhko effect., Comment: 35 pages, 26 figures, and a 19-page appendix. Accepted for publication in ApJS. Light curve data temporarily available: https://konkoly.hu/staff/lmolnar/tess_rrl_firstlight_all_lcs.txt

In this work, we aimed to derive the $gri$-band period-luminosity (PL) and period-luminosity-color (PLC) relations for late-type contact binaries, for the first time, located in the globular clusters, using the homogeneous light curves collected by the Zwicky Transient Factory (ZTF). We started with 79 contact binaries in 15 globular clusters, and retained 30 contact binaries in 10 globular clusters that have adequate number of data points in the ZTF light curves and unaffected by blending. Magnitudes at mean and maximum light of these contact binaries were determined using a fourth-order Fourier expansion, while extinction corrections were done using the {\tt Bayerstar2019} 3D reddening map together with adopting the homogeneous distances to their host globular clusters. After removing early-type and "anomaly" contact binaries, our derived $gri$-band PL and period-Wesenheit (PW) relations exhibit a much larger dispersion with large errors on the fitted coefficients. Nevertheless, the $gr$-band PL and PW relations based on this small sample of contact binaries in globular clusters were consistent with those based on a larger sample of nearby contact binaries. Good agreements of the PL and PW relations suggested both samples of contact binaries in the local Solar neighborhood and in the distant globular clusters can be combined and used to derive and calibrate the PL, PW and PLC relations. The final derived $gr$-band PL, PW and PLC relations were much improved than those based on the limited sample of contact binaries in the globular clusters., Comment: 26 pages, 4 Tables and 16 Figures; AJ Accepted

We observed Episodically Active Asteroid (6478) Gault in 2020 with multiple telescopes in Asia and North America and have found that it is no longer active after its recent outbursts at the end of 2018 and start of 2019. The inactivity during this apparation allowed us to measure the absolute magnitude of Gault of H_r = 14.63 +/- 0.02, G_r = 0.21 +/- 0.02 from our secular phasecurve observations. In addition, we were able to constrain Gault's rotation period using time-series photometric lightcurves taken over 17 hours on multiple days in 2020 August, September and October. The photometric lightcurves have a repeating $\lesssim$0.05 magnitude feature suggesting that (6478) Gault has a rotation period of ~2.5 hours and may have a semi-spherical or top-like shape, much like Near-Earth Asteroids Ryugu and Bennu. The rotation period of ~2.5 hours is near to the expected critical rotation period for an asteroid with the physical properties of (6478) Gault suggesting that its activity observed over multiple epochs is due to surface mass shedding from its fast rotation spun up by the Yarkovsky-O'Keefe-Radzievskii-Paddack effect., Comment: 23 pages, 18 figures; Accepted by ApJ Letters

Nova eruptions, thermonuclear explosions on the surfaces of white dwarfs (WDs), are now recognized to be among the most common shock-powered astrophysical transients. We present the early discovery and rapid ultraviolet (UV), optical, and infrared (IR) temporal development of AT 2019qyl, a recent nova in the nearby Sculptor Group galaxy NGC 300. The light curve shows a rapid rise lasting $\lesssim 1$ day, reaching a peak absolute magnitude of $M_V = -9.2$ mag, and a very fast decline, fading by 2 mag over 3.5 days. A steep dropoff in the light curves after 71 days and the rapid decline timescale suggest a low-mass ejection from a massive WD with $M_{\rm WD} \gtrsim 1.2~M_{\odot}$. We present an unprecedented view of the early spectroscopic evolution of such an event. Three spectra prior to the peak reveal a complex, multicomponent outflow giving rise to internal collisions and shocks in the ejecta of an He/N-class nova. We identify a coincident IR-variable counterpart in the extensive preeruption coverage of the transient location and infer the presence of a symbiotic progenitor system with an O-rich asymptotic-giant-branch donor star, as well as evidence for an earlier UV-bright outburst in 2014. We suggest that AT 2019qyl is analogous to the subset of Galactic recurrent novae with red-giant companions such as RS Oph and other embedded nova systems like V407 Cyg. Our observations provide new evidence that internal shocks between multiple, distinct outflow components likely contribute to the generation of the shock-powered emission from such systems., Comment: 32 pages, 10 figures, accepted in ApJ

HO Puppis (HO Pup) was considered as a Be-star candidate based on its gamma-Cassiopeiae-type light curve, but lacked spectroscopic confirmation. Using distance measured from Gaia Data Release 2 and the spectral-energy-distribution (SED) fit on broadband photometry, the Be-star nature of HO Pup is ruled out. Furthermore, based on the 28,700 photometric data points collected from various time-domain surveys and dedicated intensive-monitoring observations, the light curves of HO Pup closely resemble IW And-type stars (as pointed out in Kimura et al. 2020a), exhibiting characteristics such as quasi-standstill phase, brightening, and dips. The light curve of HO Pup displays various variability timescales, including brightening cycles ranging from 23 to 61 days, variations with periods between 3.9 days and 50 minutes during the quasi-standstill phase, and a semi-regular ~14-day period for the dip events. We have also collected time-series spectra (with various spectral resolutions), in which Balmer emission lines and other expected spectral lines for an IW And-type star were detected (even though some of these lines were also expected to be present for Be stars). We detect Bowen fluorescence near the brightening phase, and that can be used to discriminate between IW And-type stars and Be stars. Finally, despite only observing for four nights, the polarization variation was detected, indicating that HO Pup has significant intrinsic polarization., Comment: 24 pages, 7 tables and 12 figures; ApJ accepted

We present visible and near-infrared (NIR) photometric and spectroscopic observations of interstellar object (ISO) 2I/Borisov taken from 2019 September 10 to 2019 December 20 using the GROWTH, the Apache Point Observatory Astrophysical Research Consortium 3.5 m, and the NASA Infrared Telescope Facility 3.0 m combined with pre- and postdiscovery observations of 2I obtained by the Zwicky Transient Facility from 2019 March 17 to 2019 May 5. Comparison with imaging of distant solar system comets shows an object very similar to mildly active solar system comets with an outgassing rate of similar to 10(27)mol s(-1). The photometry, taken in filters spanning the visible and NIR range, shows a gradual brightening trend of similar to 0.03 mag day(-1)since 2019 September 10 UTC for a reddish object becoming neutral in the NIR. The light curve from recent and prediscovery data reveals a brightness trend suggesting the recent onset of significant H2O sublimation with the comet being active with super volatiles such as CO at heliocentric distances >6 au consistent with its extended morphology. Using the advanced capability to significantly reduce the scattered light from the coma enabled by high-resolution NIR images from Keck adaptive optics taken on 2019 October 4, we estimate a diameter for 2I's nucleus of less than or similar to 1.4 km. We use the size estimates of 1I/'Oumuamua and 2I/Borisov to roughly estimate the slope of the ISO size distribution, resulting in a slope of similar to 3.4 1.2, similar to solar system comets and bodies produced from collisional equilibrium.

We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo's third observing run. We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration. The GW triggers had a median localization of 4480 deg^2, median distance of 267 Mpc and false alarm rates ranging from 1.5 to 1e-25 per yr. The ZTF coverage had a median enclosed probability of 39%, median depth of 20.8mag, and median response time of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UVOIR photometric points, 64 OIR spectra, and 3 radio. We find no promising kilonova (radioactivity-powered counterpart) and we convert the upper limits to constrain the underlying kilonova luminosity function. Assuming that all kilonovae are at least as luminous as GW170817 at discovery (-16.1mag), we calculate our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than -16.6mag (extrapolated peak magnitude of GW170817) and fade at 1 mag/day (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations, the joint probability of zero detections, assuming all kilonovae are brighter than -16.6mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, <57% (<89%) of putative kilonovae could be brighter than -16.6mag assuming flat (fading) evolution, at 90% confidence. If we further account for the online terrestrial probability for each GW trigger, we find that <68% of putative kilonovae could be brighter than -16.6mag. Comparing to model grids, we find that some kilonovae must have Mej < 0.03 Msun or Xlan>1e-4 or phi>30deg to be consistent with our limits. (Abridged), Comment: Submitted to ApJ

We present a preliminary analysis of comet P/2019 LM$_4$ (Palomar) as observed by the Zwicky Transient Facility (ZTF) survey in 2019 and 2020. We find that the discovery of the comet in 2019 and the recovery in 2020 is largely attributed to two separate outbursts that are $\gtrsim2$ and $\gtrsim3.9$ mag in strength. The outbursts occurred around the end of April to early May of 2019 as well as between 2020 May 8.31 and 9.52 UTC, respectively., Comment: Submitted to RNAAS

Recently, Kundu et al (2019) reported that the globular cluster NGC 5024 (M53) possesses five extra-tidal RR Lyrae. In fact, four of them were instead known members of a nearby globular cluster NGC 5053. The status of the remaining extra-tidal RR Lyrae is controversial depending on the adopted tidal radius of NGC 5024. We have also searched for additional extra-tidal RR Lyrae within an area of $\sim8$~deg$^2$ covering both globular clusters. This includes other known RR Lyrae within the search area, as well as stars that fall within the expected range of magnitudes and colors for RR Lyrae (and yet outside the cutoff of 2/3 of the tidal radii of each globular clusters for something to be called "extra-tidal") if they were extra-tidal RR Lyrae candidates for NGC 5024 or NGC 5053. Based on the the proper motion information and their locations on the color-magnitude diagram, none of the known RR Lyrae belong to the extra-tidal RR Lyrae of either globular clusters. In the cases where the stars satisfied the magnitude and color ranges of RR Lyrae, analysis of time series data taken from the Zwicky Transient Facility do not reveal periodicities, suggesting that none of these stars are RR Lyrae. We conclude that there are no extra-tidal RR Lyrae associated with either NGC 5024 or NGC 5053 located within our search area., Comment: 17 pages, 2 Tables, 12 Figures, AJ accepted

Small planets are common around late-M dwarfs and can be detected through highly precise photometry by the transit method. Planets orbiting nearby stars are particularly important as they are often the best-suited for future follow-up studies. We present observations of three nearby M-dwarfs referred to as EIC-1, EIC-2, and EIC-3, and use them to search for transits and set limits on the presence of planets. On most nights our observations are sensitive to Earth-sized transiting planets, and photometric precision is similar to or better than TESS for faint late-M dwarfs of the same magnitude (I=15 mag). We present our photometry and transit search pipeline, which utilizes simple median detrending in combination with transit least squares based transit detection (Hippke & Heller 2019).For these targets, and transiting planets between one and two Earth radii, we achieve an average transit detection probability of 60% between periods of 0.5 and 2 days, 30% between 2 and 5 days,and 10% between 5 and 10 days. These sensitivities are conservative compared to visual searches., Comment: Accepted to AJ

We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo's third observing run (O3). We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration. The GW triggers had a median localization area of 4480 deg(2), a median distance of 267 Mpc, and false-alarm rates ranging from 1.5 to 10(-25) yr(-1). The ZTF coverage in the g and r bands had a median enclosed probability of 39%, median depth of 20.8 mag, and median time lag between merger and the start of observations of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UltraViolet/Optical/InfraRed (UVOIR) photometric points, 64 OIR spectra, and three radio images using 17 different telescopes. We find no promising kilonovae (radioactivity-powered counterparts), and we show how to convert the upper limits to constrain the underlying kilonova luminosity function. Initially, we assume that all GW triggers are bona fide astrophysical events regardless of false-alarm rate and that kilonovae accompanying BNS and NSBH mergers are drawn from a common population; later, we relax these assumptions. Assuming that all kilonovae are at least as luminous as the discovery magnitude of GW170817 (-16.1 mag), we calculate that our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than -16.6 mag (the extrapolated peak magnitude of GW170817) and fade at a rate of 1 mag day(-1) (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations based on the online classifications, the joint probability of zero detections, assuming all kilonovae are brighter than -16.6 mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, no more than <57% (<89%) of

We present new near-infrared, $JHK_s$, Period--Luminosity relations (PLRs) for RR Lyrae variables in the Messier 53 (M53 or NGC 5024) globular cluster. Multi-epoch $JHK_s$ observations, obtained with the WIRCam instrument on the 3.6-m Canada France Hawaii Telescope, are used for the first time to estimate precise mean-magnitudes for 63 RR Lyrae stars in M53 including 29 fundamental-mode (RRab) and 34 first-overtone modes (RRc) variables. The $JHK_s$-band PLRs for RR Lyrae stars are best constrained for RRab types with a minimal scatter of 22, 23, and 19 mmag, respectively. The combined sample of RR Lyrae is used to derive the $K_s$-band PLR, $K_s = -2.303 (0.063) \log P + 15.212 (0.016)$ exhibiting a $1\sigma$ dispersion of only $0.027$ mag. Theoretical Period--Luminosity--Metallicity (PLZ) relations are used to predict parallaxes for 400 Galactic RR Lyrae resulting in a median parallax zero-point offset of $-7\pm3~\mu$as in {\it Gaia} Early Data Release 3 (EDR3), which increases to $22\pm2~\mu$as if the parallax corrections are applied. We also estimate a robust distance modulus, $\mu_\textrm{M53} = 16.403 \pm 0.024$ (statistical) $\pm 0.033$ (systematic) mag, to M53 based on theoretical calibrations. Homogeneous and precise mean-magnitudes for RR Lyrae in M53 together with similar literature data for M3, M4, M5, and $\omega$ Cen are used to empirically calibrate a new RR Lyrae PLZ$_{K_s}$ relation, $K_s = -0.848 (0.007) -2.320 (0.006) \log P + 0.166 (0.011) {\rm[Fe/H]}$, anchored with {\it Gaia} EDR3 distances and theoretically predicted relations, and simultaneously estimate precise RR Lyrae based distances to these globular clusters., Comment: 18 pages, 10 figures, 4 Tables, Final Version, Accepted in The Astrophysical Journal

We present visible and mid-infrared imagery and photometry of temporary Jovian co-orbital comet P/2019 LD$_2$ taken with HST/WFC3, Spitzer/IRAC, the GROWTH telescope network, visible spectroscopy from Keck/LRIS and archival ZTF observations taken between 2019 April and 2020 August. Our observations indicate that the nucleus of LD$_2$ has a radius between 0.2-1.8 km assuming a 0.08 albedo and a coma dominated by $\sim$100$\mu$ m-scale dust ejected at $\sim$1 m/s speeds with a $\sim$1'' jet pointing in the SW direction. LD$_2$ experienced a total dust mass loss of $\sim$10$^8$ kg at a loss rate of $\sim$6 kg/s with Af$\rho$/cross-section varying between $\sim$85 cm/125 km$^2$ and $\sim$200 cm/310 km$^2$ from 2019 April 9 to 2019 Nov 8. If the increase in Af$\rho$/cross-section remained constant, it implies LD$_2$'s activity began $\sim$2018 November when within 4.8 au of the Sun, implying the onset of H$_2$O sublimation. We measure CO/CO$_2$ gas production of $\lesssim$10$^{27}$ mol/s /$\lesssim$10$^{26}$ mol/s from our 4.5 $\mu$m Spitzer observations, $g$-$r$ = 0.59$\pm$0.03, $r$-$i$ = 0.18$\pm$0.05, $i$-$z$ = 0.01$\pm$0.07 from GROWTH observations, H$_2$O gas production of $\lesssim$80 kg/s scaling from our estimated $C_2$ production of $Q_{C_2}\lesssim$7.5$\times10^{24}$ mol/s from Keck/LRIS spectroscopy. We determine that the long-term orbit of LD$_2$ is similar to Jupiter family comets having close encounters with Jupiter within $\sim$0.5 Hill radius in the last $\sim$3 y, within 0.8 Hill radius in $\sim$9 y. Additionally, 78.8$\%$ of our orbital clones are ejected from the Solar System within $1 \times 10^{6}$ years having a dynamical half-life of 3.4 $\times 10^5$ years., Comment: 28 pages, 8 figures, 2 Tables, accepted for publication in AJ

We present new near-infrared ($JHK_s$) time-series observations of RR Lyrae variables in the Messier 3 (NGC 5272) globular cluster using the WIRCam instrument at the 3.6-m Canada France Hawaii Telescope. Our observations cover a sky area of $\sim 21'\times 21'$ around the cluster center and provide an average of twenty epochs of homogeneous $JHK_s$-band photometry. New homogeneous photometry is used to estimate robust mean magnitudes for 175 fundamental-mode (RRab), 47 overtone-mode (RRc), and 11 mixed-mode (RRd) variables. Our sample of 233 RR Lyrae variables is the largest thus far obtained in a single cluster with time-resolved, multi-band near-infrared photometry. Near-infrared to optical amplitude ratios for RR Lyrae in Messier 3 exhibit a systematic increase moving from RRc to short-period ($P < 0.6$~days) and long-period ($P \gtrsim 0.6$~days) RRab variables. We derive $JHK_s$-band Period--Luminosity relations for RRab, RRc, and the combined sample of variables. Absolute calibrations based on the theoretically predicted Period--Luminosity--Metallicity relations for RR Lyrae stars yield a distance modulus, $\mu = 15.041 \pm 0.017~(\textrm{statistical}) \pm 0.036~(\textrm{systematic})$~mag, to Messier 3. When anchored to trigonometric parallaxes for nearby RR Lyrae stars from the {\it Hubble Space Telescope} and the {\it Gaia} mission, our distance estimates are consistent with those resulting from the theoretical calibrations, albeit with relatively larger systematic uncertainties., Comment: 22 pages, 14 figures, Accepted for publication in The Astronomical Journal

We present time-resolved visible spectrophotometry of minimoon 2020 CD$_3$, the second asteroid known to become temporarily captured by the Earth-Moon system's gravitational field. The spectrophotometry was taken with Keck I/LRIS between wavelengths 434 nm and 912 nm in $B$, $g$, $V$, $R$, $I$ and RG850 filters as it was leaving the Earth-Moon system on 2020 March 23 UTC. The spectrophotometry of 2020 CD$_3$ most closely resembles the spectra of V-type asteroids and some Lunar rock samples with a reddish slope of ~18$\%$/100 nm between 434 nm and 761 nm corresponding to colors of $g$-$r$ = 0.62$\pm$0.08, $r$-$i$ = 0.21 $\pm$ 0.06 and an absorption band at ~900 nm corresponding to $i$-$z$ = -0.54$\pm$0.10. Combining our measured 31.9$\pm$0.1 absolute magnitude with an albedo of 0.35 typical for V-type asteroids, we determine 2020 CD$_3$'s diameter to be ~0.9$\pm$0.1 m making it the first minimoon and one of the smallest asteroids to be spectrally studied. We use our time-series photometry to detect periodic lightcurve variations with a $<$10$^{-4}$ false alarm probability corresponding to a lightcurve period of ~573 s and a lightcurve amplitude of ~1 mag implying 2020 CD$_3$ possesses a $b/a$ axial ratio of ~2.5. In addition, we extend the observational arc of 2020 CD$_3$ to 37 days between 2020 February 15 UTC and 2020 March 23 UTC. From the improved orbital solution for 2020 CD$_3$, we estimate its likely duration of its capture to be ~2 y, and we measure the non-gravitation perturbation on its orbit due to radiation pressure with an area-to-mass ratio of 6.9$\pm$2.4$\times$10$^{-4}$ m$^2$/kg implying a density of 2.3$\pm$0.8 g/cm$^3$, broadly compatible with the densities of other meter-scale asteroids and Lunar rock. We searched for pre-discovery detections of 2020 CD$_3$ in the ZTF archive as far back as 2018 October, but were unable to locate any positive detections., Comment: 26 pages, 8 figures, 3 Tables, accepted for publication in ApJL

We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo's third observing run. We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration. The GW triggers had a median localization of 4480 deg^2, median distance of 267 Mpc and false alarm rates ranging from 1.5 to 1e-25 per yr. The ZTF coverage had a median enclosed probability of 39%, median depth of 20.8mag, and median response time of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UVOIR photometric points, 64 OIR spectra, and 3 radio. We find no promising kilonova (radioactivity-powered counterpart) and we convert the upper limits to constrain the underlying kilonova luminosity function. Assuming that all kilonovae are at least as luminous as GW170817 at discovery (-16.1mag), we calculate our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than -16.6mag (extrapolated peak magnitude of GW170817) and fade at 1 mag/day (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations, the joint probability of zero detections, assuming all kilonovae are brighter than -16.6mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, <57% (<89%) of putative kilonovae could be brighter than -16.6mag assuming flat (fading) evolution, at 90% confidence. If we further account for the online terrestrial probability for each GW trigger, we find that <68% of putative kilonovae could be brighter than -16.6mag. Comparing to model grids, we find that some kilonovae must have Mej < 0.03 Msun or Xlan>1e-4 or phi>30deg to be consistent with our limits. (Abridged), Comment: Submitted to ApJ

We present a preliminary analysis of comet P/2019 LM$_4$ (Palomar) as observed by the Zwicky Transient Facility (ZTF) survey in 2019 and 2020. We find that the discovery of the comet in 2019 and the recovery in 2020 is largely attributed to two separate outbursts that are $\gtrsim2$ and $\gtrsim3.9$ mag in strength. The outbursts occurred around the end of April to early May of 2019 as well as between 2020 May 8.31 and 9.52 UTC, respectively., Comment: Submitted to RNAAS

Recently, Kundu et al (2019) reported that the globular cluster NGC 5024 (M53) possesses five extra-tidal RR Lyrae. In fact, four of them were instead known members of a nearby globular cluster NGC 5053. The status of the remaining extra-tidal RR Lyrae is controversial depending on the adopted tidal radius of NGC 5024. We have also searched for additional extra-tidal RR Lyrae within an area of $\sim8$~deg$^2$ covering both globular clusters. This includes other known RR Lyrae within the search area, as well as stars that fall within the expected range of magnitudes and colors for RR Lyrae (and yet outside the cutoff of 2/3 of the tidal radii of each globular clusters for something to be called "extra-tidal") if they were extra-tidal RR Lyrae candidates for NGC 5024 or NGC 5053. Based on the the proper motion information and their locations on the color-magnitude diagram, none of the known RR Lyrae belong to the extra-tidal RR Lyrae of either globular clusters. In the cases where the stars satisfied the magnitude and color ranges of RR Lyrae, analysis of time series data taken from the Zwicky Transient Facility do not reveal periodicities, suggesting that none of these stars are RR Lyrae. We conclude that there are no extra-tidal RR Lyrae associated with either NGC 5024 or NGC 5053 located within our search area., Comment: 17 pages, 2 Tables, 12 Figures, AJ accepted

Small planets are common around late-M dwarfs and can be detected through highly precise photometry by the transit method. Planets orbiting nearby stars are particularly important as they are often the best-suited for future follow-up studies. We present observations of three nearby M-dwarfs referred to as EIC-1, EIC-2, and EIC-3, and use them to search for transits and set limits on the presence of planets. On most nights our observations are sensitive to Earth-sized transiting planets, and photometric precision is similar to or better than TESS for faint late-M dwarfs of the same magnitude (I=15 mag). We present our photometry and transit search pipeline, which utilizes simple median detrending in combination with transit least squares based transit detection (Hippke & Heller 2019).For these targets, and transiting planets between one and two Earth radii, we achieve an average transit detection probability of 60% between periods of 0.5 and 2 days, 30% between 2 and 5 days,and 10% between 5 and 10 days. These sensitivities are conservative compared to visual searches., Comment: Accepted to AJ

We present detailed observations of ZTF18abukavn (SN2018gep), discovered in high-cadence data from the Zwicky Transient Facility as a rapidly rising (1.3 mag/hr) and luminous (M_(g,peak) = −20 mag) transient. It is spectroscopically classified as a broad-lined stripped-envelope supernova (Ic-BL SN). The rapid rise to peak bolometric luminosity and blue colors at peak (t_(rise)∼0.5-3 days, L_(bol)≳3×10^(44) erg sec^(−1), g−r = −0.3) resemble the high-redshift Ic-BL iPTF16asu, as well as several other unclassified fast transients. The early discovery of SN2018gep (within an hour of shock breakout) enabled an intensive spectroscopic campaign, including the highest-temperature (T_(eff) ≳ 40,000K) spectra of a stripped-envelope SN. A retrospective search revealed luminous (M_g ∼ M_r ≈ −14mag) emission in the days to weeks before explosion, the first definitive detection of precursor emission for a Ic-BL. We find a limit on the isotropic gamma-ray energy release E_(γ,iso) < 4.9×10^(48) erg, a limit on X-ray emission L_X < 10^(40) erg sec^(−1), and a limit on radio emission νL_ν ≲ 10^(37) erg sec^(−1). Taken together, we find that the data are best explained by shock breakout in a massive shell of dense circumstellar material (0.02 M⊙) at large radii (3×10^(14)cm) that was ejected in eruptive pre-explosion mass-loss episodes.