Issue 
A&A
Volume 600, April 2017



Article Number  A17  
Number of page(s)  6  
Section  Extragalactic astronomy  
DOI  https://doi.org/10.1051/00046361/201630071  
Published online  21 March 2017 
Offcenter dark matter halo leading to strong central disk lopsidedness
^{1} Department of Physics, Indian Institute of Technology Bombay, Mumbai 400076, India
email: prasadchatty@phys.iitb.ac.in
^{2} Department of Physics, Indian Institute of Science, Bangalore 560012, India
email: cjjog@physics.iisc.ernet.in
Received: 16 November 2016
Accepted: 6 January 2017
There is increasing evidence from simulations and observations that the center of the dark matter halo in a Milky Waytype galaxy could be offcenter by a few 100 pc with respect to the galactic disk. We study the effect of such an offset halo on the orbits and kinematics in the central few kpc of the disk via a simple, analytical model. The equations of motion in the disk plane can be written in terms of the potentials of the disk and halo when they are concentric, and a perturbation term arising due to the offset halo. We show that this perturbation potential has an m = 1 azimuthal variation, or is lopsided, and its magnitude increases at small radii. On solving these equations, we find that the perturbed orbit shows a large deviation, ~ 40% in radius at R = 1.5 kpc, and also strong kinematical lopsidedness. Thus, even a small halo offset of 350 pc can induce surprisingly strong spatial and kinematical lopsidedness in the central region within a ~3 kpc radius. The disk lopsidedness would have important implications for the evolution of this region; for example, it could help fuel the central active galactic nucleus.
Key words: galaxies: halos / galaxies: kinematics and dynamics / galaxies: spiral / galaxies: structure / Galaxy: center
© ESO, 2017
1. Introduction
It is well known that a galactic disk is located within an extended dark matter halo, whose existence is deduced from a nearly flat rotation curve (Rubin 1983; Binney & Tremaine 1987). The dark matter halo dominates the galaxy mass and hence the disk dynamics in the outer regions of a spiral galaxy. In the inner central region of ~ 5 kpc radius, however, the disk dominates the dynamics as it does in the Milky Way (Sackett 1997; Klypin et al. 2002). In theoretical as well as observational work, the dark matter halo is generally assumed to be concentric with respect to the disk for simplicity.
There is increasing evidence that the density peak of the dark matter halo of a Milky Waytype galaxy may be located offcenter with respect to the disk center, as seen in the simulations (Kuhlen et al. 2013) where the typical separation is 300−400 pc. Once it is generated, the offset feature seems to be longlived, lasting for several Gyr. Such an offset is of interest because it could provide a dark matter explanation of the peak in the 130 GeV in the gammaray line emission that is displaced by about 200 pc from Sgr A* in the Galactic plane (Kuhlen et al. 2013). There is also observational evidence for the halo being offcenter with respect to the disk in a field galaxy, M 99 (Chemin et al. 2016), and also for a galaxy in the galaxy cluster Abell 3827 (Massey et al. 2015). The dynamical origin of this offset and its longlived nature are not wellunderstood (see Kuhlen et al. 2013). Massey et al. (2015) claim that the offset seen in a cluster galaxy could be due to an astrophysical effect or could be evidence for selfinteracting dark matter. A lower resolution may only give an upper limit to the offset; for example, Schaller et al. (2015) find that more than 95% of the simulated galaxies in their sample have an offset between the stars and dark matter halo that is smaller than the gravitational softening length of 700 pc they use. The dynamics of an offset halo in a disk could well have some similarities with the unsettled behaviour seen in central regions of galaxies (Miller & Smith 1992). Such “sloshing” or unrelaxed central region is seen in advanced mergers as well as in isolated galaxies (Jog & Maybhate 2006), though this phenomenon has received little attention.
In this paper we study how the spatial distribution and the kinematics in the central disk are affected by the offcenter dark matter halo, in a region where the disk is the dominant mass component. We do this in terms of a simple, idealized analytical model. The halo potential can be written in terms of its value with respect to the disk center and a perturbation which takes a lopsided form. Next, the equations of motion in the disk plane are solved under this joint disk and halo system. We find that the offset halo results in a surprisingly large disturbance over the central few kpc. This effect is strongly dependent on the radius. For typical exponential disk and pseudoisothermal halo parameters for the Milky Way, and assuming an offset of 350 pc as seen by Kuhlen et al. (2013) in their simulations, the resulting radial variation in the orbit is ~ 40% at R = 1.5 kpc which drops down to ~ 3% at 3 kpc. The kinematics is also strongly perturbed. Thus, even a small offset of 350 pc can induce significant spatial and kinematical lopsidedness in the central ≤ 3 kpc region.
Section 2 contains the formulation of equations and Sect. 3 contains the results. Sections 4 and 5 contain the discussion and conclusion, respectively.
2. Formulation of equations
We use cylindrical coordinates (R, φ, z). The disk is taken to have an exponential surface density distribution (Freeman 1970), and the dark matter halo is taken to have a screened or pseudoisothermal spherical density profile (Binney & Tremaine 1987). The disk is assumed to be razorthin, and the dark matter halo center is assumed to be in the same plane as the galactic disk plane, for simplicity. Thus, we only treat the planar disturbance arising due to the halo offset.
2.1. Potential due to a disk and offcenter halo
The potential at a point P in the galactic disk plane due to the disk and the offcenter halo can be written in terms of the disk and halo contributions with respect to the disk center, plus the additional perturbation term which arises due to the offset between the disk and the halo centers. We show later that this additional term is indeed a perturbation. This follows from the fact that in a typical spiral galaxy like our Galaxy, the disk dominates the halo dynamics in the central region.
The disk center and the halo center are denoted by C_{D} and C_{H}, and the offset between the two is denoted by a distance d. The schematic diagram is given in Fig. 1. Here φ is the angle between the radius vector to the point P from the disk center and the direction between the disk center and the halo center. We note that the figure is not drawn to scale, where for the sake of clarity we have shown d ~ R, whereas in our calculation we assume d ≪ R.
Fig. 1 Schematic diagram showing the centers of disk and an offcenter dark matter halo, separated by a distance d, and their effect at point P in the disk plane. 
The density distribution of an exponential disk and a pseudoisothermal halo are given respectively as (1)and (2)where r is given in spherical coordinates. Here Σ_{0} is the central surface density of the disk and R_{D} is the disk scalelength; ρ_{0} is the central density of the halo and R_{C} is the halo core radius.
Solving the Poisson equation, the potential in the midplane at a point P due to the disk is given by (see, e.g., Binney & Tremaine 1987) (3)where y = R/ 2R_{D}; and I_{n} and K_{n} (n = 0, 1) are modified Bessel functions of the first and second kind, respectively.
The halo midplane potential at the point P when the halo is concentric with respect to the disk is obtained as (4)The halo midplane potential at the point P when the halo is offcenter is given by (5)We consider small values of d (≪ R), so that (6)Using this expression, the first two terms in Eq. (5) are simplified separately as Therefore, the first terms from both expressions are simply terms from a halo that is concentric with the disk, and the other two terms arise from the nonzero value of d. Thus, the total halo potential, Φ_{H} is given by (9)where the second term on the righthand side defines the perturbation term Φ_{pert}. Thus, the total potential can be written as (10)where the first two terms give the potential due to disk and halo with the same center, C_{D}, and the last term is the perturbation potential due to the halo being offset with respect to the disk center.
2.2. Disk response to lopsided perturbation potential
The perturbation potential (see Eq. (9)) is thus given by (11)where (12)The perturbation potential depends linearly on cos φ, thus it has a lopsided form corresponding to an m = 1 azimuthal perturbation (e.g., Jog & Combes 2009). This is as expected from the geometry of the mass distribution (see the schematic plot in Fig. 1). Its magnitude is linearly proportional to the separation d between the disk and the halo centers. Thus, as the offset d → 0, the perturbation goes to zero and the net potential is then simply the sum of the potentials due to the disk and a concentric halo.
Next we define ϵ_{pert}, the dimensionless perturbation potential, to be the ratio of the perturbation potential over the total undisturbed contributions of concentric disk and halo: (13)Using Ψ(R) as the perturbation potential, and following the procedure for linear perturbation calculations as in Jog (1997, 2000), we next derive the perturbed radius and the azimuthal and radial velocity components for a particle in a closed circular orbit of radius R_{0}.
The perturbation in R is thus obtained as (14)where κ and Ω_{0} denote the epicyclic frequency and the angular speed, respectively, at R_{0}. Thus, (15)The perturbation in the azimuthal velocity v_{φ} is governed by (16)Integrating with respect to time, this gives (17)The net velocity along the azimuthal direction is given by (18)Here v_{c} is the unperturbed rotational velocity. The last two terms give δv_{φ}, the perturbed velocity along the azimuthal direction, and using Eq. (17) this is (19)In full this is given as (20)The perturbed radial velocity is given by (21)Using Eq. (15), this is given as (22)While discussing the results in Sect. 3 we refer to the radius as R for convenience (without the subscript 0).
3. Results
3.1. Input parameters
The above results are general, and are valid for a typical spiral galaxy where the disk dominates the dynamics in the central region. We now consider the input parameters for the Galaxy and study the effect of the offcenter halo on the central disk region. We take the central disk surface density, Σ_{0} = 640.9M_{⊙} pc^{2}, and the exponential disk scalelength, R_{D} = 3.2 kpc; and the halo central density, ρ_{0} = 0.035 M_{⊙} pc^{3}, and the core radius, R_{c} = 5 kpc; these values are taken from the Galaxy mass model by Mera et al. (1998).
3.2. Resulting disk lopsidedness
First, we plot the value of the dimensionless perturbation potential, ϵ_{pert}(R), defined in Eq. (13) for a choice of d = 350 pc (see Fig. 2).
Fig. 2 Perturbation potential ϵ_{pert} arising due to the halo that is offcenter by 350 pc with respect to the disk versus the disk radius R. 
We check that ϵ_{pert}(R) is indeed a perturbation, justifying the assumption made in Sect. 2.1. Its value increases as R decreases, which is expected physically, since the perturbative effect of the halo would be stronger at a smaller radius. Thus, the magnitude of the perturbation potential falls with increasing radius R. By R = 1 kpc, the fractional value of the perturbation potential is less than 1%. Yet, as we show below, this still results in strong spatial and kinematical lopsidedness of the closed orbits in the disk, because the perturbed equations have a strong dependence on radius.
Fig. 3 Resulting orbit for an initial circular orbit at R = 1.5 kpc, shown as R vs. φ in panel a), and as a polar plot in panel b). The net azimuthal and radial velocity components vs. φ are shown in panels c) and d), respectively. The dashed line in panel c) denotes v_{c}. Panels a) and b) show that the orbit is strongly distorted or lopsided spatially: it is shortened along φ = 0° and elongated in the opposite direction, i.e., along φ = 180°. 
In Fig. 3a we plot the net radius of the orbit vs. φ for an initial circular orbit of radius R = 1.5 kpc. The values for the other half of the azimuthal range are not shown as these would give a mirror image of this plot, since these quantities depend simply on cos φ. Figure 3b shows a polar plot of the perturbed orbit drawn against unperturbed, original circular orbits at R_{0} = 1, 2 kpc. Thus the perturbed orbit shows a typical lopsided form: it is shortened along φ = 0° by ~ 40% and elongated along the opposite direction (i.e., along φ = 180°). The net azimuthal and radial velocity components vs. φ are shown in Figs. 3c and d, respectively. We note that the change in azimuthal velocity is huge: its maximum magnitude (which occurs at φ = 0°) is 72.2 km s^{1}. The resulting radial or streaming velocity value, which occurs at φ = 90°, is also significant =35.2 km s^{1}. Thus, the radius and the velocity components of the perturbed orbit are strongly lopsided due to the perturbation potential resulting from the offcenter halo.
The resulting strong lopsidedness is further illustrated in Fig. 4 where the azimuthal velocity due to concentric disk and halo, as well as as the perturbation value, all obtained at φ = 0° are shown as a function of R. The net unperturbed velocity is given by the square root of the disk and halo contributions added in quadrature. The perturbation in azimuthal velocity has a negative sign with respect to the rotation, and its magnitude increases rapidly below 1.5 kpc, and can be comparable to the original rotational velocity. These features indicate that the perturbed orbit is no longer a nearly circular orbit, rather it takes on the shape of a radial orbit.
Fig. 4 Azimuthal velocity components vs. the radius R for the diskalone and the halo, both with the same center; and the perturbed azimuthal velocity component, along φ = 0°, arising when the halo is taken to be offset with respect to the disk center by 350 pc. We note that the disk dominates the dynamics in the central few kpc, and the perturbed velocity component is large. 
The fractional values of the change in the radius and the azimuthal velocity for the perturbed motion with respect to the undisturbed values are shown in Figs. 5a and b, respectively. Clearly a significant radial range up to R = 3 kpc is strongly disturbed, as can be seen from the spatial and kinematical lopsidedness. In fact, the assumption of linear perturbation in the calculations in Sect. 2 is not strictly valid below about R< 1.2 kpc. Hence, in Figs. 4 and 5 we have only shown results beyond this radius. We note that at larger radii (beyond 3.2 kpc), the sign of the radial change is opposite so that it is positive at φ = 0°; however, the magnitude of the fractional change is very small (see Fig. 5a).
Fig. 5 Fractional change in the radius of the orbit (top panel) and fractional change in the azimuthal velocity component (bottom panel) both obtained at φ = 0° vs. the radius R. 
We next check the dependence of the results on the input parameters. First we vary the halo concentration by taking R_{c} = 4 kpc and 6 kpc (instead of 5 kpc as was used above); these correspond respectively to a more centrally concentrated and low mass halo, and a more extended and massive halo. As expected, a lower core radius leads to a lower magnitude of the perturbation potential, ϵ_{pert}. The resulting fractional variation in the orbital radius at R = 1.5 kpc is now 20% for R_{c} = 4 kpc as opposed to 40% obtained for R_{c} = 5 kpc, while a higher R_{c} = 6 kpc results in a higher relative fractional change in radius of 62%. Since the perturbation potential as well as the resulting radial and velocity changes are linearly proportional to the offset value d, a value of d smaller than 350 pc that was used above will lead to correspondingly smaller perturbed values.
As shown above, an offcenter dark matter halo leads to strong central disk lopsidedness. This is in contrast to the fairly moderate lopsidedness seen in the outer disks of galaxies (Rix & Zaritsky 1995; Jog 1997). The outer disks of galaxies are often observed to be lopsided where a radially constant perturbation potential gives a good fit to the observations (Rix & Zaritsky 1995; Angiras et al. 2007; Zaritsky et al. 2013). In that case, at φ = 0°, the typical fractional change in the radius is equal to twice the fractional perturbation potential (i.e., about 10%; Jog 2000); and the fractional change in the azimuthal velocity is equal to the fractional perturbation potential (i.e., about 5%; Jog 2002). In contrast, for the offcenter halo considered here, the resulting spatial and velocity perturbations have much higher values and show a strong radial dependence. This occurs because here the perturbed equations have a strong dependence on radius and do not give a simple linear dependence on the dimensionless perturbation potential ϵ_{pert}.
4. Discussion
We next discuss a few issues related to this work.

