Adventures with Maffei 1 & 2

T. Jarrett, IPAC
(971028)

Our fearless leader M. Skrutskie was kind enough to authorize early observation of the notable galaxy pair discovered by Maffei over 30 years ago. These galaxies, known as Maffei 1 and 2, are unique in the sense that they are nearby (thus big and bright) but are vailed by the disk of our Galaxy. On the optical POSS images they appear as faint smudges -- reminiscent of low surface brightness galaxies. This is what Maffei saw way back in the days of the Beatles.

Now with our infrared eyes, we see something quite different. Roc has composed a nice 3-band JHK image of the pair. Check out Maffei 1 and Maffei 2 . Since 2MASS will survey the entire Galactic plane, we now have the means to find additional (if any) galaxies like the Maffei twins. So for now, the Maffei galaxies represent a sort of benchmark -- Can GALWORKS detect them? The answer is Yes. But the story is rather interesting. Stay tuned.

Maffei 1 had a very different detection path from that of Maffei 2, primarily because Maffei 2 is in my input catalog of known galaxies but Maffei 1 is not. The reason is that the catalog has a size limit of about 1 arcmin in diameter. Maffei 1 does not satisfy this limit in the optical (in the IR it is booming bright, as we shall see). In retrospect the catalog should probably be expanded to include smaller galaxies. I will pursue this angle at another time. In any case, Maffei 2, being in the catalog, was automatically detected and carried along in GALWORKS as a known entity. It was processed and parameterized in the usual fashion, with the information and stamp images written to disk. The stamp images are shown later in this memo. In contrast, Maffei 1 had a more unusual path.

At the onset, Maffei 1 was already different from any other galaxy encountered by GALWORKS thus far (note: relatively speaking, this is a miniscule sample): It is so bright in the core (K < 8) that it passes the flux thresholds for bright star blanking. That is to say, GALWORKS initially thought it was a bright star. The object, accordingly, was blanked from the coadd images during the bright star masking procedure. This was not a good turn of events. It was also not anticipated. Why? Because any galaxy that is bright enough to look like a bright star (to the point source processor) should in principle be big enough to be included in the catalog of known galaxies (like Maffei 2; see above). But Maffei 1 exposes the weakness in this assumption. It is big and bright, but the catalog does not know this (it is primarily an optical sample) -- thanks to Galactic extinction. There are undoubtedly other "gotchas" that I have not anticipated.

Maffei 1 is not a bright star, it is an extended object. It exhibits a "confusion radius" signature quite different from that of bright stars. GALWORKS applies a set of algorithms to find the optimum image masking of bright stars; see Automated Bright Star Masking . The first component of bright stars that is computed is the confusion radius -- corresponding to the PSF of the bright star plus residual (halo) emission symmetrically distributed about the star. The method is to compute the median value in annuli extended from the star, with the confusion radius set at about the 0.5 to 1-sigma level of the sky background. The confusion radius, however, cannot grow beyond 2.5 times the initial guess of the confusion radius (i.e., we apply a boundary limit). This limit should leave plenty of room for the confusion radius to converge. The exception occurs when either (1) the initial flux of the bright star is so far off that the initial guess at the confusion radius is severely underestimated (this can happen with very bright stars that saturate leaving a very poorly determined R1 mag), or (2) the object is not really a bright star (i.e., the initial flux of the object -- a pt. source estimated flux -- bears no resemblance to the true flux of the object). The latter case applies to Maffei 1. The confusion radius did not converge for this object (there was plenty of galaxy beyond 2.5 * estimated confusion radius). So, the obvious modification to the Automated Bright Star Masking routine was to set three confusion radii flags (one for each band), 1 for convergence and 0 for non-convergence. We then do not mask any object that does not converge. Maffei 1 satisfies this last conditional -- it is not blanked as a bright star.

