## A closer look at dwarf galaxies exhibiting MIR variability: AGN confirmation and comparison with non-variable dwarf galaxies

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00666v1-figure-313b6fe3-8859-46ba-9820-4d9c463099ce.jpg)

1 Department of Physics and Astronomy, University of California, Riverside, 900 University Ave, Riverside CA 92521, USA 2 U.S. Naval Observatory, 3450 Massachusetts Avenue NW, Washington, DC 20392-5420, USA 3 George Mason University, Department of Physics and Astronomy, MS3F3, 4400 University Drive, Fairfax, VA 22030, USA 4 Hiroshima Astrophysical Science Center, Hiroshima University, 1-3-1 Kagamiyama, Higashi-Hiroshima, Hiroshima 739-8526, Japan

(Received May 21, 2024; Revised July 18, 2024; Accepted July 31, 2024)

## ABSTRACT

Detecting active black holes in dwarf galaxies has proven to be a challenge due to their small size and weak electromagnetic signatures. Mid-infrared variability has emerged as a promising tool that can be used to detect active low-mass black holes in dwarf galaxies. We analyzed 10.4 years of photometry from the ALL WISE /NEO WISE multi-epoch catalogs, identifying 25 objects with AGNlike variability. Independent confirmation of AGN activity was found in 68% of these objects using optical and near-infrared diagnostics. Notably, we discovered a near-infrared coronal line [S IX ] λ 1.252 µ m in J1205, the galaxy with the lowest stellar mass (log M ∗ = 7.5 M ⊙ ) and low metallicity (12 + log(O/H) = 7.46) in our sample. Additionally, we found broad Pa α potentially from the BLR in two targets, and their implied black hole masses are consistent with black hole-stellar mass relations. Comparing non-variable galaxies with similar stellar masses and WISE W 1 -W 2 colors, we found no clear trends between variability and large-scale galaxy properties. However, we found that AGN activity likely causes redder W 1 -W 2 colors in variable targets, while for the non-variable galaxies, the contribution stems from strong star formation activity. A high incidence of optical broad lines was also observed in variable targets. Our results suggest that mid-infrared variability is an effective method for detecting AGN activity in low-mass galaxies and can help uncover a larger sample of active low-mass ( &lt; 10 6 M ) black holes in the universe. ⊙

Keywords: Active galactic nuclei (16), Dwarf galaxies (416), Infrared galaxies (790)

## 1. INTRODUCTION

The formation of the first black holes (BHs) in the universe is still not well understood (Volonteri &amp; Begelman 2010; Natarajan 2014; Inayoshi et al. 2020). There are two most commonly discussed formation scenarios for black hole seeds, which eventually evolved into supermassive black holes (SMBHs) at the centers of nearly every massive galaxy. In one scenario, the first BHs were formed by the collapse of early Population III stars. The BHs formed by this method are likely to have masses in the range 100 1000 M ⊙ (Madau &amp; Rees 2001) and should undergo phases of accretion at super-Eddington rates in order to reach the masses of the SMBH that are currently present in the universe (Volonteri &amp; Rees 2005; Lupi et al. 2016). In an alternate scenario, the first BHs were formed by the direct collapse of gas in the early universe (Silk &amp; Rees 1998; Bromm &amp; Loeb 2003; Begelman et al. 2006), and have heavier masses of about 10 4 -10 5 M ⊙ (Lodato &amp; Natarajan 2006; Inayoshi et al. 2020; Volonteri et al. 2021; Bogd´n et al. 2024). a Other commonly suggested formation mechanisms include runaway collisions in dense star clusters (Portegies Zwart et al. 2004; Devecchi et al. 2012; Natarajan 2021), which can produce BH seeds with masses 10 3 - 10 5 M . ⊙

aarav005@ucr.edu

Direct observational evidence for the nature of the earliest BHs is challenging to obtain. However, observing

BHs in dwarf galaxies in the nearby universe could help narrow down the formation scenarios (Mezcua 2017; Reines 2022; Haidar et al. 2022; Beckmann et al. 2023). If dwarf galaxies have a higher occupation fraction (percentage of galaxies with BHs) of low-mass ( &lt; 10 3 M ) ⊙ BHs, it is likely that the BH seeds were formed from the collapse of Population III stars (van Wassenhove et al. 2010; Greene 2012). Conversely, if they have a lower occupation fraction of high-mass ( &gt; 10 5 M ) BHs, the ⊙ direct collapse scenario is favored (Miller et al. 2015; Nguyen et al. 2019; Greene et al. 2020).

Detecting active black holes, or active galactic nuclei (AGN), which are actively accreting material, is one of the best methods to determine the presence of BHs in galaxies. Although AGN detections cannot provide a complete census of the BH occupation fraction in dwarf galaxies (since only a small fraction of BHs are likely to be actively accreting and visible as AGN), identifying AGN can still offer constraints on the occupation fraction. Moreover, AGN detections can guide followup spectroscopic observations to determine BH masses using methods such as broad-line measurements, further refining our understanding of BH formation scenarios.

With advancements in telescope technology and the advent of large-scale surveys, it has become clear that a sizeable number of dwarf galaxies are likely to host AGN (Reines et al. 2013; Moran et al. 2014; Chilingarian et al. 2018; Mezcua et al. 2018; Mezcua &amp; Dom´ ınguez S´nchez 2020; Salehirad et al. 2022; Mezcua a &amp; Dom´ ınguez S´nchez 2024). a Results from JWST indicate the presence of BH in low-metallicity galaxies in the early universe (Kocevski et al. 2023; Kokorev et al. 2023; Ubler et al. 2023; Harikane et al. 2023; Maiolino ¨ et al. 2023; Goulding et al. 2023) and local dwarfs can be used as laboratories in which to study similar BH, and their effects on the host galaxy in much greater detail than is possible at high z (Bose et al. 2018). Additionally, the detection of AGN in dwarf galaxies can be used to constrain the BH masses and occupation fractions (Chadayammuri et al. 2023).

It is challenging to identify active BH in dwarf galaxies using traditional diagnostic methods such as the Baldwin, Phillips &amp; Terlevich (BPT) diagram (Baldwin et al. 1981) or Veilleux-Osterbrock (VO) diagrams (Veilleux &amp; Osterbrock 1987) as dwarf galaxies are more likely to have low gas-phase metallicities (Lequeux et al. 1979). Cann et al. (2019) show that, as the BH mass decreases, the hardening of the spectral energy distribution from the accretion disk changes the ionization structure of the nebula. The enhanced high-energy emission from IMBHs results in a more extended partially ionized zone compared with models for higher-mass BH. This effect produces a net decrease in the predicted [O III ]/H β and [N II ]/H α emission line ratios, and the positions of galaxies can shift from the AGN region to the starforming (SF) region. While diagrams that use [S II ] and [O I ] fluxes are useful for recovering AGN across a range of BH masses and galaxy metallicities (Cann et al. 2019; Polimera et al. 2022), these diagrams still miss a large fraction of AGN (Wasleske &amp; Baldassare 2024), most likely due to the strong contribution from star formation to the line ratios. The low-mass BH in dwarf galaxies also have a small gravitational sphere of influence compared to more massive BH, making them impossible to detect dynamically in galaxies apart from the nearby ones (Do et al. 2014; den Brok et al. 2015; Nguyen et al. 2018, 2019). A low-mass BH, accreting at the Eddington luminosity, will be of comparable luminosity to that produced by stars in the galaxy, also making it challenging to accurately detect AGNs in dwarf galaxies.

A promising method to locate AGNs in dwarf galaxies is by making use of variability. AGNs are known to be variable across the entire electromagnetic spectrum and in different time scales, ranging from months to days depending on the mass of the galaxy (Ulrich et al. 1997; Vanden Berk et al. 2004; Sesar et al. 2007; MacLeod et al. 2010; Burke et al. 2020, 2021). This variability is often attributed to instabilities in the accretion disk or surface temperature fluctuations (Ruan et al. 2014). Baldassare et al. (2020) made use of optical variability using the Sloan Digital Sky Survey (SDSS, York et al. 2000) and Palomar Transient Factory (PTF, Law et al. 2009) to identify AGNs in dwarf galaxies. They found that nearly 75% of low-mass galaxies in their sample with AGN-like variability have narrow emission lines dominated by star formation and would be completely missed by the BPT diagram. Similarly Mart´ ınezPalomera et al. (2020) used the High Cadence Transit Survey (HiTS, Mart´ ınez-Palomera et al. 2018) to look for optical variability and found that nearly 95% of variable AGNs in their sample do not exhibit any AGN-like optical emission line ratios. Ward et al. (2022) found 44 out of 25,714 dwarf galaxies had optically variable AGN candidates using optical photometry from the Zwicky Transient Facility (ZTF, Bellm et al. 2019). A large sample of optically variable AGN in dwarf galaxies was also obtained from the photometric analysis of various deep field surveys such as Dark Energy Survey (DES, Dark Energy Survey Collaboration et al. 2016; Burke et al. 2022) and COSMOS (Kimura et al. 2020), as well as predictions for possibly uncovering a larger sample from future surveys such as the Rubin Observatory (Ivezi´ c et al. 2019; Burke et al. 2023, 2024).

However, high levels of extinction can severely reduce the effectiveness of using optical variability to uncover AGN in low-mass galaxies. Mid-infrared (MIR) variability, on the other hand, is less sensitive to optically obscured and Compton thick AGN, which form a significant fraction of the low luminosity AGN population (Donley et al. 2008; Annuar et al. 2015, 2017; Ricci et al. 2016). Ultraviolet and optical emission from the accretion disk is reprocessed by the dust surrounding it and is re-radiated in the MIR. Thus, any fluctuations from the disk would be manifested in the MIR as well, leading to variability in the MIR (Kawaguchi &amp; Mori 2011; Son et al. 2022; Li &amp; Shen 2023; Son et al. 2023).

Secrest &amp; Satyapal (2020) found a low incidence of MIR variable dwarf galaxies (0.02%) using 8.4 years of photometry from the All WISE /NEO WISE multiepoch catalogs and concluded that high cadence data is necessary for MIR selection of AGN in dwarf galaxies. Ward et al. (2022) used forward-modeled MIR photometry of time-resolved WISE coadded images to uncover 148 MIR variable galaxies out of a sample of 79,879 dwarf galaxies. They also found that spectroscopic approaches to AGN identification would have missed 69% of their WISE MIR variable dwarfs.

While AGNs are variable in the MIR, it is important to realize that there are several other transients that are capable of varying in the infrared as well. Although MIR variability selection is less affected by supernovae (SN) contamination than optical variability selection due to the fainter SEDs of SN (Smitka 2016), there has been evidence of infrared variability caused by SN and variable stars (Karambelkar et al. 2019; Sun et al. 2022; Soraisam et al. 2023). There have also been instances of infrared transients from interactions between stars and planets (De et al. 2023), or from mergers between a pair of compact objects such as neutron stars or white dwarfs (Rosswog 2005; Arcavi et al. 2017). Although analyzing the evolution of the light curves with time can eliminate non-AGN sources such as rapid, one-off transient events, studies have shown that intermediate-mass black holes are capable of varying over timescales ranging from a few hours to a few days, making it difficult to distinguish them from more rapid transients (Shin et al. 2022).

