Search

We present a detailed study of the MASTER OT J030227.28+191754.5 outburst in 2021-2022, reaching an amplitude of 10.2 mag and a duration of 60 d. The detections of (1) the double-peaked optical emission lines, and (2) the early and ordinary superhumps, established that MASTER OT J030227.28+191754.5 is an extremely energetic WZ Sge-type dwarf nova (DN). Based on the superhump observations, we obtained its orbital period and mass ratio as 0.05986(1) d and 0.063(1), respectively. These are within a typical range of low-mass-ratio DNe. According to the binary parameters derived based on the thermal-tidal instability model, our analyses showed that (1) the standard disk model requires an accretion rate $\simeq$ 10$^{20}$ g s$^{-1}$ to explain its peak optical luminosity and (2) large mass was stored in the disk at the outburst onset. These cannot be explained solely by the impact of its massive ($\gtrsim$ 1.15 M$_\odot$) primary white dwarf implied by Kimura et al. (2023). Instead, we propose that the probable origin of this enormously energetic DN outburst is the even lower quiescence viscosity than other WZ Sge-type DNe. This discussion is qualitatively valid for most possible binary parameter spaces unless the inclination is low ($\lesssim 40^\circ$) enough for the disk to be bright explaining the outburst amplitude. Such low inclinations, however, would not allow detectable amplitude of early superhumps in the current thermal-tidal instability model. The optical spectra at outburst maximum showed the strong emission lines of Balmer, He I, and He II series whose core is narrower than $\sim 800$ km s$^{-1}$. Considering its binary parameters, a Keplerian disk cannot explain this narrow component, but the presumable origin is disk winds., Comment: 23 pages, 13 figures, 2 tables. Accepted by PASJ. Part of the online supplemental information is included

A small number of pulsars are known to emit giant pulses, single pulses much brighter than average. Among these is PSR J0534+2200, also known as the Crab pulsar, a young pulsar with high giant pulse rates. Long-term monitoring of the Crab pulsar presents an excellent opportunity to perform statistical studies of its giant pulses and the processes affecting them, potentially providing insight into the behavior of other neutron stars that emit bright single pulses. Here, we present an analysis of a set of 24,985 Crab giant pulses obtained from 88 hours of daily observations at a center frequency of 1.55 GHz by the 20-meter telescope at the Green Bank Observatory, spread over 461 days. We study the effects of refractive scintillation at higher frequencies than previous studies and compare methods of correcting for this effect. We also search for deterministic patterns seen in other single-pulse sources, possible periodicities seen in several rotating radio transients and fast radio bursts, and clustering of giant pulses like that seen in the repeating fast radio burst FRB121102., Comment: 21 pages, 14 figures, accepted to ApJ

We report on findings from scintillation analyses using high-cadence observations of eight canonical pulsars with observing baselines ranging from one to three years. We obtain scintillation bandwidth and timescale measurements for all pulsars in our survey, scintillation arc curvature measurements for four, and detect multiple arcs for two. We find evidence of a previously undocumented scattering screen along the line of sight (LOS) to PSR J1645$-$0317, as well as evidence that a scattering screen along the LOS to PSR J2313$+$4253 may reside somewhere within the Milky Way's Orion-Cygnus arm. We report evidence of a significant change in the scintillation pattern in PSR J2022$+$5154 from the previous two decades of literature, wherein both the scintillation bandwidth and timescale decreased by an order of magnitude relative to earlier observations at the same frequencies, potentially as a result of a different screen dominating the observed scattering. By augmenting the results of previous studies, we find general agreement with estimations of scattering delays from pulsar observations and predictions by the NE2001 electron density model but not for the newest data we have collected, providing some evidence of changes in the ISM along various LOSs over the timespans considered. In a similar manner, we find additional evidence of a correlation between a pulsar's dispersion measure and the overall variability of its scattering delays over time. The plethora of interesting science obtained through these observations demonstrates the capabilities of the Green Bank Observatory's 20m telescope to contribute to pulsar-based studies of the interstellar medium., Comment: Accepted to ApJ

We report the study of a huge optical intraday flare on November 12, 2021, at 2 am UT, in the blazar OJ287. In the binary black hole model it is associated with an impact of the secondary black hole on the accretion disk of the primary. Our multifrequency observing campaign was set up to search for such a signature of the impact, based on a prediction made eight years earlier. The first I-band results of the flare have already been reported by \cite{2024ApJ...960...11K}. Here we combine these data with our monitoring in the R-band. There is a big change in the R-I spectral index by $1.0\pm0.1$ between the normal background and the flare, suggesting a new component of radiation. The polarization variation during the rise of the flare suggests the same. The limits on the source size place it most reasonably in the jet of the secondary black hole. We then ask why we have not seen this phenomenon before. We show that OJ287 was never before observed with sufficient sensitivity on the night when the flare should have happened according to the binary model. We also study the probability that this flare is just an oversized example of intraday variability, using the Krakow-dataset of intense monitoring between 2015 and 2023. We find that the occurrence of a flare of this size and rapidity is unlikely. In the Appendix, we give the full orbit-linked historical light curve of OJ287 as well as the dense monitoring sample of Krakow., Comment: to Appear in Astrophysical Journal Letters

We present photometric and spectroscopic observations of SN 2023fyq, a type Ibn supernova in the nearby galaxy NGC 4388 (D$\simeq$18~Mpc). In addition, we trace long-standing precursor emission at the position of SN 2023fyq using data from DLT40, ATLAS, ZTF, ASAS-SN, Swift, and amateur astronomer Koichi Itagaki. Precursor activity is observed up to nearly three years before the supernova explosion, with a relatively rapid rise in the final 100 days. The double-peaked post-explosion light curve reaches a luminosity of $\sim10^{43}~\rm erg\,s^{-1}$. The strong intermediate-width He lines observed in the nebular spectrum of SN 2023fyq imply the interaction is still active at late phases. We found that the precursor activity in SN 2023fyq is best explained by the mass transfer in a binary system involving a low-mass He star and a compact companion. An equatorial disk is likely formed in this process ($\sim$0.6$\rm M_{\odot}$), and the interaction of SN ejecta with this disk powers the main peak of the supernova. The early SN light curve reveals the presence of dense extended material ($\sim$0.3$\rm M_{\odot}$) at $\sim$3000$\rm R_{\odot}$ ejected weeks before the SN explosion, likely due to final-stage core silicon burning or runaway mass transfer resulting from binary orbital shrinking, leading to rapid rising precursor emission within $\sim$30 days prior to explosion. The final explosion could be triggered either by the core-collapse of the He star or by the merger of the He star with a compact object. SN 2023fyq, along with SN 2018gjx and SN 2015G, forms a unique class of Type Ibn SNe which originate in binary systems and are likely to exhibit detectable long-lasting pre-explosion outbursts with magnitudes ranging from $-$10 to $-$13., Comment: submitted to ApJ

