SciELO - Scientific Electronic Library Online

 
vol.23 número2Análisis de inyección directa en gran volumen de hidrocarburos aromáticos policíclicos en aguaDisminución de la infección por rotavirus en células MA104 por la unión de proteínas de extractos probióticos a receptores celulares integrina β3 y Hsc70 índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • En proceso de indezaciónCitado por Google
  • No hay articulos similaresSimilares en SciELO
  • En proceso de indezaciónSimilares en Google

Compartir


Universitas Scientiarum

versión impresa ISSN 0122-7483

Univ. Sci. vol.23 no.2 Bogotá mayo/ago. 2018

https://doi.org/….11144/javeriana.sc23-2.mosg 

Artículos

Mass of Spiral Galaxies by Means of a Maximum Disc Model

Masa de galaxias espirales por medio de un modelo de disco máximo

Massa de galáxias espirais por meio de um modelo de disco máximo

Jerson I. Reina1  2 

Framsol López-Suspes1  2 

Guillermo A. González2  * 

Juan Carlos Salcedo-Reyes, Edited by

1 Departamento de Ciencias Básicas, Universidad Santo Tomás, Bucaramanga, Colombia. jerson.reina@ustabuca.edu.co

2 Escuela de Física, Universidad Industrial de Santander, Bucaramanga, Colombia. framsol.lopez@ustabuca.edu.co


Abstract

Maximum disc mass models for a set of spiral galaxies from the Ursa Major Cluster are presented. We construct the models using the Hunther method and the particular solutions are chosen in such a way that the circular velocities are adjusted very accurately to the observed rotation curves of some specific spiral galaxies. Under the maximum disc hypothesis, we consider that the rotation curves of the analyzed galaxies can be modeled with only the contribution of the disc. This implies that it is not necessary to consider the contribution ot the dark matter halo in the inner part of the spiral. In this way, the models reproduce the global behavior of the rotation curves in the great majority of galaxies. Producing good adjustments to calculate the total mass of these galaxies, and yielding values of the order of 1O10 M . Based on the vertical stability criterion presented by Viera & Ramos-caro (2016), we find that all the galaxies analyzed present a vertically stable behavior. On the other hand, from the analysis of the epicyclic frequency we find that all the models exhibit mainly a radial stable behaviour except at the edge of the disc.

Keywords: Potential Theory; Stellar Dynamics; Galactic Mass

Resumen

Presentamos modelos de masa de disco máximo para un conjunto de galaxias espirales del Cluster Ursa Major. Los modelos se obtienen por medio del método de Hunter y las soluciones particulares se eligen de tal manera que las velocidades circulares se ajustan muy exactamente a las curvas de rotación observadas de algunas galaxias espirales específicas. Bajo la hipótesis del disco máximo, suponemos que la masa del disco es lo más grande posible, en consonancia con la curva de rotación de la galaxia. Por lo tanto, la contribución de la masa del halo de la materia oscura se considera insignificante en las partes internas de las espirales. Los modelos reproducen la estructura general de las curvas de rotación en la mayoría de las galaxias, proporcionando buenos ajustes para calcular la masa total de estas galaxias obteniendo valores del orden de 1010 M . Basados en el criterio de estabilidad vertical presentado por Vieira and Ramos-Caro (2016), encontramos que todas las galaxias analizadas presentan un comportamiento verticalmente estable. Por otro lado, a partir del análisis de la frecuencia epicíclica se observa que todos los modelos presentaron mayormente un comportamiento estable radial excepto en el borde del disco.

Palabras clave: Teoría del potencial; Dinámica Estelar; Masa de Galaxias

Resumo

Apresentamos modelos de massa de disco máximo para um conjunto de galáxias espirais do Cluster Ursa Major. Os modelos são obtidos por meio do método Hunter e as soluções particulares são escolhidas de tal forma que as velocidades circulares são ajustadas com muita precisão às curvas de rotação observadas de algumas galáxias espirais específicas. Sob a hipótese de disco máximo, supomos que a massa do disco é tão grande quanto possível, consistente com a curva de rotação da galáxia. A contribuição de massa do halo da matéria escura é, portanto, assumida como insignificante nas partes internas das espirais. Os modelos reproduzem a estrutura geral das curvas de rotação na maioria das galáxias, proporcionando bons ajustes para calcular a massa total dessas galáxias obtendo valores da ordem de 1010 M . Com base no critério de estabilidade vertical apresentado por Vieira and Ramos-Caro (2016), descobrimos que todas as galáxias analisadas apresentam um comportamento verticalmente estável. Por outro lado, a partir da análise da frequência epicíclica, descobrimos que todos os modelos apresentaram principalmente um comportamento estável radial, exceto na borda do disco.

