Rochester Institute of Technology RIT Scholar Works
Theses
7-30-2020
Galaxies as Probes of the Universe: From Dynamical Models to the Inference of Cosmological Parameters
Peter A. Craig [email protected]
Follow this and additional works at: https://scholarworks.rit.edu/theses
Recommended Citation Craig, Peter A., "Galaxies as Probes of the Universe: From Dynamical Models to the Inference of Cosmological Parameters" (2020). Thesis. Rochester Institute of Technology. Accessed from
This Thesis is brought to you for free and open access by RIT Scholar Works. It has been accepted for inclusion in Theses by an authorized administrator of RIT Scholar Works. For more information, please contact [email protected]. Galaxies as Probes of the Universe: From Dynamical Models to the Inference of Cosmological Parameters Peter A. Craig
AThesisSubmittedinPartialFulfillmentofthe Requirements for the Degree of Master of Science in Astrophysical Sciences & Technology
School of Physics and Astronomy College of Science Rochester Institute of Technology Rochester, NY July 30, 2020
Approved by: Andrew Robinson, Ph.D. Date Director, Astrophysical Sciences and Technology
Approval Committee: Sukanya Chakrabarti , Stefi Baum , Zoran Ninkov ASTROPHYSICAL SCIENCES AND TECHNOLOGY COLLEGE OF SCIENCE ROCHESTER INSTITUTE OF TECHNOLOGY ROCHESTER, NEW YORK
CERTIFICATE OF APPROVAL
MASTERS DEGREE THESIS
The Masters Degree Thesis of Peter Craig has been examined and approved by the thesis committee as satisfactory for the thesis requirement for the Masters degree in Astrophysical Sciences and Technology.
Dr. Sukanya Chakrabarti, Thesis Advisor
Dr. Stefi Baum
Dr. Zoran Ninkov
Date______P.A.Craig Master’sThesis
Contents
1 Introduction 1
2 Moving Groups 2 2.1 Introduction ...... 2 2.2 Simulations ...... 4 2.3 DetectionAlgorithm ...... 8 2.4 Discussion and Results ...... 12 2.4.1 GroupFindingintheMilkyWay ...... 20 2.4.2 Comparison of Simulations with the Milky Way ...... 22 2.5 Conclusions and Future Directions ...... 25
3 Magellanic Stream 27 3.1 Introduction ...... 27 3.2 Methods ...... 29 3.2.1 Codes ...... 29 3.2.2 Resolution Study ...... 30 3.2.3 MeasuringtheMassoftheMW ...... 31 3.3 Discussion and Results ...... 33 3.3.1 Galaxy Clusters ...... 40 3.4 Conclusions and Future Directions ...... 42
4LCO 42 4.1 Introduction ...... 42 4.2 Methods ...... 44 4.2.1 Di↵erenceImages...... 46 4.3 Discussion and Results ...... 47
i P.A.Craig Master’sThesis
4.4 Conclusions and Future Directions ...... 49
ii P.A.Craig Master’sThesis
Abstract
Presented here are a number of projects on topics relating to the dynamics and prop- erties of the Milky Way (MW) and measuring the Hubble constant. First is a project studying the possibility of dwarf galaxy interactions forming moving groups in the solar neighborhood. This is done using several simulations of the MW with dwarf galaxies on di↵erent orbits, and searching for velocity substructures in these simulations. The resulting structures are often similar to those observed in the MW, while such structures appear much less frequently in simulations without a perturber. The second project uses adi↵erentsetofsimulationstostudytheformationoftheMagellanicStream.These simulations make some indications about how the stream most likely formed. They also can be used to look at which models for the MW produce a stream most like that which is observed, from which we can estimate the mass of the MW. Another project aims at resolving the tension over the value of the Hubble constant. The goal is to provide a measurement of this parameter that is independent of the methods that produce this tension, and thus avoids any errors in those methods.
iii P.A.Craig Master’sThesis
1 Introduction
Presented here is my work on three di↵erent projects that I have worked on over the last two years. Most of the work I have done focuses on the Milky Way (MW) and the dynamics of the Milky Way, except for one project. This project is discussed in section
4, and it is an e↵ort to measure the Hubble constant (H0). This project is the only one here that does not involve the MW in some fashion. The first project, discussed in section 2, is an analysis of structures in the MW known as moving groups. These groups are structures found in velocity space for stars in the solar neighborhood. Each group is comprised of a number of stars that have very similar velocities. These groups are spread across the velocity space of local stars. We analyze a possible formation method for these structures as a result of interactions between the MW and a dwarf galaxy. In section 3 is my work on a project simulating the formation of the Magellanic Stream (MS). This is a stream of gas extending between two satellites of the MW, the Large and Small Magellanic Clouds (LMC and SMC). We use hydrodynamic and mhd simulations to model the formation of this structure, in an attempt to understand how the stream was formed. We also hoped to use the results of this to estimate the mass of the MW based on the observed properties of the stream. Section 4, as mentioned above, is a project designed to measure the Hubble constant. This is a result of the recent tension in the value of this parameter between the two main methods for measuring H0.Ourgoalistoprovideanindependentmeasure which can escape any possible biases that have been introduced into other measurements of H0.Thisprojectissimilartoe↵ortsbyanumberofothersworkingondi↵erent independent measurements, none of which have resolved the tension so far.
1 P.A.Craig Master’sThesis
2 Moving Groups
2.1 Introduction
Moving groups have been observed in the local region of the Milky Way for many years. Many of these groups are quite well known, including the Hyades, Pleiades, Coma Ber- nices and Sirius moving groups (Eggen, 1965). Parallaxes and proper motions from the Hipparcos satellite (ESA, 1997) were particularly useful for identifying disk moving groups (Dehnen, 1998). Gaia has significantly increased our ability to search for moving groups near the sun (Gaia Collaboration et al., 2018). Recent surveys have also added the line-of-sight velocity to the proper motions to provide full 3D velocity information for stars in the Solar neighborhood (Gontcharov, 2006), which is needed to search for moving groups. Several mechanisms for the formation of moving groups have been suggested. A common explanation is that these moving groups are the remnants of open clusters, or are formed by interactions with a bar (Dehnen, 2000a). One problem with the cluster formation idea is that stars in moving groups can have a variety of di↵erent ages and compositions, so it is unlikely that they all came from the same cluster. Analysis of GALAH data (Quillen et al., 2018b) indicates that some moving groups, such as the Hercules moving group, may be due to a resonant bar. It has also been suggested that moving groups could have been formed from perturbations due to the Magellanic Clouds via gravitational interactions (Dehnen, 1998). Recent work also finds that transient spiral structure (Hunt et al., 2018) may lead to the formation of moving groups. It is likely that there may be multiple causal mechanisms at play in the formation of moving groups in the Milky Way, including dynamical interactions from perturbing satellites. The analysis of Gaia DR-2 data has revealed many facets of a Galaxy that is clearly out of equilibrium, including the so-called phase-space spiral (Antoja et al., 2018), the
2 P.A.Craig Master’sThesis
Enceladus merger (Helmi et al., 2018), as well as the discovery of a new dwarf galaxy (Torrealba et al., 2019) that likely interacted with the Milky Way (Chakrabarti et al., 2019b). Motivated by these findings, we focus our study here in trying to understand if some of the moving groups in the Milky Way may have arisen from dwarf galaxy interactions. Here, we explore the formation of moving groups in Smoothed Particle Hydro- dynamics (SPH) simulations of the Milky Way interacting with satellites, and contrast this case with isolated simulations to understand the di↵erences in their evolution and structural appearance. By carrying out full N-body simulations of Milky Way like galax- ies that include stars, gas, and a live dark matter halo, which are perturbed by dwarf galaxies, we go beyond prior work that has used either equilibrium models to analyze moving groups, isolated simulations, or test particle calculations. Earlier work (Quillen et al., 2009) used test particle calculations to study the formation of moving groups in an interacting simulation with a perturbing satellite. Hybrid techniques (incorporat- ing both massive and tracer particles) used to simulate the Galactic disk also produce moving groups that resemble the Hercules stream (Quillen et al., 2011). Our focus on analyzing moving groups in the U-V plane is complementary to other recent work, including Khanna et al. (2019b) and Laporte et al. (2019) who have recently carried out analysis of Gaia DR-2 data, with a focus on understanding the structures in the Z V plane. It is important to emphasize that in contrast to works such as Ramos h zi et al. (2018) that have identified a comprehensive list of moving groups in Gaia DR-2, we are not attempting to present an exhaustive list of all the moving groups, but rather are focusing on identifying the large and prominent moving groups in an attempt to understand their formation scenarios. The co-planar dwarf galaxy interaction that we examine here is adopted from our prior work (Chakrabarti and Blitz, 2009), which led to the prediction of a new dark-
3 P.A.Craig Master’sThesis matter dominated dwarf galaxy on an orbit with a close pericenter ( < 10 kpc) to explain the observed large planar disturbances in the outer HI disk of the Milky Way (Levine et al., 2006). This interaction is similar to the one in Chakrabarti et al. (2019b) for the recently discovered Antila2 dwarf galaxy, and Antila2’s current radial location matches our earlier prediction (Chakrabarti and Blitz, 2009). If Antila2 is indeed the dwarf galaxy predicted earlier, some of the dynamically created velocity sub-structures in the solar neighborhood may be the result of this interaction. The simulations we analyze here are similar to ones that successfully reproduce the location and mass of a dwarf galaxy based on analysis of disturbances in the outer HI disks of galaxies (Chakrabarti et al., 2011). We also analyze a simulation that includes the three main tidal players of the Milky Way – the Large Magellanic Cloud (LMC), the Small Magellanic Cloud (SMC), as well as the Sagittarius (Sgr) dwarf galaxy. In all the interacting simulations, moving group structures form in the disk of the galaxy shortly after the dwarf galaxy begins to perturb the disk, and as we discuss below, their distribution and number di↵er from those in isolated simulations.
2.2 Simulations
We analyze three Milky Way-like simulations, two that include dwarf galaxy interactions and one without a satellite. One simulation includes the LMC and SMC, as well as the Sagittarius (Sgr) dwarf galaxy. The initial conditions and orbits for these three dwarf galaxies have been derived using HST proper motions (Lipnicky and Chakrabarti, 2017). Starting with the observed positions and velocities at present day, we carry out a backward integration in a galactic potential (the Hernquist profile matched to NFW as described in (Springel et al., 2005)) using an orbit integration code (Chang and Chakrabarti, 2011). We then evolve the simulation forwards with these initial conditions for timescales of several Gyr. We have listed our derived initial conditions (the starting
4 P.A.Craig Master’sThesis
Object x(kpc) y(kpc) z(kpc) vx (km/s) vy (km/s) vz (km/s) Sagittarius 29.6 0.58 -15.8 -66.6 30.1 164.6 LMC 48.8 259.2 -67.0 -31.7 -235.4 -2.80 SMC -10.7 184.2 4.6 40.3 -136.4 81.7 Coplanar 97.5 21.2 1.2 -251.8 -23.8 0.2
Table 1: Initial positions and velocities for the satellites in our simulations. positions and velocities) for the satellites in Table 1, and we describe our simulation setup below. In the simulation that includes the Sgr dwarf, LMC, and SMC, the Sgr dwarf experiences two pericenters over the last three Gyr, with the most recent pericenter occurring 0.3 Gyr before present day. The LMC and SMC do not approach closer than 50 kpc relative to the Galactic disk over the last few Gyr, as found also in ⇠ earlier work (Besla et al., 2012). The progenitor mass of the Sagittarius dwarf in this simulation is 2.89 1010M , and the progenitor masses for the LMC and SMC are ⇥ 1.90 1010M and 0.19 1010M respectively. The co-planar simulation that we analyze ⇥ ⇥ here is the same one that is described in Chakrabarti & Blitz (2009), i.e., it includes a
10 massive dwarf (Mdwarf 10 M )progenitorwithanorbitalplanethatisnearlyinthe ⇠ disk (Chakrabarti and Blitz, 2009), and a close pericenter approach ( 8kpc).This ⇠ close pericenter approach occurs 0.3Gyrbeforepresentday.Thethirdsimulationdoes not have a perturber, and we use this ”null” simulation as a point of contrast to the interacting galaxy simulations. The co-planar simulation is evolved forward in time for a total of 2 Gyr, the Sagittarius simulation for 5.8Gyr,andthenullsimulationfor3Gyr.Weanalyzethe simulation output files at intervals of 50 Myr. The parameters of the co-planar simulation performed by Chakrabarti & Blitz (2009) are similar to (but are not exactly the same as) those done recently by Chakrabarti et al. (2019), which describe a dynamical model for the newly discovered Antlia2 dwarf galaxy.
5 P.A.Craig Master’sThesis
As in prior work (Chakrabarti and Blitz, 2009), we use a Hernquist profile matched to the NFW profile (Springel et al., 2005) to simulate the dark matter halo of the Milky-Way like galaxy. We use a concentration of 9.4, a spin parameter of =0.036
1 and a circular velocity V200 =160kms .Wealsoincludeanexponentialdiskofgas and stars, as well as an extended HI disk. The disk mass fraction is 4.6% of the total mass, which results in a radial scale length for the exponential disk of 4.1 kpc. The mass fraction of the extended HI disk relative to the total gas mass is 0.3, with a scale length that is three times that of the exponential disk Chakrabarti and Blitz (2009). We have utilized the parallel TreeSPH code GADGET-2 (Springel, 2005) to create simulations of a Milky Way-like galaxy interacting with satellites. GADGET uses an N-body method to simulate the collisionless components (the dark matter and stars) and Smoothed Particle Hydrodynamics (SPH) to follow the evolution of the gaseous component. The gravitational softening lengths are 100 pc for the gas and stars and 200 pc for the dark matter halo particles. In the primary galaxy, we use 4 105 gas and ⇥ stellar particles and 1.2 106 dark matter halo particles. The model is initialized with ⇥ astellarmassof3.5 1010M ,agasmassof1.25 1010M ,andadarkmattermassof ⇥ ⇥ 1.44 1012M .Thesegalacticcomponentsarerepresentedinthesimulationbydi↵erent ⇥ particle types, which are gas, disk, halo and new star particles. Disk particles represent the stellar component of the galaxy which is present at the start of the simulation; these stars do not evolve over time. As the simulation progresses, the simulation gas can be converted into new stars following a Kennicutt-Schmidt prescription for star formation. Our version of GADGET-2 uses an e↵ective equation of state (Springel and Hern- quist, 2003) where a subresolution model is adopted for energy injection from supernovae. Star formation is modeled using the Kennicutt-Schmidt algorithm. We adopt a value of 0.75 for the bulk artificial viscosity parameter and the details of the implementation of the viscosity is described in Springel (2005).
6 P.A.Craig Master’sThesis
Face-on renderings of the simulated Milky Way at present day for the gas and stellar density distribution are shown in Figure 1 for the co-planar simulation and in Figure 2 for the simulation with the LMC, SMC, and the Sgr dwarf. The Fourier amplitudes for the spiral structure in this simulation can be found in Chakrabarti and Blitz (2009). Typical values for the Fourier amplitudes of the gas surface density range of the low-order modes in the outer HI disk are 0.2, with comparable values for the ⇠ stellar disk. It is likely that over the time-scales explored in this paper, where moving group formation correlates closely with the interaction with the satellite, the role of gas dynamics (i.e., dissipation and angular momentum transport to the stars) is not significant.
Figure 1: (Left) Gas density map at present day in the co-planar simulation. (Right) Stellar density map at present day in the co-planar simulation
Figure 2: (Left) Gas density map at present day in the simulation with the three satellites (LMC, SMC, and a Sgr interaction). (Right) Stellar density map at present day.
7 P.A.Craig Master’sThesis
2.3 Detection Algorithm
We developed a method to find and characterize groups in various data sets in an automated fashion given the large quantity of data. This method is designed to be aquickwaytoidentifythelargestvelocitystructuresinadataset.Itisnotmeantto detect every group that exists, and in fact is likely to miss the smallest structures. It will however reliably detect large moving groups, which is what we require for this work. Studying smaller structures is di cult without running very high resolution simulations, which we do not use here. Given that our goal is to carry out a global analysis of velocity structures in SPH simulations to see if prominent moving groups are formed by dwarf galaxy interactions, we do not require a more finely-tuned algorithm for our purposes. Several other methods have been used for the analysis of velocity structures, most notably a wavelet analysis, as in Ramos et al. (2018). These methods are quite e↵ective in large data sets like the GAIA RVs sample (Gaia Collaboration, 2016; Gaia Collaboration et al., 2018; Soubiran et al., 2018; Katz et al., 2019). Our method is distinctly di↵erent from this method, and is less well suited to analyze structure in the MW since it would miss some of the smaller groups. For our purposes, we simply require an algorithm that can detect the largest velocity structures in our simulations. For the smallest structures present within the Milky Way data it is likely the case that we would require an appreciably higher resolution to detect similar structures in our simulations. For each identified group, the average velocity, velocity dispersion, and stellar fraction is determined. We have applied the same algorithm on data for observed stars in the Solar neighborhood, and for stars selected from regions of each simulation that are similar to the Solar neighborhood in distance from the Galactic center. We fit a series of Gaussian distributions to the (U, V )velocities,whereUpoints towards the center of the Milky Way and V points in the direction of motion of a typical star on a circular orbit, and W points out of the plane. U, V and W are defined as
8 P.A.Craig Master’sThesis
U = VR, V = V Vcirc, ,andW = Vz,whereVR,V , and Vz are Galactocentric cylindrical coordinates. Throughout this work we have used the values for the solar
1 motion given in (Sch¨onrich, 2012). At the solar radius, Vcirc, 218 km s ,asgiven ⇡ in (Bovy et al., 2012). The algorithm begins by binning the 2 dimensional data into a 2D histogram. Then we fit one large, axis-aligned Gaussian, with a center and width determined by the average and standard deviation of the data set in each direction. We then can search for regions in the U V plane that have an appreciably higher density than this one Gaussian profile. This is done simply by evaluating the Gaussian profile in each bin, and then subtracting that amount from the observed star count. Bins with large residuals are most likely part of a moving group. We then fit additional Gaussian distributions, usually containing fewer stars than the original Gaussian, and with smaller standard deviations. Each of these Gaussian distributions describes one of the groups within the data set. Each of the Gaussian distributions added here has a total of five parameters: a center location in the U and V directions, standard deviations in these directions, and one parameter for the amplitude of the Gaussian. The parameters for these distributions are found by removing the background stars and then looking for over-densities in the remaining distribution. The center parameters are simply set by finding the center of an over-density. This allows us to fit a distribution on top of each group in the data set. This usually produces some false-positives due to over-densities caused by random noise. The algorithm uses a 2d Kolmogorov-Smirnov test implemented in python to evaluate the goodness of fit (Peacock, 1983; Fasano and Franceschini, 1987). Most of the parameters are reasonably close to the correct values already, however the standard deviations are sometimes inaccurate, which is remedied by the optimization routine. We can also determine the fraction by subtracting out our final distribution for this group until there is no longer an over-density at that location. These optimizations are done
9 P.A.Craig Master’sThesis for the groups one at a time, where the parameters for all of the other distributions are held constant. Ideally we would be able to fit them all at once, but this is very computationally expensive. The exception is that if two groups are close enough to e↵ect each other, then we optimize them together. This is usually not the case for any of the groups, but this condition is triggered if the two groups have estimated centers that are within 10 kilometers per second. At this point the algorithm will have identified the locations of every large group in the data set, but it also will have false positives, which need to be removed. We do this is by comparing the number of stars in the central region of the group to the number we expect in that region from the background. If there are considerably more stars present than what we expect based on the background, we determine that it is a group. The background is the number of stars in the first, large Gaussian that was fit to the data. The error in the back ground is the square root of the number of stars or bodies in that pixel. The residual is the number of stars or bodies in the pixel, in excess of the background counts. A residual value more than twice this standard deviation means the pixel deviates significantly from the model without the new group. The 2 requirement is only relevant in the higher velocity density region of a moving group. In the velocity wings of the background distribution, the primary criterion is the requirement on the minimum number of stars, since out there the expected error on the background model is small. In this region, the 2 requirement can be satisfied by a single star, but we enforce a minimum number of stars in order to classify a velocity structure as a group. If any pixel in a given group satisfies the 2 criterion and the minimum object count requirement, then we accept the group as real. Thus, even though the criterion is applied to every pixel, not every pixel in an extended structure must satisfy this requirement. The result tends to be that in groups spanning a wide range in U and V, only one of the central pixels of the group must satisfy our criteria, and pixels farther
10 P.A.Craig Master’sThesis out in the velocity tails of the group do not have to satisfy any of these criteria. This is done because it is unlikely that the remaining overdensities outside of the central pixels are large enough to satisfy any such criteria. This was selected as it is strong enough to eliminate many of the false positives, while still identifying known structures in the Hipparcos data set. This eliminates all of the cases where the algorithm has identified a velocity structure that is most likely random noise. This algorithm also requires that any group have at least 3 particles, which prevents any one particle with a high velocity from being classified as a group. A potential di culty faced by this algorithm is that there is a dependence on the number of stars, which can make the comparison of the Milky Way and the simulations di cult, as the number of particles in the simulations within the solar neighborhood is less than the number of stars in the solar neighborhood with data in Gaia DR-2. However, as we discuss below, we find in practice that this is only a problem for very small moving groups. One potential concern with a method like ours arises from asymmetric drift. It is well known that this makes a significant modification to the local velocity distribution (Golubov et al., 2013). To test whether the e↵ects in the simulations are su ciently small that it does not influence the results of our algorithm, we have used a Kolmogorov- Smirnov (K-S) test. Specifically, we check that in regions of the simulation that do not contain any large velocity structures (i.e. that are reflective of the background), atwosidedK-Stestyieldsp-valuesthattypicallyrangefrom0.2to0.65.Thus,we cannot reject the hypothesis that the background is symmetric. Asymmetric drift will be important in a search for small groups with this method, especially in the Milky Way, but any group of the size in which we are interested is still easily detectable. In the simulations the e↵ects of asymmetric drift are relatively small. The e↵ect is still present, but it is not significant enough to make our background distribution a poor fit or to impact our ability to detect large scale velocity structures.
11 P.A.Craig Master’sThesis
2.4 Discussion and Results
Figure 3 shows the U and V velocities for the disk stars in our co-planar simulation. The velocity dispersion in the z direction is much lower than in the others, which makes the velocity substructures relatively uninteresting. The stars included in this figure are only the stars that have existed since the beginning of the simulation, i.e., none of the newly generated stars are included. The data here is selected from a spherical region with a radius of 1 kpc centered at as location about where the Sun would be in the Milky Way. Figures 4 and 5 show similar plots for the remaining two simulations, with Figure 4 showing the U and V velocities for the Sagittarius simulation. We note that we recover similar structures, (but not as detailed due to our lower resolution) in the Z V plane h zi as Laporte et al. (2018) and Khanna et al. (2019a). The presence of the ”vertical waves” in the stellar disk discussed in these and earlier papers (Xu et al., 2015) manifests itself in the Z V representation, but not in the U-W or V-W representations. h zi Figures 3 , 4 and 5 show the distribution of disk particles in the U-V plane. The panels are labeled with the time from present day, with negative times occurring before the present day in the simulations. Blue circles mark the central locations of moving groups detected in the simulations. Early in the simulations there has not been su cient time for structures to form, so the earliest panels do not show any interesting velocity structures. At later times, we can clearly see interesting velocity structures form in both of the interacting simulations. There is much less prominent velocity structure formation in the null simulation. The lack of significant velocity structures in the null simulation while being present in the interacting simulations is consistent with the idea that interactions cause the formation of the moving groups that we see here. In our analysis here, we employ the old stars to identify moving groups in our simulations. We do not see significant di↵erences in moving group formation when we include the new stars in addition to the old stars in our simulations. Carrying out this
12 P.A.Craig Master’sThesis
100 -0.55 Gyr 0.0 Gyr 9 10 8 7 50 8
1 1 6 5 0 6
km s km s 4 4
V V 3 50 2 2 1
100 1 0 1 0 U km s U km s 100 0.1 Gyr 9 0.9 Gyr 8 8 7 50 7 6
1 6 1