EX Lupi, a low-mass young stellar object, went into an accretion-driven outburst in March of 2022. The outburst caused a sudden phase change of ~ 112$^{\circ}$ $\pm$ 5$^{\circ}$ in periodically oscillating multiband lightcurves. Our high resolution spectra obtained with HRS on SALT also revealed a consistent phase change in the periodically varying radial velocities, along with an increase in the radial velocity amplitude of various emission lines. The phase change and increase of radial velocity amplitude morphologically translates to a change in the azimuthal and latitudinal location of the accretion hotspot over the stellar surface, which indicates a reconfiguration of the accretion funnel geometry. Our 3D MHD simulations reproduce the phase change for EX Lupi. To explain the observations we explored the possibility of forward shifting of the dipolar accretion funnel as well as the possibility of an emergence of a new accretion funnel. During the outburst, we also found evidence of the hotspot's morphology extending azimuthally, asymmetrically with a leading hot edge and cold tail along the stellar rotation. Our high cadence photometry showed that the accretion flow has clumps. We also detected possible clumpy accretion events in the HRS spectra, that showed episodically highly blue-shifted wings in the Ca II IRT and Balmer H lines., Comment: Accepted for publication in The Astrophysical Journal

We present comprehensive optical observations of SN~2021gmj, a Type II supernova (SN~II) discovered within a day of explosion by the Distance Less Than 40~Mpc (DLT40) survey. Follow-up observations show that SN~2021gmj is a low-luminosity SN~II (LL~SN~II), with a peak magnitude $M_V = -15.45$ and Fe~II velocity of $\sim 1800 \ \mathrm{km} \ \mathrm{s}^{-1}$ at 50 days past explosion. Using the expanding photosphere method, we derive a distance of $17.8^{+0.6}_{-0.4}$~Mpc. From the tail of the light curve we obtain a radioactive nickel mass of $0.014 \pm 0.001$ M$_{\odot}$. The presence of circumstellar material (CSM) is suggested by the early-time light curve, early spectra, and high-velocity H$\alpha$ in absorption. Analytical shock-cooling models of the light curve cannot reproduce the fast rise, supporting the idea that the early-time emission is partially powered by the interaction of the SN ejecta and CSM. The inferred low CSM mass of 0.025 M$_{\odot}$ in our hydrodynamic-modeling light curve analysis is also consistent with our spectroscopy. We observe a broad feature near 4600 \AA, which may be high-ionization lines of C, N, or/and He~II. This feature is reproduced by radiation-hydrodynamic simulations of red supergiants with extended atmospheres. Several LL~SNe~II show similar spectral features, implying that high-density material around the progenitor may be common among them., Comment: Accepted version at ApJ

We present optical, infrared, ultraviolet, and radio observations of SN 2022xkq, an underluminous fast-declining type Ia supernova (SN Ia) in NGC 1784 ($\mathrm{D}\approx31$ Mpc), from $<1$ to 180 days after explosion. The high-cadence observations of SN 2022xkq, a photometrically transitional and spectroscopically 91bg-like SN Ia, cover the first days and weeks following explosion which are critical to distinguishing between explosion scenarios. The early light curve of SN 2022xkq has a red early color and exhibits a flux excess which is more prominent in redder bands; this is the first time such a feature has been seen in a transitional/91bg-like SN Ia. We also present 92 optical and 19 near-infrared (NIR) spectra, beginning 0.4 days after explosion in the optical and 2.6 days after explosion in the NIR. SN 2022xkq exhibits a long-lived C I 1.0693 $\mu$m feature which persists until 5 days post-maximum. We also detect C II $\lambda$6580 in the pre-maximum optical spectra. These lines are evidence for unburnt carbon that is difficult to reconcile with the double detonation of a sub-Chandrasekhar mass white dwarf. No existing explosion model can fully explain the photometric and spectroscopic dataset of SN 2022xkq, but the considerable breadth of the observations is ideal for furthering our understanding of the processes which produce faint SNe Ia., Comment: 38 pages, 16 figures, accepted for publication in ApJ, the figure 15 input models and synthetic spectra are now available at https://zenodo.org/record/8379254

We present optical and near-infrared observations of SN~2022crv, a stripped envelope supernova in NGC~3054, discovered within 12 hrs of explosion by the Distance Less Than 40 Mpc Survey. We suggest SN~2022crv is a transitional object on the continuum between SNe Ib and SNe IIb. A high-velocity hydrogen feature ($\sim$$-$20,000 -- $-$16,000 $\rm km\,s^{-1}$) was conspicuous in SN~2022crv at early phases, and then quickly disappeared around maximum light. By comparing with hydrodynamic modeling, we find that a hydrogen envelope of $\sim 10^{-3}$ \msun{} can reproduce the behaviour of the hydrogen feature observed in SN~2022crv. The early light curve of SN~2022crv did not show envelope cooling emission, implying that SN~2022crv had a compact progenitor with extremely low amount of hydrogen. The analysis of the nebular spectra shows that SN~2022crv is consistent with the explosion of a He star with a final mass of $\sim$4.5 -- 5.6 \msun{} that has evolved from a $\sim$16 -- 22 \msun{} zero-age main sequence star in a binary system with about 1.0 -- 1.7 \msun{} of oxygen finally synthesized in the core. The high metallicity at the supernova site indicates that the progenitor experienced a strong stellar wind mass loss. In order to retain a small amount of residual hydrogen at such a high metallicity, the initial orbital separation of the binary system is likely larger than $\sim$1000~$\rm R_{\odot}$. The near-infrared spectra of SN~2022crv show a unique absorption feature on the blue side of He I line at $\sim$1.005~$\mu$m. This is the first time that such a feature has been observed in a Type Ib/IIb, and could be due to \ion{Sr}{2}. Further detailed modelling on SN~2022crv can shed light on the progenitor and the origin of the mysterious absorption feature in the near infrared., Comment: accepted for publication in ApJ

