## Galaxies with Supermassive Binary Black Holes: (II) A Model with Cuspy Galactic Density Profiles

Ing-Guey Jiang 1 and Li-Chin Yeh 2

- 1

Department of Physics and Institute of Astronomy, National Tsing-Hua University, Hsin-Chu, Taiwan 2 Department of Applied Mathematics, National Hsinchu University of Education, Hsin-Chu, Taiwan

jiang@phys.nthu.edu.tw

## ABSTRACT

The existence and uniqueness of equilibrium points, including Lagrange Points and Jiang-Yeh Points, of a galactic system with supermassive binary black holes embedded in a centrally cuspy galactic halo are investigated herein. Differing from the previous results of non-cuspy galactic profiles that Jiang-Yeh Points only exist under a particular condition, it is found here that the Lagrange Points, L2, L3, L4 and L5, Jiang-Yeh Points, JY1 and JY2, exist under general conditions. The stability analysis shows that L2, L3, JY1 and JY2 are unstable. However, L4 and L5 are only unstable when the galactic total mass is smaller than a critical mass; otherwise they become neutrally stable centers. These results will be important for further studies on the cores of early-type galaxies.

## 1. Introduction

The general structure of galaxies has been one of the most interesting astronomical subjects due to the beauty and diversity of their morphological shapes. It is also an important research field in that galaxies are building blocks of the universe; so the formation and evolution of galactic structures, which are imprints left by the interactions of galaxies, can trace the history of the universe.

Through astronomical observations, the interaction between galaxies is confirmed to be taking place frequently, for example, the mergers of mice galaxies (NGC 4674 A&amp;B), NGC 2207/IC 2163 and NGC 2623. These results have triggered many further investigations related to mergers. For example, Wu &amp; Jiang (2009) showed that the inter-galactic population

could be ejected from galaxies during merging events. Furthermore, the ring galaxies are likely to be the outcome of major mergers, as shown in Wu &amp; Jiang (2012).

On the other hand, it is generally believed that there is a supermassive black hole at the center of galaxies. A merger of two galaxies with supermassive black holes is likely to form a supermassive binary black hole (SBBH) near the center of the newly formed galactic merging system. In fact, the major mergers of two disk galaxies are believed to be the main mechanism forming elliptical galaxies. This is partially due to the fact that elliptical galaxies are generally massive and gas-poor, and major mergers can produce a larger galactic system and reduce the gaseous component by rapid star formation or gas diffusion during merging processes. Thus, the investigations on the dynamics of elliptical galaxies are usually based on the gas-poor assumption.

In order to study the dynamic evolution of SBBH, Quinlan (1996) carried out numerical experiments on the scattering processes for the restricted three-body problem. The main issue addressed in that work was the SBBH hardening. In fact, modifications of the restricted three body problem were employed to model various problems related to planetary systems (Chermnykh 1987; Papadakis 2004, 2005a, 2005b, Jiang &amp; Yeh 2006; Yeh &amp; Jiang 2006; Kushvah 2008a, 2008b, 2009, 2011a, 2011b, 2012). In addition, the discoveries of extra-solar planets have triggered many investigations in the field of planetary systems (see Jiang &amp; Ip 2001; Ji et al. 2002; Jiang &amp; Yeh 2003, 2004a, 2004b, 2004c, 2007, 2009, 2011; Jiang et al. 2003, 2006, 2007, 2009, 2010, 2013; Chatterjee et al. 2008).

Milosavljevic and Merritt (2001) used N-body simulations to investigate the orbital decay of black holes and the formation of SBBHs during galactic mergers. They found that it only takes about a million years for black holes to sink to the center and become a hard binary. However, they claimed that black holes stall at separations of sub-parsec scales. Then, Yu (2002) semi-analytically estimated the possible time-scales to form SBBHs and found that the time-scale of dynamical friction is too long for small black holes to sink into the center, so it is difficult to form SBBHs with very small mass ratio. Thus, the SBBHs with moderate mass ratios are most likely to form and survive with semi-major axes around 10 -3 to 10 pc in spherical or nearly spherical galaxies.

Moreover, in order to explain the surface brightness of NGC 3706, Kandrup et al. (2003) considered the stellar dynamics under the gravitational fields of SBBH and a fixed galactic potential. They found that the transition between the inner and outer power-law profiles predicted by a Dehnen potential is too gradual to represent real galaxies, so the Nuker law was used in their model.

Motivated by the above work, as in Kandrup et al. (2003), a model of modified restricted

three body problems was used in Jiang &amp; Yeh (2014) to study a galactic system with SBBH. Both Kandrup et al. (2003) and Jiang &amp; Yeh (2014) employed the Nuker law as their three dimensional spherical galactic density profiles because the Nuker law is a simple and neat way to present a smooth broken power-law with two power-indexes, and the transition is located at a well defined radius, i.e. the break radius. Note that the Nuker law was originally introduced in Lauer et al. (1992) to fit the two dimensional surface brightness of M32, and later for other galaxies in Lauer et al. (1995). In order to make it clear, we will use the term spherical Nuker law to name the galactic density profiles in this paper.

Jiang &amp; Yeh (2014) discovered that when the galactic density profile is a spherical Nuker law with γ = 0 (i.e. no cusp), α = 2, β = 4, the usual five Lagrange Points always exist, and the new equilibrium points, i.e. Jiang-Yeh Points, exist if and only if the total galactic mass is larger than the critical mass. The analytic expression of this critical mass was provided and the stability analysis was performed. These results give important implications for the orbital evolution near the centers of galaxies with SBBH, and lead to a possible mechanism to form the cores of early-type galaxies through the existence of unstable Jiang-Yeh Points near supermassive black holes.

