The Gravity Spy project aims to uncover the origins of glitches, transient bursts of noise that hamper analysis of gravitational-wave data. By using both the work of citizen-science volunteers and machine-learning algorithms, the Gravity Spy project enables reliable classification of glitches. Citizen science and machine learning are intrinsically coupled within the Gravity Spy framework, with machine-learning classifications providing a rapid first-pass classification of the dataset and enabling tiered volunteer training, and volunteer-based classifications verifying the machine classifications, bolstering the machine-learning training set and identifying new morphological classes of glitches. These classifications are now routinely used in studies characterizing the performance of the LIGO gravitational-wave detectors. Providing the volunteers with a training framework that teaches them to classify a wide range of glitches, as well as additional tools to aid their investigations of interesting glitches, empowers them to make discoveries of new classes of glitches. This demonstrates that, when giving suitable support, volunteers can go beyond simple classification tasks to identify new features in data at a level comparable to domain experts. The Gravity Spy project is now providing volunteers with more complicated data that includes auxiliary monitors of the detector to identify the root cause of glitches.

Polarization observations of the Milky Way and many other spiral galaxies have found a close correspondence between the orientation of spiral arms and magnetic field lines on scales of hundreds of parsecs. This paper presents polarization measurements at 214 $\mu$m toward ten filamentary candidate ``bones" in the Milky Way using the High-resolution Airborne Wide-band Camera (HAWC+) on the Stratospheric Observatory for Infrared Astronomy (SOFIA). These data were taken as part of the Filaments Extremely Long and Dark: A Magnetic Polarization Survey (FIELDMAPS) and represent the first study to resolve the magnetic field in spiral arms at parsec scales. We describe the complex yet well-defined polarization structure of all ten candidate bones, and we find a mean difference and standard deviation of $-74^{\circ} \pm 32^{\circ}$ between their filament axis and the plane-of-sky magnetic field, closer to a field perpendicular to their length rather than parallel. By contrast, the 850 $\mu$m polarization data from \textit{Planck} on scales greater than 10 pc show a nearly parallel mean difference of $3^{\circ} \pm 21^{\circ}$. These findings provide further evidence that magnetic fields can change orientation at the scale of dense molecular clouds, even along spiral arms. Finally, we use a power law to fit the dust polarization fraction as a function of total intensity on a cloud-by-cloud basis and find indices between $-0.6$ and $-0.9$, with a mean and standard deviation of $-0.7 \pm 0.1$. The polarization, dust temperature, and column density data presented in this work are publicly available online.

Time-domain surveys such as the Zwicky Transient Facility (ZTF) have opened a new frontier in the discovery and characterization of transients. While photometric light curves provide broad temporal coverage, spectroscopic observations remain crucial for physical interpretation and source classification. However, existing spectral analysis methods -- often reliant on template fitting or parametric models -- are limited in their ability to capture the complex and evolving spectra characteristic of such sources, which are sometimes only available at low resolution. In this work, we introduce SpectraNet, a deep convolutional neural network designed to learn robust representations of optical spectra from transients. Our model combines multi-scale convolution kernels and multi-scale pooling to extract features from preprocessed spectra in a hierarchical and interpretable manner. We train and validate SpectraNet on low-resolution time-series spectra obtained from the Spectral Energy Distribution Machine (SEDM) and other instruments, demonstrating state-of-the-art performance in classification. Furthermore, in redshift prediction tasks, SpectraNet achieves a root mean squared relative redshift error of 0.02, highlighting its effectiveness in precise regression tasks as well.

