arXiv:astro-ph/0508217v1 9 Aug 2005 edee ta.20;Aurie l 05,sgetv fa of suggestive 2005), al. et ICL. Aguerri mixed poorly 2004; intracl al. 200 of et al. et distribution Feldmeier also (Durrell the variations field-to-field in marked structure: show ICL PNe irregular Bland- the an (Weil, of Studies at M87 hint 1997). around Malin envelope & well stellar as Hawthorn, extended , 1998) Virgo an many al. around as et photo features Katsiyannis tidal Deep 1994; show also Malin interactions. (e.g., tidal imaging ongoing graphic in- kine strong, 1999), are Kenney of galaxies & cluster dicative Waterman, the (Rubin, of disturbed Many matically 1987). galaxi al. of et mixes (Bingelli morphological sub- different different (Bingelli, possessing with groups kinematic substructure, (Bin- and 1993) Tammann spatial 1987) & Tammann both Popescu, signs the & possesses many for cluster Sandage, shows ideal The gelli, it environment, ICL. as complex of interest dynamically production particular a of being is of Mpc) 16 of galax a of history assembly cluster. and state dynamical ab information detailed important the holds thus phase-sp ICL The the stripp of clusters. distribution tidal of to assembly linked hierarchical 2004; intimately the al. is during ICL et the Murante of that formation showing 2004; and the al. 2005), Portinari et distri & Willman Romeo, kinematic Sommer-Larsen, (e.g., and ICL spatial hig of yield the tion to of N-body starting predictions now front, are theoretical resolution clusters the al. On et of (Arnaboldi 2005). open simulations kinematics al. latter ICL et the studying Gerhard field 2004; therein), new references a and up ing 2005 a al. 2004 al. et et intra Aguerri Feldmeier and see 2002) (IPNe; al. nebulae et planetary Durrell cluster (Ferg 1998; stars Hippel individual von in 1995; & imaging, Tanvir, both al. son, detected broadband been deep et also using has Bernstein ICL 2004) the 1991; 2002, quantified al. Kuhn and et & Feldmeier review) Boughn, for Uson, 1999 (e.g., i Vílchez-Gómez First history. (see dynamical tified their studying for tool valuable lcrncades io@aeeu [email protected] [email protected], address: Electronic 1 nti otx,teVrocutr(ta dpe distance adopted an (at cluster Virgo the context, this In a is clusters galaxy in (ICL) light intracluster diffuse The rpittpstuigL using typeset Preprint D S srnm n srpyisFlo,NA Tucson NOAO Fellow, Astrophysics and Astronomy NSF ATVERSION RAFT ih IL.Oriaerahsa1 a reaches image Our (ICL). light ujc headings: being Subject su than complex rather assembly, The cluster galaxy. subgroups. of i physical nature galaxies true Several irregular hierarchical very are pair. appearing M84/M86 they kpc, the suggesting a around 200 tails detected nearly tidal also to smaller-scale is intraclust out of diffuse traced myriad is of a M87 web as intricate well as an streamers, reveals and — surveys epeetde pia mgn fteinner the of imaging optical deep present We N A OVEMBER T E tl mltajv 6/22/04 v. emulateapj style X .C J. 1. eateto srnm,Cs etr eev University, Reserve Western Case Astronomy, of Department INTRODUCTION HRISTOPHER ,2018 7, aais lses niiul(ig)—glxe:intera galaxies: — (Virgo) individual clusters: galaxies: M IHOS IFS IH NTEVROCLUSTER VIRGO THE IN LIGHT DIFFUSE P , σ et of depth AUL s.d,jhfna.d,[email protected] [email protected], ase.edu, H rf eso oebr7 2018 7, November version Draft ARDING µ ∼ uster den- v ABSTRACT ace bu- out ing 1.5 2. mag/arcsec =28.5 nd es h- u- 2; y J , - - - - ◦ x1.5 OHN ◦ ih1.45 with ayiae hc utn nyaprino h lse it- image cluster individual the the from of galaxies portion and stars a remove To only self. subtend which images many magnitudes. in V taken Johnson exposures into 900s transformed all and are M Washington skies T blank position. and and point images time sky object in offset exposures 127 object were the from the built images bracketing covering flat ings Individual taken, sky night were Virgo. a using of images flattened degrees object square 72 2.25 of inner total A west. n t Ie24x06CDi prxmtl 1.5 approximately Schmidt the is of CCD view 2048x4096 of SITe field its a The and observing the prepara strategies). of in details reduction for 2005 data 2002 al. al. et et Feldmeier Harding and see tion, recently refurbishment the the technical using of (for tails 2004 telescope Schmidt April Burrell and 0.6m refurbished March in nights tometric ICL. reveale diffuse have of which observations web of intricate Virgo’s season first w the Here on photometry. opt report surface is wide-field Schmidt deep, Burrell for the suited mally such, compare As problems telescopes. b light open-tube optical to scattered extensive reduce and significantly design closed-tube fling fielding Its flat accurate highly possible. Burrell chip, is view CCD Case’s single of a using field onto wide imaged core Schmidt’s is Virgo’s the Because of telescope. study Schmidt imaging deep a measurements intraclus anci hot and resolution galaxies spatial the high of medium. of properties the wealth yet on a determined, data is lary poorly there still that is here Virgo the is in Accordingly sample ICL cluster. uniformly Virgo of to the distribution enough of deep depth line-of-sight gone survey full need not IPNe the have while by scales, date degree hindered to over is field angular flat also large accurately photometry to of surface clusters Deep nearby survey to size. sys- hard telescope/CCD it most makes of difficu view tems of quite field proved small actually relatively has The Virgo like clusters nearby F fteVrocutrt erhfrdfueintracluster diffuse for search to cluster Virgo the of h rdc fsot crto rudacentral a around accretion smooth of product the ELDMEIER dbigsbtenglxe.Tedfuehl of halo diffuse The galaxies. between bridges nd k utato sacalnewe oacigtogether mosaicking when challenge a is subtraction Sky epiae fteVrocutrwr ae uigpho- during taken were cluster Virgo the of images Deep make and disparity this address To i ICL the of structure large-scale the quantifying However, rlgt esesvrlln ( long several see We light. er h oeaeebde ncmo envelopes, common in embedded are core the n nteesae,wiesgicn ifs light diffuse significant while scales, these on 00 uldAe lvln,O 44106 OH Cleveland, Ave, Euclid 10900 srcueo ig’ ifs C eet the reflects ICL diffuse Virgo’s of bstructure 2 ′′ . mag/arcsec 1.5 — iesadteln xso h hpoine east- oriented chip the of axis long the and pixels 1 , AND ftedfuelgti ig,w aebegun have we Virgo, in light diffuse the of ctions H 2. EATHER MGN TECHNIQUE IMAGING M 2 ORRISON eprta previous than deeper > uniaie large-scale quantitative, 0 p)tidal kpc) 100 ◦ x0.75 the , de- af- ter nd he lt. ◦ s, i- l- it n d d e s - - , 2

FIG. 1.— Diffuse light in the core. North is up; east is to the left. The white levels saturate at µv∼26.5, while the faintest features visible have a surface brightness of µv∼28.5. we first mask stars out to where the PSF falls below 3 ADU 5 ADU, or ∼ 0.35% of the DC sky level.2 (µv=28), then subtract their extended wings out to 0.5 ADU After sky subtraction of the individual images, the images (µv=30). IRAF’s OBJMASK task is then used to mask out all are registered and medianed together to form the final mo- extended sources 2.5σ above sky, which essentially removes saic. At this point one last plane is also fit and subtracted all bright galaxies. We then remove sky in a multi-step pro- from the mosaic using the pure sky regions; the edge-to- cess. We first rebin the masked images into smaller images edge gradient of this final plane is 2 ADU (∼ 0.1% of sky). by calculating the mode in 32x32 pixel bins, and then identify We then improve the signal-to-noise by medianing the mo- regions on the images which are located far away from any saic in a 9x9 pixel boxes (∼ 1 kpc2), which brings out the bright galaxy and considered “pure sky”; on average, about faintest features. In the raw image, the 1σ level in the sky 5% of each image is covered by pure sky. Then for each im- is 4.5 ADU (µv=27.5), while in the binned, medianed image age, we calculate a DC sky level from the mode of the counts 1σ=1.7 ADU (µv=28.6). This limiting surface brightness is in these regions and subtract it from the image. This sky approximately 1.5 mag/arcsec2 deeper than the deep photo- level ranges from 1100–1500 ADU (µv=21.6–21.3), depend- graphic imaging of Malin (1994). ing on a variety of factors such as hour angle, airmass, time of As a final note, one possible source of confusion in our night, and air glow variability. Next, we perform an iterative study is contamination due to the presence of any galactic cir- plane-fitting process where the pure sky regions are used to si- rus in the field, which will be visible as large-scale diffuse multaneously fit individual residual sky planes to each binned image, constrained to minimize frame-to-frame deviations in 2 It is important to recognize that fitting and subtracting sky planes can impose a systematic underestimate of the Virgo ICL, since any degree-scale the pure-sky regions of the final sky-subtracted images. The diffuse ICL component will also be subtracted at this step. edge-to-edge gradient of the final fitted sky planes is typically 3

has total length of 38′ (178 kpc) a characteristic width of 3.′5 9 (16 kpc), and a total luminosity of ∼ 10 L⊙. Because NGC 4458 and NGC 4461 have a very high velocity relative to one another (∆v = 1296 km/s) it is unlikely the streamer comes from any strong interaction between the pair; instead, strip- ping from one of the galaxies individually is more likely. The second streamer (B) is perhaps even more remarkable. Also appearingto emanate from M87 to the northwest, it mea- sures approximately 28′ (130 kpc) long and only 1.2′ (6 kpc) wide, with a peak surface brightness of µv=27.0. The low surface brightness dwarf VCC 1149 is projected along the stream, just before the position at which the stream blends into M87’s halo; we may well be seeing the destruction of this dwarf in the tidal field of the cluster — a process which may contribute to the variation in the faint end slope of the lu- FIG. 2.— Diffuse light in the Virgo cluster, color coded by surface bright- minosity function in clusters. The thinness of the stream sup- ness: Blue: µV =25-26; Yellow: µV =26-27; Green: µV =27-28; Red: µV =28- 29. North is up; east is to the left. ports this conjecture, as the low velocity dispersion of dwarf galaxies ensures that stripped tidal streams will be dynami- cally cold, and not diffuse significantly in the tangential di- rection. The thinness and linearity of the streamer also argue that the dwarf must be on a highly radial orbit, or one that is seen projected along its orbital plane. Not all detected streams are linear, however. A broad, curved stream (C) can be seen surrounding a low surface brightness dwarf; the curvature of this stream, arcing away from M87, suggests that the galaxy is orbiting within or around the M86 group (or perhaps the NGC4435/8 pair) rather than around M87 itself. Aside from long ICL streams, many galaxies show tidal features on smaller scales. Perhaps the most obvious is the tidal tail to the north of the interacting pair NGC 4435/8 (D). Originally identified by Malin (1994) in photographic imag- ing, our imaging traces it to fainter surface brightnesses (µv> ◦ FIG. 3.— Schematic diagram showing location of diffuse features discussed 28), whereit takes a sudden 90 bendto the west. Thiskind of in the text. tidal “dogleg” is expected during close and slow encounters in a , where the tidal forces of both the individual galaxies and the cluster itself conspire to create rather com- reflected light. Examining the IRAS 100µm maps of Virgo, it plex tidal morphologies (Mihos 2005, in preparation). Other is clear that a ring of cirrus (∼ 1.5◦ in diameter) surrounds the notable features include NGC 4388’s tidal plumes to the west Virgo core and may contaminate fields to the southwest and and northeast (E), a southern tail in NGC 4425 (F; both noted north of our field. Fortuitously, however, the Virgo core itself by Malin 1994) which connects to VCC 987, another tail ex- is quite clear of galactic cirrus, arguing that features we detect tending from IC 3349 (G) and passing through VCC 942, and are bona fide structures in Virgo’s diffuse light. a loop to the east of IC 3475 (H). At low surface brightness, the giant elliptical M87 shows 3. DIFFUSE LIGHT IN VIRGO’S CORE a very extended and irregular stellar envelope. From pho- A greyscale representation of the final binned mosaic is tographic images, Weil et al. (1997) and Katsiyannis et al. shown in Figure 1, while Figure 2 shows the image colored (1998) show M87’s luminous halo extending to ∼ 100 kpc according to surface brightness. The faintest features visible scales; on our image it can be traced further out to 175 kpc are at µv∼ 28.5 or about 0.1% of sky. A wealth of diffuse (I). Here the outer isophotes are not simply the radial exten- features can be seen in Virgo’s ICL, ranging from extended sion of M87’s inner regular luminosity profile; traced to the low surface brightness envelopes to long, thin streamers, as northwest, the isophotes become boxier and show irregular well as smaller-scale tidal features associated with many of loops and shells which may be remnants from earlier strip- the Virgo galaxies. Our goal here is to focus on the qualitative ping from, or cannibalization of, smaller galaxies. The lop- morphology of the ICL; a more detailed quantitative analysis sidedness and irregularity of M87’s halo is again indicative of will be addressed in later papers. Figure 3 shows a finding a diffuse envelope which is not fully relaxed. chart for the features described below. Because of the close proximity of M86, M84, and a number In the image, two long streamers can be seen extending to of smaller galaxies, uniquely tracing out any individual halo the northwest from M87 at the lower left. One streamer (A) in that complex of galaxies proves difficult. The isophotes projects through the pair of galaxies NGC 4458/61 and ex- of both M86 and M84 appear fairly regular out to a surface tends beyond, toward a group of galaxies to the north. There brightness of µv∼ 26.5 (J), as noted by Malin (1994). At is a surface brightness gradient along its length; near M87 it fainter surface brightnesses, the isophotes blend together and has a surface brightness of µv=27, dropping to µv=27.8 near become quite irregular, and we see a “crown” of diffuse light the pair. From the point at which it emerges from M87’s stel- (at µv= 27–28)to the northofthe M84/M86system (K), anda lar halo to where it fades into the background, the streamer diffuse extension to the south of the pair, blending into the net- 4 work of galaxies to the south. Whether this extended envelope and Aguerri et al. (2005). Because IPNe detection depends is truly physical or simply a projection effect remains unclear. strongly on the underlying surface brightness, IPNe studies Also visible around M86 are several low surface brightness have complicated selection and completeness functions, par- filaments (L1–L3), one to the northeast (µv∼ 25.5), one to ticularly in fields which cover both bright galaxies and diffuse the west (µv∼ 26.5), and one to the north (µv∼ 27.5). intracluster space. With this in mind, we will restrict the dis- We also see a common envelope of light (M) surround- cussion here to a qualitative comparison, and leave a more ing the complex of galaxies to the south of M86, including quantitative analysis to a future paper. NGC 4413, IC 3363, and IC 3349, along with a tidal bridge Both IPNe studies show that the IPNe density is highly vari- connecting NGC 4413 and IC 3363, and perhaps extending able in the Virgo core; our image confirms this complexity to NGC 4388 to the east. This envelope has fairly irregular and variation in the underlying ICL. The highest IPNe densi- isophotes, arguing it is more than just simply the extended ties are found in Feldmeier et al.’s “Field 3,” which our image profiles of galaxies overlapping in projection. This group of shows is completely embedded in M87’s diffuse halo. Much galaxies may be a physically coherent subgroup within the of the IPNe detected in this field are likely to be associated Virgo cluster. with M87 itself; radial velocity measurements of the IPNe in Finally, Oosterloo & van Gorkom (2005) have recently this field agree with this assessment (Arnaboldi et al. 2005). identified an extended (110x25 kpc) HI cloud streaming from Aguerri et al.’s “SUB” field also shows a relatively high PNe NGC 4388 which they attribute to ram-pressure stripping due density; aside from containing the luminous ellipticals M84 to the lack of any stellar counterpart in previous photographic and M86, this field also partly covers the diffuse envelope imaging. We also find no extended counterpart on our im- surrounding the group of galaxies to the south of M86. The age down to a surface brightness limit of µv=28.5, save for lowest IPNe densities are reported in Aguerri et al.’s “LPC” faint streamers close to the main body of NGC 4388 coinci- field, which also corresponds to the lowest luminosity density dent with the extended emission line region noted by Yoshida in our image; indeed, over most of this field, we detect no et al. (2002). light above the background. The qualitative match between the luminosity density 4. DISCUSSION traced by IPNe and broadbandlight argues that IPne are effec- Our deep imaging of the Virgo cluster has revealed a com- tive tracers of the ICL in galaxy clusters, and vice-versa. Our plex network of extended tidal features, indicative of the on- deep imaging provides a “finding chart” for IPNe searches; going stripping and tidal evolution of galaxies in Virgo. The followup spectroscopy can then be used to study the kine- ICL is not radially symmetric around M87, the nominal cen- matics of discrete structures such as tidal tails and diffuse tral galaxy of Virgo; much of the diffuse light is centered upon envelopes around galaxies. Our survey of diffuse light in the M84/M86 complex. The impact of hierarchical assembly Virgo continues: we have follow-up imaging of additional of clusters on the ICL is clear here — rather than ICL growing portions of the Virgo cluster, and subsequent papers will focus simply via smooth accretion around a central galaxy, it distri- on quantitative descriptions of Virgo’s ICL, as well as a sur- bution reflects the substructure inherent in the cluster, as ev- vey of extremely low surface brightness galaxies in the cluster idenced in cosmologically-motivated N-body simulations of core. ICL (e.g., Dubinski 1998; Willman et al. 2004, Rudick et al. in preparation). Characterizing the ICL in terms of a simple radial profile will prove misleading in all but the most regular, relaxed clusters. This work is supported by NSF grants AST-9876143 It is also instructive to compare our map of the ICL in (JCM), AST-0098435 (HLM), and AST-0302030 (JF), and Virgo’s core to estimates of the luminosity density derived by Research Corporation Cottrell Scholarships to JCM and from studies of planetary nebulae by Feldmeier et al. (2004) HLM.

REFERENCES Aguerri, J.A.L., Gerhard, O.E., Arnaboldi, M., Napolitano, N.R., Castro- Katsiyannis, A. C., Kemp, S. N., Berry, D. S., & Meaburn, J. 1998, A&AS, Rodriguez, N., & Freeman, K.C. 2005, AJ, 129, 2585 132, 387 Arnaboldi, M. et al. 2004, ApJ, 614, L33 Malin, D. 1994, IAU Symp. 161: Astronomy from Wide-Field Imaging, 161, Bernstein, G. M., Nichol, R. C., Tyson, J. A., Ulmer, M. P., & Wittman, D. 567 1995, AJ, 110, 1507 Murante, G., et al. 2004, ApJ, 607, L83 Binggeli, B.,Sandage, A., & Tammann, G.A. 1987, AJ, 94, 251 Oosterloo, T., & van Gorkom, J. 2005, A&A, 437, 190 Binggeli, B., Popescu, C. C., & Tammann, G. A. 1993, A&AS, 98, 275 Rubin, V. C., Waterman, A. H., & Kenney, J. D. P. 1999, AJ, 118, 236 Dubinski, J. 1998, ApJ, 502, 141 Sommer-Larsen, J., Romeo, A. D., & Portinari, L. 2005, MNRAS, 357, 478 Durrell, P., Ciardullo, R., Feldmeier, J. J., Jacoby, G. H., & Sigurdsson, S. Uson, J. M., Boughn, S. P., & Kuhn, J. R. 1991, ApJ, 369, 46 2002, ApJ, 570, 119 Vílchez-Gómez, R. 1999, ASP Conf. Ser. 170: The Low Surface Brightness Feldmeier, J. J., Ciardullo, R., Jacoby, G. H., & Durrell, P.R. 2004, ApJ, 615, Universe, 349 196 Weil, M.L., Bland-Hawthorn, J., & Malin, D.F. 1997, ApJ, 490, 664 Feldmeier, J.J., Mihos, J.C., Morrison, H.L., Rodney, S.A., & Harding, P.H. Willman, B., Governato, F., Wadsley, J., & Quinn, T. 2004, MNRAS, 355, 2002, ApJ, 575, 779 159 Feldmeier, J.J., Mihos, J.C., Morrison, H.L., Harding, P.H., Kaib, N., & Yoshida, M. et al. 2002, ApJ, 567, 118 Dubinski, J. 2004, ApJ, 609, 617 Ferguson, H. C., Tanvir, N. R., & von Hippel, T. 1998, Nature, 391, 461 Gerhard, O., Arnaboldi, M., Freeman, K.C., Kashikawa, N., Okamura, S., & Yasuda, N. 2005, ApJ, 621, L93