Thus it is important to determine if MIR variability is accurately detecting AGN, specifically in dwarf galaxies. In this work, we aim to determine the effectiveness of MIR variability in tracing AGN activity by selecting a sample of dwarf galaxies that are variable in the MIR from Wide-field Infrared Survey Explorer ( WISE , Wright et al. 2010) light curves. We then obtained NIR spectroscopic observations of the subset of targets that do not show signs of AGN activity in the optical in order to look for other AGN indicators. These include the presence of NIR coronal lines which have proved to be a useful method to identify the presence of AGN in dwarf galaxies (Bohn et al. 2021) or broad lines from the broad line region (BLR) that might be too obscured to detect in the optical. We also analyzed how galaxies with variability compare to non-variable galaxies with similar WISE colors to study the differences, if any, between the two samples of galaxies. We do this to understand the effects of any biases in host galaxy properties on the observed MIR variability to indicate if MIR variability is preferentially detected in certain types of galaxies.

The paper is organized as follows: In Section 2, we describe the variable sample selection, present details of the targets, and describe the observations and data reduction. In Section 3, we present results of evidence of AGN activity in a subsample of MIR variable dwarfs. In Section 4, we compare the sample of variable dwarfs to a control sample of non-variable dwarf galaxies and analyze the differences between the properties of the two samples. We then summarize our conclusions in Section 5. Throughout the paper, we assume ΛCDM cosmology, with h = 0.73.

## 2. DATA AND OBSERVATIONS

## 2.1. Sample Selection

As in Secrest &amp; Satyapal (2020), we created a sample of dwarf galaxies from the NASA-Sloan Atlas (NSA), version nsa v1 0 1 . We started by matching the full NSA catalog to All WISE with 10 ′′ , which returned a unique match for 637168 out of 641409 (99%) of objects. We required that the elliptical Petrosian masses and S´ ersic masses differ by no more than 0.5 dex, which Secrest &amp; Satyapal (2020) found to effectively remove galaxies cut off in their SDSS images, leaving out 8739 objects.

Because we use catalog All WISE Multiepoch and NEO WISE -R PSF photometry for the light curves, we require that the galaxies not be extended or convolved with other sources in WISE , respectively, requiring ext flg&lt;=1 , nb==1 . This left 459000 galaxies. Finally, to place our sample on the WISE W 1 -W 2, W 2 -W 3 color-color diagram, we require that all objects have valid All WISE photometry in the first three bands, with corresponding cc flags equal to zero, leaving 316715 galaxies. Secrest &amp; Satyapal (2020) found that making a cut of All WISE W 2 &lt; 14 5 . mag removes the Eddington bias from the light curves, leaving 224969 systems. Of these, 6226 galaxies have elliptical Petrosian stellar masses less than 2 × 10 9 h -2 M ⊙ , which corresponds to log ( M /M ⋆ ⊙ ) &lt; . 9 6 using h = 0 73. .

4

## Table 1 . Summary of AGN confirmations in variable dwarf galaxies

## Observed with NIRES

AGN indicators other than MIR variability

Stellar Mass

Target

NIR Diagnostics

NIR BL

NIR Coronal lines

WISE

Optical variability

He II

BPT

Optical BL

(11)

(10)

(9)

(8)

(7)

(6)

(5)

(4)

(3)

(2)

(1)

· · ·

· · ·

· · ·

· · ·

✓

· · ·

· · ·

· · ·

(a)

✓

9.29

J000011.72+052317.4

· · ·

· · ·

· · ·

· · ·

(i)

✓

(g)

✓

· · ·

(e)

✓

(b,c)

✓

9.53

J082912.67+500652.3

✓

✓

· · ·

· · ·

(h)

✓

· · ·

· · ·

· · ·

· · ·

9.57

⋆

J084702.07+035202.0

· · ·

· · ·

· · ·

· · ·

· · ·

· · ·

· · ·

· · ·

(d)

✓

9.29

J085431.18+173730.5

· · ·

· · ·

· · ·

(j)

✓

✓

· · ·

✓

(e)

✓

· · ·

9.30

J090613.76+561015.3

✓

✓

✓

· · ·

(h)

✓

· · ·

· · ·

· · ·

· · ·

9.47

⋆

J092810.51+150228.0

†

· · ·

· · ·

· · ·

· · ·

· · ·

· · ·

· · ·

· · ·

8.29

J093608.61+061525.9

· · ·

· · ·

· · ·

· · ·

✓

· · ·

· · ·

· · ·

(a, b)

✓

9.52

J095438.79+403204.4

†

· · ·

· · ·

· · ·

· · ·

· · ·

✓

· · ·

· · ·

9.01

J102315.19+122227.9

✓

✓

· · ·

(h)

✓

· · ·

· · ·

· · ·

· · ·

8.87

⋆

J103206.01+225921.6

✓

✓

✓

· · ·

✓

· · ·

· · ·

· · ·

· · ·

9.51

J114433.55+072643.0

· · ·

· · ·

· · ·

· · ·

(i)

✓

· · ·

✓

(e)

✓

· · ·

9.10

J120325.67+330846.1

✓

· · ·

· · ·

✓

(h)

✓

· · ·

· · ·

· · ·

· · ·

7.52

⋆

J120503.54+455151.0

†

· · ·

· · ·

· · ·

· · ·

· · ·

· · ·

· · ·

· · ·

8.79

J122608.78-025226.2

· · ·

· · ·

· · ·

· · ·

· · ·

· · ·

(f)

✓

· · ·

· · ·

8.61

J125305.96-031258.7

· · ·

· · ·

· · ·

· · ·

✓

· · ·

(f)

✓

✓

· · ·

9.57

J130819.11+434525.6

✓

· · ·

· · ·

· · ·

✓

· · ·

· · ·

(e)

✓

· · ·

8.62

∗

J132053.67+215510.8

✓

· · ·

· · ·

· · ·

· · ·

· · ·

✓

· · ·

· · ·

9.34

J135844.73+123546.9

· · ·

· · ·

· · ·

· · ·

✓

· · ·

· · ·

· · ·

(c)

✓

8.93

J140630.09-001939.3

✓

· · ·

· · ·

· · ·

· · ·

· · ·

· · ·

· · ·

· · ·

9.42

J142835.84+570834.3

†

· · ·

· · ·

· · ·

· · ·

· · ·

(f)

✓

· · ·

· · ·

8.51

J143744.58+524333.3

· · ·

· · ·

· · ·

· · ·

✓

· · ·

(f)

✓

· · ·

· · ·

9.39

J152637.36+065941.6

· · ·

· · ·

· · ·

· · ·

· · ·

(g)

✓

· · ·

· · ·

✓

9.14

J154404.29+275334.0

✓

· · ·

· · ·

· · ·

· · ·

· · ·

· · ·

· · ·

· · ·

9.26

J161243.20+124505.2

· · ·

· · ·

· · ·

· · ·

✓

· · ·

· · ·

· · ·

(a, b)

✓

9.45

J233245.03-005845.8

## Table 1 continued

## Table 1 (continued)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00666v1-figure-47bf913d-c672-401f-81b8-f7d922972511.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00666v1-figure-8d01798d-42f6-4793-8e41-78b78614d006.jpg)

Figure 1. WISE light curves of J1205, one of the variable dwarf galaxies identified in this work, in different WISE bands W1 (3.4 µ m), W2 (4.6 µ m) and the change in color is depicted by W1W2. The single epoch photometry of the light curve is given by the grey points, while the blue points and the error bars represent the weighted mean and the standard error per 10-day observing period. The subtle redder-when-brighter behavior, consistent with the brightening being due to increasing AGN dominance over the host galaxy, and the gradual nature of the light curves indicate that the variability is likely due to AGN.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00666v1-figure-0dcb9636-ed8a-4296-bd3a-fff32a27eff6.jpg)

Single-epoch photometry from the All WISE Multiepoch Photometry (MEP) and NEO WISE tables for these 6226 dwarf galaxies were subsequently queried using the All WISE source positions, using zero position offset for MEP and a 3 ′′ tolerance for the NEO WISE positions. The single-epoch photometry was likewise cleaned, requiring valid measurements in W 1 and W 2, allowing a maximum deviation of ± 1 mag between the single-epoch W 2 and All WISE W 2 corresponding to the limit set by the Eddington bias, na==0 and nb==1 for uncomplicated PSF photometry, W 1 and W 2 cc flags and moon masked = 0, qi fact &gt; 0 and saa sep &gt;= 5 . This removed 11% of single-epoch measurements, retaining cleaned data for all 6226 dwarf galaxies. On average, there are 276 single-epoch photometric measurements for each source, from which variability metrics were calculated as in Secrest &amp; Satyapal (2020). In that work, we found that Pearson r between W 1 and W 2 is the most useful initial discriminant for genuine variability, with r &gt; 0 4 being optimal for . WISE data. This produced an initial sample of 60 objects, allowing manual inspection of their light curves. We further classified the light curves that exhibit clear variability into those with rapid AGN-like variability and those with supernova-like flares or slow declines in apparent magnitude.

Although we started with a sample of dwarf galaxies, cutting on objects with AGN-like photometric variability creates a strong selection bias on the otherwise small fraction of contaminating background quasars and galaxies with incorrect spectroscopic redshifts. To mitigate this, we uploaded the WISE coordinates of the variable dwarf galaxy candidates to the SDSS thumbnail server and removed clear instances of nearby background quasars or objects with spectra not matching their catalog redshifts. After removing these contaminants, we were left with a final sample of 25 dwarf galaxies with bona fide, AGN-like, MIR variability (see Fig. 1 for an example). As a check, we estimated the number of background quasars that could be contaminating our sample of 25 objects. Using the SDSS/BOSS quasar catalog (Lyke et al. 2020), 98% of which have WISE photometry (738088 objects), we find that 114119 have W 2 &lt; 14 5 mag like our sample. . As this catalog covers ∼ 28% of the sky, the corresponding sky density of SDSS/BOSS quasars is ∼ 10 deg -2 . However, because our AGN selection criteria is mid-IR based, we correct for the obscured population that may have been largely missed by the SDSS/BOSS surveys. We conservatively assume that about half of the mid-IR AGN-selected sources are heavily obscured (e.g., Petter et al. 2023), giving a sky density of 20 deg -2 . The final sample of 25 objects has a maximum NSA/All WISE position offset of 4 . ′′ 4, of which 6023 of the 6226 dwarf galaxies satisfy. Then, the expected number of contaminant background QSOs is ∼ 0 6. . Thus, even if mid-IR variability selects all background QSOs, it is likely that there is at most one contaminating background QSO in our final sample of 25 dwarf galaxies. The final sample is presented in Table 1.

The stellar masses listed in Table 1 were obtained from the NSA catalog (Blanton &amp; Roweis 2007). In order to estimate the uncertainty in these masses, we compared the masses listed by NSA with masses reported in other catalogs. We found an average deviation of 0.01 dex between the NSA masses and those reported in the MPA-JHU catalog (Kauffmann et al. 2003a; Salim et al. 2007). This falls in the range of the expected deviations (-0.2 dex - 0.2 dex) between the MPA-JHU and NSA stellar masses determined by Pace et al. (2019). We also determined an average deviation of 0.3 dex between the NSA masses and those reported by the Portsmouth group (Maraston et al. 2006) for galaxies with masses determined in all three catalogs. These values provide

Table 2. Observing Log

