W1 Bands Research Articles

A Data-driven Search For Mid-infrared Excesses Among Five Million Main-sequence FGK Stars

Stellar infrared excesses can indicate various phenomena of interest, from protoplanetary disks to debris disks, or (more speculatively) techno-signatures along the lines of Dyson spheres. In this paper, we conduct a large search for “extreme” infrared excesses, designed as a data-driven contextual anomaly detection pipeline. We focus our search on FGK stars close to the main sequence to favor nonyoung host stars. We look for excess in the mid-infrared, unlocking a large sample to search in while favoring extreme IR excess akin to the ones produced by extreme debris disks (EDDs) and/or planetary collision events. We combine observations from ESA Gaia Data Release 3, the Two Micron All-Sky Survey, and the unWISE version of NASA’s Wide-field Infrared Survey Explorer (WISE), and create a catalog of 4,898,812 stars with G < 16 mag. We consider a star to have an excess if it is substantially brighter in the W1 and W2 bands than what is predicted from an ensemble of machine learning models trained on the data, taking optical and near-infrared information as input features. We apply a set of additional cuts (derived from the machine learning models and the objects’ astronomical features) to avoid false positives and identify a set of 53 objects, including one previously identified EDD candidate. The typical infrared-excess fractional luminosities we find are in the range 0.005–0.1, consistent with previous EDD candidates and potential planetary collision events.

Transit timing variation of K2-237b: hints toward planet disc migration

ABSTRACT Hot Jupiters should initially form at considerable distances from host stars and subsequently migrate towards inner regions, supported directly by transit timing variation (TTV). We report the TTV of K2-237b, using reproduced timings fitted from Kepler K2 and TESS data. The timings span from 2016 to 2021, leading to an observational baseline of 5 yr. The timing evolution presents a significant bias to a constant period scenario. The model evidence is evaluated utilizing the Bayesian Information Criterion (BIC), which favours the scenario of period decay with a $\Delta$BIC of 14.1. The detected TTV induces a period decay rate ($\dot{P}$) of −1.14$\pm 0.28\times 10^{-8}$ d per day ($-$0.36 s/yr). Fitting the spectral energy distribution, we find infrared excess at the significance level of 1.5 $\sigma$ for WISE W1 and W2 bands, and 2 $\sigma$ level for W3 and W4 bands. This potentially reveals the existence of a stellar disc, consisting of hot dust at 800 $\pm$ 300 K, showing a $L_{\mathrm{ dust}}/L_{\ast }$ of 5 $\pm\ 3\times 10^{-3}$. We obtain a stellar age of 1.0$^{+1.4}_{-0.7}\, \times 10^{9}$ yr from isochrone fitting. The properties of K2-237b potentially serve as a direct observational support to the planet disc migration though more observation are needed.

A Comprehensive Study on the Mid-Infrared Variability of Blazars

We present a comprehensive investigation of mid-infrared (MIR) flux variability at 3.4 μm (W1 band) for a large sample of 3816 blazars, using Wide-field Infrared Survey Explorer (WISE) data through December 2022. The sample consists of 1740 flat-spectrum radio quasars (FSRQs), 1281 BL Lac objects (BL Lacs), and 795 blazars of uncertain type (BCUs). Considering Fermi Large Area Telescope detection, we classify 2331 as Fermi blazars and 1485 as non-Fermi blazars. Additionally, based on synchrotron peak frequency, the sample includes 2264 low-synchrotron peaked (LSP), 512 intermediate-synchrotron peaked (ISP), and 655 high-synchrotron peaked (HSP) sources. We conduct a comparative analysis of short- and long-term intrinsic variability amplitude (σm), duty cycle (DC), and ensemble structure function (ESF) across blazar subclasses. The median short-term σm values were 0.181−0.106+0.153, 0.104−0.054+0.101, 0.135−0.076+0.154, 0.173−0.097+0.158, 0.177−0.100+0.156, 0.096−0.050+0.109, and 0.106−0.058+0.100 mag for FSRQs, BL Lacs, Fermi blazars, non-Fermi blazars, LSPs, ISPs, and HSPs, respectively. The median DC values were 71.03−22.48+14.17, 64.02−22.86+16.97, 68.96−25.52+15.66, 69.40−22.17+14.42, 71.24−21.36+14.25, 63.03−33.19+16.93, and 64.63−24.26+15.88 percent for the same subclasses. The median long-term σm values were 0.137−0.105+0.408, 0.171−0.132+0.206, 0.282−0.184+0.332, 0.071−0.062+0.143, 0.218−0.174+0.386, 0.173−0.132+0.208, and 0.101−0.077+0.161 mag for the same subclasses, respectively. Our results reveal significant differences in 3.4 μm flux variability among these subclasses. FSRQs (LSPs) exhibit larger σm and DC values compared to BL Lacs (ISPs and HSPs). Fermi blazars display higher long-term σm but lower short-term σm relative to non-Fermi blazars, while DC distributions between the two groups are similar. ESF analysis further confirms the greater variability of FSRQs, LSPs, and Fermi blazars across a wide range of time scales compared to BL Lacs, ISPs/HSPs, and non-Fermi blazars. These findings highlight a close correlation between MIR variability and blazar properties, providing valuable insights into the underlying physical mechanisms responsible for their emission.