We have happily solved the problem of Maffei 1 being too bright (and landing amongst the bright star blanking set). This is not to say that all galaxies that have initial K < 9 (or J < 10) will fail to have a confusion radius that converges. We can imagine core bright galaxies (like seyfert I types) that will satisfy the bright star masking algorithm thresholds. The solution for these galaxies is unclear at this time (but we note that they must be rare because they have yet to be encountered).

With Maffei safely protected from brutal masking, the object was successfully detected and parameterized by GALWORKS. Let the sunshine in. The images below show the postage stamps of the galaxy pair. Note the large number of stars subtracted from the images (red circles) -- the stellar number density is very high in this field. Also notice that the integrated flux for Maffei 1 is about 5th mag!! Maffei 2 is about 6th mag. Still, these isophotal fluxes severely underestimate the true total flux of the galaxies because they measure only the photons coming from the inner 80 pixels (in radius) -- the GALWORKS working limit. The galaxy pair is considerably bigger than 2 arc min in radius. The third image in the gif below is a star that is located in the outer halo of Maffei 2 (the bigger of the pair).

Another interesting thing that happened to Maffei in its journey through GALWORKS was to land in the Algorithm 3 grinder -- EXTENDED BRIGHT STAR DETECTION. Algor 3 (which has not been discussed since the CDR, over 1 year ago) was designed to isolate bright stars that exhibit emission beyond the confusion radius limit (e.g., H II regions, planetary nebulae, etc). The algor is discussed briefly in the Galworks SDS. The relevant section is reproduced here. It requires some modification, however, since I have added information with regard to the confusion radius operation. (additions are coming soon)


    EXTENDED BRIGHT STAR DETECTION (as of march, 1996)

    It is possible to detect emission surrounding bright stars if the star and its associated reflected light (e.g., diffraction spikes) are properly subtracted from the coadd first. The algorithm described here (a.k.a algorithm 3) is designed to find bright emission around (but beyond the confusion radius of) bright stars.

      1. Sort bright list according to flux
      2. Select star fainter than K = 6th mag.
      3. subimage: 160 X 160
      4. blank all nearby brighter objects; blank radius = predicted confusion radius (for J, blank a larger region to annex the "ghost" seen near the core); also, blank their diffract spikes (used predicted values)
      5. blank all nearby fainter stars (but still bright stars); blank radius = predicted confusion radius
      6. blank area centered on bright source, radius = 90% of the predicted confusion radius
      7. blank all sources with a local maxima and peak flux > 3*sigma; blank radius = 6
      8. blank diffract spikes (use predicted spike lengths)
      9. compute stats for surviving pixels in subimage
        mean and standard deviation of the mean,
        standard deviation of mean from zero (mean - 0)
        >> (submean, subsigma_1, subsigma_2)
      10. parameters:
        sub_ratio1 = 100 * [ (subsigma_1 / sigma) - 1]
        sub_ratio2 = 100 * [ (subsigma_2 / sigma) - 1]

        subd = max ( sub_ratio1,sub_ratio2)

      10. Apply criteria
        subd > 25
        confusion radius flags (no convergence == pass threshold)

Maffei 1 is clearly a good candidate for this algorithm. It is bright (estimated point source mag, just less than 9th at K), but extending well beyond any predicted confusion radius. And indeed, Algor 3 easily finds Maffei 1, as well as three additional bright stars located within its massive halo. The parameters computed by algor 3, the sub_ratios, all have values well beyond the threshold for candidacy (in fact, the values "max out"). In algor 3 phase space, Maffei 1 is exceptional. The following gif image shows the postage stamp of the field "extracted" by algor 3 -- it is Maffei 1. The photometry corresponds to the R1 mags (i.e., the initial flux estimates of this object; notice how far off the flux is to the true flux). The sub ratios are called "Jsubd", "Hsubd" and "Ksubd", all maxed out in the case of Maffei 1. The string "0 0 0" refers to the confusion radius flags, 0 means that the confusion radius algorithm did not converge.