| Galaxy (SDSS Name)    | UT Date Obs. (YYYY MM DD)   |   Redshift |   Stellar Mass (log M ⊙ ) | Exp. Time a (s)   |   Slit PA (degrees) |   Airmass | Telluric   | Flux variation b (%)   |
|-----------------------|-----------------------------|------------|---------------------------|-------------------|---------------------|-----------|------------|------------------------|
| J084702.07+035202.0   | 2023 04 07                  |      0.059 |                      9.57 | 8 × 240           |                 306 |      1.05 | HD 85504   | 15                     |
| J092810.51+150228.0   | 2023 04 08                  |      0.078 |                      9.47 | 4 × 240           |                  45 |      1.5  | HD 89239   | 25                     |
| J103206.01+225921.6   | 2023 04 08                  |      0.058 |                      8.87 | 4 × 240           |                   0 |      1.5  | HD 101060  | 50                     |
| J114433.55+072643.0   | 2023 04 07                  |      0.066 |                      9.51 | 4 × 240           |                 322 |      1.3  | HD 95126   | 15                     |
| J120503.54+455151.0 c | 2023 04 07                  |      0.065 |                      7.52 | 6 × 240           |                   0 |      1.17 | HD 95126   | 15                     |
| J132053.67+215510.8   | 2023 05 08                  |      0.022 |                      8.62 | 3 × 240           |                  30 |      1.28 | HD 116960  | < 10                   |
| J135844.73+123546.9   | 2023 05 08                  |      0.024 |                      9.34 | 4 × 240           |                  90 |      1.24 | HD 122945  | < 10                   |
| J142835.84+570834.3   | 2023 05 08                  |      0.059 |                      9.42 | 6 × 240           |                 344 |      1.28 | HD 116405  | < 10                   |
| J161243.20+124505.2   | 2023 05 08                  |      0.034 |                      9.26 | 3 × 240           |                   4 |      1.46 | HD 210501  | < 10                   |

a Exposures were typically done in ABBA nodding.

b Variation in flux between individual exposures.

c Indicates galaxies with coronal line detections. See Section 3.1

a measure of the average uncertainty in the masses of our sample. Thus, the stellar masses that we use fall within the range of expected masses of dwarf galaxies even with the predicted uncertainty. Also see Fig. 2 for their positions on the optical BPT diagram.

We then extensively searched to see if any variable dwarf galaxies had already been identified as AGN in the literature. Twelve of the 25 (48%) dwarf galaxies were already identified as AGN by a combination of other optical methods, such as the presence of broad Balmer lines, strong He II emission, optical variability, and emission line ratios indicative of AGN from the BPT diagram (see Table 1 for references).

Although the MIR variable galaxy J1320 has BPT line ratios indicative of AGN, it was also classified as a potential SN candidate in the literature (Izotov &amp; Thuan 2009; Reines et al. 2013), so we do not include it among the targets for which AGN activity was previously confirmed. Additionally, only one of the 25 (J1526) variable galaxies was classified as AGN from X-ray observations (Baldassare et al. 2017). Apart from the references mentioned in the caption of Table 1, we also cross-referenced the variable targets in this work with Moran et al. (2014); Polimera et al. (2022) for BPT AGN, Chilingarian et al. (2018) for optical broad lines, and Bi et al. (2020); Ohlson et al. (2024) for X-ray data, but we did not find any matches.

## 2.2. Observations and Reductions

We selected a subsample of variable dwarf galaxies for observations in the NIR. The 13 galaxies in the subsample include nine galaxies with no confirmations of

AGNby any method (including J1320) and four galaxies that only had WISE colors indicative of AGN activity. These were the only galaxies in the sample exhibiting MIR colors indicative of AGN without any other indicators across the criteria examined. We included the latter four galaxies because young starburst galaxies can often mimic the mid-infrared colors of AGN galaxies (Hainline et al. 2016) and may not truly be AGNs. Probing the NIR spectra of the variable targets could uncover evidence of AGN activity, such as broad Paschen lines from the BLR or coronal lines indicative of a hard radiation field.

NIR spectroscopy was obtained on three separate dates (UT 2023 April 7, UT 2023 April 8 and UT 2023 May 8) with the Keck II Near-Infrared Echellette Spectrometer (NIRES; Wilson et al. 2004). NIRES is a NIR echellette spectrograph with a wavelength coverage across five bands (Z, Y, J, H, K) from 0.94 to 2.45 µ m. There is a small gap in coverage between 1.85 to 1.88 µ m, but this region is dominated by water absorption from the atmosphere and thus has low atmospheric transmission. The average spectral resolution of the five orders ranges between 84 and 89 km s -1 (R = 3400). Observations on UT 2023 April 7 and UT 2023 April 8 were taken under variable and heavy cloud cover, while observations on UT 2023 May 8 were taken through light cirrus. Due to unfavorable weather conditions, we could get data for only nine of the 13 targets. Individual exposures for each target were taken for 4 minutes using the standard ABBA nodding, amounting to total exposure times ranging from 12 to 32 minutes, depending on the brightness of the target. A0 spectral class stars were

Figure 2. Optical BPT and VO diagrams showing the positions of the variable targets observed with NIRES (blue) and the variable galaxies not observed with NIRES (cyan). The emission line ratios were derived from fluxes obtained from the SDSS MPA-JHU catalog. Most of the observed variable targets lie primarily in the SF regions, and there is no indicator of AGN activity based on these diagnostic diagrams.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00666v1-figure-7d380fd6-9295-4664-ab87-1b07e0c08db3.jpg)

used as telluric standard stars, with measured magnitudes in J, H, and K bands. Observations of the telluric stars were taken before or after each science target to correct for the atmospheric absorption features. A summary of the observations is shown in Table 2.

The data were reduced using two modified pipelines. The flat fielding and background subtraction were done by the first pipeline using techniques described in Kelson (2003) and Becker et al. (2009). Rectification, telluric correction, wavelength calibration, and extraction were done with a slightly modified version of REDSPEC 1 . The telluric correction was done by dividing each science target spectrum by the spectrum of its corresponding telluric star and multiplying it by a blackbody curve of the same temperature. Wavelength calibration was done using strong OH lines, and the 1D spectra were then median combined. Flux calibration of individual targets was done by converting the magnitude of the telluric star to the associated flux in that band 2 . A minor corrective factor ( &lt; 5%) was introduced due to the differences between the center of the NIR bands and the center of the wavelength coverage. Because we observed through clouds, we found 15 - 50 % variations in the flux levels between individual exposures for the first two nights and less than 10 % variation in flux on the final night (Flux variations for individual objects

1 https://www2.keck.hawaii.edu/inst/nirspec/redspec.html

2 https://irsa.ipac.caltech.edu/data/SPITZER/docs/ dataanalysistools/tools/pet/magtojy/

are listed in Table 2). This introduced uncertainty in the flux calibration, as the latter was done relative to telluric standards observed immediately before or after each science target. We estimated that the uncertainties are of the order of the flux variations, and we included these added uncertainties in the measurements of fluxes.

## 2.3. Spectral fitting

All NIR spectra were fit using the Bayesian AGN Decomposition Analysis for SDSS Spectra (BADASS, Sexton et al. 2021) code, which was modified to cover the NIR wavelength regime. A narrow component for each emission line was fit simultaneously with a secondorder polynomial for the continuum. Absorption features were fit using the penalized pixel fitting method (pPXF, Cappellari &amp; Emsellem (2004)) with templates from the eMILES (Vazdekis et al. 2016) stellar library.

The error (noise) in the fit was determined by calculating the standard deviation of the continuum flux, represented by the grey-shaded region in the lower panel of Figure 3 (b). The σ residuals , also indicated in the residual plot (lower panel) of Figure 3(b) denotes the standard deviation of the residuals obtained by subtracting the modeled spectrum from the observed data.

## 3. EVIDENCE OF AGN ACTIVITY

## 3.1. Detection of coronal lines in J1205

Coronal lines (CL) are high ionization fine-structure lines that arise from collisionally excited forbidden finestructure transitions in highly ionized species, with ionization potentials that extend well beyond the Lyman

)

b

Figure 3. (a) NIR spectrum of J1205 in order 5 of NIRES, shifted to rest-frame wavelength. The blue-dashed line indicates the coronal line [S IX ] detected at a 5 σ level. The grey-shaded area represents the region of significant telluric absorption. (b) The zoomed-in spectrum with spectral fits to the [S IX ] (1.252 µ m) line and the [Fe II ] (1.257 µ m)line. The lower panel depicts the residuals to the fit, with the 1 σ error in the spectrum shaded in grey around the dashed line indicating 0. The σ residuals indicates the standard deviation of the residuals determined by subtracting the fit to the spectrum from the data.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00666v1-figure-479065ce-95df-43bd-83a4-9af8127fd90e.jpg)

limit to several hundred electron volts (eV). They are typically attributed to AGN activity and have been successfully detected in dwarf galaxies hosting AGN (e.g., Bohn et al. 2021) We looked for NIR coronal lines in the nine galaxies observed with NIRES.

The NIR coronal line [S IX ] λ 1.252 µ m (I.P. 328.2 eV and critical density log (n e )=9.4 cm -3 ) was detected in J1205. The fit to the spectrum showing the [S IX ] is shown in Fig. 3. The [S IX ] line is a 5 σ detection, and the measured flux of the line is 2.11 ± 0.52 × 10 -16 erg cm -2 s -1 . The full width at half maximum (FWHM) of the line is 138 ± 10 km s -1 , which is greater than the spectral resolution of NIRES in that order (85 km s -1 ) and thus is unlikely to be a noise spike or a cosmic ray signal. The errors in the flux and FWHM are given by the standard deviation of the posterior distribution from the MCMC fitting of the lines. Based on the commonly observed NIR emission lines in galaxies (Riffel et al. 2006; Landt et al. 2008; Martins et al. 2013), no other emission lines are likely to be found at that wavelength.

The detection of [S IX ] in J1205 is tantalizing evidence that it likely hosts an AGN. This is the first instance of a NIR coronal line being discovered in J1205, providing strong evidence that J1205 has an active black hole. J1205 has the lowest stellar mass and metallicity among all the targets (log(M ∗ ) = 7.5 M ⊙ and 12 + log(O/H) = 7.46; Izotov et al. 2021). It was classified as a compact starforming galaxy with extremely high [O III ]/[O II ] flux ratios by Izotov et al. (2017) and shown to have a

[Ne V ] λ 3426 ˚ A detection and evidence of hard ionization by Izotov et al. (2021). Spectroscopic fits to the H α emission line in spectra obtained from the Large Binocular Telescope (LBT) indicate the presence of a broad line feature, which is attributed to expansion motions of supernova remnants (Izotov et al. 2017). Harish et al. (2023) classified J1205 as a green pea galaxy possibly hosting an AGN based on MIR colors and variability from WISE . The low metallicity and stellar mass of J1205 would have made it impossible to detect the presence of the AGN by traditional methods such as the BPT diagram. This indicates that MIR variability can be a promising tool and, sometimes, the only method to detect AGN in low-mass galaxies.

)

Figure 4. Keck/NIRES spectrum of J1205 in the wavelength region around [Si VI ]. The small emission feature found at the wavelength of [Si VI ] (indicated by the dashed blue line) has a flux compatible with the expected flux for [Si VI ], however, it is below the 3 σ detection limit and may not be real.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00666v1-figure-5dd66187-2867-472f-9cb5-1f6b91b5253c.jpg)

While the presence of [S IX ] in J1205 is likely to indicate AGN activity, it could also possibly be the result of a tidal disruption event (TDE; Hinkle et al. 2023; Clark et al. 2024). This can be determined by follow-up observations to check if the coronal line persists or declines after an interval of time. However, a TDE would still indicate the presence of a black hole in J1205.

We did not detect any other coronal lines apart from [S IX ] in J1205. It is possible that circumnuclear stellar populations dominate the NIR continuum and thus drown out the CL emission, making fainter lines more difficult to detect (Rodr´ ıguez-Ardila et al. 2011), especially since J1205 has been classified as a compact starforming galaxy with very young ages ( &lt; 3 Myr) for star formation bursts. There is a feature near the noise limit that is suggestive of [Si VI ] (see Fig. 4), but the feature did not have sufficient signal to be fit using BADASS. We independently fit a Gaussian to this feature and determine the flux to be 1.78 ± 0.88 × 10 -17 erg cm -2

s -1 , which is below the 3 σ limit (see Fig. 5). Thus, this feature does not appear to be a 3 σ detection, but it is possible that with deeper observations, this feature can be detected with greater confidence.