Intense outbursts in blazars are among the most extreme phenomena seen in extragalactic objects. Studying these events can offer important information about the energetic physical processes taking place within the innermost regions of blazars, which are beyond the resolution of current instruments. This work presents some of the largest and most rapid flares detected in the optical band from the sources 3C 279, OJ 49, S4 0954+658, Ton 599, and PG 1553+113, which are mostly TeV blazars. The source flux increased by nearly ten times within a few weeks, indicating the violent nature of these events. Such energetic events might originate from magnetohydrodynamical instabilities near the base of the jets, triggered by processes modulated by the magnetic field of the accretion disc. We explain the emergence of flares owing to the injection of high-energy particles by the shock wave passing along the relativistic jets. Alternatively, the flares may have also arisen due to geometrical effects related to the jets. We discuss both source-intrinsic and source-extrinsic scenarios as possible explanations for the observed large amplitude flux changes., Comment: 8 pages, 2 figures, 38th International Cosmic Ray Conference (ICRC2023) proceedings

The highly variable blazar OJ~287 is commonly discussed as an example of a binary black hole system. The 130 year long optical light curve is well explained by a model where the central body is a massive black hole of 18.35$\times$10$^9$ solar mass that supports a thin accretion disc. The secondary black hole of 0.15$\times$10$^9$ solar mass impacts the disc twice during its 12 year orbit, and causes observable flares. Recently, it has been argued that an accretion disc with a typical AGN accretion rate and above mentioned central body mass should be at least six magnitudes brighter than OJ~287's host galaxy and would therefore be observationally excluded. Based on the observations of OJ~287's radio jet, detailed in Marscher and Jorstad (2011), and up-to-date accretion disc models of Azadi et al. (2022), we show that the V-band magnitude of the accretion disc is unlikely to exceed the host galaxy brightness by more than one magnitude, and could well be fainter than the host. This is because accretion power is necessary to launch the jet as well as to create electromagnetic radiation, distributed across many wavelengths, and not concentrated especially on the optical V-band. Further, we note that the claimed V-band concentration of accretion power leads to serious problems while interpreting observations of other Active Galactic Nuclei. Therefore, we infer that the mass of the primary black hole and its accretion rate do not need to be smaller than what is determined in the standard model for OJ~287.

We present a summary of the results of the OJ 287 observational campaign, which was carried out during the 2021/2022 observational season. This season is special in the binary model because the major axis of the precessing binary happens to lie almost exactly in the plane of the accretion disc of the primary. This leads to pairs of almost identical impacts between the secondary black hole and the accretion disk in 2005 and 2022. In 2005, a special flare called "blue flash" was observed 35 days after the disk impact, which should have also been verifiable in 2022. We did observe a similar flash and were able to obtain more details of its properties. We describe this in the framework of expanding cloud models. In addition, we were able to identify the flare arising exactly at the time of the disc crossing from its photo-polarimetric and gamma-ray properties. This is an important identification, as it directly confirms the orbit model. Moreover, we saw a huge flare that lasted only one day. We may understand this as the lighting up of the jet of the secondary black hole when its Roche lobe is suddenly flooded by the gas from the primary disk. Therefore, this may be the first time we directly observed the secondary black hole in the OJ 287 binary system., Comment: 24 pages, 13 figures, published in Galaxies as part of the Special Issue Distant Glowing Objects: Quest for Quasars

We perform a comprehensive search for optical precursor emission at the position of SN~2023ixf using data from the DLT40, ZTF and ATLAS surveys. By comparing the current data set with precursor outburst hydrodynamical model light curves, we find that the probability of a significant outburst within five years of explosion is low, and the circumstellar material (CSM) ejected during any possible precursor outburst is likely smaller than $\sim$0.015\msun. By comparing to a set of toy models, we find that, if there was a precursor outburst, the duration must have been shorter than $\sim$100 days for a typical brightness of $M_{r}\simeq-9$ mag or shorter than 200 days for $M_{r}\simeq-8$ mag; brighter, longer outbursts would have been discovered. Precursor activity like that observed in the normal type II SN~2020tlf ($M_{r}\simeq-11.5$) can be excluded in SN~2023ixf. If the dense CSM inferred by early flash spectroscopy and other studies is related to one or more precursor outbursts, then our observations indicate that any such outburst would have to be faint and only last for days to months, or it occurred more than five years prior to the explosion. Alternatively, any dense, confined CSM may not be due to eruptive mass loss from a single red supergiant (RSG) progenitor. Taken together, the results of SN~2023ixf and SN~2020tlf indicate that there may be more than one physical mechanism behind the dense CSM inferred around some normal type II SNe., Comment: 10 pages, 4 figures

We present the optical spectroscopic evolution of SN~2023ixf seen in sub-night cadence spectra from 1.18 to 14 days after explosion. We identify high-ionization emission features, signatures of interaction with material surrounding the progenitor star, that fade over the first 7 days, with rapid evolution between spectra observed within the same night. We compare the emission lines present and their relative strength to those of other supernovae with early interaction, finding a close match to SN~2020pni and SN~2017ahn in the first spectrum and SN~2014G at later epochs. To physically interpret our observations we compare them to CMFGEN models with confined, dense circumstellar material around a red supergiant progenitor from the literature. We find that very few models reproduce the blended \NC{} emission lines observed in the first few spectra and their rapid disappearance thereafter, making this a unique diagnostic. From the best models, we find a mass-loss rate of $10^{-3}-10^{-2}$ \mlunit{}, which far exceeds the mass-loss rate for any steady wind, especially for a red supergiant in the initial mass range of the detected progenitor. These mass-loss rates are, however, similar to rates inferred for other supernovae with early circumstellar interaction. Using the phase when the narrow emission features disappear, we calculate an outer dense radius of circumstellar material $R_\mathrm{CSM, out}\sim5\times10^{14}~\mathrm{cm}$ and a mean circumstellar material density of $\rho=5.6\times10^{-14}~\mathrm{g\,cm^{-3}}$. This is consistent with the lower limit on the outer radius of the circumstellar material we calculate from the peak \Halpha{} emission flux, $R_\text{CSM, out}\gtrsim9\times10^{13}~\mathrm{cm}$., Comment: Published in ApJL

