SciELO - Scientific Electronic Library Online

 
vol.13 issue24Anticorrosive Properties of Chromium Coatings on AISI H13 Steel Produced via Gaseous Nitriding in a VacuumTechnical and Economic Evaluation of a Small-Scale Wind Power System Located in Berlin, Colombia. author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

  • On index processCited by Google
  • Have no similar articlesSimilars in SciELO
  • On index processSimilars in Google

Share


Tecciencia

Print version ISSN 1909-3667

Tecciencia vol.13 no.24 Bogotá Jan./June 2018

https://doi.org/10.18180/tecciencia.2018.24.6 

Artículos

Bifurcation in the Tidal Streams of Sagittarius: Numerical Simulations

Bifurcación en las Corrientes de Marea de Sagitario: Simulaciones Numéricas

Yeimy Camargo Camargo1  2  * 

Rigoberto A. Casas-Miranda2 

1 Universidad de Cundinamarca, Facatativá, Colombia

2 Universidad Nacional de Colombia, Bogotá, Colombia


Abstract

N-body simulations of the interaction between the Sagittarius dwarf galaxy and the Milky Way are performed to study whether it is possible to reproduce both the bifurcations in the tidal streams and the physical properties of Sagittarius. Both galaxies are simulated as sets of particles (live potentials). The initial structures are formed using GALIC and the gravitational interaction between these galaxies is simulated with GADGET-2. Sagittarius initially is simulated as a rotating disk-type galaxy with two components: dark matter halo and stellar disk and the Milky Way is a disk galaxy with three components: dark matter halo, stellar disk and central spheroid. Several numerical simulations of the interaction of this galaxy with the Milky Way with angles between the orbital angular momentum and the disc angular momentum in the range 180° ≤θ ≤ 0° were performed, showing that, although the position and physical properties of Sagittarius are not reproduced entirely by the remnant, the bifurcations appear clearly when 180° ≤θ < 90°.

Keywords: Galaxy: evolution; Galaxy: formation; Galaxies: individual (Sagittarius dSph)

Resumen

Se realizan simulaciones de N-cuerpos de la interacción entre la galaxia enana Sagitario y la Vía Láctea para estudiar si es posible reproducir tanto las bifurcaciones en las corrientes de marea como las propiedades físicas de Sagitario. Ambas galaxias son simuladas como conjuntos de partículas (potencial vivo). Las estructuras iniciales son conformadas en GALIC y la interacción gravitacional entre estas galaxias es simulada con GADGET-2. Sagitario inicialmente es una galaxia de tipo disco rotante con dos componentes: Halo de materia oscura y disco estelar y la Vía Láctea como es usual, es una galaxia de disco con tres componentes: halo de materia oscura, disco estelar y esferoide. Se realizan varias simulaciones de la interacción del progenitor de esta galaxia enana con la Vía Láctea, donde el ángulo entre el momento angular orbital y el momento angular del disco es 180° ≤ θ ≤ 0° mostrando que aunque no se reproducen en su totalidad la posición y propiedades físicas del remanente de Sagitario, las bifurcaciones aparecen de modo evidente cuando 180° ≤θ ≤ 90°.

Palabras clave: Galaxia: evolución; Galaxia: formación; Galaxias: individual (enana Sagitario)

1. Introduction

The study of the satellite galaxies of the Milky Way has been of great interest because they may help answer two fundamental questions about the evolution of galaxies [1]: the actual content of dark matter inside these dwarf galaxies and their implications in the hierarchical formation of the galactic halo.

One of the most interesting members of the Local Group is the Sagittarius dwarf galaxy, discovered in 1995 when R. Ibata et al. [2] did a study on the spheroid of the Milky Way. It is the closest known satellite galaxy to the Milky Way, subtends an angle of ~10° In the celestial vault and is at a galacto-centric distance of 24 ± 2 kpc [2] [3] [4].

Because of this closeness, the Sagittarius progenitor has been exposed to enormous tidal forces. Consequently, Sagittarius has been and continues to be dismembered, leaving as evidence of this interaction a trail of great tidal tails. The study of these tidal tails can reveal the shape of the gravitational potential of the Milky Way and the orbit of Sagittarius. In general, it can be stated that the study of these tidal tails of the satellite galaxies in the Local Group can give important clues about the formation scenarios of galaxies such as the Milky Way [5].