A population of mid-infrared large-amplitude variable young stellar objects from unTimely

ABSTRACT Utilizing a decade-long unTimely data set, supplemented by multiband data from archives, we search for young stellar objects (YSOs) with variations larger than one magnitude in W1 band within a region of 110 square degrees in the Galactic plane, covered by VISTA Variables in the Via Lactea. A total of 641 candidate YSOs have been identified. We classified them into bursts, dips, faders, seculars, and irregulars. Within the burst category, 18 sources were identified as FUor candidates and 1 as an EXor candidate. Irregulars are the most prevalent in the sample. In both bursts and faders, the redder sources tend to show a pattern of bluer when brighter, whereas the bluer sources display the opposite trend, possibly related to the accretion structure of YSOs at different stages. Finally, we obtained the recurrence time-scale for FUor eruptions at various stages of YSO evolution. Our findings indicate that younger YSOs generally experience more frequent eruptions compared to older ones.

Mid-infrared Period–Luminosity Relations of Gaia DR3 Long Period Variables

Long period variable (LPV) stars are very promising distance indicators in the infrared bands. We selected asymptotic giant branch (AGB) stars in the Large and Small Magellanic Cloud (LMC and SMC) from the Gaia Data Release 3 LPV catalog, and classified them into oxygen-rich (O-rich) and carbon-rich (C-rich) AGB stars. Using the Wide-field Infrared Survey Explorer database, we determined the W1- and W2-band period–luminosity relations (PLRs) for each pulsation-mode sequence of AGB stars. The dispersion of the PLRs of O-rich AGB stars in sequences C′ and C is relatively small, around 0.14 mag. The PLRs of LMC and SMC are consistent in each sequence. In the W2 band, the PLR of large-amplitude C-rich AGB stars is steeper than that of small-amplitude C-rich AGB stars, due to their more circumstellar dust. By two methods, we find that some PLR sequences of O-rich AGB stars in the LMC are dependent on metallicity. The coefficients of the metallicity effect are β = –0.533 ± 0.213 mag dex1 and β = –0.767 ± 0.158 mag dex1 for sequence C in W1 and W2 bands, respectively. The significance of the metallicity effect in W1 band for the four sequences is 2.2−3.5σ. Both of these imply that distance measurements using O-rich Mira may need to take the metallicity effect into account.

Revisiting the Dust Torus Size–Luminosity Relation Based on a Uniform Reverberation-mapping Analysis

