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)
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.
subd = max ( sub_ratio1,sub_ratio2)
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.
EXTENDED BRIGHT STAR DETECTION (as of march, 1996)
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,
10. parameters:
standard deviation of mean from zero (mean - 0)
>> (submean, subsigma_1, subsigma_2)
sub_ratio1 = 100 * [ (subsigma_1 / sigma) - 1]
10. Apply criteria
sub_ratio2 = 100 * [ (subsigma_2 / sigma) - 1]
subd > 25
confusion radius flags (no convergence == pass threshold)