PG 1553+113 is one of the few blazars with a convincing quasi-periodic emission in the gamma-ray band. The source is also a very high-energy (VHE; >100 GeV) gamma-ray emitter. To better understand its properties and identify the underlying physical processes driving its variability, the MAGIC Collaboration initiated a multiyear, multiwavelength monitoring campaign in 2015 involving the OVRO 40-m and Medicina radio telescopes, REM, KVA, and the MAGIC telescopes, Swift and Fermi satellites, and the WEBT network. The analysis presented in this paper uses data until 2017 and focuses on the characterization of the variability. The gamma-ray data show a (hint of a) periodic signal compatible with literature, but the X-ray and VHE gamma-ray data do not show statistical evidence for a periodic signal. In other bands, the data are compatible with the gamma-ray period, but with a relatively high p-value. The complex connection between the low and high-energy emission and the non-monochromatic modulation and changes in flux suggests that a simple one-zone model is unable to explain all the variability. Instead, a model including a periodic component along with multiple emission zones is required.

The Zwicky Transient Facility (ZTF) is a new optical time-domain survey that uses the Palomar 48-inch Schmidt telescope. A custom-built wide-field camera provides a 47 deg$^2$ field of view and 8 second readout time, yielding more than an order of magnitude improvement in survey speed relative to its predecessor survey, the Palomar Transient Factory (PTF). We describe the design and implementation of the camera and observing system. The ZTF data system at the Infrared Processing and Analysis Center provides near-real-time reduction to identify moving and varying objects. We outline the analysis pipelines, data products, and associated archive. Finally, we present on-sky performance analysis and first scientific results from commissioning and the early survey. ZTF's public alert stream will serve as a useful precursor for that of the Large Synoptic Survey Telescope.

The InfraRed Imaging Spectrograph (IRIS) will be a first-light client instrument for the Narrow Field Infrared Adaptive Optics System (NFIRAOS) on the Thirty Meter Telescope. IRIS includes three configurable tip/tilt (TT) or tip/tilt/focus (TTF) On-Instrument Wavefront Sensors (OIWFS). These sensors are positioned over natural guide star (NGS) asterisms using movable polar-coordinate pick-off arms (POA) that patrol an approximately 2-arcminute circular field-of-view (FOV). The POAs are capable of colliding with one another, so an algorithm for coordinated motion that avoids contact is required. We have adopted an approach in which arm motion is evaluated using the gradient descent of a scalar potential field that includes an attractive component towards the goal configuration (locations of target stars), and repulsive components to avoid obstacles (proximity to adjacent arms). The resulting vector field is further modified by adding a component transverse to the repulsive gradient to avoid problematic local minima in the potential. We present path planning simulations using this computationally inexpensive technique, which exhibit smooth and efficient trajectories.

We present results from a systematic selection of tidal disruption events (TDEs) in a wide-area (4800~deg$^2$), $g+R$ band, Intermediate Palomar Transient Factory (iPTF) experiment. Our selection targets typical optically-selected TDEs: bright (>60\% flux increase) and blue transients residing in the center of red galaxies. Using photometric selection criteria to down-select from a total of 493 nuclear transients to a sample of 26 sources, we then use follow-up UV imaging with the Neil Gehrels Swift Telescope, ground-based optical spectroscopy, and light curve fitting to classify them as 14 Type Ia supernovae (SNe Ia), 9 highly variable active galactic nuclei (AGNs), 2 confirmed TDEs, and 1 potential core-collapse supernova. We find it possible to filter AGNs by employing a more stringent transient color cut (g-r < $-$0.2 mag); further, UV imaging is the best discriminator for filtering SNe, since SNe Ia can appear as blue, optically, as TDEs in their early phases. However, when UV-optical color is unavailable, higher precision astrometry can also effectively reduce SNe contamination in the optical. Our most stringent optical photometric selection criteria yields a 4.5:1 contamination rate, allowing for a manageable number of TDE candidates for complete spectroscopic follow-up and real-time classification in the ZTF era. We measure a TDE per galaxy rate of 1.7$^{+2.9}_{-1.3}$ $\times$10$^{-4}$ gal$^{-1}$ yr$^{-1}$ (90\% CL in Poisson statistics). This does not account for TDEs outside our selection criteria, thus may not reflect the total TDE population, which is yet to be fully mapped.