We investigate the torus size–luminosity relation of Type 1 active galactic nuclei (AGNs) based on the reverberation-mapping analysis using the light curves of the optical continuum and the IR continuum obtained with the W1 and W2 bands of the Wide-field Infrared Survey Explorer survey. The final sample consists of 446 and 416 AGNs, respectively, for W1- and W2-band light curves, covering a large dynamic range of bolometric luminosity from 1043.4 to 1047.6 erg s−1, which show reliable lag measurements based on our quality assessment analysis. After correcting for the accretion disk contamination in the observed IR flux, we constrain the torus size (R dust) and AGN bolometric luminosity (L bol) relationship with the best-fit slope of 0.39 (0.33) for the W1 (W2) band, which is shallower than expected from the dust radiation equilibrium model. By combining the previous K-band lag measurements, we find that the measured torus size depends on the observed wavelength of the dust radiation, as R dust,K : R dust,W1: R dust,W2 = 1.0:1.5:1.8 (R dust ∝ λ 0.80) at L bol = 1046 erg s−1, confirming a stratified structure of the torus, where wavelength-dependent emissions originate from distinct regions of the torus. By investigating the deviation from the best-fit torus size–luminosity relation, we find a moderate correlation between the offset from the R dust–L bol relation and Eddington ratio. This suggests a possible influence of the Eddington ratio on the observed flattening of the R dust–L bol relationship.

Background Contamination of the Project Hephaistos Dyson Spheres Candidates

Project Hephaistos recently identified seven M-dwarfs as possible Dyson Spheres (DS) candidates. We have cross-matched three of these candidates (A, B and G) with radio sources detected in various all-sky surveys. The radio sources are offset from the Gaia stellar positions by ∼4.9, ∼0.4 and ∼5.″0 for candidates A, B, and G respectively. We propose that Dust obscured galaxies (DOGs) lying close to the line-of-sight of these M-dwarf stars significantly contribute to the measured WISE mid-IR flux densities in the WISE W3 and W4 wave bands. These three stars have therefore been misidentified as DS candidates. We also note that with an areal sky density of 9 × 10−6 per square arcsecond, Hot DOGs can probably account for the contamination of all 7 DS candidates drawn from an original sample of 5 million stars.

A Reverberation Mapping Study of a Highly Variable Active Galactic Nucleus 6dFGS gJ022550.0-060145

We use LCOGT observations (MJD 59434−59600) with a total exposure time of ≃50 hr and a median cadence of 0.5 day to measure the interband time delays (with respect to u) in the g, r, and i continua of a highly variable active galactic nucleus (AGN), 6dFGS gJ022550.0-060145. We also calculate the expected time delays of the X-ray reprocessing of a static Shakura–Sunyaev disk according to the source's luminosity and virial black hole mass; the two parameters are measured from the optical spectrum of our spectroscopic observation via the Lijiang 2.4 m telescope. It is found that the ratio of the measured time delays to the predicted ones is 2.6−1.3+1.3 . With optical light curves (MJD 53650–59880) from our new LCOGT and archival Zwicky Transient Facility, Pan-STARRS, Catalina Sky Survey, and ATLAS observations, and Wide-field Infrared Survey Explorer (WISE) data (MJD 55214−59055), we also measured time delays between WISE W1/W2 and the optical emission. W1 and W2 have time delays (with respect to V), 9.6−1.6+2.9×102 days and 1.18−0.10+0.13×103 days in the rest frame, respectively; hence, the dusty torus of 6dFGS gJ022550.0-060145 should be compact. The time delays of the W1 and W2 bands are higher than the dusty torus size–luminosity relationship of Lyu et al. By comparing the infrared and optical variability amplitude, we find that the dust-covering factors of the W1 and W2 emission regions are 0.7 and 0.6, respectively. Future broad emission-line reverberation mapping of this target and the results of this work enable us to determine the sizes of the AGN main components simultaneously.

Stellar Atmospheric Parameters for Cool Dwarfs in Gaia Data Release 3

We provide a catalog of atmospheric parameters for 1,806,921 cool dwarfs from Gaia Data Release 3 (DR3) that lie within the range covered by LAMOST cool dwarf spectroscopic parameters: 3200 K < T eff < 4300 K, −0.8 < [M/H] < 0.2 dex, and 4.5 < log g < 5.5 dex. Our values are derived based on machine-learning models trained with multiband photometry corrected for dust. The photometric data comprise optical data from the Sloan Digital Sky Survey r, i, and z bands, near-infrared data from the Two Micron All Sky Survey J, H, and K bands, and mid-infrared data from the ALLWISE W1 and W2 bands. We used both random forest and light gradient boosting machine machine-learning models and found similar results from both, with an error dispersion of 68 K, 0.22 dex, and 0.05 dex for T eff, [M/H], and log g, respectively. Assessment of the relative feature importance of different photometric colors indicated W1 − W2 as most sensitive to both T eff and log g, with J − H being most sensitive to [M/H]. We find that our values show a good agreement with the Apache Point Observatory Galactic Evolution Experiment, but are significantly different to those provided as part of Gaia DR3.