However, because the cuspy density profile is expected for the dark matter halo through cosmological simulations, it will be interesting to investigate the existence of Jiang-Yeh Points in a system with SBBH and a central cusp. For spherical Nuker law, i.e. Eq.(6) of Jiang &amp; Yeh (2014), the density profile has a cuspy center when γ ≥ 1. We find that when α = 2, β = 5, γ = 1, the corresponding gravitational potential of the density profile has an analytic form. Therefore, in this paper, we study the dynamics of a system wherein a test particle moves under the influence of SBBH and a galactic potential, which is dominated by a cuspy dark matter halo following a spherical Nuker law with α = 2, β = 5, γ = 1. We will investigate the existence of equilibrium points and also perform the stability analysis of these points.

We present our model in Section 2, and the analytical results for the existence of equilibrium points are in Section 3. The bifurcation diagrams, locations of equilibrium points and zero-velocity curves are shown in Section 4. The stability analysis is presented in Section 5, and the concluding remarks are offered in Section 6.

## 2. The Model

We consider the motion of test particles influenced by the gravitational force from the central binary black holes and the galaxy. The test particles are considered to move on the

same plane of the orbital plane of binary black holes.

For the description below, those parts which are exactly the same as the equations in Jiang &amp; Yeh (2014) will be skipped (refer to Jiang &amp; Yeh 2014). After the procedure of non-dimensionalization as done in Yeh et al. (2012) and Jiang &amp; Yeh (2014), in which the scale length is set to be the break radius, i.e. L 0 = r b , the density profile of the galaxy is expressed as:

where ρ c is a constant; r is the distance (in the unit of break radius r b ) from the origin in this spherical distribution. This profile gives an NFW central cusp (Navarro, Frenk, White 1995) and also a realistic sharper outer edge, i.e. when r → 0, the density ρ ∝ r -1 is the cusp in NFW profile; when r &gt;&gt; 1, the density decays as ρ ∝ r -5 which is sharper than the NFW's profile ( ρ ∝ r -3 for r &gt;&gt; 1). Because NFW profile is not for an equilibrium galaxy (Binney &amp; Trmaine 2008), our profile is a more realistic model for the system considered here. Note that all variables, functions and parameters, such as mass, potential, radius, normalization constants etc., are all written as dimensionless quantities in this paper.

$$\rho = \rho _ { c } r ^ { - 1 } \left \{ 1 + r ^ { 2 } \right \} ^ { - 2 } \,. & & ( 1 ) \\ \text{is the distance} \left ( \text{in the unit of break radius } r _ { b } \right ) from the origin in this$$

The mass of the galaxy up to r is thus:

$$M ( r ) = 2 \pi \rho _ { c } \left \{ 1 - \frac { 1 } { 1 + r ^ { 2 } } \right \}.$$

If the total galactic mass is M g , we have ρ c = M g 2 π . The corresponding potential is:

$$V ( r ) = - 4 \pi \left [ \frac { 1 } { r } \int _ { 0 } ^ { r } \rho ( r ^ { \prime } ) r ^ { \prime 2 } d r ^ { \prime } + \int _ { r } ^ { \infty } \rho ( r ^ { \prime } ) r ^ { \prime } d r ^ { \prime } \right ] = - M _ { g } \left \{ \frac { \pi } { 2 } - \tan ^ { - 1 } r \right \},$$

Moreover, from the potential, we can obtain the gravitational force as:

$$f _ { g } ( r ) \equiv - \frac { \partial V } { \partial r } = - M _ { g } \left ( \frac { 1 } { 1 + r ^ { 2 } } \right ).$$

Considering the motion on the x -y plane of the rotating frame, binary black holes, each with mass m at ( R, 0), ( -R, 0), and a test particle located on x, y with velocity u , v , the equations of motion would be as follows:

$$\cdot \\ \begin{cases} & \frac { d x } { d t } = u \\ & \frac { d y } { d t } = v \\ & \frac { d u } { d t } = 2 n v + n ^ { 2 } x - \frac { m ( x + R ) } { r _ { 1 } ^ { 3 } } - \frac { m ( x - R ) } { r _ { 2 } ^ { 3 } } - \frac { M _ { g } x } { r } \left ( \frac { 1 } { 1 + r ^ { 2 } } \right ), \\ & \frac { d v } { d t } = - 2 n u + n ^ { 2 } y - \frac { m y } { r _ { 1 } ^ { 3 } } - \frac { m y } { r _ { 2 } ^ { 3 } } - \frac { M _ { g } y } { r } \left ( \frac { 1 } { 1 + r ^ { 2 } } \right ), \end{cases}$$

where r 1 = √ ( x + R ) 2 + y 2 , r 2 = √ ( x -R ) 2 + y 2 , r = √ x 2 + y 2 and n is the angular velocity of black holes. That is, in the inertial frame, each black hole moves along the circular orbit at r = R with an angular velocity, i.e. mean motion:

$$n = \left \{ \frac { m } { 4 R ^ { 3 } } + \frac { 1 } { R } | f _ { g } ( R ) | \right \} ^ { 1 / 2 } = \left \{ \frac { m } { 4 R ^ { 3 } } + \frac { M _ { g } } { R ( 1 + R ^ { 2 } ) } \right \} ^ { 1 / 2 }$$

The corresponding Jacobi integral of this system is similar as the one in Jiang &amp; Yeh (2006, 2014):

$$C _ { J } = - u ^ { 2 } - v ^ { 2 } + n ^ { 2 } ( x ^ { 2 } + y ^ { 2 } ) + \frac { 2 m } { r _ { 1 } } + \frac { 2 m } { r _ { 2 } } + 2 M _ { g } \left \{ \frac { \pi } { 2 } - \tan ^ { - 1 } ( \sqrt { x ^ { 2 } + y ^ { 2 } } ) \right \}.$$

## 3. The Equilibrium Points

The existence and uniqueness of equilibrium points are investigated in this section. We follow the convention introduced in Jiang &amp; Yeh (2014) in naming the equilibrium points except that due to the singular property of the central point of the system, the origin will be named as the singular point . Because the singular point is not considered as an equilibrium point in our system, there will be no stability analysis of this point. The existence and uniqueness of Lagrange Points L4 and L5 will be presented in Theorem 1. The existence and uniqueness of Lagrange Points L2 and L3, and the singularity of the origin of the system will be presented in Theorem 2. The existence and uniqueness of Jiang-Yeh Points JY1 and JY2, will be presented in Theorem 3.