We report the discovery of six active galactic nuclei (AGN) caught "turning on" during the first nine months of the Zwicky Transient Facility (ZTF) survey. The host galaxies were classified as LINERs by weak narrow forbidden line emission in their archival SDSS spectra, and detected by ZTF as nuclear transients. In five of the cases, we found via follow-up spectroscopy that they had transformed into broad-line AGN, reminiscent of the changing-look LINER iPTF 16bco. In one case, ZTF18aajupnt/AT2018dyk, follow-up HST UV and ground-based optical spectra revealed the transformation into a narrow-line Seyfert 1 (NLS1) with strong [Fe VII, X, XIV] and He II 4686 coronal lines. Swift monitoring observations of this source reveal bright UV emission that tracks the optical flare, accompanied by a luminous soft X-ray flare that peaks ~60 days later. Spitzer follow-up observations also detect a luminous mid-infrared flare implying a large covering fraction of dust. Archival light curves of the entire sample from CRTS, ATLAS, and ASAS-SN constrain the onset of the optical nuclear flaring from a prolonged quiescent state. Here we present the systematic selection and follow-up of this new class of changing-look LINERs, compare their properties to previously reported changing-look Seyfert galaxies, and conclude that they are a unique class of transients well-suited to test the uncertain physical processes associated with the LINER accretion state.

IRIS (InfraRed Imaging Spectrograph) is the diffraction-limited first light instrument for the Thirty Meter Telescope (TMT) that consists of a near-infrared (0.84 to 2.4 $\mu$m) imager and integral field spectrograph (IFS). The IFS makes use of a lenslet array and slicer for spatial sampling, which will be able to operate in 100's of different modes, including a combination of four plate scales from 4 milliarcseconds (mas) to 50 mas with a large range of filters and gratings. The imager will have a field of view of 34$\times$34 arcsec$^{2}$ with a plate scale of 4 mas with many selectable filters. We present the preliminary design of the data reduction system (DRS) for IRIS that need to address all of these observing modes. Reduction of IRIS data will have unique challenges since it will provide real-time reduction and analysis of the imaging and spectroscopic data during observational sequences, as well as advanced post-processing algorithms. The DRS will support three basic modes of operation of IRIS; reducing data from the imager, the lenslet IFS, and slicer IFS. The DRS will be written in Python, making use of open-source astronomical packages available. In addition to real-time data reduction, the DRS will utilize real-time visualization tools, providing astronomers with up-to-date evaluation of the target acquisition and data quality. The quicklook suite will include visualization tools for 1D, 2D, and 3D raw and reduced images. We discuss the overall requirements of the DRS and visualization tools, as well as necessary calibration data to achieve optimal data quality in order to exploit science cases across all cosmic distance scales.

We present the current optical design for the IRIS Atmospheric Dispersion Corrector (ADC). The ADC is designed for residual dispersions less than ~1 mas across a given passband at elevations of 25 degrees. Since the last report, the area of the IRIS Imager has increased by a factor of four, and the pupil size has increased from 75 to 90mm, both of which contribute to challenges with the design. Several considerations have led to the current design: residual dispersion, amount of introduced distortion, glass transmission, glass availability, and pupil displacement. In particular it was found that there are significant distortions that appear (two different components) that can lead to image blur over long exposures. Also, pupil displacement increases the wavefront error at the imager focus. We discuss these considerations, discuss the compromises, and present the final design choice and expected performance.