(1).
We have treated a kinematic response of the disk to the offset halo for simplicity. A more realistic model would involve taking account of the backreaction of the induced disk lopsidedness on the halo. We defer such a complex, selfconsistent calculation to a future paper.

(2).
We have shown that in a typical spiral galaxy like our Galaxy, the dark matter halo that is offset by 350 pc with respect to the galactic disk center disturbs the orbits and kinematics substantially within the central few kpc of the disk. Since the offset is known to be longlived from simulation studies (Kuhlen et al. 2013), the resulting disk lopsidedness will also be longlived. This can have farreaching implications for the evolution of the central region. For example, the resulting lopsided orbits may be important in feeding the active galactic nucleus (AGN) at the center. The highly disturbed, lopsided orbits in the central region could lead to an outward angular momentum transport, thus helping gas infall, similar to what was shown for lopsidedness in the outer disk (Saha & Jog 2014), and in agreement with the infall proposed for the central region (Combes 2016). We note that this process could be important in fueling a central AGN even in an isolated galaxy, which is a longstanding problem.

(3).
We have assumed the initial orbit to be circular. However, when the random velocity dispersion is taken into account, the particle executes an epicyclic motion around this radius, hence the true undisturbed motion is not circular, and also not closed (since κ/ Ω is in general an irrational number). Incidentally, this situation also occurs in the case of lopsidedness in the outer disk (Rix & Zaritsky 1995; Jog 1997). The net disturbed motion then can be thought of as the circular motion, plus the superposition of the epicyclic motion and the perturbed motion due to the lopsided potential. The random velocity dispersion is known to be high in the galactic center region. The observations by Lewis & Freeman (1989) show that the stellar radial velocity in the Galactic disk has a central value of 95 km s^{1} and falls exponentially with radius with a scalelength of 8.7 kpc. This gives the disk velocity dispersion of 80 km s^{1} at R = 1.5 kpc. This is comparable to the rotational velocity at this radius of 91 km s^{1}, obtained by adding the disk and halo contributions in quadrature, taken from Fig. 4. This figure also gives κ = 111 km s^{1} kpc^{1}. Thus at this radius, the radial amplitude of the lopsided motion is 0.6 kpc (Fig. 5a), while the epicyclic amplitude is =80 km s^{1}/111 km s^{1} kpc^{1}(~0.7 kpc), so the two are comparable. A similar case occurs in the outer galactic disk where the two perturbations lead to comparable radial amplitudes, and yet the epicyclic amplitude is ignored for the sake of simplicity (Rix & Zaritsky 1995; Jog 1997). A similar simplification is made in the treatment of a stellar disk perturbed by the triaxial halo (Bailin et al. 2007). In the present paper also, we only take into account the lopsided perturbed motion with respect to the circular motion. Thus, our model is simple and idealized but it does still bring out the result of strong disk lopsidedness arising due to the dark matter halo offset.