Traditionally it is assumed that the Sagittarius progenitor belongs to the spheroidal type (dSph) [6] [7] [8] [9] [10][11] [12] [13], that is, there are research works trying to reproduce the characteristics observed in the tidal currents of Sagittarius (e.g. [14] [15] [16] [17] [18]), where the progenitor of the dwarf galaxy is a spheroidal galaxy (a distribution of stars and a halo of dark matter), without a subjacent theoretical model that fully explains the data provided by the 2MASS [16] and the SDSS [19] galaxy surveys. These two catalogs have allowed studying in great detail the tidal tails of Sagittarius, showing that these tails have bifurcations, besides constituting a percentage of the galactic halo of the Milky Way. Evidence of a bifurcation in the northern tidal stream of Sagittarius consisting of two branches observed from data from the SDSS DR5 [19].

Here they clearly show how the bifurcation begins around the Galactic North Pole with RA≈180°. The tidal tails are located in (l, b, R) = (350°, 50°, 46 kpc) to the north and (157°, -58°, 33 kpc) to the south [19] [20] [21], where (l, b, R) correspond to the galactic coordinates: galactic longitude, galactic latitude and distance to the Sun, respectively. Branch A is in the lowest declination while branch B is in the highest declination (See Figure 1 in [19]). In [9] it is shown that the trace of the north current of Sagittarius is between RA = 270° and RA = 190°, where the separation distance of branches A and B is ≈ 15 kpc and that there is no evidence that any part of this tidal tail passes through the solar neighborhood. A prominent Sagittarius debris arm extends from the center of this galaxy to the north galactic plane at a distance of ~45 kpc from the Sun and turns toward the top of the galactic north (NGC), from there descends back to the galactic plane.

The distribution of the debris of the Sagittarius tidal tails traces their orbit [9]. The stellar tidal tails present in the halo of the Milky Way produced by accretion of smaller galaxies are a prediction of the hierarchical model of structure formation in the Universe. It has been mentioned that branches A and B were probably formed in different stages according to [18] [19] [22] [23]. In 2006 Fellhauer et al. [23] performed several numerical N-Body simulations of the evolution of a disc-like progenitor of Sagittarius orbiting the Milky Way, which was modeled with a rigid potential.

In their simulations the progenitor of Sagittarius was a Plummer sphere. They were able to reproduce the bifurcations in the tidal streams of Sagittarius only if the potential of the Milky Way is near to spherical. It has been found that these branches have similar distances, velocity distributions and metallicities [24] [25], which would refute the hypothesis that the bifurcations are of different ages or that they are tidal streams of different progenitors [26]. There is no definitive evidence yet about that.

In 2010 Peñarubia et al. [18] considered the possibility that the progenitor of Sagittarius had not been a dSph but a rotating disk galaxy. They sought to examine how tidal currents are affected by rotation without trying to reproduce the observed characteristics of neither the tidal stream nor the current state of the galaxy. For that they used N-body simulations with a rigid potential to model the Milky Way, and found bifurcations in the simulated tidal stream similar to those observed. Table 1 shows the observed physical properties for this satellite galaxy, in addition, it is estimated that the current mass of Sagittarius is Msgr =2-5x108 M⊙ [3].

Based on the work done by [18] and [23], in this paper, we study the possibility that the progenitor of Sagittarius was a rotating disk galaxy and try to reproduce the bifurcations in its tidal tails and the observed physical properties (present) of Sagittarius modeling the Milky Way as a live potential rather than with a rigid potential, as done by [18] and [23]. Some preliminary parts of the research work presented in this paper have been published in proceedings [27]. An extended and detailed documentation about this research can be found in the master's thesis (in Spanish) [28].

This work is organized as follows: in Section 2 we show the structural parameters for each galaxy and the initial orbital conditions in the simulations. In Section 3 we present the results from the simulations, comparing the observed values of position, bifurcations and physical properties found in the cases of interest.

2. Structures and initial conditions