In general, for System (5), equilibrium points ( x , y e e ) satisfy A x , y ( e e ) = 0 and B x , y ( e e ) = 0, where:

$$A ( x, y ) = n ^ { 2 } x - \frac { m ( x + R ) } { r _ { 1 } ^ { 3 } } - \frac { m ( x - R ) } { r _ { 2 } ^ { 3 } } - \frac { M _ { g } x } { r } \frac { 1 } { ( 1 + r ^ { 2 } ) },$$

$$B ( x, y ) = n ^ { 2 } y - \frac { m y } { r _ { 1 } ^ { 3 } } - \frac { m y } { r _ { 2 } ^ { 3 } } - \frac { M _ { g } y } { r } \frac { \ell } { ( 1 + r ^ { 2 } ) }.$$

For convenience, for y = 0, we define:

/negationslash

$$h ( y ) \equiv \frac { B ( 0, y ) } { y } = n ^ { 2 } - \frac { 2 m } { [ R ^ { 2 } + y ^ { 2 } ] ^ { 3 / 2 } } - \frac { M _ { g } } { | y | } \frac { 1 } { ( 1 + y ^ { 2 } ) }$$

and

$$k ( x ) \equiv A ( x, 0 ) = n ^ { 2 } x - \frac { m ( x + R ) } { | x + R | ^ { 3 } } - \frac { m ( x - R ) } { | x - R | ^ { 3 } } - \frac { M _ { g } x } { | x | } \frac { 1 } { ( 1 + x ^ { 2 } ) } \quad \ ( 1 1 )$$

From the results in Jiang &amp; Yeh (2006), we have the following remarks:

## Remark A:

/negationslash

For y e = 0, y e satisfies h y ( ) = 0, if and only if (0 , y e Remark B:

) is the equilibrium point of System (5).

x e satisfies k x ( ) = 0, if and only if ( x , e 0) is the equilibrium point of System (5). Then, Remark (A) will be used to study the equilibrium points L4 and L5 in Theorem 1. Remark (B) will be used for all other points in Theorem 2 and 3.

Theorem 1: The Existence and Uniqueness of Lagrange Points L4 and L5 There is one and only one ¯ y 1 &gt; 0 such that h y ( ¯ ) 1 = 0 , and only one ¯ y 2 &lt; 0 such that h y ( ¯ ) = 0 2 . That is, excluding the origin (0,0) of x -y plane, there are two equilibrium points on the y -axis, i.e. L4 and L5, for System (5).

Proof:

We define

$$P ( y ) \equiv \frac { 2 m } { [ R ^ { 2 } + y ^ { 2 } ] ^ { 3 / 2 } } - n ^ { 2 }$$

$$& \text{and} \quad Q ( y ) \equiv \frac { - M _ { g } } { | y | ( 1 + y ^ { 2 } ) } \\ \text{ave } h ( u ) = - P ( u ) + O ( u )$$

so from Eq.(10), we have h y ( ) = -P y ( ) + Q y ( ).

At first, we consider the case when y &gt; 0. Because lim y → 0 Q y ( ) = -∞ , we have:

$$\lim _ { y \to 0 } h ( y ) = \lim _ { y \to 0 } - P ( y ) + Q ( y ) = - \infty ;$$

and since lim y →∞ P y ( ) = -n 2 , lim y →∞ Q y ( ) = 0, we have lim y →∞ h y ( ) = n 2 &gt; 0. Moreover, from Eqs.(12)-(13), we have:

$$P ^ { \prime } ( y ) = - 6 m y \left [ R ^ { 2 } + y ^ { 2 } \right ] ^ { - 5 / 2 },$$

and

$$Q ^ { \prime } ( y ) = M _ { g } \frac { 1 + 3 y ^ { 2 } } { y ^ { 2 } ( 1 + y ^ { 2 } ) ^ { 2 } }.$$

Since P ′ ( y ) &lt; 0 and Q y ′ ( ) &gt; 0, h ′ ( y ) = -P ′ ( y ) + Q y ′ ( ) &gt; 0 for any y &gt; 0. Thus, h y ( ) is a monotonically increasing function for any y ∈ (0 , ∞ ). With lim y → 0 h y ( ) = -∞ and lim y →∞ h y ( ) &gt; 0, we conclude there is a unique point ¯ y 1 &gt; 0 such that h y ( ¯ ) = 0. 1

For the case when y &lt; 0, from Eqs.(15)-(16), we have P ′ ( y ) &gt; 0, Q y ′ ( ) &lt; 0, so h y ( ) is a monotonically decreasing function. We also have lim y →∞ h y ( ) = n 2 &gt; 0 and lim y → 0 h y ( ) = -∞ . Thus, we find that there is a unique point ¯ y 2 &lt; 0 such that h y ( ¯ ) = 0. 2 ✷

Now we investigate the existence of equilibrium points on the x -axis. For convenience, we define and therefore,

$$S ( x ) \equiv \frac { m ( x + R ) } { | x + R | ^ { 3 } } + \frac { m ( x - R ) } { | x - R | ^ { 3 } } - n ^ { 2 } x,$$

We also define:

$$S ( x ) = \begin{cases} & \frac { m } { ( x + R ) ^ { 2 } } + \frac { m } { ( x - R ) ^ { 2 } } - n ^ { 2 } x, \quad & \text{for $x>R$,} \\ & \frac { m } { ( x + R ) ^ { 2 } } - \frac { m } { ( x - R ) ^ { 2 } } - n ^ { 2 } x, \quad & \text{for $-R<x<R$,} \\ & - \frac { m } { ( x + R ) ^ { 2 } } - \frac { m } { ( x - R ) ^ { 2 } } - n ^ { 2 } x, \quad & \text{for $x<-R$.} \end{cases} \\ \text{te also define:} & M _ { \kappa }$$

$$T ( x ) \equiv - \frac { M _ { g } x } { | x | ( 1 + x ^ { 2 } ) }. \\ ( 1 9 ) \text{, we have} \colon$$

