DOCTORAL THESIS Arxiv:1312.1643V1 [Astro-Ph.GA] 5
Total Page:16
File Type:pdf, Size:1020Kb
Charles University in Prague Faculty of Mathematics and Physics DOCTORAL THESIS Ivana Ebrov´a Shell galaxies: kinematical signature of shells, satellite galaxy disruption and dynamical friction Astronomical Institute of the Academy of Sciences of the Czech Republic Supervisor of the doctoral thesis: RNDr. Bruno Jungwiert, Ph.D. arXiv:1312.1643v1 [astro-ph.GA] 5 Dec 2013 Study program: Physics Specialization: Theoretical Physics, Astronomy and Astrophysics Prague 2013 This research has made use of NASA's Astrophysics Data System, micronised purified flavonoid fraction, and a lot of iso-butyl-propanoic-phenolic acid. Typeset in LYX, an open source document processor. For graphical presentation, we used Gnuplot, the PGPLOT (a graphics subroutine library written by Tim Pearson) and scripts and programs written by Miroslav Kˇr´ıˇzekusing Python and matplotlib. Calculations and simulations have been carried out using Maple 10, Wolfram Mathematica 7.0, and own software written in pro- gramming language FORTRAN 77, Fortran 90 and Fortran 95. The software for simulation of shell galaxy formation using test particles are based on the source code of the MERGE 9 (written by Bruno Jungwiert, 2006; unpublished); kinematics of shell galaxies in the framework of the model of radial oscillations has been studied using the smove software (written by Lucie J´ılkov´a,2011; unpublished); self-consistent simulations have been done by Kateˇrina Bartoˇskov´awith GADGET-2 (Springel, 2005). We acknowledge support from the following sources: grant No. 205/08/H005 by Czech Science Foundation; research plan AV0Z10030501 by Academy of Sciences of the Czech Republic; and the project SVV-267301 by Charles University in Prague. This work has been done with the support for a long-term development of the research institution RVO67985815. 2 I declare that I carried out this doctoral thesis independently, and only with the cited sources, literature and other professional sources. I understand that my work relates to the rights and obligations under the Act No. 121/2000 Coll., the Copyright Act, as amended, in particular the fact that the Charles University in Prague has the right to conclude a license agreement on the use of this work as a school work pursuant to Section 60 paragraph 1 of the Copyright Act. In Prague, 19. 8. 2013 Ivana Ebrov´a 3 Title: Shell galaxies: kinematical signature of shells, satellite galaxy disruption and dy- namical friction Author: Ivana Ebrov´a Department / Institute: Astronomical Institute of the Academy of Sciences of the Czech Republic Supervisor of the doctoral thesis: RNDr. Bruno Jungwiert, Ph.D., Astronomical Institute of the Academy of Sciences of the Czech Republic Abstract: Stellar shells observed in many giant elliptical and lenticular as well as a few spi- ral and dwarf galaxies presumably result from radial minor mergers of galaxies. We show that the line-of-sight velocity distribution of the shells has a quadruple-peaked shape. We found simple analytical expressions that connect the positions of the four peaks of the line profile with the mass distribution of the galaxy, namely, the circular velocity at the given shell radius and the propagation velocity of the shell. The analytical expressions were applied to a test-particle simulation of a radial minor merger, and the potential of the simulated host galaxy was successfully recovered. Shell kinematics can thus become an independent tool to determine the content and distribution of dark matter in shell galaxies up to ∼100 kpc from the center of the host galaxy. Moreover we investigate the dynamical friction and gradual disruption of the cannibalized galaxy during the shell formation in the framework of a simulation with test particles. The coupling of both effects can considerably redistribute positions and luminosities of shells. Neglecting them can lead to significant errors in attempts to date the merger in observed shell galaxies. Keywords: galaxies: kinematics and dynamics, galaxies: interactions, galaxies: evolution, methods: analytical and numerical 4 Contents 1 Objectives and motivation 9 I Introduction 11 2 Shell galaxies in brief 11 3 Observational knowledge of shell galaxies 12 3.1 Observational history . 12 3.2 Occurrence of shell galaxies . 12 3.3 Appearance of the shells . 14 3.4 Colors . 18 3.5 Gas and dust . 19 3.6 Radio and infrared emission . 20 3.7 Other features of host galaxies . 21 4 Summary of shell characteristics 22 5 Scenarios of shells' origin 23 5.1 Gas dynamical theories . 23 5.2 Weak Interaction Model (WIM) . 24 6 Merger model 28 6.1 Phase wrapping . 28 6.2 Cannibalized galaxy . 29 6.3 Ellipticity of the host galaxy . 30 6.4 Radial distribution of shells . 30 6.5 Radiality of the merger . 32 6.6 Major mergers . 33 6.7 Simulations with gas . 34 6.8 Merger model and observations . 34 7 Measurements of gravitational potential in galaxies 36 7.1 Insight into methods . 36 7.2 Use of shells . 36 II Shell kinematics 39 8 Preliminary provisions 39 8.1 Host galaxy potential model . 39 8.2 Terminology . 40 8.3 Quantities . 41 5 9 Model of radial oscillations 43 9.1 Turning point positions and their velocities . 43 9.2 Real shell positions and velocities . 44 9.3 Appearance of the shells . 46 9.4 Kinematics of shell stars . 47 9.5 Characteristics of spectral peaks . 49 9.6 Equations of LOSVD . 52 9.7 Shell-edge density distribution and LOSVD . 53 9.8 Nature of the quadruple-peaked profile . 54 10 Stationary shell 56 10.1 Motion of stars in a shell system . 57 10.2 Constant acceleration . 58 10.3 LOSVD . 59 10.4 Comparison with the model of radial oscillations . 60 11 Constant acceleration and shell velocity 61 11.1 Motion of a star in a shell system . 62 11.2 Approximative LOSVD . 63 11.3 Radius of maximal LOS velocity . 63 11.4 Approximative maximal LOS velocity . 64 11.5 Slope of the LOSVD intensity maxima . 64 11.6 Comparison of approaches . 65 11.7 Projection factor approximation . 68 12 Higher order approximation 70 12.1 Motion of a star in a shell system . 70 12.2 Comparison of approximations . 72 12.3 a1 ......................................... 74 13 Test-particle simulation 75 13.1 Parameters of the simulation . 75 13.2 Comparison of the simulation with models . 76 13.3 Recovering the potential from the simulated data . 80 13.4 Notes about observation . 84 14 Shell density 86 14.1 Projected surface density of the shell edge . 86 14.2 Time evolution . 87 14.3 Volume density . 89 14.4 Projected surface density . 90 15 Discussion 92 III Dynamical friction and gradual disruption 95 16 Motivation 95 6 17 Description of simulation 96 17.1 Configuration . 96 17.2 Plummer sphere . 96 17.3 Velocity dispersion in Plummer potential . 98 17.4 Velocity dispersion in a double Plummer sphere . 99 17.5 Standard set of parameters . 100 18 Dynamical friction 101 19 Multiple Three-Body Algorithm (MTBA) 101 19.1 Principle and characteristics . 101 19.2 Merger parameters . 102 19.3 Results of simulations . 103 20 Comparison with self-consistent simulations 104 20.1 Altering GADGET-2 computational setting . 105 20.2 Comparison of methods . 106 21 Tidal disruption 108 21.1 Massloss of the secondary . 108 21.2 Deformation of the secondary galaxy . 111 22 Simulations of shell structure 113 22.1 Dynamical friction and tidal disruption . 113 22.2 Dark halo . 115 22.3 Self-consistent versus test-particle simulations . 119 23 Discussion 121 IV Conclusions 124 V Appendix 127 A Units and conversions 127 B List of abbreviations 128 C Initial velocity distribution 129 D Introduction to dynamical friction 130 D.1 A thermodynamic meditation . 130 D.2 Chandrasekhar formula . 130 D.3 What a wonderful universe . 132 D.4 Whydoesitwork? ................................ 133 7 E Our method 135 E.1 Avoiding some approximations . 135 E.2 Back to Chandrasekhar formula . 137 E.3 Incorporation of the friction in the simulation . 139 F Tidal radius 140 G Expressions for the tidal radius 143 H Videos 144 References 146 8 1 Objectives and motivation The most successful theory of the evolution of the Universe so far seems to be the theory of the hierarchical formation based on the assumption of the existence of cold dark matter, significantly dominating the baryonic one. In such a universe, large galaxies are formed by merging of small galaxies, protogalaxies and diffuse accretion of surrounding matter. Galactic interaction and dark matter play thus a crucial role in the life of every galaxy. But the determination of both the dark matter content and the merger history of a galaxy is difficult. Firstly, the cold dark matter interacts only gravitationally (and possibly via the weak interaction) and thus the mapping of its distribution in galaxies is tricky. Secondly, the nature disallows us to see individual galaxies from different angles, thus our knowledge of their spatial properties is degenerate. Thirdly, it is non-trivial to determine anything about the history of a given galaxy as the whole existence of humanity presents only a snapshot in the evolution of the Universe. Yet this knowledge is important to confirm or disprove theories of the creation and evolution of the Universe, improve their accuracy and to understand how the Universe we live in actually looks. The deal of the galactic astronomy is to try to circumvent these obstacles. One of the possibilities is to use tidal features left by the galactic interactions. They act as dynamical tracers of the potential of their host galaxies and as hints left behind by the accreted galaxies in the past. The special case is that of arc-like fine structures found in shell galaxies. Their unique kinematics carries both qualitative and quantitative information on the distribution of the dark matter, the shape of the potential of the host galaxy and its merger history.