In this work, the Milky Way is modeled as a disk-type galaxy with three components: dark matter halo, stellar disc and stellar spheroid. Sagittarius is considered a disk-type galaxy with two components: a halo of dark matter and a stellar disk. The initial conditions of the possible orbits of Sagittarius are determined by simulating the orbit of a test particle in a rigid potential representing the Milky Way.

2.1 Initial Structures

To generate the initial galactic structures for the Milky Way and Sagittarius, the GALIC software [32] has been used.

Table 1 Physical properties of Sagittarius. References are listed to the right of each value. 

GALIC is an open source code for the construction of N-body galaxy models using an iterative method, starting with some initial particle's velocities from an initial mass density distribution. The difference of the time-averaged density response produced by the particle orbits with respect to the initial density configuration is characterized through a merit function, and a stationary solution of the collision-less Boltzmann equation is found by minimizing this merit function directly by iteratively adjusting the initial velocities. The velocity adjustment is done with a stochastic process from an approximate solution of the distribution function, seeking for an improvement of the fit.

The GALIC program requires an input file with the values of the parameters of the galaxy to be created. The values used to create the Milky Way and Sagittarius with GALIC are listed in Tables 2 and 3 respectively.

2.1.1 The Milky Way

For the modeling of the Milky Way galaxy, the parameters are shown Table. Here a h is the radial scale length of the dark matter halo, C is the concentration parameter, ad is the radial scale length of the disk, z d is the vertical scale length of the disk, N h is the number of particles in the dark matter halo, N d is the number of particles in the disc and N b is the number of particles in the bulge. In Figure 1 the initial structure of the Milky Way, projected on the planes xy, yz, xz respectively, is shown.

2.1.2 Progenitor of Sagittarius

For the model of the progenitor of the Sagittarius dwarf galaxy the parameters are shown in Table 3. Figure 2 shows the initial structure of the progenitor of Sagittarius, projected on the planes xy, yz, xz.

Table 2 Parameters used to model the Milky Way. 

Table 3 Parameters used to model the progenitor of Sagittarius [5]

Table 4 Softening length for each of the galactic structures used in the simulations [5]

The numerical structures obtained with GALIC [32] were simulated under their own gravity with GADGET-2 [33] over a period of time of 2.0 Gyr to ensure virial balance, where the values for the softening length (ε) are shown in Table 4, where εh, εd, εb, are the softening lengths for the particles in the halo, disc and bulge, respectively.

GADGET-2 is an open source code to perform N-body simulations including gas dynamics on massively parallel computers with distributed memory. The gravitational forces are computed with a hierarchical tree algorithm in combination with a particle-mesh scheme for long-range gravitational forces (optional), while the gas dynamics are simulated using smoothed particle hydrodynamics (SPH). GADGET-2 can be used to perform either Newtonian or cosmological simulations over a wide range of scales.

Figure 1 Initial structure of the Milky Way. Projections on the planes xy, yz and xz, respectively, are shown. The colors correspond to the different components of the Galaxy: blue for the bulge, green for the disk and red for the dark matter halo. 

Figure 2 Initial structure of the Sagittarius dwarf. Projections on the planes xy, yz and xz, respectively, are shown. The colors correspond to the different components of the simulated galaxy: green for the disk and red for the dark matter halo. 

Figure 3 Schematic representation of the arrangement of the orbit and the orientation and rotation direction of the dwarf galaxy. 

2.2 Initial Conditions: Orbit

The orbit of a test particle in the rigid potential of the Milky Way was simulated back in time, with an initial radial velocity of 171 km/s and a tangential velocity varying between -325 km/s and 325 km/s, at intervals of 5.0 km/s. For each simulated orbit the pericentric and apocentric distance are calculated taking into account that the observational values for the apocenter are between 50 and 70 kpc and for the pericenter are between 10 and 15 kpc. Initial conditions for the orbit are shown in the table 5.

2.3 Numerical Simulations

For the set of initial conditions in the Table 5 the evolution of the Sagittarius progenitor is simulated using GADGET-2 during 7 Gyr with several angles of inclination θ between the angular momentum of the disc and the angular momentum of the orbit in the range 0° ≤ θ ≤ 180° (see Figure 3). For the simulations the softening lengths shown in Table 4 were used. We constructed a program to set the initial conditions of the orbit of the progenitor and the inclination of the disk respect to the plane of the orbit (angle θ) in the format required by GADGET-2.