We study the properties of the galaxies hosting the first 19 tidal disruption events (TDEs) detected with the Zwicky Transient Facility (ZTF) within the context of a carefully constructed, representative host galaxy sample. We find that the ZTF sample of TDE hosts is dominated by compact "green valley" galaxies. After we restrict the comparison sample to galaxies with a similar concentration, as measured by Sersic index, we find this green valley over representation is even larger. That is, concentrated red sequence galaxies are not producing TDEs at elevated levels. We present host galaxy spectra which show that E+A galaxies are overrepresented in the ZTF sample by a factor of $\approx$22, which is lower than previous TDE host galaxy studies have found. We find that this overrepresentation can be fully accounted for when taking into account the masses, colors, and Sérsic indices of the ZTF TDE hosts. The combination of both green colors and high Sérsic index of the typical TDE host galaxy could be explained if the TDE rate is temporarily enhanced following a merger that leads to a higher central concentration of stars.

The Zwicky Transient Facility (ZTF) is conducting a wide-field survey of the northern sky in three optical bands and the collaboration cosmology working group has released 3628 spectroscopically confirmed Type Ia supernovae (SNe Ia) discovered during its first 2.5 years of operation. This "ZTF SN Ia DR2" sample is the largest SN Ia dataset to date. Fully exploiting this dataset to improve understanding of the properties of dark energy requires a photometric accuracy of O(0.1%). This can be achieved using Scene Modeling Photometry (SMP), which is optimal to extract a transient signal (SN) from a complex background (its host), while ensuring a common flux estimator with nearby stars used as calibration reference. In this paper, we present the status of the SMP development and use it to assess the precision and accuracy of the ZTF SN Ia DR2 force photometry light curves. We reach a repeatability of the star observations better than 1%. However, we have identified a new sensor effect, dubbed "pocket-effect", which distorts the Point Spread Function (PSF) in a flux-dependent manner leading to non-linearities in the photometry of a few percent. Correcting for this effect requires time- and sensor-dependent corrections to be applied at the pixel level, which is currently under development. This effects affects all light curve releases to date -- both from forced photometry and scene modelling preventing ZTF SN Ia DR2 to be used for accurate cosmological inference. Comparing the SMP and forced photometry measurements, we find that stretch and color estimated from both processings are consistent, aside from a 10 mmag shift in color. This assess the robustness of results presented as part of the the ZTF SN Ia DR2 release. The absolute calibration however shifts by 90 mmag. A reprocessing of the full ZTF SN Ia DR2 dataset using the SMP method is currently in progress.

IRIS is a near-infrared (0.84 to 2.4 microns) integral field spectrograph and wide-field imager being developed for first light with the Thirty Meter (TMT). It mounts to the advanced optics (AO) system NFIRAOS and has integrated on-instrument wavefront sensors (OIWFS) to achieve diffraction-limited spatial resolution at wavelengths longer than 1 micron. With moderate spectral resolution (R ~4,000 - 8,000) and large bandpass over a continuous field of view, IRIS will open new opportunities in virtually every area of astrophysical science. It will be able to resolve surface features tens of kilometers across Titan, while also mapping the distant galaxies at the scale of an individual star forming region. This paper summarizes the entire design and capabilities, and includes the results from the nearly completed preliminary design phase.

The InfraRed Imaging Spectrograph (IRIS) will be the first light adaptive optics instrument on the Thirty Meter Telescope (TMT). IRIS is being built by a collaboration between Caltech, the University of California, NAOJ and NRC Herzberg. We present novel aspects of the Support Structure, Rotator and On-Instrument Wavefront Sensor systems being developed at NRC Herzberg. IRIS is suspended from the bottom port of the Narrow Field Infrared Adaptive Optics System (NFIRAOS), and provides its own image de-rotation to compensate for sidereal rotation of the focal plane. This arrangement is a challenge because NFIRAOS is designed to host two other science instruments, which imposes strict mass requirements on IRIS. We have been tasked with keeping the instrument mass under seven tonnes which has resulted in a mass reduction of 30 percent for the support structure and rotator compared to the most recent IRIS designs. To accomplish this goal, while still being able to withstand earthquakes, we developed a new design with composite materials. As IRIS is a client instrument of NFIRAOS, it benefits from NFIRAOS's superior AO correction. IRIS assists this correction by sensing low-order aberrations with an On-Instrument Wavefront Sensor (OIWFS). The OIWFS consists of three independently positioned natural guide star wavefront sensors that patrol a 2-arcminute field of view. We expect tip-tilt measurements from faint stars within the IRIS imager focal plane will further stabilize the delivered image quality. We describe how the use of On-Detector Guide Windows (ODGWs) in the IRIS imager can be incorporated into the AO correction. Finally, we present our strategies for acquiring and tracking sources with this complex AO system, and for mitigating and measuring the various potential sources of image blur and misalignment due to properties of the mechanical structure and interfaces. (Abridged)