We present the densely sampled early light curve of the Type II supernova (SN) 2023ixf, first observed within hours of explosion in the nearby Pinwheel Galaxy (Messier 101; 6.7 Mpc). Comparing these data to recently updated models of shock-cooling emission, we find that the progenitor likely had a radius of $410 \pm 10\ R_\odot$. Our estimate is model dependent but consistent with a red supergiant. These models provide a good fit to the data starting about 1 day after the explosion, despite the fact that the classification spectrum shows signatures of circumstellar material around SN 2023ixf during that time. Photometry during the first day after the explosion, provided almost entirely by amateur astronomers, does not agree with the shock-cooling models or a simple power-law rise fit to data after 1 day. We consider the possible causes of this discrepancy, including precursor activity from the progenitor star, circumstellar interaction, and emission from the shock before or after it breaks out of the stellar surface. The very low luminosity ($-11\mathrm{\ mag} > M > -14\mathrm{\ mag}$) and short duration of the initial excess lead us to prefer a scenario related to prolonged emission from the SN shock traveling through the progenitor system., Comment: updated to match accepted version

We present very early photometric and spectroscopic observations of the Type Ia supernova (SN Ia) 2023bee, starting about 8 hr after the explosion, which reveal a strong excess in the optical and nearest UV (U and UVW1) bands during the first several days of explosion. This data set allows us to probe the nature of the binary companion of the exploding white dwarf and the conditions leading to its ignition. We find a good match to the Kasen model in which a main-sequence companion star stings the ejecta with a shock as they buzz past. Models of double detonations, shells of radioactive nickel near the surface, interaction with circumstellar material, and pulsational delayed detonations do not provide good matches to our light curves. We also observe signatures of unburned material, in the form of carbon absorption, in our earliest spectra. Our radio nondetections place a limit on the mass-loss rate from the putative companion that rules out a red giant but allows a main-sequence star. We discuss our results in the context of other similar SNe Ia in the literature., Comment: updated to match accepted version

We present five far- and near-ultraviolet spectra of the Type II plateau supernova, SN 2022acko, obtained 5, 6, 7, 19, and 21 days after explosion, all observed with the Hubble Space Telescope/Space Telescope Imaging Spectrograph. The first three epochs are earlier than any Type II plateau supernova has been observed in the far-ultraviolet revealing unprecedented characteristics. These three spectra are dominated by strong lines, primarily from metals, which contrasts with the relatively featureless early optical spectra. The flux decreases over the initial time series as the ejecta cools and line-blanketing takes effect. We model this unique dataset with the non-local thermodynamic equilibrium radiation transport code CMFGEN, finding a good match to the explosion of a low mass red supergiant with energy Ekin = 6 x 10^50 erg. With these models we identify, for the first time, the ions that dominate the early UV spectra. We also present optical photometry and spectroscopy, showing that SN 2022acko has a peak absolute magnitude of V = -15.4 mag and plateau length of ~115d. The spectra closely resemble those of SN 2005cs and SN 2012A. Using the combined optical and UV spectra, we report the fraction of flux redwards of the uvw2, U, B, and V filters on days 5, 7, and 19. We also create a spectral time-series of Type II supernovae in the ultraviolet, demonstrating the rapid decline of UV flux over the first few weeks of evolution. Future observations of Type II supernovae will continue to explore the diversity seen in the limited set of high-quality UV spectra., Comment: Published in ApJL

This study provides the results of simultaneous multicolor observations for the first Visorsat (STARLINK-1436) and the ordinary Starlink satellite, STARLINK-1113 in the $U$, $B$, $V$, $g'$, $r$, $i$, $R_{\rm C}$, $I_{\rm C}$, $z$, $J$, $H$, and $K_s$ bands to quantitatively investigate the extent to which Visorsat reduces its reflected light. Our results are as follows: (1) in most cases, Virorsat is fainter than STARLINK-1113, and the sunshade on Visorsat, therefore, contributes to the reduction of the reflected sunlight; (2) the magnitude at 550 km altitude (normalized magnitude) of both satellites often reaches the naked-eye limiting magnitude ($<$ 6.0); (3) from a blackbody radiation model of the reflected flux, the peak of the reflected components of both satellites is around the $z$ band; and (4) the albedo of the near infrared range is larger than that of the optical range. Under the assumption that Visorsat and STARLINK-1113 have the same reflectivity, we estimate the covering factor, $C_{\rm f}$, of the sunshade on Visorsat, using the blackbody radiation model: the covering factor ranges from $0.18 \leq C_{\rm f} \leq 0.92$. From the multivariable analysis of the solar phase angle (Sun-target-observer), the normalized magnitude, and the covering factor, the phase angle versus covering factor distribution presents a moderate anti-correlation between them, suggesting that the magnitudes of Visorsat depend not only on the phase angle but also on the orientation of the sunshade along our line of sight. However, the impact on astronomical observations from Visorsat-designed satellites remains serious. Thus, new countermeasures are necessary for the Starlink satellites to further reduce reflected sunlight., Comment: 31 pages, 9 figures, published in PASJ

Funded by a $3M Department of Defense (DoD) National Defense Education Program (NDEP) award, we are developing and deploying on a national scale a follow-up curriculum to "Our Place In Space!", or OPIS!, in which approx. 3,500 survey-level astronomy students are using our global network of "Skynet" robotic telescopes each year. The goal of this new curriculum, called "Astrophotography of the Multi-Wavelength Universe!", or MWU!, is to boost the number of these students who choose STEM majors. During Y1, our participating educators have developed MWU!'s (now renumbered) 2nd and 4th modules, and are in the process of developing its 3rd and 7th modules (out of 7). Solid progress has also been made on the software front, (1) where we have developed new graphing/analysis/modeling interfaces in support of Modules 2 and 4, and in response to feedback from the participating educators; and (2) where we are in the process of developing and adding astrophotography capabilities to Afterglow Access (AgA), our student-level, web-based, image processing and analysis application, in support of Modules 1 - 3 and 5 - 7. On the hardware front, development of our first four signal-processing units proceeds on schedule; these are key to Skynet's integration of a global network of radio telescopes, capable of exploring the invisible universe. Preparations have also been made on the evaluation and accessibility fronts, for when the first MWU! modules are deployed in Spring 2023., Comment: 10 pages, 9 figures, ASP2022 Conference Proceedings, December 2022