The simulations are performed in Cartesian coordinates with the center of the Milky Way at the origin. The positions and velocities of all particles in the simulation are stored in a snapshot file each 0.1 Gyr.

The physical properties of the remnant of Sagittarius to be determined are the central surface brightness, the mass-to-light ratio, the line-of-sight central velocity dispersion and the half-light radius.

To carry out the theoretical calculation of these quantities we perform the following sequence: first we tag the most bound particles of the progenitor at the initial time of the simulation and use them iteratively to compute the position of the center of mass of the object at the following snapshot times, taking care of those particles that are not inside the core anymore.

After this we projected the stars (disc) particles on the sky, as observed from the position of the Sun, that lie inside a circle of a given radius around the center of mass of the remnant; and computed the corresponding quantities following the method described in [34] [35] [36] [37].

Figure 4 Tidal tails of the dwarf galaxy, projected on the xy plane, as it orbits the Milky Way. The horizontal line shows the position of the disk of the Milky Way. The red dots correspond to the disk particles of the progenitor of Sagittarius. 

In order to search for the bifurcations in the tidal tails of Sagittarius we projected the positions of the stars particles (disc) of the dwarf galaxy on the sky using equatorial coordinates in declination-right ascension maps. This was done for each snapshot in each of the simulations we have performed.

3. Results

Figure 5 shows the tidal tails obtained for different angles between the angular momentum of the disc of the dwarf galaxy and the angular momentum of the orbit. It can be observed that the bifurcations are more evident when θ = 130° and θ = 160°.

In addition, since for θ= 130° the bifurcation is more similar to the features observed by the SDSS DR5, we compute the physical properties of the remnant of Sagittarius in this simulation. When θ= 0° and θ= 90° no bifurcation is generated. Our results about the bifurcation for θ= 160° are similar with those obtained by [18] (θ = -20° in their paper).

Along with the loss of particles of the dwarf galaxy due to the tidal forces produced during its interaction with the Milky Way and the deformation of the initial stellar disk, the dwarf suffers a succession of quenching and stretching events.

Table 6 Physical properties of the remnant of Sagittarius for θ = 130°. 

Table 7 Physical properties of the remnant of Sagittarius with θ = 160º 

All these processes combined produce a change in the properties an observer would measure as the remnant of the disk changes its shape and stellar content. Moreover, the observed properties could change with the relative orientation of the remnant respect to the position of the Sun.

Figure 6 shows the time evolution during t = 7.0 Gyr of the central surface brightness, the half-light radius, the central velocity dispersion along the line of sight and the mass-to-light ratio for the progenitor of Sagittarius with an angle θ = 130°. These quantities were computed as described in [36].

Figure 6 Time evolution of the physical properties of Sagittarius for θ = 130°. The horizontal lines show the corresponding value(s) reported in the literature for each of the properties of Sagittarius (see Table 3). Error bars are of the size of the symbols in the plots. 

As we already mentioned, the simulation of interest will be θ = 130° due to the great similarity between the bifurcations found in our simulations and the observational data. For this case, the bifurcation arises in t = 5.0 Gyr and it starts in RA ~170°, the separation between branches is ≈ 15 kpc.

Table 6 shows the results obtained for the physical properties of the remnant of Sagittarius. Recall that the current position of Sagittarius is (RA, dec, R) = (238.7°, -30.5°, 24 kpc). For this simulation although the bifurcation is found it is not possible to locate the remnant in a position close to the present one.

Some physical properties like the mass-to-light ratio and the of central velocity dispersion have a value close to the observed one, while quantities such as the half-light radius and the central surface brightness have values significantly different from those observed. The remnant is closer to the current position when t = 6.6 Gyr.

However, that position is significantly away from the current position and the physical properties of the remnant at that position are quite far from the observational ones.

On the other hand, in the simulation with θ = 160° it is found that the bifurcation appears in t = 5.0 Gyr and it starts at RA 160° and the separation between branches is ≈ 15kpc. Table 7, shows the results obtained for the physical properties of the remnant of Sagittarius when the bifurcation appears and when the remnant is at the position that is closer to the current one, which is t = 6.2 Gyr.