From Eqs.(11), (17) and (19), we have:

$$k ( x ) = - S ( x ) + T ( x ).$$

In the following Theorem 2, the results show that Lagrange Points L2 and L3 exist. However, due to the singular property, the origin (0,0) is a singular point, though the net force shall be physically balanced there. In Theorem 3, it is shown that another two equilibrium points, JY1 and JY2, exist. JY1 is in the region ( -R, 0) and JY2 is in the region (0 , R ) of x-axis.

## Theorem 2: The Existence and Uniqueness of Lagrange Points L2 and L3 and the Singularity of the Origin Point (0,0)

- (i) There is an unique x 1 &gt; R such that k x ( 1 ) = 0 and an unique x 2 &lt; -R such that k x ( 2 ) = 0. That is, on the x -axis, there is one and only one equilibrium point in the region ( R, ∞ ), i.e. L2, and there is one and only one equilibrium point in the region ( -∞ -, R ), i.e. L3.

/negationslash

- (ii) Due to the singularity at x = 0, the origin is not an equilibrium point (i.e. k (0) = 0), but is a singular point.

## Proof of (i):

When x &gt; R , from Eqs. (18)-(20),

$$k ( x ) = - S ( x ) + T ( x ) = - \frac { m } { ( x + R ) ^ { 2 } } - \frac { m } { ( x - R ) ^ { 2 } } + n ^ { 2 } x - \frac { M _ { g } } { 1 + x ^ { 2 } } \\ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \colon$$

$$\text{and} \ \ k ^ { \prime } ( x ) = - S ^ { \prime } ( x ) + T ^ { \prime } ( x ) = \frac { 2 m } { ( x + R ) ^ { 3 } } + \frac { 2 m } { ( x - R ) ^ { 3 } } + n ^ { 2 } + \frac { 2 x M _ { g } } { ( 1 + x ^ { 2 } ) ^ { 2 } } > 0.$$

Thus, we have lim x → R + k x ( ) = -∞ , lim x →∞ k x ( ) = ∞ . Moreover, from k ′ ( x ) &gt; 0, we know that k x ( ) is a monotonic function. Therefore, there is a unique x 1 &gt; R , such that k x ( 1 ) = 0.

When x &lt; -R , similarly, since k ′ ( x ) &gt; 0, lim x →-R -k x ( ) = ∞ , and lim x →-∞ k x ( ) = -∞ , there is a unique x 2 &lt; -R such that k x ( 2 ) = 0. ✷

## Proof of (ii):

When -R &lt; x &lt; R , from Eq. (18), we have:

$$S ( x ) & = \frac { m } { ( x + R ) ^ { 2 } } - \frac { m } { ( x - R ) ^ { 2 } } - n ^ { 2 } x, \\ \intertext { n \dots \dots \dots } \, \stackrel { n \dots } { n \dots } \, \sum \nolimits \, \dots \, \dots.$$

Thus, S (0) = 0. Moreover, from Eq.(19), we have :

$$\lim _ { x \to 0 ^ { + } } T ( x ) = \lim _ { x \to 0 ^ { + } } - \frac { M _ { g } } { ( 1 + x ^ { 2 } ) } = - M _ { g }$$

and so lim x → 0 T x ( ) does not exist. Therefore, k x ( ) is not a continuous function at x = 0, and so the origin (0,0) is a singular point. ✷

$$\lim _ { x \to 0 ^ { - } } T ( x ) = \lim _ { x \to 0 ^ { - } } - & \frac { M _ { g } x } { ( - x ) ( 1 + x ^ { 2 } ) } = M _ { g }, \\ \text{s not exist.} \text{ Therefore, $k(x)$ is not a continuous function}$$

## Theorem 3: The Existence and Uniqueness of Jiang-Yeh Points JY1 and JY2

There is a unique x 3 ∈ (0 , R ) such that k x ( 3 ) = 0 and a unique x 4 ∈ ( -R, 0) such that k x ( 4 ) = 0. That is, on the x -axis, there is one and only one equilibrium point in the region (0 , R ), i.e. JY1, and there is one and only one equilibrium point in the region ( -R, 0), i.e. JY2.

## Proof:

When -R &lt; x &lt; R , from Eq.(21), we have:

$$S ^ { \prime } ( x ) = - \frac { 2 m } { ( x + R ) ^ { 3 } } + \frac { 2 m } { ( x - R ) ^ { 3 } } - n ^ { 2 } < 0.$$

We first consider the region with 0 &lt; x &lt; R . From Eq.(19):

$$T ^ { \prime } ( x ) = - \frac { - 2 x M _ { g } } { ( 1 + x ^ { 2 } ) ^ { 2 } } = \frac { 2 x M _ { g } } { ( 1 + x ^ { 2 } ) ^ { 2 } } > 0,$$

we have k ′ ( x ) = -S ′ ( x ) + T ′ ( x ) &gt; 0, and thus k x ( ) is a monotonic function in this region. Due to the singularity at 0 for T x ( ) and at R for S(x), we consider the limits and find

lim x → 0 + k x ( ) = -S (0)+lim x → 0 + T x ( ) = -M g and lim x → R -k x ( ) = -lim x → R -S x ( )+ T R ( ) = ∞ . Thus, there is a unique x 3 ∈ (0 , R ) such that k x ( 3 ) = 0.

Similarly, for the region with -R &lt; x &lt; 0, from Eq.(19):

$$T ^ { \prime } ( x ) = \frac { - 2 x M _ { g } } { ( 1 + x ^ { 2 } ) ^ { 2 } } > 0,$$

so we have k ′ ( x ) = -S ′ ( x ) + T ′ ( x ) &gt; 0 and k x ( ) is a monotonic function in this region. Because lim x →-R + k x ( ) = -∞ and lim x → 0 -k x ( ) = M &gt; g 0, there is a unique x 4 ∈ ( -R, 0) such that k x ( 4 ) = 0. ✷