Detecting gravitationally lensed supernovae is among the biggest challenges in astronomy. It involves a combination of two very rare phenomena: catching the transient signal of a stellar explosion in a distant galaxy and observing it through a nearly perfectly aligned foreground galaxy that deflects light towards the observer. High-cadence optical observations with the Zwicky Transient Facility, with an unparalleled large field of view, led to the detection of a multiply-imaged Type Ia supernova (SN Ia), ``SN Zwicky", a.k.a. SN 2022qmx. Magnified nearly twenty-five times, the system was found thanks to the ``standard candle" nature of SNe Ia. High-spatial resolution imaging with the Keck telescope resolved four images of the supernova with very small angular separation, corresponding to an Einstein radius of only $\theta_E =0.167"$and almost identical arrival times. The small$\theta_E$ and faintness of the lensing galaxy is very unusual, highlighting the importance of supernovae to fully characterise the properties of galaxy-scale gravitational lenses, including the impact of galaxy substructures.

Water ice is a fundamental building material of comets and other bodies in the outer solar system. Yet, the properties of cometary water ice are challenging to study, due to its volatility and the typical distances at which comets are observed. Cometary outbursts, impulsive mass-loss events that can liberate large amounts of material, offer opportunities to directly observe and characterize cometary water ice. We present a study of comet 243P/NEAT, instigated by a $-3$ mag outburst that occurred in December 2018. Optical images and a 251-day lightcurve were examined to characterize the outburst and the comet's quiescent activity. Variations in the quiescent lightcurve appear to be dominated by coma asymmetries, rather than changing activity levels as the comet approached and receded from the Sun. Furthermore, the lightcurve shows evidence for 1 to 2 additional small outbursts ($-0.3$ mag) occurring in September 2018. The large December 2018 outburst likely ejected water ice grains, yet no signatures of ice were found in color photometry, a color map, nor a near-infrared spectrum. We discuss possible dynamical and thermal reasons for this non-detection. In this context, we examined the comae of comets 103P/Hartley 2 and C/2013 US$_{10}$ (Catalina), and show that a one-to-one mapping between continuum color and the presence of water ice cannot be supported. We also discuss possible causes for the large outburst, and find that there is an apparent grouping in the kinetic energy per mass estimates for the outbursts of 5 comets.

We report an observation of a transit of the hot Jupiter (HJ) KELT-23A b with the Keck Planet Finder spectrograph and a measurement of the sky-projected obliquity ($\lambda$) of its Sun-like ($T_{\rm eff} \approx 5900$ K) host star. We measured a projected stellar obliquity of $\lambda \approx 180^\circ$, indicating that the orbit of the HJ is retrograde relative to the direction of the stellar spin. Due to the slow sky-projected rotational velocity of the host star ($v \sin{i_\star} \approx 0.5$ km s$^{-1}$), the true orbit of the HJ could be closer to polar. HJs around stars with effective temperatures below the Kraft break -- such as KELT-23A -- are generally found to have prograde orbits that are well-aligned with the equatorial planes of their host stars (i.e., $\lambda \sim 0^\circ$), most likely due to spin-orbit realignment driven by stellar tidal dissipation. This system is therefore a unique outlier that strains migration and tidal theories. The fact that the HJ has a highly misaligned orbit may suggest that the planet arrived at its close-in orbit relatively recently, possibly via interactions with the wide-separation (570 AU) M-dwarf companion in the system, or that it has stalled near an antialigned or polar orientation while realigning. Using Gaia DR3, we determined the orbit of the stellar companion to be moderately face-on ($\gamma = 60 \pm 4^\circ$). We show that the distribution of observed systems in the $\gamma - \lambda$ plane can be broadly reproduced using a toy model in which the orbits of the planetary and stellar companions begin aligned with the equatorial plane of the primary star and, upon migrating inwards, the planet preferentially obtains either an aligned or polar orbit.