The bright blazar OJ~287 routinely parades high brightness bremsstrahlung flares, which are explained as being a result of a secondary supermassive black hole (SMBH) impacting the accretion disc of a more massive primary SMBH in a binary system. The accretion disc is not rigid but rather bends in a calculable way due to the tidal influence of the secondary. Below we refer to this phenomenon as a variable disc level. We begin by showing that these flares occur at times predicted by a simple analytical formula, based on general relativity inspired modified Kepler equation, which explains impact flares since 1888. The 2022 impact flare, namely flare number 26, is rather peculiar as it breaks the typical pattern of two impact flares per 12-year cycle. This is the third bremsstrahlung flare of the current cycle that follows the already observed 2015 and 2019 impact flares from OJ~287. It turns out that the arrival epoch of flare number 26 is sensitive to the level of primary SMBH's accretion disc relative to its mean level in our model. We incorporate these tidally induced changes in the level of the accretion disc to infer that the thermal flare should have occurred during July-August 2022, when it was not possible to observe it from the Earth. Thereafter, we explore possible observational evidence for certain pre-flare activity by employing spectral and polarimetric data from our campaigns in 2004/05 and 2021/22. We point out theoretical and observational implications of two observed mini-flares during January-February 2022., Comment: 29 pages, 6 figures, 1 table. arXiv admin note: text overlap with arXiv:2209.08360

Flaring episodes in blazars represent one of the most violent processes observed in extra-galactic objects. Studies of such events shed light on the energetics of the physical processes occurring in the innermost regions of blazars, which cannot otherwise be resolved by any current instruments. In this work, we present some of the largest and most rapid flares captured in the optical band in the blazars 3C 279, OJ 49, S4 0954+658, TXS 1156+295 and PG 1553+113. The source flux was observed to increase by nearly ten times within a timescale of a few weeks. We applied several methods of time series analysis and symmetry analysis. Moreover, we also performed searches for periodicity in the light curves of 3C 279, OJ 49 and PG 1553+113 using the Lomb-Scargle method and found plausible indications of quasi-periodic oscillations (QPOs). In particular, the 33- and 22-day periods found in 3C 279, i.e. a 3:2 ratio, are intriguing. These violent events might originate from magnetohydrodynamical instabilities near the base of the jets, triggered by processes modulated by the magnetic field of the accretion disc. We present a qualitative treatment as the possible explanation for the observed large amplitude flux changes in both the source-intrinsic and source-extrinsic scenarios., Comment: 11 pages, 6 figures, MNRAS accepted

Built in 2004, the Skynet robotic telescope network originally consisted of six 0.4 m telescopes located at the Cerro-Tololo Inter-American Observatory in the Chilean Andes. The network was designed to carry out simultaneous multi-wavelength observations of gamma-ray bursts (GRBs) when they are only tens of seconds old. To date, the network has been expanded to ~20 telescopes, including a 20 m radio telescope, that span four continents and five countries. The Campaign Manager (CM) is a new observing mode that has been developed for Skynet. Available to all Skynet observers, the CM semi-autonomously and indefinitely scales and schedules exposures on the observer's behalf while allowing for modification to scaling parameters in real time. The CM is useful for follow up to various transient phenomena including gravitational-wave events, GRB localizations, young supernovae, and eventually, sufficiently bright Argus Optical Array and Large Synoptic Survey Telescope events., Comment: 7 pages, 4 figures

We characterize the optical variability properties of eight lobe-dominated radio quasars (QSOs): B2 0709$+$37, FBQS J095206.3$+$235245, PG 1004$+$130, [HB89] 1156$+$631, [HB89] 1425$+$267, [HB89] 1503$+$691, [HB89] 1721$+$343, 4C $+$74.26, systematically monitored for a duration of 13 years since 2009. The quasars are radio-loud objects with extended radio lobes that indicate their orientation close to the sky plane. Five of the eight QSOs are classified as giant radio quasars. All quasars showed variability during our monitoring, with magnitude variations between 0.3 and 1 mag for the least variable and the most variable QSO, respectively. We performed both structure function (SF) analysis and power spectrum density (PSD) analysis for the variability characterization and search for characteristic timescales and periodicities. As a result of our analysis, we obtained relatively steep SF slopes ($\alpha$ ranging from 0.49 to 0.75) that are consistent with the derived PSD slopes ($\sim$2--3). All the PSDs show a good fit to single power law forms, indicating a red-noise character of variability between $\sim$13 years and weeks timescales. We did not measure reliable characteristic timescales of variability from the SF analysis which indicates that the duration of the gathered data is too short to reveal them. The absence of bends in the PSDs (change of slope from $\geq$1 to $\sim$0) on longer timescales indicates that optical variations are most likely caused by thermal instabilities in the accretion disk., Comment: Accepted for publication in ApJS; 17 pages, 5 figures, 5 tables

We present time-series CCD photometry in the $BVRI$ passbands of the recently identified symbiotic nova V1835 Aquilae (NSV 11749) over an interval of 5.1 years with 7-14 day cadence, observed during its quiescence. We find slow light variations with a range of $\sim$0.9 mag in $V$ and $\sim$0.3 mag in $I$. Analysis of these data show strong periodicity at $419 \pm 10$ days, which we interpret to be the system's orbital period. A dip in the otherwise-sinusoidal phased light curve suggests a weak ellipsoidal effect due to tidal distortion of the giant star, which in turn opens the possibility that V1835 Aql transfers some of its mass to the hot component via Roche lobe overflow rather than via a stellar wind. We also find evidence that V1835 Aql is an S-type symbiotic star, relatively free of circumstellar dust, and include it among the nuclear burning group of symbiotics. Finally, we provide photometry, periods, and light curve classifications for 22 variable stars in the field around V1835 Aql, about half of which are newly identified., Comment: Main Paper: 28 pages, 5 figures, 5 tables. Supplement: 15 pages, 4 figures, 1 table. To be published in Publications of the Astronomical Society of the Pacific

The bright blazar OJ~287 routinely parades high brightness bremsstrahlung flares which are explained as being a result of a secondary supermassive black hole (SMBH) impacting the accretion disk of a primary SMBH in a binary system. We begin by showing that these flares occur at times predicted by a simple analytical formula, based on the Kepler equation, which explains flares since 1888. The next impact flare, namely the flare number 26, is rather peculiar as it breaks the typical pattern of two impact flares per 12 year cycle. This will be the third bremsstrahlung flare of the current cycle that follows the already observed 2015 and 2019 impact flares from OJ~287. Unfortunately, astrophysical considerations make it difficult to predict the exact arrival epoch of the flare number 26. In the second part of the paper, we describe our recent OJ~287 observations. They show that the pre-flare light curve of flare number 22, observed in 2005, exhibits similar activity as the pre-flare light curve in 2022, preceding the expected flare number 26 in our model. We argue that the pre-flare activity most likely arises in the primary jet whose activity is modulated by the transit of the secondary SMBH through the accretion disk of the primary. Observing the next impact flare of OJ~287 in October 2022 will substantiate the theory of disk impacts in binary black hole systems., Comment: 16 pages, 2 figures