## 4. The Bifurcations and Zero-Velocity Curves

In order to demonstrate the analytic results proved above, we numerically determine the locations of equilibrium points by solving k x ( ) = 0. The locations of the equilibrium points on the x-axis with m = 1 as a function of M g are shown in Figs.1-2. Fig. 1(a) is for R = 0 25, Fig. . 1(b) is for R = 0 5, Fig. . 2(a) is for R = 1 and Fig. 2(b) is for R = 2.

It is clear that when M g = 0, there are three equilibrium points, L1, L2 and L3 on x -axis; when M &gt; g 0, there are four equilibrium points, L2, L3, JY1 and JY2 on x -axis. The separations between these equilibrium points are larger for larger R . Moreover, it is interesting that when M g increases from 0 to 20, the separation between L2 and L3 decreases, but the separation between JY1 and JY2 increases.

In Fig. 3, the zero-velocity curves, which are obtained through Eq.(7), of the case with m = 1 and R = 1 are presented. Fig. 3(a) is for M g = 1, Fig. 3(b) is for M g = 10, Fig. 3(c) is for M g = 30 and Fig. 3(d) is for M g = 100. The + signs are the locations for Lagrange Points, and the squares indicate the Jiang-Yeh Points. These points are numerically determined by solving k x ( ) = 0 and h y ( ) = 0. It is shown that they are completely consistent with the locations of equilibrium points implied by zero-velocity curves. The locations and the values of Jacobi integral C J of all equilibrium points shown in Fig. 3 are summarized in Table 1.

Table 1. The Locations and C J of Equilibrium Points

|           |               | L2       | L3        | L4       | L5        | JY1      | JY2       |
|-----------|---------------|----------|-----------|----------|-----------|----------|-----------|
|           | ( x e , y e ) | (1.94,0) | (-1.94,0) | (0,1.29) | (0,-1.29) | (0.19,0) | (-0.19,0) |
|           | C J           | 6.58     | 6.58      | 5.02     | 5.02      | 6.94     | 6.94      |
|           |               | L2       | L3        | L4       | L5        | JY1      | JY2       |
|           | ( x e , y e ) | (1.47,0) | (-1.47,0) | (0,1.04) | (0,-1.04) | (0.56,0) | (-0.56,0) |
|           | C J           | 28.36    | 28.36     | 23.77    | 23.77     | 28.67    | 28.67     |
|           |               | L2       | L3        | L4       | L5        | JY1      | JY2       |
|           | ( x e , y e ) | (1.33,0) | (-1.33,0) | (0,1.01) | (0,-1.01) | (0.69,0) | (-0.69,0) |
|           | C J           | 72.58    | 72.58     | 65.20    | 65.20     | 72.90    | 72.90     |
| M g = 100 |               | L2       | L3        | L4       | L5        | JY1      | JY2       |
| M g = 100 | ( x e , y e ) | (1.22,0) | (-1.22,0) | (0,1.0)  | (0,-1.0)  | (0.79,0) | (-0.79,0) |
| M g = 100 | C J           | 222.11   | 222.11    | 210.16   | 210.16    | 222.43   | 222.43    |

## 5. The Stability of Equilibrium Points

After the existence of equilibrium points is confirmed, and the locations of equilibrium points are determined, it would be interesting to understand the stability around these points. We now consider the following system:

$$\begin{cases} \quad \frac { d x } { d t } = u, \\ \quad \frac { d y } { d t } = v, \\ \quad \frac { d u } { d t } = 2 n v + A ( x, y ), \\ \quad \frac { d v } { d t } = - 2 n u + B ( x, y ), \end{cases} \\ ) \, \text{are defined in Eqs.} ( 8 ) - ( 9 ).$$

where A x, y ( ), B x, y ( ) are defined in Eqs.(8)-(9).

To study the stability of equilibrium points, we need to know the properties of the eigenvalues of equilibrium points. The characteristic equation of the eigenvalue λ is

$$\lambda ^ { 4 } + ( 4 n ^ { 2 } - A _ { x } - B _ { y } ) \lambda ^ { 2 } + 2 n ( A _ { y } - B _ { x } ) \lambda + A _ { x } B _ { y } - B _ { x } A _ { y } = 0,$$

where A x ≡ ∂A x, y /∂x ( ) , A y ≡ ∂A x, y /∂y ( ) , B x ≡ ∂B x, y /∂x ( ) , and B y ≡ ∂B x, y /∂y ( ) . Thus,

$$A _ { x } = \sum _ { \infty } ^ { 2 } - \frac { m } { r _ { 1 } ^ { 3 } } - \frac { m } { r _ { 2 } ^ { 3 } } + \frac { 3 m ( x + R ) ^ { 2 } } { r _ { 1 } ^ { 5 } } + \frac { 3 m ( x - R ) ^ { 2 } } { r _ { 2 } ^ { 5 } } - \frac { M _ { g } } { r ( 1 + r ^ { 2 } ) } + \frac { M _ { g } x ^ { 2 } ( 1 + 3 r ^ { 2 } ) } { r ( r + r ^ { 3 } ) ^ { 2 } }, ( 2 7 )$$

$$A _ { y } = \frac { 3 m ( x + R ) y ^ { 2 } } { r _ { 1 } ^ { 5 } } + \frac { 3 m ( x - R ) y } { r _ { 2 } ^ { 5 } } + \frac { M _ { g } x y ( 1 + 3 r ^ { 2 } ) } { r ( r + r ^ { 3 } ) ^ { 2 } }, \quad \cdots \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \occ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \OC },$$

$$B _ { x } = \frac { 3 m y ( \stackrel { \ell } { x } + R ) } { r _ { 1 } ^ { 5 } } + \frac { 3 m y ( \stackrel { \ell } { x } - R ) } { r _ { 2 } ^ { 5 } } + \frac { \stackrel { \ell } { M _ { g } x y ( 1 + 3 r ^ { 2 } ) } } { r ( r + r ^ { 3 } ) ^ { 2 } },$$