While the detection of the coronal line [S IX ] alone and the absence of [Si VI ], the most commonly observed coronal line in galaxies (Rodr´ ıguez-Ardila et al. 2011; Lamperti et al. 2017; M¨ller-S´nchez et al. 2018; Bohn u a et al. 2021), may seem surprising, such occurrences are not unprecedented. There have been a few instances of galaxies that have a [S IX ] detection, but only upper limits to the [Si VI] flux are reported (Lamperti et al. 2017; den Brok et al. 2022). Comparing the [Si VI ] to [S IX ] flux ratios recorded in the literature (Rodr´ ıguez-Ardila et al. 2011; Lamperti et al. 2017; Bohn et al. 2021; den Brok et al. 2022), which range from 0.1 to 8 (more than half the galaxies have ratios less than 3.5), our measured ratio of 0.25 (assuming the 3 σ flux limit for [Si VI ]) falls within this observed range. It is possible that we simply did not obtain sufficiently deep observations to detect the [Si VI ] line in J1205.

No coronal lines were detected in the remaining eight targets observed with NIRES, potentially due to the relatively low S/N of the data. To investigate this, we plotted the upper limits to the [Si VI ] flux for each target versus the WISE W2 (4.6 µ m) flux in Fig. 5. We determined the upper limits to the [Si VI ] fluxes in the non-detections by estimating the 1 σ noise levels in the region around [Si VI ] flux and multiplying it by three to obtain the 3 σ upper limit. The solid line in Fig. 5indicates the best-fit relation, which was derived based on observations from M¨ller-S´nchez et al. (2018) (inu a dicated as black dots) by Cann et al. (2021). The plot indicates that J1428 is the only target for which [Si VI ] emission is highly unlikely to be present, as its upper limits are comfortably below the scatter of the relation, and it would have been detected if it were present. For the rest of the targets, we would likely require deeper observations in order to detect [Si VI ]. The absence of CLs in these targets is unlikely to be due to high extinction, as J1205 (the target with CL detection) does not have significantly lower extinction (calculated from the H α over H β Balmer decrement) than the rest of the observed targets.

While the presence of CLs in a galaxy's spectrum confirms the existence of a black hole, the absence of CLs does not mean the contrary. Bohn et al. (2021) analyzed a sample of 9 dwarf galaxies with optical and NIR line ratios indicative of AGN and found NIR CLs only in 5 of them (55 %). Several other works in the literature have found less than 100% CL detections in samples of Sy1 and Sy2 galaxies with similar spectral windows: 66%

Figure 5. Plot of [Si VI ] vs. W2 flux for the [Si VI ] nondetections in our sample. The solid line indicating the correlation is from Cann et al. (2021), which was derived from observations from M¨ uller-S´nchez et al. (2018) and the shaded region indicates a the scatter in the data. The targets are marked as J0847 (1), J0928 (2), J1032 (3), J1144 (4), J1205 (5), J1320 (6), J1358 (7), J1428 (8), J1612 (9). Targets whose upper limits are above the relation would require deeper observations in order to detect the [Si VI ] line, while for the targets with upper limits below the relation, it is likely that the [Si VI ] line is not present. For more details regarding the upper limits, refer to the appendix

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00666v1-figure-5b3cbeda-514f-493b-8bd6-efd4d3f950ad.jpg)

(36 out of 54, Rodr´ ıguez-Ardila et al. (2011)), 25%, (5 out of 20, Mason et al. (2015)), and 43% (44 out of 102, Lamperti et al. (2017)).

The relatively low detection rates in all of these samples could be an observational effect, as ground-based observations are always hampered by telluric absorption from the atmosphere. Depending on a target's redshift, the CLs can be shifted into regions of strong telluric absorption and/or in regions dominated by strong OH skylines, making it difficult to detect weak CLs. Other factors include a dominant contribution from the stellar population to the NIR continuum and the presence of stellar absorption features.

## 3.2. Detection of broad Pa α lines

Two of the variable targets (J0928 and J1144) show a broad line component in Pa α at 1.876 µ m. We justify the need for a second component based on the Ftest. The F-test determines the significance of the fit between a complex (higher-order) and a simple (lower-

Table 3. Black hole masses

|       | L Pa α (erg s - 1 )   | FWHM (km s - 1 )   | log 10 (M BH ) (M ⊙ )   |   log 10 (M stellar ) (M ⊙ ) |
|-------|-----------------------|--------------------|-------------------------|------------------------------|
| J0928 | 10 42 . 5 ± 0 . 2     | 325 ± 5            | 6.24 ± 0.27             |                         9.47 |
| J1144 | 10 41 . 2 ± 0 . 12    | 457 ± 5            | 5.96 ± 0.14             |                         9.51 |

order) model in order to justify the number of Gaussian components we include in our fit. To do this, we calculated the standard deviation of the residuals of the fit (given by σ ). We performed the F- test, given as F = ( σ lower -order ) 2 /( σ higher -order ) 2 , for one versus two components for the Gaussian fits to the Pa α line. The Ftest was calculated over a wavelength range of 1.86-1.89 µ m, which includes the Pa α λ 1.876 µ m. An F value greater than 3 indicates that the higher-order fit (i.e., the fit using more components) is justifiable. Using this method, we determined that the Pa α emission lines in J0928 and J1144 both required a second Gaussian component, which we allowed to freely vary to fit the broad feature present in the emission lines using BADASS. The fits to the spectrum can be seen in Fig. 6.

Both J0928 and J1144 show a broad component to Pa α while not showing a broad component in the optical Balmer lines (at least at the spatial resolution of the SDSS spectrum, observed with a fiber that includes significant contaminating/diluting stellar light). In galaxies with strong dust obscuration, it is possible that the broad line region (BLR) does not manifest itself in optical emission lines. However, NIR wavelengths are less affected by dust extinction, making it easier to observe broad components that would be obscured at optical wavelengths. This method has been used to uncover optically hidden AGN in massive galaxies (Veilleux et al. 1997; Lamperti et al. 2017; Bohn et al. 2020). There were no indicators of a broad component being present in any of the forbidden lines in the NIR for the two targets. Combined with the fact that the SDSS spectrum of the two targets also did not have a significant broad component in the forbidden [O III ] lines, it is unlikely that the broad component is produced by an outflow.

Assuming the broad lines measured are from the BLR, we used the FWHM and luminosity of the broad Pa α component to estimate black hole masses. Using virial mass estimators from Kim et al. (2018), we determined the black hole mass (Table 3) from the following equation:

$$\stackrel { \vdash } { \vdash } & & M & = 1 0 ^ { 7. 0 7 } \left ( \frac { L _ { P a \alpha } } { 1 0 ^ { 4 2 } \, \text{erg} \, s ^ { - 1 } } \right ) ^ { 0. 4 9 } \times \left ( \frac { F \text{WHM} } { 1 0 ^ { 3 } \, \text{km} \, s ^ { - 1 } } \right ) ^ { 2 } \\ \stackrel { \vdash } { \vdash } & &$$

Figure 6. Broad-line fits to the Pa α line in two variable targets, J0928 (a) and J1144 (b), which could indicate a broad component from the BLR of the galaxies. There was an unidentified line next to He I line at 1.869 µ m, which was fit but left unmarked. The lower panels in both figures depict the residuals to the fit, with the 1 σ error from the noise in the spectrum shaded in grey around the dashed line indicating 0.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00666v1-figure-2a97dc19-f24b-4907-9578-57f5e3a0a40e.jpg)

The errors given in Table 3 are the standard deviation of the posterior distribution from the MCMC fitting of the lines. The measured black hole masses agree with the predicted black hole masses based on the stellar masses of J1144 and J0928 (Fig. 7) from Reines &amp; Volonteri (2015).

## 3.3. NIR diagnostic diagrams

Optical diagnostic diagrams like the BPT diagram have been well studied and used extensively to distin- guish between AGN and SF regions in galaxies. However, similar diagnostic diagrams in the NIR are scarce, chiefly due to a lack of large spectroscopic surveys in the infrared. Nevertheless, there have been several attempts to recreate similar diagnostic diagrams in the infrared that can be used to distinguish between AGN and SF regions (Larkin et al. 1998; Smith et al. 2014; Colina et al. 2015). In particular, Calabr` et al. (2023) used o CLOUDY photoionization models to identify AGN - SF diagnostics based on the ratio of bright near-infrared

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00666v1-figure-4fcf84f4-998a-44d6-9126-12f2eca8124c.jpg)

∗

/circledot

Figure 7. Black hole vs. stellar mass for a sample of galaxies at low-redshift. The yellow points are galaxies at a redshift &lt; 0.06 and include dwarf galaxies with broad lines from Reines &amp; Volonteri (2015) (RV15). The solid line is the expected linear regression from Equations 4 and 5 in RV15, along with the 0.5 dex error, indicated by the gray shaded region. The two blue points indicate the calculated BH masses based on the broad Pa α for J0928 and J1144, respectively, which fits well with the expected black hole mass for the two galaxies.

emission lines such as [S III ] 9530 ˚ , [C A I ] 9850 ˚ , [P A II ] 11880 ˚ , A [Fe II ] 12570 ˚ , A and [Fe II ] 16470 ˚ to the A Paschen lines. They used a sample of 130 AGN and SF galaxies with redshifts 0 &lt; z &lt; 3 observed using JWSTNIRSpec (Jakobsen et al. 2022) as part of the Cosmic Evolution Early Release Science Survey (CEERS; PID 1345, PI Finkelstein, Finkelstein et al. 2023) and found that this diagram is able to successfully distinguish between AGN and SF galaxies. Furthermore, they found 60% more AGN using NIR diagnostic diagrams than when using optical BPT diagrams. Thus NIR diagnostic diagrams are able to uncover optically 'hidden' AGNs, thus highlighting the importance of using NIR wavelengths to potentially locate AGN in dusty galaxies. Since optical emission lines are significantly affected by extinction, it is possible to overlook AGN indicators from optical emission lines due to extinction. On the other hand, NIR emission lines are less sensitive to dust and thus could prove to be more helpful in identifying AGN activity in dusty galaxies.

Most of the MIR variable targets in our sample fall in the SF region of the optical BPT diagram as seen in Fig. 2. We plotted them in the various NIR diagnostic diagrams in Fig. 8. We primarily used the ratios of [S III ] λ 9530 ˚ /Pa A γ and [S III ] λ 9530 ˚ /Pa A β to [Fe II ] λ 16470 ˚ /Pa A α and [Fe II ] λ 12570 ˚ /Pa A β , obtained from NIRES. For four of the targets (marked as filled circles in Fig. 8), these were the lines with signal-to-noise ratios (S/N) &gt; 3 in our sample, and their flux values could be estimated with minimum uncertainties. For the remaining targets, we could only determine one of the ratios as at least one of the emission lines on the x-axis ([Fe II ] 1.64 µ m, [Fe II ] 1.25 µ m, Pa α or Pa β ) fell in the telluric region or were in a noisy part of the spectrum, and thus we could not measure the fluxes. We instead find the upper limits (if the [Fe II ] lines were undetermined) or lower limits (if the Paschen lines were undetermined) to the emission line ratios and showed the targets as open circles with arrows to indicate upper or lower limits (both if neither of the line fluxes used in the ratio could be determined) to the ratios in Fig. 8. The upper limits for J1032, J1428, and J1320 did not exclude the possibility that these targets are AGN.

We found that four of the variable targets (J0847, J0928, J1144, and J1358) fell in the composite and AGN regions across all four NIR diagnostic diagrams. J0928, which has a broad Pa α line, moved from the SF region in the optical to the AGN region of two of the NIR diagnostic diagrams. The other broad line target, J1144, also moved into the composite region of all four NIR diagnostic diagrams.

While NIR diagnostic diagrams appear to be more effective in some cases at identifying obscured AGN in dwarf galaxies, they are still limited by the gas-phase metallicities of the galaxies. J1205, with a coronal line indicating the presence of a black hole, still lies firmly in the SF region of all the NIR diagnostic diagrams.

## 3.4. Ruling out supernova contribution to variability

Although the nine galaxies that were observed using NIRES are determined to exhibit AGN-like variability from the light curves, our selection method does not completely rule out contributions from SN. Five of the galaxies have alternate evidence of AGN activity based on the presence of coronal lines, broad Pa α lines, and emission line ratios indicative of AGN, which makes it unlikely that there is a significant contribution from SN to the variability. However, Baldassare et al. (2016) showed that in the optical, broad H α lines could be caused by SN activity. The ideal method to rule out the possibility of broad lines being caused by SN is to obtain multi-epoch spectroscopy that would show whether

AGN min limit(Calabro+23)

Figure 8. NIR diagnostic diagrams using AGN and SF limits from Calabr` et al. (2023). The targets are marked as J0847 (1), J0928 (2), o J1032 (3), J1144 (4), J1205 (5), J1320 (6), J1358 (7), J1428 (8), J1612 (9). For the targets where some of the line fluxes were not available, we plotted them as open circles with arrows, indicating upper and lower limits. Dwarf Galaxies with AGN from Bohn et al. (2021) and Type 2 AGN from Riffel et al. (2006) are plotted in red and green respectively. The gas-phase metallicity (Z) increases from bottom to top, and the ionization parameter (log (U)) increases from right to left. Four of the targets moved to the composite/AGN region across all four NIR diagnostic diagrams, as compared to their positions in the SF region in the optical diagnostic diagrams.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00666v1-figure-de0afcd7-36c6-46ce-b2aa-070257a03bbb.jpg)