Palavras-chave: Teoria do Potencial; Dinâmica estelar; Massa Galáctica.

Introduction

One of the greatest challenges of modern astronomy is the determination of the origin, nature and distribution of matter in its visible and dark forms. Regardless of whether the content of matter in a spiral galaxy is luminous or not, the gravity is the only force that affects all forms of matter, whatever its properties and conditions. So that by probing the effects of gravity through its dynamic influence, we can reveal the whole material content of a system. Currently one of the best tools available to analyze the material content of a spiral galaxy is the rotation curve (Binney & Tremaine, 2008; Leblanc, 2010; Rubin et al. 1985), which shows how the stars in a galaxy rotate around the center in approximately circular orbits (Oort, 1927), with a velocity that varies only under the effects of the galaxy’s own gravity.

In order to determine the mass distribution of a particular spiral galaxy, it is necessary to propose a model in such a way that it adjusts with the data of the rotation curve of the galaxy analyzed. Taking into account that each of the components of a spiral galaxy contributes to the rotation curve, there are in the literature several mass models that include one to several components depending on the particular situation to be modeled (Sofue et al. 2009). For example, Bajkova & Bobylev (2016) fitted a three component (bulge, disc, halo) model to the rotation velocities of Galactic objects known from observations in the range of Galactocentric distances R from 0 to 200 kpc. Recently, Gonzales & Reina (2016) presented a family of analytical potentials for flat galaxies with spheroidal halos characterized by superposition of two independent components: a potential associated with the thin disc and a potential associated with spheroidal halo. Kalnajs (1983) and Van Albada et al. (1985), among other authors, showed that the combination of a central bulb surrounded by an exponential disk produces a flat velocity curve with a reasonable stellar mass-to-light ratio.

From the maximum disc hypothesis (Palunas & Williams, 2000), the mass of the galactic disc in a spiral galaxy is accepted to be large enough, so that it is consistent with the rotation curve of the galaxy. In this way, it is not necessary to consider the contribution of the dark matter halo in the inner part to the spiral. Thus, some authors have shown that it is possible to obtain galactic models taking into account only the disc component, without having to include a halo of matter, modeling with high precision rotation curves in the optical region and according to physically reasonable mass density profiles (Kalnajs, 1983; Kent, 1986; Gonzales et al. 2010).

In particular, Gonzales et al. (2010) obtained finite thin disc models for four galaxies in the Ursa Major cluster: NGC3877, NGC3917, NGC3949 and NGC4010. The models were obtained through the Hunter method (Hunter, 1963) and the particular solutions were chosen in such a way that the circular velocities were adjusted very precisely with the data of some rotation curves taken from Verheijen & Sancisi (2001). The obtained models are vertically unstable under the criterion of vertical frequency stability (Binney & Tremaine, 2008). However, recently Viera & Ramos-caro (2016) found that in each thin disc galactic model, all circular equatorial orbits within the matter distribution are vertically stable, which reverses the instability argument of the models of Gonzales et al. (2010).

In order to illustrate the applicability of the procedure proposed by Viera & Ramos-Caro (2016), we are going to apply the methodology proposed in Gonzales et al. (2010), this time analyzing the extension of the model to all galaxies in the sample. Verheijen & Sancisi,(2001). In this way we pretend to remark that the formalism proposed by Gonzales et al. (2010) is a powerful method for obtaining stable maximum discs that describe the optical region of many galaxies.

For this purpose, we initially summarize the most remarkable aspects of the construction of the maximum disc models proposed by Gonzales et al. (2010). We then analyze the formalism developed by Viera & Ramos-caro (2016) and the way in which the analysis of the stability of the equatorial orbits must take into account the discontinuity in the partial derivative of the potential with respect to z. Later, we adjust the model to the data of the observed rotation curves of all the galaxies reported in Verheijen & Sancisi (2001) and, finally, we present a discussion of the obtained results.