$$B _ { y } = n ^ { 2 } - \frac { m } { r _ { 1 } ^ { 3 } } - \frac { m } { r _ { 2 } ^ { 3 } } + \frac { 3 m y ^ { 2 } } { r _ { 1 } ^ { 5 } } + \frac { 3 m y ^ { 2 } } { r _ { 2 } ^ { 5 } } - \frac { M _ { g } } { r ( 1 + r ^ { 2 } ) } + \frac { M _ { g } y ^ { 2 } ( 1 + 3 r ^ { 2 } ) } { r ( r + r ^ { 3 } ) ^ { 2 } }.$$

From Eqs.(28) and (29), for any x e we have A y ( x , e 0) = B x ( x , e 0) = 0, and for any y e we have A y (0 , y e ) = B x (0 , y e ) = 0. In order to do further investigation, some parameters need to be specified: we set m = 1 and R = 1 for all the results in this section. Thus, from Eq.(6) we have:

$$n ^ { 2 } = \frac { 1 } { 4 } + \frac { M _ { g } } { 2 }.$$

At first, we consider the equilibrium point L2 and JY1, ( x , y e e ), which satisfies k x ( e ) = 0 with x e &gt; 0 and y e = 0. Due to A y ( x , e 0) = 0 and B x ( x , e 0) = 0, Eq.(26) becomes:

$$\lambda ^ { 4 } + ( 4 n ^ { 2 } - A _ { x } - B _ { y } ) \lambda ^ { 2 } + A _ { x } B _ { y } = 0.$$

For convenience, we define Ω = A B x y and Π ≡ A x + B y -4 n 2 . Therefore, we have roots :

$$\max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \\ \lambda _ { + } ^ { 2 } = \frac { \Pi - \sqrt { \Pi ^ { 2 } - 4 \Omega } } } { 2 }.$$

Moreover, A x ( x , e 0) and B y ( x , e 0) can be expressed as:

$$& A _ { x } ( x _ { e }, 0 ) = n ^ { 2 } + \frac { 2 } { | x _ { e } + 1 | ^ { 3 } } + \frac { 2 } { | x _ { e } - 1 | ^ { 3 } } + \frac { 2 x _ { e } M _ { g } } { ( 1 + x _ { e } ^ { 2 } ) ^ { 2 } } > 0, \\ & \sim \varphi \quad \sim \quad 1 \, / \quad 1 \, 1 \, \bigwedge \quad \sim \quad \varphi \, \quad \, \varphi \, \quad \, \lambda \, \lambda$$

$$B _ { y } ( x _ { e }, 0 ) & = \frac { 1 } { x _ { e } } \left ( \frac { \lambda \lambda _ { 1 } ^ { 2 } } { | x _ { e } + 1 | ^ { 3 } } - \frac { \lambda \lambda _ { 1 } ^ { 2 } } { | x _ { e } - 1 | ^ { 3 } } \right ) < 0 \\ \intertext { m a l i w e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e }$$

(see Appendix A for details).

For L2, x e &gt; R = 1, from Eqs. (33)-(34), since A x ( x , e 0) &gt; 0 and B y ( x , e 0) &lt; 0, we have Ω = A x ( x , e 0) B y ( x , e 0) &lt; 0. Thus, Π 2 -4Ω &gt; 0, and we have λ 2 + &gt; 0 and λ 2 -&lt; 0. As in Szebehely (1967); this indicates that it is an unstable equilibrium point.

For JY1, it has 0 &lt; x e &lt; 1 and y e = 0. From Eqs. (33)-(34), A x ( x , e 0) &gt; 0 and B y ( x , e 0) &lt; 0, so Ω = A B x y &lt; 0. Thus, Π 2 -4Ω &gt; 0, we have λ 2 + &gt; 0 and λ 2 -&lt; 0. Therefore, JY1 is also an unstable equilibrium point.

Because our system is symmetric with respect to the y -axis, the above results are also valid for L3 and JY2. Thus, the equilibrium points L3 and JY2 are unstable.

Secondly, we study the equilibrium point L4, which can be written as (0 , y e ) with y e &gt; 0. As mentioned previously, we have A y (0 , y e ) = B x (0 , y e ) = 0. Thus, Eqs.(31) and (32) are also valid here. We find that:

$$A _ { x } ( 0, y _ { e } ) = \frac { 6 } { ( 1 + y _ { e } ^ { 2 } ) ^ { 5 / 2 } } > 0,$$

but cannot determine the sign of B y (0 , y e ) analytically here. Thus, for each given M g , the location of L4, (0 , y e ), and then the corresponding value of B y (0 , y e ), Π 2 -4Ω, and Π are determined numerically, as shown in Fig. 4(a)-(c). From Fig. 4(a), we know B y (0 , y e ) &gt; 0 and so Ω = A B x y &gt; 0. Fig. 4(b) shows the value of Π 2 -4Ω as a function of M g . It is clear that Π 2 -4Ω is not a monotonic function of M g . There is a critical value M cr ∼ 4 213, such . that when 0 &lt; M g &lt; M cr , we have Π 2 -4Ω &lt; 0; when M g ≥ M cr , we have Π 2 -4Ω ≥ 0. For the case when Π 2 -4Ω &lt; 0, both λ 2 + and λ 2 -are complex numbers. This leads to both λ + and λ -having a root which contains a positive real part, so that L4 is unstable. For the case when Π 2 -4Ω ≥ 0, we need to know the value of Π. As shown in Fig. 4(c), we find that Π &lt; 0 for the considered value of M g . We thus have λ 2 + &lt; 0 and λ 2 -&lt; 0. This leads to both λ + and λ -being pure imaginary numbers, so that L4 is a center.

Because our system is symmetric with respect to the x -axis, the above results are also valid for L5. Thus, the equilibrium point L5 is either an unstable point or a center.

## 6. Concluding Remarks