Table 7 shows that, as in the previous case, when the bifurcation appears (t = 5.0 Gyr) the remnant is not close to the actual position and the physical properties are not similar to those observed for this galaxy. Nevertheless, at t = 6.2 Gyr the position of this remnant is close to the observed position and quantities such as the half-light radius and the central velocity dispersion have values close to the current ones.

4. Conclusions

In this study we find the bifurcations in the tidal streams of Sagittarius similar to those observed by SDSS DR5 when θ = 130° and θ = 160° at t = 5.0 Gyr, However, when they appear, the main body of this galaxy is not located at the actual position of Sagittarius; for θ = 130°, close values are found in the mass-to-light ratio and the central velocity dispersion of the main body when the bifurcation becomes evident.

The position of the remnant closest to the observed one occurs when θ = 160° and t = 6.2 Gyr and values close to those observed for the half-light radius and the central velocity dispersion are found, i.e., although the bifurcations and the position of the remnant close to the present one exist, in no case they reproduce simultaneously the position of the remnant and the bifurcations nor is there a set of initial conditions that reproduce all the observational physical properties of this galaxy.

Our results regarding the reproduction of the bifurcation in the tidal streams of Sagittarius are similar to those obtained by [18 and [23], but modeling the Milky Way as a live potential. Thus, confirming their findings. It appears that under the right parameters it is possible to reproduce the bifurcation in the tidal streams of Sagittarius both with a spherical and with a rotating disc galaxy. The problem is still open and further theoretical and observational studies are required. A possible extension of this work would be to perform numerical simulations with other stellar to dark matter ratios in the progenitor and changing the orbital parameters.

References

[1] M. A. Gómez-Flechoso, R. Fux, and L. Martinet, “Sagittarius, a dwarf spheroidal galaxy without dark matter?,” A&A, vol. 347, pp. 77-91, Jul. 1999. [ Links ]

[2] R. A. Ibata, G. Gilmore, and M. J. Irwin, “Sagittarius: the nearest dwarf galaxy,” Mon. Not. R. Astron. Soc., vol. 277, no. 3, pp. 781- 800, Dec. 1995. [ Links ]

[3] D. R. Law, K. V. Johnston, and S. R. Majewski, “A Two Micron All‐ Sky Survey View of the Sagittarius Dwarf Galaxy. IV. Modeling the Sagittarius Tidal Tails,” Astrophys. J., vol. 619, no. 2, pp. 807-823, Feb. 2005. [ Links ]

[4] K. V. Johnston, D. N. Spergel, and L. Hernquist, “The Disruption of the Sagittarius Dwarf Galaxy,” Astrophys. J., vol. 451, p. 598, Oct. 1995. [ Links ]

[5] E. L. Łokas, G. Gajda, and S. Kazantzidis, “Tidal tails of dwarf galaxies on different orbits around the Milky Way,” Mon. Not. R. Astron. Soc., vol. 433, no. 1, pp. 878-888, Jul. 2013. [ Links ]

[6] W. ~B. Burton and F. J. Lockman, “HI observations towards the Sagittarius dwarf spheroidal galaxy,” A&A , vol. 349, pp. 7-10, Sep. 1999. [ Links ]

[7] R. A. Ibata, R. F. G. Wyse, G. Gilmore, M. J. Irwin, and N. B. Suntzeff, “The Kinematics, Orbit, and Survival of the Sagittarius Dwarf Spheroidal Galaxy,” Astron. J., vol. 113, p. 634, Feb. 1997. [ Links ]

[8] T. Goerdt, B. Moore, S. Kazantzidis, T. Kaufmann, A. V. Macciò, and J. Stadel, “The formation of ultra-compact dwarf galaxies and nucleated dwarf galaxies,” Mon. Not. R. Astron. Soc., vol. 385, no. 4, pp. 2136-2142, Apr. 2008. [ Links ]

[9] S. R. Majewski, M. F. Skrutskie, M. D. Weinberg, and J. C. Ostheimer, “A Two Micron All Sky Survey View of the Sagittarius Dwarf Galaxy. I. Morphology of the Sagittarius Core and Tidal Arms,” Astrophys. J., vol. 599, no. 2, pp. 1082-1115, Dec. 2003. [ Links ]

[10] B. Koribalski, S. Johnston, and R. Otrupcek, “H I observations of the new Sagittarius dwarf,” Mon. Not. R. Astron. Soc., vol. 270, no. 1, pp. L43-L45, Sep. 1994. [ Links ]

[11] M. Mateo, A. Udalski, M. Szymanski, J. Kaluzny, M. Kubiak, and W. Krzeminski, “The optical gravitational lensing experiment: Deep photometry of the Sagittarius dwarf spheroidal galaxy,” Astron. J ., vol. 109, p. 588, Feb. 1995. [ Links ]

[12] A. C. Layden and A. Sarajedini, “Photometry of the Globular Cluster M54 and the Sagittarius Dwarf Galaxy: The Age-Metallicity Relation,” Astron. J ., vol. 119, no. 4, pp. 1760-1792, Apr. 2000. [ Links ]

[13] M. Bellazzini, F. R. Ferraro, and R. Ibata, “Building Up the Globular Cluster System of the Milky Way: The Contribution of the Sagittarius Galaxy,” Astron. J ., vol. 125, no. 1, pp. 188-196, Jan. 2003. [ Links ]

[14] C. A. Martínez-Barbosa and R. A. Casas-Miranda, “Is the Sgr dSph a dark matter dominated system?,” New Astron., vol. 20, no. 1, pp. 7-14, 2013. [ Links ]

[15] A. Helmi, “Velocity Trends in the Debris of Sagittarius and the Shape of the Dark Matter Halo of Our Galaxy,” Astrophys. J., vol. 610, no. 2, pp. L97-L100, Aug. 2004. [ Links ]

[16] K. V. Johnston, D. R. Law, and S. R. Majewski, “A Two Micron All Sky Survey View of the Sagittarius Dwarf Galaxy. III. Constraints on the Flattening of the Galactic Halo,” Astrophys. J., vol. 619, no. 2, pp. 800-806, Feb. 2005. [ Links ]

[17] D. R. Law, S. R. Majewski, and K. V. Johnston, “EVIDENCE FOR A TRIAXIAL MILKY WAY DARK MATTER HALO FROM THE SAGITTARIUS STELLAR TIDAL STREAM,” Astrophys. J., vol. 703, no. 1, pp. L67-L71, Sep. 2009. [ Links ]

[18] J. Peñarrubia et al., “Was the progenitor of the Sagittarius stream a disc galaxy?,” Mon. Not. R. Astron. Soc Lett., . vol. 408, no. 1, pp. L26-L30, Oct. 2010. [ Links ]

[19] V. Belokurov et al., “The Field of Streams: Sagittarius and Its Siblings,” Astrophys. J., vol. 642, no. 2, pp. L137-L140, 2006. [ Links ]

[20] B. Yanny et al., “Identification of A‐ colored Stars and Structure in the Halo of the Milky Way from Sloan Digital Sky Survey Commissioning Data,” Astrophys. J., vol. 540, no. 2, pp. 825-841, Sep. 2000. [ Links ]

[21] R. Ibata, M. Irwin, G. F. Lewis, and A. Stolte, “Galactic Halo Substructure in the Sloan Digital Sky Survey: The Ancient Tidal Stream from the Sagittarius Dwarf Galaxy,” Astrophys. J., vol. 547, no. 2, pp. L133-L136, Feb. 2001. [ Links ]

[22] A. Helmi, “Is the dark halo of our Galaxy spherical?,” Mon. Not. R. Astron. Soc ., vol. 351, no. 2, pp. 643-648, Jun. 2004. [ Links ]

[23] M. Fellhauer et al., “The Origin of the Bifurcation in the Sagittarius Stream,” Astrophys. J., vol. 651, no. 1, pp. 167-173, Nov. 2006. [ Links ]

[24] B. Yanny et al., “TRACING SAGITTARIUS STRUCTURE WITH SDSS AND SEGUE IMAGING AND SPECTROSCOPY,” Astrophys. J., vol. 700, no. 2, pp. 1282-1298, Aug. 2009. [ Links ]

[25] A. K. Vivas et al., “The QUEST RR Lyrae Survey: Confirmation of the Clump at 50 Kiloparsecs and Other Overdensities in the Outer Halo,” Astrophys. J., vol. 554, no. 1, pp. L33-L36, Jun. 2001. [ Links ]

[26] M. Newby et al., “A SPATIAL CHARACTERIZATION OF THE SAGITTARIUS DWARF GALAXY TIDAL TAILS,” Astron. J ., vol. 145, no. 6, p. 163, May 2013. [ Links ]

[27] Y. Camargo and R. Casas-Miranda , “Simulations of the bifurcation of the tidal streams of the Sagittarius dwarf galaxy,” in Revista Mexicana de Astronomia y Astrofisica Conference Series, 2017, vol. 49, p. 113. [ Links ]

[28] Y. Camargo, “Simulaciones Numéricas de N-Cuerpos de la Formación de la Galaxia Enana Sagitario y sus Corrientes de Marea a partir de un Progenitor con forma de disco Estelar,” Universidad Nacional de Colombia, 2016. [ Links ]

[29] M. Mateo, et al., “Discovery of a Tidal Extension of the Sagittarius Dwarf Spheroidal Galaxy,” Astrophys. J., vol. 458, no. 1, Feb. 1996. [ Links ]

[30] M. Mateo, E. W. Olszewski, and H. L. Morrison, “Tracing the Outer Structure of the Sagittarius Dwarf Galaxy:Detections at Angular Distances between 10° and 34°,” Astrophys. J., vol. 508, no. 1, pp. L55-L59, Nov. 1998. [ Links ]

[31] S. E. Koposov, V. Belokurov, and N. W. Evans, “SAGITTARIUS STREAM THREE-DIMENSIONAL KINEMATICS FROM SLOAN DIGITAL SKY SURVEY STRIPE 82,” Astrophys. J., vol. 766, no. 2, p. 79, Mar. 2013. [ Links ]

[32] D. Yurin and V. Springel, “An iterative method for the construction of N-body galaxy models in collisionless equilibrium,” Mon. Not. R. Astron. Soc ., vol. 444, no. 1, pp. 62-79, Oct. 2014. (http://www.h-its.org/tap-software/galic-code/) . [ Links ]

[33] V. Springel, “The cosmological simulation code gadget-2,” Mon. Not. R. Astron. Soc ., vol. 364, no. 4, pp. 1105-1134, Dec. 2005. (https://www.h-its.org/tap-software-en/gadget-code/) [ Links ]

[34] R. S. Klessen and P. Kroupa, “Dwarf Spheroidal Satellite Galaxies without Dark Matter: Results from Two Different Numerical Techniques,” Astrophys. J., vol. 498, no. 1, pp. 143-155, May 1998. [ Links ]

[35] P. Kroupa, “Dwarf spheroidal satellite galaxies without dark matter,” New Astron ., vol. 2, no. 2, pp. 139-164, Jul. 1997. [ Links ]

[36] R. A. Casas, V. Arias, K. Peña Ramírez, and P. Kroupa, “Dwarf spheroidal satellites of the Milky Way from dark matter free tidal dwarf galaxy progenitors: Maps of orbits,” Mon. Not. R. Astron. Soc ., vol. 424, no. 3, pp. 1941-1951, 2012. [ Links ]

[37] Bohórquez Pacheco, O. A., Casas Miranda, R. A., Origin of the Milky Way Disk of Satellites: Collision of Two Disk Galaxies, TECCIENCIA, Vol. 12 No. 21, 33-38, 2016, DOI: http://dx.doi.org/10.18180/tecciencia.2016.21.6Links ]

How to cite:Casas-Miranda, R.A., Camargo, Y., Bifurcation in the Tidal Streams of Sagittarius: Numerical Simulations, TECCIENCIA, Vol. 13 No. 24, 53-62, 2018 DOI: http://dx.doi.org/10.18180/tecciencia.2018.24.6

Received: August 23, 2017; Accepted: February 13, 2018

*Corresponding Author. E-mail: ydcamarsoc@unal.edu.co

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License