An Episode of Occultation Events in Gaia21bcv

A previously unremarkable star near the Canis Major OB1/R1 association underwent an episode of multiple deep brightness minima. Light curves based on archival Gaia, Zwicky Transient Facility (ZTF), and NEOWISE data and additional observations from the Las Cumbres Observatory and UKIRT show that the star was not variable prior to 2019 August 18 (MJD 58700), and on that date started showing brightness dips of up to 3 mag in the Gaia G and ZTF r bandpasses. After MJD 59500, ≈800 days after the onset of these dipping events, the star returned to its previous brightness, and no significant dipping events have been recorded since. Compared to the stable phase, NEOWISE IR photometry in the W1 and W2 bands indicates a generally redder color, and both decreases and increases in brightness at different times during the dipping episode. The spectrum of Gaia21bcv taken after the end of the dipping episode shows several neutral and ionized metal absorption lines, including Li, indicating a spectral type of ≈K5. Variable emission from [O i] was observed. The Hα absorption in Gaia21bcv is too faint and irregular for this spectral type, indicating that the line is partly filled in by variable emission, a signature of weak episodic accretion. Gaia21bcv lies above the zero-age main sequence, but is much fainter than typical R CrB stars. We interpret the light curve of Gaia21bcv as being similar to the occultation events in ϵ Aurigae, i.e., occultation by a disk around a companion object orbiting the primary star.

Initial Results of our Spectro-photometric Monitoring of XZ Tau

We present here initial results of our spectro-photometric monitoring of XZ Tau. During our monitoring period, XZ Tau exhibited several episodes of brightness variations in timescales of months at optical wavelengths in contrast to the mid-infrared wavelengths. The color evolution of XZ Tau during this period suggests that the brightness variations are driven by changes in accretion from the disc. The mid-infrared light curve shows an overall decline in brightness by ∼ 0.5 and 0.7 magnitude respectively in WISE W1 (3.4 μm) and W2 (4.6 μm) bands. The emission profile of the hydrogen recombination lines along with that of Ca II IRT lines points towards magnetospheric accretion of XZ Tau. We have detected P Cygni profile in Hβ indicating of out-flowing winds from regions close to accretion. Forbidden transitions of oxygen are also detected, likely indicating jets originating around the central pre-main sequence star.

3 μm Phase Curves of Main-belt Asteroids from NEOWISE Photometry

The shapes of asteroid phase curves are influenced by the physical properties of asteroid surfaces. The variation of an asteroid’s brightness as a function of the solar phase angle can tell us about surface properties such as grain size distribution, roughness, porosity, and composition. Phase curves are traditionally derived from photometric observations at visible wavelengths, but phase curves using infrared data can also provide useful information about an asteroid surface. Using photometric observations centered near ∼3.4 μm from the W1 band of the Near-Earth Object Wide-field Infrared Survey Explorer mission, we construct thermally and rotationally corrected infrared phase curves for a sample of main-belt asteroids, which includes asteroids observed by the AKARI satellite, as well as subsets of the Themis and Flora dynamical families. We calculate the linear slope of the phase curves as a measure of their shape and compare W1 phase slopes to band depths of absorption features associated with hydrated materials, spectral slopes, visible albedos, W1 albedos, and diameters. We observe a steepening of the W1 phase slope of C-type asteroids with increasing 2.7 μm band depth but little correlation between the phase slope and 3 μm band depth or 3 μm spectral slope. The C-types in our sample exhibit steeper average W1 phase slopes than M- or S-types, similar to visible-light phase slopes. We also observe steeper W1 phase slopes for smaller-diameter objects within the Themis family and explore comparisons to Jupiter-family comets in phase slope versus albedo space.