Although many cases of stellar spin-orbit misalignment are known, it is usually unclear whether a single planet's orbit was tilted or if the entire protoplanetary disk was misaligned. Measuring stellar obliquities in multi-transiting planetary systems helps to distinguish these possibilities. Here, we present a measurement of the sky-projected spin-orbit angle for TOI-880 c (TOI-880.01), a member of a system of three transiting planets, using the Keck Planet Finder (KPF). We found that the host star is a K-type star ($T_{\rm eff}=5050 \pm 100$ K). Planet b (TOI-880.02) has a radius of $2.19\pm0.11\mathrm{R_{\oplus}}$ and an orbital period of $2.6$ days; planet c (TOI-880.01) is a Neptune-sized planet with $4.95\pm0.20\mathrm{R_{\oplus}}$ on a $6.4$-day orbit; and planet d (TOI-880.03) has a radius of $3.40_{-0.21}^{+0.22}\mathrm{R_{\oplus}}$ and a period of $14.3$ days. By modeling the Rossiter-McLaughlin (RM) effect, we found the sky-projected obliquity to be $|\lambda_c| = 7.4_{-7.2}^{+6.8}$$^{\circ}$, consistent with a prograde, well-aligned orbit. The lack of detectable rotational modulation of the flux of the host star and a low $\rm v\sin{i_\star}$ (1.6~km/s) imply slow rotation and correspondingly slow nodal precession of the planetary orbits and the expectation that the system will remain in this coplanar configuration. TOI-880 joins a growing sample of well-aligned, coplanar, multi-transiting systems. Additionally, TOI-880 c is a promising target for JWST follow-up, with a transmission spectroscopy metric (TSM) of $\sim 170$. We could not detect clear signs of atmospheric erosion in the H$\alpha$ line from TOI-880 c, as photoevaporation might have diminished for this mature planet.

The nature of the progenitor systems and explosion mechanisms that give rise to Type Ia supernovae (SNe Ia) are still debated. The interaction signature of circumstellar material (CSM) being swept up by expanding ejecta can constrain the type of system from which it was ejected. Most previous studies have focused on finding CSM ejected shortly before the SN Ia explosion still residing close to the explosion site, resulting in short delay times until the interaction starts. We use a sample of 3627 SNe Ia from the Zwicky Transient Facility discovered between 2018 and 2020 and search for interaction signatures over 100 days after peak brightness. By binning the late-time light curve data to push the detection limit as deep as possible, we identify potential late-time rebrightening in 3 SNe Ia (SN 2018grt, SN 2019dlf, SN 2020tfc). The late-time detections occur between 550 and 1450 d after peak brightness, have mean absolute $r$-band magnitudes of -16.4 to -16.8 mag and last up to a few hundred days, significantly brighter than the late-time CSM interaction discovered in the prototype SN 2015cp. The late-time detections all occur within 0.8 kpc of the host nucleus and are not easily explained by nuclear activity, another transient at a similar sky position, or data quality issues. This suggests environment or specific progenitor characteristics playing a role in producing potential CSM signatures in these SNe Ia. By simulating the ZTF survey we estimate that <0.5 per cent of normal SNe Ia display late-time strong H $\alpha$-dominated CSM interaction. This is equivalent to an absolute rate of $8_{-4}^{+20}$ to $54_{-26}^{+91}$ Gpc$^{-3}$ yr$^{-1}$ assuming a constant SN Ia rate of $2.4\times10^{-5}$ Mpc$^{-3}$ yr$^{-1}$ for $z \leq 0.1$. Weaker interaction signatures, more similar to the strength seen in SN 2015cp, could be more common but are difficult to constrain with our survey depth.

