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)

The optical spectroscopic observations of ASASSN-18fv observed from 2018 March 24 to 2019 Jan 26 are presented. The optical spectra are obtained from Mirranook, Vainu Bappu and South African Astronomical observatories. The spectra are dominated by hydrogen Balmer, Fe II and O I lines with P-Cygni profiles in the early phase, typical of an Fe II class nova. The spectra show He I lines along with H I and O I emission lines in the decline phase placing the nova in the hybrid class of novae. The spectra show rapid development in high ionization lines in this phase. Analysis of the light curve indicate t$_2$ and t$_3$ values of about 50 and 70 days respectively placing the nova in the category of moderately fast nova. The ejecta geometry, inclination and position angle are estimated using morpho-kinematic analysis. The geometry of the ejecta is found to be an asymmetric bipolar structure with an inclination angle of about 53$^{\circ}$. The ejected mass using photo-ionization analysis is found to be 6.07 $\times$ 10$^{-4}$ M$_{\odot}$.

The third observing run by LVC has brought the discovery of many compact binary coalescences. Following the detection of the first binary neutron star merger in this run (LIGO/Virgo S190425z), we performed a dedicated follow-up campaign with the Zwicky Transient Facility (ZTF) and Palomar Gattini-IR telescopes. The initial skymap of this single-detector gravitational wave (GW) trigger spanned most of the sky observable from Palomar Observatory. Covering 8000 deg$^2$ of the initial skymap over the next two nights, corresponding to 46\% integrated probability, ZTF system achieved a depth of $\approx$\,21 $m_\textrm{AB}$ in $g$- and $r$-bands. Palomar Gattini-IR covered 2200 square degrees in $J$-band to a depth of 15.5\,mag, including 32\% integrated probability based on the initial sky map. The revised skymap issued the following day reduced these numbers to 21\% for the Zwicky Transient Facility and 19\% for Palomar Gattini-IR. We narrowed 338,646 ZTF transient "alerts" over the first two nights of observations to 15 candidate counterparts. Two candidates, ZTF19aarykkb and ZTF19aarzaod, were particularly compelling given that their location, distance, and age were consistent with the GW event, and their early optical lightcurves were photometrically consistent with that of kilonovae. These two candidates were spectroscopically classified as young core-collapse supernovae. The remaining candidates were ruled-out as supernovae. Palomar Gattini-IR did not identify any viable candidates with multiple detections only after merger time. We demonstrate that even with single-detector GW events localized to thousands of square degrees, systematic kilonova discovery is feasible.

We study the spectroscopic and ionized structural evolution of T Pyx during its 2011 outburst, and also study the variation in degree of polarization during its early phase. Optical spectroscopic data of this system obtained from day 1.28--2415.62 since discovery, and optical, broadband imaging polarimetric observations obtained from day 1.36--29.33 during the early phases of the outburst are used in the study. The physical conditions and the geometry of the ionized structure of the nova ejecta has been modelled for a few epochs using the photo-ionization code, CLOUDY in 1D and pyCloudy in 3D. The spectral evolution of the nova ejecta during its 2011 outburst is similar to that of the previous outbursts. The variation in the line profiles is seen very clearly in the early stages due to good coverage during this period. The line profiles vary from P-Cygni (narrower, deeper and sharper) to emission profiles which are broader and structured, which later become narrower and sharper in the late post-outburst phase. The average ejected mass is estimated to be $7.03\, \times\, 10^{-6}\, M_{\odot}$. The ionized structure of the ejecta is found to be a bipolar conical structure with equatorial rings, with a low inclination angle of $14.75^\circ\, \pm\, 0.65^\circ$.

Using our previously reported observations, we derive some physical parameters of the moderately fast nova V2676 Ophiuchi 2012 # 1. The best-fit CLOUDY model of the nebular spectrum obtained on 2015 May 8 shows a hot white dwarf source with Tbb = 1.0 x 10^{5} K having a luminosity of 1.0 x 10^{38} ergs/s. Our abundance analysis shows that the ejecta are significantly enhanced relative to solar, He/H = 2.14, O/H = 2.37, S/H = 6.62 and Ar/H = 3.25. The ejecta mass is estimated to be 1.42 x 10^{-5} Msun. The nova showed a pronounced dust formation phase after 90 days from discovery. The J-H and H-K colors were very large as compared to other molecule- and dust-forming novae in recent years. The dust temperature and mass at two epochs have been estimated from spectral energy distribution (SED) fits to infrared photometry.

Since its discovery in 2008, the Andromeda galaxy nova M31N 2008-12a has been observed in eruption every single year. This unprecedented frequency indicates an extreme object, with a massive white dwarf and a high accretion rate, which is the most promising candidate for the single-degenerate progenitor of a type-Ia supernova known to date. The previous three eruptions of M31N 2008-12a have displayed remarkably homogeneous multi-wavelength properties: (i) From a faint peak, the optical light curve declined rapidly by two magnitudes in less than two days; (ii) Early spectra showed initial high velocities that slowed down significantly within days and displayed clear He/N lines throughout; (iii) The supersoft X-ray source (SSS) phase of the nova began extremely early, six days after eruption, and only lasted for about two weeks. In contrast, the peculiar 2016 eruption was clearly different. Here we report (i) the considerable delay in the 2016 eruption date, (ii) the significantly shorter SSS phase, and (iii) the brighter optical peak magnitude (with a hitherto unobserved cusp shape). Early theoretical models suggest that these three different effects can be consistently understood as caused by a lower quiescence mass-accretion rate. The corresponding higher ignition mass caused a brighter peak in the free-free emission model. The less-massive accretion disk experienced greater disruption, consequently delaying re-establishment of effective accretion. Without the early refueling, the SSS phase was shortened. Observing the next few eruptions will determine whether the properties of the 2016 outburst make it a genuine outlier in the evolution of M31N 2008-12a.