arXiv:1111.7042v2 [astro-ph.GA] 28 Feb 2012 rp rudteGlx tlatoc,ie edn and over leading covering seen dwarf i.e. each once, the be is least can from it at tails Currently, torn trailing debris the 2003). around stellar al. wraps New- et the of Majewski (e.g. that amount apparent 2002; established the really al. by that became et imaged Sgr surveys was berg to sky SDSS et the done the Johnston of damage and al. most (e.g. 2MASS et when Mateo only rate the was 1995; high It Lynden-Bell a its & 1996). at losing Lynden-Bell was 1995; was tides it al. much dwarf Galactic (1994), for Sgr the al. survive the et to to that Ibata going by realized massive, discovery not quickly its however, most After is et the Niederste-Ostholt It longer. of (e.g. 2010). SMC one and al. is LMC the galaxy after third dwarf (Sgr) tarius under- Galaxy’s the stars study distribution. and can matter gas we Way, lying stripped Milky of the streams onto frag- – in galactic manifested of observables. – accretion local cosmological of ments the signatures our the to Galactic bootstrap measuring formation to By the structure chance in of a theories us continues gives infall it halo, . stochastic the into coalescing As are Clouds Magellanic two the Niederste-Ostholt nvriy 90N hre tet atmr,M 11,US 21218, MD Baltimore, Street, Charles N. 3900 University, P University USA Laboratory, 16802, Davey PA 525 University, State vania China 100012, Beijing Sciences, Obs of Academy Astronomical Chinese National toires, China; 100871, Beijing versity, Australia 2006, NSW Sydney, of University Russia 119991, Moscow 13, pr. Universitetskiy H,UK 0HA, fteGlci aeltssriigt ae h Sagit- the date, to surviving satellites Galactic the Of as assembling, finished not clearly has Way Milky The rpittpstuigL using 2018 typeset 6, Preprint July version Draft egyE Koposov E. Sergey 7 6 5 4 3 2 1 eateto hsc srnm,TeJhsHopkins Johns The , & Physics of Department Pennsyl- The Astrophysics, and USA Astronomy 88349, of NM, Department Sunspot, Observatory, Point Apache Uni Peking Astrophysics, and Astronomy for Institute The Kavli A28, Physics, of School Astronomy, for Institute Sydney trbr srnmclIsiue ocwSaeUniversit State Moscow CB3 Institute, Astronomical Cambridge Sternberg Road, Madingley Astronomy, of Institute rnhsasbln anyt h eu temi hscnue loca confused this in stream Cetus Sa the to bo headings: mainly to Subject Using belongs mainly population belong stream. straggler stars blue Sagittarius branch the the that t show crosses Cap, we Galactic which st data, South the Stream, the mix, of Cetus region population the the stellar compared an In the as stars “bifurcation”. and clump Northern red metal-rich properties prominent e structural more a the provide and contains on We turn-offs prop broader, multiple morphological with is stream. distinct diagram stream have metal-poor brighter but more gradients The and distance similar fainter follow another is there la iia k uvyDt ees .Prle oteSgtaist Sagittarius the to Parallel 8. Release Data Survey Sky Digital Sloan h tutr fteSgtaissra nteSuhr aatch Galactic Southern the in stream Sagittarius the of structure The H AITRU TEM NTESUHR AATCHEMISPHER GALACTIC SOUTHERN THE IN STREAMS SAGITTARIUS THE 1. INTRODUCTION A 1 T .Pe J. , 1,2 E tl mltajv 5/2/11 v. emulateapj style X .Belokurov V. , aais wr aais niiul(aitru)–LclGroup Local – (Sagittarius) individual : – dwarf galaxies: narrubia ˜ 1 ..Smith M.C. , 1 ..Evans N.W. , π Schneider ain on radians rf eso uy6 2018 6, July version Draft erva- ABSTRACT 4 ark, .Bizyaev D. , A y, 6 1 - ...Wyse R.F.G. , .Gilmore G. , a k uvy(DS Fkgt ta.19;Gn et Gunn 1996; al. et (Fukugita Digi- (SDSS) pho- Sloan the the Survey the revisit of Sky from we motion available tal Here, previously proper 2012). data the al. tometric measure et Carlin to (cf. 82 stream Sgr Stripe of observations in multi-epoch the stars com- use a stream we In Sgr paper, the substructures. panion halo of stellar formation neighbouring its the and into Sgr. un- insights the an new as report to well crucial as metallicity. stream, is and Cetus South space the the of in in kinematical derstanding coincides different debris a it the has Untangling Sgr, this to Although Newberg Sgr signature by with (2009). hemisphere overlap al. Galactic to et Southern the announced in been debris recently has orbit little that made. seems been now has Mart´ınez-Delgado it et dis- progress 2010), 2006; Sgr al. Pe˜narrubia of al. et 2007; picture et al. the Fellhauer in (e.g these there the ruption Although include to to obscure. close attempts remain lie were origins that 2006) whose substructures al. plane of et Sgr (Belokurov examples tail both leading are the par- the of and In 2008) “bifurcation” (Juri´c al. remain. Over-density et difficulties Virgo (e.g. the but ticular, data 2010a), this Majewski explaining & their have Law success plane, Models considerable orbital some velocities. Sgr had radial the such and to distances as proximity heliocentric mem- classified their Sgr typically on past are based their objects to Instead, according bership. tagged simply not galaxy are Way disrupting Milky 2010b). the the Majewski from in & torn Law free-floating been (e.g. left having now originated after are have halo to 2010). and believed al. Sgr are et clusters Chou in 2009; 2007; of al. al. number et et Yanny A Chou 2007; al. 2004; et chem- al. Monaco reliable et 2006). have veloc- (Majewski also al. radial abundances these their of ical et have some tails Belokurov and Sgr measured 2003; the ities in al. stars et of Hundreds (Majewski sky the hsi h rto w bevtoa aesi hc we which in papers observational two of first the is This polar a on stream) Cetus (the stream another Finally, clusters) star (and stars unbound the Unfortunately, 5,2 .Malanushenko E. , 1 .Gieles M. , 7 eei vrapn ia ersfrom debris tidal overlapping is here emcngrto ssmlrt the to similar is configuration ream hpooercadspectroscopic and photometric th a ihrcolour-magnitude richer a has d itru n h lehorizontal blue the and gittarius ini h halo. the in tion mshr saaye ihthe with analysed is emisphere iec httetostreams two the that vidence risadselrpopulations. stellar and erties dltak but track, idal otefitrsra.Based stream. fainter the to 1 ..Irwin M.J. , 5 .Malanushenko V. , 1 ..Lewis G.F. , ∼ 10 E ◦ away, 1,3 5 D.P. , M. , 2

Fig. 1.— Density of MSTO stars with 0

Λ Λ Λ Λ Λ

− − − − − −

Fig. 2.— Left: Density of MSTO stars (with the same color-magnitude selection as for Figure 1) across the Sgr stream in the south (90◦ < Λ < 120◦) (solid line) and the density of 2MASS M giants in the same region (dotted line). A constant background (around 7000 stars per bin) has been subtracted from the histograms. Right panel: The density of MSTO stars in different slices across the stream from 70◦ < Λ < 80◦ to 110◦ < Λ < 120◦. Although the secondary stream seems to have the same offset from the main stream at different Λ, this is not actually the case as shown in Figure 3.

TABLE 1 Locations and widths of the streams

Λ Bcen,bright Bcen,faint Wbright Wfaint deg deg deg deg deg 75 -1.3 9.5 4.6 2.7 85 -2.5 10.3 5.3 4.1 95 -1.3 10.3 5.8 2.8 105 -0.9 10.7 5.6 2.9 115 -1.4 11.8 5.3 3.1 TABLE 2 − Distances to the bright stream

Λ m − M σ(m − M) − deg mag mag 92.5 16.69 0.06 97.5 16.72 0.06 102.5 16.86 0.04 − 107.5 17.01 0.07 − 112.5 17.17 0.09 Λ 117.5 17.31 0.06 122.5 17.27 0.05 127.5 17.28 0.10 Fig. 3.— Full lines show the centers of the bright and faint streams as a function of longitude Λ as determined by a double Gaussian fit to the profile of Figure 2. The dashed lines show the 1 σ widths of the streams, whilst the dotted lines show the extrap- Z . 0.02 and absolute magnitude 3 . Mi . 6, occupy- olation of the streams’ centroids on approach to the Sgr remnant, ing the range of heliocentric distances 10 . D(kpc) . 60. which is marked by the red star. The density of ∼ 13,000,000 stars that passed the above color and magnitude cuts in the SDSS DR8 dataset is shown in Figure 1 in equatorial and Galactic coordinates 19.5

Red Clump −5

16 16

15 0 15 0 120 114 108 102 96 90 120 114 108 102 96 90 Λ [deg] Λ [deg]

Fig. 5.— Measurement of distances and distance gradients along the stream using two different tracers: sub-giant branch stars and red clump stars. Left: 2D-histogram of red clump star counts as a function of longitude along the stream Λ and i-band mag- nitude. The distance gradient is clearly visible. Right: similar − 2D-histogram for sub-giant stars, selected using a combination of − g−i and g−r colors. The exact same gradient as for the red clump is clearly visible. The red line on both panels shows the distance Fig. 4.— Background subtracted Hess diagram of the Sgr stream gradient of 0.023 mag deg−1 (offset for clarity). in the area defined by 100◦ < Λ < 110◦ and −5◦ 0.85, structure. Watkins et al. (2009) showed that the density 0.22 < J − H − 0.561(J − Ks) < 0.36, 10 < Ks < 12. slice through the Sgr stream in the southern Stripe 82 Both profiles show clear evidence for bimodality, though contains at least two maxima. it is unclear whether the two structures are distinct or It is useful to define a heliocentric coordinate sys- overlapping. The right panel of Figure 2 shows cross- tem aligned with the Sgr stream. Such coordinate sys- sections across the stream in different slices. As we march tems, whose equator aligns with the stream, have already along the stream, at least in the region 90◦ < Λ < 120◦, proved useful in similar studies (e.g. Majewski et al. 2003; the cross-section remains quite invariant, but the offset Koposov et al. 2010). Using the notation of Majewski of the secondary stream from the main stream changes et al. (2003), we introduce coordinates (Λ,B) given by gently with longitude. their eqn (9). The equator of this spherical cooordinate Let us assume that the one-dimensional profile of each system coincides with the Sgr debris midplane. The bot- stream is a Gaussian whose centroid and full width at tom panel of Figure 1 shows the data in this coordinate half maximum may vary with longitude. This simple system, with the Sgr debris now straddling the equator. model of the stellar density in a tidal stream is of course This reveals that, in the new projection, rather than not completely physical, but nevertheless should be suf- a “bifurcation” of the Sgr stream (the somewhat mis- ficient to describe the pieces of the streams in the south. leading term introduced in Belokurov et al. 2006), the We now extract the centroids and widths by fitting two stellar density appears resolved into two, sometimes su- Gaussians to the starcount data. Figure 3 shows their perposed, independent streams with seemingly different behavior as a function of longitude Λ along the stream. angular widths and density profiles. Perhaps unsurprisingly, the tracks of the centers of two We are led to the conclusion that everywhere where the streams are not exactly parallel, but they are slightly Sgr tidal debris can be detected in SDSS there exists at converging. Nonetheless, it is surprising that the conver- The Sgr Streams in the South 5

− − −

Fig. 6.— Hess diagrams determined separately for the main Southern (−5◦

−− −− −− −− −− − −

1 1 − −

− − − − − − − − −

Fig. 7.— Background-subtracted distributions of stars of different metallicities across the streams. In the left panel, the stars have been selected from the Stripe 82 coadded dataset, and in the right panel from the ‘Field of Streams’ (205◦ < Λ < 240◦). The shaded histograms show the density distribution of all MS stars without splitting in metallicity bins. Note that the two streams have different metallicity distribution functions, with the brighter stream containing substantial numbers of metal-rich stars as compared to the fainter stream. This is true both for the streams in the south, and those in the north. The metallicities have been measured using the formula from Bond et al. (2010). different tracer populations in Figure 4. We select red (as used for Fig. 2). With 95 % confidence, this ratio in clump stars by the colour cut (cf. Correnti et al. 2010) the fainter stream is smaller than that in the brighter. − This suggests the existence of a simpler and more metal- 0.8

21 tight structures, whereas the BSs are poorer distance in- dicators and produce more diffuse structures (see e.g., Figure 4 of Deason et al. 2011). 20 BS Figure 8 shows the density of stars satisfying our colour cut as a function of apparent magnitude and longitude. There are indeed two evident structures present, but they 19 do not correspond to two roughly parallel density peaks which might be naturally interpreted as BS and BHBs from the same structure. Rather the figure forces upon 18 BHB us the interpretation that there are two classes of objects arising from two distinct structures: the thick structure g [mag] possesses the same longitudinal distance gradient as the 17 subgiant and red clump stars, and this can be identified with BSs from the Sgr stream. The other, thinner struc- ture has the opposite distance gradient, and we shall see 16 that these are BHB stars lying in an entirely different structure, namely the Cetus stream. The discovery of the Cetus stream was announced by 15 Newberg et al. (2009). They noticed a stream-like over- 160 140 120 100 80 60 density in low metallicity stars in SDSS DR7 that crosses ◦ Λ [deg] the Sgr stream in the south at b ∼ −70 . The Cetus stream can be distinguished from the Sgr on the ba- Fig. 8.— Magnitude distribution of BHB/BS-like stars as a func- tion of the angle along the stream. The diagram shows two classes sis of its markedly lower ratio of BS to BHB stars and of objects coming from two structures: BS from Sgr stream with its different kinematics. They also suggested that some the same distance gradient as the subgiants/red clump (shown as BHB stars previously attributed to the Sgr stream in- an offset red line), and BHBs from Cetus stream with an opposite distance gradient at g ≈ 18 and 80◦ . Λ . 120◦. stead properly belong to the Cetus stream. In the middle panel of Figure 9, we display the density distribution of stars on the sky satisfying the color cuts −0.3 < g −r < 0 and 0.9 < u−g < 1.35 together with TABLE 3 18.5

4 4 4 4 4 Λ Λ

Fig. 9.— Left: Density of high Galactic latitude (|b| > 20◦) MSTO stars selected using the same color-magnitude cuts as for Figure 1. The arrow shows the direction of motion of Sgr stars. Middle: Density of high Galactic latitude BS stars −0.3 < g−r < 0, 0.9

4 BHB BHB − −

4 4 −−− − − − − Λ −1 GSR Fig. 10.— Left: Density of Cetus stream BHB stars on the sky. The stars were selected to be BHB-like according to standard g −r, u−g cuts, Galactic latitude cut (|b| > 20◦), together with a position-dependant distance cut selecting predominantly Cetus stream stars. The arrow shows our measured distance gradient of Cetus stream BHBs. Middle: BHB stars in a coordinate system aligned with the Cetus stream, which is now clearly visible as an extended, narrow structure. The red line shows the distance gradient of 0.013 mag deg−1 offset for clarity. Right: The kinematical separation of the BS+BHB stars belonging to Sgr and Cetus stream for 80◦ < Λ < 120◦ and −10◦

We used red clump and subgiant stars to measure the model seems to be ruled out as it has been difficult to find distance gradients along the streams. This enables us strong evidence for residual rotation in the remnant of to construct composite Hess diagrams to study the stel- the Sgr dwarf (Pe˜narrubia et al. 2011). Moreover, in the lar populations in the streams. The brighter stream scenario of Pe˜narrubia et al. (2010) the two streams are shows evidence for multiple turn-offs and a prominent not expected to differ significantly in their stellar popu- red clump, whereas the secondary stream does not. This lations content. This prediction is somewhat difficult to suggests that the brighter stream is composed of more reconcile with the new data on the stellar populations of than one stellar population, and contains a significant the streams. In some other models although, such as Law number of metal-rich stars, much like the Sgr remnant & Majewski (2010a) the streams that form at consecu- itself. By contrast, the fainter stream is dominated by a tive pericentric passages could have somewhat different metal-poor population. This conclusion is also supported metallicity distributions, as metal-rich stars could have by our photometric metallicities computed for the region been torn predominantly later, from deeper within the where the Sgr streams cross the ultra-deep SDSS Stripe gravitational potential of the satellite. 82 coadded data. Finally, it is tempting to suggest that the two streams Our analysis of the new data allows us to untangle a with different properties have actually originated from complicated mix of tidal debris around the South Galac- two different progenitors. The infall of satellites in tic Cap, where the Sgr streams are crossed by the Cetus groups is not particular but general, as best illustrated Stream at an angle of ∼ 30◦. The Sgr and Cetus streams by the recent arrival of the Large and Small Magellanic have similar distances, though their stellar populations Clouds. Cosmological simulations of structure formation and kinematics are different. The structures are not part also find plenty of evidence for group infall (e.g. Li & of the same disruption event, as the Cetus stream is coun- Helmi 2008). The picture painted by the data looks more terrotating with respect to the Sgr. On the basis of their complex than this and remains a challenge to understand. density distribution and kinematics, we have shown that Although there has been substantial progress in mod- most of the BSs belong to the Sgr stream, whereas most elling the Sgr (see e.g., Fellhauer et al. 2006; Pe˜narrubia of the BHBs belong to the Cetus Stream, in this part of et al. 2010; Law & Majewski 2010a), it remains true the sky. The BS to BHB ratio in the two streams is strik- that there is no explanation of the nature of the two ingly different, as already pointed out by Newberg et al. Sgr streams. (2009). Good spectroscopic coverage of the area is essen- tial to disentangle the multiple overlapping streams with different chemical properties, following different distance gradients. SK acknowledges financial support from the Science The work in this paper has extended the ‘Field of and Technology Funding Council of the United King- Streams’ (Belokurov et al. 2006) to the south. The new dom, whilst VB and MG thanks the Royal Society for the imaging data show that, just as in the north, the Sgr award of a University Research Fellowship. GFL thanks stream is accompanied by a fainter stream. The right- the Australian Research Council for support through his hand panel in Fig 7 shows that the population mix in Future Fellowship (FT100100268) and Discovery project the fainter stream of the “bifurcation” around the North (DP110100678). Most of the data processing has been Galactic Cap does not contain as many metal-rich stars done using the python programming language and the as the brighter stream. The simplest explanation is that following open source modules: numpy9, scipy10, mat- the southern faint stream is part of the same structure plotlib11. The SDSS data was accessed using local SDSS as the northern faint stream. copy stored in PostgreSQL database and Q3C module These results raise the question: Is it possible to pro- (Koposov & Bartunov 2006). Funding for SDSS-III has duce the streams with the properties described in this been provided by the Alfred P. Sloan Foundation, the paper in a disruption of one galaxy or is more than one Participating Institutions, the National Science Founda- progenitor necessary? tion, and the U.S. Department of Energy Office of Sci- Recently, two possible scenarios, both in context of the ence. The SDSS-III web site is http://www.sdss3.org/. Sgr dSph disruption, have shown how to form two al- SDSS-III is managed by the Astrophysical Research most parallel tidal streams from the debris of one parent Consortium for the Participating Institutions of the galaxy. Although Fellhauer et al. (2006) did not fully SDSS-III Collaboration including the University of Ari- match the available data, they did suggest that branches zona, the Brazilian Participation Group, Brookhaven Na- A and B (see Belokurov et al. 2006) could be reproduced tional Laboratory, University of Cambridge, University by multiple wraps of the same stream offset on the plane of Florida, the French Participation Group, the German of the sky by small amount of differential precession. In Participation Group, the Instituto de Astrofisica de Ca- this picture, the fainter stream in the North Galactic Cap narias, the Michigan State/Notre Dame/JINA Participa- area corresponds to the dynamically old tidal debris in tion Group, Johns Hopkins University, Lawrence Berke- the Sgr trailing arm. However, this mechanism does not ley National Laboratory, Max Planck Institute for As- produce two distinct streams in the south and now seems trophysics, New Mexico State University, New York Uni- to be ruled out. versity, Ohio State University, Pennsylvania State Uni- Alternatively, Pe˜narrubia et al. (2010) point out that if versity, University of Portsmouth, Princeton University, the Sgr progenitor had a rotating stellar disk misaligned the Spanish Participation Group, University of Tokyo, with respect to its orbit, then stripping naturally pro- University of Utah, Vanderbilt University, University of duces bifurcated debris tracks on the sky. However, the Virginia, University of Washington, and Yale University.

9 http://numpy.scipy.org 10 http://www.scipy.org 10

REFERENCES Aihara, H., Allende Prieto, C., An, D., et al. 2011, ApJS, 193, 29 Koposov, S., & Bartunov, O. 2006, Astronomical Data Analysis An, D., Johnson, J. A., Clem, J. L., et al. 2008, ApJS, 179, 326 Software and Systems XV, 351, 735 Annis, J., Soares-Santos, M., Strauss, M. A., et al. 2011, Koposov, S. E., Rix, H.-W., & Hogg, D. W. 2010, ApJ, 712, 260 arXiv:1111.6619 Law, D. R., & Majewski, S. R. 2010a, ApJ, 714, 229 Bellazzini, M., Newberg, H. J., Correnti, M., Ferraro, F. R., & Law, D. R., & Majewski, S. R. 2010b, ApJ, 718, 1128 Monaco, L. 2006, A&A, 457, L21 Li, Y.-S., & Helmi, A. 2008, MNRAS, 385, 1365 Belokurov, V., Zucker, D. B., Evans, N. W., et al. 2006, ApJ, 642, Lynden-Bell, D., & Lynden-Bell, R. M. 1995, MNRAS, 275, 429 L137 Majewski, S. R., Skrutskie, M. F., Weinberg, M. D., & Bond, N. A., Ivezi´c, Z.,ˇ Sesar, B., et al. 2010, ApJ, 716, 1 Ostheimer, J. C. 2003, ApJ, 599, 1082 Bovy, J., Hogg, D. W., & Rix, H.-W. 2009, ApJ, 704, 1704 Majewski, S. R., Kunkel, W. E., Law, D. R., et al. 2004, AJ, 128, Carlin, J. L., Majewski, S. R., Casetti-Dinescu, D. I., et al. 2012, 245 ApJ, 744, 25 Mart´ınez-Delgado, D., Pe˜narrubia, J., Juri´c, M., Alfaro, E. J., & Chou, M.-Y., Majewski, S. R., Cunha, K., et al. 2007, ApJ, 670, Ivezi´c, Z. 2007, ApJ, 660, 1264 346 Mateo, M., Mirabal, N., Udalski, A., et al. 1996, ApJ, 458, L13 Chou, M.-Y., Cunha, K., Majewski, S. R., et al. 2010, ApJ, 708, Monaco, L., Bellazzini, M., Bonifacio, P., et al. 2007, A&A, 464, 1290 201 Correnti, M., Bellazzini, M., Ibata, R. A., Ferraro, F. R., & Newberg, H. J., Yanny, B., Rockosi, C., et al. 2002, ApJ, 569, 245 Varghese, A. 2010, ApJ, 721, 329 Newberg, H. J., Yanny, B., Grebel, E. K., et al. 2003, ApJ, 596, Co¸skunoˇglu, B., Ak, S., Bilir, S., et al. 2011, MNRAS, 412, 1237 L191 Deason, A.J., Belokurov V., Evans N.W., 2011, MNRAS, 416, Newberg, H. J., Yanny, B., & Willett, B. A. 2009, ApJ, 700, L61 2903 Niederste-Ostholt, M., Belokurov, V., Evans, N. W., & Eisenstein, D. J., Weinberg, D. H., Agol, E., et al. 2011, AJ, 142, Pe˜narrubia, J. 2010, ApJ, 712, 516 72 Pe˜narrubia, J., Belokurov, V., Evans, N. W., et al. 2010, Fellhauer, M., Belokurov, V., Evans, N. W., et al. 2006, ApJ, 651, MNRAS, 408, L26 167 Pe˜narrubia, J., Zucker, D. B., Irwin, M. J., et al. 2011, ApJ, 727, Fukugita, M., Ichikawa, T., Gunn, J. E., et al. 1996, AJ, 111, 1748 L2 Girardi, L., Grebel, E. K., Odenkirchen, M., & Chiosi, C. 2004, Prior, S. L., Da Costa, G. S., & Keller, S. C. 2009, ApJ, 704, 1327 A&A, 422, 205 Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, Gunn, J. E., Carr, M., Rockosi, C., et al. 1998, AJ, 116, 3040 525 Gunn, J. E., Siegmund, W. A., Mannery, E. J., et al. 2006, AJ, Totten, E. J., & Irwin, M. J. 1998, MNRAS, 294, 1 131, 2332 Watkins, L. L., Evans, N. W., Belokurov, V., et al. 2009, Hess, R., ”Die Verteilungsfunktion der absol. Helligkeiten etc.”. MNRAS, 398, 1757 Probleme der Astronomie. Festschrift fur Hugo v. Seeliger. Yanny, B., Newberg, H. J., Kent, S., et al. 2000, ApJ, 540, 825 Springer, Berlin. p. 265. Yanny, B., Newberg, H. J., Johnson, J. A., et al. 2009, ApJ, 700, Ibata, R. A., Gilmore, G., & Irwin, M. J. 1994, Nature, 370, 194 1282 Johnston, K. V., Spergel, D. N., & Hernquist, L. 1995, ApJ, 451, York, D. G., Adelman, J., Anderson, J. E., Jr., et al. 2000, AJ, 598 120, 1579 Juri´c, M., Ivezi´c, Z.,ˇ Brooks, A., et al. 2008, ApJ, 673, 864

11 http://matplotlib.sf.net