We have studied the existence and uniqueness of equilibrium points, including Lagrange Points and Jiang-Yeh Points, of a galactic system with supermassive binary black holes embedded in a central cuspy galactic halo. Due to the cuspy density profile and focusing on the case with an equal mass binary black hole, we found that the central origin is a singular point. We also found that the Lagrange Points, L2, L3, L4 and L5, and Jiang-Yeh Points, JY1 and JY2, always exist, i.e. there are six equilibrium points in the considered system. This differs from the previous results of non-cuspy galactic profiles that Jiang-Yeh Points only exist under a particular condition (Jiang &amp; Yeh 2014).

The stability analysis was performed for these equilibrium points. It is found that the equilibrium points L2, L3, JY1 and JY2 are unstable. The equilibrium point L4 (L5) is unstable when the galactic total mass M &lt; M g cr , and is a neutrally stable center when M g ≥ M cr . This critical mass M cr , which is about 4.213, is therefore an important condition for the stability of L4 and L5. These new results will be employed to investigate the cores of early-type galaxies in the near future.

## Acknowledgment

We thank the referee for very helpful suggestions. We are grateful to the National Center for High-performance Computing for computer time and facilities. This work is

supported in part by the National Science Council, Taiwan, under Ing-Guey Jiang's Grants NSC 100-2112-M-007-003-MY3 and Li-Chin Yeh's Grants NSC 100-2115-M-134-004.

## REFERENCES

Binney J., Tremaine S., 2008, Galactic Dynamics, Second Edition (Princeton University Press)

Chatterjee, S., Ford, E. B., Matsumura, S., Rasio, F. A., 2008, ApJ, 686, 580

Chermnykh, S. V., 1987, Vest. Leningrad Univ. 2, 10

- Ji, J., Li, G., Liu, L., 2002, ApJ, 572, 1041
- Jiang, I.-G., Ip, W.-H., 2001, A&amp;A, 367, 943
- Jiang, I.-G., Ip, W.-H., Yeh, L.-C., 2003, ApJ, 582, 449
- Jiang, I.-G., Yeh, L.-C., 2003, Int. J. Bifurcation and Chaos, 13, 617
- Jiang, I.-G., Yeh, L.-C., 2004a, AJ, 128, 923
- Jiang, I.-G., Yeh, L.-C., 2004b, Int. J. Bifurcation and Chaos, 14, 3153
- Jiang, I.-G., Yeh, L.-C., 2004c, MNRAS, 355, L29
- Jiang, I.-G., Yeh, L.-C., 2006, Astrophysics and Space Science, 305, 341
- Jiang, I.-G., Yeh, L.-C., 2007, ApJ, 656, 534
- Jiang, I.-G., Yeh, L.-C., 2009, AJ, 137, 4169
- Jiang, I.-G., Yeh, L.-C., 2011, MNRAS, 415, 2859
- Jiang, I.-G., Yeh, L.-C., Hung, W.-L., Yang, M.-S., 2006, MNRAS, 370, 1379
- Jiang, I.-G., Yeh, L.-C., Chang, Y.-C., Hung, W.-L., 2007, AJ, 134, 2061
- Jiang, I.-G., Yeh, L.-C., Chang, Y.-C., Hung, W.-L., 2009, AJ, 137, 329
- Jiang, I.-G., Yeh, L.-C., Chang, Y.-C., Hung, W.-L., 2010, ApJS, 186, 48
- Jiang, I.-G. et al., 2013, AJ, 145, 68
- Jiang, I.-G., Yeh, L.-C., 2014, Astrophysics and Space Science, 349, 881
- Kandrup, H. E., Sideris, I. V., Terzic, B. Bohn, C. L., 2003, ApJ, 597, 111

Kushvah, B. S., 2008a, Astrophysics and Space Science, 315, 231

Kushvah, B. S., 2008b, Astrophysics and Space Science, 318, 41 Kushvah, B. S., 2009, Astrophysics and Space Science, 323, 57 Kushvah, B. S., 2011a, Astrophysics and Space Science, 332, 99 Kushvah, B. S., 2011b, Astrophysics and Space Science, 333, 49 Kushvah, B. S., Kishor, R., Dolas, U., 2012, Astrophysics and Space Science, 337, 115 Lauer, T. R. et al., 1992, AJ, 104, 552 Lauer, T. R. et al., 1995, AJ, 110, 2622 Milosavljevic, M., Merritt, D., 2001, ApJ, 563, 34 Navarro, J. F., Frenk, C. S., White, S. D. M., 1995, MNRAS, 275, 720 Papadakis, K. E., 2004, A&amp;A, 425, 1133 Papadakis, K. E., 2005a, Astrophysics and Space Science, 299, 67 Papadakis, K. E., 2005b, Astrophysics and Space Science, 299, 129 Quinlan, G. D., 1996, New Astronomy, 1, 35 Szebehely, V., 1967, Theory of Orbits: The Restricted Problem of Three Bodies (London: Academic Press, Inc.) Wu, Y.-T., Jiang, I.-G., 2009, MNRAS, 399, 628 Wu, Y.-T., Jiang, I.-G., 2012, ApJ, 745, 105 Yeh, L.-C., Chen, Y.-C., Jiang, I.-G., 2012, Int. J. Bifurcation and Chaos, 22, 1230040 Yeh, L.-C., Jiang, I.-G., 2006, Astrophysics and Space Science, 306, 189

Yu, Q., 2002, MNRAS, 331, 935

This preprint was prepared with the AAS LT E X macros v5.0. A

## Appendix A

The details of the calculations of A x ( x , e 0) and B y ( x , e 0) in Section 5 are presented here. If m = R = 1, then r 2 1 = ( x +1) + 2 y 2 and r 2 2 = ( x -1) 2 + y 2 . From Eq. (27),