Stripped-envelope supernovae (SESNe) originate from massive stars that lose their envelopes through binary interactions or stellar winds. The connection between SESN subtypes and their progenitors remains poorly understood, as does the influence of initial mass, binarity, explodability, and metallicity on their evolutionary pathways, relative rates, ejecta masses, and progenitor ages. Here, we investigate these properties across a wide metallicity range (0.01-2 $Z_{\odot}$) using POSYDON, a state-of-the-art population synthesis code that incorporates detailed single- and binary-star model grids. We find that the common-envelope channel contributes less than 6% of SESNe, since unstable mass transfer is found less frequent than previously thought and rarely leads to CE survival when envelope binding energies are computed from detailed stellar models. The secondary channel accounts for less than 11%, while the vast majority of SESNe originate from primary stars in binaries undergoing stable mass-transfer episodes. These interactions maintain a largely metallicity-independent SESN parameter space, making the overall SESN rate almost insensitive to metallicity. In contrast, subtype fractions exhibit strong metallicity dependence, though their exact values remain affected by classification thresholds. The age distributions and therefore the progenitor masses of different SESN types also vary significantly with metallicity, revealing metallicity-dependent trends that can be tested observationally. Predicted SESN ejecta masses remain nearly constant across metallicity, in contrast to single-star models, and fall within observed ranges. Future transient surveys, combined with statistical environmental studies that constrain metallicity dependence, will provide decisive tests of these predictions and of the dominant role of binary interactions in shaping SESNe.

Low mass (sub)stellar objects represent the low end of the initial mass function, the transition to free-floating planets and a prominent interloper population in the search for high-redshift galaxies. Without proper motions or spectroscopy, can one identify these objects photometrically? JWST/NIRCam has several advantages over HST/WFC3 NIR: more filters, a greater wavelength range, and greater spatial resolution. Here, we present a catalogue of (sub)stellar dwarfs identified in the Cosmic Evolution Early Release Science Survey (CEERS). We identify 518 stellar objects down to $m_F200W \sim 28$ using half-light radius, a full three magnitudes deeper than typical HST/WFC3 images. A kNN nearest neighbour algorithm identifies and types these sources, using four HST/WFC3 and four NIRCam filters, trained on SpeX spectra of nearby brown dwarfs. The kNN with four neighbors classifies well within two subtypes: e.g M2$\pm$2 or T4$\pm$2, achieving $\sim$95% precision and recall. More granular typing results in worse metrics. In CEERS, we find 9 M8$\pm$2, 2 L6$\pm$2, 1 T4$\pm$2, and 15 T8$\pm$2. We compare the observed long wavelength NIRCam colours -- not used in the kNN -- to those expected for brown dwarf atmospheric models. The NIRCam F356W-F444W and F410M-F444W colours are redder by a magnitude for the type assigned by the kNN, hinting at a wider variety of atmospheres for these objects. We find a 300-350pc scale-height for M6$\pm$2 dwarfs plus a second structural component and a 150-200pc scale-height for T6$\pm$2 type dwarfs, consistent with literature values.

AT2022rze is a luminous, ambiguous transient located South-East of the geometric center of its host galaxy at redshift z = 0.08. The host appears to be formed by a merging galaxy system. The observed characteristics of AT2022rze are reminiscent of active galactic nuclei (AGN), tidal disruption events (TDEs), and superluminous supernovae (SLSNe). The transient reached a peak absolute magnitude of -20.2 +- 0.2 mag, showing a sharp rise (trise,1/e = 27.5 +- 0.6 days) followed by a slow decline (tdec,1/e = 382.9 +- 0.6). Its bumpy light curve and narrow Balmer lines indicate the presence of gas (and dust). Its light curve shows rather red colors, indicating that the transient could be affected by significant host extinction. The spectra reveal coronal lines, indicative of high-energy (X-ray/UV) emission. Archival data reveal no prior activity at this location, disfavoring a steady-state AGN, although an optical spectrum obtained prior to the transient is consistent with an AGN classification of the host. Based on this, we conclude that the transient most likely represents a Changing-look AGN at the center of the smallest component of the merging system.