Thin Disc Mass Models

As it was developed by Gonzales & Reina (2006), the starting point to construct finite, axially symmetric disc models is to solve the Laplace equation, which in cylindrical coordinates is written as

As we will present in the next section, in order to apply the vertical stability criterion for circular orbits, recently developed by Viera & Ramos-caro (2016), we must consider that the reflection symmetry

implies that the gravitational potential depends only on the absolute value | z | of z. so we can rewrite the normal derivative as

where Θ is the Heaviside step-function. Thus, the potential $ in these models is a smooth function of | z | with a discontinuity in the z derivative at the plane z = 0.

Given a gravitational potential Φ(R, z) that satisfies the above considerations, the circular velocity v c (R) can be obtained from the gravitational potential through the relation

Besides, by means of the Gauss Theorem and the equation (3), we can express the surface distribution of matter as

which must be adjusted to the boundary conditions

in order to obtain a discoidal distribution of matter restricted to the plane z = 0 and with a radius R in the interval 0 < R < a.

As it was developed by Gonzáles & Reina (2006) and Gonzáles et al. (2010), a gravitational potential satisfying the conditions mentioned above and describing a discoidal distribution of axially symmetric matter has the form (Bateman, 1944)

Here P n (η) are the usual Legendre polynomials and q n (ξ) = i n+1 Q n (iξ), being Qn (z) the Legendre functions of the second kind. The spheroidal oblate coordinates , n), with 0 ≤ ξ ∞ and -1 η 1, are related to the usual cylindrical coordinates (R, z) through the relations

in such a way that the disc are located at . Finally the C2n are arbitrary constants, which are chosen properly in order to adjust the circular velocity of the model and the rotation curve of some particular galaxy.

With this gravitational potential, and using the equations (4) and (5) we have that the circular velocity and the surface mass density can be written respectively as

and

where R = R/α. By integrating over the total area of the disc, we find the expression

which we will use to determine the total mass of each galaxy M. In this way, we can see from (7), (9) and (10) that the whole model is completely defined in terms of constants C2n, which will be determined by fitting (9) with the observational data associated with the rotation curve of a particular galaxy.

A Family of Particular Models

In this section we show how from (9) we can adjust the circular velocity of the model with the observed data for a given galaxy in order to construct particular models of thin disks. To this end, we must take into account that the sum must be limited to a finite number of terms. This implies that C2n = 0 for n > m, with m a positive integer, which represents the number of parameters necessary to obtain the best fit. Replacing the derivatives of the Legendre polynomials in (9) we can rewrite the circular velocity of the model as

Equation (9) and (12) and using the orthogonality properties of the associated Legendre functions (Arfken & Weber, 2005), we find that the A2n constants in (12) and the C2n constants are related through the expression

for n ≠ 0, where

This implies that if the constants A2n are determined by the adjustment of the observational data of the rotation curve, then using (13) we can obtain the values of the constants C2n which allows us to define the particular thin disc model by means of (7).

From (13) it is easy to see that the value of the constant C0 cannot be obtained through this equation. Additionally, from (10) we can see that the surface mass density diverges at the edge of the disc when q = 0, unless the following condition is imposed (Hunter, 1963)

After a couple of operations using the properties of Legendre functions, we arrive at the expression

by means of which we can obtain the value of C0, and furthermore, through (11), calculate the total mass in terms of A2n.

Previously in Gonzáles et al. (2010), it was shown that, in a first approximation, the condition for the stability of the discoidal structure, composed mainly by particles that describe nearly circular motions, is that both epicyclic and vertical frequencies are positive, (Binney & Tremaine, 2008; Danby, 1988), i.e

and

where

is the effective potential and l = Rvc is the specific axial angular momentum.

