the giant galactic globular cluster @xmath3 centauri ( ngc 5139 , hereafter @xmath3 cen ) is not only the most massive globular cluster in the milky way , but its stars have also been shown to have a large spread in both age and metallicity ( see for instance , many contributions in van leeuwen et al .
2002 ) . in light of these discoveries ,
it has been proposed ( lee et al .
1999 , majewski et al . 2000 )  and now dynamically modeled ( tsujimoto & shigeyama 2003 )  that @xmath3 cen may be the nucleus of a disrupting dsph , similar to the galactic globular cluster m54 . the second largest globular cluster , m54 is currently believed to be the core of the sgr dsph ( da costa & armandroff 1995 , bassino & muzzio 1995 , sarajedini & layden 1995 , layden & sarajedini 2000 ) and lies at a nucleation point for the several different metallicity populations of sgr ( sarajedini & layden 1995 , majewski et al . 2003 ) .
the recent observation of apparent tidal tails around @xmath3 cen by leon , meylan & combes ( 2000 , hereafter lmc00 ) bolsters the notion of this cluster being the remnant of a larger , now disrupted parent object , and , with @xmath3 cen and m54 as models , ree et al .
( 2002 ) , lee et al .
( 2001 , 2002 ) , and yoon et al .
( 2000 ) have proposed that other clusters with evidence of multiple populations , such as clusters with bimodal horizontal branches , may also be the nuclei of disrupting dsph satellite galaxies .
two clusters proposed to be m54 analogues ( ree et al .
2002 ) are ngc 6388 and ngc 6441 , the third and fifth largest globular clusters in the milky way respectively , since both have observed characteristics suggesting internal metallicity spreads . moreover , high mean metallicity clusters like these ( see table 1 ) should have only a red stub of a horizontal branch ( hb ) , but ngc 6388 and ngc 6441 also have extended blue hb tails ( piotto et al .
1997 , hereafter p97 ) .
a _ bimodal _ ( or extended ) metallicity distribution may explain both this curious hb morphology ( p97 ) as well as the unusually wide red giant branches ( rgb ) observed in these clusters .
because it is difficult to understand how multiple generations of stars might be created in a stellar system of initial mass of order @xmath4 ( which is not massive enough to retain significant amounts of gas after an initial starburst ) , it is possible that the multi - population clusters are the remnants of _ formerly _ larger stellar systems , for example dsph galaxies , where multi - epoch star formation histories are commonly found .
the verification of the existence of additional multiple - population clusters would increase the need to formulate a general process  e.g. , the disruption of initially larger systems  that leads to a more substantial _ class _ of these objects .
these scientific questions are complicated by the presence of non - negligible differential reddening in the fields of some the most interesting examples of potential multi - modal clusters . in this paper , we concentrate on two effects of differential reddening relevant to the globular cluster evolution picture outlined above : the potential creation of false signatures of both tidal tails around , and of multiple metallicities within , globular clusters . these issues are explored in case studies of the three massive globular clusters @xmath3 cen , ngc 6388 and ngc 6441 .
tidal tail searches around galactic satellites can be particularly vulnerable to differential reddening effects because search techniques typically rely on the detection of star count overdensities
. lines of sight with a lower dust column density will naturally reveal more stars to a given magnitude limit .
this may create an apparent overdensity of stars in regions with lower extinction that could be mistaken for , or overemphasize , density enhancements from tidal tails . in the case of the reported tidal tails around @xmath3 centauri by lmc00
, the contours of the detected tails curiously and closely correspond to contours of the prominent iras @xmath5 dust emission ( and therefore , presumably , optical extinction by dust ) in this field , as these authors have observed . because lmc00 did not correct for differential extinction in their optical starcounts ( s. leon , private communication ) ,
at minimum there must be some dust - shaping of the morphology of mapped @xmath3 cen tidal tails .
we re - evaluate the possible impact of differential reddening on large field studies of @xmath3 cen by comparing a customized tidal tail search on raw 2mass data to one performed on 2mass data differentially dereddened according to the dust maps of schlegel et al .
( 1998 ) .
the cause for the observed width of the rgb in ngc 6388 and ngc 6441 has been variously explained as differential reddening and/or the presence of multiple populations of different metallicity , with no clear consensus emerging .
p97 originally concluded that differential reddening on scales greater than eight arcseconds was unlikely for ngc 6388 but possible for ngc 6441 , although only to a few hundredths of a magnitude .
later evidence indicating the presence of strong differential reddening on small scales across both clusters was found by raimondo et al .
( 2002 ) using hst - wfpc2 small area surveys similar to those used by p97 .
pritzl et al .
( 2001 , 2002 ) also find differential reddening in ngc 6441 and ngc 6388 using ctio @xmath6 and @xmath7 band photometry , although they note that in ngc 6388 the differential reddening probably covers a smaller magnitude range .
if differential reddening were significant , the reddening vector could reproduce the sloping red hbs and spread the rgb stars perpendicular to the unreddened rgbs , inflating their apparent width .
numerical experiments performed by raimondo et al .
( 2002 ) indicate that differential reddening of order @xmath8e@xmath9 0.1 would be substantial enough to cause an unreddened hb of similar metallicity ( e.g. , that of 47 tucanae ) to resemble a sloping red stub similar to that of ngc 6388 or ngc 6441 .
nevertheless , raimondo et al .
conclude that reddening is not the reason for the strange hb s since their sloping nature is present even in small area subsamples exhibiting a thin rgb , which suggests that a second effect is at work .
pritzl et al .
( 2001 , 2002 ) and ree et al .
( 2002 ) argue that variable metallicity may adequately explain the observed features of the ngc 6388 and ngc 6441 cmds , and recent stellar evolutionary models ( e.g. , sweigart & catelan 1998 ) confirm that a metallicity spread can naturally create a sloping hb .
specifically , ree et al .
( 2002 ) find that two distinct populations separated by 1.2 gyr and 0.15 dex in [ fe / h ] can reproduce both the observed wide rgbs and tilted hbs in the clusters .
however , pritzl et al . ( 2001 , 2002 ) argue that the simple two - population model fails to reproduce the observed hb slope and that a more general metallicity _ spread _ is necessary .    because selective extinction at the 2mass wavelengths @xmath0 ( 1.25 @xmath10 ) , @xmath1 ( 1.65 @xmath10 ) , & @xmath11 ( 2.17 @xmath10 ,
hereafter simply referred to as @xmath12 ) is roughly one tenth that in visual bands , the magnitude of differential reddening at these wavelengths is commensurately reduced .
thus , we apply the 2mass point source catalog to test whether : ( 1 ) differential reddening affects interpretations of an _ infrared _ search for tidal tails around @xmath3 cen , and ( 2 ) the ngc 6388 and ngc 6441 rgbs are broadened at near infrared wavelengths where differential reddening will have less impact on the width of the cluster rgb s .
@xmath0 , @xmath1 , and @xmath2 band data were drawn from the 2mass all - sky point - source working survey database as of july , 2002 .
initial values of the mean field extinction , @xmath13 , and distance modulus , @xmath14 , for each cluster were obtained from the harris compilation as posted on the world wide web ( http://physun.physics.mcmaster.ca/globular.html ) on june 22 , 1999 ( see harris 1996 ) .
these were converted to the required infrared extinction parameters using the standard extinction models of cardelli , clayton , and mathis ( 1989 ) .
table 2 gives a complete list of the adopted extinction parameters .
the lmc00 study of optical starcounts in the @xmath3 cen field makes no corrections for differential reddening because magnitude calibrations were not available for the majority of clusters in their study ( s. leon , private communication ) . as figures 1a and 1b illustrate ,
the proposed tidal tails of @xmath3 cen are from enhanced star counts whose pattern appears to follow regions of lower dust absorption on the iras @xmath5 map .
this effect is particularly noticeable on the western edge of the cluster .
lmc00 note and discuss the possibility that the star counts are affected by reddening , but determine that the tidal tails are of a sufficiently large s / n that differential reddening should not significantly affect their tidal tail search .    by comparing starcount patterns obtained using raw 2mass data with those from 2mass data
differentially dereddened according to the dust maps of schlegel et al .
( 1998 ) we intend to re - evaluate whether uncorrected differential extinction can produce a false tidal tail signal .
the 2mass data has a @xmath15 completeness limit at @xmath16 mag in the direction of @xmath3 cen , yielding almost 500,000 stars within a @xmath17 box centered on the cluster .
because a photometric search for diffuse tidal tails presents considerable challenges , not least of which is substantial field star contamination in this low latitude field , a two stage method is employed to minimize field star contamination and probe for starcount excesses as deeply into the surrounding field as possible : first , a coarse filtering mask in color - color - magnitude ( ccm ) space is used to extract only those stars whose ccm location is similar to that of stars observed at radii less than the cluster core radius .
although this coarse mask does not take into account foreground / background contamination of the cluster core , it filters out a significant number of stars whose ccm locations are significantly different from those in the cluster .
second , we then apply a signal - to - noise filter similar to that used by odenkirchen et al .
( 2001 ) to separate cluster stars from field stars and emphasize starcount overdensities in the potential tidal tails .    in the first stage of filtering
, a filtering mask is applied in ccm space to pre - select those ccm locations containing stars within 100 arcseconds of the center of the cluster ( approximately @xmath18 the core radius @xmath19 , see table 1 ) .
we divide this irregularly shaped region in ( @xmath0 , ( @xmath20 ) , ( @xmath21 ) ) space into a @xmath22 grid in ccm space bounded by the limits of the ccm region ; this grid size is selected to achieve a balance between coarse binning and low number statistics . with the individual cells described by indices ( @xmath23 , @xmath24 , @xmath25 ) respectively
, we define @xmath26 as the number of stars within 100 arcseconds of the cluster center that fall within each grid cell in ccm space . to obtain an estimate of the number of field stars included in @xmath26
, we represent the field star distribution , @xmath27 , by the stars within an annulus @xmath28 ( roughly three times the tidal radius of the cluster ) .
we assume that the distribution of field stars is relatively even across the cluster and scale the field mask by the ratio @xmath29 of the areas over which the field and core regions were summed .
the selection mask @xmath30 is defined as the ratio of the core mask to the sum of the core and scaled field masks .
@xmath31 @xmath32 will have a value near one for regions of ccm space where the core sample @xmath26 is not substantially contaminated with field stars , and will decrease as contamination becomes more severe .
the complete data set is filtered using this ccm selection mask : a star is rejected from the data set if the corresponding @xmath32 is less than some minimum value @xmath33 . this algorithm is a slightly more sophisticated means to eliminate stars unrelated to @xmath3 cen than simple preliminary cuts in the color - magnitude ( cm ) plane . we find that a value of @xmath34 allows us to select cluster stars over field stars most optimally and appears to preserve the cmd of the cluster center best .
this preliminary filtering rejects roughly 90% of the stars in the original data set .    in the second stage of filtering , we apply a modified version of grillmair et al.s ( 1995 ) method in the cm plane as described by odenkirchen et al .
this second stage eliminates field stars more effectively than the first stage , but we find that the method does not perform well unless a coarser pre - selection of stars is performed first to reduce the data set to a size that the more computationally complex second - stage routine can handle . with some recycling of notation
, we define a grid in @xmath35 space , with indices ( @xmath24 , @xmath25 ) respectively , and find the signal- to - noise ratio ( @xmath36 ) of probable cluster stars to field stars in each cm grid cell .
we count the number of stars @xmath37 that fall within each grid cell and lie within @xmath38 of the cluster center ( chosen to correspond to the region in which lmc00 observe circular cluster surface density contours , see fig .
1a ) , and the number of stars @xmath39 that fall within the cm grid cell and lie between @xmath40 and @xmath41 arcminutes from the cluster center . the field star population is assumed to be uniform , and the expected number of field stars of a given ( @xmath24 , @xmath25 ) that fall within @xmath38 of the cluster core is subtracted from @xmath37 . with @xmath29 representing the ratio of the area of sky on which @xmath42 and @xmath43 respectively have been defined , the local s / n ratio of cluster to field stars for each grid cell ( @xmath24 , @xmath25 ) is , by standard poissonian statistics ,    @xmath44    it is possible to use @xmath45 as a mask for selection in the cm plane by eliminating all @xmath46 grid cells with @xmath45 less than some limiting value @xmath47 , which is chosen such that the number of stars in the tidal tail region ( as reported by lmc00 ) is a maximum relative to the field star population .
the new matrix @xmath48 represents the total number of stars in each @xmath46 grid cell for the region and outer radius @xmath49 ) located directly north and south of the circle described by the inner radius .
] of the cluster in which tidal tails have been reported by lmc00 ( see fig .
1a ) , as shown by the dashed lines in fig .
starting from @xmath50 ( where @xmath51 is the maximum value of all @xmath45 ) and iterating downward to @xmath52 , we calculate @xmath53 and @xmath54 for all ( @xmath24 , @xmath25 ) such that @xmath55 .
hence , @xmath56 and @xmath57 respectively are the total number of field stars and cluster stars in lmc00 s tidal tail regions whose locations in @xmath58 space satisfy @xmath55 .
the optimal value of @xmath47 is that which maximizes @xmath59 ; hence this final selection mask filters the data set for all stars which fall within cells in cm space such that @xmath60 .
the positions of stars from this filtered set are then converted to a smooth density distribution by application of a two - dimensional gaussian kernel with an adaptive half maximum width ( e.g. , silverman 1986 ) for each star set to the distance of the 30@xmath61 nearest neighbor in the plane of the sky .
each of these gaussians is evaluated at finely spaced grid points in the plane of the sky , and these values are summed to create an image of relative stellar densities at each location on the grid ( fig . 1c and 1d ) .
we find that this @xmath36 optimization method worked efficiently to reduce the 2mass data set by a further 80 - 90% and brought out starcount overdensities around the cluster that have the appearance of tidal tails .
focusing the definition of @xmath62 on the regions found by lmc00 to possess tidal tails is considerably more effective at bringing these apparent tails out of the background than defining @xmath62 more generally ; however , this particular approach obviously is not useful for a general survey for previously unknown tidal tail searches because it requires _ a priori _ knowledge of where tails may be found .
we verified that the overdensities found by this method are not merely artifacts of the definition of @xmath62 in the region of the lmc00 expected tails by redefining @xmath62 in locations rotated around the cluster core by 90 degrees ( that is , protruding east - west as opposed to north - south ) . using this definition of @xmath62
, no starcount overdensities were detected around the cluster , which supports that the apparent tails found with the first @xmath62 definition had not been an artificial product of the processing method and that the method of sampling @xmath62 in the tail region is a crucial step in bringing those tails out of the background noise .
in the analysis of the raw 2mass data described above a starcount excess to the north and south of the cluster is found . as fig .
1c shows , to within the limits imposed by the smoothing kernel , the location of the 2mass tails corresponds to that of the tails found by lmc00 .
the same analysis has been repeated with 2mass data dereddened according to the extinction maps of schlegel et al .
the schlegel et al .
maps have a resolution of approximately 6.1 arcminutes , which corresponds to about 1500 cells of reddening data over a @xmath63 field of view , and show differential reddening of order @xmath64 mag over that field of view .
in contrast to what is found for the raw 2mass data , the results of the analysis on the dereddened 2mass data do not reveal any elongation of the circular cluster contour lines or any overdensities resembling tidal tails ( fig .
2 plots radial surface density profiles for both original and dereddened 2mass data , and shows the surface density of the dereddened cluster to fall off smoothly with radius , like a king ( 1966 ) profile .
in contrast , the surface density of stars in the non - dereddened cluster data appears to drop at @xmath65 ( @xmath66 60 arcmin )  corresponding to the low starcount regions at this radius to the east and west of the cluster  and rise again at @xmath67 ( @xmath68 arcmin )  which corresponds to high starcount regions to the north and south of the cluster ( fig . 1c ) .
given the fidelity with which the apparent infrared and optical tails follow the regions of lower reddening on the iras 100 micron map ( fig .
1b , particularly on the western and southern edges of the cluster ) , and the disappearance of the infrared tails in the dereddened cluster data , we speculate that by analogy the @xmath3 cen tidal tails previously detected with the optical starcounts , which are more affected by dust extinction , could have been artifacts of differential reddening .
we must also point out , however , that our test of the 2mass database is not completely analagous to the lmc00 analysis : due to the @xmath16 mag completeness limit of 2mass we are unable to probe to the same effective depth as the lmc00 study .
the differences in relative depth are demonstrated by the fact that while lmc00 observe reasonably circular contours out to a radius of about 45 arcminutes ( roughly the tidal radius , see table 1 ) , we only find circular contours out to about 40 arcminutes before we begin to observe distinct elongation of the contours into potential tidal tails ( fig .
1 ) . thus , the leon et al .
counts reach fainter surface brightness limits than does 2mass and are sensitive to more diffuse , lower density tails than are we .
nevertheless , we regard the 2mass test presented here as a cautionary example demonstrating that a proper accounting of differential reddening is an important step in mapping any tidal tails in the @xmath3 cen ( or any other heavily obscured ) field .
to investigate the effects of differential reddening on the width of the rgb in the clusters ngc 6388 and ngc 6441 we use the uniform photometry of 2mass to look for spatial patterns in the distribution of red giant star colors within each cluster .
47 tucanae ( ngc 104 , hereafter 47 tuc ) is adopted as a control case for a cluster of similar mean metallicity but with little differential reddening ( average @xmath69 0.02 ) .
we also compare our results to the @xmath3 cen data  although this cluster has an overall [ fe / h ] much lower than the other clusters , it nonetheless provides a useful example of the effects of _ variable metallicity _ within a cluster over spatial scales with relatively little differential reddening on the scale of arcminutes ( @xmath70 mag , note that the previous section explored extinction variations in @xmath3 cen on the scale of _ degrees _ ) .
analysis of rgb color _ spreads _ at different magnitudes in one cluster as well as differences between clusters is simplified if we remove the mean rgb color - magnitude relation for each cluster ( i.e. make the center of the rgb locus for each cluster appear vertical in the color - magnitude diagram ) .
since the [ fe / h ] of the unreddened cluster 47 tuc is similar to that of ngc 6388 and ngc 6441 to within about 0.25 dex ( see table 1 ) , the shape , height and position of the rgbs for the three clusters in a reddening corrected , absolute color - magnitude space should be similar .
hence , we use the well - defined 2mass cmd of 47 tuc to `` straighten '' the rgb of the clusters in our sample ( figure 3a ) .
we convert apparent to absolute magnitudes using the distance moduli in table 2 and consider only that part of the rgb with absolute magnitude @xmath71 to minimize contamination from cluster horizontal branch ( hb ) as well as field stars .
a legendre polynomial is fit to the 47 tuc rgb , with a rejection threshhold of 2@xmath72 iterated until convergence .
the final polynomial is converted to power series coefficients up to third order , which are used to straighten the rgbs of all four clusters .
we define an adjusted @xmath73 color given by : @xmath74 where @xmath75 values are given in table 2 , and the coefficients @xmath76 , @xmath77 , and @xmath78 are power series coefficients given in table 3 .
similar equations are derived for the @xmath79 and @xmath80 cmds .    using the adjusted values ,
histograms of colors readily show that the rgbs of ngc 6388 and ngc 6441 ( figure 4 ) are wider than those of 47 tuc and even that of @xmath3 cen , which is well established to have a metal spread from [ fe / h ] @xmath81 -2.2 to -0.7 ( vanture , wallerstein , & suntzeff 2002 ) .
the latter result is somewhat surprising , given the fact that the metal - rich rgb of 47 tuc provides a poor approximation to the metal - poor rgb of @xmath3 cen and produces a clear tilt in the `` adjusted '' @xmath3 cen rgb ( fig .
4 ) . the fig .
4 histograms are divided into five @xmath82 colors bins in order to study the spatial distribution of rgb stars by color .
the five color bins from the lowest value of @xmath82 ( most blue ) to the highest ( most red ) are referred to as bins a through e , respectively .
because globular clusters are dynamically mixed on timescales that are short compared to chemical enrichment , there should be no significant difference in the spatial distribution of the stars in these five bins .
for example , we show below ( section 4.2.1 ) that neither 47 tuc nor @xmath3 cen exhibit any such spatial variation .
on the other hand , if differential extinction were significant , the distributions of stars by color would be uneven across the face of the cluster as a reflection of reddening variations .
this is a test we now apply to ngc 6388 and ngc 6441 .    various methods to divide the sample population into color bins were explored .
the results are not very sensitive to which method is chosen .
the following technique ensures similar treatment of each cluster : bin ranges are assigned to each cluster using a @xmath82 histogram of a subsample of the purest cluster stars , defined by an angular radius where the density of stars decreased to a limiting density given by : @xmath83 ( where the outer density is determined as the density of stars at the edge of the field of view ) . in all cases ,
this resulted in a sample of stars within a radius of between roughly one and three arcminutes , which is outside the core radius but well within the tidal radius of each cluster ( see table 1 ) . using a color histogram with binning width 0.02 mag , the central three bins are defined within the @xmath82 range from the central peak to the colors where counts decrease to @xmath84 of the central peak .
bin c is centered at the peak of the histogram and has 1/3 of the width of b+c+d .
the range of bin a is determined as the color range blue - ward of bin b up to where the counts first drop to one ( to eliminate extreme colored stars from the sample ) .
likewise , bin e is defined as the range red - ward of bin d up to where the counts first drop to one ( fig .
field star contamination ( or contamination by non - rgb stars in the cluster ) , not differential reddening , is assumed to be the reason for stars that lie very far from the fiducial rgb , and the described culling of the data set promises to retain only those stars that make up the bulk of the rgb and that could potentially have been scattered to their current locations in the cmd by differential reddening .
these color bin definitions are applied to the full 300 arcsecond data set .
similar binnings were created for @xmath85 and @xmath86 histograms , although only the @xmath82 and @xmath85 distributions were found to be useful in our analysis ( see below ) .      before proceeding further
, we apply one additional filter designed to preserve reddening effects and reduce noise from astrophysical and photometric scatter .
this filter takes advantage of the fact that the degree of scatter due to reddening in one color should be correlated to that in another color .
a legendre polynomial with an iterative @xmath88 rejection threshhold ( similar to that used earlier in section 4.1.1 ) was fit to the cluster @xmath82 vs. @xmath85 color - color diagram ( fig .
6 ) , and color - color cuts were performed to eliminate stars at greater than @xmath88 from the legendre fit . this independently derived fit for each cluster is very close to the fiducial reddening vector for standard extinction models with @xmath89 , lending further evidence to support the presence of differential reddening across each cluster . figs .
6 and 7 show color - color plots of @xmath82 versus @xmath85 and @xmath86 color , respectively , with the adopted outlier cuts shown .
note that the dynamic range in color of @xmath86 is small compared to the reddening effect sought , and it is difficult to determine which stars appear redder than average as a result of differential reddening instead of photometric errors or astrophysical reasons .
we therefore reject @xmath86 color from our analysis , since it does not have the discriminating power required to pick out the required distribution from the noise of points in figure 7 .
spatial plots of the stars in each of the @xmath82 bins for 47 tuc , @xmath3 cen , ngc 6388 and ngc 6441 are shown in figures 8 - 11 respectively .
as expected for a monometallic population with little reddening , 47 tuc shows no spatial segregation of the stars in the five bins .
although giant stars in @xmath3 cen have been found to exhibit a spatial variation in metallicity ( jurcsik 1998 ) , this variation is on scales of roughly 20 arcminutes ; in contrast , the mixed metallicity population in our smaller field of view shows no spatial segregation by rgb color .
both 47 tuc and @xmath3 cen have few stars in bins a and e relative to their total population .
in contrast , ngc 6388 exhibits a pronounced abundance of redder stars ( bins d and e ) in the northern and northwestern regions of the cluster compared to southern regions ( fig .
10 ) , and the extreme bins a and e each contain significantly more stars relative to the total population of the cluster than these bins contained for 47 tuc and @xmath3 cen .
comparison to the 47 tuc and @xmath3 cen templates for monometallic and multi - metallic clusters ( respectively ) having comparatively little differential reddening over small scales leads one to the conclusion that it is unlikely that the ngc 6388 color spread is due primarily to metallicity effects .
rather , fig .
10 suggests that differential reddening on small scales is causing ( 1 ) a greater number of stars to appear in bins a and e , and ( 2 ) large spatial variations by observed stellar color .    likewise , ngc 6441 ( fig .
11 ) shows an enhancement of redder stars ( bins d and e ) in the southeastern quadrant of the field of view , supporting the view that it too is subject to substantial differential reddening . using isochrone fitting in subfields ,
heitsch & richtler ( 1999 ) similarly find that the southeast region of ngc 6441 has higher @xmath13 than the northwestern region by about 0.15 mag .
layden et al .
( 1999 ) concur with heitsch & richtler to within a few hundredths of a magnitude using rr lyrae stars and cmds in sectors of an annulus around the cluster center .
however , in contrast to the distinct dichotomy in ngc 6388 ( redder stars in the northwest , bluer stars in the southeast ) , the spatial distribution of stars in bins a through c in ngc 6441 appears relatively constant ( with the possible exception of decreased counts in the southeastern quadrant of bin a ) , and the cluster presents the appearance of an inhomogeneously distributed redder population superposed on a centrally concentrated bluer population background .
this effect may be an artifact of the severe contamination of the cluster by field stars since ngc 6441 is located in an extremely crowded field near the galactic center ( as vividly demonstrated by the washed out starcount map shown in figure 12 ) .
indeed , heitsch & richtler ( 1999 ) note the strong contamination of their own cmd of ngc 6441 by field stars that cover the lower part of the cluster rgb . in this case ,
strong contamination by blue field stars may contribute to the apparent width of the cluster red giant branch .
figure 12 shows 2mass starcount density ( panel a ) and mean stellar color ( panel b ) maps for all four of our clusters . in this figure , higher starcounts and redder
mean colors appear whiter ( table 5 provides the numerical values for each greyscale ) .
the field of view of these maps is more than 100 times larger than that used for the analyses in the rest of this section since this method only provides coarse resolution maps which , however , are useful for understanding qualitatively the global character and degree of patchiness of the dust obscuration around each cluster .
for example , assuming that the average intrinsic star color in a given direction is relatively uniform over degree scales , the mean color maps can be interpreted as showing the general character of the reddening in the direction of each cluster . however , these maps have too coarse a greyscale to show color variations of order @xmath90 .
12 shows that the mean stellar color in the direction of ngc 6388 is not only patchier than that in the direction of 47 tuc or @xmath3 cen but that the number of redder stars appears to be greater to the north of the cluster than to the south , in agreement with the trend found above .
distinct regions of higher reddening and lower reddening are demarcated by a boundary that passes close to the center of the cluster .
differential reddening across this apparent dust boundary may produce the observed overabundance of redder rgb stars in the north - northwestern regions of the cluster compared to the south - southeastern regions .
the ngc 6441 field mean color map also shows evidence of patchy extinction , although it is not obvious whether any of these patches could produce the observed reddening patterns near the cluster that are likely the case for ngc 6388 .
however , the presence of large scale patchiness is sufficient to suggest the possibility of patchiness on the angular scales explored in section 4.2.1 .
histograms of star colors within regions of greater reddening should be offset from those of subregions with lesser reddening .
to study the color distributions in subsections of the cluster fields , which will have larger statistical noise , we apply a uniform width - normalized gaussian kernel ( e.g. , silverman 1986 ) to the @xmath82 color distribution of stars in ngc 6388 and ngc 6441 .
figure 13 shows these smoothed distributions for the northwest and southeast quadrants of each cluster .
the half maximum width of the kernel is the minimum value that produces a clearly defined central peak in the distribution for a given cluster .
ngc 6388 required this smoothing width to be only 0.02 mag , while ngc 6441 required a wider gaussian of width 0.06 mag to smooth out secondary peaks adequately .
the central peak of the distribution of the northwestern quadrant of ngc 6388 is offset from that of the southeastern quadrant by roughly @xmath91 0.061 mag , and the peaks of the northwestern and southeastern quadrants of ngc 6441 are offset by roughly @xmath92 0.072 mag . because the width of the kernel used for ngc 6441 is of order the offset between the peaks , the precision of the derived @xmath93 is less for ngc 6441 than for ngc 6388 .
nevertheless , the @xmath93 results for both clusters are in agreement with the qualitative arguments for differential reddening presented earlier in sections 4.2.1 and 4.2.2 , and our value for the differential reddening across ngc 6441 is consistent with the more detailed optical extinction maps produced by layden et al .
( 1999 ) and heitsch & richtler ( 1999 ) to within 0.02 mag in @xmath75 .
if such an offset exists between the color distribution of stars in different regions of each cluster , then the artificial width produced by the superposition of these offset distributions may contribute to the thicker rgb of the entire cluster .
14 simulates the effect of superimposing two offset color distributions by artificially reddening the northern half of the monometallic cluster 47 tucanae by 0.061 magnitudes in @xmath82 ; it is observed that while such a superposition can artificially inflate the width of the color distribution ( fig .
13 , lower panels ) , it is only by about half the amount required to re - create the width of the distributions observed in ngc 6388 and ngc 6441 . the inability of this simulation to entirely explain the rgb widths of these clusters may be due to the presence of a large reddening _ spread _
( as opposed to the simplistic two value model used here ) , or it may be indicative of a contributing effect of differential metallicity .      to test the significance of the differences in @xmath82 color distribution in different quadrants of each cluster , we applied the kolmogorov - smirnov ( ks ) test to compare spatial color distributions numerically . as a calibration of the method , we first compared the color distribution across the clusters 47 tuc and @xmath3 cen .
as tables 6 and 7 show , the @xmath82 color distribution of stars in , for example , the northwestern quadrant of the clusters correlates with the distribution of stars in the southeastern quadrant with significance levels of about 0.5 and 0.25 for 47 tuc and @xmath3 cen respectively .
these significance levels are consistent with a probability that stars within the two regions could have been randomly drawn from one parent distribution .
tables 6 and 7 demonstrate that this relation significance is roughly similar for comparisons between all quadrants for both 47 tuc and @xmath3 cen .
correlations between the southeastern and southwestern , and northeastern and northwestern quadrants of ngc 6388 have significance levels of roughly 0.9 and 0.5 respectively . in contrast , the ks test yields a probability of only 0.002 that the northwestern and southeastern quadrants share a similar color distribution ( table 8) . compare this to the variable metallicity cluster @xmath3 cen , which maintains a relatively high correlation significance of color distribution between all quadrants in spite of its metallicity spread .
it appears unlikely that a metallicity spread is the dominant cause of the rgb width in ngc 6388 ( we expect mixing in clusters in any case ) .
instead , the stars in these regions almost certainly come from statistically distinct color populations , pointing again to the presence of strong differential reddening .    application of the ks test to the @xmath82 color distribution in the southeastern and northwestern quadrants of ngc 6441 yields a near zero likelihood of correlation ( table 9 ) . qualitatively , this is not surprising because there are very few stars from bins d and e in the northwestern quadrant while there are a considerably higher number of redder stars that lie in the southeastern quadrant .
however , while these two quadrants are by orders of magnitude the least correlated regions of this globular cluster , correlations between all other combinations of quadrants are only of order 0.01 to 0.1 , which suggests ( as we have concluded in section 4.2.2 ) that reddening in ngc 6441 is patchy on smaller angular scales than in ngc 6388 .
if differential reddening is producing the observed width in these clusters , the rgb histograms in sub - regions of each cluster should also be narrower than that of each cluster as a whole .
using hst - wfpc2 visual band data for ngc 6388 and ngc 6441 , raimondo et al .
( 2002 ) find this to be the case for both clusters . in summary , it appears very probable that differential reddening is artificially inflating the apparent thickness of the rgb in these two clusters , although the thickness may not be fully explained by this effect alone .
differential extinction complicates photometric studies of galactic stellar populations , and even nearby globular clusters with @xmath13 as low as about 0.1 mag may be visibly affected . in our search for tidal tails around the giant globular cluster @xmath3 cen
we have found that differentially dereddening 2mass photometry largely eliminates apparent tidal tail candidates present in data that have not been corrected for extinction .
given the similarity of the false 2mass tidal tail signature to the optical tails reported by leon et al , it would be beneficial to verify the diffuse optical starcount distribution in the @xmath3 cen field with deep visual photometry that accounts for the @xmath94 variations .
we have also investigated the effects of differential reddening on the cmds of the clusters ngc 6388 and ngc 6441 , which exhibit abnormally large spreads in their rgbs .
the analyses presented make use of spatial maps by stars in color bins , maps of mean stellar color , kernel - smoothed histograms of stellar color in sub - regions , and ks statistical analyses of the distribution of ( @xmath21 ) and ( @xmath20 ) colors across the clusters .
these tests indicate that the northwestern quadrant of ngc 6388 is subject to about @xmath95 0.12 mag more reddening than the southeastern quadrant , and that the southeastern quadrant of ngc 6441 is more heavily reddened by about @xmath95 0.14 mag than the northwestern quadrant .
these findings are consistent with previous results from heitsch & richtler ( 1999 ) and layden et al .
( 1999 ) . as numerical experiments performed here and by raimondo et al .
( 2002 ) have shown , this amount of differential reddening could produce anomalously wider red giant branches for these clusters ( in addition to the tilt observed in their horizontal branches , as demonstrated by raimondo et al . ) .
raimondo et al . have shown that the rgb in sub - regions of these clusters is thin ( as expected if differential reddening were the cause of the width of the rgb across the entire cluster ) .
athough we find that possible width inflation due to the superposition of thin rgbs from two _ quadrants _ subject to different amounts of differential extinction may not be sufficient to completely explain the observed width of the cluster rgbs , it may be possible that the rgb width could be explained more fully by a more realistic superposition of rgbs from smaller subregions .
therefore , if not accounting for the entire rgb - widening effect , at minimum differential reddening may be exaggerating metallicity - induced rgb spreads and the actual spread in age or metalliticy within the clusters is likely not as large as has previously been suggested ( e.g. , pritzl et al .
2001 , 2002 ) .
the number of known globular clusters showing clearly anomalous metallicity spreads remains small . with few bona fide members , this `` class '' of globular cluster
can still comfortably be accomodated by origin pictures that are exceptional , e.g. , along the lines of the m54 - sgr paradigm .
however , even in the obvious case of the multi - population @xmath3 cen , the m54 analogy may be less compelling than previously thought , and the picture that some globular clusters are the `` residue '' of dsph disruption more ambiguous .
the authors would like to thank s. leon and g. meylan for helpful discussion regarding the effects of reddening on @xmath3 centauri .
this publication makes use of data products from the two micron all sky survey , which is a joint project of the university of massachusetts and the infrared processing and analysis center , funded by the national aeronautics and space administration and the national science foundation .
bassino , l.p . and
muzzio , j.c .
1995 , the observatory , 115 , 256 .
cardelli , j.a . ,
clayton , g.c . and mathis , j.s . 1989 , , 345 , 245 . da costa , g.s . and armandroff , t.e .
1995 , , 109 , 2533 .
grillmair , c.j . ,
freeman , k.c . , irwin , m. , & quinn , p.j .
1995 , , 109 , 2553 .
harris , w.e .
1996 , , 112 , 1487 .
heitsch , f. and richtler , t. 1999 , , 347 , 455 .
jurcsik , j. 1998 , , 506l , 113 .
king , i.r .
1966 , , 71 , 64 .
layden , a.c . and sarajedini , a. 2000 , , 119 , 1760 .
layden , a.c . , ritter , l.a .
, welch , d.l . , and webb , t.m.a . 1999 , , 117 , 1313 .
van leeuwen , f. , hughes , j.d .
, & piotto , g. ( eds . ) 2002 , asp conf .
265 , ( san francisco : asp ) .
lee , y.w . ,
rey , s.c . ,
ree , c.h . ,
joo , j.m . , sohn , y.j . , and yoon , s.j .
2002 , asp conf .
265 : omega centauri , a unique window into astrophysics , 305 .
lee , y.w .
, joo , j.m . ,
sohn , y.j . ,
rey , s.c . ,
lee , h.c . , and
walker , a.r .
1999 , nature , 402 , 55l .
leon , s. , meylan , g. , & combes , f. 2000 , , 359 , 907 .
majewski , s.r . ,
patterson , r.j . ,
dinescu , d.i . ,
johnson , w.y . ,
ostheimer , j.c . , kunkel , w.e . , and palma , c. 2000 , proc .
35th liege international astrophysics colloq . , the galactic halo : > from globular clusters to field stars , ed .
a. noels et al .
( belgium : institut dastrophysique et de geophysique ) , 619 .
odenkirchen , m. , grebel , e.k . ,
rockosi , c.m . ,
dehnen , w. , ibata , r. , rix , h.w . ,
stolte , a. , wolf , c. , anderson , j.e .
bahcall , n.a .
, brinkmann , j. , csabai , i. , hennessy , g. , hindsley , r.b . , ivezi , z. , lupton , r.h . ,
munn , j.a . ,
pier , j.r . ,
stoughton , c. , and york , d.g .
2001 , , 548 , 165 .
piotto , g. , sosin , c. , king , i.r . ,
djorgovski , s.g . , rich , r.m . ,
dorman , b. , renzini , a. , phinney , s. , and liebert , j. 1997 , in advances in stellar evolution , ed .
rood & a. renzini ( cambridge : cambridge university press ) , 84 .
pritzl , b.j . ,
smith , h.a . , catelan , m. , & sweigart , a.v .
2001 , , 122 , 2600 .
pritzl , b.j . ,
smith , h.a .
, catelan , m. , & sweigart , a.v .
2002 , , 124 , 949 .
raimondo , g. , castellani , v. , cassisi , s. , brocato , e. and piotto , g. 2002 , , 569 , 975 .
ree , c.h . ,
yoon , s.j . ,
rey , s.c . , and lee , y.w .
2002 , asp conf .
265 : omega centauri , a unique window into astrophysics , 101 .
sarajedini , a. and layden , a.c .
1995 , , 109 , 1086 .
schlegel , d. , finkbeiner , d. , & davis , m. 1998 , , 500 525 .
silverman , b.w .
1986 , density estimation for statistics and data analysis ( london : chapman & hall ) .
sweigart , a.v .
, & catelan , m. 1998 , , 501 , l63 .
trager , s.c . , king , i.r . , & djorgovski , s. 1995 , , 109 , 218 .
tsujimoto , t. and shigeyama , t. 2003 , astro - ph/0302300 vanture , a.d . , wallerstein , g. , & suntzeff , n.b .
2002 , , 569 , 984 .
yoon , s.j . ,
lee , c.h . , and lee , y.w . 2000 , american astronomical society meeting , 196 .
cccccccc 47 tucanae & 0.124 & 13.37 & 13.26 & 0.04 & 0.021 & 0.011 & 0.009@xmath3 centauri & 0.299 & 13.92 & 13.63 & 0.12 & 0.060 & 0.033 & 0.027ngc 6388 & 1.240 & 16.54 & 15.44 & 0.40 & 0.209 & 0.109 & 0.090ngc 6441 & 1.395 & 16.62 & 15.38 & 0.45 & 0.234 & 0.123 & 0.102


in order to gain information about a quantum state of the radiation field one has to measure some observable .
the measurement process unavoidably involves some kind of interaction , thus coupling the mode under examination to one or more other modes of the field .
one has to admit that , in general , the measured observable is not defined on the sole hilbert space of the signal mode @xcite . on the contrary ,
it reflects properties of the global state , possibly entangled , coming from the interaction among the signal mode and the set of the probe modes .
sometimes , it is possible to get rid of these probe modes , so that the statistics of the output is related only to the quantum statistics of the signal mode .
this , as an example , is the case of balanced homodyne detectors @xcite . in that case , in fact , an appropriate rescaling of the output photocurrent allows to completely neglect the local oscillator in the definition of the measured observable , which simply results to be proportional to the signal field quadrature @xmath0
.    however , in more general cases this procedure can not be pursued and the output statistics always remind us the way we are probing the signal under examination .
this is a common feature of three relevant detection schemes in quantum optics , which provide different setup for jointly measuring a couple of photocurrents .
they are the heterodyne @xcite , the eight - port homodyne @xcite and the recently introduced six - port homodyne detectors @xcite .
it is a purpose of this paper to show that , although they involve very different ways of coupling signal to probe modes , they provide the same information on the signal under examination .
actually , we will prove their full equivalence by demonstrating that their output photocurrents have the same operatorial form .
more surprisingly , this remains true also when taking into account the inefficiency of photodetection , even though the latter occurs at very different stages in the three detection schemes .
the paper will be organized as follows . in the next section we briefly
review how to describe inefficient photodetection in terms of beam splitters and ideal detectors . in the three sections [ s : eig],[s : het ] and [ s : tri ]
we will examine the three detection schemes respectively , in order to show that all of them jointly measure the real and the imaginary part of the complex photocurrent @xcite @xmath1 being @xmath2 the signal mode and @xmath3 a probe mode .
the inefficiency of realistic photodetectors will be also taken into account , in order to show it does not affect the equivalence of the considered detection schemes .    in section refs
: pom we will derive the output statistics of considered schemes as members of the more general class of two - photocurrent devices .
section [ s : out ] will close the paper with some concluding remarks .
the final stage of any detection scheme is represented by photodetection , namely counting photon through their conversion to electric pulses .
let us consider a light beam @xmath4 entering a phototube which converts to electric pulses a fraction @xmath5 of the incoming photons . by keeping open the detector window for a time interval @xmath6 ,
the probability @xmath7 of counting @xmath8 photons is given by @xcite @xmath9^m } \over{m ! } } \exp [ -\zeta\hat i ( t)t]\nr\right\ }          \label{pc - gen}\;,\end{aligned}\ ] ] where @xmath10 denotes normal ordering of operator and @xmath11 is the beam intensity @xmath12 @xmath13 denotes the positive ( negative ) frequency part of the field .
for a single - mode field excited in a stationary state eq .
[ pc - gen ] can be rewritten as @xmath14 being @xmath15=1 $ ] the single mode field operator and @xmath16 the global quantum efficiency of the photodetectors . for unit quantum efficiency
this coincides with the actual photon number distribution of the state under examination @xmath17 whereas , in the realistic case of non - unit quantum efficiency eq .
( [ pc - singlemode ] ) becomes a binomial convolution @xcite . in formula @xmath18
let us now consider the scheme in fig .
[ f : eta ] .
the signal mode @xmath2 is impinged in a beam splitter of transmissivity @xmath19 whose second port @xmath3 is placed in the vacuum .
a perfect ( @xmath20 ) photodetection on the exiting mode reveals @xmath8 photons with a probability @xmath21 being @xmath22 the identity operator on the second input of the beam splitter and @xmath23 the evolution operator of the beam splitter .
a straightforward calculation shows that @xmath24 eq .
( [ conv_n2 ] ) coincides with eq .
( [ conv_n ] ) for @xmath25 .
this means that a photodetection process with efficiency @xmath26 is equivalent to a perfect photodetection process performed on the original signal mixed with vacuum by a beam splitter with a value of the transmissivity equal to the quantum efficiency @xcite . in the following
we will adopt this equivalent scheme .
eight - port homodyne detector is known for a long time for the joint determination of phase and amplitude of the field in the microwave domain .
it was introduced in the optical domain by walker and carrol @xcite and successively analyzed by different authors @xcite . a schematic diagram of the experimental setup is reported in fig .
[ f : eig ] . there are four balanced beam splitters whereas a @xmath27 phase shifter is inserted in one arm .
the four input modes are denoted by @xmath28 , @xmath29 whereas the detected output modes are denoted by @xmath30 , @xmath29 .
there are four identical photodetectors whose quantum efficiency is given by @xmath26 .
the _ noise _ modes used to take into account inefficiency , according to the scheme of the previous section , are denoted by @xmath31 , @xmath32 .
we consider @xmath33 as the signal mode , whereas @xmath34 is referred to be the idler of the device .
the mode @xmath35 is unexcited , whereas @xmath36 is placed in a highly excited coherent state @xmath37 provided by an intense laser beam ( local oscillator ) .
the detected photocurrents are @xmath38 , @xmath29 which form the eight - port homodyne observables @xmath39 the latter are derived by rescaling the difference photocurrent , each of them obtained in an homodyne scheme .
for this reason eight - port homodyne is known also as double - homodyne detection . in eq .
( [ 8phot ] ) @xmath26 denotes the quantum efficiency of the photodetectors whereas @xmath40 is the intensity of the local oscillator . in order to obtain @xmath41 and @xmath42 in terms of the input modes we first note that the input - output mode transformation is necessarily linear , as only passive components are involved in the detection scheme @xcite
thus , we can write @xmath43 where the transformation matrix can be computed starting from the corresponding transformations for the beam splitters and the phase shifter @xcite @xmath44 \label{matrix}\:.\end{aligned}\ ] ] eqs .
( [ linear ] ) and ( [ matrix ] ) together with the equivalent scheme for the inefficient detector leads to the following expression for the output modes @xmath45&+\sqrt{1-\eta}\ ; u_1    \nonumber \\
b_2 = & \sqrt{\eta}\;\left[a_1 + ia_2 -a_3-i a_4\right]&+\sqrt{1-\eta}\ ; u_2   \nonumber   \\
b_3 = & \sqrt{\eta}\;\left[a_1- a_2 + ia_3- a_4\right]&+\sqrt{1-\eta}\ ; u_3   \nonumber \\
b_4 = & \sqrt{\eta}\;\left[a_1- ia_2 -a_3 + i a_4\right]&+\sqrt{1-\eta}\ ; u_4 \label{output}\:.\end{aligned}\ ] ] upon inserting eqs .
( [ output ] ) in eq .
( [ 8phot ] ) and considering the limit of highly excited local oscillator we obtain the eight - port photocurrents in terms of the input modes @xmath46 + o[\frac{1}{|z| } ] \nonumber \\
\fl \widehat { \cal z}_2&= & -\hat a_{1 } ( \pi/2 ) + \hat a_{2 } ( \pi/2 ) + \sqrt{\frac{1-\eta}{\eta } } \left [ \hat u_{3 } ( \pi/2 ) - \hat u_{4 } ( \pi/2 ) \right ] + o[\frac{1}{|z| } ] \label{zphot}\:.\end{aligned}\ ] ] in eq . ( [ zphot ] )
@xmath47 denotes a quadrature of the field .
the complex photocurrent @xmath48 is given by @xmath49 for unit quantum efficiency , whereas for non unit quantum efficiency it becomes a gaussian convolution of eq .
( [ 8zed ] ) , we will consider this point in detail in [ s : pom ] .
it is worth noticing here that the mode transformation defined by eqs .
( [ linear ] ) and ( [ matrix ] ) is distinctive for a canonical @xmath50-port linear coupler as defined in refs .
it has been rigorously shown @xcite that a @xmath51-port linear coupler can always be realized in terms of a number of beam splitters and phase - shifters . however
, this implementation is , in general , not unique .
the interest of eight - port homodyne scheme lies in the fact it provides the minimal scheme for realizing a @xmath50-multiport @xcite .
if one considers the multiport as a given _ black - box _ device the eight - port homodyne scheme can be depicted as in fig .
[ f : rough ] .
this will facilitate the comparison with the six - port homodyne detection presented in section [ s : tri ] .
heterodyne detection scheme is known for a long time in radiophysics .
it has been introduced in the domain of optics @xcite in order to describe the joint measurement of two conjugated quadratures of the field .
the term heterodyne is used as the involved field modes are excited on different frequencies .    in fig .
[ f : het ] we show a schematic diagram of the detector .
we denote by @xmath52 the signal field , whereas @xmath53 describes the local oscillator .
the field @xmath54 accounts for the losses due to inefficient photodetection .
the input signal is excited in a single mode ( say @xmath2 ) at the frequency @xmath55 , whereas also the local oscillator is excited in only one mode at the frequency @xmath56 .
this local oscillator mode is placed in a strong coherent state @xmath37 by means of an intense laser beam .
the beam splitter has a transmissivity given by @xmath19 , whereas the photodetectors shows quantum efficiency @xmath26 .
the heterodyne output photocurrents are given by the real @xmath57 and the imaginary @xmath58 part of the complex photocurrent @xmath59 .
the latter is obtained after the rescaling of the output photocurrent @xmath60 , which is measured at the intermediate frequency @xmath61 . by fourier
transform of eq .
( [ intensity ] ) we have @xmath62 being @xmath63 positive and negative part of the output field . in terms of the input fields
( [ hetero1 ] ) can be written as @xmath64 \nonumber \\ & & \left[\sqrt{\eta\tau } \hat{\e}^{(+)}_s ( \omega ' ) + \sqrt{\eta ( 1-\tau)}\hat{\e}^{(+)}_{lo } ( \omega ' ) + \sqrt{1-\eta}\hat{\e}^{(+)}_l ( \omega ' ) \right ]   \;.\label{hetero2}\end{aligned}\ ] ] heterodyne photocurrent is obtained by the following rescaling @xmath65 physically this definition corresponds to require a very intense local oscillator , which however is allowed only for a little mixing with the signal mode @xcite . in this limit
only terms containing the local oscillator field @xmath66 at the frequency @xmath56 can survive in eq .
( [ hetero2 ] ) , so that we have @xmath67 where @xmath68 in writing eq .
( [ hetero5 ] ) we have substituted @xmath69 for the relevant modes involved . provided @xmath70 and @xmath71 to be noise modes placed in the vacuum the expression ( [ hetero5 ] ) for the heterodyne photocurrents leads to the identical output statistics of eight - port homodyne photocurrents , the mode @xmath72 playing the role of the idler of the device . the full equivalence of the two detection schemes has been thus proved .
a linear , symmetric three - port optical coupler is a straightforward generalization of the customary lossless symmetric beam splitter .
the three input modes @xmath73 , @xmath74 are combined to form 3 output modes @xmath75 , @xmath76 . in analogy
to lossless beam splitters , which are described by unitary 2@xmath772 matrices  @xcite , any lossless triple coupler is characterized by a unitary @xmath78 matrix  @xcite . for the symmetric case we have the form    @xmath79    where each matrix element @xmath80 represents the transmission amplitude from the @xmath81-th input port to the @xmath82-th output port , that is @xmath83 .
such devices have already been implemented in single - mode optical fiber technology and commercial triple coupler have been available for some time  @xcite .
any triple coupler can be also implemented by discrete optical components using symmetric beam splitters and phase shifters only @xcite . as it has already mentioned in section [ s : eig ] , this is due to remarkable mathematical fact that that any unitary @xmath84-dimensional matrix can be factorized into a sequence of 2-dimensional transformations plus phase - shift @xcite .
moreover , this decomposition is not , in general , unique . in fig .
[ f : ttt ] a possible implementation of a triple coupler is schematically reported .
experimental realizations of triple couplers has been reported for both cases , the passive elements case and the optical fiber one @xcite .
let us now consider the measurement scheme of fig .
[ f : thd ] .
the three input modes are mixed by a triple coupler and the resulting output modes are subsequently detected by three identical photodetectors .
the measured photocurrents are proportional to @xmath85 , @xmath86 given by @xmath87 after photodetection a fourier transform ( ft ) on the photocurrents is performed @xmath88 this procedure is a straightforward generalization of the customary two - mode balanced homodyning technique . in that case , in fact ,
the sum and the difference of the two output photocurrents are considered , which actually represent the fourier transform in a two - dimensional space . by means of the identity @xmath89 for the periodic ( modulus 3 )
kronecker delta @xmath90 , we obtain our final expressions for the fourier transformed photocurrents @xmath91 @xmath92 gives no relevant information as it is insensitive to the phase of the signal field , whereas @xmath93 and @xmath94 are hermitian conjugates of each other and contain the relevant information in their real and imaginary part . in the following let us assume @xmath33 is the signal mode and @xmath34 is fed by a highly excited coherent state @xmath95 representing the local oscillator . for large @xmath40
the output photocurrents are intense enough to be easily detected .
they can be combined to give the reduced photocurrents @xmath96 \nonumber \\ \widehat { \cal z}_2 & = & \sqrt{3 } \frac{\hat { \cal i}_2 - \hat { \cal i}_3 } { 2 i|z|}=-\hat a_{1 } ( \pi/2 ) + \hat a_{3 } ( \pi/2 ) + o[\frac{1}{|z| } ] \label{thp}\;,\end{aligned}\ ] ] which we refer to as the _ triple homodyne photocurrents_. again the complex photocurrent @xmath97 has the form @xmath98 , being @xmath33 the signal mode and @xmath35 the idler of the device .    when accounting for the non unit quantum efficiency @xmath26 of the photodetectors the output modes are written as @xmath99 so that the reduced photocurrents are now given by @xmath100 + o[\frac{1}{|z| } ] \nonumber \\
\fl \widehat { \cal z}_2= \sqrt{3 } \frac{\hat { \cal i}_2 - \hat { \cal i}_3 } { 2 i\eta |z|}= -\hat a_{1 } ( \pi/2 ) + \hat a_{3 } ( \pi/2 ) + \sqrt{\frac{1-\eta}{\eta } } \left [ \hat u_{1 } ( \pi/2 ) - \hat u_{3 } ( \pi/2 ) \right ] + o[\frac{1}{|z| } ] \label{6zphot}\:.\end{aligned}\ ] ] when , as it is the case , the modes @xmath101 are placed in the vacuum the six - port photocurrents in eq .
( [ 6zphot ] ) leads to the same statistics of the eight - port photocurrents in eq .
( [ zphot ] ) .
indeed , they describe different devices leading to the same amount of information on the signal mode @xmath33 .
some comments are , however , in order . by comparison of fig .
[ f : rough ] and fig .
[ f : thd ] it appears obvious that six - port homodyne is an optimized scheme relative to the eight - port one .
one mode less is needed to reach the same amount of information , thus decreasing the possible sources of noise .
the reason for this lies in the final stage of the two schemes .
the fourier transform of the six - port photocurrents , in fact , is a more effective procedure relative to the double - homodyne final stage of the eight - port scheme .
this is related to the noise suppression mechanism of homodyne detection .
the latter , in fact , is generalized to the multi - mode case by the fourier transform rather than duplication of the original two - mode scheme .
in this section we analyze with some detail the output statistics of an abstract two - photocurrent device .
the latter is characterized by the joint measurement of the real @xmath57 and the imaginary @xmath58 part of the complex photocurrent @xmath102 when equipped with perfect photodetectors . on the other hand , in the realistic case of inefficient photodetection
the photocurrent is given by @xmath103 in eqs .
( [ zed ] ) and ( [ zeta ] ) @xmath2 and @xmath3 denote single mode radiation field which can be excited in any quantum state .
the operators @xmath70 and @xmath71 also denote single mode radiation field , however strictly placed in the vacuum states .
they are used to simulate losses due inefficient photodetection , according to the equivalent detection scheme of section [ s : eta ] . in the following
we will refer to the mode @xmath2 as the signal mode which is under examination , whereas the mode @xmath3 is in a known and fixed state , playing the role of the probe mode of the device .
this is only for sake of convenience .
it is obvious that the roles of the two modes can be interchanged and any argument can be reversed considering the mode @xmath3 as a signal mode . according to this scheme
we denote the outcome probability density distribution by @xmath104 .
the latter describes , in the complex plane , the state of the mode @xmath2 as probed by the mode @xmath3 . indeed , different choices for the probe mode lead to very different features in the probability distribution .
each experimental random outcome @xmath105 from the joint measurement of @xmath57 and @xmath58 can be considered as a point @xmath106 in the complex plane . on the other hand ,
the output photocurrent @xmath59 is expressed as a sum of different contributions coming from different modes .
therefore , it appears intuitively rather obvious that the resulting probability distribution will be given by a convolution . to be more specific
let us start by considering the ideal case of unit quantum efficiency @xmath20 .
we write the probability distribution @xmath104 as the fourier transform @xmath107 of the characteristic function @xmath108\right\}\label{xi}\:,\end{aligned}\ ] ] being @xmath109 the global density matrix describing both modes @xmath2 and @xmath3 .
we consider probe mode to be independent on the signal mode , so that the input mode is factorized as @xmath110 upon substituting eqs .
( [ zed ] ) and ( [ rho ] ) in eq .
( [ xi ] ) we are able to write the characteristic function @xmath111 as a product @xmath112 being @xmath113 \right\}\ ] ] the displacement operator and @xmath114 the single mode characteristics function , the latter entering in the definition of the wigner function of a single mode radiation field @xcite @xmath115 we now insert eq .
( [ pxi ] ) in eq .
( [ ft ] ) . by means of eq .
( [ wdf ] ) and using the convolution theorem we arrive at the final result @xmath116 the symbol @xmath117 denoting convolution .
( [ con ] ) provides a justification for referring to the mode @xmath3 as the probe of the device .
in fact , it acts as a filter on the signal wigner function .
the latter is not a genuine probability distribution , as it can be negative when describing quantum interference .
thus , it can not be directly sampled by experiments .
however , the convolution ( [ con ] ) make it more regular , leading to @xmath104 which is a measurable distribution .    phase space density as in eq .
( [ con ] ) have been already introduced by wodkiewicz @xcite to account for the effect of the measuring apparatus in a joint measurement of conjugated variables .
originally , they have been termed phase - space _ propensities_. more recently , they also have been used in entropic description of quantum mechanical states @xcite .
two - photocurrent devices appear as the natural setup to start from , in order to experimentally access such kind of phase - space distribution .
wigner function , though it contains a complete description of the quantum state , can not be directly measured .
the phase distributions @xmath104 are smoothed version of it , corresponding to the occurrence of additional noise of purely quantum origin .
indeed , two - photocurrent devices provide the generalized joint measurement of position and momentum , thus unavoidably introducing additional noise by first principles @xcite .
the crucial point to notice here is that they offer the remarkable possibility to manipulate this quantum noise . as it emerges from eq .
( [ con ] ) it can be redirected in the desidered region of the complex plane by suitable choice of the probe mode , according to the kind of information of which is of interest .    in order to incorporate the effects of inefficient detection
we first note that eq .
( [ zeta ] ) differs from eq .
( [ zed ] ) by two additional _ additive _ terms .
therefore , further convolutions to the original wigner function are expected .
indeed , the characteristic function of the whole device is now expressed as @xmath118 where characteristic function for the noise modes can be easily evaluated as @xmath119 \qquad j=1,2 \label{noisexi}\:.\end{aligned}\ ] ] upon inserting eqs .
( [ pxieta ] ) and ( [ noisexi ] ) in eq .
( [ ft ] ) we obtain our final result , namely the output probability distribution @xmath120 of a two - photocurrent device equipped with photodetectors of quantum efficiency @xmath26 @xmath121 \label{noisecon}\;,\end{aligned}\ ] ] being @xmath122 the probability distribution obtained for ideal photodetection and @xmath123 conclusions -----------    in this paper we have introduced the class of two - photocurrent devices .
this kind of detectors are characterized by the fact that they jointly measure the real and the imaginary part of the complex photocurrent @xmath124 , where @xmath2 and @xmath3 describe two single mode of the radiation field .
eight - port homodyne , heterodyne and six - port homodyne detectors belong to this general class , thus their fully equivalence in probing radiation field has been proved .
it has also been proved that this equivalence still holds when the inefficiency of the photodetection process is taken into account .
this is an interesting and unexpected result , as photodetection takes place at very different stages in the three detection schemes .
the three schemes analyzed in this paper are equivalent from the point of view of provided information on the measured signal .
nevertheless , they have different physical implementations which have to be compared .
we pointed out the advantages of the six - port homodyne scheme in comparison with the customary eight - port one .
actually , it provides the minimal scheme to access generalized phase space distribution @xcite .    in a two - photocurrent device
a generalized joint measurement of position and momentum is performed on the signal mode .
this results in a smoothing of the signal wigner function to a measurable distribution , which represents the output probability distribution of the measurement .
some additional noise is unavoidably introduced , according to the heisenberg principle for joint measurement .
however , the filtering process has been shown to be a convolution with the wigner function of the probe mode . therefore , it is possible to manipulate and redirect the noise .
a suitable choice of the probe mode enhances different features of the signal phase space distribution , according to the kind of desidered information .
i would thank prof .
harry paul for his kind hospitality in the group nichtklassiche strahlung of max - planck society and prof .
mauro dariano for introducing me to the subject of phase - space measurements .
i would also thank valentina de renzi , ole steuernagel , gabriel drobny , alexei chizhov and alfred wnsche for profitable discussions .
this work has been partially supported by the university of milan by a scholarship for postgraduate studies in foreign countries .
99 busch p , lahti p j 1995 _ riv .
nuovo cim . _ * 18 * 1 yuen h p , chan v w s 1983 _ opt .
* 8 * 177 yuen h p , shapiro j h 1980 _ ieee trans .
inf . theory _ * it-26 * 78 shapiro j h , wagner s s 1984 _ ieee j. quant .
_ * qe-20 * 803 shapiro j h 1985 _ ieee j. quant .
_ * qe-21 * 237 walker n g , carrol j e 1984 _ electron
. lett . _ * 20 * 981 walker n g , carrol j e 1986 _ opt .
_ * 18 * 355 walker n g 1987 _ j. mod .
opt . _ * 34 * 15 lai y , haus h a 1989 _ quantum opt .
_ * 1 * 99 paris m g a , chizhov a , steuernagel o 1996 _ opt .
comm._*134 * 117 zucchetti a , vogel w , welsch d g 1996 _ phys .
_ a*54 * 856 dariano g m , sacchi m f 1995 _ phys .
* r4309 kelley p l , kleiner w h , 1964 _ phys .
rev . _ * 136 * 316 mandel l , wolf e 1995 _ optical coherence and quantum optics _ , ( cambridge university press ) leonhardt u , paul h 1993 _ phys .
_ a*48 * 4598 dariano g m , macchiavello c , paris m g a 1995 _ phys .
_ a*198 * 286 noh j w , fougeres a , mandel l 1991 _ phys . rev
. lett . _ * 67 * 1426 noh j w , fougeres a , mandel l 1992 _ phys .
_ a*45 * 424 noh j w , fougeres a , mandel l 1992 _ phys .
_ a*46 * 2840 bandilla a , ritze h h 1993 _ quantum opt . _ * 5 * 213 freyberger m , schleich w p 1993 _ phys .
_ a*47 * r30 leonhardt u , paul h 1993 _ phys . rev .
_ a*47 * 2460 dariano g m , paris m g a 1993 _ phys .
_ a*48 * r4039 dariano g m , paris m g a 1994 _ phys .
_ a*49 * 3022 freyberger m , heni m , schleich w p 1995 _ quantum opt . _
* 7 * 187 paris m g a 1996 _ phys .
_ a*53 * 2658 paris m g a 1996 _ opt .
_ * 124 * 277 stenholm s 1994 _ j. mod
_ * 41 * 2483 stenholm s 1995 _ appl .
_ b*60 * 243 trm p , stenholm s , jex i 1995 _ phys . rev . _
a*52 * 4853 jex i , stenholm s , zeilinger a 1995 _ opt . comm . _ * 117 * 95 reck m , zeilinger a , bernstein h j , bertani p 1994 _ phys .
lett . _ * 73 * 58 trm p , jex i , stenholm s 1996 _ j. mod . opt . _
* 43 * 245 yuen h p 1982 _ phys .
_ a*91 * 101 paris m g a 1996 _ phys .
_ a*217 * 78 campos r a , saleh b e a , teich m c 1989 _ phys .
_ a*40 * 1371 mattle k , michler m , weinfurther h , zeilinger a , zukowski m 1995 _ appl .
_ b*60 * s111 weihs g , reck m , weinfurther h , zeilinger a 1996 _ opt .
* 21 * 302 sheem s k 1981 _ j. appl .
* 52 * 3865 hariharan p , sen d 1959 _ j. sci .
_ _ 36 _ 70 zernike f 1950 _ j. opt .
* 40 * 326 wigner e p 1932 _ phys .
_ * 40 * 749 cahill k e , glauber r j 1969 _ phys .
_ * 177 * 1857 cahill k e , glauber r j 1969 _ phys .
rev . _ * 177 * 1882 wodkiewicz k 1984 _ phys .
lett . _ * 52 * 1064 wodkiewicz k 1986 _ phys .
_ a*115 * 304 wodkiewicz k 1988 _ phys .
_ a*129 * 1 buzek v , keitel c h , knight p l 1995 phys .
a*51 * 2575 buzek v , keitel c h , knight p l 1995 phys .
rev . a*51 * 2594 she c y , heffner h 1966 _ phys .
rev . _ * 152 * 1103 arthurs e , goodman m s 1988 _ phys
_ * 60 * 2447 paris m g a , steuernagel o 1996 , in preparation .
( 10,10)(-5,-5)(-2,-4.5)(0,1)0.6(-2,-3.5)(0,1)7.5 ( -4,2)(1,0)2(-4.5,-2)(1,0)0.6 ( -1.2,1.2)(-1,1)1.65(1.2,-1.2)(1,-1)1.65 ( -1.2,-1.2)(-1,-1)1.65(1.2,1.2)(1,1)1.65 ( -3.5,-2)(1,0)7.5(-4,-4.5)(1,0)2 ( -4.5,-4)(0,1)2(2,-4)(0,1)3.5 ( 2,.5)(0,1)3(2,3.9)(0,1)0.6 ( -4,2)(1,0)7.5(3.9,2)(1,0)0.6 ( 2,4.5)(1,0)2(4.5,2)(0,1)2 ( 2,3.5)(1,0).4(2,3.5)(-1,0).4 ( -2,-3.5)(1,0).4(-2,-3.5)(-1,0).4 ( -3.5,-2)(0,1).4(-3.5,-2)(0,-1).4 ( 3.5,2)(0,1).4(3.5,2)(0,-1).4 ( -4.5,-4.5)(2,0)(4.5,4.5 ) ( 2,3.5)(-2,-3.5 ) ( -3.5,-2)(3.5,2 ) ( -2,-1.3)(0,-1)0(-2,-3)(0,-1)0 ( 1,-2)(-1,0)0(2.7,-2)(-1,0)0 ( -1.3,-2)(-1,0)0(-3,-2)(-1,0)0 ( -2,1)(0,-1)0(-2,2.7)(0,-1)0 ( 2,1.3)(0,1)0(2,3)(0,1)0 ( 1.3,2)(1,0)0(3,2)(1,0)0 ( -1,2)(1,0)0(-2.7,2)(1,0)0 ( 2,-1)(0,1)0(2,-2.7)(0,1)0 ( 4.25,4.35 ) ( -4.75,-4.65 ) ( 1.7,-.1 ) ( 4.2,-2.15)@xmath35(-4.6,2)@xmath33 ( -1.5,-3.9)@xmath125(-4.1,-1.3)@xmath126 ( 2.5,3.7)@xmath127(3.9,1.2)@xmath128 ( 1.7,-4.5)@xmath36(-2.35,4.35)@xmath34        ( 15.00,20.00)(25.00,20.00)(15.00,30.00)(25.00,30.00 ) ( 15.00,40.00)(25.00,40.00)(15.00,50.00)(25.00,50.00 ) ( 14.00,50.00)(14.00,40.00 ) ( 14.00,30.00)(14.00,20.00 ) ( 52.00,50.00)(65.00,50.00)(65.00,40.00)(52.00,40.00 ) ( 52.00,30.00)(65.00,30.00)(65.00,20.00)(52.00,20.00 ) ( 65.00,47.00)(65.00,53.00)(65.00,43.00)(65.00,37.00 ) ( 65.00,33.00)(65.00,27.00)(65.00,23.00)(65.00,17.00 ) ( 65,20)(65,30 ) ( 65,40)(65,50 ) ( 68.00,50.00)(75.00,50.00)(75.00,40.00)(68.00,40.00 ) ( 68.00,30.00)(75.00,30.00)(75.00,20.00)(68.00,20.00 ) ( 89.00,45.00)(100.00,45.00)(81.00,23.00)(75.00,20.00 ) ( 85.00,45.00)(85.00,25.00 ) ( 100.00,45.00)(100.00,25.00 ) ( 100.00,25.00)(89.00,25.00)(75.00,50.00)(81.00,47.00 ) ( 75.00,40.00)(81.00,43.00)(75.00,30.00)(81.00,27.00 ) ( 9.,53.)(0,0)[cc]@xmath33(9.,43.)(0,0)[cc]@xmath34 ( 9.,33.)(0,0)[cc]@xmath35(9.,23.)(0,0)[cc]@xmath36 ( 57,53.)(0,0)[cc]@xmath129(57,43.)(0,0)[cc]@xmath130 ( 57,33.)(0,0)[cc]@xmath131(57,23.)(0,0)[cc]@xmath132 ( 72,53.)(0,0)[cc]@xmath128 ( 72,43.)(0,0)[cc]@xmath127 ( 72,33.)(0,0)[cc]@xmath126 ( 72,23.)(0,0)[cc]@xmath125 ( 105.00,45.00)(0,0)[lc]@xmath133 ( 105.00,25.00)(0,0)[lc]@xmath134 ( 85,45)(0,0)[cc]**(85,25)(0,0)[cc]* *    ( 14,7)(-4.5,-3.5)(-2,0)(1,0)6.5(-2,-0.2)(1,0 ) ( -.8,-.2)(1,0)0(0,-2)(0,1)4 ( -.3,-2)(0,1)(.3,-2)(0,1 ) ( -.3,-.8)(0,1)0(.3,-.8)(0,1)0 ( -1,-1)(1,1)2(-.8,0)(1,0)0 ( 0,-.8)(0,1)0(1,0)(1,0)0 ( 0,1)(0,1)0(4.5,0 ) ( 4.5,-.5)(0,1)1(5,0)(1,0).5 ( 5.5,0)(1,1).8(5.5,0)(1,-1).8 ( 6.5,1)(6.5,-1 ) ( 6.7,1)(1,0)1.8(6.7,-1)(1,0)1.8 ( 8.5,1)(1,0)0(8.5,-1)(1,0)0 ( -3.,.4)@xmath135 ( -4.49,-.7)@xmath136 ( -.75,-3)@xmath137 ( 3,-1)(0,1)(3,0)(0,1)0 ( 3,-1.3)@xmath138 ( 2.8,0)(0,1)(2.8,0)(0,-1)0 ( 2,1.3)@xmath139 ( 6,1.5)@xmath140(6,-2)@xmath141 ( 8.65,1)@xmath142(8.65,-1)@xmath143    ( 143.33,53.00)(25.00,25.00 ) ( 50.33,50.00)(55.33,50.00)(50.33,40.00)(55.33,40.00 ) ( 50.33,30.00)(55.33,30.00)(55.33,30.00)(65.33,40.00 ) ( 55.33,40.00)(65.33,30.00)(65.33,40.00)(70.33,40.00 ) ( 65.33,30.00)(70.33,30.00)(55.33,50.00)(70.33,50.00 ) ( 78.33,30.00)(80.33,30.00)(80.33,30.00)(80.33,30.00 ) ( 80.33,30.00)(85.33,30.00)(80.33,40.00)(85.33,40.00 ) ( 80.33,50.00)(85.33,50.00)(80.33,50.00)(70.33,40.00 ) ( 70.33,50.00)(80.33,40.00)(70.33,30.00)(72.33,30.00 ) ( 85.33,50.00)(87.33,50.00)(85.33,40.00)(95.33,40.00 ) ( 85.33,30.00)(95.33,30.00)(93.33,50.00)(95.33,50.00 ) ( 95.33,50.00)(100.33,50.00)(95.33,40.00)(100.33,40.00 ) ( 95.33,30.00)(100.33,30.00)(100.33,40.00)(110.33,50.00 ) ( 100.33,50.00)(110.33,40.00)(110.33,50.00)(120.33,50.00 ) ( 110.33,40.00)(120.33,40.00)(100.33,30.00)(120.33,30.00 ) ( 120.33,30.00)(130.33,40.00)(120.33,40.00)(130.33,30.00 ) ( 130.33,30.00)(140.33,30.00)(130.33,40.00)(140.33,40.00 ) ( 120.33,50.00)(140.33,50.00)(90.33,50.00 ) ( 75.33,30.00)(50.33,50.00 ) ( 50.33,40.00)(50.33,30.00 ) ( 140.33,50.00)(140.33,40.00 ) ( 140.33,30.00 ) ( 57.33,34.00 )                      ( 75.66,30.00)(0,0)[cc]@xmath144 ( 90.66,50.00)(0,0)[cc]@xmath145 ( 145,50.00)(0,0)[cc]@xmath146 ( 145,40.00)(0,0)[cc]@xmath147 ( 145,30.00)(0,0)[cc]@xmath148 ( 46.33,50.00)(0,0)[cc]@xmath149 ( 46.33,40.00)(0,0)[cc]@xmath150 ( 46.33,30.00)(0,0)[cc]@xmath151 ( 67.00,36.00)(0,0)[cc]bs ( 82.00,46.00)(0,0)[cc]bs ( 112.00,46.00)(0,0)[cc]bs ( 132.00,36.00)(0,0)[cc]bs        ( 9.34,36.33)(1,0)10.00(64.67,36.00)(5,6)9.00 ( 73.67,46.66)(1,0)10.00(83.67,46.66)(5,-6)8.67 ( 92.34,36.33)(-4,-5)8.33(84.00,26.00)(-1,0)10.33 ( 73.67,26.00)(-5,6)8.33(83.67,46.66)(4,3)5.67 ( 89.34,51.00)(1,0)15.33(92.00,36.00)(1,0)12.33 ( 84.00,26.00)(5,-4)5.33(89.34,21.66)(1,0)15.33 ( 73.34,26.00)(-1,-1)4.33(9.34,51.33 ) ( 9.34,36.33)(9.34,21.99 ) ( 9.34,36.33)(9.34,51.33 ) ( 105.00,51.00)(105.00,36.00 ) ( 105.00,21.66)(73.34,46.66 ) ( 84.00,46.66)(92.34,36.00 ) ( 84.34,26.33)(65.34,36.33 ) ( 73.34,26.00 ) ( 78.00,36.00)(0,0)[cc]ft ( 6.67,53.83)(0,0)[cc]@xmath149 ( 6.67,38.83)(0,0)[cc]@xmath150 ( 6.67,24.44)(0,0)[cc]@xmath151 ( 37.67,53.83)(0,0)[cc]@xmath146 ( 37.67,38.83)(0,0)[cc]@xmath147 ( 37.67,23.83)(0,0)[cc]@xmath148 ( 57.67,54.33)(0,0)[cc]@xmath152 ( 57.67,39.33)(0,0)[cc]@xmath153 ( 57.67,24.99)(0,0)[cc]@xmath154 ( 107.67,54.33)(0,0)[cc]@xmath155 ( 107.67,39.33)(0,0)[cc]@xmath156 ( 107.67,24.99)(0,0)[cc]@xmath157 ( 68.67,51.00)(1,-1)4.67(9.34,51.33)(1,0)10.33 ( 9.34,22.00)(1,0)10.00(44.67,36.33)(-1,0)12.00 ( 41.34,36.33)(1,0)6.30 ( 47.67,36.33 ) ( 47.67,22.00 ) ( 47.67,51.08 ) ( 47.67,47.58)(0,1)7.00(47.67,32.83)(0,1)7.00 ( 47.67,18.50)(0,1)7.00(51.34,36.33)(1,0)14.00 ( 68.67,51.00)(-1,0)17.67(69.00,21.66)(-1,0)17.89 ( 51.67,21.66)(0,1)0.00(47.67,51.00)(-1,0)13.67 ( 47.67,21.66)(-1,0)13.00(45.67,51.00 ) ( 45.67,21.66)(45.67,36.33 ) ( 53.67,51.00)(53.67,36.33 ) ( 53.67,21.66 )