The Structure Function of Mid-infrared Variability in Low-redshift Active Galactic Nuclei

Using the multi-epoch mid-infrared (MIR) photometry from the Wide-field Infrared Survey Explorer spanning a baseline of ∼10 yr, we extensively investigate the MIR variability of nearby active galactic nuclei (AGNs) at 0.15 < z < 0.4. We find that the ensemble structure function in the W1 band (3.4 μm) can be modeled with a broken power law. Type 1 AGNs tend to exhibit larger variability amplitudes than type 2 AGNs, possibly due to the extinction by the torus. The variability amplitude is inversely correlated with the AGN luminosity, consistent with a similar relation known in the optical. Meanwhile, the slope of the power law increases with AGN luminosity. This trend can be attributed to the fact that the inner radius of the torus is proportional to the AGN luminosity, as expected from the size−luminosity relation of the torus. Interestingly, low-luminosity type 2 AGNs, unlike low-luminosity type 1 AGNs, tend to exhibit smaller variability amplitude than do high-luminosity AGNs. We argue that either low-luminosity type 2 AGNs have distinctive central structures due to their low luminosity or their MIR brightness is contaminated by emission from the cold dust in the host galaxy. Our findings suggest that the AGN unification scheme may need to be revised. We find that the variability amplitude of dust-deficient AGNs is systematically larger than that of normal AGNs, supporting the notion that the hot and warm dust in dust-deficient AGNs may be destroyed and reformed according to the strength of the ultraviolet radiation from the accretion disk.

Compact object probed by microlensing in the universe: MACHOs, stars, planets?

Nearly 30% of the components in the universe are dark matter, hitherto, astronomers are still uncertain about their properties. This project attempts to constrain whether dark matter is MACHOs or diffused particles via statistics of microlensing events observed toward quasars, which are the brightest objects in the Universe. The identification of isolated microlensing events enables us to study stellar and planetary objects in distant universes that would be otherwise unobservable. During the research two isolated microlensing events toward quasar J1821 is discovered in the first 25 quasars in our sample. The WISE light curves in W1 and W2 bands are almost identical in magnification, consistent with achromatic variability due to microlensing. The crossing times for the two events are about 2 days and 1.3 days, thus the lenses are probably a star or sub-stellar object. By assuming the lenses are in galaxies of the galactic cluster with a redshift of 0.81 lying in our sight of view to the quasar, calculations of the mass of the lenses as a function of its velocity can be made, yielding a planetary mass or less in the relative velocity dominated by the movements of galaxies. If the superluminous motion of jets in a quasar dominates the relative velocity, stellar masses are derived. The data analysis yielded from the microlensing light curves of J1821 suggests microlensing as a method for astronomers to further study stellar and sub-stellar objects in our Universe.

Optical properties of metal-poor T dwarf candidates

Context. Metal-poor brown dwarfs are poorly understood because they are extremely faint and rare. Only a few candidates have been identified as T-type subdwarfs in infrared surveys and their optical properties remain unconstrained. Aims. We aim to improve the knowledge of the optical properties of T subdwarf candidates to break the degeneracy between metallicity and temperature and to investigate their atmospheric properties. Methods. Deep z-band images of ten known T subdwarf candidates were collected with the 10.4-m Gran Telescopio Canarias. Low-resolution optical spectra for two of them were obtained with the same telescope. Photometric measurements of the z-band flux were performed for all the targets and they were combined with infrared photometry in J, H, K, W1, and W2 bands from the literature to obtain the colours. The spectra were compared with solar-metallicity T dwarf templates and with laboratory spectra. Results. We found that the targets segregate into three distinct groups in the W1 − W2 versus z − W1 colour-colour diagram. Group I objects are mixed with solar-metallicity T dwarfs. Group III objects have W1 − W2 colours similar to T dwarfs but very red z − W1 colours. Group II objects lie between Group I and III. The two targets for which we obtained spectra are located in Group I and their spectroscopic properties resemble normal T dwarfs but with water features that are deeper and have a shape akin to pure water. Conclusions. We conclude that the W1 − W2 versus z − W1 colour-colour diagram is excellent to break the metallicity-temperature degeneracy for objects cooler than L-type ones. A revision of the spectral classification of a T subdwarf might be needed in the future, according to the photometric and spectroscopic properties of WISE1810 and WISE0414 in Group III discussed in this work.