The relationship between Type Ia supernovae (SNe Ia) and their host galaxy stellar masses is well documented. In particular, Hubble residuals display a luminosity shift based on host mass, known as the mass step, which is often used as an extra correction in the standardisation of SN Ia luminosities. Here we investigate Hubble residuals and the mass step in the context of Si II $\lambda 6355$ velocities, using 277 near-peak SNe Ia from ZTF DR2. We divide the sample into high-velocity (HV) and normal-velocity (NV) SNe Ia, separated at 12,000 km/s, resulting in 70 HV and 207 NV objects. We then examine links between Si II $\lambda$6355 velocities, light-curve stretch $x_{1}$, colour $c$, and host properties to explore potential environmental and/or progenitor-related effects. Although we only find a marginal difference between the Hubble residuals of HV and NV SNe Ia, the NV mass step is $0.149 \pm 0.024$ mag ($6.3\sigma$), while HV SNe Ia show $0.046 \pm 0.041$ mag ($1.1\sigma$), consistent with zero. The NV-HV mass-step difference is $\sim 2.2\sigma$. The clearest subtype difference is seen in central regions (d_{DLR} < 1), where NV SNe Ia show a strong mass step but HV SNe Ia none, yielding a $3.1-3.6\sigma$ difference. A host-colour step appears for both: NV $0.142 \pm 0.024$ mag ($5.9\sigma$) and HV $0.158 \pm 0.042$ mag ($3.8\sigma$). Overall, NV and HV colour steps are consistent. HV SNe Ia show modest ($\sim 2.5$-$3\sigma$) steps in outer regions (d_{DLR} > 1), while NV SNe show stronger environmental trends. Thus, NV SNe Ia appear more environmentally sensitive, especially in central, likely metal-rich and older regions, while HV SNe Ia show weaker, subset-dependent trends, and applying a universal mass-step correction could introduce biases. Refined classifications or environment-dependent factors may improve future cosmological analyses beyond standard $x_{1}$ and $c$ cuts.

In this paper, we continue our study on the evolution of black holes (BHs) that receive velocity kicks at the origin of their host star cluster potential. We now focus on BHs in rotating clusters that receive a range of kick velocities in different directions with respect to the rotation axis. We perform N-body simulations to calculate the trajectories of the kicked BHs and develop an analytic framework to study their motion as a function of the host cluster and the kick itself. Our simulations indicate that for a BH that is kicked outside of the cluster's core, as its orbit decays in a rotating cluster the BH will quickly gain angular momentum as it interacts with stars with high rotational frequencies. Once the BH decays to the point where its orbital frequency equals that of local stars, its orbit will be circular and dynamical friction becomes ineffective since local stars will have low relative velocities. After circularization, the BH's orbit decays on a longer timescale than if the host cluster was not rotating. Hence BHs in rotating clusters will have longer orbital decay times. The timescale for orbit circularization depends strongly on the cluster's rotation rate and the initial kick velocity, with kicked BHs in slowly rotating clusters being able to decay into the core before circularization occurs. The implication of the circularization phase is that the probability of a BH undergoing a tidal capture event increases, possibly aiding in the formation of binaries and high-mass BHs.