The Zwicky Transient Facility (ZTF) performs a systematic neutrino follow-up program, searching for optical counterparts to high-energy neutrinos with dedicated Target-of-Opportunity (ToO) observations. Since first light in March 2018, ZTF has taken prompt observations for 24 high-quality neutrino alerts from the IceCube Neutrino Observatory, with a median latency of 12.2 hours from initial neutrino detection. From two of these campaigns, we have already reported tidal disruption event (TDE) AT 2019dsg and likely TDE AT 2019fdr as probable counterparts, suggesting that TDEs contribute >7.8% of the astrophysical neutrino flux. We here present the full results of our program through to December 2021. No additional candidate neutrino sources were identified by our program, allowing us to place the first constraints on the underlying optical luminosity function of astrophysical neutrino sources. Transients with optical absolutes magnitudes brighter that $-21$ can contribute no more than 87% of the total, while transients brighter than $-22$ can contribute no more than 58% of the total, neglecting the effect of extinction and assuming they follow the star formation rate. These are the first observational constraints on the neutrino emission of bright populations such as superluminous supernovae. None of the neutrinos were coincident with bright optical AGN flares comparable to that observed for TXS 0506+056/IC170922A, with such optical blazar flares producing no more than 26% of the total neutrino flux. We highlight the outlook for electromagnetic neutrino follow-up programs, including the expected potential for the Rubin Observatory.

We present the goals, strategy and first results of the high-cadence Galactic plane survey using the Zwicky Transient Facility (ZTF). The goal of the survey is to unveil the Galactic population of short-period variable stars, including short period binaries and stellar pulsators with periods less than a few hours. Between June 2018 and January 2019, we observed 64 ZTF fields resulting in 2990 deg$^2$ of high stellar density in ZTF-$r$ band along the Galactic Plane. Each field was observed continuously for 1.5 to 6 hrs with a cadence of 40 sec. Most fields have between 200 and 400 observations obtained over 2-3 continuous nights. As part of this survey we extract a total of $\approx$230 million individual objects with at least 80 epochs obtained during the high-cadence Galactic Plane survey reaching an average depth of ZTF-$r$ $\approx$20.5 mag. For four selected fields with 2 million to 10 million individual objects per field we calculate different variability statistics and find that $\approx$1-2% of the objects are astrophysically variable over the observed period. We present a progress report on recent discoveries, including a new class of compact pulsators, the first members of a new class of Roche Lobe filling hot subdwarf binaries as well as new ultracompact double white dwarfs and flaring stars. Finally we present a sample of 12 new single-mode hot subdwarf B-star pulsators with pulsation amplitudes between ZTF-$r$ = 20-76 mmag and pulsation periods between $P$ = 5.8-16 min with a strong cluster of systems with periods $\approx$ 6 min. All of the data have now been released in either ZTF Data Release 3 or data release 4.

In over a thousand known cataclysmic variables (CVs), where a white dwarf is accreting from a hydrogen-rich star, only a dozen have orbital periods below 75 minutes. One way to achieve these short periods requires the donor star to have undergone substantial nuclear evolution prior to interacting with the white dwarf, and it is expected that these objects will transition to helium accretion. These transitional CVs have been proposed as progenitors of helium CVs. However, no known transitional CV is expected to reach an orbital period short enough to account for most of the helium CV population, leaving the role of this evolutionary pathway unclear. Here we report observations of ZTF J1813+4251, a 51-minute orbital period, fully eclipsing binary system consisting of a star with a temperature comparable to that of the Sun but a density 100 times greater due to its helium-rich composition, accreting onto a white dwarf. Phase-resolved spectra, multi-band light curves and the broadband spectral energy distribution allow us to obtain precise and robust constraints on the masses, radii and temperatures of both components. Evolutionary modeling shows that ZTF J1813+4251 is destined to become a helium CV binary, reaching an orbital period under 20 minutes, rendering ZTF J1813+4251 a previously missing link between helium CV binaries and hydrogen-rich CVs.