The procedure used to arrive at the analysis of the epicyclic and vertical frequencies around circular orbits (Binney & Tremaine, 2008), is valid only for potentials that are smooth or that have at least continuous second derivatives. This means that the analysis of the vertical frequency of disc potentials, such as those analyzed in Gonzáles et al. (2010), is not a reliable indicator of the vertical stability of the corresponding circular orbit. The analysis must take into account the discontinuity in the partial derivative of the potential with respect to z due to the surface distribution of matter in the equatorial plane Viera & Ramos-caro (2016).

In order to analyze the behavior of the z coordinate of a vertically perturbed orbit in the vicinity of a discoidal distribution, as it was recently pointed out by Viera & Ramos-caro (2016), we consider an equatorial circular orbit of radius R under the action of a small disturbance, which can be seen as an instantaneous vertical increment v0z, in the neighborhood of the particle, which does not affect the value of l. The equation of movement that allows us to study the temporal evolution of the vertical perturbation is

From this equation it can be noted that the vertical disturbance will remain small if and only if . This guarantees that if v0z is small enough, the perturbed path will oscillate around the original one for a long time. The particle will always traverse the disc with a velocity whose vertical component has an absolute value |v0z |. Then, just after crossing the disc, the particle will also have a velocity with vertical component of magnitude |v0z |, which means that the movement after crossing the disc will have the same behavior as the movement before crossing it. Therefore, the disturbed orbit will remain oscillating around the original one for a sufficiently small initial vertical velocity.

According to (5), the condition can be written in terms of the surface mass density Σ. Then it can be established that in a Newtonian model of a thin disc, a sufficient and necessary condition for vertical stability of a circular orbit of radius R is

Since we know that all disc models which attempt to represent realistic distributions of matter satisfy such a condition, then it can be established that equatorial circular orbits in thin disc models are always stable under vertical perturbations.

Therefore, by using then (19) in (17), we can easily obtain the relation

in such a way that, from (12), the epiciclic frequency can be cast as

where X = αk.

A brief analysis of (22) shows that the obtained models can have radial instabilities near the edge of the disc, because the models are adjusted with the observed rotation curves, which implies that the circular velocity will increase rapidly with the radius and then remain approximately constant. Then, the largest positive value of dv 2 c /dR occurs at R = 0 and hence so x 2 will have a positive maximum at the center of the disc, and then decrease towards the edge.

Fitting the Rotation Curve

We will start this section by pointing out the most relevant aspects of the adjustment, which we have previously published in (Gonzales & Reina, 2016). In order to adjust the previous model to the real observed data, we have taken a sample of spiral galaxies from the Ursa Major cluster. We pick the corresponding data out from Table 4 of the paper by Verheijen & Sancisi (2001), which presents the results of an extensive 21 cm-line synthesis imaging survey of 41 galaxies in the neighborhood of the Ursa Major cluster using the Westerbork Synthesis Radio Telescope. The mean distance between this telescope and the cluster is 18.6 Mpc. At this distance, 1 arcmin corresponds to 5.4 kpc.

For each rotation curve data, we take as the value of α the value given by the last tabulated radius, i.e. we are assuming that the radius of each galaxy is defined by its corresponding last observed value. Thereafter we take the radii normalized in units of α to fit the rotation curve of every galaxy by means of the model (12). The fits are made through a non-linear least squares fitting using the Levenberg Marquardt algorithm, implemented internally by ROOT (Brun & Rademakers, 1997) version 5.28, which minimizes the weighted sum of squares of deviations between the fit and the data. We assigned weights to the data points inversely proportional to the square of their errors. These errors correspond to 2vΔv, being Δv the galaxy velocity measurement error.

For each galaxy, we look initially for all the possible fits starting at m = 1 up to m = N - 1, with N the number of measured data pairs (R, v 2), and we take as value for m those for which we get the minimum reduced chi squared xr 2 (the best fit). Now we can discard the galaxies that do not pass the reduced chi squared test (Bevington & Keith, 2003) with a confidence level of 95 %. From the initial 41 galaxies, 26 have survived the reduced chi squared test. These 26 galaxies with their respective information of the fit done are tabulated in Table 1, where the columns 1, 2, 3, 4, 5, 6, 7, and 8 are the name, type of morphology, radius, number of parameters of the fitting function? mass (calculated according to our model), minimum reduced chi squared number of degrees of freedom (equal to N - 1), and the confidence level of the fit, respectively.