We present panchromatic observations and modeling of the Calcium-rich supernova 2019ehk in the star-forming galaxy M100 (d$\approx$16.2 Mpc) starting 10 hours after explosion and continuing for ~300 days. SN 2019ehk shows a double-peaked optical light curve peaking at $t = 3$ and $15$ days. The first peak is coincident with luminous, rapidly decaying $\textit{Swift}$-XRT discovered X-ray emission ($L_x\approx10^{41}~\rm{erg~s^{-1}}$ at 3 days; $L_x \propto t^{-3}$), and a Shane/Kast spectral detection of narrow H$\alpha$ and He II emission lines ($v \approx 500$ km/s) originating from pre-existent circumstellar material. We attribute this phenomenology to radiation from shock interaction with extended, dense material surrounding the progenitor star at r<10^{15} cm and the resulting cooling emission. We calculate a total CSM mass of $\sim$ $7\times10^{-3}$ $\rm{M_{\odot}}$ with particle density $n\approx10^{9}\,\rm{cm^{-3}}$. Radio observations indicate a significantly lower density n < 10^{4}\,\rm{cm^{-3}} at larger radii. The photometric and spectroscopic properties during the second light curve peak are consistent with those of Ca-rich transients (rise-time of $t_r =13.4\pm0.210$ days and a peak B-band magnitude of $M_B =-15.1\pm0.200$ mag). We find that SN 2019ehk synthesized $(3.1\pm0.11)\times10^{-2} ~ \rm{M_{\odot}}$ of ${}^{56}\textrm{Ni}$ and ejected $M_{\rm ej} = (0.72\pm 0.040)~\rm{M_{\odot}}$ total with a kinetic energy $E_{\rm k}=(1.8\pm0.10)\times10^{50}~\rm{erg}$. Finally, deep $\textit{HST}$ pre-explosion imaging at the SN site constrains the parameter space of viable stellar progenitors to massive stars in the lowest mass bin (~10 $\rm{M_{\odot}}$) in binaries that lost most of their He envelope or white dwarfs. The explosion and environment properties of SN 2019ehk further restrict the potential WD progenitor systems to low-mass hybrid HeCO WD + CO WD binaries.

Young star clusters are the most common birth-place of massive stars and are dynamically active environments. Here, we study the formation of black holes (BHs) and binary black holes (BBHs) in young star clusters, by means of 6000 N-body simulations coupled with binary population synthesis. We probe three different stellar metallicities (Z=0.02, 0.002 and 0.0002) and two initial density regimes (density at the half-mass radius $\rho_{\rm h}\ge{}3.4\times10^4$ and $\ge{1.5\times10^2}$ M$_\odot$ pc$^{-3}$ in dense and loose star clusters, respectively). Metal-poor clusters tend to form more massive BHs than metal-rich ones. We find $\sim{}6$, $\sim{}2$, and <1% of BHs with mass m_{\rm BH}>60 M$_\odot$ at Z=0.0002, 0.002 and 0.02, respectively. In metal-poor clusters, we form intermediate-mass BHs with mass up to $\sim{}320$ M$_\odot$. BBH mergers born via dynamical exchanges (exchanged BBHs) can be more massive than BBH mergers formed from binary evolution: the former (latter) reach total mass up to $\sim{}140$ M$_\odot$ ($\sim{}80$ M$_\odot$). The most massive BBH merger in our simulations has primary mass $\sim{}88$ M$_\odot$, inside the pair-instability mass gap, and a mass ratio of $\sim{}0.5$. Only BBHs born in young star clusters from metal-poor progenitors can match the masses of GW170729, the most massive event in O1 and O2, and those of GW190412, the first unequal-mass merger. We estimate a local BBH merger rate density $\sim{}110$ and $\sim{}55$ Gpc$^{-3}$ yr$^{-1}$, if we assume that all stars form in loose and dense star clusters, respectively.

(Abridged) Using the Zwicky Transient Facility alert stream, we are conducting a large campaign to spectroscopically classify all transients occurring in galaxies in the Census of the Local Universe (CLU) catalog. The aim of the experiment is to construct a spectroscopically complete, volume-limited sample of transients coincident within 100" of CLU galaxies out to 200 Mpc, and to a depth of 20 mag. We describe the survey design and spectroscopic completeness from the first 16 months of operations. We present results from a systematic search for Calcium rich gap transients in the sample of 22 low luminosity (peak absolute magnitude $M > -17$), hydrogen poor events found in the experiment (out of 754 spectroscopically classified SNe). We report the detection of eight Calcium rich gap transients, and constrain their volumetric rate to be at least $\approx 15\pm5$% of the SN Ia rate. Combining this sample with ten events from the literature, we find a likely continuum of spectroscopic properties ranging from events with SN Ia-like features (Ca-Ia objects) to SN Ib/c-like features (Ca-Ib/c objects) at peak light. Within the Ca-Ib/c events, we find two populations of events distinguished by their red ($g - r \approx 1.5$ mag) or green ($g - r \approx 0.5$ mag) spectral colors at $r$-band peak, wherein redder events show strong line blanketing signatures, slower light curves, weaker He lines and lower [Ca II]/[O I] in the nebular phase. Together, we find that the spectroscopic continuum, volumetric rates and striking old environments are consistent with the explosive burning of He shells on low mass white dwarfs. We posit that Ca-Ia and red Ca-Ib/c objects are consistent with the double detonation of He shells with high He burning efficiency, while green Ca-Ib/c objects could arise from less efficient He burning scenarios such as detonations in low density He shells or He shell deflagrations.