$$& \, = R = 1, \text{ then } r _ { 1 } ^ { 2 } = ( x + 1 ) ^ { 2 } + y ^ { 2 } \text{ and } r _ { 2 } ^ { 2 } = ( x - 1 ) ^ { 2 } + y ^ { 2 }. \text{ From Eq. (27),} \\ & \, A _ { x } ( x _ { e }, 0 ) = n ^ { 2 } - \frac { 1 } { | x _ { e } + 1 | ^ { 3 } } - \frac { 1 } { | x _ { e } - 1 | ^ { 3 } } + \frac { 3 ( x _ { e } + 1 ) ^ { 2 } } { | x _ { e } + 1 | ^ { 5 } } - \frac { 3 ( x _ { e } - 1 ) ^ { 2 } } { | x _ { e } - 1 | ^ { 5 } } \\ & \quad \, \quad - \frac { M _ { g } } { | x _ { e } | ( 1 + x _ { e } ^ { 2 } ) } + \frac { M _ { g } x _ { e } ^ { 2 } ( 1 + 3 x _ { e } ^ { 2 } ) } { | x _ { e } | ^ { 3 } ( 1 + x _ { e } ^ { 2 } ) ^ { 2 } } \\ & \, = n ^ { 2 } + \frac { 2 } { | x _ { e } + 1 | ^ { 3 } } + \frac { 2 } { | x _ { e } - 1 | ^ { 3 } } - \frac { M _ { g } ( 1 + x _ { e } ^ { 2 } ) } { | x _ { e } | ( 1 + x _ { e } ^ { 2 } ) ^ { 2 } } + \frac { M _ { g } ( 1 + 3 x _ { e } ^ { 2 } ) } { | x _ { e } | ( 1 + x _ { e } ^ { 2 } ) ^ { 2 } } \\ & \, = n ^ { 2 } + \frac { 2 } { | x _ { e } + 1 | ^ { 3 } } + \frac { 2 } { | x _ { e } - 1 | ^ { 3 } } + \frac { 2 | x _ { e } | M _ { g } } { ( 1 + x _ { e } ^ { 2 } ) ^ { 2 } } > 0. \\ \\ & \text{ On the other hand, from Remark } B, \text{ since } ( x _ { e }, 0 ) \text{ is an equilibrium point, } k ( x _ { e } ) = 0$$

On the other hand, from Remark B , since ( x , e 0) is an equilibrium point, k x ( e ) = 0. By Eq.(11) and m = R = 1, we have

$$- \frac { M _ { g } } { | x _ { e } | ( 1 + x _ { e } ^ { 2 } ) } & = \frac { ( x _ { e } + 1 ) } { x _ { e } | x _ { e } + 1 | ^ { 3 } } + \frac { ( x _ { e } - 1 ) } { x _ { e } | x _ { e } - 1 | ^ { 3 } } - n ^ { 2 }. \\ \tilde { \wedge \wedge } \quad \cdot = \tilde { \wedge \wedge } \quad \cdot$$

From Eq. (30) and Eq. (36), we have:

$$\Pi \, \mathcal { L } _ { \ell } \, \langle \mathfrak { O } _ { \ell } \, \text{and} \, \mathfrak { L } _ { \ell } \, \langle \mathfrak { O } _ { \ell } \, \text{,} \, \text{we move.} \\ & B _ { y } ( x _ { e }, 0 ) = n ^ { 2 } - \frac { 1 } { | x _ { e } + 1 | ^ { 3 } } - \frac { 1 } { | x _ { e } - 1 | ^ { 3 } } - \frac { M _ { g } } { | x _ { e } | ( 1 + x _ { e } ^ { 2 } ) } \\ & = n ^ { 2 } - \frac { x _ { e } } { x _ { e } | x _ { e } + 1 | ^ { 3 } } - \frac { x _ { e } } { x _ { e } | x _ { e } - 1 | ^ { 3 } } + \frac { ( x _ { e } + 1 ) } { x _ { e } | x _ { e } + 1 | ^ { 3 } } + \frac { ( x _ { e } - 1 ) } { x _ { e } | x _ { e } - 1 | ^ { 3 } } - n ^ { 2 } \\ & = \frac { 1 } { x _ { e } } \left ( \frac { 1 } { | x _ { e } + 1 | ^ { 3 } } - \frac { 1 } { | x _ { e } - 1 | ^ { 3 } } \right ) < 0.$$

Fig. 1.- The locations of the equilibrium points on the x -axis with m = 1 as a function of M g . (a) is for R = 0 25 and (b) is for . R = 0 5. .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0070v1-figure-cba18bf8-c0bb-4587-9652-6699c7bb6939.jpg)

Fig. 2.- The locations of the equilibrium points on the x -axis with m = 1 as a function of M g . (a) is for R = 1 and (b) is for R = 2.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0070v1-figure-8d644d97-5524-4902-9928-19e4329983c1.jpg)

Fig. 3.- The zero-velocity curves of the system when R = 1 and m = 1, on which the corresponding values of the Jacobi integral C J are labeled. (a) is for M g = 1, (b) is for M g = 10, (c) is for M g = 30, (d) is for M g = 100. The + signs indicate the locations of Lagrange Points and the squares indicate the locations of Jiang-Yeh Points.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0070v1-figure-12da4265-051c-43b8-8b22-15152dcb7632.jpg)

Fig. 4.(a) The value of B y (0 , y e ) as a function of M g for L4. (b) Π 2 -4Ω as a function of M g for L4. (c) Π as a function of M g for L4.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0070v1-figure-6ffba1a3-3cd6-43e2-a552-fea97330c926.jpg)
Ing-Guey Jiang Li-Chin Yeh
Galaxies with Supermassive Binary Black Holes: (II) A Model with Cuspy Galactic Density Profiles
The existence and uniqueness of equilibrium points, including Lagrange Points
and Jiang-Yeh Points, of a galactic system with supermassive binary black holes
embedded in a centrally cuspy galactic halo are investigated herein. Differing
from the previous results of non-cuspy galactic profiles that Jiang-Yeh Points
only exist under a particular condition, it is found here that the Lagrange
Points, L2, L3, L4 and L5, Jiang-Yeh Points, JY1 and JY2, exist under general
conditions. The stability analysis shows that L2, L3, JY1 and JY2 are unstable.
However, L4 and L5 are only unstable when the galactic total mass is smaller
than a critical mass; otherwise they become neutrally stable centers. These
results will be important for further studies on the cores of early-type
galaxies.