the broad emission lines persist. Although we do not have multi-epoch spectra and cannot check for line variability, we did inspect all the light curves by eye, and determined that they do not show the rapid increase and decline in their curves that is characteristic of supernovae. Furthermore, we inspected the NIR spectrum for the common NIR emission lines associated with SN (Rho et al. 2018) and did not find any features that indicated the presence of SN.

3.5. Summary of results

Figure 9. Flowchart summarizing the results characterizing the effectiveness of MIR variability in tracing AGN activity in low mass galaxies.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00666v1-figure-63e5a9ca-d121-4b1c-9422-b264c7323e93.jpg)

Dwarf galaxies with active black holes have been found to not display AGN signatures uniformly across different detection methods. If we have a dwarf galaxy that is variable in the MIR, how reliably can we determine that it is actually hosting an active black hole? In our work, we found that out of 25 dwarf galaxies that exhibit MIR variability, we have independent evidence that 17 of them (68%) host an AGN, strongly suggesting that MIR variability is due to the presence of an AGN in these targets. Twelve of these were classified as optical AGN, and the remaining five had NIR signatures of hosting an AGN. Four did not have any optical indicators of AGN, but we were not able to obtain NIR observations for them, so the true fraction of independently confirmed AGN in the sample of MIR variable AGN can be higher than 68% (see Fig. 9 for a summary).

Our results showed that if a galaxy is variable in the MIR, there is a high probability that it is actually hosting an active black hole. Using TESS light curves, Treiber et al. (2023) identified a similar fraction (62%) of optically variable galaxies that were either previously classified as AGNs in the literature or confirmed as AGNs based on a combination of emission-line diagnostics, mid-IR colors, or X-ray luminosity in their analysis. Thus, variability is a good indicator of AGN activity in dwarf galaxies and can be used to recover a larger fraction of AGN that is missed by other methods.

Signs of active black holes in dwarf galaxies are often difficult to distinguish from those originating from SF

activity. The low luminosity of the AGN can lead to the possibility that SF activity capable of producing similar luminosities often act as interlopers. In the infrared in particular, Hainline et al. (2016) found that most of the BPT-confirmed AGNs from Reines et al. (2013) were not selected as WISE AGN, which could be due to the host galaxy light dominating in the mid-infrared color space. Conversely, star-forming dwarf galaxies are capable of heating dust in such a way that mimics the infrared colors of more luminous AGNs. That is why W 1 -W 2 colors alone should not be used to select for AGN in dwarf galaxies.

## 4. COMPARISON WITH NON-VARIABLE DWARF GALAXIES

Given the potential contamination from SF processes at infrared wavelengths, it is crucial to test whether MIR variability is accurately tracing AGN activity or if it is being caused by SF interlopers in dwarf galaxies. We have found that MIR variability is frequently accompanied by alternate indicators of AGN (see Sec. 2.1 and Sec. 4), making it less likely that it is caused by stellar processes. Additionally, we compared dwarf galaxies that exhibit MIR variability with similar dwarf galaxies that do not exhibit MIR variability. We did this to recover differences in optical and infrared galaxy properties, if any, between the two samples of galaxies and understand the effects of any potential biases caused by host galaxy properties on the observed MIR variability.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00666v1-figure-5cc50eaa-9acc-4855-9527-154e748abd75.jpg)

∗

Figure 10. Histograms showing the distribution of the WISE W 1 -W 2 colors and stellar masses of the target sample of variable dwarf galaxies and the control sample of non-variable dwarf galaxies.

## 4.1. Control sample selection

We selected a control sample of non-variable dwarf galaxies, which were matched in WISE W 1 -W 2 colors and stellar mass with the sample of variable dwarf galaxies. We obtained five control galaxies on average for each of the 25 target galaxies that we initially classified as variable and thus created a control sample of 136 non-variable dwarf galaxies. A two-sample KolmogorovSmirnov (K-S) test between the W 1 -W 2 colors and stellar masses of the variable galaxies and non-variable control galaxies gives p=0.99 and p=0.56, respectively, indicating that the control galaxies do not have systematically different WISE colors or stellar masses (see Fig. 10). None of the galaxies in the control sample passed the variability tests based on our initial constraints (see Section 2.1), and thus can be safely assumed to be non-variable. The two samples have similar redshift distributions, ranging from 0.01 to 0.1.

vidual galaxies from the H α /H β Balmer decrement. We then obtained extinction corrected [O III ] λ 5007 ˚ luA minosity values for all the galaxies in the target and the control samples. We also measured the star-formation rates (SFRs) for the galaxies using a number of methods. We directly obtained the SFRs for the galaxies from the MPA-JHU catalog, which were calculated based on the methods used in Brinchmann et al. (2004). We independently estimated SFRs from the extinction-corrected H α luminosity values based on the relations by Kennicutt (1998) as well as from the extinction-corrected [O II ] fluxes using Equation 11 in Kewley et al. (2004). The SFRs calculated by these different optical methods agreed with each other. We determined the optical colors (u-r) from SDSS for the different galaxies.

We then obtained emission line flux values for [O II ] λλ 3727,3729 ˚, He A II λ 4686, H β λ 4861 A, [O ˚ III ] λ 5007 ˚ , [O A I ] λ 6300 ˚ , [N A II ] λ 6584 ˚ , H A α λ 6563 ˚ and A [S II ] λ 6731 ˚ from the MPA-JHU catalog (Kauffmann A et al. 2003a). We measured extinction values for indi-

## 4.2. Large-scale properties of variable and non-variable dwarf galaxies

Comparing various properties of variable and nonvariable targets with similar W 1 -W 2 colors, we did not find a significant difference between the extinction values, [O III ] luminosities, optical colors, and SFRs between both samples (Fig. 11).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00666v1-figure-179821f9-8c6e-4174-bcd5-119a42e1505e.jpg)

-

Figure 11. Histograms showing the comparison optical properties of the variable and the non-variable dwarf galaxies. The dashed lines indicate the average values for the two samples. We did not find a significant difference in several key optical properties, such as extinction, [O III ] luminosity, optical colors and SFRs.

Figure 12. The BPT diagram of the variable and non-variable samples. Galaxies with log [N II ]/H α &lt; -1 (indicative of low metallicity) are in the blue shaded region). Using log [N II ]/H α as a proxy for metallicity, we found that a larger fraction of variable dwarfs have higher metallicities than non-variable dwarfs.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00666v1-figure-ce51a5d4-8cc8-4b2f-8a43-75d81bab93a4.jpg)

Using the emission line values obtained from the MPAJHU catalog, we plotted the BPT diagrams of both samples (see Fig. 12). We found that 3/25 (12%) of the variable targets lie above the Kauffmann et al. (2003b) demarcation line for the AGN region, while 8/136 (0.05%) of the non-variable controls lie in the AGN region.

We obtained He II λ 4686 emission line fluxes for eight of the variable galaxies and eighty-two of the nonvariable galaxies from the MPA-JHU catalog. No significant difference was found between the positions of the variable and non-variable galaxies in the He II λ 4686/H β versus [N II ] λ 6584/H α diagnostic diagram presented in Shirazi &amp; Brinchmann (2012). While some of the variable galaxies were previously identified as AGNs using this diagnostic, we identified four additional variable galaxies that meet the AGN criteria. These galaxies are marked with blue checkmarks in Column 5 in Table 1.

Determining gas-phase metallicities for the targets is challenging as the traditional methods that depend on various emission lines such as [O III ] and [O II ] often do not account for the contamination due to AGN activity. Moreover, studies have shown that there is an ambiguity in determining metallicities using methods such as the ratio of [O II ] λλ 3726,3729 + [O III ] λλ 4959,5007 over

Figure 13. Plot of the H α FWHM vs the [OIII] FWHM for the variable and non-variable targets. A larger fraction of variable targets have very broad (FWHM &gt; 500 km s -1 ) Balmer lines than non-variable targets.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00666v1-figure-e8cc1e2e-3361-4f86-85e8-dedfd9d70ccf.jpg)

H β (R23, Pagel et al. 1979), particularly in the range of metallicities (7.5 &lt; 12 + log(O/H) &lt; 8.5) that we would expect these dwarf galaxies to have (Maiolino et al. 2008; Jones et al. 2015). Thus, to comment on the metallicities of the targets, we looked solely at the [N II ]/H α line ratios from the BPT diagram, which can be used as a proxy for metallicity (Kewley &amp; Dopita 2002; Pettini &amp; Pagel 2004; Kewley &amp; Ellison 2008).

We found that, on average, variable targets have higher metallicities (log [N II ]/H α &gt; -1) as compared to the non-variable targets. Only 2/25 (8%) targets have [N II ]/H α lower than -1 as compared to 59/136 (43%) for the control sample, as seen in Fig. 12. BPTselected AGNs in dwarf galaxies are expected to be biased towards higher metallicities (Reines et al. 2013). However, we adopted a completely independent method of selecting AGN in dwarf galaxies, and we still found that variable dwarf galaxies (tracing the AGN) have higher metallicities than non-variable dwarfs. Since AGN contribution can also lead to an increase in the [N II ]/H α ratio, a large fraction of variable targets with high [N II ]/H α ratios may be a direct consequence of the presence of AGNs in the variable targets.

## 4.3. Optical broad lines and W 2 -W 3 colors

We obtained the SDSS spectra of all the galaxies in the target and the control samples. We then performed spectral fits for the entire sample in the regions that include the prominent Balmer lines (4400 - 7000 ˚ ) and fit A simultaneously for H β λ 4861 ˚ , [O A III ] λλ 4963, 5007 ˚ , A H α λ 4861 ˚ , and [N A II ] λλ 6548, 6583 ˚ . We found that A variable targets have broader Balmer lines compared to non-variable targets (Fig. 13): 6/25 (24%) targets have H α FWHM &gt; 500 km s -1 while only 2/136 (0.015%) of the controls have H α FWHM &gt; 500 km s -1 . Broad Balmer lines could indicate the presence of AGN activity, but there could be potential contamination from star formation. We did not find any broad outflowing component in the forbidden [O III ] lines of both samples, indicating that neither sample has strong outflows in ionized gas.