Modern time-domain surveys like the Zwicky Transient Facility (ZTF) and the Legacy Survey of Space and Time (LSST) generate hundreds of thousands to millions of alerts, demanding automatic, unified classification of transients and variable stars for efficient follow-up. We present AppleCiDEr (Applying Multimodal Learning to Classify Transient Detections Early), a novel framework that integrates four key data modalities (photometry, image cutouts, metadata, and spectra) to overcome limitations of single-modality classification approaches. Our architecture introduces (i) two transformer encoders for photometry, (ii) a multimodal convolutional neural network (CNN) with domain-specialized metadata towers and Mixture-of-Experts fusion for combining metadata and images, and (iii) a CNN for spectra classification. Training on ~ 30,000 real ZTF alerts, AppleCiDEr achieves high accuracy, allowing early identification and suggesting follow-up for rare transient spectra. The system provides the first unified framework for both transient and variable star classification using real observational data, with seamless integration into brokering pipelines, demonstrating readiness for the LSST era.

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.

Young star clusters are likely the most common birthplace of massive stars across cosmic time and influence the formation of compact binaries in several ways. Here, we simulate the formation of black hole -- neutron star binaries (BHNSs) in young star clusters, by means of the binary population synthesis code \texttt{MOBSE} interfaced with the $N$-body code \texttt{NBODY6++GPU}. BHNSs formed in young star clusters (dynamical BHNSs) are significantly more massive than BHNSs formed from isolated binaries (isolated BHNSs): $\sim{}40$~\% of the dynamical BHNS mergers have total mass $>15$ M$_\odot$, while only $\sim{}0.01$~\% of the isolated BHNS mergers have mass in excess of this value. Hence, our models strongly support a dynamical formation scenario for GW190814, given its total mass $\sim{}26$ M$_\odot$, if this event is a BHNS merger. All our dynamical BHNSs are ejected from their parent star cluster before they reach coalescence. Thus, a significant fraction of BHNS mergers occurring in the field might have originated in a young star cluster. The mass spectrum of BHNS mergers from gravitational-wave detections will provide a clue to differentiate between dynamical and isolated formation of BHNSs.

We present JWST/MIRI and complementary ground-based near-infrared observations of the Type II SN 2017eaw taken 6 years post-explosion. SN 2017eaw is still detected out to 25 $\mu$m and there is minimal evolution in the mid-infrared spectral energy distribution (SED) between the newly acquired JWST/MIRI observations and those taken a year earlier. Modeling of the mid-infrared SED reveals a cool $\sim$160 K dust component of $5.5\times10^{-4}\ \mathrm{M}_\odot$ and a hot $\sim$1700 K component of $5.4\times10^{-8}\ \mathrm{M}_\odot$ both composed of silicate dust. Notably there is no evidence of temperature or mass evolution in the cool dust component in the year between JWST observations. We also present new and archival HST and ground-based ultraviolet (UV) and optical observations which reveal reduced but continued circumstellar medium (CSM)-ejecta interaction at $>$2000 days post-explosion. The UV and mid-infrared emission show similar decline rates, suggesting both probe the interface between the ejecta and CSM. Given this, the continued existence of boxy H$\alpha$ emission in the nebular spectra, the low inferred optical depth of the dust, and the lack of temperature and mass evolution, we suggest that the cool dust component in SN 2017eaw may be primarily due to pre-existing dust rather than newly-formed dust in the ejecta or cold dense shell.