(4).
Generally, once an m = 1 feature develops, it is longlived, lasting for several Gyr, as has been shown for the lopsidedness in the outer galactic disk (Ideta 2002; Saha et al. 2007; also see Jog & Combes 2009). The selfgravity of the disk reduces the differential precession, and that leads to this long life. A similar process would operate in the central region, thus the m = 1 configuration in the disk arising due to halo offset could be longlived. This issue needs to be addressed by future Nbody simulations.
5. Conclusion
We study the effect of an offcenter dark matter halo on the central galactic disk, where the disk dominates the dynamics. To do this, we develop a simple, analytical model where the equations of motion are formulated in terms of the potential contributions of the disk and halo when these two are concentric, plus a perturbation term that arises due to the halo offset. We show that the perturbation potential has a lopsided form, i.e., it shows an m = 1 azimuthal dependence. The halo offset thus acts as a lopsided perturbation potential, and we study the disk response to this. The magnitude of the perturbation potential depends on the halo parameters and also the offset distance d between the disk and the halo centers. Though the fractional perturbation potential is small (<1%), we find that it has a surprisingly strong effect on the orbits and the kinematics in the central region within a radius of about 3 kpc. The resulting strong disk lopsidedness can have important dynamical consequences, for example it can help fuel the central AGN via an outward angular momentum transfer. An offcenter dark matter halo thus results in complex and interesting dynamics in the central region that should be explored in future work.
Acknowledgments
We would like to thank the anonymous referee for thoughtful comments which have greatly improved the paper. C.J. would like to thank the DST, Government of India for support via a J.C. Bose fellowship (SB/S2/JCB31/2014).
References
 Angiras, R. A., Jog, C. J., Dwarakanath, K. S., & Verheijen, M. A. W. 2007, MNRAS, 378, 276 [NASA ADS] [CrossRef] [Google Scholar]
 Bailin, J., Simon, J. D., Bolatto, A. D., Gibson, B. K., & Power, C. 2007, ApJ, 667, 191 [Google Scholar]
 Binney, J., & Tremaine, S. 1987, Galactic Dynamics (Princeton: Princeton University Press) [Google Scholar]
 Chemin, L., Hure, J.M., Soubiran, C., et al. 2016, A&A, 588, A48 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Combes, F. 2016, in IAU Symp. 322, eds. S. Longmore, G. Bicknell, & R. Crocker, 245 [Google Scholar]
 Freeman, K. C. 1970, ApJ, 160, 811 [NASA ADS] [CrossRef] [Google Scholar]
 Ideta, M. 2002, ApJ, 568, 190 [NASA ADS] [CrossRef] [Google Scholar]
 Jog, C. J. 1997, ApJ, 488, 642 [NASA ADS] [CrossRef] [Google Scholar]
 Jog, C. J. 2000, ApJ, 542, 216 [NASA ADS] [CrossRef] [Google Scholar]
 Jog, C. J. 2002, A&A, 391, 471 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Jog, C. J., & Combes 2009, Phys. Rep., 471, 75 [NASA ADS] [CrossRef] [Google Scholar]
 Jog, C. J., & Maybhate, A. 2006, MNRAS, 370, 891 [NASA ADS] [Google Scholar]
 Klypin, A., Zhao, H., & Somerville, R. S. 2002, ApJ, 573, 597 [Google Scholar]
 Kuhlen, M., Guedes, J., Pillepich, A., Madau, P., & Mayer, L. 2013, ApJ, 765, 10 [NASA ADS] [CrossRef] [Google Scholar]
 Lewis, J. R., & Freeman, K. C. 1989, AJ, 97, 139 [NASA ADS] [CrossRef] [Google Scholar]
 Massey, R., Williams, L., Smit, R., et al. 2015, MNRAS, 449, 3393 [Google Scholar]
 Mera, D., Chabrier, G., & Schaeffer, R. 1998, A&A, 330, 953 [NASA ADS] [Google Scholar]
 Miller, R. H., & Smith, B. F. 1992, ApJ, 393, 508 [NASA ADS] [CrossRef] [Google Scholar]
 Rix, H.W., & Zaritsky, D. 1995, ApJ, 447, 82 [NASA ADS] [CrossRef] [Google Scholar]
 Rubin, V. 1983, Science, 220, 1339 [NASA ADS] [CrossRef] [Google Scholar]
 Sackett, P. 1997, ApJ, 483, 103 [NASA ADS] [CrossRef] [Google Scholar]
 Saha, K., & Jog, C. J. 2014, MNRAS, 444, 352 [NASA ADS] [CrossRef] [Google Scholar]
 Saha, K., Combes, F., & Jog, C. J. 2007, MNRAS, 382, 419 [NASA ADS] [CrossRef] [Google Scholar]
 Schaller, M., Robertson, A., Massey, R., Bower, R. G., & Eke, V. R. 2015, MNRAS, 453, L58 [NASA ADS] [CrossRef] [Google Scholar]
 Zaritsky, D., Salo, H., & Laurikainen, E., et al. 2013, ApJ, 772, 135 [NASA ADS] [CrossRef] [Google Scholar]