Constraining AGN Torus Sizes with Optical and Mid-infrared Ensemble Structure Functions

We propose a new method to constrain the size of the dusty torus in broad-line active galactic nuclei (AGNs) using optical and mid-infrared (MIR) ensemble structure functions (SFs). Because of the geometric dilution of the torus, the MIR response to optical continuum variations has suppressed variability with respect to the optical that depends on the geometry (e.g., size, orientation, opening angle) of the torus. More extended tori have steeper MIR SFs with respect to the optical SFs. We demonstrate the feasibility of this SF approach using simulated AGN light curves and a geometric torus model. While it is difficult to use SFs to constrain the orientation and opening angle due to the insensitivity of the SF on these parameters, the size of the torus can be well determined. Applying this method to the ensemble SFs measured for 587 SDSS quasars, we measure a torus R–L relation of in the WISE W1 band and sizes ∼1.4 times larger in the W2 band, which are in good agreement with dust reverberation mapping measurements. Compared with the reverberation mapping technique, the SF method is much less demanding in data quality and can be applied to any optical+MIR light curves for which a lag measurement may not be possible, as long as the variability process and torus structure are stationary. While this SF method does not extract all information contained in the light curves (i.e., the transfer function), it provides an intuitive interpretation for the observed trends of AGN MIR SFs compared with optical SFs.

Mid-infrared dusty torus sizes in active galactic nuclei with Hβ reverberation mapping

ABSTRACT We compile optical and mid-infrared light curves from the time-domain surveys (i.e. CRTS, PTF, ZTF, and ASAS-SN) and Wide-field Infrared Survey Explorer (WISE) archive for a selected sample of active galactic nuclei (AGNs) with Hβ reverberation mapping (RM) measurements. We measure the time lags (and thus torus sizes) of W1 (∼3.4 $\mu$ m) and W2 (∼4.6 $\mu$ m) band light curves relative to the optical one using the MICA method. Through Hβ RM, the sample has well-measured AGN properties, therefore allowing us to reliably constrain the relations between torus sizes and AGN properties. We perform linear regressions for the relations between torus sizes and 5100 Å luminosities ($R\propto L_{5100}^{\beta }$ ) in two cases: β = 0.5 and β set free. The latter case yields β ≈ 0.37 ± 0.028 for both W1 and W2 bands, shallower than the expected value of 0.5, possibly due to the dependence of torus size on accretion rate. For β = 0.5, by combining with the previous K band RM measurements, we obtain the characteristic broad-line region (BLR) and tours sizes following RBLR:RK :RW1:RW2 = 1.0:6.2:9.2:11.2. We investigate the deviations of the W1 and W2 band observed torus sizes from the corresponding best-fitting relations (with β = 0.5) and find that they both are correlated with accretion rate. As the accretion rate increases, the torus sizes tend to be shortened compared to the anticipated sizes from the best-fitting relations, similar to the behaviour found in BLRs. Such behaviours can be explained by the self-shadowing effect of slim discs. This is further supported by ratios of the W1 and W2 band torus sizes to BLR sizes, which do not show significant correlations with AGN properties.

RR Lyrae Mid-infrared Period–Luminosity–Metallicity and Period–Wesenheit–Metallicity Relations Based on Gaia DR3 Parallaxes