Table 1 Best Fit results of the UMa Cluster. 

Fig. 1 - 4 shown the twenty six adjusted rotation curves of the surviving galaxies. The points with error bars are the observations, as reported in (Verheijen & Sancisi, 2001), while the solid black line is the fitted rotation curve by using (12), where the values for the A 2n are given by the best fit. As we can see, relation (12) fits quite accurately the observed data of the referred galaxies. Starting from the values obtained for A2n, the values of the constants C2n are calculated by using the relations (13) and (16). The C2n constants are tabulated in Table 3. In a similar way, we can use (11) and (16) in order to calculate the mass for each galaxy with its respective error; these values are shown together in the fifth column of Table 1.

Figure 1 Circular Velocity in [km/s] of galaxies NGC3726, NGC3877, NGC3917, NGC3949, NGC3953, NGC3972, NGC4010 and NGC4085. The solid lines represent the rotation curve determined from (12) and the parameters of the best fit listed in Table 1. The dashed lines correspond to the rotation curve for the galaxies radially stable by using the parameters m tabulated in Table 2. The points with error bars are the observations, as reported in (Verheijen & Sancisi, 2001

Figure 2 Circular Velocity in [km/s] of galaxies NGC4088, NGC4100, NGC4157, NGC4183, NGC4217, NGC4389, UGC6399 and UGC6446. The solid lines represent the rotation curve determined from (12) and the parameters of the best fit listed in Table 1. The dashed lines correspond to the rotation curve for the galaxies radially stable by using the parameters m tabulated in Table 2. The points with error bars are the observations, as reported in (Verheijen & Sancisi, 2001

Figure 3 Circular Velocity in [km/s] of galaxies UGC6667, UGC6818, UGC6917, UGC6923, UGC6930, UGC6962, UGC6969 and UGC6973. The solid lines represent the rotation curve determined from (12) and the parameters of the best fit listed in Table 1. The dashed lines correspond to the rotation curve for the galaxies radially stable by using the parameters m tabulated in Table 2. The points with error bars are the observations, as reported in (Verheijen & Sancisi, 2001

Figure 4 Circular Velocity in [km/s] of galaxies UGC6983 and UGC7089. The solid lines represent the rotation curve determined from (12) and the parameters of the best fit listed in Table 1. The dashed lines correspond to the rotation curve for the galaxies radially stable by using the parameters m tabulated in Table 2. The points with error bars are the observations, as reported in (Verheijen & Sancisi, 2001

Given that the set of constants C2n defines completely each particular thin disc model, it is straightforward to compute all the physical quantities characterizing each galaxy. Explicit expressions for the gravitational potential ϕ (R, z) and the surface mass density Σ(R) can be easily obtained by using constants C2n at expressions (7) and (10), they are not presented here. Instead of that, we plot in Fig. 5-8 the surface densities for the twenty six galaxies, as functions of the dimensionless radial coordinate R = R/a. The solid black line represents the behavior of the surface mass density that makes use of the parameters obtained by means of the best fit. For the twenty six galaxies we obtain a well behaved surface mass density, having a maximum at the disc center, and hence, decreasing until totally vanishing at the disc’s edge. This implies, according to (21), that all galaxies are stable against vertical perturbations.

Figure 5 Surface density mass in [M / pc2] for galaxies NGC3726, NGC3877, NGC3917, NGC3949, NGC3953, NGC3972, NGC4010 and NGC4085. The solid lines represent the behavior of the surface density mass determined from (10) and the best fit parameters m listed in Table (1). The dashed lines correspond to the surface density mass for the galaxies radially stable by using the parameters m tabulated in Table (2). For all the galaxies of the sample the surface density has a physically reasonable behavior: maximum at the center of the disc and then decreasing monotonically with the radius. 

Figure 6 Surface density mass in [M / pc2] for galaxies NGC4088, NGC4100, NGC4157, NGC4183, NGC4217, NGC4389, UGC6399 and UGC6446. The solid lines represent the behavior of the surface density mass determined from (10) and the best fit parameters m listed in Table (1). The dashed lines correspond to the surface density mass for the galaxies radially stable by using the parameters m tabulated in Table (2). For all the galaxies of the sample the surface density has a physically reasonable behavior: maximum at the center of the disc and then decreasing monotonically with the radius. 

Figure 7 Surface density mass in [M / pc2] for galaxies UGC6667, UGC6818, UGC6917, UGC6923, UGC6930, UGC6962, UGC6969 and UGC6973. The solid lines represent the behavior of the surface density mass determined from (10) and the best fit parameters m listed in Table (1). The dashed lines correspond to the surface density mass for the galaxies radially stable by using the parameters m tabulated in Table (2). For all the galaxies of the sample the surface density has a physically reasonable behavior: maximum at the center of the disc and then decreasing monotonically with the radius. 

Figure 8 Surface density mass in [M / pc2] for galaxies UGC6983 and UGC7089. The solid lines represent the behavior of the surface density mass determined from (10) and the best fit parameters m listed in Table (1). The dashed lines correspond to the surface density mass for the galaxies radially stable by using the parameters m tabulated in Table (2). For all the galaxies of the sample the surface density has a physically reasonable behavior: maximum at the center of the disc and then decreasing monotonically with the radius. 

In a similar way, we can compute the epicyclic frequencies by using (23), and the values of the constants A 2n obtained from the best fit; however, as with the surface mass densities, we will not present the explicit expressions here, and we only show the corresponding plots. Fig. 9-12 show with a solid black line the plots of the epicyclic frequencies for the set of galaxies considered. As we can see, seventeen galaxies of the sample have small regions of radial instability close to the disc’s edge. For each of these radially unstable galaxies, we consider the possibility that there exists a set of parameters m, different from those of the best fit (column 4, Table 1), so that the x2 test continues to be satisfied with a confidence level of 95 %, and which in turn can stabilize some of these galaxies, without altering the physical behavior of the rotation curve and the surface mass density.

In order to guarantee a confidence level of 95 %, for each of these galaxies we take a set of parameters m = m ± 1. For each of these two values we perform the adjustment, and check whether the instability continues.

If the instability disappears, we have obtained the value of m and the parameters A 2n needed to calculate the mass of the galaxy that stabilizes the model. From this new fit, we manage to stabilize 7 galaxies, keeping the rotation curves and surface mass density with a good behavior. In Fig. 1 - 4, 5 - 8, and 9 - 12 we show with dashed lines, the rotation curve, the surface mass density and the epiciclic frequency, respectively, for these 7 galaxies. In Table 2 we show the NGC/UGC name of those galaxies stable after having done the new adjustment.

Table 2 Stable Galaxies by using the others parameters (different to the best fit ones). 

Table 3 Constants C 2n [km2s-2] and values of m

Figure 9 Epicyclic frequency in [Km 2 /s 2] for 26 galaxies NGC3726, NGC3877, NGC3917, NGC3949, NGC3953, NGC3972, NGC4010 and NGC4085. The solid lines represent the epicyclic frequency determined from (23) and the best fit parameters m listed in Table 1. The regions where the frequency are negative represent zones of the radial instability. The dashed lines correspond to the epicyclic frequency for the galaxies radially stable by using the parameters m tabulated in Table 2

Figure10. Epicyclic frequency in [Km2/s2] for 26 galaxies NGC4088, NGC4100, NGC4157, NGC4183, NGC4217, NGC4389, UGC6399 and UGC6446. The solid lines represent the epicyclic frequency determined from (23) and the best fit parameters m listed in Table 1. The regions wherethe frequency are negative represent zones of the radial instability. The dashed lines correspond to the epicyclic frequency for the galaxies radially stable by using the parameters m tabulated in Table 2

Figure 11 Epicyclic frequency in [Xm 2 / s 2] for 26 galaxies UGC6667, UGC6818, UGC6917, UGC6923, UGC6930, UGC6962, UGC6969 and UGC6973. The solid lines represent the epicyclic frequency determined from (23) and the best fit parameters m listed in Table 1. The regions where the frequency are negative represent zones of the radial instability. The dashed lines correspond to the epicyclic frequency for the galaxies radially stable by using the parameters m tabulated in Table 2

Figure 12 Epicyclic frequency in [Km 2 /s 2] for 26 galaxies for galaxies UGC6983 and UGC7089. The solid line represent the epicyclic frequency determined from (23) and the best fit parameters m listed in Table 1. The regions where the frequency are negative represent zones of the radial instability. The dashed line correspond to the epicyclic frequency for the galaxies radially stable by using the parameters m tabulated in Table 2

The number of parameters m used in this new fit are shown in column (2). In columns (4) and (5) we show the reduced chi square and the number of degrees freedom, respectively. In (3) we present the mass values in order to obtain galaxies radially stable.

Concluding Remarks

We have presented a mass model for a set of twenty six spiral galaxies of the Ursa Major Cluster. The model is adjusted very precisely with the data of the rotation curve of the galaxies of the sample, considering that all the mass is concentrated only in the stellar disc. In this way the model implies that it is not always necessary to introduce a halo of dark matter or MOND theories in order to describe properly the rotation curve of the galaxies. In all cases the best statistical fit between the data and the model is fixed to a confidence level of 95 %. As was mentioned by Gonzales et al. (2010), the equation (12) can be considered as a kind of “universal rotation curve” for flat galaxies, which can easily be adjusted to the observed data of the rotation curve of any particular spiral galaxy.

The surface mass density of each galaxy present a good behavior: maximum at the center and vanishing at the edge, which is in concordance with the observed luminosity profile of many spiral galaxies, and the obtained values for the corresponding total mass M it agrees with the expected order of magnitude. Besides, considering the stability criterion constructed by Vieira & Ramos-Caro (2016), and taking into account that all of the mass density profiles obtained here are positive, we can conclude that all the models are vertically stable, contrary to the statements shown in Gonzales et al. (2010). Thus, we remark that our formalism is a powerful method to obtain stable maximum disc describing the optical region of many HSB galaxies.

On the other hand, the analysis of the epicyclic frequency reveals that 17 galaxies of the sample present small regions of radial instability near the disc edge. From the equation (23) was expected that this unstable behavior depends only of the fit parameters A2n. We have analyzed the behavior of the radial frequency for all the possible values of m, different to those given by the best fit, but so that the confidence level remains at 95 %, and thus, we have achieved that 7 of the 17 unstable galaxies present a stable behavior.

According to the results obtained, we consider that the models presented can be considered as a first approximation for obtaining more realistic models. We particularly believe that the mass values M of the 26 galaxies modeled can be taken as an estimate of the upper limit for the mass of the disc of these galaxies, since in the models it was considered that the whole mass was concentrated in the galactic disc.

Acknowledgements

This research was partially supported by VIE-UIS, under grants number 1838 and 2314, and COLCIENCIAS, Colombia, under grant number 8840. JIR and FLS wants to thank the support from Centro de Investigaciones, Universidad Santo Tomás, Bucaramanga.

References

Arfken G, Weber H. Mathematical Methods for Physicists, 6th Ed., Academic Press, 2005. [ Links ]

Bajkova AT, Bobylev VV. Rotation curve and mass distribution in the galaxy from the velocities of objects at distances up to 200 kpc, Astronomy Letters, 42(49): 567-582, 2016. doi: 10.1134/S1063773716090012 [ Links ]

Bateman H. Partial Differential Equations, First American Edition, New York, Dover, 1944. [ Links ]

Bevington P, Keith D. Data Reduction and Error Analysis for the Physical Sciences, 3rd Ed., Mc Graw Hill. 2003. [ Links ]

Binney J, Tremaine S. Galactic Dynamics, 2nd Ed., Princeton University Press. 2008. [ Links ]

Brun R, Rademakers F. ROOT - An object oriented data analysis framework, Nuclear Instruments and Methods in Physics Research A, 389(1): 81-86, 1997. doi: 10.1016/S0168-9002(97)00048-X [ Links ]

Danby JMA. An Introduction to Stellar Astrophysics, Wiley. 1988. [ Links ]

González GA, Plata-Plata S, Ramos-Caro J. Finite thin disc models of four galaxies in the Ursa Major cluster: NGC3877, NGC3917, NGC3949 and NGC4010, Monthly Notices of the Royal Astronomical Society, 404(1): 468-474, 2010. doi: 10.1111/j.1365-2966.2010.16303.x [ Links ]

González GA, Reina JI. An infinite family of generalized Kalnajs disks, Monthly Notices of the Royal Astronomical Society, 371(4): 1873-1876, 2006. doi: 10.1111/j.1365-2966.2006.10819.x [ Links ]

González GA, Reina JI. Analytical Potentials for Flat Galaxies with Spheroidal Halos. Revista de la Academia Colombiana de Ciencias Exactas Físicas y Naturales, 40(156): 402-411, 2016. doi: 10.18257/raccefyn.376 [ Links ]

Hunter C. The structure and stability of self-gravitating disks, Monthly Notices of the Royal Astronomical Society, 126(4): 299-315, 1963. doi: 10.1093/mnras/126.4.299 [ Links ]

Kalnajs AJ. Mass distribution and dark halos: Dussion, in: Internal kinematics and dynamics of galaxies, Proceedings of the IAU Symposium, Dordrecht, D. Reidel Publishing Co., 87-88, 1983. [ Links ]

Kent SM. Dark matter in spiral galaxies. I - Galaxies with optical rotation curves, The Astronomical Journal, 91(6): 1301-1327, 1986. doi: 10.1086/114106 [ Links ]

LeBlanc F. Galactic Dynamics, 2nd Ed., Princeton University Press. 2010. [ Links ]

Oort JH. Observational evidence confirming Lindblad’s hypothesis of a rotation of the galactic system, Bulletin of the Astronomical Institutes of the Netherlands, 3(120): 275-282, 1927. [ Links ]

Palunas P, Williams TB. Maximum disc mass models for spiral galaxies, The Astronomical Journal, 120(6): 2884-2903, 2000. doi: 0.1086/316878 [ Links ]

Rubin VC, Burstein D, Ford WK, Thonnard N. Rotation Velocities of 16 SA Galaxies and a Comparison of Sa, Sb, and SC Rotation Properties, The Astrophysical Journal, 289(1): 81-104, 1985. doi: 10.1086/162866 [ Links ]

Sofue Y, Honma M, Omodaka T. Unified Rotation Curve of the Galaxy - Decomposition into de Vaucouleurs Bulge, disc, Dark Halo, and the 9-kpc Rotation Dip, Publications of the Astronomical Society of Japan, 61(2): 227-236, 2009. doi: 10.1093/pasj/61.2.227 [ Links ]

Van Albada TS, Bahcall JN, Begeman K, Sancisi R. Distribution of dark matter in the spiral galaxy NGC 3198, The Astrophysical Journal., 295(1): 305-313, 1985. doi: 10.1086/163375 [ Links ]

Verheijen MAW, Sancisi R. The Ursa Major cluster of galaxies. IV. HI synthesis observations, Astronomy and Astrophysics, 370(3): 765-867, 2001. doi: 10.1051/0004-6361:20010090 [ Links ]

Vieira RSS, Ramos-Caro J. Integrability of motion around galactic razor- thin disks, Celestial Mechanics and Dynamical Astronomy, 126(4): 483-500, 2016. doi: 10.1007/s10569-016-9705-0 [ Links ]

Funding: N.A.

Electronic supplementary material: N.A.

Jerson I. Reina He is a professor full time of mathematics and physics at the Universidad Santo Tomás. The author obtained his PhD degree in Physics at Universidad Industrial de Santander. His main research areas are astrophysics and galactic dynamics.

Framsol López-Suspes His is PhD in Natural Sciences - Physics of the Universidad Industrial de Santander. He is currently a professor at the Universidad Santo Tomás in the Department of Basic Sciences. His subjects of interest are dynamical systems, dynamics of galaxies, and general relativity.

Guillermo A. González He is Professor of Physics at Universidad Industrial de Santander. He obtained his PhD degree in Applied Mathematics at Universidade Estadual de Campinas. His research areas are Astrophysics, General Relativity and Mathematical Physics.

Citation: Reina JI, Lopez-Suspes F, Gonzalez GA. Mass of Spiral Galaxies by Means of a Maximum Disc Model, Universitas Scientiarum, 23 (2): 191-218, 2018. doi: 10.11144/Javeriana.SC23-2.mosg

Received: September 22, 2017; Accepted: May 31, 2018

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