We present photometric and spectroscopic data of SN 2018lab, a low luminosity type IIP supernova (LLSN) with a V-band peak luminosity of $-15.1\pm0.1$ mag. SN 2018lab was discovered by the Distance Less Than 40 Mpc (DLT40) SNe survey only 0.73 days post-explosion, as determined by observations from the Transiting Exoplanet Survey Satellite (TESS). TESS observations of SN 2018lab yield a densely sampled, fast-rising, early time light curve likely powered by circumstellar medium (CSM) interaction. The blue-shifted, broadened flash feature in the earliest spectra ($<$2 days) of SN 2018lab provide further evidence for ejecta-CSM interaction. The early emission features in the spectra of SN 2018lab are well described by models of a red supergiant progenitor with an extended envelope and close-in CSM. As one of the few LLSNe with observed flash features, SN 2018lab highlights the need for more early spectra to explain the diversity of flash feature morphology in type II SNe.

SN~2018aoz is a Type Ia SN with a $B$-band plateau and excess emission in the infant-phase light curves $\lesssim$ 1 day after first light, evidencing an over-density of surface iron-peak elements as shown in our previous study. Here, we advance the constraints on the nature and origin of SN~2018aoz based on its evolution until the nebular phase. Near-peak spectroscopic features show the SN is intermediate between two subtypes of normal Type Ia: Core-Normal and Broad-Line. The excess emission could have contributions from the radioactive decay of surface iron-peak elements as well as ejecta interaction with either the binary companion or a small torus of circumstellar material. Nebular-phase limits on H$\alpha$ and He~I favour a white dwarf companion, consistent with the small companion size constrained by the low early SN luminosity, while the absence of [O~I] and He~I disfavours a violent merger of the progenitor. Of the two main explosion mechanisms proposed to explain the distribution of surface iron-peak elements in SN~2018aoz, the asymmetric Chandrasekhar-mass explosion is less consistent with the progenitor constraints and the observed blueshifts of nebular-phase [Fe~II] and [Ni~II]. The helium-shell double-detonation explosion is compatible with the observed lack of C spectral features, but current 1-D models are incompatible with the infant-phase excess emission, $B_{\rm max}-V_{\rm max}$ color, and absence of nebular-phase [Ca~II]. Although the explosion processes of SN~2018aoz still need to be more precisely understood, the same processes could produce a significant fraction of Type Ia SNe that appear normal after $\sim$ 1 day., Comment: Submitted for publication in ApJ. 35 pages, 16 figures, 7 tables

We present high-cadence optical and ultraviolet light curves of the normal Type Ia supernova (SN) 2021aefx, which shows an early bump during the first two days of observation. This bump may be a signature of interaction between the exploding white dwarf and a nondegenerate binary companion, or it may be intrinsic to the white dwarf explosion mechanism. In the case of the former, the short duration of the bump implies a relatively compact main-sequence companion star, although this conclusion is viewing-angle dependent. Our best-fit companion-shocking and double-detonation models both overpredict the UV luminosity during the bump, and existing nickel-shell models do not match the strength and timescale of the bump. We also present nebular spectra of SN 2021aefx, which do not show the hydrogen or helium emission expected from a nondegenerate companion, as well as a radio nondetection that rules out all symbiotic progenitor systems and most accretion disk winds. Our analysis places strong but conflicting constraints on the progenitor of SN 2021aefx; no current model can explain all of our observations., Comment: updated to match accepted version

We present high-cadence optical, ultraviolet (UV), and near-infrared data of the nearby ($D\approx23$ Mpc) Type II supernova (SN) 2021yja. Many Type II SNe show signs of interaction with circumstellar material (CSM) during the first few days after explosion, implying that their red supergiant (RSG) progenitors experience episodic or eruptive mass loss. However, because it is difficult to discover SNe early, the diversity of CSM configurations in RSGs has not been fully mapped. SN 2021yja, first detected within ${\approx}5.4$ hours of explosion, shows some signatures of CSM interaction (high UV luminosity, radio and x-ray emission) but without the narrow emission lines or early light curve peak that can accompany CSM. Here we analyze the densely sampled early light curve and spectral series of this nearby SN to infer the properties of its progenitor and CSM. We find that the most likely progenitor was an RSG with an extended envelope, encompassed by low-density CSM. We also present archival Hubble Space Telescope imaging of the host galaxy of SN 2021yja, which allows us to place a stringent upper limit of ${\lesssim}9\ M_\odot$ on the progenitor mass. However, this is in tension with some aspects of the SN evolution, which point to a more massive progenitor. Our analysis highlights the need to consider progenitor structure when making inferences about CSM properties, and that a comprehensive view of CSM tracers should be made to give a fuller view of the last years of RSG evolution., Comment: updated to match accepted version

Type Ia Supernovae are thermonuclear explosions of white dwarf stars. They play a central role in the chemical evolution of the Universe and are an important measure of cosmological distances. However, outstanding questions remain about their origins. Despite extensive efforts to obtain natal information from their earliest signals, observations have thus far failed to identify how the majority of them explode. Here, we present infant-phase detections of SN 2018aoz from a brightness of -10.5 absolute AB magnitudes -- the lowest luminosity early Type Ia signals ever detected -- revealing a hitherto unseen plateau in the $B$-band that results in a rapid redward color evolution between 1.0 and 12.4 hours after the estimated epoch of first light. The missing $B$-band flux is best-explained by line-blanket absorption from Fe-peak elements in the outer 1% of the ejected mass. The observed $B-V$ color evolution of the SN also matches the prediction from an over-density of Fe-peak elements in the same outer 1% of the ejected mass, whereas bluer colors are expected from a purely monotonic distribution of Fe-peak elements. The presence of excess nucleosynthetic material in the extreme outer layers of the ejecta points to enhanced surface nuclear burning or extended sub-sonic mixing processes in some normal Type Ia Supernova explosions., Comment: Accepted for publication in Nature Astronomy. Main text = 8 pages, 4 figures, 1 table; Full document = 46 pages, with Methods, Supplementary Information, 7 Supplementary figures, 2 Supplementary tables and references. Nat Astron (2022)