We present new empirical infrared period–luminosity–metallicity (PLZ) and period–Wesenheit–metallicity (PWZ) relations for RR Lyae based on the latest Gaia Early Data Release 3 (EDR3) parallaxes. The relations are provided in the Wide-field Infrared Survey Explorer (WISE) W1 and W2 bands, as well as in the W(W1, V − W1) and W(W2, V − W2) Wesenheit magnitudes. The relations are calibrated using a very large sample of Galactic halo field RR Lyrae stars with homogeneous spectroscopic [Fe/H] abundances (over 1000 stars in the W1 band), covering a broad range of metallicities (−2.5 ≲ [Fe/H] ≲ 0.0). We test the performance of our PLZ and PWZ relations by determining the distance moduli of both galactic and extragalactic stellar associations: the Sculptor dwarf spheroidal galaxy in the Local Group (finding ), the Galactic globular clusters M4 (), and the Reticulum globular cluster in the Large Magellanic Cloud (). The distance moduli determined through all our relations are internally self-consistent (within ≲0.05 mag) but are systematically smaller (by ∼2–3σ) than previous literature measurements taken from a variety of methods/anchors. However, a comparison with similar recent RR Lyrae empirical relations anchored with EDR3 likewise shows, to varying extents, a systematically smaller distance modulus for PLZ/PWZ RR Lyrae relations.

Evidence for Black Holes in Green Peas from WISE Colors and Variability

We explore the presence of active galactic nuclei (AGNs)/black holes in Green Pea galaxies (GPs), motivated by the presence of high-ionization emission lines such as He ii and [Ne iii] in their optical spectra. In order to identify AGN candidates, we used mid-infrared (MIR) photometric observations from the all-sky Wide-field Infrared Survey Explorer (WISE) mission for a sample of 1004 GPs. Considering only >5σ detections with no contamination from neighboring sources in AllWISE, we select 31 GPs out of 134 as candidate AGNs based on a stringent three-band WISE color diagnostic. Using multi-epoch photometry in W1 and W2 bands based on time-resolved unWISE coadd images, we find two sources exhibiting variability in both the WISE bands among 112 GPs with W1 ≤16 mag and no contamination from neighboring sources in unWISE. These two variable sources were selected as AGNs by the WISE three-band color diagnostic as well. Compared to variable AGN fractions observed among low-mass galaxy samples in previous studies, we find a higher fraction (∼1.8%) of MIR variable sources among GPs, which demonstrates the uniqueness and importance of studying these extreme objects. Through this work, we demonstrate that MIR diagnostics are promising tools to select AGNs that may be missed by other selection techniques (including optical emission-line ratios and X-ray emission) in star formation-dominated, low-mass, low-metallicity galaxies.

Galaxy and Mass Assembly (GAMA)

Aims. We investigate how different mid-infrared (mid-IR) properties of galaxies are correlated with the environment in which the galaxies are located. For this purpose, we first study the dependence of galaxy clustering on the absolute magnitude at 3.4 μm and redshift. Then, we look into the environmental dependence of mid-IR luminosities and the galaxy properties derived from these luminosities. We also explore how various IR galaxy luminosity selections influence the galaxy clustering measurements. Methods. We used a set of W1 (3.4 μm) absolute magnitude (MW1) selected samples from the Galaxy and Mass Assembly (GAMA) survey matched with mid-IR properties from the Wide-field Infrared Survey Explorer (WISE) in the redshift range 0.07 ≤ z &lt; 0.43. We computed the galaxy two-point correlation function (2pCF) and compared the clustering lengths between subsamples binned in MW1 and in redshift. We also measured the marked correlation function (MCF), in which the galaxies are weighted by marks when measuring clustering statistics, using the luminosities in the WISE W1 to W4 (3.4 to 22 μm) bands as marks. Additionally, we compared the measurements of MCFs with different estimates of stellar mass and star formation rate (SFR) used as marks. Finally, we checked how different selections applied to the sample affect the clustering measurements. Results. We show strong clustering dependence on the W1 absolute magnitude: galaxies brighter in the W1 band are more strongly clustered than their fainter counterparts. We also observe a lack of significant redshift dependence of clustering in the redshift range 0.07 ≤ z &lt; 0.43. We show that although the W1 and W2 bands are direct indicators of stellar mass, a galaxy sample selected based on W1 or W2 bands does not perfectly show the clustering behaviour of a stellar mass-selected sample. The proxy relation between W3 and W4 bands and SFR is similar. We also demonstrate the influence of estimation techniques of stellar mass and SFR on the clustering measurements.