NEXUS: A Search for Nuclear Variability with the First Two JWST NIRCam Epochs

The multi-cycle JWST Treasury program NEXUS will obtain cadenced imaging and spectroscopic observations around the North Ecliptic Pole during 2024-2028. Here we report a systematic search for nuclear variability among $\sim 25\,$k sources covered by NIRCam (F200W+F444W) imaging using the first two NEXUS epochs separated by 9 months in the observed frame. Difference imaging techniques reach $1\sigma$ variability sensitivity of 0.18~mag (F200W) and 0.15~mag (F444W) at 28th magnitude (within 0".2 diameter aperture), improved to $0.01$~mag and $0.02$~mag at $<25$th magnitude, demonstrating the superb performance of NIRCam photometry. The difference imaging results represent significant improvement over aperture photometry on individual epochs (by $>30\%$). We identify 465 high-confidence variable sources among the parent sample, with 2-epoch flux difference at $>3\sigma$ from the fiducial variability sensitivity. Essentially all these variable sources are of extragalactic origin based on preliminary photometric classifications, and follow a similar photometric redshift distribution as the parent sample up to $z_{\rm phot}>10$. While the majority of these variability candidates are likely normal unobscured AGNs, some of them may be rare nuclear stellar transients and tidal disruption events that await confirmation with spectroscopy and continued photometric monitoring. We also constrain the photometric variability of ten spectroscopically confirmed broad-line Little Red Dots (LRDs) at $3\lesssim z \lesssim 7$, and find none of them show detectable variability in either band. We derive stringent $3\sigma$ upper limits on the F444W variability of $\sim 3-10\%$ for these LRDs, with a median value of $\sim 5\%$. These constraints imply weak variability in the rest-frame optical continuum of LRDs.