We present observations and modeling of SN 2016hnk, a Ca-rich supernova (SN) that is consistent with being the result of a He-shell double-detonation explosion of a C/O white dwarf. We find that SN 2016hnk is intrinsically red relative to typical thermonuclear SNe and has a relatively low peak luminosity ($M_B = -15.4$ mag), setting it apart from low-luminosity Type Ia supernovae (SNe Ia). SN 2016hnk has a fast-rising light curve that is consistent with other Ca-rich transients ($t_r = 15$ d). We determine that SN 2016hnk produced $0.03 \pm 0.01 M_{\odot}$ of ${}^{56}\textrm{Ni}$ and $0.9 \pm 0.3 M_{\odot}$ of ejecta. The photospheric spectra show strong, high-velocity Ca II absorption and significant line blanketing at \lambda < 5000 Angstroms, making it distinct from typical (SN 2005E-like) Ca-rich SNe. SN 2016hnk is remarkably similar to SN 2018byg, which was modeled as a He-shell double-detonation explosion. We demonstrate that the spectra and light curves of SN 2016hnk are well modeled by the detonation of a $0.02 \ M_{\odot}$ helium shell on the surface of a $0.85 \ M_{\odot}$ C/O white dwarf. This analysis highlights the second observed case of a He-shell double-detonation and suggests a specific thermonuclear explosion that is physically distinct from SNe that are defined simply by their low luminosities and strong [Ca II] emission.

Stellar streams -- formed from tidally stripped globular clusters or dwarf galaxies -- are sensitive tracers of a galaxy's accretion history and gravitational potential. While numerous streams are known in the Milky Way (MW), the formation and evolution of stellar streams have been primarily studied in isolated settings. The impact of subsequent galaxy interactions on stellar streams remains largely unexplored. Understanding merger-induced effects is however, crucial given the accretion of the Large Magellanic Cloud (LMC) onto the MW, and the fact that for example M31 and Cen A have experienced recent mergers. We analyze the detailed evolution of 1024 stellar streams during a complete MW--LMC-like merger, systematically varying initial stream properties and considering various orbits for the infalling perturber. We find that an MW--LMC mass-ratio merger significantly alters stellar stream properties, including energy, angular momentum, orbit, and morphology. Some streams exhibit dramatic morphological changes or develop complex substructures, while others see substantial shifts in energy and/or angular momentum, or re-orient their orbital plane. Interestingly, strong morphological alterations do not necessarily correlate with large changes in energy or orbit. A few streams split apart with parts moving to different orbits, appearing disconnected in position and kinematics despite their common origin. Strong effects correlate with close encounters between stream particles and the infalling perturber at various times during the merger. Our findings highlight the considerable impact of significant accretion events on the properties of stellar streams, and the challenge to recover the initial orbits of streams from their appearance at the present-day. Visualizations of the detailed evolution of all 1024 stellar streams are available at this https URL

Globular clusters may host intermediate mass black holes (IMBHs) at their centres. Here we propose a new method for their identification using millisecond pulsars (MSPs) as probes. We show that measuring the first (jerk) and second (jounce) derivatives of the accelerations of an ensemble of MSPs will let us infer the presence of an IMBH in a globular cluster better than measuring the sole accelerations. We test this concept by simulating a set of star clusters with and without a central IMBH to extract the distributions of the stellar jerks and jounces. We then apply this technique to the ensemble of MSPs in the Galactic globular cluster 47 Tucanae. Current timing observations are insufficient to constrain the presence of an IMBH and can only be used to pose upper limits on its mass. But, with few more years of observations it will be possible to test for the presence of a central IMBH with mass smaller than $\sim$ 1000 M$_{\odot}$. We conclude that jerks and jounces help significantly in reducing the upper limit of the mass of IMBHs in Galactic globular clusters.