We present the results of photometric reverberation mapping observations on the changing look active galactic nucleus Mrk 590 at z = 0.026. The observations were carried out from July to December, 2018 using broad band B, R and narrow band H{\alpha} and S II filters. B-band traces the continuum emission from the accretion disk, R-band encompasses both the continuum emission from the accretion disk and the redshifted H{\alpha} line from the broad line region (BLR), S II band contains the redshifted H{\alpha} emission and the H{\alpha} band traces the continuum emission underneath the S II band. All the light curves showed strong variation with a fractional root-mean-square variation of 0.132 (+/-) 0.001 in B-band and 0.321 (+/-) 0.001 in H{\alpha} line. From cross-correlation function analysis, we obtained a delayed response of H{\alpha} line emission to the optical B-band continuum emission of 21.44(+1.49/-2.11) days in the rest-frame of the source, corresponding to a linear size of the BLR of 0.018 pc. This is consistent with previous estimates using H{\beta}. By combining the BLR size with the H{\alpha} line full width at half maximum of 6478 (+/-) 240 km/s measured from a single-epoch spectrum obtained with the Subaru telescope, we derived a black hole mass of 1.96 (+0.15 / -0.21) X 10^8 Msun., Comment: 16 pages, 12 figures, 4 tables, Accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS), 5 October 2021

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

We present multi-band photometry and spectroscopy of SN 2018cuf, a Type IIP ("P" for plateau) supernova (SN) discovered by the Distance Less Than 40 Mpc survey (DLT40) within 24 hours of explosion. SN 2018cuf appears to be a typical Type IIP SN, with an absolute $V$-band magnitude of $-$16.73 $\pm$ 0.32 at maximum and a decline rate of 0.21 $\pm$ 0.05 mag/50d during the plateau phase. The distance of the object is constrained to be 41.8 $\pm$ 5.7 Mpc by using the expanding photosphere method. We use spectroscopic and photometric observations from the first year after the explosion to constrain the progenitor of SN 2018cuf using both hydrodynamic light curve modelling and late-time spectroscopic modelling. The progenitor of SN 2018cuf was most likely a red supergiant of about 14.5 $\rm M_{\odot}$ that produced 0.04 $\pm$ 0.01 $\rm M_{\odot}$ $\rm ^{56}Ni$ during the explosion. We also found $\sim$ 0.07 $\rm M_{\odot}$ of circumstellar material (CSM) around the progenitor is needed to fit the early light curves, where the CSM may originate from pre-supernova outbursts. During the plateau phase, high velocity features at $\rm \sim 11000\ km~s^{-1}$ are detected both in the optical and near-infrared spectra, supporting the possibility that the ejecta were interacting with some CSM. A very shallow slope during the post-plateau phase is also observed and it is likely due to a low degree of nickel mixing or the relatively high nickel mass in the SN., Comment: 22 pages, 15 figures, submitted to ApJ

We describe the GW follow-up strategy and subsequent results of the DLT40 during the second science run (O2) of the LVC. Depending on the information provided in the GW alert together with the localization map sent by the LVC, DLT40 would respond promptly to image the corresponding galaxies selected by our ranking algorithm in order to search for possible EM counterparts in real time. During the LVC O2 run, DLT40 followed ten GW triggers, observing between $\sim$20-100 galaxies within the GW localization area of each event. From this campaign, we identified two real transient sources within the GW localizations with an appropriate on-source time -- one was an unrelated type Ia supernova (SN~2017cbv), and the other was the optical kilonova, AT 2017fgo/SSS17a/DLT17ck, associated with the binary neutron star coalescence GW170817 (a.k.a gamma-ray burst GRB170817A). We conclude with a discussion of the DLT40 survey's plans for the upcoming LVC O3 run, which include expanding our galaxy search fields out to $D\approx$65 Mpc to match the LVC's planned three-detector sensitivity for binary neutron star mergers., Comment: 38 pages, 8 figures and 14 tables, Submitted to ApJ

Using a global network of small telescopes, we have obtained light curves of Proxima Centauri at 329 observation epochs from 2006 - 2017. The planet Proxima b discovered by Anglada-Escude et al. (2016) with an orbital period of 11.186 d has an a priori transit probability of ~1.5%; if it transits, the predicted transit depth is about 5 millimagnitudes. In Blank et al. (2018), we analyzed 96 of our light curves that overlapped with predicted transit ephemerides from previously published tentative transit detections, and found no evidence in our data that would corroborate claims of transits with a period of 11.186 d. Here we broaden our analysis, using 262 high-quality light curves from our data set to search for any periodic transit-like events over a range of periods from 1 - 30 d. We also inject a series of simulated planet transits and find that our data are sufficiently sensitive to have detected transits of 5 millimagnitude depth, with recoverability ranging from ~100% for an orbital period of 1 d to ~20% for an orbital period of 20 d for the parameter spaces tested. Specifically at the 11.186 d period and 5 millimagnitude transit depth, we rule out transits in our data with high confidence. We are able to rule out virtually all transits of other planets at periods shorter than 5 d and depths greater than 3 millimagnitudes; however, we cannot confidently rule out transits at the period of Proxima b due to incomplete orbital phase coverage and a lack of sensitivity to transits shallower than 4 millimagnitudes., Comment: 11 pages, 10 figures, 4 Tables, accepted to ApJ

We present photometric and spectroscopic analysis of the peculiar core-collapse SN 2013gc, spanning seven years of observations. The light curve shows an early maximum followed by a fast decline and a phase of almost constant luminosity. At +200 days from maximum, a brightening of 1 mag is observed in all bands, followed by a steep linear luminosity decline after +300 d. In archival images taken between 1.5 and 2.5 years before the explosion, a weak source is visible at the supernova location, with mag$\approx$20. The early supernova spectra show Balmer lines, with a narrow ($\sim$560 km s$^{-1}$) P-Cygni absorption superimposed on a broad ($\sim$3400 km s$^{-1}$) component, typical of type IIn events. Through a comparison of colour curves, absolute light curves and spectra of SN 2013gc with a sample of supernovae IIn, we conclude that SN 2013gc is a member of the so-called type IId subgroup. The complex profile of the H$\alpha$ line suggests a composite circumstellar medium geometry, with a combination of lower velocity, spherically symmetric gas and a more rapidly expanding bilobed feature. This circumstellar medium distribution has been likely formed through major mass-loss events, that we directly observed from 3 years before the explosion. The modest luminosity ($M_I\sim-16.5$ near maximum) of SN 2013gc at all phases, the very small amount of ejected $^{56}$Ni (of the order of $10^{-3}$ M$_\odot$), the major pre-supernova stellar activity and the lack of prominent [O I] lines in late-time spectra support a fall-back core-collapse scenario for the massive progenitor of SN~2013gc., Comment: 20 pages, 11 figures, 8 tables, accepted by MNRAS