This higher incidence of optical broad lines in the sample of galaxies exhibiting MIR variability as compared to non-variable galaxies is surprising. The presence of broad lines in the spectrum of AGN is typically attributed to an orientation effect. In type 1 AGNs, the nucleus is oriented so that we can observe the BLR and narrow line regions (NLR). MIR variability, on the other hand, should have no dependence on orientation, as the variability is believed to originate from the dust surrounding the accretion disk and is expected to be isotropic. However, Son et al. (2023) found that lowluminosity type 2 AGNs tend to have lower variability amplitudes than their type 1 counterparts. They attribute this to the type 2 AGN having distinctive central structures due to their low luminosity or their MIR brightness being contaminated by emission from the cold dust in the host galaxy. On the other hand, recent observations of extended (tens to hundreds of parsecs) MIR emission from the polar regions of AGNs suggest that MIR emission may be dominated by dust from the polar regions rather than from the torus (Asmus et al. 2016; Asmus 2019). This interpretation possibly suggests that we can simultaneously observe the BLR as well as MIR emission. Our observations of a large number of Type 1 AGN that are variable in MIR can be explained with either of these scenarios.

We also obtained WISE W 1, W 2, W 3, and W 4 fluxes for all the galaxies. Examining the infrared colors of both samples (Fig. 14), we found that non-variable targets have higher values of W 2 -W 3 (redder W 2 -W 3 colors) on average than variable targets (3.92 vs 3.18). Redder W 2 -W 3 colors indicate activity caused by starburst events (Jarrett et al. 2013). A p-value of 0.027 ( &lt; 0.05) indicates that there is a statistically significant difference in W 2 -W 3 between the variable and control samples. This suggests that variable targets have

Figure 14. The WISE color-color diagram for the variable and non-variable galaxies, with histograms plotting the distribution of the values in W 1 -W 2 (right) and W 2 -W 3 (top). Although both samples have similar W 1 -W 2 colors, non-variable galaxies have redder W 2 -W 3 colors.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00666v1-figure-cd23234e-2580-4e60-90ff-35b2ffe4e9d6.jpg)

bluer W 2 -W 3 colors, indicating that there is likely lesser contribution from dust heated by star formation in the MIR for the variable targets. This difference in the contribution of SF activity is not apparent from the optical, where both samples have similar SFRs (see Fig. 11). Even with similar W 1 -W 2 colors, variable dwarfs still have bluer W 2 -W 3 colors. Thus, it is likely that AGN activity causes redder W 1 -W 2 colors in variable targets, while for the non-variable galaxies, the contribution stems from strong star formation activity.

## 5. CONCLUSIONS

In this work, we selected a sample of twenty-five dwarf galaxies that are identified to be variable in the MIR from the light curves from the most recent release of the All WISE /NEO WISE catalog. We determined that 12/25 (48%) of these galaxies had been confirmed as AGNin the literature by several optical diagnostic methods, such as the BPT diagram, presence of broad lines or strong He II emission, indicating that MIR variability accurately traces black hole activity in these small dwarf galaxies. For a subsample of the variable dwarf galaxies that did not have any optical confirmations of AGN activity, we obtained NIR observations using Keck/NIRES. Our main results can be summarized as follows (also see Fig. 9):

- 1. We present the first detection of a NIR coronal line ([S IX]) in J1205, which has the lowest mass (log(M ) = 7.5 M ∗ ⊙ ) among all the objects in our

- sample. Most commonly used AGN diagnostics, such as BPT diagrams and He II emission lines, fail to detect the presence of an AGN in J1205, thus highlighting that MIR variability can be a powerful and sometimes the only tool to detect AGNs in low-mass galaxies.
- 2. We find evidence of a broad Pa α line potentially from the BLR in two of the variable targets (J0928 and J1144). These targets lie in the SF region of the optical BPT diagram, but they lie in the composite or AGN regions of the NIR diagnostic diagrams. The black hole masses implied by their broad lines are consistent with the black hole to stellar mass relation from RV15.
- 3. Using NIR diagnostic diagrams, which make use of NIR emission line ratios in order to differentiate between AGN and SF regions in galaxies, we find that several of the variable targets (4/9) that we observed moved into the AGN/composite region of the diagnostic diagram compared to their location in the SF regions of the optical BPT diagram.
- 4. Comparing a sample of variable and non-variable dwarf galaxies with similar stellar mass and WISE W 1 -W 2 colors, we find that variable dwarfs have bluer W 2 -W 3 colors than nonvariable dwarfs. The redder W 1 -W 2 colors of non-variable dwarfs are likely because of high star formation. Thus, variable dwarfs accurately trace AGN activity independent of mid-infrared colors.
- 5. Variable dwarfs have a higher fraction of galaxies that have broad ( &gt; 500 km s -1 ) Balmer (H β and H α ) lines than the non-variable dwarfs (24% compared to 0.015%). While this could simply indicate a higher incidence of AGN, it is remarkable that such a large fraction would be type 1 AGN.

We find independent confirmation of AGN activity in at least 68 % out of 25 MIR variable dwarf galaxies, indicating that MIR variability can be used to accurately detect AGN in dwarf galaxies. Improving methods to detect low-mass active black holes is the first step in constraining the black hole occupation fraction in dwarf galaxies, thus determining the evolutionary pathways for black hole seeds in the early universe. Our study shows that using indicators in the infrared is crucial to obtaining a more complete census of black holes in dwarf galaxies. With the increasing amount of infrared data being available from JWST and the upcoming infrared transient survey telescopes such as the Palomar GattiniIR (Moore &amp; Kasliwal 2019), harnessing the power of MIR variability will likely be fruitful in uncovering large numbers of active black holes in dwarf galaxies.

We thank the anonymous referee for their thoughtful feedback and constructive comments that greatly helped to improve this paper. Partial support for this project was provided by the National Science Foundation under grant No. AST 1817233. The data presented herein were obtained at the W. M. Keck Observatory, which is operated as a scientific partnership among the California Institute of Technology, the University of California, and the National Aeronautics and Space Administration. The Observatory was made possible by the generous financial support of the W. M. Keck Foundation. The authors wish to recognize and acknowledge the very significant cultural role and reverence that the summit of Maunakea has always had within the indigenous Hawaiian community. We are most fortunate to have the opportunity to conduct observations from this mountain.

This publication makes use of data products from the Wide-field Infrared Survey Explorer, which is a joint project of the University of California, Los Angeles, and the Jet Propulsion Laboratory/California Institute of Technology, and NEO WISE , which is a project of the Jet Propulsion Laboratory/ California Institute of Technology. WISE and NEO WISE are funded by the National Aeronautics and Space Administration.

SDSS-IV is managed by the Astrophysical Research Consortium for the Participating Institutions of the SDSS Collaboration. Funding for the Sloan Digital Sky Survey IV has been provided by the Alfred P. Sloan Foundation, the U.S. Department of Energy Office of Science, and the Participating Institutions. SDSS-IV acknowledges support and resources from the Center for High-Performance Computing at the University of Utah.

## Facilities: WISE , NEO WISE , KCWI/Keck, SDSS

