Arxiv:1710.08935V2 [Astro-Ph.CO] 27 Oct 2017
Total Page:16
File Type:pdf, Size:1020Kb
Submitted to ApJ Preprint typeset using LATEX style AASTeX6 v. 1.0 ACTION DYNAMICS OF THE LOCAL SUPERCLUSTER Edward J. Shaya Astronomy Department, University of Maryland, College Park, MD 20742, USA R. Brent Tully Institute for Astronomy, University of Hawaii, Honolulu, HI 96822, USA Yehuda Hoffman Racah Institute of Physics, Hebrew University, Jerusalem 91904, Israel Daniel Pomarede` Institut de Recherche sur les Lois Fondamentales de l’Univers, CEA, Universit´eParis-Saclay, F-91191 Gif-sur-Yvette, France (Received 2017, August 16) ABSTRACT The fully nonlinear gravitationally induced trajectories of a nearly complete set of galaxies, groups, and clusters in the Local Supercluster are constructed in a Numerical Action Method (NAM) model constrained by data from the CosmicFlows survey and various distance indicators. We add the gravity field due to inhomogeneities external to the sample sphere by making use of larger scale peculiar flow measurements. Assignments of total masses were made to find the best overall set of mutual attractions, as determined by a goodness criterion based on present day radial velocities, individually for the Virgo Cluster, M31, and the Milky Way (MW), and via a mass-to-light ratio relationship for other masses. The low median chi-square found indicates the model fits the present day velocity flow well, but a slightly high mean chi-square may indicate that some masses underwent complex orbits. The best fit, when setting the value of H0 to the CosmicFlows value of 75 km/s/Mpc and the WMAP value for Ωm = 0.244 consistent with that H0, occurs with the following parameters: 0.15 10 Ωorphan =0.077 0.016, M/LK = 40 2L10 M⊙/L⊙ (L10 is K-band luminosity in units of 10 L⊙), ± 14 ± a Virgo mass of 6.3 0.8 10 M⊙ (M/LK = 113 15M⊙/L⊙), and a mass for the MW plus M31 12± × ± of 5.15 0.35 10 M⊙. The best constant mass-to-light ratio is M/L = 58 3M⊙/L⊙. The ± × K ± Virgocentric turnaround radius is 7.3 0.3 Mpc. We explain several interesting trends in peculiar ± motions for various regions now that we can construct the 3D orbital histories. Keywords: large-scale structure arXiv:1710.08935v2 [astro-ph.CO] 27 Oct 2017 1. INTRODUCTION cluster dynamics date back 36 years (Tonry & Davis The large scale overdensity region of the universe that 1981; Aaronson et al. 1982). The Virgo Cluster is the we reside in was called the Virgo Supergalaxy and later main condensation in this region and its collective mass the Local Supercluster by de Vaucouleurs (1953, 1956). causes substantial peculiar motions. The Centaurus and We focus on a domain of 3000 km s−1 40 Mpc. Hydra clusters, just beyond our sample volume, and sev- ∼ ∼ This article, on the dynamics of the region, is comple- eral other rich clusters collectively make up the “Great mentary to an article on the groups and larger scale as- Attractor” that is a significant source of a large scale sociations found in the same region (Kourkchi & Tully flow enveloping the study region (Dressler et al. 1987a; 2017). The first extensive discussions of Local Super- Lynden-Bell et al. 1988); see Hoffman et al. (2017) for an accounting of attractors and repellers that affect lo- cal bulk motion. It has become evident that the so-called [email protected] Local Supercluster is only an appendage of the Great At- 2 Shaya, Tully, Hoffman & Pomarede` tractor complex which, in the fullness of its basin of at- used in the past for the nonlinear regime near clusters, traction is called the Laniakea Supercluster (Tully et al. has been to assume spherical infall (Hoffman et al. 2014). Nearby components that call for our attention in- 1980; Tully & Shaya 1984; Ekholm et al. 1999; clude the Fornax Cluster and adjacent Eridanus Cloud Karachentsev et al. 2014) in which each shell of that de Vaucouleurs (1956) referred to as the Southern matter follows the Friedmann equations of a standard Supergalaxy, and the Local Void (Tully & Fisher 1987) universe of the same age, density, and ΩΛ. The actual that is actually part of an extensive network of voids. situation is generally more complex than spherical Our interest is in galaxy flow patterns, departures infall, therefore realistic models need to be more elab- from cosmic expansion. The radial component of these orate. Peebles (1989) demonstrated a way forward for motions are decoded from galaxy distance measures and the more general case with the approach initially called the degree to which they deviate from the mean Hubble ‘least action’ but now called the ‘numerical action law expectations. method’ (NAM) (Peebles et al. 2001; Phelps et al. Distance estimates have become plentiful with the re- 2006; Shaya & Tully 2013). The solutions are paths lease of the Cosmicflows-3 compendium of 18,000 galax- at either extrema or saddle points of the action, the ies (Tully et al. 2016). The almost 400 nearby, high pre- integral of the Lagrangian over time, as test paths cision, measures based on the tip of the red giant branch are varied arbitrarily. NAM works with the comoving method (Jacobs et al. 2009) are particularly important Lagrangian for collisionless particles, ie, their kinematic for the present discussion. We are able to recover the minus potential energies transformed into comoving mass distribution from the observed motions of galaxies coordinates. even though orbits in high density regions can be com- As in the usual celestial mechanics problem, six phase- plex. Orbits manifest curvature and, across much of the space constraints are needed per particle to solve for par- local volume, dynamics are in the non-linear regime, i.e. ticle motions influenced by their mutual gravity. Since, peculiar velocities are not linearly related to the over- at early times, peculiar velocities were small (v a3/2), p ∝ densities. consistent with a smooth homogeneous universe seen in At still higher densities, there are shell-crossings the cosmic microwave background, NAM presumes for (Shandarin et al. 2012), however such situations are too the first time step that peculiar velocities are consistent complicated to be disentangled on the scales of current with linear theory. Since we now have fairly accurate interest. Therefore, this work considers collapsed knots distances to many galaxies, the other three constraints as indivisible units and refers to them as “masses.” can be fulfilled by the present 3D positions. Before discussing our methodology, we mention other In NAM, each path is the path of the center of mass approaches to the mapping of mass in large-scale of the atoms of each final assemblage which, according regions. Weak gravitational lensing ((Kaiser 1992; to the center of mass theorem, moves independently of Kaiser & Squires 1993) provide mass projected onto the internal motions, including inelastic collisions such the plane of the sky, and studies of cluster caustics as mergers. Hence, knowledge of the detailed history (Diaferio & Geller 1997; Rines et al. 2003) provide to- of assembly of each galaxy or group is, to first order, tal masses within some rich cluster. not needed. By assuming point or spherical particles, On large scales, where the fractional overdensities are we loose high order terms in the gravity field associated low, with non-spherical shapes. These are only important ρ ρ¯ δ − << 1 (1) during very close approaches or at very early times when ≡ ρ¯ overdensities were adjacent. Our paths begin at redshift a simple linear relation holds between the divergence of z = 4, before the present constituents were fully formed peculiar velocities, vp, and the distribution of matter but after significant mass collapse had taken place. (Dekel et al. 1990) It is now evident that a significant component of pe- culiar velocities can be generated by large scale struc- ∇ vp = H0fδ (2) · − ture out to great distances. Nearby, the repulsion whereρ ¯ is the mean global density, and the dimension- of the Local Void causes flow perturbations that are less velocity factor f is defined in terms of the density as important as Virgo attraction (Tully et al. 2008; growth factor D, approximated by Lahav et al. (1991) Karachentsev et al. 2015). While Virgo infall might be as: our primary focus, a holistic approach is needed to attain a realistic history of any part of the Local Supercluster. ˙ a D 0.6 ΩΛ 1 As a first step, we construct paths over a region of 38 f = Ω + (1 + Ωm). (3) a˙ D ≈ m 70 2 Mpc in radius centered on ourselves, including what we However, within a supercluster, fractional overden- anticipate to be the full region of infall around the Virgo sities are not typically well below unity. One option, Cluster and all pertinent nearby attractors. This region Action Dynamics of the Local Supercluster 3 roughly coincides with what has historically been called turnaround radius today? Where is the gravitationally the Local Supercluster. bounded limit, and what fraction of the total super- The Local Void extends to 5, 000 km s−1 and only cluster volume is associated with these domains? What ∼ if our model includes the absence of matter on this fraction of matter is presently unassociated with galax- scale will the expansion of the Local Void boundary ies, either still in the primordial intergalactic medium have the correct magnitude. In addition, the align- or tossed out by gravitational encounters, ram pressure ment of major masses, from the Shapley Concentration or internal explosive activity? Then, what are the dom- (Raychaudhury 1989; Scaramella et al. 1989) through inant flow patterns besides Virgo infall? Streaming to- the Norma Centaurus Hydra complex (Great Attrac- ward the Great Attractor and away from the Local Void − − tor) (Lynden-Bell et al.