Proxima Centauri has become the subject of intense study since the radial-velocity discovery by Anglada-Escud\'e et al. 2016 of a planet orbiting this nearby M-dwarf every ~ 11.2 days. If Proxima Centauri b transits its host star, independent confirmation of its existence is possible, and its mass and radius can be measured in units of the stellar host mass and radius. To date, there have been three independent claims of possible transit-like event detections in light curve observations obtained by the MOST satellite (in 2014-15), the BSST telescope in Antarctica (in 2016), and the Las Campanas Observatory (in 2016). The claimed possible detections are tentative, due in part to the variability intrinsic to the host star, and in the case of the ground-based observations, also due to the limited duration of the light curve observations. Here, we present preliminary results from an extensive photometric monitoring campaign of Proxima Centauri, using telescopes around the globe and spanning from 2006 to 2017, comprising a total of 329 observations. Considering our data that coincide directly and/or phased with the previously published tentative transit detections, we are unable to independently verify those claims. We do, however, verify the previously reported ubiquitous and complex variability of the host star. We discuss possible interpretations of the data in light of the previous claims, and we discuss future analyses of these data that could more definitively verify or refute the presence of transits associated with the radial-velocity discovered planet., Comment: Published in The Astronomical Journal, Volume 155, Number 6. 15 pages, 13 figures, and 2 tables. Updated Table 2 FOVs

CS 1246 is a hot subdwarf B star discovered in 2009 to exhibit a single, large-amplitude radial pulsation. An O-C diagram constructed from this mode revealed reflex motion due to the presence of a low-mass M dwarf, as well as a long-term trend consistent with a decrease in the pulsational period. The orbital reflex motion was later confirmed with radial velocity measurements. Using eight years of data collected with the Skynet Robotic Telescope Network, we show that the pulsation amplitude of CS 1246 is decaying nonlinearly. We also present an updated O-C diagram, which might now indicate a positive $\dot P$ and a new $2.09 \pm 0.05$ yr oscillation consistent with orbital reflex motion of the entire inner sdB+dM binary, possibly due to the gravitational influence of a circumbinary planet with minimum mass $m\sin i = 3.3 \, \pm \, ^{2.1}_{1.2}$ $M_{\rm Jup}$. However, unlike the presence of the M dwarf, we hesitate to claim this object as a definitive detection since intrinsic variability of the pulsation phase could theoretically produce a similar effect., Comment: 6 pages, 4 figures, 2 tables; to appear in the Proceedings of the Eighth Meeting on Hot Subdwarfs and Related Objects, 9 to 15 July 2017, Cracow, Poland. Eds. A. Baran, A. E. Lynas-Gray, Open Astronomy, in press

During the second observing run of the Laser Interferometer gravitational- wave Observatory (LIGO) and Virgo Interferometer, a gravitational-wave signal consistent with a binary neutron star coalescence was detected on 2017 August 17th (GW170817), quickly followed by a coincident short gamma-ray burst trigger by the Fermi satellite. The Distance Less Than 40 (DLT40) Mpc supernova search performed pointed follow-up observations of a sample of galaxies regularly monitored by the survey which fell within the combined LIGO+Virgo localization region, and the larger Fermi gamma ray burst error box. Here we report the discovery of a new optical transient (DLT17ck, also known as SSS17a; it has also been registered as AT 2017gfo) spatially and temporally coincident with GW170817. The photometric and spectroscopic evolution of DLT17ck are unique, with an absolute peak magnitude of Mr = -15.8 \pm 0.1 and an r-band decline rate of 1.1mag/d. This fast evolution is generically consistent with kilonova models, which have been predicted as the optical counterpart to binary neutron star coalescences. Analysis of archival DLT40 data do not show any sign of transient activity at the location of DLT17ck down to r~19 mag in the time period between 8 months and 21 days prior to GW170817. This discovery represents the beginning of a new era for multi-messenger astronomy opening a new path to study and understand binary neutron star coalescences, short gamma-ray bursts and their optical counterparts., Comment: ApJL in press, 4 figures

Binary neutron star mergers are important to understand stellar evolution, the chemical enrichment of the universe via the r-process, the physics of short gamma-ray bursts, gravitational waves and pulsars. The rates at which these coalescences happen is uncertain, but it can be constrained in different ways. One of those is to search for the optical transients produced at the moment of the merging, called a kilonova, in ongoing SN searches. However, until now, only theoretical models for kilonovae light curve were available to estimate their rates. The recent kilonova discovery AT~2017gfo/DLT17ck gives us the opportunity to constrain the rate of kilonovae using the light curve of a real event. We constrain the rate of binary neutron star mergers using the DLT40 Supernova search, and the native AT~2017gfo/DLT17ck light curve obtained with the same telescope and software system. Excluding AT~2017gfo/DLT17ck due to visibility issues, which was only discovered thanks to the aLIGO/aVirgo trigger, no other similar transients detected during 13 months of daily cadence observations of $\sim$ 2200 nearby ($<$40 Mpc) galaxies. We find that the rate of BNS mergers is lower than 0.47 - 0.55 kilonovae per 100 years per $10^{10}$ $L_{B_{\odot}}$ (depending on the adopted extinction distribution). In volume, this translates to $<0.99\times 10^{-4}\,_{-0.15}^{+0.19},\rm{Mpc^{-3}}\,\rm{yr^{-1}}$(SNe Ia-like extinction distribution), consistent with previous BNS coalescence rates. Based on our rate limit, and the sensitivity of aLIGO/aVirgo during O2, it is very unlikely that kilonova events are lurking in old pointed galaxy SN search datasets., Comment: 3 figures, 2 tables