Vineet Ojha, Xue-Bing Wu, Luis C. Ho

We performed a comprehensive analysis of flux and color variability in a redshift-matched sample of Seyfert galaxies, comprising 23 gamma-ray-detected narrow-line Seyfert 1 galaxies (gNLS1s), 190 non-gamma-ray-detected narrow-line Seyfert 1 galaxies (ngNLS1s), and 10 gamma-ray-detected broad-line Seyfert 1 galaxies (gBLS1s). Utilizing multi-band light curves from the Zwicky Transient Facility (ZTF) in g, r, and i bands, along with mid-infrared (MIR) observations in W1 and W2 bands from the Wide-Field Infrared Survey Explorer (WISE), we observed that gBLS1s exhibit more significant variability than gNLS1s, while ngNLS1s display minimal variability across both optical and MIR wavelengths. The pronounced variability in gBLS1s may be attributed to a more closely aligned jet relative to the observer's line of sight or their comparatively lower accretion rates. In contrast, the subdued variability in ngNLS1s suggests that their flux changes are primarily driven by accretion disk instabilities. A significant correlation between optical and MIR variability amplitudes found here supports the reprocessing scenario, wherein variations in the accretion disk emission are re-emitted by surrounding dust. Furthermore, our long-term color variability analysis revealed a stronger redder-when-brighter (RWB) trend in approximately 35%, 61%, and 80% of gNLS1s, ngNLS1s, and gBLS1s, respectively, in optical wavelength, strengthens the reprocessing scenario with the observed trend of stronger RWB for approximately 68%, 96%, and 30%, respectively, in MIR wavelength. The prevalent RWB trend observed in both optical and MIR wavelengths from the current sample of Seyfert galaxies on the longer time scales is likely associated with accretion disk instabilities.