Software: astropy (Astropy Collaboration et al. 2013, 2018), BADASS (Sexton et al. 2021), pPXF, (Cappellari &amp; Emsellem 2004), REDSPEC (https:// www2.keck.hawaii.edu/inst/nirspec/redspec.html), SciencePlots (Garrett 2021)

Table 4. 3 σ Upper limits to the [Si VI ] flux measurements for the non-detections of coronal lines in this galaxies observed with NIRES (See Figure 5)

| Target   |   [Si VI ] flux upper limits (10 - 17 erg cm - 2 s - 1 ) |
|----------|----------------------------------------------------------|
| J0847    |                                                     8.12 |
| J0928    |                                                   100.97 |
| J1032    |                                                    81.28 |
| J1144    |                                                    19.47 |
| J1205    |                                                     5.55 |
| J1320    |                                                     3.31 |
| J1358    |                                                    27.54 |
| J1428    |                                                     0.55 |
| J1612    |                                                     4.57 |

The non-uniform detection of coronal lines across galaxies is likely due to insufficient observation times in order to detect the typically faint coronal lines. We present the upper limits for the [Si VI ] flux for the nine galaxies observed with NIRES in this paper (Table 4). The upper limits to the [Si VI ] fluxes were determined by estimating the 1 σ noise levels in the region around [Si VI ] flux and multiplying it by three to obtain the 3 σ upper limit indicated in Table 4. None of the observed galaxies had [Si VI ] detections.

## REFERENCES

Annuar, A., Gandhi, P., Alexander, D. M., et al. 2015, ApJ, 815, 36, doi: 10.1088/0004-637X/815/1/36 Annuar, A., Alexander, D. M., Gandhi, P., et al. 2017, ApJ, 836, 165, doi: 10.3847/1538-4357/836/2/165 Arcavi, I., Hosseinzadeh, G., Howell, D. A., et al. 2017, Nature, 551, 64, doi: 10.1038/nature24291 Asmus, D. 2019, MNRAS, 489, 2177, doi: 10.1093/mnras/stz2289 Asmus, D., H¨nig, S. F., &amp; Gandhi, P. 2016, ApJ, 822, 109, o doi: 10.3847/0004-637X/822/2/109 Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&amp;A, 558, A33, doi: 10.1051/0004-6361/201322068 Astropy Collaboration, Price-Whelan, A. M., Sip˝cz, B. M., o et al. 2018, AJ, 156, 123, doi: 10.3847/1538-3881/aabc4f Baldassare, V. F., Geha, M., &amp; Greene, J. 2020, ApJ, 896, 10, doi: 10.3847/1538-4357/ab8936 Baldassare, V. F., Reines, A. E., Gallo, E., &amp; Greene, J. E. 2017, ApJ, 836, 20, doi: 10.3847/1538-4357/836/1/20 Baldassare, V. F., Reines, A. E., Gallo, E., et al. 2016, ApJ, 829, 57, doi: 10.3847/0004-637X/829/1/57 Baldwin, J. A., Phillips, M. M., &amp; Terlevich, R. 1981, PASP, 93, 5, doi: 10.1086/130766 Becker, G. D., Rauch, M., &amp; Sargent, W. L. W. 2009, ApJ, 698, 1010, doi: 10.1088/0004-637X/698/2/1010

Begelman, M. C., Volonteri, M., &amp; Rees, M. J. 2006, MNRAS, 370, 289, doi: 10.1111/j.1365-2966.2006.10467.x Bellm, E. C., Kulkarni, S. R., Graham, M. J., et al. 2019, PASP, 131, 018002, doi: 10.1088/1538-3873/aaecbe Bi, S., Feng, H., &amp; Ho, L. C. 2020, ApJ, 900, 124, doi: 10.3847/1538-4357/aba761 Blanton, M. R., &amp; Roweis, S. 2007, AJ, 133, 734, doi: 10.1086/510127 Bogd´ an, A., Goulding, A. D., Natarajan, P., et al. 2024, ´ Nature Astronomy, 8, 126, doi: 10.1038/s41550-023-02111-9 Bohn, T., Canalizo, G., Satyapal, S., &amp; Pfeifle, R. W. 2020, ApJ, 899, 82, doi: 10.3847/1538-4357/aba52c Bohn, T., Canalizo, G., Veilleux, S., &amp; Liu, W. 2021, ApJ, 911, 70, doi: 10.3847/1538-4357/abe70c Bose, S., Deason, A. J., &amp; Frenk, C. S. 2018, ApJ, 863, 123, doi: 10.3847/1538-4357/aacbc4 Brinchmann, J., Charlot, S., White, S. D. M., et al. 2004, MNRAS, 351, 1151, doi: 10.1111/j.1365-2966.2004.07881.x Bromm, V., &amp; Loeb, A. 2003, ApJ, 596, 34, doi: 10.1086/377529 Burke, C. J., Liu, Y., Ward, C. A., et al. 2024, arXiv e-prints, arXiv:2402.06882, doi: 10.48550/arXiv.2402.06882

Beckmann, R. S., Dubois, Y., Volonteri, M., et al. 2023, MNRAS, 523, 5610, doi: 10.1093/mnras/stad1544

Burke, C. J., Shen, Y., Chen, Y.-C., et al. 2020, ApJ, 899, 136, doi: 10.3847/1538-4357/aba3ce

| Burke, C. J., Shen, Y., Liu, X., et al. 2023, MNRAS, 518, 1880, doi: 10.1093/mnras/stac2478                                               |
|-------------------------------------------------------------------------------------------------------------------------------------------|
| Burke, C. J., Shen, Y., Blaes, O., et al. 2021, Science, 373, 789, doi: 10.1126/science.abg9933                                           |
| Burke, C. J., Liu, X., Shen, Y., et al. 2022, MNRAS, 516, 2736, doi: 10.1093/mnras/stac2262                                               |
| Calabr`, o A., Pentericci, L., Feltre, A., et al. 2023, A&A, 679, A80, doi: 10.1051/0004-6361/202347190                                   |
| Cann, J. M., Satyapal, S., Abel, N. P., et al. 2019, ApJL, 870, L2, doi: 10.3847/2041-8213/aaf88d                                         |
| Cann, J. M., Satyapal, S., Rothberg, B., et al. 2021, ApJL, 912, L2, doi: 10.3847/2041-8213/abf56d                                        |
| Cappellari, M., & Emsellem, E. 2004, PASP, 116, 138, doi: 10.1086/381875                                                                  |
| Chadayammuri, U., Bogd´n, a ´ A., Ricarte, A., & Natarajan, P. 2023, ApJ, 946, 51, doi: 10.3847/1538-4357/acbea6                          |
| Chilingarian, I. V., Katkov, I. Y., Zolotukhin, I. Y., et al. 2018, ApJ, 863, 1, doi: 10.3847/1538-4357/aad184                            |
| Clark, P., Graur, O., Callow, J., et al. 2024, MNRAS, 528, 7076, doi: 10.1093/mnras/stae460                                               |
| Colina, L., Piqueras L´pez, o J., Arribas, S., et al. 2015, A&A, 578, A48, doi: 10.1051/0004-6361/201425567                               |
| Dark Energy Survey Collaboration, Abbott, T., Abdalla, F. B., et al. 2016, MNRAS, 460, 1270, doi: 10.1093/mnras/stw641                    |
| De, K., MacLeod, M., Karambelkar, V., et al. 2023, Nature, 617, 55, doi: 10.1038/s41586-023-05842-x                                       |
| den Brok, J. S., Koss, M. J., Trakhtenbrot, B., et al. 2022, ApJS, 261, 7, doi: 10.3847/1538-4365/ac5b66                                  |
| den Brok, M., Seth, A. C., Barth, A. J., et al. 2015, ApJ, 809, 101, doi: 10.1088/0004-637X/809/1/101                                     |
| Devecchi, B., Volonteri, M., Rossi, E. M., Colpi, M., & Portegies Zwart, S. 2012, MNRAS, 421, 1465, doi: 10.1111/j.1365-2966.2012.20406.x |
| Do, T., Wright, S. A., Barth, A. J., et al. 2014, AJ, 147, 93, doi: 10.1088/0004-6256/147/4/93                                            |
| Donley, J. L., Rieke, G. H., P´rez-Gonz´lez, e a P. G., & Barro, G. 2008, ApJ, 687, 111, doi: 10.1086/591510                              |
| Finkelstein, S. L., Bagley, M. B., Ferguson, H. C., et al. 2023, ApJL, 946, L13, doi: 10.3847/2041-8213/acade4                            |
| Garrett, J. 2021, doi: 10.5281/zenodo.5512926                                                                                             |
| Goulding, A. D., Greene, J. E., Setton, D. J., et al. 2023, ApJL, 955, L24, doi: 10.3847/2041-8213/acf7c5                                 |
| Greene, J. E. 2012, Nature Communications, 3, 1304, doi: 10.1038/ncomms2314                                                               |
| Greene, J. E., Strader, J., & Ho, L. C. 2020, ARA&A, 58, 257, doi: 10.1146/annurev-astro-032620-021835                                    |

| Haidar, H., Habouzit, M., Volonteri, M., et al. 2022, MNRAS, 514, 4912, doi: 10.1093/mnras/stac1659                                                   |
|-------------------------------------------------------------------------------------------------------------------------------------------------------|
| Hainline, K. N., Reines, A. E., Greene, J. E., & Stern, D. 2016, ApJ, 832, 119, doi: 10.3847/0004-637X/832/2/119                                      |
| Harikane, Y., Ouchi, M., Oguri, M., et al. 2023, ApJS, 265, 5, doi: 10.3847/1538-4365/acaaa9                                                          |
| Harish, S., Malhotra, S., Rhoads, J. E., et al. 2023, ApJ, 945, 157, doi: 10.3847/1538-4357/acb99c                                                    |
| Hinkle, J. T., Tucker, M. A., Shappee, B. J., et al. 2023, MNRAS, 519, 2035, doi: 10.1093/mnras/stac3659                                              |
| Inayoshi, K., Visbal, E., & Haiman, Z. 2020, ARA&A, 58, 27, doi: 10.1146/annurev-astro-120419-014455                                                  |
| Ivezi´, c ˇ Z., Kahn, S. M., Tyson, J. A., et al. 2019, ApJ, 873, 111, doi: 10.3847/1538-4357/ab042c                                                  |
| Izotov, Y. I., & Thuan, T. X. 2009, ApJ, 707, 1560, doi: 10.1088/0004-637X/707/2/1560                                                                 |
| Izotov, Y. I., Thuan, T. X., & Guseva, N. G. 2017, MNRAS, 471, 548, doi: 10.1093/mnras/stx1629 -. 2021, MNRAS, 508, 2556, doi: 10.1093/mnras/stab2798 |
| A&A, 661, A80, doi: 10.1051/0004-6361/202142663 Jarrett, T. H., Masci, F., Tsai, C. W., et al. 2013, AJ, 145,                                         |
| 6, doi: 10.1088/0004-6256/145/1/6                                                                                                                     |
| Jones, T., Martin, C., & Cooper, M. C. 2015, ApJ, 813, 126, doi: 10.1088/0004-637X/813/2/126                                                          |
| Karambelkar, V. R., Adams, S. M., Whitelock, P. A., et al. 2019, ApJ, 877, 110, doi: 10.3847/1538-4357/ab1a41                                         |
| Kauffmann, G., Heckman, T. M., White, S. D. M., et al. 2003a, MNRAS, 341, 33, doi: 10.1046/j.1365-8711.2003.06291.x                                   |
| Kauffmann, G., Heckman, T. M., Tremonti, C., et al. 2003b, MNRAS, 346, 1055, doi: 10.1111/j.1365-2966.2003.07154.x                                    |
| Kawaguchi, T., & Mori, M. 2011, ApJ, 737, 105, doi: 10.1088/0004-637X/737/2/105                                                                       |
| Kelson, D. D. 2003, PASP, 115, 688, doi: 10.1086/375502                                                                                               |
| Kennicutt, Robert C., J. 1998, ARA&A, 36, 189, doi: 10.1146/annurev.astro.36.1.189                                                                    |
| Kewley, L. J., & Dopita, M. A. 2002, ApJS, 142, 35, doi: 10.1086/341326                                                                               |
| Kewley, L. J., & Ellison, S. L. 2008, ApJ, 681, 1183, doi: 10.1086/587500                                                                             |
| Kewley, L. J., Geller, M. J., & Jansen, R. A. 2004, AJ, 127, 2002, doi: 10.1086/382723                                                                |
| Kim, D., Im, M., Canalizo, G., et al. 2018, ApJS, 238, 37, doi: 10.3847/1538-4365/aadfd5                                                              |
| Kimura, Y., Yamada, T., Kokubo, M., et al. 2020, ApJ, 894, 24, doi: 10.3847/1538-4357/ab83f3 Kocevski, D. D., Onoue, M., Inayoshi, K., et al. 2023,   |
| ApJL, 954, L4, doi: 10.3847/2041-8213/ace5a0                                                                                                          |

| Kokorev, V., Fujimoto, S., Labbe, I., et al. 2023, ApJL, 957, L7, doi: 10.3847/2041-8213/ad037a                                                   |
|---------------------------------------------------------------------------------------------------------------------------------------------------|
| Lamperti, I., Koss, M., Trakhtenbrot, B., et al. 2017, MNRAS, 467, 540, doi: 10.1093/mnras/stx055                                                 |
| Landt, H., Bentz, M. C., Ward, M. J., et al. 2008, ApJS, 174, 282, doi: 10.1086/522373                                                            |
| Larkin, J. E., Armus, L., Knop, R. A., Soifer, B. T., & Matthews, K. 1998, ApJS, 114, 59, doi: 10.1086/313063                                     |
| Law, N. M., Kulkarni, S. R., Dekany, R. G., et al. 2009, PASP, 121, 1395, doi: 10.1086/648598                                                     |
| Lequeux, J., Peimbert, M., Rayo, J. F., Serrano, A., & Torres-Peimbert, S. 1979, A&A, 80, 155                                                     |
| Li, J., & Shen, Y. 2023, ApJ, 950, 122, doi: 10.3847/1538-4357/accade                                                                             |
| Lodato, G., & Natarajan, P. 2006, MNRAS, 371, 1813, doi: 10.1111/j.1365-2966.2006.10801.x                                                         |
| Lupi, A., Haardt, F., Dotti, M., et al. 2016, MNRAS, 456, 2993, doi: 10.1093/mnras/stv2877                                                        |
| Lyke, B. W., Higley, A. N., McLane, J. N., et al. 2020, ApJS, 250, 8, doi: 10.3847/1538-4365/aba623                                               |
| MacLeod, C. L., Ivezi´, c ˇ Z., Kochanek, C. S., et al. 2010, ApJ, 721, 1014, doi: 10.1088/0004-637X/721/2/1014                                   |
| Madau, P., & Rees, M. J. 2001, ApJL, 551, L27, doi: 10.1086/319848                                                                                |
| Maiolino, R., Nagao, T., Grazian, A., et al. 2008, A&A, 488, 463, doi: 10.1051/0004-6361:200809678                                                |
| Maiolino, R., Scholtz, J., Curtis-Lake, E., et al. 2023, arXiv e-prints, arXiv:2308.01230, doi: 10.48550/arXiv.2308.01230                         |
| Maraston, C., Daddi, E., Renzini, A., et al. 2006, ApJ, 652, 85, doi: 10.1086/508143                                                              |
| Marleau, F. R., Clancy, D., Habas, R., & Bianconi, M. 2017, A&A, 602, A28, doi: 10.1051/0004-6361/201629832                                       |
| Mart´nez-Palomera, ı J., Lira, P., Bhalla-Ladd, I., F¨rster, o F., & Plotkin, R. M. 2020, ApJ, 889, 113, doi: 10.3847/1538-4357/ab5f5b            |
| Mart´nez-Palomera, ı J., F¨rster, o F., Protopapas, P., et al. 2018, AJ, 156, 186, doi: 10.3847/1538-3881/aadfd8                                  |
| Martins, L. P., Rodr´guez-Ardila, ı A., Diniz, S., Gruenwald, R., & de Souza, R. 2013, MNRAS, 431, 1823, doi: 10.1093/mnras/stt296                |
| Mason, R. E., Rodr´guez-Ardila, ı A., Martins, L., et al. 2015, ApJS, 217, 13, doi: 10.1088/0067-0049/217/1/13                                    |
| Mezcua, M. 2017, International Journal of Modern Physics D, 26, 1730021, doi: 10.1142/S021827181730021X                                           |
| Mezcua, M., Civano, F., Marchesi, S., et al. 2018, MNRAS, 478, 2576, doi: 10.1093/mnras/sty1163                                                   |
| Mezcua, M., & Dom´nguez ı S´nchez, a H. 2020, ApJL, 898, L30, doi: 10.3847/2041-8213/aba199 -. 2024, MNRAS, 528, 5252, doi: 10.1093/mnras/stae292 |

| Miller, B. P., Gallo, E., Greene, J. E., et al. 2015, ApJ, 799, 98, doi: 10.1088/0004-637X/799/1/98                                                         |
|-------------------------------------------------------------------------------------------------------------------------------------------------------------|
| Moore, A. M., & Kasliwal, M. M. 2019, Nature Astronomy, 3, 109, doi: 10.1038/s41550-018-0675-x                                                              |
| Moran, E. C., Shahinyan, K., Sugarman, H. R., V´lez, e D. O., & Eracleous, M. 2014, AJ, 148, 136, doi: 10.1088/0004-6256/148/6/136                          |
| M¨ller-S´nchez, u a F., Hicks, E. K. S., Malkan, M., et al. 2018, ApJ, 858, 48, doi: 10.3847/1538-4357/aab9ad                                               |
| Natarajan, P. 2014, General Relativity and Gravitation, 46, 1702, doi: 10.1007/s10714-014-1702-6 -. 2021, MNRAS, 501, 1413, doi: 10.1093/mnras/staa3724     |
| Nguyen, D. D., Seth, A. C., Neumayer, N., et al. 2018, ApJ, 858, 118, doi: 10.3847/1538-4357/aabe28 -. 2019, ApJ, 872, 104, doi: 10.3847/1538-4357/aafe7a   |
| Oh, K., Yi, S. K., Schawinski, K., et al. 2015, ApJS, 219, 1, doi: 10.1088/0067-0049/219/1/1                                                                |
| Ohlson, D., Seth, A. C., Gallo, E., Baldassare, V. F., & Greene, J. E. 2024, AJ, 167, 31, doi: 10.3847/1538-3881/acf7bc                                     |
| Pace, Z. J., Tremonti, C., Chen, Y., et al. 2019, ApJ, 883, 83, doi: 10.3847/1538-4357/ab3726                                                               |
| Pagel, B. E. J., Edmunds, M. G., Blackwell, D. E., Chun, M. S., & Smith, G. 1979, MNRAS, 189, 95, doi: 10.1093/mnras/189.1.95                               |
| Paudel, S., Hilker, M., Ree, C. H., & Kim, M. 2016, ApJL, 820, L19, doi: 10.3847/2041-8205/820/1/L19                                                        |
| Petter, G. C., Hickox, R. C., Alexander, D. M., et al. 2023, ApJ, 946, 27, doi: 10.3847/1538-4357/acb7ef                                                    |
| Pettini, M., & Pagel, B. E. J. 2004, MNRAS, 348, L59, doi: 10.1111/j.1365-2966.2004.07591.x                                                                 |
| Polimera, M. S., Kannappan, S. J., Richardson, C. T., et al. 2022, ApJ, 931, 44, doi: 10.3847/1538-4357/ac6595                                              |
| Portegies Zwart, S. F., Baumgardt, H., Hut, P., Makino, J., & McMillan, S. L. W. 2004, Nature, 428, 724, doi: 10.1038/nature02448                           |
| Reines, A. E. 2022, Nature Astronomy, 6, 26, doi: 10.1038/s41550-021-01556-0                                                                                |
| Reines, A. E., Greene, J. E., & Geha, M. 2013, ApJ, 775, 116, doi: 10.1088/0004-637X/775/2/116                                                              |
| Reines, A. E., & Volonteri, M. 2015, ApJ, 813, 82, doi: 10.1088/0004-637X/813/2/82                                                                          |
| Rho, J., Geballe, T. R., Banerjee, D. P. K., et al. 2018, ApJL, 864, L20, doi: 10.3847/2041-8213/aad77f                                                     |
| Ricci, C., Bauer, F. E., Arevalo, P., et al. 2016, ApJ, 820, 5, doi: 10.3847/0004-637X/820/1/5 Riffel, R., Rodr´guez-Ardila, ı A., & Pastoriza, M. G. 2006, |

| Rodr´guez-Ardila, ı A., Prieto, M. A., Portilla, J. G., & Tejeiro, J. M. 2011, ApJ, 743, 100, doi: 10.1088/0004-637X/743/2/100                        |
|-------------------------------------------------------------------------------------------------------------------------------------------------------|
| Rosswog, S. 2005, ApJ, 634, 1202, doi: 10.1086/497062                                                                                                 |
| Ruan, J. J., Anderson, S. F., Plotkin, R. M., et al. 2014, ApJ, 797, 19, doi: 10.1088/0004-637X/797/1/19                                              |
| Salehirad, S., Reines, A. E., & Molina, M. 2022, ApJ, 937, 7, doi: 10.3847/1538-4357/ac8876                                                           |
| Salim, S., Rich, R. M., Charlot, S., et al. 2007, ApJS, 173,                                                                                          |
| 267, doi: 10.1086/519218 Sartori, L. F., Schawinski, K., Treister, E., et al. 2015,                                                                   |
| MNRAS, 454, 3722, doi: 10.1093/mnras/stv2238                                                                                                          |
| Sesar, B., Ivezi´, c ˇ Z., Lupton, R. H., et al. 2007, AJ,                                                                                            |
| 134, 2236, doi: 10.1086/521819                                                                                                                        |
| Sexton, R. O., Matzko, W., Darden, N., Canalizo, G., & Gorjian, V. 2021, MNRAS, 500, 2871, doi: 10.1093/mnras/staa3278                                |
| Shin, L., Woo, J.-H., Son, D., et al. 2022, AJ, 163, 73, doi: 10.3847/1538-3881/ac4038                                                                |
| Shirazi, M., & Brinchmann, J. 2012, MNRAS, 421, 1043, doi: 10.1111/j.1365-2966.2012.20439.x                                                           |
| Silk, J., & Rees, M. J. 1998, A&A, 331, L1, doi: 10.48550/arXiv.astro-ph/9801013                                                                      |
| Smith, K. L., Koss, M., & Mushotzky, R. F. 2014, ApJ, 794, 112, doi: 10.1088/0004-637X/794/2/112 Smitka, M. T. 2016, PhD thesis, Texas A&M University |
| System Son, S., Kim, M., & Ho, L. C. 2022, ApJ, 927, 107, doi: 10.3847/1538-4357/ac4dfc                                                               |
| -. 2023, ApJ, 958, 135, doi: 10.3847/1538-4357/ad01bc                                                                                                 |
| Sun, L., Xiao, L., & Li, G. 2022, MNRAS, 513, 4057,                                                                                                   |
| doi: 10.1093/mnras/stac1121                                                                                                                           |
| Treiber, H. P., Hinkle, J. T., Fausnaugh, M. M., et al. MNRAS, 525, 5795, doi: 10.1093/mnras/stad2530                                                 |

- Ulrich, M.-H., Maraschi, L., &amp; Urry, C. M. 1997, ARA&amp;A, 35, 445, doi: 10.1146/annurev.astro.35.1.445
- van Wassenhove, S., Volonteri, M., Walker, M. G., &amp; Gair, J. R. 2010, MNRAS, 408, 1139,
- doi: 10.1111/j.1365-2966.2010.17189.x
- Vanden Berk, D. E., Wilhite, B. C., Kron, R. G., et al. 2004, ApJ, 601, 692, doi: 10.1086/380563
- Vazdekis, A., Koleva, M., Ricciardelli, E., R¨ck, B., &amp; o Falc´n-Barroso, J. 2016, MNRAS, 463, 3409, o doi: 10.1093/mnras/stw2231
- Veilleux, S., Goodrich, R. W., &amp; Hill, G. J. 1997, ApJ, 477, 631, doi: 10.1086/303735
- Veilleux, S., &amp; Osterbrock, D. E. 1987, ApJS, 63, 295, doi: 10.1086/191166
- Veron-Cetty, M. P., &amp; Veron, P. 2006, VizieR Online Data Catalog: Quasars and Active Galactic Nuclei (12th Ed.), CDS/ADC Collection of Electronic Catalogues, Vol. 7248, p. 0-1 (2006)
- Volonteri, M., &amp; Begelman, M. C. 2010, MNRAS, 409, 1022, doi: 10.1111/j.1365-2966.2010.17359.x
- Volonteri, M., Habouzit, M., &amp; Colpi, M. 2021, Nature Reviews Physics, 3, 732, doi: 10.1038/s42254-021-00364-9

Volonteri, M., &amp; Rees, M. J. 2005, ApJ, 633, 624, doi: 10.1086/466521

- Ward, C., Gezari, S., Nugent, P., et al. 2022, ApJ, 936, 104, doi: 10.3847/1538-4357/ac8666
- Wasleske, E. J., &amp; Baldassare, V. F. 2024, arXiv e-prints, arXiv:2405.20312, doi: 10.48550/arXiv.2405.20312

Wilson, J. C., Henderson, C. P., Herter, T. L., et al. 2004, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 5492, Ground-based Instrumentation for Astronomy, ed. A. F. M. Moorwood &amp; M. Iye, 1295-1305, doi: 10.1117/12.550925

Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868, doi: 10.1088/0004-6256/140/6/1868

- York, D. G., Adelman, J., Anderson, John E., J., et al. 2000, AJ, 120, 1579, doi: 10.1086/301513
Archana Aravindan Gabriela Canalizo Nathan Secrest Shobita Satyapal Thomas Bohn
A closer look at dwarf galaxies exhibiting MIR variability: AGN confirmation and comparison with non-variable dwarf galaxies
Detecting active black holes in dwarf galaxies has proven to be a challenge
due to their small size and weak electromagnetic signatures. Mid-infrared
variability has emerged as a promising tool that can be used to detect active
low-mass black holes in dwarf galaxies. We analyzed 10.4 years of photometry
from the ALL$WISE$/NEO$WISE$ multi-epoch catalogs, identifying 25 objects with
AGN-like variability. Independent confirmation of AGN activity was found in 68%
of these objects using optical and near-infrared diagnostics. Notably, we
discovered a near-infrared coronal line [S IX] $\lambda$ 1.252 $\mu$m in J1205,
the galaxy with the lowest stellar mass (log M$_{*}$ = 7.5 M$_{\odot}$) and low
metallicity (12 + log(O/H) = 7.46) in our sample. Additionally, we found broad
Pa$\alpha$ potentially from the BLR in two targets, and their implied black
hole masses are consistent with black hole-stellar mass relations. Comparing
non-variable galaxies with similar stellar masses and $WISE$ $W1-W2$ colors, we
found no clear trends between variability and large-scale galaxy properties.
However, we found that AGN activity likely causes redder $W1-W2$ colors in
variable targets, while for the non-variable galaxies, the contribution stems
from strong star formation activity. A high incidence of optical broad lines
was also observed in variable targets. Our results suggest that mid-infrared
variability is an effective method for detecting AGN activity in low-mass
galaxies and can help uncover a larger sample of active low-mass ($<$ 10$^{6}$
M$_{\odot}$) black holes in the universe.