All Figures
Fig. 1 Schematic diagram showing the centers of disk and an offcenter dark matter halo, separated by a distance d, and their effect at point P in the disk plane. 

In the text 
Fig. 2 Perturbation potential ϵ_{pert} arising due to the halo that is offcenter by 350 pc with respect to the disk versus the disk radius R. 

In the text 
Fig. 3 Resulting orbit for an initial circular orbit at R = 1.5 kpc, shown as R vs. φ in panel a), and as a polar plot in panel b). The net azimuthal and radial velocity components vs. φ are shown in panels c) and d), respectively. The dashed line in panel c) denotes v_{c}. Panels a) and b) show that the orbit is strongly distorted or lopsided spatially: it is shortened along φ = 0° and elongated in the opposite direction, i.e., along φ = 180°. 

In the text 
Fig. 4 Azimuthal velocity components vs. the radius R for the diskalone and the halo, both with the same center; and the perturbed azimuthal velocity component, along φ = 0°, arising when the halo is taken to be offset with respect to the disk center by 350 pc. We note that the disk dominates the dynamics in the central few kpc, and the perturbed velocity component is large. 

In the text 
Fig. 5 Fractional change in the radius of the orbit (top panel) and fractional change in the azimuthal velocity component (bottom panel) both obtained at φ = 0° vs. the radius R. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.